Dynamic response to climate change in the radial growth of Picea schrenkiana in western Tien Shan, China

Liang Jiao · Shengjie Wang · Ke Chen · Xiaoping Liu

Abstract Forests are important ecosystems for economic and social development.However, the response of tree radial growth to climate has produced ‘divergent problems’ at high latitudes under global warming.In this study, the response stability and trend of Picea schrenkiana radial growth to variability in climate factors were analyzed in the mid-latitudes of the western Tien Shan Mountains.Radial growth of P.schrenkiana was mainly limited by minimum and mean temperatures.The divergent responses of radial growth occurred in response to the minimum and mean temperatures at the beginning of the growing season (April-May) of the current year, but responses to drought occurred in July-September of the previous year.And the mean and minimum temperatures in June-September of the current year were both stable.Radial growth first increased and then decreased according to the basal area increment, with a gradual increase in temperature.Therefore, forest ecosystems in mountainous arid areas will be increasingly affected by future climate warming.

Keywords Dendroecology · Divergent response · Climate change · Forest degradation · Western Tien Shan mountains

Introduction

Forests are important natural resources covering approximately 3.45 billion ha or 30% of the global land area.Forest play important roles in the formation and maintenance of biodiversity, carbon sequestration, the carbon balance, the hydrological cycles and paleoclimatic reconstructions (FAO 2006; Anderegg et al.2020).Global warming has a significant effect on the sensitivity of tree species and the composition of forest ecosystems (Williams et al.2010).The Intergovernmental Panel on Climate Change (IPCC) fifth report reported that the global mean temperature at the Earth’s surface increased by 0.85 °C ± 0.20 °C between 1880 and 2012, the hottest 30-year-period in the past 800 years, and the future climate was forecast to warm further (IPCC 2013).Therefore, the sustainable development of forest ecosystems has been integrated with the core global ecological problems as global warming continues to intensify (Jiao et al.2017).

As a database of climate change, tree rings are widely selected to study regional historical climate reconstructions and growth-climate relationships due to accurate dating, strong continuity, high sensitivity and easy access (Wang et al.2020).However, several scholars have found that climate warming led to a decrease in the sensitivity of trees to temperatures at high latitudes in the Northern Hemisphere (Wang et al.2017; Zhang et al.2019).The phenomenon of an unstable response between radial growth and climatic factors was termed the ‘divergence problem’, which challenged the accuracy and scientificity of tree-ring climate reconstruction work based on the ‘uniformitarian principle’ (D’Arrigo et al.1992; Briffa et al.1998; Wilmking et al.2020).At present, the results of divergent responses differ between regions, tree species, elevation gradients and age groups.For example, Siberian Larch (Larix sibiricaLedeb.) did not show a distinct divergent response to temperature in the Altai Mountains, but divergent response ofLarix sibiricaand the relatively stable responses by Schrenk spruce (Picea schrenkianaFisch.& C.A.Mey.) to climate factors were shown in the eastern Tien Shan (Shang et al.2010; Jiao et al.2019).Global warming is believed to have positive effects on forests at high elevations but negative effects at low elevations reflected by radial growth to climate along elevation gradients according to dendroclimatology, as observed in the western mountains of North America (Salzer et al.2009).Radial growth response of Chinese white pine (Pinus armandiiFranch.) of different age groups to climate factors were different in the Baiyunshan National Nature Reserve, central China (Peng et al.2018).Moreover, the mechanisms of the ‘divergence problem’ are also controversial, mainly focusing on the heterogeneous age effects of biological species, the threshold changes of the growth responses to climate, inappropriate detrended methods, and low-quality climate data (Cook et al.1997; D’Arrigo et al.2008; Allen et al.2018; Wilmking et al.2020).Therefore, it is necessary to further assess the stability of the responses to climate changes by different tree species and in diverse areas, and to explore the causes of the divergent responses.

The Tien Shan Mountains or simply Tien Shan, as a wet island in an arid region, play an important role in maintaining the stability of oases and the sustainability of water resources in northwest China (Liang et al.2016).Moreover, the Tien Shan adequately reflects the fluctuations in regional climate change as there is limited interference from areas of human activities.Thus, the Tien Shan has become a “hot spot” for dendrochronological research (Zhang et al.2016a).Although numerous studies on paleoclimate reconstruction have been performed in the Tien Shan since the 1960s, few have focused on the ‘divergence problem’ (Zhang et al.2016b).At the same time, climate transformations from warm-dry to warm-wet areas since the end of the 1980s also provide a possibility for further analysis of the impacts of climate change on forest ecosystems (Shi et al.2007).However, conclusions of climatic impacts on trees and forests have been inconsistent.For example, the basal area increment (BAI) ofP.schrenkianahas shown n an increase in the Tien Shan with the extension of the growing season caused by global warming (Qi et al.2015).However, other studies found thatP.schrenkianagrowth declined trend with intensified water deficits in inland Asia (Liu et al.2013).Therefore, the impacts of climate change on ecosystems in the Tien Shan need further analysis.

A virginP.schrenkianaforest in the western Tien Shan was selected and the stability of the responses to climate changes were analyzed and changes to the growth patterns of dominant coniferous trees with climate transformation were described.The research questions were: (1) What are the principal climatic factors limiting the growth ofP.schrenkianain the western Tien Shan? (2) Is there a ‘divergence problem’ in the dynamic growth-climate relationship? and (3) What was the trend in growth ofP.schrenkianaunder climate change?

Materials and methods

Study site and tree ring acquisition

Tree rings were collected from species in theP.schrenkianaNational Nature Reserve in western Tien Shan (XTS) (Fig.1).The research region has a typical temperate continental climate with rainy and hot conditions in the same period.However, warm, humid air from the Atlantic enters the western Tien Shan under the influence of the westerly winds, resulting in relatively warmer and more humid climates than eastern Tien Shan.Annual total precipitation is approximately 508 mm, mainly concentrated in April-September, accounting for 84.0% of the annual rainfall.Annual mean temperature is 3.4 °C, with January the coldest month at-11.2 °C, and July the hottest month at 15.1 °C (Figs.2, 3).The annual mean temperature in the study region has been increasing since 1959 (mean maximum temperature: 0.19 °C per decade, mean temperature: 0.33 °C per decade, mean minimum temperature: 0.53 °C per decade,P＜ 0.001).The interannual variation in precipitation was not significant (3.82 mm per decade,P=0.574).We further studied the rate of temperature increase (slope/10a) and found that it had a higher rate of increase in the 1980s and the decade before the twenty-first century (P＜ 0.01).The above meteorological data are from Zhaosu (ZS) Meteorological station.P.schrenkianais the dominant conifer species in the region and mainly distributed at an altitude of 1500-2800 m a.s.l.with a canopy density of approximately 40%.

Fig.1 Locations of the sampling site and the nearest meteorological station (Western Tien Shan)

To minimize the effects of non-climate factors, cores from healthy trees in undisturbed open stands with similar slopes (10°-30°), aspects (north-facing) and soil conditions (grey cinnamon soil) were selected.Two cores from each tree were taken at the 120° position with 5-mm increment borers.Sixty cores from 30 live trees were obtained in June 2017.The characteristics of the forest communities are shown (Table 1).

Fig.2 Monthly variation of mean temperature and total precipitation (1959-2016)

Tree-ring chronology development

In the laboratory, the cores were mounted on woody supports, air-dried, sanded to produce clear boundaries for identification, and the widths measured using LINTAB (RINNTECH e.K.Heidelberg, Germany) with 0.001 mm resolution (Stokes and Smiley 1968; Schweingruber 1996).Two cores with ambiguous rings and large growth-curve variations were removed.The accuracy of all measurement series was verified by the COFFECHA program (Holmes 1983).Ring widths were detrended using negative exponential or straight line curves, and single chronologies were combined with the individual series by calculating the biweight robust means chronology using the ARSTAN program (Cook et al.2013).The variance in chronologies was stabilized with the Briffa Rbar-weighted method (Cook and Kairiukstis 1990).Finally, a standard chronology (STD) was developed.

Statistical parameters of the chronology reflect the strength and quality of the climate information (Fritts 1976).Mean sensitivity (MS) and standard deviation (SD) were used to evaluate chronology sensitivity to climate factors.The interpretation of the variance in the first principal component (PC1) and signal-to-noise ratio (SNR) characterized the common climate information contained in the chronology.The average correlation coefficient (R1, R2, R3) indicated the consistency of the radial growth fluctuations.An expressed population signal (EPS) greater than 0.85 indicated that the signal contained in the tree-ring series represented the overall characteristics (Wigley et al.1984).A high first-order autocorrelation coefficient (AC1) indicated a high growth response in the current year to the climate in the previous year (Fig.4; Table 2).

Table 2 Comparison of the statistical characteristic parameters of tree-ring width chronologies in study area with those in the near study areas

Fig.4 Interannual variability of the tree-ring standard chronology, its running expressed population signal (EPS), sample depth, and mean inter-series correlation (Rbar)

Table 1 The information of tree-ring sampling site

Meteorological data

Meteorological data were acquired from the ZS meteorological station at 43.15°N, 83.13°E, and at 1854 m (Fig.1).The climatic data included monthly maximum, mean, and minimum temperatures and total precipitation (1959-2016).These temperature data were adjusted for the altitude differences between the sampling sites and the meteorologicalstation (the temperature decreases by 0.6 °C for every 100 m increase).Radial growth exhibited a delayed response to climate, showing that trees are affected not only by the current year’s climate but also by the previous year’s climate.Therefore, the monthly climate data were selected from the previous July to the current September for the analysis of the growth-climate relationship (Fritts 1976).

The standardized precipitation and evapotranspiration index (SPEI) were also calculated using monthly climate data (Vicente-Serrano et al.2010).The SPEI is an ideal tool for monitoring and studying the process of drought because it organically integrates the Palmer drought severity index (PDSI) with evaporation demand changes and standardized precipitation index (SPI) properties over multiple time scales (Diego Galván et al.2015).The SPEI data on a 24-month scale were calculated using R (R Core Team R 2014) with the ‘SPEI’ R package (Beguería and Vincente-Serrano 2013).

Statistical methods

Climatic factors controlling radial growth were determined by calculating the Pearson correlation coefficient between monthly climate data and chronology with SPSS 22.0 (Kirkpatrick and Feeny 2014).Two significance levels (0.05 and 0.01) were set and a two-tailed test was performed.

whereYis the ecological security index,Xitheilandscape pattern index,C1(Y,Xi) the Pearson correlation coefficient betweenYandXi,C2(Y,Xi) the Pearson correlation coefficient betweenYand lnXi,C(Y,Xi)the Pearson correlation coefficient betweenYandXi.The sign (positive or negative) of the coefficient is consistent with that of the original value.

The moving correlation analysis was used to test the dynamic responses of radial growth to climate using the Dendroclim 2002 program (Biondi and Waikul 2004).We employed a fixed 30-year period starting from 1960 that slid across time to compute the correlation coefficients between tree-ring width chronologies and seasonal climatic factors with a bootstrap procedure using 1000 replications.

Radial growth pattern

Basal area index (BAI) is a more realistic indicator than tree-ring chronology to quantify trends in tree radial growth.In fact, the BAI overcomes undetectable increasing growth trends caused by a decrease in ring width with age in mature trees (Weiner and Thomas 2001; Gomez-Guerrero et al.2013).The BAI eliminates any data-transformation biases resulting from the use of the original tree-ring width data; in particular, downward trends of BAI indicate that growth was restricted by the environment (Schuster and Oberhuber 2013).Therefore, BAI is considered a better means for assessing radial growth (LeBlanc 1990; Pedersen 1998).The formula is:

whereRis the tree radius in yearst-1, and TRW is the raw width of the tree-ring in yeart.

Results

Correlation between chronology and climatic factors

The standard chronology (STD) forP.schrenkianacontained extensive climate information of good quality and sensitive to climate change.Thus, it was suitable for the analysis of radial growth-climate relationships (Table 2).Correlation analysis indicated that temperature was the controlling factor for tree growth (Fig.5).On a monthly scale, the STD chronology was positively correlated to total precipitation (previous July:r=0.382; previous September:r=0.290), with a high level of significance (P＜ 0.01,P＜ 0.05), negatively correlated to temperatures (maximum, mean and minimum temperature in the previous July to October and the current April to September), with a high level of significance (P＜ 0.05), and positively correlated to the SPEI (previous July:r=0.399) with a high level of significance (P＜ 0.01).On the growing season scale, STD was correlated to the maximum temperature (r=-0.669), mean temperature (r=-0.813), minimum temperature (r=-0.806), total precipitation (r=0.413) and SPEI (r=0.361) with a high level of significance (P＜ 0.01) late in the growing season of the previous year (July-September).The STD was correlated to the maximum, mean and minimum temperatures in the early growing season (April-May) and the middle-late growing season (June-September) of the current year with a high level of significance (P＜ 0.01,P＜ 0.05).

Fig.5 Correlation between chronology and climatic factors

Dynamic growth-climate relationship

Based on the results of the growth-climate relationships, the moving correlations of radial growth and climatic factors in the growing season were tested during 1960-2016 (Fig.6).The tree-ring width chronology to mean and minimum temperatures, precipitation and SPEI were stable during the middle-late growing season of the previous year (July-September), and the significant negative correlation between the STD and mean and minimum temperatures was also stable during the middle-late growth season of the current year (June-September).However, the divergent responses of the STD to mean and minimum temperatures in the early growing season of the current year showed gradually increasing trends (mean temperature:y=-0.015x+ 30.158,R2=0.693,P＜ 0.001; minimum temperature:y=-0.013x+ 26.553,R2=0.692,P＜ 0.001).

Fig.6 Moving correlations between chronology and climatic factors; (a-e represents the mean maximum temperature, mean temperature, mean minimum temperature, total precipitation and SPEI, T temperature, p previous year, c current year)

Radial growth trend

Interannual variation in the BAI was evaluated during the period 1816 to 2016 (Fig.7).The radial growth trended to increase and then decrease, specifically an upward trend from 1816 to 1958 (1.05 cm2per decade,P＜ 0.001), but downward trend from 1959 to 2016 (-2.16 cm2per decade,P＜ 0.001).To avoid influence of the age effect, the trend of radial growth of trees of different age classes (＞ 200 years old; ＜ 200 years old) was analyzed.The results show that the radial growth of the two age classes trended upwards (1.60 and 2.22 cm2per decade, respectively) from 1816 to 1958, and downward from 1959 to 2016 (-3.01 and-1.60 cm2per decade, respectively,P＜ 0.001).

Fig.7 Interannual trends of BAI from 1816 to 2016

Discussion

Ecological response of radial growth to the main limiting climatic factors

To verify that our tree-ring chronology contained sufficient climate information, we compared MS, SD, EPS and SNR with those near the research area and found that our chronology had a high quality and could meet the requirements of growth-climate research (Cook and Kairiukstis 1990).Tree rings not only record climate information but also represent the physiological response of vegetation to climate change (Liang et al.2009).P.schrenkianagrowth in western Tien Shan is mainly restricted by previous drought conditions and current temperatures in the growing season.

Due to the influence of elevation-dependent responses in plants, trees at high altitudes are restricted by low temperatures, a positive correlation with temperature.Trees at low altitudes are mainly restricted by drought (Shao et al.2005).In our study, STD was negatively correlated with temperature and positively with precipitation and the SPEI from the previous July to September period, indicating that drought in the middle to late growing seasons of previous years had a profound impact on tree growth (Fig.5).This is because the water storage capacity was poor forP.schrenkianadistributed on steep slopes with shallow roots at low elevations, so the limited water barely met the needs of trees during the growing season (Table 1).Moreover, high temperatures increased the transpiration of leaves and evaporation of surface moisture, increasing the scarcity of available water resources and reducing the accumulation of photosynthetic nutrients, easily resulting in the formation of narrow tree rings (King et al.2013; Gou et al.2015; Zhang et al.2017).There was a significant negative response reflected in tree-ring width to temperature in the growing season of the current year, especially the mean minimum temperature (Fig.5).This indicated that radial growth ofP.schrenkianawas mainly limited by the mean minimum temperature.The same result has been found in the Tien Shan and other regions (Wu et al.2013; Liu et al.2016).As global warming further intensifies, temperatures in the western Tien Shan will show a significant upward trend, with the mean minimum temperature rising at the highest rate (Fig.3).In fact, the mean minimum temperature effects plant respiration, net productivity, nutrient absorption and other physiological processes (Palombo et al.2014).Therefore, a significant increase in the mean minimum temperature would result in a decrease in the diurnal temperature differences during the growing season and an increase in respiration and photosynthetic consumption, thereby restricting the differentiation and growth ofP.schrenkianacambium cells (Yu et al.2013).

Fig.3 Interannual variations of temperature (mean maximum, mean, and mean minimum) and total precipitation (1959-2016)

Response stability of radial growth to climate change in time

The results of growth-climate dynamic relationships exhibit complex characteristics due to differences in habitat, tree species, age, and climate change patterns (D’Arrigo et al.2008; Leonelli et al.2009; Fonti et al.2018; Jiao et al.2019).In our study,P.schrenkianawas affected by the persistence of drought conditions during the mid to late growing season of the previous year (Fig.5).Temperature increase resulted in increased evaporation which increased the adsorption of water by leaf mesophyll cells and slowed the movement of water, resulting in closing of the stomata and weakening photosynthesis (Berninger et al.2004).In addition,decreased precipitation did not satisfy tree water demands, so the trees were under water shortage stress for considerable time (Jiang et al.2017).Therefore, radial growth exhibited consistent responses to the temperature, precipitation and SPEI of the previous year.A similar study took place in the mid-western Tien Shan, showing that Schrenk spruce exhibited positive responses to the standardized precipitation and evapotranspiration index (SPEI) from July-October of the previous growing season (Chen et al.2016).

The limiting effects of the mean and minimum temperatures on annual rings were also relatively stable during the mid to late growing season of the current year (Fig.6).This was because the relatively high amount of precipitation would address the demands of tree growth in this period, so the radial growth exhibited no correlation with precipitation or drought.However, the continuously high temperatures during this period increased tree respiration and carbohydrates consumption, leading to the maintenance of radial growth sensitivity of radial growth to temperature (Schaberg et al.2000).Therefore, the ring chronology showed a stable, negative response to temperature in the mid to late growing season of the current year.

The response by Schrenk spruce showed divergence to mean and minimum temperatures at the beginning of the growing season of the current year, mainly supporting the range from uncorrelated to significantly negatively correlated (Fig.6).Carbohydrates produced by photosynthesis are a major source of energy and nutrients for tree growth.When increased temperatures were still within the physiological threshold range of Schrenk spruce, relatively high temperatures increase photosynthetic efficiency and promote radial growth (Jiang et al.2014).However, Franceschini et al.(2013) found that the physiological adaptive threshold of Norway spruce (Picea abiesFranch.) was gradually exceeded with an abrupt increase in mean temperature in the early growing season, and the species showed a divergent response to the effect of temperature as reflected by the annual rings changing from promotion to reduction.At the same time, the increase in precipitation from rain and melting snow increased soil moisture, largely satisfying moisture demands for radial growth in the early growing season (Dang et al.2007).Therefore, species sensitivity to total precipitation did not change significantly in the early growing season.The same result was also found in a study on the divergent response ofLarix sibiricaLedeb.and climatic factors at high altitudes of the eastern Tien Shan (Jiao et al.2015).

Radial growth under climate change

The effects of global climate change on tree growth at different temporal and spatial scales remain highly uncertain (IPCC 2013).Our research results show BAI increased initially and then decrease (Fig.7).In particular, BAI of trees ＞ 200 years old and ＜ 200 years old showed the same growth patterns, indicating that tree-ring width was mainly controlled by climate and is not related to tree age.

Species are often limited by primary ecological factors, and a change beyond the threshold of ecological adaptation will lead to growth reduction and often death (Wiens 1989).At the same time, studies of tree physiology have also concluded that the ‘threshold effect’ occurs, supporting the inverted U-shaped relationships of tree growth to climate (Kramer and Kozlowski 1979).This was also verified by simulating the effect of temperature rise on the size and number of stem cells of trees along the northern edge of a forest (Tuovinen 2005).Precipitation in the western parts of the Tien Shan is more abundant than in the eastern.The length of the growing season is prolonged and the rate of cambium cell division is accelerated with temperature increase.Therefore, radial growth was promoted in the western Tien Shan.However, with further increases in temperature, the tree-ring width chronology showed a significant negative response, and the BAI also showed a significant downward trend, indicating temperatures rose beyond the optimal threshold value in recent decades, and that the effects of temperature on growth have not been restricted with increased respiration, worsening drought and decreased photosynthetic nutrient accumulation.The same result was also found in the eastern Tien Shan (Jiao et al.2019).

Drought due to climate warming has led to an increase in tree mortality, a decline in forest productivity and the instability of forest ecosystems worldwide (Sarris et al.2010; Anderegg et al.2012), especially in arid inland Asia (Poulter et al.2013; Anderegg et al.2015).Forest decline caused by climate warming represents a significant threat to forest resources and, to a certain extent, ecological problems such as substantial reductions in available water resources, adverse changes to forest carbon storage, and weakened forest climate regulation (Liang et al.2006).In addition, increasing temperatures will not only cause a decline in forests in arid regions but also impact humid area and highaltitude forests in the future (Jiao et al.2015).

Conclusion

In the mid-latitude western Tien Shan of the Northern Hemisphere, the minimum and mean temperatures in the growing season are the main factors limiting the radial growth ofP.schrenkiana.Moreover, the dynamic responses of radial growth to the different climate factors differ, showing relatively stable response to drought in the mid to late growing season of the previous year and to the minimum and mean temperatures in same period of the middle-late growing seasons of the current year.However, responses to minimum and mean temperatures in the early growing season of the current year diverged.In addition, radial growth of Schrenk spruce increased and then decreased with climate warming.If temperatures continue to rise in the future, radial growth of Schrenk spruce will become more restricted.Therefore, the influence of climatic change on mountain forest ecosystems in arid regions is important.According to the regular responses by radial growth to climate factors, the monitoring of dominant tree species, targeted forest conservation management and scientific implementation of forest resource protection measures should be strengthened to mitigate the threat of future climate warming.