Spatial biases reduce the ability of Earth system models to simulate soil heterotrophic respiration fluxes

Heterotrophic respiration (Rh) is, at a global scale, one of the largest CO2 fluxes between the Earth's surface and atmosphere and may increase in the future. The previous generation of Earth system models (ESMs) was able to reproduce global fluxes relatively well, but at that, time no gridded products were available to perform an in-depth evaluation. The capacity of the new generation of ESMs used within the Coupled Model Intercomparison Project Phase 6 (CMIP6) to reproduce this flux has not been evaluated, meaning that the realism of resulting CO2 flux estimates is unclear. In this study, we combine recently released observational data on Rh and ESM simulations to evaluate the ability of 13 ESMs from CMIP6 to reproduce Rh. Only 4 of the 13 tested ESMs were able to reproduce the total Rh flux, but spatial analysis underlined important bias compensation for most of the ESMs, which generally showed an overestimation in tropical regions and an underestimation in arid regions. To identify the main drivers of the bias, we performed an analysis of the residuals and found that mean annual precipitation was the most important driver explaining the difference between ESM simulations and observation-derived products of Rh, with a higher bias between ESM simulations and Rh products where precipitation was high. Based on our results, next-generation ESMs should focus on improving the response of Rh to soil moisture.</p

Soil organic carbon quantity, chemistry and thermal stability in a mountainous landscape : a Rock-Eval pyrolysis survey

International audienceMountain soils store huge amounts of carbon which may be highly vulnerable to the strong land use and climate changes that mountain areas currently experience worldwide. Here, we tested the Rock-Eval (RE) pyrolysis as a proxy technique to (i) quantify soil organic carbon (SOC) stocks, (ii) bring insights into SOC bulk chemistry and (iii) investigate biogeochemical stability at the landscape scale in a mountain area of the French calcareous Prealps. A total of 109 soils from 11 eco-units representing the variety of ecosystems of the study area were analyzed with RE pyrolysis. RE pyrolysis showed an excellent predictive performance (R2 = 0.99) for SOC content even in calcareous soils. The technique revealed specific chemical fingerprints for some eco-units and soil types, with decreasing hydrogen index values from Anthroposols (425 ± 62 mg HC/g SOC) to Umbrisols, Leptosols (311 ± 49 mg HC/g SOC) and to Cambisols (278 ± 35 mg HC/g SOC), associated with an increase in SOC maturation. Newly developed RE pyrolysis indices revealed the high stability of SOC in most eco-units developed on Cambisols (acidic grasslands, alpine meadows, bushy facies) and a significantly lower stability of SOC in mountain ridges, sheepfold areas and coniferous forest soils. The persistence of SOC in this mosaic of ecosystems may depend not only on its chemistry or thermal stability, but also on local environmental factors such as climatic conditions or pH, especially for high altitude soils. Overall, RE pyrolysis appears as an appropriate tool for landscape scale carbon inventories and could become a standardized proxy for assessing the vulnerability of SOC stocks

A model based on Rock-Eval thermal analysis to quantify the size of the centennially persistent organic carbon pool in temperate soils

Changes in global soil carbon stocks have considerable potential to influence the course of future climate change. However, a portion of soil organic carbon (SOC) has a very long residence time ( &gt;  100 years) and may not contribute significantly to terrestrial greenhouse gas emissions during the next century. The size of this persistent SOC reservoir is presumed to be large. Consequently, it is a key parameter required for the initialization of SOC dynamics in ecosystem and Earth system models, but there is considerable uncertainty in the methods used to quantify it. Thermal analysis methods provide cost-effective information on SOC thermal stability that has been shown to be qualitatively related to SOC biogeochemical stability. The objective of this work was to build the first quantitative model of the size of the centennially persistent SOC pool based on thermal analysis. We used a unique set of 118 archived soil samples from four agronomic experiments in northwestern Europe with long-term bare fallow and non-bare fallow treatments (e.g., manure amendment, cropland and grassland) as a sample set for which estimating the size of the centennially persistent SOC pool is relatively straightforward. At each experimental site, we estimated the average concentration of centennially persistent SOC and its uncertainty by applying a Bayesian curve-fitting method to the observed declining SOC concentration over the duration of the long-term bare fallow treatment. Overall, the estimated concentrations of centennially persistent SOC ranged from 5 to 11 g C kg−1 of soil (lowest and highest boundaries of four 95 % confidence intervals). Then, by dividing the site-specific concentrations of persistent SOC by the total SOC concentration, we could estimate the proportion of centennially persistent SOC in the 118 archived soil samples and the associated uncertainty. The proportion of centennially persistent SOC ranged from 0.14 (standard deviation of 0.01) to 1 (standard deviation of 0.15). Samples were subjected to thermal analysis by Rock-Eval 6 that generated a series of 30 parameters reflecting their SOC thermal stability and bulk chemistry. We trained a nonparametric machine-learning algorithm (random forests multivariate regression model) to predict the proportion of centennially persistent SOC in new soils using Rock-Eval 6 thermal parameters as predictors. We evaluated the model predictive performance with two different strategies. We first used a calibration set (n  =  88) and a validation set (n  =  30) with soils from all sites. Second, to test the sensitivity of the model to pedoclimate, we built a calibration set with soil samples from three out of the four sites (n  =  84). The multivariate regression model accurately predicted the proportion of centennially persistent SOC in the validation set composed of soils from all sites (R2  =  0.92, RMSEP  =  0.07, n  =  30). The uncertainty of the model predictions was quantified by a Monte Carlo approach that produced conservative 95 % prediction intervals across the validation set. The predictive performance of the model decreased when predicting the proportion of centennially persistent SOC in soils from one fully independent site with a different pedoclimate, yet the mean error of prediction only slightly increased (R2  =  0.53, RMSEP  =  0.10, n  =  34). This model based on Rock-Eval 6 thermal analysis can thus be used to predict the proportion of centennially persistent SOC with known uncertainty in new soil samples from different pedoclimates, at least for sites that have similar Rock-Eval 6 thermal characteristics to those included in the calibration set. Our study reinforces the evidence that there is a link between the thermal and biogeochemical stability of soil organic matter and demonstrates that Rock-Eval 6 thermal analysis can be used to quantify the size of the centennially persistent organic carbon pool in temperate soils

Dynamics of soil organic matter based on new Rock-Eval indices

This paper aims to develop a new approach to interpret soil organic matter (SOM) dynamics from Rock-Eval pyrolysis. Rock-Eval standard parameters (TpS2, HI, OI) are limited when applied to SOM, as they were defined for tracking the origin of sedimentary organic matter (i.e. terrestrial vs aquatic and lacustrine vs marine). This study proposes new Rock Eval based indices, projected on a new diagram plotting > 1300 samples, including litter and soil horizons, pure compounds and organic material. These new parameters reflect the thermal stability of SOM rather than its bulk chemistry. Their calculations are based on the contribution of four different areas (A1 to A4) integrated below the S2 pyrogram (amounts of released hydrocarbon compounds during the pyrolysis step). Results demonstrate that the relative values of A1 to A4 parameters can be used to survey the evolution of SOM during pedogenesis. These Rock-Eval parameters revealed a consistent thermal differentiation of SOM with depth, from thermally labile biogenic SOM in soil organic horizons to thermally stable mineral-associated SOM compounds in organo-mineral and mineral soil horizons, indicating a pedogenetic stabilization of SOM. Finally, newly defined I- and R-indices integrate the respective parts of A1 to A4 parameters into SOM dynamics, the I-index emphasizing the degree of transformation of the immature organic fraction (related to SOM stabilization), the R-index highlighting the contribution of the most refractory fraction or persistent SOM (related to pedogenic and inherited contributions). Analyses of a wide range of soils and reference samples (cellulose, lignin, humic substances, lignite, charcoal, coal, etc.) allow end-members as well as particular trends (humic, spodic, inherited) to be drawn. Specific examples are provided in order to illustrate some applications and uses for this new Rock-Eval based I/R diagram, from the study of grain size fractions to the evolution of SOM in soil chronosequences

Data from: Peatland vascular plant functional types affect methane dynamics by altering microbial community structure

1. Peatlands are natural sources of atmospheric methane (CH4), an important greenhouse gas. It is established that peatland methane dynamics are controlled by both biotic and abiotic conditions, yet the interactive effect of these drivers is less studied and consequently poorly understood. 2. Climate change affects the distribution of vascular plant functional types (PFTs) in peatlands. By removing specific PFTs, we assessed their effects on peat organic matter chemistry, microbial community composition and on potential methane production (PMP) and oxidation (PMO) in two microhabitats (lawns and hummocks). 3. Whilst PFT removal only marginally altered the peat organic matter chemistry, we observed considerable changes in microbial community structure. This resulted in altered PMP and PMO. PMP was slightly lower when graminoids were removed, whilst PMO was highest in the absence of both vascular PFTs (graminoids and ericoids), but only in the hummocks. 4. Path analyses demonstrate that different plant–soil interactions drive PMP and PMO in peatlands and that changes in biotic and abiotic factors can have auto-amplifying effects on current CH4 dynamics. 5. Synthesis. Changing environmental conditions will, both directly and indirectly, affect peatland processes, causing unforeseen changes in CH4 dynamics. The resilience of peatland CH4 dynamics to environmental change therefore depends on the interaction between plant community composition and microbial communities

Ecosystem respiration

To asses the effect of vascular plant removal in aforementioned experiment, one year after the start of the experiment, we measured carbon dioxide respiration rates on 10 cm diameter collars using an automated soil CO2 flux system (LI-8100, LI-COR Biosciences, USA). We compare the fluxes from the vascular plant removal plots with respiration rates from comparable areas in the control plots with little or no vascular plant cover, and assessed the relationship with the amount of removed biomass