Improving predictions of soil organic carbon (SOC) dynamics by multi-compartment models requires validation of turnover times of different SOC pools. Techniques such as laboratory incubation and isotope analysis have been adopted to estimate C turnover times, yet no studies have systematically compared these techniques and assessed the uncertainties associated with them. Here, we tested whether C turnover times of soil fractions were biased by methodology, and how this changed across soil particle sizes and ecosystems. We identified 52 studies that quantified C turnover times in different soil particles fractionated either according to aggregate size (e.g., macro- versus micro-aggregates) or according to soil texture (e.g., sand versus silt versus clay). C turnover times of these soil fractions were estimated by one of three methods: laboratory incubation (16 studies), δ13C shift due to C3–C4 vegetation change (25 studies), and 14C dating (19 studies). All methods showed that C turnover times of soil fractions generally increase with decreasing soil particle size. However, estimates of C turnover times within soil fractions differed significantly among methods, with incubation estimating the shortest turnover times and 14C the longest. The short C turnover times estimated by incubation are likely due to optimal environmental conditions for microbial decomposition existing in these studies, which is often a poor representation of field conditions. The 13C method can only be used when documenting a successive C3 versus C4 vegetation shift. C turnover times estimated by 14C were systematically higher than those estimated by 13C, especially for fine soil fractions (i.e., silt and clay). Overall, our findings highlight methodological uncertainties in estimating C turnover times of soil fractions, and correction factors should be explored to account for methodological bias when C turnover times estimated from different methods are used to parameterize soil C models.

Soil properties, such as clay content, are hypothesized to control decomposition of soil organic carbon (SOC). However, these hypotheses of soil property-C decomposition relationships have not been explicitly tested at large spatial scales. Here, we used a data-assimilation approach to evaluate the roles of soil properties and environmental factors in regulating decomposition of SOC. A three-pool (active, slow, and passive) C-cycling model was optimally fitted with 376 published laboratory incubation data from soils acquired from 73 sites with mean annual temperature ranging from -. 15 to 26. °C. Our results showed that soil physical and chemical properties regulated decomposition rates of the active and the slow C pools. Decomposition rates were lower for soils with high clay content, high field water holding capacity (WHC), and high C:N ratio. Multifactor regression and structural equation modeling (SEM) analyses showed that clay content was the most important variable in regulating decomposition of SOC. In contrast to the active and slow C pools, soil properties or environmental factors had little effect on the decomposition of the passive C pool. Our results show inverse soil property-C decomposition relationships and quantitatively evaluate the essential roles of soil texture (clay content) in controlling decomposition of SOC at a large spatial scale. The results may help model development and projection of changes in terrestrial C sequestration in the future.

Quantifying soil organic carbon (SOC) decomposition under warming is critical to predict carbon-climate feedbacks. According to the substrate regulating principle, SOC decomposition would decrease as labile SOC declines under field warming, but observations of SOC decomposition under warming do not always support this prediction. This discrepancy could result from varying changes in SOC components and soil microbial communities under warming. This study aimed to determine the decomposition of SOC components with different turnover times after subjected to long-term field warming and/or root exclusion to limit C input, and to test whether SOC decomposition is driven by substrate lability under warming. Taking advantage of a 12-year field warming experiment in a prairie, we assessed the decomposition of SOC components by incubating soils from control and warmed plots, with and without root exclusion for 3 years. We assayed SOC decomposition from these incubations by combining inverse modeling and microbial functional genes during decomposition with a metagenomic technique (GeoChip). The decomposition of SOC components with turnover times of years and decades, which contributed to 95% of total cumulative CO₂ respiration, was greater in soils from warmed plots. But the decomposition of labile SOC was similar in warmed plots compared to the control. The diversity of C-degradation microbial genes generally declined with time during the incubation in all treatments, suggesting shifts of microbial functional groups as substrate composition was changing. Compared to the control, soils from warmed plots showed significant increase in the signal intensities of microbial genes involved in degrading complex organic compounds, implying enhanced potential abilities of microbial catabolism. These are likely responsible for accelerated decomposition of SOC components with slow turnover rates. Overall, the shifted microbial community induced by long-term warming accelerates the decomposition of SOC components with slow turnover rates and thus amplify the positive feedback to climate change.

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.

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.

Atmospheric CO₂ concentration is increasing, largely due to anthropogenic activities. Previous studies of individual free-air CO₂ enrichment (FACE) experimental sites have shown significant impacts of elevated CO₂ (eCO₂) on soil microbial communities; however, no common microbial response patterns have yet emerged, challenging our ability to predict ecosystem functioning and sustainability in the future eCO₂ environment. Here we analyzed 66 soil microbial communities from five FACE sites, and showed common microbial response patterns to eCO₂, especially for key functional genes involved in carbon and nitrogen fixation (e.g., pcc/acc for carbon fixation, nifH for nitrogen fixation), carbon decomposition (e.g., amyA and pulA for labile carbon decomposition, mnp and lcc for recalcitrant carbon decomposition), and greenhouse gas emissions (e.g., mcrA for methane production, norB for nitrous oxide production) across five FACE sites. Also, the relative abundance of those key genes was generally increased and directionally associated with increased biomass, soil carbon decomposition, and soil moisture. In addition, a further literature survey of more disparate FACE experimental sites indicated increased biomass, soil carbon decay, nitrogen fixation, methane and nitrous oxide emissions, plant and soil carbon and nitrogen under eCO₂. A conceptual framework was developed to link commonly responsive functional genes with ecosystem processes, such as pcc/acc vs. soil carbon storage, amyA/pulA/mnp/lcc vs. soil carbon decomposition, and nifH vs. nitrogen availability, suggesting that such common responses of microbial functional genes may have the potential to predict ecosystem functioning and sustainability in the future eCO₂ environment.

Climate warming can result in both abiotic (e.g., permafrost thaw) and biotic (e.g., microbial functional genes) changes in Arctic tundra. Recent research has incorporated dynamic permafrost thaw in Earth system models (ESMs) and indicates that Arctic tundra could be a significant future carbon (C) source due to the enhanced decomposition of thawed deep soil C. However, warming-induced biotic changes may influence biologically related parameters and the consequent projections in ESMs. How model parameters associated with biotic responses will change under warming and to what extent these changes affect projected C budgets have not been carefully examined. In this study, we synthesized six data sets over 5 years from a soil warming experiment at the Eight Mile Lake, Alaska, into the Terrestrial ECOsystem (TECO) model with a probabilistic inversion approach. The TECO model used multiple soil layers to track dynamics of thawed soil under different treatments. Our results show that warming increased light use efficiency of vegetation photosynthesis but decreased baseline (i.e., environment-corrected) turnover rates of SOC in both the fast and slow pools in comparison with those under control. Moreover, the parameter changes generally amplified over time, suggesting processes of gradual physiological acclimation and functional gene shifts of both plants and microbes. The TECO model predicted that field warming from 2009 to 2013 resulted in cumulative C losses of 224 or 87 g/m² , respectively, without or with changes in those parameters. Thus, warming-induced parameter changes reduced predicted soil C loss by 61%. Our study suggests that it is critical to incorporate biotic changes in ESMs to improve the model performance in predicting C dynamics in permafrost regions.

Increases in carbon (C) inputs to soil can replenish soil organic C (SOC) through various mechanisms. However, recent studies have suggested that the increased C input can also stimulate the decomposition of old SOC via priming. Whether the loss of old SOC by priming can override C replenishment has not been rigorously examined. Here we show, through data-model synthesis, that the magnitude of replenishment is greater than that of priming, resulting in a net increase in SOC by a mean of 32% of the added new C. The magnitude of the net increase in SOC is positively correlated with the nitrogen-to-C ratio of the added substrates. Additionally, model evaluation indicates that a two-pool interactive model is a parsimonious model to represent the SOC decomposition with priming and replenishment. Our findings suggest that increasing C input to soils likely promote SOC accumulation despite the enhanced decomposition of old C via priming.

Realistic projection of future climate-carbon (C) cycle feedbacks requires better understanding and an improved representation of the C cycle in permafrost regions in the current generation of Earth system models. Here we evaluated 10 terrestrial ecosystem models for their estimates of net primary productivity (NPP) and responses to historical climate change in permafrost regions in the Northern Hemisphere. In comparison with the satellite estimate from the Moderate Resolution Imaging Spectroradiometer (MODIS; 246 ± 6 g C m−2 yr−1), most models produced higher NPP (309 ± 12 g C m−2 yr−1) over the permafrost region during 2000–2009. By comparing the simulated gross primary productivity (GPP) with a flux tower-based database, we found that although mean GPP among the models was only overestimated by 10% over 1982–2009, there was a twofold discrepancy among models (380 to 800 g C m−2 yr−1), which mainly resulted from differences in simulated maximum monthly GPP (GPPmax). Most models overestimated C use efficiency (CUE) as compared to observations at both regional and site levels. Further analysis shows that model variability of GPP and CUE are nonlinearly correlated to variability in specific leaf area and the maximum rate of carboxylation by the enzyme Rubisco at 25°C (Vcmax_25), respectively. The models also varied in their sensitivities of NPP, GPP, and CUE to historical changes in climate and atmospheric CO2 concentration. These results indicate that model predictive ability of the C cycle in permafrost regions can be improved by better representation of the processes controlling CUE and GPPmax as well as their sensitivity to climate change.

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.