Explore the vast range of titles available.

Surface soil moisture (mrsos) and vertically integrated soil moisture (mrsol) over the top 10 cm should, by definition, be physically consistent in Earth System Models (ESMs). However, an evaluation of nine CMIP6 models reveals substantial inconsistencies: in some models, mrsos and integrated mrsol agree globally; in others, they align only in specific regions; and in a few, they diverge across all grid cells. These discrepancies arise from a combination of factors, including metadata errors, inconsistent variable definitions, or diagnostic sequencing within the model. We demonstrate how such issues can lead to significant biases, even when both variables are present and seemingly well-defined. As model complexity increases and multi-model comparisons become more common, assumptions about variable equivalence may lead to flawed conclusions. This study highlights the need for routine consistency checks, improved metadata standards, and community-wide practices that ensure reliability of derived variables across ESM outputs, particularly in preparation for CMIP7.

The increasing complexity of Earth system models has inspired efforts to quantitatively assess model fidelity through rigorous comparison with best available measurements and observational data products. Earth system models exhibit a high degree of spread in predictions of land biogeochemistry, biogeophysics, and hydrology, which are sensitive to forcing from other model components. Based on insights from prior land model evaluation studies and community workshops, the authors developed an open source model benchmarking software package that generates graphical diagnostics and scores model performance in support of the International Land Model Benchmarking (ILAMB) project. Employing a suite of in situ, remote sensing, and reanalysis data sets, the ILAMB package performs comprehensive model assessment across a wide range of land variables and generates a hierarchical set of web pages containing statistical analyses and figures designed to provide the user insights into strengths and weaknesses of multiple models or model versions. Described here is the benchmarking philosophy and mathematical methodology embodied in the most recent implementation of the ILAMB package. Comparison methods unique to a few specific data sets are presented, and guidelines for configuring an ILAMB analysis and interpreting resulting model performance scores are discussed. ILAMB is being adopted by modeling teams and centers during model development and for model intercomparison projects, and community engagement is sought for extending evaluation metrics and adding new observational data sets to the benchmarking framework. Key Point The ILAMB benchmarking system broadly compares models to observational data sets and provides a synthesis of overall performance

Over the past several decades, the land modeling community has recognized the importance of nutrient regulation on the global terrestrial carbon cycle. Implementations of nutrient limitation in land models are diverse, varying from applying simple empirical down‐regulation of potential gross primary productivity under nutrient deficit conditions to more mechanistic treatments. In this study, we introduce a new approach to model multinutrient (nitrogen [N] and phosphorus [P]) limitations in the Energy Exascale Earth System Model (E3SM) Land Model version 1 (ELMv1‐ECA). The development is grounded on (1) advances in representing multiple‐consumer, multiple‐nutrient competition; (2) a generic dynamic allocation scheme based on water, N, P, and light availability; (3) flexible plant CNP stoichiometry; (4) prognostic treatment of N and P constraints on several carbon cycle processes; and (5) global data sets of plant physiological traits. Through benchmarking the model against best knowledge of global plant and soil carbon pools and fluxes, we show that our implementation of nutrient constraints on the present‐day carbon cycle is robust at the global scale. Compared with predecessor versions, ELMv1‐ECA better predicts global‐scale gross primary productivity, ecosystem respiration, leaf area index, vegetation biomass, soil carbon stocks, evapotranspiration, N2O emissions, and NO3‐ leaching. Factorial experiments indicate that representing the phosphorus cycle improves modeled carbon fluxes, while considering dynamic allocation improves modeled carbon stock density. We also highlight the value of using the International Land Model Benchmarking (ILAMB) package to evaluate and document performance during model development. Key Points We developed a new carbon‐nitrogen‐phosphorus land model (ELMv1‐ECA) integrated in the E3SM Earth system model We benchmarked the simulated present‐day carbon cycle using the International Land Model Benchmarking package (ILAMB) We documented model performance and identified necessary future improvements

Global estimates of the land carbon sink are often based on simulations by terrestrial biosphere models (TBMs). The use of a large number of models that differ in their underlying hypotheses, structure and parameters is one way to assess the uncertainty in the historical land carbon sink. Here we show that the atmospheric forcing datasets used to drive these TBMs represent a significant source of uncertainty that is currently not systematically accounted for in land carbon cycle evaluations. We present results from three TBMs each forced with three different historical atmospheric forcing reconstructions over the period 1850–2015. We perform an analysis of variance to quantify the relative uncertainty in carbon fluxes arising from the models themselves, atmospheric forcing, and model-forcing interactions. We find that atmospheric forcing in this set of simulations plays a dominant role on uncertainties in global gross primary productivity (GPP) (75% of variability) and autotrophic respiration (90%), and a significant but reduced role on net primary productivity and heterotrophic respiration (30%). Atmospheric forcing is the dominant driver (52%) of variability for the net ecosystem exchange flux, defined as the difference between GPP and respiration (both autotrophic and heterotrophic respiration). In contrast, for wildfire-driven carbon emissions model uncertainties dominate and, as a result, model uncertainties dominate for net ecosystem productivity. At regional scales, the contribution of atmospheric forcing to uncertainty shows a very heterogeneous pattern and is smaller on average than at the global scale. We find that this difference in the relative importance of forcing uncertainty between global and regional scales is related to large differences in regional model flux estimates, which partially offset each other when integrated globally, while the flux differences driven by forcing are mainly consistent across the world and therefore add up to a larger fractional contribution to global uncertainty.

Accurate assessment of anthropogenic carbon dioxide (CO2) emissions and their redistribution among the atmosphere, ocean, and terrestrial biosphere in a changing climate is critical to better understand the global carbon cycle, support the development of climate policies, and project future climate change. Here we describe and synthesise datasets and methodologies to quantify the five major components of the global carbon budget and their uncertainties. Fossil CO2 emissions (EFOS) are based on energy and cement production data. Emissions from land-use change (ELUC) are estimated by bookkeeping models based on land-use data. The global atmospheric CO2 growth rate (GATM) is computed from changes in concentration measured at surface stations. The global net uptake of CO2 by the ocean (SOCEAN) is estimated with global ocean biogeochemistry models and observation-based fCO2-products. The global net uptake of CO2 by the land (SLAND) is estimated with dynamic global vegetation models. Additional lines of evidence are provided by atmospheric inversions, atmospheric oxygen measurements, ocean interior observation-based estimates, and Earth System Models. This year, we introduced corrections on the ELUC, SOCEAN and SLAND estimates. The sum of all sources and sinks results in the carbon budget imbalance (BIM), a measure of imperfect data and incomplete understanding of the contemporary carbon cycle. All uncertainties are reported as ± 1σ. For the year 2024, EFOS increased by 1.1 % relative to 2023, with fossil emissions at 10.3 ± 0.5 GtC yr−1 (including the cement carbonation sink, 0.2 GtC yr−1), ELUC was 1.3 ± 0.7 GtC yr−1, for total anthropogenic CO2 emissions of 11.6 ± 0.9 GtC yr−1 (42.4 ± 3.2 GtCO2 yr−1). Also, for 2024, GATM was 7.9 ± 0.2 GtC yr−1 (3.73 ± 0.1 ppm yr−1), 2.2 GtC above the 2023 growth rate. SOCEAN was 3.4 ± 0.4 GtC yr−1 and SLAND was 1.9 ± 1.1 GtC yr−1, leaving a large negative BIM (−1.7 GtC yr−1), suggesting that the total sink or GATM is strongly overestimated in 2024. The global atmospheric CO2 concentration averaged over 2024 reached 422.8 ± 0.1 ppm. Preliminary data for 2025 suggest an increase in EFOS relative to 2024 of +1.0 % (0.2 % to 1.7 %) globally, and atmospheric CO2 concentration increasing by 2.1 ppm reaching 425.6 ppm, 53 % above the pre-industrial level (around 278 ppm in 1750). Overall, the mean and trend in the components of the global carbon budget are consistently estimated over the period 1959–2024, with a near-zero overall budget imbalance, although discrepancies of up to around 1 GtC yr−1 persist for the representation of annual to decadal variability in CO2 fluxes. Comparison of estimates from multiple approaches and observations shows: (1) a persistent large uncertainty in the estimate of land-use change emissions, (2) a low agreement between the different methods on the magnitude of the land CO2 flux in the northern extra-tropics, and (3) a discrepancy between the different methods on the mean ocean sink. This living data update documents changes in methods and datasets applied to this most-recent global carbon budget as well as evolving community understanding of the global carbon cycle. The data presented in this work are available at https://doi.org/10.18160/GCP-2025 (Friedlingstein et al., 2025c).

Land surface models require continuous validation against observations to improve and reduce simulation uncertainty. However, inferred model performance can be heavily influenced by subjective choices made in the selection and application of observational data products. A key area often misrepresented by models is the Arctic–Boreal region, which is a potential tipping point region in Earth’s climate system due to large permafrost carbon stocks that are vulnerable to release with climate warming. We use the International Land Model Benchmarking (ILAMB) framework to evaluate how the model skill of TRENDY-v9 models varies based on the choice of observational-based benchmark and how benchmarks are applied in model evaluation. This analysis uses global datasets integrated into ILAMB and new, regionally-specific observational products from the Arctic–Boreal Vulnerability Experiment. Our results cover the overall time period of 1979–2019 and show that model scores can vary substantially depending on the data product applied, with higher model scores indicating better model performance against observations. The lowest model scores occur when benchmarked against regional, compared to global, datasets. We also evaluate observed and modeled functional relationships between ecosystem respiration and air temperature and between gross primary production and precipitation. Here, we find that the magnitude and shape of the responses are strongly impacted by the choice of observational dataset and the approach used to construct the functional relationship benchmark. These results suggest that model evaluation studies could conclude a false sense of model skill if only using a single benchmark data product or if not applying regional data products when performing a regional model analysis. Collectively, our findings highlight the influence of benchmarking choices on model evaluation and point to the need for benchmarking guidelines when assessing model skill.

Earth system models (ESMs) have been rapidly developed in recent decades to advance our understanding of climate change‐carbon cycle feedback. However, those models are massive in coding, require expensive computational resources, and have difficulty in diagnosing their performance. It is highly desirable to develop ESMs with modularity and effective diagnostics. Toward these goals, we implemented a matrix approach to the Community Land Model version 5 (CLM5) to represent carbon and nitrogen cycles. Specifically, we reorganized 18 balance equations each for carbon and nitrogen cycles among the 18 vegetation pools in the original CLM5 into two matrix equations. Similarly, 140 balance equations each for carbon and nitrogen cycles among the 140 soil pools were reorganized into two additional matrix equations. The vegetation carbon and nitrogen matrix equations are connected to soil matrix equations via litterfall. The matrix equations fully reproduce simulations of carbon and nitrogen dynamics by the original model. The computational cost for forwarding simulation of the CLM5 matrix model was 26% more expensive than the original model, largely due to calculation of additional diagnostic variables, but the spin‐up computational cost was significantly saved. We showed a case study on modeled soil carbon storage under two forcing data sets to illustrate the diagnostic capability that the matrix approach uniquely offers to understand simulation results of global carbon and nitrogen dynamics. The successful implementation of the matrix approach to CLM5, one of the most complex land models, demonstrates that most, if not all, the biogeochemical models can be reorganized into the matrix form to gain high modularity, effective diagnostics, and accelerated spin‐up. Plain Language Summary Land models are widely used in climate change research. Due to the complex system, the model is not easily comprehended, nor the results are easily interpreted, even by a specialist. Enhancing the model tractability is imperative to make climate change prediction more effective, especially as models become more and more complex. In this study, we developed a matrix model by reorganizing six carbon and nitrogen modules of Community Land Model version 5 (CLM5) into four matrix equations to represent vegetation carbon, vegetation nitrogen, soil carbon, and soil nitrogen balance equations, respectively. The CLM5 matrix model gains high modularity, effective diagnostics, and accelerated spin‐up. The success of applying the matrix approach to CLM5, one of the most complex land models, support the theoretical analysis that the matrix approach is applicable to almost all land biogeochemical models. Key Points CLM5 matrix model reproduces carbon and nitrogen dynamics of the original model and offers more modularized code structure The diagnostic capability uniquely offered by the matrix approach make it easy to understand results of carbon and nitrogen dynamics The successful implementation of the matrix approach in CLM5 demonstrates its applicability to other biogeochemistry models

Emissions‐driven (prognostic CO2) simulations are essential for representing two‐way carbon‐climate feedback in Earth System Models. We present an emissions‐driven land–atmosphere coupled biogeochemistry (BGC) configuration (BGCLNDATM_progCO2) in version 2.1 of the Energy Exascale Earth System Model (E3SMv2.1). This is the first E3SM configuration that performs land‐atmosphere emission‐hindcasts. Here, we document its implementation, evaluate the model's performance against observations and other models, and propose a structured evaluation protocol for such emissions‐driven simulations. We conducted transient historical simulations (1850–2014) with BGCLNDATM_progCO2 and compare them to reference simulations—a land‐atmosphere coupled simulation without BGC and a standalone land simulation with BGC, both using prescribed CO2 concentrations—and to observations. BGCLNDATM_progCO2 overestimates atmospheric CO2 concentrations by 11–23 ppm yet stays within the 40‐ppm spread CMIP6 emission‐driven models and retains physical climate properties comparable to the reference runs. The CO2 biases are partly attributed to underrepresented oceanic CO2 uptake and inadequate representations of some terrestrial processes. In general, introducing prognostic CO2 did not change physical climate metrics at the global scale but had larger regional effects, particularly over land where spatially heterogeneous CO2 and prognostic leaf area index influenced surface energy balance. Finally, we propose a general evaluation protocol including spin‐up assessment, atmospheric CO2 benchmarking, physical climate evaluation, and land biogeochemical analysis to support scientific rigor and facilitate inter‐model comparisons. The new configuration lays the groundwork for future enhancements, including improved terrestrial biogeochemical processes, integrated marine biogeochemistry, and additional human–Earth system interactions. These developments advance E3SM toward fully coupled emissions‐driven simulations, enabling more accurate carbon–climate feedback projections and informing mitigation policy by providing physically consistent carbon‐budget metrics for mitigation scenarios. Plain Language Summary Understanding the impact of carbon dioxide (CO2) emissions on climate is vital for predicting future changes and crafting effective policies. Earth System Models (ESMs) are essential tools for simulating Earth's climate and assessing various influencing factors. In this study, we extended the Energy Exascale Earth System Model (E3SM)'s capabilities so that CO2 levels are calculated directly from human and natural emissions instead of being prescribed as a single global value. This extension allows for a more realistic representation of CO2 exchange between the atmosphere and land. We conducted historical simulations from 1850 to 2014 using this new development and compared results with observations and other models. Our model slightly overestimates atmospheric CO2 levels compared to measurements but is comparable to other models in capturing key climate features. To help other researchers build and test similar “emission‐driven” models, we created a step‐by‐step evaluation framework that checks CO2 behavior, climate variables, and land‐atmosphere interactions. Our work advances E3SM modeling by accurately representing how CO2 emissions affect Earth's systems. This enhancement lays the groundwork for modeling interactions between human‐Earth interactions, thereby enabling future studies that can inform mitigation and adaption. Key Points Implemented emissions‐driven land–atmosphere biogeochemistry in E3SMv2.1 (BGCLNDATM_progCO2), enabling prognostic CO2 simulations Established a structured evaluation protocol ensuring scientific rigor and facilitating inter‐model comparisons of model performance Emissions‐driven BGCLNDATM_progCO2 simulations maintain a physical climate similar to reference runs with prescribed CO2 concentrations

As regionally-relevant high-resolution Earth system data is increasingly relied upon across scientific, policy, and practitioner communities, there is an urgent need for coordinated and federated infrastructure to store, manage, standardize, and distribute decision-relevant community data products. Substantial effort is required to ensure that these products, which are often critical for regional impact assessments and decision-making, are findable, accessible, interoperable, and reusable. The Earth System Grid Federation US project (ESGF2-US) is addressing this challenge by expanding its open-source, distributed platform to support the hosting and dissemination of downscaled Earth system datasets. This expansion includes aligning new downscaled datasets with developing community standards for metadata and file structure, consistent with existing ESGF archives. This includes ensuring CF-compliance, applying CMORization where appropriate, and developing tools to streamline user access. In this paper, we highlight the technical and coordination work required to bring downscaled data into ESGF2-US and aim to inform the broader Earth system data user community about the growing availability and utility of these curated resources.

Vast carbon stocks stored in permafrost soils of Arctic tundra are under risk of release to the atmosphere under warming climate scenarios. Ice-wedge polygons in the low-gradient polygonal tundra create a complex mosaic of microtopographic features. This microtopography plays a critical role in regulating the fine-scale variability in thermal and hydrological regimes in the polygonal tundra landscape underlain by continuous permafrost. Modeling of thermal regimes of this sensitive ecosystem is essential for understanding the landscape behavior under the current as well as changing climate. We present here an end-to-end effort for high-resolution numerical modeling of thermal hydrology at real-world field sites, utilizing the best available data to characterize and parameterize the models. We develop approaches to model the thermal hydrology of polygonal tundra and apply them at four study sites near Barrow, Alaska, spanning across low to transitional to high-centered polygons, representing a broad polygonal tundra landscape. A multiphase subsurface thermal hydrology model (PFLOTRAN) was developed and applied to study the thermal regimes at four sites. Using a high-resolution lidar digital elevation model (DEM), microtopographic features of the landscape were characterized and represented in the high-resolution model mesh. The best available soil data from field observations and literature were utilized to represent the complex heterogeneous subsurface in the numerical model. Simulation results demonstrate the ability of the developed modeling approach to capture – without recourse to model calibration – several aspects of the complex thermal regimes across the sites, and provide insights into the critical role of polygonal tundra microtopography in regulating the thermal dynamics of the carbon-rich permafrost soils. Areas of significant disagreement between model results and observations highlight the importance of field-based observations of soil thermal and hydraulic properties for modeling-based studies of permafrost thermal dynamics, and provide motivation and guidance for future observations that will help address model and data gaps affecting our current understanding of the system.