Introduction

Severe drought events in Central Europe in 1976, 2003, and 2015 triggered multiple studies on the effects of episodic drought on the growth and mortality of forest tree species (Allen et al. 2015; Bréda et al. 2006; Ciais et al. 2005). The findings suggest that tree species cultivated at or beyond the border of their natural range, such as Norway spruce (Picea abies [L.] Karst) and European larch (Larix decidua Mill), in Central Europe show severe growth reduction and mortality (Kölling et al. 2009; Lévesque et al. 2013) during extreme drought events. To mitigate the effects of drought on tree productivity and survival, silviculture practices aim to select better acclimated species and provenances (Atzmon et al. 2004; Arend et al. 2011; Zang et al. 2011). Scots pine (Pinus silvestris L.) and sessile oak (Quercus petraea L.), for instance, are less susceptible to drought (Walentowski et al. 2007; Zang et al. 2011, 2012) than Norway spruce and becoming more suitable for forestry in Central Europe under climate change scenarios that predict future warm and dry conditions. Possible silvicultural practices in view of climate change include down regulating stand density (D’Amato et al. 2013; Sohn et al. 2016), modified thinning practices (Gebhardt et al. 2014; Pretzsch et al. 2018; Rodríguez-Calcerrada et al. 2011), and the promotion of mixed tree species plantings (Ammer 2017). The latter’s efficacy, however, has yet to be assessed for drought mitigation (Grossiord 2018).

Most current knowledge on tree responses to drought is derived from the analyses of episodic drought events like those in 1976, 2003, and 2015. However, the effects of extended drought periods on tree growth, as expected under future climate scenarios, are still unknown. It is currently thought that the ability of trees to acclimate to drought is underestimated (Lapenis et al. 2005; Reich et al. 2016). Forests may acclimate to extended drought by physiological, morphological, and allometric adjustment at the tree level (Aasamaa et al. 2004; Pretzsch et al. 2013; Schuldt et al. 2016), and by density reduction, structural, and species compositional changes at the stand level (Lapenis et al. 2005).

Here, we analysed and compared the growth responses of Norway spruce (Picea abies [L.] Karst.) and European beech (Fagus sylvatica [L.]) to natural episodic and experimentally extended drought in mature monospecific and mixed-species stands of Norway spruce and European beech in the Kranzberg Forest. This study utilized the throughfall exclusion experiment KROOF in the Kranzberg Forest (Pretzsch et al. 2014; Tomasella et al. 2018; Hesse et al. 2019) and additional long-term tree measurements nearby (Pretzsch et al. 1998). To better understand the long-term effects of drought on tree growth in intra- and inter-specific environments, we concentrated on the following questions and hypotheses:

Q1: How do species react to natural drought events (represented by the years 2003 and 2015) in intra- versus inter-specific environments?

H1: The growth adjustments of Norway spruce and European beech do not differ and are equal in intra-specific and inter-specific environments.

Q2: How do species respond to extended (5-year-long) experimentally induced drought? What drives adjustments in growth with a focus on intra- versus inter-specific environments?

H2: The growth of Norway spruce is equal to European beech; intra-specific responses do not differ from inter-specific responses; and all trees in a stand react similarly.

Q3: How does the extended (5-year-long) experimentally induced drought affect the tree mortality?

H3: Tree mortality does not differ between the treatment and control plots.

We also further discuss the ecological and practical silvicultural implications of growth responses to episodic and extended drought.

Materials and methods

Description of the study sites

Kranzberg Forest (longitude: 11° 39′ 42″ E, latitude: 48° 25′ 12″ N, altitude 490 m a.s.l) is located in Southern Germany, approximately 35 km Northeast of Munich. Average annual precipitation is 750–800 mm yr−1 with 460–500 mm during the growing season (May–September). The average annual air temperature is 7.8 °C and 13.8 °C on a seasonal basis. At the site, monospecific and mixed-species stands of Norway spruce and European beech stock grow on luvisol originating from loess over Tertiary sediments that provide a high nutrient and water supply (Göttlein et al. 2012; Pretzsch et al. 1998). Depending on soil depth, the water holding capacity for plant available water ranges between 17 and 28% of volumetric soil water content, while soil pHH2O varied between 4.1 and 5.1.

We characterized the water supply for each year by calculating the index of de Martonne (1926) (M = precipitation/(temperature + 10)) on the basis of the precipitation (in mm) and temperature (in °C) for the whole year (My) and for the growing season from April to September (Mgs). Because of its minimal data requirement, this index has been widely used in the recent studies to describe the drought conditions or aridity for a given region (Rötzer et al. 2012; Pretzsch et al. 2013; Quan et al. 2013). The water supply for plant growth improves with increasing M index. Within our study, Mgs varied between 12 and 24 and My between 30 and 65.

Throughfall exclusion experiment and control plots

We established 12 experimental plots in Kranzberg, i.e., 6 throughfall exclusion (TE) plots and 6 control plots (CO). Plot sizes varied between 110 and 200 m2. Summed over all plots, the total area was 868 m2 and 862 m2 for the CO and TE plots, respectively.

Before starting the throughfall exclusion experiment, soil and root trenching was performed in spring 2010. Soil was trenched to about 1 m deep and 15 cm wide and lined with a heavy-duty plastic tarp, impermeable to water and root growth, and refilled with the original soil material (Pretzsch et al. 2016). At about 1 m depth, a dense clay layer of tertiary sediments limits further downward-rooting (Häberle et al. 2015). In the six TE plots, roofs were installed about 3 m from the ground, completely underneath the stand canopy, to exclude all forms of precipitation. Roofs were first closed in 2014.

Roofs closed automatically in response to precipitation, and only stayed closed during precipitation events to prevent micro-meteorological and physiological effects (Pretzsch et al. 2014). Because the aim of the experiment was to induce summer drought, the roofs were kept open in the winter months. This resulted in small annual precipitation amounts for the throughfall exclusion plots in the years 2014–2018. The winter precipitation amounts for the 5 years of the experimental drought were below 150 mm (Fig. 1).

Fig. 1
figure 1

Temperature (left), precipitation (middle), and Martonne index (right) since 1998 by years (black line) and growing seasons (May–September, grey line) for the Kranzberg Forest site. The thin lines represent to conditions on the plots with throughfall exclusion from 2014 to 2018 (TE plots). The years 2003 and 2015 are indicated by broken vertical lines and represent years with extremely low water availability

Due to the natural drought in 2015 a bark beetle infestation was observed across the entire Kranzberg Forest. Therefore, starting within the year 2015, bark beetle damage was confined through annually spraying the spruce crowns and stem surfaces with the contact insecticide Karate Forst liquid using the canopy crane.

Stand water was variable in the study years 1998–2018. Extreme dry years in 2003 and 2015 had significantly lower precipitation amounts compared to the rather moist years of 2001–2002 and 2005–2013. Accordingly, the Martonne index varied from 30 (2003) to 65 (2001) for the whole year and from 15 (2003) to 25 (2005) for individual growing seasons. Meteorological data were acquired from a nearby forest weather station “Freising”, which is part of the Bavarian Environmental Monitoring System (LWF 2017). For further information about the Kranzberg Forest, see Göttlein et al. (2012); Häberle et al. (2012), and for more details about the KROOF experiment, see Pretzsch et al. (2014, 2018).

Dendrometric survey

A full survey of the Kranzberg Forest experimental plots in 2016 determined Norway spruce was 65 and European beech was 85 years old. Mean and dominant tree sizes were similar between the plots. The tallest trees (as used for calculating height of the dominant trees by Assmann and Franz 1963) had heights of 34.3 m (spruce) and 33.0 m (beech) indicating optimal growing conditions, i.e., site indexes of O40 according to the yield table of Assmann and Franz (1963, 1965) for Norway spruce and I. site class according to Schober (1975) for European beech.

The quadratic mean stem diameters at breast height were 27.1–36.4 cm, with mean heights of 27.2–36.4 m. Dominant tree diameters measured 41.4–44.9 cm. The stem diameters were the lowest in the monospecific stands; in the beech by 20% lower than in spruce. Stem diameters were the highest in mixed-species stands; with beech again by about 20% lower than in spruce. The tree heights were similar in monospecific and mixed species stands; on average, beech is by 5 m lower than spruce. Collectively, there were 639–926 trees per hectare with a stand basal area of 54.0–60.1 m2 ha−1, standing stem volume of 802–981 m3 ha−1, and a mean periodic volume growth (1998–2016) of 19.4–26.3 m3 ha−1 yr−1. The lower values of the given ranges for tree number, stand basal area, standing volume, and volume growth of the monospecific beech stands, the upper values of the monospecific spruce stands, and the mixed species stands lie in between (for more stand information, see Pretzsch et al. 2014, 2018).

We utilized two data sources to evaluate tree diameter.

Since 1998, all trees of the Kranzberg Forest site were equipped with permanent diameter tapes with Vernier scales for circumferential recording to a 1 mm resolution (UMS, Germany). These data, excluding those from the throughfall exclusion experiment, were used to analyse the natural episodic drought in 2003 and 2015 to answer question Q1. Thus, diameter and circumferential stem growth at breast height were recorded for 268 spruce and 141 beech trees for 2003 and for 214 spruce and 108 beech trees for 2015 (Table 1).

Table 1 Characteristics of Norway spruce and European beech in response to episodic droughts in 2003 and 2015

To analyse species response to extended drought (Question Q2), another 51 Norway spruce and 51 European beech were equipped with girth tapes and first measured in 2011. Half of the trees were under the throughfall exclusion roofs to mimic extreme summer drought conditions from 2014 to 2018, the other half served as controls (Pretzsch et al. 2016) (Table 2).

Table 2 Characteristics of the Norway spruce and European beech trees included in the throughfall exclusion experiment from 2011 to 2018 (throughfall exclusion from 2014 to 2018)

To compare the mortality of Norway spruce and European beech under episodic and extended drought (Question 3), we utilized both data sets, the long-term records from 1998–2018 (episodic droughts) and the time series from 2014 to 2018 (experimentally extended drought) (Table 4).

Based on the stem diameter, di at the beginning of each year i and the annual circumferential growth ici, equal to the annual diameter growth idi = ici/π within the year i, we calculated the annual basal area growth \({\text{iba}}_{i} = \pi /4 \times ({\text{d}}_{i} + i{\text{d}}_{i} )^{2} - \pi /4 \times {\text{d}}_{i}^{2} = \pi /4 \times (2 \times {\text{d}}_{i} \times i{\text{d}}_{i} + i{\text{d}}_{i} 2)\) (Assmann 1961, p. 52).

Methods

Quantification of intra- and inter-specific environments

Species composition within each tree’s environment was quantified via an algorithm that counted the species identity of its six nearest neighbours (Fig. 2). The neighbours were chosen irrespectively of their size and social position; due to its advanced development state, the stand mainly consisted of codominant and dominant trees. Based on the results, we assigned each tree to one of four groups, ss = spruce surrounded by spruce, sb = spruce surrounded by beech, bb = beech surrounded by beech, and bs = beech surrounded by spruce.

Fig. 2
figure 2

Method of characterizing a central tree’s (tree in the middle of the respective hexagons) intra- or inter-specific environment; this approach was developed as an algorithm for automatic sorting into the groups ss, sb, bs, and bb. a In mature stands, trees grow in a hexagonal distribution pattern and have 6 direct neighbours (no. 1…6) on average. b In this case, Norway spruce is in the centre and the environment is by 100% (6 trees/6 trees × 100), the category is ss. c Here, Norway spruce is surrounded by three Norway spruce and three European beech, the admixture of another species apart from Norway spruce is > 0 and the group is sb. d Here, the grouping is bs, as beech is surrounded by Norway spruces (ss Norway spruces in neighbourhood of Norway spruces, sb Norway spruces in neighbourhood of European beeches, European beeches in neighbourhood of European beeches, bs European beeches in neighbourhood of Norway spruces)

In an advanced stand development phase like Kranzberg Forest, the trees in even-aged stands grow more or less in a hexagonal pattern (Prodan 1968a, b), i.e., each tree has on average six direct neighbours (n = 1…6) (Fig. 2a). Figure 2b shows an ss constellation where the central tree is a spruce tree surrounded by spruces. The proportion of other tree species in its environment is 0% as. The constellation in Fig. 2c results analogously in a group membership of sb and \(m_{{{\text{other}}}} = 3/6 \times 100 = 50\;\%\). Figure 2d shows a constellation of bs where beech is surrounded just by spruces, so that \(m_{{{\text{other}}}} = 6/6 \times 100 = 100\;\%\). We choose a rather strict separation between monospecific and mixed environments. As soon as the environment included another species other than the species of the central tree, it was re-characterized. In other words, only completely pure tree groups were characterized as ss or bb.

Calculation of resistance and resilience

The response of tree basal area increment, iba (cm2 yr−1), to the natural drought stress events in the years 2003 and 2015 was characterized by three different phases: (a) the growth PreDr in the 3-year-periods before the drought years 2003 and 2015, respectively, (b) the growth Dr during the drought years 2003 and 2015, respectively, and (c) the 3-year growth PostDr after the two drought years 2003 and 2015, respectively (Lloret et al. 2011). Indices for resistance, Rt = Dr/PreDr, recovery, Rc = PostDr/Dr, and resilience, Rs = postDr/PreDr, were used for the characterization of the stress response patterns. Resistance quantifies the growth decrease from the pre-drought period to the drought period. Rt = 1 indicates complete resistance; the further the value decreases below Rt = 1, the lower the resistance. Recovery describes the tree growth response after the drought period. Rc = 1 indicates a persistence at the low-growth level even after the drought period, values of Rc < 1 indicate a further decline, and values Rc > 1 represent a recovery from the drought period. Resilience is defined as tree growth after the drought period compared to the tree growth before the drought period. Rs values >  = 1 indicate high resilience with growth levels that are equal to or above the level before the drought event; Rs values < 1 indicate low resilience with growth levels below the one before the drought period. For a more detailed description of these indices, see Lloret et al. (2011).

Indexing, trend elimination, and smoothing

To evaluate the individual tree growth response to drought, we used the original annual iba data from the permanent girth tapes. We used the original data without trend elimination, smoothing, etc. due to the following reasons. (i) in contrast to the annual diameter or tree ring width growth, the trends of the annual iba growth rates tracked more or less parallel to the x-axis, except near drought years (2003, 2015) and the throughfall exclusion period (2014–2018). Therefore, no significant up- or down-age trends would bias the resistance or resilience analyses. (ii) The time span from 1998–2018 was too short to smooth or eliminate any trend, since in this time span, there were two natural drought events (2003, 2015) and one experimentally induced growth decline resulting from water limitation. Any attempt to fit a smooth curve through the 20-year-period would be questionable as the period was too short and more than a quarter of the period would have been overlayed by non-age-related disturbances; (iii) because the stands are even-aged and the trees all show more or less the same age trend. This applies especially for the trees of the precipitation exclusion experiment, as they were all dominant and even more homogeneous in the growing conditions and trends than the full data set. (iv) We compared the results only between groups with the same general age trend (Norway spruce vs. European beech, intra- inter-specific growth, TE vs. CO), so any influence of the age on the resistance or resilience indices should be eliminated as the trends in both groups were similar.

Estimation of mortality rates

Mortality rate calculations were based on the tree numbers, N, at the beginning Nb, and end Ne, of the observation periods. Using the compound interest formula, \(N_{{\text{e}}} = N_{{\text{b}}} \times 1.0{\text{m}}^{n}\), the mortality rates, m, and percent of mortality, \(m_{\% } = m \times 100\), were calculated for defined groups of trees (e.g., CO and TE). Hereby, n represents the length of the period in years. For our purpose, the basic equation \(N_{{\text{e}}} = N_{{\text{b}}} \times 1.0m^{n}\) was transformed to \(1.0m = \sqrt[n]{{N_{{\text{e}}} /N_{{\text{b}}} }}\) and \({\text{m}}\, = \,{1} - \sqrt[n]{{N_{e} /N_{b} }}\) to arrive at the mortality rate m. Note that the term 1.0 m is the convention of writing 1.0 + m in financial mathematics.

Mortality rates were calculated separately for the tree groups under episodic and experimentally extended drought and separately within these groups for Norway spruce and European beech.

The statistical software R 3.4.1 (R Core Team 2018) was used for all calculations, in particular the glht and t test functions for group comparison and lme function for regression analyses from the nlme package (Pinheiro et al. 2017).

Results

Growth response to natural episodic drought

Trees exposed to throughfall exclusion were excluded from the results presented in this section (= Tnt). For the drought events of 2003 and 2015, we show the periodic basal area increment in the 3 years before, during, and after the drought years (PreDr, Dr, PostDr) (Table 1). The long-term trend in annual basal area (± se) growth from 1998 to 2018 decreased slightly for Norway spruce and remained stable for European beech (Fig. 3). This long-term trend, however, was interrupted by dips in annual growth in 2003 and 2015, especially for Norway spruce. In general, European beech was much more resistant to the drought years.

Fig. 3
figure 3

Mean annual basal area increment (± SE) of (a) N. spruce (n = 268) and (b) E. beech (n = 141) from 1998 to 2018 excluding the trees under the throughfall exclusion roofs since 2014 (tree numbers refer to the year 2003). The broken vertical lines mark the drought event years 2003 and 2015

Norway spruce had reduced growth in both drought years (2003 and 2015), while European beech was only slightly reduced in 2003 and even increased in 2015 (see bold printed ratios of resistance Rt in Table 1). These data clearly show that the growth of Norway spruce is severely negatively impacted (50–60%) by drought, while European beech trees are much less effected.

The basal area increment of Norway spruce, in general, grew twice as much as European beech. In the drought year 2003, Norway spruce’s growth decreased to 41% and European beech to Rt = 76% compared to the 3-year-period before. After 2003, spruce was slower to recover than beech. In 2015, Norway spruce was also less resistant than European beech: spruce decreased to Rt = 51% of the initial increment level in response to the drought event, while beech increased incremental growth, surpassing the rate of growth in the 3-year-period before the drought (see bold numbers in Table 1).

To reveal any intra- and inter species-specific response pattern to drought, we analysed the growth response in the drought years 2003 and 2015 (Dr) compared with the 3-year-period before (PreDr) and after (PostDr) the events. Drought had a much stronger effect on Norway spruce growth compared to European beech despite their intra- and inter-specific environments (Fig. 4a). Since the relationships between the species, and between the intra- and inter-specific differences were similar, we show the results for 2003 only (Fig. 4).

Fig. 4
figure 4

Visualization of the growth response to the 2003 drought based on the annual basal area increment (± SE). The pre-drought growth in the period 2000–2002 is set to 1.0 (1.0-line). The growth in the drought year 2003 and in the post-drought period 2004–2006 was sketched in relation to this reference level. a On average, growth of Norway spruce (N. sp.) dropped steeply and recovered slowly, European beech (E. be.) was hardly affected by the 2003 drought. b When Norway spruce grew in inter-specific environment with beech (sb, broken lines), it was 20% less affected by drought compared with intra-specific constellations (ss). c When European beech grew in inter-specific environment with spruce (bs, broken lines), it was more affected by drought by app. 20% compared with intra-specific constellation (bb). However, the recovery and resilience was much faster in inter-specific compared with intra-specific constellation. (ss growth of Norway spruces in neighbourhood of Norway spruces, sb growth of Norway spruces in neighbourhood of European beeches, bb growth of European beeches in neighbourhood of European beeches, bs growth of European beeches in neighbourhood of Norway spruces)

Interestingly, Norway spruce was 10–20% less effected by drought when growing in the environment of beech trees (see sb in Fig. 4b), while reductions in spruce’s growth was greater in intra-specific spruce environments in 2003. Intra-specific competition (group ss, n = 192, mean 0.43 ± 0.02) had significantly reduced growth (p < 0.05) compared to inter-specific competition (group sb, n = 62, mean 0.56 ± 0.05).

European beech, on the other hand, were significantly more effected under drought in inter-specific environments but recovered quickly (see bs in Fig. 4c). Beech growing in the environment of other beech trees, in contrast, were much less affected by drought (see bb in Fig. 4c). In 2003, beech trees growing intra-specifically (group bb, n = 93, mean 0.87 ± 0.07) grew significantly more (p < 0.05) than trees in inter-specific environments (group bs, n = 23, mean 0.67 ± 0.05). This means that when beech grew in an inter-specific environment, it was by 23% more affected by drought compared to intra-specific conditions. This implies that in dry years, Norway spruce benefited from growing in mixed stands, obviously at the expense of European beech, as trees of the latter species significantly reduced their growth in inter-specific neighbourhood. Growth responses to experimentally extended drought by throughfall exclusion

Before the start of the throughfall exclusion experiment in 2014, we measured tree growth on the 6 CO and 6 TE plots for the years 2011–2013 to have an initial growth level reference. Compared to spruce, beech had less than half the mean basal area increment in the reference period 2011–2013 with some variation between the CO and TE plots of each species (Table 2). On the 6 CO and 6 TE plots, we recorded the course of growth of in total 102 dominant trees. The following analyses are based on 51 trees for each of the two species with 25 trees on the control plots and 26 on the treatment plots. Trees that suffered mortality were excluded from the analyses of growth reactions.

Norway spruce grew less in the period of 2014–2018 compared to the prior years, most likely due to the dry year in 2015. Trees in the CO plots exhibited a slight growth decrease from 2011 to 2013 compared to 2014–2018 from 17.4 to 13.4 cm2yr−1. However, on the TE plots, the decrease was much more severe from 19.9 to 4.6 cm2yr−1.

Using the relative growth (0.77) of the CO plot as a reference for the relative growth on the TE plots (0.23), the growth level was 0.30 (see Table 2, in bold and italic numbers), i.e., a loss of 70% in annual growth.

On the CO plots, European beech grew more in the period 2014–2018 than in the years 2011–2013, maybe because of late frost event in spring 2011 (Bayerische Forstverwaltung 2015). On the TE plots, we found a medium decrease from 7.1 to 4.9 cm2yr−1, i.e., a relative decrease by 31%. Using, analogously to Norway spruce, the relative growth (1.25) of the CO plot as a reference for the relative growth on the TE plots was 0.69; the growth level was 0.55 (see Table 2, in bold and italic numbers), i.e., a loss of 45% of basal area growth.

The basal area increment of Norway spruce, in both CO and TE plots decreased over time, mainly a result of the dry year in 2015 (Fig. 5, a, b). However, in the TE plots (Fig. 5b), the decrease was drastically more pronounced. In 2016 and 2017, a few of the trees had an upward trend, i.e., demonstrated recovery.

Fig. 5
figure 5

Annual basal area increment of individual trees from 2010 to 2018 shown for N. spruce and E. beech on the CO and TE plots. The vertical lines show the beginning of the throughfall exclusion (2014) on the TE plots. On the TE plots (2014–2018), Norway spruce reduced its growth (b), while European beech responded minimally to drought (d)

Most European beeches had a positive growth trend on the CO plots (Fig. 5c) and a negative trend on the TE plots (Fig. 5c) during the treatment period. However, some beech trees acclimated or even recovered during the throughfall exclusion period, i.e., in the years 2016–2018 (Fig. 5d).

In summary, we found clear negative responses to the experimentally induced drought in Norway spruce; the average loss in annual basal area growth amounted to 70% (Table 2). We found medium drought-induced negative responses in European beech; the average loss in annual basal area growth amounted to 45% (Table 2). Finally, we found some indications of acclimation and recovery for both tree species (see Fig. 6).

Fig. 6
figure 6

Visualization of the resistance to the 2014–2018 throughfall exclusion based on the annual basal area increment (± SE). Pre-drought growth in the period 2011–2013 is set to 1.0 (1.0-line, solid black), the growth in the years of the throughfall exclusion is shown in relation to the pre-drought level. a On average, growth of Norway spruce (N. sp.) dropped steeply and recovered slowly; growth of European beech (E. be.) dropped less strongly or even increased above the level of the pre-drought period after four years. b When Norway spruce grew in inter-specific environments with European beech (sb, broken lines), it was 20–30% less affected during the throughfall exclusion and recovered remarkably in subsequent years compared to spruce in an intra-specific constellation (ss). c When European beech grew in an inter-specific environment with spruce (bs, broken lines), the beech trees grew better in the first years of drought, but then fell behind the growth of beech in intra-specific environment (bb) (ss growth of Norway spruces in neighbourhood of Norway spruces, sb growth of Norway spruces in neighbourhood of European beeches, bb growth of European beeches in neighbourhood of European beeches, bs growth of European beeches in neighbourhood of Norway spruces)

To more closely examine the stress response to the throughfall exclusion, we analysed the annual basal area increments in the years 2011–2018 (Fig. 6 and Table 3). We first compared the performance of the species (Fig. 6a). Note that, in Fig. 6, the reference period 2011–2013 was marked by a bold horizontal line at level 1.0, reaching from 2011 to 2013. The mean year of this period is 2012. To visualize the growth after this reference period, we drew a connecting line from 2012, the mean of this reference period, to the relative growth in the years 2014; from there, we continued the line to the relative growth in 2015–2018.

Table 3 Comparison of the drought resistance on the TE plots in the experimentally induced drought period 2014–2018 on annual basis

Both species strongly reduced their growth in 2014, the first year after the throughfall exclusion experiment was initiated. Norway spruce continued to decrease over time but stabilized in 2016–2018, while European beech stabilized earlier and recovered to the initial level by 2018 (Table 3).

As the average growth of European beech in the reference period 2011–2013 was probably reduced by the late frost in spring 2011, we also calculated the growth response in the years afterwards after elimination of the year 2011 from the reference. However, this hardly changed the results as the beeches quickly recovered already in 2012 from the late frost. This is visible in Fig. 5c and d, where the course of beech growth shows a strong upward trend in 2012.

For Norway spruce in particular, environment effected growth after experimentally extended drought stress (Fig. 6b). Spruces growing in the environment of other spruces exhibited decreased growth much more than spruces close to beech. There was a significantly lower stress response and greater growth in spruce growing in inter-specific environments compared with those growing in an intra-specific constellation (Table 3).

Beech trees in both inter- and intra-specific conditions, however, responded to drought similarly at first (Fig. 6c). However, from 2016, we found significant differences between the two groups, i.e., beech in inter-specific environments outgrew those with spruce as neighbours (Table 3).

Mortality of Norway spruce and European beech caused by experimentally extended drought

Mortality of trees within the plots with experimentally extended drought would certainly impact the water supply and growth of the remaining trees, and, therefore, analysed the mortality of Norway spruce and European beech in the TE plots compared with untreated reference groups, i.e., the CO group and the group of all trees of the site Kranzberg Forest without the ones of the throughfall exclusion experiment (= Tnt) to improve the interpretation of our results.

The mortality rate under natural conditions between 1998 and 2018 was 1.24% for Norway spruce, 2.00% for European beech, and 1.50% for the total stand. In the period 1998–2019, no tree thinning occurred in the plots, so the given mortality rates represent the mean dropout under self-thinning conditions (Table 4).

Table 4 Tree numbers at the beginning, Nb, and at the end, Ne of the defined observation periods for the trees of the Tnt group and for the trees of the TE and the CO plots

The throughfall exclusion experiment CO plots were also not thinned. Mortality rates of 0.00% Norway spruce, 0.83% European beech, and 0.45% for the total stand occurred between 2011 and 2018. Tree mortality under throughfall exclusion was much higher than for the Tnt group. However, since 2011–2018 does not encompasses as many years as the period from 1998–2018 and climatic conditions differed between the two periods, this comparison should be viewed within the context of its limitations.

The comparison of the CO plot’s with the TE plot’s mortality rates is more interesting, as they refer to the same time period 2011–2018. Trees in plots with extended experimental drought had mortality rates of 7.45% for Norway spruce, 1.46% for European beech, and 4.07% for the entire stand (Table 4). Therefore, both species had higher mortality in TE plots compared with the CO plots, and in the case of Norway spruce, the mortality rate was five times as high as the long-term mortality of the Tnt group. The ranking of the mortality rate between the considered groups was ss > sb > bs > bb. Although it is difficult to assess the final causes of mortality (drought, bark beetle, and competition for light), we assume that in case of Norway spruce, most of the mortality (70%) was caused by bark beetle despite of the chemical control measures, some directly by drought (20%) and the rest by self-thinning due to competition (10%). The latter assumptions are based on the annual assessment of the vitality of all individual trees on the 12 plots. An indication for mortality caused by bark beetle were boring holes in the bark, boring dust on the ground and galleries under the bark. We assumed dropout by self-thinning in case of subdominant trees that became continuously more competed by their neighbours in the previous years. In case of those trees with transparent crowns that died although showing neither bark beetle infestation nor suppression by neighbours, we assumed a dropout by drought stress.

Discussion

Many studies have tackled species-specific drought resistance outcomes in monoculture tree plantings. However, species structural and functional trait differences can result in a particular species-specific stress responses when growing in monocultures (Bréda et al. 2006; Niinemets and Valladares 2006) and a potential reduction of stress response when growing in mixture (Ammer 2019; Grossiord 2018). Norway spruce is commonly assessed as a highly drought susceptible species (Lévesque et al. 2013; Zang 2012), while European beech, although under debate (Rennenberg et al. 2004), is less affected by drought (Ammer et al. 2005; Ewald et al. 2004). Whether tree species growing in mixtures can reduce drought susceptibility may depend among other things on the species combination (Metz et al. 2016; Pretzsch et al. 2013), the site conditions (Grossiord et al. 2014; Trouvé et al. 2017), and the stand density (Bottero et al. 2017; Sohn et al. 2016).

Norway spruce and European beech represent species with different hydraulic systems (xylem anatomy). Spruce exhibits a more isohydric strategy (Lyr et al. 1992), reducing stomatal conductance at early stages of soil drought. In contrast, beech displays a more anisohydric strategy, with less stomata sensitivity to soil drought, allowing for more carbon gain, and stem and root growth during prolonged time spans under mild to moderate drought (Leuschner 2009; Nikolova et al. 2009). These differences along with the high drought susceptibility of Norway spruce at the edge of its natural range and the maximum stand density within the experimental stands used here may have contributed to the substantial and lasting decrease in spruces’ growth compared to the minor growth reduction of beech under both episodic (Fig. 4a) and extended (Fig. 6a) drought.

Whether the potential of resource use, stress reduction, and even overyielding in mixed stands can be exploited by a given species assemblage depends on the respective site conditions (Forrester et al. 2014). Under ample water supply, e.g., a spatial or temporal complementarity of water uptake may be less useful than under drought. This explains why even rather complementary tree species may change the way which they grow in mixtures from beneficial to disadvantageous along ecological gradients (Pretzsch et al. 2015).

The temporal shift in the water uptake, i.e., that the transpiration of Norway spruce starts earlier than European beech (Rötzer et al.2017a), may explain the benefit of Norway spruce when growing in inter-specific neighborhood in the analysed stands (Figs. 4b and 6b). We hypothesize that spruce in proximity to beech benefits from a better water supply in the spring when beech is still leafless (see e.g. Rötzer et al. 2017a). This pre-emptive water uptake by spruce may reduce the water availability and growth of beech in the environment of spruce as observed on the TE plots of this study. This assumption is substantiated by measurements of soil moisture and water uptake by depth (Goisser et al. 2016) and micro-dendrometer trajectories (Rötzer et al. 2017b) which show seasonal shifts during spring drought and negative (during summer drought) soil moisture effects of beech neighboured by spruce. This underlines that the time of the year in which a drought occurs in mixed-species stands determines which species may benefit or lose in inter-specific neighbourhood.

The slight basal area growth recovery of both species after the initial downtrend of under experimentally extend drought (Fig. 6) is of special interest as it suggests an ability to adapt to drought stress.

Enhanced compensation growth of fine roots upon drought (e.g., in beech, Meier and Leuschner et al. 2008), adjustment of the mycorrhiza to an increased share of long-distance exploration types (Nickel et al. 2018), and acclimation of the branches and leaves to drought (Barbeta and Penuelas, 2016; Tomasella et al. 2018) may be effective measures of drought acclimation.

An increase of mortality of Norway spruce, e.g., caused by bark beetle attacks combined with extended drought, may reduce the stand density and, in this way, may improve the water supply of the remaining trees on the TE plots and contribute to their recovery.

The analysed stands are within the range of natural occurrence of European beech but at the limit of the distribution range of Norway spruce (Bayerisches Staatsministerium für Ernährung, Landwirtschaft und Forsten (2001). The site conditions at Kranzberg Forest allow both species nearly maximum productivity, indicated by the site indexes of O40 according to the yield table of Assmann and Franz (1963) for Norway spruce and site class I, according to Schober (1975) for European beech. But growing at the edge of its ecological niche, Norway spruce achieved its maximum productivity in years with ample water supply and when disturbances [e.g., bark beetle (Ips typographus L.) or gregarious spruce sawfly (Pristiphora abietina (Christ.) (Hym., Tenthredinidae)] were controlled by forest management (Skatulla et al. 1989; Wermelinger 2004). Living at the edge of an ecological niche can have amplified deleterious effects on species when small temporal environmental changes can have strong non-linear effects on growth and fitness. For Norway spruce, this means that the trees may be more susceptible to decline in drought years (Biermayer and Tretter 2016; Kölling et al. 2009) or that the trees have a requirement for facilitated positive inter-specific interactions (Brandl and Falk 2019; del Río et al. 2014; Pretzsch et al. 2012). In addition, Norway spruce is generally characterized by rather high mortality rates also in its natural range (Synek et al. 2020).

The restrictions and risks of cultivating Norway spruce beyond its natural occurrence are important to understand for forest practice. Because of its high productivity, excellent timber, and multi-purpose use, Norway spruce is highly valued and has a long history and tradition, especially in monocultures far off its natural range in mountainous regions of Central Europe and the Boreal region. The increasing tree damage in monocultures by both biotic and abiotic (snow breakage, wind) disturbances have resulted changes to forest practices including a move away from planting Norway spruce solely in monocultures. An alternative is mixed stand plantings that support more stable tree species, e.g., European beech, silver fir, Douglas-fir, or Scots pine, while maintaining a significant population of Norway spruce.

The silvicultural tools mitigating forest damage from drought are comprised of a selection of well acclimated species and provenances (Atzmon et al. 2004; Arend et al. 2011; Bolte et al. 2010; Zang et al. 2011), downregulation of stand density (D’Amato et al. 2013; Sohn et al. 2016), and thinning (Gebhardt et al. 2014; Pretzsch et al. 2018; Rodríguez-Calcerrada et al. 2011). An additional measure may be increased tree species mixtures although not yet rated effective for drought mitigation in general (Grossiord 2018). Our study provides an example of how tree mixtures can reduce stress and allow for continued growth of Norway spruce when growing closely mixed with European beech. This required single tree mixture, whereas most common in forest practice are mixtures in groups or clusters. Cultivation of European beech in two or three groups or clusters embedded in Norway spruce stands has the economic advantages of facilitated beech establishment (Wagner et al.2010), better timber quality when growing in intra-specific environment (Höwler et al. 2019; Pretzsch and Rais 2016), and an easier harvest (Hanewinkel 2001). In the common group or cluster mixtures, Norway spruce would most likely mainly benefit when growing at the edges of the beech groups, in close environment of beech, but not in the other zones. This means that Norway spruce, growing in the warm dry limit of its natural distribution, seems to be facilitated most effectively when growing directly next to groups of European beech which is natural in this area. This suggests that the choice of a climate smart species mixing pattern might be another tool in the silvicultural package of measures mitigating drought damages.

Conclusions

Experimentally extended drought established by a 5-year throughfall exclusion experiment enabled new insights into how Norway spruce and European beech may respond to future climate change scenarios that predict longer and more intense drought periods. The extended drought caused a drastically reduced growth in the first years, followed by a less severe decline in the subsequent period. To some extent, both species were able to acclimate to the drought and recover from the initial growth collapse, after exposure to episodic droughts. Norway spruce benefited significantly from growing in the environment of European beech, while beech overcame drought slightly better in intra- specific environments.

The considered site is representative for many areas in Southern Germany where Norway spruce is cultivated beyond its natural range, and while it can achieve optimal productivity under average climatic conditions, which becomes susceptible to drought and biotic disturbances during dry years.

Many recent studies show that tree species mixing can result in overyielding compared to monospecific stands and can increase the resistance of Norway spruce against biotic disturbances (e.g. bark beetle damages). The mixture of the highly productive and economically valuable Norway spruce with stabilizing trees species such as European beech may reconcile economy with ecology. The revealed drought stress relief of Norway spruce in inter-specific environments may be a strong argument in favour of a transition to mixed species forest stands and their superior ecosystem services.

Author contributions statement

HP, TG, and KP initiated the project. TR, K.-HH, TG, KP, TB, and HP developed and established the experimental design, HP and TR evaluated the data and wrote the manuscript. TR, TG, K.-HH, KP, and TB revised the manuscript.