Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Earth System Models”

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 1,549 records · Page 86

Improved Time-Stepping Methods in Global to Regional Ocean Modeling (Annual Status Report 2020)

Time stepping algorithms are an important part of ocean models, and strongly influence both the accuracy of solution and performance. There have been a number of projects investigating various improvements for ocean time-stepping schemes in the Model for Prediction Across Scales-Ocean (MPAS-Ocean), a component of the DOE Energy Exascale Earth System Model. Ocean dynamics include fast surface gravity waves, which are two-dimensional, and slower internal waves, which are three-dimensional, so ocean models use a split time-stepping scheme that separates these barotropic and baroclinic modes for efficiency. MPAS-Ocean runs on variable-resolution horizontal meshes, and must scale to tens of thousands of cores and millions of horizontal gridcells. Ocean models require time stepping algorithms that are customized to these needs, and which are tuned for performance on various resolutions and architectures.

58 GEOSCIENCES↗

SCREAM Beijing Flood RRM Technical Note

The Beijing Convection-Permitting (CP) Regionally Refined Model (RRM) used for 2023 Jingjinji Flood event is developed based on the Simple Cloud-Resolving E3SM Atmosphere Model (SCREAM) version 0 (Fortran code) under the United States (U.S) Department of Energy (DOE) Energy Exascale Earth System Model (E3SM) project (Caldwell et al., 2021) and the regionally refined model (RRM) configuration (Tang et al., 2019).

58 GEOSCIENCES↗

EAMxx 3.25 km Western Pacific RRM Technical Note

The Western Pacific Convection-Permitting (CP) Regionally Refined Model (RRM) used for 2023 Jingjinji Flood event is developed based on the Simple Cloud-Resolving E3SM Atmosphere Model (SCREAM) version 0 (Fortran code) under the United States (U.S) Department of Energy (DOE) Energy Exascale Earth System Model (E3SM) project (Caldwell et al., 2021) and the regionally refined model (RRM) configuration (Tang et al., 2019, 2023; Zhang et al., 2024; Bogenschutz et al., 2024).

58 GEOSCIENCES↗

A multi-algorithm approach for modeling coastal wetland eco-geomorphology

Coastal wetlands play an important role in the global water and biogeochemical cycles. Climate change makes it more difficult for these ecosystems to adapt to the fluctuation in sea levels and other environmental changes. Given the importance of eco-geomorphological processes for coastal wetland resilience, many eco-geomorphology models differing in complexity and numerical schemes have been developed in recent decades. However, their divergent estimates of the response of coastal wetlands to climate change indicate that substantial structural uncertainties exist in these models. To investigate the structural uncertainty of coastal wetland eco-geomorphology models, we developed a multi-algorithm model framework of eco-geomorphological processes, such as mineral accretion and organic matter accretion, within a single hydrodynamics model. The framework is designed to explore possible ways to represent coastal wetland eco-geomorphology in Earth system models and reduce the related uncertainties in global applications. We tested this model framework at three representative coastal wetland sites: two saltmarsh wetlands (Venice Lagoon and Plum Island Estuary) and a mangrove wetland (Hunter Estuary). Through the model–data comparison, we showed the importance of using a multi-algorithm ensemble approach for more robust predictions of the evolution of coastal wetlands. We also found that more observations of mineral and organic matter accretion at different elevations of coastal wetlands and evaluation of the coastal wetland models at different sites in diverse environments can help reduce the model uncertainty.

58 GEOSCIENCES↗

Spatial heterogeneity effects on land surface modeling of water and energy partitioning

Abstract. Understanding the influence of land surface heterogeneity on surface water and energy fluxes is crucial for modeling earth system variability and change. This study investigates the effects of four dominant heterogeneity sources on land surface modeling, including atmospheric forcing (ATM), soil properties (SOIL), land use and land cover (LULC), and topography (TOPO). Our analysis focused on their impacts on the partitioning of precipitation (P) into evapotranspiration (ET) and runoff (R), partitioning of net radiation into sensible heat and latent heat, and corresponding water and energy fluxes. An initial set of 16 experiments were performed over the continental US (CONUS) using the E3SM land model (ELMv1) with different combinations of heterogeneous and homogeneous datasets. The Sobol' total and first-order sensitivity indices were utilized to quantify the relative importance of the four heterogeneity sources. Sobol' total sensitivity index measures the total heterogeneity effects induced by a given heterogeneity source, consisting of the contribution from its own heterogeneity (i.e., the first-order index) and its interactions with other heterogeneity sources. ATM and LULC are the most dominant heterogeneity sources in determining spatial variability of water and energy partitioning, mainly contributed by their own heterogeneity and slightly contributed by their interactions with other heterogeneity sources. Their heterogeneity effects are complementary, both spatially and temporally. The overall impacts of SOIL and TOPO are negligible, except TOPO dominates the spatial variability of R/P across the transitional climate zone between the arid western and humid eastern CONUS. Accounting for more heterogeneity sources improves the simulated spatial variability of water and energy fluxes when compared with ERA5-Land reanalysis dataset. An additional set of 13 experiments identified the most critical components within each heterogeneity source, which are precipitation, temperature, and longwave radiation for ATM, soil texture, and soil color for SOIL and maximum fractional saturated area parameter for TOPO.

54 ENVIRONMENTAL SCIENCES↗

A Community Atmosphere Model with Superparameterized Clouds

In 1999, National Center for Atmospheric Research (NCAR) scientists Wojciech Grabowski and Piotr Smolarkiewicz created a "multiscale" atmospheric model in which the physical processes associated with clouds were represented by running a simple high-resolution model within each grid column of a lowresolution global model. In idealized experiments, they found that the multiscale model produced promising simulations of organized tropical convection, which other models had struggled to produce. Inspired by their results, Colorado State University (CSU) scientists Marat Khairoutdinov and David Randall created a multiscale version of the Community Atmosphere Model (CAM). They removed the cloud parameterizations of the CAM, and replaced them with Khairoutdinov's high-resolution cloud model. They dubbed the embedded cloud model a "super-parameterization," and the modified CAM is now called the "SP-CAM." Over the next several years, many scientists, from many institutions, have explored the ability of the SP-CAM to simulate tropical weather systems, the day-night changes of precipitation, the Asian and African monsoons, and a number of other climate processes. Cristiana Stan of the Center for Ocean-Land-Atmosphere Interactions found that the SP-CAM gives improved results when coupled to an ocean model, and follow-on studies have explored the SP-CAM's utility when used as the atmospheric component of the Community Earth System Model. Much of this research has been performed under the auspices of the Center for Multiscale Modeling of Atmospheric Processes, a National Science Foundation (NSF) Science and Technology Center for which the lead institution is CSU.

atmosphere↗

Antarctic Ice Sheet Freshwater Discharge Drives Substantial Southern Ocean Changes Over the 21st Century

Multidecadal satellite observations indicate that the Antarctic Ice Sheet (AIS) is losing mass at an accelerating rate, which has the potential to impact many aspects of the climate system. While previous studies demonstrated the importance of AIS freshwater (FW) discharge for regional and global climate processes using climate model experiments, many have applied unrealistic FW forcings. Here, we explore potential Southern Ocean (SO) impacts of realistic AIS mass loss over the 21st century in the Community Earth System Model version 2 by applying observation-based historical and ice sheet model-based future AIS FW forcing. The added FW reduces wintertime deep convective area by 72% while retaining 83% more sea ice. Congruent with other studies, we find the increased FW discharge extensively impacts local and remote SO surface and subsurface temperature and stratification. These results demonstrate the necessity of accounting for AIS mass loss in global climate models for projecting future climate.

Geology↗

A Scientific Data Provenance Harvester for Distributed Applications

Data provenance provides a way for scientists to observe how experimental data originates, conveys process history, and explains influential factors such as experimental rationale and associated environmental factors from system metrics measured at runtime. The US Department of Energy Office of Science Integrated end-to-end Performance Prediction and Diagnosis for Extreme Scientific Workflows (IPPD) project has developed a provenance harvester that is capable of collecting observations from file based evidence typically produced by distributed applications. To achieve this, file based evidence is extracted and transformed into an intermediate data format inspired in part by W3C CSV on the Web recommendations, called the Harvester Provenance Application Interface (HAPI) syntax. This syntax provides a general means to pre-stage provenance into messages that are both human readable and capable of being written to a provenance store, Provenance Environment (ProvEn). HAPI is being applied to harvest provenance from climate ensemble runs for Accelerated Climate Modeling for Energy (ACME) project funded under the U.S. Department of Energy’s Office of Biological and Environmental Research (BER) Earth System Modeling (ESM) program. ACME informally provides provenance in a native form through configuration files, directory structures, and log files that contain success/failure indicators, code traces, and performance measurements. Because of its generic format, HAPI is also being applied to harvest tabular job management provenance from Belle II DIRAC scheduler relational database tables as well as other scientific applications that log provenance related information.

Stephan, Eric G.↗

Drivers of coupled climate model biases in representing Labrador Sea convection

Abstract This study investigates the representation of ocean convection in the Labrador Sea in seven Earth System Models (ESMs) from the Coupled Model Intercomparison Project Phase 5 and 6 datasets. The relative role of the oceanic and atmospheric biases in the subpolar North Atlantic gyre are explored using regional ocean simulations where the atmospheric forcing or the ocean initial and boundary conditions are replaced by reanalysis data in the absence of interactive air-sea coupling. Commonalities and differences among model behaviors are discussed with the objective of finding a pathway forward to improve the representation of the ocean mean state and variability in a region of fundamental importance for climate variability and change. Results highlight that an improved representation of ocean stratification in the North Atlantic subpolar gyre is urgently needed to constrain future climate change projections. While improving the ocean model resolution in the North Atlantic alone may contribute a better representation of both boundary currents and propagation of heat and freshwater anomalies into the Labrador Sea, it may not be sufficient. Addressing the atmospheric heat flux bias with better resolution in the atmosphere and land topography may allow for deep convection to occur in the Labrador Sea in some of the models that miss it entirely, but the greatest priority remains improving the representation of ocean stratification.

Liu, Guangpeng↗

Enhancing Extreme Precipitation Predictions With Dynamical Downscaling: A Convection‐Permitting Modeling Study in Texas and Oklahoma

Precipitation in the Southern Plains of the United States is relatively well depicted by the Community Earth System Model (CESM). However, despite its ability to capture seasonal mean precipitation anomalies, CESM consistently underestimates extreme pluvial and drought events, rendering it an insufficient tool for extending simulation lead times for exceptional events, such as the abnormally dry May 2011, which helped drive Texas into its worst period of drought in more than a century, and the abnormally wet May 2015, which led to widespread flooding in that state. Ensemble‐based regional climate experiments are completed for the two extreme years using Weather Research and Forecasting model (WRF) and downscaled from CESM. WRF simulations are at convection‐permitting grid resolution for improved physical representation of simulated precipitation over the Southern Great Plains. By integrating convection‐permitting models (CPMs) into each individual member of a CESM ocean data assimilation ensemble, this study demonstrates that high‐resolution dynamical downscaling can improve model skillfulness at capturing these two events and is thus a potentially useful tool for forecasting extremely high and extremely low precipitation events at subseasonal or even seasonal lead times.

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↗

A Bayesian Approach to Regional Decadal Predictability: Sparse Parameter Estimation in High-Dimensional Linear Inverse Models of High-Latitude Sea Surface Temperature Variability

Stochastic reduced models are an important tool in climate systems whose many spatial and temporal scales cannot be fully discretized or underlying physics may not be fully accounted for. One form of reduced model, the linear inverse model (LIM), has been widely used for regional climate predictability studies—typically focusing more on tropical or midlatitude studies. However, most LIM fitting techniques rely on point estimation techniques deriving from fluctuation–dissipation theory. In this methodological study we explore the use of Bayesian inference techniques for LIM parameter estimation of sea surface temperature (SST), to quantify the skillful decadal predictability of Bayesian LIM models at high latitudes. We show that Bayesian methods, when compared to traditional point estimation methods for LIM-type models, provide better calibrated probabilistic skill, while simultaneously providing better point estimates due to the regularization effect of the prior distribution in high-dimensional problems. We compare the effect of several priors, as well as maximum likelihood estimates, on 1) estimating parameter values on a perfect model experiment and 2) producing calibrated 1-yr SST anomaly forecast distributions using a preindustrial control run of the Community Earth System Model (CESM). Finally, we employ a host of probabilistic skill metrics to determine the extent to which an LIM can forecast SST anomalies at high latitudes. We find that the choice of prior distribution has an appreciable impact on estimation outcomes, and priors that emphasize physically relevant properties enhance the model’s ability to capture variability of SST anomalies.

54 ENVIRONMENTAL SCIENCES↗

Impact of Earth Greening on the Terrestrial Water Cycle

Leaf area index (LAI) is increasing throughout the globe, implying Earth greening. Global modeling studies support this contention, yet satellite observations and model simulations have never been directly compared. Here, for the first time, a coupled land–climate model was used to quantify the potential impact of the satellite-observed Earth greening over the past 30 years on the terrestrial water cycle. The global LAI enhancement of 8% between the early 1980s and the early 2010s is modeled to have caused increases of 12.0 ± 2.4 mm yr −1 in evapotranspiration and 12.1 ± 2.7 mm yr −1 in precipitation—about 55% ± 25% and 28% ± 6% of the observed increases in land evapotranspiration and precipitation, respectively. In wet regions, the greening did not significantly decrease runoff and soil moisture because it intensified moisture recycling through a coincident increase of evapotranspiration and precipitation. But in dry regions, including the Sahel, west Asia, northern India, the western United States, and the Mediterranean coast, the greening was modeled to significantly decrease soil moisture through its coupling with the atmospheric water cycle. This modeled soil moisture response, however, might have biases resulting from the precipitation biases in the model. For example, the model dry bias might have underestimated the soil moisture response in the observed dry area (e.g., the Sahel and northern India) given that the modeled soil moisture is near the wilting point. Thus, an accurate representation of precipitation and its feedbacks in Earth system models is essential for simulations and predictions of how soil moisture responds to LAI changes, and therefore how the terrestrial water cycle responds to climate change.

54 ENVIRONMENTAL SCIENCES↗

Carbon stored in live ponderosa pines in the Sierra Nevada will not return to pre-drought (2012) levels during the 21st century due to bark beetle outbreaks

Outbreaks of several bark beetle species can develop rapidly in response to drought and may result in large transfers of carbon (C) stored in live trees to C stored in dead trees (10s of Tg C yr -1 in the western U.S. alone), which over time will be released back to the atmosphere. The western pine beetle (WPB) outbreak incited by the 2012–2015 mega-drought in the Sierra Nevada, California, U.S., could portend more frequent and/or severe bark beetle outbreaks as the temperature warms and drought frequency and intensity increase in the future. However, changes in the frequency and/or severity (resultant levels of host tree mortality) of beetle outbreaks are difficult to predict as outbreaks are complex with non-linear and eruptive processes primarily driven by interactions among beetle populations, the demography of hosts and other tree species, and climate and weather. Using an insect phenology and tree defense model, we projected the future likelihood of WPB outbreaks in the Sierra Nevada with climate drivers from different Earth System Models. Our goal was to understand how host (ponderosa pine, PIPO) recovery and future warming and drought affect the frequency and severity of WPB outbreaks and their C consequences. Our projections suggested that by 2100 the C stored in live PIPO (mean: 1.98 kg C m -2 , 95% CI: 1.74–2.21 kg C m -2 ) will not return to levels that occurred before the 2012–2015 drought (2012: ~2.30 kg C m –2 ) due to future WPB outbreaks. However, differences in climate models indicate a wide range of possible WPB outbreak frequencies and severities. Our results suggest that total plot basal area is the most significant factor in the mortality rate of PIPO by WPB in any given year, followed by drought severity and temperature. High levels of host basal area, higher temperature, and extreme drought all contribute to the frequency and severity of future WPB outbreaks. While PIPO basal area may decline under increased drought and warming, limiting high-stand basal area (>60 m 2 ha -1 ) may reduce the severity of future WPB outbreaks in the Sierra Nevada.

54 ENVIRONMENTAL SCIENCES↗

Global Modeling, Field Campaigns, Upscaling and Ray Desjardins

In the early 1980's, it became apparent that land surface radiation and energy budgets were unrealistically represented in Global Circulation models (GCM's), Shortly thereafter, it became clear that the land carbon budget was also poorly represented in Earth System Models (ESM's), A number of scientific communities, including GCM/ESM modelers, micrometeorologists, satellite data specialists and plant physiologists, came together to design field experiments that could be used to develop and validate the contemporary prototype land surface models. These experiments were designed to measure land surface fluxes of radiation, heat, water vapor and CO2 using a network of flux towers and other plot-scale techniques, coincident with satellite measurements of related state variables, The interdisciplinary teams involved in these experiments quickly became aware of the scale gap between plot-scale measurements (approx 10 - 100m), satellite measurements (100m - 10 km), and GCM grid areas (l0 - 200km). At the time, there was no established flux measurement capability to bridge these scale gaps. Then, a Canadian science learn led by Ray Desjardins started to actively participate in the design and execution of the experiments, with airborne eddy correlation providing the radically innovative bridge across the scale gaps, In a succession of brilliantly executed field campaigns followed up by convincing scientific analyses, they demonstrated that airborne eddy correlation allied with satellite data was the most powerful upscaling tool available to the community, The rest is history: the realism and credibility of weather and climate models has been enormously improved enormously over the last 25 years with immense benefits to the public and policymakers.

Sellers, P. J.↗

CMIP7 data request: impacts and adaptation priorities and opportunities

The Coupled Model Intercomparison Project Phase 7 (CMIP7) undertook an extensive process to gather community input and refine data requests related to impacts and adaptation applications of Earth System Model (ESM) outputs. The Impacts and Adaptation (I&A) Data Request Team worked with CMIP7 leadership to distribute an open solicitation across many communities that use climate model outputs requesting inputs for new and existing variables, the most applicable temporal characteristics, and groupings of variables that together allow for specific application opportunities. This input was then collated and translated into CMIP7 standard templates for inclusion in the broader data request, leading to 13 I&A data request opportunities, 60 variable groups and 539 unique variables sought by vulnerability, impacts, adaptation, and climate services user communities. Here, we describe these opportunities and variable groups, as well as new insights into how ESM groups can prioritize outputs that set off a chain of further analyses, ultimately informing decisions impacting society and natural systems. These include an emphasis on high-resolution outputs to allow further modeling of climate impacts at regional and local scales, improved representation of extreme weather events, enhanced accuracy of downscaling and bias-adjustment techniques, and support for more detailed assessments for decision-making in adaptation and mitigation strategies. There is also broad interest in more extensive provisioning of two-dimensional variables at the Earth's surface, prioritizing experiments that enhance our understanding of both the recent past and future scenarios, and providing outputs that allow further downscaling and bias adjustment. We emphasize that variable groups are the fundamental level at which to engage with the I&A data request, matching the scale of input and the way output provision enables specific I&A applications. Given resource constraints, we applaud CMIP7 efforts to foster strong engagement and communication between ESM groups and the I&A team to build consensus around prudent compromises in priority variables, temporal resolutions, simulation experiments, time subsets, and ensemble members.

Ruane, Alex C. [NASA Goddard Inst. for Space Studi↗

Climate change impacts on net primary production (NPP) and export production (EP) regulated by increasing stratification and phytoplankton community structure in the CMIP5 models

Abstract. We examine climate change impacts on net primary production (NPP) and export production (sinking particulate flux; EP) with simulations from nine Earth system models (ESMs) performed in the framework of the fifth phase of the Coupled Model Intercomparison Project (CMIP5). Global NPP and EP are reduced by the end of the century for the intense warming scenario of Representative Concentration Pathway (RCP) 8.5. Relative to the 1990s, NPP in the 2090s is reduced by 2–16 % and EP by 7–18 %. The models with the largest increases in stratification (and largest relative declines in NPP and EP) also show the largest positive biases in stratification for the contemporary period, suggesting overestimation of climate change impacts on NPP and EP. All of the CMIP5 models show an increase in stratification in response to surface–ocean warming and freshening, which is accompanied by decreases in surface nutrients, NPP and EP. There is considerable variability across the models in the magnitudes of NPP, EP, surface nutrient concentrations and their perturbations by climate change. The negative response of NPP and EP to increasing stratification reflects primarily a bottom-up control, as upward nutrient flux declines at the global scale. Models with dynamic phytoplankton community structure show larger declines in EP than in NPP. This pattern is driven by phytoplankton community composition shifts, with reductions in productivity by large phytoplankton as smaller phytoplankton (which export less efficiently) are favored under the increasing nutrient stress. Thus, the projections of the NPP response to climate change are critically dependent on the simulated phytoplankton community structure, the efficiency of the biological pump and the resulting levels of regenerated production, which vary widely across the models. Community structure is represented simply in the CMIP5 models, and should be expanded to better capture the spatial patterns and climate-driven changes in export efficiency.

54 ENVIRONMENTAL SCIENCES↗

Improved Groundwater Table and L-Band Brightness Temperature Estimates for Northern Hemisphere Peatlands Using New Model Physics and SMOS Observations in a Global Data Assimilation Framework

There is an urgent need to include northern peatland hydrology in global Earth system models to better understand land-atmosphere interactions and sensitivities of peatland functions to climate change, and, ultimately, to improve climate change predictions. In this study, we introduced for the first time peatland-specific model physics into an assimilation scheme for L-band brightness temperature (Tb) data from the Soil Moisture Ocean Salinity (SMOS) mission to improve groundwater table estimates. We conducted two sets of model-only and data assimilation experiments using the Catchment Land Surface Model (CLSM), applying (over peatlands only) in one of them a peatland-specific adaptation (PEATCLSM). The evaluation against in-situ measurements of peatland groundwater table depth indicates the superiority of PEATCLSM model physics and additionally improved performance after assimilating SMOS Tb observations. The better performance of PEATCLSM over nearly all Northern Hemisphere peatlands is further supported by the better agreement between SMOS Tb observations and Tb estimates from the model-only and data assimilation runs. Within the data assimilation scheme, PEATCLSM reduces Tb observation-minus-forecast residuals and leads to reduced data assimilation updates of water storage components and, thus, reduced water budget imbalances in the assimilation system.

SMOS↗