Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Longitude”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 127 records · Page 7

Multiscale ACI Satellite Database

The SATELLITE_EAGLES_PNNL NetCDF dataset contains a suite of satellite- and reanalysis-derived atmospheric and surface parameters on a regular latitude–longitude grid. The dataset includes core geophysical fields such as land fraction, aerosol optical depth at multiple wavelengths (465, 550, 667, and 865 nm), sea surface temperature, estimated inversion strength, and various thermodynamic and dynamic quantities (e.g., relative humidity, vertical velocity, boundary-layer height, and surface fluxes) from both MERRA and ERA reanalysis products, provided as daily-mean and instantaneous values. A major component of the dataset consists of MODIS-retrieved cloud microphysical properties, including cloud droplet number concentration, cloud effective radius, optical thickness, and liquid water path, provided for three compositing regimes (“All,” “Q06,” and “G18”). Corresponding cloud-top parameters—temperature, height, and pressure—along with total and domain-mean cloud fraction fields are also included. The file further integrates additional satellite data from AMSR-E (for cloud water, rain water, and surface precipitation retrievals) and CERES (for top-of-atmosphere radiative fluxes, cloud fractions, and albedo). This dataset is designed to evaluate aerosol–cloud interactions in warm clouds, emphasizing the use of MODIS for deriving cloud droplet number concentration and liquid water path statistics. The complementary satellite and reanalysis fields are co-located and time-matched to the same instantaneous MODIS observations, enabling consistent comparisons between cloud properties, aerosol loading, and large-scale meteorological conditions. The dataset is recently featured in Christensen et al. (2025), Machine Learning Reveals Strong Grid-Scale Dependence in the Satellite Nd–LWP Relationship, EGUsphere [preprint], https://doi.org/10.5194/egusphere-2025-3850, 2025.

Christensen, Matthew [Pacific Northwest National L↗

Multiscale ACI Satellite Database

The SATELLITE_EAGLES_PNNL NetCDF dataset contains a suite of satellite- and reanalysis-derived atmospheric and surface parameters on a regular latitude–longitude grid. The dataset includes core geophysical fields such as land fraction, aerosol optical depth at multiple wavelengths (465, 550, 667, and 865 nm), sea surface temperature, estimated inversion strength, and various thermodynamic and dynamic quantities (e.g., relative humidity, vertical velocity, boundary-layer height, and surface fluxes) from both MERRA and ERA reanalysis products, provided as daily-mean and instantaneous values. A major component of the dataset consists of MODIS-retrieved cloud microphysical properties, including cloud droplet number concentration, cloud effective radius, optical thickness, and liquid water path, provided for three compositing regimes (“All,” “Q06,” and “G18”). Corresponding cloud-top parameters—temperature, height, and pressure—along with total and domain-mean cloud fraction fields are also included. The file further integrates additional satellite data from AMSR-E (for cloud water, rain water, and surface precipitation retrievals) and CERES (for top-of-atmosphere radiative fluxes, cloud fractions, and albedo). This dataset is designed to evaluate aerosol–cloud interactions in warm clouds, emphasizing the use of MODIS for deriving cloud droplet number concentration and liquid water path statistics. The complementary satellite and reanalysis fields are co-located and time-matched to the same instantaneous MODIS observations, enabling consistent comparisons between cloud properties, aerosol loading, and large-scale meteorological conditions. The dataset is recently featured in Christensen et al. (2025), Machine Learning Reveals Strong Grid-Scale Dependence in the Satellite Nd–LWP Relationship, EGUsphere [preprint], https://doi.org/10.5194/egusphere-2025-3850, 2025.

54 ENVIRONMENTAL SCIENCES↗

Classification Analysis of Southwest Pacific Tropical Cyclone Intensity Changes Prior to Landfall

This study evaluates the ability of a random forest classifier to identify tropical cyclone (TC) intensification or weakening prior to landfall over the western region of the Southwest Pacific Ocean (SWPO) basin. For both Australia mainland and SWPO island cases, when a TC first crosses land after spending ≥24 h over the ocean, the closest hour prior to the intersection is considered as the landfall hour. If the maximum wind speed (V max ) at the landfall hour increased or remained the same from the 24-h mark prior to landfall, the TC is labeled as intensifying and if the V max at the landfall hour decreases, the TC is labeled as weakening. Geophysical and aerosol variables closest to the 24 h before landfall hour were collected for each sample. The random forest model with leave-one-out cross validation and the random oversampling example technique was identified as the best-performing classifier for both mainland and island cases. The model identified longitude, initial intensity, and sea skin temperature as the most important variables for the mainland and island landfall classification decisions. Incorrectly classified cases from the test data were analyzed by sorting the cases by their initial intensity hour, landfall hour, monthly distribution, and 24-h intensity changes. TC intensity changes near land strongly impact coastal preparations such as wind damage and flood damage mitigations; hence, this study will contribute to improve identifying and prioritizing prediction of important variables contributing to TC intensity change before landfall.

54 ENVIRONMENTAL SCIENCES↗

Spatial Interpolation of Air Pollutant and Meteorological Variables in Central Amazonia

The Amazon Rainforest is highlighted by the global community both for its extensive vegetation cover that constantly suffers the effects of anthropic action and for its substantial biodiversity. This dataset presents data of meteorological variables from the Amazon Rainforest region with a spatial resolution of 0.001° in latitude and longitude, resulting from an interpolation process. The original data were obtained from the GoAmazon 2014/5 project, in the Atmospheric Radiation Measurement (ARM) repository, and then processed through mathematical and statistical methods. The dataset presented here can be used in experiments in the field of Data Science, such as training models for predicting climate variables or modeling the distribution of species.

54 ENVIRONMENTAL SCIENCES↗

AATTENUATION—The Atmospheric Attenuation Model for CSP Tower Plants: A Look-Up Table for Operational Implementation

Attenuation of solar radiation between the receiver and the heliostat field in concentrated solar power (CSP) tower plants can reduce the overall system performance significantly. The attenuation varies strongly with time and the average attenuation at different sites might also vary strongly from each other. If no site specific attenuation data is available, the optimal plant design cannot be determined and rough estimations of the attenuation effect are required leading to high uncertainties of yield analysis calculations. The attenuation is caused mainly by water vapor content and aerosol particles in the lower atmospheric layer above ground. Although several on-site measurement systems have been developed during recent years, attenuation data sets are usually not available to be included during the plant project development. An Atmospheric Attenuation (AATTENUATION) model to derive the atmospheric transmittance between a heliostat and receiver on the basis of common direct normal irradiance (DNI), temperature, relative humidity, and barometric pressure measurements was developed and validated by the authors earlier. The model allows the accurate estimation of attenuation for sites with low attenuation and gives an estimation of the attenuation for less clear sites. However, the site-dependent coefficients of the AATTENUATION model had to be developed individually for each site of interest, which required time-consuming radiative transfer simulations, considering the exact location and altitude, as well as the pre-dominant aerosol type at the location. This strongly limited the application of the model despite its typically available input data. In this manuscript, a look-up table (LUT) is presented which enables the application of the AATTENUATION model at the site of interest without the necessity to perform the according complex radiative transfer calculations for each site individually. This enables the application of the AATTENUATION model for virtually all resource assessments for tower plants and in an operational mode in real time within plant monitoring systems around the world. The LUT also facilitates the generation of solar attenuation maps on the basis of long-term meteorological data sets which can be considered during resource assessment for CSP tower plant projects. The LUTs are provided together with this manuscript as supplementary files. The LUT for the AATTENUATION model was developed for a solar zenith angle (SZA) grid of 1°, an altitude grid of 100 m, 7 different standard aerosol types and the standard AFGL atmospheres for mid-latitudes and the tropics. The LUT was tested against the original version of the AATTENUATION model at 4 sites in Morocco and Spain, and it was found that the additional uncertainty introduced by the application of the LUT is negligible. With the information of latitude, longitude, altitude above mean sea level, DNI, relative humidity (RH), ambient temperature (Tair), and barometric pressure (bp), the attenuation can be now derived easily for each site of interest.

14 SOLAR ENERGY↗

Open Cluster IC 1369 and Its Vicinity: Multicolor Photometry and Gaia DR2 Astrometry

The open cluster IC 1369, located in Cygnus close to the galactic equator at the 89.°6 longitude, and its vicinity are investigated by applying two-dimensional photometric classification of stars down to V = 19 mag measured in the Vilnius seven-color photometric system. Employing the Gaia DR2 coordinates, parallaxes, and proper motions, 143 cluster members are found having cluster membership probabilities ≥0.75. The hottest stars of the evolved main sequence are of spectral class B9. Distances to the cluster, determined by two methods based on the Gaia parallaxes, are 3280 and 3370 pc with an uncertainty of 400–500 pc. Employing the physical HR diagram with the coordinates logL/L{sub ⊙} versus logT{sub eff} and the Padova isochrones, the cluster age is found to be between 300 and 350 Myr. Among the cluster members we identified six red giants of spectral classes G5–G8 and luminosities II–III, 17 Am stars, one Ap star, and one B5 blue straggler. The extinction shows a steep rise up to ∼2.0 mag at 300–700 pc, caused by the extensions of the dust cloud LDN 970 and the Great Cygnus Rift, and up to ∼2.5 mag at the edge of the Perseus arm at 3.0 kpc. The extinction is almost constant between 0.9 and 2.5 kpc from the Sun. The average extinction of the cluster members is 2.54 mag. A systematic error of color excesses found in the Argonaut calculator results is discussed.

79 ASTRONOMY AND ASTROPHYSICS↗

Further Investigation on Chromospheric and Prominence Activity of the RS Canum Venaticorum Star SZ Piscium

To continue our study on chromospheric activity and detection for possible prominence events of the very active RS Canum Venaticorum star SZ Piscium (SZ Psc), long-term high-resolution spectroscopic observations were obtained during several observing runs from 2014 to 2018. Based on the spectral subtraction technique, the chromospheric emission of the Ca ii IRT (λ8662, λ8542, and λ8498), H{sub α}, Na i D{sub 1}, D{sub 2} doublet, H{sub β}, and Ca ii H & K lines is mainly associated with the K1 IV primary star of the SZ Psc system, in good agreement with the previous studies, and the F8 V secondary star also shows some chromospheric emission, implying its active chromosphere. Moreover, an optical flare characterized by the He i D{sub 3} line emission together with stronger emission in the other indicators was detected. Furthermore, two chromospheric active longitudes around the two quadratures of the system were identified for most of the time, and the chromospheric activity shows significant changes during a few orbital cycles. The chromospheric activity level seems to show a long-term variation during our observations. There were some excess absorption features in the subtracted H{sub α} line and the other activity indicators, which would be caused by prominence-like materials associated with the K1 IV primary star of the system. Prominence materials could absorb the chromospheric emission and continuum from the K1 IV primary star and even the F8 V secondary one.

79 ASTRONOMY AND ASTROPHYSICS↗

A Census of Star Formation in the Outer Galaxy. II. The GLIMPSE360 Field

We have conducted a study of star formation in the outer Galaxy from 65°< l < 265° in the region observed by the GLIMPSE360 program. This Spitzer warm mission program mapped the plane of the outer Milky Way with IRAC at 3.6 and 4.5 μm. We combine the IRAC, Wide-field Infrared Survey Explorer (WISE), and Two Micron All Sky Survey catalogs and our previous results from another outer Galaxy survey and identify a total of 47,338 young stellar objects (YSOs) across the field spanning >180° in Galactic longitude. Using the DBSCAN method on the combined catalog, we identify 618 clusters or aggregations of YSOs having five or more members. We identify 10,476 class I, 29,604 class II, and 7325 anemic class II/class III YSOs. The ratio of YSOs identified as members of clusters was 25,528/47,338, or 54%. We found that 100 of the clusters identified have previously measured distances in the WISE H ii survey. We used these distances in our spectral energy distribution (SED) fitting of the YSOs in these clusters, of which 96 had YSOs with <3σ fits. We used the derived masses from the SED model fits to estimate the initial mass function (IMF) in the inner and outer Galaxy clusters; dividing the clusters by galactocentric distances, the slopes were Γ = 1.87 ± 0.31 above 3 M {sub ⊙} for R {sub Gal} < 11.5 kpc and Γ = 1.15 ± 0.24 above 3 M {sub ⊙} for R {sub Gal} > 11.5 kpc. The slope of the combined IMF was found to be Γ = 1.92 ± 0.42 above 3 M {sub ⊙}. These values are consistent with each other within the uncertainties and with literature values in the inner Galaxy high-mass star formation regions. The slopes are likely also consistent with a universal Salpeter IMF.

79 ASTRONOMY AND ASTROPHYSICS↗

An Investigation of Poorly Studied Open Cluster NGC 4337 Using Multicolor Photometric and Gaia DR2 Astrometric Data

We present a comprehensive analysis (photometric and kinematical) of the poorly studied open cluster NGC 4337 using 2MASS, WISE, APASS, and Gaia DR2 databases. By determining the membership probabilities of stars, we identified the 624 most probable members with membership probability higher than 50% by using proper motion and parallax data taken from Gaia DR2. The mean proper motion of the cluster is obtained as μ{sub x}=−8.83±0.01 and μ {sub y} = 1.49 ± 0.006 mas yr{sup −1}. We find the normal interstellar extinction toward the cluster region. The radial distribution of members provides a cluster radius of 7.′75 (5.63 pc). The estimated age of 1600 ± 180 Myr indicates that NGC 4337 is an old open cluster with a bunch of red giant stars. The overall mass function slope for main-sequence stars is found as 1.46±0.18 within the mass range 0.75−2.0 M{sub ⊙}, which is in fair agreement with Salpeter’s value (x = 1.35) within uncertainty. The present study demonstrates that NGC 4337 is a dynamically relaxed open cluster. Using the Galactic potential model, Galactic orbits are obtained for NGC 4337. We found that this object follows a circular path around the Galactic center. Under the kinematical analysis, we compute the apex coordinates (A, D) by using two methods: (i) the classical convergent point method and (ii) the AD-diagram method. The obtained coordinates are (A {sub conv}, D {sub conv}) = (96.°27 ± 0.°10, 13.°14 ± 0.°27) and (A {sub ◦}, D {sub ◦}) = (100.°282 ± 0.°10, 9.°577 ± 0.°323) respectively. We also computed the Velocity Ellipsoid Parameters, matrix elements (μ {sub ij}), direction cosines (l {sub j}, m {sub j}, n {sub j}), and the Galactic longitude of the vertex (l {sub 2}).

79 ASTRONOMY AND ASTROPHYSICS↗

BVR{sub C}I{sub C} Observations and Analysis of the UV LEO-type Binary, the Pre-WUMa, V1023 Persei

V1023 Per is a polar spotted, but well-detached dwarf binary, very likely a Pre-WUMa eclipsing binary. It was observed on 11 nights in 2015 October, November, and December at Dark Sky Observatory in North Carolina with the 0.81 m reflector of Appalachian State University. The period behavior is complex and may be increasing with a cubic or quadratic ephemeris. Its odd light curves are of high amplitude but have no totality and reach maximum light just before and following the secondary eclipse, indicating that it has polar spots similar to UV Leo. The 28 yr orbital period study calls the attention of observers to further monitor this binary to understand the complex nature of the period evolution. Its presently fixed polar spot does indicate that it must have a strong magnetic field and that it is synchronously rotating. The BVR{sub c}I{sub c} simultaneous 2016 Wilson–Devinney program (W-D) solution gives fillouts of 68% and 75% for the primary and secondary components, respectively. The polar position of the spot, its radius (24°), and T-fact (0.75) indicate that a strong magnetic field is present. The small ΔT in the components (∼289 K) show that the stars are similar in spectral type despite them being well detached. The inclination is high, ∼85.°19 ± 0.°04, but there is not a total eclipse due to the high mass ratio (∼0.739 ± 0.001). Due to the low galactic position (longitude = 150.°390, latitude = −1.°023), its reddening is addressed. The primary components temperature is estimated to be 5250 ± 250 K.

79 ASTRONOMY AND ASTROPHYSICS↗

Signatures of a Tidally Induced Spiral Arm at the Anticenter of the Milky Way and a Kinematically Extended Anticenter Stream Using DESI Data Release 2

Using the Dark Energy Spectroscopic Instrument (DESI) Milky Way Survey, we examine the six-dimensional space of the anticenter region of the Milky Way stellar disk (150° < Galactic longitude < 220°) using 61,883 main-sequence turnoff stars. We focus on two well-known stellar overdensities in the anticenter: the Monoceros Ring (MRi) and Anticenter Stream (ACS). We find that the MRi overdensity has kinematic signatures consistent with a tidally induced spiral arm, a type of dynamic spiral arm created by an interaction with a satellite galaxy, most likely the Sagittarius dwarf spheroidal galaxy (Sgr). We use the kinematics of the MRi to calculate the two most recent passage times of Sgr, finding 0.25 ± 0.09 Gyr and 1.10 ± 0.23 Gyr from the present day. We validate that the ACS is kinematically decoupled from the MRi because they are moving in opposite radial and vertical directions. We find that the kinematics associated with the ACS extends beyond our defined overdensity. The features we see in the ACS region are likely part of a broader distribution of stars with the same kinematic signature as detected in other places, like the vertical wave in the outer disk and phase spiral.

Lambert, Mika [University of California, Santa Cru↗

Fermi-LAT Detection of the Supernova Remnant G312.4-0.4 in the Vicinity of 4FGL J1409.1-6121e

Abstract Gamma-ray emission provides constraints on the nonthermal radiation processes at play in astrophysical particle accelerators. This allows both the nature of accelerated particles and the maximum energy that they can reach to be determined. Notably, it remains an open question to what extent supernova remnants contribute to the sea of Galactic cosmic rays. In the Galactic plane, at around 312° of Galactic longitude, Fermi-LAT observations show an extended source (4FGL J1409.1−6121e) around five powerful pulsars. This source is described by one large disk of 0.°7 radius with a high significance of 45σin the 4FGL-DR3 catalog. Using 14 yr of Fermi-LAT observations, we revisited this region with a detailed spectro-morphological analysis in order to disentangle its underlying structure. Three sources have been distinguished, including the supernova remnant G312.4−0.4 whose gamma-ray emission correlates well with the shell observed at radio energies. The hard spectrum detected by the LAT, extending up to 100 GeV without any sign of a cutoff, is well reproduced by a purely hadronic model.

Astronomy & Astrophysics↗

A Data-driven, Physics-based Transport Model of Solar Energetic Particles Accelerated by Coronal Mass Ejection Shocks Propagating through the Solar Coronal and Heliospheric Magnetic Fields

Abstract In an effort to develop computational tools for predicting radiation hazards from solar energetic particles (SEPs), we have created a data-driven physics-based particle transport model to calculate the injection, acceleration, and propagation of SEPs from coronal mass ejection (CME) shocks traversing through the solar corona and interplanetary magnetic fields. The model runs on an input of corona and heliospheric plasma and magnetic field configuration from a magnetohydrodynamic model driven by solar photospheric magnetic field measurements superposed with observed CME shocks determined from coronagraph images. SEP source particles are injected at the shock using the result of diffusive shock acceleration formulation from a characteristic obliquity-dependent injection from a heated solar wind thermal tail population. With several advanced computation techniques involving stochastic simulation and integration, the model obtains the particle intensity at any location in interplanetary space through the rigorous solution to the time-dependent 5D focus transport equation in the phase space that includes perpendicular diffusion. We apply the model to the 2011 November 3 CME event. The calculation results reproduce multispacecraft SEP observations at Earth and STEREO-B reasonably well without normalization of particle flux. The observations at STEREO-A can be reproduced by rescaling particle energy or modified energy dependence of particle diffusion coefficients. This circumsolar SEP event seen by spacecraft at Earth, STEREO-A, and STEREO-B at widely separated longitudes can be explained by diffusive shock acceleration by a single CME shock with a moderate speed.

79 ASTRONOMY AND ASTROPHYSICS↗

Global ocean dimethyl sulfide climatology estimated from observations and an artificial neural network

Marine dimethyl sulfide (DMS) is important to climate due to the ability of DMS to alter Earth's radiation budget. Knowledge of the global-scale distribution, seasonal variability, and sea-to-air flux of DMS is needed in order to improve understanding of atmospheric sulfur, aerosol/cloud dynamics, and albedo. Here we examine the use of an artificial neural network (ANN) to extrapolate available DMS measurements to the global ocean and produce a global climatology with monthly temporal resolution. A global database of 82 996 ship-based DMS measurements in surface waters was used along with a suite of environmental parameters consisting of latitude–longitude coordinates, time of day, time of year, solar radiation, mixed layer depth, sea surface temperature, salinity, nitrate, phosphate, and silicate. Linear regressions of DMS against the environmental parameters show that on a global-scale mixed layer depth and solar radiation are the strongest predictors of DMS. These parameters capture ~9 % and ~7 % of the raw DMS data variance, respectively. Multilinear regression can capture more of the raw data variance (~39 %) but strongly underestimates DMS in high-concentration regions. In contrast, the artificial neural network captures ~66 % of the raw data variance in our database. Like prior climatologies our results show a strong seasonal cycle in surface ocean DMS with the highest concentrations and sea-to-air fluxes in the high-latitude summertime oceans. We estimate a lower global sea-to-air DMS flux (20.12±0.43 Tg S yr –1 ) than the prior estimate based on a map interpolation method when the same gas transfer velocity parameterization is used. Our sensitivity test results show that DMS concentration does not change unidirectionally with each of the environmental parameters, which emphasizes the interactions among these parameters. The ANN model suggests that the flux of DMS from the ocean to the atmosphere will increase with global warming. Given that larger DMS fluxes induce greater cloud albedo, this corresponds to a negative climate feedback.

54 ENVIRONMENTAL SCIENCES↗

Organic matter transformations are disconnected between surface water and the hyporheic zone

Abstract. Biochemical transformations of organic matter (OM) are a primary driver of river corridor biogeochemistry, thereby modulating ecosystem processes at local to global scales. OM transformations are driven by diverse biotic and abiotic processes, but we lack knowledge of how the diversity of those processes varies across river corridors and across surface and subsurface components of river corridors. To fill this gap we quantified the number of putative biotic and abiotic transformations of organic molecules across diverse river corridors using ultra-high-resolution mass spectrometry. The number of unique transformations is used here as a proxy for the diversity of biochemical processes underlying observed profiles of organic molecules. For this, we use public data spanning the contiguous United States (ConUS) from the Worldwide Hydrobiogeochemical Observation Network for Dynamic River Systems (WHONDRS) consortium. Our results show that surface water OM had more biotic and abiotic transformations than OM from shallow hyporheic zone sediments (1–3 cm depth). We observed substantially more biotic than abiotic transformations, and the numbers of biotic and abiotic transformations were highly correlated with each other. We found no relationship between the number of transformations in surface water and sediments and no meaningful relationships with latitude, longitude, or climate. We also found that the composition of transformations in sediments was not linked with transformation composition in adjacent surface waters. We infer that OM transformations represented in surface water are an integrated signal of diverse processes occurring throughout the upstream catchment. In contrast, OM transformations in sediments likely reflect a narrower range of processes within the sampled volume. This indicates decoupling between the processes influencing surface water and sediment OM, despite the potential for hydrologic exchange to homogenize OM. We infer that the processes influencing OM transformations and the scales at which they operate diverge between surface water and sediments.

58 GEOSCIENCES↗

Global anthropogenic emissions (CAMS-GLOB-ANT) for the Copernicus Atmosphere Monitoring Service simulations of air quality forecasts and reanalyses

Anthropogenic emissions are the result of many different economic sectors, including transportation, power generation, industrial, residential and commercial activities, waste treatment and agricultural practices. Air quality models are used to forecast the atmospheric composition, analyze observations and reconstruct the chemical composition of the atmosphere during the previous decades. In order to drive these models, gridded emissions of all compounds need to be provided. This paper describes a new global inventory of emissions called CAMS-GLOB-ANT, developed as part of the Copernicus Atmosphere Monitoring Service (CAMS; https://doi.org/10.24380/eets-qd81, Soulie et al., 2023). The inventory provides monthly averages of the global emissions of 36 compounds, including the main air pollutants and greenhouse gases, at a spatial resolution of 0.1° × 0.1° in latitude and longitude, for 17 emission sectors. The methodology to generate the emissions for the 2000–2023 period is explained, and the datasets are analyzed and compared with publicly available global and regional inventories for selected world regions. Depending on the species and regions, good agreements as well as significant differences are highlighted, which can be further explained through an analysis of different sectors as shown in the figures in the Supplement.

54 ENVIRONMENTAL SCIENCES↗

The E3SM Diagnostics Package (E3SM Diags v2.7): a Python-based diagnostics package for Earth system model evaluation

Abstract. The E3SM Diagnostics Package (E3SM Diags) is a modern, Python-based Earth system model (ESM) evaluation tool (with Python module name e3sm_diags), developed to support the Department of Energy (DOE) Energy Exascale Earth System Model (E3SM). E3SM Diags provides a wide suite of tools for evaluating native E3SM output, as well as ESM data on regular latitude–longitude grids, including output from Coupled Model Intercomparison Project (CMIP) class models. E3SM Diags is modeled after the National Center for Atmospheric Research (NCAR) Atmosphere Model Working Group (AMWG, 2022) diagnostics package. In its version 1 release, E3SM Diags included a set of core essential diagnostics to evaluate the mean physical climate from model simulations. As of version 2.7, more process-oriented and phenomenon-based evaluation diagnostics have been implemented, such as analysis of the quasi-biennial oscillation (QBO), the El Niño–Southern Oscillation (ENSO), streamflow, the diurnal cycle of precipitation, tropical cyclones, ozone and aerosol properties. An in situ dataset from DOE's Atmospheric Radiation Measurement (ARM) program has been integrated into the package for evaluating the representation of simulated cloud and precipitation processes. This tool is designed with enough flexibility to allow for the addition of new observational datasets and new diagnostic algorithms. Additional features include customizable figures; streamlined installation, configuration and execution; and multiprocessing for fast computation. The package uses an up-to-date observational data repository maintained by its developers, where recent datasets are added to the repository as they become available. Finally, several applications for the E3SM Diags module were introduced to fit a diverse set of use cases from the scientific community.

54 ENVIRONMENTAL SCIENCES↗

Evaluation of native Earth system model output with ESMValTool v2.6.0

Earth system models (ESMs) are state-of-the-art climate models that allow numerical simulations of the past, present-day, and future climate. To extend our understanding of the Earth system and improve climate change projections, the complexity of ESMs heavily increased over the last decades. As a consequence, the amount and volume of data provided by ESMs has increased considerably. Innovative tools for a comprehensive model evaluation and analysis are required to assess the performance of these increasingly complex ESMs against observations or reanalyses. One of these tools is the Earth System Model Evaluation Tool (ESMValTool), a community diagnostic and performance metrics tool for the evaluation of ESMs. Input data for ESMValTool needs to be formatted according to the CMOR (Climate Model Output Rewriter) standard, a process that is usually referred to as “CMORization”. While this is a quasi-standard for large model intercomparison projects like the Coupled Model Intercomparison Project (CMIP), this complicates the application of ESMValTool to non-CMOR-compliant climate model output. In this paper, we describe an extension of ESMValTool introduced in v2.6.0 that allows seamless reading and processing of “native” climate model output, i.e., operational output produced by running the climate model through the standard workflow of the corresponding modeling institute. This is achieved by an extension of ESMValTool's preprocessing pipeline that performs a CMOR-like reformatting of the native model output during runtime. Thus, the rich collection of diagnostics provided by ESMValTool is now fully available for these models. For models that use unstructured grids, a further preprocessing step required to apply many common diagnostics is regridding to a regular latitude–longitude grid. Extensions to ESMValTool's regridding functions described here allow for more flexible interpolation schemes that can be used on unstructured grids. Currently, ESMValTool supports nearest-neighbor, bilinear, and first-order conservative regridding from unstructured grids to regular grids. Example applications of this new native model support are the evaluation of new model setups against predecessor versions, assessing of the performance of different simulations against observations, CMORization of native model data for contributions to model intercomparison projects, and monitoring of running climate model simulations. For the latter, new general-purpose diagnostics have been added to ESMValTool that are able to plot a wide range of variable types. Currently, five climate models are supported: CESM2 (experimental; at the moment, only surface variables are available), EC-Earth3, EMAC, ICON, and IPSL-CM6. As the framework for the CMOR-like reformatting of native model output described here is implemented in a general way, support for other climate models can be easily added.

58 GEOSCIENCES↗