Terrestrial carbon cycle models have incorporated increasingly more processes as a means to achieve more-realistic representations of ecosystem carbon cycling. Despite this, there are large across-model variations in the simulation and projection of carbon cycling. Several model intercomparison projects (MIPs), for example, the fifth phase of the Coupled Model Intercomparison Project (CMIP5) (historical simulations), Trends in Net Land-Atmosphere Carbon Exchange (TRENDY), and Multiscale Synthesis and Terrestrial Model Intercomparison Project (MsTMIP), have sought to understand intermodel differences. In this study, the authors developed a suite of new techniques to conduct post-MIP analysis to gain insights into uncertainty sources across 25 models in the three MIPs. First, terrestrial carbon storage dynamics were characterized by a three-dimensional (3D) model output space with coordinates of carbon residence time, net primary productivity (NPP), and carbon storage potential. The latter represents the potential of an ecosystem to lose or gain carbon. This space can be used to measure how and why model output differs. Models with a nitrogen cycle generally exhibit lower annual NPP in comparison with other models, and mostly negative carbon storage potential. Second, a transient traceability framework was used to decompose any given carbon cycle model into traceable components and identify the sources of model differences. The carbon residence time (or NPP) was traced to baseline carbon residence time (or baseline NPP related to the maximum carbon input), environmental scalars, and climate forcing. Third, by applying a variance decomposition method, the authors show that the intermodel differences in carbon storage can be mainly attributed to the baseline carbon residence time and baseline NPP (>90% in the three MIPs). The three techniques developed in this study offer a novel approach to gain more insight from existing MIPs and can point out directions for future MIPs. Since this study is conducted at the global scale for an overview on intermodel differences, future studies should focus more on regional analysis to identify the sources of uncertainties and improve models at the specified mechanism level.

Soil is the largest terrestrial reservoir of organic carbon and is central for climate change mitigation and carbon-climate feedbacks. Chemical and physical associations of soil carbon with minerals play a critical role in carbon storage, but the amount and global capacity for storage in this form remain unquantified. Here, we produce spatially-resolved global estimates of mineral-associated organic carbon stocks and carbon-storage capacity by analyzing 1144 globally-distributed soil profiles. We show that current stocks total 899 Pg C to a depth of 1 m in non-permafrost mineral soils. Although this constitutes 66% and 70% of soil carbon in surface and deeper layers, respectively, it is only 42% and 21% of the mineralogical capacity. Regions under agricultural management and deeper soil layers show the largest undersaturation of mineral-associated carbon. Critically, the degree of undersaturation indicates sequestration efficiency over years to decades. We show that, across 103 carbon-accrual measurements spanning management interventions globally, soils furthest from their mineralogical capacity are more effective at accruing carbon; sequestration rates average 3-times higher in soils at one tenth of their capacity compared to soils at one half of their capacity. Our findings provide insights into the world's soils, their capacity to store carbon, and priority regions and actions for soil carbon management.

The enhanced vegetation growth by climate warming plays a pivotal role in amplifying the seasonal cycle of atmospheric CO 2 at northern lands (>50° N) since 1960s. However, the correlation between vegetation growth, temperature and seasonal amplitude of atmospheric CO 2 concentration have become elusive with the slowed increasing trend of vegetation growth and weakened temperature control on CO 2 uptake since late 1990s. Here, based on in situ atmospheric CO 2 concentration records from the Barrow observatory site, we found a slowdown in the increasing trend of the atmospheric CO 2 amplitude from 1990s to mid-2000s. This phenomenon was associated with the paused decrease in the minimum CO 2 concentration ([CO 2 ] min ), which was significantly correlated with the slowdown of vegetation greening and growing-season length extension. We then showed that both the vegetation greenness and growing-season length were positively correlated with spring but not autumn temperature over the northern lands. Furthermore, such asymmetric dependences of vegetation growth upon spring and autumn temperature cannot be captured by the state-of-art terrestrial biosphere models. These findings indicate that the responses of vegetation growth to spring and autumn warming are asymmetric, and highlight the need of improving autumn phenology in the models for predicting seasonal cycle of atmospheric CO 2 concentration.

Soil organic matter (SOM) supports the Earth's ability to sustain terrestrial ecosystems, provide food and fiber, and retains the largest pool of actively cycling carbon. Over 75% of the soil organic carbon (SOC) in the top meter of soil is directly affected by human land use. Large land areas have lost SOC as a result of land use practices, yet there are compensatory opportunities to enhance productivity and SOC storage in degraded lands through improved management practices. Large areas with and without intentional management are also being subjected to rapid changes in climate, making many SOC stocks vulnerable to losses by decomposition or disturbance. In order to quantify potential SOC losses or sequestration at field, regional, and global scales, measurements for detecting changes in SOC are needed. Such measurements and soil-management best practices should be based on well established and emerging scientific understanding of processes of C stabilization and destabilization over various timescales, soil types, and spatial scales. As newly engaged members of the International Soil Carbon Network, we have identified gaps in data, modeling, and communication that underscore the need for an open, shared network to frame and guide the study of SOM and SOC and their management for sustained production and climate regulation.

Soil carbon (C) is a critical component of Earth system models (ESMs), and its diverse representations are a major source of the large spread across models in the terrestrial C sink from the third to fifth assessment reports of the Intergovernmental Panel on Climate Change (IPCC). Improving soil C projections is of a high priority for Earth system modeling in the future IPCC and other assessments. To achieve this goal, we suggest that (1) model structures should reflect real-world processes, (2) parameters should be calibrated to match model outputs with observations, and (3) external forcing variables should accurately prescribe the environmental conditions that soils experience. First, most soil C cycle models simulate C input from litter production and C release through decomposition. The latter process has traditionally been represented by first-order decay functions, regulated primarily by temperature, moisture, litter quality, and soil texture. While this formulation well captures macroscopic soil organic C (SOC) dynamics, better understanding is needed of their underlying mechanisms as related to microbial processes, depth-dependent environmental controls, and other processes that strongly affect soil C dynamics. Second, incomplete use of observations in model parameterization is a major cause of bias in soil C projections from ESMs. Optimal parameter calibration with both pool- and flux-based data sets through data assimilation is among the highest priorities for near-term research to reduce biases among ESMs. Third, external variables are represented inconsistently among ESMs, leading to differences in modeled soil C dynamics. We recommend the implementation of traceability analyses to identify how external variables and model parameterizations influence SOC dynamics in different ESMs. Overall, projections of the terrestrial C sink can be substantially improved when reliable data sets are available to select the most representative model structure, constrain parameters, and prescribe forcing fields.