Increasing radial and latewood growth rates of Larix cajanderi Mayr. and Pinus sylvestris L. in the continuous permafrost zone in Central Yakutia (Russia)

Abstract

Key message

Annual radial growth and latewood formation increase during the period 1966–2015 in Larix cajanderi Mayr. and Pinus sylvestris L. along a latitudinal gradient on Siberian permafrost. Warmer temperatures and precipitation at the onset of the growth season favor growth in both species.

Context

A vast area of Siberian boreal forest is covered by continuous permafrost, where tree radial growth is strongly limited by low temperatures. Understanding how climate controls Siberian conifer growth is critical in order to predict their response to ongoing climate change.

Aims

To explore variations in the radial growth and earlywood/latewood formation of Pinus sylvestris L. and Larix cajanderi Mayr. trees along a climate gradient within the continuous permafrost zone in Central Yakutia (Russia), a region experiencing accelerated warming.

Methods

We evaluated spatiotemporal differences in tree-ring and latewood width among sites by using generalized additive mixed models. Pearson’s correlations were used to analyze the effect of local climate (temperature, precipitation) and drought (SPEI) on the formation of earlywood and latewood.

Results

Both conifers showed increased radial growth in recent decades. Radial growth was controlled by temperature and precipitation in both species along the gradient, with earlier precipitation events being critical for the formation of wide rings.

Conclusion

Although temperature is considered as the main limiting factor triggering tree growth in the continuous permafrost zone, our results showed that water availability also plays a key role in radial tree growth even in the coldest sites.

1 Introduction

Boreal forests play a leading role in regional and global climate systems, affecting the global carbon balance and water cycle (Melnikov and Drozdov 2006; Falkowski et al. 2010; Schlesinger and Jasechko 2014). However, up to 80% of these forests occur in the permafrost zone (Helbig et al. 2016) and are thus highly sensitive to changes in climate (Gartenr 1995; Sperry et al. 2006). A vast area of Siberian boreal forest in high latitudes is covered by continuous permafrost (Duchkov 2006), where the active soil layer remains frozen for most of the year (Lopez et al. 2010; Nikolaev et al. 2011) and tree growth is strongly limited by spring–summer temperatures (Briffa et al. 1998, 2008; Vaganov et al. 1999; Hellmann et al. 2016).

The accelerated warming observed in northeastern Siberia during recent decades (Gorokhov and Fedorov 2018) has induced an increase in the active soil layer depth, facilitating tree growth (Bosikov 1998; Fedorov and Konstantinov 2003). The impact of warming (IPCC 2014) might, therefore, be positive for tree growth in this environment due to the lengthening of the growing season (Nikolaev et al. 2011). However, the increase of temperature over a certain threshold might also lead to an excess of soil moisture, triggering the degradation of permafrost and leading to forest fragmentation and forest loss (Lopez et al. 2008; Baltzer et al. 2014; Jorgenson and Osterkamp 2005; Carpino et al. 2018), potentially compromising the boreal forest role in global carbon and water cycles. Therefore, understanding tree responses to a wide range of environmental conditions within the continuous permafrost zone might provide a clearer insight into the factors controlling tree performance and ultimately forest functioning.

Tracheid cells make up more than 90% of xylem in conifers (Vaganov et al. 2006), combining hydraulic (earlywood) and structural (latewood) functions (Tyree et al. 2002). Since tracheid formation, anatomy, and functioning are strongly related to climate conditions (Antonova and Stasova 1993, 1997; Zweifel et al. 2006; Olano et al. 2012; Babushkina et al. 2018), they encode different environmental information (Vaganov 1990; Fonti et al. 2010; Arzac et al. 2018a, b; Popkova et al. 2018). Thus, an independent analysis of the climatic factors controlling the formation of different tree-ring sectors (i.e., earlywood and latewood) is critical to determine the extent to which climate control on tree growth varies along the growing season.

The intra-seasonal change in xylem morphology (from earlywood to latewood cells) might be driven by the need for structural reinforcement by latewood (e.g., Gartenr 1995; Sperry et al. 2006), photoperiod (e.g., Gyllenstrand et al. 2007), or changes in water availability (e.g., Antonova and Stasova 1997; Olano et al. 2014) throughout the growing season. Since earlywood and latewood cells are related towards different functions (Björklund et al. 2017), the ratio between the two tree-ring sectors plays a critical role in tree water–carbon interactions (Hacke et al. 2001; Domec and Gartner 2002). In the permafrost zone, temperature determines growing season length due to the timing of permafrost thawing (Lopez et al. 2010; Nikolaev et al. 2011), controlling the cell production and radial cell expansion (Vaganov et al. 1999), and, also, in the secondary wall thickening and the formation of latewood cells, as evidenced by the high correlations between maximum latewood density and summer temperature in the permafrost area (Briffa et al. 1998; Vaganov et al. 1999).

In this work, we analyzed the climatic factors controlling the radial growth (tree-ring width) and earlywood/latewood formation of two dominant conifers of Siberian boreal forests, the deciduous and frost-resistant Larix cajanderi Mayr and the evergreen Pinus sylvestris L. along a southwest–northeast gradient on the continuous permafrost zone in Central Yakutia. P. sylvestris grows mostly on dry sands or sandy loams, with high permafrost thawing (Nikolaev 2003; Nikolaev et al. 2011), whereas L. cajanderi grows on moderately moist loams (Nikolaev et al. 2009), and is the dominant conifer in the permafrost zone (Abaimov 2010). Also, the two conifers develop shallow root systems that enable them to thrive despite the occurrence of permafrost in the deeper soil horizons (Nikolaev et al. 2009, 2011).

Central Yakutia is experiencing an accelerated increase in temperatures in recent years (Gorokhov and Fedorov 2018), which may have extended the tree-growing season. Along an altitudinal/latitudinal gradient in Central Yakutia, we explored (a) whether radial growth has increased in response to warmer conditions. Additionally, we evaluated the relative roles of temperature and water availability in the two conifers. We hypothesized that (b) temperature and soil moisture availability would be critical for the formation of earlywood, while temperature would be the main factor triggering the formation of latewood along the gradient. Finally, we hypothesized that (c) higher growth rates associated with wide rings would be linked with an earlier timing of precipitation, expanding the period of soil moisture availability.

2 Materials and methods

2.1 Study area and sampling design

We selected four sampling sites along a 1500-km southwest–northeast climatic gradient of increasing altitude within the continuous permafrost zone in Central Yakutia, from Yakutsk to Oymyakon (northeastern Siberia, Russia; Fig. 1). The region is characterized by extreme cold and dry continental climate (with an annual thermal amplitude of 60 °C). The first site (S1), the most southwestern, is located at 206 m asl near the locality of Yakutsk. The second site (S2) is located at 245 m asl close to the locality of Krest-Khaldzhay. The third site (S3) is located at 465 m asl near the locality of Tompo. Finally, the fourth site (S4), the most northeastern, is located at 626 m asl close to the locality of Oymyakon. The mean annual temperature along the gradient ranged from − 9.9 to − 18.6 °C and the total annual precipitation from 269 to 377 mm, with S3 being the wettest site (see Table 1 for further details). Climate data was obtained from a high-resolution (0.5° grids) CRU data set for a 50-year period (1966–2015), available at https://crudata.uea.ac.uk/cru/data/hrg/cru_ts_3.23/ge/.

Fig. 1

Study area and comparison of mean yearly temperature and total yearly precipitation among sampling sites. The meteorological data covers the common period 1966–2015 and was from 0.5°-resolution CRU grids

Table 1 Characterization of the sampling sites. Meteorological data was obtained by 0.5°-resolution CRU grids for the common period 1966–2015. Lc, Larix cajanderi; Ps, Pinus sylvestris

At each site (except in S4), we sampled mono-specific forest stands of L. cajanderi and P. sylvestris during summer 2015. Forest stands at each site were around 10–20 km apart. Larch trees were sampled at all the sites, whereas pine trees were sampled in all except the coolest S4. A total of 125 dominant or codominant trees were sampled (74 larches and 51 pines; Table 1), and two cores per tree were taken at breast height with a 5-mm diameter increment borer. Cores were air-dried and polished with successively finer grades of sandpaper until the xylem structure was visible under magnification.

2.2 Tree-ring measurements and chronologies

After visual cross-dating, tree-ring width (RW), earlywood width (EW), and latewood width (LW) were measured to the nearest 0.01 mm by using a LINTAB-5 sliding stage micrometer interfaced with the software TSAP Win (RINNTECH, Heidelberg, Germany) in a total of 12,156 annual rings. Cross-dating accuracy was checked using the software COFECHA (Grissino-Mayer 2001), and each raw series was standardized using the software ARSTAN (Cook and Holmes 1996). Spline functions were fitted to the raw series to a spline function with a 50% frequency cutoff at 32 years, which was flexible enough to reduce the non-climatic variance preserving high-frequency climatic information (Cook and Peters 1981). In addition, dendrochronological statistical measures (e.g., mean sensitivity and expressed population signal) were obtained from ARSTAN. Due to the high correlation between EW and LW chronologies (mean correlations r = 0.7 for larch and r = 0.5 for pine), a linear regression was used to remove the dependence of LW on EW (Babst et al. 2016), obtaining an adjusted latewood index (LW_adj) (Meko and Baisan 2001) non-correlated with EW chronology (Stahle et al. 2009).

2.3 Pointer years

We analyzed the climatic factors driving the formation of wide and narrow tree-rings (pointer years) for each species and site during a 50-year period (1966–2015). The definition of pointer years was based on the deviation from mean values of the tree-ring width index chronologies. Thus, years with RW index 1 SD above the mean RW chronology were considered as wide rings, whereas years with RW index 1 SD below the mean RW chronology were considered as narrow rings (Shi et al. 2008; Tychkov et al. 2019). Once pointer years were defined, monthly precipitation and temperature were averaged for pointer years to compare the environmental conditions triggering the occurrence of wide and narrow rings.

2.4 Statistical analysis

We performed generalized additive mixed models (GAMM, Zuur et al. 2009) to assess inter-site differences in whole ring and latewood width (absolute values) as well as temporal trends in tree growth. We evaluated trends in growth rate using “year” as a linear factor; if growth has increased/reduced along the studied period, “year” would have a positive/negative effect and no significant effect of “year” would mean no temporal trend. Decrease in growth with tree age may obscure temporal growth trend effect. Therefore, to remove this confounding factor, we modeled age-related radial growth changes with a non-linear function using a smoother. Additionally, we considered measurements nested within tree and temporally autocorrelated. Thus, each tree-ring trait was evaluated as a function of the site (nominal variable), year (linear variable representing the temporal trend), tree age (modeled using a smoother), and nested within tree with a temporal autocorrelation. Analyses were performed in R, version 3.3.1 (R Development Core Team 2016) with the “nmle” package.

A climate analysis was performed to evaluate the trend in the mean annual temperature and the total annual precipitation for the 1966–2015 period in all the sites. Additionally, we also evaluated seasonal trends in climate. We used Pearson’s correlations to explore the relationship between residual RW, EW, and LW_adj chronologies and monthly climatic factors for the common period (1966–2015) at each site, for a time window from August of the previous year to September of the growth year (July in the case of EW). Climate data series were standardized following the same procedure as for the tree-ring chronologies.

Finally, we used the standardized precipitation–evapotranspiration index (SPEI, Vicente-Serrano et al. 2010), available at http://spei.csic.es/database.html, to evaluate the differential effect of soil water balance on the analyzed chronologies at different temporal scales in which the water deficit is accumulated. We calculated Pearson’s correlation between RW, EW, and LW_adj chronologies and SPEI from September prior to the growth year to the current (growth year) August over different time scales (1–18 months), which are also available at the SPEI database.

3 Results

3.1 Climate trends

During the period 1966–2015, mean annual temperature has increased at a rate between 0.47 °C/decade (R² = 0.45; P < 0.001) at site S4 to 0.65 °C/decade (R² = 0.55; P < 0.001) at site S1 (Fig. 2). The tree growth season, defined by daily mean temperature above 5 °C, extended from May to September in S1 and S2 and from June to August in S3 and S4. Interestingly, the lower southwestern sites (S1 and S2) have experienced a greater increase in summer and winter temperatures, while the higher northeastern sites (S3 and S4) experienced the highest rates of temperature increase in spring and autumn (data not shown). No significant trend in precipitation was detected in any of the sites for the study period.

Fig. 2

Mean annual temperature trends for the four sampling sites during the common period 1966–2015. Mean annual temperature trends (a). Variation in temperature by seasons (b)

3.2 Variation of radial growth and latewood width

Tree age ranged from 75 to 148 years (103 ± 36; mean ± SD) for larch to 84–128 years (105 ± 37) for pine. RW was 0.57–1.11 mm (0.78 ± 0.23; mean ± SD) in larch and 0.56–1.29 mm (0.94 ± 0.37) in pine (see Annex Table 4 for further details). Ring width chronologies showed values of mean correlation between trees (r_bt) from 0.32 to 0.43 for larch and 0.18 to 0.40 for pine, while RW residual chronologies mean sensitivity (m_sx) ranged from 0.17 to 0.21 for larch and 0.09 to 0.16 for pine. Expressed population signal (EPS) was over 0.85 for all the sites except pines in S3 (Annex Table 4), suggesting adequate replication and a high common signal in each locality except for pine chronology in S3. Within-species RW correlations showed non-significant values for larch sites and moderate but significant values for pine sites (r = 0.34–0.48; P < 0.001, Annex Table 5). Between-species correlation was only significant in the lowest site (S1, r = 0.56; P < 0.001, Annex Table 6).

GAMM produced a satisfactory model for tree age in both species and parameters, suggesting a negative effect of tree age on radial and latewood growth (Table 2; Annex Fig. 6). At the same time, “year” (as linear variable representing the temporal trend) had a positive effect in all the models indicating that RW and LW growth rates had increased throughout the life of sampled trees for both species. RW had increased at a rate of 1.7 μm/decade in pine and 3.5 μm/decade in larch, whereas LW increased at a rate of 0.4 and 1.2 μm/decade in pine and larch, respectively (Table 2). Pairwise comparisons also indicated growth rate differences between sites. Pine RW showed the lowest growth rate in S1, being marginally higher at S3 than S2. Pine latewood growth rate was highest at S3, with no differences between S1 and S2. Larch RW growth was higher in S3 than S1 and S2, with S4 being intermediate (marginally lower than S3 and higher than S1, S2). Finally, larch LW growth was also higher in S3, but there was no difference between the other sites (Table 2; Fig. 3).

Table 2 Estimated parameters of the generalized additive model for radial growth and latewood in Larix cajanderi and Pinus sylvestris

Fig. 3

Descriptive diagram of the main results from generalized additive mixed model, showing the components of the model and its effects on ring width (RW) and latewood width (LW) on Pinus sylvestris and Larix cajanderi

3.3 Climate–growth response

Pearson’s correlations between standardized climatic variables and tree-ring trait chronologies were significant in both species for the common period 1966–2015. Since RW and EW were highly correlated in all the cases (r > 0.86; P < 0.001), climate–growth relationships are only shown for EW and LW_adj (Figs. 4 and 5). The effect of temperature on larch EW was site-dependent, with a positive response to temperature in previous November in S2 and March and July of the current year in S3, whereas a negative effect of May temperature appeared in S4. In contrast, temperatures from March to May had a positive effect in LW_adj in all the sites except S1 (March), while a negative effect was seen for temperature in previous year October and current year June in S4, July in S1, and August in S3. The effect of temperature on pine growth was clearer (Fig. 4). April temperature positively affected EW, with an increase of the climatic signal from the warmer to the coldest sites. However, a negative effect was observed in previous December and current January temperatures in S2 and S3, respectively. The positive effect of temperature on LW_adj extended for two months (April and May) in all cases except S3, where temperature did not show a significant signal (Fig. 5).

Fig. 4

Correlations (Pearson’s coefficient) between earlywood (EW) and latewood adjusted (LW_adj) chronologies of Larix cajanderi with total precipitation and minimum temperature for the period 1966–2015. Correlations were calculated from August of the previous year (uppercase letters) to September of the current growth year (lowercase letters). Dashed lines indicated P < 0.05, and dotted lines indicate P < 0.01. Bar colors indicate the four sampled sites

Fig. 5

Correlations (Pearson’s coefficient) between earlywood (EW) and latewood adjusted (LW_adj) chronologies of Pinus sylvestris with total precipitation and minimum temperature for the period 1966–2015. Correlations were calculated from August of the previous year (uppercase letters) to September of the current growth year (lowercase letters). Dashed lines indicated P < 0.05, and dotted lines indicate P < 0.01. Bar colors indicate the four sampled sites

Precipitation also showed a tight correlation with the two tree-ring sectors, although patterns were less clear in larch than in pine. EW was favored by June precipitation in the S1, S3, and S4 sites, February precipitation in S2 (but with a negative effect in S3), and April precipitation in S1. May precipitation also negatively affected EW in S2. Previous September precipitation was positively correlated with LW_adj in S1 and S3. March, June, and July precipitation were correlated with higher LW_adj values in S2, S3, and S4, respectively. Previous December (S3), March (S1), and April (S2) precipitation showed a negative effect on LW_adj. In contrast, the effect of precipitation was clearer in pine, with May precipitation favoring growth in S2 and S3, and June precipitation favoring growth in S1. LW_adj responded to June (S1), July (S2), and August (S3) precipitation (Figs. 4 and 5).

In agreement with the significant impact of precipitation on residual chronologies, SPEI was significantly correlated in both species (Table 3). Larch EW residual chronologies showed moderate correlations with SPEI (r = 0.39–0.53; P < 0.01), with maximum intensity occurring in June in all sites except S2 (February), with a temporal scale ranging from 1 month in S2 to 12 months in S4. Pine EW residual chronologies were correlated with SPEI (r = 0.31–0.38; P < 0.05) with maximum signals occurring earlier than in larch (May) in all but S1 (June) with a temporal scale from 8 months in S3 to 15 months in S1 (Table 3). In both species, the strongest signals were found in the drier and warmer site S1. LW_adj larch chronologies also showed significant correlations with SPEI (r = 0.33–0.54; P < 0.01), with maximum signals occurring from March in S2 to July in S3 (Table 3). The highest signal intensity was found in the driest site S1 in June. Temporal scales ranged from 1 month in S1 to 11 months in S3. In Scots pine, correlations between SPEI and LW_adj chronologies were lower but still significant in all sites except S3 (r = 0.33–0.38; P < 0.01). The timing of the signal was earlier in S1 (January) than in S2 (July), with shorter temporal scales (1 to 2 months) than in larch sites.

Table 3 Correlation coefficients between residual chronologies and monthly standardized precipitation–evaporation index (SPEI) at different temporal scales for the period 1966–2015. Correlations were calculated from September of previous year to August of the current year. Values indicate maximum and minimum correlation indicating the month with highest signal (previous year in capital letters) and time lag in months (m)

3.4 Climate effect on pointer year occurrence

In each site, seven to nine years met the requirements to be considered pointer years (i.e., 1 SD above or below mean RW index chronologies), revealing differences between sites and species (see Annex Fig. 7). Although some pointer years were shared between sites, we also found unique pointer years (both narrow and wide rings) in each of the sites, highlighting local differences in the environmental conditions. The effect of precipitation triggering the formation of narrow and wide rings (Annex Fig. 8) showed a common pattern between sites and species (except in Larch S2), with wider rings occurring in years with higher precipitation events in the late spring and early summer (May–June), while narrow rings occurred in years with lower spring precipitation. Higher precipitation events that occurred later in the growing season (July–August) seem not to have affected the formation of wider rings, as shown by pine in sites S2 and S3.

4 Discussion

Temperature and radial growth of L. cajanderi and P. sylvestris have increased from 1966 to 2015 along the studied gradient in Central Yakutia (northeastern Siberia). Spring–summer temperatures were critical for tree radial growth, although radial growth was also affected by water availability, since EW and LW_adj were correlated with precipitation. In addition, the occurrence of wide rings was related to an earlier occurrence of precipitation during the pointer years.

Tree growth is strongly limited by low temperatures in the Yakutia territory (Tchebakova et al. 2006) due to the supercooling of the root area during the winter (Sugimoto et al. 2003; Nikolaev et al. 2009, 2011), which can extend to May and June. The onset of cambial activity requires a minimum temperature threshold (Rossi et al. 2008, 2013), which in the continuous permafrost zone is related to permafrost thawing (Vaganov et al. 1999; Bryukhanova et al. 2013). This effect was noticeable in pines, which showed greater earlywood growth under warm March–April conditions and more latewood growth when it was warm in April–May. Therefore, the observed warming in recent decades may have promoted an earlier permafrost thawing, triggering an earlier onset of growth and lengthening the total growing season (Nikolaev et al. 2011), resulting in higher growth rates and supporting the higher growth showed by the GAMM over time, after age effect being removed.

Temperature was not the only factor determining tree growth in the Yakutia territory. Precipitation showed high correlations with growth: pines responded to late spring rainfall, with later timing of signals in cooler sites. Larch also showed strong correlations with precipitation but without a clear pattern in the climate signal. This lack of a common signal may reflect the effect of local factors, such as topography on water availability, and may explain the lack of regional synchrony for this species. However, the use of SPEI revealed a common effect of spring water availability on ring width. This result was reinforced by pointer year analysis. Wide rings occurred in years with earlier precipitation events (May–June) in six out of seven sites. The seasonal thaw layer is a critical factor for water availability (Lopez et al. 2007), and any extra water supply from precipitation at the beginning of the growing season may enhance tree growth. In addition, liquid water from precipitation may also induce a deeper active soil layer as a result of heat transfer, and thus increase the nutrients available for tree growth. This result is in agreement with a previous experimental study, which showed that irrigation conducted in the permafrost area increased tree transpiration rates at the beginning of treatment, but this effect decreased later on when rainfall supplied the forest demand (Lopez et al. 2010). Thus, the association of earlier precipitation events with wide rings may be a consequence of the positive effect of extra water availability. We observed the same pattern in both species, meaning that although larch trees are better adapted to grow in the continuous permafrost zone, its radial growth may also respond to the earlier availability of liquid water. On the other hand, narrow rings occurred in years with later precipitation events (July), when the growing season would have been reduced in length before temperature again became the limiting factor (Yang et al. 2017).

The removal of the EW effect on LW reinforced the climatic signal in the last part of the ring, enhancing the latewood climate signal (Arzac et al. 2018b). Temperature had a positive effect on LW_adj in agreement with the high correlations found between maximum latewood density and summer temperature in higher latitudes (Briffa et al. 1998; Vaganov et al. 1999; Kirdyanov et al. 2008). Interestingly, precipitation also showed correlations with growth: pines responded to late spring rainfall, with later timing of the signals in cooler sites, and a month delay in the signal from EW to LW_adj. Again, the responses in larch were strong but showed idiosyncratic results for each site. Overall, these results indicate that water availability during summer may also limit growth at the end of the growing season in these cool environments.

Ongoing climate change will significantly impact the continuous permafrost zone in northern regions (Groisman et al. 2017). Warming is inducing an accelerated permafrost thawing (Groisman et al. 2013; Romanovsky et al. 2017), affecting water and carbon cycles (Groisman et al. 2017) as well as vegetation distribution (Tchebakova et al. 2006). Our results found significant growth increases in both species, despite previous growth projections in Yakutia territory predicting that an increase in P. sylvestris would be to the detriment of L. cajanderi distribution, due to the increase in the active soil layer depth which will negatively impact the growth of the better adapted to permafrost L. cajanderi (Tchebakova et al. 2006). Future projections should take water availability into account; although during warm and dry periods the permafrost may provide an extra water source for tree growth (Sugimoto et al. 2002; Saurer et al. 2016), if warming trends are persistent over time, melting of permafrost and increases in evapotranspiration may worsen the water balance during the growing season, potentially constraining radial growth (Groisman et al. 2013).

5 Conclusions

The study of environmental factors controlling tree growth and the formation of earlywood and latewood in conifers growing on permafrost in Central Yakutia revealed increased growth across sites in two species with different life strategies (i.e., evergreen and deciduous). The temperature was a critical growth factor in both conifers, but water availability also affected the radial growth. The formation of wide rings was triggered by earlier precipitation events in the two, highlighting the importance of soil water availability to the onset of the growing season at northern latitudes.

Annexes

Table 4 Mean tree-ring, earlywood, and latewood width for measured Larix cajanderi and Pinus sylvestris trees, and summary of dendrochronological statistics for residual chronologies for each site and parameter. The statistics shown correspond to the common period 1966–2015 in which climate analyses were performed

Table 5 Inter-site correlations (Pearson’s r) between Larix cajanderi and Pinus sylvestris RW chronologies.

Table 6 Inter-species correlations (Pearson’s r) between Larix cajanderi and Pinus sylvestris RW chronologies by species

Fig. 6

Effect of tree age modeled by the GAMM on tree tree-ring width (a and b, pine and larch respectively) and latewood width (c and d, pine and larch respectively).

Fig. 7

Pointer years in ring width (RW) chronologies (narrow and wide rings). Red lines show Larix cajanderi RW chronologies, and blue lines Pinus sylvestris RW chronologies. Pointer years were calculated as 1 SD above (wide rings) and below (narrow rings) average

Fig. 8

Effect of precipitation on the formation of pointer years. Lines with upward-pointing triangles and downward-pointing triangles indicate wide and narrow rings, respectively