SEARCH · Engineering Papers

Results for “Earth Model”

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.

SourceResearch connectionFrom yearThrough year

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

At least 55 records · Page 3

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Disentangling atmospheric, hydrological, and coupling uncertainties in compound flood modeling within a coupled Earth system model

Compound riverine and coastal flooding is usually driven by complex interactions among meteorological, hydrological, and ocean extremes. However, existing efforts to model this phenomenon often do not integrate hydrological processes across atmosphere–land–river–ocean systems, leading to substantial uncertainties that have not been fully examined. To bridge this gap, we leverage the new capabilities of the Energy Exascale Earth System Model (E3SM) that enable a multi-component framework that integrates coastally refined atmospheric, terrestrial, and oceanic components. We evaluate compound uncertainties arising from two-way land–river–ocean coupling in E3SM and track the cascading meteorological and hydrological uncertainties through ensemble simulations over the Delaware River basin and estuary during Hurricane Irene (2011). Our findings highlight the importance of two-way river–ocean coupling to compound flood modeling and demonstrate E3SM's capability in capturing compound flood extent near the coast, with a hit rate over 0.75. Our study shows the growing uncertainties that transition from atmospheric forcings to flood distribution and severity. Furthermore, an analysis based on artificial neural networks is used to assess the roles of hydrological drivers, such as infiltration and soil moisture, in the generation of compound flooding. The response of compound floods to tropical cyclones (TCs) is found to be susceptible to these often overlooked drivers. For instance, the flooded area could increase more than 2-fold (∼2.4) if Hurricane Irene were preceded by an extreme antecedent soil moisture condition (AMC). The results not only support the use of a multi-component framework for interactive flooding processes, but also underscore the necessity of broader definitions of compound flooding that encompass the simultaneous occurrence of intense precipitation, storm surge, and high AMC during TCs.

Feng, Dongyu [Pacific Northwest National Laborator↗

DOE OSTI2020

An improved mechanistic model for ammonia volatilization in Earth system models: Flow of Agricultural Nitrogen version 2 (FANv2)

Abstract. Volatilization of ammonia (NH 3 ) from fertilizers and livestock wastes forms a significant pathway of nitrogen losses in agricultural ecosystems and constitutes the largest source of atmospheric emissions of NH 3 . This paper describes a major update to the process model FAN (Flow of Agricultural Nitrogen), which evaluates NH 3 emissions interactively within an Earth system model; in this work, the Community Earth System Model (CESM) is used. The updated version (FANv2) includes a more detailed treatment of both physical and agricultural processes, which allows the model to differentiate between the volatilization losses from animal housings, manure storage, grazed pastures, and the application of manure and different types of mineral fertilizers. The modeled ammonia emissions are first evaluated at a local scale against experimental data for various types of fertilizers and manure, and they are subsequently run globally to evaluate NH3 emissions for 2010–2015 based on gridded datasets of fertilizer use and livestock populations. Comparison of regional emissions shows that FANv2 agrees with previous inventories for North America and Europe and is within the range of previous inventories for China. However, due to higher NH 3 emissions in Africa, India, and Latin America, the global emissions simulated by FANv2 (48 Tg N) are 30 %–40 % higher than in the existing inventories.

Vira, Julius↗

DOE OSTI2022

More Realistic Intermediate Depth Dry Firn Densification in the Energy Exascale Earth System Model (E3SM)

Abstract Earth system models account for seasonal snow cover, but many do not accommodate the deeper snowpack on ice sheets (aka firn) that slowly transforms to ice under accumulating snowfall. To accommodate and resolve firn depths of up to 60 m in the Energy Exascale Earth System Model's land surface model (ELM), we add 11 layers to its snowpack and evaluate three dry snow compaction equations in multi‐century simulations. After comparing results from ELM simulations (forced with atmospheric reanalysis) with empirical data, we find that implementing into ELM a two‐stage firn densification model produces more accurate dry firn densities at intermediate depths of 20–60 m. Compared to modeling firn using the equations in the (12 layer) Community Land Model (version 5), switching to the two‐stage firn densification model (with 16 layers) significantly decreases root‐mean‐square errors in upper 60 m dry firn densities by an average of 41 kg m −3 (31%). Simulations with three different firn density parameterizations show that the two‐stage firn densification model should be used for applications that prioritize accurate upper 60 m firn air content (FAC) in regions where the mean annual surface temperature is greater than roughly −31°C. Because snow metamorphism, firn density, and FAC are major components in modeling ice sheet surface albedo, melt water retention, and climatic mass balance, these developments advance broader efforts to simulate the response of land ice to atmospheric forcing in Earth system models.

58 GEOSCIENCES↗

DOE OSTI2023

Ensemble Spread Behavior in Coupled Climate Models: Insights From the Energy Exascale Earth System Model Version 1 Large Ensemble

AbstractAssessing uncertainty in future climate projections requires understanding both internal climate variability and external forcing. For this reason, single‐model initial condition large ensembles (SMILEs) run with Earth System Models (ESMs) have recently become popular. Here we present a new 20‐member SMILE with the Energy Exascale Earth System Model version 1 (E3SMv1‐LE), which uses a “macro” initialization strategy choosing coupled atmosphere/ocean states based on inter‐basin contrasts in ocean heat content (OHC). The E3SMv1‐LE simulates tropical climate variability well, albeit with a muted warming trend over the twentieth century due to overly strong aerosol forcing. The E3SMv1‐LE's initial climate spread is comparable to other (larger) SMILEs, suggesting that maximizing inter‐basin ocean heat contrasts may be an efficient method of generating ensemble spread. We also compare different ensemble spread across multiple SMILEs, using surface air temperature and OHC. The Community Earth system Model version 1, the only ensemble which utilizes a “micro” initialization approach perturbing only atmospheric initial conditions, yields lower spread in the first ∼30 years. The E3SMv1‐LE exhibits a relatively large spread, with some evidence for anthropogenic forcing influencing spread in the late twentieth century. However, systematic effects of differing “macro” initialization strategies are difficult to detect, possibly resulting from differing model physics or responses to external forcing. Notably, the method of standardizing results affects ensemble spread: control simulations for most models have either large background trends or multi‐centennial variability in OHC. This spurious disequlibrium behavior is a substantial roadblock to understanding both internal climate variability and its response to forcing.

Stevenson, Samantha↗

DOE OSTI2022

Building a machine learning surrogate model for wildfire activities within a global Earth system model

Abstract. Wildfire is an important ecosystem process, influencing land biogeophysical and biogeochemical dynamics and atmospheric composition. Fire-driven loss of vegetation cover, for example, directly modifies the surface energy budget as a consequence of changing albedo, surface roughness, and partitioning of sensible and latent heat fluxes. Carbon dioxide and methane emitted by fires contribute to a positive atmospheric forcing, whereas emissions of carbonaceous aerosols may contribute to surface cooling. Process-based modeling of wildfires in Earth system land models is challenging due to limited understanding of human, climate, and ecosystem controls on fire counts, fire size, and burned area. Integration of mechanistic wildfire models within Earth system models requires careful parameter calibration, which is computationally expensive and subject to equifinality. To explore alternative approaches, we present a deep neural network (DNN) scheme that surrogates the process-based wildfire model with the Energy Exascale Earth System Model (E3SM) interface. The DNN wildfire model accurately simulates observed burned area with over 90 % higher accuracy with a large reduction in parameterization time compared with the current process-based wildfire model. The surrogate wildfire model successfully captured the observed monthly regional burned area during validation period 2011 to 2015 (coefficient of determination, R2=0.93). Since the DNN wildfire model has the same input and output requirements as the E3SM process-based wildfire model, our results demonstrate the applicability of machine learning for high accuracy and efficient large-scale land model development and predictions.

58 GEOSCIENCES↗

DOE OSTI2020

Initial Land Use/Cover Distribution Substantially Affects Global Carbon and Local Temperature Projections in the Integrated Earth System Model

Initial land cover distribution varies among Earth system models, an uncertainty in initial conditions that can substantially affect carbon and climate projections. We use the integrated Earth System Model to show that a 3.9 M km2 difference in 2005 global forest area (9–14% of total forest area) generates uncertainties in initial atmospheric CO2 concentration, terrestrial carbon, and local temperature that propagate through a future simulation following the Representative Concentration Pathway 4.5. By 2095, the initial 6 ppmv uncertainty range increases to 9 ppmv and the initial 26 PgC uncertainty range in terrestrial carbon increases to 33 PgC. The initial uncertainty range in annual average local temperature of -0.74 to 0.96 °C persists throughout the future simulation, with a seasonal maximum in Dec-Jan-Feb. These results highlight the importance of accurately characterizing historical land use and land cover to reduce overall initial condition uncertainty.

Di Vittorio, Alan↗

DOE OSTI2025

Implementation of a turbine farm model into the Energy Exascale Earth System Model for investigation and quantification of global climate impacts

Although there has been widespread deployment of wind farms in the United States, and plans to continue deployment into the future, the complete effects of wind farms on Earth systems are not well understood. The work performed here has incorporated wind farm models into the Energy Exascale Earth System Model (E3SM) capable of simulating the effects of extracting momentum from the atmospheric flow field using power generating wind farms. This new capability will allow scientists to quantify the impacts of wind farm induced changes on Earth systems by exploiting E3SM’s ability to couple atmospheric, oceanic, and biogeochemical (BGC) models on a global scale and monitor precipitation levels, extreme weather events, soil moisture content and jet stream location over decades-long time periods. This tool will be used to inform decision making on wind farm citing and will contribute to the Lab’s ability to assess energy technology impacts on the environment and evaluate the trade-offs between energy infrastructure investments and their impacts on natural systems.

17 WIND ENERGY↗

DOE OSTI2025

Improved Representations of Longwave Surface Emissivity to Reduce Surface and Atmospheric Heating Biases in Earth System Models

Abstract Many Earth system models (ESMs) approximate surface emissivity as a broadband constant. This approximation reduces the computational burden, yet omits the spectral structure of emissivity and atmospheric absorption. Neglecting spectral variation in surface emission introduces biases in longwave (LW) atmospheric fluxes and heating. Biases are strongest over surfaces with strongly varying emissivity and minimal atmospheric opacity. We examine these biases over water, ice, and snow surfaces. We partition spectral emissivity into the 16 spectral bands utilized by a single‐column atmospheric radiative transfer model (RRTMG_LW) commonly used in ESMs. We quantify flux and heating biases introduced by broadband assumptions relative to the spectrally resolved case for standard atmospheric profiles over each surface type. Current assumptions tend to overestimate upwelling surface fluxes; for example, the greybody assumption overestimates flux by 1.6 W/m 2 (0.52%) at the bottom of a mid‐latitude winter atmosphere over ice, and by 2.33 (1.0%) at the top of atmosphere. The blackbody assumption tends to artificially cool Earth's surface, stabilizing the lower troposphere. Interestingly, the optimal broadband emissivity can deviate from the Planck‐weighted mean by up to 3% depending on surface type and atmospheric profile. We investigate bias sensitivity to surface temperature, cloud water path, and atmospheric water vapor. Bias is most sensitive to water vapor content, and least sensitive to cloud water path. Lastly, we show that a modified greybody method with updated broadband values can reduce total surface flux bias up to 1.69 , comparable to a five‐band approach and at a fraction of the computational cost.

Manzo, L. [Department of Earth System Science Univ↗

DOE OSTI2020

Exploring topography-based methods for downscaling subgrid precipitation for use in Earth System Models

Topography exerts major control on land surface processes. To improve representation of topographic impacts on land surface processes, a new topography-based subgrid structure has been introduced to the Energy Exascale Earth System Model to represent the subgrid heterogeneity of surface elevation. To take advantage of the new subgrid structure for improving land surface modeling, this study explores four variations of topography-based methods for downscaling grid precipitation to the corresponding subgrids. In the first three methods, the deviation of the subgrid precipitation from the grid’s is equal to the grid precipitation multiplied by the ratio of the elevation difference between the subgrid and grid mean to a specified elevation equals to the grid elevation, the difference between the maximum and minimum subgrid elevation, and the maximum subgrid elevation, respectively. The second method limits the ratio to 0.5 to avoid extreme values on mountains and the third method accounts for the slope effect. The fourth method is similar to the third method except that the Froude Number is used to limit the numerator in a blocking regime of the ambient flow. The downscaled precipitation is evaluated using the PRISM precipitation data over the U.S. using statistical metrics. Results show that by accounting for topographic slope besides elevation, the third and fourth methods show clear advantages over the first and second methods. Furthermore, introducing the Froude Number in the fourth method improves downscaling skill and shows consistent advantages over the third method in areas with larger subgrid heterogeneity across different grids sizes.

Tesfa, Teklu K.↗

DOE OSTI2025

A Fortran-Python Interface for Integrating Machine Learning Parameterization into Earth System Models

Parameterizations in Earth System Models (ESMs) are subject to biases and uncertainties arising from subjective empirical assumptions and incomplete understanding of the underlying physical processes. Recently, the growing representational capability of machine learning (ML) in solving complex problems has spawned immense interests in climate science applications. Specifically, ML-based parameterizations have been developed to represent convection, radiation and microphysics processes in ESMs by learning from observations or high-resolution simulations, which have the potential to improve the accuracies and alleviate the uncertainties. Previous works have developed some surrogate models for these processes using ML. These surrogate models need to be coupled with the dynamical core of ESMs to investigate the effectiveness and their performance in a coupled system. In this study, we present a novel Fortran-Python interface designed to seamlessly integrate ML parameterizations into ESMs. This interface showcases high versatility by supporting popular ML frameworks like PyTorch, TensorFlow, and Scikit-learn. We demonstrate the interface's modularity and reusability through two cases: a ML trigger function for convection parameterization and a ML wildfire model. We conduct a comprehensive evaluation of memory usage and computational overhead resulting from the integration of Python codes into the Fortran ESMs. By leveraging this flexible interface, ML parameterizations can be effectively developed, tested, and integrated into ESMs.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Report on the LLNL Global Seismic Waveform Tomography Modeling Project

Earth models have important applications for seismic event monitoring including the location and characterization of potential underground nuclear tests. Seismic waves generated by earthquakes and man-made events are altered by variations in Earth’s subsurface properties; and these distortions effect the arrival times of energy packets (phases) and obscure the characteristics of the original energy source (event mechanism). Three-dimensional models of Earth’s seismic properties can help explain and predict the distortions to the seismic wavefield and reveal the properties of a source that generated the waves, including the event location using the model-predicted timing of the phases and source mechanism using the modelpredicted waveform characteristics.

58 GEOSCIENCES↗

DOE OSTI2021

Description and Demonstration of the Coupled Community Earth System Model v2 – Community Ice Sheet Model v2 (CESM2‐CISM2)

Abstract Earth system/ice‐sheet coupling is an area of recent, major Earth System Model (ESM) development. This work occurs at the intersection of glaciology and climate science and is motivated by a need for robust projections of sea‐level rise. The Community Ice Sheet Model version 2 (CISM2) is the newest component model of the Community Earth System Model version 2 (CESM2). This study describes the coupling and novel capabilities of the model, including: (1) an advanced energy‐balance‐based surface mass balance calculation in the land component with downscaling via elevation classes; (2) a closed freshwater budget from ice sheet to the ocean from surface runoff, basal melting, and ice discharge; (3) dynamic land surface types; and (4) dynamic atmospheric topography. The Earth system/ice‐sheet coupling is demonstrated in a simulation with an evolving Greenland Ice Sheet (GrIS) under an idealized high CO 2 scenario. The model simulates a large expansion of ablation areas (where surface ablation exceeds snow accumulation) and a large increase in surface runoff. This results in an elevated freshwater flux to the ocean, as well as thinning of the ice sheet and area retreat. These GrIS changes result in reduced Greenland surface albedo, changes in the sign and magnitude of sensible and latent heat fluxes, and modified surface roughness and overall ice sheet topography. Representation of these couplings between climate and ice sheets is key for the simulation of ice and climate interactions.

Muntjewerf, Laura↗

DOE OSTI2026

Knowledge Graph for End-to-End Traceability of an Integrated Human-Earth System Model

Integrated human-Earth system models inform energy-water-land system dynamics and policies, yet their results are difficult to trace through input-data, model structure, scenario configurations, and solved outputs. Because this information is siloed across disconnected artifacts, process-based IAMs have historically lacked a unified, queryable representation. Such lack of traceability prevents researchers from systematically isolating the multi-sector drivers of complex outcomes (such as tracing water-scarcity results back to distant energy-system dynamics) or conducting holistic uncertainty attribution across hundreds of interacting parameters. To address this concern, our work documents the software engineering process of a knowledge graph that unifies these four layers for the Global Change Analysis Model (GCAM-USA_Reference scenario, GCAM v9.1). The graph was built as a relational property graph in DuckDB from the run’s own artifacts: the input-preparation dependency map (gcamdata chunk map), the model’s XML input files, the run configuration, and the results database (BaseX), successfully mapping the model’s declared structure. The resulting graph comprises 204,321 nodes and 1,687,814 edges across 16 node types and 15 edge types, with approximately 16.3 million time-series values stored separately to maintain structural efficiency. To ensure representation fidelity, every edge carries an epistemic-status annotation recording the warrant for the relationship (structural, provenance, dependency, or model-derived), and a machine-readable provenance ledger classifying the origin of every schema element. Evaluation against a fixed five-benchmark suite with locked baselines reports zero structural orphans, zero dangling edge endpoints, and 100% of output-producing technologies traceable to raw input files. Two interactive interfaces present the graph, including a serverless browser application built on DuckDB-Wasm. By establishing the first end-to-end provenance framework for an IAM, this work enables researchers and scientists to systematically audit complex policy scenarios, debug model structures, and trace policy-relevant outputs to their data origins in real time.

Artifical Intelligence↗

DOE OSTI2021

Ocean Surface Flux Algorithm Effects on Earth System Model Energy and Water Cycles

Earth system models parameterize ocean surface fluxes of heat, moisture, and momentum with empirical bulk flux algorithms, which introduce biases and uncertainties into simulations. We investigate the atmosphere and ocean model sensitivity to algorithm choice in the Energy Exascale Earth System Model (E3SM). Flux differences between algorithms are larger in atmosphere simulations (where wind speeds can vary) than ocean simulations (where wind speeds are fixed by forcing data). Surface flux changes lead to global scale changes in the energy and water cycles, notably including ocean heat uptake and global mean precipitation rates. Compared to the control algorithm, both COARE and University of Arizona (UA) algorithms reduce global mean precipitation and top of atmosphere radiative biases. Further, UA may slightly reduce biases in ocean meridional heat transport. We speculate that changes seen here, especially in the ocean, could be even larger in coupled simulations.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2026

Modeling Large Dust Aerosols in the Community Earth System Model Version 2 (CESM2)

Dust aerosols have a wide size distribution from less than 0.1 to over 100 μm and dominate Earth's atmospheric aerosol mass. However, most Earth system models (ESMs) inadequately represent dust aerosols larger than 10 μm in diameter, limiting the accuracy of the simulated dust cycle and climate impacts. Here, we introduce a new modeling framework that captures the full observed size distribution of dust aerosols, incorporating recent advances into a mineral-resolved version of the Community ESM, while addressing known issues in previous versions. Comprehensive evaluation against diverse observations of bulk dust and component minerals demonstrates that the model reproduces the observed dust cycle across particle sizes. Incorporating the previously unrepresented large-dust fractions substantially alters dust budget estimates, highlighting potential changes in simulated climate impacts and underscoring the importance of comprehensive size-resolved dust modeling. Despite these advancements, uncertainties persist. Our results indicate that a size-dependent reduction in settling velocity is required to reproduce the observed dust size distribution downwind of source regions. Specifically, in the new model, the gravitational settling velocity of dust particles larger than 10 μm in diameter must be reduced by as much as 85% to achieve agreement with observations. This empirical reduction serves as a constraint on physics-based models of dust settling. Future developments should address misrepresented physical processes that hinder accurate modeling of the large dust aerosol transport. Expanding observational data sets covering the full-size distribution is also essential to better constrain the dust cycle and improve the representation of dust optical properties and climate effects.

Li, Longlei [Cornell Univ., Ithaca, NY (United Sta↗

DOE OSTI2023

Incorporating an Interactive Fire Plume-Rise Model in the DOE's Energy Exascale Earth System Model Version 1 (E3SMv1) and Examining Aerosol Radiative Effect

The vertical distribution of biomass burning aerosol (BBA) is important in regulating their impacts on weather and climate. The plume-rise process affects the injection height of BBA and interacts with the air parcel lifting and cloud processes. However, these processes are not represented in most global climate models. In this study, we replaced the fixed vertical profiles of monthly BBA emissions in the Department of Energy's Energy Exascale Earth System Model version 1 (E3SMv1) with an interactive fire plume-rise model. The vertical distribution of BBA emissions was calculated as a function of ambient thermodynamic conditions from the host E3SMv1, with distributions of fire sizes and sensible heat fluxes derived from the observations. The maximum fire radiative power (FRP) technique was used to determine the fire size. Scaling-FRP technique is used to calculate the wildfire heat release. Daily BBA emission, superimposed with a fire diurnal cycle retrieved from the satellite observation, was included in model simulations. The model shows improved agreement with satellite retrievals and in situ measurement during the National Oceanic and Atmospheric Administration Wildfire Experiment for Cloud chemistry, Aerosol absorption, and Nitrogen campaign. The model-observation comparison demonstrates the importance of the plume-rise model and fire diurnal cycle assumption in determining the BBA fields. We also find that E3SMv1 with new features produces a larger carbonaceous aerosol burden, leading to 0.13 W m –2 warming at the top of atmosphere compared to the default E3SMv1. This highlights the importance of accurately representing the BBA injection height and indicates a no-linear nature in the BBA-induced radiative effect.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Modeling Noise: Paths toward AI-Enabled Stochastic Earth System Models and Parameterizations

Some of the key challenges in Earth system prediction arise from an uncertain representation of unpredictable natural variability, across time and space scales. This variability, hereafter referred to as noise, represents the quantity against which the strength of a signal of interest is measured to assess predictability. Underestimating or overestimating this noise in Earth system models can lead to issues ranging between overconfidence in an erroneous prediction and a lack thereof in an otherwise accurate one. Noise is part of the physical processes in the Earth system and can amplify or damp a signal of interest in a complex manner that is not systematically characterized. Specifically, while there has been significant progress in our understanding of multiscale interactions among known modes of variability (potentially predictable signals), how these are impacted by the noise is not clear. Recent effort toward stochastic parametrization schemes that provide a representation of noise due to uncertain sub-grid processes in climate models (Berner et al., 2017 and references therein) has shown promise not only in reducing biases and improving probabilistic prediction but also in improving the representation of natural variability in the model toward what is observed. However, today’s standard schemes for stochastic parametrization are still simplistic as they do not take the state of the atmosphere into account (other than by the use of multiplicative noise for model tendencies). Furthermore they require hand-tuning and trial-and-error testing for the magnitude of the added noise, which is difficult in a chaotic system that is as complex as the Earth system.

58 GEOSCIENCES↗

DOE OSTI2022

Nitrification, denitrification, and competition for soil N : Evaluation of two Earth System Models against observations

Abstract Earth System Models (ESMs) have implemented nitrogen (N) cycles to account for N limitation on terrestrial carbon uptake. However, representing inputs, losses, and recycling of N in ESMs is challenging. Here, we use global rates and ratios of key soil N fluxes, including nitrification, denitrification, mineralization, leaching, immobilization, and plant uptake (both NH 4 + and NO 3 − ), from the literature to evaluate the N cycles in the land model components of two ESMs. The two land models evaluated here, E3SM Land Model version 1 (ELMv1)‐ECA and CLM5.0, originated from a common model but have diverged in their representation of plant–microbe competition for soil N. The models predict similar global rates of gross primary productivity (GPP) but have approximately two‐fold to three‐fold differences in their underlying global mineralization, immobilization, plant N uptake, nitrification, and denitrification fluxes. Both models dramatically underestimate the immobilization of NO 3 − by soil bacteria compared with literature values and predict dominance of plant uptake by a single form of mineral nitrogen (NO 3 − for ELM, with regional exceptions, and NH 4 + for CLM5.0). CLM5.0 strongly underestimates the global ratio of gross nitrification:gross mineralization and both models are likely to substantially underestimate the ratio of nitrification:denitrification. Few experimental data exist to evaluate this last ratio, in part because nitrification and denitrification are quantified using different techniques and because denitrification fluxes are difficult to measure at all. More observational constraints on soil nitrogen fluxes such as nitrification and denitrification, as well as greater scrutiny of the functional impact of introducing separate NH 4 + and NO 3 − pools into ESMs, could help to improve confidence in present and future simulations of N limitation on the carbon cycle.

Nevison, Cynthia↗