Soil Organic Matter Persistence as a Stochastic Process: Age and Transit Time Distributions of Carbon in Soils

The effect of repeated hurricanes on the age of organic carbon in humid tropical forest soil

Increasing hurricane frequency and intensity with climate change is likely to affect soil organic carbon (C) stocks in tropical forests. We examined the cycling of C between soil pools and with depth at the Luquillo Experimental Forest in Puerto Rico in soils over a 30-year period that spanned repeated hurricanes. We used a nonlinear matrix model of soil C pools and fluxes ("soilR") and constrained the parameters with soil and litter survey data. Soil chemistry and stable and radiocarbon isotopes were measured from three soil depths across a topographic gradient in 1988 and 2018. Our results suggest that pulses and subsequent reduction of inputs caused by severe hurricanes in 1989, 1998, and two in 2017 led to faster mean transit times of soil C in 0-10 cm and 35-60 cm depths relative to a modeled control soil with constant inputs over the 30-year period. Between 1988 and 2018, the occluded C stock increased and δ13C in all pools decreased, while changes in particulate and mineral-associated C were undetectable. The differences between 1988 and 2018 suggest that hurricane disturbance results in a dilution of the occluded light C pool with an influx of young, debris-deposited C, and possible microbial scavenging of old and young C in the particulate and mineral-associated pools. These effects led to a younger total soil C pool with faster mean transit times. Our results suggest that the increasing frequency of intense hurricanes will speed up rates of C cycling in tropical forests, making soil C more sensitive to future tropical forest stressors.

Controls on timescales of soil organic carbon persistence across sub-Saharan Africa

Given the importance of soil for the global carbon cycle, it is essential to understand not only how much carbon soil stores but also how long this carbon persists. Previous studies have shown that the amount and age of soil carbon are strongly affected by the interaction of climate, vegetation, and mineralogy. However, these findings are primarily based on studies from temperate regions and from fine-scale studies, leaving large knowledge gaps for soils from understudied regions such as sub-Saharan Africa. In addition, there is a lack of data to validate modeled soil C dynamics at broad scales. Here, we present insights into organic carbon cycling, based on a new broad-scale radiocarbon and mineral dataset for sub-Saharan Africa. We found that in moderately weathered soils in seasonal climate zones with poorly crystalline and reactive clay minerals, organic carbon persists longer on average (topsoil: 201 ± 130 years; subsoil: 645 ± 385 years) than in highly weathered soils in humid regions (topsoil: 140 ± 46 years; subsoil: 454 ± 247 years) with less reactive minerals. Soils in arid climate zones (topsoil: 396 ± 339 years; subsoil: 963 ± 669 years) store organic carbon for periods more similar to those in seasonal climate zones, likely reflecting climatic constraints on weathering, carbon inputs and microbial decomposition. These insights into the timescales of organic carbon persistence in soils of sub-Saharan Africa suggest that a process-oriented grouping of soils based on pedo-climatic conditions may be useful to improve predictions of soil responses to climate change at broader scales.

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.

Increased atmospheric CO2 and the transit time of carbon in terrestrial ecosystems

The response of terrestrial ecosystems to increased atmospheric CO2 concentrations is controversial and not yet fully understood, with previous large-scale forest manipulation experiments exhibiting contrasting responses. Although there is consensus that increased CO2 has a relevant effect on instantaneous processes such as photosynthesis and transpiration, there are large uncertainties regarding the fate of extra assimilated carbon in ecosystems. Filling this research gap is challenging because tracing the movement of new carbon across ecosystem compartments involves the study of multiple processes occurring over a wide range of timescales, from hours to millennia. We posit that a comprehensive quantification of the effect of increased CO2 must answer two interconnected questions: How much and for how long is newly assimilated carbon stored in ecosystems? Therefore, we propose that the transit time distribution of carbon is the key concept needed to effectively address these questions. Here, we show how the transit time distribution of carbon can be used to assess the fate of newly assimilated carbon and the timescales at which it is cycled in ecosystems. We use as an example a transit time distribution obtained from a tropical forest and show that most of the 60% of fixed carbon is respired in less than 1 year; therefore, we infer that under increased CO2 , most of the new carbon would follow a similar fate unless increased CO2 would cause changes in the rates at which carbon is cycled and transferred among ecosystem compartments. We call for a more frequent adoption of the transit time concept in studies seeking to quantify the ecosystem response to increased CO2 .

Relating mineral-organic matter stabilization mechanisms to carbon quality and age distributions using ramped thermal analysis

Organic carbon (OC) association with soil minerals stabilizes OC on timescales reflecting the strength of mineral-C interactions. We applied ramped thermal oxidation to subsoil B horizons with different mineral-C associations to separate OC according to increasing temperature of oxidation, i.e. thermal activation energy. Generally, OC released at lower temperatures was richer in bioavailable forms like polysaccharides, while OC released at higher temperatures was more aromatic. Organic carbon associated with pedogenic oxides was released at lower temperatures and had a narrow range of 14C content. By contrast, N-rich compounds were released at higher temperatures from samples with 2 : 1 clays and short-range ordered (SRO) amorphous minerals. Temperatures of release overlapped for SRO minerals and crystalline oxides, although the mean age of OC released was older for the SRO. In soils with more mixed mineralogy, the added presence of older OC released at temperatures greater than 450°C from clays resulted in a broader distribution of OC ages within the sample, especially for soils rich in 2 : 1 layer expandable clays such as smectite. While pedogenic setting affects mineral stability and absolute OC age, mineralogy controls the structure of OC age distribution within a sample, which may provide insight into model structures and OC dynamics under changing conditions. This article is part of the Theo Murphy meeting issue 'Radiocarbon in the Anthropocene'.

A decrease in the age of respired carbon from the terrestrial biosphere and increase in the asymmetry of its distribution

We provide here a model-based estimate of the transit time of carbon through the terrestrial biosphere, since the time of carbon uptake through photosynthesis until its release through respiration. We explored the consequences of increasing productivity versus increasing respiration rates on the transit time distribution and found that while higher respiration rates induced by higher temperature increase the transit time because older carbon is respired, increases in productivity cause a decline in transit times because more young carbon is available to supply increased metabolism. The combined effect of increases in temperature and productivity results in a decrease in transit times, with the productivity effect dominating over the respiration effect. By using an ensemble of simulation trajectories from the Carbon Data Model Framework (CARDAMOM), we obtained time-dependent transit time distributions incorporating the twentieth century global change. In these simulations, transit time declined over the twentieth century, suggesting an increased productivity effect that augmented the amount of respired young carbon, but also increasing the release of old carbon from high latitudes. The transit time distribution of carbon becomes more asymmetric over time, with more carbon transiting faster through tropical and temperate regions, and older carbon being respired from high latitude regions. This article is part of the Theo Murphy meeting issue 'Radiocarbon in the Anthropocene'.

Ecosystem and soil respiration radiocarbon detects old carbon release as a fingerprint of warming and permafrost destabilization with climate change

The permafrost region has accumulated organic carbon in cold and waterlogged soils over thousands of years and now contains three times as much carbon as the atmosphere. Global warming is degrading permafrost with the potential to accelerate climate change as increased microbial decomposition releases soil carbon as greenhouse gases. A 19-year time series of soil and ecosystem respiration radiocarbon from Alaska provides long-term insight into changing permafrost soil carbon dynamics in a warmer world. Nine per cent of ecosystem respiration and 23% of soil respiration observations had radiocarbon values more than 50‰ lower than the atmospheric value. Furthermore, the overall trend of ecosystem and soil respiration radiocarbon values through time decreased more than atmospheric radiocarbon values did, indicating that old carbon degradation was enhanced. Boosted regression tree analyses showed that temperature and moisture environmental variables had the largest relative influence on lower radiocarbon values. This suggested that old carbon degradation was controlled by warming/permafrost thaw and soil drying together, as waterlogged soil conditions could protect soil carbon from microbial decomposition even when thawed. Overall, changing conditions increasingly favoured the release of old carbon, which is a definitive fingerprint of an accelerating feedback to climate change as a consequence of warming and permafrost destabilization. This article is part of the Theo Murphy meeting issue 'Radiocarbon in the Anthropocene'.

Moisture-driven divergence in mineral-associated soil carbon persistence

Mineral stabilization of soil organic matter is an important regulator of the global carbon (C) cycle. However, the vulnerability of mineral-stabilized organic matter (OM) to climate change is currently unknown. We examined soil profiles from 34 sites across the conterminous USA to investigate how the abundance and persistence of mineral-associated organic C varied with climate at the continental scale. Using a novel combination of radiocarbon and molecular composition measurements, we show that the relationship between the abundance and persistence of mineral-associated organic matter (MAOM) appears to be driven by moisture availability. In wetter climates where precipitation exceeds evapotranspiration, excess moisture leads to deeper and more prolonged periods of wetness, creating conditions which favor greater root abundance and also allow for greater diffusion and interaction of inputs with MAOM. In these humid soils, mineral-associated soil organic C concentration and persistence are strongly linked, whereas this relationship is absent in drier climates. In arid soils, root abundance is lower, and interaction of inputs with mineral surfaces is limited by shallower and briefer periods of moisture, resulting in a disconnect between concentration and persistence. Data suggest a tipping point in the cycling of mineral-associated C at a climate threshold where precipitation equals evaporation. As climate patterns shift, our findings emphasize that divergence in the mechanisms of OM persistence associated with historical climate legacies need to be considered in process-based models.

Younger carbon dominates global soil carbon efflux

Soil carbon (C) is comprised of a continuum of organic compounds with distinct ages (i.e., the time a C atom has experienced in soil since the C atom entered soil). The contribution of different age groups to soil C efflux is critical for understanding soil C stability and persistence, but is poorly understood due to the complexity of soil C pool age structure and potential distinct turnover behaviors of age groups. Here, we build upon the quantification of soil C transit times to infer the age of C atoms in soil C efflux (a_efflux ) from seven sequential soil layer depths down to 2 m at a global scale, and compare this age with radiocarbon-inferred ages of C retained in corresponding soil layers (a_soil ). In the whole 0-2 m soil profile, the mean a_efflux is 194 21 1021 (mean with 5%-95% quantiles) year and is just about one-eighth of a_soil ( 1476 717 2547 year), demonstrating that younger C dominates soil C efflux. With increasing soil depth, both a_efflux and a_soil are increased, but their disparities are markedly narrowed. That is, the proportional contribution of relatively younger soil C to efflux is decreased in deeper layers, demonstrating that C inputs (new and young) stay longer in deeper layers. Across the globe, we find large spatial variability of the contribution of soil C age groups to C efflux. Especially, in deep soil layers of cold regions (e.g., boreal forests and tundra), a_efflux may be older than a_soil , suggesting that older C dominates C efflux only under a limited range of conditions. These results imply that most C inputs may not contribute to long-term soil C storage, particularly in upper layers that hold the majority of new C inputs.

Beyond bulk: Density fractions explain heterogeneity in global soil carbon abundance and persistence

Understanding the controls on the amount and persistence of soil organic carbon (C) is essential for predicting its sensitivity to global change. The response may depend on whether C is unprotected, isolated within aggregates, or protected from decomposition by mineral associations. Here, we present a global synthesis of the relative influence of environmental factors on soil organic C partitioning among pools, abundance in each pool (mg C g^-1  soil), and persistence (as approximated by radiocarbon abundance) in relatively unprotected particulate and protected mineral-bound pools. We show that C within particulate and mineral-associated pools consistently differed from one another in degree of persistence and relationship to environmental factors. Soil depth was the best predictor of C abundance and persistence, though it accounted for more variance in persistence. Persistence of all C pools decreased with increasing mean annual temperature (MAT) throughout the soil profile, whereas persistence increased with increasing wetness index (MAP/PET) in subsurface soils (30-176 cm). The relationship of C abundance (mg C g^-1  soil) to climate varied among pools and with depth. Mineral-associated C in surface soils (<30 cm) increased more strongly with increasing wetness index than the free particulate C, but both pools showed attenuated responses to the wetness index at depth. Overall, these relationships suggest a strong influence of climate on soil C properties, and a potential loss of soil C from protected pools in areas with decreasing wetness. Relative persistence and abundance of C pools varied significantly among land cover types and soil parent material lithologies. This variability in each pool's relationship to environmental factors suggests that not all soil organic C is equally vulnerable to global change. Therefore, projections of future soil organic C based on patterns and responses of bulk soil organic C may be misleading.

Land-Use Change Depletes Quantity and Quality of Soil Organic Matter Fractions in Ethiopian Highlands

The depletion of soil organic matter (SOM) reserve after deforestation and subsequent management practices are well documented, but the impacts of land-use change on the persistence and vulnerability of storage C and N remain uncertain. We investigated soil organic C (SOC) and N stocks in a landscape of chrono-sequence natural forest, grazing/crop lands and plantation forest in the highlands of North-West Ethiopia. We hypothesized that in addition to depleting total C and N pools, multiple conversions of natural forest significantly change the relative proportion of labile and recalcitrant C and N fractions in soils, and thus affect SOM quality. To examine this hypothesis, we estimated depletion of SOC and N stocks and labile (1 & 2) and recalcitrant (fraction 3) C and N pools in soil organic matter following the acid hydrolysis technique. Our studies showed the highest loss of C stock was in grazing land (58%) followed by cropland (50%) and eucalyptus plantation (47%), while on average ca. 57% N stock was depleted. Eucalyptus plantation exhibited potential for soil C recovery, although not for N, after 30 years. The fractionation of SOM revealed that depletions of labile 1 C stocks were similar in grazing and crop lands (36%), and loss of recalcitrant C was highest in grazing soil (56%). However, increases in relative concentrations of labile fraction 1 in grazing land and recalcitrant C and N in cropland suggest the quality of these pools might be influenced by management activities. Also, the C:N ratio of C fractions and recalcitrant indices (RIC and RIN) clearly demonstrated that land conversion from natural forest to managed systems changes the inherent quality of the fractions, which was obscured in whole soil analysis. These findings underscore the importance of considering the quality of SOM when evaluating disturbance impacts on SOC and N stocks.

Stochastic and deterministic interpretation of pool models

Carbon and element cycling models can be expressed in terms of the dynamics of individual particles or collection of them in aggregated pools. In both cases, the models represent the same dynamics and provide similar predictions. The time required for individual particles to pass through a system, that is, the transit time, can be obtained from both approaches. Pool models can be analyzed from a stochastic or a deterministic point of view.

Assessing the diverse environmental effects of biochar systems: An evaluation framework

Biochar has been recognised as a carbon dioxide removal (CDR) technology. Unlike other CDR technologies, biochar is expected to deliver various valuable effects in e.g. agriculture, animal husbandry, industrial processes, remediation activities and waste management. The diversity of biochar side effects to CDR makes the systematic environmental assessment of biochar projects challenging, and to date, there is no common framework for evaluating them. Our aim is to bridge the methodology gap for evaluating biochar systems from a life-cycle perspective. Using life cycle theory, actual biochar projects, and reviews of biochar research, we propose a general description of biochar systems, an overview of biochar effects, and an evaluation framework for biochar effects. The evaluation framework was applied to a case study, the Stockholm Biochar Project. In the framework, biochar effects are classified according to life cycle stage and life cycle effect type; and the biochar's end-of-life and the reference situations are made explicit. Three types of effects are easily included in life cycle theory: changes in biosphere exchanges, technosphere inputs, and technosphere outputs. For other effects, analysing the cause-effect chain may be helpful. Several biochar effects in agroecosystems can be modelled as future productivity increases against a reference situation. In practice, the complexity of agroecosystems can be bypassed by using empirical models. Existing biochar life cycle studies are often limited to carbon footprint calculations and quantify a limited amount of biochar effects, mainly carbon sequestration, energy displacements and fertiliser-related emissions. The methodological development in this study can be of benefit to the biochar and CDR research communities, as well as decision-makers in biochar practice and policy.

Long-term fate of rapidly eroding carbon stock soil profiles in coastal wetlands

Marsh edge erosion is one of the major causes of land and associated carbon loss in wetland-dominated coastlines. Assessing carbon stocks and understanding fate of eroding carbon is an essential component of wetland carbon budget. This study aims to understand the vertical soil carbon profile of an eroding marsh and potential mineralization of carbon in estuaries. Eleven soil cores (~2 m deep) were collected from the edge of four highly eroding marsh sites and three cores from the estuarine bottom (~50 cm deep). Cores were sectioned into 10-cm intervals and analyzed for total, labile and refractory carbon, carbon density, select enzyme and microbial activities, and organic and inorganic phosphorus forms. The total carbon, labile carbon, and carbon density increased with depth at all sites. The carbon density at 1-1.5 m deep (0.04 ± 0.003 g cm^-3) was significantly higher (p < 0.0001) than the top 1 m soil (0.032 ± 0.002 g cm^-3), indicating the need for considering deeper carbon profile for blue carbon stock assessment. The age of the carbon at the estuarine bottom was 388 ± 84 years before present (ybp) indicating the recently eroded wetland carbon is not reburied in the estuary. Significant anaerobic microbial activity was present at all the soil depths suggesting high potential of mineralization of eroded carbon in the aerobic estuarine water. The coastlines experiencing high relative sea-level rise at present or coastlines that are projecting high sea-level rise in the near future are susceptible to losing an enormous amount of previously sequestered carbon over a relatively short period of time.

From pools to flow: The PROMISE framework for new insights on soil carbon cycling in a changing world

Soils represent the largest terrestrial reservoir of organic carbon, and the balance between soil organic carbon (SOC) formation and loss will drive powerful carbon-climate feedbacks over the coming century. To date, efforts to predict SOC dynamics have rested on pool-based models, which assume classes of SOC with internally homogenous physicochemical properties. However, emerging evidence suggests that soil carbon turnover is not dominantly controlled by the chemistry of carbon inputs, but rather by restrictions on microbial access to organic matter in the spatially heterogeneous soil environment. The dynamic processes that control the physicochemical protection of carbon translate poorly to pool-based SOC models; as a result, we are challenged to mechanistically predict how environmental change will impact movement of carbon between soils and the atmosphere. Here, we propose a novel conceptual framework to explore controls on belowground carbon cycling: Probabilistic Representation of Organic Matter Interactions within the Soil Environment (PROMISE). In contrast to traditional model frameworks, PROMISE does not attempt to define carbon pools united by common thermodynamic or functional attributes. Rather, the PROMISE concept considers how SOC cycling rates are governed by the stochastic processes that influence the proximity between microbial decomposers and organic matter, with emphasis on their physical location in the soil matrix. We illustrate the applications of this framework with a new biogeochemical simulation model that traces the fate of individual carbon atoms as they interact with their environment, undergoing biochemical transformations and moving through the soil pore space. We also discuss how the PROMISE framework reshapes dialogue around issues related to SOC management in a changing world. We intend the PROMISE framework to spur the development of new hypotheses, analytical tools, and model structures across disciplines that will illuminate mechanistic controls on the flow of carbon between plant, soil, and atmospheric pools.

Dark microbial CO2 fixation in temperate forest soils increases with CO2 concentration

Dark, that is, nonphototrophic, microbial CO₂ fixation occurs in a large range of soils. However, it is still not known whether dark microbial CO₂ fixation substantially contributes to the C balance of soils and what factors control this process. Therefore, the objective of this study was to quantitate dark microbial CO₂ fixation in temperate forest soils, to determine the relationship between the soil CO₂ concentration and dark microbial CO₂ fixation, and to estimate the relative contribution of different microbial groups to dark CO₂ fixation. For this purpose, we conducted a ¹³ C-CO₂ labeling experiment. We found that the rates of dark microbial CO₂ fixation were positively correlated with the CO₂ concentration in all soils. Dark microbial CO₂ fixation amounted to up to 320 µg C kg^-1  soil day^-1 in the Ah horizon. The fixation rates were 2.8-8.9 times higher in the Ah horizon than in the Bw1 horizon. Although the rates of dark microbial fixation were small compared to the respiration rate (1.2%-3.9% of the respiration rate), our findings suggest that organic matter formed by microorganisms from CO₂ contributes to the soil organic matter pool, especially given that microbial detritus is more stable in soil than plant detritus. Phospholipid fatty acid analyses indicated that CO₂ was mostly fixed by gram-positive bacteria, and not by fungi. In conclusion, our study shows that the dark microbial CO₂ fixation rate in temperate forest soils increases in periods of high CO₂ concentrations, that dark microbial CO₂ fixation is mostly accomplished by gram-positive bacteria, and that dark microbial CO₂ fixation contributes to the formation of soil organic matter.