SEARCH · Engineering Papers
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.
We review the most current knowledge pertaining to carbon cycle components in an Earth System model and discuss data assimilation applications to inform and improve process-based representations in these models.
Abstract The sea ice component of the Community Earth System Model version 2 (CESM2) contains new “mushy‐layer” physics that simulates prognostic salinity in the sea ice, with consequent modifications to sea ice thermodynamics and the treatment of melt ponds. The changes to the sea ice model and their influence on coupled model simulations are described here. Two simulations were performed to assess the changes in the vertical thermodynamics formulation with prognostic salinity compared to a constant salinity profile. Inclusion of the mushy layer thermodynamics of Turner et al. (2013, https://doi.org/10.1002/jgrc.20171 ) in a fully coupled Earth system model produces thicker and more extensive sea ice in the Arctic, with relatively unchanged sea ice in the Antarctic compared to simulations using a constant salinity profile. While this is consistent with the findings of uncoupled ice‐ocean model studies, the role of the frazil and congelation growth is more important in fully coupled simulations. Melt pond drainage is also an important contribution to simulated ice thickness differences as also found in the uncoupled simulations of Turner and Hunke (2015; https://doi.org/10.1002/2014JC010358 ). However, it is an interaction of the ponds and the snow fraction that impacts the surface albedo and hence the top melt. The changes in the thermodynamics and resulting ice state modify the ice‐ocean‐atmosphere fluxes with impacts on the atmosphere and ocean states, particularly temperature.
Understanding terrestrial carbon cycling has relied primarily on studies of topsoils that are typically characterized to depths shallower than 0.5 m. At a semiarid site instrumented down to 7 m, here we measured seasonal- and depth-resolved carbon inventories and fluxes and groundwater and unsaturated zone flow rates. Measurements showed that ~30% of the CO 2 efflux to the atmosphere (60% in winter) originates from below 1 m, contrary to predictions of less than 1% by Earth System Model land modules. Respiration from deeper roots and deeper microbial communities is supported by favorable subsurface temperatures, moisture, and oxygen availability. Below 1 m, dissolved organic carbon fluxes from the overlying soil and C from deep roots and exudates are expected to be important in sustaining microbial respiration. Finally, because these conditions are characteristic of semiarid climate regions, we contend that Earth System Model land modules should incorporate such deeper soil processes to improve CO 2 flux predictions.
Abstract Composted manure and green waste amendments have been shown to increase net carbon (C) sequestration in rangeland soils and have been proposed as a means to help lower atmospheric CO 2 concentrations. However, the effect of climate change on soil organic C (SOC) stocks and greenhouse gas emissions in rangelands is not well understood, and the viability of climate change mitigation strategies under future conditions is even less certain. We used a process‐based biogeochemical model (DayCent) at a daily time step to explore the long‐term effects of potential future climate changes on C and greenhouse gas dynamics in annual grassland ecosystems. We then used the model to explore how the same ecosystems might respond to climate change following compost amendments to soils and determined the long‐term viability of net SOC sequestration under changing climates. We simulated net primary productivity (NPP), SOC, and greenhouse gas fluxes across seven California annual grasslands with and without compost amendments. We drove the DayCent simulations with field data and with site‐specific daily climate data from two Earth system models (CanESM2 and HadGEM‐ES) and two representative concentration pathways (RCP4.5 and RCP8.5) through 2100. NPP and SOC stocks in unamended and amended ecosystems were surprisingly insensitive to projected climate changes. A one‐time amendment of compost to rangeland acted as a slow‐release organic fertilizer and increased NPP by up to 390–814 kg C ha −1 year −1 across sites. The amendment effect on NPP was not sensitive to Earth system model or emissions scenario and endured through the end of the century. Net SOC sequestration amounted to 1.96 ± 0.02 Mg C ha −1 relative to unamended soils at the maximum amendment effect. Averaged across sites and scenarios, SOC sequestration peaked 22 ± 1 years after amendment and declined but remained positive throughout the century. Though compost stimulated nitrous oxide (N 2 O) emissions, the cumulative net emissions (in CO 2 equivalents) due to compost were far less than the amount of SOC sequestered. Compost amendments resulted in a net climate benefit of 69.6 ± 0.5 Tg CO 2 e 20 ± 1 years after amendment if applied to similar ecosystems across the state, amounting to 39% of California's rangeland. These results suggest that the biogeochemical benefits of a single amendment of compost to rangelands in California are insensitive to climate change and could contribute to decadal‐scale climate change mitigation goals alongside emissions reductions.
The absorption of nearly half of all anthropogenic carbon dioxide (CO2) emissions by terrestrial and marine ecosystems has played a critical role in mitigating climate change. However, a persistent lack of understanding of these uptake processes impedes attribution, and thereby introduces large uncertainties into projections of the trajectory of carbon and climate in coming decades. A better understanding and quantification of the seasonal and interannual variability of carbon flux between the atmosphere and ocean and terrestrial biospheres is critically needed to improve predictions from Earth system models. We propose to bring together a diverse suite of remote sensing observations, mechanistic models, and NASA's Goddard Earth Observing System, version 5 (GEOS-5) model to address two objective 1) quantify the processes controlling the temporal variability of biosphere-atmosphere CO2 flux on seasonal and interannual timescales from 1989 to 2018 and 2) evaluating the predictability of the components of carbon flux on seasonal timescales. Because GEOS-5 has been developed to incorporate a variety of Earth system observations and includes a complete and physically consistent depiction of the atmosphere, land, and ocean carbon cycles, it represents a unique and ideal synthesis framework for this effort. By expanding the realism of processes that could be included in GEOS-5 seasonal forecasts, this effort supports NASA's goals toward an integrated Earth system modeling and prediction framework.We propose to integrate a number of observations into the GEOS-5 modeling system to refine retrospective estimates of both terrestrial and marine carbon flux and to test hypotheses regarding the effects of climate drivers on terrestrial carbon flux. This will include assimilating several satellite ocean color products in the NASA Ocean Biogeochemistry Model (NOBM) to improve the representation of ocean biology. We will also utilize remote sensing datasets to increase the temporal resolution of global land-use change estimates and fire emissions to better constrain the temporal variability of the terrestrial carbon cycle. This information will be incorporated into two NASA terrestrial biosphere models Catchment-CN, a dynamic vegetation model that is integrated within GEOS-5 and the Ecosystem Demography (ED) model, a core modeling resource supporting the Global Ecosystem Dynamics Investigation (GEDI) mission to map biomass. Fluxes will be evaluated through 1) comparison to top-down flux estimates inferred from surface CO2 observations and 2) comparison against available aircraft and satellite observations.By combining predictions of net ecosystem exchange and air-sea CO2 flux from land and ocean models with statistical predictions of land use change, fire, and fossil fuel emissions, we propose to create the first dynamical CO2 forecast on seasonal timescales. Because of the experimental nature of this endeavor, an important project objective will be to evaluate the limits to the predictability of different carbon cycle processes over lead times of 1-12 months using different satellite and surface data streams.
Abstract Northern peatlands are a terrestrial carbon store, with an annual sink of 0.1 Pg C yr −1 and a total storage estimate of 547 Pg C. Northern peatlands are also major contributors of atmospheric methane. Most land surface models do not accurately represent peatland carbon emissions, partly because they do not represent the hydrologic cycle and/or microtopography adequately. Interactions between water table depth and microtopography in peatlands influence decomposition and modulate CO 2 and CH 4 fluxes. A modified version of the land surface component of the Energy Exascale Earth System Model, was recently created to represent the microtopography and hydrology of a raised dome bog in northern Minnesota, USA. In this study, three microtopographic parameters are analyzed in the modified version: hummock height, hummock‐hollow spacing, and percent hollow. Terrestrial laser scanning observations are used to set uncertainty bounds for these parameters. Our model experiment results suggest that carbon‐related quantities of interest (QOI) were typically the most sensitive to hummock height, and those QOI (especially net ecosystem exchange, NEE) were sensitive to interactions between parameters. Furthermore, NEE was most relatively influenced by microtopographic parameters in the model, varying by 35%. We found that increasing hummock height resulted in more C being stored in plant tissue and less in soil organic matter. This coincided with decreases in Sphagnum and increases in Picea and shrub net primary production. These results suggest that future studies may consider extending prognostic capabilities of carbon cycling by incorporating hummock hollow microtopography into earth system models.
The parameterization of key photosynthesis parameters is one of the key uncertain sources in modeling ecosystem gross primary productivity (GPP). Solar-induced chlorophyll fluorescence (SIF) offers a good proxy for GPP since it marks the actual process of photosynthesis; while machine learning (ML) provides a robust approach to model the GPP-SIF relationship. Here, we trained the boosted regressing tree (BRT) and the Random Forest ML models with Greenhouse Gases Observing Satellite SIF data and in situ GPP observations from 49 eddy covariance towers. These trained ML GPP-SIF models were fed into the Energy Exascale Earth System Model (E3SM) Land Model (ELM) to generate ELM-simulated global SIF estimates, which were then benchmarked against satellite SIF observations with a surrogate modeling approach. Our results indicated good modeling performance of the ML-based GPP-SIF relationship. The ELM model when fed with the ML GPP-SIF models also can well predict the spatial-temporal variations in SIF. We also found high model accuracy for the surrogate modeling. Model parameter sensitivity analysis suggested that the fraction of leaf nitrogen in RuBisCO (flnr) is the most sensitive parameter to the SIF; other sensitive parameters include the Ball-Berry stomatal conductance slope (mbbopt) and the vcmax entropy (vcmaxse). The posterior uncertainty in simulated GPP was greatly reduced after benchmarking, and the model produced improved spatial patterns of mean GPP relative to FLUXCOM GPP. Our integrated approach provides a new avenue for improving land models and using remote-sensing SIF, which can be further improved in the future with more ground- and satellite-based observations.
Climate change is expected to increase the frequency of flooding events and, thus, the risks of flood-related mortality and infrastructure damage. Global-scale assessments of future flooding from Earth system models based only on precipitation changes neglect important processes that occur within the land surface, particularly plant physiological responses to rising CO 2 . Higher CO 2 can reduce stomatal conductance and transpiration, which may lead to increased soil moisture and runoff in some regions, promoting flooding even without changes in precipitation. Here we assess the relative impacts of plant physiological and radiative greenhouse effects on changes in daily runoff intensity over tropical continents using the Community Earth System Model. We find that extreme percentile rates increase significantly more than mean runoff in response to higher CO 2 . Plant physiological effects have a small impact on precipitation intensity but are a dominant driver of runoff intensification, contributing to one half of the 99th and one third of the 99.9th percentile runoff intensity changes.
The Coupled Model Intercomparison Project (CMIP) is a flagship of the World Climate Research Programme (WCRP). CMIP has become a recognised ‘brand’ in climate circles evolving over the last thirty years from a targeted research activity by a small number of climate modelling centres intercomparing their Earth System Model (ESM) simulations to a broad international coordinated research effort (Durack et al, 2025). CMIP is organized as a research activity leveraging funded and in-kind contributions from experts within modelling centres and the broader scientific community supported more recently by a fully-funded International Project Office. Within CMIP, Model Intercomparison Projects (MIPs) are community-designed to understand past, present and future climate. CMIP data provides a valuable resource for climate research and is routinely used to assess model representation of climate processes and test scientific hypotheses in the context of model uncertainty and (forced and internal) variability as evident from its prolific use in scientific publications1 . The impact relies on enabling infrastructure (most prominently via the Earth System Grid Federation (ESGF)), which allows sharing of simulation output, provision of the boundary conditions used in each simulation, and definition of the data standards that are essential to facilitating wide use of the data. The impact is supplemented by the wide-ranging scrutiny to which model simulations are subjected. Beyond its use in research, CMIP data is a key resource for communities producing derived climate information from downscaling and impact studies, such as the Coordinated Regional Downscaling Experiment (CORDEX; Gutowski et al., 2016) and the Intersectoral Impacts MIP (ISIMIP; Frieler et al., 2024). Government, academic and commercial entities also increasingly rely on CMIP and its downstream data for climate risk assessments and climate services (for example, Copernicus Climate Change Service and World Bank portal). This means that, although CMIP is a research activity, it increasingly serves a secondary and very relevant role as a provider of climate data – a long-recognised dichotomy (Stevens, 2024). Research and applications have distinct needs, with the former requiring flexibility and generality and the latter consistency. Here we explain how the design of the research activity has been adapted to reduce the burdens imposed by applications and how the research infrastructure might evolve to further enable scientific inquiry. We propose one possible approach to consistently providing model information and projections for applications in the future.
Plant water-use efficiency (WUE), which is the ratio of the uptake of carbon dioxide through photosynthesis to the loss of water through transpiration, is a very useful metric of the functioning of the land biosphere. WUE is expected to increase with atmospheric CO 2 , but to decline with increasing atmospheric evaporative demand – which can arise from increases in near-surface temperature or decreases in relative humidity. We have used Δ 13 C measurements from tree rings, along with eddy covariance measurements from Fluxnet sites, to estimate the sensitivities of WUE to changes in CO 2 and atmospheric humidity deficit. This enables us to reconstruct fractional changes in WUE, based on changes in atmospheric climate and CO 2 , for the entire period of the instrumental global climate record. We estimate that overall WUE increased from 1900 to 2010 by 48 ± 22 %, which is more than double that simulated by the latest Earth System Models. This long-term trend is largely driven by increases in CO 2 , but significant inter-annual variability and regional differences are evident due to variations in temperature and relative humidity. Here, there are several highly populated regions, such as western Europe and East Asia, where the rate of increase of WUE has declined sharply in the last 2 decades. Our data-based analysis indicates increases in WUE that typically exceed those simulated by Earth System Models – implying that these models are either underestimating increases in photosynthesis or underestimating reductions in transpiration.
Earth System Model (ESM) development benefits from attention to the downstream applications for nature and society. The evaluation of ESMs often focuses on a limited set of climatic impact-drivers (CIDs)1 and numerous large-scale parameters that are not directly relevant for impact analysis. To address this limitation, we propose a comprehensive framework for impact-relevant model version comparison. This framework explores both the intended and unintended consequences of shifts in model components, parameters, and resolutions, emphasizing the perspectives of the climate impacts community. We have selected a core set of CIDs as a basis for comparing model versions.
The Functionally-Assembled Terrestrial Ecosystem Simulator (FATES) is a vegetation model for use in Earth system models (ESMs). The model includes a size- and age-structured representation of tree dynamics, competition between functionally diverse plant functional types, and the biophysics underpinning plant growth, competition, mortality, as well as the carbon, water, and energy exchange with the atmosphere. The FATES model is designed as a modular vegetation model that can be integrated within a host land model for inclusion in ESMs. The model is designed for use in global change studies to understand and project the responses and feedbacks between terrestrial ecosystems and the Earth system under changing climate and other forcings.
The Functionally-Assembled Terrestrial Ecosystem Simulator (FATES) is a vegetation model for use in Earth system models (ESMs). The model includes a size- and age-structured representation of tree dynamics, competition between functionally diverse plant functional types, and the biophysics underpinning plant growth, competition, mortality, as well as the carbon, water, and energy exchange with the atmosphere. The FATES model is designed as a modular vegetation model that can be integrated within a host land model for inclusion in ESMs. The model is designed for use in global change studies to understand and project the responses and feedbacks between terrestrial ecosystems and the Earth system under changing climate and other forcings.
The EC-Earth climate model is a seamless Earth System Model (ESM) used to carry out climate research in 24 academic institutions and meteorological services from 11 countries in Europe. This model couples several components and it is continuously under development. In this work we present a study regarding the impact of the I/O and data management when using EC-Earth in well-known supercomputing environments. Most large-scale and long-term climate simulators have been developed bearing in mind the paramount importance of its scalability. However, the computational capabilities of the High Performance Computing (HPC) environments increase at so great speed that it is almost impossible to re-implement the whole models so that they are able to exploit efficiently the new features. Therefore, it is necessary to design different strategies to take advantage of them. In this work we present an operational framework to run ensemble simulations in HPC platforms. A set of experiments are presented in order to validate the suitability of this technique. Moreover, the derived impact regarding the I/O and data management aspects is analyzed.
We report that earth system models are complex and represent a large number of processes, resulting in a persistent spread across climate projections for a given future scenario. Owing to different model performances against observations and the lack of independence among models, there is now evidence that giving equal weight to each available model projection is suboptimal. This Perspective discusses newly developed tools that facilitate a more rapid and comprehensive evaluation of model simulations with observations, process-based emergent constraints that are a promising way to focus evaluation on the observations most relevant to climate projections, and advanced methods for model weighting. Lastly, these approaches are needed to distill the most credible information on regional climate changes, impacts, and risks for stakeholders and policy-makers.
Genomes reveal nutrient stress patterns Within the surface ocean, nitrogen, iron, and phosphorous can all be limiting nutrients for phytoplankton depending on location. Ustick et al. used the prevalence of Prochlorococcus genes involved in nutrient acquisition to develop maps of inferred nutrient stress across the global ocean (see the Perspective by Coleman). They found broad patterns of limitation consistent with an Earth system model and nutrient addition experiments. Leveraging metagenomic data in this manner is an appealing approach that will help to expand our understanding of the biogeochemistry in the vast open ocean. Science , this issue p. 287 ; see also p. 239
The projected loss of soil carbon to the atmosphere resulting from climate change is a potentially large but highly uncertain feedback to warming. The magnitude of this feedback is poorly constrained by observations and theory, and is disparately represented in Earth system models (ESMs). In this study, to assess the climatological temperature sensitivity of soil carbon, we calculate apparent soil carbon turnover times that reflect long-term and broad-scale rates of decomposition. We show that the climatological temperature control on carbon turnover in the top metre of global soils is more sensitive in cold climates than in warm climates and argue that it is critical to capture this emergent ecosystem property in global-scale models. We present a simplified model that explains the observed high cold-climate sensitivity using only the physical scaling of soil freeze-thaw state across climate gradients. Current ESMs fail to capture this pattern, except in an ESM that explicitly resolves vertical gradients in soil climate and carbon turnover. An observed weak tropical temperature sensitivity emerges in a different model that explicitly resolves mineralogical control on decomposition. Lastly, these results support projections of strong carbon-climate feedbacks from northern soils and demonstrate a method for ESMs to capture this emergent behaviour.
Riverine dissolved organic carbon (DOC) plays a vital role in regional and global carbon cycles. However, the processes of DOC conversion from soil organic carbon (SOC) and leaching into rivers are insufficiently understood, inconsistently represented, and poorly parameterized, particularly in land surface and Earth system models. As a first attempt to fill this gap, we propose a generic formula that directly connects SOC concentration with DOC concentration in headwater streams, where a single parameter, the transformation rate from SOC in the soil to DOC leaching flux (P r ), accounts for the overall processes governing SOC conversion to DOC and leaching from soils (along with runoff) into headwater streams. We then derive high-resolution P r maps over the contiguous US (CONUS) using SOC data from two different sources: the Harmonized World Soil Database v1.2 (HWSD) and SoilGrids 2.0. Both maps are developed following the same five major steps: (1) selecting independent catchments where observed riverine DOC data are available with reasonable quality; (2) estimating catchment-average SOC for the independent catchments; (3) estimating the P r values for these catchments based on the generic formula and catchment-average SOC; (4) developing a predictive model of P r with machine learning (ML) techniques and catchment-scale climate, hydrology, geology, and other attributes; and (5) deriving a national map of P r based on the ML model. For evaluation, we compare the DOC concentration derived using the P r map and the observed DOC concentration values at evaluation catchments. The resulting mean absolute scaled error and coefficient of determination are 0.73 and 0.47 for the HWSD-based model and 0.58 and 0.72 for the SoilGrids-based model, respectively, suggesting the effectiveness of the overall methodology. Efforts to constrain uncertainty and evaluate sensitivity of P r to different factors are discussed. To illustrate the use of such maps, we derive a riverine DOC concentration reanalysis dataset over CONUS. The two P r maps, robustly derived and empirically validated, lay a critical cornerstone for better simulating the terrestrial carbon cycle in land surface and Earth system models. Our findings not only set a foundation for improving our predictive understanding of the terrestrial carbon cycle at the regional and global scales, but also hold promises for informing policy decisions related to decarbonization and climate change mitigation. The data presented in this study are publicly available at https://doi.org/10.5281/zenodo.14563816 (Li et al., 2024).