Abstract
Global ecosystem models may require microbial components to accurately predict feedbacks between climate warming and soil decomposition, but it is unclear what parameters and levels of complexity are ideal for scaling up to the globe. Here we conducted a model comparison using a conventional model with first-order decay and three microbial models of increasing complexity that simulate short- to long-term soil carbon dynamics. We focused on soil carbon responses to microbial carbon use efficiency (CUE) and temperature. Three scenarios were implemented in all models: constant CUE (held at 0.31), varied CUE (−0.016 °C−1), and 50 % acclimated CUE (−0.008 °C−1). Whereas the conventional model always showed soil carbon losses with increasing temperature, the microbial models each predicted a temperature threshold above which warming led to soil carbon gain. The location of this threshold depended on CUE scenario, with higher temperature thresholds under the acclimated and constant scenarios. This result suggests that the temperature sensitivity of CUE and the structure of the soil carbon model together regulate the long-term soil carbon response to warming. Equilibrium soil carbon stocks predicted by the microbial models were much less sensitive to changing inputs compared to the conventional model. Although many soil carbon dynamics were similar across microbial models, the most complex model showed less pronounced oscillations. Thus, adding model complexity (i.e. including enzyme pools) could improve the mechanistic representation of soil carbon dynamics during the transient phase in certain ecosystems. This study suggests that model structure and CUE parameterization should be carefully evaluated when scaling up microbial models to ecosystems and the globe.
Similar content being viewed by others
Explore related subjects
Discover the latest articles, news and stories from top researchers in related subjects.References
Ågren GI, Bosatta E (2002) Reconciling differences in predictions of temperature response of soil organic matter. Soil Biol Biochem 34:129–132
Ågren GI, Wetterstedt JÅM (2007) What determines the temperature response of soil organic matter decomposition? Soil Biol Biochem 39:1794–1798
Allison SD, Wallenstein MD, Bradford MA (2010) Soil-carbon response to warming dependent on microbial physiology. Nat Geosci 3:336–340
Bolker BM, Pacala SW, Parton WJ (1998) Linear analysis of soil decomposition: insights from the century model. Ecol Appl 8:425–439
Bradford MA, Davies CA, Frey SD, Maddox TR, Melillo JM, Mohan JE, Reynolds JF, Treseder KK, Wallenstein MD (2008) Thermal adaptation of soil microbial respiration to elevated temperature. Ecol Lett 11:1316–1327
del Giorgio PA, Cole JJ (1998) Bacterial growth efficiency in natural aquatic systems. Annu Rev Ecol Syst 29:503–541
Devevre OC, Horwath WR (2000) Decomposition of rice straw and microbial carbon use efficiency under different soil temperatures and moistures. Soil Biol Biochem 32:1773–1785
Fontaine S, Mariotti A, Abbadie L (2003) The priming effect of organic matter: a question of microbial competition? Soil Biol Biochem 35:837–843
Frey SD, Lee J, Melillo JM, Six J (2013) The temperature response of soil microbial efficiency and its feedback to climate. Nat Clim Change 3:395–398
Friedlingstein P, Cox P, Betts R, Bopp L, Von Bloh W, Brovkin V, Cadule P, Doney S, Eby M, Fung I, Bala G, John J, Jones C, Joos F, Kato T, Kawamiya M, Knorr W, Lindsay K, Matthews HD, Raddatz T, Rayner P, Reick C, Roeckner E, Schnitzler KG, Schnur R, Strassmann K, Weaver AJ, Yoshikawa C, Zeng N (2006) Climate-carbon cycle feedback analysis: results from the (CMIP)-M-4 model intercomparison. J Clim 19:3337–3353
German DP, Marcelo KRB, Stone MM, Allison SD (2012) The Michaelis–Menten kinetics of soil extracellular enzymes in response to temperature: a cross-latitudinal study. Glob Change Biol 18:1468–1479
Hartley IP, Heinemeyer A, Ineson P (2007) Effects of three years of soil warming and shading on the rate of soil respiration: substrate availability and not thermal acclimation mediates observed response. Glob Change Biol 13:1761–1770
Hartley IP, Hopkins DW, Garnett MH, Sommerkorn M, Wookey PA (2008) Soil microbial respiration in arctic soil does not acclimate to temperature. Ecol Lett 11:1092–1100
Jobbagy EG, Jackson RB (2000) The vertical distribution of soil organic carbon and its relation to climate and vegetation. Ecol Appl 10:423–436
Knorr W, Prentice IC, House JI, Holland EA (2005) Long-term sensitivity of soil carbon turnover to warming. Nature 433:298–301
Kuzyakov Y (2010) Priming effects: interactions between living and dead organic matter. Soil Biol Biochem 42:1363–1371
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
Li J, Ziegler SE, Lane CS, Billings SA (2013) Legacies of native climate regime govern responses of boreal soil microbes to litter stoichiometry and temperature. Soil Biol Biochem 66:204–213
Luo YQ, Wan SQ, Hui DF, Wallace LL (2001) Acclimatization of soil respiration to warming in a tall grass prairie. Nature 413:622–625
Malcolm GM, Lopez-Gutierrez JC, Koide RT, Eissenstat DM (2008) Acclimation to temperature and temperature sensitivity of metabolism by ectomycorrhizal fungi. Glob Change Biol 14:1169–1180
Manzoni S, Porporato A (2007) A theoretical analysis of nonlinearities and feedbacks in soil carbon and nitrogen cycles. Soil Biol Biochem 39:1542–1556
Manzoni S, Taylor P, Richter A, Porporato A, Agren GI (2012) Environmental and stoichiometric controls on microbial carbon-use efficiency in soils. New Phytol 196:79–91
McGuire AD, Melillo JM, Kicklighter DW, Joyce LA (1995) Equilibrium responses of soil carbon to climate change: empirical and process-based estimates. J Biogeogr 22:785–796
Melillo JM, Steudler PA, Aber JD, Newkirk K, Lux H, Bowles FP, Catricala C, Magill A, Ahrens T, Morrisseau S (2002) Soil warming and carbon-cycle feedbacks to the climate system. Science 298:2173–2176
Melillo JM, Butler S, Johnson J, Mohan J, Steudler P, Lux H, Burrows E, Bowles F, Smith R, Scott L, Vario C, Hill T, Burton A, Zhou YM, Tang J (2011) Soil warming, carbon–nitrogen interactions, and forest carbon budgets. Proc Natl Acad Sci USA 108:9508–9512
Parton WJ, Schimel DS, Cole CV, Ojima DS (1987) Analysis of factors controlling soil organic-matter levels in great-plains grasslands. Soil Sci Soc Am J 51:1173–1179
Parton WJ, Stewart JWB, Cole CV (1988) Dynamics of C, N, P and S in grassland soils—a model. Biogeochemistry 5:109–131
Rustad LE, Campbell JL, Marion GM, Norby RJ, Mitchell MJ, Hartley AE, Cornelissen JHC, Gurevitch J, Gcte-News (2001) A meta-analysis of the response of soil respiration, net nitrogen mineralization, and aboveground plant growth to experimental ecosystem warming. Oecologia 126:543–562
Schimel J (2013) Soil carbon: microbes and global carbon. Nat Clim Change 3:867–868
Schimel JP, Weintraub MN (2003) The implications of exoenzyme activity on microbial carbon and nitrogen limitation in soil: a theoretical model. Soil Biol Biochem 35:549–563
Sinsabaugh RL, Antibus RK, Linkins AE (1991) An enzymatic approach to the analysis of microbial activity during plant litter decomposition. Agric Ecosyst Environ 34:43–54
Sinsabaugh RL, Manzoni S, Moorhead DL, Richter A (2013) Carbon use efficiency of microbial communities: stoichiometry, methodology and modelling. Ecol Lett 16:930–939
Six J, Frey SD, Thiet RK, Batten KM (2006) Bacterial and fungal contributions to carbon sequestration in agroecosystems. Soil Sci Soc Am J 70:555–569
Thornley JHM, Cannell MGR (2001) Soil carbon storage response to temperature: an hypothesis. Ann Bot 87:591–598
Thornton PE, Doney SC, Lindsay K, Moore JK, Mahowald N, Randerson JT, Fung I, Lamarque JF, Feddema JJ, Lee YH (2009) Carbon-nitrogen interactions regulate climate-carbon cycle feedbacks: results from an atmosphere-ocean general circulation model. Biogeosciences 6:2099–2120
Todd-Brown KEO, Hopkins FM, Kivlin SN, Talbot JM, Allison SD (2012) A framework for representing microbial decomposition in coupled climate models. Biogeochemistry 109:19–33
Todd-Brown KEO, Randerson JT, Post WM, Hoffman FM, Tarnocai C, Schuur EAG, Allison SD (2013) Causes of variation in soil carbon simulations from CMIP5 Earth system models and comparison with observations. Biogeosciences 10:1717–1736
Treseder KK, Balser TC, Bradford MA, Brodie EL, Dubinsky EA, Eviner VT, Hofmockel KS, Lennon JT, Levine UY, MacGregor BJ, Pett-Ridge J, Waldrop MP (2012) Integrating microbial ecology into ecosystem models: challenges and priorities. Biogeochemistry 109:7–18
Tucker CL, Bell J, Pendall E, Ogle K (2013) Does declining carbon-use efficiency explain thermal acclimation of soil respiration with warming? Glob Change Biol 19:252–263
Wang G, Post WM (2012) A theoretical reassessment of microbial maintenance and implications for microbial ecology modeling. FEMS Microbiol Ecol 81:610–617
Wang GS, Post WM, Mayes MA (2013a) Development of microbial-enzyme-mediated decomposition model parameters through steady-state and dynamic analyses. Ecol Appl 23:255–272
Wang YP, Chen BC, Wieder WR, Luo YQ, Leite M, Medlyn BE, Rasmussen M, Smith MJ, Agusto FB, Hoffman F (2013b) Oscillatory behavior of two nonlinear microbial models of soil carbon decomposition. Biogeosci Discuss 10:19661–19700
Wieder WR, Bonan GB, Allison SD (2013) Global soil carbon projections are improved by modelling microbial processes. Nat Clim Change 3:909–912
Zelenev VV, van Bruggen AHC, Semenov AM (2005) Short-term wavelike dynamics of bacterial populations in response to nutrient input from fresh plant residues. Microb Ecol 49:83–93
Zhou JZ, Xue K, Xie JP, Deng Y, Wu LY, Cheng XH, Fei SF, Deng SP, He ZL, Van Nostrand JD, Luo YQ (2012) Microbial mediation of carbon-cycle feedbacks to climate warming. Nat Clim Change 2:106–110
Acknowledgments
We thank two anonymous reviewers for their valuable and insightful comments. This research was funded by US National Science Foundation (NSF) grants DBI 0850290, EPS 0919466, DEB 0743778, DEB 0840964, EF 1137293, and EF 0928388 and was also funded in part by the Laboratory Directed Research and Development (LDRD) Program of the Oak Ridge National Laboratory (ORNL) and by the U.S. Department of Energy Biological and Environmental Research program. ORNL is managed by UT-Battelle, LLC, for the U.S. Department of Energy under contract DE-AC05-00OR22725. Part of the model runs were performed at the Supercomputing Center for Education & Research (OSCER), University of Oklahoma. This manuscript has been authored by UT-Battelle, LLC, under Contract No. DE-AC05-00OR22725 with the U.S. Department of Energy. The United States Government retains and the publisher, by accepting the article for publication, acknowledges that the United States Government retains a non-exclusive, paid-up, irrevocable, world-wide license to publish or reproduce the published form of this manuscript, or allow others to do so, for United States Government purposes.
Author information
Authors and Affiliations
Corresponding author
Additional information
Responsible Editor: W. Troy Baisden.
Electronic supplementary material
Below is the link to the electronic supplementary material.
Appendix 1
Appendix 1
Conventional model (CON)
The conventional model is representative of first-order models of soil organic carbon (SOC) dynamics. This model includes SOC, dissolved organic C (DOC), and microbial biomass C (MBC) pools with the decomposition rate of each pool represented as a first-order process. The decay constant k i increases exponentially with temperature T according to the Arrhenius relationship:
where k i,ref is the decay constant at the reference temperature T ref (K), and E ai is the activation energy with i = D, S, or B representing DOC, SOC, and MBC pools, respectively. R is the ideal gas constant, 8.314 J mol−1 K−1. Decomposition of each pool is represented as:
The change in the SOC pool is proportional to external inputs (I S ), transfers from the other pools, and losses due to first-order decomposition:
where a DS is the transfer coefficient from the DOC to the SOC pool, a B is the transfer coefficient from the MBC to the DOC and SOC pools, and a BS is the partition coefficient for dead microbial biomass between the SOC and DOC pools. Transfer coefficients can range from 0.0 to 1.0, with lower values indicating a larger fraction of C respired as CO2. The change in the DOC pool is represented similarly, but includes a transfer from SOC to DOC in proportion to a SD and a loss due to microbial uptake, u * D:
The change in the microbial biomass pool is the difference between uptake and turnover, where u represents the fraction h−1 of the DOC pool taken up by microbial biomass:
The CO2 respiration rate is the sum of the proportion of fluxes that do not enter soil pools:
Steady state analytical solution
The steady state analytical solutions for the DOC, SOC, and MBC pools in CON are:
GER
The GER microbial model represents SOC change as a function of input rate I S , microbial turnover r B , MBC, and extracellular enzyme V max and K m :
C inputs and dead biomass enter the SOC pool, and SOC is lost through decomposition, which is assumed to be a Michaelis–Menten process represented by the last term in Eq. 15. MBC change is a function of microbial turnover and assimilation of decomposed soil organic C, which occurs with C use efficiency E C :
where E C is a linear function of temperature with slope m:
The CO2 respiration rate (C R ) is then the fraction of decomposition not assimilated by microbial biomass:
V max and K m (Y) have an Arrhenius dependence on temperature, similar to Eq. 4 in the conventional model:
Steady state analytical solution
The steady state analytical solutions for the SOC and MBC pools in GER are:
where E C must be larger than r B /V, otherwise microbes cannot assimilate enough C to compensate for microbial turnover; if E C = 1, then microbes respire no C, all C is assimilated, and biomass grows indefinitely.
AWB
AWB is a more complex version of GER that includes explicit DOC and ENZC pools. Microbial biomass increases with DOC uptake (F U ) times C use efficiency and declines with death (F B ) and enzyme production (F E ):
where assimilation is a Michaelis–Menten function scaled to the size of the microbial biomass pool:
Microbial biomass death is modeled as a first-order process with a rate constant r B :
Enzyme production is modeled as a constant fraction (r E ) of microbial biomass:
Temperature sensitivities for V, V U , K, and K U follow the Arrhenius relationship as in Eq. 19. Note that this relationship differs from the published version of AWB that used a linear relationship for K and K U temperature sensitivity. We used the Arrhenius relationship here to facilitate comparison with the other models and used the parameter values from the linear relationship at 20 °C as the reference values in Eq. 19. CO2 respiration is the fraction of DOC uptake that is not assimilated into MBC:
The enzyme pool increases with enzyme production and decreases with enzyme turnover:
where enzyme turnover is modeled as a first-order process with a rate constant r L :
The SOC pool increases with external inputs and a fraction (a BS ) of microbial biomass death and decreases due to decomposition losses:
where decomposition of SOC is catalyzed according to Michaelis–Menten kinetics by the enzyme pool:
The DOC pool receives external inputs, the remaining fraction of microbial biomass death, the decomposition flux, and dead enzymes, while assimilation of DOC by microbial biomass is subtracted:
Steady state analytical solution
The steady state analytical solutions for SOC, DOC, MBC, and ENZC in AWB are:
which simplifies to the following if I D = I S :
MEND
Five C pools are considered in MEND: (1) particulate organic carbon (POC, represented by the variable P in model equations), (2) mineral-associated organic carbon (MOC, M), (3) active layer of MOC (Q) interacting with dissolved organic carbon through adsorption and desorption, (4) dissolved organic carbon (DOC, D), (5) microbial biomass carbon (MBC, B), and (6) extracellular enzymes (EP and EM). The component fluxes are DOC uptake by microbes (denoted by the flux F 1), POC decomposition (F 2), MOC decomposition (F 3), microbial growth respiration (F 4) and maintenance respiration (F 5), adsorption (F 6) and desorption (F 7), microbial mortality (F 8), enzyme production (F 9), and enzyme turnover (F 10). Model equations for each component are listed as follows:
where V i and K i represent the V max and K m for enzymatic degradation of pool i, m R is the maintenance respiration rate, Q max is the maximum DOC sorption capacity, K des and K ads are the specific adsorption and desorption rates, p i is the fraction of m R associated with production of enzyme i, and r i is the turnover rate of enzyme pool i. V i , K i , m R , K des , and K ads follow Arrhenius temperature sensitivity similar to Eq. 19, and E C is linearly dependent on temperature as in Eq. 17. The differential equations are as follows for the pools:
and the CO2 respiration rate is calculated as:
MEND represents microbial respiration as a fraction of assimilation (Eqs. 40, 41) whereas GER and AWB represent respiration as a fraction of microbial uptake (Eqs. 18, 26); note that these representations are algebraically identical with respect to CUE.
Steady state analytical solution
The steady state analytical solutions to the MEND differential equations are as follows:
where
Equations 55–56 simplify to the following if I D ≪ I P :
See Table 2 for all model parameter values.
Rights and permissions
About this article
Cite this article
Li, J., Wang, G., Allison, S.D. et al. Soil carbon sensitivity to temperature and carbon use efficiency compared across microbial-ecosystem models of varying complexity. Biogeochemistry 119, 67–84 (2014). https://doi.org/10.1007/s10533-013-9948-8
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s10533-013-9948-8