Next Article in Journal
Broad-Leaved Tree Growth Modulated by Industrial Air Pollution in the Northern Romania (Baia Mare Region)
Previous Article in Journal
Differed Adaptive Strategies to Nutrient Status between Native and Exotic Mangrove Species
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Growth and Adaptive Capacity of Douglas Fir Genetic Resources from Western Romania under Climate Change

1
Department of Forest Genetics and Tree Breeding, “Marin Dracea” National Institute for Research and Development in Forestry, 077190 Bucharest, Romania
2
Faculty of Silviculture and Forest Engineering, Transilvania University of Brasov, 500036 Brasov, Romania
3
Department of Research and Meteo Infrastructure Projects, Meteo Romania (National Meteorological Administration), 013686 Bucharest, Romania
*
Author to whom correspondence should be addressed.
Submission received: 25 April 2022 / Revised: 13 May 2022 / Accepted: 18 May 2022 / Published: 20 May 2022
(This article belongs to the Special Issue Introduced Forest Tree Species and Their Genetic Diversity)

Abstract

:
The most recent climate change scenarios show that Southern and Eastern Europe will be affected by a significant increase in temperature and drought frequency by the end of the 21st century. Romania has already recorded very high temperatures and long periods of drought over recent decades, the most affected regions being the south, west and east of the country. Considering that successful forest management requires suitable species and high-quality reproductive material for reforestation, the aim of this study was to evaluate and compare the growth and drought response of Douglas fir (Pseudotsuga menziessi var. menziesii) and Norway spruce (Picea abies) in two stands installed at the end of the 19th century in western Romania. The growth traits, wood characteristics and drought parameters (resistance, recovery, resilience and relative resilience) of Douglas fir and Norway spruce trees have been analyzed and compared. The climate–growth relationship was determined using growth response functions over the period 1938–2017. Additionally, to simulate the potential impact of climate change on Douglas fir in this region, the RCP4.5 scenario was used over two periods: 2041–2070 and 2071–2100. The results reveal that Douglas fir has an exceptional growth capacity, overcoming the Norway spruce since the early ages in both site conditions. The highest growth performances were seen in the low-productivity site. From analyzing the responses to drought events, considerable differences were found between species. The results highlight the high resistance and relative resilience to extreme droughts of Douglas fir compared to Norway spruce. However, autumn–winter temperatures play an important role in the adaptation of Douglas fir to site conditions in Romania. The use of appropriate provenances of Douglas fir in mixed stands with native broadleaved species may be an option for climatically exposed sites, thus increasing the value of these stands.

1. Introduction

The Douglas fir (Pseudotsuga menziessi) is one of the most important forest tree species in North America. Since the middle of the 19th century, it has been successfully introduced into many regions with a temperate climate in both hemispheres [1]. In Europe, Douglas fir was introduced in 1827 and currently covers the largest area outside its natural range, being present in 35 countries [2]. In Romania, the first Douglas fir plantations were established in 1887 in eastern Romania (Moldavia) and one year later, in 1888, in the western part of the country (Banat) [3,4]. Although it is one of the most important non-native tree species, it covers only 0.40% of the total forest area in Europe (approximately 823,534 ha) [2] and 0.12% (7300 ha) in Romania (NFI 2016).
Douglas fir has a wide distribution area in North America, where it is represented by two geographically distinct varieties: coastal (P. menziessi var. menziesii (Mirbel) Franco) and interior (P. menziessi var. glauca (Beissn.) Franco). The coastal variety of Douglas fir is the most widespread in Europe because it was more suitable to climate conditions than the interior variety. The costal variety ranges along the North American Pacific coast (California, Oregon, Washington, and British Columbia) from the sea level to about 1500 m and is adapted to an oceanic climate characterized by mild, wet winters and cool, relatively dry summers, and long frost-free growing seasons. Compared to the climatic conditions in Europe, most of the precipitation occurs in winter, with the annual precipitation rates ranging from 760 to 3000 mm [1,5].
Results from long-term provenance trials in the natural range indicate high adaptive genetic variation among Douglas fir populations and clinal variation patterns in growth, phenology, and cold hardiness [6,7,8,9]. Douglas fir can be referred to as an “adaptive specialist” because its populations are closely adapted to their ecological niches [10,11]. The coastal Douglas fir populations reach an increased height, diameter, and volume compared with the interior variety, and they are more tolerant to needle diseases but less tolerant to fall frost, winter cold, and drought. Winter temperatures and frost date are the major factors involved in the adaptation of the coastal Douglas fir, whereas summer drought is of less importance [12]. Furthermore, there is a large variability for wood characteristics within the natural range of the species. The wood of the northern populations is considered more valuable (denser and stronger) than that of the southern ones [13,14].
Since the 20th century, a large number of Douglas fir provenance trials have been established at the national and European levels, indicating a substantial level of interest for the species. Thus, after the Second World War, it became a major tree species used in reforestation in Western Europe. In most European countries, it showed higher productivity compared with the native conifers and good resistance against fungal pathogens and low numbers of pests and diseases [15,16].
Subsequently, the results of International Union of Forest Research Organizations (IUFRO) experiments on Douglas fir provenances in Europe show its sensitivity to winter drought and frosts and the importance of the geographic origin of forest reproductive material [17,18]. Thus, the interior populations (particularly in British Columbia) are the most resistant to early and winter frosts but susceptible to late frosts [19,20]. Coastal populations, with late flushing, are the most resistant to late frost but rather sensitive to early and winter frosts [15]. Further, Breidenstein [21], in a synthesis of the results from 108 Douglas fir provenance tests established by IUFRO (1967) in 15 European countries, has shown that low-elevation coastal and Cascade provenances from Washington State underwent the most rapid growth over most sites and broad adaptability, although mortality was high on several colder sites [22]. Some northern Oregon provenances from west of the Cascades as well as a few southwestern British Columbia sources also had high productivity across much of Europe. Only in continental climates (e.g., in Sweden, Finland, and the Czech Republic), the interior variety outperforms the coastal one. The highest mortality rates were found for the southern coastal Oregon provenances on the coldest sites and the interior of British Columbia provenances on the mildest sites. The success in reforestation work is mostly dependent on the genetic selection of the appropriate forest reproductive material for planting sites.
In Romania, Douglas fir breeding activities started in the 1970s and focused on assessing the genetic variability in provenance trials [23], establishing first-generation seed orchards [24], and preserving valuable stands as genetic resources [25]. Unfortunately, the breeding program was abandoned in the last two decades, and recommendations on the use of appropriate provenances in Romania still need scientific validation. There is also a lack of knowledge regarding the growth potential of Douglas fir in Romania, relationships between growth and climatic variables, the selection of suitable seed sources, and species response to drought and rapid climate change.
Climate change scenarios show that Southern and Eastern Europe will be affected by a significant increase in temperature and drought frequency by the end of the 21st century [26]. Over recent decades, Romania has experienced record temperatures [27], decreasing wind speed [28,29] and snow depth [30] changes in the frequency of hydrometeors [31], as well as drought periods, which were associated with global climate change [32,33,34,35] or changes in atmospheric circulation [36,37,38,39]. Climate change will certainly have detrimental effects on forest ecosystems in these regions [40,41,42,43]. Norway spruce, one of the most planted conifer species outside the natural distribution area, showed a high sensitivity to decreases in water availability and weak adaptive capacity at lower elevations [44,45]. However, recent studies showed that Douglas fir might be an alternative species because some provenances are well adapted to drought and warmer conditions [46,47,48]. Therefore, interest in using Douglas fir at lower altitudes or in response to climate change has increased in Europe in recent years.
The aim of this study was to evaluate the growth response and adaptive capacity of Douglas fir genetic resources from western Romania in the context of climate change. These populations were established at the end of the 19th century and are mixed stands with Douglas fir, Norway spruce, and broadleaved species. They are among the oldest Douglas fir plantations in Romania and were designated as forest genetic resources for Douglas fir to be included in the National Catalog of Forest Genetic Resources [49]. Based on the assumption that global warming will reduce the growth and adaptive capacity of Norway spruce at lower elevations in Southeastern Europe, the objectives were to (1) comparatively assess the growth and adaptive capacity of Douglas fir and Norway spruce, (2) evaluate the Douglas fir response to drought, (3) determine climate–growth relationships, (4) evaluate the Douglas fir suitability to future climate conditions, and (5) provide recommendations for the use of Douglas fir.
Considering that sustainable forest management requires productive species and high-quality reproductive material for reforestation, the need to identify valuable seed sources has emerged. Furthermore, knowing the provenance of seed sources or plant material is particularly important for the success of reforestation in response to climate change. This study aimed to provide further information and offer solutions for the sustainable management of forest ecosystems and forest adaptation strategies to climate change.

2. Materials and Methods

2.1. Study Sites

This study was conducted in two of the most valuable forest genetic resources (FGR) of Douglas fir located in the Ana Lugojana Forest District (Banat region). These stands are 105 and 130 years old, respectively, and are mixed stands with Douglas fir, Norway spruce, and broadleaved species. The two stands are located in highly contrasting site conditions: FGR 1 in high-productivity site conditions for European beech at an altitude of 880 m (45°35′ N, 22°15′ E), whereas FGR 2 is in low-productivity site conditions at an altitude of 460 m (45°35′ N, 22°25′ E). The characteristics of the two stands are presented in Table 1.

2.2. Field Measurements and Analyses

In each stand, 20 dominant or (co-)dominant trees per species of 105 and 130 years old, respectively, were measured for total height and diameter at 1.30 m (DBH) and cored at breast height using 5 mm increment borers (Haglof, Sweden). Cores were prepared using standard dendrochronological methods [50]. Then, each core sample was scanned at 1200 dpi, using an Epson Expression 10,000 XL, and the wood characteristics were measured using the Ligno Vision software package to the nearest 0.001 mm. The assessed characteristics were: ring width (RW), earlywood width (EW), latewood width (LW), and latewood percentage (LWP), as an indication of wood quality.
For each study site, a master series was constructed and cross-dated using COFECHA [51] to avoid dating errors due to missing or false rings, which could be present in an increment radial core. Only dendrochronological series that presented intercorrelation values >0.328 (p < 0.01) were included in final tree ring data. All tree-ring time series were standardized to a mean value of one to obtain a dimensionless ring-width index (RWI), thus preserving a large portion of low-frequency variability due to the influence of climatic events [52,53]. The negative exponential regression in the R package (dplR) [54] was applied for each raw measurement series because it is deterministic, meaning that it follows a model of tree growth. The final data set consisted of data from 38 Douglas fir trees and 32 Norway spruce trees. The analyzed period was 1938–2017 for each tree-ring series.

2.3. Climate Data

Climatic data series for the study sites were extracted from the CRU TS (Climatic Research Unit gridded Time Series) dataset v.4.04, over the period 1901–2017. CRU TS is a widely used dataset made by interpolating the monthly climate anomalies from weather station measurements at a spatial resolution of 0.5° × 0.5°, covering all continents except Antarctica [55].
High-resolution gridded climate data over Romania (1961–2020) were derived from daily precipitation and temperature grids over the period 1961–2020. The dataset consists of a spatial resolution of 1 km × 1 km and was made using state-of-the-art interpolation techniques for an improved reproduction of the spatial climatic variability (e.g., [56,57]).
Climate projections were derived from the RoCliB dataset [58,59], which includes air temperature and precipitation data from 10 general circulation models dynamically downscaled by several regional climate models and bias corrected (adjusted) over Romania for the 1971–2100 period. We used the 10-model ensemble data for the scenarios RCP 4.5 and RCP 8.5 for the periods 2041–2070 and 1971–2100.
The following climatic variables of FGR location were calculated: mean annual temperature (MAT); mean temperature of the growing season (April to September) (MTVEG); mean temperatures for January (MTJAN) and July (MTJUL) (i.e., the coldest and warmest months, respectively); mean temperature from October to December of the previous year (MTOCT-DEC); mean temperature from October of the previous year to March of the current year (MTOCT-MAR); mean temperature from January to March of the current year (MTJAN-MAR); sum annual precipitation (SAP); sum precipitation during the growing season (SPVEG); sum precipitation of the coldest (SPJAN) and warmest (SPJUL) months; sum precipitation from October to December of the previous year (SPOCT-DEC); sum precipitation from January to March of the current year (SPJAN-MAR); and sum precipitation from October of the previous year to March of the current year (SPOCT-MAR).
According to the Köppen and Geiger climate classification [60], the site’s climate is warm and temperate (Cfb), with slight Mediterranean influences being representative of vegetation conditions in the west of the country. The climate is characterized by mild winters and a richer rainfall regime. Comparative climatic data for Douglas fir populations from western Romania and for the natural distribution range are presented in Table 2.

2.4. Determination of Drought Events and Drought Response Parameters

As an indicator for meteorological droughts, we calculated the Standardized Precipitation Evapotranspiration Index (SPEI) [61], which is based on precipitation and the potential evapotranspiration (PET) over the period 1901–2019. The PET was calculated according to the Thornthwaite equation [62]. Given that an extreme drought event must last for a minimum of 2 to 3 months, we calculated the SPEI at a time scale of 3 months (SPEI-3), with the R package “SPEI” [63]. This approach allowed us to detect both seasonal and annual variations of drought events during the analyzed period. The drought years were classified as follows: SPEI ≤ −2 extreme drought year; SPEI between −1.99 and −1.50: severe drought year; SPEI between −1.49 and −1.0: moderate drought year; SPEI between −1.0 and + 1.0: normal precipitation year.
The species response to drought events was evaluated using four drought parameters [64]: resistance (Res), recovery (Rec), resilience (Rsl), and relative resilience (rRsl). Resistance (Res) was calculated as the ratio between ring width during drought (Dr) and before the drought event (preDr)—Res = Dr/preDr—and indicates how much the radial growth decreased during drought (Res = 1 means high tolerance; Res < 1 means low tolerance).
Recovery (Rec) was calculated as the ratio between the ring width after the drought event (postDr) and during drought (Rc = postDr/Dr) and indicates the revitalization capacity after a drought period.
Resilience (Rsl) represents the ratio of the ring width after drought (postDr) and pre-drought (preDr)—Rsl = postDr/preDr—and describes the species capacity to reach pre-drought increment after a drought event (Rsl = 1 means full restoration; Rs < 1 means long-term growth reductions). Relative resilience (rRsl) was calculated by rRsl = (postDr-Dr)/preDr. Pre- and post-drought ring widths were calculated as average values for a three-year period before or after the drought year.

2.5. Data Analysis

Analysis of variance were performed using the Generalized Linear Mixed Models (GLM) procedure (SPSS v19). The total amount of variation was divided into tree species, years, FGR (sites), and the species site interaction. Apart from the FGR, which was considered fixed, all effects were considered random. The assumptions of the model were checked by a Shapiro and Wilk test for normality and by Levene’s test for homogeneity.
ANOVA was performed as described in the following mixed model:
Yijk = μ + Si + Rj + Yl + BFij + eijlk
where: Xijk = performance of kth tree in ith species in jth RGF; μ = overall mean; Si = effect of ith species; Rj = effect of jth RGF; Yl = effect of lst year; SYil = interaction of ith species and lst year; eijlk = random error associated with ijlkth trees.
The relationships between climatic variables of the plantation sites and ring-width indices and latewood percentage of Douglas fir and Norway spruce were evaluated through response function analysis. The quadratic models based on both temperature and precipitation as predictor variables were applied in the SPSS program (stepwise selection method), considering them more suitable [64,65,66]. We used seven temperature and seven precipitation variables, and the best models were chosen based on the R2 coefficient. Response functions analysis was performed for the period 1938–2017 using the formula of Wang:
Yij = β0 + β1Tnj + β2Pnj3Tnj2 + β4Pnj2 + ej
where Yij is the observation of the population i at the site j; βs are the intercept and regression coefficients; Tnj and Pnj are the temperature and precipitation variables, respectively, at the site j; and ej is the residual.
To assess the potential distribution area of Douglas fir in Romania, the species climate envelope was developed based on climate predictors that reflect species ecological requirements [67]. The predictors included the yearly (annual) precipitation amount, the sum of precipitation of the coldest month and the mean annual temperature for the period 1961–2020. The same ecological indicators were also used for modelling the species suitability to future environmental conditions in Romania under the RCP 4.5 and RCP 8.5 scenarios by 2100—for the intervals 2041–2070 and 2071–2100.

3. Results

3.1. Variation in Growth and Wood Characteristics

The analysis of variance for each FGR and trait are presented in Table 3. There were significant differences for all traits (p < 0.001) within and between species for each site (FGR) and among sites as well. There was also a highly significant year effect for traits and species x year interaction in both study sites.
In the case of Douglas fir, the total height ranged from 45.70 m to 54.70 m in FGR 1 and 41.80 m to 57.60 m in FGR 2. DBH ranged from 70.38 cm to 100.34 cm in RGF 1 and 54.78 cm to 129.62 cm in FGR 2. The average of volume/tree was 11.9 m3 in FGR 1 and 15.6 m3 in FGR 2. The best growth performances were obtained by both conifer species in FGR 2, located on low-productivity site conditions for European beech. Overall, Douglas fir exceeded the Norway spruce for all growth characteristics, with percentages on average between 24 and 27% for total height, 56% for DBH, and between 75 and 77% for volume per tree (Figure 1).
The ring width and latewood percentage varied significantly across species, site, and year. The ring width for Douglas fir ranged from 0.14 mm to 8.21 mm along all years and sites, whereas for Norway spruce, it ranged from 0.11 mm to 6.20 mm. Averaged by site and species, the average ring width for Douglas fir ranged from 2.74 mm in FGR 1 to 2.02 mm in FGR 2 and from 1.64 mm in FGR 1 to 1.20 mm in FGR 2 for Norway spruce. The overall average values per species were 2.39 mm for Douglas fir and 1.42 mm for Norway spruce.
The average latewood percentage ranged from 6 to 96% for Douglas fir and from 7 to 43% for Norway spruce. The overall average values per species were 53% for Douglas fir and 43% for Norway spruce. The average ring widths and latewood percentage were about 41% and 19% higher, respectively, for Douglas fir than those of Norway spruce for the period of 1938–2017. In both study sites, ring width decreased with increasing age, but growing patterns differed between species (Figure 2 and Figure 3). Thus, at the same age, 105 and 130 years, respectively, the ring width of Norway spruce was 50% lower in FGR 1 and 31% in FGR 2 compared with Douglas fir.

3.2. Climate Variation and Identification of Drought Years

The two FGRs are located very close to each other and, therefore, a single climate database was used. The gridded climate data analysis revealed considerable changes in climate conditions of the FGR sites over the analyzed period (Figure 4). MAT varied between 6.8 °C in 1940 and 10.9 °C in 2014. SAP ranged from 467 mm in 2000 to 1055 mm in 2010, whereas SPVEG varied from 284 mm in 2000 to 686 mm in 2010. The most significant changes occurred over the last two decades; MAT increased by 0.9 °C, TMVEG by 1.0 °C, and TMJAN-MAR by 1.1 °C. Surprisingly, in the last 20 years, the precipitation amount did not decrease; in fact, there was a slight increase—SAP by 30 mm, SPVEG by 17 mm, and SPOCT-MAR by 14 mm—whereas MPJAN remained constant.
Based on SPEI values, we identified the years of severe and extreme drought for the period 1901–2017. Figure 5 presents the frequency of drought events in these sites. There were six extreme drought years (1958, 1968, 1986, 2000, 2003, and 2012) in this period, whereas the years of severe drought were a total of 32. Most of the drought events (13 out of 38 extreme and severe drought years) occurred in the period 2000–2017. Further, during this period, the highest number of consecutive years of drought was recorded (2006–2009 and 2011–2015). Other such drought periods were 1961–1963 and 1972–1975.
The most significant drought event occurred in 2000, when two negative peaks were observed. The 2000 drought had the highest intensity and longest duration over the last century. In that year, eight months of drought were observed, of which six extreme droughts included the growing season and fall, until the end of December. The same intensity of drought was observed in 1934 only.
Furthermore, the years characterized by a long duration of drought but lower intensity were 2003 (five months of drought, of which four were of extreme/severe drought), 2011 (seven months of drought, of which three were of extreme/severe drought), and 2012 (five months of drought, of which three were of extreme/severe drought).

3.3. Species Response to Drought

The analysis of variance revealed significant differences between species for all drought parameters in FGR 2, located on a low-productivity site (Table 4). The drought reaction parameters of the two species differed significantly depending on the extreme drought years in both FGRs. Significant differences between study sites regarding the species drought reaction were found for Norway spruce only. Within-species differences were non-significant for both species and FGR. The variation of drought parameters for all extreme drought years on species and FGR is presented in Figure 6.
In FGR 2, differences in species reaction were obtained for all extreme drought years, except the year 2000 (Table 5). In FGR 1, only the extreme droughts during 2000 and 2003 produced a significant interspecific variation for recovery, resilience, and relative resilience.
In FGR 1, Douglas fir exceeded Norway spruce regarding resistance capacity during extreme drought in 2000 and resistance, recovery, resilience, and relative resilience during extreme droughts in 2003 and 2012. In FGR 2, Douglas fir displayed a higher recovery during the extreme drought in 2003 and higher resistance, recovery, resilience, and relative resilience during the extreme drought in 2012. It should be noted that we used a three-year period before and after the extreme event to calculate drought parameters. For the year 2000, the analyzed period after the drought included the years 2002 (with moderate drought) and 2003 (with extreme drought).
Norway spruce was slightly better than Douglas fir in terms of recovery and resilience capacity during the extreme event in 2000 in both sites, and resistance and resilience in FGR 2 during the summer drought in 2003. Norway spruce also showed a higher capacity of adapting to the spring drought in 1968 and autumn drought in 1986 than Douglas fir.
Establishing the ranking of drought parameters in both study sites for all years of extreme drought, Norway spruce showed a higher resilience to extreme drought events, whereas Douglas fir showed higher resistance and relative resilience. Both tree species recorded similar recovery capacities (Figure 7). Correlations among growth traits and drought parameters demonstrated that Douglas fir had the highest resistance to drought and the widest growth rings (Table 6).

3.4. Growth Response Functions

The climate–growth response relationships varied between conifer species and study sites (Table 7). The main climate variables linked to tree radial growth of Douglas fir located on the high-productivity site (FGR 1) were MTOCT-MAR and SPVEG, whereas for FGR 2 on the low-productivity site, these variables were MTJUL, MTOCT-MAR, and SPOCT-MAR. The main climatic drivers explaining radial growth of Norway spruce were MTJUL and SPOCT-MAR in the first study site and MTVEG and SPJAN-MAR in the second study site.
Douglas fir appears to be less sensitive to precipitation than to temperature. The autumn and winter temperatures were positively correlated with Douglas fir growth, but the July temperatures negatively influenced both species. A negative relationship between radial growth and MTVEG for Norway spruce was also found. Increasing precipitation in the growing season and in the autumn–winter interval positively influenced Douglas fir radial growth, but this influence was negative on Norway spruce.
The growth–climate relationships were weak. R2 ranged between 0.10 and 0.12 for Douglas fir and between 0.06 and 0.10 for Norway spruce, which suggest that only a small portion of the radial growth variation can be explained by climatic factors.
Similarly, a modest response was found for latewood percentage. Thus, significant correlations between climate and LWP of Douglas fir were found in FGR 2, whereas for Norway spruce in FGR 1, MTJUL, MTVEG, and SPOCT-MAR negatively influenced the LWP.

3.5. Potential Impact of Climate Change Projections

The regions suitable for growing Douglas fir in Romania based on climate data over 1961–2020 are shown in Figure 8. For climatic envelope modeling, we used climatic variables, considering them to be the most important in the context of adaptation to climate change and because it is well known that Douglas fir grows on different soil types that are moderately acidic.
The projection suitability maps under RCP4.5 and the RCP8.5 scenarios over the periods 2041–2070 and 2071–2100 are presented in Figure 9. The maps show that the species suitability will increase along altitudes and in other Carpathians regions with moderate climates. Our results suggest that this species may be cultivated on a larger scale than has been considered so far.

4. Discussion

In this study, we analyzed the growth and adaptive capacity to climate change of Douglas fir compared to Norway spruce in two forest genetic resources, which are among the few oldest plantations with Douglas fir in Romania.
Substantial differences in growth traits and wood characteristics of the two stands were found. Douglas fir exceeded Norway spruce for all studied characteristics in both locations, but the best growth performances were obtained on the low-productivity site. The average volume per tree, ring width, and latewood percentage were approximately 76%, 41%, and 19%, respectively, higher for Douglas fir than for Norway spruce.
The analysis of climate data showed significant changes in climate conditions in the last century in western Romania, indicating that the growing conditions for forest tree species have changed. Significant variations occurred in the last two decades, when the mean annual temperature increased by 0.9 °C, mean temperature in January by 1.3 °C, and mean temperature during the growing period by 1.0 °C. It is surprising that the precipitation amount did not decrease; in fact, a slight increase was observed.
The growth response functions revealed that the climatic variables of the study sites were a major driver of the growth performance of the two conifer species, but the climate variables depend on species. There was a significant negative response to the temperature in July and a positive influence of previous autumn–winter temperature for Douglas fir. Provenance studies show that winter air temperatures are of the utmost importance for populations’ adaptation, limiting Douglas fir growth in Europe [15,21]. Our study confirms that autumn–winter temperatures are the most important factor in determining the tree-ring width in Douglas fir. Furthermore, increasing precipitation during the growing season or in the previous autumn–winter period had a positive effect on trees’ radial growth but negatively influenced latewood percentage. These growth–climate relationships are similar to those observed in a number of regional studies in North America [68,69,70].
For Norway spruce growing outside its natural range in the study locations, the July temperatures and temperatures during the growing season were the limiting factors in determining tree growth. The observed significant negative influence of SPOCT-MAR and SPIAN-MAR point to Norway spruce being less sensitive to precipitation during the period of vegetative rest compared to Douglas fir. The latewood percentage was sensitive to increasing the July and growing season temperatures and autumn–winter precipitation for both species.
The growth–climate relationships were, in general, weak. According to Fritts [52], the strength of the climate–growth relationship depends on how strongly growth is limited by the climate of the study site. At Ana Lugojana Forest District, the mean values of SAP for the 1901–2017 period was 754 mm, and MAT was 8.97 °C, which indicates that Norway spruce is placed in a climatic suboptimum and that decreased precipitation and increased temperatures can have negative effects on its growth and survival. The coastal Douglas fir is a variety adapted to highly different site conditions; in western Oregon and Washington, it occurs from sea level to over 1700 m. Its adaptation is a consequence of trade-offs between selection for traits to avoid exposure to frost and traits that confer high vigor in mild environments [71].
The climate of Ana Lugojana Forest District is at the lower limit of its climatic optimum concerning SAP, winter temperature, and frost-free days. This may explain the exceptional growth performance of Douglas fir in the two study sites from western Romania. The difference between the climate within the native distribution range and western Romania consists of the precipitation distribution over the year. Although precipitation is nearly evenly distributed over the year in Romania, precipitation at the place of origin has a pronounced maximum in winter and a minimum in summer [1].
Considerable differences in response to drought events were found between the two tree species. The species reaction depended on the timing and duration of the drought event. Generally, Douglas fir had higher resistance, recovery, resilience, and relative resilience to the summer droughts from 2000, 2003, and 2012 than Norway spruce, which showed a higher tolerance to spring drought in 1968 and autumn drought in 1986 than Douglas fir. Regarding the ranking of drought parameters for both sites and all extreme drought years, Norway spruce displayed a higher resilience to extreme droughts, whereas Douglas fir showed a higher resistance and relative resilience. Both species recorded a similar recovery capacity.
Norway spruce is particularly vulnerable to drought [72,73,74,75]. In contrast, Douglas fir is more drought resistant than Norway spruce [2,47,48]. Rising temperatures and decreasing precipitation in the near future may increase mortality risk for Norway spruce and other native coniferous species, particularly at low elevation sites [44,76,77,78]. Douglas fir is a fast-growing tree and can be a potential species for biomass production. In Europe, Douglas fir produces high-quality timber, which generally equals or exceeds the value of indigenous softwoods species. Thus, the Douglas fir dry wood density is on average 0.45 t/m−3 [2], while of Norway spruce is on average 0.40 t/m−3 (wood-database.com/Norway spruce). One of the reasons for this potential is that under the global warming scenario, the likelihood of cold damage to low-elevation sources will decrease with time. Recent studies have shown that planting Douglas fir with broadleaf species had a positive effect on the survival of this species [2,78]. In this context, the selection of suitable Douglas fir forest reproductive material is essential because it affects the growth, frost sensitivity, and tolerance to diseases [15,22,71,79].
Our results based on growth response functions and climate models by 2100 suggest that Douglas fir has a high growth potential in many Carpathians regions, not only in the west of Romania as has been considered so far. Our results confirm that climate change can increase Douglas fir productivity at higher elevations as a consequence of improving growth conditions. Further, the developed models used for the assessment of Douglas fir suitability under future climate in Romania showed good survivability.

5. Conclusions

Selecting suitable tree species and provenances adapted to ongoing climate change is of great significance in forest management. The cultivation of Douglas fir, a non-native tree species, has often been associated with higher risk and uncertainty. The developed models used for the assessment of Douglas fir suitability under the predicted future climate in Romania showed good survivability in many Carpathian regions. Douglas fir has an exceptional growth capacity, overcoming the Norway spruce in both high- and low-productivity sites, currently as well as in the past. Considerable differences were found in drought tolerance as well, with our results demonstrating the high resistance and relative resilience to extreme droughts of Douglas fir compared to Norway spruce. However, autumn–winter temperatures play an important role in the adaptation of Douglas fir to site conditions in Romania.
The use of appropriate provenances of Douglas fir in mixed stands with native broadleaved species may be an option for climatically exposed sites, thus increasing the value of these stands. Therefore, the conservation of the most valuable genetic resources of Douglas fir should have priority because these stands can be potential seed sources for ecosystem restoration. This study reveals the importance of improving our knowledge about the growth, ecology, and adaptive capacity of this non-native species in the context of climate change.

Author Contributions

Conceptualization: G.M., A.-L.C. and M.-V.B.; methodology: G.M. and A.-L.C., E.C. and M.-V.B.; field measurements and dendrochronological analyzes: A.-M.A.; climate data analysis and processing: I.-A.N. and M.-V.B.; maps: I.-A.N. All authors have read and agreed to the published version of the manuscript.

Funding

This study was carried out within the framework of the project PN 19070303 (revision of the provenance regions for the production and deployment of the forest reproductive materials in Romania in order to increase the adaptability of forest ecosystems to climate change). This work was supported by the Ministry of Higher Education, Research and Innovation in Romania, in BIOSERV Nucleu Program.

Data Availability Statement

The RoCliB climate data is publicly available.

Acknowledgments

We thank Alexandru Dumitrescu and Vlad Amihaesei (Meteo Romania) for providing the ensemble data of RoCliB. The CRU-TS dataset was made available by the Climatic Research Unit of the University of East Anglia.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Lavender, D.P.; Hermann, R.K.; OSU Forest Research Laboratory. Douglas-Fir: The Genus Pseudotsuga; OSU College of Forestry: Corvallis, OR, USA, 2014; ISBN 9780615979953. [Google Scholar]
  2. Spiecker, H.; Lindner, M.; Schuler, J. Douglas Fir—An Option for Europe; European Forest Institute: Joensuu, Finland, 2019. [Google Scholar]
  3. Ivanovici, D.C. Relativ La Bradul Duglas (Regarding the Douglas Fir). Rev. Padur. 1952, 12, 73. [Google Scholar]
  4. Lazarescu, C.; Ionescu, C. Cultura Duglasului Verde Si a Pinului Strob (Culture of Douglas Fir and Weymouth Pine); Agrosilvica: Bucharest, Romania, 1964. [Google Scholar]
  5. Lacaze, J.F.; Tomassone, R. Contribution a l’etude de la Variabilite du Douglas (Pseudotsuga menziesii Mirb.). Ann. For. Sci. 1967, 24, 85–106. [Google Scholar] [CrossRef] [Green Version]
  6. Sorensen, F. Provenance Variation in Pseudotsuga menziesii Seedlings from the Var. menziesii—Var. glauca Transition Zone in Oregon. Silvae Genet. 1979, 28, 96–103. [Google Scholar]
  7. Campbell, R.K. Mapped Genetic Variation of Douglas-Fir to Guide Seed Transfer in Southwest Oregon. Silvae Genet. 1986, 35, 85–96. [Google Scholar]
  8. Rehfeldt, G.E. Ecological Adaptations in Douglas-Fir (Pseudotsuga menziesii Var. glauca): A Synthesis. For. Ecol. Manag. 1989, 28, 203–215. [Google Scholar] [CrossRef] [Green Version]
  9. Silen, R.; Mandel, N. Clinal Genetic Growth Variation within Two Douglas Fir Breeding Zones. J. For. 1993, 81, 216–220. [Google Scholar]
  10. Konnert, M.; Alizoti, P.; Bastien, J. European Provenance Recommendations for Selected Non-Native Species 1—Introduction. In Proceedings of the WG2 Technical Report; Institute of Silviculture, University of Natural Resources and Life Sciences Vienna: Vienna, Austria, 2019; pp. 1–50. [Google Scholar]
  11. Alizoti, P.; Bastien, J.C.; Chakraborty, D.; Klisz, M.M.; Kroon, J.; Neophytou, C.; Schueler, S.; van Loo, M.; Westergren, M.; Konnert, M.; et al. Non-Native Forest Tree Species in Europe: The Question of Seed Origin in Afforestation. Forests 2022, 13, 273. [Google Scholar] [CrossRef]
  12. St Clair, J.B.; Mandel, N.L.; Vance-Borland, K.W. Genecology of Douglas Fir in Western Oregon and Washington. Ann. Bot. 2005, 96, 1199–1214. [Google Scholar] [CrossRef] [Green Version]
  13. Nepveu, G.; Blachon, J. Growth Ring Width and Suitability for Use as Construction Timber of Conifers: Douglas Fir, Scots Pine, Maritime Pine, Sitka Spruce, Norway Spruce and Silver Fir. Rev. For. Français 1989, 41, 497–506. [Google Scholar] [CrossRef] [Green Version]
  14. Vargas-Hernandez, J.; Adams, W.T. Genetic Variation of Wood Density Components in Young Coastal Douglas-Fir: Implications for Tree Breeding. Can. J. For. Res. 1991, 21, 1801–1807. [Google Scholar] [CrossRef]
  15. Bastien, J.C.; Sanchez, L.; Michaud, D. Douglas Fir (Pseudotsuga menziesii (Mirb.) Franco). In Forest Tree Breeding in Europe: Current State-of-the-Art; Pâques, L.E., Ed.; Managing Forest Ecosystems; Springer: Dordrecht, The Netherlands, 2013; Volume 25, pp. 325–373. ISBN 9789400761452. [Google Scholar]
  16. Vejpustková, M.; Čihák, T. Climate Response of Douglas Fir Reveals Recently Increased Sensitivity to Drought Stress in Central Europe. Forests 2019, 10, 97. [Google Scholar] [CrossRef] [Green Version]
  17. Schober, R. Experiences with the Douglas Fir in Europa; FAO/FORGEN: Stockholm, Sweeden, 1963. [Google Scholar]
  18. Silen, R. Genetics of Douglas Fir; U.S. Department of Agriculture: Washington, DC, USA, 1978.
  19. Burzyński, G. Resistance to the Frost Damages and Breeding Value of Different Douglas Fir Provenances. Ph.D. Thesis, Instytut Badawczy Leśnictwa, Sękocin Stary, Poland, 1999. [Google Scholar]
  20. Braun, H.; Wolf, H. Untersuchungen Zum Wachstum und Zur Frosthärte von Douglasien-Populationen in Ostdeutschland. Beiträge Forstwirtsch. Landsch. 2001, 35, 211–214. [Google Scholar]
  21. Breidenstein, I.; Bastien, J.C.; Roman-Amat, B. Douglas Fir Range Wide Variation. Results from the IUFRO Database. In Proceedings of the Joint Meeting of Western Forest Genetics Association and IUFRO Working Parties, Olympia, Greece, 20–24 August 1990. [Google Scholar]
  22. Kleinschmit, J.; Bastien, J.C. IUFRO’s Role in Douglas-Fir (Pseudotsuga menziesii (Mirb.) Franco) Tree Improvement. Silvae Genet. 1992, 4, 161–173. [Google Scholar]
  23. Enescu, V. Cercetări de Provenienţă la Duglas Şi Larice. Teste Timpurii (Researches on Douglas Fir and Larch Provenances. Early Tests), 2nd ed.; Series I; Institutul de Cercetari si Amenajări Silvice: Bucharest, Romania, 1984. [Google Scholar]
  24. Enescu, V. Ameliorarea Arborilor (Breeding of the Tree Species); CERES: Bucharest, Romania, 1973. [Google Scholar]
  25. Enescu, V.; Chereches, D.; Bandiu, C. Conservation of Biodiversity and of Genetic Resources; Agris: Bucharest, Romania, 1998. [Google Scholar]
  26. Meinshausen, M.; Smith, S.J.; Calvin, K.; Thomson, A.; Daniel, J.S.; Kainuma, M.L.T.; Matsumoto, K.; Lamarque, J.; Raper, S.C.B.; Riahi, K.; et al. The RCP Greenhouse Gas Concentrations and Their Extensions from 1765 to 2300. Clim. Chang. 2011, 109, 213–241. [Google Scholar] [CrossRef] [Green Version]
  27. Birsan, M.-V.; Micu, D.-M.; Niţă, I.-A.; Mateescu, E.; Szép, R.; Keresztesi, Á. Spatio-Temporal Changes in Annual Temperature Extremes over Romania (1961–2013). Rom. J. Phys. 2019, 64, 816. [Google Scholar]
  28. Birsan, M.V.; Marin, L.; Dumitrescu, A. Seasonal Changes in Wind Speed in Romania. Rom. Rep. Phys. 2013, 65, 1479–1484. [Google Scholar]
  29. Birsan, M.V.; Nita, I.-A.; Craciun, A.; Sfica, L.; Keresztesi, Á.; Szep, R.; Micheu, M. Observed Changes in Mean and Maximum Monthly Wind Speed over Romania since AD 1961. Rom. Rep. Phys. 2020, 72, 702. [Google Scholar]
  30. Micu, D.M.; Dumitrescu, A.; Cheval, S.; Nita, I.A.; Birsan, M.V. Temperature Changes and Elevation-Warming Relationships in the Carpathian Mountains. Int. J. Climatol. 2021, 41, 2154–2172. [Google Scholar] [CrossRef]
  31. Manea, A.; Birsan, M.V.; Tudorache, G.; Cărbunaru, F. Changes in the Type of Precipitation and Associated Cloud Types in Eastern Romania (1961–2008). Atmos. Res. 2016, 169, 357–365. [Google Scholar] [CrossRef]
  32. Bojariu, R.; Bîrsan, M.-V.; Cică, R.; Velea, L.; Burcea, S.; Dumitrescu, A.; Dascălu, S.I.; Gothard, M.; Dobrinescu, A.; Cărbunaru, F.; et al. Schimbările Climatice—De la Bazele Fizice la Riscuri Și Adaptare; Bojariu, R., Ed.; Printech: Bucharest, Romania, 2015; ISBN 9786062303631. [Google Scholar]
  33. Dumitrescu, A.; Birsan, M.V. ROCADA: A Gridded Daily Climatic Dataset over Romania (1961–2013) for Nine Meteorological Variables. Nat. Hazards 2015, 78, 1045–1063. [Google Scholar] [CrossRef]
  34. Cheval, S.; Busuioc, A.; Dumitrescu, A.; Birsan, M.V. Spatiotemporal Variability of Meteorological Drought in Romania Using the Standardized Precipitation Index (SPI). Clim. Res. 2014, 60, 235–248. [Google Scholar] [CrossRef]
  35. Cheval, S.; Dumitrescu, A.; Birsan, M.V. Variability of the Aridity in the South-Eastern Europe over 1961–2050. Catena 2017, 151, 74–86. [Google Scholar] [CrossRef]
  36. Nita, I.-A.; Sfica, L.; Apostol, L.; Radu, C.; Keresztesi, Á.; Szep, R. Changes in Cyclone Intensity over Romania According to 12 Tracking Methods. Rom. Rep. Phys. 2020, 72, 706. [Google Scholar]
  37. Nita, I.-A.; Apostol, L.; Patriche, C.; Sfica, L.; Bojariu, R.; Birsan, M.-V. Frequency of Atmospheric Circulation Types over Romania According to Jenkinson-Collison Method Based on Two Long-Term Reanalysis Datasets. Rom. J. Phys. 2022. In press. [Google Scholar]
  38. Țîmpu, S.; Sfîcă, L.; Dobri, R.-V.; Cazacu, M.-M.; Nita, A.-I.; Birsan, M.-V. Tropospheric Dust and Associated Atmospheric Circulations over the Mediterranean Region with Focus on Romania’s Territory. Atmosphere 2020, 11, 349. [Google Scholar] [CrossRef] [Green Version]
  39. Sfîcă, L.; Beck, C.; Nita, A.I.; Voiculescu, M.; Birsan, M.V.; Philipp, A. Cloud Cover Changes Driven by Atmospheric Circulation in Europe during the Last Decades. Int. J. Climatol. 2021, 41, E2211–E2230. [Google Scholar] [CrossRef]
  40. Ciais, P.; Sabine, C.; Bala, G.; Bopp, L.; Brovkin, V.; Canadell, J.; Chhabra, A.; DeFries, R.; Galloway, J.; Heimann, M.; et al. Carbon and Other Biogeochemical Cycles. In Climate Change 2013: The Physical Science Basis. Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change; Cambridge University Press: Cambridge, UK, 2014; pp. 465–570. [Google Scholar]
  41. Bolte, A.; Ammer, C.; Löf, M.; Madsen, P.; Nabuurs, G.J.; Schall, P.; Spathelf, P.; Rock, J. Adaptive Forest Management in Central Europe: Climate Change Impacts, Strategies and Integrative Concept. Scand. J. For. Res. 2009, 24, 473–482. [Google Scholar] [CrossRef]
  42. Lindner, M.; Maroschek, M.; Netherer, S.; Kremer, A.; Barbati, A.; Garcia-Gonzalo, J.; Seidl, R.; Delzon, S.; Corona, P.; Kolström, M.; et al. Climate Change Impacts, Adaptive Capacity, and Vulnerability of European Forest Ecosystems. For. Ecol. Manag. 2010, 259, 698–709. [Google Scholar] [CrossRef]
  43. Thom, D.; Seidl, R.; Steyrer, G.; Krehan, H.; Formayer, H. Slow and Fast Drivers of the Natural Disturbance Regime in Central European Forest Ecosystems. For. Ecol. Manag. 2013, 307, 293–302. [Google Scholar] [CrossRef]
  44. Cochard, H. Vulnerability of Several Conifers to Air Embolism. Tree Physiol. 1992, 11, 73–83. [Google Scholar] [CrossRef]
  45. Bouriaud, O.; Popa, I. Comparative Dendroclimatic Study of Scots Pine, Norway Spruce, and Silver Fir in the Vrancea Range, Eastern Carpathian Mountains. Trees—Struct. Funct. 2009, 23, 95–106. [Google Scholar] [CrossRef]
  46. Eilmann, B.; de Vries, S.M.G.; den Ouden, J.; Mohren, G.M.J.; Sauren, P.; Sass-Klaassen, U. Origin Matters! Difference in Drought Tolerance and Productivity of Coastal Douglas-Fir (Pseudotsuga menziesii (Mirb.)) Provenances. For. Ecol. Manag. 2013, 302, 133–143. [Google Scholar] [CrossRef]
  47. Podrázský, V. Potential of Douglas—Fir as a Partial Substitute for Norway Spruce—Review of the Newest Czech Literature. Beskydy 2015, 8, 55–58. [Google Scholar] [CrossRef]
  48. Vitali, V.; Büntgen, U.; Bauhus, J. Silver Fir and Douglas Fir Are More Tolerant to Extreme Droughts than Norway Spruce in South-Western Germany. Glob. Chang. Biol. 2017, 23, 5108–5119. [Google Scholar] [CrossRef] [PubMed]
  49. Parnauta, G.; Stuparu, E.; Budeanu, M. Catalogul National Al Resurselor Genetice Forestiere (National Catalog of Forest Genetics Resources); Editura Silvica: Bucharest, Romania, 2011. [Google Scholar]
  50. Stokes, M.A.; Smiley, T.L. An Introduction to Tree-Ring Dating; University of Arizona Press: Tucson, AZ, USA, 1996; ISBN 0816516804. [Google Scholar]
  51. Grissino-Mayer, H. Evaluating Crossdating Accuracy: A Manual and Tutorial for the Computer Program COFECHA. Tree-Ring Res. 2001, 57, 205–221. [Google Scholar]
  52. Frits, H.C. Tree Rings and Climate; Academic Press: Cambridge, MA, USA, 1976; ISBN 978-0-12-268450-0. [Google Scholar]
  53. Briffa, K.R.; Jones, P.D. Basic Chronology Statistics and Assesment; Springer: Dordrecht, The Netherlands, 1990; ISBN 978-90-481-4060-2. [Google Scholar]
  54. Bunn, A. A Dendrochronology Program Library in R (DplR). Dendrochronologia 2008, 26, 115–124. [Google Scholar] [CrossRef]
  55. 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] [Green Version]
  56. Dumitrescu, A.; Birsan, M.V.; Manea, A. Spatio-Temporal Interpolation of Sub-Daily (6 h) Precipitation over Romania for the Period 1975–2010. Int. J. Climatol. 2016, 36, 1331–1343. [Google Scholar] [CrossRef] [Green Version]
  57. Dumitrescu, A.; Birsan, M.-V.; Nita, I.-A. A Romanian Daily High-Resolution Gridded Dataset of Snow Depth (2005–2015). Geofizika 2017, 34, 2. [Google Scholar] [CrossRef]
  58. Dumitrescu, A.; Amihaesei, V.-A.; Cheval, S. RoCliB—bias-corrected CORDEX RCM dataset over Romania. Geosci. Data J. 2022. [Google Scholar] [CrossRef]
  59. Dumitrescu, A.; Amihaesei, V. RoCliB—Bias Corrected CORDEX RCM Dataset over Romania (2.0) [Data Set]. Available online: https://zenodo.org/record/6336837#.YoduUKhByUl (accessed on 1 April 2022).
  60. Kottek, M.; Grieser, J.; Beck, C.; Rudolf, B.; Rubel, F. World Map of the Köppen-Geiger Climate Classification Updated. Meteorol. Z. 2006, 15, 259–263. [Google Scholar] [CrossRef]
  61. 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] [Green Version]
  62. Thornthwaite, C.W. An Approach toward a Rational Classification of Climate. Soil Sci. 1948, 66, 77. [Google Scholar] [CrossRef] [Green Version]
  63. Beguería, S.; Vicente-Serrano, S.M. Calculation of the Standardised Precipitation-Evapotranspiration Index. R-Package. CRAN. 2017. Available online: https://cran.r-project.org/web/packages/SPEI/index.html (accessed on 1 April 2022).
  64. 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]
  65. Wang, T.; O’Neill, G.; Aitken, S.N. Integrating Environmental and Genetic Effects to Predict Responses of Tree Populations to Climate. Ecol. Appl. 2010, 20, 153–163. [Google Scholar] [CrossRef] [PubMed]
  66. Mihai, G.; Birsan, M.V.; Dumitrescu, A.; Alexandru, A.; Mirancea, I.; Ivanov, P.; Stuparu, E.; Teodosiu, M.; Daia, M. Adaptive Genetic Potential of European Silver Fir in Romania in the Context of Climate Change. Ann. For. Res. 2018, 61, 95–108. [Google Scholar] [CrossRef]
  67. Sofletea, N.; Curtu, L. Dendrologie; Editura Pentru Viata: Brasov, Romania, 2001. [Google Scholar]
  68. González-Elizondo, M.; Jurado, E.; Návar, J.; González-Elizondo, M.S.; Villanueva, J.; Aguirre, O.; Jiménez, J. Tree-Rings and Climate Relationships for Douglas-Fir Chronologies from the Sierra Madre Occidental, Mexico: A 1681–2001 Rain Reconstruction. For. Ecol. Manag. 2005, 213, 39–53. [Google Scholar] [CrossRef]
  69. Littel, J.; Peterson, D.; Tjoelker, M. Douglas Fir Growth in Mountain Ecosystems: Water Limits Tree Growth from Stand to Region. Ecol. Monogr. 2008, 28, 349–368. [Google Scholar] [CrossRef]
  70. Chen, P.Y.; Welsh, C.; Hamann, A. Geographic Variation in Growth Response of Douglas-Fir to Interannual Climate Variability and Projected Climate Change. Glob. Chang. Biol. 2010, 16, 3374–3385. [Google Scholar] [CrossRef]
  71. St. Clair, J.B. Genetic Variation in Fall Cold Hardiness in Coastal Douglas-Fir in Western Oregon and Washington. Can. J. Bot. 2006, 84, 1110–1121. [Google Scholar] [CrossRef]
  72. Pretzsch, H.; Dursky, J. Growth Reaction of Norway Spruce (Picea abies) and European Beechn (Fagus silvatica) to Possible Climatic Changes in Germany. A Sensitiy Study. Blackwell Verl. 2002, 121, 145–154. [Google Scholar]
  73. Pichler, P.; Oberhuber, W. Radial Growth Response of Coniferous Forest Trees in an Inner Alpine Environment to Heat-Wave in 2003. For. Ecol. Manag. 2007, 242, 688–699. [Google Scholar] [CrossRef]
  74. Boden, S.; Kahle, H.P.; von Wilpert, K.; Spiecker, H. Resilience of Norway Spruce (Picea abies (L.) Karst) Growth to Changing Climatic Conditions in Southwest Germany. For. Ecol. Manag. 2014, 315, 12–21. [Google Scholar] [CrossRef]
  75. Zang, C.; Hartl-Meier, C.; Dittmar, C.; Rothe, A.; Menzel, A. Patterns of Drought Tolerance in Major European Temperate Forest Trees: Climatic Drivers and Levels of Variability. Glob. Chang. Biol. 2014, 20, 3767–3779. [Google Scholar] [CrossRef]
  76. 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]
  77. Lebourgeois, F.; Rathgeber, C.B.K.; Ulrich, E. Sensitivity of French Temperate Coniferous Forests to Climate Variability and Extreme Events (Abies alba, Picea abies and Pinus sylvestris). J. Veg. Sci. 2010, 21, 364–376. [Google Scholar] [CrossRef]
  78. Brandl, S.; Paul, C.; Knoke, T.; Falk, W. The Influence of Climate and Management on Survival Probability for Germany’s Most Important Tree Species. For. Ecol. Manag. 2020, 458, 117652. [Google Scholar] [CrossRef]
  79. Konnert, M.; Fady, B.; Gömöry, D.; A’Hara, S.; Wolter, F.; Ducci, F.; Koskela, J.; Bozzano, M.; Maaten, T.; Kowalczyk, J. Use and Transfer of Forest Reproductive Material: In Europe in the Context of Climate Change; European Forest Genetic Resources Programme (EUFORGEN): Barcelona, Spain, 2015; 75p. [Google Scholar]
Figure 1. Total height (a) and volume/tree (b) of Douglas fir (DU) and Norway spruce (NS) in forest genetics resources (FGR) located in western Romania.
Figure 1. Total height (a) and volume/tree (b) of Douglas fir (DU) and Norway spruce (NS) in forest genetics resources (FGR) located in western Romania.
Forests 13 00805 g001
Figure 2. Variation of the ring width index (RWI) and latewood percentage (LWI) of Douglas fir (DU) and Norway spruce (NS) in FGR 1.
Figure 2. Variation of the ring width index (RWI) and latewood percentage (LWI) of Douglas fir (DU) and Norway spruce (NS) in FGR 1.
Forests 13 00805 g002
Figure 3. Variation of the ring width index (RWI) and latewood percentage (LWP) of Douglas fir (DU) and Norway spruce (NS) in FGR 2.
Figure 3. Variation of the ring width index (RWI) and latewood percentage (LWP) of Douglas fir (DU) and Norway spruce (NS) in FGR 2.
Forests 13 00805 g003
Figure 4. Variation of the mean annual temperature (MAT) and sum annual precipitation (SAP) in 1901–2017 period.
Figure 4. Variation of the mean annual temperature (MAT) and sum annual precipitation (SAP) in 1901–2017 period.
Forests 13 00805 g004
Figure 5. Variation of standardized precipitation evapotranspiration index (SPEI) across the analyzed period.
Figure 5. Variation of standardized precipitation evapotranspiration index (SPEI) across the analyzed period.
Forests 13 00805 g005
Figure 6. Variation of drought parameters for Douglas fir (DU) and Norway spruce (NS) in forest genetic resources (FGR 1 and FGR 2).
Figure 6. Variation of drought parameters for Douglas fir (DU) and Norway spruce (NS) in forest genetic resources (FGR 1 and FGR 2).
Forests 13 00805 g006
Figure 7. The ranking of drought parameters for Douglas fir (DU) and Norway spruce (NS) in forest genetic resources (FGR 1 and FGR 2).
Figure 7. The ranking of drought parameters for Douglas fir (DU) and Norway spruce (NS) in forest genetic resources (FGR 1 and FGR 2).
Forests 13 00805 g007
Figure 8. Suitability areas for the Douglas fir in Romania according to the high-resolution climatological data over 1961–2020. The polygons denote the regions of provenance (as in [49]): A1: Eastern Carpathians, western cline; A2: Eastern Carpathians, eastern cline; A3: Giurgeu–Ciuc depression; B1: Brasov depression; B2: Curvature Carpathians, outer cline; C1: Southern Carpathians, northern cline; C2: Southern Carpathians, southern cline; D1: Mehedinti/Cerna/Semenic Mountains; D2: Tarcu/Poiana Rusca Mountains; E1: Zarand/Metaliferi Mountains; E2: Western Apuseni Mountains; E3: Eastern Apuseni Mountains; F1: Transylvania Plain; F2: Transylvania Plateau; G1: Suceva/Siret/Iasi Hills; G2: Jijia Plain; G3: Barlad Plateau; H1: Covur Plateau; H2: Siret and Baragan Plains; H3: Danube water holes; I1: Danube Delta; I2: Dobrogea Plateau; J1: Bucharest Plain; J2: Oltenia Plain; K1: Timis and Arad Plain; K2: Cris/Carei/Somes Plain.
Figure 8. Suitability areas for the Douglas fir in Romania according to the high-resolution climatological data over 1961–2020. The polygons denote the regions of provenance (as in [49]): A1: Eastern Carpathians, western cline; A2: Eastern Carpathians, eastern cline; A3: Giurgeu–Ciuc depression; B1: Brasov depression; B2: Curvature Carpathians, outer cline; C1: Southern Carpathians, northern cline; C2: Southern Carpathians, southern cline; D1: Mehedinti/Cerna/Semenic Mountains; D2: Tarcu/Poiana Rusca Mountains; E1: Zarand/Metaliferi Mountains; E2: Western Apuseni Mountains; E3: Eastern Apuseni Mountains; F1: Transylvania Plain; F2: Transylvania Plateau; G1: Suceva/Siret/Iasi Hills; G2: Jijia Plain; G3: Barlad Plateau; H1: Covur Plateau; H2: Siret and Baragan Plains; H3: Danube water holes; I1: Danube Delta; I2: Dobrogea Plateau; J1: Bucharest Plain; J2: Oltenia Plain; K1: Timis and Arad Plain; K2: Cris/Carei/Somes Plain.
Forests 13 00805 g008
Figure 9. Projected suitability areas for the Douglas fir in Romania under the RCP4.5 (a,b) and RCP8.5 (c,d) scenarios, over the periods 2041–2070 (a,c) and 2071–2100 (b,d). The polygons denote various regions of provenance (as defined by [49]). A1: Eastern Carpathians, western cline; A2: Eastern Carpathians, eastern cline; A3: Giurgeu–Ciuc depression; B1: Brasov depression; B2: Curvature Carpathians, outer cline; C1: Southern Carpathians, northern cline; C2: Southern Carpathians, southern cline; D1: Mehedinti/Cerna/Semenic Mountains; D2: Tarcu/Poiana Rusca Mountains; E1: Zarand/Metaliferi Mountains; E2: Western Apuseni Mountains; E3: Eastern Apuseni Mountains; F1: Transylvania Plain; F2: Transylvania Plateau; G1: Suceva/Siret/Iasi Hills; G2: Jijia Plain; G3: Barlad Plateau; H1: Covur Plateau; H2: Siret and Baragan Plains; H3: Danube water holes; I1: Danube Delta; I2: Dobrogea Plateau; J1: Bucharest Plain; J2: Oltenia Plain; K1: Timis and Arad Plain; K2: Cris/Carei/Somes Plain.
Figure 9. Projected suitability areas for the Douglas fir in Romania under the RCP4.5 (a,b) and RCP8.5 (c,d) scenarios, over the periods 2041–2070 (a,c) and 2071–2100 (b,d). The polygons denote various regions of provenance (as defined by [49]). A1: Eastern Carpathians, western cline; A2: Eastern Carpathians, eastern cline; A3: Giurgeu–Ciuc depression; B1: Brasov depression; B2: Curvature Carpathians, outer cline; C1: Southern Carpathians, northern cline; C2: Southern Carpathians, southern cline; D1: Mehedinti/Cerna/Semenic Mountains; D2: Tarcu/Poiana Rusca Mountains; E1: Zarand/Metaliferi Mountains; E2: Western Apuseni Mountains; E3: Eastern Apuseni Mountains; F1: Transylvania Plain; F2: Transylvania Plateau; G1: Suceva/Siret/Iasi Hills; G2: Jijia Plain; G3: Barlad Plateau; H1: Covur Plateau; H2: Siret and Baragan Plains; H3: Danube water holes; I1: Danube Delta; I2: Dobrogea Plateau; J1: Bucharest Plain; J2: Oltenia Plain; K1: Timis and Arad Plain; K2: Cris/Carei/Somes Plain.
Forests 13 00805 g009
Table 1. The sites and stands characteristics of the Douglas fir forest genetic resources.
Table 1. The sites and stands characteristics of the Douglas fir forest genetic resources.
Stand CharacteristicsFGR 1FGR 2
Species composition80% Douglas fir; 20% Norway spruce30% European beech; 50% Douglas fir; 10% Norway spruce; 10% other species
Age105130
Class of productionII
SoilEutricambosoil mollicEutricambosoil lithic
Site conditionsHigh productivity Low productivity
Vegetation layerEuropean beech layerEuropean beech layer
Altitude880 m450 m
Climatic province
(by Köppen–Geiger)
Cfb-warm and temperate climate with slight Mediterranean influencesCfb-warm and temperate climate with slight Mediterranean influences
Table 2. Comparative climatic data for Douglas fir forest genetic resources from western Romania and at the natural distribution range in Pacific Northwest.
Table 2. Comparative climatic data for Douglas fir forest genetic resources from western Romania and at the natural distribution range in Pacific Northwest.
SiteMAT
(°C)
MTJAN
(°C)
MTJAN–MAR
(°C)
MTVEG
(°C)
SAP
(mm)
MPJAN
(mm)
SPJAN–MAR
(mm)
SPVEG
(mm)
Frost-Free Days
Western Romania8.97−2.560.0115.6275443183477248
Pacific Northwest3.5 to 14.4–2.0 to 3.0–1.5 to 57.4 to 23760–300015 to 52443 to 123390 to 750195–260
MAT—the mean annual temperature; MTJAN—the mean temperature of the coldest month (January); MTJAN–MAR—the mean temperature from January to March of the current year; MTVEG—the mean temperature of the growing season; SAP—the sum annual precipitation; SPJAN—the sum precipitation of the coldest month (January); SPVEG—the sum precipitation during the growing season.
Table 3. Analysis of variance of growth and wood traits for the studied FGRs.
Table 3. Analysis of variance of growth and wood traits for the studied FGRs.
Source of VariationDFVariance (s2)
Total
Height (m)
Diam.
1.30 m (cm)
Volume/
Tree (m3)
Ring
Width (mm)
Latewood
%
FGR 1Douglas fir (DU)185.917 ***102.894 ***9.030 ***34.859 ***2104.688 ***
Norway spruce (NS)159.772 ***29.857 ***0.622 ***13.413 ***2489.699 ***
Species (S)11282.276 ***12613.321 ***694.107 ***841.898 ***5330.840 ***
Year (Y)79---1.322 ***388.354 ***
Interaction S x Y79---0.618 ***297.571 ***
Error26407.66969.6955.2080.080164.581
FGR 2Douglas fir (DU)1814.618 ***238.937 ***24.728 ***27.718 ***1592.266 ***
Norway spruce (NS)157.947 ***52.250 ***1.007 ***13.223***2137.951 ***
Species (S)11686.527 ***15754.121 ***1260.260 ***451.099228.388 ***
Year (Y)79---2.371 ***375.961 ***
Interaction S x Y79---0.543 ***387.809 ***
Error256012.160181.35819.5890.098182.471
FGR 1 and FGR 2Douglas fir (DU)14.995974.623 *127.097 *383.776 ***55.155 ***
Norway spruce (NS)112.005220.500 *1.620121.424 ***30.932 ***
Species (S)12953.924 ***28115.320 ***1898.461 ***121.424 ***149.833
Year (Y) for DU 79---1.144 ***335.557
Interaction DU x Y79---0.721 ***311.100
Year (Y) for NS79---2.429 ***408.470 ***
Interaction NS x Y79---0.495 ***384.963 ***
Year (Y)79---2.911 ***424.742 ***
Interaction S x Y79---0.742 ***320.916 ***
Error53609.839138.78413.8420.092208.275
The level of significance is represented as follows: * p < 0.05; ** p < 0.01; *** p < 0.001.
Table 4. Analysis of variance of Douglas fir and Norway spruce drought parameters for years of extreme drought.
Table 4. Analysis of variance of Douglas fir and Norway spruce drought parameters for years of extreme drought.
TrialSource of
Variation
DFVariance (s2)
ResistanceRecoveryResilienceRelative
Resilience
FGR 1Douglas fir (DU)180.0640.2520.1060.152
Norway spruce (NS)150.0840.0830.3600.131
Between species10.0270.0560.0010.035
Extreme drought
year for DU
50.658 ***3.130 ***1.289 ***1.706 ***
Extreme drought
year for NS
50.1201.652 ***1.414 **1.208 ***
FGR 2Douglas fir (DU)180.0450.2410.0820.132
Norway spruce (NS)150.0900.1160.1180.089
Between species10.595 **1.454 *0.1470.156
Extreme drought
year for DU
50.638 ***3.303 ***0.865 ***1.947 ***
Extreme drought
year for NS
50.305 **1.109 ***1.977 ***1.187 ***
FGR 1 and FGR 2Douglas fir (DU)10.0760.0290.4420.143
Norway spruce (NS)10.423 *1.221 **0.0670.835 *
Between Species10.1831.055 *0.0720.025
Between sites for DU10.0760.0290.4420.143
Between sites for NS10.423 *1.221 **0.0670.835 *
Extreme drought
year for DU
50.420 ***5.353 ***1.589 ***3.282 ***
Extreme drought
year for NS
50.164 *2.669 ***2.688 ***2.278 ***
The level of significance is represented as follows: * p < 0.05; ** p < 0.01; *** p < 0.001.
Table 5. Differences between Douglas fir and Norway spruce for the drought parameters during the years of extreme drought.
Table 5. Differences between Douglas fir and Norway spruce for the drought parameters during the years of extreme drought.
FGRDrought
Parameters
Variance (s2)
195819681986200020032012
FGR 1Resistance0.0360.198 *0.0420.2840.0440.014
Recovery0.0200.1200.0081.373 *0.840 *0.510
Resilience0.0260.3890.0241.1051.355 *0.224
Rel. resilience0.0010.0340.0022.504 *0.893 *0.122
FGR 2Resistance0.552 *0.2840.1070.0010.972 ***0.004
Recovery1.791 *0.782 **0.112 *0.2383.773 **1.903 ***
Resilience0.0412.609 **0.0170.2240.0051.326 **
Rel. resilience0.893 *1.159 **0.209 *0.2040.845 **1.183 **
The level of significance is represented as follows: * p < 0.05; ** p < 0.01; *** p < 0.001.
Table 6. Phenotypic correlations between wood characters and drought parameters of Douglas fir (DU) and Norway spruce (NS).
Table 6. Phenotypic correlations between wood characters and drought parameters of Douglas fir (DU) and Norway spruce (NS).
SpeciesTraitRWLWEW
DUResistance0.594*0.769 **0.524
Recovery−0.888 ***−0.692 *−0.643 *
Resilience−0.476−0.105−0.224
Rel. resilience−0.811 **−0.594 *−0.580*
NSResistance−0.046−0.2690.151
Recovery−0.790 **−0.173−0.476
Resilience−0.776 **−0.363−0.462
Rel. resilience−0.797 **−0.159−0.483
RW—ring width; LW—late wood width; EW—early wood width; The level of significance is represented as follows: * p < 0.05; ** p < 0.01; *** p < 0.001.
Table 7. Climatic response models for radial growth (RW) and latewood percentage (LWP) of Douglas fir and Norway spruce.
Table 7. Climatic response models for radial growth (RW) and latewood percentage (LWP) of Douglas fir and Norway spruce.
SpeciesFGRGrowth Response ModelSignf.R2
Douglas firFGR 1RW = 0.864 + 0.008 MT2OCT-MAR + 0.001 SPVEG*0.12
Douglas firFGR 2RW = 1.686 − 0.002 MT2JUL + 0.005 MT2OCT-MAR + 0.001 SPOCT-MAR*0.10
Norway spruceFGR 1RW = 1.527 − 0.001 MT2JUL −0.0001 SP2OCT-MAR*0.06
Norway spruceFGR 2RW = 2.813 − 0.003 MT2VEG − 0.001 SP2IAN-MAR*0.10
Douglas firFGR 2LWP = 74.941 − 0.044 TM2JUL − 0.001 SP2OCT-MAR***0.23
Norway spruceFGR 1LWP = 56.419 − 0.038 MT2VEG − 0.014 SPOCT-MAR*0.07
MTJUL—the mean temperature for July; MTVEG—the mean temperature of the growing season; MTOCT-MAR—the mean temperature from October of the previous year to March of the current year; SPOCT-MAR—the sum precipitation from October of the previous year to March of the current year; SPJAN-MAR—the sum precipitation from January to March of the current year. The level of significance is represented as follows: * p < 0.05; ** p < 0.01; *** p < 0.001.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Mihai, G.; Curtu, A.-L.; Alexandru, A.-M.; Nita, I.-A.; Ciocîrlan, E.; Birsan, M.-V. Growth and Adaptive Capacity of Douglas Fir Genetic Resources from Western Romania under Climate Change. Forests 2022, 13, 805. https://0-doi-org.brum.beds.ac.uk/10.3390/f13050805

AMA Style

Mihai G, Curtu A-L, Alexandru A-M, Nita I-A, Ciocîrlan E, Birsan M-V. Growth and Adaptive Capacity of Douglas Fir Genetic Resources from Western Romania under Climate Change. Forests. 2022; 13(5):805. https://0-doi-org.brum.beds.ac.uk/10.3390/f13050805

Chicago/Turabian Style

Mihai, Georgeta, Alexandru-Lucian Curtu, Alin-Madalin Alexandru, Ion-Andrei Nita, Elena Ciocîrlan, and Marius-Victor Birsan. 2022. "Growth and Adaptive Capacity of Douglas Fir Genetic Resources from Western Romania under Climate Change" Forests 13, no. 5: 805. https://0-doi-org.brum.beds.ac.uk/10.3390/f13050805

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