Growth and drought reaction of European hornbeam, European white elm, field maple and wild service tree

Considering their drought tolerance and growth characteristics, rare native tree species are well-suited admixed species for the development of climate-stable forests in Central Europe. In our study, we assessed the growth and drought reaction of the four rare native tree species European hornbeam (Carpinus betulus L.), European white elm (Ulmus laevis Pall.), field maple (Acer campestre L.), and wild service tree (Sorbus torminalis (L.) Crantz). Based on tree-ring data, we (I) evaluated their species-specific growth characteristics and variability and examined the influencing site and tree characteristics on annual growth. (II) We quantified their reaction to single drought events, also depending on site and tree variables. (III) We compared our results to oak (Quercus robur L., Quercus petraea (Matt.) Liebl.) and European beech (Fagus sylvatica L.). As they are well-known Central European tree species, there is a broad knowledge about their growth and drought response across wide geographical ranges available. Bringing the results of European beech and oak in relation with the rare native species, it allows to categorise their growth and drought reaction and to contextualise their performance. Our results show, that besides European white elm, the rare species showed an overall lower annual growth with a higher variability than European beech and oak. However, especially field maple and wild service tree were better adapted to drought than European beech and partially even recovered better than oak. Combining the aspects of growth stability and drought tolerance, we conclude that rare native tree species are well suited as admixed species in future forest stands. European hornbeam is a suitable match for European beech on wetter sites, while field maple and wild service tree are a sensible complement for the climate stable oak on drier sites.


Introduction
In recent years, Central Europe has experienced several severe drought summers with high temperatures and low precipitation . Climate change projections for Central and Northern Europe predict even drier and hotter climate conditions in the future, with more frequent, severe and prolonged droughts (Cook et al. 2020;Eyring et al. 2016;Ionita and Nagavciuc 2021;Spinoni et al. 2019;Zscheischler and Seneviratne 2017). The associated consequences for forest ecosystems are substantial. Persistent climatic changes enhance the risk of growth decline, forest mortality (Allen et al. 2010(Allen et al. , 2015Senf et al. 2020), and diminish forest productivity (Rita et al. 2020), diversity, and carbon storage (Pilli et al. 2022). Recent droughts in Central European forests have caused canopy damages (Beloiu et al. 2022;Buras et al. 2020;Philipp et al. 2021;Sturm et al. 2022; Thonfeld et al. 2022) resulting in reduced crown size and height growth (Jacobs et al. 2021). The composition of Central European forests is expected to change in future due to different drought resistance of tree species. The introduction and promotion of alternative admixed species is a viable option to increase forest stability (Sturm et al. 2022). In this context, non-native tree species, as western hemlock (Tsuga heterophylla (Raf.) Sarg), northern red oak (Quercus rubra L.) or Communicated by van der Maaten. sweet chestnut (Castanea sativa Mill.), that are expected to be well adapted to hot summer temperatures and drought seem to be promising and recommended for an increased cultivation (Bolte et al. 2009;Hoffmann et al. 2018;Pötzelsberger et al. 2020). However, non-native species might be maladapted to colder winter temperatures and late frost events, which are also expected to increase in the future (Kreyling et al. 2015;Vitasse et al. 2019;Zohner et al. 2020). Therefore, native and drought resilient species may be a more reasonable alternative to complement present prevailing species. Other than for non-native species, they exhibit less uncertainties concerning the exposure to pests, adoption by local fauna, the potential use of their wood and their influence on other ecosystem services (Castro-Díez et al. 2019;Matevski and Schuldt 2021;Sapsford et al. 2020).
Among native species, especially European hornbeam (Carpinus betulus L.), European white elm (Ulmus laevis Pall.), field maple (Acer campestre L.), and wild service tree (Sorbus torminalis (L.) Crantz) seem to be promising options for admixture in Central European forests. The distribution range of European hornbeam is mainly determined by annual mean and winter temperature of the region (Varol et al. 2022). It tolerates annual mean temperatures up to 15 °C, which qualifies it for the use in Central European silviculture, where its future distribution range is considered to still be congruent to its present one (Koch et al. 2022). In the strictly preserved zone of Bialowieza national park, it even increased its share in the tree species composition in the last decades (Brzeziecki et al. 2018(Brzeziecki et al. , 2020. Severe and short droughts can be withstood by European hornbeam due to a relatively stable sap flow under stress conditions (Leuzinger et al. 2005). Seedlings, however, were found to show a high drought sensitivity, but also a high capacity to recover (Beloiu et al. 2020). In addition, Scharnweber et al. (2020) found a strong growth decline of mature European hornbeams after the drought of 2019. Overall, detailed studies about the drought tolerance of adult trees of European hornbeam are currently scarce.
European white elm is a tree species of riparian forests with a current distribution range that covers climates from dry and cool to warm and moist (Collin and Bozzano 2015). It is expected to keep (Koch et al. 2022) or even expand (Thurm et al. 2018) this distribution range in central Europe under predicted climate change scenarios, especially where a decline of current dominant species is observed. At the same time, European white elm is able to maintain a high growth performance under different climatic scenarios (Thurm et al. 2018). Mild drought stress does not have an influence on morphological traits of its saplings (Black-Samuelsson et al. 2003). However, drought reactions may highly vary between different provenances (Black-Samuelsson et al. 2003;Venturas et al. 2015). Leonova et al. (2022) found that European white elm, similar to oak, can maintain a constant amount of fine root biomass under drought conditions and can even enlarge it during drought.
Field maple favours a mild climate with low autumn precipitation and occurs within the entire temperate climate zone, except high mountain areas, regions with high precipitation (e.g. the British Isles) and the Mediterranean (Coudun et al. 2006). In a previous study, field maple showed a high drought tolerance, even exceeding that of oak (Kunz et al. 2018). In a greenhouse experiment, seedlings were found to have a high drought resistance (Kunz et al. 2016), which could, however, not be confirmed under field conditions where a high mortality and low recovery was noted (Beloiu et al. 2020). In general, field maple is expected to be well adapted to future climate change scenarios (Walentowski et al. 2014).
Wild service tree is currently growing in warm and dry parts of Europe (Rasmussen and Kollmann 2004) which are regions that correspond to future climate conditions in Central Europe. This adaptation makes the wild service tree a promising tree species under climate change conditions (Walentowski et al. 2014). A recent study noticed a higher drought tolerance than European beech, but not than oak (Kunz et al. 2018). Seedlings are sensitive to drought, but recover quickly afterwards (Kunz et al. 2016).
General growth patterns of the species are widely unknown. Only for European hornbeam, a yield table for Northern Germany exists (Lockow 2009). For wild service tree, Pyttel et al. (2013) used stem analysis to determine annual increments and deduced age-height and heightdiameter curves. In other studies, rare species are often aggregated in categories as "other hard wood species". For weighing out different species against each other and choose sensible species for cultivation, it is, however, important for forest managers to assess the potential growth of the species. It is also crucial to further evaluate the potential of rare domestic tree species to cope with climate change to develop adapted management strategies for Central European forests. Hence, more information on drought response patterns of these species is urgently needed. Most studies addressing the drought sensitivity of rare species are based on experiments with seedling or young plants or simulations (e.g. Kunz et al. 2016;Thurm et al. 2018;Varol et al. 2022). Studies using data from older trees are underrepresented. However, not only drought tolerance is an important characteristic of tree species for forest managers, but also a reliable and stable growth with economic potential.
To make our results for the rare species easier to access for forest managers and forest scientists, we compare them to European beech (Fagus sylvatica L.) and oak (Quercus 1 3 robur L., Quercus petraea (Matt.) Liebl.). As these are the dominant deciduous tree species in central Europe (Meyer et al. 2020), their growth and drought reactions are already well studied and widely known. European beech shows pronounced growth declines in large parts of is distribution range as a reaction to recent drought events (Leuschner 2020;Schuldt et al. 2020). Although it has been found to have the potential to genetically adapt to increasing drought stress (Petrik et al. 2022;Pfenninger et al. 2021), a declined growth and distribution range are anticipated, especially at its drier and warmer distribution limits (Del Martinez Castillo et al. 2022;Muffler et al. 2020). Pedunculate oak and sessile oak show a higher drought tolerance than European beech (Mette et al. 2013;Meyer et al. 2020) and are expected to show stable growth even under drier and hotter conditions. However, a decrease in sap flow rates and canopy conductance were observed in response to extreme drought events (Süßel and Brüggemann 2021). These negative effects of extreme drought on the ecophysiological processes of oaks can be mitigated by mixing with other species, e.g. European beech or Scot's pine (Steckel et al. 2020;Stimm et al. 2021). In general, oak and European beech show different drought reactions (Scharnweber et al. 2011;van der Werf et al. 2007), with higher stability of oak (Kasper et al. 2022) and also a higher resistance and resilience towards drought (Meyer et al. 2020).
In our study, we used tree-ring data to: (I) Assess the dendrometric growth characteristics of European hornbeam, European white elm, field maple, and wild service tree and the overall variability of growth, depending on site and tree variables. We hypothesise that there are species-specific growth characteristics and differences in the influence of site conditions. (II) Quantify and compare the species' reaction to single drought events. We hypothesise that the species show differences in their resistance and resilience towards drought and their recovery after drought events. We further hypothesise that the speciesspecific drought reactions are influenced by the individual tree size and site-specific climate conditions. (III) Compare the growth characteristics and drought reaction of the four rare native tree species with the well-studied European beech and oak which themselves differ in their growth and reaction to drought. We hypothesise that rare species show generally a smaller annual growth than European beech and oak and that they are at least as drought tolerant as European beech.

Sites and sample tree selection
For each species, we took tree cores from trees in two stands on two different sites in South-Eastern Germany in winter 2020/2021. Except for wild service tree for which on the second site, we only took cores in one stand. For European hornbeam, we additionally sampled trees of two stands located in eastern and central Poland, respectively. In total, we sampled the selected rare species in 17 stands on 9 sites (Fig. 1). We selected stands where the species grew in large proportions to reflect mono-specific or only slightly mixed conditions. The sites per species ideally covered different site conditions. However, due to the rarity of the species the desired variation in growing conditions could not be always realised. Therefore, the two sites of field maple and European white elm did not differ strongly in site conditions. All stands were undergoing regular silvicultural treatment in the past. As all stands were even-aged, the two stands per site covered a dbhrange per species (Table 1). The two stands per site were chosen as close to each other as possible to ensure similar site conditions. In each stand, we sampled 15-20 trees, covering the whole diameter spectrum of the stand. Trees were randomly selected, however, not standing directly next to each other. Depending on the occurrence, additional 10-15 European beech or oak trees or both were sampled (in total 105 European beeches and 107 oaks). They were chosen to grow in the direct vicinity of each stand with the sampled rare species and within a similar species composition, past treatment, age and site conditions (Table 2). When selecting oaks, we did not differentiate between sessile oak and pedunculate oak. We additionally included tree-ring data from two existing experimental plots of European beech that were in close vicinity to the plots with rare species (Schmied et al. 2023). The values of basal area included in Tables 1 and 2 were determined doing an angle count sample at each cored tree.
The selected sites for European beech, field maple, oak and wild service tree were located in their core distribution area in Europe. The two European hornbeam sites in Poland showed lower precipitation rates than the German sites. Both sites of European white elm were located on sites with a rather high precipitation compared to the rest of the European white elm's distribution range (Fig. 2).

Climate data
For the sites located in Germany, we used data of monthly mean, maximum and minimum temperature and precipitation derived from a 1 × 1 km grid of the German weather service (DWD) (DWD Climate Data Center 2022a; b). For the sites in Poland, we used climate data of the Climate Research Unit (CRU) (Harris et al. 2020).
Using the temperature and precipitation data, we calculated the standardised precipitation evaporation index (SPEI) for each stand (Vicente-Serrano et al. 2010). The SPEI is a multiscalar drought index which combines precipitation and temperature data and is a well performing and frequently used index in studies which evaluate the impact of drought on forest growth (Bhuyan et al. 2017;Ionita and Nagavciuc 2021;Skiadaresis et al. 2019;Vicente-Serrano et al. 2012). For the potential evapotranspiration, that is included in the index, we were using the Hargreaves equation (Droogers and Allen 2002;Hargreaves 1994). We calculated the SPEI for periods of 3, 6, 10 and 12 months to also cover lagged climate-induced growth reactions.
To characterise the long-term climatic conditions on a site, we used the De Martonne aridity index (DMI) (Martonne 1926), based on the reference period of 1991-2020. The index is calculated as DMI = P∕(T + 10) with P being the sum of annual precipitation and T the annual mean temperature. For a general characterisation of sites, the DMI is more suitable than the SPEI, as the SPEI is a standardised index that approaches zero for long-term means.
For the sites located in Germany, we also used the soil moisture index (SMI) (Helmholtz Centre for Environmental Research 2021; Samaniego et al. 2013;Zink et al. 2016) obtained from the Helmholtz Centre for Environmental Research (Helmholtz Centre for Environmental Research 2021; Zink et al. 2016). The SMI is an index scaled between 0 and 1 describing the soil moisture in comparison to a longterm expected value. It uses interpolated climate data and implements it into the hydrological model system mHM Samaniego et al. 2010) to simulate soil moisture at a resolution of 4 × 4 km. In this paper, we used  the data of the total soil column (1.8 m) from 1951 to 2020 in a monthly resolution.

Sample and data preparation
For the collection of tree cores, we used a HAGLÖF increment corer with a diameter of 5 mm. We extracted two cores per tree at breast height (dbh, 1.30 m), one from northern and one from eastern cardinal direction. This minimised the influence of reaction wood induced by the main wind direction from SW and allowed a better representation of overall growth Speer 2010). Since wild service tree and European white elm are very rare species with valuable timber, we were only able to extract cores at a height of 30 cm due to forest owner preferences. For an easy handling and measurement, we glued the increment cores on wooden boards and subsequently sanded with increasingly finer abrasive paper from 400 to 800 grit to enhance the visibility of tree-ring borders. Measurements were taken to the nearest 1/100 mm using a digital positioning table (Kutschenreiter and Johann; Digitalpositiometer, Biritz and Hatzl GmbH, Austria). Visual crossdating was performed based on common matching patterns of wide and narrow rings (Stokes and Simley 1996;Schweingruber et al. 1990;Speer 2010). Afterwards, we statistically verified crossdating using the dplR package in R (Bunn 2010).
Due to the low visibility of the tree rings of European hornbeam, we used an alternative method for measurements of this species. After sanding, the cores were photographed using the microscopic camera of a Lignostation (Rinntech) with a resolution of 1/10 mm. The resulting photographs were imported into Adobe Photoshop (version 22). Using a high pass filter with a radius of 40 pixels, a linear light blending mode and changing the displayed colour spectrum, tree-ring borders were easier to detect. More compact wood sections appeared more yellow, less compact section more blueish. The resulting photos were then imported into the Lignovision software (Rinntech, version 1.37). Here, tree rings were marked, measured, and visually crossdated.
During the measurement and crossdating process, we had to reject several cores due to an extremely weak visibility of tree rings that made it impossible to measure and crossdate the cores. The final number of cores per plot used in this study can be found in Tables 1 and 2.  Table 1 can be found in the supplementary material

Selection of drought years
For the selection of single drought years, we identified one single drought indicator. Requirements of this index were a reliable identification of drought, a high spatial resolution, and a high correlation with tree growth over all species and stands. Therefore, we calculated bootstrapped Pearson's correlation between all climate variables and drought indices and our species-specific site chronologies. For all variables and indices, we used monthly values as well as seasonal means for spring (March-May), summer (June-August), autumn (September-November), winter (December-February of the following calendar year) and the vegetation period (April-September). For precipitation, we summed up the values for the corresponding season.
To calculate average site-overarching coefficients for each tree species, we transformed the correlation coefficients using Fisher's z scores to avoid underestimations due to skewed distributions (Silver and Dunlap 1987) and averaged them. After the calculation, they were back-transformed. Chronologies of young European white elm and oak of stand OFF1 were excluded from the analysis of climate growth relationships, as well of the analysis of drought reaction due to the short time-period covered. Coefficients of correlations for all sites and species can be obtained from the supplements (Figs. S1 and S2).
Of the climate data available for this study, SMI and SPEI were the most reliable drought indices (Schwarz et al. 2020). As they combine both temperature and precipitation, they are more suitable to determine drought years than indices only considering one of the variables (Zscheischler and Seneviratne 2017). While SMI also includes soil data, SPEI was available in a higher spatial resolution. The seasonal SPEI3 for summer (SPEI3 summer ) showed the highest correlation and lowest inter-stand variation for all species and stands and was used in further analysis. Drought years selected with this index are also resembling future climatic conditions, where drought events will mainly happen in summer, but also in spring months (Cook et al. 2020;Ionita et al. 2020;Ionita and Nagavciuc 2021), as it includes mainly summer data but due to its 3-month frame also spring months. Furthermore, SPEI3 is a suitable index as it is targeting shortterm extreme events and seasonal precipitation changes and is widely used to assess drought in forest ecosystems (Bachmair et al. 2018;Dell'Oro et al. 2020;Obladen et al. 2021;Spinoni et al. 2017).
Drought year identification was based on classifications by Slette et al. (2019). We considered all years with SPEI3 summer ≤ -1 as potential drought years. From these years, we chose the three driest years of the time covered by the chronology of each stand separately (Supplements Table S2).

Analysis of tree-ring data
Tree-ring data were detrended using a 30-year spline with a 50% frequency cutoff (Cook 1992;Cook and Peters 1997;Klesse 2021). This spline smoothes low-frequency variation in the tree-ring series which are associated with management (e.g. due to thinnings) and age trends. Inter-annual variability and high-frequency variation, however, are still preserved. The detrending procedure resulted in dimensionless ring-width indices (RWI) (see Fig. 3).
For the two cores of each tree, the mean value per year was calculated using a Tukey's biweight robust mean. Finally, we transformed the RWI of each site into respective average chronologies. Years with measurements from less than 5 trees were truncated. In addition, we calculated the subsample signal strength (SSS) (Wigley et al. 1984) for each plot and year covered by the data. The SSS (Wigley et al. 1984) is a value quantifying the correlation between a subsample of time series and a larger sample and is an indicator for the strength of representation of the larger sample by the smaller subsample (Buras 2017). Years in a chronology with a SSS smaller than 0.85 were truncated, following the recommendations of Wigley et al. (1984). The quality of the chronologies was furthermore assessed by using common dendroecological key figures. The values can be obtained from Supplementary Table S1.

Variability of growth and influence of site and tree variables
To assess the dispersion of tree-ring series, we used the Gini coefficient, calculated from the detrended year ring widths, as proposed by Biondi and Qeadan (2008). The Gini coefficient is a frequently used index of data variability, covering all lags in a tree-ring series. To test differences between the species-specific Gini coefficients, we pooled the detrended tree-ring data for each species and calculated bootstrapped 95% confidence intervals (Dixon et al. 1987).
To assess the sensitivity of annual growth on tree and site variables, we fitted a mixed effects model for each tree species. For the evaluation of annual drought conditions, we used the SPEI3 summer . The global model was fitted using the following formula: where t referred to the year and i to the random intercept on tree level and a 0 , a 1 , a 2 , a 3 , a 4 , a 5 and a 6 were regression coefficients and it the normally distributed error term. The b i variable respected the nested data structure on tree level. To consider the presence of an autocorrelation between consecutive tree-ring measurements, we included an autoregressive correlation structure (Pinheiro and Bates 2000;Venables and Ripley 2002).
The model assumptions of homoscedasticity and normality of residuals were checked visually using qq-plots and by plotting residuals vs. fitted values. The plots showed no violation of model assumptions.
The global model was then stepwise reduced using the Akaike Information Criterion (AIC). For European white (1) log(annual growth) =a 0 + a 1 * log(dbh t ) + a 2 * SPEI3 t + a 3 * DMI + a 4 * log(dbh t ) * SPEI3 t + a 5 * log(dbh t ) * DMI + a 6 * SPEI3 t * DMI + b i + it , elm, we additionally removed the DMI from the model, as both sites showed very similar values (36.11 and 38.68). We scaled and centred the predictive variables to enhance the comparability and comprehensibility of variables and their influence on the response variable (Schielzeth 2010). A summary of the unstandardised input variables can be found in supplementary table S3.

Reaction to single drought events
For the selected drought years, we subsequently calculated indices of resistance (R t ), recovery (R c ) and resilience (R s ) according to Lloret et al. (2011) to quantify growth responses towards droughts. These indices are commonly  (Schwarz et al 2020). They are calculated as ratios between the growth in a period before (PreDr), during (Dr), and after drought (PostDr): We chose variable pre-and post-drought periods of 2, 3 and 5 years as recommended by Schwarz et al. (2020). However, we could not find an evident difference between different periods. We decided to use a pre-and post-drought period of 2 years, as this period restricts the influence of other effects like mast years or defoliation by insects (Bottero et al. 2021;Schwarz et al. 2020) and also allowed us to include the recovery values of 2018 in our study.
To detect significant differences between the median values of the species for all Lloret indices, we used a Wilcoxon rank-sum test with Bonferroni correction and a p value of 0.05, as the assumption of normal distribution of the samples could not be accepted for all species.
As proposed by Schwarz et al. (2020), we also compared our species-specific relationship between resistance and recovery to the "line of full resilience". This line is derived by the relation between recovery, resilience and resistance (5), when the resilience values is set to 1 (6): The species-specific relationship follows a power-function with the following formula: The line of full resilience shows the hypothetical recovery values a tree needs to reach to fully obtain its pre-drought growth level (full resilience) for each resistance value. By observing the deviation and progression of the species-specific curve in comparison to the line of full resilience, we can rank and summarise the growth responses to drought and assess the recovery potential of the species.
To detect differences in the relationships of recovery and resistance between the species, we linearized formula (7) by using a logarithmic transformation and included the species as a linear term: We performed post hoc comparisons of coefficient combinations using the glht function from the multcomp package (Hothorn et al. 2015) to evaluate differences in the progression of the curve of full resilience between species.

Effect of site and tree variables on drought response
To assess the dependency of resilience, resistance and recovery on site variables and diameter, we fitted mixed effect models for every species, following the same approach as for the model of annual growth sensitivity.
This resulted in the following global models: where t referred to the year and i to the random intercept on tree level and a 0 , a 1 , a 2 , a 3 , a 4 , a 5 and a 6 were regression coefficients and it the normally distributed error term. SPEI3 t stands for the SPEI3 of summer months of the year t. The b i variable represented the random intercept on tree level. The model's explanatory variables were stepwise reduced, resulting in the model with the lowest AIC. The model assumptions of homoscedasticity and normality of residuals were again checked visually using qq-plots and by plotting residuals vs. fitted values. The plots showed no violation of model assumptions.

Statistical software
All analyses in our study were carried out using R, version 4.2.1 (R Core Team 2022). For calculating SPEI values, we used the SPEI package (Beguería and Vicente-Serrano 2017), described in Beguería et al. (2014). For detrending, the calculation of descriptive statistics and chronology building of tree-ring data, we used the R package dplR (Bunn 2010;Bunn et al. 2021) and pointRes (van der Maaten-Theunissen et al. 2015, 2021. Climate growth relationships were calculated using the treeclim package (Zang and Biondi 2015). Linear mixed effect model were fitted using lme4 (Bates et al. 2015) and lmerTest (Kuznetsova et al. 2017). For general linear hypothesis testing, we used the package multcomp (Hothorn et al. 2015), and for the calculation of the Gini coefficient, the package DescTools (Signorell 2022).

Growth rates and growth variability
The basic tree-ring values for all species can be obtained from Table 3. European white elm showed the highest maximal year ring widths of all species, with also the highest standard deviation. This resulted in a strong variability of diameters at a given age (Fig. 4). Compared to the other tree species, the tree-ring widths of European hornbeam did not vary much between different sites and stands. The dbh-age curves of the trees were very similar. This applied to both the plots located in Germany and the plots located in Poland. While mean tree-ring widths of wild service tree remained under the level of European beech and oak, its maximum tree-ring width showed a higher value than oak and European beech.
The growth variability of detrended tree rings was higher for the rare species than for European beech and oak (Fig. 5).
Here, European white elm showed the highest Gini coefficient, European beech the lowest.

Effect of site and tree variables on annual growth
The fit and reduction of the global model (Eq. 1) resulted in the species-specific models displayed in Supplementary  Table S4.
The annual growth of all species was significantly affected by dbh, annual SPEI and, was applicable, by the DMI. For European beech, European hornbeam, field maple, wild service tree and oak we could observe significant differences in annual growth between sites with poor and better water supply. On very dry sites, the annual growth decreased for bigger trees. On better sites, however, this decrease was weaker for field maple and wild service tree or did not happen at all for European beech, European hornbeam, and oak.
Here, we could even observe an increase in annual growth for bigger dbhs. For European beech and field maple also, the reaction on annual weather conditions, displayed by the SPEI, differed between sites. For these species, the increase of annual growth as a reaction to high SPEI values was higher on good sites than on very dry sites. All other species did not show a difference in their reaction between sites. However, higher SPEI values generally led to a higher annual growth. European white elm and European hornbeam showed a higher sensitivity and increase in annual growth towards wetter years for thinner trees.

Resistance, resilience, and recovery towards drought
We could sort the species by their median values for the index of resistance in the following order: Field maple (0.61) < European hornbeam (0.73) < wild service tree (0.74) < oak (0.76) < European beech (0.79) < European white elm (0.92) (Fig. 6). After testing the hypothesis of equal medians, we could differ between three groups: field maple with low resistance values, oak and European hornbeam with a middle position and European white elm with the highest resistance values. European beech and wild service tree took an intermediate position between European hornbeam and European white elm.
The median recovery values were ranked as follows: European white elm (1.05) < oak (1.13) < European hornbeam and wild service tree (1.18) < European beech (1.20) < field maple (1.63). Overall, after testing for equal means, we could not find significant differences between medians. Only field maple showed a significant higher recovery value than the other species.
Median resilience values followed the subsequent order: oak (0.86) < wild service tree, European whit elm and European hornbeam (0.90) < European beech (0.94) < field maple (1.02). Between the median of oak and field maple, we could observe a significant difference. For all other species, the hypothesis of equal medians could not be rejected. Figure 7 shows the relationship between recovery and resistance of the species and the theoretical line of full resilience. European hornbeam, field maple and European beech already showed a very close fit to the line of full resilience for very low values of resistance. Field maple and wild service tree intersected the full-resilience line even for lower resistance values than the other species. The model intercept of European beech was significantly higher (p < 0.05) than the one of wild service tree, European hornbeam, and field maple. In addition, field maple showed a significant higher intercept than oak. For resistance values above 1, field maple, oak and wild service tree showed the highest positive deviation from the line of all species.
In general, we could observe three major curve shapes. European hornbeam and European beech on the one hand showed a good fit to the line of full resilience for all resistance values. On the other hand, European white elm, wild service tree and oak showed a negative deviation from the line for low resistance values but, especially for wild service tree, with a an early intersection and a subsequent clear positive deviation. Field maple, as the only species, shows both a good fit on the line of full resilience for low resistance values and its exceeding for high resistance values.

Influence of site and tree size on drought reaction
Concerning the drought reaction, we could observe significant influences of tree size, drought severity and overall site conditions. However, the expression of the effects strongly differed between species.
The resilience of European hornbeam, European white elm and field maple was significantly affected by the severity of the drought year, displayed by the SPEI values. Field maple and European white elm were more resilient after milder drought years, while European hornbeam showed an opposite reaction. For field maple, this effect was especially visible on sites with a lower DMI, while on sites with more favourable climatic conditions, the resilience was similar over different drought intensities. In addition, for European beech and oak, we could observe an influence of the long-term  6 Index values of the Lloret indices of recovery, resilience and resistance for European beech, European hornbeam, European white elm, field maple, oak, and wild service tree. Significant differences (p < 0.05) between the medians of the index value per species, determined by a Wilcoxon rank-sum test with Bonferroni correction, are displayed with letters climatic conditions of the site, displayed by the DMI, and a significant influence of the dbh. For European beech, the resilience values were higher on sites with higher water supply and for thinner trees. For oak, we could observe the higher resilience of thin trees only on sites with a high DMI. On sites with a lower DMI trees of all diameters showed similar resilience values. Notable is, that for wild service tree, we could not find a significant effect of any variable.
Concerning the resistance, we could observe a significant influence of the SPEI values of the drought year for all species beside wild service tree. For oak, European white elm and field maple, a less severe drought led to higher resistance values. European beech showed low resistance values in very dry years on sites with a low DMI, while the resistance of European hornbeams increased with more severe droughts on those sites. Concerning the effect of tree size on the resistance European beech showed a higher resistance in bigger trees on sites with a high DMI and thinner trees on sites with a low DMI. Again, it is notable that for wild service tree, we could not find a significant effect of any variable.
The tree size did not have a significant influence on the recovery of European white elm, field maple, and wild service tree. For oak, we could observe a better recovery for thinner trees. Thinner European hornbeams showed a better recovery on favourable sites than bigger trees, while on sites with a low DMI no such size effect could be observed. For European beech, we could observe the opposite, with a better recovery of thin trees on sites with a low DMI and of bigger trees on sites with a high DMI. This significant influence of the climate conditions on the sites could also be observed for oak and field maple. Field maple and European beech showed generally higher recovery values on sites with low water supply and less severe drought years. On more favourable site, this effect was reversed. The recovery values of oak generally increased with a higher DMI. Wild service tree and European white did not show a significant effect of diameter, site and SPEI on their recovery.
In summary, we could observe a high susceptibility to severe drought for beech on sites with a low DMI, while oak was more stable, especially on good sites and over different tree sizes. Among the rare native species, field maple again stuck out with a high robustness towards drought, even on very dry sites and towards severe drought. European hornbeam reacted similarly as European beech, however, with Fig. 7 Species-specific relationship between resistance and recovery (black line (see Eq. 7) with confidence bands calculated with a bootstrap with 1000 iterations), the red line (see Eq. 6) represents the theoretical line of full resilience a higher tolerance towards severe droughts and more constant over different tree sizes. For wild service tree, our data showed a stable drought reaction without any influence of site and tree variables. European white elm showed a particularly high tolerance towards mild droughts.

Growth rates and growth variability
Concerning the age-growth relationships, we could observe that the growth of the rare species is within the growth range of European beech and oak on the same sites. When comparing the growth of European hornbeam with existing yield tables for northern Germany (Lockow 2009), sites in Germany reached a site index of II.5 (site WAL) and II (site EBR). Both sites in Poland reached a site index of II.25, despite their lower DMI. For all stands, the mean diameters were higher than the diameters specified by the yield table. This indicates that the validity of the yield table of Lockow (2009) might be regionally limited to stands in northern Germany with their specific climatic and soil conditions. On other sites, the yield table may underestimate the actual growth of European hornbeam. The similar growth patterns of European hornbeam in Poland and Germany underline its high ecological amplitude. Even in eastern Poland, under climate conditions that are limiting the distribution range of European beech, the European hornbeam shows high growth rates.
For wild service tree, the annual increments were similar to the ones observed by Pyttel et al. (2013). They also observed very low annual increments which they explained by the high shade tolerance of wild service tree and its ability to even survive under strong suppression. This indicates that some of the wild service trees used in our study might have been suppressed in the past or growing under a lot of competition. This also explains the high variation in the dbh-age curves showed in Fig. 4, where some trees had much higher diameters for a given age than others.
In general, European beech and oak grew more stable over all sites. The missing significant difference in annual growth between sites with different DMI of oak in the linear mixed model is another indication of the lower climatic sensitivity of oak in comparison to European beech (Scharnweber et al. 2011). The maximum tree-ring widths of the rarer species show, that under good conditions and with proper silvicultural management, they can reach high year ring widths. Combined with the results of the model assessing the sensitivity on site variables, we can see that the high Gini coefficients of rare species are also connected to a high susceptibility of the species to favourable years. Apparently, the species can easily implement good annual growing conditions in annual growth. The high Gini coefficients may also be related to a different silvicultural treatment of the species in the past. European hornbeams are often used as a so called "serving" tree species in oak stands for the shading of stems and the prevention of the development of secondary branches. After a release of suppression by the extraction of the canopy, a higher growth can be provoked. However, this again underlines that with enough growing space high year ring widths of the rare species are possible.
The lower sensitivity to drier years of trees with a smaller dbh that we observed could be related to their smaller height and thus the lower gravimetric potential they have to overcome to transport water in the tree trunk and the higher stomatal conductance and sensitivity that make smaller trees less prone to transpirational stress (Grote et al. 2016;Ryan et al. 2006;Ryan and Yoder 1997). The decline in growth with an increasing tree size on dry sites can be traced back to a typical age effect and the culmination of growth already at lower diameters on poorer sites. The decrease of annual growth in field maple and wild service tree even on better sites indicates that both species culminate earlier in their growth than European beech and oak. Again, it also indicates a different silvicultural treatment of the species in the past. The increase of growth with higher diameters on better sites for European beech, oak and European hornbeam points to a release from suppression or competition in higher ages. Outstanding is the high potential of annual growth of European white elm. Especially on good sites and in young stands, this species can easily obtain year ring widths of more than 1 cm.

Drought reaction
When comparing the Lloret indices of all species, not considering site and tree variables, we could not find many significant differences between species. Nevertheless, it was noticeable that field maple showed especially high values of recovery and resilience, however, with low resistance values. The resistance values of European beech being higher than the ones of field maple is contradicting the results of Kunz et al. (2016) who had reverse findings. This could be due to a different selection of drought years and stands in both studies and also the missing consideration between resistance and recovery values (Schwarz et al. 2020). However, when combining resistance and recovery values and comparing them with the line of full resilience, field maple and wild service tree stick out as drought tolerant species. Wild service tree due to its capability to surpass the line of full resilience for milder droughts and field maple for both its high recovery values in years with low resistance and its good performance after milder drought. This again matches the findings of Kunz et al. (2016). It once more shows the importance of rather combining both resistance and recovery than looking at the single indices when evaluating the drought tolerance of different species. Thereby, it can be considered that species with high resistance values were often showing low recovery values and vice versa, as after a strong decline, there is naturally a stronger potential for recovery. Species with higher values of recovery for lower resistance values, as wild service tree and field maple, can recover more easily even in years with a strong growth decline. Oak showed better overall resistance to drought than the rare species, however, in severe drought years field maple and wild service tree seem to be better adapted.
The good recovery of oak and European white elm in years with milder drought can be attributed to their wood anatomy. In contrast to the other species considered in this study, oak and elm are ring porous species. They are less sensitive to climate and drought than diffuse-porous species (Elliott et al. 2015) as they produce larger early wood vessels in which then the majority of the hydraulic transport takes place (Hacke and Sauter 1996;Zimmermann 1983). After a drought event, they can, therefore, recover more quickly as larger vessels are already produced at the start of the next growing season. Diffuse porous species must use their smaller vessels, created throughout the previous growing period. Embolised conduits can be repaired by refilling them (Hacke and Sauter 1996). Furthermore, our results correspond to the findings of Leonova et al. (2022), who found an increased fine root growth of European white elm during drought. This also leads to an enhanced growth during the next growing season.

Influence of size and site variables
Regarding the drought reaction, we could observe a higher drought tolerance of European hornbeam towards severe drought and on dry sites, as already observed by Leuzinger et al. (2005). This could explained by differences in their water consumption and stomatal regulation strategies, where European hornbeam follows a more isohydric behaviour (Köcher et al. 2009;Leuschner et al. 2019). In contrast to our results, Scharnweber et al. (2020) observed a strong growth decline for European hornbeam after two consecutive drought years. The drought tolerance of European beech and European hornbeam was similar in their study, with even a higher growth decline for European hornbeam. However, they suppose that a part of this growth decline can be attributed to an increased growth of roots in drought years as an adaptation strategy. This again matches our results where the higher resistance of European hornbeam towards severe drought could be observed for dry sites, where the European hornbeam might already be adapted to drought events. For field maple on very dry sites, we could observe a high resilience and very high recovery towards milder droughts. On more favourable sites field maple also recovered very well after severe drought events. This one more time points out the high drought tolerance of the species and its high potential on dry sites. As already studied by Kunz et al. (2018), field maple is a highly drought tolerant species with a high potential under climate change. Beloiu et al. (2022) found a low recovery and high mortality after drought for field maple saplings, while Kunz et al. (2016) found much younger seedlings to be highly drought resistant. This indicates that the initial drought resistance of the species may decrease with age. However, the high mortality after drought events may lead to a natural selection of trees with a higher drought resistance and, therefore, an adaptation to drought. At the same time, it may result in a higher drought tolerance of surviving trees due to an acclimation effect after the drought event (Kozlowski and Pallardy 2002). This again results in a higher drought tolerance of mature trees. For European white elm, we could not differentiate between different climatic site conditions, as the two sites had very similar DMI values. However, we could observe a high resistance against milder drought. Although more detailed research is needed, our study confirms the assumptions of Walentowski et al. (2014), that European white elm could be a suitable species under drier climate conditions. Site variables did not have any influence on drought resilience, recovery, and resistance of wild service tree. This is also a sign of the high stability of wild service tree under different site conditions. On both sites, we could find a good recovery after mild droughts, which was already observed by Kunz et al. (2016) for seedlings.
Oak could cope better with drought on wetter sites, however, showing higher resistance values for large trees and a faster recovery of small trees on dry sites. This can be attributed to the weaker stomatal control on photosynthesis of small oaks trees (Zang et al. 2012). In our study, we could observe the dependence of drought tolerance of European beech on climatic site conditions. On wetter sights, the tolerance of severe droughts was higher than on dry sites, while the tolerance towards mild droughts was lower. This higher susceptibility to mild drought on sites with a good water supply was already observed in other studies dealing with the growth of European beech towards its distribution margins (Cavin and Jump 2017;Muffler et al. 2020;Weber et al. 2013). It might be the result of a phenotypical reaction (Cavin and Jump 2017;Leuschner 2020) to site conditions and an acclimation after multiple drought years (Petrik et al. 2022). This could point to a more isohydric behaviour of trees on drier sites (Nguyen et al. 2017). The good recovery of beech in milder droughts on dry sites can be connected to the rapid recovery of leaf metabolism (Leuschner 2020).
High recovery values on wet sites after droughts might be related to the stimulation of fine root growth, as observed for saplings of European beech ).

Study limitations
Our study results show some limitation as tree cores were partially taken of different stem heights. The expression of climate sensitivity decreases on lower parts of the stem (Hoffmann et al. 2018) which may lead to some inaccuracy when comparing values of stands with different retrieval height. Furthermore, we did not sample dead trees as proposed by Schwarz et al. (2020). Therefore, we were not able to assess climatic factors that might have a very negative, even lethal, effect on tree growth. The low mean-Gleichläufigkeit values of European white elm indicate that the growth of trees on the plots was not primarily limited by precipitation and evapotranspiration, but rather by other site factors. As both sampled sites of European white elm were in riparian forests, the flooding regime and groundwater access could be factors determining tree growth. This follows the critique of Schwarz et al. (2020) and Zang et al. (2020), who suggest including more soil parameters in studies. For our study this was not possible, as, besides the SMI, no data on past soil moisture were available for those stands. The SMI, however, only has a resolution of 4 × 4 km, which again leads to imprecise results for our forest stands.
By sampling two sites per species, we tried to cover a variety of different site conditions. However, this was not always possible, as the occurrence of our species was often limited to sites where European beech or coniferous species were not dominating. These were, in case of field maple and wild service tree, very dry sites, in case of European white elm riparian forests. However, as seen in Fig. 2, our selected stands were located mostly in the core distribution range of the species and thereby represent typical climatic conditions. Nevertheless, an expansion of the study area with more plots on the edges of the distribution range could be interesting for future studies.
In our study, we did not account for the fact that repeated drought events may influence the drought reaction of the species. Especially for younger stands, the selected drought years might be the first droughts that the trees were experiencing. However, as all species covered younger and older trees with a different number of past drought events, we decided to neglect that effect in our analysis.

Conclusion: implications for silviculture
Even if European beech shows a certain drought resistance on sites with a good water supply and during mild droughts, it is generally displaying a high sensitivity towards drought. This points out a possible limitation on its occurrence and cultivation. Even on more favourable sites, the risk of drought events will increase in the future. As other rare species were more resistant towards drought on these sites, they could increase their shares in the species composition while the share of European beech might decrease.
Oak turned out to be a stable species, even under dry conditions. Comparably high increments combined with a high drought resistance make it a solid basis for future forest stands. For sites were in the future many very dry years are expected, wild service tree and field maple are good options to enrich these stands. These two species are very well adapted to drought and can keep up with oak with their ability in transforming good weather conditions into annual growth. The proportions of mixture should increase with increasing drought potential.
Concerning the drought sensitivity and growth, European hornbeam showed overall a similar behaviour as European beech, however, with a slightly smaller climate sensitivity and a stable growth even over different climate zones in Poland and Germany. Furthermore, our results show a good resistance against severe drought years. We, therefore, recommend European hornbeam as a complementary species for the enrichment of European beech stands. As for both species drought sensitivity increased with diameter, an adjustment of target diameters might be considered. In addition, oak stands with slightly better water supply could be suitable for an admixture of European hornbeam.
The higher drought tolerance of younger trees of wild service tree, field maple and European white elm is especially important concerning the relative uncertainty of weather development in the next years. Our results suggest that young stands with those species, that were recently established or will be established in the next years, can cope with a variety of different weather conditions and are, therefore, a low-risk option to sustain unstable stands on dry sites.
European white elm can be highly recommended as a species for riparian forest systems. Even if we could only get a limited assertion of its behaviour under drought conditions, our results generally indicate a high drought tolerance which implies a suitability also on dry sites. European white elm is prone to the invasive Dutch elm disease (DED), but much less than other native elm species like wych elm (Ulmus glabra Huds.) (Jürisoo et al. 2019(Jürisoo et al. , 2021. As common ash (Fraxinus excelsior L.) is also suffering from ash dieback, European white elm remains as a relatively stable species on floodplain sites, fulfilling both ecological and economical functions. The high growth rates also make it an interesting species concerning the carbon sequestration in forests.
In general, rare native species are a promising option to enrich the species composition on dry sites and to contribute to the establishment of climate tolerant forest stands. However, with ongoing climate change, even those species will eventually reach the limits of their drought tolerance. Therefore, the fight against climate change stays the most important factor to guarantee the stability of future forests.
Author contribution statement JS collected and analysed data and wrote the manuscript; GS supervised the tree-ring analysis and manuscript writing process; EU supervised data collection, analysis and writing process; HP initiated the study and supervised the whole data collection, analysis and writing process.
Funding Open Access funding enabled and organized by Projekt DEAL. This study was financed by the Bavarian State Ministry for Food, Agriculture and Forestry within the framework of the project W048 'Growth dynamics and climate sensitivity of rare domestic tree species'. In addition, this research has received funding from the European Union's Horizon 2020 research and innovation programme under the Marie Skłodowska-Curie grant agreement No. 778322 (CARE4C).

Data availability
The materials described in the manuscript including all relevant raw data will be freely available upon request from the corresponding author.

Conflict of interest
The authors have no relevant financial or non-financial interests to disclose.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/.