skip to main content

Title: Divergent controls of soil organic carbon between observations and process-based models

The storage and cycling of soil organic carbon (SOC) are governed by multiple co-varying factors, including climate, plant productivity, edaphic properties, and disturbance history. Yet, it remains unclear which of these factors are the dominant predictors of observed SOC stocks, globally and within biomes, and how the role of these predictors varies between observations and process-based models. Here we use global observations and an ensemble of soil biogeochemical models to quantify the emergent importance of key state factors – namely, mean annual temperature, net primary productivity, and soil mineralogy – in explaining biome- to global-scale variation in SOC stocks. We use a machine-learning approach to disentangle the role of covariates and elucidate individual relationships with SOC, without imposing expected relationshipsa priori. While we observe qualitatively similar relationships between SOC and covariates in observations and models, the magnitude and degree of non-linearity vary substantially among the models and observations. Models appear to overemphasize the importance of temperature and primary productivity (especially in forests and herbaceous biomes, respectively), while observations suggest a greater relative importance of soil minerals. This mismatch is also evident globally. However, we observe agreement between observations and model outputs in select individual biomes – namely, temperate deciduous forests and grasslands, which both show stronger relationships of SOC stocks with temperature and productivity, respectively. This approach highlights biomes with the largest uncertainty and mismatch with observations for targeted model improvements. Understanding the role of dominant SOC controls, and the discrepancies between models and observations, globally and across biomes, is essential for improving and validating process representations in soil and ecosystem models for projections under novel future conditions.

more » « less
Award ID(s):
1929393 1637686
Author(s) / Creator(s):
; ; ; ; ; ; ; ; ; ; ; ;
Publisher / Repository:
Springer Science + Business Media
Date Published:
Journal Name:
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Estimates of soil organic carbon (SOC) stocks are essential for many environmental applications. However, significant inconsistencies exist in SOC stock estimates for the U.S. across current SOC maps. We propose a framework that combines unsupervised multivariate geographic clustering (MGC) and supervised Random Forests regression, improving SOC maps by capturing heterogeneous relationships with SOC drivers. We first used MGC to divide the U.S. into 20 SOC regions based on the similarity of covariates (soil biogeochemical, bioclimatic, biological, and physiographic variables). Subsequently, separate Random Forests models were trained for each SOC region, utilizing environmental covariates and SOC observations. Our estimated SOC stocks for the U.S. (52.6 ± 3.2 Pg for 0–30 cm and 108.3 ± 8.2 Pg for 0–100 cm depth) were within the range estimated by existing products like Harmonized World Soil Database, HWSD (46.7 Pg for 0–30 cm and 90.7 Pg for 0–100 cm depth) and SoilGrids 2.0 (45.7 Pg for 0–30 cm and 133.0 Pg for 0–100 cm depth). However, independent validation with soil profile data from the National Ecological Observatory Network showed that our approach (R2 = 0.51) outperformed the estimates obtained from Harmonized World Soil Database (R2 = 0.23) and SoilGrids 2.0 (R2 = 0.39) for the topsoil (0–30 cm). Uncertainty analysis (e.g., low representativeness and high coefficients of variation) identified regions requiring more measurements, such as Alaska and the deserts of the U.S. Southwest. Our approach effectively captures the heterogeneous relationships between widely available predictors and the current SOC baseline across regions, offering reliable SOC estimates at 1 km resolution for benchmarking Earth system models.

    more » « less
  2. null (Ed.)
    Wood formation consumes around 15% of the anthropogenic CO 2 emissions per year and plays a critical role in long-term sequestration of carbon on Earth. However, the exogenous factors driving wood formation onset and the underlying cellular mechanisms are still poorly understood and quantified, and this hampers an effective assessment of terrestrial forest productivity and carbon budget under global warming. Here, we used an extensive collection of unique datasets of weekly xylem tissue formation (wood formation) from 21 coniferous species across the Northern Hemisphere (latitudes 23 to 67°N) to present a quantitative demonstration that the onset of wood formation in Northern Hemisphere conifers is primarily driven by photoperiod and mean annual temperature (MAT), and only secondarily by spring forcing, winter chilling, and moisture availability. Photoperiod interacts with MAT and plays the dominant role in regulating the onset of secondary meristem growth, contrary to its as-yet-unquantified role in affecting the springtime phenology of primary meristems. The unique relationships between exogenous factors and wood formation could help to predict how forest ecosystems respond and adapt to climate warming and could provide a better understanding of the feedback occurring between vegetation and climate that is mediated by phenology. Our study quantifies the role of major environmental drivers for incorporation into state-of-the-art Earth system models (ESMs), thereby providing an improved assessment of long-term and high-resolution observations of biogeochemical cycles across terrestrial biomes. 
    more » « less
  3. Abstract

    The “hierarchy of factors” hypothesis states that decomposition rates are controlled primarily by climatic, followed by biological and soil variables. Tropical montane forests (TMF) are globally important ecosystems, yet there have been limited efforts to provide a biome‐scale characterization of litter decomposition. We designed a common litter decomposition experiment replicated in 23 tropical montane sites across the Americas, Asia, and Africa and combined these results with a previous study of 23 sites in tropical lowland forests (TLF). Specifically, we investigated (1) spatial heterogeneity in decomposition, (2) the relative importance of biological factors that affect leaf and wood decomposition in TMF, and (3) the role of climate in determining leaf litter decomposition rates within and across the TMF and TLF biomes. Litterbags of two mesh sizes containingLaurus nobilisleaves or birchwood popsicle sticks were spatially dispersed and incubated in TMF sites, for 3 and 7 months on the soil surface and at 10–15 cm depth. The within‐site replication demonstrated spatial variability in mass loss. Within TMF, litter type was the predominant biological factor influencing decomposition (leaves > wood), with mesh and burial effects playing a minor role. When comparing across TMF and TLF, climate was the predominant control over decomposition, but the Yasso07 global model (based on mean annual temperature and precipitation) only modestly predicted decomposition rate. Differences in controlling factors between biomes suggest that TMF, with their high rates of carbon storage, must be explicitly considered when developing theory and models to elucidate carbon cycling rates in the tropics.

    Abstract in Spanish is available with online material.

    more » « less
  4. Abstract

    Solar‐induced chlorophyll fluorescence (SIF) has been increasingly used as a proxy for terrestrial gross primary productivity (GPP). Previous work mainly evaluated the relationship between satellite‐observed SIF and gridded GPP products both based on coarse spatial resolutions. Finer resolution SIF (1.3 km × 2.25 km) measured from the Orbiting Carbon Observatory‐2 (OCO‐2) provides the first opportunity to examine the SIF–GPP relationship at the ecosystem scale using flux tower GPP data. However, it remains unclear how strong the relationship is for each biome and whether a robust, universal relationship exists across a variety of biomes. Here we conducted the first global analysis of the relationship between OCO‐2 SIF and tower GPP for a total of 64 flux sites across the globe encompassing eight major biomes. OCO‐2 SIF showed strong correlations with tower GPP at both midday and daily timescales, with the strongest relationship observed for daily SIF at the 757 nm (R2 = 0.72,p < 0.0001). Strong linear relationships between SIF and GPP were consistently found for all biomes (R2 = 0.57–0.79,p < 0.0001) except evergreen broadleaf forests (R2 = 0.16,p < 0.05) at the daily timescale. A higher slope was found for C4grasslands and croplands than for C3ecosystems. The generally consistent slope of the relationship among biomes suggests a nearly universal rather than biome‐specific SIF–GPP relationship, and this finding is an important distinction and simplification compared to previous results. SIF was mainly driven by absorbed photosynthetically active radiation and was also influenced by environmental stresses (temperature and water stresses) that determine photosynthetic light use efficiency. OCO‐2 SIF generally had a better performance for predicting GPP than satellite‐derived vegetation indices and a light use efficiency model. The universal SIF–GPP relationship can potentially lead to more accurate GPP estimates regionally or globally. Our findings revealed the remarkable ability of finer resolution SIF observations from OCO‐2 and other new or future missions (e.g., TROPOMI, FLEX) for estimating terrestrial photosynthesis across a wide variety of biomes and identified their potential and limitations for ecosystem functioning and carbon cycle studies.

    more » « less
  5. Abstract

    Wetland soils are a key global sink for organic carbon (C) and a focal point for C management and accounting efforts. The ongoing push for wetland restoration presents an opportunity for climate mitigation, but C storage expectations are poorly defined due to a lack of reference information and an incomplete understanding of what drives natural variability among wetlands. We sought to address these shortcomings by (1) quantifying the range of variability in wetland soil organic C (SOC) stocks on a depressional landscape (Delmarva Peninsula, USA) and (2) investigating the role of hydrology and relative topography in explaining variability among wetlands. We found a high degree of variability within individual wetlands and among wetlands with similar vegetation and hydrogeomorphic characteristics. This suggests that uncertainty should be presented explicitly when inferring ecosystem processes from wetland types or land cover classes. Differences in hydrologic regimes, particularly the rate of water level recession, explained some of the variability among wetlands, but relationships between SOC stocks and some hydrologic metrics were eclipsed by factors associated with separate study sites. Relative topography accounted for a similar portion of SOC stock variability as hydrology, indicating that it could be an effective substitute in large-scale analyses. As wetlands worldwide are restored and focus increases on quantifying C benefits, the importance of appropriately defining and assessing reference systems is paramount. Our results highlight the current uncertainty in this process, but suggest that incorporating landscape heterogeneity and drivers of natural variability into reference information may improve how wetland restoration is implemented and evaluated.

    more » « less