Next Article in Journal
Banana Classification Using Sanger Sequencing of the Ribosomal DNA Internal Transcribed Spacer (ITS) Region
Next Article in Special Issue
Integrated Effects of Soil Moisture on Wheat Hydraulic Properties and Stomatal Regulation
Previous Article in Journal
Enhancing Withanolide Production in the Withania Species: Advances in In Vitro Culture and Synthetic Biology Approaches
Previous Article in Special Issue
Recent Advances in Studies of Genomic DNA Methylation and Its Involvement in Regulating Drought Stress Response in Crops
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Plants
Volume 13
Issue 16
10.3390/plants13162172
plants-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Tree Resilience Indices of Norway Spruce Provenances Tested in Long-Term Common Garden Experiments in the Romanian Carpathians

by
Alin Madalin Alexandru
1,2,*,
Georgeta Mihai
1,
Emanuel Stoica
1,2 and
Alexandru Lucian Curtu
2,*
1
Department of Forest Genetics and Tree Breeding, “Marin Dracea” National Institute for Research and Development in Forestry, 077190 Voluntari, Romania
2
Faculty of Silviculture and Forest Engineering, Transilvania University of Brasov, 500123 Brasov, Romania
*
Authors to whom correspondence should be addressed.
Plants 2024, 13(16), 2172; https://doi.org/10.3390/plants13162172
Submission received: 13 June 2024 / Revised: 25 July 2024 / Accepted: 4 August 2024 / Published: 6 August 2024
(This article belongs to the Special Issue Drought Responses and Adaptation Mechanisms in Plants)
Browse Figures
Versions Notes

Abstract

:
Provenance trials provide a valuable opportunity to evaluate the impact of extreme events on growth and wood properties. In this study, we have evaluated 81 Norway spruce provenances, tested in three provenance trials established in the Romanian Carpathians in 1972. The response to drought of the Norway spruce provenances has been examined using the following tree resilience indices: resistance, recovery, resilience, and relative resilience. The relationship between climate and growth, the correlations between wood traits, and the coordinates of the origin and tree resilience indices were also analysed. In each provenance trial, there were significant differences between provenances and years regarding wood widths and latewood percentage (LWP). Regarding drought extreme events, the years when they occurred in all three provenance trials were 2000 and 2003. Significant differences between provenances for at least one tree resilience index have been found in all provenance trials, for the year 2000. By using subperiods of 25 years, changes in the relationship between climate and growth have been observed. Several provenances with high radial growth and good resistance and/or recovery have been identified. Provenances that performed better in common garden experiments could be used in assisted migration, even in the proximity of the current natural range.

1. Introduction

The 1.5 °C limit of global warming might be reached in the first half of the 2030s [1]. As a result, severe and extreme weather events will become more frequent. Among them, drought and heat stress are considered to have the greatest adverse impact on forest ecosystems. The regions most vulnerable to increasing drought risk will be Southern Europe, with the Mediterranean region as a hotspot, and Southeastern Europe [2,3,4].
Although it is considered that forests have a large potential for climate change mitigation, droughts and other extreme weather events can threaten their growth, productivity, genetic diversity, and distribution [5,6,7,8,9,10]. In Europe, drought events have caused significant impacts on forests in the last decades, as the century passes [11]. Between 1987 and 2016, drought was the cause of excess forest mortality for almost half a million ha [12].
In Romania, from 1901, every decade has had from one to four years with extreme droughts/rain events. The number of identified droughts has been increasing since 1981 [13]. Annual precipitation is expected to decrease, particularly in the southeastern part of the country [14]. Analysing the predictions made for Romania for the short term, 2021–2050, and the long-term, 2071–2100 versus 1971–2000, the eastern and northeastern parts will have large changes in winter temperature, and the eastern and southern parts will have larger changes for the other seasons, more than 2 or 4 °C, for RCP4.5 and RCP8.5, respectively. Regarding precipitation, it is expected to increase, more so in winter, in the central part of the country. At the same time, for the summer, a decline is seen in most cases, the most affected being the southeastern and southwestern areas [15].
With a distribution area of around 30 million ha, Norway spruce (Picea abies) is an economically important coniferous species in Europe. The vulnerability of Norway spruce to drought has been expressed in many studies [16,17]. It is highly sensitive to rising temperatures and water shortfall [18] and to declining water availability, and it has a weak adaptive capacity at lower elevations [19,20]. Drought causes a decrease in its vitality [21], which leads to a lower resistance against bark beetles, especially Ips typographus [22]. Therefore, the risk of mortality for Norway spruce and other native coniferous species might be greater with increased temperatures and decreased precipitation amounts, especially for populations that are established outside the natural range [23].
In Europe, according to RCP 8.5 projections, the frequency of droughts is projected to increase over the whole of Europe [11]. The Norway spruce might retreat in the Alps, the Carpathians, and Scandinavia above 60° latitude, disappearing completely from the Central European lowlands [24]. In Romania, the suitable area lost by Norway spruce until 2100 is estimated to be 8%, while the gained area is around 2% [25]. Mihai et al. [26] have revealed that the climate envelopes for the main forest species in Romania have already shifted to another ecosystem climate. Norway spruce will expand to higher altitudes but will significantly decrease in frequency and lose its habitat, particularly in the Eastern Carpathians.
Although the climate will change more rapidly than some populations can adapt, forest species hold different adaptive capacities to mitigate the impacts of climate change. Resilience to drought is a trait with environmental and genetic components, and different populations of the same species might respond differently to the same climate [27]. Recent studies have shown that there is a genetic variation in the sensitivity/tolerance of forest species against drought events, both within and among populations [28,29,30,31]. Using genome analyses, significant associations were found between SNPs (single-nucleotide polymorphisms) and traits related to the drought adaptation of Norway spruce [32].
Therefore, an approach to ensure the genetic adaptation of forest species to the negative impacts of drought is to select the most resilient provenances [28,33]. Using seed sources that will be climatically adapted during the rotation cycle will increase the resilience of forest ecosystems and will maintain the services provided by them. Selection of the most suitable provenances, together with assisted migration, can be considered the most effective strategy for the management and conservation of forest genetic resources in the face of environmental change.
Field tests, such as provenance trials, provide a valuable opportunity to evaluate the impact of extreme events on growth and wood properties. By testing different provenances in homogenous conditions, the most valuable and adapted provenances can be selected. Therefore, provenance trials can provide important data concerning intraspecific adaptive capacity and the selection of suitable populations for reforestation programs. Several studies have analysed the intraspecific variation in climate growth of Norway spruce in provenance trials [34,35,36,37,38,39]. However, the trials from the previous studies were not comprised of so many provenances, and those provenances were studied at younger ages, than the trials from our study, which are comprised of 81 provenances that cover well the distribution area of the species in Europe and are at half the rotation age of this species in Romania.
To quantify trees’ resilience to climate change, tree-ring width measurements are frequently used. Ring width is linked with conifers’ wood structure; hence, changes in radial growth patterns will alter forest productivity, wood properties, and timber quality [40]. In the first half of the growing season, large and thin-walled earlywood tracheids are formed that are more efficient in water transport but are also more vulnerable to cavitation. In the second half of the growing season, narrow and thick-walled latewood tracheids are formed, which provide mechanical support [41].
This study aimed to assess the growth response and adaptive capacity of Norway spruce provenances originating from thirteen European countries to extreme drought events which have occurred in this region in the last 49 years. The objectives were to (i) determine the genetic variation of wood traits among the Norway spruce provenances, (ii) evaluate the response to drought of the provenances, (iii) determine the relationship between growth and the climate of the trial sites, (iv) compute correlations between wood traits, tree resilience indices, and the geographical coordinates of the provenances’ origin, and (v) provide information for the implementation of sustainable forest management and forest genetic resource conservation.
This study will improve knowledge regarding the adaptive capacity of Norway spruce populations of different geographical origins to the cumulative effects of drought and heat stress associated with climate change.

2. Materials and Methods

2.1. The Provenance Trials

The three provenance trials from Romania, Dorna Candrenilor, Zarnesti, and Turda, analysed in this study, are in the Eastern, Southern, and Western Romanian Carpathians, respectively. The Zarnesti and Dorna Candrenilor provenance trials are in the mixed beech and coniferous species zone, while Turda is in the Norway spruce vegetation zone. The climate conditions are continental, with Scandinavian–Baltic influences in the north (Dorna Candrenilor) and temperate–continental with oceanic influence in the west (Turda) and central parts (Zarnesti) [42]. The three Norway spruce provenance trials were established in 1972 [43]. Eighty-one provenances were tested in these trials, ten from Romania and the rest from twelve other European countries (Figure 1). The geographic coordinates and details about the tested provenances are presented in Supplementary Table S1.
Provenances were planted in a randomised complete block design, each plot with 16 (4 × 4) individuals per provenance at a 2 × 2 m spacing and three blocks. No thinning or artificial pruning was performed in these provenance trials.
The wood cores were extracted in the fall of 2020. Given the advanced age of the provenance trials, where possible, four trees were randomly selected for each provenance in each repetition. With the Haglöf increment borer (Haglöf Sweden AB, Långsele, Sweden), one core/tree was extracted at breast height, on the slope line, to avoid compression and tension wood. Only one core from each tree was extracted to avoid tree damage.
The cores were then dried and progressively sanded to obtain a smooth surface. For the image to be of good quality (1200 dpi), an Epson Expression 12000XL scanner (Seiko Epson Corporation, Nagano, Japan) was used. Ring width (RW), earlywood (EW), and latewood (LW) were measured with CooRecorder software version 7.4 [44]. The latewood percentage (LWP) was calculated as the ratio between LW and RW. For each trial, the series were checked and cross-dated using the detrendeR version 1.0.5 [45] and dplR version 1.7.3 [46] packages in the R free software environment version 4.0.0 [47]. Dendrochronological series with intercorrelation values below 0.328 were excluded from the analysis. The intercorrelation is the correlation coefficient between different cores. The final number of wood cores was 2709: 947 for Dorna Candrenilor, 906 for Zarnesti, and 856 for the Turda trial.
The spline method has been used to detrend the ring series, where the frequency response is 0.50 at a wavelength of 0.67 × “series length in years”. By detrending, the tree’s natural biological growth trend is estimated and removed, to produce the dimensionless ring width index (RWI).
Climatic data for the period 1972–2020 were obtained using the Climate Downscaling Tool (ClimateDT) [48], at a resolution of 1 km. For the period 1961–1990, it uses a high-resolution grid as the baseline, while for the historical period, it uses monthly anomalies derived from Climatic Research Unit gridded Time Series (CRU-TS) [49]. For future scenarios, UKCP18 surfaces [50] are being used. Except for the minimum and maximum monthly temperatures, all other climatic variables were used to evaluate the relation between climate and growth.
To determine the years with meteorological droughts, the Standardised Precipitation Evaporation Index (SPEI) [51] was calculated based on precipitation and the potential evapotranspiration (PET) for the 1972–2020 period, using the SPEI R package version 1.7 [52]. The Thornthwaite equation [53] was used for calculating PET, which can be computed only with temperature data. Using the SPEI at a time scale of 3 months (SPEI-3) (which uses the moving average of precipitation for three months, e.g., the value of SPEI-3 for March is calculated based on the precipitation and PET from January, February, and March), the drought years that recorded values below −1 were classified as follows: from −1 to −1.49, moderate drought; from −1.5 to −1.99, severe drought; less than −2, extreme drought.
The response to drought events was evaluated using the tree resilience indices proposed by Lloret et al.: resistance (Res), recovery (Rec), resilience (Rsl), and relative resilience (relRsl) [54]. Resistance indicates the reversal of the decrease in radial growth during the drought. It is the ratio of RW during the drought (Dr) and before the drought (preDr). Values above 1 indicate high tolerance, while values below 1 mean low tolerance. Recovery indicates the capacity to revitalise after a drought. It is the ratio of RW after drought (postDR) and RW during drought. Resilience indicates the capability of reaching a pre-drought RW after a drought event. It is the ratio of RW after drought and RW before the drought. Values below 1 indicate long-term growth reductions. Relative resilience was calculated using the formula relRsl = (postDr − Dr)/preDr. It is the resilience accounting for the impact that the drought year had. The before and after RW (preDR and post DR) were calculated as average values for a period of three years before and after a year with a drought event.
To verify that the years with extreme drought events that we identified were actually years when the trees were stressed, pointer years were calculated in the dplR R package version 1.7.3 [46], based on the Becker algorithm [55]. The term ‘pointer years’ denote years in which the greater part of the trees show very narrow (negative pointer year) or wide ring widths (positive pointer year) [56,57]. A threshold of 60% of influenced trees was used to identify the pointer years.

2.2. Statistical Analysis

An analysis of the genetic variation in wood traits and tree resilience indices was performed for each trial site, considering that droughts and post-drought effects vary between ecoregions [58]. Using the lmerTest R package version 3.1-3 [59], a mixed linear model was applied and tested, where provenance and year x provenance were considered random effects and the block as a fixed effect.
Xijk = μ + Pj + Bk + Yl + Pj × Yl + ejklm,
where μ is the general mean, Pj is the effect of the jth provenance, Bk is the effect of the kth block, Yl is the effect of the lth year, Pj × Yl is the provenance-by-year interaction, and ejklm is the error term associated with jklm trees [60].
Pearson coefficients were used to assess the correlations between the geographical coordinates (latitude, longitude, altitude) of the provenances’ origin and the wood traits and mean tree resilience indices, to examine the degree to which trait variation is influenced by local adaptation to climatic conditions from the origin location. The coordinates were not transformed prior to the correlation analyses.
The response and correlation function analysis from treeclim R package version 2.0.6.0 [61] was used to assess the non-stationarity of climate–growth correlations in each provenance trial. Non-stationarity was defined as a modification in the relationship’s magnitude or direction between variables over time [62]. The climatic variables used for this analysis were monthly values of temperature and precipitation from each year, from March to September (growing season), for the 1981–2020 interval. Subperiods of 25 years, each following one starting a year later, were analysed, resulting in 16 subperiods.

3. Results

3.1. Drought Years Identification

In each study site, large variations in annual mean temperature and annual precipitation amounts were recorded (Figure 2). The year 2000 had the lowest amount of precipitation for the analysed period in all three sites, ranging from 511 mm at the Turda trial, to 646 mm at the Zarnesti trial. The highest MAT was recorded for the year 2019, in all three trials, with values ranging from 6.6 °C at the Turda trial to 7.3 °C at the Dorna Candrenilor and Zarnesti trials, respectively.
Large variations in temperature and precipitation can be seen for the analysed period (1972–2020). For temperature, both the annual mean and the one from the warmest quarter, there is a clear trend of increasing values. Regarding precipitation, although the annual amount has decreased in the last decade, the most remarkable fact is the difference between 2010 and the 2011 APs. The year 2010 was the one with the highest amount of annual precipitation for the Dorna Candrenilor and Turda trials and the second highest for the Zarnesti one. In contrast, the year 2011, even though it was not classified as extreme according to SPEI-3 in our study, has one of the lowest APs from the analysed period: the second, third, and fourth ones for the Turda, Dorna Candrenilor, and Zarnesti trials, respectively.
The year 2000 had the lowest AP for the analysed period, the values dropping to 511, 553, and 646 mm in the Turda, Dorna Candrenilor, and Zarnesti trials, respectively. The optimum climate for Norway spruce in Romania is represented by an MAT between +4 and +7 °C and an AP between 800 and 1200 mm [63]. At Dorna Candrenilor, the AP dropped in seven years even below the limit value of 600 mm (1986, 1990, 1994, 2000, 2003, 2011, 2015) and only 16 out of the 49 years had optimum AP values (above 800 mm). At the Turda trial, regarding AP, four years had values below the limit value. At the Zarnesti trial, ten years showed APs below the optimum value, but none were below the limit.
For the period 1972–2020, moderate, severe, and extreme drought years were identified based on the SPEI-3 values (Figure 3). In the Dorna Candrenilor and Turda trials, 23 years with severe and extreme droughts have been identified for the analysed period, while for Zarnesti, the number was 24 (Table 1).
The number of years with extreme drought varied among the sites and ranged between three at the Turda trial to five at the Dorna Candrenilor trial. The years when extreme drought occurred in all three provenance trials were 2000 and 2003, as confirmed by the pointer years analysis (Figure 4).

3.2. Norway Spruce Series

For the analysed period, the mean RW ranged from 2.55 mm, at the Zarnesti trial, to 2.81 mm, at the Turda trial.
The mean series intercorrelation (rbar) was calculated for each provenance trial. It is the mean correlation coefficient among the tree-ring series, and it ranged from 0.876 to 0.908 (Table 2).

3.3. The Effect of Provenance and Year on Radial Width and Latewood Percentage

In each provenance trial, there were significant differences between provenances and years regarding wood width and LWP (Table 3). The block effect was significant for all traits, except RW at the Turda trial. Regarding year x provenance interaction, it was significant for RW and EW at the Zarnesti trial, and for RW, EW, and LWP at the Dorna Candrenilor trial.
The wood trait variation over the years for each provenance trial can be seen in Figure 5.

3.4. Genetic Variation in Drought Response

Analysing the provenances’ different responses to the drought from the year 2000 (Figure 6), significant differences were found between the provenances for all tree resilience indices in the Zarnesti trial. At the Turda trial, the differences were significant for all indices, except resilience, while at Dorna Candrenilor, resilience was the only factor with significant differences between provenances (Table 4).
The best results regarding resistance in the year 2000 were obtained at the Turda and Dorna Candrenilor provenance trials, with a mean resistance of 0.89 and 0.88, respectively, higher by 23.61 and 22.22%, respectively, than the mean resistance from the Zarnesti provenance trial (0.72).
The mean value for recovery of the year 2000 was above one in all three provenance trials, the highest being achieved in the Zarnesti trial, 1.15 ± 0.12 (SD).
The provenance trials Dorna Candrenilor and Turda had the greatest mean resilience for the year 2000, 0.90 and 0.89, respectively.
The mean of the relative resilience was between 0.0, in the Turda trial, and 0.07, in the Zarnesti trial.
For the drought event from 2003, significant differences between provenances for all tree resilience indices were found only in the Zarnesti trial. No significant differences between provenances were found for any tree resilience indices in the other two sites (Table 4).
The drought event from 2000 impacted the three provenance trials differently. Its highest impact was at the Zarnesti trial, where no Norway spruce provenance had a mean resistance above 1; in other words, all provenances were affected (Figure 6).
The highest amplitude for the mean value of the provenances’ resistance was at the Zarnesti trial, ranging from 0.55, provenance 26-Winterthur, to 0.97, provenance 83-Bramarv. For this trial, the site mean was 0.72, the lowest of the three trial sites. For the Turda trial, it ranged from 0.67, provenance 26-Winterthur to 1.07, provenance 38-Val Di Fiemme. And for the Dorna Candrenilor trial, the mean provenance resistance ranged from 0.75, provenance 34-Borki, to 1.02, provenance 2-Branstad.
Only five provenances had values below 1 regarding recovery at the Zarnesti trial: 88-Pualanka, 85-Heinola, 82-Sund, 101-Valke Karlovice, and 83-Bramarv. At the Dorna Candrenilor trial, the number was 20 provenances, and at the Turda trial, it was 28.
No provenances had a mean resilience above one at the Zarnesti trial, with values ranging from 0.65, provenance 5-Seljord, to 0.95, provenance 37-Latemar. The Northern provenances 2-Branstad, 3-Sandar, 92-Padasjoki, and 55-Munkahus were the ones that had a resilience above 1 at the Dorna Candrenilor trial, for the year 2000; 37-Latemar, 40-Wietersdf and 20-Le Brassus had a mean resilience of 1; the rest had mean values ranging from 0.76, provenance 60-Keletbukki Allami, to 0.99, provenance 20-Le Brassus. At the Turda trial, values above 1 for resilience were shown by provenances 83-Bramarv (1.06) and 5-Seljord (1.05). Provenances 93-Urjala and 85-Heinola had a resilience of 1, and the rest had values ranging from 0.75, provenance 1-Senum, to 0.99, provenance 45-Hollenburg and 67-Frasin.

3.5. Phenotypic Correlations

Significant and positive correlations between elevation and wood traits were found only at the Dorna Candrenilor trial, except LWP (Table 5). Also, at this trial, the longitude was negatively and significantly correlated with LW and LWP. At the Turda trial, longitude was negatively correlated with EW and positively correlated with LWP.
Regarding latitude, significant and negative correlations were found at the Dorna Candrenilor trial for all traits except LWP. At the Zarnesti trial, latitude was negatively correlated with LWP only. No significant correlations were found between wood traits and latitude at the Turda trial.
The correlation between EW and LW ranged from 0.309 at Dorna Candrenilor to 0.753 at Turda. Correlations between RW and EW were between 0.922 in the Dorna Candrenilor trial and 0.983 in the Turda trial. Negative correlations were found between RW and LWP in the Dorna Candrenilor and Zarnesti trials.
Regarding the correlations between RW and LW, the highest value was obtained at the Turda trial, 0.861.
Regarding the mean of tree resilience indices from 2000 and 2003, at the Dorna Candrenilor trial, elevation was positively correlated with resistance (r = 0.128) and negatively correlated with recovery (r = −0.132) (Table 6). Latitude was positively correlated with recovery.
At the Zarnesti trial, however, no significant correlations were found between elevation and tree resilience indices. Latitude was negatively correlated with all tree resilience indices. The same for longitude, except the correlation with resistance was not significant (r = 0.120, p = 0.061).
At the Turda trial, no significant correlations were found regarding elevation and latitude. Longitude was positively correlated with resistance (r = 0.167), and negatively correlated with recovery (−0.151) and relative resilience (r = −0.152).
Regarding the correlations between the mean of the 2000 and 2003 tree resilience indices and mean wood traits (Table 7), resistance was negatively and significantly correlated with EW and RW at the Dorna Candrenilor trial, and with EW, LW, and RW at the Turda trial. In the same two trials, the correlations were positive with LWP.
There was a positive correlation between recovery and EW and a negative one with LWP at the Dorna Candrenilor trial. At the Zarnesti trial, there was a low correlation between recovery and LW (r = 137*). In contrast, recovery was highly correlated with EW, LW, and RW at the Turda trial.

3.6. The Response and Correlation Function Analysis

Analysing the response and correlation function, precipitation in July always had a positive effect on ring width index at the Dorna Candrenilor trial (Figure 7), but it became significant after 1991. Contrarily, the effect of the precipitation from April slowly decreased, having even negative values for the last period, 1996–2020. The temperature of March also had a positive effect, more so after 1990. The temperature of May, although positive in the first periods, from 1981–2010, had only a negative impact later on, and it seems it is only getting stronger.
At the Zarnesti trial, precipitation in June was significant for the first three periods of 1981–2007 (Figure 8); the value of the correlation coefficient, while still positive, was decreasing. Precipitation from September and the temperature in May, while not significant, had a negative impact but became stronger in the last periods (three and four, respectively).
At the Turda trial, precipitation in June and July always had a positive effect on ring growth, but more so for the last three and four periods, respectively (Figure 9). The negative effect of March temperature for the first period was analysed, 1981–2005; it became positive starting with the 1983–2007 period, even being significant in the 1993–2017 period. The temperature of May, which was significant for the 1987–2011 and 1988–2012 periods, had a negative effect on ring growth for all the analysed periods, except the last three (1994–2020), when the correlation coefficients had values close to 0.

4. Discussion

The radial width and early- and latewood, as well as latewood percentage and the response to drought, of 81 Norway spruce provenances have been analysed in three long-term provenance trials located in three different geographic regions of the Romanian Carpathian Mountains.
The plant response to drought stress is very intricate, and various strategies can be adopted to endure drought events. Baldi and La Porta classify these strategies into three main groups: drought avoidance, drought resistance, and drought resilience [64]. The drought avoidance strategy (e.g., deep roots or stomatal adjustments) represents a tree’s capacity to minimise the amount of water lost and maximise water uptake [65,66]. Drought resistance was defined as the capacity of a tree to endure stress, and drought resilience was defined as how swiftly a tree can resume normal growth, once the stress stopped [54,67]. The traits and processes that are involved in the drought avoidance, resistance, and resilience of conifers have been reviewed in other studies [68].
The negative effect of rising temperatures on conifer growth was also observed in the study of Gómez-Aparicio et al. [69], which can be connected to a higher atmospheric water demand, which causes greater drought stress. Positive correlations between Norway spruce RW and the climate of July and August of the current year but also October of the previous year, and a negative effect of the temperature of September from the previous year, have been reported in the Czech Republic [70].
The precipitation amounts from June and July also became important for pine growth in Lithuania [71], a change from earlier studies [72]. In the study of Lévesque [16], soil water availability from December to July was also important. Roibu [73] found that the correlations with temperature and precipitation also changed over time for ash and oak from the Republic of Moldova, using moving window correlations of 25 years.
The differences between provenances regarding wood traits and drought response were very significant in each trial. This can be explained by the broad geographic amplitude of the tested provenances, which comprised provenances of the Alpine range, the Bohemian Massive, the Carpathian Mountains, the Rila and Rhodope Mountains, the Vosges and Jura Mountains, and provenances from Northern Europe (Sweden, Norway, and Finland), as well as provenances from NE Germany and NE Poland. The high adaptive variation of Norway spruce has been observed in Central and Southeastern Europe [32].
Recovery was positively correlated with latitude in Dorna Candrenilor but negatively correlated with it in the Zarnesti trial. In the Dorna Candrenilor trial, which is a forest site of superior quality, the Southern provenances had a greater resistance and thus a lower recovery. On the other hand, in the Zarnesti trial, which is a forest site of medium–inferior quality, the Southern provenances had a lower resistance and, thus, a greater recovery (see Figure 6).
The provenances with a high radial growth and good resistance and/or recovery in all three trials are the following: 53-Neustift, 49-Redl-Zipf-Fuchsberg, and 40-Wietersdf, from the Eastern Alps; 19-Kerns and 25-Wassen, from the Central Alps; 31-Bremenhagen, from NE Germany; 99-Zelesna Ruda, from the Bohemian Forest Mountains (Czech Republic); and 67-Frasin, 70-Coşna, 71-Moldovița, and 75-Broşteni, from the Eastern Romanian Carpathians.
Provenance 53-Neustift, from the Eastern Alps, had a great mean RW in all three trials. The resistance for the year 2000 had values between 0.65, at the Zarnesti trial, to 0.86, at the Turda trial. Regarding recovery, it ranged between 0.97 and 1.35 at the Dorna Candrenilor and Zarnesti trials, respectively. And the resilience was between 0.80 and 0.89 at the Dorna Candrenilor and Turda trials, respectively.
Another valuable provenance is 25-Wassen, from the Central Alps; with a mean RW of 3.11 and 3.52 mm at the Dorna Candrenilor and Turda trial, respectively, it was the first and third ranking in these provenance trials. In the Zarnesti trial, it occupied rank 32, with a mean RW of 2.58 mm. It was not very resistant, with values for the resistance index ranging from 0.66, at the Zarnesti trial, to 0.85, at the Turda trial. The resilience had values between 0.86 and 0.93 at the Zarnesti and Turda trial, respectively. This provenance’s recovery had great values: from 1.12 to 1.38 at the Dorna Candrenilor and Zarnesti trials, respectively
Provenance 19-Kerns, again from the Central Alps, presented good growth in all three trials and a good recovery (from 0.99 to 1.42 in the Dorna Candrenilor and Zarnesti trials, respectively). It had a great resistance at the Dorna Candrenilor trial (0.98), but less so at the Turda and Zarnesti trials, 0.86 and 0.62, respectively.
Romanian provenance 75-Broşteni had a mean RW above average (2.92) at Dorna Candrenilor and Turda trials, and below average (2.46) at the Zarnesti trial. The resistance was between 0.81 and 0.88, at the Zarnesti and Dorna Candrenilor trials, respectively. Even though the resistance was low, it was above average at the Zarnesti trial, but below average in the other two trials. This provenance managed to obtain a recovery above one in all three trials, and a resilience close to one, 0.97 in the Dorna Candrenilor trial. In the other two trials, it was below average, but close to it, 0.88 compared with 0.89, and 0.77 compared with 0.80, in Turda and Zarnesti trial, respectively.
Provenance 70-Coşna, likewise from Romania, had a great RW in all three trials too. It had a high resistance in the Dorna Candrenilor trial, 0.92, but lower and below average in the Turda trial, 0.85, and lower but above average in the Zarnesti trial, 0.73. In contrast, the recovery of this provenance was high in the Turda and Zarnesti trial, but below one in the Dorna Candrenilor trial (0.92). The resilience was close to the mean value of the trials, but higher only at the Turda trial.
It should be noted that the extreme drought from 2003 and the severe drought from 2002 affected the value of the calculated recovery, resilience, and relative resilience for the extreme event from 2000. Also, the values of resistance, resilience, and relative resilience of the year 2003 were influenced by the extreme event from 2000 and the severe drought from 2002. Further studies should investigate extreme events that are more distant, a fact that is becoming difficult, considering recent events from the 2018–2022 period in Europe [74], or the projections for Romania [75].
Some provenances may represent artificial stands of unknown origin, as is clearly the case for the provenances from Hungary and Northeastern Germany, which are outside the range. Molecular analyses should be conducted in future research to confirm the origin of thsee provenances, especially for the ones selected for being the most resilient. It is important to know the exact origin of the valuable provenances, so that they can be used in breeding programs and assisted migration and for the conservation of the species’ genetic resources.

5. Conclusions

The number of years with extreme drought events has increased in the last two decades in the Romanian Carpathian Mountains. We found high genetic variation in the wood characteristics and drought response of the Norway spruce provenances tested in this study. These results could be further exploited in breeding, seed selection, and genetic conservation programs but also represent a stage in the development of the adaptation strategy of the species.
The drought year and the environmental conditions of the planting site influenced the provenances’ drought response. Norway spruce provenances with high growth rates are the most affected by extreme droughts, a fact highlighted by the negative correlations between RW and resistance in all three trials. Moreover, a prolonged period for recovery was needed after the drought, indicating a lower fitness for drought tolerance of Norway spruce. However, several provenances with high radial growth and good resistance and/or recovery have been identified.
Provenances that performed better in common garden experiments could be used in assisted migration, even in the proximity of the current natural range. The adaptation process could be assisted by relying on assisted migration, especially at the edges of the species distribution, and the genetic diversity would be easier to conserve and increase.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/plants13162172/s1, Table S1: Details about the tested Norway spruce provenances and the provenance trials. Table S2: Details about the Norway spruce provenance trials.

Author Contributions

Conceptualisation, A.M.A.; methodology, A.M.A.; maps, A.M.A.; supervision, G.M. and A.L.C.; project administration, G.M.; funding acquisition, G.M.; writing—all authors. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Romanian Ministry of Research, Innovation, and Digitalisation, in the BIOSERV Nucleu Program, within the framework of project PN 19070303 (Revision of the provenance regions for the production and deployment of forest reproductive materials in Romania to increase the adaptability of forest ecosystems to climate change) and “Programul 1—Dezvoltarea sistemului naţional de cercetare—dezvoltare, Subprogram 1.2—Performanţă instituţională—Proiecte de finanţare a excelenţei în CDI”—project “Creșterea capacității și performanței instituționale a INCDS “Marin Drăcea”) în activitatea de CDI—CresPerfInst” (Contract No. 34PFE./30.12.2021).

Data Availability Statement

Data presented in this study are available from the corresponding author upon request.

Acknowledgments

The authors express their gratitude to the anonymous reviewers for their valuable comments and suggestions, which greatly enhanced this paper.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. AR6 Synthesis Report: Climate Change 2023. Available online: https://www.ipcc.ch/report/ar6/syr/ (accessed on 4 January 2024).
  2. Good, P.; Bärring, L.; Giannakopoulos, C.; Holt, T.; Palutikof, J. Non-Linear Regional Relationships between Climate Extremes and Annual Mean Temperatures in Model Projections for 1961–2099 over Europe. Clim. Res. 2006, 31, 19–34. [Google Scholar] [CrossRef]
  3. Briffa, K.R.; van der Schrier, G.; Jones, P.D. Wet and Dry Summers in Europe since 1750: Evidence of Increasing Drought. Int. J. Climatol. 2009, 29, 1894–1905. [Google Scholar] [CrossRef]
  4. Gudmundsson, L.; Seneviratne, S.I. European Drought Trends. Proc. IAHS 2015, 369, 75–79. [Google Scholar] [CrossRef]
  5. Allen, C.D.; Macalady, A.K.; Chenchouni, H.; Bachelet, D.; McDowell, N.; Vennetier, M.; Kitzberger, T.; Rigling, A.; Breshears, D.D.; Hogg, E.H.; et al. A Global Overview of Drought and Heat-Induced Tree Mortality Reveals Emerging Climate Change Risks for Forests. For. Ecol. Manag. 2010, 259, 660–684. [Google Scholar] [CrossRef]
  6. Anderegg, W.R.L.; Trugman, A.T.; Badgley, G.; Anderson, C.M.; Bartuska, A.; Ciais, P.; Cullenward, D.; Field, C.B.; Freeman, J.; Goetz, S.J.; et al. Climate-Driven Risks to the Climate Mitigation Potential of Forests. Science 2020, 368, eaaz7005. [Google Scholar] [CrossRef]
  7. Ciais, P.; Reichstein, M.; Viovy, N.; Granier, A.; Ogée, J.; Allard, V.; Aubinet, M.; Buchmann, N.; Bernhofer, C.; Carrara, A.; et al. Europe-Wide Reduction in Primary Productivity Caused by the Heat and Drought in 2003. Nature 2005, 437, 529–533. [Google Scholar] [CrossRef] [PubMed]
  8. Lloret, F.; Escudero, A.; Iriondo, J.M.; Martínez-Vilalta, J.; Valladares, F. Extreme Climatic Events and Vegetation: The Role of Stabilizing Processes. Glob. Change Biol. 2012, 18, 797–805. [Google Scholar] [CrossRef]
  9. Reichstein, M.; Bahn, M.; Ciais, P.; Frank, D.; Mahecha, M.D.; Seneviratne, S.I.; Zscheischler, J.; Beer, C.; Buchmann, N.; Frank, D.C.; et al. Climate Extremes and the Carbon Cycle. Nature 2013, 500, 287–295. [Google Scholar] [CrossRef] [PubMed]
  10. Kastridis, A.; Kamperidou, V.; Stathis, D. Dendroclimatological Analysis of Fir (A. Borisii-Regis) in Greece in the Frame of Climate Change Investigation. Forests 2022, 13, 879. [Google Scholar] [CrossRef]
  11. Spinoni, J.; Vogt, J.V.; Naumann, G.; Barbosa, P.; Dosio, A. Will Drought Events Become More Frequent and Severe in Europe? Int. J. Climatol. 2018, 38, 1718–1736. [Google Scholar] [CrossRef]
  12. Senf, C.; Buras, A.; Zang, C.S.; Rammig, A.; Seidl, R. Excess Forest Mortality Is Consistently Linked to Drought across Europe. Nat. Commun. 2020, 11, 6200. [Google Scholar] [CrossRef] [PubMed]
  13. Mateescu, E.; Smarandache, M.; Jeler, N.; Apostol, V. Drought Conditions and Management Strategies in Romania. In Initiative on “Capacity Development to Support. National Drought Management Policy”(WMO, UNCCD, FAO and UNW-DPC). Country Report. 2013, Volume 600. Available online: https://procon.bg/article/drought-conditions-and-management-strategies-romania (accessed on 13 June 2024).
  14. Busuioc, A.; von Storch, H. Changes in the Winter Precipitation in Romania and Its Relation to the Large Scale Circulation, Tellus, 48A, p 538–552. 1996.
  15. Dumitrescu, A.; Amihaesei, V.-A.; Cheval, S. RoCliB– Bias-Corrected CORDEX RCMdataset over Romania. Geosci. Data J. 2022, 10, 262–275. [Google Scholar] [CrossRef]
  16. Lévesque, M.; Saurer, M.; Siegwolf, R.; Eilmann, B.; Brang, P.; Bugmann, H.; Rigling, A. Drought Response of Five Conifer Species under Contrasting Water Availability Suggests High Vulnerability of Norway Spruce and European Larch. Glob. Change Biol. 2013, 19, 3184–3199. [Google Scholar] [CrossRef]
  17. Zang, C.; Pretzsch, H.; Rothe, A. Size-Dependent Responses to Summer Drought in Scots Pine, Norway Spruce and Common Oak. Trees 2012, 26, 557–569. [Google Scholar] [CrossRef]
  18. Jansson, G.; Danusevičius, D.; Grotehusman, H.; Kowalczyk, J.; Krajmerova, D.; Skrøppa, T.; Wolf, H. Norway Spruce (Picea Abies (L.) H.Karst.). In Forest Tree Breeding in Europe: Current State-of-the-Art and Perspectives; Pâques, L.E., Ed.; Managing Forest Ecosystems; Springer: Dordrecht, The Netherlands, 2013; pp. 123–176. ISBN 978-94-007-6146-9. [Google Scholar]
  19. Bouriaud, O.; Popa, I. Comparative Dendroclimatic Study of Scots Pine, Norway Spruce, and Silver Fir in the Vrancea Range, Eastern Carpathian Mountains. Trees 2009, 23, 95–106. [Google Scholar] [CrossRef]
  20. Cochard, H. Vulnerability of Several Conifers to Air Embolism. Tree Physiol. 1992, 11, 73–83. [Google Scholar] [CrossRef] [PubMed]
  21. Schuldt, B.; Buras, A.; Arend, M.; Vitasse, Y.; Beierkuhnlein, C.; Damm, A.; Gharun, M.; Grams, T.E.E.; Hauck, M.; Hajek, P.; et al. A First Assessment of the Impact of the Extreme 2018 Summer Drought on Central European Forests. Basic Appl. Ecol. 2020, 45, 86–103. [Google Scholar] [CrossRef]
  22. Marini, L.; Økland, B.; Jönsson, A.M.; Bentz, B.; Carroll, A.; Forster, B.; Grégoire, J.-C.; Hurling, R.; Nageleisen, L.M.; Netherer, S.; et al. Climate Drivers of Bark Beetle Outbreak Dynamics in Norway Spruce Forests. Ecography 2017, 40, 1426–1435. [Google Scholar] [CrossRef]
  23. Liepe, K.J.; van der Maaten, E.; van der Maaten-Theunissen, M.; Liesebach, M. High Phenotypic Plasticity, but Low Signals of Local Adaptation to Climate in a Large-Scale Transplant Experiment of Picea abies (L.) Karst. in Europe. Front. For. Glob. Change 2022, 5, 804857. [Google Scholar] [CrossRef]
  24. Buras, A.; Menzel, A. Projecting Tree Species Composition Changes of European Forests for 2061–2090 Under RCP 4.5 and RCP 8.5 Scenarios. Front. Plant Sci. 2019, 9, 1986. [Google Scholar] [CrossRef]
  25. Mihai, G.; Teodosiu, M.; Birsan, M.-V.; Alexandru, A.-M.; Mirancea, I.; Apostol, E.-N.; Garbacea, P.; Ionita, L. Impact of Climate Change and Adaptive Genetic Potential of Norway Spruce at the South–Eastern Range of Species Distribution. Agric. For. Meteorol. 2020, 291, 108040. [Google Scholar] [CrossRef]
  26. Mihai, G.; Alexandru, A.-M.; Nita, I.-A.; Birsan, M.-V. Climate Change in the Provenance Regions of Romania over the Last 70 Years: Implications for Forest Management. Forests 2022, 13, 1203. [Google Scholar] [CrossRef]
  27. Rehfeldt, G.E.; Leites, L.P.; Bradley St Clair, J.; Jaquish, B.C.; Sáenz-Romero, C.; López-Upton, J.; Joyce, D.G. Comparative Genetic Responses to Climate in the Varieties of Pinus Ponderosa and Pseudotsuga Menziesii: Clines in Growth Potential. For. Ecol. Manag. 2014, 324, 138–146. [Google Scholar] [CrossRef]
  28. Arend, M.; Kuster, T.; Günthardt-Goerg, M.S.; Dobbertin, M. Provenance-Specific Growth Responses to Drought and Air Warming in Three European Oak Species (Quercus robur, Q. petraea and Q. pubescens). Tree Physiol. 2011, 31, 287–297. [Google Scholar] [CrossRef] [PubMed]
  29. Burczyk, J.; Giertych, M. Response of Norway Spruce (Picea abies [L.] Karst.) Annual Increments to Drought for Various Provenances and Locations. Silvae Genet. 1991, 40, 146–152. [Google Scholar]
  30. George, J.-P.; Grabner, M.; Campelo, F.; Karanitsch-Ackerl, S.; Mayer, K.; Klumpp, R.T.; Schüler, S. Intra-Specific Variation in Growth and Wood Density Traits under Water-Limited Conditions: Long-Term-, Short-Term-, and Sudden Responses of Four Conifer Tree Species. Sci. Total Environ. 2019, 660, 631–643. [Google Scholar] [CrossRef] [PubMed]
  31. Mihai, G.; Alexandru, A.M.; Stoica, E.; Birsan, M.V. Intraspecific Growth Response to Drought of Abies Alba in the Southeastern Carpathians. Forests 2021, 12, 387. [Google Scholar] [CrossRef]
  32. Trujillo-Moya, C.; George, J.-P.; Fluch, S.; Geburek, T.; Grabner, M.; Karanitsch-Ackerl, S.; Konrad, H.; Mayer, K.; Sehr, E.M.; Wischnitzki, E.; et al. Drought Sensitivity of Norway Spruce at the Species’ Warmest Fringe: Quantitative and Molecular Analysis Reveals High Genetic Variation among and within Provenances. G3 Genes Genomes Genet. 2018, 8, 1225–1245. [Google Scholar] [CrossRef] [PubMed]
  33. Atzmon, N.; Moshe, Y.; Schiller, G. Ecophysiological Response to Severe Drought in Pinus Halepensis Mill. Trees of Two Provenances. Plant Ecol. 2004, 171, 15–22. [Google Scholar] [CrossRef]
  34. Schmidtling, R.C. Use of Provenance Tests to Predict Response to Climate Change: Loblolly Pine and Norway Spruce. Tree Physiol. 1994, 14, 805–817. [Google Scholar] [CrossRef]
  35. Vitasse, Y.; Delzon, S.; Bresson, C.C.; Michalet, R.; Kremer, A. Altitudinal Differentiation in Growth and Phenology among Populations of Temperate-Zone Tree Species Growing in a Common Garden. Can. J. For. Res. 2009, 39, 1259–1269. [Google Scholar] [CrossRef]
  36. Gömöry, D.; Longauer, R.; Hlásny, T.; Pacalaj, M.; Strmeň, S.; Krajmerová, D. Adaptation to Common Optimum in Different Populations of Norway Spruce (Picea abies Karst.). Eur. J. For. Res. 2012, 131, 401–411. [Google Scholar] [CrossRef]
  37. Kapeller, S.; Lexer, M.J.; Geburek, T.; Hiebl, J.; Schueler, S. Intraspecific Variation in Climate Response of Norway Spruce in the Eastern Alpine Range: Selecting Appropriate Provenances for Future Climate. For. Ecol. Manag. 2012, 271, 46–57. [Google Scholar] [CrossRef]
  38. Klisz, M.; Ukalska, J.; Koprowski, M.; Tereba, A.; Puchałka, R.; Przybylski, P.; Jastrzębowski, S.; Nabais, C. Effect of Provenance and Climate on Intra-Annual Density Fluctuations of Norway Spruce Picea abies (L.) Karst. in Poland. Agric. For. Meteorol. 2019, 269–270, 145–156. [Google Scholar] [CrossRef]
  39. Sáenz-Romero, C.; Kremer, A.; Nagy, L.; Újvári-Jármay, É.; Ducousso, A.; Kóczán-Horváth, A.; Hansen, J.K.; Mátyás, C. Common Garden Comparisons Confirm Inherited Differences in Sensitivity to Climate Change between Forest Tree Species. PeerJ 2019, 7, e6213. [Google Scholar] [CrossRef] [PubMed]
  40. Zhang, S.; Belien, E.; Ren, H.; Rossi, S.; Huang, J.-G. Wood Anatomy of Boreal Species in a Warming World: A Review. Iforest–Biogeosci. For. 2020, 13, 130. [Google Scholar] [CrossRef]
  41. Denne, M.P. Definition of Latewood According to Mork (1928). IAWA J. 1989, 10, 59–62. [Google Scholar] [CrossRef]
  42. Badea, L. Geografia României; Editura Academiei Republicii Socialiste România: Bucharest, Romania, 1983. [Google Scholar]
  43. Nițu, C.; Benea, V.; Duran, V.; Florescu, I.; Gruescu, A.; Marcu, A.; Răescu, V. Aspecte Privind Variabilitatea Genetică a Unor Proveniențe de Molid. An. Institutul De Cercet. Si Amenaj. Silvice 1974, 31, 49–58. [Google Scholar]
  44. Cybis.Se: Technical Writing, Software Development, Dendrochronology. Available online: https://www.cybis.se/ (accessed on 30 May 2023).
  45. Campelo, F. detrendeR: A Graphical User Interface (GUI) to Visualize and Analyze Dendrochronological Data. 2022. Available online: https://cran.r-project.org/web/packages/detrendeR/index.html (accessed on 23 June 2022).
  46. Bunn, A.; Korpela, M.; Biondi, F.; Campelo, F.; Mérian, P.; Qeadan, F.; Zang, C.; Buras, A.; Cecile, A.; Mudelsee, M.; et al. dplR: Dendrochronology Program Library in R. 2022. Available online: https://CRAN.R-project.org/package=dplR (accessed on 23 June 2022).
  47. R Foundation for Statistical Computing. R: A Language and Environment for Statistical Computing; R core team: Vienna, Austria, 2020; Available online: https://www.R-project.org (accessed on 24 June 2021).
  48. Bucci, G. ClimateDT—Climate Downscaling Tool. Available online: https://www.ibbr.cnr.it/climate-dt/ (accessed on 30 May 2023).
  49. Harris, I.; Osborn, T.J.; Jones, P.; Lister, D. Version 4 of the CRU TS Monthly High-Resolution Gridded Multivariate Climate Dataset. Sci. Data 2020, 7, 109. [Google Scholar] [CrossRef]
  50. Lowe, J.; Bernie, D.; Bett, P.; Bricheno, L.; Brown, S.; Calvert, D.; Clark, R.D.; Karen; Eagle; Edwards, T.; et al. UKCP 18 Science Overview Report November 2018 (Updated March 2019); Met Office Hadley Centre: Exeter, UK, 2019.
  51. Vicente-Serrano, S.M.; Beguería, S.; López-Moreno, J.I. A Multiscalar Drought Index Sensitive to Global Warming: The Standardized Precipitation Evapotranspiration Index. J. Clim. 2010, 23, 1696–1718. [Google Scholar] [CrossRef]
  52. Beguería, S.; Vicente-Serrano, S.M. SPEI: Calculation of the Standardised Precipitation-Evapotranspiration Index. 2017. Available online: https://CRAN.R-project.org/package=SPEI (accessed on 9 May 2022).
  53. Thornthwaite, C.W. An Approach toward a Rational Classification of Climate. Geogr. Rev. 1948, 38, 55–94. [Google Scholar] [CrossRef]
  54. Lloret, F.; Keeling, E.G.; Sala, A. Components of Tree Resilience: Effects of Successive Low-Growth Episodes in Old Ponderosa Pine Forests. Oikos 2011, 120, 1909–1920. [Google Scholar] [CrossRef]
  55. Becker, M.; Nieminen, T.M.; Gérémia, F. Short-Term Variations and Long-Term Changes in Oak Productivity in Northeastern France. The Role of Climate and Atmospheric CO2. Ann. For. Sci. 1994, 51, 477–492. [Google Scholar] [CrossRef]
  56. Jetschke, G.; van der Maaten, E.; van der Maaten-Theunissen, M. Towards the Extremes: A Critical Analysis of Pointer Year Detection Methods. Dendrochronologia 2019, 53, 55–62. [Google Scholar] [CrossRef]
  57. Popa, I. Fundamente Metodologice Şi Aplicaţii de Dendrocronologie; Editura Tehnicǎ Silvicǎ: Bucharest, Romania, 2004. [Google Scholar]
  58. Wu, X.; Liu, H.; Hartmann, H.; Ciais, P.; Kimball, J.S.; Schwalm, C.R.; Camarero, J.J.; Chen, A.; Gentine, P.; Yang, Y.; et al. Timing and Order of Extreme Drought and Wetness Determine Bioclimatic Sensitivity of Tree Growth. Earth’s Future 2022, 10, e2021EF002530. [Google Scholar] [CrossRef]
  59. Kuznetsova, A.; Brockhoff, P.B.; Christensen, R.H.B.; Jensen, S.P. lmerTest: Tests in Linear Mixed Effects Models. 2020. Available online: https://CRAN.R-project.org/package=lmerTest (accessed on 8 December 2022).
  60. Nanson, A. Génétique et Amélioration des Arbres Forestiers; Presses Agronomiques de Gembloux: Gembloux, Belgium, 2004; ISBN 978-2-87016-070-1. [Google Scholar]
  61. Zang, C.; Biondi, F. Treeclim: An R Package for the Numerical Calibration of Proxy-Climate Relationships. Ecography 2015, 38, 431–436. [Google Scholar] [CrossRef]
  62. Bueno de Mesquita, C.P.; White, C.T.; Farrer, E.C.; Hallett, L.M.; Suding, K.N. Taking Climate Change into Account: Non-Stationarity in Climate Drivers of Ecological Response. J. Ecol. 2021, 109, 1491–1500. [Google Scholar] [CrossRef]
  63. Stănescu, V.; Şofletea, N.; Popescu, O.C. Flora Forestiera Lemnoasa a Romaniei; Ceres: Bucharest, Romania, 1997; ISBN 978-973-40-0383-9. [Google Scholar]
  64. Baldi, P.; La Porta, N. Toward the Genetic Improvement of Drought Tolerance in Conifers: An Integrated Approach. Forests 2022, 13, 2016. [Google Scholar] [CrossRef]
  65. Forner, A.; Valladares, F.; Aranda, I. Mediterranean Trees Coping with Severe Drought: Avoidance Might Not Be Safe. Environ. Exp. Bot. 2018, 155, 529–540. [Google Scholar] [CrossRef]
  66. McDowell, N.; Pockman, W.T.; Allen, C.D.; Breshears, D.D.; Cobb, N.; Kolb, T.; Plaut, J.; Sperry, J.; West, A.; Williams, D.G.; et al. Mechanisms of Plant Survival and Mortality during Drought: Why Do Some Plants Survive While Others Succumb to Drought? New Phytol. 2008, 178, 719–739. [Google Scholar] [CrossRef]
  67. Eilmann, B.; Rigling, A. Tree-Growth Analyses to Estimate Tree Species’ Drought Tolerance. Tree Physiol. 2012, 32, 178–187. [Google Scholar] [CrossRef] [PubMed]
  68. Moran, E.; Lauder, J.; Musser, C.; Stathos, A.; Shu, M. The Genetics of Drought Tolerance in Conifers. New Phytol. 2017, 216, 1034–1048. [Google Scholar] [CrossRef] [PubMed]
  69. Gómez-Aparicio, L.; García-Valdés, R.; Ruíz-Benito, P.; Zavala, M.A. Disentangling the Relative Importance of Climate, Size and Competition on Tree Growth in Iberian Forests: Implications for Forest Management under Global Change. Glob. Change Biol. 2011, 17, 2400–2414. [Google Scholar] [CrossRef]
  70. Čermák, P.; Kolar, T.; Žid, T.; Trnka, M.; Rybníček, M. Norway Spruce Responses to Drought Forcing in Area Affected by Forest Decline. For. Syst. 2019, 28, e016. [Google Scholar] [CrossRef]
  71. Augustaitis, A.; Augustaitienė, I.; Baugarten, M.; Bičenkienė, S.; Girgždienė, R.; Kulbokas, G.; Linkevičius, E.; Marozas, V.; Mikalajūnas, M.; Mordas, G.; et al. Tree-Ring Formation as an Indicator of Forest Capacity to Adapt to the Main Threats of Environmental Changes in Lithuania. Sci. Total Environ. 2018, 615, 1247–1261. [Google Scholar] [CrossRef] [PubMed]
  72. Linkevičius, E.; Kliučius, A.; Šidlauskas, G.; Augustaitis, A. Variability in Growth Patterns and Tree-Ring Formation of East European Scots Pine (Pinus sylvestris L.) Provenances to Changing Climatic Conditions in Lithuania. Forests 2022, 13, 743. [Google Scholar] [CrossRef]
  73. Roibu, C.-C.; Sfeclă, V.; Mursa, A.; Ionita, M.; Nagavciuc, V.; Chiriloaei, F.; Leșan, I.; Popa, I. The Climatic Response of Tree Ring Width Components of Ash (Fraxinus excelsior L.) and Common Oak (Quercus robur L.) from Eastern Europe. Forests 2020, 11, 600. [Google Scholar] [CrossRef]
  74. Knutzen, F.; Averbeck, P.; Barrasso, C.; Bouwer, L.M.; Gardiner, B.; Grünzweig, J.M.; Hänel, S.; Haustein, K.; Johannessen, M.R.; Kollet, S.; et al. Impacts and Damages of the European Multi-Year Drought and Heat Event 2018–2022 on Forests, a Review. EGUsphere 2023, 1–56. [Google Scholar] [CrossRef]
  75. Antonescu, B.; Ene, D.; Boldeanu, M.; Andrei, S.; Mărmureanu, L.; Marin, C.; Pîrloagă, R. Future Changes in Heatwaves Characteristics in Romania. Theor. Appl. Clim. 2023, 153, 525–538. [Google Scholar] [CrossRef]
Plants 13 02172 g001
Figure 1. The locations of the Norway spruce provenances (dots and numbers) and of the provenance trials (red triangles).
Figure 1. The locations of the Norway spruce provenances (dots and numbers) and of the provenance trials (red triangles).
Plants 13 02172 g001
Plants 13 02172 g002
Figure 2. The variation in mean annual temperature (MAT) and annual precipitations (AP) for the Norway spruce provenance trials for the 1972–2020 period; (A)—Dorna Candrenilor, (B)—Turda, (C)—Zarnesti.
Figure 2. The variation in mean annual temperature (MAT) and annual precipitations (AP) for the Norway spruce provenance trials for the 1972–2020 period; (A)—Dorna Candrenilor, (B)—Turda, (C)—Zarnesti.
Plants 13 02172 g002
Plants 13 02172 g003
Figure 3. Variation in SPEI-3 for the Norway spruce provenance trials, for the period 1972–2020; (A)—Dorna Candrenilor, (B)—Turda, (C)—Zarnesti.
Figure 3. Variation in SPEI-3 for the Norway spruce provenance trials, for the period 1972–2020; (A)—Dorna Candrenilor, (B)—Turda, (C)—Zarnesti.
Plants 13 02172 g003
Plants 13 02172 g004
Figure 4. The pointer years in each Norway spruce provenance trial. A nature value of 1 means the year had a positive influence on the trees’ ring growth, while a value of −1 means the year had a negative influence on the trees’ ring growth.
Figure 4. The pointer years in each Norway spruce provenance trial. A nature value of 1 means the year had a positive influence on the trees’ ring growth, while a value of −1 means the year had a negative influence on the trees’ ring growth.
Plants 13 02172 g004
Plants 13 02172 g005
Figure 5. Variation in the wood traits of Norway spruce over the years for the period 1981–2020; (A)—Dorna Candrenilor, (B)—Turda, (C)—Zarnesti. RWI –ring width index; EWI –earlywood width index; LWI—latewood width index; LWP—latewood percentage.
Figure 5. Variation in the wood traits of Norway spruce over the years for the period 1981–2020; (A)—Dorna Candrenilor, (B)—Turda, (C)—Zarnesti. RWI –ring width index; EWI –earlywood width index; LWI—latewood width index; LWP—latewood percentage.
Plants 13 02172 g005
Plants 13 02172 g006
Figure 6. Variation in tree resilience indices of Norway spruce provenances for the year 2000.
Figure 6. Variation in tree resilience indices of Norway spruce provenances for the year 2000.
Plants 13 02172 g006
Plants 13 02172 g007
Figure 7. The response and correlation function for the Dorna Candrenilor trial of Norway spruce provenances. * significant at 5%.
Figure 7. The response and correlation function for the Dorna Candrenilor trial of Norway spruce provenances. * significant at 5%.
Plants 13 02172 g007
Plants 13 02172 g008
Figure 8. The response and correlation function for the Zarnesti trial of Norway spruce provenances. * significant at 5%.
Figure 8. The response and correlation function for the Zarnesti trial of Norway spruce provenances. * significant at 5%.
Plants 13 02172 g008
Plants 13 02172 g009
Figure 9. The response and correlation function for the Turda trial of Norway spruce provenances. * significant at 5%.
Figure 9. The response and correlation function for the Turda trial of Norway spruce provenances. * significant at 5%.
Plants 13 02172 g009
Table 1. Years and months with severe and extreme drought events for each Norway spruce trial. The bold represents the extreme events, based on SPEI-3.
Table 1. Years and months with severe and extreme drought events for each Norway spruce trial. The bold represents the extreme events, based on SPEI-3.
YearMonths with Severe and Extreme Droughts
Dorna CandrenilorTurdaZarnesti
1972MarchMarchMarch, April
1973JanJanuary, February-
1974Febr.–AprilFebr.–AprilFebr.–April
1975MarchMarchMarch
1977--December
1978--January
1982NovemberNovemberNovember
1983DecemberDecemberApril
1986May, NovemberNovemberMay, October
1989MarchMarchMarch
1990Febr., MarchMayMarch
1992-May, SeptemberSeptember
1994JulyJuly-
1996--July
2000June, DecemberJune–August, October–DecemberJune, July, August, December
2001January-January
2002Febr., May, JuneFebr., May, Junefebr., May, June
2003May, JuneMay, JuneMay, June, July, August
2006-DecemberDecember
2007June-June, July
2008-February-
2009May, SeptemberMay-
2011Oct., NovemberOctober, NovemberOctober, November
2012Aug, SeptemberAugust–OctoberAugust, September
2013December--
2015July, August, September-July
2017-February-
2018May, OctoberSeptember–NovemberMay, October
2019August-DecemberOctoberAugust, September, November, December
2020JanuaryJanuaryJanuary
Table 2. Descriptive statistics of the Norway spruce series at 49 years after planting.
Table 2. Descriptive statistics of the Norway spruce series at 49 years after planting.
rbarProvenance rbar minProvenance rbar maxMean RW (mm)Min RW
(mm)		Max RW
(mm)
Dorna Candrenilor0.9080.802, prov 900.961, prov 542.732.15, prov 833.11, prov 25
Turda0.8760.729, prov 670.941, prov 752.812.31, prov 833.92, prov 10
Zarnesti0.8760.705, prov 180.942, prov 422.552.13, prov 943.18, prov 55
rbar—the mean series intercorrelation; RW—ring width.
Table 3. The results for the random and fixed effects of the wood width and LWP in each Norway spruce provenance trial evaluated at age 49.
Table 3. The results for the random and fixed effects of the wood width and LWP in each Norway spruce provenance trial evaluated at age 49.
TraitLRTpLRT Year x provMS BMS YearMean ± SD
ZarnestiRW 711.55 ***66.98 ***91.56 ***1116.66 ***2.55 ± 0.19
EW677.85 ***84.63 ***38.94 ***739.12 ***1.90 ± 0.17
LW376.41 ***0.0011.12 ***46.88 ***0.65 ± 0.05
LWP510.66 ***0.004583.1 ***13,944.9 ***27.78 ± 1.65
Dorna CandrenilorRW712.85 ***12.42 ***47.88 ***1996.86 ***2.73 ± 0.18
EW579.62 ***30.60 ***12.07 ***1310.91 ***2.08 ± 0.14
LW639.68 ***1.29 16.981 ***65.855 ***0.66 ± 0.07
LWP491.92 ***26.37 ***3207.5 ***20,586.8 ***27.06 ± 1.54
TurdaRW1259.9 ***0.002.091315.16 ***2.81 ± 0.29
EW1003.0 ***0.002.28 *835.49 ***2.03 ± 0.22
LW996.03 ***0.000.66 **56.16 ***0.78 ± 0.08
LWP326.7 ***0.001389.3 ***15,916 ***30.76 ± 1.62
*, **, *** Significant at 5%, 1%, and 0.1%, respectively; LRTp—likelihood ratio test for provenance random effect; LRT year x prov—likelihood ratio test for provenance x year random effect; MS B—mean squares for block effect; MS Year—mean squares for year effect; RW—ring width; EW—earlywood width; LW—latewood width; LWP—latewood percentage.
Table 4. The results for the random and fixed effects of the tree resilience indices in each Norway spruce provenance trial evaluated in their common extreme years.
Table 4. The results for the random and fixed effects of the tree resilience indices in each Norway spruce provenance trial evaluated in their common extreme years.
YearProvenance TrialTraitLRTpMS B
2000Dorna CandrenilorResistance0.084 ns0.118 *
Recovery2.593 ns0.011 ns
Resilience5.418 *0.063 ns
Rel. resilience1.966 ns0.009 ns
TurdaResistance7.859 **0.124 *
Recovery8.014 **0.0009 ns
Resilience1.850 ns0.120 *
Rel. resilience9.308 **0.0002 ns
ZarnestiResistance25.268 ***0.254 ***
Recovery38.035 ***0.666 ***
Resilience7.015 **0.048 ns
Rel. resilience32.521 ***0.458 ***
2003Dorna CandrenilorResistance2.604 ns0.192 ***
Recovery0 ns0.073 ns
Resilience0.429 ns0.234 ***
Rel. resilience5.12 × 10−13 ns0.029 ns
TurdaResistance2.297 ns0.087 ns
Recovery0.021 ns0.371 *
Resilience0.879 ns0.123 *
Rel. resilience1.42 × 10−12 ns0.308 **
ZarnestiResistance14.843 ***0.219 **
Recovery24.834 ***0.576 ***
Resilience14.084 ***0.064 ns
Rel. resilience22.538 ***0.383 **
*, **, *** Significant at 5%, 1%, and 0.1%, respectively; ns—non-significant; LRTp—likelihood ratio test for provenance random effect; MS B—mean squares for block effect.
Table 5. Correlation coefficients between wood traits and geographic coordinates of Norway spruce provenances at 49 years after planting.
Table 5. Correlation coefficients between wood traits and geographic coordinates of Norway spruce provenances at 49 years after planting.
Provenance Trial LWLWPRWLatitudeLongitudeElevation
Dorna CandrenilorEW0.309 ***−0.445 ***0.922 ***−0.185 **−0.0430.179 ***
LW 0.625 ***0.653 ***−0.159 *−0.162 *0.220 **
LWP −0.099−0.005−0.135 *0.034
RW −0.212 *−0.1000.232 *
ZarnestiEW0.452 ***−0.468 ***0.963 ***0.048−0.010−0.065
LW 0.502 ***0.676 ***−0.048−0.042−0.016
LWP −0.235 **−0.131 *−0.0200.057
RW 0.025−0.021−0.059
TurdaEW0.753 ***−0.468 ***0.983 ***−0.068−0.130 *0.010
LW 0.134 *0.861 ***−0.036−0.0580.065
LWP −0.324 ***0.0700.155 *0.029
RW −0.062−0.1170.026
*, **, *** Significant at 5%, 1%, and 0.1%, respectively; RW—ring width; EW—earlywood width; LW—latewood width; LWP—latewood percentage.
Table 6. Correlation coefficients between mean tree resilience indices and geographic coordinates of Norway spruce provenances at 49 years after planting.
Table 6. Correlation coefficients between mean tree resilience indices and geographic coordinates of Norway spruce provenances at 49 years after planting.
Provenance Trial LatitudeLongi-TudeElevation
Dorna CandrenilorResistance−0.037−0.1140.128 *
Recovery0.131 *0.038−0.132 *
Resilience0.074−0.0990.043
Rel. resilience0.1230.025−0.103
ZarnestiResistance0.215 ***0.120−0.121
Recovery−0.236 ***−0.317 ***0.111
Resilience−0.223 ***−0.277 ***0.124
Rel. resilience−0.239 ***−0.312 ***0.117
TurdaResistance0.1030.167 **0.047
Recovery−0.099−0.151 *0.048
Resilience0.0780.0190.049
Rel. resilience−0.047−0.152 *−0.012
*, **, *** Significant at 5%, 1%, and 0.1%, respectively.
Table 7. Correlation coefficients between mean tree resilience indices and wood traits of Norway spruce provenances at 49 years after planting.
Table 7. Correlation coefficients between mean tree resilience indices and wood traits of Norway spruce provenances at 49 years after planting.
Provenance Trial ResistanceRecoveryResilienceRel. Resilience
Dorna CandrenilorEW−0.255 ***0.153 *−0.175 **0.106
LW0.108−0.0970.018−0.106
LWP0.289 ***−0.207 **0.150 *−0.173 **
RW−0.159 *0.082−0.132 *0.042
ZarnestiEW−0.0070.063−0.0430.077
LW−0.0350.137 *0.0840.141 *
LWP−0.0350.0520.0890.034
RW−0.0160.094−0.0110.106
TurdaEW−0.370 ***0.347 ***−0.0500.330 ***
LW−0.267 ***0.277 ***−0.0040.260 ***
LWP0.192 **−0.0650.141 *−0.090
RW−0.360 ***0.345 ***−0.0390.327 ***
*, **, *** Significant at 5%, 1%, and 0.1%, respectively; RW—ring width; EW—earlywood width; LW—latewood width; LWP—latewood percentage.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Alexandru, A.M.; Mihai, G.; Stoica, E.; Curtu, A.L. Tree Resilience Indices of Norway Spruce Provenances Tested in Long-Term Common Garden Experiments in the Romanian Carpathians. Plants 2024, 13, 2172. https://doi.org/10.3390/plants13162172

AMA Style

Alexandru AM, Mihai G, Stoica E, Curtu AL. Tree Resilience Indices of Norway Spruce Provenances Tested in Long-Term Common Garden Experiments in the Romanian Carpathians. Plants. 2024; 13(16):2172. https://doi.org/10.3390/plants13162172

Chicago/Turabian Style

Alexandru, Alin Madalin, Georgeta Mihai, Emanuel Stoica, and Alexandru Lucian Curtu. 2024. "Tree Resilience Indices of Norway Spruce Provenances Tested in Long-Term Common Garden Experiments in the Romanian Carpathians" Plants 13, no. 16: 2172. https://doi.org/10.3390/plants13162172

APA Style

Alexandru, A. M., Mihai, G., Stoica, E., & Curtu, A. L. (2024). Tree Resilience Indices of Norway Spruce Provenances Tested in Long-Term Common Garden Experiments in the Romanian Carpathians. Plants, 13(16), 2172. https://doi.org/10.3390/plants13162172

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop