Impact of Climate Trends and Drought Events on the Growth of Oaks ( Quercus robur L. and Quercus petraea (Matt.) Liebl.) within and beyond Their Natural Range

: Due to predicted climate change, it is important to know to what extent trees and forests will be impacted by chronic and episodic drought stress. As oaks play an important role in European forestry, this study focuses on the growth response of sessile oak ( Quercus petraea (Matt.) Liebl.) and pedunculate oak ( Quercus robur (L.)) under contrasting climatic conditions. Analyses cover both site conditions of their natural occurrence (Southern Germany and Northeast Italy) and site conditions beyond their natural range (South Africa). The sites beyond their natural range represent possible future climate conditions. Tree-ring series from three different sites were compared and analysed using dendrochronological methods. The long-term growth development of oak trees appears to be similar across the sites, yet the growth level over time is higher in the drier and warmer climate than in the temperate zone. When compared with previous growth periods, growth models reveal that oak trees grew more than expected during the last decades. A recent setback in growth can be observed, although growth is still higher than the model predicts. By focusing on the short-term reactions of the trees, distinct drought events and periods were discovered. In each climatic region, similar growth reactions developed after drought periods. A decline in growth rate occurred in the second or third year after the drought event. Oaks in South Africa are currently exposed to a warmer climate with more frequent drought events. This climatic condition is a future prediction also for Europe. In view of this climate change, we discuss the consequences of the long- and short- term growth behaviour of oaks grown in the climate of South Africa for a tree species selection that naturally occurs in Europe.


Introduction
Knowledge of the resilience and resistance of tree species is crucial for forestry in a changing climate. In boreal and temperate climates where the rotation may last a century or more, species choice is a far-reaching decision, which should anticipate the species' performance under climate change. In order to generate respective knowledge, the growth and vitality of trees is often studied along and resilience of less adapted species like spruce and beech. Pretzsch et al. [18] found a ranking pertaining to the growth in the drought periods of 1976 and 2003. When compared to spruce and beech, oak had the top ranking position for growth. This underpins the relative high resilience of oak to drought in its natural range. The considered oak, like other tree species, showed remarkable growth trends in Europe [19,20]. Data was collected from 14 long-term experimental plots for almost one and a half centuries. From this data, it was concluded that tree size and stand parameters of oak trees currently develop much faster than in the past. This is highly relevant for forestry in Central Europe. Thus, from the management perspective, particular threshold sizes are reached decades earlier than when compared with the past. Due to the accelerated stand development, considering the same stand age, stem numbers per unit area are presently lower than in the past, while at the same time, stand density in terms of volume is higher [19,21]. In addition, the level of the tree growth rate vs. tree size allometry increased significantly.
Within the natural distribution range, a broad range of parameters regarding the growth of oak trees have already been investigated in several studies [22][23][24][25][26]. To identify the climatic potential oak growth has and where the geographical limit of oak trees lies, we extended the site spectrum of the species beyond its natural range. We pooled growth data from various sites within Central Europe (Germany and Italy) and compared it to the growth data of oaks from various sites in the Western Cape, South Africa. This reflects the growth occurring outside the natural range of oaks [27]. To gain knowledge on how oak performs within and outside its natural distribution area, we formulated the following research questions: 1.
How do oaks perform in regard to their general growth, both within and beyond their natural range, over a long-term basis? 2.
Has the oak growth revealed a change in growth dynamic during the last decades? 3.
Are there site-specific differences in the growth reaction of oaks to drought events, particularly in terms of resistance and resilience?

Study Areas and Climatic Conditions
The study sites are located in Southern Germany, Northeast Italy, and Western Cape, South Africa. Each study site (hereafter termed "experiment") consists of different research areas within each country. For Germany and Italy, four sites were selected, while nine sites were designated in South Africa.
The four sites in Southern Germany are situated between 9 • 26 E-11 • 48 E and 48 • 11 N-49 • 57 N, covering North to South Bavaria (Figure 1a). A detailed description of these experiments can be found in Pretzsch et al. [28]. The German sites represent the temperate climate zone, where the annual mean temperature is within a range of 7.0-8.5 • C and the mean precipitation within 660-960 mm. The climate data for this data set was obtained from weather stations close to each of the four sites ( Figure 2a) [29].
The chosen Italian sites were located in the Veneto region (Northeast Italy), with coordinates ranging from 45 • 52 N-45 • 38 N to 11 • 55 E-12 • 07 E (Figure 1b). This experiment represents the transition zone from warm-temperate to subtropical climate. Mean annual temperature ranges between 12.9-14.1 • C and mean annual precipitation varies around 590-800 mm/year ( Figure 2b). Climate data for the Italian research sites were available in a monthly resolution from climate stations close to the research plot.
Both, the German and Italian study areas are located within the natural range of the species Quercus robur L./Qu. petraea Liebl. [27,30].
The study areas in the Western Cape, South Africa, represent growing conditions that are not typified by the natural occurrence of the chosen oak species. Research locations were selected within the Western Cape Province, in the districts of Cape Town and the Cape Winelands (coordinates 33 • 56 S-34 • 08 S/18 • 28 E-18 • 97 E). These sites represent oak growth in the Mediterranean climate zone of the Southern Hemisphere ( Figure 1c). Mean annual temperature is substantially higher than in the natural distribution area of these    Both, the German and Italian study areas are located within the natural range of the species Quercus robur L./Qu. petraea Liebl. [27,30].
The study areas in the Western Cape, South Africa, represent growing conditions that are not typified by the natural occurrence of the chosen oak species. Research locations were selected within the Western Cape Province, in the districts of Cape Town and the Cape Winelands (coordinates 33°56′ S-34°08′ S/18°28′ E-18°97′ E). These sites represent oak growth in the Mediterranean climate

Tree-Ring Data
For the experiments in Germany and Italy, pure oak stands were chosen. From each research site, we sampled solely dominant trees. The sampled trees in South Africa are primarily solitary grown oaks in rural areas around Cape Town, Stellenbosch, and Paarl. In total, 139 trees were sampled; 66 oaks from Germany, 20 from Italy, and 53 from South Africa were used for dendrochronological measurements and analyses.
Standard measurements of tree height and diameter at breast height (DBH) were recorded for all investigated trees. Trees were cored from two directions with an angle of 90 • using a 5 mm increment corer (Haglöf, Langsele, Sweden). Thereafter, the cores were polished with progressively finer sandpaper from grit size 180 up to 800. Annual ring widths were measured with the Lintab measuring system using the TsapWin software (Rinntech, Heidelberg, Germany) [31]. Crossdating and synchronization of ring-width data was accomplished by the software TsapWin [31] and COFECHA [32,33] using standard dendrochronological methods. Chronologies of each experiment were built from a minimum sample depth of five trees. We removed the biological age trend in the ring-width data. A double de-trending procedure was applied to all series by using the Dendrochronology Program Library in R (dplR-package) [34,35]. The resulting index series contains only year-to-year variability associated with fluctuations in climate [36][37][38]. In this procedure, the raw ring-width series were conservatively standardized by adjusting a negative exponential or a linear regression function with negative slope. Thereafter, the resulting index values were de-trended by a smoothing cubic spline with a 50% frequency-response cut-off equal to 67% of the series length to calculate dimensionless ring-width indices (RWIs). To produce trend-eliminated time-series of the ring-width indices, in a final step, autocorrelation was removed using autoregressive (AR) models and the series were averaged using Tukey's biweight robust mean. This reduces bias caused by extreme values. Through this procedure, the long-term growth trends in the raw tree-ring series were removed. Mean growth rate average and standard deviation, as well as expressed population signal (EPS), mean sensitivity (MS), and mean inter-series correlation (Rbar) were calculated for the chronologies using the dplR package in R [34,39,40]. EPS was used to test the population strength's representation of the resulting tree-ring chronologies [41]. This was calculated on a 50-year running window with an overlap of 25 years [36,42] for the common period of each experiments' tree-ring series. Mean sensitivity (MS) depends on the year-to-year variability and was employed as a measure for variability in the tree-ring series [38,40,43]. The mean inter-series correlation (Rbar) is a metric for agreement between all tree ring series to be averaged to a site chronology. It was calculated to assess the signal strength of the site chronologies [36,42,44] over the common period of each experiments' tree-ring series.

Performance of Growth Level within and beyond the Natural Range
Many studies have reported significant changes of tree growth due to environmental changes during the last decades [19,20,[45][46][47][48]. Consequently, we analysed our data for recent growth trends by fitting a biologically plausible age-increment model to all tree ring data up to 1980. We extrapolated the fitted model beyond this cut-off and compared its predictions to the observed tree increments after 1980. We chose this year as a cut-off due to the length of the tree-ring data sets Therefore, we achieved a considerable length of tree-ring data from both periods (before and after 1980) for our long-term analysis. Assuming, in the null hypothesis, that variables have no influence on the age-increment relationship, we chose the Hugershoff equation for the age-increment model. The Hugershoff equation is, in addition to being biologically meaningful [49], easy to linearize and can therefore be handled with standard regression methods. We chose the individual tree basal area growth iba (m 2 /year) as our measure for tree increment because, in contrast to diameter increment, it directly represents information about biological production [50]. Applied to iba, the Hugershoff model reads as follows (Equation (1)): where t is tree age and a 0 , . . . , a 2 are constants (a 0 , a 1 > 0, a 2 < 0). In linearized form it can be written as (Equation (2)): ln(iba) = ln(a 0 ) + a 1 · ln(t) + a 2 ·t, This linearized form was fitted to all observations of tree age and basal area increment up to 1980 using the following linear mixed effects model (Equation (3)): where the indexes i, j, k, and l represent the country, the study site, the individual tree, and an individual observation, respectively. The fixed effect parameters are β 0 , β 1 , and β 2 . b i , b ij , b ijk , c i , and d i are normally distributed random effects with an expected value of 0, and ε ijkl are I.I.D. errors (independent and identically distributed errors). By adding to the fixed effect parameter β 0 , the random effects b i , b ij , and b ijk allow for a country-, study-site-in-country-, and a tree-in-study-site-specific scaling of the Hugershoff model. The random effects c i and d i in addition permit a general country-specific curve shape by adding to the fixed effect parameters β 1 and β 2 . This model structure enabled us to perform biologically plausible model fits and predictions even for trees with short tree ring series, where simple individual tree model fits would not lead to useful results.

Long-Term Growth Behaviour and Recent Growth Trends
For all observations after 1980, we calculated the expected basal area incrementsîba ijkl by applying the fitted model from Equation (3) including all fixed and random effects. Dividing the actual basal area increments riba ijkl by the expected ones, we obtained Equation (4): This quantifies, for each single tree and for each year in the extrapolation period, how much the actual increment deviates from the one that would be expected if growth conditions were the same as in 1980 and earlier. In order to test for significant country-specific patterns in riba, we applied the following generalized additive mixed effects model (Equation (5)): The indexes in Equation (5) have the same meaning as in Equation (3); the same is true for the fixed effect β 0 and the random effects b ij and b ijk , and for the I.I.D. errors ε ijkl . The functions f 1 , f 2 , and f 3 are non-parametric smoothers as functions of the calendar year (year); each one applies to a specific country-this is indicated by the binary coded variables a 1 , a 2 , a 3 , which are defined as follows (Equations (6a)-(6c)): The fixed effects β 1 and β 2 allow for estimating country-specific base levels of riba in combination with β 0 , because they are associated with the binary variables a 2 and a 3 , which represents a dummy coding. Therefore, β 0 is the base-level prediction for Germany while β 0 + β 1 and β 0 + β 2 is the same prediction for Italy and South Africa, respectively.
All statistical computations were performed with the software R version 3.2.2 [51], specifically the packages lme4 [52] and mgcv [53]. For a tree growing under constant conditions, an allometric relationship (see below, Equation (7)) ln(y) = a + b· ln(x), (7) can be expected between its size x and its size growth y [54]. Consequently, any change in the parameters a and b of this relationship would indicate a change in the tree's growing conditions. We utilized this concept for indirect testing of whether the growth conditions of the oaks were different in the years before and after 1980. Taking the tree basal area (ba) at a specific point in time as the size variable and the basal area increment in the subsequent year (iba) as the size growth variable, we formulated the following linear mixed-effect model (Equation (8)): The indexes i, j, and t represent the study site, the individual tree, and a specific point in time (calendar year), respectively. The variable period indicates whether an observation is from a year before (period = 0) or after 1980 (period = 1). The model's fixed effect parameters are α 0 , α 1 , β 0 and β 1 , while a i , a ij , b i and b ij are random effects on the study site level (index i) and on the tree-in-study-site level (indexes ij). These random effects are assumed to be normally distributed with an expectancy of 0. The same assumption can be made for ε ijt which represents I.I.D. errors. If α 1 and β 1 or both significantly differ from zero, this would indicate different allometric relationships before and after 1980.
The random effects in the model have important functions: a i and b i cover specific effects that result from the given site level combinations and thus allow for significant overarching relationships to be identified. The role of a ij and b ij is to cover autocorrelation effects due to repeated measurements from each tree. We considered it important to not only include random effects on the intercept of the allometric relationship (random effects a i , a ij ) but also on the slope (random effects b i , b ij ), which might differ considerably between study sites [55].
We fitted the mixed linear model (Equation (8)) for the whole dataset as well as separately for each climate zone (Temperate Germany, Warm-Temperate/Subtropical Italy, and Mediterranean South Africa). This allowed us to test for a general period effect (before and after 1980) on the size-growth-relation of oaks on the one hand. On the other hand, the zone-specific fits informed us in more detail without ending up with an over-complicated overall model. The range of tree basal areas in both periods, before and after 1980, largely overlaps. Thus, biased results due to unbalanced data had not to be suspected.

Reactions of Oak Growth on Drought Events
To ascertain whether extreme event years are appearing in all experiments and the differentiation in strength of these events, we applied pointer year analyses to the data sets [56,57].
In a first step, all raw time-series were standardized, for this analysis, with their 13-year moving average [57]. In a second step, these time-series were normalized to a standard deviation of one and a mean of zero over the period 1956-2004 for Germany, 1946-2009 for Italy, and 1903-2005 for South Africa [56]. This procedure enables an identification of growth anomalies to the probability density function of the standardized normal distribution and facilitates a clear assignment of the normalized Cropper values (C) as follows: weak for (C) > +1.000/−1.000, strong for (C) > +1.280/−1.280, and extreme for (C) > +1.645/−1.645 [56]. Secondary to climatic factors, additional matters exist that can cause pointer years, including fruiting, masting behaviour, and insect outbreaks.
This study employed the standardized precipitation evapotranspiration index (SPEI) as a drought indicator to determine the occurrence and extent of drought events and periods in South Germany, Northeast Italy, and the Western Cape, South Africa. To analyse the reaction on drought events comparing oak growth within and outside their natural distribution area, we calculated a drought index for the three experiments (Germany, Italy, and South Africa) using the available climate data from each country/experiment. Since drought is determined by temperature, rainfall, and evaporative changes, it is important to consider each variable individually in order to accurately analyse its impact. Vicente-Serrano et al. [58] proposed a drought index: SPEI, which depends on the potential evapotranspiration (PET). SPEI accounts for the effect of temperature variability in drought monitoring, and it can be computed at different time-scales. The process of calculating SPEI is based on a monthly balance between precipitation and PET, which is adjusted using a three-parameter log-logistic distribution to take into account common negative values [59]. PET was calculated using the Thornthwaite [60] equation. It is a simple method used to compute monthly potential evapotranspiration, when, as is the case here, only temperature data are available. For this study, we calculated the SPEI with a time scale of 6 months for Germany from 1947-2011, for Italy from 1955-2012, and for South Africa from 1945-2011.
The resulting drought years attained from the SPEI evaluations for each experiment were further analysed. For this, we applied superposed epoch analyses (SEAs) [61,62] to the ring-width data of the measured oaks in order to examine whether drought events have an influence on the growth of oak and how this varies in different climates. One point of interest was the recovery time after drought events with different intensities. In a first step, we removed the biological age trend in the ring-width data, following the same method as described under Section 2.2. In a second step, we averaged the difference between the mean ring-width index of each tree and the corresponding average value for the, in total, 9 years (4 years before the event year, the event year itself, and 4 years after the event year), respectively. Bootstrap resampling was used to determine whether RWIs for these years were significantly different to a random set of other years (lag: +1 year) from the data set. Bootstrapping conduces in the significance test for the derivation of confidence intervals (p < 0.05), for which the dendrochronological package dplR was used [39].

Tree-Ring Data and Their Basic Statistics
The selected trees at the four German, the four Italian, and the nine South African research sites covered a representative scope of age, height, and diameter at breast height (DBH) ( Table 1). Tree-ring series from all sites of each country were successfully cross-dated, resulting in the establishment of three chronologies. For each experiment's chronology, the basic tree ring statistics are shown in Table 1, according to their equivalent periods.  Matching the extent of increment growth (mean growth rate, Table 1) between the three research areas, we observed higher ring-width variability for the trees in Italy and South Africa compared to the trees from Germany ( Figure 3). Compared with trees in Germany, which had a maximum growth rate of 4.75 mm a −1 , oaks in Italy showed a growth rate of 10.62 mm a −1 and oaks in South Africa a rate of 12.48 mm a −1 . Differences in ring-width variability could also be seen when matching the level of raw ring-width data ( Figure 3) between the three experiments.
Mean sensitivity (MS) and mean inter-series correlation (Rbar) values were similar at the German and South African research sites and range within the threshold values given by Speer [63], whereas at the Italian research sites, MS showed the highest value but Rbar the lowest. EPS values, calculated for the common periods, respectively, revealed the highest value for the German sites. The EPS value for South African oak trees was just a little bit lower, whereas Italian oaks once again showed the lowest value (Table 1). The dark grey graph represents the averaged raw value chronologies of each experiment's sites. On the second Y-axis, the sample depth is pictured as an area-plot.

Evaluation of Long-Term Growth Trends
Performance of Growth Level within and beyond the Natural Range A first overview comparing the growth level of oaks within and outside their natural distribution area is shown in Figure 4 (derived from model Equation (3); model parameters can be found in the Appendix A Table A1). It clarifies whether long-term growth trends have different characteristics.
Comparing the three experiments, the basal area increment development over age is at the lowest level for the oaks in Germany. The Italian oak trees showed the fastest increase in basal area development over age. South African oaks start at the highest level but then increase only slightly with age, i.e., they have a similar development pattern with increasing age as the German oaks. The dotted lines in Figure 4 represent the model Equation (3) results for the regions within the experiments, respectively. South Africa and Italy revealed higher site-specific variability than Germany did. To study the long-term growth of oaks, we proposed a cut-off of the data set in 1980. Thus, we looked at a basal area growth development for the period before and after 1980 for all experiments ( Figure 5, derived from model Equation (8); model parameters found in the Appendix A Table A2). The growth trend follows a similar pattern at all sites but with different magnitudes. It appears small trees are growing faster than historical trends until they reach a certain size threshold, and larger trees tend to grow more slowly than historical trends ( Figure 5). At all three experiments the slopes in basal area increment are significantly different when the periods before and after 1980 are contrasted. German oaks have the least growth change before and after 1980. On the Italian sites, the incremental gap between growth before and after 1980 is extensive. Oak growth has the steepest slope before 1980 starting at an individual tree basal area (m 2 ) under 0.01 m 2 until 0.12 m 2 and flattens in the period after 1980. Looking at the growth of oak trees in South Africa, the slope between basal area increment, before and after 1980, differs considerably. The crossover point (before/after 1980) is at about 0.06 m 2 in Europe but at 0.20 m 2 in South Africa. As such, moderate-sized trees (ba = 0.10 m 2 ) in South Africa currently grow faster, as opposed to the same size of trees in Europe, which grow slower than historic rates.
Significant differences between the experiments are found by observing whether oak growth is of a variable pattern by focusing on the long-term development. Figure 6 (derived from model Equation (5); model parameters in the Appendix A Table A3 and Figure A1) shows the ratio values (after 1980) of the individual tree basal area increment growth, comparing real to expected growth, between the three experiments Germany-Italy-South Africa.   Compared to the other study sites, growth starts here on the highest level but decreases after 1995. Trees in South Africa indicate a slightly higher growth after 1980 than expected from the historic model. An increase in growth is visible from around the year 1990 until 1997-1998.

Reactions of Oak Growth on Drought Events
Pointer year analyses were applied to all data sets to examine the impact of extreme events on the growth of oaks grown in different climates. The pointer year totals for Germany   Oaks at the German research sites had the strongest and most extreme negative pointer years in the time span 1900-1930. Two extreme positive and three extreme negative pointer years could be identified (Figure 7). On average, there are 1-2 negative pointer years and 2-3 positive pointer years (from weak to extreme per definition) that occur every century. On Italian research sites, no extreme pointer year was found in our data. However, from 1957 to 2012 there are two strong negative and three strong positive pointer years (Figure 7). South African oaks showed an equal ratio of negative to positive pointer years. Just one extreme positive pointer year was identified, in addition to six strong positive and five strong negative ones (Figure 7). On all sites, negative and positive pointer years are evenly balanced.
Identified pointer years were matched with their congruence to a drought index with annual resolution. Pointer years did not reveal an exact correlation to the identified drought years.
To identify drought years at each experiment, we chose the standardized precipitation evapotranspiration index (SPEI). According to Potop et al. [59], resulting negative and/or positive index values can be classified according to seven classes of SPEI categories, which range from moderate wet/dry over normal to severe and extreme wet and/or dry. In our approach, we chose two of these categories to differentiate extreme from severe droughts. A SPEI value of −2 to −1.5 is classified as severe drought and values ≤ −2 as extreme drought events ( Figure 8). Oaks at the German research sites had the strongest and most extreme negative pointer years in the time span 1900-1930. Two extreme positive and three extreme negative pointer years could be identified (Figure 7). On average, there are 1-2 negative pointer years and 2-3 positive pointer years (from weak to extreme per definition) that occur every century. On Italian research sites, no extreme pointer year was found in our data. However, from 1957 to 2012 there are two strong negative and three strong positive pointer years (Figure 7). South African oaks showed an equal ratio of negative to positive pointer years. Just one extreme positive pointer year was identified, in addition to six strong positive and five strong negative ones (Figure 7). On all sites, negative and positive pointer years are evenly balanced.
Identified pointer years were matched with their congruence to a drought index with annual resolution. Pointer years did not reveal an exact correlation to the identified drought years.
To identify drought years at each experiment, we chose the standardized precipitation evapotranspiration index (SPEI). According to Potop et al. [59], resulting negative and/or positive index values can be classified according to seven classes of SPEI categories, which range from moderate wet/dry over normal to severe and extreme wet and/or dry. In our approach, we chose two of these categories to differentiate extreme from severe droughts. A SPEI value of −2 to −1.5 is classified as severe drought and values ≤ −2 as extreme drought events ( Figure 8). Climate data showed no extreme drought years according to SPEI at the German research areas. Severe dry years were identified in 1953, 1972, 1997, and 2003. Of all the experiments, the Italian study sites had the most severe and extreme dry years. These were 1982,1986,1990,1993,2003, as well as 2012. Looking at the selected South African study areas, the years 1960, 1973, and 1994 could be identified as extreme drought events. Furthermore, several severe drought events were evident from 1960-2010 ( Figure 8). Besides the drought events, we also identified drought periods for our study sites in South Africa (e.g., 1997-2000) whose effects were also reflected in the long-term analysis. For all experiment sites, the amount of drought events and/or periods increased from 1980 onwards, which reinforced our decision to set 1980 as a cut-off date for our long-term analysis ( Figure 5 and 6). Calculations of the SPEI indicated several drought years per research area, which could be used in the superposed epoch analyses (SEA). For calculating SEA, we chose drought years that appeared in the SPEI analyses as either severe or extreme drought years. Furthermore, to minimize bias, we excluded years that were less than four years apart. Drought years occurring outside the common period (Table 1) of the research site's tree-ring chronology were also excluded. In total, three to four drought years per experiment were incorporated into the analyses. For the oak chronologies from Germany, we chose the drought years 1972, 1997, and 2003; for the Italian sites, we selected the years 1982, 1986, 1990, and 2003; for South African oak trees, the drought years 1960, 1973, and 1994 were selected.
Oaks at each site showed a reduction and recovery phase after drought, but intensity and time were different between the experiments (Figure 9). All sites displayed the highest growth reduction two years after the drought event. This fact is most pronounced in the German experiment. Oaks growing on South African sites had the longest recovery time. Climate data showed no extreme drought years according to SPEI at the German research areas. Severe dry years were identified in 1953, 1972, 1997, and 2003. Of all the experiments, the Italian study sites had the most severe and extreme dry years. These were 1982,1986,1990,1993,2003, as well as 2012. Looking at the selected South African study areas, the years 1960, 1973, and 1994 could be identified as extreme drought events. Furthermore, several severe drought events were evident from 1960-2010 ( Figure 8). Besides the drought events, we also identified drought periods for our study sites in South Africa (e.g., 1997-2000) whose effects were also reflected in the long-term analysis. For all experiment sites, the amount of drought events and/or periods increased from 1980 onwards, which reinforced our decision to set 1980 as a cut-off date for our long-term analysis ( Figure 5 and 6). Calculations of the SPEI indicated several drought years per research area, which could be used in the superposed epoch analyses (SEA). For calculating SEA, we chose drought years that appeared in the SPEI analyses as either severe or extreme drought years. Furthermore, to minimize bias, we excluded years that were less than four years apart. Drought years occurring outside the common period (Table 1) of the research site's tree-ring chronology were also excluded. In total, three to four drought years per experiment were incorporated into the analyses. For the oak chronologies from Germany, we chose the drought years 1972, 1997, and 2003; for the Italian sites, we selected the years 1982, 1986, 1990, and 2003; for South African oak trees, the drought years 1960, 1973, and 1994 were selected.
Oaks at each site showed a reduction and recovery phase after drought, but intensity and time were different between the experiments (Figure 9). All sites displayed the highest growth reduction two years after the drought event. This fact is most pronounced in the German experiment. Oaks growing on South African sites had the longest recovery time.
When comparing the three experiments, oak chronologies from Germany showed the strongest growth setback after a drought event. In the drought year itself, decline was minimal on all sites. In the first and second year after the drought, the majority of the trees experienced a setback in growth, whereas oaks on the Italian research sites showed the least reaction after a drought event. Based on the chosen drought years, oaks on the South African study sites indicated the strongest decline in the second year after a drought (Figure 9). Results revealed a high tolerance of water deficits for oaks on all research sites, as none of the experiments' trees showed significant growth decline in the drought year itself nor in the following years. Mean RWI departures were also calculated for each site within the three experiments, these results are plotted in the Appendix A (Table A4). When comparing the three experiments, oak chronologies from Germany showed the strongest growth setback after a drought event. In the drought year itself, decline was minimal on all sites. In the first and second year after the drought, the majority of the trees experienced a setback in growth, whereas oaks on the Italian research sites showed the least reaction after a drought event. Based on the chosen drought years, oaks on the South African study sites indicated the strongest decline in the second year after a drought (Figure 9). Results revealed a high tolerance of water deficits for oaks on all research sites, as none of the experiments' trees showed significant growth decline in the drought year itself nor in the following years. Mean RWI departures were also calculated for each site within the three experiments, these results are plotted in the Appendix A (Table A4).

Discussion
Changing climate conditions will likely affect tree growth and have motivated forest science to identify possible pathways of growth reactions to facilitate adaptation. During the last decades, a wide range of studies have focused on the growth dynamics of trees within their natural habitat [23,65,66]. Other studies investigated tree growth response to climate change by looking at tree growth in regions that currently exhibit climate conditions that are expected for the area of interest under future climate [67][68][69][70]. Both approaches have the limitation that the real flexibility of a tree species against climate conditions is not fully covered. So far, species growing at conditions that are distinct from their natural range conditions rarely have been investigated. However, such studies could provide an idea of how native trees grow if climatic conditions are changing. Pretzsch et al. [71] showed, for trees growing under urban conditions, that growth could even be fostered due to higher temperatures. In our study, we included oaks growing at sites that exhibit natural range growing conditions as well as strongly diverging conditions. We have the advantage to revert to a data set of more than 100-years tree growth within and outside the trees natural range.

Discussion
Changing climate conditions will likely affect tree growth and have motivated forest science to identify possible pathways of growth reactions to facilitate adaptation. During the last decades, a wide range of studies have focused on the growth dynamics of trees within their natural habitat [23,65,66]. Other studies investigated tree growth response to climate change by looking at tree growth in regions that currently exhibit climate conditions that are expected for the area of interest under future climate [67][68][69][70]. Both approaches have the limitation that the real flexibility of a tree species against climate conditions is not fully covered. So far, species growing at conditions that are distinct from their natural range conditions rarely have been investigated. However, such studies could provide an idea of how native trees grow if climatic conditions are changing. Pretzsch et al. [71] showed, for trees growing under urban conditions, that growth could even be fostered due to higher temperatures. In our study, we included oaks growing at sites that exhibit natural range growing conditions as well as strongly diverging conditions. We have the advantage to revert to a data set of more than 100-years tree growth within and outside the trees natural range.

General Growth Level and Growth Trends
For the growth of oaks, the current study points out a lower variability in ring-width data for German oaks compared to the oak trees growing in the examined study sites in Italy and South Africa. Hereby, the year-to-year variability ( Figure 3) and the general level of growth appears to be accelerated (Figure 4) for the trees that are growing in a warmer climate with possibly longer drought periods. The growth level of oaks growing at average temperatures of 11 • C (Italy) or 16 • C (South Africa) but at similar precipitation regimes as in Germany show a two and four times higher productivity respectively than those in Germany. Growth trajectories for all sites are far beyond those predicted by the existing yield table. Considering that all our results are based on a single tree level, we cannot provide any statements for stand level, where smaller trees would also be included. Higher temperatures, either through climate change or through generally effective climate conditions (Italy and South Africa), might explain the departure of growth trajectories from the yield table, which was developed for German growing conditions ( Figure 10). Figure 10. Observed diameter at breast height (DBH) over tree age (black lines) for study sites (a) Germany; (b) Italy; and (c) South Africa respectively. DBH is compared with the mean diameter range (grey shaded area), according to the yield table for oak (heavy thinning, site classes I-III) by Jüttner [72]. Yield table values are pictured on different scales for each country.
Even if we assume that for the study sites in South Africa the trees had more favourable conditions through possible tending activities like irrigation and reduction of competing vegetation, the long-term growth development seems to be more related to prevailing climate conditions. Relatively warm spring seasons, higher temperatures throughout the year, as well as an extended vegetation period potentially foster growth.
At all sites, a more rapid growth during the last decades (after 1980) was observed. Changes in size-growth relations ( Figure 5) suggest that younger trees are currently growing quicker in diameter, but increment rates for older trees seem to be damped compared to before 1980. Consequently, for the trees investigated in our study, the general growth curve seems to be modified by one that has a steeper slope and an earlier inflection point. (Figure 11). This change in growth pattern appears to be most pronounced within oak trees in South Africa.
The accelerated growth can be traced by the comparison of recent growth level with that before 1980. Trees are currently performing better than they did before 1980 ( Figure 6). Decreasing relation over time between modelled growth rates (before 1980) and observed growth rates, as visible from all three countries from 2009-2010 onwards, might be related to approaching growth curves through aging of the sampled trees. For the sampled oak trees from South Africa, the drop in the curve found around 1997-1998 could be exaggerated by the coinciding Super El Niño event in 1997-1998. That Super Nino event had a huge negative impact on the rainfall and temperatures in the Western Cape, South Africa [73]. A protracted period of below average levels in the SPEI data in South Africa approximately occurred for the period 1997-2007. This pattern has not occurred elsewhere in our data. The same goes for Germany (2003)(2004)(2005)(2006)(2007)(2008). Drought has never persisted that long in our data set. In Italy, there does not appear to be a protracted dry period that coincided with the drop in modelled trends (2000)(2001)(2002)(2003)(2004)(2005)(2006)(2007)(2008)(2009)(2010). The deferred change in growth at the German sites might be influenced by high atmospheric deposition during the 1960's to 1980's, which limited growth all over Central Europe.
Similar global change triggered changes in growth patterns have also been described by Pretzsch et al. [19] for Central Europe; Kauppi et al. [74] for boreal forests; and Fang et al. [75] for forests in Japan. trends (2000-2010). The deferred change in growth at the German sites might be influenced by high atmospheric deposition during the 1960's to 1980's, which limited growth all over Central Europe.
Similar global change triggered changes in growth patterns have also been described by Pretzsch et al. [19] for Central Europe; Kauppi et al. [74] for boreal forests; and Fang et al. [75] for forests in Japan.

Reactions of Oak Growth on Drought Events/Periods
Oak trees reveal sustainable growth reduction after severe drought events. This general growth behaviour is observed in all analysed climate zones. Although growth is not limited much during the drought year itself, growth reduction is even stronger in the subsequent years. The heaviest growth reduction over time was observed at German sites and least at Italian sites. Oaks growing in South Africa showed the longest period for recovery from drought. Being a ring-porous tree species, oak is most likely categorised as an anisohydric type. Anisohydric plants strive to regulate stomatal conductance to keep photosynthetic processes alive, even through drought periods [76]. However, anisohydric traits may lead to a longer drought memory due to high investments to sustain photosynthetic activity, thus leading to a hangover of drought effects. We already know from several studies from Middle Europe [22,25,26,77,78] that oaks react with a delay on drought events. The relative high tolerance of oaks [79] towards water deficits is confirmed by the fact that trees investigated in our study showed no significant growth decline during drought years against average growth. Even for much warmer regions growth reaction of oak on drought does not differ significantly in intensity and duration from that in German climates. These results lead to the assumption that growth of oak is much more determined by temperature than by water supply.

Conclusions
Currently, in Central European forests, the two main oak species are decreasing in proportion due to the tendency towards close-to-nature forestry favouring shade tolerant species. Light demanding tree species, such as oak or pine, are often outcompeted in continuously covered forests because of their need for light, especially in the early development state. However, oak seems to have the potential to cope with warmer climates, which are likely to arrive in Central Europe in the future [69,80,81].
This study shows high and accelerated growth of oak (Figure 12a), even beyond its natural range, which is probably related to changing environments. However, prevailing climates in the region of the natural distribution of a tree species are not necessarily the limiting factor for its growth. Rather, a higher growth level is possible and therefore generally assumed climate-growth relations may need to be reconsidered (Figure 12b). Even assumptions about stress-related growth reaction might be fragmentary when drawn on natural range conditions (Figure 12c).

Reactions of Oak Growth on Drought Events/Periods
Oak trees reveal sustainable growth reduction after severe drought events. This general growth behaviour is observed in all analysed climate zones. Although growth is not limited much during the drought year itself, growth reduction is even stronger in the subsequent years. The heaviest growth reduction over time was observed at German sites and least at Italian sites. Oaks growing in South Africa showed the longest period for recovery from drought. Being a ring-porous tree species, oak is most likely categorised as an anisohydric type. Anisohydric plants strive to regulate stomatal conductance to keep photosynthetic processes alive, even through drought periods [76]. However, anisohydric traits may lead to a longer drought memory due to high investments to sustain photosynthetic activity, thus leading to a hangover of drought effects. We already know from several studies from Middle Europe [22,25,26,77,78] that oaks react with a delay on drought events. The relative high tolerance of oaks [79] towards water deficits is confirmed by the fact that trees investigated in our study showed no significant growth decline during drought years against average growth. Even for much warmer regions growth reaction of oak on drought does not differ significantly in intensity and duration from that in German climates. These results lead to the assumption that growth of oak is much more determined by temperature than by water supply.

Conclusions
Currently, in Central European forests, the two main oak species are decreasing in proportion due to the tendency towards close-to-nature forestry favouring shade tolerant species. Light demanding tree species, such as oak or pine, are often outcompeted in continuously covered forests because of their need for light, especially in the early development state. However, oak seems to have the potential to cope with warmer climates, which are likely to arrive in Central Europe in the future [69,80,81].
This study shows high and accelerated growth of oak (Figure 12a), even beyond its natural range, which is probably related to changing environments. However, prevailing climates in the region of the natural distribution of a tree species are not necessarily the limiting factor for its growth. Rather, a higher growth level is possible and therefore generally assumed climate-growth relations may need to be reconsidered (Figure 12b). Even assumptions about stress-related growth reaction might be fragmentary when drawn on natural range conditions (Figure 12c).
Considering these results, oak can play a crucial role when forests are to be adapted to climate change in Central European regions. Therefore, silvicultural approaches towards climate smart forestry should tap these climate stability features of oak by establishing the required light regime for this species in order to increase its share in European forests.
The study further underpins the potential of dendrochronological analyses of tree species growing beyond their natural range. Further studies (e.g., in arboreta, urban environments, or common garden experiments) may additionally exploit the growth potential of trees. Such a tree provides a so far neglected potential for climate change research. Considering these results, oak can play a crucial role when forests are to be adapted to climate change in Central European regions. Therefore, silvicultural approaches towards climate smart forestry should tap these climate stability features of oak by establishing the required light regime for this species in order to increase its share in European forests.
The study further underpins the potential of dendrochronological analyses of tree species growing beyond their natural range. Further studies (e.g., in arboreta, urban environments, or common garden experiments) may additionally exploit the growth potential of trees. Such a tree provides a so far neglected potential for climate change research.
Author Contributions: H.P. and E.U. conceived, designed and performed the experiments; D.P. and P.B. analysed the data; V.C. contributed to materials; D.P. wrote the paper and performed the experiments; H.P., P.B., and E.U. contributed to writing; B.D.T. and T.R. revised the manuscript and performed the experiments.

Conflicts of Interest:
The authors declare no conflicts of interest. Table A1. Model parameters for the linear mixed effect model Equation (3) displaying the performance of growth level (Figure 4). Parameter estimates for each experiment's fixed effects with standard errors (Std. Error) and significances (p) and the standard deviation (Std. Dev.) of random effects. Estimated country-specific values were added to the random effects, as they are needed for the curve display in Figure 4.

Conflicts of Interest:
The authors declare no conflicts of interest. Table A1. Model parameters for the linear mixed effect model Equation (3) displaying the performance of growth level ( Figure 4). Parameter estimates for each experiment's fixed effects with standard errors (Std. Error) and significances (p) and the standard deviation (Std. Dev.) of random effects. Estimated country-specific values were added to the random effects, as they are needed for the curve display in Figure 4.  Table A2. Model parameters for the linear mixed effect model Equation (8) displaying the change of size-growth relation ( Figure 5). Parameter estimates for each experiment's fixed effects with standard errors (Std. Error) and significances (p) and the standard deviation (Std. Dev.) of random effects.  Table A3. Model parameters for the generalized additive mixed effects model Equation (5) displaying the long-term growth behaviour and recent growth trends ( Figure 6). Parameter estimates for each experiment's fixed effects with standard errors (Std. Error) and significances (p) and the standard deviation (Std. Dev.) of random effects. The non-parametric smoothers as function of the calendar year are added to the fixed effects and are displayed in Figure A1 in the Appendix A.  Figure A1 0.0000 f2(year) see Figure A1 0.0000 f3(year) see Figure A1 0.0000