Main and interactive effects of warming and canopy or rainfall on soil respiration
Warming generally increased long-term rates of soil respiration (P < 0.001, Fig. 1, Table 1, and Extended Data Fig. 1 and Fig. 2). Specifically, on average across all observations in all years in both experiments, soil respiration increased by 6.5% and 17.2% under + 1.7°C and + 3.3°C, respectively (Extended Data Fig. 3a). The stimulation varied among years from lows of near zero to a high of 17.7% and 34.3% in 2019 for + 1.7°C and + 3.3°C, respectively (Extended Data Fig. 4). On average across all years, there was no difference in soil respiration in open and closed canopies (Table 1); soil respiration was lower by 4.7% in rainfall reduction than ambient rainfall treatment (Extended Data Fig. 3b).
Table 1
Repeated measures linear mixed-effects model analyses of interactive effects of warming, canopy condition, and rainfall reduction on soil respiration.
Warming and canopy experiment | Warming and rainfall experiment |
Effect | DF | F | P > F | Effect | DF | F | P > F |
Warming | 1 | 159.2 | < 0.001 | Warming | 1 | 40.7 | < 0.001 |
Canopy | 1 | 0.9 | 0.37 | Rainfall | 1 | 12.8 | < 0.001 |
Year (yr) | 12 | 9.9 | < 0.001 | Year (yr) | 9 | 19.9 | < 0.001 |
Site | 1 | 4.1 | 0.07 | Site | 1 | 3.2 | 0.14 |
Warming × Canopy | 1 | 10.1 | < 0.01 | Warming × Rainfall | 1 | 4.2 | 0.04 |
Warming × yr | 12 | 5.4 | < 0.001 | Warming × yr | 9 | 4.1 | < 0.001 |
Canopy × yr | 11 | 8.4 | < 0.001 | Rainfall × yr | 9 | 3.6 | < 0.001 |
Warming × Site | 1 | 1.0 | 0.32 | Warming × Site | 1 | 16.4 | < 0.001 |
Canopy × Site | 1 | 0.0 | 0.88 | Rainfall × Site | 1 | 0.1 | 0.72 |
yr × Site | 12 | 3.3 | < 0.001 | yr × Site | 9 | 2.9 | < 0.01 |
Warming × Canopy × yr | 11 | 4.0 | < 0.001 | Warming × Rainfall × yr | 9 | 3.3 | < 0.001 |
Warming × Canopy × Site | 1 | 12.9 | < 0.001 | Warming × Rainfall × Site | 1 | 5.1 | 0.02 |
Warming × yr × Site | 12 | 1.9 | 0.03 | Warming × yr × Site | 9 | 3.9 | < 0.001 |
Canopy × yr × Site | 11 | 1.0 | 0.41 | Rainfall × yr × Site | 9 | 2.1 | 0.03 |
Warming × Canopy × yr × Site | 11 | 6.0 | < 0.001 | Warming × Rainfall × yr × Site | 9 | 0.9 | 0.50 |
Warming: experimental warming level; Canopy: canopy condition; Rainfall: experimental rainfall amount; Site: field site; F: F ratio. Significant effects (P < 0.05) are highlighted in bold. The measurements of soil respiration from all campaigns were used in the analyses. The values of increased temperature were used to represent the experimental warming level in the analysis (0, 1.7, and 3.3 for ambient, + 1.7°C, and + 3.3°C treatment, respectively); and year, canopy condition, and rainfall amount were nominal factors in the mixed model. |
In the “warming and canopy condition” experiment, warming was a main factor affecting soil respiration (Table 1), and the effect of warming on soil respiration varied among years, canopies, and sites (significant interactions, Table 1 and Fig. 1). Similarly, in the “warming and rainfall amount” experiment, warming and rainfall reduction individually and interactively affected soil respiration, and their effects varied with site and year (significant interactions, Table 1).
Whether warming increases or decreases soil respiration may depend in part on the net effect of the direct stimulatory effect and the potential negative effect of soil drying. Our long-term data showed that warming treatment increased soil temperature (Fig. 2 and Extended Data Fig. 5 and decreased soil moisture (Fig. 2 and Extended Data Fig. 6) regardless of canopy condition, rainfall reduction treatment, and site (Extended Data Fig. 7). Since soil respiration showed a greater response to realized soil temperature than soil moisture (Fig. 2), we found overall positive effects of warming on soil respiration (Extended Data Fig. 3a).
As two of the only other long-term studies of soil respiration response to warming in forests, the 26-year warming experiment in the Harvard Forest and 16-year Flakaliden warming experiment in spruce forest in northern Sweden, provide useful context. Those studies, which elevated soil temperature to 5°C above ambient, but did not warm aboveground, found that warming increased soil respiration on average at Harvard Forest by 8.5%21 and by 2% in years 14–16 of the experiment at Flakaliden26. In contrast, combined plant and soil warming of + 3.3°C in the present B4WarmED study over all years stimulated approximately two and eight times more soil respiration than + 5.0°C warming of only soils in the Harvard Forest and long-term results from Flakaliden. Moreover, warming of + 1.7°C in B4WarmED resulted in almost as large, and much larger, respectively, average stimulation of soil respiration as + 5.0°C soil warming at Harvard Forest and in years 14–16 at Flakaliden. Moreover, we observed no tendency for warming stimulation of soil respiration averaged across sites and sub-experiment to diminish over time (Extended Data Fig. 4).
Two potential explanations for the much greater effect per unit C warming of combined above and belowground warming are as follows: first, when compared to belowground warming alone, simultaneous above- and belowground warming might further increase plant growth including root biomass, and thus enhance autotrophic root respiration, an important component of soil respiration27; and second, combined above- and belowground warming could increase plant aboveground and belowground litter and exudate inputs, increasing heterotrophic use of plant-derived C. Additionally, although respiration at our site being more responsive per degree warming than at the warmer Harvard Forest site28 is consistent with past evidence that soil organic C decomposition is more sensitive to warming in colder regions, the lesser sensitivity at colder Flakiliden26 is not.
Mid-term (2–6 years) studies in forests also found warming generally increased soil respiration in those early years by more than the longer-term studies after 10 years. In studies in tropical forest in Panama (+ 4°C for 2 years)29, deciduous forest in Japan (+ 3°C for 4 years)30, boreal pine forest in Finland (+ 2 to 4°C for 4 years)31, temperate oak forest in China (+ 1.6 to 2.1°C for 5 years)32, mixed coniferous forest in California (+ 4°C for 5 years)33, and Norway spruce forest in Austria (+ 4°C for 6 years)16, warming increased soil respiration by 55%, 16%, 35%, 29%, 30%, and 37%, respectively. It should be noted that the differences in climate zone (e.g. boreal, temperate, and tropical zones), warming method (e.g. soil surface warming, whole-soil warming, only aboveground warming, and both above- and below-ground warming) and level, and longevity among these mid-term studies make it difficult to compare their results, and/or their results in relation to the three longer-term studies.
Pathways regulating effects of warming, canopy condition, and rainfall amount on soil respiration
The effects of warming, rainfall reduction, and canopy removal logically might occur at least in part through changes in soil temperature and soil moisture. Thus we used structural equation modeling to assess direct and indirect effects of treatments including via these two pathways (Fig. 2). Warming and open canopy conditions increased soil temperature, ultimately increasing soil respiration (Fig. 2a and 2b); warming and rainfall reduction decreased soil moisture, ultimately decreasing soil respiration (Fig. 2c and 2d). Significant three-way interactions between warming, site, and canopy condition or rainfall reduction (Table 1) were mainly explained by the change in soil moisture (Fig. 2), which was partly consistent with Hypothesis 1. It should be noted, however, that the interactions between warming and canopy condition or rainfall reduction were more complex than proposed in Hypothesis 1, in that warming and canopy condition interacted with site in affecting soil respiration for the “warming and canopy condition” experiment (P < 0.001, Table 1). Generally speaking, soil respiration increased with increasing temperature except for in the open canopy condition at HWRC site, which showed insignificant differences in soil respiration between + 1.7°C and + 3.3°C (Extended Data Fig. 8a). This lack of temperature effect is likely because plots in open canopy condition at HWRC had significantly lower soil moisture (Extended Data Fig. 7a). The negative effects of both warming and open canopy condition on soil moisture led to the lowest soil moisture under + 3.3°C, and thus, the negative effects of soil drying offset the positive effect of warming on soil respiration. However, at CFC, canopy condition did not significantly affect soil moisture (Extended Data Fig. 7a). Therefore, soil respiration was stimulated at both levels of warming in both canopy treatments at CFC. Overall, because of the different effects of canopy condition on soil moisture at the two sites (insignificant difference in soil moisture between open and closed canopy condition at CFC but significantly higher soil moisture under closed than open canopy condition at HWRC, Extended Data Fig. 7a), soil respiration increased with increasing warming level regardless of canopy condition at CFC, but not at HWRC (Extended Data Fig. 8a), which was partly consistent with Hypothesis 1.
A similar mechanism (e.g. warming effects on soil respiration mediated by rainfall reduction effects on soil moisture) appeared to be at work for the “warming and rainfall amount” experiment. For example, significant interactions between warming, rainfall reduction, and site were found (P < 0.05, Table 1). Specifically, when compared to CFC, smaller differences in soil respiration between warming levels were found at HWRC (Extended Data Fig. 8b). The absolute decline in soil moisture caused by rainfall reduction was similar between the two sites (Extended Data Fig. 7b). However, because HWRC had lower ambient soil moisture overall than CFC, rainfall reduction led to much lower levels of soil moisture at HWRC, likely inhibiting soil respiration response to warming more at that site. This interpretation is further supported by the results of SEM analysis, which showed a stronger correlation between soil moisture and respiration at HWRC than CFC (Fig. 2c and 2d). Overall, our findings suggest that canopy condition and rainfall reduction can significantly change warming effects on soil respiration by altering soil moisture.
Role of ambient soil moisture in warming effects on soil respiration
When compared to how warming affects soil respiration by altering soil temperature and moisture, little is known about how fluctuations in ambient soil temperature and moisture (e.g. soil temperature and moisture under ambient warming treatment) modulate warming effects on soil respiration, partly due to the few measurements from long-term studies. Across 13 years of our study, ambient soil moisture and temperature from 22,386 observations ranged from 0.05 m3 m− 3 to 0.44 m3 m− 3 and from − 0.14°C to 27.74°C during times when soil respiration was measured, respectively. The large variation in ambient soil moisture and temperature could help us advance our understanding of their roles in mediating warming effects on soil respiration in boreal forests. Generally speaking, when the soil was extremely dry (e.g. 0th − 10th percentile ambient soil moisture), warming effects on soil respiration were minor or even negative, but became positive when ambient soil moisture was higher (Fig. 3), which supports Hypothesis 2. In addition, warming effects on soil respiration remained the same or decreased when the soil was extremely moist (e.g. 90th − 100th percentile) compared to intermediate moisture levels (e.g. 10th − 90th percentile). When soils are very dry, the growth of both plants and soil microorganisms is significantly inhibited13,14,22. The decrease in soil moisture induced by warming can further exacerbate the water limitation of plants and microorganisms. Therefore, insignificant or even negative effects of warming on soil respiration occurred when ambient soil moisture was extremely low. When soils are very wet, soil microorganisms are oxygen limited4,34, thus weakening the positive effect of warming on soil respiration. Given that extreme precipitation events (e.g. flooding and prolonged drought) are projected to become more frequent and severe in the future35,36, our results suggest that warming effects on soil respiration may fluctuate from positive to negative as ambient soil moisture levels vary. By overlooking the negative effect of warming on soil respiration in extremely dry soils, we would overestimate the overall positive effect of warming on soil respiration and thus the positive climate feedback, especially in regions that are projected to experience more drought events in the future.
Long-term interannual temporal pattern of warming effects on soil respiration
Significant interannual variability in soil respiration was found for all treatments at both sites (Figs. 1 and 4). Overall, the response of soil respiration to warming in our experiment exhibited either two- (decrease and then increase) or three-phase (decrease-increase-decrease) patterns (Fig. 4 and Extended Data Fig. 9), similar to the results from the study at Harvard Forest, MA, USA21. For example, in closed canopy plots, the effects of both warming levels on soil respiration decreased for the first four years at both sites, then increased in years five and six, and then either increased continuously (at CFC) or showed another decline (at HWRC). Under open canopy condition and ambient rainfall at both sites, the effects of both warming levels were initially positive, generally peaked in year 2, then decreased to neutral or negative with a trough after 4–6 years, then increased to positive values again. For open canopy plots with rainfall reduction, warming effects started out positive and declined at CFC but were small and did not change a lot over time at HWRC. Overall, our findings indicate that the temporal pattern of long-term warming effects on soil respiration in boreal forests is complex, varying significantly with canopy condition, rainfall amount, and site.
To test Hypothesis 3, we determined the relationships between annual ambient soil moisture, temperature, precipitation, and the changes in soil respiration caused by warming. Unlike short-term fluctuations in soil moisture, which were significantly related to the effects of warming on soil respiration across individual measurements, the annual average change in soil respiration with warming treatment was not significantly correlated with annual precipitation, ambient soil moisture, or temperature, under almost all situations. These findings indicate that the long-term interannual patterns of warming effects on soil respiration in boreal forests were not explained by the interannual variations of either precipitation, soil moisture, or temperature, which is inconsistent with Hypothesis 3. Thus, interannual changes in soil biochemical or microbial variables, rather than in weather, likely contributed to the interannual pattern of warming effects on soil respiration21.
For the “warming and canopy condition” experiment, warming effects on soil respiration declined for all treatment combinations over the first few years (Fig. 4). This decline might be driven by the depletion of more easily decomposable C under warming over time21. The duration of the decrease in warming effects on soil respiration was much shorter under closed than open canopy conditions especially at HWRC, possibly due to the relatively higher soil moisture but lower soil temperature under closed than open canopy conditions (Extended Data Fig. 7a and 7c). The same increase in temperature has been shown to stimulate decomposition rates more in cold than warm regions28; furthermore, higher soil moisture may promote labile C decomposition. Therefore, colder and wetter soils under closed than open canopy conditions might have resulted in a faster labile C decomposition rate, shortening the duration of declining warming effects on soil respiration. Despite no direct evidence from the present study, changes in the soil microbial community as were found in the Harvard Forest21 might explain why warming effects on soil respiration oscillated over time in our study. Lower litterfall amount and sparse grass growth under closed canopy condition at HWRC might explain why warming effects on soil respiration decreased again after the decrease-increase pattern under this condition. Overall, our results indicate that the interannual pattern of warming effects on soil respiration varies with canopy condition.
For the “warming and rainfall amount” experiment, the relatively modest temporal changes in warming effects on soil respiration under rainfall reduction than ambient rainfall treatment (Fig. 4) suggest that by causing soil water deficits, rainfall reduction might slow the depletion of soil labile C and changes in soil microbial community caused by long-term warming. Since global climate change is projected to encompass both warming and changes in rainfall amount in the future8, the role of rainfall reduction in mediating the temporal pattern of warming effects on soil C cycling should be considered in future experimental and modeling studies.