Allometry and Post-Drought Growth Resilience of Pedunculate Oak ( Quercus robur L.) Varieties

: This paper presents an analysis of the radial growth, tree dimensions, and allometry of three phenological pedunculate oak ( Quercus robur L.; QURO) varieties (early (E-QURO), typical (T-QURO), and late (L-QURO)), from a common garden experiment. We focused on the resistance and resilience of each variety to drought events, which occurred in 2012 and 2017, as well as their recovery potential during juvenile and mature growth phases, with the goal of clarifying how QURO drought sensitivity is inﬂuenced by tree phenology and growth stage. Our results indicate that E-QURO is more drought resistant, while T-QURO and L-QURO exhibit greater recovery potential after a drought event. Hence, typical and late QURO varieties are better prepared to withstand climate change. We also noted differences in the physical dimensions and the allometry of the studied QURO varieties. On average, 21-year-old QURO specimens from the analyzed stand are 9.35 m tall, have a crown width (CW) of 8.05 m, and a diameter at breast height (DBH) of 23.71 cm. Although T-QURO varieties had the greatest DBH and CW, they were shorter than E- and L-QURO, which are similar in height. T-QURO is also shorter relative to DBH, while L-QURO has a wider crown relative to tree height (TH). Intra-variety variations are higher than variations among half-sib (open-pollinated) families of each variety. Moreover, the adopted regression model provided a better ﬁt to the CW/DBH ratio than to TH/DBH and CW/TH.


Introduction
The selection of tree specimens with traits appropriate to the local environment is an important aspect of "climate-smart" forest management [1][2][3][4]. Due to its high intra-species and intra-provenance variations [5], pedunculate oak (Quercus robur L.; QURO) provides a valuable gene pool for the selection of superior, climate-resistant genetic material. In European-scale provenance trials, QURO radial growth performances showed a high level of heredity [5]. Furthermore, phenological varieties in a common garden experiment kept the phenological dynamics of the mother tree, with more than a month's difference between the leafing and flowering phenological phases [6] of the early and late varieties. Empirical evidence further indicates that wood anatomical characteristics are inherited from mother trees [7] and influence the water use efficiency (WUE) [8], which is crucial for tree fitness. These phenomena provide a sustainability in heredity of superior traits across an appropriately managed QURO forest gene pool.
QURO radial growth variations are also affected by climate conditions [9,10] and site characteristics [11,12]. Hence, intra-site comparisons can provide valuable evidence regarding historical environmental conditions. Recent drought events (21st century) have caused dieback in most European oak forests, which are also adversely affected by pest attacks [13,14]. As QURO is highly water-sensitive, both above-ground (e.g., precipitation

Site Description
The experimental trial that is the focus of our investigation is located in the lowland region of Serbia (southeastern Europe), close to the town of Sremska Mitrovica (45 • 0 19.076 N; 19 • 36 51.137 E) and within the Sava River catchment area. The plot is located at 85 m a.s.l., with a slope of <2%. The stand conditions and terrain characteristics are uniform. According to the SoilGrids (https://soilgrids.org, accessed on 15 May 2021) database, which is based on the World Reference Base classification, the soil in the area can be characterized as a gleysol. As the plot is located~10 km away from the river, it is not directly influenced by the anomalies in discharge. However, climate conditions significantly affect the analyzed trees.

Climate Conditions and Drought Events
Climate data for this study (spanning the 1950−2020 period) were sourced from the FORESEE database [34]. The obtained data were interpreted on annual and monthly timescales. Average temperature (TEMP; • C), sum of precipitation (PRCP; mm m −2 ), and the forest aridity index (FAI) served as climate indicators. Average temperature was interpreted at an annual scale (TEMP) and as averages for spring (March to May, TEMP MAM ) and summer (June to August, TEMP JJA ). Likewise, precipitation was interpreted as the annual sum of precipitation (PRCP) and calculated separately for spring (March to May, PRCP MAM ) and summer (June to August. PRCP JJA ). All meteorological data were interpreted as averages with standard deviations (SD).
Drought events were interpreted using the Standardized Precipitation Index (SPI) [35] for the 2005−2020 period. In the analyses, 6-and 12-month August SPIs were adopted as indicators of annual drought (SPI.12.Aug) and drought during the intensive vegetation period (SPI.6.Aug).

Trial Setup
During the second half of the 20th century, intensive phenotyping observations were made in different types of flooded and non-flooded QURO forests within the Danube and Sava River catchment areas, which are located in the northern lowland region of Serbia [6]. Based on those findings, 871 representatives of phenotypically determinate 21-year-old QURO specimens exhibiting excellent performance and fitness were included in our dataset. In this phenologically controlled trial, significant variations (up to 28 days) in leafing and flowering among early and late QURO varieties were noted [6].
Within the analyzed QURO gene pool, four QURO phenological varieties were determined and denoted as (i) early (Quercus robur var. praecox; E-QURO); (ii) typical (Q. robur var. typica; T-QURO); (iii) late (Q. robur var. tardiflora; L-QURO); and (iv) extremely late (Q. robur var. tardissima). Extremely late QURO is rare and its characteristics are not comparable with those of other analyzed QURO varieties, and thus was excluded from our study. A trial involving selected phenologically different QURO specimens was established in 2000 by Vojvodinašume Public Enterprise in collaboration with the Institute of Lowland Forestry and Environment (Novi Sad, Serbia). All of the planted one year-old QURO seedlings that were included in the trial originated from open-pollinated (half-sib) families taken from superior mother trees (plus trees) growing in natural forests across Serbia. During the 21-year period, only basic forest management techniques were implemented on the plot, namely dry tree removal and the clearing of lower branches. Planted QUROs were also phenologically checked. All selected QURO specimens were planted in a 6 m × 8 m raster and should be interpreted as individual trees in this study. Phenological varieties were organized in three blocks to form a single plot.

Dendrochronological and Biometric Analyses
Wood cores and biometric data were taken during the summer of 2020. Two cores per tree were taken (west and east orientation) at an angle of 90 • at breast height (~1.3 m), for a total of 120 cores from 60 trees (20 trees per variety). Only healthy, symmetrical QURO specimens free from any visible anomalies were included in the sample. To minimalize the influence of any variations in site conditions (soil properties, terrain, etc.), the sample was randomized and the representative trees were selected evenly across the entire plot, excluding trees located at the perimeter. After sampling, the collected wood cores were air-dried, fixed to wooden mounts, and sanded until tree rings were clearly visible. The prepared cores were scanned, and tree ring width (TRW), early wood width (EW), and late wood width (LW) were measured with an accuracy of 0.01 mm using CDendro 9.5 and CooRecorder 9.5 software (Cybis Elektronik & Data AB).
Biometric analyses (tree dimension measurements and allometric relations) were performed on all representative E-QURO, T-QURO, and L-QURO specimens, thus involving 871 trees. Tree dimension measurements included tree height (TH) and crown width (CW), both expressed in meters, as well as diameter at breast height (DBH; 1.3 m), expressed in centimeters, and tree trunk height (TTH) were uniform and counted at 2.4 m. We measured two radii at an angle of 90 • to the trunk of CW and DBH and their averages were interpreted in further analysis. Trees with DBH < 5 cm and damaged trunks or asymmetric or damaged crowns were excluded from the sample. Allometry analyses included TH/DBH, CW/DBH, and CW/TH ratios.

Statistical Analyses
Radial growth was interpreted using the cumulative basal area increment (BAI; cm 2 ), TRW, EW, LW, and their detrended values, i.e., tree-ring width index (TRWi), early wood index (EWi), and late wood index (LWi), respectively. In line with the approach adopted by Fritts [36], single detrending was adopted to remove the age trend using a modified negative exponential curve. Smoothed BAI, TRW, and TRWi chronologies were obtained from a generalized additive model (GAM) smoother using the equation (y~s(x)). Intraand inter-variety deviations in BAI, TRW, EW, and LW were interpreted using principal component analysis (PCA). The juvenile and mature tree growth phases were defined by applying the optimal segmentation model following radial growth chronology trends [37].
Variations in drought response across QURO varieties were interpreted using the resistance, resilience, and recovery indices [38][39][40]. For this purpose, drought years were defined with respect to SPI values. The influence of variety and half-sib families were interpreted using two factorial analysis of variance (ANOVA) to the unbalanced samples. The responses of different varieties to drought events were interpreted using the ANOVA honestly significant difference (HSD) Tukey post-hoc test, with p < 0.05 signifying statistical significance. Resistance, resilience, and recovery were calculated using the following equations: Resiliance = X post event year mean X 2 pre−event years (2) where X denotes the TRW value and 2012 and 2017 are considered drought event years. Tree biometry data were interpreted via descriptive statistics, while significance of difference was tested using two-way ANOVA with variety and half-sib as dependent variables and p < 0.05 indicating statistical significance. Allometry findings were interpreted by conducting ordinary least squares (OLS) linear regression analysis using the equation below: where parameters a and b denote intercept and regression slope, respectively, and se represents standard error of the regression model. All data processing was performed using R [41]. Chronology data were detrended using the "dplR" package [42]. The findings were presented graphically using the "ggplot2" package [43], while the "Car" package was used to calculate the descriptive statistics and conduct ANOVA with Tukey post-hoc test.

Climate Conditions, Stand Characteristics, and Drought Events
According to the World Climate Classification [44], climate at the study site is temperate continental to modified continental. August and January are the warmest and coldest months, respectively. Over the 70-year period from 1950 to 2020, TEMP and PRCP were 14.47 • C and 684.0 mm m −2 , respectively ( Table 1). During spring (TEMP MAM ) and summer (TEMP JJA ), the average temperatures were 14.75 • C and 24.47 • C, respectively, whereas the sums of precipitation during spring (PRCP MAM ) and summer (PRCP JJA ) were 165.1 and 207.8 mm m −2 , respectively. As shown in Table 1 Figure 1e), according to the Bjerknes's scale [45]. It is also worth noting that, the FAI values (16.81 and 12.88, respectively) in 2012 and 2017 were more than 200% higher than the optimal oak stand aridity due to the extreme droughts in these years, and thus, were extremely stressful and not conducive to QURO growth.
Several extremely wet and dry years were noted throughout the 2005−2020 period ( Figure 2). Based on 6-and 12-month accumulation periods, 2010, 2014, and 2019 were more humid years (SPI ≥ 1), whereas 2012 and 2017 were more arid years (SPI ≤ −1). According to the drought indices, 2007 and 2009 could also be considered as drought years.
However, because both the annual drought indices and those based on vegetation periods were the lowest for 2012 and 2017, these years are defined as the representative cases with regard to the growth performance of different oak varieties during drought. Likewise, in 2012 and 2017 were noted as different tree growth phases (juvenile and mature), which are included in this study.  Figure 1e) according to the Bjerknes's scale [45]. It is also worth noting that, the FAI values (16.81 and 12.88, respectively) in 2012 and 2017 were more than 200% higher than the optima oak stand aridity due to the extreme droughts in these years, and thus, were extremely stressful and not conducive to QURO growth.
Several extremely wet and dry years were noted throughout the 2005−2020 period ( Figure 2). Based on 6-and 12-month accumulation periods, 2010, 2014, and 2019 were more humid years (SPI ≥ 1), whereas 2012 and 2017 were more arid years (SPI ≤ −1). According to the drought indices, 2007 and 2009 could also be considered as drought years However, because both the annual drought indices and those based on vegetation periods were the lowest for 2012 and 2017, these years are defined as the representative cases with regard to the growth performance of different oak varieties during drought. Likewise, in 2012 and 2017 were noted as different tree growth phases (juvenile and mature), which are included in this study.   [45]. It is also worth noting that, the FAI valu and 12.88, respectively) in 2012 and 2017 were more than 200% higher than th oak stand aridity due to the extreme droughts in these years, and thus, were e stressful and not conducive to QURO growth.
Several extremely wet and dry years were noted throughout the 2005−20 (

Radial Growth Performance
Significant differences in radial growth performance were noted among QURO varieties and across the analyzed timespan (Figure 3a-c). Firstly, varieties differed with respect to the number of tree rings at breast height, with 14, 13, and 12 years recorded for T-QURO, L-QURO, and E-QURO, respectively. and are likely the results of the extreme droughts in 2012 and 2017.
The decrease in the radial increment after ~10th year of growth could be attributed to the end of the juvenile growth phase, which is the most productive phase in a tree's life Based only on the radial growth dynamics in the juvenile growth phase, L-QURO can be distinguished from T-QURO and E-QURO. The faster growth observed in L-QURO during its first years is important for its survival and competitiveness, as well as for faster nursery production of oak seedlings.
Inter-variety differences were examined by subjecting the data for the 2010-2019 period to PCA analysis, which revealed stronger deviations in BAI and TRW (Figure 4a,b) relative to EW and LW (Figure 4c,d). Based on the BAI and TRW chronologies, more pronounced differences were noted between L-QURO and E-QURO relative to T-QURO.  The decrease in the radial increment after~10th year of growth could be attributed to the end of the juvenile growth phase, which is the most productive phase in a tree's life. Based only on the radial growth dynamics in the juvenile growth phase, L-QURO can be distinguished from T-QURO and E-QURO. The faster growth observed in L-QURO during its first years is important for its survival and competitiveness, as well as for faster nursery production of oak seedlings.
Inter-variety differences were examined by subjecting the data for the 2010-2019 period to PCA analysis, which revealed stronger deviations in BAI and TRW (Figure 4a,b) relative to EW and LW (Figure 4c,d). Based on the BAI and TRW chronologies, more pronounced differences were noted between L-QURO and E-QURO relative to T-QURO.

Early and Late Wood Characteristics
Throughout the studied period, detrended EW and LW chronologies (Figure 3d showed some differences among analyzed QURO varieties. Although their patterns w similar, EW > LW was observed for L-QURO during its juvenile growth phase. Deviatio between the EWi and LWi chronologies were less pronounced in T-QURO and E-QUR The strongest EWi/LWi deviations were observed in L-QURO. According to the PCA results, EW and LW inter-variety differences were lower th those noted for TRW and BAI ( Figure 4). As shown in Figure 4c,d, the PC1 values for E and LW were nearly identical at 36.21% and 36.28%, while the PC2 values were sign cantly different (32.94% vs. 17.39%). It was not possible to define a clear pattern in ratios between early wood and late wood among phenological varieties.

Drought Resistance, Resilience, and Recovery Potential
The responses of QURO varieties to drought events (with 2012 and 2017 as test cas were interpreted via resistance, resilience, and recovery potential based on the droug indices and TRW data. The first drought event (2012) occurred during the QURO juven growth phase, whereas the second drought event (2017) occurred during the matu (adult and reproductive) growth phase.
Phenological variety and tree growth phase were significant factors in determini the drought sensitivity of QURO radial growth (Figure 5a-c). Higher resistance and bet

Early and Late Wood Characteristics
Throughout the studied period, detrended EW and LW chronologies (Figure 3d-f) showed some differences among analyzed QURO varieties. Although their patterns were similar, EW > LW was observed for L-QURO during its juvenile growth phase. Deviations between the EWi and LWi chronologies were less pronounced in T-QURO and E-QURO. The strongest EWi/LWi deviations were observed in L-QURO.
According to the PCA results, EW and LW inter-variety differences were lower than those noted for TRW and BAI ( Figure 4). As shown in Figure 4c,d, the PC1 values for EW and LW were nearly identical at 36.21% and 36.28%, while the PC2 values were significantly different (32.94% vs. 17.39%). It was not possible to define a clear pattern in the ratios between early wood and late wood among phenological varieties.

Drought Resistance, Resilience, and Recovery Potential
The responses of QURO varieties to drought events (with 2012 and 2017 as test cases) were interpreted via resistance, resilience, and recovery potential based on the drought indices and TRW data. The first drought event (2012) occurred during the QURO juvenile growth phase, whereas the second drought event (2017) occurred during the mature (adult and reproductive) growth phase.
Phenological variety and tree growth phase were significant factors in determining the drought sensitivity of QURO radial growth (Figure 5a-c). Higher resistance and better  Likewise, statistically significant intra-variety variations during both drought events were noted with regard to resistance, resilience, and recovery (Figure 5a-c). Although E-QURO exhibited the greatest resistance and resilience, these responses were significantly different from those of T-and L-QURO only in 2017.
Recovery values also differed across varieties and between drought events ( Figure  5c). In 2012, T-QURO exhibited the greatest recovery, followed by L-QURO and finally E-QURO. During 2017, which coincided with the mature growth phase, lower recovery performances were detected, with L-QURO demonstrating the best recovery in radial growth but lower resistance and resilience relative to the other two varieties.

Tree Dimensions and Allometry
On average, 21-year-old QURO specimens from the analyzed stand were 9.35 m tall and had crown widths of 8.05 m and DBH of 23.71 cm; however, significant variations were observed among QURO varieties. In total, 24, 25, and 27 half-sib families of early, typical, and late QURO varieties were included in this study. Based on two-way ANOVA, Likewise, statistically significant intra-variety variations during both drought events were noted with regard to resistance, resilience, and recovery (Figure 5a-c). Although E-QURO exhibited the greatest resistance and resilience, these responses were significantly different from those of T-and L-QURO only in 2017.
Recovery values also differed across varieties and between drought events (Figure 5c). In 2012, T-QURO exhibited the greatest recovery, followed by L-QURO and finally E-QURO. During 2017, which coincided with the mature growth phase, lower recovery performances were detected, with L-QURO demonstrating the best recovery in radial growth but lower resistance and resilience relative to the other two varieties.

Tree Dimensions and Allometry
On average, 21-year-old QURO specimens from the analyzed stand were 9.35 m tall and had crown widths of 8.05 m and DBH of 23.71 cm; however, significant variations were observed among QURO varieties. In total, 24, 25, and 27 half-sib families of early, typical, and late QURO varieties were included in this study. Based on two-way ANOVA, half-sib families have a lower influence on DBH, TH, and CW deviations among QURO varieties ( Table 2). On the other hand, varieties as well as varieties x half-sib statistically significantly separated all QUROs' DBH, TH, and CW. Otherwise, half-sib samples did not exhibit statistically significant deviations. Based on F-test results, varieties are more strongly differentiated in terms of DBH (F 261.40) than TH (F 2.48) or CW (25.89). However, half-sib has the strongest effect on DBH (F 26.24), followed by TH (F 2.91), and finally CW (F 1.98). Variety x half-sib interaction exerts the strongest influence on CW (F 51.60) followed by TH (F 48.47). The greatest trunk diameter at breast height was measured for T-QURO (25.63 cm), followed by L-QURO (24.20 cm), and E-QURO (21.35 cm), with each variety forming a statistically significant homogenous group according to Tukey's HSD rank test results ( Figure 6). The same pattern was noted for CW, but not for tree height. E-QURO (9.46 m) and L-QURO (9.38 m) performed better than T-QURO (9.19 m) in terms of height increment. Analyzed varieties vary by 1.50 m in CW in favour of T-QURO, which had wider average crown radii (8.73 m) than did E-QURO (7.23 m).  The greatest trunk diameter at breast height was measured for T-QURO (25.63 cm), followed by L-QURO (24.20 cm), and E-QURO (21.35 cm), with each variety forming a statistically significant homogenous group according to Tukey's HSD rank test results ( Figure 6). The same pattern was noted for CW, but not for tree height. E-QURO (9.46 m) and L-QURO (9.38 m) performed better than T-QURO (9.19 m) in terms of height increment. Analyzed varieties vary by 1.50 m in CW in favour of T-QURO, which had wider average crown radii (8.73 m) than did E-QURO (7.23 m). Tree dimension proportions (allometric relations) for each QURO variety are shown in Figure 7, while the OLS linear regression findings are summarized in Table 3. We analyzed three allometric pairs (TH/DBH, CW/DBH, and CW/TH) to explain tree dimension proportions. In contrast to measured tree dimensions, allometric proportions vary less among QURO varieties though some important differences were observed. As shown in Figure 6 and confirmed via regression model outputs (mostly Adj. r 2 ), the model provided Tree dimension proportions (allometric relations) for each QURO variety are shown in Figure 7, while the OLS linear regression findings are summarized in Table 3. We analyzed three allometric pairs (TH/DBH, CW/DBH, and CW/TH) to explain tree dimension proportions. In contrast to measured tree dimensions, allometric proportions vary less among QURO varieties though some important differences were observed. As shown in Figure 6 and confirmed via regression model outputs (mostly Adj. r 2 ), the model provided a better fit for CW/DBH relative to TH/DBH and CW/TH. All analyzed allometric relationships were highly statistically significant (p < 0.00 Fisher test (F-test) results ranged from 80.4 (CW/TH; T-QURO) to 629.1 (CW/DBH; QURO). Then, adjusted coefficient of determination (Adj. r 2 ) values for the regressi model were similar across the varieties and ranged from weak (CW/TH 0.22; T-QURO) strong (CW/DBH 0.67; E-QURO). The model also provided a better fit to the CW and DB data relative to the TH/DBH and CW/TH pairs. No patterns were observed among QUR varieties with regard to Adj. r 2 strength. Although T-QURO's TH/DBH and CW/DBH lometric pairs showed weaker fitting, this pattern was not observed in the CW/TH allom tric relationship. However, differences between L-QURO and E-QURO lowered t strength of their relationships with the Adj. r 2 of T-QURO.   All analyzed allometric relationships were highly statistically significant (p < 0.001). Fisher test (F-test) results ranged from 80.4 (CW/TH; T-QURO) to 629.1 (CW/DBH; E-QURO). Then, adjusted coefficient of determination (Adj. r 2 ) values for the regression model were similar across the varieties and ranged from weak (CW/TH 0.22; T-QURO) to strong (CW/DBH 0.67; E-QURO). The model also provided a better fit to the CW and DBH data relative to the TH/DBH and CW/TH pairs. No patterns were observed among QURO varieties with regard to Adj. r 2 strength. Although T-QURO's TH/DBH and CW/DBH allometric pairs showed weaker fitting, this pattern was not observed in the CW/TH allometric relationship. However, differences between L-QURO and E-QURO lowered the strength of their relationships with the Adj. r 2 of T-QURO.
Regression line slope (parameter b) defines the allometry of the analyzed tree dimension parameters. Among QURO varieties, stronger variations were noted with respect to CW/TH (0.43-0.50) and, to a lesser extent, to TH/DBH (0. 13-0.19). The b coefficient was nearly identical for CW/DBH among all three analyzed QURO varieties, ranging from 0.19 to 0.22. T-QURO was shorter relative to DBH, while L-QURO had a wider crown relative to tree height. All three QURO varieties had an identical relationship between crown width and DBH.

Intra-Variety Deviations in Radial Growth Dynamics
The study findings reveal variations among QURO phenological varieties as well as across time. Differences in radial growth can be attributed to high QURO genetic diversity [21] across the species' distribution range [18]. In the provenance trial conducted by Bert and colleagues, spatial arrangement (i.e., provenance-separated QURO populations) was shown to contribute to significant deviations in radial growth and resistance to extreme climate conditions [46]. These authors observed different radial growth in two provenance trials, but also the same pattern across provenances. The southern and xeromorphic QURO provenances are better adapted to the changes expected in climate conditions over the next decades.
The sensitivity of QURO to stand conditions has been demonstrated in previous studies [9,12,13,47], as well as in the current study. Nonetheless, clear spatial patterns among QURO populations across Europe cannot easily be defined based on wood properties alone. Above-and belowground environmental factors and the genetic potential for biomass production also influence QURO radial growth. An earlier study noted that soil water availability and above-ground water sources (e.g., precipitation, soil water, ice, etc.) equally affect radial growth [12]. The same authors observed that soil properties significantly influence radial growth and the tree's sensitivity to drought. Radial growth is also influenced by site-specific geomorphological elements. For instance, rivers are significant modifiers of QURO radial growth performance, as flooded QURO populations exhibit greater radial increments as compared to non-flooded forests [9]. Similarly, micro-site stand conditions have been shown to shape the radial growth of Scots pine in temperate climate zones [48].
If many phenological QURO varieties were selected, the relationship between phenology and radial growth would not be fully understood. On the European level, studies showed that sessile oak vegetation dynamics and their adaptive mechanism were confirmed [49]. A study by Puchałka et al. [25] confirmed a link at the anatomical level between QURO phenology and tree-ring formation dynamics. Following the cambial phenology of 12 QUROs, the study found higher cambial activity up to the middle of April, and EW development until the end of May [25].

Drought Sensitivity of Different QURO Varieties during Juvenile and Mature Tree Growth Phases
Tree performance varies throughout the life cycle [50]. In particular, significantly different growth dynamics are observed in juvenile and mature trees, thus affecting wood properties [37]. Younger trees grow faster, exhibit more variable radial growth, and have smaller vessels with thinner walls than more mature trees. As a result, younger trees have a lower wood density [37,51] and a different physiological response than older trees [52]. The length of the transition period from the juvenile to the mature growth phase can vary considerably (by up to few decades) depending on the species-specific characteristics and stand conditions [51,52]. In the climate conditions at the study site, QURO transitions from juvenile to mature tree grow phase around the 10th year [19].
Since our common garden experiment was established, two extreme drought events were recorded. The first extreme drought occurred in 2012, which coincided with QURO's juvenile growth phase. The second was recorded in 2017, after the analyzed trees had reached the mature growth phase (according to the radial increment dynamics). All QURO specimens exhibited better performance, resistance, resilience, and recovery potential during the juvenile growth phase than during the mature phase. These findings are in line with the results reported by McGowran et al. [52], who confirmed that trees exhibit superior morpho-physiological performance during the juvenile growth phase as compared to the mature growth phase. Because the analyzed QUROs exhibit good resistance, resilience, and recovery characteristics, these varieties should be considered when establishing "climate-smart" forests [53].
However, the responses of the analyzed phenological varieties to drought events varied. E-QURO had a greater resistance and resilience but slightly lower recovery potential, which could be attributed to vegetation dynamics and climate synchronization. Under optimal climate conditions, early QURO finishes its intensive growth phase near the end of April [25]. Later varieties can therefore experience optimal conditions that are up to one month shorter. We also found that trees that were growing in stressed conditions showed better recovery potential than trees growing under optimal conditions. This is because they have already been adapting to stress-a situation known as stress priming [54,55]and have therefore already adapted the hydraulic anatomy of their wood tissue to more stressful conditions [56].
Because the flowering phenophases of early and late QURO varieties do not coincide [6], their crossing in natural populations is not possible. These two varieties therefore show significant genetic deviations from one another. Vessel properties, which define plant hydraulic capacity and shape species-specific responses to drought events [57], provide a high level of heredity [7]. Hence, L-QURO's better resistance and resilience should be considered as a genotypic characteristic rather than a response to the surrounding environment.

QURO Inter-Variety Deviations in Tree Dimensions and Allometry
Tree dimensions and their proportions (i.e., allometry) are influenced by the surrounding environment [58,59], as well as by a tree's fitness [60] and genetic potential [61]. The differences in tree dimensions among QURO varieties observed in this study are in line with the high variation in the biomass production of QURO across the different provenances [46]. Trunk cambial activity is linked to tree phenology [25].
Tree allometry is more than tree dimension shaped via genetics. If provenance tests showed significant and stronger provenance than site influence of QURO growth [46], surrounded effects could not be ignored in studies. Dai et al. [60] showed that stressors, such as drought, significantly affect tree allometry. Moreover, the same authors emphasize that mentioned relations are not simple and are caused by a lot of surrounded influences/stressors, the same as a tree's survival potential in stress periods. Light availability [58], stand mixture, and density [59], for instance, also significantly shape tree allometry. As shown by this study and other allometric studies, trunk diameter is recognized as a more reliable variable than crown dimensions [62].
Previous studies have confirmed inter-site variations in QURO tree allometry [63] and radial growth [46], which are closely linked [64]. Because our study focuses on even-aged trees from a unique common garden experiment, the influence of external factors was minimized. Thus, any observed differences can be attributed to genetic differences.

Perspective of Implications Phenological QURO Varieties for Climate-Resistant Forestry
The ability of oak forests to adapt to climate change is key to their survival [65], especially for populations growing in lowlands and arid regions that are extremely sensitive to upcoming climate conditions. Surrounding environment impacts on tree performances are complex and not completely clear. The phenological diversity of oaks is a valuable resource for the cultivation of more climate-resistant QURO specimens/populations.
According to the established RCP scenarios, climate changes in the near future will cause significant changes in forest dynamics [50]. Trees must be able to adapt to the new conditions, but this can be difficult to do for long-lived tree species like QURO. If endangered QURO forests are to survive, forest managers must focus on adaptationoriented management strategies. Our study shows that phenological variety exists with regard to the resistance and recovery potential of QURO to drought events. Observed deviations labelled this approach correctly based on phenology dynamics in the context of overcoming problems with intensive climate changes.

Conclusions
Coetaneous (21 years) QURO phenological varieties differ from one another with regard to radial growth, biometry, and allometry performances. Their response to drought events varies across both variety and tree growth stage (juvenile and mature). Likewise, oak varieties differ with regard to growth dynamics. Late QURO shows faster radial growth in the first years, whereas typical and early QURO varieties have similar radial growth chronologies. Although early QURO experiences a longer vegetation period (by one month) under optimal conditions, differences between early and latewood are smaller than differences noted between whole tree-ring measurements.
The analyzed QURO varieties exhibit different resistance, resilience, and recovery responses to drought. These responses also differ across growth phases. The early QURO variety shows better resistance, but late and typical varieties exhibit better recovery after drought events. In both juvenile and mature tree growth stages, all QURO varieties show sensitivity. All QURO varieties exhibited better performance during the juvenile growth phase.
In conclusion, the phenotypical differences observed in this common garden experiment should be of great relevance to both forestry practitioners and scientists. The observed differences show that the analyzed gene pool has the potential to produce drought-resistant but highly productive QURO specimens. An understanding of the range of phenotypical variety available may play an important role in helping forest managers shift from "classic" to "climate smart" forest management.

Data Availability Statement:
The data presented in this study are available from the corresponding author upon request.