ORIGINAL RESEARCH article
One Century of Forest Monitoring Data in Switzerland Reveals Species- and Site-Specific Trends of Climate-Induced Tree Mortality
- 1Swiss Federal Institute for Forest, Snow and Landscape Research WSL, Birmensdorf, Switzerland
- 2SwissForestLab, Birmensdorf, Switzerland
- 3Institute of Meteorology and Climate Research – Atmospheric Environmental Research, Karlsruhe Institute of Technology, Garmisch-Partenkirchen, Germany
- 4Institute of Terrestrial Ecosystems, ETH Zurich, Zurich, Switzerland
Climate-induced tree mortality became a global phenomenon during the last century and it is expected to increase in many regions in the future along with a further increase in the frequency of drought and heat events. However, tree mortality at the ecosystem level remains challenging to quantify since long-term, tree-individual, reliable observations are scarce. Here, we present a unique data set of monitoring records from 276 permanent plots located in 95 forest stands across Switzerland, which include five major European tree species (Norway spruce, Scots pine, silver fir, European beech, and sessile and common oak) and cover a time span of over one century (1898–2013), with inventory periods of 5–10 years. The long-term average annual mortality rate of the investigated forest stands was 1.5%. In general, species-specific annual mortality rates did not consistently increase over the last decades, except for Scots pine forests at lower altitudes, which exhibited a clear increase of mortality since the 1960s. Temporal trends of tree mortality varied also depending on diameter at breast height (DBH), with large trees generally experiencing an increase in mortality, while mortality of small trees tended to decrease. Normalized mortality rates were remarkably similar between species and a modest, but a consistent and steady increasing trend was apparent throughout the study period. Mixed effects models revealed that gradually changing stand parameters (stand basal area and stand age) had the strongest impact on mortality rates, modulated by climate, which had increasing importance during the last decades. Hereby, recent climatic changes had highly variable effects on tree mortality rates, depending on the species in combination with abiotic and biotic stand and site conditions. This suggests that forest species composition and species ranges may change under future climate conditions. Our data set highlights the complexity of forest dynamical processes such as long-term, gradual changes of forest structure, demography and species composition, which together with climate determine mortality rates.
Forest ecosystems play a crucial role in maintaining the balance of land-atmosphere cycles, sequestering carbon, fostering biodiversity and providing esthetic value to humanity (Bonan, 2008; Pan et al., 2011; Millar and Stephenson, 2015). Hence, understanding forest mortality is a pressing task of our times in order to predict and regulate the future development of forests. Although forest mortality is a natural process within forest stand dynamics (Franklin et al., 1987), mortality rates have been reported to increase in recent years due to climatic changes (Allen et al., 2010; McDowell et al., 2018), and are predicted to increase further in the future (Steinkamp and Hickler, 2015).
Drivers of tree mortality are diverse and intertwined with one another. The decline-disease theory is a conceptual framework relating tree mortality to a sequence of abiotic and biotic factors, which can be of predisposing, inciting or contributing nature (Manion, 1981; Houston, 1984). Predisposing factors (e.g., age or air pollutants) devitalize a tree over a long time period. Inciting factors, such as drought or insect defoliation, are short-term events that considerably reduce tree vigor. Finally, the occurrence of contributing factors (e.g., additional drought events, fungi) determines the fate of the weakened tree. Contributing factors can act in a long term, but can also cause a fast die-off when a certain threshold is reached (e.g., bark beetle infestations). Forest stand characteristics play an important role in forest mortality. For example, higher stand basal area and subsequent competition have been shown to increase mortality rates (Bradford and Bell, 2017; Young et al., 2017). Mortality risk also depends on tree or stand age (van Mantgem et al., 2009; Cailleret et al., 2017; Neumann et al., 2017) with usually a U-shaped curve of mortality probability in relation to age and/or tree size (Coomes and Allen, 2007; Lines et al., 2010). In recent years, changing climatic conditions, such as increasing drought intensity and frequency as well as rising temperatures, have been found to strongly influence forest mortality rates (van Mantgem et al., 2009; Allen et al., 2010; Anderegg et al., 2013). Disturbances that are often a consequence of changing climate, such as more numerous and severe fires (Brando et al., 2014; Stephens et al., 2018) or damaging insect outbreaks associated with storms and droughts (Weed et al., 2013; Kolb et al., 2016; Wood et al., 2018) have also been linked to increased forest mortality. Given this diversity of factors, mortality rates vary depending on the geographical location, its climate and consequent vegetation structure. Nevertheless, an increase in forest mortality due to climatic changes has been observed across continents, under different climatic conditions and across diverse forest types, from gymnosperm to angiosperm dominated ones, and from tropical to boreal regions (Allen et al., 2010; Anderegg et al., 2015; McDowell et al., 2018).
Drought-induced tree mortality can develop after substantial long periods of suboptimal water supply. The underlying physiological mechanisms include substantial damages to the tree hydraulic system through air entering the water conducting tissues of trees (Tyree and Sperry, 1989) under high xylem tensions induced by low soil moisture and high evaporative demand during drought. If these so called embolisms are substantial, then trees are not able to recover (Hartmann et al., 2018) following drought release and may subsequently die (Anderegg et al., 2016; Adam et al., 2017). Trees have developed multiple physiological strategies, at root, stem and leaf level, which allow them to cope with drought stress and to avoid hydraulic failure to a certain degree. For example, they can grow deep roots for better access to water, build xylem that is robust to embolism and/or control water loss via stomatal closure or leaf shedding (Martin-StPaul et al., 2017; Choat et al., 2018).
Growing recognition of the impacts of extreme droughts on forest ecosystems in Europe has spurred scientific attention on how major tree species including Scots pine (Pinus sylvestris L.), Norway spruce [Picea abies (L.) H. Karst], silver fir (Abies alba Mill.), European beech (Fagus sylvatica L.) and oak [Quercus petraea (Matt.) Liebl., Quercus robur L.] will cope with the expected changes in climate. These species follow different physiological strategies to cope with drought stress and exhibit different levels of drought tolerance (Breda et al., 2006). Among these species, Scots pine usually occurs on the driest sites used for commercial forestry in Central Europe. Compared to Norway spruce and silver fir, Scots pine is considered to be more drought tolerant. This is presumably due to its early stomatal closure to avoid xylem embolism under mild-moderate drought (Martinez-Vilalta et al., 2004) in combination with a deep rooting system (Richardson, 2000). Despite these strategies, Scots pine may not be well adapted to the combined effects of drought and heat (Dobbertin et al., 2005; Giuggiola et al., 2010) or to several consecutive drought events (Bigler et al., 2006) and is further vulnerable to diverse drought-related pests and diseases (Wermelinger et al., 2008). Norway spruce has repeatedly been described as particularly susceptible to drought (Zang et al., 2012, 2014; Pretzsch et al., 2013), which is likely due to its shallow root system and its vulnerability to drought-related insect outbreaks. There is an agreement that Norway spruce is already negatively affected by the impacts of climate change, which will be more pronounced in the future (Hanewinkel et al., 2013; Zang et al., 2014). Silver fir, however, has been shown to be more resilient to drought stress and associated phenomena like bark beetle outbreak or storm damages (Zang et al., 2014; George et al., 2015; Vitali et al., 2017). European beech is another drought sensitive species (Leuschner et al., 2001; Gessler et al., 2007; Charru et al., 2010). Several studies show a strong decline in radial growth of beech trees in response to decreasing water availability particularly for the Mediterranean region (Jump et al., 2006; Piovesan et al., 2008). However, beech has also shown a high physiological plasticity and adaptability to changing growing conditions (Vitasse et al., 2014; Cocozza et al., 2016; Stojnic et al., 2018), and the acclimation of beech to future climatic changes is highly uncertain. Oaks are generally considered to be well adapted to drought, due to their high resistance to xylem embolism (Robert et al., 2017; Lobo et al., 2018), their high capacity to withdraw water from the soil also under dry conditions by deep root systems and their ability to maintain water potential gradients along the soil-stem-crown continuum (Zweifel et al., 2007), as well as their xeromorphic leaf structure. Moreover, they are able to rapidly resume assimilation after periods of water deficiency (Kubiske and Abrams, 1993; Kuster et al., 2013).
Across Europe, a considerable number of studies found increasing mortality rates due to drought and heatwave events in Spain (Peñuelas et al., 2001; Martínez-Vilalta and Piñol, 2002), France (Bréda et al., 2006; Vennetier et al., 2007), Switzerland (Bigler et al., 2006), Poland (Siwecki and Ufnalski, 1998), Norway (Solberg et al., 2004), and Russia (Kauhanen et al., 2008; Ogibin and Demidova, 2009). However, opposite trends in mortality rates have also been reported. When comparing the period of 1900–1960 with 1960–2000, mortality rates of Norway spruce in Germany did not change, while for beech, mortality rates decreased by 17% (Pretzsch et al., 2014b). Similarly, annual mortality rates of four angiosperm species in Sweden were also found to be smaller in recent decades (1988–2013) than earlier (1912–1988, Hytteborn et al., 2017). Our current understanding on forest mortality trends is strongly influenced by the time span of available data. A long-term perspective is necessary to obtain a comprehensive picture of forest dynamics, especially when dealing with long-living organisms, such as trees.
It is also important to make a distinction between mortality rates of individual species and mortality rates of entire forest ecosystems. Studies investigating physiological mechanisms driving mortality usually focus on species-specific mortality, and often cover a relatively short time span, limited sample size and geographical range (e.g., Chao et al., 2008; Mitchell et al., 2013). In contrast, ecosystem-scale mortality rates are typically assessed across multi-year or multi-decadal time span, covering a large sample size and geographical range (e.g., Williams et al., 2013; Neumann et al., 2017; Huelsmann et al., 2018). The two perspectives complement each other and are essential for improving our current understanding and future predictions of tree and forest mortality across the globe.
The majority of ecosystem-scale studies assess mortality rates across one or two decades (e.g., Monserud and Sterba, 1999; Bradford and Bell, 2017; Pillet et al., 2018; Rogers et al., 2018) or after a specific disturbance event (e.g., Breshears et al., 2005; Michaelian et al., 2011; Brodrick and Asner, 2017; Reed et al., 2018). However, long-term mortality trends can add a unique perspective and significantly improve our understanding of forest dynamics. Long-term monitoring requires extensive effort, resources and commitment; hence, the data are extremely scarce. We are aware of only three studies reporting tree mortality rates of seven European species and spanning a time frame of at least 100 years (Pretzsch et al., 2014a,b; Hytteborn et al., 2017). An additional challenge is to assess mortality across broad climatic gradients and across a more diverse suite of species. For example, a recent geographically comprehensive study encompassed entire Europe, but included only two species and covered one decade (Neumann et al., 2017).
In this study, we address these knowledge gaps by compiling a dataset spanning ca. 120 years of inventory data for five dominant European tree species in Switzerland: Scots pine, Norway spruce, silver fir, European beech, sessile and common oak. Despite the fact that our study is geographically restricted to Switzerland, the large heterogeneity of the Swiss landscape enables an assessment of forest mortality under a wide range of environmental conditions. For example, the 276 sampled plots spanned a wide altitudinal (∼1800 m) and precipitation gradient (∼800 mm during growth season). We combined long-term inventory data of forest mortality, stand characteristics such as basal area and mean stand diameter at breast height (used as proxy for stand age), and high-resolution climate data to: (i) assess annual tree mortality rates and temporal trends of five dominant tree species across Switzerland over the last century, and (ii) identify the main drivers of mortality, out of a set of 12 predictors, including climate, stand characteristics and topography. We hereby focused on the impact of summer drought conditions on mortality, being aware that also other climatic factors, such as winter conditions, might have an impact on forest mortality at certain sites. We hypothesized that drought-induced mortality is species-specific with high drought-induced mortality rates in spruce and beech and lowest mortality rates in more drought-resistant oak and pine. Within species we expected different responses to drought along environmental gradients, due to different limiting factors that depend on site conditions (e.g., temperature vs. water availability). For instance, pine or spruce growing at dry (lowland) sites and showing already climate related declines should exhibit higher drought induced mortality rates compared to wet or high-altitude sites. On the other hand, it has been shown that, e.g., beech trees that are adapted to a generally lower water availability are more resistant to occurring droughts compared to beech trees growing under generally good water supply (Grossiord et al., 2014; Kunz et al., 2018). Therefore, we would expect a higher drought-induced mortality on wet sites compared to drier sites for species with a high adaptive capacity, e.g., beech and oak.
Materials and Methods
Forest Inventory Data
The study includes data from three long-term data sources of Swiss forest growth monitoring networks (Supplementary Table S1): Swiss Experimental Forest Management plots network (EFM), Swiss Long-Term Forest Ecosystem monitoring network (LWF), and the Swiss Nature Reserves Network (SNR). Within all three networks, the diameter at breast height (DBH, 1.3 m height) and the status (dead/alive) of all marked trees larger than a network-specific DBH threshold were recorded during repeated inventory campaigns (Supplementary Table S1). For consistency, only trees with DBH > 5 cm were included in this study.
The EFM plots were established between 1887 and the early 1900s, and currently comprise ca. 390 yield plots across Switzerland in order to study the development of forests in a changing environment under the influence of forest management (cf., Schütz and Zingg, 2010). Inventory intervals ranged from 3 to 22 years (mean 7 years). The LWF plots were established in 1995 and are part of the intensively monitored Level II plot network of the International Co-operative Programme on Assessment and Monitoring of Air Pollution Effects on Forests (ICP Forests) in Europe (Ferretti et al., 2010). The LWF plots are assessed every 3–6 years (mean 5 years) since 1995 according to a harmonized sampling protocol (Dobbertin and Neumann, 2010; Hug et al., 2011). The SNR network currently consists of 49 unmanaged forest reserves throughout Switzerland (Brang et al., 2006, 2011; Wunder et al., 2007). Repeated inventory campaigns of core plots within the nature reserves have been conducted every 5–20 years (mean 12 years) since the late 1940s. It needs to be noted that the permanent plots of the three monitoring networks were assessed at different intervals (3–22 years, on average 9 years), for different inventory years and time periods, and that the number of plots increased with time.
A total of 276 permanent plots within 95 forest sites (Supplementary Table S1 and Figure 1) were selected according to the following criteria: (1) at least two inventories per plot were available; (2) the plots were not affected – as far as documented – by major disturbances or management interventions; for the EFM, only plots with natural thinning (A grade according to IUFRO rules, Versuchsanstalten, 1902) or plots with low management intensity, where suppressed, fallen or dying stems are removed from below but without impact on crown layer (B grade), were included, so that the measured mortality rates should reflect as far as possible the natural mortality rates; (3) plots consisted of at least 20% of the basal area of one of the following species: Scots pine (Pinus sylvestris L.), Norway spruce [Picea abies (L.) H. Karst], silver fir (Abies alba Mill.), European beech (Fagus sylvatica L.) or oak [Quercus petraea (Matt.) Liebl., Quercus robur L.]. The two oak species were treated as one entity and are referred throughout the text as ‘oak’, due to their morphological similarities and tendency to hybridize, leading to their disputed taxonomical identity (Kissling, 1980; Aas, 1998; Muir et al., 2000).
Figure 1. Location of plots from EFM (red symbols), LWF (green symbols) and SNR (blue symbols) networks included in this study. In a number of locations, several plots are located in close proximity to each other. See Supplementary Material for maps showing locations of each species separately (Supplementary Figure S1).
The plots within the SNR and the EFM forest sites were considered as separate plots, despite their geographical proximity (on average 100 m apart, ranging from 3 to 423 m), as their stand structure and abiotic characteristics were very heterogeneous (cf. Brang et al., 2011). The individual plots covered an area of 0.02–3.47 ha (mean 0.5 ha). Out of the total of 276 plots, 175 were pure stands (basal area of a given species > 70%) and 101 mixed stands. Among the pure stands, 47 were dominated by pine, 23 by spruce, 10 by fir, 42 by beech, 50 by oak, and 3 by other species (for these plots, species with remaining 20–30% basal area were included in analyses). According to annual changes in tree numbers and basal area as described in Brang et al. (2011), 40 stands were classified to the juvenile forest phase, 17 stands to the young forest phase, 178 to the optimal phase and 41 to the old growth phase. The mean stand age of the 100 forest stands with age information was 130 years (ranging from 45 to 360 years). In total, the data series used in this study consisted of 2–17 inventories per plot (on average 9) from 1898 to 2013. Average data series length was 45 years and the longest series comprised 113 years from 1898 to 2011. In total, data of 247,529 trees were included.
All site-specific climate data (temperature, precipitation, drought stress indicators) for the time period from 1901 to 2012 were obtained with a daily resolution from Remund et al. (2016). Climate data from 1901 to 1930 were interpolated from CRU grid data “CRU TS Version 1.2” (Mitchell et al., 2004) based on the change factor method, and since 1931 from observational plots (MeteoSchweiz) based on the Shepard’s Gravity interpolation method (Zelenka et al., 1992).
For all study plots within one forest site identical climate data were used. All climate variables were calculated for the growth period from May to September. The minimal site water balance (SWBmin) during the growth period was used as drought indicator as it was found to be one of the best predictors of mortality in a comparison of various multivariate models for beech and spruce in Switzerland (Braun, 2016). In addition, SWBmin resulted in the best agreement with soil water content measured in 10–20 cm soil depth for a subset of LWF sites used in this study, relative to other drought indices (R = 0.48, data not shown). SWBmin was calculated according to Grier and Running (1977) as sum of the daily differences between precipitation sum and potential evapotranspiration with field capacity as starting value. The detailed methodology is described here: http://www.wsl.ch/staff/niklaus.zimmermann/programs/amls/swb.aml. Lower SWBmin, denotes higher drought stress. In this study, the minimal value for the growth season was applied.
Within each species, plots were categorized into two distinct groups, referred to as “ecoregions” (Table 1 and Supplementary Figure S1). The grouping was intended to emphasize the diverse mortality patterns in forests growing under contrasting climate conditions. Recent studies have shown that temperature and drought have divergent effects on high elevation vs. low elevation forests, as well as on temperature-limited vs. water-limited forests (Jolly et al., 2005; Lindner et al., 2010; Sarris et al., 2011). For Scots pine in the inner-Alps in Switzerland, mortality was found to be highest below 1000 m a.s.l., which was also related to drought (Rigling et al., 2013). We divided plots so that both ecoregions per species contained approximately equal number of plots. The resulting ecoregions encompassed low and high-altitude plots for pine (threshold at 1000 m) and spruce (1300 m). The grouping of pine and spruce by elevation resulted also in significantly different climate conditions with warmer temperatures in the lowlands, and for spruce also drier conditions, compared to the highlands (Table 1). For beech, oak and fir, which grow predominantly at lower elevations, we distinguished dry and wet plots based on the SWBmin (beech, fir SWBmin = 100 mm, oak SWBmin = 50 mm).
Table 1. Stand characteristics for 1960–2012 period summarizing mean values and the 0.05–0.95% quantile range (in parentheses) for plots within a given species-ecoregion.
Based on the individual tree data, annual mortality rates m were calculated at a population level per inventory period (interval) as follows:
where N0 and Nt are the numbers of living trees at the beginning and end of the interval, respectively, and t is the inventory interval in years (Sheil and May, 1996). Mortality rates m were only calculated for populations with N0 > 10. In populations with different mortality probabilities, mortality rates decline with the length of the inventory interval, because the fraction of trees with a higher mortality probability declines faster than the fraction of trees with a low mortality probability. Therefore, mortality rates calculated from diverse interval lengths were compared for subpopulations for which homogeneous mortality probabilities may be assumed (Sheil and May, 1996). To account for the different interval length in our data set, we followed the approach adopted by Rohner et al. (2012) and calculated mortality rates for subpopulations, i.e., for each tree species and for three DBH-classes (small, medium, large) separately, assuming homogeneous mortality probability for these subpopulations. DBH-classes for each plot and each time period were defined so that all three DBH-classes contained the same number of trees.
To allow for comparison among plots, we standardized the mortality rate m per plot and inventory period by dividing it by the mean mortality rate of the plot during the total observation period:
where rm is the normalized mortality (%), m is the mortality rate as in Equation (1) and i is the inventory period.
We distinguished two time periods: the full time period covered by the inventory data set of 1898 to 2013, and a shorter time period including only the recent decades since 1960. This was done because (i) many plots (especially within the LWF and SNR networks) were established later than 1960, and (ii) climate in Switzerland started to change around 1960 (Figure 2). Hereby, inventory data were available for the period 1898–2013 and climate data for the period 1901–2012. Therefore, for analyses based only on inventory data, we considered the whole length of the inventory time series, whereas for analyses including also climate data we used inventory data only for 1901–2012 to match with the time-period of climate data availability.
Figure 2. Mean trends of normalized annual temperature (A), precipitation (B), and SWBmin (C) during growth period for the study sites during the years 1901–2012. Yearly data were normalized by their long-term average. Black lines with gray areas indicate the mean ± SD. For a better trend visualization loess smoother was added.
To assess the temporal trend of mortality, generalized linear mixed effects models (GLMM) for proportions with annual mortality rate as response variable were used, i.e., logistic regression models assuming binomial distribution (Zuur et al., 2009). Random effects were defined for the intercept with study plots nested in forest sites (see “Study sites” section) as grouping factor. We used two different model approaches: (1) We fitted one model per species and ecoregion with year at the beginning of the inventory period and DBH-class as fixed effects (Equation 3). (2) In order to detect whether mortality within the three DBH-classes followed different (e.g., counteracting) trends, we fitted additional GLMMs for the three DBH-classes separately (Equation 4).
Where Yijk is the annual mortality rate per observation i, study site j and study plot k, α is the intercept, Year is the start year of inventory, and DBH-class is a factor with three levels (small, medium, large). Fixed effects parameters are given as β, with β1 indicating the temporal trend of the mortality time series and β2 indicating the additive effect of DBH-class. aj and ajjk are site and plot random effects, assuming normal distribution with mean = 0 and σj2 and σjjk 2, respectively.
where Yijk,DBH-class is the annual mortality rate per observation i, study site j and study plot k of a given DBH-class (small, medium or large). All other parameters are the same as described for Equation 3.
Multivariate relationships between explanatory variables and temporal changes in mortality rates were analyzed using GLMMs as described for Equations 3 and 4, but considering a pool of 8 initial predictors (stand characteristics, plot properties, and climate variables) and 4 interactions as fixed effects (Equation 5). All possible combinations of models were calculated per species and ecoregion. To avoid oversaturation of a model, a maximum of 5 (fir) and 7 (other species) explanatory variables were included simultaneously (resulting in 310 and 560 potential models, respectively). Preliminary analysis indicated that an additional error structure to account for plot spatial autocorrelation did not improve model performance and was not incorporated into the final model.
where x1…xn are the included explanatory variables as described below. All other parameters are the same as described for Equation 3.
We considered the following initial predictors: Basal area (BA, m2 ha-1) was recorded at the beginning of each inventory period and indicated the competitive state in the stand. As stand age was not available for each plot, we used the mean stand DBH as a proxy for stand age (mDBH, cm). Plot-wise regressions of mDBH against stand age for 100 plots with available stand age information resulted in a median adj. R2 of 0.96. Further, we included the DBH-classes and topographic properties, such as altitude (m), slope (degree) and aspect (as factor with four levels) as fixed effects. Mean temperature (Temp, °C) and SWBmin (mm) as drought indicator during the growth period were calculated for all inventory periods and then related to the long-term average (ΔTemp, ΔSWBmin). All numerical variables (BA, mDBH, slope, altitude, ΔTemp, ΔSWBmin) were standardized. We tested for interactions between BA and/or mDBH and the climate variables (ΔTemp and ΔSWBmin). Only plots with at least two inventory periods available were included in the GLMMs (n = 224). To account for collinearity of explanatory variables only models with a maximum correlation coefficient between explanatory variables smaller than 0.7 were considered. Models were ranked according to their corrected AIC, AICc (Hurvich and Tsay, 1989). In case of several comparable models (Akaike weight of the best model < 0.9) an average model of all models with a ΔAICc (calculated from the model with lowest AICc) < 4 was calculated (Burnham and Anderson, 2002). Model residuals were visually checked for heterogeneity vs. fitted values and included variables.
All statistical analyses were performed in R statistical software (R Development Core Team, 2011) with the function glmer (package ‘lme4’) for GLMMs and dredge and model.avg (package ‘MuMin’) for model selection and model averaging, respectively.
Overview of Annual Mortality Rates
Average annual mortality rate for all species across the ∼120 years was 1.5%, and a considerable variation in mortality rates across species and ecoregions (dry and wet, high and low altitudes) was observed (Table 1 and Figure 3). Oak and fir had highest mean annual mortality rates of 2 and 1.8%, respectively, pine and beech had 1.3%, and spruce had the lowest mortality rate of 1% (Table 2). However, mortality rates of individual plots ranged from 0% to over 13%, with majority of plots (95%) averaging below 6%.
Figure 3. Annual mortality rates (%) for the five studied species (A, pine; B, spruce; C, fir; D, beech; E, oak) and ecoregions across the entire study period (1898–2013). Each data point represents mean annual mortality rate ± SE calculated at 10-year intervals and averaged across all measured plots. Numbers in the gray band (‘n obs/plots’) indicates the number of observations included/number of plots monitored during the respective inventory interval. Inventory intervals ranged from 3 to 22 years and were carried out in different years for different plots. Because of the different inventory intervals and dates, the temporal trajectories can differ from results shown in Table 2.
Table 2. Mean annual mortality rates per species ± standard error (SE) calculated for different study periods.
Driving Factors of Annual Mortality Rates
Out of an initial set of 8 predictor variables and 4 interactions, GLMMs were calculated for each species and ecoregion to identify the drivers of temporal variation in tree mortality (Table 3), which are presented in detail in the following.
Table 3. Coefficients of the averaged generalized mixed effects models (GLMM) with mortality rate as response variable and stand characteristics, climate and topography as explanatory variables, as well as interactions of stand characteristics with climate variables (indicated by ‘:’).
Basal area (BA) was the most important and consistent predictor of annual mortality rates for all species and ecoregions. Increasing annual mortality rates were correlated with higher BA (Table 3 and Figure 4). The strongest effect was observed for fir and oak growing on dry sites and beech from wet sites. BA was not significant for only two species-ecoregions, pine at low altitudes and spruce at high altitudes. In the latter case, however, interaction of BA with climate variables was significant (Table 3).
Figure 4. Curves illustrating predicted response of annual mortality rates to stand and climate variables as derived from the GLMMs (Table 3). (A–E) Correspond to the five species studied. (A) pine, (B) spruce, (C) fir, (D) beech, and (E) oak. Response curves were calculated by holding all predictors in the model constant at their mean, except for the selected displayed variable. Curves represent a mean response over the three DBH-classes, and four aspect levels (if significant in the model). ΔSWBmin: deviation from the long-term mean SWBmin (minimum Site Water Balance) during the growth season, ΔTemp: deviation from the long-term mean temperature during the growth season.
Mean DBH (mDBH, proxy for stand age) had a significant effect in six out of the ten species-ecoregion combinations and both negative and positive relationships were observed (Table 3). Increased mortality rates were associated with higher mDBH in pine from low altitudes, spruce from high altitudes, and fir and oak both from dry regions, while for spruce at low altitudes, and oak and beech in wet regions, increased mortality rates were related to lower mDBH. Annual mortality rates varied considerably between the three DBH-classes. Overall, smallest trees had the largest annual mortality rates, while largest trees had, in most cases, lower mortality rates (Table 3).
Changes in the minimum site water balance (ΔSWBmin) and in temperature (ΔTemp) had variable effects on annual mortality rates, and were significant in four species-ecoregions for ΔSWBmin and for ΔTemp (Table 3 and Figure 4). Increasing temperatures were associated with increasing mortality only in oak from wet regions, while mortality decreased in pine at high altitudes, fir and oak both in dry regions. The strongest negative impact of ΔSWBmin on annual mortality rates was recorded for pine in the lowlands and a weaker effect for spruce in lowlands, thus mortality increased under drier conditions. In contrast, for pine and spruce at high altitudes no significant impact of ΔSWBmin on annual mortality rates was found. Fir mortality was not related to ΔSWBmin in neither of the two ecoregions. For beech and oak, mortality rates were positively associated with ΔSWBmin, indicating an increase in mortality under wet conditions. For beech at dry sites the relationship was only marginally significant, however, ΔSWBmin was included in 24 of the best 25 models. Although the relationship of ΔSWBmin and mortality was not significant in the averaged model for beech at wet sites, ΔSWBmin was included in the best ranked models either as significant single parameter or in a significant interaction (Supplementary Table S2). For both, beech and oak, mortality rates at a given ΔSWBmin appeared to be slightly higher on wet sites compared to dry sites (Figure 4), for beech especially in denser stands (cf. Supplementary Figure S3), but differences were generally small. Larger differences were observed for ΔTemp. For beech, only wet sites showed a negative temperature effect on mortality in interaction with stand age, whereas mortality was not related to ΔTemp on dry sites. For oak, mortality was increased at higher ΔTemp on wet sites, whereas the sign was reversed at dry sites (Figure 4).
The impact of ΔSWBmin and/or ΔTemp on mortality was dependent on mDBH for all species, particularly at low and dry sites (Supplementary Figures S2, S3). In general, the effect of ΔSWBmin was stronger in old stands, except for pine at low altitudes with a reversed interaction. The effect of ΔTemp negatively interacted with mDBH, indicating a greater vulnerability of young stands relative to old stands to increased temperature. At high-altitude and wet sites, the magnitude of the climatic effect (ΔSWBmin and/or ΔTemp) on mortality was stronger for dense forests stands (i.e., higher stand BA) (Table 3 and Supplementary Figure S3). This was consistent for all species.
Altitude and Site Topography
Altitude, slope and aspect were only occasionally related to mortality (Table 3) with no clear pattern. For instance, for beech growing on wet sites, mortality rates were negatively correlated with altitude, while for pine from low altitude regions, the pattern was opposite. The steepness of the slope was significantly related to mortality for beech from wet and fir from dry regions. Aspect was negatively related to mortality rates for beech from both dry and wet regions and positively in pine from high altitudes. For all other species-ecoregions the relationships of mortality with topography were not significant.
Temporal Trends of Annual Mortality Rates
Temporal trends of annual mortality rates varied considerably across species and ecoregions throughout the entire study period (starting in ∼1900s, 1910s, 1920s, and 1930s depending on the species and ecoregion). Generalized mixed effect models indicated that over the last ∼120 years, pine and spruce from high altitudes and fir from dry regions exhibited the highest increase of annual mortality rates by 1.6, 1.0, and 1.3% per year, respectively (Table 4). For other species and ecoregions either no significant or very small shifts were observed (<±0.5%). For pine in lowlands and beech in dry regions, time series started later than 1960 and were therefore not included in this analysis.
Table 4. Results of the GLMMs with annual mortality rate as response variable, with first year of the inventory period as fixed effect, and plot nested in forest site (see Materials and Methods, Equations 3 and 4) as random effect.
Throughout 1960–2013, pine in lowlands showed the highest increase of mortality rates by 2.3% per year, followed by spruce in both ecoregions, with an increase by 1.8 and 1.6% per year at low and high altitudes, respectively (Table 4). Fir from dry regions, beech from wet regions and oak from both dry and wet regions showed decline in mortality rates with most pronounced shift in oak from dry regions (-4.2% per year). In the remaining three species-ecoregions the change was not significant.
Over the last ∼120 years, mortality of trees within the large and small DBH-classes exhibited opposite trends: mortality of large trees increased in seven species-ecoregions, while mortality of small trees decreased, increased or showed no significant trend (Table 4). Since 1960, this tendency became even more pronounced with an increase of mortality in large trees for six species-ecoregions, and a decrease in small trees for five species-ecoregions or no significant changes. Exceptions were fir at dry sites, where all DBH-classes showed increases in mortality, and in oak where a general decrease of mortality across the DBH-classes was observed. These opposing tendencies also emerged when comparing mortality rates of DBH-classes between two periods before and after 1960, although changes were not always significant due to large data heterogeneity (Figure 5).
Figure 5. Annual mortality rate by species and DBH-class for time periods before and after 1960. Asterisks indicate statistically significant differences in mortality rates per DBH-class and species between the two time periods (∗∗∗P < 0.001, ∗∗P < 0.01, and ∗P < <0.05, unpaired t-test on normally distributed data). Number of plots included is the same as in Table 2.
Despite notable differences between species in temporal trends of annual mortality (Figure 3 and Table 4) and heterogeneous data (different number of plots, inventory dates and intervals between the inventories), the five species displayed similar fluctuations in normalized mortality rates (Figure 6 and Supplementary Table S3). Normalized mortality rates express annual mortality rate relative to the long-term average annual mortality rate. Although annual mortality rates changed only mildly during the last ∼120 years (Table 4), it is noteworthy that in the most recent three decades (1980–2010) the cross-species average mortality rate did not drop below the 100-year average. In addition, the frequency of mortality peaks in the last three decades was also higher than earlier. The magnitude of the peaks, however, decreased probably due to the higher plot number in recent years. There were also distinct differences in the temporal mortality trends (cf. Table 4). For instance, the temporal trajectory of beech mortality correlated well with all four other species (Figure 6 and Supplementary Table S3), while trajectory of oak mortality was well associated with only two other species: beech and pine.
Figure 6. Normalized mortality rates averaged across all species (bold line) and individual species (thin lines) over the entire study period. The lines are smoothed with a spline function. Normalized mortality rate is annual mortality rate divided by average mortality rate across the entire study period. The misalignment of trajectories may be partly due to different inventory dates and intervals. The number of plots is indicated by the gray polygon (see also Figure 3).
In this study, we assessed mortality rates of five dominant species in Swiss forests throughout the last century (1898–2013) and examined factors driving mortality change. The complex study design with different points in time of forest inventories as well as differing length of the inventory periods within and across the three monitoring networks challenges statistical analysis and generates limitations on data interpretation. Nevertheless, some very clear patterns emerged from the one century-long data set. Since 1960, mortality rates of pine in lowlands and of spruce increased significantly, whereas those of oak trees decreased (Table 4). Fir and beech showed only minor or not significant changes in mortality rates. Stand properties, especially BA, were the most important factors influencing variability of mortality rates, while climate impact (temporal changes in minimum site water balance and temperature) was inconsistent across species and ecoregions, and depended on interactions with mDBH (i.e., an indicator of stand age) and BA (indicating competition) (Table 3). Thus, we could not conclusively confirm our initial hypothesis, that oak and pine would be most tolerant, and beech and spruce would be most vulnerable to drought, at least not in absence of other stand-related factors. The impact of drought and heat on mortality appeared rather to result from a combination of species and site effects, and being further modulated by competition and stand age. The long-term data presented here highlight the complex nature of the processes shaping forest dynamics. Throughout one century, unmanaged forests gradually alter their structure, demography and in some cases also their species composition. Consequently, structural and demographical changes in combination with climate affected mortality rates in our study, implying that a variety of factors need to be accounted for in mortality studies at the ecosystem scale. Below, we discuss the annual mortality rates, the effect of competition and climate variability on mortality, and its impact on the future development of forests in Switzerland.
Annual Mortality Rates
Average annual mortality rate through the entire study period was 1.5% but varied considerably across species and plots (Figure 3 and Tables 1, 2). On average, mortality rates of various tree species in Europe range between 0.1 and 2.9 %, with a median of 0.5% (Supplementary Table S5), which is lower compared to mortality rates in this study. Beyond the overall inherent, high variability of mortality rates across forests, this may be due to several methodological and conceptual reasons: (1) the 5 cm DBH threshold applied here was lower than in other studies, and this might have increased the average mortality rate of the investigated populations, since mortality tends to be higher in small trees (Rohner et al., 2012; Rigling et al., 2013); (2) The data set included a high proportion of young stands, for which mortality rates are usually higher compared to forests in the optimum phase (Oliver and Larson, 1996); (3) The long study period of over 100 years may have increased the probability of mortality events, although fluctuations of mortality rates across different decadal time periods were small (Supplementary Table S4); (4) Our study encompassed a large number of plots across broad altitudinal and climate gradients, including also plots at species distribution margins with exceptionally high mortality, rising the overall average mortality rates; (5) Finally, the absence of management in our study sites might explain the relatively high mortality rates, especially in small and medium size tree classes. Thus, a direct comparison of mortality rates across different studies should be interpreted with caution. These and other factors challenge a direct comparison of absolute mortality rates across different studies, while relative changes of mortality rates as well as their driving factors should be comparable across studies.
Drivers of Mortality
Mortality rates were best explained by a combination of stand properties and climate variables, whereas altitude and topography were less significant. In concordance with other studies (Condés and del Río, 2015; Pillet et al., 2018) our results showed that when climate and competition considered together a larger portion of variation in mortality rates can be explained, as when considered separately.
Effects of Stand Characteristics on Mortality Rates
Basal area, and thus competition, was the strongest predictor variable of mortality rates. For all species and ecoregions a consistent effect of increasing mortality with higher BA was found (Table 3 and Figure 4), as high stand density leads to increased competition for resources, which consequently may result in higher mortality (Giuggiola et al., 2013; Ruiz-Benito et al., 2013; Condés and del Río, 2015; Bradford and Bell, 2017; Young et al., 2017).
In all species, mortality rates were largest among the smallest diameter trees, while the largest trees had lowest mortality rates, as reported also in many other studies (Monserud and Sterba, 1999; Mantgem and Stephenson, 2007; van Mantgem et al., 2009; Peng et al., 2011; Neumann et al., 2017) and could be explained by smaller trees being outcompeted by larger trees (Coomes and Allen, 2007) or by trees of the same size class. The impact of stand age (implied by mean DBH) on mortality rates was significant for most species and ecoregions, but with variable sign and strength, which might reflect differences in stand development of the included plots.
Effects of Climate on Mortality Rates
Climate variables, i.e., changes in temperature and SWBmin, had contrasting influence on mortality rates, depending on species and ecoregion. For spruce and pine, mortality of lowland forests increased with increasing dryness (negative ΔSWBmin). Contrary to our initial hypothesis, this effect was most conspicuous for pine, for which SWBmin was the overall strongest predictor of annual mortality rates (Table 3 and Figure 4). Increased mortality of pine due to drought became, however, an increasingly observed pattern in the dry valleys of the Valais, in Grisons and the Churer Rheintal in Switzerland (Dobbertin et al., 2005; Bigler et al., 2006; Rigling et al., 2006, 2013; Wohlgemuth and Rigling, 2014). Although pine can recover after a few incidents of drought (Eilmann et al., 2010; Eilmann and Rigling, 2012), several consecutive dry years can lead to a reduction in vitality and associated higher susceptibility of trees to pests, pathogens and parasites (Bigler et al., 2006; Rigling et al., 2010). Although Scots pine is generally drought-tolerant, but - growing already on dry sites – an increase in frequency and severity of drought events may force it beyond its physiological limits. Given that pine has been introduced by humans outside its native range and was promoted until 1950 (Gimmi et al., 2010), it is now being outcompeted by indigenous or invasive species, or species which are better adapted to dry conditions, like pubescent oak (Rigling et al., 2013). Spruce in the lowlands has also been found to be susceptible to drought (Vanoni et al., 2016a,b) and is expected to be substituted by more drought-tolerant species in the future (Hanewinkel et al., 2013). However, we did not find a very strong impact of drought on spruce mortality. Instead, BA and mDBH were more influential. This might be due to the fact that we excluded plots that were substantially affected by storm damages or bark beetle attacks. These plots might have been also experienced severe drought beforehand since the intensity of the attack is often amplified by preceding drought events (Meier et al., 2014; Kolb et al., 2016). Thus, secondary factors, such as bark beetles, affect spruce mortality rates often to a larger extent than drought itself. Additionally, in contrast to the lowland pine plots, which were restricted to very dry regions, lowland spruce plots exhibited on average a relatively high soil water availability (Table 1). Contrary to the lowlands, mortality rates of pine and spruce at higher altitudes were not related to drought, but increased under low temperatures, probably due to direct low-temperature limitation (Körner, 1998, 2003; Vanoni et al., 2016a) combined with indirect temperature stress such as frost damage, winter desiccation, or low-temperature photo-inhibition (Barbeito et al., 2012). Similar patterns were also observed for Scots pine growing at higher altitudes in Spain (Ruiz-Benito et al., 2013).
For silver fir trees, drought did not show any significant influence on annual mortality (Table 3). This pattern confirms other studies from the Churer Rheintal in Switzerland, where the detrimental impact of drought on growth of silver fir trees is mitigated by deep soils with higher water holding capacity and northern exposition (Wohlgemuth and Rigling, 2014). On drier sites, the strong impact of low temperatures on mortality rates, could presumably be explained by the sensitivity of silver fir to frost (Lebourgeois et al., 2010). In summary, silver fir in Switzerland seems to be able to cope with the current level of experienced droughts and is still within its physiological boundaries, as also observed in other parts of Europe (George et al., 2015).
Increasing mortality with higher water availability was found for beech and, even more noticeably, for oak (Table 3 and Figure 4). This pattern is counter-intuitive, especially for beech on dry sites, although it was also observed in other studies (Rohner et al., 2012; Nothdurft, 2013; Ruiz-Benito et al., 2013), and might be related to other concurrent abiotic conditions. Presumably, higher water availability in moist years can lead to an increase of total leaf area and subsequent overshadowing of short, young trees limiting their access to light and water, consequently resulting in their increased mortality. Beech is often considered to be drought sensitive (Gessler et al., 2007) and, in general, performs better on wetter sites (Aranda et al., 2000; Gessler et al., 2007). A recent study also found that beech growth appears to be very location specific, with growth decline at lower altitudes and growth increase at higher altitudes (Dulamsuren et al., 2017). In concordance with this, we observed higher beech mortality at lower altitudes compared to higher altitudes.
It has been observed that beech (Vitasse et al., 2014) and also oak (Jensen and Hansen, 2010; Grossiord et al., 2014) trees can adapt very well to dry conditions, possibly due to their deep rooting penetration, efficient stomatal control (Nothdurft, 2013), and high degree of evolutionary adaptability (Roloff and Grundmann, 2008). Thus, higher mortality rates per given ΔSWBmin at wet sites compared to dry sites for beech and oak partly confirm the hypothesis that drought resistance of trees increases with lower site water availability (Kunz et al., 2018). This might be due to the high degree of genetic adaptability. However, differences related to SWB were small and might also be associated with generally higher turn-over at sites with higher water availability. The temperature effect on mortality of the two species clearly differed between dry and wet sites of the two species, indicating that trees that are adapted to dry and warm conditions might be able to cope better with occurring heat waves as projected for the future. For oak species, the impact of pathogens related to wet conditions and frost (Führer, 1998; Gaertig et al., 2005) but also to drought (Wood et al., 2018) is an additional relevant factor driving mortality (Haavik et al., 2015).
Effects of Stand Characteristics and Climate on Mortality Rates
The inclusion of an interaction term of BA or mDBH (i.e., stand age) with climate variables improved the model parsimony for pine, spruce, fir and oak, indicating that the effect of climatic variables varies across stand structures and developmental stages (Table 3). These effects were more pronounced on dry and lowland sites compared to wet and highland sites (Supplementary Figures S2, S3). Our results indicate that the impact of drought and/or temperature on tree mortality increased with competition (i.e., higher BA), which has also been reported for many tree species in Spain (Ruiz-Benito et al., 2013). In addition, pine has been found to perform better after stand density reduction on xeric sites in Switzerland (Giuggiola et al., 2013, 2018), further indicating that competition is an important factor in driving mortality and demography (Kohler et al., 2010; Sohn et al., 2016). For age, the pattern reversed along the temperature gradient and young stands were generally more strongly affected by increased temperatures compared to old stands. This might be due to the fact that young stands are usually denser and prone to self-thinning. Our results indicate that heat and/or drought may trigger the self-thinning process in young stands.
Temporal Trends in Annual Mortality
Temporal trends in annual mortality rates varied depending on the data grouping factor (e.g., species, ecoregions, DBH-classes, Figures 3, 5, 6 and Table 4) partly due to methodological reasons (i.e., data averaging and integration procedures). For example, discrepancies between temporal trajectories in GLMMs (Table 4) and figures (Figures 3, 5) originated from the fact that models used actual dates of inventories and interval lengths and accounted for plots as grouping factor, while figures show mortality estimates averaged per decade or even larger time periods and/or per species and ecoregions. Nevertheless, some clear patterns emerged. Species-specific trends of annual mortality rates across the ∼120 years were modest, but became larger during the second half of the century (Figure 3 and Table 4). Spruce and pine exhibited increasing mortality rates, while fir on dry sites, oak on dry and wet sites, and, to a lesser degree, beech on wet sites showed decreasing mortality rates. Thus, the normalized mortality rates (Figure 6) averaged over all species showed a slight, but consistent, increasing trend in mortality across all species, but with a large uncertainty at the beginning of the study period due to low plot number. The steady increase in the normalized mortality rates since 1900, rather than an obvious break point at the start of pronounced climatic changes around 1960 (Figure 2), suggests that the increase in mortality is more strongly driven by gradually changing stand parameters (such as density and age, as shown in the previous paragraphs) and to a lesser extent by climate. However, it is also noteworthy that fluctuations of normalized mortality rates of all species are remarkably similar, especially since the 1960s (Figure 6 and Supplementary Table S3). Thus, mortality peaks were for all species relatively (given the heterogeneity of the data structure) synchronous every 10–20 years, especially during the recent decades. Additionally, since the 1980s, annual mortality rates were consistently, although only slightly, above the long-term average. Temporal disturbance patterns strongly synchronized across landscapes have been also observed in temperate forests in Europe during 1986–2016 and were related to preceding drought and storm events (Senf and Seidl, 2018). Thus, forest stand properties, such as BA, might drive the overall long-term trend of forest mortality, while climate modulates these trends with an assumed increasing importance during the recent decades.
Our results also demonstrate that trends in annual mortality rates differed depending on tree size and even run in opposing directions. Mortality rates of small trees tended to decrease or did not change, while mortality rates of large trees generally increased over time (Figure 5 and Table 4). This might be a stand aging effect, since observed forests aged for up to 120 years throughout the study period and run through different stand developmental stages. Moreover, recent studies in temperate forests indicate that susceptibility to drought increases with tree age or size (Carrer and Urbinati, 2004; Lloret et al., 2011; Ding et al., 2017), probably due to higher risk of hydraulic failure (Bennett et al., 2015; McDowell and Allen, 2015) or higher costs involved in growing new functional xylem after the drought has passed (Trugman et al., 2018). As a consequence, an increased mortality of large trees may have allowed for higher survival of young trees. However, also contradictory observations have been made, showing higher susceptibility to drought of young trees (Mantgem and Stephenson, 2007; Colangelo et al., 2017). More research is needed to clarify this context dependency (Ding et al., 2017).
Contrarily to our initial hypothesis, we found a significant drought-induced mortality only in pine at low altitudes and, to a lesser extent, in spruce from low altitudes (Table 3). It appears, that drought-tolerant pine, growing already on dry sites, was severely impacted by drought events, but that fir, beech and oak, growing mainly on soils with a good water holding capacity, seemed not yet to be affected by the occurrence of drought in Switzerland. They appear to be quite well adapted to the current conditions within their ecological niche. This pattern is contrary to many studies reporting rising mortality rates due to drought and heat (Breshears et al., 2005; Gitlin et al., 2006; Allen et al., 2010, 2015; Anderegg et al., 2013; Gustafson and Sturtevant, 2013; Holmgren et al., 2013; Reichstein et al., 2013; Bradford and Bell, 2017; Chen et al., 2018; McDowell et al., 2018). One reason for this discrepancy may stem from the fact that many studies reported mortality after a specific drought or heatwave event and/or were restricted to relatively short time intervals of several years, including a large climatic disturbance event (for example, Martínez-Vilalta and Piñol, 2002; Dobbertin, 2005; Landmann and Dreyer, 2006; Ogibin and Demidova, 2009). The inventory intervals of five or more years also may dilute the climate signal on mortality (Dobbertin et al., 2005; Huelsmann et al., 2018). Nevertheless, our results show that mortality rates of five dominant tree species in Switzerland increased only slightly over the last ∼120 years, which could mainly be related to changes in stand structure. This suggests that Swiss forests have been resilient to recent climates change so far (with the exception of some hotspots) and that instead of an abrupt transition of forests, changes in species composition might occur more gradually and subtly. Nevertheless, the amplified effect of drought and heat under competition might indicate potentially strong changes in the demographic structure and species composition of forests in the future (Ruiz-Benito et al., 2013).
It is predicted that the future will bring longer and more intense droughts in many regions across Europe (Fischer and Schär, 2010; Dai, 2013; CH2014-Impacts, 2014). The resilience of trees to climatic changes requires further investigation, for instance, how quickly and effectively trees recover from episodes of drought and heatwaves. Some studies suggest that drought tolerant species may shift to sites that are becoming more arid, as these species are already acclimated to drier conditions. This shift may eliminate species that are already at their ecophysiological limits (Lévesque et al., 2013, 2014; Rigling et al., 2013). Moreover, species specific responses to insect attacks, pest and pathogens, may shape the structure of future forests. For example, in Germany and Austria, an increase in mortality risk is predicted especially for spruce at low altitudes, while beech and oak may be more robust to climate change (Lexer et al., 2002; Nothdurft, 2013; Ding et al., 2017), which could be supported by the observations from Switzerland reported in this study.
The datasets generated for this study are available on request to the corresponding author.
AR, SE, and AZ conceived the study. SE analyzed the data with support from BR, AB, AKB, NKR, and AR. SE and KZ wrote the manuscript with contribution from all other co-authors.
This study was funded by the Swiss Federal Office for the Environment FOEN and the Swiss Federal Institute for Forest, Snow and Landscape Research WSL as part of the research program ‘Forests and Climate Change.’ NKR received funding by the German Research Foundation through its Emmy Noether Programme (RU 1657/2-1).
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
We would like to thank Jan Remund from meteotest for providing the climate data. We also would like to thank the networks EFM, LWF, and NKF for providing the data, and all the field teams, who did the extensive monitoring work since the last 100 years. This study is dedicated to Matthias Dobbertin (†2012) who initiated this project.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2019.00307/full#supplementary-material
Adam, H. D., Zeppel, M. J. B., Anderegg, W. R. L., Hartmann, H., Landhäusser, S. M., Tissue David, T., et al. (2017). A multi-species synthesis of physiological mechanisms in drought-induced tree mortality. Nat. Ecol. Evol. 1, 1285–1291. doi: 10.1038/s41559-017-0248-x
Allen, C. D., Breshears, D. D., and McDowell, N. G. (2015). On underestimation of global vulnerability to tree mortality and forest die-off from hotter drought in the Anthropocene. Ecosphere 6:art129. doi: 10.1890/ES15-00203.1
Allen, C. D., Macalady, A. K., Chenchouni, H., Bachelet, D., McDowell, N., Vennetier, M., et al. (2010). A global overview of drought and heat-induced tree mortality reveals emerging climate change risks for forests. For. Ecol. Manag. 259, 660–684. doi: 10.1016/j.foreco.2009.09.001
Anderegg, W. R. L., Flint, A., Huang, C.-Y., Flint, L., Berry, J. A., Davis, F. W., et al. (2015). Tree mortality predicted from drought-induced vascular damage. Nat. Geosci. 8, 367–371. doi: 10.1038/ngeo2400
Anderegg, W. R. L., Kane, J. M., and Anderegg, L. D. L. (2013). Consequences of widespread tree mortality triggered by drought and temperature stress. Nat. Clim. Chang. 3, 30–36. doi: 10.1038/nclimate1635
Anderegg, W. R. L., Klein, T., Bartlett, M. K., Sack, L., Pellegrini, A. F. A., Choat, B., et al. (2016). Meta-analysis reveals that hydraulic traits explain cross-species patterns of drought-induced tree mortality across the globe. Proc. Natl. Acad. Sci. U.S.A. 113, 5024–5029. doi: 10.1073/pnas.1525678113
Aranda, I., Gil, L., and Pardos, J. A. (2000). Water relations and gas exchange in Fagus sylvatica L. and Quercus petraea (Mattuschka) Liebl. in a mixed stand at their southern limit of distribution in Europe. Trees 14, 344–352. doi: 10.1007/s004680050229
Barbeito, I., Dawes, M. A., Rixen, C., Senn, J., and Beri, P. (2012). Factors driving mortality and growth at treeline: a 30-year experiment of 92 000 conifers. Ecology 93, 389–401. doi: 10.1890/11-0384.1
Bigler, C., Braeker, O. U., Bugmann, H., Dobbertin, M., and Rigling, A. (2006). Drought as inciting mortality factor in Scots pine stands of the Valais, Switzerland. Ecosystems 9, 330–343. doi: 10.1007/s10021-005-0126-2
Bradford, J. B., and Bell, D. M. (2017). A window of opportunity for climate-change adaptation: easing tree mortality by reducing forest basal area. Front. Ecol. Environ. 15, 11–17. doi: 10.1002/fee.1445
Brando, P. M., Balch, J. K., Nepstad, D. C., Morton, D. C., Putz, F. E., Coe, M. T., et al. (2014). Abrupt increases in Amazonian tree mortality due to drought–fire interactions. Proc. Natl. Acad. Sci. U.S.A. 111, 6347–6352. doi: 10.1073/pnas.1305499111
Braun, S. (2016). Die Eignung Verschiedener Trockenheitsindices zur Abschätzung des Trockenheitsrisikos für den Wald und ihre standortspezifische Entwicklung bis 2050/2100: Schlussbericht des Projektes im Forschungsprogramm “Wald und Klimawandel”. Zürich: ETH Zurich.
Breda, N., Huc, R., Granier, A., and Dreyer, E. (2006). Temperate forest trees and stands under severe drought: a review of ecophysiological responses, adaptation processes and long-term consequences. Ann. For. Sci. 63, 625–644. doi: 10.1051/forest:2006042
Bréda, N., Huc, R., Granier, A., and Dreyer, E. (2006). Temperate forest trees and stands under severe drought: a review of ecophysiological responses, adaptation processes and long-term consequences. Ann. For. Sci. 63, 625–644. doi: 10.1051/forest:2006042
Breshears, D. D., Cobb, N. S., Rich, P. M., Price, K. P., Allen, C. D., Balice, R. G., et al. (2005). Regional vegetation die-off in response to global-change-type drought. Proc. Natl. Acad. Sci. U.S.A. 102, 15144–15148. doi: 10.1073/pnas.0505734102
Cailleret, M., Jansen, S., Robert, E. M. R., Desoto, L., Aakala, T., Antos, J. A., et al. (2017). A synthesis of radial growth patterns preceding tree mortality. Glob. Chang. Biol. 23, 1675–1690. doi: 10.1111/gcb.13535
Chao, K.-J., Phillips, O. L., Gloor, E., Monteagudo, A., Torres-Lezama, A., and Martínez, R. V. (2008). Growth and wood density predict tree mortality in Amazon forests. J. Ecol. 96, 281–292. doi: 10.1111/j.1365-2745.2007.01343.x
Charru, M., Seynave, I., Morneau, F., and Bontemps, J. D. (2010). Recent changes in forest productivity: an analysis of national forest inventory data for common beech (Fagus sylvatica L.) in north-eastern France. For. Ecol. Manag. 260, 864–874. doi: 10.1016/j.foreco.2010.06.005
Chen, L., Huang, J.-G., Dawson, A., Zhai, L., Stadt, K. J., Comeau, P. G., et al. (2018). Contributions of insects and droughts to growth decline of trembling aspen mixed boreal forest of western Canada. Glob. Chang. Biol. 24, 655–667. doi: 10.1111/gcb.13855
Cocozza, C., de Miguel, M., Psidova, E., Ditmarova, L., Marino, S., Maiuro, L., et al. (2016). Variation in ecophysiological traits and drought tolerance of beech (Fagus sylvatica L.) seedlings from different populations. Front. Plant Sci. 7:886. doi: 10.3389/fpls.2016.00886
Colangelo, M., Camarero, J. J., Borghetti, M., Gazol, A., Gentilesca, T., and Ripullone, F. (2017). Size matters a lot: drought-affected Italian oaks are smaller and show lower growth prior to tree death. Front. Plant Sci. 8:135. doi: 10.3389/fpls.2017.00135
Condés, S., and del Río, M. (2015). Climate modifies tree interactions in terms of basal area growth and mortality in monospecific and mixed Fagus sylvatica and Pinus sylvestris forests. Eur. J. For. Res. 134, 1095–1108. doi: 10.1007/s10342-015-0912-0
Ding, H., Pretzsch, H., Schuetze, G., and Roetzer, T. (2017). Size-dependence of tree growth response to drought for Norway spruce and European beech individuals in monospecific and mixed-species stands. Plant Biol. 19, 709–719. doi: 10.1111/plb.12596
Dobbertin, M., Mayer, P., Wohlgemuth, T., Feldmeyer-Christe, E., Graf, U., Zimmermann, N. E., et al. (2005). The decline of Pinus sylvestris L. forests in the swiss Rhone Valley - a result of drought stress? Phyton Ann. Rei Bot. A 45, 153–156.
Dobbertin, M., and Neumann, M. (2010). “Part V: Tree growth,” in Manual on Methods and for Harmonized Sampling, Assessment, Monitoring and Analysis of the Effects of Air Pollution on Forests, ed. UNECE ICP-Forests Programme Co-ordinating Centre (Hamburg: UNECE ICP-Forests Programme Co-ordinating Centre), 1–29.
Dulamsuren, C., Hauck, M., Kopp, G., Ruff, M., and Leuschner, C. (2017). European beech responds to climate change with growth decline at lower, and growth increase at higher elevations in the center of its distribution range (SW Germany). Trees 31, 673–686. doi: 10.1007/s00468-016-1499-x
Eilmann, B., Buchmann, N., Siegwolf, R., Saurer, M., Cherubini, P., and Rigling, A. (2010). Fast response of Scots pine to improved water availability reflected in tree-ring width and δ13C. Plant Cell Environ. 33, 1351–1360. doi: 10.1111/j.1365-3040.2010.02153.x
Ferretti, M., Fischer, R., Mues, V., Granke, O., and Lorenz, M. (2010). “Basic design principles for the ICP forests monitoring networks. Manual part II,” in Manual on Methods and Criteria for Harmonized Sampling, Assessment, Monitoring and Analysis of the Effects of Air Pollution on Forests, ed. UNECE ICP Forests Programme Co-ordinating Centre (Hamburg: UNECE ICP Forests Programme Coordinating Centre).
Führer, E. (1998). “Oak decline in central Europe: a synopsis of hypotheses,” in Proceedings of Population Dynamics, Impacts, and Integrated Management of Forest Defoliating Insects (General Technical Report NE-247), eds M. L. McManus and A. M. Liebhold (Radnor, PA: USDA Forest Service), 7–24.
Gaertig, T., von Wilpert, K., and Seemann, D. (2005). Differentialdiagnostische Untersuchungen zu Eichenschäden in Baden-Württemberg. Freiburg: Forstliche Versuchs- und Forschungsanstalt Baden-Württemberg, 61.
George, J.-P., Schueler, S., Karanitsch-Ackerl, S., Mayer, K., Klumpp, R. T., and Grabner, M. (2015). Inter- and intra-specific variation in drought sensitivity in Abies spec. and its relation to wood density and growth traits. Agric. For. Meteorol. 214–215, 430–443. doi: 10.1016/j.agrformet.2015.08.268
Gessler, A., Keitel, C., Kreuzwieser, J., Matyssek, R., Seiler, W., and Rennenberg, H. (2007). Potential risks for European beech (Fagus sylvatica L.) in a changing climate. Trees 21, 1–11. doi: 10.1007/s00468-006-0107-x
Gimmi, U., Wohlgemuth, T., Rigling, A., Hoffmann, C. W., and Bürgi, M. (2010). Land-use and climate change effects in forest compositional trajectories in a dry Central-Alpine valley. Ann. For. Sci. 67:701. doi: 10.1051/forest/2010025
Gitlin, A. R., Sthultz, C. M., Bowker, M. A., Stumpf, S., Paxton, K. L., Kennedy, K., et al. (2006). Mortality gradients within and among dominant plant populations as barometers of ecosystem change during extreme drought. Conserv. Biol. 20, 1477–1486. doi: 10.1111/j.1523-1739.2006.00424.x
Giuggiola, A., Bugmann, H., Zingg, A., Dobbertin, M., and Rigling, A. (2013). Reduction of stand density increases drought resistance in xeric Scots pine forests. For. Ecol. Manag. 310, 827–835. doi: 10.1016/j.foreco.2013.09.030
Giuggiola, A., Kuster, T. M., and Saha, S. (2010). Drought-induced mortality of Scots pines at the southern limits of its distribution in Europe: causes and consequences. iForest Biogeosci. For. 3, 95–97. doi: 10.3832/ifor0542-003
Giuggiola, A., Zweifel, R., Feichtinger, L., Vollenweider, P., Bugmann, H., Haeni, M., et al. (2018). Competition for water in a xeric forest ecosystem – effects of understory removal on soil micro-climate, growth and physiology of dominant Scots pine trees. For. Ecol. Manag. 409, 241–249. doi: 10.1016/j.foreco.2017.11.002
Grossiord, C., Gessler, A., Granier, A., Pollastrini, M., and Bussotti, F. (2014). Interspecific competition influences the response of oak transpiration to increasing drought stress in a mixed mediterranean forest. For. Ecol. Manag. 318, 54–61. doi: 10.1016/j.foreco.2014.01.004
Haavik, L. J., Billings, S. A., Guldin, J. M., and Stephen, F. M. (2015). Emergent insects, pathogens and drought shape changing patterns in oak decline in North America and Europe. For. Ecol. Manag. 354, 190–205. doi: 10.1016/j.foreco.2015.06.019
Hanewinkel, M., Cullmann, D. A., Schelhaas, M. J., Nabuurs, G. J., and Zimmermann, N. E. (2013). Climate change may cause severe loss in the economic value of European forest land. Nat. Clim. Chang. 3, 203–207. doi: 10.1038/nclimate1687
Hartmann, H., Moura, C. F., Anderegg, W. R. L., Ruehr, N. K., Salmon, Y., Allen, C. D., et al. (2018). Research frontiers for improving our understanding of drought-induced tree and forest mortality. New Phytol. 218, 15–28. doi: 10.1111/nph.15048
Huelsmann, L., Bugmann, H., Cailleret, M., and Brang, P. (2018). How to kill a tree: empirical mortality models for 18 species and their performance in a dynamic forest model. Ecol. Appl. 28, 522–540. doi: 10.1002/eap.1668
Hug, C., Dobbertin, M., Nussbaumer, C., and Stettler, Y. (2011). Aufnahmeanleitung für die Brusthöhenumfang- und Höheninventur auf den Flächen der Langfristige Waldökosystem-Forschbung (LWF). Birmensdorf: Swiss Federal Institute for Forest, Snow and Landscape Research (WSL).
Hytteborn, H., Svensson, B. M., Kempe, K., Press, A., and Rydin, H. (2017). Century-long tree population dynamics in a deciduous forest stand in central Sweden. J. Veg. Sci. 28, 1057–1069. doi: 10.1111/jvs.12556
Jolly, W. M., Dobbertin, M., Zimmermann, N. E., and Reichstein, M. (2005). Divergent vegetation growth responses to the 2003 heat wave in the Swiss Alps. Geophys. Res. Lett. 32:L18409. doi: 10.1029/2005gl023252
Jump, A. S., Hunt, J. M., and Penuelas, J. (2006). Rapid climate change-related growth decline at the southern range edge of Fagus sylvatica. Glob. Chang. Biol. 12, 2163–2174. doi: 10.1111/j.1365-2486.2006.01250.x
Kauhanen, H., Wallenius, T., Kuuluvainen, T., Aakala, T., and Mikkola, K. (2008). “Extensive mortality of spruce forests in Arkhangelsk Region: satellite image analysis,” in Poster Presentation At: International Conference “Adaptation of Forests and Forest Management to Changing Climate with Emphasis on Forest Health: A Review of Science, Policies, and Practices”, Umeå.
Kohler, M., Sohn, J., Nägele, G., and Bauhus, J. (2010). Can drought tolerance of Norway spruce (Picea abies (L.) Karst.) be increased through thinning? Eur. J. For. Res. 129, 1109–1118. doi: 10.1007/s10342-010-0397-9
Kolb, T. E., Fettig, C. J., Ayres, M. P., Bentz, B. J., Hicke, J. A., Mathiasen, R., et al. (2016). Observed and anticipated impacts of drought on forest insects and diseases in the United States. For. Ecol. Manag. 380, 321–334. doi: 10.1016/j.foreco.2016.04.051
Kubiske, M. E., and Abrams, M. D. (1993). Stomatal and nonstomatal limitations of photosynthesis in 19 temperate tree species on contrasting sites during wet and dry years. Plant Cell Environ. 16, 123–129. doi: 10.1111/j.1365-3040.1996.tb02070.x
Kunz, J., Löffler, G., and Bauhus, J. (2018). Minor European broadleaved tree species are more drought-tolerant than Fagus sylvatica but not more tolerant than Quercus petraea. For. Ecol. Manag. 414, 15–27. doi: 10.1016/j.foreco.2018.02.016
Kuster, T. M., Arend, M., Bleuler, P., Günthardt-Goerg, M. S., and Schulin, R. (2013). Water regime and growth of young oak stands subjected to air-warming and drought on two different forest soils in a model ecosystem experiment. Plant Biol. 15, 138–147. doi: 10.1111/j.1438-8677.2011.00552.x
Landmann, G., and Dreyer, E. (2006). Impacts of drought and heat on forest. Synthesis of available knowledge, with emphasis on the 2003 event in Europe. Ann. For. Sci. 3, 567–652. doi: 10.1051/forest:2006062
Lebourgeois, F., Rathgeber, C. B. K., and Ulrich, E. (2010). Sensitivity of French temperate coniferous forests to climate variability and extreme events (Abies alba, Picea abies and Pinus sylvestris). J. Veg. Sci. 21, 364–376. doi: 10.1111/j.1654-1103.2009.01148.x
Leuschner, C., Backes, K., Hertel, D., Schipka, F., Schmitt, U., Terborg, O., et al. (2001). Drought response at leaf, stem and fine root levels of competitive Fagus sylvatica L. and Quercus petraea (Matt.) Liebl. trees in dry and wet years. For. Ecol. Manag. 149, 33–46. doi: 10.1016/S0378-1127(00)00543-0
Lévesque, M., Rigling, A., Bugmann, H., Weber, P., and Brang, P. (2014). Growth response of five co-occurring conifers to drought across a wide climatic gradient in Central Europe. Agric. For. Meteorol. 197, 1–12. doi: 10.1016/j.agrformet.2014.06.001
Lévesque, M., Saurer, M., Siegwolf, R., Eilmann, B., Brang, P., Bugmann, H., et al. (2013). Drought response of five conifer species under contrasting water availability suggests high vulnerability of Norway spruce and European larch. Glob. Chang. Biol. 19, 3184–3199. doi: 10.1111/gcb.12268
Lexer, M. J., Hönninger, K., Scheifinger, H., Matulla, C., Groll, N., Kromp-Kolb, H., et al. (2002). The sensitivity of Austrian forests to scenarios of climatic change: a large-scale risk assessment based on a modified gap model and forest inventory data. For. Ecol. Manag. 162, 53–72. doi: 10.1016/S0378-1127(02)00050-6
Lindner, M., Maroschek, M., Netherer, S., Kremer, A., Barbati, A., Garcia-Gonzalo, J., et al. (2010). Climate change impacts, adaptive capacity, and vulnerability of European forest ecosystems. For. Ecol. Manag. 259, 698–709. doi: 10.1016/j.foreco.2009.09.023
Lines, E. R., Coomes, D. A., and Purves, D. W. (2010). Influences of forest structure, climate and species composition on tree mortality across the eastern US. PLoS One 5:e13212. doi: 10.1371/journal.pone.0013212
Lloret, F., Keeling, E. G., and Sala, A. (2011). Components of tree resilience: effects of successive low-growth episodes in old ponderosa pine forests. Oikos 120, 1909–1920. doi: 10.1111/j.1600-0706.2011.19372.x
Lobo, A., Torres-Ruiz, J. M., Burlett, R., Lemaire, C., Parise, C., Francioni, C., et al. (2018). Assessing inter- and intraspecific variability of xylem vulnerability to embolism in oaks. For. Ecol. Manag. 424, 53–61. doi: 10.1016/j.foreco.2018.04.031
Martínez-Vilalta, J., and Piñol, J. (2002). Drought-induced mortality and hydraulic architecture in pine populations of the NE Iberian Peninsula. For. Ecol. Manag. 161, 247–256. doi: 10.1016/S0378-1127(01)00495-9
McDowell, N., Allen Craig, D., Anderson-Teixeira, K., Brando, P., Brienen, R., Chambers, J., et al. (2018). Drivers and mechanisms of tree mortality in moist tropical forests. New Phytol. 219, 851–869. doi: 10.1111/nph.15027
Michaelian, M., Hogg, E. H., Hall, R. J., and Arsenault, E. (2011). Massive mortality of aspen following severe drought along the southern edge of the Canadian boreal forest. Glob. Chang. Biol. 17, 2084–2094. doi: 10.1111/j.1365-2486.2010.02357.x
Mitchell, P. J., O’Grady, A. P., Tissue, D. T., White, D. A., Ottenschlaeger, M. L., and Pinkard, E. A. (2013). Drought response strategies define the relative contributions of hydraulic dysfunction and carbohydrate depletion during tree mortality. New Phytol. 197, 862–872. doi: 10.1111/nph.12064
Mitchell, T. D., Carter, T. R., Jones, P. D., Hulme, M., and New, M. (2004). A Comprehensive Set of Climate Scenarios for Europe and the Globe: The Observed Record (1900–2000) and 16 Scenarios (2000–2100). Norwich: University of East Anglia, 30.
Nothdurft, A. (2013). Spatio-temporal prediction of tree mortality based on long-term sample plots, climate change scenarios and parameter frailty modeling. For. Ecol. Manag. 291, 43–54. doi: 10.1016/j.foreco.2012.11.028
Ogibin, B. N., and Demidova, N. A. (2009). “Successional dynamics of old-growth spruce forests in the watersheds of the rivers Northern Dvina-Pinega in the Arkhangelsk Region,” in . Northern Coniferous Forests–Tools Through Research for the Sustainable Use of Forests in the Barents Region, eds H. Kauhanen, V. Neshataev, and M. Vuopio (Helsinki: Finnish Forest Research Institute), 108–113.
Pan, Y., Birdsey, R. A., Fang, J., Houghton, R., Kauppi, P. E., Kurz, W. A., et al. (2011). A large and persistent carbon sink in the world’s forests. Science 333, 988–993. doi: 10.1126/science.1201609
Peng, C., Ma, Z., Lei, X., Zhu, Q., Chen, H., Wang, W., et al. (2011). A drought-induced pervasive increase in tree mortality across Canada’s boreal forests. Nat. Clim. Chang. 1, 467–471. doi: 10.1038/nclimate1293
Pillet, M., Joetzjer, E., Belmin, C., Chave, J., Ciais, P., Dourdain, A., et al. (2018). Disentangling competitive vs. climatic drivers of tropical forest mortality. J. Ecol. 106, 1165–1179. doi: 10.1111/1365-2745.12876
Piovesan, G., Biondi, F., Di Fillipo, A., and Maugeri, M. (2008). Drought-driven growth reduction in old beech (Fagus sylvatica L.) forests of the central Apennines, Italy. Glob. Chang. Biol. 14, 1265–1281. doi: 10.1111/j.1365-2486.2008.01570.x
Pretzsch, H., Biber, P., Schütze, G., and Bielak, K. (2014a). Changes of forest stand dynamics in Europe. Facts from long-term observational plots and their relevance for forest ecology and management. For. Ecol. Manag. 316, 65–77. doi: 10.1016/j.foreco.2013.07.050
Pretzsch, H., Schütze, G., and Uhl, E. (2013). Resistance of European tree species to drought stress in mixed versus pure forests: evidence of stress release by inter-specific facilitation. Plant Biol. 15, 483–495. doi: 10.1111/j.1438-8677.2012.00670.x
Reed, D. E., Ewers, B. E., Pendall, E., Frank, J., and Kelly, R. (2018). Bark beetle-induced tree mortality alters stand energy budgets due to water budget changes. Theor. Appl. Climatol. 131, 153–165. doi: 10.1007/s00704-016-1965-9
Rigling, A., Bigler, C., Eilmann, B., Feldmeyer-Christe, E., Gimmi, U., Ginzler, C., et al. (2013). Driving factors of a vegetation shift from Scots pine to pubescent oak in dry Alpine forests. Glob. Chang. Biol. 19, 229–240. doi: 10.1111/gcb.12038
Rigling, A., Dobbertin, M., Bürgi, M., Feldmeier-Christe, E., Gimmi, U., Ginzler, C., et al. (2006). “Baumartenwechsel in den walliser waldföhrenwäldern,” in Wald und Klimawandel. Forum für Wissen, ed. T. Wohlgemuth (Birmensdorf: Eidgenössische Forschungsanstalt WSL), 71.
Robert, E. M. R., Mencuccini, M., and Martinez-Vilalta, J. (2017). “The anatomy and functioning of the xylem in oaks,” in Tree Physiology. Oaks Physiological Ecology. Exploring the Functional Diversity of Genus Quercus L, eds E. Gil-Pelegrín, J. J. Peguero-Pina, and D. Sancho-Knapik (Cham: Springer International Publishing), 261–302.
Rogers, B. M., Solvik, K., Hogg, E. H., Ju, J., Masek, J. G., Michaelian, M., et al. (2018). Detecting early warning signals of tree mortality in boreal North America using multiscale satellite data. Glob. Chang. Biol. 24, 2284–2304. doi: 10.1111/gcb.14107
Rohner, B., Bigler, C., Wunder, J., Brang, P., and Bugmann, H. (2012). Fifty years of natural succession in Swiss forest reserves: changes in stand structure and mortality rates of oak and beech. J. Veg. Sci. 23, 892–905. doi: 10.1111/j.1654-1103.2012.01408.x
Ruiz-Benito, P., Lines, E. R., Gómez-Aparicio, L., Zavala, M. A., and Coomes, D. A. (2013). Patterns and drivers of tree mortality in iberian forests: climatic effects are modified by competition. PLoS One 8:e56843. doi: 10.1371/journal.pone.0056843
Sarris, D., Christodoulakis, D., and Körner, C. (2011). Impact of recent climatic change on growth of low elevation eastern Mediterranean forest trees. Clim. Chang. 106, 203–223. doi: 10.1007/s10584-010-9901-y
Schütz, J. P., and Zingg, A. (2010). Improving estimations of maximal stand density by combining Reineke’s size-density rule and yield level, using the example of spruce (Picea abies (L.) Karst.) and European Beech (Fagus sylvatica L.). Ann. For. Sci. 67:507. doi: 10.1051/forest/2010009
Senf, C., and Seidl, R. (2018). Natural disturbances are spatially diverse but temporally synchronized across temperate forest landscapes in Europe. Glob. Chang. Biol. 24, 1201–1211. doi: 10.1111/gcb.13897
Sohn, J. A., Hartig, F., Kohler, M., Huss, J., and Bauhus, J. (2016). Heavy and frequent thinning promotes drought adaptation in Pinus sylvestris forests. Ecol. Appl. 26, 2190–2205. doi: 10.1002/eap.1373
Solberg, S., Andreassen, K., Clarke, N., Tørseth, K., Tveito, O. E., Strand, G. H., et al. (2004). The possible influence of nitrogen and acid deposition on forest growth in Norway. For. Ecol. Manage. 192, 241–249. doi: 10.1016/j.foreco.2004.01.036
Stephens, S. L., Collins, B. M., Fettig, C. J., Finney, M. A., Hoffman, C. M., Knapp, E. E., et al. (2018). Drought, tree mortality, and wildfire in forests adapted to frequent fire. BioScience 68, 77–88. doi: 10.1093/biosci/bix146
Stojnic, S., Suchocka, M., Benito-Garzón, M., Torres-Ruiz, J. M., Cochard, H., Bolte, A., et al. (2018). Variation in xylem vulnerability to embolism in European beech from geographically marginal populations. Tree Physiol. 38, 173–185. doi: 10.1093/treephys/tpx128
Trugman, A. T., Detto, M., Bartlett, M. K., Medvigy, D., Anderegg, W. R. L., Schwalm, C., et al. (2018). Tree carbon allocation explains forest drought-kill and recovery patterns. Ecol. Lett. 21, 1552–1560. doi: 10.1111/ele.13136
van Mantgem, P. J., Stephenson, N. L., Byrne, J. C., Daniels, L. D., Franklin, J. F., Fulé, P. Z., et al. (2009). Widespread increase of tree mortality rates in the western United States. Science 323, 521–524. doi: 10.1126/science.1165000
Vanoni, M., Bugmann, H., Noetzli, M., and Bigler, C. (2016a). Drought and frost contribute to abrupt growth decreases before tree mortality in nine temperate tree species. For. Ecol. Manag. 382, 51–63. doi: 10.1016/j.foreco.2016.10.001
Vanoni, M., Bugmann, H., Noetzli, M., and Bigler, C. (2016b). Quantifying the effects of drought on abrupt growth decreases of major tree species in Switzerland. Ecol. Evol. 6, 3555–3570. doi: 10.1002/ece3.2146
Vennetier, M., Vila, B., Liang, E., Guibal, F., Taahbet, A., and Gadbin-Henry, C. (2007). Impact of climate change on pine forest productivity and on the shift of a bioclimatic limit in a Mediterranean area. Option Médit. Séries A 75, 189–197.
Versuchsanstalten, V. D. F. (1902). Beratungen der vom Vereine deutscher Forstlicher Versuchsanstalten eingesetzten Kommission zur Feststellung des neuen Arbeitsplanes für Durchforstungs-und Lichtungsversuche. Allg. For. Jadgztg 78, 180–184.
Vitali, V., Büntgen, U., and Bauhus, J. (2017). Silver fir and Douglas fir are more tolerant to extreme droughts than Norway spruce in south-western Germany. Glob. Chang. Biol. 23, 5108–5119. doi: 10.1111/gcb.13774
Vitasse, Y., Lenz, A., Kollas, C., Randin, C., Hoch, G., and Körner, C. (2014). Genetic vs. non-genetic responses of leaf morphology and growth to elevation in temperate tree species. Funt. Ecol. 28, 243–252. doi: 10.1111/1365-2435.12161
Wermelinger, B., Rigling, A., Schneider-Mathis, D., and Dobbertin, M. (2008). Infestation preferences of bark and wood boring insects in declining Scots pine (Pinus sylvestris) forests in the Swiss Rhone valley. Ecol. Entomol. 33, 239–249. doi: 10.1111/j.1365-2311.2007.00960.x
Williams, A. P., Allen, C. D., Macalady, A. K., Griffin, D., Woodhouse, C. A., Meko, D. M., et al. (2013). Temperature as a potent driver of regional forest drought stress and tree mortality. Nat. Clim. Chang. 3, 292–297. doi: 10.1038/nclimate1693
Wohlgemuth, T., and Rigling, A. (2014). Kurz-und Langfristige Auswirkungen des Klimas auf die Wälder im Churer Rheintal: Schlussbericht Projekt Bündner Wald im Klimawandel. Birmensdorf: Eidg. Forschungsanstalt für Wald, Schnee und Landschaft WSL.
Wood, J. D., Knapp, B. O., Muzika, R.-M., Stambaugh, M. C., and Gu, L. (2018). The importance of drought–pathogen interactions in driving oak mortality events in the Ozark Border Region. Environ. Res. Lett. 13:015004. doi: 10.1088/1748-9326/aa94fa
Wunder, J., Reineking, B., Matter, J.-F., Bigler, C., and Bugmann, H. (2007). Predicting tree death for Fagus sylvatica and Abies alba using permanent plot data. J. Veg. Sci. 18, 525–534. doi: 10.1111/j.1654-1103.2007.tb02567.x
Young, D. J. N., Stevens, J. T., Earles, J. M., Moore, J., Ellis, A., Jirka, A. L., et al. (2017). Long-term climate and competition explain forest mortality patterns under extreme drought. Ecol. Lett. 20, 78–86. doi: 10.1111/ele.12711
Zang, C., Hartl-Meier, C., Dittmar, C., Rothe, A., and Menzel, A. (2014). Patterns of drought tolerance in major European temperate forest trees: climatic drivers and levels of variability. Glob. Chang. Biol. 20, 3767–3779. doi: 10.1111/gcb.12637
Keywords: drought, competition, stand basal area, climate change, tree size, mortality
Citation: Etzold S, Ziemińska K, Rohner B, Bottero A, Bose AK, Ruehr NK, Zingg A and Rigling A (2019) One Century of Forest Monitoring Data in Switzerland Reveals Species- and Site-Specific Trends of Climate-Induced Tree Mortality. Front. Plant Sci. 10:307. doi: 10.3389/fpls.2019.00307
Received: 12 November 2018; Accepted: 26 February 2019;
Published: 22 March 2019.
Edited by:Sebastian Leuzinger, Auckland University of Technology, New Zealand
Reviewed by:Martina Pollastrini, Università degli Studi di Firenze, Italy
Cyrille Barthélémy Karl Rathgeber, INRA Centre Nancy-Lorraine, France
Henrik Hartmann, Max-Planck-Institut für Biogeochemie, Germany
Copyright © 2019 Etzold, Ziemińska, Rohner, Bottero, Bose, Ruehr, Zingg and Rigling. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Sophia Etzold, email@example.com