Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Earth Systems”

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 973 records · Page 54

High-Resolution Modeling of the Gulf of Mexico using E3SM

Coastal ocean modeling is a high priority in the DOE‘s Energy Exascale Earth System Model (E3SM). The goal is to accurately predict the risk of damage to coastal resources and infrastructure due to a changing climate in the coming decades. North American coastal communities are areas of particular interest, as this fits under the topic of US national security and planning assessments in a changing climate. LANL Institutional Computing time for the Tier 1 allocation ”Coastal Ocean and Sea Ice Modeling” have been used for development and testing of numerical methods needed for E3SM coastal applications.

54 ENVIRONMENTAL SCIENCES↗

Can Biogeochemical Tracer Observations Constrain Southern Ocean Diapycnal Mixing Rates?

Abstract Direct observations of background diapycnal mixing rates in the Southern Ocean (SO) are limited spatially and temporally, making the choice of an appropriate value to parameterize this mixing in Earth system models a challenge. However, the deployment of Argo floats throughout the SO has provided an extensive range of observations of both physical and biogeochemical parameters. We use an ocean state estimate run with various background diapycnal mixing coefficients to assess if biogeochemical tracer observations can be used to better constrain SO diapycnal mixing rates. We find that vertical tracer distributions in the SO are highly sensitive to the rate of background diapycnal mixing and can provide an upper limit on background mixing rates. This demonstrates the importance of biogeochemical tracer observations throughout the full depth of the water column to validate ocean models.

Ellison, Elizabeth [The Australian National Univer↗

Air‐Ice‐Ocean Coupling During a Strong Mid‐Winter Cyclone: Observing Coupled Dynamic Interactions Across Scales

Abstract Arctic cyclones are key drivers of sea ice and ocean variability. During the 2019–2020 Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) expedition, joint observations of the coupled air‐ice‐ocean system were collected at multiple spatial scales. Here, we present observations of a strong mid‐winter cyclone that impacted the MOSAiC site as it drifted in the central Arctic pack ice. The sea ice dynamical response showed spatial structure at the scale of the evolving and translating cyclonic wind field. Internal ice stress and ocean stress play significant roles, resulting in timing offsets between the atmospheric forcing and the ice response and post‐cyclone inertial ringing in the ice and ocean. Ice motion in response to the wind field then forces the upper ocean currents through frictional drag. The strongest impacts to the sea ice and ocean from the passing cyclone occur as a result of the surface impacts of a strong atmospheric low‐level jet (LLJ) behind the trailing cold front and changing wind directions between the warm‐sector LLJ and post cold‐frontal LLJ. Impacts of the cyclone are prolonged through the coupled ice‐ocean inertial response. Local impacts of the approximately 120 km wide LLJ occur over a 12 hr period or less and at scales of a kilometer to a few tens of kilometers, meaning that these impacts occur at combined smaller spatial scales and faster time scales than most satellite observations and coupled Earth system models can resolve.

Meteorology & Atmospheric Sciences↗

The GFDL Variable-Resolution Global Chemistry-Climate Model for Research at the Nexus of US Climate and Air Quality Extremes

We present a variable-resolution global chemistry-climate model (AM4VR) developed at NOAA's Geophysical Fluid Dynamics Laboratory (GFDL) for research at the nexus of US climate and air quality extremes. AM4VR has a horizontal resolution of 13 km over the US, allowing it to resolve urban-to-rural chemical regimes, mesoscale convective systems, and land-surface heterogeneity. With the resolution gradually reducing to 100 km over the Indian Ocean, we achieve multi-decadal simulations driven by observed sea surface temperatures at 50% of the computational cost for a 25-km uniform-resolution grid. In contrast with GFDL's AM4.1 contributing to the sixth Coupled Model Intercomparison Project at 100 km resolution, AM4VR features much improved US climate mean patterns and variability. In particular, AM4VR shows improved representation of: precipitation seasonal-to-diurnal cycles and extremes, notably reducing the central US dry-and-warm bias; western US snowpack and summer drought, with implications for wildfires; and the North American monsoon, affecting dust storms. AM4VR exhibits excellent representation of winter precipitation, summer drought, and air pollution meteorology in California with complex terrain, enabling skillful prediction of both extreme summer ozone pollution and winter haze events in the Central Valley. AM4VR also provides vast improvements in the process-level representations of biogenic volatile organic compound emissions, interactive dust emissions from land, and removal of air pollutants by terrestrial ecosystems. We highlight the value of increased model resolution in representing climate–air quality interactions through land-biosphere feedbacks. AM4VR offers a novel opportunity to study global dimensions to US air quality, especially the role of Earth system feedbacks in a changing climate.

54 ENVIRONMENTAL SCIENCES↗

Contrasting Carbon–Water–Energy Dynamics in Perennial and Annual Bioenergy Agroecosystems Using Eddy Covariance and Interpretable Machine Learning

Understanding how agroecosystems respond to environmental variability is fundamental to predicting productivity and sustainability under a changing climate. We analyzed 55 site-years of high-frequency eddy covariance observations from five agroecosystems—two perennial grasses (miscanthus and switchgrass), two annual rotation systems (maize–soybean and sorghum–soybean), and a restored native prairie—to examine ecosystem-scale carbon, water, and energy fluxes. Using an interpretable machine-learning framework with regression tree ensembles, Shapley Additive Explanations, and Accumulated Local Effects, we quantified how environmental and temporal factors regulate gross primary productivity (GPP), evapotranspiration (ET), water-use efficiency, and the Bowen ratio. Perennials exhibited stronger physiological buffering and maintained fluxes across a broader range of temperature and moisture conditions, reflecting deeper rooting and persistent canopy cover. Annuals, in contrast, showed greater short-term variability and stronger coupling to atmospheric demand, with GPP and ET declining rapidly under low humidity or soil moisture. Differences in temperature sensitivity of Bowen ratio further revealed that perennials sustained proportionally greater sensible heat flux under cool conditions, whereas annuals exhibited constrained energy exchange when evaporative demand was low. Together, these results demonstrate that crop life cycle and canopy structure are fundamental determinants of ecosystem-scale carbon–water–energy coupling. By integrating long-term flux observations with interpretable machine learning, this study identifies the environmental drivers that shape agroecosystem function and highlights how conversion from annual to perennial feedstocks can enhance climatic resilience and alter land–atmosphere energy feedbacks. These findings provide a data-driven basis for improving crop and Earth-system models and for guiding bioenergy landscape design under future climate scenarios.

Accumulated Local Effects↗

Evapotranspiration depletes groundwater under warming over the contiguous United States

A warmer climate increases evaporative demand. However, response to warming depends on water availability. Existing earth system models represent soil moisture but simplify groundwater connections, a primary control on soil moisture. Here we apply an integrated surface-groundwater hydrologic model to evaluate the sensitivity of shallow groundwater to warming across the majority of the US. We show that as warming shifts the balance between water supply and demand, shallow groundwater storage can buffer plant water stress; but only where shallow groundwater connections are present, and not indefinitely. As warming persists, storage can be depleted and connections lost. Similarly, in the arid western US warming does not result in significant groundwater changes because this area is already largely water limited. The direct response of shallow groundwater storage to warming demonstrates the strong and early effect that low to moderate warming may have on groundwater storage and evapotranspiration.

54 ENVIRONMENTAL SCIENCES↗

Irrigation-induced land water depletion aggravated by climate change

Agricultural irrigation has experienced rapid expansion, and its growing freshwater consumption is potentially exacerbating water scarcity issues. Previous studies predominantly relied on observations or land-only simulations, often neglecting land–atmosphere interactions or failing to capture long-term evolution. We therefore analyse the effects of historical irrigation expansion on water fluxes and resources using seven Earth system models. Here we show that irrigation expansion in many regions substantially decreases the net water influx from the atmosphere to land, further aggravating the existing drying trends caused by climate change. For example, irrigation expansion changed the trend of this net influx from −0.664 ( ± 0.283) to −1.461 ( ± 0.261) mm yr −2 in South Asia after 1960. Consequently, the local terrestrial water storage depletion rate is substantially enlarged by irrigation expansion (for example, from −2.559 ( ± 0.094) to −16.008 ( ± 0.557) mm yr −1 ). Our results attribute the land water loss to irrigation expansion and climate change, calling for immediate solutions to tackle the negative trends.

irrigation↗

Properties and Controlling Processes of Aerosol and Cloud Condensation Nuclei in Marine Boundary Layer Over Eastern North Atlantic (Final Report)

Remote marine low cloud systems have large spatial and temporal coverages. Because of their relatively low optical thickness and background aerosol concentrations, marine low clouds are particularly susceptible to perturbations in aerosols associated with anthropogenic emissions. Recent studies show that a large fraction of the global aerosol indirect forcing can be attributed to changes in marine low clouds, despite their relatively long distance from most anthropogenic sources. Currently, there remain large uncertainties in the magnitude of the global aerosol indirect radiative forcing. One major contribution to this large uncertainty derives from poor understanding of the aerosols under natural conditions, the perturbation by anthropogenic emissions, and the processes that drive them. In this project, we took advantage of the comprehensive, multi-dimensional characterization of meteorological parameters, trace gases, aerosol, cloud, and drizzle fields during the Aerosol and Cloud Experiments in Eastern North Atlantic (ACE-ENA) and the long-term measurements at the Atmospheric Radiation Measurement (ARM) Eastern North Atlantic (ENA) site on Graciosa Island. By synergistically combining in situ measurements, remote sensing retrievals, and model simulations, we quantified aerosol properties in the ENA, their vertical profiles, and variations with season, synoptic state, and airmass. We investigated key aerosol processes, including entrainment of free troposphere aerosol into marine boundary layer, new particle formation, growth of Aitken mode particles, and the influences of these processes on the cloud condensation nuclei (CCN) population in marine boundary layer. The representation of key aerosol processes in the Energy Exascale Earth System Model (E3SM) was evaluated using the observations. This project has led to 26 publications, which are listed at the end of this report.

54 ENVIRONMENTAL SCIENCES↗

Concurrent Changes in Extreme Hydroclimate Events in the Colorado River Basin

Extreme events resulting in catastrophic damage have more than doubled in the last five years, costing hundreds of lives and thousands of homes, and heavily undermining regional economic stability. At present, most of these hydroclimatic extreme events are documented by the media as individual events; however, in scientific terms, many are better understood as concurrent events—concurrent extremes of both temperature and precipitation (e.g., drought, floods). This paper considers concurrent changes in hydroclimate extremes, including heatwaves, drought, flooding, and low flows, in six historical-to-future (1970–1999, 2070–2099) Earth System Model (ESM) climate scenarios for the Colorado River basin. Results indicate that temperature-driven Impacts (heatwaves, drought) have the strongest responses while precipitation-driven Impacts have weaker responses. All Impacts exhibit an increase in magnitude from synoptic to annual time scales, with heatwaves increasing in strength about three times at the annual time scale versus the synoptic, while low flows only increase slightly. Critical watersheds in the Colorado were identified, highlighting the Blue River basin, Uncompahgre, East Taylor, Salt/Verde watersheds, locations of important water infrastructures, water resources, and hydrological research. Our results indicate that concurrent extreme hydroclimate events are projected to increase in the future and intensify within critical regions of the Colorado River basin. Considering extreme hydroclimate events concurrently is an important step towards linking economic and social effects of these events and their associated instabilities on a regional scale.

54 ENVIRONMENTAL SCIENCES↗

Challenges in studying water fluxes within the soil-plant-atmosphere continuum: A tracer-based perspective on pathways to progress

Tracing and quantifying water fluxes in the hydrological cycle is crucial for understanding the current state of ecohydrological systems and their vulnerability to environmental change. Especially the interface between ecosystems and the atmosphere that is strongly mediated by plants is important to meaningfully describe ecohydrological system functioning. Many of the dynamic interactions generated by water fluxes between soil, plant and the atmosphere are not well understood, which is partly due to a lack of interdisciplinary research. This opinion paper reflects the outcome of a discussion among hydrologists, plant ecophysiologists and soil scientists on open questions and new opportunities for collaborative research on the topic "water fluxes in the soil-plant-atmosphere continuum" especially focusing on environmental and artificial tracers. We emphasize the need for a multi-scale experimental approach, where a hypothesis is tested at multiple spatial scales and under diverse environmental conditions to better describe the small-scale processes (i.e., causes) that lead to large-scale patterns of ecosystem functioning (i.e., consequences). Novel in-situ, high-frequency measurement techniques offer the opportunity to sample data at a high spatial and temporal resolution needed to understand the underlying processes. Here we advocate for a combination of long-term natural abundance measurements and event-based approaches. Multiple environmental and artificial tracers, such as stable isotopes, and a suite of experimental and analytical approaches should be combined to complement information gained by different methods. Virtual experiments using process-based models should be used to inform sampling campaigns and field experiments, e.g., to improve experimental designs and to simulate experimental outcomes. On the other hand, experimental data are a pre-requisite to improve our currently incomplete models. Interdisciplinary collaboration will help to overcome research gaps that overlap across different earth system science fields and help to generate a more holistic view of water fluxes between soil, plant and atmosphere in diverse ecosystems.

54 ENVIRONMENTAL SCIENCES↗

Atmosphere-Biosphere Interaction Study (Field Campaign Report)

The objective of this campaign was to collect hyperspectral reflectance measurements and complementary optical and infrared data of land surface (i.e., vegetation and soil) throughout the growing season (i.e., from March to September) at high temporal frequency (e.g., < 30 minutes). To collect the measurements, the EcoSpec system, a tower-based sensor system housing seven optical and infrared sensors that was developed at Argonne National Laboratory, was deployed to the U.S. Department of Energy Atmospheric Radiation Measurement (ARM) user facility Southern Great Plains (SGP) atmospheric observatory (within the crop, south of the eddy correlation flux measurement system [ECOR] tower). The EcoSpec system consists of a spectroradiometer having 2,151 channels (350–2,500 nm), red-greenblue (RGB) camera, diffuse radiometer, thermal infrared (TIR) sensors, albedometers, Photochemical Reflectance Index (PRI) sensors, Normalized Difference Vegetation Index (NDVI) sensors, and wetness sensor. The system is designed to collect an array of land surface measurements, including hyperspectral reflectance of plants and soils, RGB photos capturing contextual information of land surface, radiant temperature of sky and land surface, incoming and reflected shortwave radiation, and incoming radiation components (direct and diffuse) as frequent as every minute from up to 12 discrete positions around the tower. This would allow us to capture inherent heterogeneity of the terrestrial ecosystem and to obtain representative spectral signatures of the plant and soil surfaces via averaging those measurements. The SGP site provides advanced atmospheric and other physical measurements that are essential for our investigation of ecosystem functions using hyperspectral reflectance measurements of land surfaces. Located in the mid-latitudes, the SGP site exhibits dynamic seasonality and a wide range of meteorological variability and flux properties that make it ideal for studying interactions between near-surface atmosphere and biosphere across a range of environmental conditions. The atmosphere, plants, and soil control terrestrial carbon and water cycles. To make more accurate climate forecasts, researchers need to understand ecosystem dynamics at the biosphere–atmosphere interface. Even though our ability to forecast ecosystem dynamics and climate at both the regional and global scales has advanced significantly, we still do not know how local phenomena affecting water and carbon fluxes—such as diurnal variations—interrelate with large-scale atmospheric/climate phenomena, and vice versa. The lack of a more thorough understanding of ecosystem dynamics and climate interactions is a key knowledge gap for the improvement of earth system models (ESMs).

54 ENVIRONMENTAL SCIENCES↗

Litter Production and Foliar Nutrient Resorption in Pioneer and Non-Pioneer Species in a Selective Logging Experiment in the Central Amazon, BIONTE, ZF-2, Manaus, 2022-23

This dataset was collected near the city of Manaus, Brazil, at the Experimental Station of Tropical Forestry (EEST, aka “ZF2”), inside the BIONTE (BIOmass and NuTrient Experiment). The experiment included three levels of increasing selective logging intensity, along with control, with 1-hectare permanent plots (12 total) located at the center of 4-hectare treatment plots. The vegetation has a high floristic diversity, the soils of the region are poor in nutrients, and the topography is characterized by plateaus (where BIONTE is located), and also valley bottoms and slopes. Three treatments of differing logging intensities were applied in the BIONTE experiment (T1, T2 and T3). The study was conducted in Treatment 3 (Block I – permanent plot), which represents the most intensive logging treatment, with 69% of the basal area (m²∙ha⁻¹) removed in 1988. The present dataset spans the period from May 1, 2022, to May 1, 2023. The data package includes leaf_nutrient_data, litterfall_total_data, leaf_litterfall_species_specific_data, and species_info, all provided in .csv format. These formats allow users to process and analyze the data in various software applications and programming languages, such as Python and R. This dataset was collected to advance knowledge on nutrient cycling in Amazonian forests, specifically distinguishing between species with two distinct functional traits: fast-growing and slow-growing. It also aims to improve Earth System Models, such as the E3SM Functionally Assembled Terrestrial Ecosystem Simulator (FATES). Additionally, it was used in a paper currently in preparation (Carvalho et al., in prep.), which aims to quantify seasonal litter production and foliar nutrient resorption in pioneer (fast-growing) and non-pioneer (slow-growing) tree species in the central Amazon. Specifically, it seeks to answer two key questions: 1) Is there a difference in leaf litter production, leaf nutrient flux and leaf nutrient concentration between pioneers and non-pioneers species? Is there a difference in the efficiency of foliar nutrient resorption between pioneers and non-pioneers species?

54 ENVIRONMENTAL SCIENCES↗

Wildfire Smoke Demonstrates Significant and Predictable Black Carbon Light Absorption Enhancements

Abstract Black carbon (BC) is estimated to have the second largest anthropogenic radiative forcing in earth‐systems models (ESMs), but there is significant uncertainty in its impact due to complex mixing with organics. Laboratory‐generated particles show that co‐mixed non‐absorbing material enhances absorption by BC by a factor of 2–3.5 as predicted by optical models. However, weak or no enhancements are often reported for field studies. The cause of lower‐than‐expected absorption is not well understood and implies a lower radiative impact of BC compared to how many ESMs currently treat aerosols. By analyzing BC aerosol particle‐by‐particle we reconcile observed and expected absorption for ambient smoke plumes varying in geographic origin, fuel types, burn conditions, atmospheric age and transport. Although particle‐by‐particle tracking is computationally prohibitive for sophisticated ESMs we show that realistic BC absorption is reliably estimated by bulk properties of the plume providing a suitable parameterization to constrain black carbon radiative forcing.

54 ENVIRONMENTAL SCIENCES↗

A More La Niña–Like Response to Radiative Forcing after Flux Adjustment in CESM2

Abstract In response to greenhouse gas forcing, most coupled global climate models project the tropical Pacific SST trend toward an “El Niño–like” state, with a reduced zonal SST gradient and a weakened Walker circulation. However, observations over the last five decades reveal a trend toward a more “La Niña–like” state with a strengthening zonal SST gradient. Recent research indicates that the identified trend differences are unlikely to be entirely due to internal variability and probably result, at least in part, from systematic model biases. In this study, Community Earth System Model, version 2 (CESM2), is used to explore how mean-state biases within the model may influence its forced response to radiative forcing in the tropical Pacific. The results show that using flux adjustment to reduce the mean-state bias in CESM2 over the tropical regions results in a more La Niña–like trend pattern in the tropical Pacific, with a strengthening of the tropical Pacific zonal SST gradient and a relatively enhanced Walker circulation, as hypothesized to occur if the ocean thermostat mechanism is stronger than the atmospheric mechanisms which by themselves would weaken the Walker circulation. We also find that the historical strengthening of the tropical Pacific zonal gradient is transient but persists into the near term in a high-emissions future warming scenario. These results suggest the potential of flux adjustment as a method for developing alternative projections that represent a wider range of possible future tropical Pacific warming scenarios, especially for a better understanding of regional patterns of climate risk in the near term.

Zhuo, Jing-Yi [Lamont-Doherty Earth Observatory, C↗

Shrub Expansion Simulations at Trail Valley Creek Tundra site using E3SM Land Model (ELM) Arctic-focused Version

The warming of the Arctic is causing substantial compositional, structural, and functional changes in tundra vegetation including shrub and densification in parts of the Arctic. Assessing the impact of these changes in vegetation composition on the Arctic’s carbon and energy budgets is important to constrain projected local and global surface-atmosphere exchanges. We conduct a sensitivity analysis of the projected surface energy fluxes, soil carbon pools, and carbon dioxide fluxes (net ecosystem exchange, gross primary production, and ecosystem respiration) between present day and 2100 to different shrub expansion rates and air temperature increases under future emission scenarios (intermediate – RCP4.5, and high – RCP8.5) using the Arctic-focused version of the Energy Exascale Earth System Model (E3SM) Land Model (ELM). We focus on Trail Valley Creek (TVC), a mineral upland tundra site located in the western Canadian Arctic, which is experiencing tall shrub densification and expansion. In this study, we run TVC under two different warming scenarios RCP4.5 and RCP 8.5 and simulate different shrubification rates projected until year 2100. In this repository, we include all the forcing, input, parameters, and output data corresponding to all the simulations performed. flmd.csv includes a detailed description of the datasets files.

54 ENVIRONMENTAL SCIENCES↗

CO2 and CH4 leaf-level fluxes and soil porewater concentrations from common vegetation patches in Louisiana’s coastal wetlands

This dataset contains leaf-level flux and soil porewater concentration measurements of carbon dioxide (CO2) and methane (CH4 ) in plots in the footprint of Ameriflux sites US-LA2 and US-LA3. Leaf fluxes in US-LA2 were measured on patches dominated by Sagittaria lancifolia and co-dominated by Sagittaria lancifolia and Typha latifolia. In US-LA3, fluxes were measured from distinct Juncus roemerianus and Spartina alterniflora patches. The porewater concentrations were collected across a vertical profile (~50 cm depth) at centric locations within 25 m2 plots where we measured the leaf fluxes. US-LA3 included an additional set of measurements in open water spots. We aimed to evaluate differences in leaf fluxes and porewater pools of CO2 and CH4 of representative ecohydrological patches across a salinity gradient. We also used this dataset to help develop ELM-Wet, a more realistic representation of wetland carbon biogeochemical processes within the U.S. Department of Energy’s Energy Exascale Earth System Model (E3SM) Land Model version 1 (ELM v.1). The files can be opened with regular text editors or spreadsheet programs. Version 2.0 (8/26/2025): This is the latest version of this dataset. The update includes additional samples of soil porewater CH4/CO2 concentrations from June-2021 to November-2022, as well as minor adjustments made to V1 samples via changing Henry's solubility to account for porewater salinity. Additionally leaf-level measurments of spectral indices, PSRI, NDVI, and PRI have been added to complement Leaf-level flux measurements. All V1 data sets have been integrated into V2 sheets, ensuring data from the previous version is contained with the additional samples and consistent with V2 metadata.

54 ENVIRONMENTAL SCIENCES↗

Global gross primary production simulated by three terrestrial biosphere models to analyse smoothing of limiting rates

Processed global terrestrial gross primary production (GPP) from three terrestrial biosphere models (TBMs). Two simulations were run with each model, one that selects the rate limiting process of photosynthesis based on the minimum of the potential rates (following Farquhar et al. 1980) and the other that uses quadratic smoothing between the potential rates (following Collatz et al. 1991). The three TBMs are: 1) the Energy Exascale Earth System Model (E3SM) land model (ELM) (release: v1.1.0) (Burrows et al., 2020), a coupled carbon, nitrogen, and phosphorus model with sun/shade big-leaf canopy photosynthesis scaling. 2) The Functionally Assembled Terrestrial Ecosystem Simulator (FATES, tag: sci.1.30.0_api.8.0.0) (Koven et al., 2019), coupled with the Community Land Model (CLM, version 5; Lawrence et al., 2019), a carbon-only vegetation demography model with multi-leaf and multi-canopy layers for scaling photosynthesis, and with a leaf area index (LAI) optimisation scheme. 3) The Sheffield Dynamic Global Vegetation Model (SDGVM, tag: Walkeretal2020_GCB), a carbon-only model with multi-leaf-layer canopy photosynthesis scaling and also with an LAI optimisation scheme (Walker et al., 2017; Woodward & Lomas, 2004). These data were generated to evaluate the effect of Collatz smoothing of limiting rates on simulated GPP across a range of terrestrial biosphere models. Analysis of results can be found in Walker et al. (2020). This dataset is comprised of R scripts to process data and produce figures and .RDS files, a binary format specific to R and read with the 'readRDS' function.

54 ENVIRONMENTAL SCIENCES↗

Performance Results on CPU/GPU Exascale Architectures for OMEGA: The Ocean Model for E3SM Global Applications

The US Department of Energy (DOE) conducts climate simulations on some of the world’s largest supercomputers. These exascale machines use heterogeneous architectures with both CPUs and GPUs, and scientific codes must adapt to make full use of this computing power. Los Alamos National Lab is developing Omega: The Ocean Model for E3SM Global Applications, which is specifically designed for modern exascale computers. It uses external libraries that have been optimized for a variety of architectures to run on different supercomputers. Omega is an unstructured-mesh ocean model based on TRiSK numerical methods. It will be the new ocean component of the DOE’s Energy Exascale Earth System Model (E3SM). The algorithms in Omega follow those of the current ocean component, MPAS-Ocean, but it will be written in C++ rather than Fortran to take advantage of the Kokkos performance portability library. Omega spatial operators are written as Kokkos kernels to run efficiently on both CPUs and GPUs. Work on Omega began in 2023 with a new C++ framework for unstructured mesh partitioning, halo exchanges, parallel IO, and Kokkos interfaces. The current version, Omega-0, is being developed to solve the shallow water equations and at present includes all of the tendency terms but not time stepping. Here we share the results of Omega-0 verification and performance testing. Verification includes unit tests implemented with CTest as well as convergence tests in Polaris, an in-house python package with a large suite of test problems. Performance tests compare simulations conducted on CPUs versus GPUs and across different architectures: tests are run on Frontier, which has AMD “Optimized 3rd Gen EPYC” CPUs and AMD MI250X GPUs, as well as Perlmutter, which is composed of AMD EPYC 7763 CPUs and NVIDIA A100 GPUs.

58 GEOSCIENCES↗