A large and poorly understood component of global warming is the terrestrial carbon cycle feedback to the climate system (
1). Simulation experiments with fully coupled, three-dimensional carbon-climate models suggest that carbon cycle feedbacks could substantially accelerate or slow climate change over the 21st century (
2–
4). Both the sign and magnitude of these feedbacks in the real Earth system are still highly uncertain because of gaps in basic understanding of terrestrial ecosystem processes. For example, the potential switch of the terrestrial biosphere from its current role as a carbon sink to a carbon source is critically dependent on the long-term temperature sensitivity of soil organic matter (SOM) decay (
5–
7) and complex carbon-nitrogen interactions that will likely occur in a warmer world (
8–
12). However, without long-term field-based experiments, the sign of the feedback cannot be determined, the complex mechanisms regulating that feedback cannot be quantified, and models that incorporate the soil’s role in carbon feedbacks to the climate system cannot be tested. Here, we present results from a long-term (26-year) soil-warming experiment designed to explore these feedback issues in an ecosystem context.
We started our soil warming study in 1991 in an even-aged mixed hardwood forest stand at the Harvard Forest in central Massachusetts (42.54°N, 72.18°W), where the dominant tree species are red maple (Acer rubrum L.) and black oak (Quercus velutina Lam.). The soil is a stony loam with a distinct organic matter–rich forest floor. (See the supplementary materials for more information on the site’s soils, climate, and land-use history.)
The field manipulation contains 18 plots, each 6 × 6 m, that are grouped into six blocks. The three plots within each block are randomly assigned to one of three treatments: (i) heated plots in which the average soil temperature is continuously elevated 5°C above ambient by the use of buried heating cables; (ii) disturbance control plots that are identical to the heated plots except that they receive no electrical power; and (iii) undisturbed control plots that have been left in their natural state (no cables). The heating method works well under a variety of moisture and temperature conditions (
13). Here, we compared carbon dynamics measured in the heated plots to those measured in the disturbance control plots, so as to isolate heating effects from the effects of cable installation (e.g., root cutting and soil compaction) (
14).
We used a static chamber technique (
14) to measure soil CO
2 emission rates in the study plots monthly between April and November each year for 26 years (see supplementary materials). An ephemeral but significant soil respiration response to warming occurred over the experiment’s first decade (
Fig. 1), with soil respiration greater in the heated plots than in the controls (phase I). In the next phase of the response, phase II (years 11 to 17), soil respiration rates in the heated plots were generally equal to or less than those in the control plots. During years 18 to 23 of the study (phase III), we observed a second trend reversal, with soil respiration once again higher in the heated plots. In the three most recent years (phase IV) of the study, 2014 to 2016, soil respiration rates in the heated plots were again equal to or less than those in the control plots.
We partitioned soil respiration into its two components, root and microbial respiration. On the basis of field measurements of root respiration made between June and November 2009, we developed a temperature-driven root respiration model that we used to estimate root respiration over the course of the experiment (see supplementary materials). We calculated microbial respiration by difference, with microbial respiration equaling soil respiration minus root respiration (
10). Our estimate is that two-thirds of the cumulative CO
2-C emitted from the plots over the 26-year study has been microbial. By this analysis, we calculate a warming-induced soil carbon loss from the full soil profile over the 26-year study of 1510 ± 160 g C m
−2, which is equivalent to a 17% loss of the soil carbon found in the top 60 cm of the soil at the start of the experiment. With respect to timing, we estimate that about three-quarters of this soil carbon loss occurred during phase I, and the remaining quarter during phase III. No measurable carbon loss occurred during either phase II or phase IV (
Fig. 2).
In both the control and heated plots, we made direct measurements of carbon stocks in the upper horizons of the soil profile—the distinct, organic matter–rich surface horizon or forest floor and the top 30 cm of the mineral soil just beneath the forest floor (fig. S1). We measured a carbon loss from the forest floor in response to soil warming of 800 ± 300 g C m
−2, which represents a 31% reduction in forest floor carbon stock over the 26-year study. With our direct measurements, we did not detect any statistically significant changes in the carbon stocks across the top 30 cm of the mineral soil horizon. However, combining our estimate of carbon loss from the full profile based on the respiration measurements with the carbon loss measured from the forest floor, we estimate that the warming-induced carbon loss from the mineral soil in the full soil profile over the study period was 710 g C m
−2. An in situ soil warming experiment in a California forest also shows that warming increases the decay of subsoil organic matter (
15).
We explored possible relationships between the multiphase pattern in soil respiration and multiyear variations in climate, and found none (see supplementary materials and fig. S2, A and B). We propose that several biogeochemical factors combine to affect the timing, magnitude, and thermal acclimation of soil carbon loss. These include depletion of microbially accessible carbon pools, reductions in microbial biomass, a shift in microbial carbon use efficiency, and changes in microbial community composition. Thermal acclimation, a phenomenon observed in response to soil warming in a grassland ecosystem (
16), was evident in all four phases of our long-term soil warming experiment, such that at a given temperature there was less respiration from the heated plots than from the control plots (
Fig. 3). During phases I and III, the acclimation was insufficient to compensate for the 5°C increase in soil temperature, so that CO
2 emissions from heated plots were greater than from control plots. During phases II and IV, however, acclimation was large enough to compensate for the 5°C increase in soil temperature.
In this study, the importance of thermal acclimation of the soil respiration response became evident when power to the heated plots was off for part of the year in 1995 and 2005 and throughout 2010. The power shutdowns in 1995 and 2005 resulted when summer lightning strikes damaged the system that controlled experimental heating. The power shutdown in 2010 was planned and carried out to increase our understanding of the thermal acclimation response. In the three instances, soil respiration rates in the powered-down heated plots dropped below those of the controls after the power was off, and they returned (within weeks) to rates above the controls once the power was restored (fig. S3).
Our biogeochemical and molecular observations suggest that warming causes cycles of soil carbon decay punctuated by periods of structural and functional changes in the microbial community. Sustained reductions of microbial biomass over the course of the experiment [(
17–
19) and this study] have been accompanied by several other changes, including (i) altered respiratory and lipid profiles (
17); (ii) changes in microbial community structure and function as determined using small subunit ribosomal RNA analysis (
20,
21), metagenomics (
22), and enzyme assays (
19,
21); (iii) characterization of substrate utilization profiles of bacterial isolates (
22); and (iv) measurements of microbial carbon use efficiency (
23). We have integrated these observations (
Table 1) to develop a conceptual model of the time-varying (four-phase) effects of soil warming on feedback to climate.
Phase I was a period of substantial soil carbon loss, especially from the surface organic horizon. The rate of carbon loss essentially followed an exponential decay pattern, rapid at first, slowing to near zero over the experiment’s first decade. In phase II, soil respiration rates in the heated plots were generally equal to or less than those in the control plots. The transition from phase I to phase II was characterized by a depletion of a labile C pool (
18), which is considered to be the driver of reduced microbial biomass (
24).
Phase II appears to have been a period of microbial community reorganization, leading to changes in structure and function. During this time, soil heating reduced the abundance of fungal biomarkers and also caused a shift toward Gram-positive bacteria and especially actinobacteria (
17). Along with these structural changes, heating reduced the capacity of microbial community to utilize simple C substrates (
17).
The transition from phase II to phase III was characterized by a continued shift toward a more oligotrophic microbial community with increased diversity due to increased evenness (Pielou’s J statistic), reduced microbial biomass, and reduced fungal dominance, as evidenced by the same population size of bacteria and narrowed fungal/bacterial ratios for C-degrading genes (
20).
In phase III, soil respiration rates were higher in the heated plots than in the controls. This third phase appears to have been a period when recalcitrant substrates such as lignin became an important source of carbon for the microbial community. This is consistent with the evidence that during phase III, relative to controls, there was an increase by a factor of 4 in potential lignin-degrading enzyme (lignase) activity in surface soils from the heated plots (
21). As phase III transitioned toward phase IV, we observed a reduction in the relative abundance of lignin in the soil C pool in the surface mineral horizon (
19). At the same time, the magnitude of the carbon loss to the atmosphere through decomposition during this period may have been attenuated somewhat by a shift toward higher microbial carbon use efficiency for recalcitrant substrates in warmed soils relative to control soils (
23).
Our study is just 3 years into phase IV. This new phase may turn out to be another period of microbial community reorganization that will eventually transition to yet another phase of further carbon loss from decay of recalcitrant forms of SOM. Because recalcitrant SOM pools make up a substantial fraction of global soil carbon stocks (
25), small changes in the decay rates of these pools could result in a large self-reinforcing feedback to the climate system over multiple decades (
26). As a preliminary test of global significance, extrapolating our results to the world’s forests, we estimate a global aggregate soil carbon loss from the upper 1 m of soil (
27) over the 21st century of ~190 Pg C. This does not account for possible future climate-driven changes in plant-soil interactions that could affect the long-term balance between the formation and decomposition of SOM. Critical to this balance will be changes in the amount of fresh carbon transferred from plants to the soils as the world warms. Inputs of this fresh carbon can contribute to soil carbon sequestration, but they can also accelerate the decomposition of more recalcitrant forms of SOM through biological priming mechanisms (
28).
Our first-order estimate of a warming-induced loss of 190 Pg of soil carbon over the 21st century is equivalent to the past two decades of carbon emissions from fossil fuel burning (
29) and is comparable in magnitude to the cumulative carbon losses to the atmosphere due to human-driven land use change during the past two centuries (
30). A transfer of carbon of this magnitude from forest soils to the atmosphere in response to warming would amplify the mitigation challenge already faced by society. It is also important to recognize that a global-scale, microbially mediated feedback could be very difficult, if not impossible, to halt.
Acknowledgments
Supported by U.S. Department of Energy grants DE-SC0010740 and DE-SC0016590 and by NSF grants DEB 1237491 (Long-Term Ecological Research) and DEB 1456528 (Long-Term Research in Environmental Biology). We thank the team of scientists who worked with us on the study over the past 26 years, including J. Aber, T. Ahrens, C. Baldino, J. Blanchard, M. Bradford, L. Burrows, A. Burton, S. Butler, C. Catricala, R. Hanifin, T. Hill, J. Johnson, H. Lux, J. Mohan, S. Morrisseau, K. Newkirk, W. Peterjohn, A. Ricca, L. Scott, R. Smith, P. Steudler, and C. Vario. Data sets pertaining to this study are available on the Harvard Forest Online Data Archive (
http://harvardforest.fas.harvard.edu/harvard-forest-data-archive), data sets hf005 and hf171.
Nonnegligible role of warming-induced soil drying in regulating warming effect on soil respiration
Nonnegligible role of warming-induced soil drying in regulating warming effect on soil respiration
Enzai Du
State Key Laboratory of Earth Surface Processes and Resource Ecology, Faculty of Geographical Science, Beijing Normal University, Beijing, 100875, China.
Corresponding author: Enzai Du, Tel: +86-10-58808085, Email: [email protected]
Abstract
Based on results of a 26-year soil warming experiment (soil temperature being elevated by 5 °C) in a Harvard hardwood forest, Melillo et al. demonstrated a four-phase pattern of long-term warming effect on soil respiration, while the mechanisms were not fully elucidated because they neglected the indirect effect due to warming-induced soil drying. By showing a significant correlation between precipitation anomaly and inter-annual variation of warming effect on soil respiration, we suggest a nonnegligible role of warming-induced soil drying in regulating the long-term warming effect on soil respiration. Our analysis recommends further efforts to consider both the direct and indirect (i.e., warming-induced soil drying) warming effects to gain more in-depth understanding of the long-term soil C dynamics.
Understanding the long-term effect of climate warming on soil respiration is a prerequisite for the projection of future change in global soil carbon (C) pools and their consequent feedbacks to climate change. Based on results of a 26-year soil warming experiment (soil temperature being elevated by 5 °C) in a Harvard hardwood forest, Melillo et al. (1) proposed a hypothesis of four-phase pattern of soil warming effect on soil respiration (ΔRs = RsHeated−RsControl) and attributed it to soil substrate and microbial changes. However, the mechanisms were not fully elucidated because the authors improperly neglected the indirect effect due to warming-induced decline in soil moisture.
Experimental results have evidenced an essential role of soil moisture in regulating the warming effect on soil respiration (2-5). Lower soil moisture can reduce ΔRs, and even leads to negative ΔRs, when soil moisture becomes extremely limiting to plant roots and soil microbial organisms (3-6). In addition to the direct warming effect (1), warming-induced soil drying (2-6) may indirectly reduce ΔRs, especially in drier years (2-5). Moreover, soil moisture stress has been evidenced to reduce temperature sensitivity of soil respiration (7, 8). Warming-induced soil drying may have also contributed to a decrease in temperature sensitivity of soil respiration in heated plots, being observed to occur consistently across all four phases during the long-term soil warming experiment in Harvard forest [see Fig. 3 in Melillo et al. (1)]. Unfortunately, we cannot re-analyze the role of warming-induced soil drying due to unavailable soil moisture data (1).
As a primary source of soil moisture, precipitation is expected to buffer the soil moisture stress and result in a positive effect on inter-annual variation of ΔRs . By retrieving data on total annual precipitation (TAP) anomaly (showing no significant temporal trend, p=0.94) and ΔRs from Melillo et al. (1), we tested the effect of TAP anomaly on inter-annual variation of ΔRs. As ΔRs showed a significant decrease over time (ΔRs = –6.1547*Year + 12402, adjusted R² = 0.29, p=0.005) due to long-term microbial C substrate depletion (1, 9, 10), we thus used detrended (i.e., residual) ΔRs [being calculated as, measured ΔRsi – (–6.1547*Yeari +12402), where i indicates a specific year] to indicate the inter-annual variation of ΔRs. As expected, we found a significantly positive effect of precipitation anomaly on the inter-annual variation of ΔRs (slope = 0.15 g C m−2 year−1 per mm; adjusted R2=0.17, p=0.03). Our result suggests a significant role of soil moisture in regulating the warming effect on soil respiration. It further implies that warming-induced soil drying stress has very likely exerted a negative effect on the response of soil respiration to long-term warming treatment in Harvard forest.
The long-term experiment, with each plot covering a small area of 6 × 6 m2 (1), generally represented a scenario of soil C depletion under soil warming, while it failed to demonstrate the story in a whole C cycle perspective. The small plot size was not able to track the direct and indirect (i.e., warming-induced soil drying) warming effects on forest primary productivity, which could consequently provide fresh C substrates for soil respiration (11). In that case, negative feedbacks between C substrate limitation and respirational C loss likely dominate a long-term reduction of ΔRs (9, 10). By inducing a reduction of soil moisture (2-6, 12, 13), soil warming also indirectly regulates ΔRs via inducing changes in the composition and activity of microbial community (5, 14). For instance, warming-induced soil moisture stress may have contributed substantially to the declining trends of ΔRs during the relatively drier periods 1997-2002 and 2015-2016 [see Fig. S2 in Melillo et al. (1)].
As the indirect effect due to warming-induced soil drying stress has not been addressed properly, the mechanisms of the four-stage pattern of long-term warming effect on soil respiration are likely not fully elucidated (1). Moreover, Melillo et al. (1) assumed that microbial respiration accounted for a constant proportion (approximately 2/3) of soil respiration, while they ignored temporal changes in the microbial proportion over the four-phase pattern of soil respiration. Therefore, in-depth understanding of the long-term warming effect on soil C dynamics calls for further efforts to consider both the direct and indirect (i.e., warming-induced soil drying) effects in a whole C cycle perspective at an ecosystem scale.
Acknowledgements: This work was supported by the National Natural Science Foundation of China (Nos. 41630750 & 31400381) and State Key Laboratory of Earth Surface and Resource Ecology (No. 2017-ZY-07).
References:
1. J. M. Melillo et al., Science 358, 101-105 (2017).
2. J.C. Carey et al., PNAS 113, 13797-13802 (2016).
3. G. Li et al., Forests 8, 56 (2017).
4. W. Liu et al., Global Change Biol. 15, 184-195 (2009).
5. S. D. Allison et al., Global Change Biol. 14, 2898-2909 (2008).
6. J. Harte et al., Ecol. Appl. 5, 132–150 (1995).
7. R.S. Jassal et al., Global Change Biol. 14, 1305–1318 (2008).
8. A. Gershenson et al., Global Change Biol. 15, 176-183 (2009).
9. M. U. Kirschbaum Global Change Biol. 10, 1870–1877 (2004).
10. A. L. Romero-Olivares et al., Soil Biol. Biochem. 107, 32–40 (2017).
11. A. Hursh et al., Global Change Biol. 23, 2090–2103 (2017).
12. J. C. Yuste et al., Tree physiol. 23, 1263–1270 (2003).
13. V. Suseela et al., Global Change Biol. 18, 336–348 (2012).
14. S. Manzoni et al., Ecology 93, 930–938 (2012).