Convergence in simulating global soil organic carbon by structurally different models after data assimilation.

Current biogeochemical models produce carbon-climate feedback projections with large uncertainties, often attributed to their structural differences when simulating soil organic carbon (SOC) dynamics worldwide. However, choices of model parameter values that quantify the strength and represent properties of different soil carbon cycle processes could also contribute to model simulation uncertainties. Here, we demonstrate the critical role of using common observational data in reducing model uncertainty in estimates of global SOC storage. Two structurally different models featuring distinctive carbon pools, decomposition kinetics, and carbon transfer pathways simulate opposite global SOC distributions with their customary parameter values yet converge to similar results after being informed by the same global SOC database using a data assimilation approach. The converged spatial SOC simulations result from similar simulations in key model components such as carbon transfer efficiency, baseline decomposition rate, and environmental effects on carbon fluxes by these two models after data assimilation. Moreover, data assimilation results suggest equally effective simulations of SOC using models following either first-order or Michaelis-Menten kinetics at the global scale. Nevertheless, a wider range of data with high-quality control and assurance are needed to further constrain SOC dynamics simulations and reduce unconstrained parameters. New sets of data, such as microbial genomics-function relationships, may also suggest novel structures to account for in future model development. Overall, our results highlight the importance of observational data in informing model development and constraining model predictions.

Carbon sequestration in the subsoil and the time required to stabilize carbon for climate change mitigation.

Soils store large quantities of carbon in the subsoil (below 0.2 m depth) that is generally old and believed to be stabilized over centuries to millennia, which suggests that subsoil carbon sequestration (CS) can be used as a strategy for climate change mitigation. In this article, we review the main biophysical processes that contribute to carbon storage in subsoil and the main mathematical models used to represent these processes. Our guiding objective is to review whether a process understanding of soil carbon movement in the vertical profile can help us to assess carbon storage and persistence at timescales relevant for climate change mitigation. Bioturbation, liquid phase transport, belowground carbon inputs, mineral association, and microbial activity are the main processes contributing to the formation of soil carbon profiles, and these processes are represented in models using the diffusion-advection-reaction paradigm. Based on simulation examples and measurements from carbon and radiocarbon profiles across biomes, we found that advective and diffusive transport may only play a secondary role in the formation of soil carbon profiles. The difference between vertical root inputs and decomposition seems to play a primary role in determining the shape of carbon change with depth. Using the transit time of carbon to assess the timescales of carbon storage of new inputs, we show that only small quantities of new carbon inputs travel through the profile and can be stabilized for time horizons longer than 50 years, implying that activities that promote CS in the subsoil must take into consideration the very small quantities that can be stabilized in the long term.

Microbial carbon use efficiency promotes global soil carbon storage.

Soils store more carbon than other terrestrial ecosystems1,2. How soil organic carbon (SOC) forms and persists remains uncertain1,3, which makes it challenging to understand how it will respond to climatic change3,4. It has been suggested that soil microorganisms play an important role in SOC formation, preservation and loss5-7. Although microorganisms affect the accumulation and loss of soil organic matter through many pathways4,6,8-11, microbial carbon use efficiency (CUE) is an integrative metric that can capture the balance of these processes12,13. Although CUE has the potential to act as a predictor of variation in SOC storage, the role of CUE in SOC persistence remains unresolved7,14,15. Here we examine the relationship between CUE and the preservation of SOC, and interactions with climate, vegetation and edaphic properties, using a combination of global-scale datasets, a microbial-process explicit model, data assimilation, deep learning and meta-analysis. We find that CUE is at least four times as important as other evaluated factors, such as carbon input, decomposition or vertical transport, in determining SOC storage and its spatial variation across the globe. In addition, CUE shows a positive correlation with SOC content. Our findings point to microbial CUE as a major determinant of global SOC storage. Understanding the microbial processes underlying CUE and their environmental dependence may help the prediction of SOC feedback to a changing climate.

Soil respiration-driven CO2 pulses dominate Australia's flux variability.

The Australian continent contributes substantially to the year-to-year variability of the global terrestrial carbon dioxide (CO2) sink. However, the scarcity of in situ observations in remote areas prevents the deciphering of processes that force the CO2 flux variability. In this study, by examining atmospheric CO2 measurements from satellites in the period 2009-2018, we find recurrent end-of-dry-season CO2 pulses over the Australian continent. These pulses largely control the year-to-year variability of Australia's CO2 balance. They cause two to three times larger seasonal variations compared with previous top-down inversions and bottom-up estimates. The pulses occur shortly after the onset of rainfall and are driven by enhanced soil respiration preceding photosynthetic uptake in Australia's semiarid regions. The suggested continental-scale relevance of soil-rewetting processes has substantial implications for our understanding and modeling of global climate-carbon cycle feedbacks.

Increase in carbon input by enhanced fine root turnover in a long-term warmed forest soil.

Fine root litter represents an important carbon input to soils, but the effect of global warming on fine root turnover (FRT) is hardly explored in forest ecosystems. Understanding tree fine roots' response to warming is crucial for predicting soil carbon dynamics and the functioning of forests as a sink for atmospheric carbon dioxide (CO2). We studied fine root production (FRP) with ingrowth cores and used radiocarbon signatures of first-order, second- to third-order, and bulk fine roots to estimate fine root turnover times after 8 and 14 years of soil warming (+4 °C) in a temperate forest. Fine root turnover times of the individual root fractions were estimated with a one-pool model. Soil warming strongly increased fine root production by up to 128 % within one year, but after two years, the production was less pronounced (+35 %). The first-year production was likely very high due to the rapid exploitation of the root-free ingrowth cores. The radiocarbon signatures of fine roots were overall variable among treatments and plots. Soil warming tended to decrease fine root turnover times of all the measured root fractions after 8 and 14 years of warming, and there was a tendency for trees to use older carbon reserves for fine root production in warmed plots. Furthermore, soil warming increased fine root turnover from 50 to 106 g C m-2 yr-1 (based on two different approaches). Our findings suggest that future climate warming may increase carbon input into soils by enhancing fine root turnover. If this increase may partly offset carbon losses by increased mineralization of soil organic matter in temperate forest soils is still unclear and should guide future research.

Human and climate global-scale imprint on sediment transfer during the Holocene.

Accelerated soil erosion has become a pervasive feature on landscapes around the world and is recognized to have substantial implications for land productivity, downstream water quality, and biogeochemical cycles. However, the scarcity of global syntheses that consider long-term processes has limited our understanding of the timing, the amplitude, and the extent of soil erosion over millennial time scales. As such, we lack the ability to make predictions about the responses of soil erosion to long-term climate and land cover changes. Here, we reconstruct sedimentation rates for 632 lakes based on chronologies constrained by 3,980 calibrated 14C ages to assess the relative changes in lake-watershed erosion rates over the last 12,000 y. Estimated soil erosion dynamics were then complemented with land cover reconstructions inferred from 43,669 pollen samples and with climate time series from the Max Planck Institute Earth System Model. Our results show that a significant portion of the Earth surface shifted to human-driven soil erosion rate already 4,000 y ago. In particular, inferred soil erosion rates increased in 35% of the watersheds, and most of these sites showed a decrease in the proportion of arboreal pollen, which would be expected with land clearance. Further analysis revealed that land cover change was the main driver of inferred soil erosion in 70% of all studied watersheds. This study suggests that soil erosion has been altering terrestrial and aquatic ecosystems for millennia, leading to carbon (C) losses that could have ultimately induced feedbacks on the climate system.

Global covariation of carbon turnover times with climate in terrestrial ecosystems.

The response of the terrestrial carbon cycle to climate change is among the largest uncertainties affecting future climate change projections. The feedback between the terrestrial carbon cycle and climate is partly determined by changes in the turnover time of carbon in land ecosystems, which in turn is an ecosystem property that emerges from the interplay between climate, soil and vegetation type. Here we present a global, spatially explicit and observation-based assessment of whole-ecosystem carbon turnover times that combines new estimates of vegetation and soil organic carbon stocks and fluxes. We find that the overall mean global carbon turnover time is 23(+7)(-4) years (95 per cent confidence interval). On average, carbon resides in the vegetation and soil near the Equator for a shorter time than at latitudes north of 75° north (mean turnover times of 15 and 255 years, respectively). We identify a clear dependence of the turnover time on temperature, as expected from our present understanding of temperature controls on ecosystem dynamics. Surprisingly, our analysis also reveals a similarly strong association between turnover time and precipitation. Moreover, we find that the ecosystem carbon turnover times simulated by state-of-the-art coupled climate/carbon-cycle models vary widely and that numerical simulations, on average, tend to underestimate the global carbon turnover time by 36 per cent. The models show stronger spatial relationships with temperature than do observation-based estimates, but generally do not reproduce the strong relationships with precipitation and predict faster carbon turnover in many semi-arid regions. Our findings suggest that future climate/carbon-cycle feedbacks may depend more strongly on changes in the hydrological cycle than is expected at present and is considered in Earth system models.

Reconcilable differences: a joint calibration of fine-root turnover times with radiocarbon and minirhizotrons.

We used bomb-radiocarbon and raw minirhizotron lifetimes of fine roots (< 0.5 mm in diameter) in the organic layer of Norway spruce (Picea abies) forests in southern Sweden to test if different models are able to reconcile the apparently contradicting turnover time estimates from both techniques. We present a framework based on survival functions that is able to jointly model bomb-radiocarbon and minirhizotron data. At the same time we integrate prior knowledge about biases of both techniques--the classification of dead roots in minirhizotrons and the use of carbon reserves to grow new roots. Two-pool models, either in parallel or in serial setting, were able to reconcile the bomb-radiocarbon and minirhizotron data. These models yielded a mean residence time of 3.80 ± 0.16 yr (mean ± SD). On average 60 ± 2% of fine roots turned over within 0.75 ± 0.10 yr, while the rest was turning over within 8.4 ± 0.2 yr. Bomb-radiocarbon and minirhizotron data alone give a biased estimate of fine-root turnover. The two-pool models allow a mechanistic interpretation for the coexistence of fast- and slow-cycling roots--suberization and branching for the serial-two-pool model and branching due to ectomycorrhizal fungi-root interactions for the parallel-two-pool model.