Fresh carbon input (above and belowground) contributes to soil carbon sequestration, but also accelerates decomposition of soil organic matter through biological priming mechanisms. Currently, poor understanding precludes the incorporation of these priming mechanisms into the global carbon models used for future projections. Here, we show that priming can be incorporated based on a simple equation calibrated from incubation and verified against independent litter manipulation experiments in the global land surface model, ORCHIDEE. When incorporated into ORCHIDEE, priming improved the model's representation of global soil carbon stocks and decreased soil carbon sequestration by 51% (12 ± 3 Pg C) during the period 1901-2010. Future projections with the same model across the range of CO₂ and climate changes defined by the IPCC-RCP scenarios reveal that priming buffers the projected changes in soil carbon stocks - both the increases due to enhanced productivity and new input to the soil, and the decreases due to warming-induced accelerated decomposition. Including priming in Earth system models leads to different projections of soil carbon changes, which are challenging to verify at large spatial scales.

The role of soil microorganisms in regulating soil organic matter (SOM) decomposition is of primary importance in the carbon cycle, in particular in the context of global change. Modeling soil microbial community dynamics to simulate its impact on soil gaseous carbon (C) emissions and nitrogen (N) mineralization at large spatial scales is a recent research field with the potential to improve predictions of SOM responses to global climate change. In this study we present a SOM model called ORCHIMIC, which utilizes input data that are consistent with those of global vegetation models. ORCHIMIC simulates the decomposition of SOM by explicitly accounting for enzyme production and distinguishing three different microbial functional groups: Fresh organic matter (FOM) specialists, SOM specialists, and generalists, while also implicitly accounting for microbes that do not produce extracellular enzymes, i.e., cheaters. ORCHIMIC and two other organic matter decomposition models, CENTURY (based on first-order kinetics and representative of the structure of most current global soil carbon models) and PRIM (with FOM accelerating the decomposition rate of SOM), were calibrated to reproduce the observed respiration fluxes of FOM and SOM from the incubation experiments of Blagodatskaya et al. (2014). Among the three models, ORCHIMIC was the only one that effectively captured both the temporal dynamics of the respiratory fluxes and the magnitude of the priming effect observed during the incubation experiment. ORCHIMIC also effectively reproduced the temporal dynamics of microbial biomass. We then applied different idealized changes to the model input data, i.e., a 5K stepwise increase of temperature and/or a doubling of plant litter inputs. Under 5K warming conditions, ORCHIMIC predicted a 0.002K-1 decrease in the C use efficiency (defined as the ratio of C allocated to microbial growth to the sum of C allocated to growth and respiration) and a 3% loss of SOC. Under the double litter input scenario, ORCHIMIC predicted a doubling of microbial biomass, while SOC stock increased by less than 1% due to the priming effect. This limited increase in SOC stock contrasted with the proportional increase in SOC stock as modeled by the conventional SOC decomposition model (CENTURY), which can not reproduce the priming effect. If temperature increased by 5K and litter input was doubled, ORCHIMIC predicted almost the same loss of SOC as when only temperature was increased. These tests suggest that the responses of SOC stock to warming and increasing input may differ considerably from those simulated by conventional SOC decomposition models when microbial dynamics are included. The next step is to incorporate the ORCHIMIC model into a global vegetation model to perform simulations for representative sites and future scenarios.

Soil nutrient availability has a strong influence on the fate of soil carbon (C) during microbial decomposition, contributing to Earth's C balance. While nutrient availability itself can impact microbial physiology and C partitioning between biomass and respiration during soil organic matter decomposition, the availability of labile C inputs may mediate the response of microorganisms to nutrient additions. As soil organic matter is decomposed, microorganisms retain or release C, nitrogen (N) or phosphorus (P) to maintain a stoichiometric balance. Although the concept of a microbial stoichiometric homeostasis has previously been proposed, microbial biomass CNP ratios are not static, and this may have very relevant implications for microbial physiological activities. Here, we tested the hypothesis that N, P and potassium (K) nutrient additions impact C cycling in a tropical soil due to microbial stoichiometric constraints to growth and respiration, and that the availability of energy-rich labile organic matter in the soil (i.e. leaf litter) mediates the response to nutrient addition. We incubated tropical soil from French Guiana with a 13C labeled leaf litter addition and with mineral nutrient additions of +K, +N, +NK, +PK and +NPK for 30 days. We found that litter additions led to a ten-fold increase in microbial respiration and a doubling of microbial biomass C, along with greater microbial N and P content. We found some evidence that P additions increased soil CO2 fluxes. Additionally, we found microbial biomass CP and NP ratios varied more widely than CN in response to nutrient and organic matter additions, with important implications for the role of microorganisms in C cycling. The addition of litter did not prime soil organic matter decomposition, except in combination with +NK fertilization, indicating possible P-mining of soil organic matter in this P-poor tropical soil. Together, these results point toward an ultimate labile organic substrate limitation of soil microorganisms in this tropical soil, but also indicate a complex interaction between C, N, P and K availability. This highlights the difference between microbial C cycling responses to N, P, or K additions in the tropics and explains why coupled C, N and P cycle modeling efforts cannot rely on strict microbial stoichiometric homeostasis as an underlying assumption.

Increasing atmospheric CO₂ stimulates photosynthesis which can increase net primary production (NPP), but at longer timescales may not necessarily increase plant biomass. Here we analyse the four decade-long CO₂-enrichment experiments in woody ecosystems that measured total NPP and biomass. CO₂ enrichment increased biomass increment by 1.05 ± 0.26 kg C m^-2 over a full decade, a 29.1 ± 11.7% stimulation of biomass gain in these early-secondary-succession temperate ecosystems. This response is predictable by combining the CO₂ response of NPP (0.16 ± 0.03 kg C m^-2 y^-1) and the CO₂-independent, linear slope between biomass increment and cumulative NPP (0.55 ± 0.17). An ensemble of terrestrial ecosystem models fail to predict both terms correctly. Allocation to wood was a driver of across-site, and across-model, response variability and together with CO₂-independence of biomass retention highlights the value of understanding drivers of wood allocation under ambient conditions to correctly interpret and predict CO₂ responses.

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.

Peatlands store substantial amounts of carbon and are vulnerable to climate change. We present a modified version of the Organising Carbon and Hydrology In Dynamic Ecosystems (ORCHIDEE) land surface model for simulating the hydrology, surface energy, and CO2 fluxes of peatlands on daily to annual timescales. The model includes a separate soil tile in each 0.5° grid cell, defined from a global peatland map and identified with peat-specific soil hydraulic properties. Runoff from non-peat vegetation within a grid cell containing a fraction of peat is routed to this peat soil tile, which maintains shallow water tables. The water table position separates oxic from anoxic decomposition. The model was evaluated against eddy-covariance (EC) observations from 30 northern peatland sites, with the maximum rate of carboxylation (Vcmax) being optimized at each site. Regarding short-term day-to-day variations, the model performance was good for gross primary production (GPP) (r2 Combining double low line 0.76; Nash-Sutcliffe modeling efficiency, MEF Combining double low line 0.76) and ecosystem respiration (ER, r2 Combining double low line 0.78, MEF Combining double low line 0.75), with lesser accuracy for latent heat fluxes (LE, r2 Combining double low line 0.42, MEF Combining double low line 0.14) and and net ecosystem CO2 exchange (NEE, r2 Combining double low line 0.38, MEF Combining double low line 0.26). Seasonal variations in GPP, ER, NEE, and energy fluxes on monthly scales showed moderate to high r2 values (0.57-0.86). For spatial across-site gradients of annual mean GPP, ER, NEE, and LE, r2 values of 0.93, 0.89, 0.27, and 0.71 were achieved, respectively. Water table (WT) variation was not well predicted (r2<0.1), likely due to the uncertain water input to the peat from surrounding areas. However, the poor performance of WT simulation did not greatly affect predictions of ER and NEE. We found a significant relationship between optimized Vcmax and latitude (temperature), which better reflects the spatial gradients of annual NEE than using an average Vcmax value.