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,405 records · Page 78

Leveraging High-resolution Molecular Composition of Soil Organic Matter to Enhance Carbon Cycling Modeling

Soils store more carbon than the atmosphere and vegetation combined, yet Earth system models still struggle to predict how this vast reservoir will respond to environmental change. A central limitation is that most soil biogeochemical models represent organic matter using bulk conceptual pools or chemically homogeneous fractions, preventing direct use of rapidly expanding molecular-scale datasets. Here we develop and test a new soil decomposition framework that explicitly integrates high-resolution information on organic matter composition. First, we construct a molecularly informed litter decomposition module in which plant inputs are partitioned into five functional compound classes—carbohydrates, proteins, lignin-like aromatics, lipids, and carbonyls—using a molecular mixing model calibrated to solid-state 13 C Nuclear Magnetic Resonance (NMR) spectra. Class-specific kinetics, lignin-dependent physical protection, and substrate-driven microbial carbon use efficiency allow the module to capture metabolic tradeoffs associated with enzyme production and nutrient limitation. We then embed this litter module within a microbially explicit whole-soil model that tracks the transformation of these compound classes through particulate organic matter, dissolved organic matter, mineral-associated organic matter, and microbial biomass. High-resolution Fourier Transform Ion Cyclotron Resonance mass spectrometry (FTICR-MS) data are used to link internal pools to measurable soil organic matter fractions and to constrain key process parameters. Applications at soil-core and ecosystem scales demonstrate that the new model reproduces observed soil respiration dynamics while providing mechanistic attribution of CO 2 fluxes to specific chemical classes and pools. Compared to existing frameworks such as the Community Land Model soil biogeochemistry module and the Millennial model, our approach maintains competitive predictive skill while substantially improving interpretability and opportunities for data–model integration. This work illustrates a viable pathway for leveraging molecular-scale observations to reduce structural uncertainty in soil carbon–climate feedback projections.

54 ENVIRONMENTAL SCIENCES↗

ELM model simulations of Plum Island Ecosystems LTER low marsh site 2018-2020

Model simulations using the Department of Energy's Energy Exascale Earth System Model (E3SM) land model (ELM) with improved capabilities to represent vegetation response to salinity and inundation. The simulations were conducted for a tidal salt marsh at Plum Island Ecosystems Long Term Ecological Research (LTER) site near Rowley, Massachusetts, USA; the site is a low marsh dominated by Spartina alterniflora. The model was forced with site-specific meteorology, salinity and tidal cycles from 2018-2020. Four sets of model simulations are included and described below:1. Parameterization of the salinity response function. These simulations tested different combinations of values for optimal salinity and salinity tolerance.2. Model evaluation. This comparison conducted simulations using the default model, the salinity function only, the submergence function only, and both the salinity and submergence functions. 3. Salinity scenarios. These simulations used the 2018 salinity input data varied by -5 to +10 ppt salinity.4. Water level scenarios. These simulations used the tide height varied by -10 to +50 cm. These simulations were used to demonstrate how the salinity and submergence functions better represent carbon uptake by tidal salt marshes.The data package includes netCDF files used as forcing files for tide height and salinity, one for each year 2018-2020 at observed salinity concentrations, and an additional three forcing files in which salinity concentrations were varied 5 ppt lower, 5 ppt higher, and 10 ppt higher than the measured 2018 time series. Also included are python scripts for creating forcing files, plain text parameter and command files for running simulations, model outputs in netCDF format, and python scripts for visualizing outputs. Code for the modified E3SM model is archived in Sulman et al 2023 at doi:10.15485/1991625. More detail about files is provided in the README.md file.

54 ENVIRONMENTAL SCIENCES↗

Forecasting the response of Earth's surface to future climatic and land use changes: A review of methods and research needs

In the future, Earth will be warmer, precipitation events will be more extreme, global mean sea level will rise, and many arid and semiarid regions will be drier. Human modifications of landscapes will also occur at an accelerated rate as developed areas increase in size and population density. We now have gridded global forecasts, being continually improved, of the climatic and land use changes (C&LUC) that are likely to occur in the coming decades. However, besides a few exceptions, consensus forecasts do not exist for how these C&LUC will likely impact Earth-surface processes and hazards. In some cases, we have the tools to forecast the geomorphic responses to likely future C&LUC. Fully exploiting these models and utilizing these tools will require close collaboration among Earth-surface scientists and Earth-system modelers. This paper assesses the state-of-the-art tools and data that are being used or could be used to forecast changes in the state of Earth's surface as a result of likely future C&LUC. We also propose strategies for filling key knowledge gaps, emphasizing where additional basic research and/or collaboration across disciplines are necessary. The main body of the paper addresses cross-cutting issues, including the importance of nonlinear/threshold-dominated interactions among topography, vegetation, and sediment transport, as well as the importance of alternate stable states and extreme, rare events for understanding and forecasting Earth-surface response to C&LUC. Five supplements delve into different scales or process zones (global-scale assessments and fluvial, aeolian, glacial/periglacial, and coastal process zones) in detail.

54 ENVIRONMENTAL SCIENCES↗

Challenging a Global Land Surface Model in a Local Socio-Environmental System

Land surface models (LSMs) predict how terrestrial fluxes of carbon, water, and energy change with abiotic drivers to inform the other components of Earth system models. Here, we focus on a single human-dominated watershed in southwestern Michigan, USA. We compare multiple processes in a commonly used LSM, the Community Land Model (CLM), to observational data at the single grid cell scale. For model inputs, we show correlations (Pearson’s R) ranging from 0.46 to 0.81 for annual temperature and precipitation, but a substantial mismatch between land cover distributions and their changes over time, with CLM correctly representing total agricultural area, but assuming large areas of natural grasslands where forests grow in reality. For CLM processes (outputs), seasonal changes in leaf area index (LAI; phenology) do not track satellite estimates well, and peak LAI in CLM is nearly double the satellite record (5.1 versus 2.8). Estimates of greenness and productivity, however, are more similar between CLM and observations. Summer soil moisture tracks in timing but not magnitude. Land surface reflectance (albedo) shows significant positive correlations in the winter, but not in the summer. Looking forward, key areas for model improvement include land cover distribution estimates, phenology algorithms, summertime radiative transfer modelling, and plant stress responses.

54 ENVIRONMENTAL SCIENCES↗

US energy system transitions under cumulative emissions budgets

Cumulative emissions budgets are increasingly being used by decision makers and analysts to understand emissions reductions and associated transitions in the context of long-term goals such as limiting global mean temperature increase over the century to 1.5 or 2°C (IPCC, 2014). While previous studies have explored the implications of cumulative emissions budgets for the global economy, few studies have conducted regional and national-level analyses (Iyer et al., 2015; Riahi et al., 2015). This paper explores cumulative emissions budgets through 2050 consistent with the 1.5 and 2°C long-term temperature goals in the context of the United States. We employ a state-level model of the United States embedded within a global human-Earth system model (GCAM-USA) to study the implications of such budgets for the U.S. energy system (Iyer et al., 2017a; Iyer et al., 2017b). Our results show that achieving the stringent cumulative emissions budgets entails accelerated deployment of energy conserving technology, almost complete decarbonization of the power sector and increased electrification of buildings and industrial end-use sectors and decarbonization of transport employing a combination of electrification and the substitution of bioenergy for fossil fuels by the mid-century. We also find substantial state-level differences in the relative roles of these decarbonization strategies. Furthermore, our results highlight that increased ambition in the near-term will be valuable in setting the stage for smoother transformations in the future to achieve stringent cumulative emissions budgets consistent with long-term temperature goals established by the international community.

Feijoo, Felipe↗

The Effects of an Improved Dynamic Vegetation Phenology Representation in a Global Land Surface Model

Evapotranspiration (ET) is a major driver of the interaction between the land surface and the atmosphere through its component mechanisms, including plant transpiration (T) and soil evaporation. To accurately capture land-atmosphere interactions in global Earth System Models, it is thus critical that the underlying land surface models accurately model both the land hydrology as well as the dynamic response of vegetation to environmental drivers. In an effort to introduce a more realistic vegetation representation, the NASA Catchment land surface model, which is part of the Goddard Earth Observing System (GEOS), has previously been merged with the carbon and nitrogen physics modules of the Community Land Model version 4, resulting in the new Catchment-CN model. Catchment-CN has inherited the advanced treatment of land surface hydrology of Catchment, but is additionally able to dynamically model the response of vegetation to environmental drivers, in contrast to the fixed vegetation climatology that was prescribed in Catchment. Recently, the parameterization of Catchment-CN vegetation has been augmented to better account for variability of vegetation responses to external forcings within existing plant functional types, and vegetation parameters have been calibrated against Moderate Resolution Imaging Spectrometer observations of the fraction of absorbed photosynthetically radiation. These efforts have led to a significant reduction in the RMSE of modeled photosynthetic activity with respect to observations.This presentation investigates the effect of the improved vegetation representation on the partitioning of ET within Catchment-CN. Specifically, we compare global maps of the T:ET ratio across different temporal scales in (1) the original Catchment model, (2) the original Catchment-CN model, and (3) the augmented and calibrated Catchment-CN model. The modeled T and ET estimates are compared against a comprehensive set of ground observations from various field studies, as well as independent global T:ET estimates from different ET algorithms provided in the context of the Water Cycle Observation Multi-mission Strategy ? Evapotranspiration (WACMOS-ET) initiative.

Kolassa, J.↗

Exploring a Variable-Resolution Approach for Simulating Regional Climate in the Rocky Mountain Region Using the VR-CESM: VR-CESM Simulation in the Rocky Mountains

The reliability of climate simulations and projections, particularly in the regions with complex terrains, is greatly limited by the model resolution. Here, we evaluate the variable-resolution Community Earth System Model (VR-CESM) with a high-resolution (0.125°) refinement over the Rocky Mountain region. The VR-CESM results are compared with observations, as well as CESM simulation at a quasi-uniform 1° resolution (UNIF) and Canadian Regional Climate Model version 5 (CRCM5) simulation at a 0.11° resolution. We find that VR-CESM is effective at capturing the observed spatial patterns of temperature, precipitation, and snowpack in the Rocky Mountains with the performance comparable to CRCM5, while UNIF is unable to do so. VR-CESM and CRCM5 simulate better the seasonal variations of precipitation than UNIF, although VR-CESM still overestimates winter precipitation whereas CRCM5 and UNIF underestimate it. All simulations distribute more winter precipitation along the windward (west) flanks of mountain ridges with the greatest overestimation in VR-CESM. VR-CESM simulates much greater snow water equivalent peaks than CRCM5 and UNIF, although the peaks are still 10–40% less than observations. Moreover, the frequency of heavy precipitation events (daily precipitation ≥ 25 mm) in VR-CESM and CRCM5 is comparable to observations, whereas the same events in UNIF are an order of magnitude less frequent. In addition, VR-CESM captures the observed occurrence frequency and seasonal variation of rain-on-snow days and performs better than UNIF and CRCM5. Lastly, these results demonstrate the VR-CESM's capability in regional climate modeling over the mountainous regions and its promising applications for climate change studies.

54 ENVIRONMENTAL SCIENCES↗

Impacts of climate change on future hurricane induced rainfall and flooding in a coastal watershed: A case study on Hurricane Harvey

The warming climate is likely to increase hurricane-associated extreme rainfall and lead to sea-level rise (SLR). Thus, how the floods induced by intense hurricanes respond to these potential changes is of great concern. This study investigates the future warmer climate impacts on hurricane-induced extreme rainfall, and—more importantly—the subsequent compound flooding at the watershed scale (from an event-based analysis perspective). To this goal, a modeling framework is designed based on the Distributed Hydrology Soil Vegetation Model (DHSVM), the Two-Dimensional fvand the Regional Community Earth System Model (R-CESM). The framework was applied to Hurricane Harvey (2017) at the Clear Creek watershed (a coastal watershed in the southern Houston) as a case study. The results show that the projected maximum rainfall totals over the watershed would be exacerbated by 17.7 % and 49.7 % in the 2050s and 2090s (respectively) under Representative Concentration Pathway 8.5 (RCP 8.5). This means a 16.1 % increase in Harvey rainfall over the watershed per degree Celsius increase in Mean Surface Temperature over the Gulf of Mexico region (18°~31° N, 77°~98° W). Meanwhile, the increases in maximum inundation extent would be 11.0 % (2050s) and 19.5 % (2090s). Furthermore, considerable increases in maximum inundation depth and duration in regions along the middle and downstream of Clear Creek (and also those around Clear Lake) are expected. The projected SLR will have little effect on the maximum inundation depth and extent if storm surge changes are not taken into account; meanwhile, it will influence the inundation duration at downstream locations. In conclusion, this modeling framework can be also applied at other coastal watersheds to evaluate the projected climate change impacts on the compound flooding induced by extreme climate events.

54 ENVIRONMENTAL SCIENCES↗

Different Spatiotemporal Patterns in Global Human Population and Built‐Up Land

Abstract Population concentration and built‐up land expansion are two prominent features of contemporary urbanization. Existing literature on the population aspect of urbanization has mostly focused on national and regional aggregates, and literature on the land development aspect has often relied on spatial case studies of individual cities or their meta‐analyses. Using newly‐available data, here we conduct the first global‐coverage, spatial analysis of the relationship between (changes in) population and built‐up land at multiple spatial scales, and compare to existing common beliefs about urbanization based on individual city studies. We find that population and built‐up land show distinctly different spatial and temporal patterns (with a global correlation coefficient around 0.6). Contrary to common impressions, our results show that during recent decades, developed and developing regions across the world experienced comparable amounts of built‐up land expansion. While meta‐analyses have reported that built‐up land in urban areas expands globally on average twice as fast as population grows, our results show the global change rates of built‐up land and population are similar. Also, most global population, including what national statistics agencies call urban population, reside in areas with low land development levels (which are frequently less than 5% built up). These changes in perspective suggest that urbanization's potential large‐scale impacts may need to be re‐evaluated, and lead to best‐practice recommendations for urbanization modeling and analysis. Especially, the common practice in large‐scale earth system modeling of assuming demographically‐defined urban population resides in areas with medium to high built‐up land development levels should change.

Gao, Jing↗

Pyrocumulonimbus Events Over British Columbia in 2017: An Ensemble Model Study of Parameter Sensitivities and Climate Impacts of Wildfire Smoke in the Stratosphere

Abstract The Pyrocumulonimbus (pyroCb) events over British Columbia in 2017 were observed in the lower stratosphere for about 8–10 months after the smoke injections. Several previous studies used global climate models to investigate the physical parameters for the 2017 pyroCb events, but the conclusions show strong model dependency. In this study, we use the Energy Exascale Earth System Model (E3SM) and complete an ensemble of runs exploring three injection parameters: smoke aerosol mass, the percentage of black carbon within the smoke aerosols, and plume injection height. Additionally, we consider the heterogeneous reaction of ozone and primary organic matter. According to the satellite daily observed aerosol optical depth (AOD), we find that the best ensemble member is the simulation with 0.4 Tg of smoke, 3% of which is black carbon, a 13.5 km smoke injection height, and a 10 −5 probability factor of the heterogeneous reaction. Besides AOD, we examine the ensemble score based on the metrics used in previous studies: the metrics of the extinction coefficient at 18 km altitude and the maximum plume height. The conclusion of the best estimate of the injection parameters for 2017 pyroCb events shows strong not only model but also evaluation metric dependency. We use the Random Forest machine learning technique to quantify the relative importance of each parameter in accurately simulating the 2017 pyroCb events and find that the injection height is the most critical feature, no matter which metric is used to score the ensemble members.

54 ENVIRONMENTAL SCIENCES↗

Using Satellite and ARM Observations to Evaluate Cold Air Outbreak Cloud Transitions in E3SM Global Storm‐Resolving Simulations

Abstract This study examines marine boundary layer cloud regime transition during a cold air outbreak (CAO) over the Norwegian Sea, simulated by a global storm‐resolving model (GSRM) known as the Simple Cloud‐Resolving Energy Exascale Earth System Model Atmosphere Model (SCREAM). By selecting observational references based on a combination of large‐scale conditions rather than strict time‐matched comparisons, this study finds that SCREAM qualitatively captures the CAO cloud transition, including boundary layer growth, cloud mesoscale structure, and phase partitioning. SCREAM also accurately locates the greatest ice and liquid in the mesoscale updrafts, however, underestimates supercooled liquid water in cumulus clouds. The model evaluation approach adopted by this study takes advantages of the existing computational‐expensive global simulations of GSRM and the available observations to understand model performance and can be applied to assessments of other cloud regimes in different regions. Such practice provides valuable guidance on the future effort to correct and improve biased model behaviors.

54 ENVIRONMENTAL SCIENCES↗

Exploring the use of novel programming models in land surface models

A wide range of programming models are currently under rapid development to meet the needs of application developers looking to work on more complex machines. These models fill a variety of roles. Some look to abstract supercomputer architecture, including both processors and memory, to present a strategy for portable performance across a wide range of machines. Others look to expose concurrency by explicitly constructing task-driven dependency graphs that allow a scheduler to find parallelism. Here we explore the implications for application codes of adopting two such programming models, Kokkos and Legion, one from each class of models. We specifically focus on the software design implications on refactoring existing applications, rather than the performance and performance tuning of these models. We identify a strategy for refactoring the Energy Exascale Earth System Model's Land Surface Model, an extremely complex code for climate applications, and prototype a series of mini-apps that explore the adoption of Kokkos and Legion. In doing this, we identify commonalities across the models, leading to a series of conclusions about application software design and refactoring for the adoption of novel programming models.

Coon, Ethan↗

From the Arctic to the tropics: multi-biome prediction of leaf mass per area using leaf reflectance

Leaf-mass-per-area (LMA) is a key plant trait, reflecting tradeoffs between leaf photosynthetic function, longevity, and structural investment. Capturing spatial and temporal variability in LMA has been a long-standing goal of ecological research and is an essential component for advancing Earth system models. Despite the substantial variation in LMA within and across Earth's biomes, an efficient, globally generalizable approach to predict LMA is still lacking. In this work, we explored the capacity to predict LMA from leaf spectra across much of the global LMA trait space, with values ranging from 17 g m -2 to 393 g m -2 . Our dataset contained leaves from a wide range of biomes from the high Arctic to the tropics, included broad- and needle-leaf species, and upper and lower canopy (i.e. sun and shade) growth environments. Here we demonstrate the capacity to rapidly estimate LMA using only spectral measurements across a wide range of species, leaf age, and canopy position, from diverse biomes. Our model captures LMA variability with high accuracy and low error (R 2 =0.89; RMSE=15.45 g m -2 ). Lastly, our finding highlights that the leaf economic spectrum is mirrored by the leaf optical spectrum paving the way for this technology to predict the diversity of LMA in ecosystems across global biomes.

54 ENVIRONMENTAL SCIENCES↗

Bayesian Calibration of the Community Land Model using Surrogates

We present results from the Bayesian calibration of hydrological parameters of the Community Land Model (CLM), which is often used in climate simulations and Earth system models. A statistical inverse problem is formulated for three hydrological parameters, conditioned on observations of latent heat surface fluxes over 48 months. Our calibration method uses polynomial and Gaussian process surrogates of the CLM, and solves the parameter estimation problem using a Markov chain Monte Carlo sampler. Posterior probability densities for the parameters are developed for two sites with different soil and vegetation covers. Our method also allows us to examine the structural error in CLM under two error models. We find that accurate surrogate models can be created for CLM in most cases. The posterior distributions lead to better prediction than the default parameter values in CLM. Climatologically averaging the observations does not modify the parameters’ distributions significantly. The structural error model reveals a correlation time-scale which can potentially be used to identify physical processes that could be contributing to it. While the calibrated CLM has a higher predictive skill, the calibration is under-dispersive.

Bayesian calibration↗

Bayesian calibration of the Community Land Model using surrogates

We present results from the Bayesian calibration of hydrological parameters of the Community Land Model (CLM), which is often used in climate simulations and Earth system models. A statistical inverse problem is formulated for three hydrological parameters, conditional on observations of latent heat surface fluxes over 48 months. Our calibration method uses polynomial and Gaussian process surrogates of the CLM, and solves the parameter estimation problem using a Markov chain Monte Carlo sampler. Posterior probability densities for the parameters are developed for two sites with different soil and vegetation covers. Our method also allows us to examine the structural error in CLM under two error models. We find that surrogate models can be created for CLM in most cases. The posterior distributions are more predictive than the default parameter values in CLM. Climatologically averaging the observations does not modify the parameters' distributions significantly. The structural error model reveals a correlation time-scale which can be used to identify the physical process that could be contributing to it. While the calibrated CLM has a higher predictive skill, the calibration is under-dispersive.

54 ENVIRONMENTAL SCIENCES↗

Applying Corrective Machine Learning in the E3SM Atmosphere Model in C++ (EAMxx)

The Simplified Cloud-Resolving E3SM Atmosphere Model (SCREAM) is the newest addition to the family of Earth System Models capable of explicitly resolving convective systems. SCREAM is a kilometer-scale configuration of the advanced E3SM Atmosphere Model (EAMxx), designed for heterogeneous systems. While the enhanced accuracy of kilometer-scale modeling offers significant benefits, it comes with a substantial computational cost, limiting feasible simulation durations to only a few years, even on the fastest supercomputers. Machine learning presents an opportunity for scientists to achieve the high accuracy of storm-resolving models at a significantly reduced cost. Building on the previous success of applying corrective machine learning (ML) to the FV3 model, this study explores the effects of implementing corrective ML in EAMxx-SCREAM. We also address the computational challenges of integrating the corrective ML, which is written in Python, with the C++/Kokkos EAMxx driver, as well as the potential pitfalls of generalizing an approach that was effective with one atmosphere model to another.

54 ENVIRONMENTAL SCIENCES↗

National Climate Database (NCDB)

The National Climate Database (NCDB) is a high resolution, bias-corrected climate dataset consisting of the three most widely used variables of solar radiation- global horizontal (GHI), direct normal (DNI), and diffuse horizontal irradiance (DHI)- as well as other meteorological data. The goal of the NCDB is to provide unbiased high temporal and spatial resolution climate data needed for renewable energy modeling. The NCDB is modeled using a statistical downscaling approach with Regional Climate Model (RCM)-based climate projections obtained from the North American Coordinated Regional Climate Downscaling Experiment (NA-CORDEX; linked below). Daily climate projections simulated by the Canadian Regional Climate Model 4 (CanRCM4) forced by the second-generation Canadian Earth System Model (CanESM2) for two Representative Concentration Pathways (RCP4.5 or moderate emissions scenario and RCP8.5 or highest baseline emission scenario) are selected as inputs to the statistical downscaling models. The National Solar Radiation Database (NSRDB) is used to build and calibrate statistical models.

Array↗

An extensible perturbed parameter ensemble for the Community Atmosphere Model version 6

This paper documents the methodology and preliminary results from a perturbed parameter ensemble (PPE) technique, where multiple parameters are varied simultaneously and the parameter values are determined with Latin hypercube sampling. This is done with the Community Atmosphere Model version 6 (CAM6), the atmospheric component of the Community Earth System Model version 2 (CESM2). We apply the PPE method to CESM2–CAM6 to understand climate sensitivity to atmospheric physics parameters. The initial simulations vary 45 parameters in the microphysics, convection, turbulence and aerosol schemes with 263 ensemble members. These atmospheric parameters are typically the most uncertain in many climate models. Control simulations and targeted simulations to understand climate forcing due to aerosols and fast climate feedbacks are analyzed. The use of various emulators is explored in the multi-dimensional space mapping input parameters to output metrics. Parameter impacts on various model outputs, such as radiation, cloud and aerosol properties, are evaluated. Machine learning is also used to probe optimal parameter values against observations. Our findings show that PPE is a valuable tool for climate uncertainty analysis. Furthermore, by varying many parameters simultaneously, we find that many different combinations of parameter values can produce results consistent with observations, and thus careful analysis of tuning is important. The CESM2–CAM6 PPE is publicly available and extensible to other configurations to address questions of other model processes in the atmosphere and other model components (e.g., coupling to the land surface).

54 ENVIRONMENTAL SCIENCES↗