Abstract:The spatial and temporal variations in terrestrial carbon storage play a pivotal role in regulating future climate change. However, Earth system models (ESMs), which have coupled the terrestrial biosphere and atmosphere, show great uncertainty in simulating the global land carbon storage. Here, based on multiple global datasets and a traceability analysis, we diagnosed the uncertainty source of terrestrial carbon storage in 22 ESMs that participated in phases 5 and 6 of the Coupled Model Intercomparison Project (CMIP5 and CMIP6). The modeled global terrestrial carbon storage has converged among ESMs from CMIP5 (1936.9 ± 739.3 PgC) to CMIP6 (1774.4 ± 439.0 PgC) but is persistently lower than the observation-based estimates (2285 ± 669 PgC). By further decomposing terrestrial carbon storage into net primary production (NPP) and ecosystem carbon residence time ( τ E ), we found that the decreased intermodel spread in land carbon storage primarily resulted from more accurate simulations on NPP among ESMs from CMIP5 to CMIP6. The persistent underestimation of land carbon storage was caused by the biased τ E . In CMIP5 and CMIP6, the modeled τ E was far shorter than the observation-based estimates. The potential reasons for the biased τ E could be the lack of or incomplete representation of nutrient limitation, vertical soil biogeochemistry, and the permafrost carbon cycle. Moreover, the modeled τ E became the key driver for the intermodel spread in global land carbon storage in CMIP6. Overall, our study indicates that CMIP6 models have greatly improved the terrestrial carbon cycle, with a decreased model spread in global terrestrial carbon storage and less uncertain productivity. However, more efforts are needed to understand and reduce the persistent data–model disagreement on carbon storage and residence time in the terrestrial biosphere. The spatial and temporal variations in terrestrial carbon storage play a pivotal role in regulating future climate change. However, Earth system models (ESMs), which have coupled the terrestrial biosphere and atmosphere, show great uncertainty in simulating the global land carbon storage. Here, based on multiple global datasets and a traceability analysis, we diagnosed the uncertainty source of terrestrial carbon storage in 22 ESMs that participated in phases 5 and 6 of the Coupled Model Intercomparison Project (CMIP5 and CMIP6). The modeled global terrestrial carbon storage has converged among ESMs from CMIP5 (1936.9 ± 739.3 PgC) to CMIP6 (1774.4 ± 439.0 PgC) but is persistently lower than the observation-based estimates (2285 ± 669 PgC). By further decomposing terrestrial carbon storage into net primary production (NPP) and ecosystem carbon residence time ( τ E ), we found that the decreased intermodel spread in land carbon storage primarily resulted from more accurate simulations on NPP among ESMs from CMIP5 to CMIP6. The persistent underestimation of land carbon storage was caused by the biased τ E . In CMIP5 and CMIP6, the modeled τ E was far shorter than the observation-based estimates. The potential reasons for the biased τ E could be the lack of or incomplete representation of nutrient limitation, vertical soil biogeochemistry, and the permafrost carbon cycle. Moreover, the modeled τ E became the key driver for the intermodel spread in global land carbon storage in CMIP6. Overall, our study indicates that CMIP6 models have greatly improved the terrestrial carbon cycle, with a decreased model spread in global terrestrial carbon storage and less uncertain productivity. However, more efforts are needed to understand and reduce the persistent data–model disagreement on carbon storage and residence time in the terrestrial biosphere.

Overestimated nitrogen loss from denitrification for natural terrestrial ecosystems in CMIP6 Earth System Models

Soil Nitrous Oxide Emissions Across the Northern High Latitudes

A global synthesis of soil denitrification: Driving factors and mitigation strategies

Feedback of carbon and nitrogen cycles enhances carbon sequestration in the terrestrial biosphere

Isotopic Elucidation of Microbial Nitrogen Transformations in Forest Soils

Warming and Redistribution of N Inputs Drive Long-Term Increase in Terrestrial N2O Emission Factor

Evaluating nitrogen cycling in terrestrial biosphere models: a disconnect between the carbon and nitrogen cycles

Microbial denitrification dominates nitrate losses from forest ecosystems

Global soil nitrous oxide emissions since the preindustrial era estimated by an ensemble of terrestrial biosphere models: Magnitude, attribution, and uncertainty

Benchmarking Carbon Fluxes of the Isimip2a Biome Models

Warming and Redistribution of Nitrogen Inputs Drive an Increase in Terrestrial Nitrous Oxide Emission Factor

Simulated Effects of Nitrogen Saturation on the Global Carbon Budget Using the IBIS Model

Marine Biogeochemical Cycling and Oceanic CO2 Uptake Simulated by the NUIST Earth System Model Version 3 (NESM V3)

Spatial Changes in Nitrogen Inputs Drive Short- and Long-Term Variability in Global Nitrous Oxide Emissions

Evolution of Uncertainty in Terrestrial Carbon Storage in Earth System Models from CMIP5 to CMIP6

Denitrification dominates dissimilatory nitrate reduction across global natural ecosystems

Representation of the Terrestrial Carbon Cycle in CMIP6

Characterizing the Importance of Denitrification for N2O Production in Soils Using Natural Abundance and Isotopic Labeling Techniques

Soil Denitrification Fluxes in a Northern Hardwood Forest: The Importance of Snowmelt and Implications for Ecosystem N Budgets

Assessment of Model Estimates of Land-Atmosphere CO2 Exchange Across Northern Eurasia

Projecting atmospheric N 2 O rise until the end of the 21st century: an Earth system model study