Corresponding author: Simon Knüsel ( simon.knuesel@wsl.ch ) Academic editor: Ingo Kowarik
© 2019 Simon Knüsel, Marco Conedera, Harald Bugmann, Jan Wunder.
This is an open access article distributed under the terms of the Creative Commons Attribution License (CC BY 4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Citation:
Knüsel S, Conedera M, Bugmann H, Wunder J (2019) Low litter cover, high light availability and rock cover favour the establishment of Ailanthus altissima in forests in southern Switzerland. NeoBiota 46: 91-116. https://doi.org/10.3897/neobiota.46.35722
|
Future forest composition is uncertain in many areas due to climate change. The spread of non-native species adds to these uncertainties, particularly in forests recently colonised by novel tree species. To anticipate future forest composition, and thus the provision of ecosystem services, a thorough understanding of the factors influencing the establishment of non-native tree species is essential. We studied the presence and abundance of regeneration of Ailanthus altissima (Mill.) Swingle in 89 plots on a regular grid in three sites in southern Switzerland to determine the main drivers of its establishment. All sites are located in abandoned, i.e. currently unmanaged stands of Castanea sativa Mill. Propagule pressure is caused by single mature A. altissima that colonised the area ca 40 years ago. We found high rock cover, low litter cover and high light availability to be the most important predictors for the presence of A. altissima regeneration, whereas its abundance was positively influenced by high light availability, low litter cover and high browsing on regeneration of competing species. However, the presence models performed much better than the abundance models. Interestingly, the most important factors favouring the establishment of A. altissima in recently undisturbed sites were found to be similar in a nearby site after a severe forest fire, which suggests a similar establishment strategy after a disturbance as in recently undisturbed forests. Based on our results we expect a further expansion of the species in lowland forests currently dominated by C. sativa, likely controlled primarily by light availability.
Ailanthus altissima (tree of heaven), browsing, Castanea sativa (sweet chestnut), Gap Light Index (GLI), invasive tree, light availability, litter cover, non-native tree, propagule pressure
Many forest ecosystems are expected to face compositional changes as a consequence of climate change, the introduction and spread of novel tree species, or interactions thereof (
The establishment of non-native tree species has mainly been studied in disturbed forest ecosystems (e.g.,
In this study, we focus on the establishment of the early-successional tree species Ailanthus altissima (Mill.) Swingle. Its spread in forest ecosystems has often been associated with disturbances induced by forest management (
i) What are the most important environmental factors that drive the presence and abundance of regeneration of A. altissima in recently undisturbed forests?
Once identified, we compared these factors with the environmental drivers of the post-fire establishment of A. altissima in the same study region (
ii) What are the key differences in the environmental drivers of the establishment of A. altissima in disturbed versus recently undisturbed forests?
Knowledge on key differences in the establishment requirements of A. altissima in forests after a disturbance versus recently undisturbed forests may facilitate predictions on the spread of the species.
Ailanthus altissima is a dioecious early successional tree species originating from China that was brought to Europe around the 1740s as an ornamental tree (
The study was conducted in the sweet chestnut (C. sativa) forests of southern Switzerland. The climate in the study area is characterized by a mild average annual temperature of ca 12.7 °C and a mean annual precipitation of 1862 mm (climate normal 1981–2010, MeteoSwiss station Locarno-Monti, 383 m a.s.l.). A large amount of precipitation falls in short and heavy spells during the growing season (1023 mm from May to September) followed by a relatively dry period in winter. The soils are classified as haplic podzol (cryptopodzol) on crystalline bedrock (
The studied forests used to be dominated by C. sativa, which was introduced by the Ancient Romans more than 2000 years ago (
We selected three study sites in southern Switzerland based on an occurrence map of A. altissima (
To get an overview of the colonization history as well as the existing propagule pressure of A. altissima in the study sites, we searched for and mapped all seed-bearing female A. altissima within the boundary of the study sites and in a buffer of approximately 300 m around the studied forests. Of all accessible female trees, an increment core was extracted at a height of ca 1.3 m to determine their age. In the site Locarno, only a single female tree was found with an age of ca 18 years. In Sementina, the oldest of the four female trees was ca 27 years old, while the other three were only ca 15 years old. In the site San Vittore, several adult and at least 40-years-old A. altissima trees that grew along the road crossing the study site (Fig.
Study area and study sites. Study area in southern Switzerland (left column) and selected study sites: Locarno (a; 46°10'36"N, 8°46'12"E), Sementina (b; 46°10'53"N, 8°58'13"E) and San Vittore (c; 46°14'19"N, 9°04'46"E; right column with scale bars representing 500 m). (d) Represents the study area analysed by
We applied a systematic sampling design that was centred on A. altissima. An adaptation was necessary to obtain a sufficient number of sampling points with regeneration of A. altissima, as the studied forests are located at the distribution frontier where the species occurs at low densities only. To this end, we placed a regular grid with a mesh size of 50 m over all three study sites and defined potential sampling points at the nodes of the grid. At each of these points we searched for A. altissima saplings (>40 cm tall and <12 cm DBH) in a circular area of 800 m2 (ca 16 m radius) around the sampling points. Plants smaller than 40 cm were not considered to avoid individuals that had germinated in the year of the field survey (
At the absence points, two nested, slope-corrected circular plots with sizes of 100 and 200 m2 were temporarily installed at the nodes of the grid. At presence points, in contrast, the nested plots were not installed at the nodes of the grid, but were moved to have the nearest A. altissima sapling in the centre.
A range of environmental variables were then measured in each plot. First, the light environment was recorded at the plot centre with a hemispherical photograph (Canon EOS 50D camera with a Sigma EX DC 4.5 mm fisheye lens). The photograph was taken horizontally at a height of 1.7 m to avoid shading by small regeneration. If the A. altissima sapling representing the plot centre (later referred to as central A. altissima) at presence points was taller than 1.7 m it was bent to the forest floor to avoid self-shading. In addition, the age of the central A. altissima was assessed by counting growth units (
In the 100 m2 plot, we surveyed ground cover as well as seedlings (10–40 cm tall) and saplings (>40 cm tall) of all tree species. Percentage cover of bare soil (later referred to as mineral soil cover), litter, herbs and shrubs (later referred to as vegetation cover), coarse woody debris and vegetation-free rock (later referred to as rock cover) were assessed visually in 10% classes. Ground cover values were assessed so that they added up to 100%. For tree seedlings, the species and the presence of recent (i.e., last three years) browsing damages were noted. For saplings, we additionally recorded height, DBH and damage by fraying (only for saplings >130 cm tall).
In the 200 m2 plot, we assessed stand structure and signs of management (i.e., presence of tree stumps). For all trees and coppice shoots or stools (DBH >12 cm), DBH, height and the vitality class following
We used two types of multiple regression models to analyse A. altissima establishment. First, we used logistic models (LOG) to analyse the presence or absence of A. altissima saplings. Second, we used negative binomial models (NB) to analyse the abundance (number) of A. altissima saplings (see Methods S1 for details).
The light environment was characterized by the Gap Light Index (GLI), which integrates diffuse and direct irradiation over a growing season. GLI was calculated from the hemispherical photographs with the program Hemisfer (www.schleppi.ch/hemisfer;
Three of the five assessed ground cover variables, mineral soil, litter, and rock cover were used as explanatory variables (Table
All continuous variables were standardized before running the models by subtracting the mean and dividing the result by the standard deviation.
Environmental parameters recorded for the establishment models.
Continuous predictors | Mean (SD) | Range |
Light environment: | ||
Gap light index (GLI) | 21 (14.6) | 4–73 |
Ground cover: | ||
Mineral soil cover (%) | 15 (12.4) | 0–52 |
Litter cover (%) | 32 (20.8) | 0–85 |
Rock cover (%) | 18 (13.9) | 0–70 |
Coarse woody debris cover (%) | 5 (6.1) | 0–32 |
Vegetation cover (%) | 29 (24.7) | 0–92 |
Propagule pressure: | ||
Distance to seed source (m) | 127 (65.1) | 15–303 |
Stand characteristics: | ||
Basal area per plot (m2 ha−1) | 36 (21.3) | 5–121 |
Browsing: | ||
Browsed competitors (%) | 40 (18.0) | 0–100 |
Categorical predictors | Levels | n |
Site | Locarno | 29 |
San Vittore | 25 | |
Sementina | 35 |
Starting with the eight selected variables (Table
For the averaging of all models within 2 AICC points, we used the zero method to compute parameter estimates (
Goodness-of-fit of the models was assessed by comparing the most parsimonious models against the null model with a Chi-squared test. For the LOG models we further calculated AUC values with the R package pROC (
To discriminate factors driving the establishment of A. altissima after a forest fire (
In particular, we aimed at quantifying differences in ground cover and light availability in presence and absence points in this study and in fire affected plots in
The presence points were on average slightly closer to seed sources (118 m, ± 61.3 m standard deviation) than the absence points (136 m, ± 68.1 m), however not significantly (P > 0.2). In the 44 presence points, a total of 243 saplings of A. altissima were recorded. They had a mean height of 224 cm (± 179.9 cm), a mean DBH of 1.0 cm, and a mean age of 4 years (± 1.8 yr), with a maximum of 11 years. When present, A. altissima saplings often dominated the regeneration at the sampling point (57% ± 31.3% of all saplings, Table
In addition, we recorded a total of 188 seedlings of A. altissima, accounting on average for 16% (± 18.1%) of all seedlings at the presence points (Table
Only A. altissima showed an increasing relative frequency of individuals with increasing height (Fig.
Fractions of stem numbers per species for three different height classes. Data were pooled from all three sites.
Presence and abundance of regeneration. Shown are seedlings (10–40 cm tall) and saplings (>40 cm tall and <12 cm DBH) of A. altissima and the four most abundant co-occurring seedlings or saplings. Values are averages per sampling point, with standard deviations in parentheses.
Class | Species | Sampling points | |||
With A. altissima saplings | Without A. altissima saplings | ||||
n*ha−1 | proportion per sampling point (%) | n*ha−1 | proportion per sampling point (%) | ||
Seedlings | A. altissima | 494 (513) | 16 (18.1) | 163 (177) | 4 (6.5) |
Fraxinus spp. | 4023 (11312) | 35 (26.1) | 1082 (1812) | 20 (18.1) | |
Quercus spp. | 1161 (1375) | 28 (22.4) | 2146 (3329) | 37 (31.0) | |
C. sativa | 1068 (1377) | 20 (18.6) | 1912 (3174) | 42 (27.5) | |
R. pseudoacacia | 563 (392) | 15 (10.5) | 500 (831) | 10 (8.7) | |
Saplings | A. altissima | 560 (603) | 57 (31.3) | – | – |
C. sativa | 391 (429) | 30 (23.3) | 465 (470) | 72 (32.0) | |
R. pseudoacacia | 492 (776) | 26 (20.3) | 300 (490) | 18 (12.4) | |
Fraxinus spp. | 257 (234) | 16 (9.6) | 275 (236) | 45 (46.1) | |
S. aria | 500 (693) | 23 (31.6) | 229 (138) | 46 (29.1) |
The occurrence of A. altissima saplings was positively (P < 0.05) affected by rock cover and negatively (P < 0.001) by litter cover in LOGpars (Table
The three most important variables according to LOGavg, i.e. rock cover, litter cover and GLI, showed significant differences (P < 0.01) between presence-and absence points (Fig.
Comparison of ground cover (a–d) and light environment (e, f) between sampling points with or without saplings (>40 cm tall and <12 cm DBH) of A. altissima. Values originate from recently undisturbed forests (this study), and from a forest analysed four years after a fire (
Results of the logistic models of the occurrence of A. altissima saplings. The null model (LOGnull), the most parsimonious model (LOGpars) as well as an average of all models (n = 8) within 2 AICC points (LOGavg) are shown. Significant predictors are written in bold, values in parentheses represent standard errors.
LOGnull | LOGpars | LOGavg | |
(Intercept) | −0.05 | −0.19 | 0.21 |
(0.21) | (0.29) | (0.65) | |
Rock cover (%) | 0.74* | 0.82* | |
(0.36) | (0.41) | ||
Litter cover (%) | −1.86*** | −1.11 | |
(0.42) | (0.71) | ||
Gap Light Index (-) | 0.92 | ||
(0.68) | |||
Distance to seed source (m) | −0.44 | ||
(0.43) | |||
Mineral soil cover (%) | 0.12 | ||
(0.25) | |||
Site - San Vittore | −0.19 | ||
(0.51) | |||
Site - Sementina | −0.52 | ||
(1.00) | |||
AICC | 123.99 | 82.41 | – |
AUC | 0.50 | 0.88 | 0.91 |
Relative importance | |||
Rock cover | 1.00 | ||
Litter cover | 0.87 | ||
Gap Light Index | 0.83 | ||
Distance to seed source | 0.69 | ||
Mineral soil cover | 0.27 | ||
Site | 0.26 |
According to NBpars, the abundance of A. altissima saplings was significantly positively affected by GLI and the percentage of browsed competitors, whereas litter cover had a significantly negative effect (Table
NB
avg identified the percentage of browsed competitors, GLI and litter cover as the most important variables as well, supplemented by mineral soil cover, total basal area and rock cover. However, the influence of total basal area and rock cover were marginal only (Table
Results of the negative binomial model of the number of A. altissima saplings. The null model (NBnull), the most parsimonious model (NBpars) as well as a model average of all models (n = 6) within 2 AICC points (NBavg) are shown. Significant predictors are written in bold, values in parentheses represent standard errors.
NBnull | NBpars | NBavg | |
(Intercept) | 1.01*** | 0.53** | 0.51** |
(0.21) | (0.19) | (0.19) | |
Gap Light Index (-) | 0.70** | 0.75** | |
(0.24) | (0.26) | ||
Browsed competitors (%) | 0.45* | 0.45* | |
(0.18) | (0.18) | ||
Litter cover (%) | −0.64* | −0.53 | |
(0.28) | (0.33) | ||
Mineral soil cover (%) | 0.23 | ||
(0.23) | |||
Total basal area (m2ha-1) | 0.05 | ||
(0.13) | |||
Rock cover (%) | 0.01 | ||
(0.07) | |||
AICC | 349.95 | 330.73 | – |
Relative importance | |||
Gap Light Index | 1.00 | ||
Browsed competitors | 1.00 | ||
Litter cover | 0.89 | ||
Mineral soil cover | 0.61 | ||
Total basal area | 0.24 | ||
Rock cover | 0.11 |
Overall, A. altissima regeneration (10–300 cm tall) experienced less (P < 0.01) browsing damage compared to the four most frequent other species (Fraxinus spp., C. sativa, Quercus spp. and R. pseudoacacia) and an average of all other species (Fig.
Compared to the two species with most abundant saplings taller than 130 cm (C. sativa and R. pseudoacacia) as well as an average over all other species, A. altissima was frayed significantly less (56% ± 33% standard deviation, P < 0.05) than C. sativa (81% ± 27%) but did not differ significantly (P > 0.4) from R. pseudoacacia (75% ± 35%) or an average of all other species (52% ± 31%, Fig. S2).
Proportion of browsed regeneration. Browsing damages were recorded on saplings 10–300 cm tall for the three most recent years. Widths of the bars are proportional to the number of individuals recorded (n > 360). Different lowercase letters indicate significant (P < 0.05) differences between proportions.
Two of the three most important variables for the occurrence of A. altissima saplings identified in this study, litter cover and light availability, also differed significantly (P < 0.01) between presence- and absence points in a study conducted after a forest fire (
The number of A. altissima saplings was significantly negatively (P < 0.05) correlated with the number of S. aria saplings only (Table S6). However, S. aria was merely present in 10 sampling points, with a maximum of 13 saplings at one sampling point (cf. Table
The number of A. altissima saplings was positively (P < 0.001) correlated with the number of seedlings of R. pseudoacacia and Fraxinus spp., while there was no significant relationship (P > 0.5) with the other two most frequent seedling species, C. sativa and Quercus spp. (Table S5).
The amount of basal area of live trees of the four most prominent tree species (C. sativa, Quercus spp., B. pendula and R. pseudoacacia) at presence points did not show a significant relationship (P > 0.05) with A. altissima sapling density (Table S7).
The occurrence of A. altissima saplings was mainly favoured by high rock cover, low litter cover and high light availability (Table
The negative effect of litter cover for the occurrence and the abundance of regeneration of A. altissima may have been caused by a delay in germination and an increased invertebrate herbivory and seedling mortality, as found in an early successional site (
As expected, light availability played a major role for the occurrence of the early successional A. altissima. Although the species is able to persist in low light levels for several years (
The models for the abundance of A. altissima saplings had a much lower predictive power than the occurrence models, which implies that the following interpretations represent tendencies only. While high light availability and low litter cover were also important for the density of A. altissima regeneration, rock cover did not influence sapling abundance. This might be related to the lower availability of suitable sites for establishment in rocky forest patches (
However, the distribution of tree species is determined not only by the availability of suitable microsites for their establishment, but also by seed production and dispersal (
Vegetation cover, which was excluded from the models due to multicollinearity with light availability and litter cover, was significantly higher at presence points (Fig.
We did not consider soil characteristics as A. altissima is known to tolerate a broad range of soil conditions (
Ailanthus altissima suffered almost no damage from ungulate browsing compared to co-occurring regeneration (Fig.
Forest fires temporarily increase the availability of light and nutrients, and reduce competition (
Due to the early successional character of both A. altissima and R. pseudoacacia (
The increasing relative abundance of A. altissima with increasing development stage (Fig.
High light availability, low litter cover and to a lesser degree high rock cover resulted as the most important factors determining the occurrence and abundance of A. altissima regeneration in forests lacking of recent large scale disturbances like fire or coppice management (research question i). Interestingly, very similar factors were also important for the establishment of A. altissima after a disturbance (research question ii), suggesting a consistent establishment strategy across sites with different disturbance intensities. Thus, light availability remains a key factor that controls the establishment and spread of the species, particularly for generative regeneration analysed in this study.
We thank Florim Sabani, Anna Katharina Calmet, Vanessa Ingold and Laura Jirjahn for field assistance. Further, we thank Boris Pezzatti, Patrick Schleppi and Janet Maringer for help with data processing and analysis. We further thank Eric Gehring and Janet Maringer for helpful discussions on earlier versions of the manuscript and Luca Plozza for information on the study sites. This study was funded by the Swiss National Science Foundation (SNF) and the French National Research Agency (ANR) (Project-Nr: 310030L_156661), the Federal Office for the Environment (FOEN) and the Swiss Federal Research Institute for Forest, Snow and Landscape Research (WSL).