Herbivory modulates soil CO2 fluxes after windthrow: a case study in temperate mountain forests

Ungulate herbivory can alter functional plant communities of early-successional forest ecosystems. The consequences of such vegetation changes on soil carbon cycling are still not fully understood. Here, we used an ungulate exclusion experiment to investigate how different levels of herbivory and associated changes in vegetation succession modulate soil CO2 efflux and its heterotrophic and autotrophic sources following windthrow in temperate mountain forests. Our results indicate that only high levels of ungulate herbivory and associated vegetation shifts from tree to rather grass dominated plant communities affect soil CO2 fluxes. We did not find evidence that a moderate herbivory level and accompanied smaller shifts in the functional plant community affect soil CO2 fluxes. A greater soil CO2 efflux under the influence of high herbivory pressure was primarily attributed to accelerated heterotrophic respiration, likely due to warmer soil conditions. Moreover, autotrophic respiration from grass roots and associated microbial communities is suggested to contribute to higher soil CO2 fluxes. We conclude that intense herbivory and accompanied successional changes in the functional plant community enhance soil carbon losses following forest windthrow. This might have negative consequences for the soil carbon stocks and for the climate system.


Introduction
Temperate forests store large amounts of carbon (C) and act thereby as important sink for atmospheric CO 2 (Goodale et al. 2002;Luyssaert et al. 2010). Natural disturbances, such as from windthrows or bark beetle attacks, can turn forest net CO 2 sinks into distinct net CO 2 sources due to reduced photosynthetic uptake and continuing respiratory losses (Amiro et al. 2010;Matthews et al. 2017;Yamanoi et al. 2015). The magnitude of net CO 2 losses after disturbance depends largely on vegetation regrowth. With tree regeneration, it takes roughly one to two decades until disturbed forests become CO 2 neutral again (Amiro et al. 2010;Matthews et al. 2017;Yamanoi et al. 2015). Across many European forests, tree regeneration is, however, delayed or completely inhibited following disturbance, often as a result of ungulate herbivory (Ammer 1996;Pröll et al. 2014;Ramirez et al. 2019;Thrippleton et al. 2018). In Austria, for example, about 60% of the regenerating forests are heavily affected by browsing from deer and chamois (Schodterer 2016). Moreover, herbaceous species such as Calamagrostis sp. can represent strong competitors for seedlings and saplings and once a dense ground vegetation layer has established, sites can remain in herbaceous dominated, non-forest states for several decades post-disturbance (Kupferschmid and Bugmann 2005;Pröll et al. 2014;Rebele and Lehmann 2001;Thrippleton et al. 2018). This raises the question of whether ungulate herbivory and accompanied shifts in the functional Communicated by Christian Ammer.

Mathias Mayer and David
Keßler have contributed equally to this work.

Electronic supplementary material
The online version of this article (https ://doi.org/10.1007/s1034 2-019-01244 -9) contains supplementary material, which is available to authorized users. 1 3 plant community affect C cycle dynamics of early-successional forests, with potential consequences for the ecosystem C balance and the CO 2 exchange with the climate system. Soil CO 2 efflux (Fs; = soil respiration) represents the largest respiratory flux across forest ecosystems (~ 60% of ecosystem respiration), and its relative contribution to the forest C balance could be shown to increase following disturbance (Janssens et al. 2001;Paul-Limoges et al. 2015). Soil CO 2 efflux is strongly mediated by vegetation type and plant community characteristics (Bond-Lamberty et al. 2004b;Metcalfe et al. 2011;Raich and Tufekciogul 2000). Plants drive Fs principally via effects on belowground allocation of photosynthetically fixed C, litter quality and quantity, and microclimatic conditions (e.g. temperature and moisture) (Metcalfe et al. 2011). Moreover, plant groups can be associated with distinct microbial communities (e.g. mycorrhizal symbionts) (Brundrett 2009;Van Der Heijden et al. 2015), thereby affecting Fs differently. Successional shifts in the functional plant community composition due to herbivory are therefore assumed to affect both, the autotrophic (i.e. CO 2 from roots and associated micro-organisms) and heterotrophic (i.e. CO 2 from decomposition of soil organic matter) source of Fs. Particularly, the response of heterotrophic respiration to herbivory-affected succession might be a key factor influencing both, the net ecosystem-atmosphere CO 2 exchange and the amount of C stored in forest soils.
Several studies have investigated forest disturbance effects on Fs and its sources (Kobler et al. 2015;Kulmala et al. 2014;Mayer et al. 2014Mayer et al. , 2017bZehetgruber et al. 2017) and herbivory effects on a variety of soil characteristics and processes, including soil C storage, litter decomposition, and Fs (Andriuzzi and Wall 2017;Bardgett and Wardle 2003;Ellis and Leroux 2017;Prietzel and Ammer 2008;Tanentzap and Coomes 2012). Yet, it remains largely unknown how ungulate herbivory and accompanied shifts in successional plant community composition modulate in situ soil CO 2 fluxes after disturbance in temperate forests. Here, we evaluate how different levels of herbivory by deer and chamois affect the heterotrophic and autotrophic sources of Fs at forest windthrow sites in the Austrian Alps using ungulate exclusion (fence) and root exclusion (trenching) treatments. We hypothesized that Fs would increase with an increasing influence of herbivores and an accompanied increase in grass cover (Fig. 1a). An increase in Fs was Fig. 1 Hypothesized effect of ungulate herbivory and succession on soil CO 2 efflux and its heterotrophic and autotrophic respiratory sources (a). Control and fence treatments at windthrow sites in Höllengebirge and Reutte, respectively; a stronger herbivory pressure and accompanied stronger vegetation changes are apparent at Reutte (b) thereby primarily attributed to accelerated heterotrophic respiration rates due to more favourable microclimatic conditions (e.g. warmer) under herbs/grasses (Mayer et al. 2014(Mayer et al. , 2017a. Additionally, higher autotrophic respiration rates from a dense herbaceous vegetation layer (Mayer et al. 2014(Mayer et al. , 2017bZehetgruber et al. 2017) were assumed to contribute to higher Fs under the influence of herbivores (Fig. 1a).

Study sites
The study took place in the Höllengebirge mountain range (henceforth 'Höllengebirge'; 47° 47′ 19″ N, 13° 38° 21″ E) and in the Tannheimer mountain area near the town of Reutte (henceforth 'Reutte'; 47° 27′ 27″ N, 10° 40′ 12″ E), located in the Austrian Calcareous Alps. Höllengebirge is a south-west exposed site at an altitude of 1000 m a.s.l, and Reutte is a south-east exposed site at an altitude of 950 m a.s.l. Average air temperature and precipitation in 2016 were 7.1 °C and 1903 mm at the Höllengebirge site and 8.6 °C and 1351 mm at the Reutte site (data provided by Zentralanstalt für Meteorologie und Geodynamik-ZAMG). The Höllengebirge and Reutte sites were similar to one another regarding bedrock which was mainly limestone in paragenesis with dolomite. Chromic Cambisols, Rendzic Leptosols, and Folic Histosols (IUSS Working Group WRB 2006) were the dominant soil types and Moder and Tangel (Zanella et al. 2019) the main humus forms. Mean soil properties of the sites are given in Table 1.
The natural woodland community at both sites is dominated by Picea abies, Abies alba, and Fagus sylvatica (Kilian et al. 1994). Both sites were affected by a storm event during which several hectares of the forest stand were either blown over or destroyed by wind snap. The Höllengebirge and Reutte site was disturbed in 2007 and 2003, respectively. At both sites, the timber (mainly the stem fraction) was removed after the disturbance event. The Höllengebirge windthrow site (henceforth 'Höllengebirge control') was dominated by herbaceous ground vegetation (Calamagrostis sp., Carex alba, and Adenostyles glabra); a mixture of naturally regenerating and planted trees (Picea abies and Larix decidua) covered roughly a quarter to a third of the site. The Reutte windthrow site (henceforth 'Reutte control') was dominated by a dense layer of Calamagrostis sp. and other grass species. Based on the tree regeneration density/cover at the control plots, the herbivory effects were classified as 'moderate herbivory level' at the Höllengebirge site and 'high herbivory level' at the Reutte site (Fig. 1b).
A fence (height 2 m) covering an area of ~ 0.7 ha was installed at each site. Year of installation was 2008 and 2010 at Reutte and Höllengebirge, respectively. This ungulate exclusion treatment allowed for studying the effect of herbivory (particularly browsing) on tree regeneration and successional plant communities after windthrow. In 2016, vegetation inside the fence was dominated by a mixture of naturally regenerating and planted trees (Picea abies, Larix decidua, Abies alba, Acer pseudoplatanus, Sorbus aucuparia, Salix sp., Fagus sylvatica, and Pinus sylvestris) which covered most of the surface (Fig. 1b). Fence treatments at Höllengebirge and Reutte are henceforth denoted as 'Höllengebirge fence' and 'Reutte fence', respectively.

Field measurements
Prior to Fs measurements, plastic collars (4 cm height, 10 cm diameter, 3 cm inserted into soil) were installed across Höllengebirge control (n = 23), Höllengebirge fence (n = 16), Reutte control (n = 16), and Reutte fence (n = 16) treatments; locations were distributed randomly. Root-exclusion plots (Subke et al. 2006) were established at each site by digging trenches down to either bedrock or a maximum depth of 80 cm, each encompassing an area of 1 × 1 m. All roots within the trenches were cut, and a plastic sheet was installed to inhibit root and mycorrhizal in-growth. For Fs measurements, one plastic collar (see above) was installed in the centre of each root-exclusion plot. In total, 14 and 10 root-exclusion plots were installed across the Höllengebirge and Reutte site, respectively (distributed equally between control and fence treatments). Ground vegetation within root-exclusion plots and collars was removed prior to measurements, and regrowth was cut regularly (Bahn et al. 2008).
During the vegetation period of 2016, Fs was measured in monthly intervals by means of the closed chamber technique, using a portable infrared gas analyser (EGM-4; PP Systems International Inc., Amesbury, MA, USA). Measurements were taken by connecting a respiration chamber (SRC-1; PP Systems International Inc.) to the plastic collars. Soil temperature at 5 cm depth (handheld thermometer) and soil moisture between 0 and 7 cm depth (Field Scout TDR Soil Moisture Meter; Spectrum Technologies Inc., Aurora, IL, USA) were measured next to the collars at the same time. To avoid a temporal sampling bias, sequential arrangement of treatments and plots was alternated between measurements. At each site, soil temperature and moisture were also monitored continuously (1 h interval) by means of five GS3 sensors (installation depth 5 cm) and EM50 data loggers (Decagon Devices, USA). To get spatially representative hourly temperature and moisture data for each site and treatment, continuous measurements were corrected by linear regression to the manually gathered soil temperature and moisture measurements during the corresponding time (Wangdi et al. 2017). At Reutte, tree and grass cover was estimated for each plot using a 1 × 1 m frame; collars represented the centre.

Statistical analysis
The effect of herbivory exclusion on Fs, soil temperature, and soil moisture were tested by means of analysis of variance with a mixed-effects model structure at each site (Pinheiro et al. 2014). To account for the repeated measurement structure within the data, the plots were assigned random effects and the treatments were assigned fixed effects in each model. Linear mixed-effects models were used to predict Fs (µmol CO 2 m − 2 sec − 1 ) as a function of soil temperature (°C) and moisture (vol%) (Reichstein and Beer 2008): Due to an exponential relationship between Fs and soil temperature, data were log-transformed prior to analysis. Fixed effects model coefficients β 0 − β 2 were estimated from the data and are given in Table 2. The terms Fs ij , T ij , M ij are the measured soil CO 2 efflux, soil temperature, and soil moisture, respectively (for observation j at plot i). Plots were assigned random effects where terms a i and b i depict random intercepts and slope coefficients for the response to soil temperature, respectively. The term ε ij is the residual error. Including a random slope to describe a plot-wise response of Fs to soil moisture and an interaction of temperature and moisture did not improve the model fit. Marginal and conditional R 2 served as indicators for model goodness of fit (Nakagawa and Schielzeth 2013). Additionally, the relationship between Fs and temperature only was tested (Table S1).
The linear mixed-effects model (Table 2) and the average continuous soil climate data of the control and fence treatments were used to predict log(Fs) for each plot (also for root-exclusion plots) and each day between 7 June until 7 November. Soil CO 2 efflux from root-exclusion plots was assumed to be from heterotrophic respiration only. Autotrophic respiration rates were calculated by subtracting the predicted daily mean heterotrophic respiration rates from log(Fs) predictions of non-root exclusion plots. The log(Fs) predictions were transformed to Fs by taking the exponent.
Structural equation modelling (Beaujean 2014;Grace 2006) was used to explore the effects of grass and tree cover on the sums of Fs. This analysis was conducted for Reutte only, as no vegetation data were available for Höllengebirge plots. An initial a priori model included pathways between Fs, tree and grass cover and a correlation between tree and grass cover. Two alternative final models were obtained by removing pathways in a stepwise procedure. Only significant pathways were kept in the final models. Model selection was based on Akaike's information criterion (AIC), chi-squared test results, and comparative fit index (CFI) (Beaujean 2014;Grace 2006).

Results and discussion
In this study, we investigated how ungulate herbivory and associated differences in successional plant communities affect soil CO 2 efflux (Fs) and its heterotrophic and autotrophic sources at two windthrow sites in the Austrian Alps. Herbivory effects were assessed by means of fence treatments, and the sources of Fs were separated by means of root exclusion treatments. Fenced treatments were characterized by a dense tree regeneration. Herbivory promoted the growth of herbal and grass species. A moderate herbivory level at Höllengebirge control plots resulted in a mixture of regenerating trees and herbaceous plants while a high herbivory level at Reutte control plots resulted in a dense layer of Calamagrostis grasses and other graminoids, respectively (Fig. 1b). This successional pattern is in line with Tremblay et al. (2006), who showed that the abundance of browsetolerant species, such as grasses, was positively related to deer density in harvested forest stands. Soil CO 2 efflux was significantly higher at herbivoryaffected control plots at Reutte when compared to fence plots, while no fence effect was present at Höllengebirge (Fig. 2a, b). Average sums of Fs of 457 ± 29, 449 ± 36, 860 ± 67, and 594 ± 40 g C m −2 (mean ± SE for a 5 months' period) were calculated for Höllengebirge control, Höllengebirge fence, Reutte control, and Reutte fence, respectively (Fig. 3). The control plots at Reutte respired on average 45% more CO 2 than the fence plots. Our hypothesis, Fs would increase with increasing herbivory pressure, was therefore only supported in the case of a high herbivory level and Fig. 2 Soil CO 2 efflux, heterotrophic respiration (= soil CO 2 efflux from root exclusion plots), soil temperature, and soil moisture at windthrow sites in Höllengebirge and Reutte (mean ± SE). Measurements were conducted in control and fence treatments. Given are p-values of ANOVA with mixed effects model structure       Fig. 3 Modelled, cumulative soil CO 2 efflux sums (7 June until 7 November 2016) from control and fence treatments at Höllengebirge and Reutte. Sums are separated into heterotrophic and autotrophic respiration. Autotrophic respiration was calculated as the difference between soil CO 2 efflux and heterotrophic respiration. Error bars represent standard error of the mean accompanied large plant community differences inside and outside the fence (i.e. tree versus grass dominated plots).
Conversely, it appears that a moderate herbivory level and smaller plant community differences inside and outside the fence (i.e. tree versus tree/grass/herb dominated plots) did not affect Fs rates. Generally, these findings are partly in line with the results of a recent meta-analysis by Andriuzzi and Wall (2017), who reported an overall positive response of Fs to herbivory across temperate ecosystems. However, their meta-analysis was biased towards temperate grasslands and herbivory by sheep and included also laboratory measurements of heterotrophic respiration. Support that Fs would increase with decreasing tree and increasing grass abundance is also given by structural equation modelling. The best model reveals a negative relation between Fs and tree cover (Fig. 4a). A clear negative correlation between grass and tree cover masked a positive relation between Fs and grass cover. This relationship was significant after removing the path between tree cover and Fs from the model (alternative model, Fig. 4b).
Heterotrophic respiration represented the dominant CO 2 flux at the studied sites, with an average contribution to Fs ranging between 56 and 77% (Fig. 3). These values are in good agreement with previous findings from a fire-chronosequence study in boreal forests (Bond-Lamberty et al. 2004a) and from studies on clear-cut sites in temperate forests (Williams et al. 2014;Zehetgruber et al. 2017). Similar to Fs, heterotrophic respiration was only affected by high herbivory pressure at Reutte (Fig. 2c, d). There, herbivory resulted in 33% higher CO 2 fluxes from heterotrophic respiration as suggested by our modelling results (Fig. 3). We therefore argue that greater Fs rates under high herbivory levels can be primarily explained by enhanced heterotrophic respiration from the decomposition of soil organic matter.
Soil CO 2 efflux and heterotrophic respiration showed distinct seasonal variations in all treatments, strongly following the patterns in soil climate (Fig. 2). Soil temperature and moisture explained 34% and 69% of the overall and plot specific variation in Fs (Table 2). However, soil temperature only explained 32% of the overall and 64% of the plot specific variation in Fs (Table S1). As temperature was the major abiotic driver of Fs, the greater heterotrophic respiration rates at Reutte control plots can, to a large degree, be explained by warmer soil conditions (Kulmala et al. 2014;Mayer et al. 2014Mayer et al. , 2017aZehetgruber et al. 2017); soil temperature was on average 1.9 ± 0.5 °C higher at the grass-dominated control plots when compared to the fence plots (Fig. 2f). In contrast, no differences in soil temperature could be seen at the Höllengebirge site where regenerating trees were present at both, fence and control plots (Fig. 2e). Thus, higher temperatures at Reutte control plots are mainly caused by a lack in crown shading from regenerating trees. These results correspond well to Mayer et al. (2017a), where the absence of regenerating trees increased soil temperature and decomposition processes following forest gap disturbance, while their presence kept temperature and decomposition even close to control stand levels. Additionally, higher heterotrophic respiration rates at Reutte control plots might be related to changes in the litter provided to the microbial community (Bardgett and Wardle 2003). As shown for clear-cut sites, particularly grasses can produce leaf biomass stocks (and thus litter) which exceed those of regenerating trees by far (Tremblay et al. 2006). Herbivores may have also directly affected microbial activity by the deposition of faeces or urine (Bardgett and Wardle 2003).
The higher Fs under the influence of intense herbivory is also likely due to a greater contribution of autotrophic respiration from roots and associated micro-organisms. Estimates for autotrophic respiration suggest ~ 30% higher plant-associated soil CO 2 fluxes at Reutte (Fig. 3). Particularly, Calamagrostis species develop large root biomass pools which can be in a similar size of the fine root pools of mature forest stands (Brunner et al. 2013;Rebele and Lehmann 2001). The supply of fresh C from roots to the microbial community might have additionally enhanced the decomposition of older soil organic matter, a process which is known as rhizosphere priming (Gavazov et al. 2018;Kuzyakov 2010). Herbaceous species are, moreover, Fig. 4 Structural equation models describing the influence of tree (a) and grass (b) cover (%) on soil CO 2 efflux. Single-and double-headed arrows represent significant relationships and correlations between variables, respectively. Values represent model coefficients and their standard errors (in brackets). Comparative fit index (CFI) is a measure for the goodness of the model fit associated with different root symbionts when compared to trees. While most tree species form symbioses with ectomycorrhizal fungi, grasses and herbs form symbioses primarily with arbuscular mycorrhizal fungi (Brundrett 2009;Van Der Heijden et al. 2015). Respiration from mycorrhizal fungi can make a great contribution to Fs (Heinemeyer et al. 2007;Subke et al. 2011). Differences in mycorrhizal symbionts are therefore likely to additionally modulate autotrophic soil CO 2 fluxes under an herbivory associated shift in the plant community composition.
It should be stressed that our results are based on pseudoreplicated measurements (Hurlbert 1984) from only one site per herbivory stage. Thus, we cannot claim for an ecological generalization. The results are representative for an intermediate stage after forest disturbance only (i.e. ~ one decade after windthrow). Moreover, the pre-fence conditions are unknown and the sites differed with regard to time since disturbance and soil properties (e.g. organic layer thickness). A thicker organic layer, for example, may be a cause for the generally higher Fs rates at Reutte (Fig. 3). It is also likely that the Reutte site is in a different phase of soil organic matter decomposition due to a later successional stage after windthrow. A conclusion regarding temporal changes of post-disturbance development of Fs is therefore not possible. A conclusion regarding herbivory levels, related plantcommunity changes, and how this affects soil CO 2 fluxes is, nonetheless, possible.
Collectively, our results indicate that only high levels of ungulate herbivory and accompanied shifts from tree to rather grass dominated plant communities affect Fs and its heterotrophic and autotrophic sources at the studied windthrow sites. We did not find evidence that a moderate herbivory level and accompanied smaller shifts in the functional plant community affect soil CO 2 fluxes. Higher Fs rates under the influence of intense herbivory were primarily attributed to accelerated heterotrophic respiration, likely due to warmer soil conditions. Moreover, autotrophic respiration from grass roots and associated microbial communities might have additionally promoted higher Fs rates. Disturbances, such as windthrows, can cause significant C losses from the forest soil and can turn former forest CO 2 sinks into distinct CO 2 sources to the atmosphere (Amiro et al. 2010;Thom and Seidl 2015). Greater heterotrophic respiration rates due to intense ungulate herbivory may further enhance net soil C losses following disturbance, with consequences for both the soil C storage and the climate system. On the other hand, higher autotrophic respiration rates from grasses might indicate higher photosynthetic fixation rates and thus greater input rates of fresh C to soil (Litton et al. 2007). This could mitigate post-disturbance net C losses at least to some degree (Zehetgruber et al. 2017). Additionally, grass roots have a higher root turnover than tree roots suggesting higher belowground litter input rates when compared to trees (Solly et al. 2014). Further studies on how ungulate herbivory and associated changes in vegetation affect soil C cycling are, therefore, urgently required to better understand plant-soil-atmosphere interactions following forest disturbance.