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.

Hydrodynamics-driven uncertainty study in simulating terrestrial carbon cycle with a coupled water-carbon cycle model

Uncertainty in Future High Flows in Qiantang River Basin, China

Uncertainty in High and Low Flows Due to Model Structure and Parameter Errors

Assessment of Extreme Flows and Uncertainty under Climate Change: Disentangling the Uncertainty Contribution of Representative Concentration Pathways, Global Climate Models and Internal Climate Variability

Identification of Key Physical Processes and Improvements for Simulating and Predicting Net Primary Production over the Tibetan Plateau

Sensitivity of Global Terrestrial Gross Primary Production to Hydrologic States Simulated by the Community Land Model Using Two Runoff Parameterizations

Overall uncertainty of climate change impacts on watershed hydrology in China

Uncertainty Analysis of Terrestrial Net Primary Productivity and Net Biome Productivity in China During 1901–2005

Uncertainties of gross primary productivity of Chinese grasslands based on multi-source estimation

The parametric uncertainty estimation of water and nitrogen transport simulation in a paddy field experiment using HYDRUS‐1D

Hydrological Projections under CMIP5 and CMIP6: Sources and Magnitudes of Uncertainty

Impact Assessment of Multiple Uncertainty Sources on High Flows under Climate Change

Evaluating the uncertainty of terrestrial water budget components over High Mountain Asia

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

Uncertainty propagation in a global biogeochemical model driven by leaf area data

Sources of Uncertainty in Modeled Land Carbon Storage Within and Across Three MIPs: Diagnosis with Three New Techniques

Decomposing Uncertainties in the Future Terrestrial Carbon Budget Associated with Emission Scenarios, Climate Projections, and Ecosystem Simulations Using the ISI-MIP Results

Vegetation Functional Properties Determine Uncertainty of Simulated Ecosystem Productivity: A Traceability Analysis in the East Asian Monsoon Region

Relative importance of uncertain model parameters driving water fluxes in a Land Surface Model

Quantifying uncertainty sources in extreme flow projections for three watersheds with different climate features in China

Estimation of Gross Primary Production in China (1982–2010) with Multiple Ecosystem Models