Effects of Soil Warming and Nitrogen Addition on Soil Respiration in a New Zealand Tussock Grassland

Soil respiration (R S) represents a large terrestrial source of CO2 to the atmosphere. Global change drivers such as climate warming and nitrogen deposition are expected to alter the terrestrial carbon cycle with likely consequences for R S and its components, autotrophic (R A) and heterotrophic respiration (R H). Here we investigate the impacts of a 3°C soil warming treatment and a 50 kg ha−1 y−1 nitrogen addition treatment on R S, R H and their respective seasonal temperature responses in an experimental tussock grassland. Average respiration in untreated soils was 0.96±0.09 μmol m−2 s−1 over the course of the experiment. Soil warming and nitrogen addition increased R S by 41% and 12% respectively. These treatment effects were additive under combined warming and nitrogen addition. Warming increased R H by 37% while nitrogen addition had no effect. Warming and nitrogen addition affected the seasonal temperature response of R S by increasing the basal rate of respiration (R 10) by 14% and 20% respectively. There was no significant interaction between treatments for R 10. The treatments had no impact on activation energy (E 0). The seasonal temperature response of R H was not affected by either warming or nitrogen addition. These results suggest that the additional CO2 emissions from New Zealand tussock grassland soils as a result of warming-enhanced R S constitute a potential positive feedback to rising atmospheric CO2 concentration.


Introduction
Soils contain a pool of carbon approximately double that stored in terrestrial biomass [1]. Soil respiration (R S ), the primary pathway for return of soil carbon to the atmosphere, is increasing globally by 0.1 Pg C y 21 at present [2]. This increase, hypothesised to be a result of global warming, is concerning as temperatures are expected to rise by as much as 6.4uC over the next century [3]. Coupled climate models indicate a likely soildriven positive feedback to climate change, although uncertainty remains in the magnitude of this feedback [4,5].
Numerous warming experiments have investigated the impacts of long-term climate warming on carbon cycling, suggesting that, on average, warming of 0.3-6.0uC increases soil respiration (R S ) by 20% [6]. However, several notable examples have shown the effect of warming on R S to be only transient [7,8]. Mechanisms for this acclimation of R S to prolonged warming include depletion of labile carbon substrates [8][9][10], changes to the microbial community structure [7,11], physiological acclimation of soil microbes [12], reduction in root biomass [13] and reduction in the specific root respiration rate [14]. Acclimation of soil respiration may limit potential soil carbon loss as a result of climate warming.
Global change scenarios also suggest that nitrogen cycling in terrestrial ecosystems will be altered. Nitrogen deposition due to crop fertilisation and fossil fuel combustion currently exceeds terrestrial nitrogen fixation and is expected to increase in the future [15]. As warming also increases nitrogen mineralization [6], there exists the possibility for synergistic effects of warming and anthropogenic nitrogen deposition on plant-available nitrogen.
While warming-induced increases in R S represent a likely positive feedback to rising atmospheric CO 2 concentration, enhanced nitrogen deposition has been suggested as a possible mitigating factor due to negative impacts of nitrogen addition on R S [16]. The findings from nitrogen addition experiments in forests suggest that reduction in R S may represent a carbon offset equivalent to the nitrogen fertilisation effect on primary production. As well, reductions in R S have been observed in grasslands as a result of nitrogen addition [17,18]. Due to feedbacks between the nitrogen and carbon cycles, nitrogen availability will likely influence the magnitude of the terrestrial feedback to rising atmospheric CO 2 concentration [19].
The net response of R S to warming and nitrogen addition depends largely on the combined response of its components, autotrophic soil respiration (R A ) and heterotrophic soil respiration (R H ), which are likely to have different responses to environmental drivers. Autotrophic respiration refers to respiratory activity of roots and associated rhizosphere microbes, while R H refers to soil organic matter decomposition by soil microbes [20]. The important distinction between R A and R H is that the former represents respiration of carbon recently assimilated by plants, whereas the latter releases carbon that may have residence times in the soil reaching millennia [21].
Heterotrophic respiration is widely expected to increase under warming scenarios [22,23]. Several field warming experiments demonstrated such increases [24,25]. While warming generally increases R H , nitrogen addition can decrease microbial biomass [26], potentially reducing R H . This reduction in R H may explain the overall decrease in R S observed in response to nitrogen addition [16].
In this study, we investigated the impacts of soil warming and nitrogen addition, as well as their interaction, on R S and its components, R A and R H . Such multifactor experiments are important to improve the predictive ability of coupled-climate models, as single factor experiments may fail to predict interactive effects of global change drivers [27,28]. Likewise, partitioning the autotrophic and heterotrophic components of R S can lead to greater mechanistic understanding of the response of R S to environmental drivers [29].
Native tussock grassland was selected as a model system because grasslands are a widespread and important store of carbon in New Zealand [30], and globally [31]. Soil respiration and R H were measured over a period of 27 months with the objective of determining the likely feedback effect that increases in R S in grasslands will have on rising atmospheric CO 2 concentration in response to soil warming and nitrogen addition.

Study site
This study was conducted at the Cass Warming Experiment at the University of Canterbury Cass Field Station in central South Island, New Zealand (43.03u S, 171.75u E, 590 m a.s.l.). The site was constructed in January 2009 in an area of tussock grassland. Soils at the site are classified as acidic allophane brown by New Zealand Soil Classification (Typic Dystrochrept by USDA) [32,33]. Prior to this study, vegetation and the top 200 mm of topsoil were removed, twenty 12.25 m 2 plots were laid out and 90 m of resistance heating cable (Argus Heating, Ltd., Christchurch, New Zealand) were arranged in rows with 140 mm spacing between cables in each plot to achieve a heating density of 76 W m 22 [34]. Dummy cables were arranged similarly in unheated plots. The cables were then covered with 200 mm of topsoil and the native New Zealand tussock grasses Festuca novae-zelandiae (50 individuals per plot), Poa cita (50 per plot), Chionochloa rigida (22 per plot), and Chionochloa flavescens (12 per plot) were planted.
Five plots were assigned to each of four treatments: control, warming only, nitrogen addition only and combined warming and nitrogen addition (Appendix S1). In each of 10 plots designated for warming, three thermocouples (Type-E, Campbell Scientific, Logan, UT, USA) were buried to a depth of 100 mm in a stratified design which captured a range of horizontal distances from heating cables (directly above, one quarter of the distance between two cables and the midpoint). In each of the control plots, one thermocouple was buried to 100 mm soil depth. The heating cables were switched on and off to maintain a 3uC difference between the average of the three thermocouples in warmed plots and the nearest un-warmed plot. Warming was controlled by a datalogger (CR1000X, Campbell Scientific, Logan, UT, USA) and hourly average plot soil temperatures were recorded. An auxiliary weather station measured hourly average air temperature, soil temperature and volumetric water content at 100 mm depth.
Nitrogen addition began in February 2009. Nitrogen was applied as calcium ammonium nitrate at a rate of 10 kg N ha 21 five times throughout the growing season to achieve a total amendment of 50 kg N ha 21 y 21 . For each plot, the fertiliser was dissolved in 4 L water and distributed using a watering can over both plants and soil. The continuous 3uC warming treatment was started in July 2009. Two plots, one each of the warming only and combined warming and nitrogen addition treatments, were subsequently dropped from analyses due to malfunction of the heating cables.

Respiration measurements
Measurements of soil respiration were carried out over a 27 month period beginning in August 2009 (winter) and continuing through October 2011 (spring). Six 100 mm diameter polyvinyl chloride measurement collars were installed to a soil depth of 70 mm in each plot. The rate of soil respiration in each collar was measured at 2-4 week intervals using a portable respiration system (SRC-1 and EGM-4, PP Systems, Amesbury, MA, USA). An additional two measurement collars were installed in each plot to a soil depth of 300 mm in order to exclude roots and provide an estimate of heterotrophic respiration (R H ). These deep collars extended into the clay subsoil, limiting potential root growth into the soil beneath the collar. In contrast, the shallow collars were inserted to a depth that would allow root infiltration beneath the collar while providing a seal with the soil surface. The collars remained in place for the duration of the experiment to avoid soil or root disturbance. Measurement using the deep collars began in January 2010. Simultaneous with each soil respiration measurement, soil temperature and soil water content at 50 mm depth were measured using a thermocouple (Type-E, Omega Engineering, Ltd, Stamford, CT, USA) and a soil moisture sensor (Theta Probe type ML1 and ML2, DeltaT Devices, Cambridge, UK), respectively.

Substrate addition
Availability of labile substrates in the soil is important in regulating R S [35]. In order to assess levels of substrate limitation induced by warming and nitrogen addition treatments and the presence of roots, a substrate addition experiment was carried out in late-October 2011. In each of 16 plots representing four replicates for each treatment, two pairs of soil respiration measurement collars were selected: one pair with roots present and another pair with roots excluded. One collar from each pair was selected randomly for substrate addition. All collars were measured immediately prior to substrate addition. Subsequent to initial measurement, the two collars from each plot selected for substrate addition were amended with 20 ml of 0.2 M sucrose solution (an amount approximately equivalent to 10 days of carbon losses from R H ). In order to ensure that the sucrose solution infiltrated beyond the soil surface, 5 ml were injected with a syringe to a depth of 25 mm at four locations within each collar. The collars designated as controls were treated similarly with water. Soil respiration was then measured in each collar at 30 min, 1 h, 2.5 h, 4 h and then at 4-8 h intervals until the substrate response was no longer evident. Substrate-induced respiration (S I ) was calculated for each pair of collars as the difference between respiration rates of the control and substrate-added collars, as a proportion of the rate for the control treatment.

Soil analyses
Soils were sampled in January 2010, March 2011 and March 2012. Three 54 mm diameter soil cores were taken to a depth of 100 mm in each plot and the soil was homogenized into a single sample. Roots were removed by 8 mm sieve and dried at 60uC. As the grass roots were very fine, in 2012 a subsample was removed Effect of Warming and Nitrogen on Soil Respiration PLOS ONE | www.plosone.org from the whole sample and washed over a 650 mm sieve to obtain root biomass. Microbial biomass was estimated using the fumigation-extraction technique adapted from Vance et al. [36]. Remaining soil was dried at 60uC, passed through a 2 mm sieve to remove remaining roots and ground in a ball mill. Samples were then analysed for total carbon and nitrogen concentration on an elemental analyser (CNS2000, Leco, St. Joseph, MO, USA).
Between 20 September 2011 and 22 October 2011, plant available nitrogen was estimated using ion exchange membranes (PRS probes, Western Ag Innovations, Saskatoon, Canada). The PRS probes were installed to a depth of 100 mm at three locations in each plot. Following a one-month burial period, probes were removed, rinsed with deionized water and returned to Western Ag Innovations for analysis of NH 4 + and NO 3 2 .

Statistical analyses
The effects of warming and nitrogen addition on seasonal measurements of R S , R H and soil water content were assessed using linear mixed-effects models conducted in the 'nlme' package [37] in R version 2.12.1 [38]. Warming, nitrogen addition and measurement date, along with their interactions were included as fixed effects, while measurement collars nested within plots were included as random effects to account for the non-independence of multiple samples through time and multiple collars per plot. Residual analyses were undertaken and log transformation was used for R S and R H , to correct for heteroscedasticity. The effects of warming and nitrogen on the proportion of R S constituted by R H (f RH ) were similarly assessed by treating plot averages of f RH on a given date as a sample and evaluating random effects at the plot level.
Temperature responses of R S and R H were fitted to an Arrhenius-type curve [39], modified with a soil water content response function [40]: where T S is soil temperature (K), R 10 is the basal respiration rate at 10uC, E 0 is the activation energy of enzymatic reactions, h is the soil volumetric water content and a and b are parameters that determine the shape of a sigmoidal response of respiration to soil water content.
Nonlinear mixed-effects models (also conducted in the 'nlme' package for R) were used to fit Equation (1) initially to measurements of R S , and subsequently to R H . First, the effect of roots on the temperature response of R S was investigated by testing how root presence (as a fixed effect) altered parameter values for R 10 and E 0 . Subsequently, temperature responses of R S and R H were investigated in separate models, with the latter substituting R H in place of R S in Equation (1). Warming and nitrogen addition, as well as their interactions, were investigated as fixed effects on R 10 and E 0 for both the R S and R H models. For all the above nonlinear models, measurement collars nested within plots were evaluated as random effects. Fixed and random effects on a and b were not evaluated, as few measurement dates occurred under water-limited conditions, and a generic water content response curve which limited respiration when soil water content was less than 0.2 m 3 m 23 was deemed appropriate based on analysis of residuals of a temperature-only model.
The final fixed effects structure was determined by first constructing a maximal model which included the presence of plant roots, warming, nitrogen and all interactions. A power variance function was fitted in order to correct for heteroscedasticity [41]. To account for autocorrelation in repeated measurements of the same collar, a first order autoregressive structure was used [42]. Fixed effects and interactions were removed iteratively based initially on their p-values and the best fit model was selected. During each step of this procedure, model comparisons were undertaken using a likelihood ratio test and selection of the bestfitting model was achieved through minimisation of Akaikie's Information Criterion (AIC). To test for potential acclimation of R S to warming, the data were bisected such that the first full year of measurement was separated from the second. The interaction between measurement year and the warming treatment was added as a fixed effect on both R 10 and E 0 parameters. This interaction was tested for significance to determine whether the warming effect was consistent across measurement years, with a significant negative interaction term for either parameter indicating acclimation.
Soil carbon content, soil nitrogen content, microbial biomass, plant available nitrogen and substrate induced respiration (S I ) were all assessed by multi-way ANOVA, with temperature and nitrogen treatments, as well as their interaction, as factors. For S I , the maximum value recorded for each pair off collars over the measurement period was tested. For those variables that were measured repeatedly (soil carbon, soil nitrogen, microbial biomass carbon), a separate ANOVA was conducted for each time point.

Seasonal variation of soil temperature and water content
Average soil temperature over the entire 27 month measurement period was 9.6uC. Warming increased soil temperature by an average (6SE) of 3.160.2uC over the course of the experiment (Fig. 1). Soil water content varied seasonally, falling below 0.

Seasonality of R S
Soil respiration showed a seasonal pattern driven primarily by seasonal temperature (Fig. 1, Fig. 2A). As such, measurement date had a significant effect in the linear mixed-effects model (p,0.0001, Table S2). Heterotrophic respiration showed a similar response to seasonal temperature and, as such, measurement date was significant in predicting R H (p,0.0001, Table S3, Fig. 2B). The average proportion of R S constituted by R H (i.e., f RH ) was 71% (Fig. 2C). No apparent seasonal pattern was observed in f RH . Both R S and R H were sensitive to soil water content, with a reduction in respiration rate observed below 0.2 m 3 m 23 soil water content. This was particularly evident on 5 March 2010 and 12 December 2010, the driest measurement dates (Fig. 2).

Effects of warming and nitrogen addition on R S
On average, R S was increased by 41% due to warming (p,0.0001) and by 12% due to nitrogen addition (p = 0.004). The treatments combined additively, as no significant interaction was observed. Warming significantly increased R H by 37% (p = 0.014), though nitrogen did not significantly affect R H (p = 0.798), nor was there a significant interaction between the warming and nitrogen treatments. The proportional contribution of R H to R S was reduced to an average of 59% by nitrogen addition, although this reduction was of marginal significance (p = 0.051, Table S4) due to high variability in f RH . Warming had no significant impact on f RH .
Responses of R S and R H to added substrate were highly variable across replicates, with values of S I ranging from a 0.8 to 4.75 fold increase in respiration. There was no effect of the presence of roots on S I (p = 0.487). For R S , maximum S I was decreased significantly by warming (p = 0.040, Table 1). In contrast, nitrogen addition did not affect S I significantly (p = 0.146), nor was there any significant interaction with warming (p = 0.835).

Modelled temperature responses of R S
The presence of roots increased R 10 significantly (p,0.0001, Table S5), but had no effect on E 0 (p = 0.340). As both R S and R H were shown to be sensitive to soil water content, addition of the soil water content response function in Equation (1) resulted in a significant improvement in model fit over a temperature-only model (DAIC = 1350, p,0.0001). The best-fit model of the temperature response of R S included warming and nitrogen addition as fixed effects on R 10 , which was increased significantly by both warming (p,0.0001, Table S6) and nitrogen addition (p,0.0001, Table 2). E 0 was unaffected by either warming or nitrogen and there were no significant interactions between warming and nitrogen for either R 10 or E 0 , so they were removed from the model. The inclusion of a warming by measurement year interaction for R 10 resulted in small (0.02 mmol m 22 s 21 ) and marginally significant decrease in R 10 during the second year of warming (p = 0.065, Table S7). Inclusion of the temperature by year interaction for E 0 resulted in a non-significant interaction term (p = 0.347), indicating little effect of treatment time on the temperature response of R S . For R H , all treatments exhibited a single temperature response curve regardless of treatment ( Table 2).

Variation of soil properties and microbial biomass
The average soil carbon concentration was 43 g kg 21 and this was not affected significantly by either warming or nitrogen addition (Table 1). Likewise, total nitrogen concentration, which averaged 3.4 g kg 21 , was unaffected by any of the treatments. The average microbial biomass carbon was 646 mg kg 21 . This did not change significantly under the main effects of warming and nitrogen addition. However, a significant negative interaction between warming and nitrogen addition was observed on the final soil sampling date (p = 0.023) indicating a reduced microbial biomass under combined warming and nitrogen addition. Cumulative exchange of plant available nitrogen, as estimated from the one month burial of PRS probes, was 8 mg N m 22 on average and was not significantly different between treatments. At the time of the final soil sample, average (6 SE) root biomass in the top 100 mm of soil was 465626 g m 22 . Root biomass was highly variable and no treatment differences were detected.

Discussion
The average (6 SE) soil respiration rate measured in the control plots over the course of this study was 0.9660.09 mmol m 22 s 21 . This value falls well within the range reported for temperate grasslands [43]. The relative contribution of R H to R S of 71% for control plots was very close to the average for non-forest ecosystems of 63% [20] and also agreed well with another temperate grassland warming experiment where R H contributed 56-72% of R S annually [25]. The 3uC warming treatment led to an increased rate of both R S and R H over the 27 month measurement period. The average 41% increase in R S due to warming falls well within the reported range of a 25% reduction to a 45% increase in R S , as a result of experimental warming [6]. Likewise, warming-induced enhancement of R H has been observed in other grassland warming experiments [25]. Warming treatment had no significant impact on f RH indicating that R A and R H were similarly sensitive to temperature. This is supported by the results of the temperature response curve fitting, in which similar values of E 0 were obtained for both R S and R H . However, the temperature response analysis did reveal a slightly higher basal respiration rate, R 10 , for warmed soils. As this increase in R 10 was only evident in R S , and not R H , we must assume that R A is responsible for the increase in R 10 . While an increase in the basal rate of R A as a result of warming would appear to contradict the previous finding that f RH was unaffected by warming, measurements of f RH were highly variable and, as they were calculated from plot averages, subject to within plot temperature variation. Thus, our results are consistent with a slight increase in root activity in the warmed plots. As no differences in root biomass were observed as a result of warming, this increase in R A may be due to increased specific root respiration. Further, long-term warming has been associated with increased root exudation [44], and labile carbon in these exudates may have stimulated Effect of Warming and Nitrogen on Soil Respiration PLOS ONE | www.plosone.org rhizosphere microbial activity leading to an increase in R A , as measured by the root exclusion approach.
The acclimation of R S to warming frequently observed in many long-running soil warming experiments [7,8] was absent in our study. We expected that acclimation would result in a significant, negative warming by measurement year interaction (i.e., a decrease in R 10 or E 0 relative to the control during the second year of warming). However, we observed only a small, marginally significant treatment by measurement year interaction effect for R 10 . One possible explanation for this lack of acclimation is the relatively short duration of this experiment. Physiological acclimation of roots and soil microbes should occur rapidly compared with the duration of our experiment, though changes to biomass and soil carbon pools may take longer. Another explanation is that acclimation has been linked to depletion of labile carbon substrates [9,10]. As the study site was recently cleared of vegetation and soil structure was disturbed, it is likely that the size of the labile carbon pool was reduced as a result of the disturbance. Our measurements would then reflect the temperature response of decomposition of more recalcitrant soil organic matter in the absence of a large labile carbon pool to which size adjustments can occur. As the system advances and labile carbon accumulates, acclimation may become evident. However, S I was significantly higher in the control treatment, indicating that labile substrates represent a greater limitation to R S in the control plots. This is consistent with observations of other grassland warming experiments which showed higher labile carbon content in warmed soils due to greater belowground allocation and turnover of roots [11].
The significant reduction in soil volumetric water content as a result of soil warming has potential implications for the effects of soil warming on R S . Both R S and R H were observed to be water limited below a soil volumetric water content value of 0.2 m 3 m 23 . Thus, warming-induced soil drying may serve to mitigate warming-enhanced carbon losses to R S , as other warming experiments have shown [45]. However, the soil-drying effects observed here were small except when water was already limiting in all treatments. We suggest that soil drying effects of warming are contributing little to the mitigation of warming effects at this site, due to the frequency of rainfall and the relatively short duration of water limited periods.
The 12% increase in R S with nitrogen addition is consistent with findings for temperate grasslands [46]. In young and severely nitrogen limited ecosystems, added nitrogen may increase the amount of photosynthate allocated belowground [16]. This may be a reasonable explanation for increased R S in our experimental tussock grassland, which was planted shortly before measurements began and has very low levels of plant-available nitrogen.
The increase in R S due to nitrogen addition can be attributed entirely to R A , as R H remained unaffected by nitrogen addition.
The analysis of f RH confirmed that autotrophic contribution to R S increased with addition of nitrogen. Likewise, nitrogen increased R 10 for R S , but had no effect on R H . Similar to warming, we found no significant increase in root biomass in the treatment with added nitrogen, though there was a trend for higher root biomass. Specific root respiration has been shown to increase with increasing root tissue nitrogen concentration in grasslands [47]. Thus, it is likely that increased specific root respiration rate as a result of nitrogen addition contributed to this increase in R S .
Contrary to expectations, plant-available nitrogen in the soil was not increased by warming or nitrogen addition. There may be several factors contributing to this result. First, nitrogen was applied to both the plant and soil. As a result, a portion of the nitrogen was likely intercepted by foliar uptake [48]. Additionally, the tussock grassland soils are subject to heavy leaching, likely decreasing the residence time of added nitrogen in soils. Further, the PRS probes used to estimate plant available nitrogen were inserted into soil with roots. Strong competition for nitrogen amongst roots may have contributed to the low level of plantavailable nitrogen in all treatments.
No interactive effects of warming and nitrogen addition were observed for R S , R H or their respective temperature responses. This indicates that the effects of these separate global change drivers are additive. It has been suggested that global change drivers may interact, resulting in smaller effect sizes than those reported for single drivers [28]. However, few such instances have been noted for R S [49,50]. The only significant interaction observed in this study was the negative interaction between warming and nitrogen addition on microbial biomass. Nitrogen addition generally results in decreased microbial biomass [26]. While we found no significant decrease of microbial biomass under nitrogen addition alone, a decrease was observed under combined warming and nitrogen addition. However, this decrease did not result in reduced R H . Absent from this study is the inclusion of the rhizosphere priming effect in our estimate of R H . This refers to the effect that living roots have on R H as a result of their impact on the physical and chemical environment within the soil [51]. Priming effects can influence both the rate and temperature response of R H [52][53][54]. This may represent a potential source of error in our determination of the contribution of R H to R S . A previous study in tussock grassland soils showed a dampening of the short-term response of R H to temperature when plants were present [54]. However, in that study, priming effects were absent when plant and soil were held at a constant temperature of 15uC. Only when the soil temperature was perturbed from the constant incubation temperature over a period of hours were priming effects observed. As such, use of the root exclusion method may be appropriate for evaluating longer-term, seasonal temperature responses of R H , as in the present study.
Our results highlight the potential impacts of warming and nitrogen addition on the global carbon cycle. Over the course of the 27 month experiment, simulated cumulative CO 2 emissions, based on measured temperature response curves of R S , were 621 g C m 22 for the control treatment (Fig. 3). Warming increased cumulative R S to 953 g C m 22 , nitrogen addition resulted in cumulative emissions of 750 g C m 22 and the combined effect resulted in emissions from R S of 1127 g C m 22 . While these represent substantial differences in emissions, the contrasting responses of autotrophic and heterotrophic respiration to the treatments must be considered. While increases in R A may have consequences for the carbon economy of plants, they are likely to be offset by increased primary production. However, increases in R H due to warming present the potential for sustained loss of stored soil carbon. Extrapolation of our results to the 4.3 Mha of tussock grassland in New Zealand suggests the additional 70 g C m 22 y 21 carbon losses to R H as a result of 3 uC warming would amount to a positive feedback to rising atmospheric CO 2 concentration equivalent to 30% of New Zealand's current annual fossil fuel emissions. Table S1 F-values for fixed effects in the best-fit linear mixed-effects model of soil volumetric water content. (DOC)    Table S7 F-values for fixed effects in a nonlinear mixedeffects model of soil respiration including the interaction between the warming treatment and measurement year as a fixed effect on R 10 and E 0 parameters.

(DOC)
Appendix S1 Experimental layout for the Cass Warming Experiment. (DOC)