Next Article in Journal
Nitrogen and Phosphorus of Plants Associated with Arbuscular and Ectomycorrhizas Are Differentially Influenced by Drought
Next Article in Special Issue
Metal- and Organ-Specific Response to Heavy Metal-Induced Stress Mediated by Antioxidant Enzymes’ Activities, Polyamines, and Plant Hormones Levels in Populus deltoides
Previous Article in Journal
Comparative Metabolomics Analysis of Stigmas and Petals in Chinese Saffron (Crocus sativus) by Widely Targeted Metabolomics
Previous Article in Special Issue
Visible Foliar Injury and Ecophysiological Responses to Ozone and Drought in Oak Seedlings
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Do Different Tree-Ring Proxies Contain Different Temperature Signals? A Case Study of Norway Spruce (Picea abies (L.) Karst) in the Eastern Carpathians

by
Andrei Popa
1,2,
Ionel Popa
1,3,*,
Cătălin-Constantin Roibu
4 and
Ovidiu Nicolae Badea
1,2
1
National Institute for Research and Development in Forestry ‘Marin Drăcea’, 007190 Bucharest, Romania
2
Faculty of Silviculture and Forest Engineering, Transilvania University of Brașov, 500036 Brașov, Romania
3
Center for Mountain Economy (CE–MONT), 725700 Vatra Dornei, Romania
4
Forest Biometrics Laboratory, Faculty of Forestry, ‘Stefan cel Mare’ University of Suceava, 720229 Suceava, Romania
*
Author to whom correspondence should be addressed.
Plants 2022, 11(18), 2428; https://doi.org/10.3390/plants11182428
Submission received: 9 August 2022 / Revised: 9 September 2022 / Accepted: 13 September 2022 / Published: 17 September 2022

Abstract

:
One of the most important proxy archives for past climate variation is tree rings. Tree-ring parameters offer valuable knowledge regarding how trees respond and adapt to environmental changes. Trees encode all environmental changes in different tree-ring parameters. In this study, we analyzed how air temperature is encoded in different Norway spruce tree-ring proxies along an altitude gradient in an intramountain valley of the Carpathians. The study area, in the Gheorgheni region, Romania (Eastern Carpathians), has a mountain climate with a frequent temperature inversion in winter. The climate–growth relationship was analyzed for two contrasting altitudes: low elevation, i.e., below 1000 m a.s.l., and high elevation, i.e., above 1500 m a.s.l. Two local weather stations, one in the valley and the other on the upper part of the mountains, provide daily temperatures (Joseni—750 m a.s.l. and Bucin—1282 m a.s.l.). The bootstrap Pearson correlation between cumulative daily temperature data and three tree-ring proxies (tree-ring width—TRW, basal area increment—BAI, and blue intensity—BI) was computed for each series. The results show that elevation modulates the climate response pattern in the case of BI, and remains relatively similar for TRW and BAI. The winter temperature’s positive influence on spruce growth was observed in both TRW and BAI chronologies. Additionally, the BAI chronology highlights a positive relationship with summer temperature. The highest correlation coefficient (r = 0.551, p < 0.05, n = 41) was recorded between BI residual chronology from high elevation series and summer/autumn temperature from the upper-part weather station for a cumulative period of 59 days (the second half of August to the beginning of October). Our results show that, for this intramountain valley of the Eastern Carpathians, different tree-ring proxies capture different climate signals.

Graphical Abstract

1. Introduction

Tree growth is driven by multiple factors, and climatic conditions represent one of the main drivers of wood accumulation [1,2]. Changes in climatic conditions (e.g., temperature rise, precipitation decrease) and higher frequencies of extreme events (e.g., heat waves, long periods of drought) induce stress on trees and forest stands, with negative ecological and economic effects [3,4]. All these environmental changes are encoded in different tree-ring parameters.
In the context of climate change, increasing efforts have been made to understand and assess the effect of environmental change on forest ecosystems [5]. Knowledge and mitigation of climate-change impacts represent a central goal for forest management systems [6,7,8,9,10].
Dendrochronological studies analyzing how trees are influenced by climate and how tree species adapt to new climate conditions also provide a record of the past climate [1,2,11]. Tree rings can be used as an important proxy to highlight annual climate variations [12]. The main tree-ring proxy used in dendrochronology is tree-ring width (TRW). However, in some cases, TRW does not provide a strong and robust climate signal compared to other tree-ring parameters such as maximum latewood density (MXD) or stable isotopes [13,14,15,16,17,18]. Usually, to ascertain MXD, expensive equipment is required, and there are logistical limitations. A relatively new parameter (image-based blue reflectance—blue intensity, BI) has been developed to respond to these limitations [19,20,21]. BI is a proxy that represents measured reflected light in specific wavelengths of the color spectrum. Studies have shown a strong correlation (over r = 0.95, p < 0.05) between BI and MXD [19,20,21,22,23,24]. BI has a stronger climate signal in temperature-limited ecosystems compared to TRW and is less sensitive to disturbances [24,25]. Based on these findings, there is potential for BI to be used as a substitute for MXD. However, the basal area increment (BAI) represents a two-dimensional measurement, specifically on the surface of the tree ring. The basal area increment is more related to the biomass increment of the tree and stand productivity [26,27]. Moreover, BAI is a suitable proxy that can preserve low and medium-frequency growth variability [28,29].
Norway spruce (Picea abies (L.) Karst) is one of the dominant coniferous species at both the Romanian and European levels [30]. On a large scale, Norway spruce is managed in monocultures due to its high productivity and ability to grow at high rates inside and outside of natural areas, with significant economic advantages [31]. Monoculture management, which focuses mainly on wood production, is prone to more issues compared to mixed-forest management, which focuses on both productivity and biodiversity, with higher stand resistance and resilience [8,32,33]. Generally, Norway spruce growth is driven by summer temperature in mountainous areas and by precipitation at lower elevations [34,35,36,37,38,39]. Spruce has demonstrated high sensitivity to extreme events such as drought or heatwaves [40,41,42].
In mountainous areas, due to alternating slopes and valleys, it is possible to have a regional climate modulated by a local topography [43]. The particular local climate specific to the intramountain valleys of the Carpathians can induce different climate responses of tree species (in our case, Norway spruce) compared with general patterns observed in mountainous regions.
In this study, we aimed to determine the climate signals captured in three tree proxies (TRW, BAI, and BI) of Norway spruce in an intramountain valley in the Eastern Carpathians along an altitudinal gradient. The specific research questions were:
  • How does air temperature modulate Norway spruce growth in an intramountain valley of the Carpathians?
  • Is the correlation between temperature and the investigated tree-ring parameters stable through time?

2. Results and Discussion

2.1. Description of Chronologies

The individual age of sampled trees varies from 68 years to 135 years with a mean chronologies age of 98 years, with insignificant differences between low- and high-elevation chronologies (Table 1). The average growth rate is 1.96 mm·year−1, and ranges from 1.77 mm·year−1 (high elevation) to 2.16 mm·year−1 (low elevation). A reduction in tree growth with increasing elevation is expected, being the consequence of temperature decrease and shortening of the growing season [34,44].
Common variance, expressed by the mean series intercorrelation (Rbar) with values around 0.3 for all tree-ring proxies except for BI at the low elevation, reflects a medium climatic control of Norway spruce growth in the study area. The mean sensitivity of the analyzed chronologies showed lower interannual variability for BI chronologies compared to TRW or BAI chronologies. Low mean-sensitivity values are typical for Norway spruce growing in optimal climatic conditions [35]. Lower values for mean sensitivity of BI chronologies have also been recorded for Norway spruce in other regions [45,46] or for other species from the Picea, Abies genus [45,47,48]. The temporal memory, expressed by serial autocorrelation of raw data, is highest in the case of TRW and lowest for BI.
The TRW, BAI, and BI index chronologies were truncated for the period 1978–2019 to overlap with the climatic data (Figure 1). The expressed population signal (EPS) for the analyzed period exceeds the threshold of 0.85 [49] for all chronologies.

2.2. Climate–Growth Relationships for Three Tree-Ring Parameters

The TRW residual index chronologies correlate positively with winter temperatures (cumulative windows width starting from 21 to 120 days) (Figure 2). The correlation coefficient between high-elevation TRW index chronology and mean temperature from the up-hill weather station (Bucin) has the highest value (r = 0.494, p < 0.05, n = 41) with the 3 December–18 January period. The low-elevation TRW index chronology has the highest correlation (r = 0.485, p < 0.05, n = 41) with the 1 November–12 February mean temperature from the up-hill weather station. Regarding the correlation between the TRW residuals index and the winter air temperature from the valley weather station, the maximum correlations were lower (r= 0.442, p < 0.05, n = 41—high-elevation chronology and r = 0.435, p < 0.05, n = 41—low-elevation chronology).
Moreover, the TRW residual index chronology from low-elevation sites shows a positive and significant correlation with spring temperature (recorded in the valley) from March to April (r = 0.374, p < 0.05, n = 41). For both chronologies, low and high elevation, a negative correlation between TRW residual index and mean temperatures is present in the previous vegetation season, in July. The maximum negative correlation between TRW residual index chronologies and previous summer temperature varies from r = −0.502 (p < 0.05, n = 41) for low-elevation series (temperature from the up-hill weather station) to r = −0.462 (p < 0.05, n = 41) for high-elevation series (temperature from the up-hill weather station). A significant negative correlation between TRW index chronology and temperature in the previous autumn is present only at low-elevation sites, regardless of the weather station.
A positive correlation between the TRW index and December temperatures has also been reported for other forests in the Eastern Carpathians [50]. In mountainous regions, and mainly at high elevation, the growth of Norway spruce is usually positively correlated with summer temperatures [35,37,38,50,51,52]. The positive correlation between TRW index chronologies and winter temperatures at an elevation above 1500 m a.s.l. is not a common dendroclimatic pattern for Norway spruce. The possible explanation for this climate–growth relationship could be related to temperature inversion with a high frequency during winter.
Generally, the air temperature is characterized by a decreasing gradient proportional to the elevation [53]. The average gradient for decreasing air temperature with the increase in elevation, in mountainous areas, is 0.6 °C/100 m [54]. Temperature inversions are phenomena caused by the abnormal variation of the radiative heat balance induced by terrain fragmentation, depth of valleys, and local topographical peculiarities [55,56,57]. Temperature inversions are a common occurrence in intramountain valleys, especially in the Eastern Carpathians [56,58]. Thermal inversion leads to cold-air stagnation at the bottom of the valleys, which favors fog formation with a significant impact on both human health and vegetation growth [59]. At the same time, the upper part of the slopes benefits from higher global radiation and warmer air mass with positive effects on vegetation and soil conditions.
Due to the temperature inversion phenomenon or mild winters, a higher air temperature in the cold season can lead to faster snow melt and higher soil temperature. These factors have a positive influence on cambium reactivation and apical growth in the spring, decreasing the risk of tissue freezing and even creating a low probability of xylem embolism in the spring [60,61,62]. These hypotheses could offer a possible explanation for the positive and significant correlation between temperatures in the winter and TRW index chronologies.
A negative correlation between TRW index series and temperature from the previous vegetation season has been observed in other Eastern Carpathian sites [50], Tatra Mountains [38,63], in the timberline in East-Central Europe [37,64], and in the Alps [65]. This correlation highlights a sensitivity of Norway spruce growth to the previous year’s temperature. Increasing summer temperatures may induce a reduction in the growth of spruce in the next year. This temporal memory of growth can be linked with structural carbohydrate dynamics and extending the growing season in previous years [51], and is supported by high values of first order autocorrelation (Table 1). The correlation of growth with the previous year temperature highlighted the importance of the carryover effects of climate variability, such as photosynthetic gain and storage of assimilates [1,66]. Moreover, a warm summer promotes the flowering of spruce, which is associated with a decrease in growth in the next year [67]. Hansen and Beck [68] highlight that the carbohydrates accumulated in previous autumn are depleted in spring. Also, the dynamics of spruce reserves involve an accumulation of lipids in summer which are metabolized during the autumn [69].
The BAI residual index chronology correlation pattern differs depending on elevation (Figure 3). The correlation between high-elevation BAI index chronology and winter temperatures recorded at the up-hill weather station is significant and has a higher value for the period of 26 November to 18 February. Interestingly, the positive correlation between winter–spring temperatures and the BAI index chronology from a low elevation is no longer significant when the temperatures are considered from the up-hill weather station. Regarding the high-elevation BAI index chronology, the correlation patterns with respect to the mean temperatures from both weather stations are similar. In the case of the winter period, the correlation with temperatures from the valley weather station is significant for longer cumulative window lengths.
Compared to the dendroclimatic pattern of TRW, a positive relationship with summer temperature was observed for BAI. The highest correlation was recorded between the BAI chronology from the high elevation and cumulative temperatures from the valley weather station from 17 May to 22 July (r = 0.468, p < 0.05, n = 41). A positive correlation with summer temperature on the current year’s wood accumulation is logical for high-altitude chronology, since temperature is a limiting factor for these habitats [35]. A negative correlation (r = −0.507, p < 0.05, n = 41) between BAI residual chronology from lower sites and previous summer temperatures is significant for cumulative periods of 29 days (second half of June to the middle of July) for both weather stations. The authors of [51] point out that higher temperatures during the summer can induce a high rate of respiration with negative effects on carbohydrate reserves used in the first phase of growth of the next year. The photosynthetic gain during the previous summer has a strong effect on current year ring width [67,70].
A clear pattern of positive and significant correlation coefficient between temperature and BI chronologies was found only for the high-elevation chronology (Figure 4). The highest correlation between the BI residual chronology from high-elevation series and summer/autumn from the up-hill station is 0.551, p < 0.05, n = 41 for cumulative windows of 59 days (second half of August to the beginning of October). An unusual correlation was found between the BI index and winter temperature from the valley weather station. This may be a false-positive correlation because it is less likely that winter temperature has a strong influence on the wood density of latewood. The correlation between low-elevation BI chronology and mean temperature from the previous summer is negative and significant (r = −0.481, p < 0.05, n = 41) for a cumulative window of 28 days; that is, 19 June–17 July. No significant correlation between BI index and previous year temperature was noted in the case of the high-elevation chronology.
In contrast, at low elevation, the BI chronology shows almost no significant correlation with current year temperatures. This suggests that at elevations below 1000 m, in this intramountain valley, the late-summer temperature is not a limiting factor in the thickening and lignification of cell walls of Norway spruce. This can be linked to a longer growing season at lower-altitude sites and a higher stand productivity [35]. The highest correlation between tree-ring parameters and temperature has been reported for BI chronologies from the high elevation. These BI correlation patterns show that thickening of the secondary cellular wall and the lignification process at high altitude are driven by the late-summer temperature [71]. The negative relationship with the previous year’s temperature can be explained by the trade-offs in the carbohydrate allocation for seed production, increment and formation of buds, with significant effects on nutrients reserves available for next year’s growing start [72,73]. It has already been reported in the literature that BI chronologies, as surrogates for maximum latewood density, express a stronger relationship with climate compared to TRW at sites where the temperature is the most limiting factor [16,20,74,75,76].
Tree-ring width or basal area increment are tree-ring proxies containing aggregated information about the climate conditions throughout the whole growth season, and about disturbances. Meanwhile, the blue intensity contains information about the second part of the growing season. The combination of these three tree-ring proxies can offer an integrated perspective on the climate–growth relationship.

2.3. Time Stability of Climate–Growth Relationship

The time stability correlation between TRW, BAI, and BI residual chronologies and temperature from the up-hill weather station was assessed using the cumulative window with the highest correlation coefficient. It was preferred to assess the time stability for windows with the highest correlations despite different periods between low and high elevation chronologies. The periods with the highest correlations may be the most relevant for assessing the climate sensitivity of Norway spruce. In the case of TRW, for the low-elevation series, the highest correlation was for the interval between 31 October and 18 February and, in the case of the high-elevation series, for the period between 5 December and 14 January. The time variation pattern of the correlation has evident temporal shifts, mainly in the second part of the analyzed period (Figure 5A). For the high-elevation chronology, the correlation has an increasing general trend with a small decrease around 1992 and stabilization after 2000. Contrary to the low-elevation chronology, the correlation intensity increased until the 1989–2003 period. After that, a slight decrease was recorded, with robust stability in time in the last years. The correlation between the low-elevation TRW chronology and temperature from the up-hill weather station is lower than the correlation recorded for high-elevation sites for most of the periods analyzed.
The stability of correlation between the BAI chronology from upper sites and growing season temperatures (24 May to 17 July) is stable in time with a slight decrease around the 2000s (Figure 5B). The correlation coefficient for the entire period is r = 0.441, with a minimum correlation of r = 0.298 for the period 1995–2009. Similarly, for lower sites, the BAI chronology moving correlation with previous summer temperatures (22 June to 15 July) records a reduction of the intensity in the middle of the analyzed periods, which is statistically non-significant. The decrease in correlation highlights a temporal instability in the climate sensitivity of BAI for both elevations, which has also been reported in other Norway spruce sites around Europe [34,44].
The time stability for BI correlation with air temperature was checked only for the high-elevation chronology (Figure 5C). Two periods were analyzed, one in the first part of the growing season (9 May to 28 July) and one in the late part (9 August to 8 October). The correlation between the temperature at the beginning of the growing season and BI index chronology has lower values in the first part of the analyzed period. The correlation is stable in time and significant only after the 1997–2010 period. The temperature from the last part of the growing season leads to a constant decrease in time of the correlation with the BI residual chronology, followed by a significant increase after 2004.
The running correlation for a 15-year period highlights the non-stationary correlation between TRW, BAI, and BI index chronologies, and air temperature. The continuous climate change leads to different responses of trees to temperature.

3. Materials and Methods

3.1. Study Area

The study was carried out in the Gheorgheni region, a large intramountain valley in the center of the Eastern Carpathians (Romania). The altitude in the study region ranges from 700 m a.s.l. to 1770 m a.s.l. The study area is located between 46°37′ N, 25°25′ E and 46°50′ N, 25°36′ E (Figure 6). The general geology of the study area is represented by a volcanogenic–sedimentary complex.
All the forests in the study area are managed forests, even-aged, and the Norway spruce is the main species. The forests belong to private and community owners and are managed by Gheorgheni Forest District manly in a high forest silvicultural system (clear-cutting system and shelterwood system).
In the study area mean annual temperature varies from 6.0 °C at the lower part of the valley to 4.8 °C at higher elevations (Figure 7). The average amount of annual precipitation is 542 mm at 750 m a.s.l. and 1004 mm in the upper part of the valley. The coldest month is January (the mean temperature at the lower part of the valley is −6.4 °C and −5.7 °C at higher elevations) and the warmest is July (the mean temperature at the lower part of the valley is 16.6 °C and 14.1 °C at up-hill weather station). The month with the most precipitation is June with 130 mm in the upper part of the valley and 90 mm at the lower part. The study area is characterized by frequent thermal inversion [59]. January is the month with the highest average of days with temperature inversion (13.8 days) followed by December (12.9 days) and November (11.9 days). During the spring/summer months, on average, less than 1–2 days with temperature inversion were recorded.
The climate diagrams according to Walter and Leith were created in R using the climatol package [78].

3.2. Sample Collection and Data Processing

To investigate the effect of altitude on the Norway spruce, climate response increment cores from 12 locations were collected: five locations at low elevation (altitude varied from 880 m a.s.l. to 1020 m a.s.l.) and seven locations at high elevation (altitude varied from 1510 m a.s.l. to 1630 m a.s.l.). In each plot, 15 to 20 dominant and healthy trees were selected to extract increment cores. By applying a sampling strategy, we selected mature trees, and the limitation of the analysis period to available temperature data ensured the exclusion of the juvenile growth part from the climate–growth relationship analysis. One core per tree was extracted using a 5-mm-diameter Pressler borer at breast height (i.e., 1.3 m). For each sampled tree, 2 perpendicular breast-height diameters (DBH) were measured with a forest caliper. The mean of these two diameters was used to compute the BAI. All three tree-ring proxies (TRW, BAI, and BI) were measured/calculated for the same cores. A subset of 50 cores for each altitude level (high and low elevations), were selected for measurements and analysis. In order to obtain reliable BI measurements, we selected only cores with no discoloration due to fungi, no gaps due to broken cores, and with parallel rings, as well as from cores containing pith or those that allowed easy determination of any missing rings.
Increment cores were dried, mounted on wooden supports, and sanded with successive sanding grits (from 80 grit to 400 grit) to ensure the flatness of the sanded surface and ring boundaries’ visibility. Sandpaper of 400 grains/mm2, according to [79], can assure the quality of BI measurements. After sanding, cores were scanned at 1200 DPI using an Epson Expression 12,000 XL scanner to measure TRW and BI. The scanning resolution (1200 DPI) used in this study represents the scanner’s hardware optical resolution, and we did not use an interpolation algorithm that could affect the BI measurements [79].
Norway spruce is a coniferous species with no visible differences between heartwood and sapwood. Based on this characteristic of spruce wood, it is possible to measure latewood blue reflectance without any chemical treatments [22,80,81,82]. To measure latewood blue reflectance, the standard protocol was followed [21,23]. Window parameter settings were adjusted according to [21]. The measured values of blue reflectance were transformed into blue intensity (BI) using Equation (1):
x i ( a d j ) = 2.56   x i 100
where xi represents the raw BI value for the year i; the constant 2.56 is used to ensure that all the values of BI are positive (all the possible values for xi are from 0 to 255). The inversion step was computed in CooRecorder 9.6 [83].
The tree-ring measurements (TRW and BI) were computed using CooRecorder software on scanned images [83]. No missing or false rings were observed.
The BAI was reconstructed for each ring individually and adjusted according to the mean DBH of each tree. The DBH of each tree (i) and for each year (t) was reconstructed by subtracting the doubled radial growth in the year t from the DBH in the year t. The calculation starts with the DBH measured in the year of coring (Equation (2)) [84].
D B H i ( t 1 ) = D B H i ( t ) 2 × T R W i ( t )
Based on the annual reconstructed DBH, the BAI for each tree was calculated using the following equation:
B A I t = π 4 × ( D B H t 2 D B H t 1 2 )
The BAI was calculated using an iterative function in R [85].
The measured individual series were cross-dated visually using TSAP-Win software [86]. COFECHA software was used for statistically cross-dating and measurement checking using correlation analysis of detrended 50-year intervals with a 25-year overlap [87,88].
A cubic smoothing spline function with a 50% frequency cutoff at 30 years was used to eliminate the age trend and any other disturbance signals in TRW and BAI chronologies [89,90]. According to [21], for the BI value, there is a possible increase in the juvenile period due to the transformation from blue reflectance to blue intensity. Because of this, the detrending of the BI series was made by fitting a Hugershoff curve to raw measurements [89].
For each analyzed tree-ring parameter, 2 chronologies were developed, one for low elevation and the second for high elevation. Indices were computed as the ratio between raw measurements and fitted values. In order to eliminate the autocorrelation that was still present in the standard index series, an autoregressive model was applied. In the analyses, we used the residual index chronologies obtained by bi-weight mean without variance stabilization. The chronologies quality and proprieties were assessed by calculating the classical dendrochronological statistics parameters as mean sensitivity, mean series intercorrelation, first-order autocorrelation, and expressed population signal. The mean sensitivity highlights the change in year-to-year tree growth [1]. The mean series intercorrelation (known as R-bar) and expressed population signal (known as EPS) express the signal strength in the analyzed chronologies and the common variability [49,89]. The carry-over effect of the previous year’s climate condition on the current year’s growth can be assessed by first-order autocorrelation [1,2]. The detrending, chronology development, and statistical parameters were computed using the dplR package on R software [91].

3.3. Climatic Dataset

Two local weather stations, Joseni (750 m a.s.l.)—valley weather station—and Bucin (1282 m a.s.l.)—up-hill weather station—provided the climatic data (daily mean temperature) (Figure 6 and Figure 7). Climatic data are available from 1963 at the Joseni weather station and from 1978 at the Bucin weather station. Therefore, the analyzed interval was limited to the common period, 1978–2019.

3.4. Climate–Growth Relationship Assessment

The availability of local daily air temperature data in the study area allowed us to analyze the climate–growth relationship on a cumulative daily scale [92]. Using the monthly temperature data, a common practice in dendroclimatology [89], to investigate the relationship between climate and tree growth conduce to a signal lost because of the artificial split of year in months, without any physiological bases [92,93]. The tree photosynthesis and wood increment are continuous processes and require more flexible time windows to investigate the relationship between tree growth and climate. The bootstrapped (n = 1000 replications) Pearson correlation between the TRW, BAI, and BI residual index chronologies and cumulative daily mean temperature was calculated [94,95]. Bootstrap samples were selected randomly with replacement from the entire interval (1978–2019) [96]. The significance of the bootstrapped correlation coefficient, at p < 0.05, was established according with [97]. The daily temperatures were aggregated in moving time windows of 21 days to 120 days starting from June in the previous year of growth to October of the current growing season. The time stability of correlations was checked for running windows with a length of 15 years, moved by one year [96]. The correlation between residual index chronologies and cumulative daily temperatures and time stability were computed using the dendroTools package in R [98].

4. Conclusions

Our results show that each tree-ring proxy contains a different climate signal. In the studied intramountain valley in the Eastern Carpathians, the growth of Norway spruce is influenced by winter temperatures, and signals demonstrating this are present in both TRW and BAI chronologies. The presence of temperature inversions can explain the unusual correlation pattern. The TRW and winter temperature correlation is unstable through time for both elevations, with evident temporal shifts after 1996. The highest correlations for TRW were obtained for cumulative windows of 45 days at high elevation and for cumulative windows of 113 days for low elevation; both cumulative windows are centered at the beginning of January. The BAI at high elevation contains a strong temperature signal from the middle of the growing season. In the case of BI, a clear climatic signal was observed only for the high-elevation chronology for cumulative windows starting in the second half of August to the beginning of October. The preliminary results obtained in this study need to be replicated for other intramountain valleys in the Carpathians to confirm the change in the general dendroclimatic pattern (growth driven by summer temperature) due to specific local climate modulated by thermal inversions.

Author Contributions

Conceptualization, A.P. and I.P.; methodology, A.P. and I.P.; software, A.P.; formal analysis, A.P.; investigation, A.P. and I.P.; resources, A.P. and I.P.; data curation, A.P. and I.P.; visualization, A.P. and C.-C.R.; writing—original draft preparation, A.P.; writing—review and editing, I.P., C.-C.R. and O.N.B.; supervision, I.P. and O.N.B.; project administration, I.P.; funding acquisition, I.P. and O.N.B. All authors have read and agreed to the published version of the manuscript.

Funding

A.P., I.P. and O.N.B. were supported by a grant of the Ministry of Research, Innovation and Digitization, CNCS-UEFISCDI, project number PN-III-P4-PCE-2021–1002, within PNCDI III and from project CresPerfInst (contract 34PFE/30.12.2021). C-C.R. was supported by a grant of the Ministry of Research, Innovation and Digitization, CNCS-UEFISCDI, project number PN-III-P1-1.1-TE-2021-1419, within PNCDI III.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on a reasonable request from the corresponding author.

Acknowledgments

Permission to work in the forests from the Gheorgheni area was granted to us by the Gheorgheni Forest District (Arnold Horvath). We want to thank to Mihai Balabașciuc for help in the field and in the laboratory. We would like to thank the reviewers for all valuable comments and suggestions, which helped us to improve the quality of the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Fritts, H. Tree Rings and Climate; Elsevier: Amsterdam, The Netherlands, 1976; ISBN 0323145280. [Google Scholar]
  2. Speer, J.H. Fundamentals of Tree-Ring Research; University of Arizona Press: Tucson, AZ, USA, 2010; ISBN 0816526842. [Google Scholar]
  3. Kramer, K.; Leinonen, I.; Loustau, D. The importance of phenology for the evaluation of impact of climate change on growth of boreal, temperate and Mediterranean forests ecosystems: An overview. Int. J. Biometeorol. 2000, 44, 67–75. [Google Scholar] [CrossRef]
  4. 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]
  5. European Environment Agency. European Forest Ecosystems: State and Trends; Publications Office of European Union: Luxembourg, 2016. [Google Scholar]
  6. Bigler, C.; Veblen, T.T. Increased early growth rates decrease longevities of conifers in subalpine forests. Oikos 2009, 118, 1130–1138. [Google Scholar] [CrossRef]
  7. Di Filippo, A.; Pederson, N.; Baliva, M.; Brunetti, M.; Dinella, A.; Kitamura, K.; Knapp, H.D.; Schirone, B.; Piovesan, G. The longevity of broadleaf deciduous trees in Northern Hemisphere temperate forests: Insights from tree-ring series. Front. Ecol. Evol. 2015, 3, 46. [Google Scholar] [CrossRef]
  8. Lindenmayer, D.B.; Franklin, J.F. Conserving Forest Biodiversity: A Comprehensive Multiscaled Approach; Island Press: Washington, DC, USA, 2002; ISBN 1559639350. [Google Scholar]
  9. Pretzsch, H.; Biber, P.; Schütze, G.; Uhl, E.; Rötzer, T. Forest stand growth dynamics in Central Europe have accelerated since 1870. Nat. Commun. 2014, 5, 4967. [Google Scholar] [CrossRef]
  10. Vitali, V.; Büntgen, U.; Bauhus, J. Seasonality matters—The effects of past and projected seasonal climate change on the growth of native and exotic conifer species in Central Europe. Dendrochronologia 2018, 48, 1–9. [Google Scholar] [CrossRef]
  11. Parobeková, Z.; Sedmáková, D.; Kucbel, S.; Pittner, J.; Jaloviar, P.; Saniga, M.; Balanda, M.; Vencurik, J. Influence of disturbances and climate on high-mountain Norway spruce forests in the Low Tatra Mts., Slovakia. For. Ecol. Manag. 2016, 380, 128–138. [Google Scholar] [CrossRef]
  12. Jones, P.D.; Briffa, K.R.; Osborn, T.J.; Lough, J.M.; van Ommen, T.D.; Vinther, B.M.; Luterbacher, J.; Wahl, E.R.; Zwiers, F.W.; Mann, M.E.; et al. High-resolution palaeoclimatology of the last millennium: A review of current status and future prospects. Holocene 2009, 19, 3–49. [Google Scholar] [CrossRef]
  13. Büntgen, U.; Frank, D.C.; Nievergelt, D.; Esper, J. Summer Temperature Variations in the European Alps, a.d. 755–2004. J. Clim. 2006, 19, 5606–5623. [Google Scholar] [CrossRef]
  14. Schweingruber, F.; Fritts, H.; Bräker, O.; Drew, L.; Schär, E. The X-ray technique as applied to dendroclimatology. Tree Ring Bull. 1978, 38, 61–91. [Google Scholar]
  15. Wilson, R.; Anchukaitis, K.; Briffa, K.R.; Büntgen, U.; Cook, E.; D’Arrigo, R.; Davi, N.; Esper, J.; Frank, D.; Gunnarson, B.; et al. Last millennium northern hemisphere summer temperatures from tree rings: Part I: The long term context. Quat. Sci. Rev. 2016, 134, 1–18. [Google Scholar] [CrossRef] [Green Version]
  16. Nagavciuc, V.; Roibu, C.-C.; Ionita, M.; Mursa, A.; Cotos, M.-G.; Popa, I. Different climate response of three tree ring proxies of Pinus sylvestris from the Eastern Carpathians, Romania. Dendrochronologia 2019, 54, 56–63. [Google Scholar] [CrossRef]
  17. Nagavciuc, V.; Kern, Z.; Ionita, M.; Hartl, C.; Konter, O.; Esper, J.; Popa, I. Climate signals in carbon and oxygen isotope ratios of Pinus cembra tree-ring cellulose from the Călimani Mountains, Romania. Int. J. Clim. 2020, 40, 2539–2556. [Google Scholar] [CrossRef]
  18. Koprowski, M.; Duncker, P. Tree ring width and wood density as the indicators of climatic factors and insect outbreaks affecting spruce growth. Ecol. Indic. 2012, 23, 332–337. [Google Scholar] [CrossRef]
  19. McCarroll, D.; Pettigrew, E.; Luckman, A.; Guibal, F.; Edouard, J.-L. Blue Reflectance Provides a Surrogate for Latewood Density of High-latitude Pine Tree Rings. Arct. Antarct. Alp. Res. 2002, 34, 450–453. [Google Scholar] [CrossRef]
  20. Campbell, R.; McCarroll, D.; Loader, N.J.; Grudd, H.; Robertson, I.; Jalkanen, R. Blue intensity in Pinus sylvestris tree-rings: Developing a new palaeoclimate proxy. Holocene 2007, 17, 821–828. [Google Scholar] [CrossRef]
  21. Rydval, M.; Larsson, L.-Å.; McGlynn, L.; Gunnarson, B.E.; Loader, N.J.; Young, G.H.F.; Wilson, R. Blue intensity for dendroclimatology: Should we have the blues? Experiments from Scotland. Dendrochronologia 2014, 32, 191–204. [Google Scholar] [CrossRef]
  22. Björklund, J.A.; Gunnarson, B.E.; Seftigen, K.; Esper, J.; Linderholm, H.W. Blue intensity and density from northern Fennoscandian tree rings, exploring the potential to improve summer temperature reconstructions with earlywood information. Clim. Past 2014, 10, 877–885. [Google Scholar] [CrossRef]
  23. Campbell, R.; McCarroll, D.; Robertson, I.; Loader, N.J.; Grudd, H.; Gunnarson, B. Blue Intensity In Pinus sylvestris Tree Rings: A Manual for A New Palaeoclimate Proxy. Tree-Ring Res. 2011, 67, 127–134. [Google Scholar] [CrossRef]
  24. Wilson, R.; Rao, R.; Rydval, M.; Wood, C.; Larsson, L.-Å.; Luckman, B.H. Blue Intensity for dendroclimatology: The BC blues: A case study from British Columbia, Canada. Holocene 2014, 24, 1428–1438. [Google Scholar] [CrossRef]
  25. Björklund, J.A.; Gunnarson, B.E.; Seftigen, K.; Esper, J.; Linderholm, H.W. Is Blue Intensity Ready to Replace Maximum Latewood Density as a Strong Temperature Proxy? A Tree-Ring Case Study on Scots Pine from Northern Sweden. Clim. Past Discuss. 2013, 9, 5227–5261. [Google Scholar] [CrossRef]
  26. Biondi, F. Comparing tree-ring chronologies and repeated timber inventories as forest monitoring tools. Ecol. Appl. 1999, 9, 12. [Google Scholar] [CrossRef]
  27. West, P. Use of diameter increment and basal area increment in tree growth studies. Can. J. For. Res. 1980, 10, 71–77. [Google Scholar] [CrossRef]
  28. Biondi, F.; Qeadan, F. A Theory-Driven Approach to Tree-Ring Standardization: Defining the Biological Trend from Expected Basal Area Increment. Tree-Ring Res. 2008, 64, 81–96. [Google Scholar] [CrossRef]
  29. Han, Y.; Wang, Y.; Liu, B.; Huang, R.; Camarero, J.J. Moisture mediates temperature-growth couplings of high-elevation shrubs in the Tibetan plateau. Trees 2022, 36, 273–281. [Google Scholar] [CrossRef]
  30. Caudullo, G.; Tinner, W.; de Rigo, D. Picea abies in Europe: Distribution, habitat, usage and threats. In European Atlas of Forest Tree Species; Publication Office of EU: Luxembourg, 2016. [Google Scholar]
  31. Klimo, E.; Hager, H.; Kulhavý, J. Spruce Monocultures in Central Europe: Problems and Prospects; European Forest Institute Joensuu: Joensuu, Finland, 2000; Volume 33. [Google Scholar]
  32. Paquette, A.; Messier, C. The effect of biodiversity on tree productivity: From temperate to boreal forests: The effect of biodiversity on the productivity. Glob. Ecol. Biogeogr. 2011, 20, 170–180. [Google Scholar] [CrossRef]
  33. Schütz, J.-P.; Götz, M.; Schmid, W.; Mandallaz, D. Vulnerability of spruce (Picea abies) and beech (Fagus sylvatica) forest stands to storms and consequences for silviculture. Eur. J. For. Res 2006, 125, 291–302. [Google Scholar] [CrossRef]
  34. Bošel’a, M.; Sedmák, R.; Sedmáková, D.; Marušák, R.; Kulla, L. Temporal shifts of climate–growth relationships of Norway spruce as an indicator of health decline in the Beskids, Slovakia. For. Ecol. Manag. 2014, 325, 108–117. [Google Scholar] [CrossRef]
  35. Levanič, T.; Gričar, J.; Gagen, M.; Jalkanen, R.; Loader, N.J.; McCarroll, D.; Oven, P.; Robertson, I. The climate sensitivity of Norway spruce [Picea abies (L.) Karst.] in the southeastern European Alps. Trees 2009, 23, 169–180. [Google Scholar] [CrossRef]
  36. Mäkinen, H.; Nöjd, P.; Kahle, H.-P.; Neumann, U.; Tveite, B.; Mielikäinen, K.; Röhle, H.; Spiecker, H. Radial growth variation of Norway spruce (Picea abies (L.) Karst.) across latitudinal and altitudinal gradients in central and northern Europe. For. Ecol. Manag. 2002, 171, 243–259. [Google Scholar] [CrossRef]
  37. Ponocná, T.; Spyt, B.; Kaczka, R.; Büntgen, U.; Treml, V. Growth trends and climate responses of Norway spruce along elevational gradients in East-Central Europe. Trees 2016, 30, 1633–1646. [Google Scholar] [CrossRef]
  38. Savva, Y.; Oleksyn, J.; Reich, P.B.; Tjoelker, M.G.; Vaganov, E.A.; Modrzynski, J. Interannual growth response of Norway spruce to climate along an altitudinal gradient in the Tatra Mountains, Poland. Trees 2006, 20, 735–746. [Google Scholar] [CrossRef]
  39. Koprowski, M.; Zielski, A. Dendrochronology of Norway spruce (Picea abies (L.) Karst.) from two range centres in lowland Poland. Trees 2006, 20, 383–390. [Google Scholar] [CrossRef]
  40. 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]
  41. 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]
  42. Van der Maaten-Theunissen, M.; Kahle, H.-P.; van der Maaten, E. Drought sensitivity of Norway spruce is higher than that of silver fir along an altitudinal gradient in southwestern Germany. Ann. For. Sci. 2013, 70, 185–193. [Google Scholar] [CrossRef]
  43. Barry, R.G. Mountain Weather and Climate; Psychology Press: London, UK, 1992; ISBN 0415071135. [Google Scholar]
  44. Kolář, T.; Čermák, P.; Trnka, M.; Žid, T.; Rybníček, M. Temporal changes in the climate sensitivity of Norway spruce and European beech along an elevation gradient in Central Europe. Agric. For. Meteorol. 2017, 239, 24–33. [Google Scholar] [CrossRef]
  45. Begović, K.; Rydval, M.; Mikac, S.; Čupić, S.; Svobodova, K.; Mikoláš, M.; Kozák, D.; Kameniar, O.; Frankovič, M.; Pavlin, J.; et al. Climate-growth relationships of Norway Spruce and silver fir in primary forests of the Croatian Dinaric mountains. Agric. For. Meteorol. 2020, 288–289, 108000. [Google Scholar] [CrossRef]
  46. Kaczka, R.J.; Spyt, B.; Janecka, K.; Beil, I.; Büntgen, U.; Scharnweber, T.; Nievergelt, D.; Wilmking, M. Different maximum latewood density and blue intensity measurements techniques reveal similar results. Dendrochronologia 2018, 49, 94–101. [Google Scholar] [CrossRef]
  47. Tsvetanov, N.; Dolgova, E.; Panayotov, M. First measurements of Blue intensity from Pinus peuce and Pinus heldreichii tree rings and potential for climate reconstructions. Dendrochronologia 2020, 60, 125681. [Google Scholar] [CrossRef]
  48. Schwab, N.; Kaczka, R.; Janecka, K.; Böhner, J.; Chaudhary, R.; Scholten, T.; Schickhoff, U. Climate Change-Induced Shift of Tree Growth Sensitivity at a Central Himalayan Treeline Ecotone. Forests 2018, 9, 267. [Google Scholar] [CrossRef] [Green Version]
  49. Wigley, T.M.; Briffa, K.R.; Jones, P.D. On the average value of correlated time series, with applications in dendroclimatology and hydrometeorology. J. Appl. Meteorol. Climatol. 1984, 23, 201–213. [Google Scholar] [CrossRef]
  50. Sidor, C.G.; Popa, I.; Vlad, R.; Cherubini, P. Different tree-ring responses of Norway spruce to air temperature across an altitudinal gradient in the Eastern Carpathians (Romania). Trees 2015, 29, 985–997. [Google Scholar] [CrossRef]
  51. 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]
  52. Leonelli, G.; Pelfini, M. Influence of climate and climate anomalies on norway spruce tree-ring growth at different altitudes and on glacier responses: Examples from the central italian alps. Geogr. Ann. Ser. A Phys. Geogr. 2008, 90, 75–86. [Google Scholar] [CrossRef]
  53. Apăvaloae, M.; Apostol, L.; Pîrvulescu, I. Inversiunile termice din culoarul Moldovei (sectorul Câmpulung Moldovenesc–Frasin) și influența lor asupra poluării atmosferei. Sci. Ann. “Stefan Cel Mare” Univ. Geogr. Ser. 1994, 5. [Google Scholar]
  54. Ahrens, C.D. Meteorology Today: An Introduction to Weather, Climate, and the Environment; Cengage Learning Canada, Inc.: Toronto, ON, Canada, 2015; ISBN 0176728333. [Google Scholar]
  55. Apăvaloae, M.; Pîrvulescu, I.; Apostol, L. Caracteristici ale inversiunilor termice din Podișul Fălticenilor. Lucr. Semin. Geogr. “Dimitrie Cantemir” 1987, 8. [Google Scholar]
  56. Ciutea, A.; Jitariu, V. Thermal inversions identification through the analysis of the vegetation inversions occurred in the forest ecosystems from the Eastern Carpathians. Present Environ. Sustain. Dev. 2020, 14, 29–42. [Google Scholar] [CrossRef]
  57. Ichim, P.; Apostol, L.; Sfîcă, L.; Kadhim-Abid, A.-L.; Istrate, V. Frequency of Thermal Inversions Between Siret and Prut Rivers in 2013. Present Environ. Sustain. Dev. 2014, 8, 267–284. [Google Scholar] [CrossRef]
  58. Sfîcă, L.; Nicuriuc, I.; Niță, A. Boundary Layer Temperature Stratification as Result of Atmospheric Circulation within the Western Side of Brașov Depression. In Proceedings of the 2019 Air and Water Components of the Environment Conference, Cluj-Napoca, Romania, 22–24 March 2019; pp. 53–64. [Google Scholar]
  59. Palffy, E. Temperature inversion in the Csik basin. Acta Clim. 1995, 28, 41–45. [Google Scholar]
  60. Mayr, S.; Schmid, P.; Beikircher, B.; Feng, F.; Badel, E. Die hard: Timberline conifers survive annual winter embolism. New Phytol. 2020, 226, 13–20. [Google Scholar] [CrossRef] [PubMed]
  61. Rossi, S.; Deslauriers, A.; Anfodillo, T.; Carraro, V. Evidence of threshold temperatures for xylogenesis in conifers at high altitudes. Oecologia 2007, 152, 1–12. [Google Scholar] [CrossRef] [PubMed]
  62. Rossi, S.; Rathgeber, C.B.K.; Deslauriers, A. Comparing needle and shoot phenology with xylem development on three conifer species in Italy. Ann. For. Sci. 2009, 66, 206. [Google Scholar] [CrossRef]
  63. Buntgen, U.; Frank, D.C.; Kaczka, R.J.; Verstege, A.; Zwijacz-Kozica, T.; Esper, J. Growth responses to climate in a multi-species tree-ring network in the Western Carpathian Tatra Mountains, Poland and Slovakia. Tree Physiol. 2007, 27, 689–702. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  64. Putalová, T.; Vacek, Z.; Vacek, S.; Štefančík, I.; Bulušek, D.; Král, J. Tree-ring widths as an indicator of air pollution stress and climate conditions in different Norway spruce forest stands in the Krkonoše Mts. Cent. Eur. For. J. 2019, 65, 21–33. [Google Scholar] [CrossRef]
  65. Leal, S.; Melvin, T.M.; Grabner, M.; Wimmer, R.; Briffa, K.R. Tree-ring growth variability in the Austrian Alps: The influence of site, altitude, tree species and climate. Boreas 2007, 36, 426–440. [Google Scholar] [CrossRef]
  66. Matisons, R.; Elferts, D.; Krišāns, O.; Schneck, V.; Gärtner, H.; Wojda, T.; Kowalczyk, J.; Jansons, Ā. Nonlinear Weather–Growth Relationships Suggest Disproportional Growth Changes of Norway Spruce in the Eastern Baltic Region. Forests 2021, 12, 661. [Google Scholar] [CrossRef]
  67. Selås, V.; Piovesan, G.; Adams, J.M.; Bernabei, M. Climatic factors controlling reproduction and growth of Norway spruce in southern Norway. Can. J. For. Res. 2002, 32, 9. [Google Scholar] [CrossRef]
  68. Hansen, J.; Beck, E. Seasonal changes in the utilization and turnover of assimilation products in 8-year-old Scots pine (Pinus sylvestris L.) trees. Trees 1994, 8, 172–182. [Google Scholar] [CrossRef]
  69. Holl, W. Seasonal Fluctuation of Reserve Materials in the Trunkwood of Spruce [Picea abies (L.) Karst.]. J. Plant Physiol. 1985, 117, 355–362. [Google Scholar] [CrossRef]
  70. Kozlowski, T.T.; Pallardy, S.G. Growth Control in Woody Plants; Elsevier: Amsterdam, The Netherlands, 1997; ISBN 0080532683. [Google Scholar]
  71. Gindl, W.; Grabner, M.; Wimmer, R. The influence of temperature on latewood lignin content in treeline Norway spruce compared with maximum density and ring width. Trees 2000, 14, 409–414. [Google Scholar] [CrossRef]
  72. Hacket-Pain, A.; Ascoli, D.; Berretti, R.; Mencuccini, M.; Motta, R.; Nola, P.; Piussi, P.; Ruffinatto, F.; Vacchiano, G. Temperature and masting control Norway spruce growth, but with high individual tree variability. For. Ecol. Manag. 2019, 438, 142–150. [Google Scholar] [CrossRef]
  73. Matisons, R.; Elferts, D.; Krišāns, O.; Schneck, V.; Gärtner, H.; Bast, A.; Wojda, T.; Kowalczyk, J.; Jansons, Ā. Non-linear regional weather-growth relationships indicate limited adaptability of the eastern Baltic Scots pine. For. Ecol. Manag. 2021, 479, 118600. [Google Scholar] [CrossRef]
  74. Akhmetzyanov, L.; Sánchez-Salguero, R.; García-González, I.; Buras, A.; Dominguez-Delmás, M.; Mohren, F.; den Ouden, J.; Sass-Klaassen, U. Towards a new approach for dendroprovenancing pines in the Mediterranean Iberian Peninsula. Dendrochronologia 2020, 60, 125688. [Google Scholar] [CrossRef]
  75. Fuentes, M.; Salo, R.; Björklund, J.; Seftigen, K.; Zhang, P.; Gunnarson, B.; Aravena, J.-C.; Linderholm, H.W. A 970-year-long summer temperature reconstruction from Rogen, west-central Sweden, based on blue intensity from tree rings. Holocene 2018, 28, 254–266. [Google Scholar] [CrossRef]
  76. Știrbu, M.-I.; Roibu, C.-C.; Carrer, M.; Mursa, A.; Unterholzner, L.; Prendin, A.L. Contrasting Climate Sensitivity of Pinus cembra Tree-Ring Traits in the Carpathians. Front. Plant Sci. 2022, 13, 855003. [Google Scholar] [CrossRef] [PubMed]
  77. Walter, H.; Lieth, H. Klimadiagramm-Weltatlas: Von Heinrich Walter Und Helmut Lieth; Gustav Fischer Verlag Jena: Jena, Germany, 1967. [Google Scholar]
  78. Guijarro, J.A.; Guijarro, M.J.A. Package ‘Climatol’. 2019. Available online: https://cran.r-project.org/web/packages/climatol/climatol.pdf (accessed on 20 April 2020).
  79. Babst, F.; Frank, D.; Büntgen, U.; Nievergelt, D.; Esper, J. Effect of sample preparation and scanning resolution on the Blue Reflectance of Picea abies. TRACE Proc. 2009, 7, 188–195. [Google Scholar]
  80. Björklund, J.; Gunnarson, B.E.; Seftigen, K.; Zhang, P.; Linderholm, H.W. Using adjusted Blue Intensity data to attain high-quality summer temperature information: A case study from Central Scandinavia. Holocene 2015, 25, 547–556. [Google Scholar] [CrossRef]
  81. Wilson, R.; D’Arrigo, R.; Andreu-Hayles, L.; Oelkers, R.; Wiles, G.; Anchukaitis, K.; Davi, N. Experiments based on blue intensity for reconstructing North Pacific temperatures along the Gulf of Alaska. Clim. Past 2017, 13, 1007–1022. [Google Scholar] [CrossRef]
  82. Wilson, R.; Anchukaitis, K.; Andreu-Hayles, L.; Cook, E.; D’Arrigo, R.; Davi, N.; Haberbauer, L.; Krusic, P.; Luckman, B.; Morimoto, D.; et al. Improved dendroclimatic calibration using blue intensity in the southern Yukon. Holocene 2019, 29, 1817–1830. [Google Scholar] [CrossRef]
  83. Maxwell, R.S.; Larsson, L.-A. Measuring tree-ring widths using the CooRecorder software application. Dendrochronologia 2021, 67, 125841. [Google Scholar] [CrossRef]
  84. Bosela, M.; Tumajer, J.; Cienciala, E.; Dobor, L.; Kulla, L.; Marčiš, P.; Popa, I.; Sedmák, R.; Sedmáková, D.; Sitko, R.; et al. Climate warming induced synchronous growth decline in Norway spruce populations across biogeographical gradients since 2000. Sci. Total Environ. 2021, 752, 141794. [Google Scholar] [CrossRef] [PubMed]
  85. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2022. [Google Scholar]
  86. Rinn Tech. TSAPWin Scientific: Time Series Analysis and Presentation for Dendrochronology and Related Applications; Rinn Tech: Heidelberg, Germany, 2012. [Google Scholar]
  87. Holmes, R. Computer assisted quality control. Tree-Ring Bull 1983, 43, 69–78. [Google Scholar]
  88. Grissino-Mayer, H.D. Evaluating crossdating accuracy: A manual and tutorial for the computer program COFECHA. Tree Ring Res. 2001, 57, 205–221. [Google Scholar]
  89. Cook, E.R.; Kairiukstis, L.A. Methods of Dendrochronology: Applications in the Environmental Sciences; Springer Science & Business Media: Berlin/Heidelberg, Germany, 1990; ISBN 9401578796. [Google Scholar]
  90. Schweingruber, F.H. Tree Rings: Basics and Applications of Dendrochronology; Springer Science & Business Media: Berlin/Heidelberg, Germany, 1989; ISBN 9400912730. [Google Scholar]
  91. Bunn, A.G. A dendrochronology program library in R (dplR). Dendrochronologia 2008, 26, 115–124. [Google Scholar] [CrossRef]
  92. Jevšenak, J. Daily climate data reveal stronger climate-growth relationships for an extended European tree-ring network. Quat. Sci. Rev. 2019, 221, 105868. [Google Scholar] [CrossRef]
  93. Beck, W.; Sanders, T.G.M.; Pofahl, U. CLIMTREG: Detecting temporal changes in climate–growth reactions—A computer program using intra-annual daily and yearly moving time intervals of variable width. Dendrochronologia 2013, 31, 232–241. [Google Scholar] [CrossRef]
  94. Efron, B.; Tibshirani, R. Bootstrap methods for standard errors, confidence intervals, and other measures of statistical accuracy. Stat. Sci. 1986, 1, 54–75. [Google Scholar] [CrossRef]
  95. Guiot, J. The bootstrapped response function. Tree-Ring Bull 1991, 51, 39–41. [Google Scholar]
  96. Biondi, F.; Waikul, K. DENDROCLIM2002: A C++ program for statistical calibration of climate signals in tree-ring chronologies. Comput. Geosci. 2004, 30, 303–311. [Google Scholar] [CrossRef]
  97. Dixon, P. Bootstrap resampling. In The Encyclopedia of Environmetrics; El-Shaarawi, A.H., Piegorsch, W.W., Eds.; Springer: Berlin/Heidelberg, Germany, 2001. [Google Scholar]
  98. Jevšenak, J.; Levanič, T. dendroTools: R package for studying linear and nonlinear responses between tree-rings and daily environmental data. Dendrochronologia 2018, 48, 32–39. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Norway spruce TRW, BAI, and BI residual index series chronologies for 1978–2019 (the shadowed line represents the individual series and the thicker line represents the master chronologies).
Figure 1. Norway spruce TRW, BAI, and BI residual index series chronologies for 1978–2019 (the shadowed line represents the individual series and the thicker line represents the master chronologies).
Plants 11 02428 g001
Figure 2. Correlation between TRW residual index chronologies and cumulative daily temperature from valley weather station (Joseni) (A) and from up-hill weather station (Bucin) (B) (vertical black line represents the limit between previous (with lowercase) and current (with uppercase) year).
Figure 2. Correlation between TRW residual index chronologies and cumulative daily temperature from valley weather station (Joseni) (A) and from up-hill weather station (Bucin) (B) (vertical black line represents the limit between previous (with lowercase) and current (with uppercase) year).
Plants 11 02428 g002
Figure 3. Correlation between BAI residual index chronologies and cumulative daily temperature from valley weather station (Joseni) (A) and from up-hill weather station (Bucin) (B) (vertical black line represents the limit between previous (with lowercase) and current (with uppercase) year).
Figure 3. Correlation between BAI residual index chronologies and cumulative daily temperature from valley weather station (Joseni) (A) and from up-hill weather station (Bucin) (B) (vertical black line represents the limit between previous (with lowercase) and current (with uppercase) year).
Plants 11 02428 g003
Figure 4. Correlation between BI residual index chronologies and cumulative daily temperature from valley weather station (Joseni) (A) and from up-hill weather station (Bucin) (B) (vertical black line represents the limit between previous (with lowercase) and current (with uppercase) year).
Figure 4. Correlation between BI residual index chronologies and cumulative daily temperature from valley weather station (Joseni) (A) and from up-hill weather station (Bucin) (B) (vertical black line represents the limit between previous (with lowercase) and current (with uppercase) year).
Plants 11 02428 g004
Figure 5. The 15-year running correlations between air temperature and tree-ring proxies ((A)—tree-ring width; (B)—basal area increment; (C)—blue intensity). The legend contains the site location and period (with the dashed line representing a significant correlation coefficient for p < 0.05, n = 15).
Figure 5. The 15-year running correlations between air temperature and tree-ring proxies ((A)—tree-ring width; (B)—basal area increment; (C)—blue intensity). The legend contains the site location and period (with the dashed line representing a significant correlation coefficient for p < 0.05, n = 15).
Plants 11 02428 g005
Figure 6. The study localization—with blue the high elevation sites, up-hill weather station, and with red the low elevation sites, valley weather station (image source: https://geoportal.ancpi.ro (accessed on 1 August 2022); https://earth.google.com (accessed on 1 August 2022)).
Figure 6. The study localization—with blue the high elevation sites, up-hill weather station, and with red the low elevation sites, valley weather station (image source: https://geoportal.ancpi.ro (accessed on 1 August 2022); https://earth.google.com (accessed on 1 August 2022)).
Plants 11 02428 g006
Figure 7. Climate diagrams according with Walter and Leith [77] for the local weather stations in the study area (valley station on the right and up-hill station on the left). The red line represents the variations of the mean monthly temperatures, the blue line represents the variations of the monthly precipitation, the blue hatched area represents the humid period, and the blue filled area represents the wet period. The bar in the lower part suggests the indication of months where frost is likely.
Figure 7. Climate diagrams according with Walter and Leith [77] for the local weather stations in the study area (valley station on the right and up-hill station on the left). The red line represents the variations of the mean monthly temperatures, the blue line represents the variations of the monthly precipitation, the blue hatched area represents the humid period, and the blue filled area represents the wet period. The bar in the lower part suggests the indication of months where frost is likely.
Plants 11 02428 g007
Table 1. Basis statistics for the six chronologies: TRW—tree-ring width; BAI—basal area increment; BI—blue intensity; AGR—average growth rate (mm year−1 for TRW, mm2 year−1 for BAI, and no unit for BI); SD—standard deviation; Rbar—series intercorrelation; Mean sens—mean sensitivity; Auto corr.—serial autocorrelation.
Table 1. Basis statistics for the six chronologies: TRW—tree-ring width; BAI—basal area increment; BI—blue intensity; AGR—average growth rate (mm year−1 for TRW, mm2 year−1 for BAI, and no unit for BI); SD—standard deviation; Rbar—series intercorrelation; Mean sens—mean sensitivity; Auto corr.—serial autocorrelation.
Tree Ring ProxySeries LengthAGR ± SDRbarMean SensAuto Corr.
MeanMinMax
Low Elevation Chronologies
TRW95791352.16 ± 0.530.3040.1760.869
BAI1445.7 ± 614.70.2990.2110.773
BI2.27 ± 0.170.1460.0630.548
High Elevation Chronologies
TRW101681331.77 ± 0.540.3350.1500.857
BAI1003.6 ± 541.60.3330.1630.844
BI2.15 ± 0.200.2870.0770.459
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Popa, A.; Popa, I.; Roibu, C.-C.; Badea, O.N. Do Different Tree-Ring Proxies Contain Different Temperature Signals? A Case Study of Norway Spruce (Picea abies (L.) Karst) in the Eastern Carpathians. Plants 2022, 11, 2428. https://doi.org/10.3390/plants11182428

AMA Style

Popa A, Popa I, Roibu C-C, Badea ON. Do Different Tree-Ring Proxies Contain Different Temperature Signals? A Case Study of Norway Spruce (Picea abies (L.) Karst) in the Eastern Carpathians. Plants. 2022; 11(18):2428. https://doi.org/10.3390/plants11182428

Chicago/Turabian Style

Popa, Andrei, Ionel Popa, Cătălin-Constantin Roibu, and Ovidiu Nicolae Badea. 2022. "Do Different Tree-Ring Proxies Contain Different Temperature Signals? A Case Study of Norway Spruce (Picea abies (L.) Karst) in the Eastern Carpathians" Plants 11, no. 18: 2428. https://doi.org/10.3390/plants11182428

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