- Research Paper
- Open access
- Published:
Elevation affects both the occurrence of ungulate browsing and its effect on tree seedling growth for four major tree species in European mountain forests
Annals of Forest Science volume 81, Article number: 13 (2024)
Abstract
Key message
In European mountain forests, the growth of silver fir (Abies alba Mill.), sycamore maple (Acer pseudoplatanus L.), European beech (Fagus sylvatica L.) and Norway spruce (Picea abies (L.) H. Karst.) seedlings is more strongly affected by ungulate browsing than by elevation. But, the constraint exerted by ungulates, in particular the probability for seedlings to be browsed, increases with elevation for most species.
Context
While concerns about mountain forest regeneration rise due to their high vulnerability to climate change, the increase in wild ungulate populations and the expansion of their range in the last decades exert an additional constraint on the survival and growth of young trees. Understanding how this constraint can vary with elevation is thus a key to assess the consequences of this population increase for the regeneration of mountain forests.
Aims
In this study, we investigate the effect of elevation on (i) the occurrence of browsing for seedlings and on (ii) the reduction in seedling growth induced by ungulate browsing.
Methods
We monitored height growth and browsing occurrence on silver fir, sycamore maple, European beech and Norway spruce seedlings across seven elevation gradients (from 400 to 2013 m) located from France to northern Sweden.
Results
Seedlings of the two most palatable species—fir and maple—were more likely to be browsed at high elevation while the opposite effect was observed for spruce. Browsing strongly reduced seedling growth for all species but Norway spruce, while elevation had no direct effect on seedling growth. This browsing-induced growth reduction was stronger at high elevation for fir seedlings.
Conclusions
Browsing is overall a stronger constraint on seedling growth than elevation for four dominant species of European mountain forests. Elevation can, however, affect both browsing probability and the effect of browsing on seedling growth. Our results highlight the importance of taking into account ungulate pressure and its interactive effect with elevation when forecasting the regeneration of mountain forests under a changing climate.
1 Introduction
Mountain forests are keystone ecosystems known to provide a wide range of services to human society, e.g. protection against natural hazards, carbon storage or timber production (Schirpke et al. 2019; Stritih et al. 2021). In addition to being particularly vulnerable to climate change (Thuiller et al. 2005), mountain forests are increasingly threatened by the current dynamics of ungulate populations. Indeed, the past decades have witnessed an important increase in wild ungulate populations in the Northern Hemisphere, in terms of both abundance and geographical range (Fuller and Gill 2001; Milner et al. 2006; Beguin et al. 2016). As a consequence, ungulate damage is a growing concern for agricultural production, forest ecosystems sustainability and conservation areas (Reimoser and Putman 2011). Ungulates are key drivers of forest regeneration leading to reduced seedling growth and survival and to changes in plant assemblages benefiting species that are less palatable or more tolerant to browsing (Boulanger et al. 2015; Laurent et al. 2017). In mountain forests, elevation can also be a driver of trees’ demographic performance, mainly through a sharp decrease in temperature with increasing elevation which reduces growth (Marcora et al. 2013; Jiang et al. 2022; Muñoz Mazón et al. 2022), but also through changes in e.g. direct solar radiations or snow pack thickness in winter (Körner 2007, 2012). Many of the environmental conditions that change along elevation gradients can also affect the browsing pressure exerted on seedlings, e.g. snow cover (Hallinger et al. 2010) and tree species composition (Champagne et al. 2016), or seedlings’ response to browsing, e.g. climate (Barton and Shiels 2020; Champagne et al. 2021). The resultant effect of elevation on the constraint exerted by deer on forest regeneration remains, however, largely unknown. As high-elevation regions are expected to experience amplified warming (Pepin et al. 2015), understanding how elevation and deer browsing interactively affect natural forest regeneration is crucial to anticipate the consequences of global changes for mountain forests.
The effect of elevation on browsing was a matter of conflicting results (Speed et al. 2013; Ameztegui and Coll 2015; Herrero et al. 2016). First, because several species of ungulates only use woody species as alternate forage during winter when their main food resources (e.g. herbs and shrubs) become scarce (Heggberget et al. 2002), browsing rates on tree seedlings could increase at higher elevation, where the winter season is harsher. This would be consistent with the attractant-decoy hypothesis, which predicts that plants avoid herbivory when surrounded by more palatable neighbours that divert browsing (Ruttan and Lortie 2013). The opposite effect could, however, also be expected: theories predicting food selection at broad spatial scales suggest that ungulates are rather expected to primarily browse vegetation in feeding patches with abundant and highly palatable forage (Bee et al. 2009; Kuijper et al. 2009), which should mostly be located at low elevation, where productivity is higher. A typical illustration is the observation that several ungulate species such as red deer (Cervus elaphus Linnaeus) (Kropil et al. 2015) or roe deer (Capreolus capreolus L.) (Cagnacci et al. 2011) can migrate downslope in winter, where conditions are milder. In addition, more frequent browsing occurrences on seedlings could also be expected at low elevations due to a thinner snow pack, since snow cover has a protective effect against browsing, particularly on small seedlings (Hallinger et al. 2010). Overall, it remains thus largely disputed how elevation might affect the browsing pressure exerted by ungulates on tree seedlings.
Environmental variations along elevation gradients also exert strong effects on seedling performance, with potentially complex interactions with browsing. For instance, theoretical studies showed that herbivory should have a stronger impact on plant performance in conditions of low resource availability (Wise and Abrahamson 2005). In mountain ecosystems, resource availability can be strongly conditioned by elevation (i.e. higher temperatures leading to reduced water availability at low elevation, higher cloud cover reducing light availability at high elevation) (Körner 2007, 2012). In practice, however, studies that investigated the interactions between climate and browsing on seedling growth produced contrasted results. At the tree line, herbivory by large mammals can either decrease (Brodie et al. 2012; Vuorinen et al. 2020) or increase climate effect on seedling growth (Speed et al. 2011; Sanders-DeMott et al. 2018). Similarly, in controlled conditions, reduced water availability can have a positive (Hawkes and Sullivan 2001) or neutral (Kullberg and Welander 2003; Barton and Shiels 2020) effect on seedling tolerance to browsing.
While it is crucial to keep in mind that environmental changes along elevation gradients are not perfectly analogous to those expected to occur in time with climate change, elevation-driven temperature gradients remain a powerful research tool to infer the ecological response of organisms to climate change (Körner 2007; Sundqvist et al. 2013). Investigating the relative effect of elevation and browsing on seedling growth could provide insights into the relative importance of changes in temperature and ungulate populations for forest regeneration. Several studies have indeed suggested that browsing might have stronger impact than climate in controlling tree regeneration (Cairns and Moen 2004; Speed et al. 2010; Herrero et al. 2012, 2016), which could have strong repercussions for future changes in tree species distribution. For instance, Fisichelli et al. (2012) showed that across North America, high browsing pressure could annihilate the positive effect of temperature on the growth of temperate tree species, thus pushing the temperate–boreal transition zone further south.
In this study, we analysed the joint response of tree seedling performance to elevation and browsing in tree species of mountainous forests throughout Europe. We focused on silver fir (Abies alba Mill.), sycamore maple (Acer pseudoplatanus L.), European beech (Fagus sylvatica L.) and Norway spruce (Picea abies (L.) H. Karst.), which are four of the most dominant and economically important tree species in European mountains. We analysed seedling growth and browsing damage along seven elevational gradients, spread across a European latitudinal gradient reaching from France to northern Sweden, through Germany and Poland. Based on these data, we explored the following questions:
-
1.
How does elevation affect the occurrence of browsing for the seedlings of these four species?
-
2.
Are there interactive effects of elevation and browsing on tree seedling growth?
-
3.
Is seedling growth more affected by elevation or by browsing?
2 Materials and methods
2.1 Study sites, species and sampling
Seven elevational gradients (sites) were sampled, in four European countries (France, Germany, Poland and Sweden) covering a large latitudinal (from 44° to 67°) and elevational (from 446 to 2013 m above sea level) range (Fig. 1 and Table 1). Sampled sites were all located in mountain forests and were overall representative of the local forest structure (i.e. even or uneven-aged low intensity management in all sites except in the Bavarian Forest National Park (FNP) which is a protected area with no management). Free-ranging populations of ungulates occurred in each site (see species list in Table 1). All sites were covered by a 30-cm-deep snow cover or more during winter (Appendix 1). The dominant tree species, the ungulate species present, local climate and site characteristics are listed in Table 1 and in Rabasa et al. (2013). Extraction of climatic data for each site is provided in Appendix 2.
We focused on silver fir (Abies alba Mill.), Norway spruce (Picea abies (L.) H. Karst.), European beech (Fagus sylvatica L.) and sycamore maple (Acer pseudoplatanus L.) (Table 1) that were the four most abundant tree species across all sites and in each of the plots selected within sites. They display different ecological strategies: sycamore maple and Norway spruce are considered post-pioneer species, while silver fir and European beech are late-successional species (Rameau et al. 1993).
On each site, we sampled seedlings at different elevation levels, which were separated by at least 100 m in elevation, except for Swedish transects (i.e. 50 m), where the temperature decreases more rapidly with increasing elevation (Stone and Carlson 1979). The total number of elevation levels sampled varied from 2 to 7 depending on the length of the elevation gradient. In all sites but Bavarian Forest National Park, we established between four to six circular plots per elevation level (plots were 20 m in radius; see details in Rabasa et al. (2013)). In the Bavarian FNP, the plots were located every 100 m along four transects, for a total of 30 plots. To minimize the effect of forest management, we selected within each site plots with no evidence trees having been harvested over the past 4 years (defined as no observation of recent stumps in the vicinity of the plot).
In each plot, seedlings (defined in this study as young trees with a height < 70 cm) were sampled in square-shaped subplots located in the centre of each plot. The density of small seedlings can be extremely variable at small spatial scale (depending on the distance to a seeding tree and local light availability). Therefore, we used a sampling protocol that allowed adjusting to this natural variability. More specifically, the size of subplots varied among species depending on the local seedling density: 4 m2 when seedling density was sufficient to sample 5 individuals per species, 16 m2 in plots of overall low seedling abundance. We sampled a minimum of 5 seedlings and up to 20 seedlings per species per plot. When more than 20 seedlings per species were present in a plot, we randomly sampled seedlings with different heights. We verified a posteriori that height distribution was consistent whether there was a random selection of seedlings (i.e. in plots with more than 20 seedlings) or not (see Appendix 3). In rare cases, it was not possible to find 5 seedlings of each species in the plots so we randomly sampled up to 5 seedlings (< 70 cm) of each dominant species at each elevation searching up to 100 m away from the plot (while avoiding the 100-m search radius of nearby plots). The height of each seedling was measured, and browsing marks were recorded (1 when the terminal shoot was browsed during the current year, 0 when not). The annual height growth of each seedling was measured as the length of the last shoot. Seedlings were sampled at the end of the growing season in each site and elevation to ensure that fully expanded shoots were measured. Measuring the last shoot length is a common proxy for seedling growth (e.g. see Canham et al. (1994) and Palmer and Truscott (2003)). Using such a metric does not capture overcompensation by new shoots replacing the apical one, but we assumed this pattern to be negligible as polycyclic growth is rare in mountain forests, especially for small seedlings. Seedlings were sampled in 2010, except in 2017 for the Bavarian Forest National Park site. In total, 2981 seedlings were measured (Bernard et al 2023).
2.2 Statistical analysis
2.2.1 Effect of climate on browsing occurrence
We tested the effect of elevation on browsing occurrence (binary variable (Br)) separately for each tree species. In addition to the effect of elevation (Z), we chose to take into account seedling height (Ht) because taller seedlings may be more apparent than smaller ones (Diaci 2002; Nichols et al. 2015). We included an interaction between elevation and seedling height because taller seedlings might be less protected by the snow pack (Visscher et al. 2006). We used a generalized mixed linear model for browsing probability (p Br) (model 1) with a binomial error term and a logit link function
Elevation and seedling height (Ht) were centred and scaled. Elevation, seedling height and their interaction were treated as fixed effects. Elevation level (ZL) and sites (S) were treated as random effects, with elevation level nested in the site (θ S/ZL). We chose to treat the site as a random effect to avoid biases related to between-sites variations in unobserved environmental conditions (e.g. in ungulate density and thus ungulate pressure). Elevation level was also treated as a random effect to account for the fact that plots located at the same elevation level are spatially clustered. α is the estimated parameter, and ε is the error term.
2.2.2 Effect of climate and browsing on seedling growth
Second, to test the relative effects of elevation and browsing (Br) on seedling growth, i.e. as measured by the last shoot length (LSL), we used a linear mixed model (2) with here again seedling height (Ht) as a covariate to account for the fact that taller seedlings tend to grow faster.
The response LSL was log-transformed to match the normality assumption. Elevation (Z), browsing (Br) and seedling height (Ht) were centred and scaled. We included an interaction between elevation and browsing to test if elevation had an influence on seedling growth response to browsing. We added the same random effects as described before (θ S/ZL), i.e. elevation level nested in the site, and ε is the residual normal error.
2.2.3 Model validation
For each model, we graphically checked that the structure of the residuals matched the assumptions of normality and homoscedasticity. Because Norway spruce was the only tree species present in the coldest site (Gällivare) which also had a very different climate from the other sites (mean annual temperature = − 8.7 °C vs. mean annual temperatures comprised between + 0.3 and + 3.2 °C in the other sites), we first ran the models without this site. For Norway spruce, we repeated each model after adding data from Gällivare to check whether the trends were driven only by this particular site. Results related to these analyses are reported in Appendix 4 (Table 5 and Fig. 6) for model 1 and Appendix 5 (Table 6 and Fig. 7) for model 2.
Mixed models were fitted using glmer and lmer functions from lme4 package, in R version 3.4.1 (R Core Team 2019), and model statistics were obtained using lmerTest package (Bates et al. 2014; Kuznetsova et al. 2017).
3 Results
The residuals of the models are analysed in Appendix 6 and show that hypotheses of the model are respected.
3.1 Effects of elevation on browsing probability
Browsing probability (p Br) increased significantly with seedling height for all species, whereas elevation had more contrasted effects (Fig. 2 and Table 2). Maple seedlings were more likely to be browsed at high elevation. In contrast, elevation had a negative effect on the browsing probability in Norway spruce (Fig. 2 and Table 2), and this effect was also observed when including the Swedish site Gällivare (Fig. 8 in Appendix 6). For fir seedlings, there was a significant interaction between height and elevation (Table 2), so that for the tallest seedlings only, browsing probability was higher at high elevation. The probability for beech seedlings to be browsed was unaffected by elevation. When significant, the effect of elevation on browsing probability was overall of similar magnitude than the effect of height (Fig. 2b).
3.2 Relative effects of climate and browsing on seedling growth
Browsing reduced seedling growth (as measured by the last shoot length) for all species but Norway spruce, with a particularly strong effect in fir and beech (Fig. 3 and Table 3). Unbrowsed seedlings had an average growth rate of 2 cm year−1 and 3 cm year−1 (for fir and maple, respectively), while browsed seedlings grew by only 0.5 cm year−1 (Fig. 3).
Elevation had an overall very weak effect on seedling growth. For silver fir, there was a significant interaction between elevation and browsing (Table 3) so that browsing only reduced seedling growth at high elevation (Fig. 3a). The interaction between elevation and browsing was also significant for beech and Norway spruce (Table 3), but its magnitude was too weak to affect the hierarchy between browsed and unbrowsed seedlings. Besides, the interaction was no longer significant for Norway spruce when the Swedish site Gällivare was included (Appendix 5). Lastly, height growth increased with seedling height for all species, although the trend was not significant for silver fir (Fig. 3c), which was the species with the smallest range of seedling height (75% of the seedlings were less than 8.0 cm tall).
4 Discussion
Based on seven elevation gradients from Central to Northern Europe, we show that elevation can affect both the probability of being browsed and the browsing-induced growth reduction of seedlings for four emblematic tree species of European mountain forests. We also show that seedling height growth is more strongly affected by ungulate browsing than by environmental changes along elevation gradients for the four studied species. Given the tight correlation between elevation and temperature (Körner 2007), this suggests that for these four species, changes in ungulate populations could exert a constraint on the regeneration that is comparable in magnitude to the expected increase in temperature under future climatic scenarios (without considering other components of climate change such as drought or other disturbances).
4.1 Palatable seedling species are more frequently browsed at high elevation
We found that for the two most palatable species (i.e. silver fir and maple), seedlings located at higher elevation were more likely to be browsed by ungulates. Our finding is consistent with the idea that higher availability of alternative forage (i.e. forbs and herbaceous species) at low elevation, where productivity is higher and winter conditions are milder, leads to a lower browsing intensity on woody species (Gordon and Prins 2008). Indeed, browsing on woody species often occurs because of a lack of more palatable resource during low food-availability periods (Redjadj et al. 2014), except for roe deer (Storms et al. 2008; Barančeková et al. 2010) and moose (Felton et al. 2021) that feed mostly on woody plants. As forbs and herbaceous species are more abundant and less protected by snow as elevation decreases, our result would be consistent with a shift of food choice towards plants with a higher and more easily accessible nutritional content (Ruttan and Lortie 2013). The adverse effect could also have been expected, with ungulates exerting a higher browsing pressure in feeding patches where forage availability is higher (Bee et al. 2009; Kuijper et al. 2009), especially since species such as red deer (Kropil et al. 2015) or roe deer (Cagnacci et al. 2011) have been reported to migrate downslope in winter, where conditions are more favourable. In that regard, the results from our correlative study tend to support the idea that food selection at local scale is more important in driving the browsing probability on a given species than food selection at broad spatial scale (Champagne et al. 2016). To draw more solid conclusions on this question, the development of methodological approaches to quantify ungulate density that are both easy to implement at continental scale and comparable between and within sites would be an important step.
For the four species, taller seedlings displayed a higher probability of being browsed, which is in agreement with the observation that ungulates barely eat close to the ground (Renaud et al. 2003) and rather browse taller seedlings that are more detectable and less energetically costly to search (Bergquist and Örlander 1998; Bergström and Bergqvist 1999; Kullberg and Bergström 2001; Diaci 2002). In general, this size preference was consistent across the full climatic gradient. For silver fir, however, we observed that height only increased the browsing probability at high elevation, which is likely related to the protective effect of snow (Visscher et al. 2006). Indeed, in most of our study sites, minimum snow cover at high elevation during winter was 20 cm (Appendix 1), so that seedlings below 20 cm are always protected by the snow pack while the protection of higher seedlings depends on snow depth. Our results thus support the idea that the presence of a thick snow cover at high elevation should further amplify the already existing effect of seedling height on the probability to be browsed. The presence of a snow cover in winter could also explain our observation that Norway spruce seedlings are preferentially browsed at low elevation, unlike the other species. Indeed, unlike maple or fir that can also be browsed in spring and summer, Norway spruce is exclusively browsed in winter (Storms et al. 2008). As we focused our sampling on seedlings below 70 cm (most being smaller than 40 cm), spruce seedlings of this size were likely only accessible to ungulates at low elevation in winter. Taken together, our results show that the browsing pressure exerted on young seedlings in mountain forest is highly driven by elevation. This correlation could arise from both changes in the availability of alternative forage and protection by the snow pack.
4.2 Stronger effects of browsing than elevation on seedling growth
For all species, we found that browsing strongly impacted seedlings’ growth, except for Norway spruce that is the least palatable species. This is consistent with the findings of Motta (2003) on rowan (Sorbus aucuparia L.) and Norway spruce: in highly browsed sites, there was no more rowan growing beyond a browsing height (100–160 cm), whereas Norway spruce was barely affected by browsing. The same pattern was found in fenced designs for silver fir and Norway spruce, with very few silver fir seedlings exceeding 20 cm in height in unfenced plots while Norway spruce growth only slightly differed between fenced and unfenced plots (Bernard et al. 2017). In lowlands, the palatability of European beech is equivalently low (Boulanger et al. 2009) or even lower (Mann 2009) than that of spruce, but browsing on beech may occur in mountainous areas when there is less alternate food resource (Heinrichs et al. 2012; Diaci et al. 2012).
In contrast, elevation had no effect on seedling growth for the four studied species, while faster growth rates could have been expected at low elevation where conditions are warmer, implying longer growing season that should promote seedling growth (Fisichelli et al. 2014). Using an elevation gradient, Callaway (1998) showed that higher productivity at low elevation could be offset by stronger competitive interactions, which would be consistent with our results. Yet, the most striking result is that the effect of elevation on height growth was negligible compared to the effect of browsing for all the studied species. This stronger browsing effect is consistent with studies reporting the ability of ungulate herbivory to affect the climatic limit of the tree line (Cairns and Moen 2004; Speed et al. 2010; Herrero et al. 2012, 2016) or hinder growth responses to changes in temperature (Post and Pedersen 2008; Fisichelli et al. 2012). While it remains important to keep in mind that the changes in environmental conditions along elevational gradients are not perfectly analogous to those expected in the context of climate change, they remain strongly temperature-driven and can thus be informative on the ecological response of organisms to climate change (Körner 2007; Sundqvist et al. 2013). Considering the important magnitude of the temperature gradient within our study sites (i.e. 3.3 °C in average; see Appendix 2), our results suggest that for the four tree species we studied, the general trend of increasing wild ungulate populations could represent a comparable threat to their natural regeneration than temperature increase expected under global warming. It is thus crucial that future studies analyse jointly the effect of all components of climate change (e.g. CO2 concentration and disturbance regimes) and ungulate browsing on forest regeneration to draw clear conclusion.
4.3 Elevation affects the growth reduction induced by browsing
We found weak and often significant interactive effects of elevation and browsing on seedling growth. In particular, we found a quite strong interaction for silver fir seedlings, whose height growth was more strongly reduced by browsing at high elevation. This most likely relates to the colder conditions, shorter growing season and reduced direct irradiance in the PAR at high elevation (Körner 2012), which leaves less time and resources for seedlings to compensate for the loss of photosynthetic tissue due to browsing (Wise and Abrahamson 2005). Previous studies had shown that lower availability of resources such as light (Baraza et al. 2010) or water (Hawkes and Sullivan 2001) can modulate the browsing-induced growth reduction. But, our study is, to our knowledge, the first to demonstrate that changing environmental conditions along elevation gradients can affect seedlings’ response to browsing. When put in perspectives with our result on browsing probability, it implies that silver fir seedlings are both more frequently browsed and more severely impacted by browsing at high elevation. As such, our results shed light on the mechanisms explaining how ungulate browsing can affect the limit of the tree line for some tree species (Cairns and Moen 2004; Speed et al. 2010; Herrero et al. 2012, 2016) and advocates for more studies on the interactive effects of climate and herbivory on forest regeneration (Champagne et al. 2021; Segar et al. 2022).
5 Conclusion
Our results reveal that for four major tree species in European mountain forests, the first stage of tree regeneration is overall more strongly impacted by browsing pressure than by environmental changes observed along elevation gradients. This suggests that the general trend of increasing native ungulate populations will likely exert a major constraint on the regeneration of European mountain forests. It would thus be important to analyse how these effects on the first regeneration stage, which is a key bottleneck in tree population dynamics (Clark et al. 1999), translate into impacts at ecosystem level. Furthermore, we show that for several species, elevation influences both the intensity of the browsing pressure and the impact of browsing on seedling growth. Therefore, forecasting future dynamics of European mountain forests should be based on models, considering the critical impacts of increasing density of ungulates, and their interactions with elevation.
Availability of data and materials
Data are archived in the following Zenodo repository: https://zenodo.org/records/10370097.
References
Ameztegui A, Coll L (2015) Herbivory and seedling establishment in Pyrenean forests: influence of micro- and meso-habitat factors on browsing pressure. For Ecol Manage 342:103–111. https://doi.org/10.1016/j.foreco.2015.01.021
Barančeková M, Krojerová-Prokešová J, Šustr P, Heurich M (2010) Annual changes in roe deer (Capreolus capreolus L.) diet in the Bohemian Forest, Czech Republic/Germany. Eur J Wildl Res 56:327–333. https://doi.org/10.1007/s10344-009-0321-0
Baraza E, Zamora R, Hódar JA (2010) Species-specific responses of tree saplings to herbivory in contrasting light environments: an experimental approach. Écoscience 17:156–165. https://doi.org/10.2980/17-2-3286
Barton KE, Shiels AB (2020) Additive and non-additive responses of seedlings to simulated herbivory and drought. Biotropica 52:1217–1228. https://doi.org/10.1111/btp.12829
Bates D, Mächler M, Bolker BM, Walker SC (2014) Fitting linear mixed-effects models using lme4. J Stat Softw 67:1–51
Bee JN, Tanentzap AJ, Lee WG et al (2009) The benefits of being in a bad neighbourhood: plant community composition influences red deer foraging decisions. Oikos 118:18–24. https://doi.org/10.1111/j.1600-0706.2008.16756.x
Beguin J, Tremblay J, Thiffault N et al (2016) Management of forest regeneration in boreal and temperate deer–forest systems: challenges, guidelines, and research gaps. Ecosphere 7:1–16
Bergquist J, Örlander G (1998) Browsing damage by roe deer on Norway spruce seedlings planted on clearcuts of different ages: 2. effect of seedling vigour. For Ecol Manage 105:283–293
Bergström R, Bergqvist G (1999) Large herbivore browsing on conifer seedlings related to seedling morphology. Scand J Res 14:361–367. https://doi.org/10.1080/02827589950152683
Bernard M, Boulanger V, Dupouey J et al (2017) Deer browsing promotes Norway spruce at the expense of silver fir in the forest regeneration phase. For Ecol Manage 400:269–277. https://doi.org/10.1016/j.foreco.2017.05.040
Bernard M, Barrere J, Granda E, Benavides R, Rabasa SG, Heurich M, Kunstler G (2023) Data from “Elevation affects both the occurrence of ungulate browsing and its effect on tree seedling growth for four major tree species in European mountain forests. Zenodo. https://doi.org/10.5281/zenodo.10370097
Boulanger V, Baltzinger C, Saïd S et al (2009) Ranking temperate woody species along a gradient of browsing by deer. For Ecol Manage 258:1397–1406. https://doi.org/10.1016/j.foreco.2009.06.055
Boulanger V, Baltzinger C, Saïd S et al (2015) Decreasing deer browsing pressure influenced understory vegetation dynamics over 30 years. Ann Sci 72:367–378. https://doi.org/10.1007/s13595-014-0431-z
Brodie J, Post E, Watson F, Berger J (2012) Climate change intensification of herbivore impacts on tree recruitment. Proc Biol Sci 279:1366–1370. https://doi.org/10.1098/rspb.2011.1501
Cagnacci F, Focardi S, Heurich M, et al (2011) Partial migration in roe deer: migratory and resident tactics are end points of a behavioural gradient determined by ecological factors
Cairns DM, Moen J (2004) Herbivory influences tree lines. J Ecol 92:1019–1024
Callaway RM (1998) Competition and facilitation on elevation gradients in subalpine forests of the Northern Rocky Mountains, USA
Canham CD, McAninch JB, Wood DM (1994) Effects of the frequency, timing and intensity of simulated browsing on growth and mortality of tree seedlings. Can J Res 24:817–825
Champagne E, Tremblay JP, Côté SD (2016) Spatial extent of neighboring plants influences the strength of associational effects on mammal herbivory. Ecosphere 7:1–13. https://doi.org/10.1002/ecs2.1371
Champagne E, Turgeon R, Munson AD, Raymond P (2021) Seedling response to simulated browsing and reduced water availability: insights for assisted migration plantations. Forests 12:1396
Clark JS, Silman M, Kern R, Macklin E, HilleRisLambers J (1999) Seed Dispersal Near and Far: Patterns Across Temperate and Tropical Forests. Ecology 80:1475–1494
Diaci J (2002) Regeneration dynamics in a Norway spruce plantation on a silver fir-beech forest site in the Slovenian Alps. For Ecol Manage 161:27–38
Diaci J, Adamic T, Rozman A (2012) Gap recruitment and partitioning in an old-growth beech forest of the Dinaric Mountains: influences of light regime, herb competition and browsing. For Ecol Manage 285:20–28. https://doi.org/10.1016/j.foreco.2012.08.010
Felton AM, Wam HK, Felton A et al (2021) Macronutrient balancing in free-ranging populations of moose. Ecol Evol 11:11223–11240. https://doi.org/10.1002/ece3.7909
Fisichelli N, Frelich LE, Reich PB (2012) Sapling growth responses to warmer temperatures “cooled” by browse pressure. Glob Chang Biol 18:3455–3463. https://doi.org/10.1111/j.1365-2486.2012.02785.x
Fisichelli N, Vor T, Ammer C (2014) Broadleaf seedling responses to warmer temperatures “chilled” by late frost that favors conifers. Eur J Res 133:587–596. https://doi.org/10.1007/s10342-014-0786-6
Fuller RJ, Gill RMA (2001) Ecological impacts of increasing numbers of deer in British woodland. Forestry 74:193–199. https://doi.org/10.1093/forestry/74.3.193
Gordon IJ, Prins HHT (2008) The ecology of browsing and grazing. Springer-Verlag, Berlin
Hallinger M, Manthey M, Wilmking M (2010) Establishing a missing link: warm summers and winter snow cover promote shrub expansion into alpine tundra in Scandinavia. New Phytol 186:890–899. https://doi.org/10.1111/j.1469-8137.2010.03223.x
Hawkes CV, Sullivan JJ (2001) The impact of herbivory on plants in different resource conditions: a meta-analysis. Ecology 82:2045–2058. https://doi.org/10.1890/0012-9658(2001)082[2045:TIOHOP]2.0.CO;2
Heggberget TM, Gaare E, Ball JP (2002) Reindeer (Rangifer tarandus) and climate change: importance of winter forage. Rangifer 22:13–31
Heinrichs S, Winterhoff W, Schmidt W (2012) Vegetation dynamics of beech forests on limestone in central Germany over half a century – effects of climate change, forest management, eutrophication or game browsing? Biodiver Ecol 4:49–61. https://doi.org/10.7809/b-e.00059
Herrero A, Zamora R, Castro J, Hódar JA (2012) Limits of pine forest distribution at the treeline: herbivory matters. Ecology 213:459–469. https://doi.org/10.1007/s11258-011-9993-0
Herrero A, Almaraz P, Zamora R et al (2016) From the individual to the landscape and back: time-varying effects of climate and herbivory on tree sapling growth at distribution limits. J Ecol 104:430–442. https://doi.org/10.1111/1365-2745.12527
Hijmans RJ, Cameron SE, Parra JL et al (2005) Very high resolution interpolated climate surfaces for global land areas. Int J Climatol 25:1965–1978. https://doi.org/10.1002/joc.1276
Jiang F, Cadotte MW, Jin G (2022) Size- and environment-driven seedling survival and growth are mediated by leaf functional traits. Proc Biol Sci 289:20221400. https://doi.org/10.1098/rspb.2022.1400
Karger DN, Conrad O, Böhner J et al (2017) Climatologies at high resolution for the earth’s land surface areas. Sci Data 4:1–20. https://doi.org/10.1038/sdata.2017.122
Karger DN, Conrad O, Böhner J, et al (2018) Data from: climatologies at high resolution for the earth’s land surface areas. Dryad, dataset
Körner C (2007) The use of “altitude” in ecological research. Trends Ecol Evol 22:569–574. https://doi.org/10.1016/j.tree.2007.09.006
Körner C (2012) Alpine treelines. Springer, Basel
Kropil R, Smolko P, Garaj P (2015) Home range and migration patterns of male red deer Cervus elaphus in Western Carpathians. Eur J Wildl Res 61:63–72. https://doi.org/10.1007/s10344-014-0874-4
Kuijper DPJ, Cromsigt JPGM, Churski M et al (2009) Do ungulates preferentially feed in forest gaps in European temperate forest? For Ecol Manage 258:1528–1535. https://doi.org/10.1016/j.foreco.2009.07.010
Kullberg Y, Bergström R (2001) Winter browsing by large herbivores on planted deciduous seedlings in southern Sweden. Scand J Res 16:371–378. https://doi.org/10.1080/02827580117385
Kullberg Y, Welander NT (2003) Effects of simulated winter browsing and drought on growth of Quercus robur L. seedlings during establishment. For Ecol Manage 173:125–133
Kunstler G, Albert CH, Courbaud B et al (2011) Effects of competition on tree radial-growth vary in importance but not in intensity along climatic gradients. J Ecol 99:300–312. https://doi.org/10.1111/j.1365-2745.2010.01751.x
Kuznetsova A, Brockhoff PB, Christensen RHB (2017) lmerTest package: tests in linear mixed effects models. J Stat Softw 82:1–26. https://doi.org/10.18637/JSS.V082.I13
Laurent L, Mårell A, Korboulewsky N et al (2017) How does disturbance affect the intensity and importance of plant competition along resource gradients ? For Ecol Manage 391:239–245. https://doi.org/10.1016/j.foreco.2017.02.003
Mann TE (2009) Vegetationsökologisches Monitoring im Nationalpark Harz unter besonderer Berücksichtigung des Schalenwild-Einflusses und der Waldstruktur. Cuvillier Verlag, Göttingen
Marcora PI, Renison D, País-Bosch AI et al (2013) The effect of altitude and grazing on seedling establishment of woody species in central Argentina. For Ecol Manage 291:300–307. https://doi.org/10.1016/j.foreco.2012.11.030
Milner JM, Bonenfant C, Mysterud A et al (2006) Temporal and spatial development of red deer harvesting in Europe: biological and cultural factors. J Appl Ecol 43:721–734. https://doi.org/10.1111/j.1365-2664.2006.01183.x
Motta R (2003) Ungulate impact on rowan (Sorbus aucuparia L.) and Norway spruce (Picea abies (L.) Karst.) height structure in mountain forests in the eastern Italian Alps. For Ecol Manage 181:139–150. https://doi.org/10.1016/S0378-1127(03)00128-2
Muñoz Mazón M, Klanderud K, Sheil D (2022) Canopy openness modifies tree seedling distributions along a tropical forest elevation gradient. Oikos 2022:e09205. https://doi.org/10.1111/oik.09205
Nichols RV, Cromsigt JPGM, Spong G (2015) DNA left on browsed twigs uncovers bite-scale resource use patterns in European ungulates. Oecologia 178:275–284. https://doi.org/10.1007/s00442-014-3196-z
Palmer SCF, Truscott AM (2003) Browsing by deer on naturally regenerating Scots pine (Pinus sylvestris L.) and its effects on sapling growth. For Ecol Manage 182:31–47. https://doi.org/10.1016/S0378-1127(03)00026-4
Pepin N, Bradley RS, Diaz HF et al (2015) Elevation-dependent warming in mountain regions of the world. Nat Clim Chang 5:424–430. https://doi.org/10.1038/nclimate2563
Post E, Pedersen C (2008) Opposing plant community responses to warming with and without herbivores. Proc Natl Acad Sci 105:12353–12358. https://doi.org/10.1073/pnas.0802421105
R Core Team (2019) R: A Language and Environment for Statistical Computing. R Foundation for Statistical Computing. Vienna
Rabasa SG, Granda E, Benavides R et al (2013) Disparity in elevational shifts of European trees in response to recent climate warming. Glob Chang Biol 19:2490–2499. https://doi.org/10.1111/gcb.12220
Rameau J-C, Mansion D, Dumé G (1993) Flore forestière française (guide écologique illustré), tome 1: Plaines et collines. Institut Pour le Développement Forestier, Paris
Redjadj C, Darmon G, Maillard D et al (2014) Intra- and interspecific differences in diet quality and composition in a large herbivore community. PLoS One 9:1–13. https://doi.org/10.1371/journal.pone.0084756
Reimoser F, Putman RJ (2011) Impacts of wild ungulates on vegetation: costs and benefits. In: Putman RJ, Apollonio M, Andersen R (eds) Ungulate management in Europe, Cambridge. Cambridge University Press, Cambridge, p 144–191
Renaud PC, Verheyden-Tixier H, Dumont B (2003) Damage to saplings by red deer (Cervus elaphus): effect of foliage height and structure. For Ecol Manage 181:31–37. https://doi.org/10.1016/S0378-1127(03)00126-9
Ruttan A, Lortie CJ (2013) A systematic review of the attractant-decoy and repellent-plant hypotheses: do plants with heterospecific neighbours escape herbivory? J Plant Ecol 8:337–346
Sanders-DeMott R, McNellis R, Jabouri M, Templer PH (2018) Snow depth, soil temperature and plant–herbivore interactions mediate plant response to climate change. J Ecol 106:1508–1519. https://doi.org/10.1111/1365-2745.12912
Schirpke U, Tappeiner U, Tasser E (2019) A transnational perspective of global and regional ecosystem service flows from and to mountain regions. Sci Rep 9:6678. https://doi.org/10.1038/s41598-019-43229-z
Segar J, Pereira HM, Baeten L et al (2022) Divergent roles of herbivory in eutrophying forests. Nat Commun 13:7837. https://doi.org/10.1038/s41467-022-35282-6
Speed JDM, Austrheim G, Hester AJ, Mysterud A (2010) Experimental evidence for herbivore limitation of the treeline. Ecology 91:3414–3420
Speed JDM, Austrheim G, Hester AJ, Mysterud A (2011) Growth limitation of mountain birch caused by sheep browsing at the altitudinal treeline. For Ecol Manage 261:1344–1352. https://doi.org/10.1016/j.foreco.2011.01.017
Speed JDM, Austrheim G, Hester AJ, Mysterud A (2013) The response of alpine Salix shrubs to long-term browsing varies with elevation and herbivore density. Arct Antarct Alp Res 45:584–593. https://doi.org/10.1657/1938-4246-45.4.584
Stone PH, Carlson JH (1979) Atmospheric lapse rate regimes and their parameterization. J Atmos Sci 36:415–423
Storms D, Aubry P, Hamann J et al (2008) Seasonal variation in diet composition and similarity of sympatric red deer Cervus elaphus and roe deer Capreolus capreolus. Wildlife Biol 14:237–250
Stritih A, Bebi P, Rossi C, Grêt-Regamey A (2021) Addressing disturbance risk to mountain forest ecosystem services. J Environ Manage 296:113188. https://doi.org/10.1016/j.jenvman.2021.113188
Sundqvist MK, Sanders NJ, Wardle DA (2013) Community and ecosystem responses to elevational gradients: processes, mechanisms, and insights for global change. Annu Rev Ecol Evol Syst 44:261–280
Szwed M, Pińskwar I, Kundzewicz ZW et al (2017) Changes of snow cover in Poland. Acta Geophys 65:65–76. https://doi.org/10.1007/s11600-017-0007-z
Thuiller W, Lavorel S, Araujo MB et al (2005) Climate change threats to plant diversity in Europe. Proc Natl Acad Sci 102:8245–8250
Visscher DR, Merrill EH, Fortin D, Frair JL (2006) Estimating woody browse availability for ungulates at increasing snow depths. For Ecol Manage 222:348–354. https://doi.org/10.1016/j.foreco.2005.10.035
Vuorinen KEM, Rao SJ, Hester AJ, Speed JDM (2020) Herbivory and climate as drivers of woody plant growth: do deer decrease the impacts of warming? Ecol Applic 30:e02119. https://doi.org/10.1002/eap.2119
Wise MJ, Abrahamson WG (2005) Beyond the compensatory continuum: environmental resource levels and plant tolerance of herbivory. Oikos 109:417–428. https://doi.org/10.1111/j.0030-1299.2005.13878.x
Zimmermann NE, Edwards TC, Moisen GG et al (2007) Remote sensing-based predictors improve distribution models of rare, early successional and broadleaf tree species in Utah. J Appl Ecol 44:1057–1067. https://doi.org/10.1111/j.1365-2664.2007.01348.x
Acknowledgements
We are grateful to Kieran Leigh-Moy and Suzanne van Beeck Calkoen for their help in collecting field data in the Bavarian Forest National Park. We also thank Peter Kjellander, Tomasz Podgorski and Christoph Dreiser for providing us data about ungulate presence and/or snow cover in Gällivare, Krościenko and Freiburg, respectively.
Code availability
The R script to reproduce the analyses is freely available on GitHub at https://github.com/jbarrere3/BaccaraPaper/tree/main.
Funding
Financial support was provided by the European Union with the project BACCARA (FP7-226299, 7FP) and the Spanish Ministry for Innovation and Science with the grant Consolider Montes (CSD2008_00040). This study was funded by the French National Forest Office (ONF), French Office for Biodiversity (OFB) and National Centre for Scientific Research (CNRS) under partnership agreement 2016/05/6171. XM received support from the project DISTIMACC (BGF, ECOFOR 2014–23). GK received support from the REFORCE - EU FP7 ERA-NET Sumforest 2016 through the call ‘Sustainable forests for the society of the future’, with the ANR as national funding agency (grant 408 ANR-16-SUMF-0002) and from the ANR DECLIC (grant ANR-20-CE32-0005-01).
Author information
Authors and Affiliations
Contributions
MB, JB, XM, SS, VB, HJ, EG, SGR, FV and GK conceived the ideas and designed methodology. EG, RB, SGR, GK and MH collected the data. MB, JB and GK analysed the data. MB, JB, XM and GK led the writing of the manuscript. All authors contributed critically to the drafts and gave final approval for publication.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
Not applicable.
Consent for publication
All authors gave their informed consent to this publication and its content.
Competing interests
The authors declare that they have no conflict of interest.
Additional information
Handling editor: Erwin Dreyer
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Appendices
Appendix 1
1.1 Snow cover per site
Appendix 2
2.1 Climate in each site
2.1.1 General approach
To derive the fine scale spatio-temporal variation experienced by seedlings along each elevation gradient, we first downscaled the long-term mean monthly temperature (over the period 1960–1990) of WorldClim database, version 1 (Hijmans et al. 2005), at the transect level using the local elevation based on a moving-window regression method (Kunstler et al. 2011; Zimmermann et al. 2007). Then, to account for inter-annual variability, we used data from CHELSA portal (T mean), version 1.2 online (Karger et al. 2017, 2018), to derive annual anomalies and correct the long-term monthly mean temperature of WorldClim. Lastly, we also extracted at each plot location the mean annual precipitation (P mean) from CHELSA portal, version 1.2 online (Karger et al. 2017, 2018).
2.1.2 Correction of anomalies and calculation of mean winter temperature
All sites but Bavarian Forest National Park
All sites but Bavarian Forest National Park (BFNP) were sampled in 2010 so we considered mean December 2009 and January and February 2010 temperature data. At each plot of all sites but BFNP, mean winter temperature (T wint) was calculated for each month as
where ΔT m_meanCHELSA is the anomaly of the monthly CHELSA temperature for a given year (computed as the difference between the monthly temperature of that year and the long term).
Then, for the winter before sampling: T wint = (T Dec09 + T Jan10 + T Feb10) / 3
BFNP
BFNP was sampled in 2017, so we considered mean December 2016, January 2017 and February 2017 temperature data. However, CHELSA version 1.2 online contains data up to 2013. Since we had data measured from weather stations at this site, we decided to correct downscaled mean monthly temperature data by the anomalies of the weather stations. Therefore, at each plot of BFNP, mean winter temperature was calculated for each month as
Then, for the winter before sampling: T wint = (T Dec16 + T Jan17 + T Feb17) / 3
2.1.3 Climate along each elevation gradient
Appendix 3
3.1 Height distribution of seedlings
In plots where we found more than 20 seedlings of a given species, we randomly selected 20 of the available seedlings. To verify the randomness of our selection, we compared the height distribution of seedlings when the number of seedlings of this species is lower than 20 (i.e. no random selection) or not (i.e. random selection).
Overall, the height distribution of each species is consistent regardless of whether seedlings were randomly selected or not.
Appendix 4
4.1 Effect of elevation on browsing probability when the Swedish site Gällivare is included
Appendix 5
Effect of browsing and elevation on seedling growth when the Swedish site Gällivare is included
Appendix 6
6.1 Structure of the residuals of the growth model
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Bernard, M., Barrere, J., Morin, X. et al. Elevation affects both the occurrence of ungulate browsing and its effect on tree seedling growth for four major tree species in European mountain forests. Annals of Forest Science 81, 13 (2024). https://doi.org/10.1186/s13595-024-01226-x
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s13595-024-01226-x