INTRODUCTION

Soil organic matter (SOM) represents the single largest actively cycling reservoir of terrestrial organic carbon, accounting for more than three times as much carbon as that present in the atmosphere or terrestrial vegetation (Schmidt et al. 2011; Lehmann and Kleber 2015). SOM is vulnerable to decomposition to either CO2 or CH4, which can increase atmospheric greenhouse gas concentrations (GHGs) and serve as a positive feedback to climate change. Conversely, the formation and stabilization of SOM within aggregates or associated with soil minerals can lead to carbon sequestration, representing a negative feedback to climate change. However, the conundrum as to why some SOM decomposes rapidly, while other thermodynamically unstable SOM can persist on centennial time scales (Hedges et al. 2000), leads to substantial uncertainty in model structures, as well as uncertainty in the predictability of the land carbon sink trajectory.

This chapter aims to tackle a part of this problem through generating recommendations for improved model structure in the representation of SOM cycling within global Earth System Models. The chapter is organized as follows: We first review the recent work contributing to the understanding of SOM stability. Then, we discuss current model structures and how these model structures deal with various processes, including data-model integration, microbial modeling, organo–mineral interactions, and SOM persistence. Finally, we conclude with recommendations for the future development of soil carbon models and the required relevant processes that determine the stability of SOM.

THE EMERGENT PICTURE OF SOIL ORGANIC MATTER STABILITY

Traditional conceptual models of SOM formation and stability have focused on the chemical structure of SOM compounds. These models consider the recalcitrance of certain SOM compounds to biological decomposition as the primary determinant of persistence and long-term stability (Schmidt et al. 2011; Lehmann and Kleber 2015). Within this perspective, plant-derived polymers are the primary source of organic matter (Fig. 1a). These polymers were thought to form through a suite of reactions, via biotic (Fig. 1a.i) or abiotic mechanisms, leaving large, dark-colored (humic) compounds rich in carbon and nitrogen (N) (Lehmann and Kleber 2015). It was thought that the formation of these compounds renders SOM resistant to further decomposition. However, the hypothesis that chemical recalcitrance (and humification) explains SOM stability has gradually fallen out of favor, due in part to a failure to demonstrate that the secondary synthesis of humic compounds is relevant in natural systems (Keiluweit et al. 2015; Kleber et al. 2015), but also because chemically complex compounds (e.g., lignin) that should theoretically persist in soil can have shorter residence times than supposedly labile molecules (Schmidt et al. 2011). In response, new theories of SOM formation and stability have begun to emerge (Schmidt et al. 2011; Riley et al. 2014; Cotrufo et al. 2015; Lehmann and Kleber 2015; Castellano et al. 2015; Dwivedi et al. 2017a).

Figure 1.

Schematic depicting the (a) traditional and (b) emergent conceptual understanding of soil organic matter (SOM) cycling. The traditional view does not explicitly represent the primary underlying agents and processes (e.g., microbes, aggregation) that drive SOM cycling. In contrast, the emergent view suggests that SOM decomposition is a function of a wide range of ecosystem properties and mechanisms. The key processes depicted in the figure include: (i) exoenzyme production, (ii) decomposition, (iii) organo-mineral interactions, (iv) microbial necromass, (v) ecology of belowground biota, and (vi) root contributions and mobilization of nutrients. Additional details are provided in the text.

Figure 1.

Schematic depicting the (a) traditional and (b) emergent conceptual understanding of soil organic matter (SOM) cycling. The traditional view does not explicitly represent the primary underlying agents and processes (e.g., microbes, aggregation) that drive SOM cycling. In contrast, the emergent view suggests that SOM decomposition is a function of a wide range of ecosystem properties and mechanisms. The key processes depicted in the figure include: (i) exoenzyme production, (ii) decomposition, (iii) organo-mineral interactions, (iv) microbial necromass, (v) ecology of belowground biota, and (vi) root contributions and mobilization of nutrients. Additional details are provided in the text.

Many of these emerging theories put interactions between physical and geochemical conditions and biology at the forefront of SOM stability. Schmidt et al. (2011), for example, concluded that SOM persistence depended on ecosystem properties; they described an array of characteristic environmental traits that contribute to the stability of SOM and sequestration of carbon in soils. These properties include climate, physical properties (e.g., the formation of preferential flow paths within soils that determine the degree of connectivity between microbes and substrates—Hunt and Manzoni 2016), chemical properties (e.g., soil pH or the chemical quality of plant exudates and litter—Castellano et al. 2015), physical protection mechanisms (e.g., the sorption of organic matter onto mineral surfaces, Fig. 1b.iiiKleber et al. 2015), soil redox status, the ecology of belowground biota (including bacteria, archaea, fungi and soil macrofauna, Fig. 1b.v—Bardgett and Van Der Putten 2014), and the metabolic constraints on OM mineralization. In such a network of factors, the diverse biotic components interact chemically and physically with each other, and with the surrounding heterogeneous matrix of minerals and soil aggregates, resulting in the emergence of SOM.

Microbial and plant processes

The majority of SOM is decomposed by microbes and respired into the atmosphere over short (< 50 year) time scales. However, a significant portion of this SOM can become chemically or physically protected—before or after microbial processing—and persist for much longer (100 to 1,000 year) time scales. Several studies have recognized the importance of diverse microorganisms in the formation of SOM (Allison et al. 2010; Knicker 2011; Cotrufo et al. 2013). The Microbial Efficiency-Matrix Stabilization framework, for example, describes SOM dynamics resulting from individual physiologies of diverse microorganisms (Cotrufo et al. 2013). The concept behind the “framework” asserts that plant exudates and litter, while being significant components of SOM themselves (Castellano et al. 2015), are decomposed (Fig. 1a/b.ii), assimilated (Fig. 1b.v), and processed by the microbial community.

Plants play an active role in soil carbon cycling through litterfall and root exudation, thereby influencing microbial growth (Fig. 1b.iv). For instance, roots excrete small-molecular-weight compounds to aid plant acquisition of nutrients and water. In addition, when plants supply small molecular-weight compounds such as sugars or reducible compounds, they prime microbial growth in soil. In general, the ‘priming’ of the soil microbial community can stimulate the production of exoenzymes near the root zone, increasing depolymerization of large biopolymers to monomers, which can become accessible to plant roots (Fig. 1b.i) (Kuzyakov et al. 2000; Gunina and Kuzyakov 2015). Moreover, roots also play a significant role in SOM stability through activities aimed at liberating and mobilizing nutrients (Fig. 1b.vi). Organic acids and chelators can disrupt the ligand bonds between mineral surfaces and organic matter or nutrients. Such actions can subsequently release previously immobilized phosphorus and other inorganic nutrients (Hocking 2001; Wu et al. 2018), and destabilize mineral-associated organic matter (MAOM) (Keiluweit et al. 2015; Jilling et al. 2018).

Overall, microbial processing of decomposed plant exudates and litter results in the microbial “exudation” of different metabolic byproducts, including enzymes, lipids, glycoproteins, extracellular polysaccharides (EPS), and the eventual accumulation of chemically heterogeneous necromass (Fig. 1b.iv) that are the major precursors to stable SOM (Miltner et al. 2012; Cotrufo et al. 2013, 2015; Kögel-Knabner 2017). Indeed, recent evidence suggests that the chemical composition of persistent SOM more closely resembles that of microbial cells and their byproducts than that of plant materials (Kögel-Knabner 2017). Decomposition-resistant SOM often exhibits a narrow carbon-to-nitrogen (carbon:nitrogen) ratio (between 8 to 10) that is more representative of microbial residues than the wide C:N ratio (upwards of 30) characteristic of plant materials (Baisden 2002; Conen et al. 2008). These microbial cells and excretions are often stabilized via organo–mineral and organo–metal–oxide interactions.

Organo–mineral interactions

Physicochemical associations between soil minerals and plant or microbial-derived organic matter can result in the protection of theoretically thermodynamically unstable (or labile) compounds within soil for decades. Soil minerals often constitute a mixture of phyllosilicates, oxides, and hydroxides, as well as short-range order and noncrystalline minerals in some soil types. These particles are generally found in the clay (< 2 μm) size-fraction of soils. The accessibility of SOM to microorganisms can be hindered through mineral interactions by adsorbing or occluding organic compounds and immobilizing exoenzymes responsible for the depolymerization reactions. Recent studies have demonstrated that different classes of microbial biochemicals show variability in the degree to which they adsorb onto mineral surfaces (Jagadamma et al. 2012, 2014). For example, Swenson et al. (2015) and Newcomb et al. (2017) have demonstrated generalizable patterns in the competitive sorption of different functional groups to the surface of ferrihydrite, an amorphous iron oxide. Phosphate-, N-, and aromatic-containing metabolites and carboxylates showed strong adsorptive behavior onto mineral surfaces (Swenson et al. 2015). Conversely, differential biochemical adsorption shapes the composition, stoichiometry (e.g., carbon:nitrogen:phosphorous ratio), and accessibility of the dissolved organic matter (DOM) pool. For example, the immobilization of exoenzymes on mineral surfaces can alter catalytic function by blocking the active site or changing enzyme structure (Quiquampoix and Ratcliffe 1992; Baron et al. 1999; Norde and Giacomelli 2000; Servagent-Noinville et al. 2000). Although adsorbed exoenzymes generally decrease in catalytic capacity, their residence times are higher than free exoenzymes (Yan et al. 2010). These mechanisms, in turn, can shape the structure, function, and activity of the microbial community, which can promote the production of enzymes or organic acids that mineralize available polymers and desorb stable SOM compounds (Kaiser et al. 2014).

Finally, empirical observations and modeling efforts have typically concentrated on surface soils (< 20 cm) for the estimation of carbon storage and turnover rates (Trumbore et al. 1995; Hicks Pries et al. 2018). However, over half of global SOM is stored in subsoils (> 20 cm), and this SOM is significantly older than surface SOM (Jobbágy and Jackson 2000; Mathieu et al. 2015). These patterns suggest a distinct balance of decomposition and preservation processes at depth relative to that at the surface. In all, this emerging picture of MAOM places a premium on understanding the interdependent impacts of soil hydrology, geology, geochemistry, and ecology across the complete soil profile. In particular, it calls for an improved understanding of the feedback between soil properties and the metabolisms of microbes (fungi, bacteria, and archaea). Insight into these emergent properties of soil carbon cycling has in particular benefited from improvements in approaches to study fine-scale MAOM dynamics.

EMERGING TECHNOLOGIES

Emerging technologies are shaping the picture of MAOM dynamics as affected by biology, chemistry, and soil structure. Advances in analytical techniques are revealing the chemistry of minerals and the composition of soil organic matter. Extremely detailed information on the composition of extractable soil organic matter can be obtained by using Fourier Transform Ion Cyclotron Resonance Mass spectrometry (FT-ICR-MS) and quantified by 1H and 13C solution-state Nuclear Magnetic Resonance (Simpson et al. 2001; Tfaily et al. 2015). For nonextractable organic matter, solid-state NMR can characterize 13C labeled organic matter sorbed on surfaces of non-paramagnetic minerals (Clemente et al. 2012). FT-Infrared spectrometry is compatible with living systems and can detect functional group characteristics of MAOM, even when associated with paramagnetic minerals such as iron oxides (Omoike and Chorover 2004; Lehmann et al. 2007). The capacity for imaging soil in nondestructive ways has grown. Advancements in X-ray computed tomography of intact soils yields information regarding the soil structure, pore network and even roots, which can be processed into three-dimensional images (Taina et al. 2011; Mairhofer et al. 2012). More importantly, when combining imaging techniques with chemical analysis, components of MAOM can be spatially resolved. Techniques such as synchrotron radiation scanning transmission X-ray microscopy (STXM) coupled with near edge X-ray absorption fine structure (NEXAFS) has been used to determine the spatial heterogeneity of MAOM in two dimensions (Lehmann et al. 2008). A more detailed literature review on the benefits and limitations of advanced analytical techniques used for MAOM has been presented previously (Kleber et al. 2015).

Microbial traits that affect SOM stability

The emerging understanding of SOM dynamics also calls for greater efforts to identify and quantify the role of microbial functional traits in determining the fate of organic compounds within soils. Microbial communities are complex adaptive systems (Levin 1999), in which heterogeneous collections of individual organisms that coexist, co-operate, and compete across wide-ranging environmental gradients shape the evolution and emergence of community structure. The fitness of individuals within a community is governed by their inherent physiological traits and the energetic costs of intracellular trade-offs (Maharjan et al. 2013). Microbial traits can be defined as physiological or morphological characteristics closely linked with ecosystem function and performance (Viollea et al. 2014). These traits can be broadly separated into response (determining how organisms respond to environmental change, such as drought) and effect (determining an organism's effect on ecosystem processes, such as nitrification) traits (Allison 2012). The expression of, and resource allocation to, different traits (e.g., exoenzyme production) is regulated by intracellular energetic trade-offs (Beardmore et al. 2011). These factors determine an individual microbe's ability to compete for resources within a wider community, and to tolerate and thrive under extraneous selective pressures, including fluctuating environmental conditions. Traits and trade-offs that promote optimal relative fitness will subsequently be selected for and fixed (at least temporarily) within a population (Razinkov et al. 2013).

Many relevant traits (e.g., osmolyte production or exoenzyme production) show deep evolutionary divergence and nonrandom distributions (Schimel and Schaeffer 2012), meaning that environmental selection for specific ecological strategies can lead to concomitant changes in a suite of genome-linked traits related to SOM formation and degradation. Reorganization of microbial community composition can also change anabolic trait distribution and select for traits that convert plant-derived carbon to metabolites (e.g., osmolytes), or necromass that may be preferentially stabilized on soil minerals (Liang et al. 2017). On a broader scale, community changes that result in an increase in the soil fungal:bacterial ratio can also result in higher SOM due to the higher biomass associated with the proliferation of fungi (Strickland and Rousk 2010). The extent of bacterial secondary metabolism (e.g., antibiotics, non-structural carbohydrates, siderophores, etc.), which plays a large role in shaping the chemical composition of microbial exudates and necromass, is beginning to be quantified through advances in genomic and metabolomic technologies. For example, Crits-Christoph et al. (2018) recently used large-scale metagenomic sequencing to reconstruct the genomes of a number of common soil bacteria. They uncovered the potential for extensive secondary metabolite production, suggesting a largely underestimated capacity of these bacteria to produce and exude biosynthetic products and contribute to the composition of SOM. However, while “omic” techniques provide greater understanding of the complexity of the microbial community and its metabolic diversity, the need for large amounts of material for nucleic acid extraction (covering multiple niches) and the destructive nature of that extraction, limit understanding of temporal dynamics. These dynamics include how community-composition interactions necessitate the development and application of microbial centric models that permit testing of hypotheses related to cell-cell interactions, response to environmental fluctuations, and the microbial feedback to carbon dynamics (Riley et al. 2014; Georgiou et al. 2017).

While these small-scale processes undoubtedly interact, giving rise to larger-scale phenomena, a tractable balance of simplicity and complexity must be sought in the development of parsimonious, yet sufficiently accurate, macro-scale models of SOM cycling. To this end, it is critical to thoroughly understand key SOM decomposition and stabilization mechanisms, where and when they are important, and how nonlinear interactions between components may compound. Changes in climate and vegetation may also cause important shifts in the relative role of key mechanisms, with broad implications for the terrestrial carbon cycle.

MODELING SOM DYNAMICS

There is a suite of SOM dynamics models incorporated into land models ranging from simple pool-based CENTURY-like representations (Parton et al. 1987, 1998, 2010; Jenkinson et al. 2008) to more complex models such as the Biotic and Abiotic Model of SOM (BAMS; Riley et al. 2014; Dwivedi et al. 2017; Tang et al. 2019) and the COntinuous representation of SOC in the organic layer and the mineral soil, Microbial Interactions and Sorptive StabilizatION (COMISSION; Ahrens et al. 2015), which explicitly represents physicochemical processes and microbially mediated reactions. Here we briefly describe the traditional and some of the emerging modeling approaches. More importantly, this section aims to describe differences in the traditional and emerging modeling approaches, and the need for a reactive transport modeling framework representing relevant mechanisms.

Traditional modeling approaches

Traditionally, mean residence times (MRT) are used to categorize soil organic carbon into different “pools”. The decomposition reaction rate is the inverse of MRT, and therefore pools with shorter MRT have higher inherent decomposition reactivity (Parton et al. 1987; Davidson and Janssens 2006; Kleber 2010). These pools have been extensively used to investigate SOM decomposition in models such as CENTURY, ROTH-C, and their progeny models. These pools are intended to represent from the most labile to the most recalcitrant carbon in CENTURY-Like models (Parton et al. 1987, 2010), and from humified organic matter to inert in ROTH-C models (Jenkinson and Coleman 2008; Jenkinson et al. 2008). Similarly, SOM decomposition rates are also empirically referred to as fast, slow, and passive for pools ranging from the most labile, or humified, to the most recalcitrant, or inert. These traditional approaches have been very widely applied, but lack representation of the underlying microbial and abiotic processes described above. They have also been shown to have very high parameter equifinality, in spite of being calibratable in many systems (Luo et al. 2015, 2017).

The traditional CENTURY-like models are formulated by assuming sufficient microbial activity to support organic matter decomposition, because litter bag and soil incubation experiments often observe consistent patterns of exponential organic matter decomposition over time (Davidson et al. 1995; Berg and McClaugherty 2014). Soil minerals are considered to exert only a multiplicative effect on the turnover rates of the different components of SOM. Other multiplicative factors are similarly applied to represent how temperature and moisture regulate SOM decomposition. However, spectrum analysis of field measurements such as concentration-gradient-based soil surface CO2 efflux time series often show that there are many dynamic signals that cannot be resolved with traditional models (e.g., Vargas et al. 2010; Manzoni et al. 2018). Notably, the traditional models fail to capture the priming effect in soil, which involves strong interactions between microbes and their bioproducts, substrates, and soil minerals (e.g., Fontaine and Barot 2005). Furthermore, when organo–mineral interactions are involved, substrate bio-availability usually decreases. This observation suggests that in soils where soil minerals are plentiful, oligotrophic microbes are favored. However, model representation of oligotrophic microbes may not be necessary, considering that microbes can live in a dormant state, making it difficult to determine whether oligotrophic microbes are the major contributors to observed microbial activity or simply whether a smaller population of more copiotrophic organisms are active (Joergensen and Wichern 2018).

Although these simplified pool-based models have been important for estimating SOM stocks, they are unable to capture SOM dynamics beyond first order effects. For example, most site- and global-scale terrestrial SOM models apply first-order relationships, and therefore cannot predict explicit turnover times for different compounds. Further, a major obstacle remains in matching the model with data in terms of understanding which components of SOM are decomposed and which are stabilized, because the pools in the traditional models are not structurally well defined (i.e., not observable; Abramoff et al. 2018). Most importantly, these traditional modeling approaches cannot answer why organic matter that is thermodynamically unstable persists in soils, sometimes for a millennium (Schmidt et al. 2011). In view of these modeling inconsistencies, new models have been developed and applied to understand the effects of microbial processes and organo–mineral interactions on SOM dynamics.

Transitional modeling approaches

Interactions between microbes and soil minerals can be explicitly represented in more mechanistic models—some of the earlier attempts can be found in Smith (1979) and Grant et al. (1993). The recently published batch (i.e., no vertical resolution) models (e.g., MEND; Wang et al. 2013, RESOM; Tang and Riley 2014, CORPSE; Sulman et al. 2014, and MIMICS; Wieder et al. 2014) use diverse approaches to formulate organo–mineral interactions. There are also several new models (described below in Section ‘Current use of reactive transport modeling of SOM dynamics’) that explicitly represent microbes, soil mineral interactions, and vertical structure.

The MEND model incorporates a mineral-adsorbed soil carbon pool affected by separate adsorption and desorption kinetics that result in a Langmuir isotherm at equilibrium. However, Wang et al. (2013) argued that by incorporating microbial dormancy into the MEND modeling framework, only a single microbial functional group is required to reproduce laboratory aerobic incubation measurements of biomass and CO2 respiration rates. A similar approach is adopted by the CORPSE model. RESOM uses the equilibrium form of the Langmuir isotherm, which is then integrated with the Equilibrium Chemistry Approximation kinetics (Tang and Riley 2013; Tang 2015; Zhu et al. 2017) to represent interactions between dissolvable carbon, mineral surfaces, exoenzymes, and microbes. The MIMICS model avoids an explicit formulation of organo–mineral interactions, and instead lumps these interactions by changing the half saturation constants used to drive Michaelis-Menten type microbial kinetics. These models were evaluated together in a recent model intercomparison (Sulman et al. 2018); that study concluded that the observational data were insufficient to invalidate any of the model structures. The authors suggested that long-term observations of soil carbon input (from plants) and perturbation experiments (such as warming) are needed to better evaluate and constrain these newly developed models.

Next-generation process representation approaches

Soil is a complex system and considering its physical and chemical evolution is critical for accurately capturing SOM dynamics (Riley et al. 2019). The following discussion is intended to review modeling needs for the following important topics: (1) microbial dynamics, (2) mineral associated organic matter (MAOM) interactions, and (3) the effects of SOM molecular structure. We subsequently make recommendations as to how these processes can be integrated into a reactive transport model.

Microbial Processes.

The structure, function, and response to perturbation of the soil microbial community plays an important role in shaping soil biogeochemistry and SOM dynamics (see Section ‘Microbial and plant processes’ above), motivating the active representation of the belowground community in next-generation model structures. Allison (2012), using the DEMENT model to explore how microbial traits would covary with the evolution of decomposing litter, found that incorporating microbes with different growth traits enabled the model to explain more than 60% of the variance in decomposition rates of 15 Hawaii litter types. Using a similar modeling strategy, Kaiser et al. (2014) found that microbial diversity is essential in explaining the stoichiometric evolution in decomposing litters. However, relatively few studies explored the influence of microbial diversity on carbon formation and decomposition in the mineral soil. When this is done, microbial community structure is often aggregated into the r versus K strategists that are adopted from the MacArthur and Wilson (1967) classification scheme. For instance, the MIMICS model (Wieder et al. 2014) used a copiotrophic group to represent the r-strategists and an oligotrophic group to represent the K-strategists. This r vs K representation of microbial community structure has been argued to be critical to accurately simulating the priming effect (Fontaine et al. 2003). However, models like MEND, CORPSE, and RESOM are able to simulate priming with only one microbial group, as long as substrate diversity is included in the model. In addition, Lawrence et al. (2009) showed that short-term SOM dynamics resulting from interactions between microbes and environmental conditions (e.g., drying–rewetting events) could not be adequately represented using a first-order representation of decomposition (as in Century or RothC). They suggested that while first-order models can capture the bulk response of SOM dynamics under steady-state conditions, the inclusion of exoenzyme and microbial controls on decomposition is necessary to appropriately simulate pulsed rewetting dynamics.

Mineral Associated Organic Matter (MAOM).

Traditional models that employ a function of soil texture to represent soil mineral effects on soil carbon dynamics, as a multiplicative modifier to decomposition rates, are likely underestimating the roles played by soil minerals. For example, soil minerals shift between different oxidation status due to particular electron orbital configurations, which enable them to act as either reductants or oxidants of organic matter. A number of microorganisms are able to take advantage of this chemical flexibility and couple these soil minerals with organic carbon to sequester energy when the most favorable oxidant, i.e., oxygen, is in short supply. Likewise, none of the models mentioned above sufficiently accounts for organo–mineral interactions by only considering sorption–desorption dynamics.

Soil minerals are important factors in determining soil aqueous chemistry. Minerals, microbial metabolic byproducts (e.g., organic acids), and inorganic compounds exert strong controls on soil pH and redox potential. However, few contemporary soil carbon models have taken such organo–mineral interactions into account. Rather, when they are needed, such as in modeling methane dynamics, pH and redox potential are formulated as multiplier functions without accounting for their mechanistic evolutions (Zhang et al. 2002; Li et al. 2004; Zhuang et al. 2004; Riley et al. 2011). A slightly more mechanistically based model is attempted in Tang et al. (2016), where they merged the Windermere Humic Aqueous Model (WHAM; Tipping 1994) with the CLM-CN model to simulate the effect of iron cycling on CO2 and CH4 production in anoxic Arctic soil microcosms. Their model also simulated the evolution of pH using the aqueous chemistry capability provided by WHAM. However, the decomposition dynamics follow the traditional century-like formulation, so the model is therefore still not mechanistically representing organo–mineral interactions. Zheng et al. (2019) recently updated the framework from Tang et al (2016) with the PHREEQC 3.0 model (Charlton and Parkhurst 2011) and applied the resultant code to incubated organic soils from an Arctic tundra site. Although their simulations indicated some favorable comparisons with CO2 and CH4 fluxes, the representation of organo–mineral interactions was not improved.

The most comprehensive model integrating soil pH, redox, and sorption-desorption aspects of organo–mineral interaction is likely the ecosys model (e.g., Grant et al. 2003, 2017a, b), which includes a comprehensive parameterization of aqueous electrolyte dynamics and their relations with soil pH and ion-exchange capacity (Grant and Heaney 2010). However, this parameterization was primarily used to model inorganic phosphorus dynamics. As such, the role of soil minerals (e.g., ferrihydrite) as alternative electron acceptors for organic carbon decomposition is not modeled.

In summary, although much of the necessary chemistry knowledge exists, a modeling framework that systematically evaluates how organo–mineral interactions will play out in various aspects of soil biogeochemistry has not been developed (but see Riley et al. 2019 for a possible path forward in this regard). Further, more knowledge on how organo–mineral interactions affect biological aspects of microbial dynamics is needed. For instance, there is not a good method applicable to large-scale models to represent how organo–mineral interactions contribute to the formation of biofilms at different places within the soil, even though such structures are of great significance in determining how microbes deal with fluctuations in soil water (Or et al. 2007).

Molecular Structures.

Most land models use CENTURY-Like SOM pools to represent SOM dynamics by applying pseudo-first-order approximations (Jenkinson et al. 2008; Parton et al. 2010). Although these simplified representations of SOM compounds are informative, some investigators have argued that the structures of diverse organic matter compounds may need to be represented to capture SOM dynamics, for several reasons. First, the Century-like pools are largely qualitative and vague, thereby leading to large uncertainty when confronting predictions of SOM dynamics with observations. In general, it is not possible to precisely categorize pools of SOM by turnover times that have such a broad range (e.g., a few days to thousands of years). In nature, the formation or decomposition of SOM compounds are properties of the ecosystem (Schmidt et al. 2011), and the traditional qualitative structural classification is a major obstacle toward providing a mechanistic basis for SOM dynamics. For instance, a group of SOM compounds may have a completely different response to soil moisture and temperature than another group of SOM compounds, even though these distinct groups may have comparable turnover times. Moreover, organo–mineral interactions may also show contrasting behaviors; for example, a group of compounds falling into comparable intrinsic decomposition rates may have different affinities to soil minerals or responses to environmental conditions (Gordon and Millero 1985; Gu et al. 1994; Mikutta et al. 2007; Kleber et al. 2011). Overall, we argue that there is a need to represent SOM compounds using their molecular structures in reactive transport models along with physical protection mechanisms (e.g., MAOM) and different functional groups of microbes (e.g., Riley et al. 2014; Dwivedi et al. 2017; Tang et al. 2019.)

As discussed in the Introduction and the sction ‘Current Use of Reactive Transport Modeling of SOM Dynamic’, SOM decomposition is governed by environmental conditions such as physical heterogeneity, physical disconnection, plant inputs, microbial diversity, and microbial activity, as well as the molecular structure of organic matter (Schmidt et al. 2011; Riley et al. 2014; Dwivedi et al. 2017a). To capture such complex physico-chemical and biological interactions, we propose that a vertically resolved reactive transport modeling framework is required.

CURRENT USE OF REACTIVE TRANSPORT MODELING OF SOM DYNAMICS

Flow and reactive transport processes influence SOM formation and decomposition. In particular, geochemical processes can affect the mobility of chemical species through different mechanisms of dissolution-precipitation, sorption-desorption, ion-exchange, redox-reactions, complexation, and colloidal interactions (Arora et al. 2016, 2018; Yabusaki et al. 2017; Dwivedi et al. 2018b, a). There are a variety of RTMs and numerical codes—such as TOUGHReact (Xu et al. 2006; Maggi et al. 2008; Gu et al. 2009), PFLOTRAN (Hammond et al. 2014), CRUNCH (Steefel et al. 2015), ecosys (Grant 2013), BAMS (Riley et al. 2014; Dwivedi et al. 2017a; Tang et al. 2019), and BeTR (Tang et al. 2013; Tang and Riley 2018)—that are available to describe and can represent the interaction of various complex and competing biogeochemical processes across spatial and time scales. These numerical codes can represent aqueous, gaseous, and sorbed phases, along with a large suite of processes—including hydrology, soil energy dynamics, advection and diffusion, aqueous speciation, minerals precipitation and dissolution, microbial dynamics, adsorption and desorption, and equilibrium and nonequilibrium chemical reactions.

Several studies have demonstrated the use of reactive transport models to address SOM decomposition and stabilization (Grant et al. 2003; Ahrens et al. 2015; Dwivedi et al. 2017a; Georgiou et al. 2018). These models explicitly represent physical protection mechanisms and different functional groups of microbes. BAMS also represented a reduced number of SOM molecular structures (Riley et al. 2014; Dwivedi et al. 2017; Tang et al. 2019). Here, we provide a brief description of how these processes are represented in reactive transport models.

Organo–mineral interactions are represented in a fashion similar to the isotherm-based approach in the single layer batch model. For instance, BAMS1 (Riley et al. 2014) adopted linear kinetics to model sorption and desorption processes, whereas in BAMS2 (Dwivedi et al. 2017a) the formulation was replaced with the Surface Complexation Model (SCM), allowing variable biogeochemical conditions within a thermodynamic framework. Both BAMS1 and BAMS2 were able to reproduce the almost exponentially decreasing profile of soil carbon content and increasing age of soil radiocarbon with depth. The COMISSION model (Ahrens et al. 2015) similarly considers interactions between microbes, minerals, and substrates, using Langmuir kinetics and reactive transport modeling, and was able to obtain analogous results. In contrast, ecosys uses Freundlich kinetics to model the sorption dynamics of soluble carbon to soil mineral surfaces and unhydrolyzed organic matter (Grant et al. 2003). Of all these models, ecosys is the only one that also considers a full range of terrestrial ecosystem properties; it has been shown in many studies across various ecosystems to be able to capture much of the observed temporal variability in net ecosystem carbon exchanges (e.g., Grant et al. 2003, 2006, 2017a, b). Overall, these modeling exercises suggest transport combined with soil–mineral-regulated SOM stabilization and microbial decomposition help explain observed soil carbon content and age profiles. Furthermore, the specific mathematical form used to represent organo–mineral interaction plays a minor role, given the limited data available to evaluate the models.

Riley et al. (2014) made the first attempt to explicitly represent molecular structures of SOM compounds and integrate processes into a reactive transport model. Because SOM compounds consist of several plant-synthesized, degraded compounds, and microbial biomass as well as necromass, it is difficult to include a full representation of each species. Therefore, they grouped SOM compounds based on different metrics such as O:C ratio, charges (positive or negative), and degree of polarity. Riley et al. (2014) include above- and belowground litter, as well as root exudates for carbon inputs. The litter was degraded into several simpler organic compound groups such as cellulose, hemicellulose, lignin, monosaccharides, phenols, amino acids, lipids, and nucleotides. Root exudates were considered to be simpler organics, and were included among monomers (i.e., monosaccharides, amino acids, organic acids, and lipids pools). Subsequently, Dwivedi et al. (2017) modified this reaction network by collapsing all the monomers into a single pool. Although observations of individual SOM compounds were not available for comparison with predictions, individual SOM compounds highlighted the processes predicted to control SOM stocks and cycling.

To a certain extent, these reactive transport models demonstrated the value of including different microbial functional groups. BAMS included fungi and aerobic bacteria functional groups by assigning affinities to decompose plant litter and SOM (Neely et al. 1991; Romaní et al. 2006; Thevenot et al. 2010; DeAngelis et al. 2013). Subsequently, dissolved organic carbon (DOC) and peptidoglycan (i.e., necromass or dead cell wall material (DCWM); Frostegard and Baath 1996) were produced upon death of these microbial groups. Even though substantially underrepresenting functional diversity in soils, these models indicated key processes and pulse responses, and the value of examining SOM dynamics with this type of approach. It is important to realize that the governing ecological functions in the ecosystem are based upon biodiversity, and it is very data-intensive to decipher linkages between ecosystem function and biodiversity (Goldfarb et al. 2011).

There are some models that have gone beyond representing microbial diversity by considering soil fauna in the decomposition dynamics of SOM. For instance, Hunt et al. (1987) modeled the detrital food-web in a grassland by representing 15 groups of organisms, ranging from bacteria and fungi to predacious mites. The simulations show that trophic interactions are essential in maintaining the ecosystem wide nitrogen, echoing the message by Kaiser et al. (2014) that biodiversity is important for enabling nutrient-based regulation of the carbon cycle. Considering that SOM decomposition always involves a spectrum of bioreactions that include nitrification, denitrification, methanogenesis, and methanotrophy, it may be important to consider microbial diversity when all these different bioreaction pathways and different gas fluxes are to be modeled. For a carbon-only model, before more comprehensive data are available to benchmark model evaluations, it is hard to justify whether microbial community structure should be included in the model formulation.

Below, we describe the most important governing equations for simulating reactive flow and transport in the continuum subsurface environment.

Variably saturated flow

Multiphase flow (e.g., liquid and gas) or the Richards equation can be used to describe variably saturated flow. The Richards equation assumes only one active phase (e.g., liquid), while the gas phase is assumed to be passive.

Richards equation.

The Richards equation (Richards 1931) has been widely used in the literature for describing partially saturated flow assuming a passive air phase present at atmospheric pressure (e.g., Neuman 1973; Panday et al. 1993; Dwivedi et al. 2016, 2017b). The saturation-based formulation of the Richards equation is given as:

SaSsht+ϕSat=[kraKh]+Qa
(1)

where ∇ denotes the divergence operator; h [m] is the hydraulic head; K [m s−1] is the hydraulic conductivity tensor; t [s] is the time; ϕ [m3 void m−3 medium] is a dimensionless quantity, which represents the porosity of the porous media; Sa [m3 H2O m−3 void] and dimensionless quantity kra [–] denote the saturation of the aqueous phase and the relative permeability, respectively; Ss [m−1] is the specific storage coefficient; Qa [m3 H2O m−3 medium s−1] is the volumetric source or sink term.

The corresponding pressure-based formulation is given as:

[Saϕρf]t=[ρfkraksatμ(Pρfgez)]+ρfQa
(2)

ρf [kg m−3 water] represents the fluid density; ∇P [kg m−1 s−2, or Pa m−1] represents the gradient of the fluid pressure; g [m s−2] denotes the acceleration due to gravity; ez is a unit vector in the vertical direction; ksat [m2] is the permeability tensor for fully saturated conditions.

Multiphase Flow.

The occurrence of components like CO2 or H2O in various phases (gas, liquid, solid) can be described using multiphase flow. Multiphase flow is described by solving the mass conservation equation of multiphase transport of each component j in phase α:

[ϕΣαραSαYjα]t=[ΣραYjαkraksatμ(Pραgez)]+ραQj
(3)

where Yja is the mass fraction of component j in phase α; m [kg s−1 m−1] is the dynamic viscosity; and all other variables (e.g., kra, etc.) are same as defined earlier.

There are several models in the literature that can be used to describe the relationship between the aqueous phase saturation and capillary pressure. However, Brooks-Corey and van Genuchten formulations are the most well-known. In addition, the Mualem and Burdine formulations are typically used to describe relative permeability functions (Brooks and Corey 1964; Mualem 1976; van Genuchten 1980). We encourage readers to explore user's manuals of numerical codes such as TOUGHReact, PFLOTRAN, or Crunch for more details on these models.

Reactive transport equations

To simulate reactive transport, the transport and biogeochemical reactions are combined in the continuity equation as follows:

(ϕSLCi)t=(ϕSLDi*Ci)(qCi)Σr=1NrνirRrΣm=1NmνimRmΣg=1NsνilRl
(4)

The left-hand side of the Equation (4) denotes the accumulation term (mol m−3 medium s−1), which is the product of the porosity (φ) and liquid saturation (SL), and the concentration (Ci). On the right-hand side, the first and second term together describe advective-diffusive transport; the third, fourth, and fifth terms represent various reactions that are partitioned between aqueous-phase reactions (Rr), mineral reactions (Rm), and gas reactions (Rl), respectively. Furthermore, a stoichiometric relationship is used to compute transformation from reactant to product species. In Equation (4)Di* denotes the diffusion coefficient, which is specific to chemical species considered as indicated by the subscript i; nij are the stoichiometric coefficients of reactant j in reaction i, whereas there are Nr, Nm, Ng reactants in aqueous, mineral, and gas phases, respectively.

These reactions can further be classified into equilibrium-based or kinetic. Equilibrium-based reactions typically require the assumption that there exists a local equilibrium between reactants and products. However, rate-based reactions (kinetic) are always the more general (Steefel and Lasaga 1994). SOM decomposition involves microbial reactions, which we describe below.

Biological reaction rates

The microbial decomposition rate of substrate Ci can be described using Michaelis-Menten (Michaelis and Menten 1913) kinetics:

dCidt=−μiCiKi+CiO2Ko2+O2BY
(5)

where μi (s−1) is the maximum specific consumption rate of substrate i (s−1); O2 (mol O2 m−3) is the aqueous O2 concentration; and B (mg C-wet-biomass L−1) is the wet biomass carbon of microbial functional groups considered in the reaction network (e.g., bacteria and fungi). The effects of pH, saturation levels of soil, and temperature (T) can be important for decomposition rates (e.g., Schimel et al. 2011). Equation (5) can be modified to using functions f1(θ), f2(pH), and f3(T) to describe these environmental effects on microbial activity (Maggi et al. 2008), as follows:

dCidt=−μiCiKi+CiO2Ko2+O2BYf1(θ)f2(pH)f3(T)
(6)

While the Monod-Michaelis-Menten kinetics, including the related multi-substrate progeny, have enjoyed a wide range of applications, Tang and Riley (2013) showed that the Monod kinetics fail to account for the substrate limitation in their approximation to the law of mass action, and misplace such limitation as a linear competition by juxtaposing many Monod terms (e.g., for Eqn. (5), the Monod representation of competition may lead to simple addition of more Monod terms without accounting for one microbe's influence on another microbe's K parameter (Tang 2015; Tang and Riley 2017). This shortcoming will result in poor model performances under conditions of substrate limitation, a situation that is characteristic of soil. Further, the Monod kinetics will oversimplify the connections between substrates and consumers when describing a multi-substrate—multi-consumer network, resulting in the modeled fluxes being oversensitive to the kinetic parameters that would be formulated in the corresponding law of action model. The Equilibrium Chemistry Approximation kinetics (ECA) and its extended SUPECA (synthesizing unit plus ECA) kinetics are both able to adequately address these two shortcomings of the Monod kinetics, and also provide a very straightforward way to link the kinetic parameters with thermodynamics and biological traits (Tang and Riley 2013, 2017; Zhu et al. 2016a, b).

In most formulations of microbial growth in reactive transport models, substrate use efficiency is the only parameter to describe the partition of an assimilated substrate into either biomass or respired product. However, the catabolic and anabolic separation is more complex than this simple representation. When one takes into account building up metabolic reserves, cell maintenance and growth, and the enzyme and polymer exudation involved in microbial physiology, the emergent CUE can have strong nonlinear variations with respect to the environmental conditions, such as temperature or moisture (Tang and Riley 2014; Allison and Goulden 2017). CUE may be further modified by regulations of nutrient availability (Manzoni et al. 2018). All these conclusions suggest a more nuanced parameterization of substrate use is needed in future modeling of microbial-activity-regulated SOM dynamics.

Mineral reactions

SOM dynamics are influenced by minerals present in the geochemical system. For example, pH affects SOM reactions, and the presence or absence of calcite minerals may therefore affect SOM dynamics. Similarly, iron minerals (e.g., ferrihydrite) can be important for SOM dynamics, as a significant fraction of the aqueous phase of SOM can undergo sorption and stay protected from microbial reactions (Dzombak and Morel 1996; Dwivedi et al. 2017a). Mineral precipitation and dissolution reactions can be described using Transition State Theory (or TST) type rate laws:

Rm=(kneutral+kH+[aH+]+ΣkjΠi[aij])(1Qm/Keq.m)
(7)

where rate constants for neutral, acid or additional (jth) reaction mechanisms are denoted by kneutral, kH+ and kj, respectively; aij indicates the activity of the ith aqueous species in the jth reaction; finally, Qm denotes the ion activity product of the mth mineral phase and Keq, m denotes its corresponding equilibrium constant.

Sorption reactions

Sorption is a complex process, one that depends upon organic molecule characteristics, surface area, and site density of minerals, as well as aqueous chemistry (Dudal and Gérard 2004). There are several approaches that can be used to describe sorption processes and MAOM. However, kinetic and equilibrium adsorption isotherms (e.g., linear, Langmuir, Freundlich; Davis and Kent 1990) have been widely used in the literature to describe SOM dynamics (Grant et al. 2011; Mayes et al. 2012; Riley et al. 2014; Ahrens et al. 2015). Typically, sorption isotherms describe a functional relationship between adsorbate and adsorbent for a constant-temperature condition. Correspondingly, empirical sorption isotherms have been developed to describe both kinetic (i.e., sorption site is slow and assumed to achieve adsorption equilibrium described using forward and reverse rates) and equilibrium (i.e., sorption site is assumed to reach adsorption equilibrium fast) systems. Because these empirical formulations cannot account for the out-of-sample effects of chemical conditions for which they are derived, thermodynamic surface complexation models (SCM) were also used (Dwivedi et al., 2017). A simple linear sorption relationship can be imposed using forward (adsorption; k (s−1)) and reverse (desorption; kr (s−1)) rates. In the absence of any competing source or sink of a specific SOM species, an effective equilibrium linear sorption relationship will become as follows:

Kd=kfkr
(8)

Nonlinear kinetic models require estimates of these rate constants (forward and reverse) and an additional exponential parameter (Goldfarb et al. 2011). On the other hand, SCMs offer a mechanistic approach for representing sorption processes and are generally more robust and applicable over variable geochemical conditions (Dzombak 1990; Goldfarb et al. 2011). The SCM models describe the sorption of solutes on solid surfaces as a chemical reaction between aqueous species and surface sites (surface complexation) (Dzombak 1990).

A+SAS
(9)

where A and S are the sorbate molecule and sorbent site (surface), respectively; and AS is the sorbed complex. The sorption reactions are described by a mass action law equation assuming a local equilibrium:

KEq=[AS][A][S]
(10)

where KEq denotes the equilibrium sorption constant; [AS] is the concentration of surface complexes (mol kg−1water); [A] denotes the concentration of sorbate molecule (mol kg−1water), and [S] denotes the sorbent site (mol kg−1water).

Although we are not aware of any study considering rate-limited and equilibrium-based models together to describe sorption processes, we argue that rate-limited and equilibrium-based sorption do co-occur in soils. We suggest future modeling approaches consider rate-limited and equilibrium-based sorption to determine the relative importance of these processes for SOM dynamics.

SUMMARY

Soils store vast amounts of terrestrial organic carbon, more than the atmosphere and terrestrial vegetation combined. This carbon is vulnerable to release to the atmosphere under a changing climate. Although the mechanisms leading to the decomposition of SOM are well documented, uncertainties persist regarding the stability of SOM. To address this critical gap, a robust predictive understanding and modeling of SOM dynamics is essential for examining short- and long-term changes in soil carbon storage and feedbacks with climate. In this chapter, we reviewed recent research that improves emergent understanding of the important factors contributing to SOM stability. While there currently exists a suite of models representing SOM dynamics that span a range of complexity, some recent mechanistic models are more consistent with this emergent understanding of SOM persistence. Yet even those recent models do not represent several processes that can be important for SOM dynamics. We conclude that the next-generation models need to represent the full spectrum of quantitatively important mechanisms for determining SOM persistence—including rate-limited and equilibrium-based sorption, the formation of soil aggregates, representative soil minerals, microbial community dynamics, and vegetation interactions—to accurately predict short- and long-term SOM dynamics. Because this recommendation is obviously challenging, we have assembled an open-source, reactive-transport based SOM model that can be used to robustly integrate many of these processes (called BeTR-S; https://github.com/BeTR-biogeochemistry-modeling/sbetr) and invite the community to experiment with it (Riley et al. 2019).

Overall, it is important to understand SOM dynamics because SOM decomposition gives rise to potential greenhouse gases. Microbial dynamics, MAOM, and the molecular structure of SOM compounds have implications for understanding and modeling the short- and long-term responses of soil carbon stocks under local and regional climatic perturbations. Finally, this chapter illustrates the need to evaluate SOM dynamics in a reactive transport modeling framework that includes ecosystem properties.

Acknowledgment

This material is based upon work supported by the U.S. Department of Energy, Office of Science, Office of Biological and Environmental Research (as part of the Watershed Function Scientific Focus Area), and Office of Science, Office of Advanced Scientific Computing (as part of the project “Deduce: Distributed Dynamic Data Analytics Infrastructure for Collaborative Environments”) under Contract No. DE-AC02-05CH11231. Jinyun Tang acknowledges support from the Next Generation Ecosystem Experiment (NGEE-Arctic). The USDA NIFA Postdoctoral Fellowship program supported Katerina Georgiou. Nicholas Bouskill acknowledges support from a DOE Early Career Research Project (# FP00005182). William Riley acknowledges support from the Terrestrial Ecosystem Science Scientific Focus Area of Berkeley Lab. We thank Diana Swantek of Berkeley Lab for assistance with preparing Figure 1.

REFERENCES

1.
Abramoff
R
,
Xu
X
,
Hartman
M
,
O'Brien
S
,
Feng
W
,
Davidson
E
,
Finzi
A
,
Moorhead
D
,
Schimel
J
,
Torn
M
,
Mayes
MA
(
2018
)
The Millennial model: in search of measurable pools and transformations for modeling soil carbon in the new century
.
Biogeochemistry
 
137
:
51
71
2.
Ahrens
B
,
Braakhekke
M
,
Guggenberger
G
,
Schrumpf
M
,
Reichstein
M
(
2015
)
Contribution of sorption, DOC transport and microbial interactions to the 14C age of a soil organic carbon profile: Insights from a calibrated process model
.
Soil Biol Biochem
 
88
:
390
402
3.
Allison
SD
(
2012
)
A trait-based approach for modelling microbial litter decomposition
.
Ecol Lett
 
15
:
1058
1070
4.
Allison
SD
,
Goulden
ML
(
2017
)
Consequences of drought tolerance traits for microbial decomposition in the DEMENT model
.
Soil Biol Biochem
 
7
:
104
113
5.
Allison
SD
,
Wallenstein
MD
,
Bradford
MA
(
2010
)
Soil-carbon response to warming dependent on microbial physiology
.
Nat Geosci
 
3
:
336
340
6.
Arora
B
,
Spycher
NF
,
Steefel
CI
,
Molins
S
,
Bill
M
,
Conrad
ME
,
Dong
W
,
Faybishenko
B
,
Tokunaga
TK
,
Wan
J
,
Williams
KH
,
Yabusaki
SB
(
2016
)
Influence of hydrological, biogeochemical and temperature transients on subsurface carbon fluxes in a flood plain environment
.
Biogeochemistry
 
127
:
367
396
7.
Arora
B
,
Davis
JA
,
Spycher
NF
,
Dong
W
,
Wainwright
HM
(
2018
)
Comparison of electrostatic and non-electrostatic models for U (VI) sorption on aquifer sediments
.
Ground Water
 
56
:
73
86
8.
Baisden
WT
(
2002
)
Turnover and storage of C and N in five density fractions from California annual grassland surface soils
.
Global Biogeochem. Cycles
 
16
:
1117
9.
Bardgett
RD
,
Van Der Putten
WH
(
2014
)
Belowground biodiversity and ecosystem functioning
.
Nature
 
515
:
505
511
10.
Baron
MH
,
Revault
M
,
Servagent-Noinville
S
,
Abadie
J
,
Quiquampoix
H
(
1999
)
Chymotrypsin adsorption on montmorillonite: Enzymatic activity and kinetic FTIR structural analysis
.
J Colloid Interface Sci
 
214
:
319
332
11.
Beardmore
RE
,
Gudelj
I
,
Lipson
DA
,
Hurst
LD
(
2011
)
Metabolic trade-offs and the maintenance of the fittest and the flattest
.
Nature
 
472
:
342
346
12.
Berg
B
,
McClaugherty
C
(
2014
)
Plant litter: Decomposition, Humus Formation, Carbon Sequestration
.
Springer-Verlag
Berlin Heidelberg
13.
Brooks
R
,
Corey
A
(
1964
)
Hydraulic properties of porous media
.
Hydrol Pap Color State Univ
 
3
:
37
pp
14.
Castellano
MJ
,
Mueller
KE
,
Olk
DC
,
Sawyer
JE
,
Six
J
(
2015
)
Integrating plant litter quality, soil organic matter stabilization and the carbon saturation concept
.
Glob Change Biol
 
21
:
3200
3209
15.
Charlton
SR
,
Parkhurst
DL
(
2011
)
Modules based on the geochemical model PHREEQC for use in scripting and programming languages
.
Comput Geosci
 
37
:
1653
1663
16.
Clemente
JS
,
Gregorich
EG
,
Simpson
AJ
,
Kumar
R
,
Courtier-Murias
D
,
Simpson
MJ
(
2012
)
Comparison of nuclear magnetic resonance methods for the analysis of organic matter composition from soil density and particle fractions
.
Environ Chem
 
9
:
97
107
17.
Conen
F
,
Zimmermann
M
,
Leifeld
J
,
Seth
B
,
Alewell
C
(
2008
)
Relative stability of soil carbon revealed by shifts in δ 15N and C:N ratio
.
Biogeosciences
 
5
:
123
128
18.
Cotrufo
MF
,
Soong
JL
,
Horton
AJ
,
Campbell
EE
,
Haddix
ML
,
Wall
DH
,
Parton
WJ
(
2015
)
Formation of soil organic matter via biochemical and physical pathways of litter mass loss
.
Nat Geosci
 
8
:
776
780
19.
Cotrufo
MF
,
Wallenstein
MD
,
Boot
CM
,
Denef
K
,
Paul
E
(
2013
)
The Microbial Efficiency-Matrix Stabilization (MEMS) framework integrates plant litter decomposition with soil organic matter stabilization: do labile plant inputs form stable soil organic matter?
Glob Change Biol
 
19
:
988
995
20.
Crits-Christoph
A
,
Diamond
S
,
Butterfield
CN
,
Thomas
BC
,
Banfield
JF
(
2018
)
Novel soil bacteria possess diverse genes for secondary metabolite biosynthesis
.
Nature
 
558
:
440
444
21.
Davidson
EA
,
Agren
G
,
Daniel
O
,
Emeis
K
,
Largeau
C
,
Lee
C
,
Mopper
K
,
Oades
JM
,
Reeburgh
WS
,
Schimel
DS
,
Zepp
RG
(
1995
)
What are the physical, chemical and biological processes that control the formation and degradation of nonliving organic matter?
In:
The Role of Nonliving Organic Matter in the Earth's Carbon Cycle
 .
John Wiley & Sons Ltd
22.
Davis
JA
,
Kent
DB
(
1990
)
Surface complexation mode+ling in aqueous geochemistry
.
Rev Mineral
 
23
:
176
260
23.
DeAngelis
KM
,
Chivian
D
,
Fortney
JL
,
Arkin
AP
,
Simmons
B
,
Hazen
TC
,
Silver
WL
(
2013
)
Changes in microbial dynamics during long-term decomposition in tropical forests
.
Soil Biol Biochem
 
66
:
60
68
24.
Dudal
Y
,
Gérard
F
(
2004
)
Accounting for natural organic matter in aqueous chemical equilibrium models: A review of the theories and applications
.
Earth Sci Rev
 
66
:
199
216
25.
Dwivedi
D
,
Mohanty
BP
,
Lesikar
BJ
(
2016
)
Impact of the linked surface water-soil water-groundwater system on transport of E. coli in the Subsurface
.
Water Air Soil Pollut
 
227
:
351
26.
Dwivedi
D
,
Riley
WJ
,
Torn
MS
,
Spycher
N
,
Maggi
F
,
Tang
JY
(
2017a
)
Mineral properties, microbes, transport, and plant-input profiles control vertical distribution and age of soil carbon stocks
.
Soil Biol Biochem
 
107
:
244
259
27.
Dwivedi
D
,
Steefel
CI
,
Arora
B
,
Bisht
G
(
2017b
)
Impact of intra-meander hyporheic flow on nitrogen cycling
.
Procedia Earth Planet Sci
 
17
:
404
407
28.
Dwivedi
D
,
Arora
B
,
Steefel
CI
,
Dafflon
B
,
Versteeg
R
(
2018a
)
Hot spots and hot moments of nitrogen in a riparian corridor
.
Water Resour Res
 
54
:
205
222
29.
Dwivedi
D
,
Steefel
CI
,
Arora
B
,
Newcomer
M
,
Moulton
JD
,
Dafflon
B
,
Faybishenko
B
,
Fox
P
,
Nico
P
,
Spycher
N
,
Carroll
R
,
Williams
KH
(
2018b
)
Geochemical exports to river from the intrameander hyporheic zone under transient hydrologic conditions: East River mountainous watershed, Colorado
.
Water Resour Res
 
4
:
205
222
30.
Dzombak
DA
(
1990
)
Surface Complexation Modeling: Hydrous Ferric Oxide
.
John Wiley & Sons
 
31.
Dzombak
DA
,
Morel
FMM
(
1996
)
Surface Complexation Modeling: Hydrous Ferric Oxide
.
John Wiley & Sons
 
32.
Fontaine
S
,
Barot
S
(
2005
)
Size and functional diversity of microbe populations control plant persistence and long-term soil carbon accumulation
.
Ecol Lett
 
8
:
1075
1087
33.
Fontaine
S
,
Mariotti
A
,
Abbadie
L
(
2003
)
The priming effect of organic matter: A question of microbial competition?
Soil Biol Biochem
 
35
:
837
843
34.
Frostegard
A
,
Baath
E
(
1996
)
The use of phospholipid fatty acid analysis to estimate bacterial and fungal biomass in soil
.
Biol Fertil Soils
 
22
:
59
65
35.
Georgiou
K
,
Abramoff
RZ
,
Harte
J
,
Riley
WJ
,
Torn
MS
(
2017
)
Microbial community-level regulation explains soil carbon responses to long-term litter manipulations
.
Nat Commun
 
8
:
1223
36.
Georgiou
K
,
Harte
J
,
Mesbah
A
,
Riley
WJ
(
2018
)
A method of alternating characteristics with application to advection-dominated environmental systems
.
Comput Geosci 2018
 
22
:
851
865
37.
Goldfarb
KC
,
Karaoz
U
,
Hanson
CA
,
Santee
CA
,
Bradford
MA
,
Treseder
KK
,
Wallenstein
MD
,
Brodie
EL
(
2011
)
Differential growth responses of soil bacterial taxa to carbon substrates of varying chemical recalcitrance
.
Front Microbiol
 
2
:
94
38.
Gordon
AS
,
Millero
FJ
(
1985
)
Adsorption mediated decrease in the biodegradation rate of organic compounds
.
Microbiol Ecol
 
11
:
289
39.
Grant
RF
(
2013
)
Modelling changes in nitrogen cycling to sustain increases in forest productivity under elevated atmospheric CO2 and contrasting site conditions
.
Biogeosciences
 
10
:
7703
7721
40.
Grant
RF
,
Heaney
DJ
(
2010
)
Inorganic phosphorus transformation and transport in soils: mathematical modeling in Ecosys
.
Soil Sci Soc Am J
 
61
:
752
764
41.
Grant
RF
,
Juma
NG
,
McGill
WB
(
1993
)
Simulation of carbon and nitrogen transformations in soil: Mineralization
.
Soil Biol Biochem
 
42.
Grant
RF
,
Oechel
WC
,
Ping
CL
(
2003
)
Modelling carbon balances of coastal arctic tundra under changing climate
.
Glob Change Biol
 
9
:
16
36
43.
Grant
RF
,
Black
TA
,
Gaumont-Guay
D
,
Klujn
N
,
Barr
AG
,
Morgenstern
K
,
Nesic
Z
(
2006
)
Net ecosystem productivity of boreal aspen forests under drought and climate change: Mathematical modelling with Ecosys
.
Agric For Meteorol
 
140
:
152
170
44.
Grant
RF
,
Humphreys
ER
,
Lafleur
PM
,
Dimitrov
DD
(
2011
)
Ecological controls on net ecosystem productivity of a mesic arctic tundra under current and future climates
.
J Geophys Res Biogeosci
 
116
:
G01031
45.
Grant
RF
,
Mekonnen
ZA
,
Riley
WJ
,
Arora
B
,
Torn
MS
(
2017a
)
Mathematical modelling of arctic polygonal tundra with ecosys: 2. microtopography determines How CO2 and CH4 exchange responds to changes in temperature and precipitation
.
J Geophys Res Biogeosci
 
122
:
3174
3187
46.
Grant
RF
,
Mekonnen
ZA
,
Riley
WJ
,
Wainwright
HM
,
Graham
D
,
Torn
MS
(
2017b
)
Mathematical modelling of arctic polygonal tundra with ecosys: 1. Microtopography determines how active layer depths respond to changes in temperature and precipitation
.
J Geophys Res Biogeosci
 
122
:
3161
3173
47.
Gu
BH
,
Schmitt
J
,
Chen
ZH
,
Liang
LY
,
McCarthy
JF
(
1994
)
Adsorption and desorption of natural organic matter on iron oxide: mechanisms and models
.
Environ Sci Technol
 
8
:
38
46
48.
Gu
CH
,
Maggi
F
,
Riley
WJ
,
Hornberger
GM
,
Xu
T
,
Oldenburg
CM
,
Spycher
N
,
Miller
NL
,
Venterea
RT
,
Steefel
C
(
2009
)
Aqueous and gaseous nitrogen losses induced by fertilizer application
,
J Geophys Res-Biogeo
 
114
:
G01006
49.
Gunina
A
,
Kuzyakov
Y
(
2015
)
Sugars in soil and sweets for microorganisms: Review of origin, content, composition and fate
.
Soil Biol Biochem
 
90
:
87
100
50.
Hammond
GE
,
Lichtner
PC
,
Mills
RT
(
2014
)
Evaluating the performance of parallel subsurface simulators: An illustrative example with PFLOTRAN
.
Water Resour Res
 
50
:
208
228
51.
Hedges
JI
,
Eglinton
G
,
Hatcher
PG
,
Kirchman
DL
,
Arnosti
C
,
Derenne
S
,
Evershed
RP
,
Kögel-Knabner
I
,
De Leeuw
JW
,
Littke
R
,
Michaelis
W
,
Rullkötter
J
(
2000
)
The molecularly-uncharacterized component of nonliving organic matter in natural environments
.
Org Geochem
 
31
:
945
958
52.
Hicks Pries
CE
,
Sulman
BN
,
West
C
,
O'Neill
C
,
Poppleton
E
,
Porras
RC
,
Castanha
C
,
Zhu
B
,
Wiedemeier
DB
,
Torn
MS
(
2018
)
Root litter decomposition slows with soil depth
.
Soil Biol Biochem
 
25
:
103
114
53.
Hocking
PJ
(
2004
)
Organic acids exuded from roots in phosphorus uptake and aluminum tolerance of plants in acid soils
.
Adv Agron
 
74
:
63
97
54.
Hunt
AG
,
Manzoni
S
(
2016
)
Networks on Networks The Physics of Geobiology and Geochemistry
.
Morgan & Claypool Publishers
 
55.
Hunt
HW
,
Coleman
DC
,
Ingham
ER
,
Ingham
RE
,
Elliott
ET
,
Moore
JC
,
Rose
SL
,
Reid
CPP
,
Morley
CR
(
1987
)
The detrital food web in a shortgrass prairie
.
Biol Fertil Soils
 
3
:
57
56.
Jagadamma
S
,
Mayes
MA
,
Phillips
JR
(
2012
)
Selective sorption of dissolved organic carbon compounds by temperate soils
.
PLoS One
 
7
:
e50434
57.
Jagadamma
S
,
Mayes
MA
,
Zinn
YL
,
Gísladóttir
G
,
Russell
AE
(
2014
)
Sorption of organic carbon compounds to the fine fraction of surface and subsurface soils
.
Geoderma
 
213
:
79
86
58.
Jenkinson
DS
,
Coleman
K
(
2008
)
The turnover of organic carbon in subsoils. Part 2. Modelling carbon turnover, The turnover of organic carbon in subsoils. Part 2. Modelling carbon turnover
.
Eur J Soil Sci
 
59
:
400
413
59.
Jenkinson
DS
,
Poulton
PR
,
Bryant
C
(
2008
)
The turnover of organic carbon in subsoils. Part 1. Natural and bomb radiocarbon in soil profiles from the Rothamsted long-term field experiments
.
Eur J Soil Sci
 
59
:
391
399
60.
Jilling
A
,
Keiluweit
M
,
Contosta
AR
,
Frey
S
,
Schimel
J
,
Schnecker
J
,
Smith
RG
,
Tiemann
L
,
Grandy
AS
(
2018
)
Minerals in the rhizosphere: overlooked mediators of soil nitrogen availability to plants and microbes
.
Biogeochemistry
 
139
:
103
61.
Jobbágy
EG
,
Jackson
RB
(
2000
)
The vertical distribution of soil organic carbon and its relation to climate and vegetation
.
Ecol Appl
 
10
:
423
436
62.
Joergensen
RG
,
Wichern
F
(
2018
)
Alive and kicking: Why dormant soil microorganisms matter
.
Soil Biol Biochem
 
116
:
419
430
63.
Kaiser
C
,
Franklin
O
,
Dieckmann
U
,
Richter
A
(
2014
)
Microbial community dynamics alleviate stoichiometric constraints during litter decay
.
Ecol Lett
 
17
:
680
690
64.
Keiluweit
M
,
Bougoure
JJ
,
Nico
PS
,
Pett-Ridge
J
,
Weber
PK
,
Kleber
M
(
2015
)
Mineral protection of soil carbon counteracted by root exudates
.
Nat Clim Change
 
5
:
588
595
65.
Kleber
M
,
Nico
PS
,
Plante
A
,
Filley
T
,
Kramer
M
,
Swanston
C
,
Sollins
P
(
2011
)
Old and stable soil organic matter is not necessarily chemically recalcitrant: Implications for modeling concepts and temperature sensitivity
.
Glob Change Biol
 
17
:
1097
1107
66.
Kleber
M
,
Eusterhues
K
,
Keiluweit
M
,
Mikutta
C
,
Mikutta
R
,
Nico
PS
(
2015
)
Mineral–organic associations: formation, properties, and relevance in soil environments
.
Adv Agron
 
130
:
1
140
67.
Knicker
H
(
2011
)
Soil organic N—An under-rated player for C sequestration in soils?
Soil Biol Biochem
 
3
:
1118
1129
68.
Kögel-Knabner
I
(
2017
)
The macromolecular organic composition of plant and microbial residues as inputs to soil organic matter: Fourteen years on
.
Soil Biol Biochem
 
105
:
A3
A8
69.
Kuzyakov
Y
,
Friedel
JK
,
Stahr
K
(
2000
)
Review of mechanisms and quantification of priming effects
.
Soil Biol Biochem
 
32
:
1485
1498
70.
Lawrence
CR
,
Neff
JC
,
Schimel
JP
(
2009
)
Does adding microbial mechanisms of decomposition improve soil organic matter models? A comparison of four models using data from a pulsed rewetting experiment
.
Soil Biol Biochem
 
41
:
1923
1934
71.
Lehmann
J
,
Kleber
M
(
2015
)
The contentious nature of soil organic matter
.
Nature
 
528
:
60
72.
Lehmann
J
,
Kinyangi
J
,
Solomon
D
(
2007
)
Organic matter stabilization in soil microaggregates: Implications from spatial heterogeneity of organic carbon contents and carbon forms
.
Biogeochemistry
 
85
:
45
57
73.
Lehmann
J
,
Solomon
D
,
Kinyangi
J
,
Dathe
L
,
Wirick
S
,
Jacobsen
C
(
2008
)
Spatial complexity of soil organic matter forms at nanometre scales
.
Nat Geosci
 
1
:
238
242
74.
Levin
SA
(
1999
)
Fragile Dominion: Complexity and the Commons
.
Basic Books. Perseus Books
75.
Li
C
,
Cui
J
,
Sun
G
,
Trettin
C
(
2004
)
Modeling impacts of management on carbon sequestration and trace gas emissions in forested wetland ecosystems
.
Environ Manage
 
33
(
Suppl 1
):
S176
S186
76.
Liang
C
,
Schimel
JP
,
Jastrow
JD
(
2017
)
The importance of anabolism in microbial control over soil carbon storage
.
Nat. Microbiol
 
2
:
17105
77.
Luo
Z
,
Wang
E
,
Zheng
H
,
Baldock
JA
,
Sun
OJ
,
Shao
Q
(
2015
)
Convergent modelling of past soil organic carbon stocks but divergent projections
.
Biogeosciences
 
12
:
4373
4383
78.
Luo
Z
,
Wang
E
,
Sun
OJ
(
2017
)
Uncertain future soil carbon dynamics under global change predicted by models constrained by total carbon measurements
.
Ecol Appl
 
27
:
1001
1009
79.
MacArthur
RH
,
Wilson
EO
(
1967
)
The Theory of Island Biogeography
.
Princeton Univ. Press
,
Princeton, NJ
80.
Maggi
F
,
Gu
C
,
Riley
WJ
,
Hornberger
GM
,
Venterea
RT
,
Xu
T
,
Spycher
N
,
Steefel
C
,
Miller
NL
,
Oldenburg
CM
(
2008
)
A mechanistic treatment of the dominant soil nitrogen cycling processes: Model development, testing, and application
.
J Geophys Res
 
113
:
G02016
81.
Maharjan
R
,
Nilsson
S
,
Sung
J
,
Haynes
K
,
Beardmore
RE
,
Hurst
LD
,
Ferenci
T
,
Gudelj
I
(
2013
)
The form of a tradeoff determines the response to competition
.
Ecol Lett
 
16
:
1267
1276
82.
Mairhofer
S
,
Zappala
S
,
Tracy
SR
,
Sturrock
C
,
Bennett
M
,
Mooney
SJ
,
Pridmore
T
(
2012
)
RooTrak: automated recovery of three-dimensional plant root architecture in soil from X-ray microcomputed tomography images using visual tracking
.
Plant Physiol
 
158
:
561
569
83.
Manzoni
S
,
Capek
P
,
Porada
P
,
Thurner
M
,
Winterdahl
M
,
Beer
C
,
Bruchert
V
,
Frouz
J
,
Herrmann
AM
,
Lindahl
BD
,
Lyon
SW
(
2018
)
Reviews and syntheses: Carbon use efficiency from organisms to ecosystems—Definitions, theories, and empirical evidence
.
Biogeosciences
 
15
:
5929
5949
84.
Mathieu
JA
,
Hatté
C
,
Balesdent
J
,
Parent
É
(
2015
)
Deep soil carbon dynamics are driven more by soil type than by climate: A worldwide meta-analysis of radiocarbon profiles
.
Glob Change Biol
 
1
:
4278
4292
85.
Mayes
MA
,
Heal
KR
,
Brandt
CC
,
Phillips
JR
,
Jardine
PM
(
2012
)
Relation between soil order and sorption of dissolved organic carbon in temperate subsoils
.
Soil Sci Soc Am J
 
76
:
1027
86.
Michaelis
L
,
Menten
ML
(
1913
)
Die Kinetik der Invertinwirkung
.
Biochem Z
 
49
:
333
369
87.
Mikutta
R
,
Mikutta
C
,
Kalbitz
K
,
Scheel
T
,
Kaiser
K
,
Jahn
R
(
2007
)
Biodegradation of forest floor organic matter bound to minerals via different binding mechanisms
.
Geochim Cosmochim Acta
 
71
:
2569
2590
88.
Miltner
A
,
Bombach
P
,
Schmidt-Brücken
B
,
Kästner
M
(
2012
)
SOM genesis: Microbial biomass as a significant source
.
Biogeochemistry
 
111
:
41
89.
Mualem
Y
(
1976
)
A new model for predicting the hydraulic conductivity of unsaturated porous media
.
Water Resour Res
 
12
:
513
522
90.
Neely
CL
,
Beare
MH
,
Hargrove
WL
,
Coleman
DC
(
1991
)
Relationships between fungal and bacterial substrate-induced respiration, biomass and plant residue decomposition
.
Soil Biol. Biochem.
 
23
:
947
954
91.
Neuman
SP
(
1973
)
Saturated-unsaturated seepage
.
J Hydraul Div Am Soc Civ Eng
 
99
(
HY12
):
2233
2250
92.
Newcomb
CJ
,
Qafoku
NP
,
Grate
JW
,
Bailey
VL
,
De Yoreo
JJ
(
2017
)
Developing a molecular picture of soil organic matter-mineral interactions by quantifying organo–mineral binding
.
Nat Commun
 
8
:
396
93.
Norde
W
,
Giacomelli
CE
(
2000
)
BSA structural changes during homomolecular exchange between the adsorbed and the dissolved states
.
J Biotechnol
 
79
:
259
268
94.
Omoike
A
,
Chorover
J
(
2004
)
Spectroscopic study of extracellular polymeric substances from Bacillus subtilis: Aqueous chemistry and adsorption effects
.
Biomacromolecules
 
5
:
1219
1230
95.
Or
D
,
Smets
BF
,
Wraith
JM
,
Dechesne
A
,
Friedman
SP
(
2007
)
Physical constraints affecting bacterial habitats and activity in unsaturated porous media—a review
.
Adv Water Resour
 
30
:
1505
1527
96.
Panday
S
,
Huyakorn
PS
,
Therrien
R
,
Nichols
RL
(
1993
)
Improved three-dimensional finite-element techniques for field simulation of variably saturated flow and transport
.
J Contam Hydrol
 
12
:
3
33
97.
Parton
WJ
,
Schimel
DS
,
Cole
CV
,
Ojima
DS
(
1987
)
Analysis of factors controlling soil organic matter levels in great plains grasslands 1
.
Soil Sci Soc Am J
 
51
:
1173
98.
Parton
WJ
,
Hartman
M
,
Ojima
D
,
Schimel
D
(
1998
)
DAYCENT and its land surface submodel: Description and testing
.
Glob Planet Change
 
19
:
35
48
99.
Parton
WJ
,
Hanson
PJ
,
Swanston
C
,
Torn
M
,
Trumbore
SE
,
Riley
W
,
Kelly
R
(
2010
)
ForCent model development and testing using the enriched background isotope study experiment
.
J Geophys Res Biogeosci
 
115
:
G04001
100.
Quiquampoix
H
,
Ratcliffe
RG
(
1992
)
A 31P NMR study of the adsorption of bovine serum albumin on montmorillonite using phosphate and the paramagnetic cation Mn2+: modification of conformation with pH
.
J Colloid Interface Sci
 
148
:
343
352
101.
Razinkov
IA
,
Baumgartner
BL
,
Bennett
MR
,
Tsimring
LS
,
Hasty
J
(
2013
)
Measuring competitive fitness in dynamic environments
.
J Phys Chem B
 
117
:
13175
13181
102.
Richards
LA
(
1931
)
Capillary conduction of liquids through porous mediums
.
J Appl Phys
 
1
:
318
333
103.
Riley
WJ
,
Subin
ZM
,
Lawrence
DM
,
Swenson
SC
,
Torn
MS
,
Meng
L
,
Mahowald
NM
,
Hess
P
(
2011
)
Barriers to predicting changes in global terrestrial methane fluxes : analyses using CLM4Me, a methane biogeochemistry model integrated in CESM
.
Biogeosciences
 
8
:
1925
1953
104.
Riley
WJ
,
Maggi
F
,
Kleber
M
,
Torn
MS
,
Tang
JY
,
Dwivedi
D
,
Guerry
N
(
2014
)
Long residence times of rapidly decomposable soil organic matter: application of a multi-phase, multi-component, and vertically resolved model (BAMS1) to soil carbon dynamics
.
Geosci Model Dev
 
7
:
1335
1355
105.
Riley
WJ
,
Sierra
C
,
Tang
JY
,
Bouskill
NJ
,
Zhu
Q
,
Abramoff
R
(
2019
)
Next generation soil biogeochemistry model representations: A proposed community open source model farm (BeTR-S)
. In:
Multi-scale Biogeochemical Processes in Soil Ecosystems: Critical Reactions and Resilience to Climate Changes
 ,
Y.
Yang
,
M.
Keiluweit
,
N.
Senesi
,
B.
Xing
(eds)
In Press Dwivedi et al. Abiotic and Biotic Controls on Soil Organo–Mineral Interactions
106.
Romaní
AM
,
Fischer
H
,
Mille-Lindblom
C
,
Tranvik
LJ
(
2006
)
Interactions of bacteria and fungi on decomposing litter: Differential extracellular enzyme activities
.
Ecology
 
87
:
2559
2569
107.
Schimel
JP
,
Schaeffer
SM
(
2012
)
Microbial control over carbon cycling in soil
.
Front Microbiol
 
3
:
348
108.
Schimel
JP
,
Wetterstedt
JÅM
,
Holden
PA
,
Trumbore
SE
(
2011
)
Drying/rewetting cycles mobilize old C from deep soils from a California annual grassland
.
Soil Biol Biochem
 
43
:
1101
1103
109.
Schmidt
MWI
,
Torn
MS
,
Abiven
S
,
Dittmar
T
,
Guggenberger
G
,
Janssens
IA
,
Kleber
M
,
Kögel-Knabner
I
,
Lehmann
J
,
Manning
DAC
,
Nannipieri
P
,
Rasse
DP
,
Weiner
S
,
Trumbore
SE
(
2011
)
Persistence of soil organic matter as an ecosystem property
.
Nature
 
478
:
49
56
110.
Servagent-Noinville
S
,
Revault
M
,
Quiquampoix
H
,
Baron
MH
(
2000
)
Conformational changes of bovine serum albumin induced by adsorption on different clay surfaces: FTIR analysis
.
J Colloid Interface Sci
 
221
:
273
283
111.
Simpson
AJ
,
Kingery
WL
,
Shaw
DR
,
Spraul
M
,
Humpfer
E
,
Dvortsak
P
(
2001
)
The application of 1H HR-MAS NMR spectroscopy for the study of structures and associations of organic components at the solid—Aqueous interface of a whole soil
.
Environ Sci Technol
 
35
:
3321
3332
112.
Smith
OL
(
1979
)
An analytical model of the decomposition
.
J Biomech
 
9
:
397
405
113.
Steefel
CI
,
Lasaga
AC
(
1994
)
A coupled model for transport of multiple chemical species and kinetic precipitation/dissolution reactions with application to reactive flow in single phase hydrothermal systems
.
Am J Sci
 
294
:
529
592
114.
Steefel
CI
,
Appelo
CAJJ
,
Arora
B
,
Jacques
D
,
Kalbacher
T
,
Kolditz
O
,
Lagneau
V
,
Lichtner
PC
,
Mayer
KU
,
Meeussen
JCLL
,
Molins
S
,
Moulton
D
,
Shao
H
,
Šimůnek
J
,
Spycher
N
,
Yabusaki
SB
,
Yeh
GT
(
2015
)
Reactive transport codes for subsurface environmental simulation
.
Comput Geosci
 
19
:
445
478
115.
Strickland
MS
,
Rousk
J
(
2010
)
Considering fungal: Bacterial dominance in soils—Methods, controls, and ecosystem implications
.
Soil Biol Biochem
 
42
:
1385
1395
116.
Sulman
BN
,
Phillips
RP
,
Oishi
AC
,
Shevliakova
E
,
Pacala
SW
(
2014
)
Microbe-driven turnover offsets mineral-mediated storage of soil carbon under elevated CO2
.
Nat Clim Change
 
4
:
1099
1102
117.
Sulman
BN
,
Moore
JAM
,
Abramoff
R
,
Averill
C
,
Kivlin
S
,
Georgiou
K
,
Sridhar
B
,
Hartman
MD
,
Wang
G
,
Wieder
WR
,
Bradford
MA
,
Luo
Y
,
Mayes
MA
,
Morrison
E
,
Riley
WJ
,
Salazar
A
,
Schimel
JP
,
Tang
J
,
Classen
AT
(
2018
)
Multiple models and experiments underscore large uncertainty in soil carbon dynamics
.
Biogeochemistry
 
141
:
109
118.
Swenson
TL
,
Bowen
BP
,
Nico
PS
,
Northen
TR
(
2015
)
Competitive sorption of microbial metabolites on an iron oxide mineral
.
Soil Biol Biochem
 
90
:
34
41
119.
Taina
IA
,
Heck
RJ
,
Elliot
TR
(
2011
)
Application of X-ray computed tomography to soil science: A literature review
.
Can J Soil Sci
 
88
:
1
19
120.
Tang
JY
(
2015
)
On the relationships between the Michaelis-Menten kinetics, reverse Michaelis-Menten kinetics, equilibrium chemistry approximation kinetics, and quadratic kinetics
.
Geosci Model Dev
 
8
:
3823
3835
121.
Tang
JY
,
Riley
WJ
(
2013
)
A total quasi-steady-state formulation of substrate uptake kinetics in complex networks and an example application to microbial litter decomposition
.
Biogeosciences
 
10
:
8329
8351
122.
Tang
J
,
Riley
WJ
(
2014
)
Weaker carbon–climate feedbacks resulting from microbial and abiotic interactions
.
Nature Clim Change
 :
2438
123.
Tang
JY
,
Riley
WJ
(
2017
)
SUPECA kinetics for scaling redox reactions in networks of mixed substrates and consumers and an example application to aerobic soil respiration
.
Geosci Model Dev
 
10
:
3277
3295
124.
Tang
JY
,
Riley
WJ
,
Koven
CD
,
Subin
ZM
(
2013
)
CLM4-BeTR, a generic biogeochemical transport and reaction module for CLM4: model development, evaluation, and application
.
Geosci Model Dev
 
6
:
127
140
125.
Tang
G
,
Zheng
J
,
Xu
X
,
Yang
Z
,
Graham
DE
,
Gu
B
,
Painter
SL
,
Thornton
PE
(
2016
)
Biogeochemical modeling of CO2 and CH4 production in anoxic Arctic soil microcosms
.
Biogeosciences
 
13
:
5021
5041
126.
Tang
J
,
Riley
WJ
(
2018
)
Predicted land carbon dynamics are strongly dependent on the numerical coupling of nitrogen mobilizing and immobilizing processes : A Demonstration with the E3SM Land Model
.
Earth Interact
 
22
:
1
8
127.
Tang
FHM
,
Riley
WJ
,
Maggi
F
(
2019
),
Hourly and daily rainfall intensification causes opposing effects on C and N emissions, storage, and leaching in dry and wet grasslands
.
Biogeochemistry
 .
In Press
.
128.
Tfaily
MM
,
Chu
RK
,
Tolić
N
,
Roscioli
KM
,
Anderton
CR
,
Paša-Tolić
L
,
Robinson
EW
,
Hess
NJ
(
2015
)
Advanced solvent based methods for molecular characterization of soil organic matter by high-resolution mass spectrometry
.
Anal Chem
 
87
:
5206
5215
129.
Thevenot
M
,
Dignac
MF
,
Rumpel
C
(
2010
)
Fate of lignins in soils: A review
.
Soil Biol Biochem
 
42
:
1200
1211
130.
Tipping
E
(
1994
)
WHAMC-A chemical equilibrium model and computer code for waters, sediments, and soils incorporating a discrete site/electrostatic model of ion-binding by humic substances
.
Comput Geosci
 
20
:
973
1023
131.
Trumbore
SE
,
Davidson
EA
,
Barbosa De Camargo
P
,
Nepstad
DC
,
Martinelli
LA
(
1995
)
Belowground cycling of carbon in forests and pastures of eastern Amazonia
.
Global Biogeochem Cycles
 
9
:
515
528
132.
van Genuchten
MT
(
1980
)
A closed-form equation for predicting the hydraulic conductivity of unsaturated soils 1
.
Soil Sci Soc Am J
 
44
:
892
133.
Vargas
R
,
Detto
M
,
Baldocchi
DD
,
Allen
MF
(
2010
)
Multiscale analysis of temporal variability of soil CO2 production as influenced by weather and vegetation
.
Glob Change Biol
 
16
:
1589
1605
134.
Viollea
C
,
Reich
PB
,
Pacala
SW
,
Enquist
BJ
,
Kattge
J
(
2014
)
The emergence and promise of functional biogeography
.
PNAS
 
111
:
13690
13696
135.
Wang
G
,
Post
WM
,
Mayes
MA
(
2013
)
Development of microbial-enzyme-mediated decomposition model parameters through steady-state and dynamic analyses
.
Ecol Appl
 
3
:
255
272
136.
Wieder
WR
,
Grandy
AS
,
Kallenbach
CM
,
Bonan
GB
(
2014
)
Integrating microbial physiology and physio-chemical principles in soils with the MIcrobial-MIneral Carbon Stabilization (MIMICS) model
.
Biogeosciences
 
11
:
3899
3917
137.
Wu
L
,
Kobayashi
Y
,
Wasaki
J
,
Koyama
H
(
2018
)
Organic acid excretion from roots: a plant mechanism for enhancing phosphorus acquisition, enhancing aluminum tolerance, and recruiting beneficial rhizobacteria
.
Soil Sci Plant Nutr
 
64
:
697
704
138.
Xu
T
,
Sonnenthal
E
,
Spycher
N
,
Pruess
K
(
2006
)
TOUGHREACT—A simulation program for non-isothermal multiphase reactive geochemical transport in variably saturated geologic media: Applications to geothermal injectivity and CO2 geological sequestration
.
Comput Geosci
 
32
:
145
165
139.
Yabusaki
SB
,
Wilkins
MJ
,
Yilin
F
,
Williams
KH
,
Arora
B
,
Bargar
JR
,
Beller
HR
,
Bouskill
NJ
,
Brodie
EL
,
Christensen
JN
,
Conrad
ME
,
Danczak
RE
,
King
E
,
Spycher
NF
,
Steefel
CI
,
Tokunaga
TK
,
Versteeg
RJ
,
Waichler
SR
,
Wainwright
HM
(
2017
)
Water table dynamics and biogeochemical cycling in a shallow, variably-saturated floodplain
.
Environ Sci Technol
 
51
:
3307
3317
140.
Yan
J
,
Pan
G
,
Li
L
,
Quan
G
,
Ding
C
,
Luo
A
(
2010
)
Adsorption, immobilization, and activity of β-glucosidase on different soil colloids
.
J Colloid Interface Sci
 
348
:
565
570
141.
Zhang
Y
,
Li
C
,
Trettin
CC
,
Li
H
,
Sun
G
(
2002
)
An integrated model of soil, hydrology, and vegetation for carbon dynamics in wetland ecosystems
.
Global Biogeochem Cycles
 
16
:
1
17
142.
Zheng
J
,
Thornton
PE
,
Painter
SL
,
Gu
B
,
Wullschleger
SD
, and
Graham
DE
(
2019
)
Modeling anaerobic soil organic carbon decomposition in Arctic polygon tundra: Insights into soil geochemical influences on carbon mineralization
.
Biogeosciences
 
16
:
663
680
143.
Zhu
Q
,
Iversen
CM
,
Riley
WJ
,
Slette
IJ
(
2016a
)
A new theory of plant–microbe nutrient competition resolves inconsistencies between observations and model predictions
.
Ecol Appl
 
27
:
875
886
144.
Zhu
Q
,
Riley
WJ
,
Tang
J
,
Koven
CD
(
2016b
)
Multiple soil nutrient competition between plants, microbes, and mineral surfaces: model development, parameterization, and example applications in several tropical forests
.
Biogeosci ences
 
13
:
341
363
145.
Zhu
Q
,
Riley
WJ
,
Tang
JY
(
2017
),
A new theory of plant and microbe nutrient competition resolves inconsistencies between observations and models
.
Ecol Appl
 
27
:
875
888
146.
Zhuang
Q
,
Melillo
JM
,
Kicklighter
DW
,
Prinn
RG
,
McGuire
AD
,
Steudler
PA
,
Felzer
BS
,
Hu
S
(
2004
)
Methane fluxes between terrestrial ecosystems and the atmosphere at northern high latitudes during the past century: A retrospective analysis with a process-based biogeochemistry model
.
Global Biogeochem Cycles
 
16
:
663
680
https://github.com/BeTR-biogeochemistry-modeling/sbetr; accessed on 06 May 2019.
Open access: Article available to all readers online.