1. Introduction
Wild ungulates, such as roe deer (
Capreolus capreolus L.), red deer (
Cervus elaphus L.) and chamois (
Rupicapra rupicapra L.), browse on plants [
1,
2]. In particular, by browsing on the leader shoot (also called the terminal shoot), these animals can influence tree regeneration. Ungulates selectively browse specific ‘palatable’ tree species [
3,
4]. In addition, not all tree species are able to respond quickly and effectively to browsing (see review by [
5]). For example,
Abies alba Mill. saplings from high-elevation provenances showed a delayed response to simulated browsing, whereas
Picea abies (L.) H. Karst. and
Fagus sylvatica L. saplings did not [
6,
7]. Both effects can lead to browsing-induced retarded growth of palatable saplings (e.g., oak) relative to competing less palatable tree saplings [
8], sometimes changing the growth rate ranking between the species [
9,
10]. Observations, experiments and model simulations in many temperate and boreal forests have shown that selective browsing by wild ungulates can affect the development of a forest stand and lead to substantial changes in the composition and structure of plant communities [
11,
12], such as a decline of
A. alba [
13].
In forest stands with few tree species, like many mountain forests, biodiversity losses are more apparent than in forests with many tree species. Nonetheless, a loss of admix tree species in, e.g., vigorous beech forests (Milio-Fagetum or Carici-Fagetum) can have important consequences. In many sub- to lower montane beech forests in Switzerland,
F. sylvatica is expected to grow less vigorously and thus become less dominant in the future due to climate change. Therefore, other tree species should be promoted, such as
Acer platanoides L.,
Acer pseudoplatanus L.,
Fraxinus excelsior L.,
Quercus petraea (Matt.) Liebl. and
Quercus robur L., as well as
Acer campestre L.,
Sorbus aria (L.) Crantz and
Tilia cordata Mill. on drier sites (see recommendations in the Tree App [
www.tree-app.ch]). Further,
A. alba,
Prunus avium L.,
Taxus baccata L.,
Ulmus glabra Huds.,
Betula pendula Roth and
Juglans regia L. are conditionally recommended (Tree App). This means a wide range of tree species, some of which demand more light than
F. sylvatica (e.g.,
Q. robur [
14]) and/or are very attractive to wild ungulates (
Table 1; [
15]). Thus, ungulate browsing may also limit the silvicultural potential with respect to climate adaptation.
In many vigorous forests, foresters and wildlife managers evaluate the situation regarding the browsing influence of wild ungulates only through estimations, often with very different outcomes. In simplified terms, the latter consider the abundance of all saplings, while the former focus on the browsed saplings. Browsing percentage, i.e., the number of saplings with a browsed leader shoot, expressed as a percentage of all observed saplings, is often used as an objective and reproducible measure of the frequency of browsing on tree saplings [
16]. However, frequent browsing does not imply a relevant browsing impact. Indeed, browsing percentage alone captures little of the effective long-term influence of browsing on tree regeneration, such as reduction in stem number or loss of tree species in the future stand [
17]. Apart from the proportion of browsed saplings, at least the following factors have been stated to be important for estimating the long-term influence of browsing: (i) the spatial distribution and quantity (e.g., density) of tree saplings of all species; (ii) the within-tree browsing intensity (i.e., if only buds or large parts of the leader shoot are eaten) per season; (iii) the site-specific height growth of the tree saplings, and thus the time needed for trees to grow out of the reach of browsers, and the browsing-induced changes in growth rate ranking between the differently selected species; (iv) the possible delay in the tree response after browsing; and (v) the tree mortality induced by browsing (e.g., [
18,
19]). Mortality caused by ungulates can only be assessed by individually marking saplings and tracking them over time (e.g., [
20]) or by comparing fenced and unfenced sites (e.g., [
21]). Both methods are very costly and time-consuming. In addition, tree mortality due to browsing is often less important than the number of tree saplings that survive [
19]. Determining browsing-induced mortality is, therefore, usually not necessary for an initial assessment of the browsing impact. In contrast, the first four of the factors listed above can be assessed easily in inventories by using the k-tree method, i.e., by measuring and assessing at least the two tree saplings nearest to the plot center per species and height class (i.e., k > 1; see k-tree sampling [
22] or inventories by [
23,
24,
25]).
Alternative methods in which all tree saplings are counted within a predefined area (i.e., fixed-area plots [
26,
27,
28] or plot count methods [
29]) have greater accuracy for tree density assessments. Theoretically, height increments and each individual tree’s response after browsing could also be evaluated with these methods. However, in stands with a large number of saplings, this would mean that many trees would have to be measured, and many trees of the same species or height class would be sampled per plot. Generally, the large number of trees in densely regenerated areas is the main reason why the selected sample area size in large plot-count inventories is very small (e.g., 0.9 m in the Swiss National Forest Inventory [
30]), at the expense of sampling of rarer tree species. The k-tree method is much faster and more efficient for collecting information on rarer species as well. On the one hand, only some saplings of common species, selected according to clear criteria, are measured (the k trees nearest to the plot center). On the other hand, large areas can be searched for rarer species (up to the maximum defined search distance [
29]). The k-tree method also has disadvantages, as density calculations based on this method are prone to biases [
27,
31]. However, if the density calculation is omitted and the focus is instead laid on the calculation of the occupied area per species [
29] or on the stocking percentage [
27], the advantages outweigh the disadvantages in terms of the detailed evaluation of trees [
25,
32] selected according to their proximity to the plot center (and thus the method is not biased by choosing the largest trees per species and plot, cf. [
33]).
As wild ungulates move in space and selectively choose their food, and as favorable areas for tree regeneration are irregularly distributed, the measured frequency of browsing depends on the spatial distribution of the sample plots. First, differences are expected between plots sampled in a systematic grid over the entire forest (i.e., in all stages of development and across all forest types) and plots sampled only in open sites where tree regeneration is more likely [
34]. Second, forest edge density has a large influence on deer [
35] and, thus, on browsing [
36]. Small ‘islands’ of forests within agricultural land may, therefore, have more browsing than larger forest complexes. Thus, an inventory adapted to local conditions and silvicultural practices is important.
The main objectives of this study were (i) to launch an objective but easily feasible survey to assess the impact of browsing in the species-rich, patchily distributed beech forests of the eastern Swiss Plateau; (ii) to estimate the browsing impact based on the first four factors listed above (i.e., excluding browsing-induced mortality); and (iii) to reduce the number of factors to be assessed in future inventories at this location or in comparable beech forests while still including the main factors needed to estimate not only the frequency but also the impact of browsing, e.g., the browsing-induced shifts in the growth rate ranking between the species.
The following hypotheses were tested:
Browsing is selective regarding preference for known palatable tree species and individual vigorously growing saplings. Species recommended for enhancing the adaptation of these beech forests to climate change are particularly selected by ungulates.
The vigorous growth of saplings leads to light browsing and thus few impacts.
Coniferous tree species are mainly lightly browsed in winter and respond in the following vegetation period by forming a new leader shoot. Thus, there are no delays in the tree response after browsing, which can be neglected in further surveys.
Deciduous tree species are mainly lightly browsed in summer and can respond within the same vegetation period. Therefore, a survey in autumn is optimal for capturing the impact of browsing on these tree species.
Browsing in forest islands is more frequent than in larger forest complexes, and this difference should be considered when systematic grids are established.
3. Results
Twenty-one different tree species were found (
Table 1). The most common species were
F. sylvatica,
A. alba,
A. pseudoplatanus,
F. excelsior,
P. abies and
P. avium. These species, however, were not equally distributed within the height classes. HC5 included 17.6% of all
F. sylvatica saplings but only 3.4% of
F. excelsior, 5.2% of
P. avium, 7.2% of
A. alba, 8.5% of
Picea and 8.6% of
A. pseudoplatanus saplings. Other relatively frequent but palatable and susceptible species, such as
T. baccata, Sorbus aucuparia L. and
Quercus spp., were not found in larger height classes (
Table 1). In contrast, 11.8% of
Ulmus spp. were found in HC5.
3.1. Stocked Stand Area and Spatial Distribution
No tree regeneration was found within the maximum search distance of 5 m in only 2 of the 133 surveyed plots. In all other plots, at least one tree sapling was found. The proportion of stocked area, with a stocking goal of 3000 saplings per hectare (artificial reduction to 1.03 m search distance), was 98.5 ± 0.8% for saplings in HC1 to HC5.
Though more
F. sylvatica saplings were measured (
Table 1), the proportion of stocked area was largest for
A. alba,
F. excelsior and
A. pseudoplatanus; (black symbols in
Figure 1).
Prunus avium was found in only ca. one-third of the plots. Nevertheless,
P. avium was quite evenly distributed throughout the Kirchberg forest district (
Figure 2). All species that were found in several plots were spread over the entire territory and not only in a small local area (
Figure 2). Just
S. aria was missing from the northern part, and
J. regia was missing from the southeastern part of the territory (
Figure 2).
The proportion of stocked area per developmental stage was also calculated. Young growth (3 plots), thicket (2 plots) and pole timber (12 plots) were too seldom within the 133 plots for this analysis. Picea was more widespread in mixed uneven-aged stands (79.2 ± 8.5%) than in young to old timber (40 ± 10% to 57.6 ± 8.7%). Fagus sylvatica was less widespread in young timber (60 ± 10%) than in medium to old timber and mixed uneven-aged stands (81.8 ± 6.8% to 83.3 ± 7.8%). Prunus avium was less frequent in medium timber (27.3 ± 7.9%) than in old timber (50.0 ± 8.7%). All other species showed no significant differences in the proportion of stocked area between the developmental stages.
The proportion of stocked area was higher in larger forested areas (total 108 plots) than in the forest islands (25 plots) for A. alba, T. baccata, Ulmus spp. and S. aria.
Notably, the very few plots with young growth and thicket had a mean canopy shading of 44 ± 18.2%. All other plots were rather shaded, i.e., canopy shading in plots with pole timber was 80 ± 9.5%, with young timber 82.4 ± 9.7%, with median timber 80.9 ± 11.8%, with old timber 82.9 ± 10.3%, and with mixed uneven-aged stands 83.3 ± 9.6%. Further, there was a negative correlation (cor = −0.55) between the percentage cover of competing vegetation and canopy shading. The competing plant species was mostly Rubus fruticosus.
There was no correlation between canopy shading and hill slope (cor = 0.08), canopy shading and elevation (cor = −0.014), or hill slope and elevation (cor = 0.05).
Figure 2.
Spatial distribution (longitude and latitude) of the tree species found in >10 plots in the Kirchberg forest district. Open black circles represent plots without saplings, blue filled circles represent plots with saplings in HC1 to HC4 (for which the proportion of stocked area was calculated), green filled circles represent plots with no HC1 to HC4 saplings but with an HC5 sapling, and orange filled circles represent plots with no HC1 to HC4 saplings but with an HC0 sapling (only possible for Abies alba and Taxus baccata).
Figure 2.
Spatial distribution (longitude and latitude) of the tree species found in >10 plots in the Kirchberg forest district. Open black circles represent plots without saplings, blue filled circles represent plots with saplings in HC1 to HC4 (for which the proportion of stocked area was calculated), green filled circles represent plots with no HC1 to HC4 saplings but with an HC5 sapling, and orange filled circles represent plots with no HC1 to HC4 saplings but with an HC0 sapling (only possible for Abies alba and Taxus baccata).
3.2. Frequency, Season and Within-Tree Browsing Intensity
The proportion of browsed area was largest for
A. alba (13.5 ± 3.0%; blue triangles in
Figure 1). The proportion of area occupied by saplings with other damage (drought, frost, insect) or with no leader shoot (e.g., after heavy browsing in the preceding year) was also highest for
A. alba (9.8 ± 2.6%), while it was only 3.0 ± 1.5% for
F. excelsior, 2.3 ± 1.3% for
A. pseudoplatanus and
P. avium, and 1.5 ± 1.1% for
P. abies and
F. sylvatica.
The proportion of browsed area out of the total stocked area was highest for
P. avium (30.6 ± 6.7%), followed by
A. alba (16.1 ± 3.6%),
F. excelsior (11.8 ± 3.1%),
A. pseudoplatanus (11.0 ± 3.0%),
F. sylvatica (8.4 ± 2.9%) and, finally,
P. abies (0%). Winter browsing was much more frequent than summer browsing for all species—both evergreen and deciduous, as the proportion of winter-browsed area out of the total stocked area was larger than the proportion of summer-browsed area out of the total stocked area (red and green x symbols in
Figure 1). The only exception was
F. sylvatica, with 4.2 ± 2.1% in both periods. The rarer the species, the higher the proportion of winter-browsed area out of the total stocked area (note that
Figure 1 is sorted according to descending species occurrence). This was also the case for some species that were classified as unpalatable, i.e., ungulate preference class 1 (
Table 1), such as
J. regia. The only outliers were
Carpinus betulus and
M. sylvestris, with 0% browsing. However, area proportions of very rare species should be interpreted with caution.
The only important differences in the proportion of browsed area out of the total stocked area regarding developmental stage were observed for F. excelsior, with less winter browsing in young timber (4.8 ± 4.8%) than in mixed uneven-aged stands (20.0 ± 9.2%). Taxus baccata, F. excelsior, P. avium and Ulmus spp. had significantly higher proportions of browsed areas out of the total stocked area in larger forests compared to forest islands.
Considering the two nearest saplings per height class, winter browsing was more frequent than summer browsing for all tree species (including
F. sylvatica), and light and heavy summer browsing were rare or even absent for some species (
Figure 3).
For some species, light winter browsing was more frequent than heavy winter browsing, including
T. baccata,
F. sylvatica,
A. pseudoplatanus and
P. avium. In contrast, heavy winter browsing was more frequent for
A. alba,
F. excelsior,
Ulmus,
S. aucuparia,
Quercus,
J. regia and
A. campestre, and all browsed
A. platanoides were heavily browsed (
Figure 3).
Few individuals of
F. sylvatica,
A. pseudoplatanus and
P. avium were browsed in winter and in summer. For
Ulmus,
S. aria and
Quercus, between 10% and 15% were browsed in both seasons. Notably,
J. regia was among the most frequently and also most heavily browsed species here (
Figure 3). Out of all assessed
A. alba of HC0, 7.3% were browsed.
Figure 3.
Percentage of browsed saplings of HC1 to HC4 per season and within-tree browsing intensity (including the nearest and second nearest saplings per height class). Only species with ≥13 observed trees in these height classes are displayed. The species are sorted into evergreen and deciduous species and according to descending occurrence, apart from Acer campestre and Acer platanoides for better comparison with Acer pseudoplatanus.
Figure 3.
Percentage of browsed saplings of HC1 to HC4 per season and within-tree browsing intensity (including the nearest and second nearest saplings per height class). Only species with ≥13 observed trees in these height classes are displayed. The species are sorted into evergreen and deciduous species and according to descending occurrence, apart from Acer campestre and Acer platanoides for better comparison with Acer pseudoplatanus.
3.3. Height Growth
Generally, RGR2021 values were rather small, as the median over all six frequent species was only 0.073 (blue line in
Figure 4), e.g., 7 cm for a tree height of 1 m or 3.5 cm for a tree height of 50 cm. RGR was larger in 2021 (mean = 0.118) than in 2022 (mean = 0.099, median = 0.066). The actual difference was probably even more pronounced given that the first and all further flushes were measured for the height increment in 2022, while only the last flush was measured in 2021. However, only six
Ulmus saplings and one single sapling of
P. abies,
A pseudoplatanus and
F. excelsior of HC1 to HC4 had a second flush in 2022. No
F. sylvatica saplings in HC1 to HC4 had a second flush, and only three saplings in HC5 had one. The overall growth rate ranking calculated based on the height increment in 2021 was
P. avium >
F. excelsior >
A. alba and
F. sylvatica ≥
A. pseudoplatanus ≥
P. abies. For all species, similar variables were important in explaining RGR2021 and height increment in 2021 (
Table 2). Canopy shading had a negative effect on all six main species and on all species together, except
P. avium and the other susceptible species (
Table 2). Tree height had a positive effect on the height increment of all species. For RGR, tree height had a negative influence on
F. sylvatica and
A. pseudoplatanus, and thus, also on all six main species or all species together. This indicates that the taller the tree, the larger the height increment, with a strictly linear relationship for
A. alba,
P. abies, F. excelsior and
P. avium but with a flattening for
F. sylvatica and
A. pseudoplatanus. Damage frequency negatively influenced tree growth for all species apart from
P. abies (
Table 2). Except for
P. abies, within-tree browsing intensity was an important variable explaining RGR2021 (
Figure 4) and height increment in 2021 (
Table 2). For
A. alba and F. excelsior, RGR2021 and height increment in 2021 were equal between saplings with no sign of browsing or damage on their leader shoot and for the lightly browsed saplings. For all the other species, lightly browsed saplings had a larger remaining terminal-shoot piece than unbrowsed trees (e.g.,
Figure 4). This suggests that the best-growing, most vital trees were browsed. Heavy browsing significantly reduced the height increment in 2021 and RGR2021 for all species (there were no
P. abies with heavy browsing;
Figure 3). There were too few saplings with a dried leader shoot or with other damage to interpret the results. Light twig browsing in winter had a negative effect on both height increment and RGR in the model with the six main species considered together, but in the separate models for single species, it only had a significant effect on RGR2021 of
F. excelsior. All the other variables, such as developmental stage, forest structure, elevation, hill slope, hunting zone and cover of competing vegetation, had no effect on height increment in 2021 or on RGR2021.
Similar results were obtained for RGR2022 as for RGR2021. However, as there was much less summer than winter browsing (
Figure 3), the summer within-tree browsing intensity dropped out of all single species models. However, ‘indirect’ browsing in earlier years had a negative effect, as damage frequency negatively affected RGR2022 of all species. Leaf browsing had no influence on leader shoot growth. Light twig browsing in winter reduced RGR2022 only in the single species model for
F. excelsior and in the model with the six main species considered together.
Table 2.
Regression model for the height increment in 2021 for saplings of height class (HC) 1 to HC4. The model with all 21 species included browsing preference by ungulate (BPU) class as a predictor, and models with the six main species included species (S) as a predictor. Susceptible species are defined in
Table 1, but the ones included as main species were excluded. N = number of species in the respective analysis, * = significant at
p ≤ 0.01 and ns = not significant. All saplings with no terminal shoot in 2021 were omitted. The random factor ‘plot’ was significant in all models except that for
Prunus avium (but was retained for comparison between models). Winter within-tree browsing intensity was highly significant in all models and thus is not shown.
Table 2.
Regression model for the height increment in 2021 for saplings of height class (HC) 1 to HC4. The model with all 21 species included browsing preference by ungulate (BPU) class as a predictor, and models with the six main species included species (S) as a predictor. Susceptible species are defined in
Table 1, but the ones included as main species were excluded. N = number of species in the respective analysis, * = significant at
p ≤ 0.01 and ns = not significant. All saplings with no terminal shoot in 2021 were omitted. The random factor ‘plot’ was significant in all models except that for
Prunus avium (but was retained for comparison between models). Winter within-tree browsing intensity was highly significant in all models and thus is not shown.
Regression Model | N | Species | Tree Height | Damage Frequency | Twig Browsing | Canopy Shading |
---|
All 21 species | 2226 | BPU * | 0.017 | −0.144 | * | −0.017 |
6 main species | 1954 | S * | 0.017 | −0.150 | * | −0.020 |
Abies alba | 334 | - | 0.020 | −0.122 | ns | −0.015 |
Picea abies | 231 | - | 0.015 | ns | ns | −0.019 |
Fagus sylvatica | 527 | - | 0.012 | −0.307 | ns | −0.021 |
Acer pseudoplatanus | 391 | - | 0.015 | −0.123 | ns | −0.022 |
Fraxinus excelsior | 363 | - | 0.027 | −0.193 | * | −0.015 |
Prunus avium | 108 | - | 0.022 | −0.239 | ns | ns |
Susceptible species | 232 | - | 0.019 | −0.164 | ns | ns |
Depending on the year, F. sylvatica, P. avium and F. excelsior required the least amount of time to grow out of the reach of roe deer, with ca. 15 years. P. abies required around 28 years, A. alba 28–34 years, and A. pseudoplatanus 27–36 years. Additionally, taking into account the growth of lightly browsed saplings reduced the time by around 1–2 years.
Figure 4.
Relative growth rate in 2021 (RGR2021) per within-tree browsing intensity category in winter 2021/2022. Saplings from HC1 to HC4 are included. The blue dotted line represents the median growth rate over all tree species. The numbers given in all panels indicate the number of saplings sampled in the respective categories for that species. Different lower-case letters indicate significant differences at p ≤ 0.05 between the within-tree browsing intensity in post-hoc tests (for cases where the variable was significant at p < 0.01 in the linear mixed effects models). Median (bold line), first and third quartile (bottom and top of box), quartile ± 1.5 × interquartile range (whiskers) and individual points more extreme in value (circles) were drawn using the boxplot function in the base R package. The width of the boxes represents the number of trees within the various categories.
Figure 4.
Relative growth rate in 2021 (RGR2021) per within-tree browsing intensity category in winter 2021/2022. Saplings from HC1 to HC4 are included. The blue dotted line represents the median growth rate over all tree species. The numbers given in all panels indicate the number of saplings sampled in the respective categories for that species. Different lower-case letters indicate significant differences at p ≤ 0.05 between the within-tree browsing intensity in post-hoc tests (for cases where the variable was significant at p < 0.01 in the linear mixed effects models). Median (bold line), first and third quartile (bottom and top of box), quartile ± 1.5 × interquartile range (whiskers) and individual points more extreme in value (circles) were drawn using the boxplot function in the base R package. The width of the boxes represents the number of trees within the various categories.
There were a considerable number of pairs of unbrowsed saplings of
A. alba,
P. abies,
A. pseudoplatanus,
F. excelsior and
P. avium to unbrowsed
F. sylvatica saplings of the same height class per plot (
Figure 5). However, there were only a few pairs with browsed saplings to unbrowsed
F. sylvatica saplings. At the sample plot level, none of these five species had a growth advantage over
F. sylvatica in the unbrowsed/undamaged condition. However, the ‘remaining’ RGR2021 of lightly browsed saplings of these four species was still larger than that of unbrowsed
F. sylvatica (winter within-tree browsing intensity
p < 0.01; no–light
p = 0.047). In contrast, the heavily browsed saplings had lower remaining RGR2021 than unbrowsed
F. sylvatica (no–heavy
p < 0.01, light–heavy
p < 0.01).
3.4. Response to Browsing
After heavy browsing in winter 2021/2022, 40% of
A. alba saplings, and a considerable number of
P. avium,
F. sylvatica, A. pseudoplatanus and
F. excelsior saplings, had no leader shoot in 2022, i.e., did not respond after heavy browsing in the dormant season (
Table 3). More
A. alba saplings responded after light than after heavy winter browsing. All
F. excelsior and
P. avium saplings with light winter browsing formed a new shoot in 2022. Surprisingly, more
F. sylvatica and
A. pseudoplatanus saplings did not respond after light than after heavy winter browsing (
Table 3). Only winter twig browsing significantly reduced the response (
p < 0.01).
Of the A. alba saplings with no shoot in 2021, 73.7% had no new leader shoot by the end of the vegetation season in 2022, i.e., had a time lag in their response by at least two years. For all the other main species, there were only one to three saplings (none for P. avium) with no leader shoot in 2021, while none of them had formed a new leader shoot in 2022.
After browsing during the vegetation period, only 10% of F. sylvatica, but at least 40% of A. pseudoplatanus, 50% of P. avium and 67% of F. excelsior saplings responded to light browsing with the formation of a new shoot within the same vegetation period.
Table 3.
Percentage of tree saplings with no leader shoot in 2022 after light or heavy winter browsing in 2021/2022 for the six main tree species of height class (HC) 1 to HC4 (Picea abies excluded, due to too few browsed saplings). Values in italics are based on only 14–24 saplings. Values in brackets include the HC0 saplings of Abies alba.
Table 3.
Percentage of tree saplings with no leader shoot in 2022 after light or heavy winter browsing in 2021/2022 for the six main tree species of height class (HC) 1 to HC4 (Picea abies excluded, due to too few browsed saplings). Values in italics are based on only 14–24 saplings. Values in brackets include the HC0 saplings of Abies alba.
Species | Light Winter Browsing [%] | Heavy Winter Browsing [%] |
---|
Abies alba | 28.6 (22.7) | 40.0 (36.7) |
Fagus sylvatica | 35.7 | 21.1 |
Acer pseudoplatanus | 22.9 | 15.6 |
Fraxinus excelsior | 0.0 | 17.9 |
Prunus avium | 0.0 | 27.8 |
3.5. Damage Frequency
The damage frequency ranking was P. abies < F. sylvatica < A. alba < A. pseudoplatanus. F. excelsior and P. avium had damage frequencies similar to those of A. alba and A. pseudoplatanus.
In the model for the damage frequency of the main six tree species, only tree species and tree height were significant. The greater the tree height, the higher the damage frequency (estimate ± std. error 0.0079 ± 0.0005). Notably, there were already very small
A. pseudoplatanus saplings with damage along the main stem axis, and there were only four
A. pseudoplatanus taller than 50 cm that did not have at least one damage point (
Figure 6). It should also be considered that it becomes more difficult to detect browsing damage from earlier years as a sapling becomes older and broader. Thus, it is likely that several of the larger saplings were damaged when they were small (10–20 cm) but that this damage was not detected when the saplings were 1–2 m tall.
4. Discussion
In the beech forests of the Kirchberg forest district in the canton of St. Gallen, Switzerland, 21 different tree species were found to regenerate naturally in the browsing impact survey conducted in 133 systematically placed plots. All tree species recommended for such beech forests, in order to cope with the changing climate, were found. However, only
Acer pseudoplatanus and
Fraxinus excelsior from the highly recommended and
Abies alba from the conditionally recommended species were very widespread.
Prunus avium was found in at least ca. one-third of the plots, and
Taxus baccata,
Quercus spp.,
Ulmus spp. and
Sorbus aucuparia in one of every five to seven plots. For these species, it can be assumed that there are enough seed trees to ensure natural regeneration. In addition, the stocking goal of 3000 tree saplings per hectare was fulfilled in all but two plots. Therefore, the natural potential to adapt to future conditions is present in these beech forests. However, browsing can jeopardize this adaptability to future climate conditions [
52]. We will discuss some aspects based on our hypotheses.
- 1.
Browsing is selective regarding preference for known palatable tree species and individual vigorously growing saplings. Species recommended for enhancing the adaptation of these beech forests to climate change are particularly selected by ungulates.
Browsing was indeed selective in terms of both species and individual trees. However, not only the species known to be palatable by wild ungulates were selected, and the differences were rather small.
Picea abies saplings were hardly browsed at all, as expected [
15]. However, the proportion of browsed area out of the total stocked area was not significantly different for
A. alba,
F. excelsior,
A. pseudoplatanus and the generally much less palatable
F. sylvatica. One possible reason for this is that these otherwise more palatable tree species were much more widespread than
F. sylvatica. It is known that wild ungulates do not select species equally everywhere (e.g., [
15,
53]) and in any season [
54] and often prefer rarer species [
55,
56], which makes these tree species even more rare [
57]. For example, in the Kirchberg forest district, the rarely admixed ‘unpalatable’
Juglans regia was also browsed frequently and heavily. This indicates that the animals in this region selected the common species (apart from
P. abies) according to their frequency but also selectively browsed rare species. Rarer species currently include
Acer platanoides and
A. campestre,
Quercus spp.,
Sorbus aria,
Tilia spp.,
T. baccata,
Ulmus glabra,
Betula pendula and
J. regia, which seem to be more suitable in terms of adaptation to future climates and are, therefore, recommended to foresters (TreeApp). Thus, ungulate browsing at least slows down the silvicultural potential concerning climate adaptation.
The most vigorously growing trees were browsed, so the length of the measured remnant after light browsing was equal to (
A. alba and
F. excelsior) or greater than the leader shoot length of the unbrowsed trees (all other species with browsing). This is in agreement with the plant vigor hypothesis [
58] and other studies in beech forests (e.g., [
40]).
- 2.
The vigorous growth of saplings leads to light browsing and thus few impacts.
The tree saplings did not grow particularly well in these supposedly vigorous forests (cf. description of the forest communities [
37]). Compared with other studies, RGR was rather low [
25]. However, it has to be considered that, with a mean shading of 80.7 ± 12.9%, the forests in the Kirchberg forest district were severely shaded. Even the shade-tolerant tree species
A. alba [
14] grows worse under >ca. 85–90% shade than under more open conditions [
20,
59]. In
A. alba plantations in the Czech Republic, canopy openness of 30% was found to be optimal [
60]. The best growth conditions for natural
A. alba seedlings were found under the relatively light-permeable crowns of light-demanding tree species (
Larix decidua,
Pinus sylvestris), e.g., ca. 43–45% canopy openness in the Karkonosze Mountains [
61]. In the case of
F. sylvatica, increased growth with increased light (relative light intensity between 0% and 35%) has been reported [
62]. The relatively shaded stands may also be one reason why
A. alba grew as well as
F. sylvatica, and
A. pseudoplatanus grew as well as (2021) or less than (2022)
F. sylvatica.
Quercus spp. definitely need more light than these more shade-tolerant species, so low light availability, in addition to browsing, is probably one of the main reasons for the absence of larger oaks (but see the possible role of large grazing herbivores, e.g., [
63]). The increment was smaller in 2022 than in 2021, probably due to the locally warm and dry summer of 2022. In addition to the shading, the warm and dry climate could be a reason why vegetation growth was completed early and practically no trees formed a second flush in 2022.
A hill slope of up to 30° was found to favor superior productivity classes for
A. alba [
64].
Abies alba saplings originating from steeper slopes had smaller height increments before and after browsing than those from flat areas [
6]. In both studies, hill slope and elevation were positively correlated. The average slope of the forests in the Kirchberg forest district was 33 ± 23%. Neither hill slope nor elevation had an influence on the height increment or damage frequency of any species. As hill slope, elevation and canopy shading were not correlated, we assume that canopy shading was indeed the main reason for the relatively poor growth.
When height increments are small, browsing is often heavier, as a single bite automatically removes a larger percentage of the terminal shoot. In fact, more heavy than light browsing was found for many species (
Figure 3). Heavy browsing usually leads to greater damage to the individual tree as (i) tree height is reduced more, (ii) reserves and meristems are lost and (iii) the new terminal shoot is, therefore, often shorter and sometimes formed with a time delay (cf. review in [
5]). For example,
F. sylvatica saplings grew equally with or without light-simulated browsing (terminal bud removal) but were negatively affected by heavy clipping [
7].
Despite the similar growth of
A. alba and
F. sylvatica, the time
A. alba needed to grow out of the reach of roe deer was double that required by
F. sylvatica. The same was true for
A. pseudoplatanus. However, apart from
F. sylvatica, there were rather few saplings in the larger height classes (
Table 1), and many or all (e.g.,
A. pseudoplatanus) of them had been browsed several times on their main stem (
Figure 6). Moreover, the linear mixed effects models revealed a negative relationship between damage frequency and height increment (
Table 2). The time needed to grow out of the reach of roe deer is, therefore, biased (prolonged) by ungulate browsing.
- 3.
Coniferous tree species are mainly lightly browsed in winter and respond in the following vegetation period. Thus, there are no delays in the tree response after browsing, which can be neglected in further surveys.
Conifers, i.e.,
A. alba, were almost exclusively browsed in winter. However, browsing was more often heavy, and many
A. alba saplings showed a delayed response even after light browsing. The increments of
A. alba were rather small, and thus it can be assumed that practically no preformed buds remained on the remaining pieces from which the
A. alba saplings could respond the following year. This is because trees that are taller prior to browsing have more buds remaining after browsing on the residual leader shoot sections [
65], and a certain tree height has to be reached to have internodal buds out of which new leader shoots can grow [
66]. In the Kirchberg forests, however, most
A. alba saplings (still) belonged to HC1. The observed other damage in
Figure 3 corresponds to
A. alba trees without terminal shoots. If we look only at the proportion of browsed area out of the total stocked area (or percentage of browsed assessed saplings, cf.
Figure 3), not many
A. alba were affected by browsing (compared with a 15.4–19% threshold of tolerable browsing intensity in beech forests in Aargau [
67] and Thurgau [
68]). By summing the browsed saplings and the ones with no leader shoot after browsing, however, every fourth
A. alba was affected. In addition, damage along the main stem axis reduced growth in the following years, i.e., the greater the damage frequency, the smaller the increment. This means that practically all
A. alba had reduced growth as a result of current or former browsing. In an experiment with planted but naturally browsed
A. alba saplings, the browsed saplings were not able to compensate partially. That is, even two years after browsing, they did not achieve the same height increment as unbrowsed saplings [
65]. Thus, browsing-induced height differences increase over time. The stands in the cited experiment were also Milio-Fagetum stands, though at about 1000 m a.s.l. and with a northeastern aspect [
65]. Therefore, it cannot simply be assumed that sites with a vigorous forest community have trees with a quick browsing response (no delay) and partial compensation.
Based on a clipping experiment in which
A. alba originating from the Swiss Plateau (lowland provenances) responded with no time lag even when planted at 1090 m a.s.l. [
6] and all
F. sylvatica provenances fully compensated for light clipping [
7], we had expected all tree species in the lowland Kirchberg forests to have a greater ability to recover from browsing. If only the browsing percentage had been measured, this lacking or at least delayed recovery would not have been detected. For this reason, it is important to carry out the browsing inventories in autumn so that the response to winter browsing can be considered in all forest stands, independent of vigor, tree species composition (coniferous and deciduous) and elevation.
- 4.
Deciduous tree species are mainly lightly browsed in summer and can respond within the same vegetation period. Therefore, a survey in autumn is optimal for capturing the impact of browsing on these species.
In contrast to our expectations, many deciduous tree species were browsed more frequently in winter than in summer, and many saplings were unable to respond to browsing during the vegetation season. Both results are probably related to the low availability of light. Due to the darkness under the dense canopy of these beech forests, the shrub and herb layers are often only weakly developed [
37]. This was also the case in the studied forests. This limits the food supply, especially in winter without herbs, and favors browsing in the dormant season [
55].
More vigorously growing saplings have longer leader shoots, and more buds remain on the browsed remnant, meaning that resprouting is possible sooner. Nonetheless, growth in the Kirchberg forests was limited by the dark stand conditions (cf. section above). Second and further flushes, especially in beech, are clearly more frequent when more light is available [
69] and are formed in fully open sites by all beech provenances [
7]. In the Kirchberg region, the stands were probably too dark for the formation of such further flushes. In addition, the growth of
F. sylvatica was found to be inversely proportional to maximum summer temperatures [
52], and 2022 was hotter than other years in the Kirchberg region. Therefore, it is not surprising that many trees were not able to form a second flush after browsing.
Moreover, of the frequent deciduous trees, some saplings also responded with a delay of at least one year to light and heavy winter browsing (
Table 3). As expected, all
F. excelsior and
P. avium saplings with light winter browsing formed a new shoot in 2022. In contrast, more
F. sylvatica and
A. pseudoplatanus saplings did not respond after light than after heavy winter browsing. We can only speculate that this result has to do with shading. However, in a different study, the sister species
Acer saccharum was not sufficiently stimulated to compensate for the loss of the old leader meristems caused by the removal of the uppermost buds [
70]. It could thus be a question of apical dominance in
A. pseudoplatanus. In
F. sylvatica, with sympodial branching according to the Troll architectural model [
71] and all buds being ‘the same’, this seems unlikely.
The larger the trees, the higher the damage frequency. No larger
A. pseudoplatanus saplings were without damage along the stem axis. The more damage there is, the smaller the height increment in the following year (in our study but also elsewhere [
72]). Thus, the poor responsiveness of
A. pseudoplatanus could also be due to a too large loss of reserves caused by former damage under already limited conditions (i.e., nonstructural carbohydrates decrease in darkness [
73]).
Twig browsing reduces the height growth of the leader shoots in
F. excelsior. It is known that when many twigs are browsed, reserve material in the twigs is lost [
74], which can reduce height growth in subsequent years [
75,
76]. Apparently, this was a sufficiently widespread phenomenon in Kirchberg for
F. excelsior but not for the other species.
Despite the less frequent summer browsing, we recommend keeping the autumn inventory in the Kirchberg forest district and carrying out inventories in autumn at other locations as well. First, an inventory should not be changed, if possible, so that long-term comparisons remain possible. Second, it is important to have information on how vigorous the ‘whole’ growth is to judge other beech forests (e.g., up to five flushes of
F. sylvatica in open sites of beech forests in Birmensdorf, Switzerland [
7], which would be recorded as separate years in inventories conducted in spring (or performed by assessing only the previous year, cf. Swiss National Forest Inventory [
39]). Third, it cannot be assumed that deciduous trees respond to browsing in the current year or in the following growing season. The response must, therefore, be checked explicitly.
- 5.
Browsing in forest islands is more frequent than in larger forest complexes, and this difference should be considered when systematic grids are established.
Unexpectedly, browsing in small forests isolated within agricultural land (forest islands) was not more frequent than browsing in larger forest complexes. The proportion of browsed area out of the total stocked area was not different or even smaller, and the variable ‘forest island’ dropped out of the models for height increment and damage frequency. This means that densification of the systematic grid in such small forests could be avoided. However, it is important to note the plot placement in this study. First, the sample plots were moved 15 m into the forest away from the forest edge area since the forest edge area in the canton of St. Gallen is upgraded and maintained specifically for wildlife and biodiversity and thus does not primarily serve for silvicultural production. However, this led to the omission of the smallest ‘forests’, less than approx. 50 m in diameter. Thus, to test whether browsing is more frequent or more intense at forest edges or in the smallest forest islands, no shift into the forest should be implemented. Second, the stocked area of four tree species was smaller in forest islands than in the larger forest complexes. Therefore, the mortality due to browsing could be higher in forest islands. Third, the form of the larger forest complexes probably plays a role. Many of these forests are intermeshed with the agricultural landscape (
Appendix A Figure A1). The distance to forest edges would, therefore, have to be explicitly considered (see, e.g., [
36]). In a very small-scale landscape, however, most plots have small distances to the forest edge.