- 1College of Forestry, Shanxi Agricultural University, Taigu, Shanxi, China
- 2Zhongtiaoshan Forestry Bureau of Shanxi Province, Houma, Shanxi, China
Introduction: Global climate change can affect the sensitivity of tree radial growth to climate factors, but the specific responses of tree radial growth to microclimate along the altitudinal gradient in the long term are still unclear.
Methods: In this study, the tree-ring width chronologies of Pinus tabuliformis Carr. in Shanxi Province of China were studied at three altitude gradients (1200-1300 m (low altitude), 1300-1400 m (medium altitude) and 1400-1500 m (high altitude)) during 1958-2017.
Results: The results showed that (1) the climate background could be divided into two periods based on the Mann-Kendall test analysis: 1958–1996 was a stable period (mean annual temperature (MAT)=10.25°C, mean annual precipitation (MAP)=614.39 mm), and 1997–2017 was a rapid change period (MAT=10.91°C, MAP=564.70 mm), indicating a warming and drying trend in the study region. (2) The radial growth of P. tabuliformis at different altitudes showed inconsistent variation patterns. The tree radial growth at low and medium altitudes (CV=27.01% for low altitude and CV=24.69% for medium altitude) showed larger variation amplitudes during the rapid change period than that in the stable period (CV=12.40% for low altitude and CV=18.42% for medium altitude). In contrast to the increasing trend, the tree radial growth rates at the high altitude showed a decreasing trend across years. (3) In the stable period, the radial growth of P. tabuliformis at the low altitude showed a significantly negative response to temperature and a positive response to precipitation in May and June. The tree radial growth at the medium altitude was positively related to precipitation in June and minimum temperature in February. The tree growth at the high altitude was mainly positively correlated with the temperature in May and August. In the rapid change period, the radial growth of P. tabuliformis at the low altitude was affected by more meteorological factors than that in the stable period. Medium-altitude trees were positively influenced by precipitation in June and minimum temperature in January, whereas high-altitude trees responded positively to wind speed in February. (4) Along altitudinal gradients, tree radial growth was more related to temperature than precipitation in the stable period. The tree radial growth at the high altitude during the rapid change period was only affected by wind speed in February, whereas the tree radial growth at low and medium altitudes was mainly affected by temperature to a similar extent during the two periods.
Discussion: The study indicated that tree growth-climate response models could help deeply understand the impact of climate change on tree growth adaptation and would be beneficial for developing sustainable management policies for forest ecosystems in the transition zone from warm-temperate to subtropical climates.
1 Introduction
The global warming trend continues (WMO, 2022), and forests will be greatly affected (Bonan, 2008). However, knowledge about how forests have responded to continuous warming since the latter half of the last century has not been well documented (Matías et al., 2017). Therefore, exploring the impact of climate change on forest ecosystems has become an important scientific issue (Babst et al., 2018). With accurate dating, high resolution, long time series and retention of rich tree growth information (Fritts, 1976), tree rings have been widely used by scholars from various countries for paleoclimate reconstruction, assessing the response of tree growth to climate change, and revealing the impact of temperature increase on regional vegetation (Zhang and Qiu, 2007; Liang et al., 2009; Fan et al., 2010; Wang et al., 2017).
Climate warming since the latter half of the 20th century has changed the environmental conditions and growth trends of trees (Rai et al., 2020), as well as the way trees interact with the climate (Babst et al., 2019). It also alters tree physiological processes and ecosystem productivity in an interactive and complex manner (Keenan et al., 2013; Reed et al., 2018). Studies have shown that the threshold effect leads to changes in the relationship between tree radial growth and climate (Jacoby and D’Arrigo, 1995). For example, under the influence of climate change, the changes in tree-ring width showed synchronization with climate factor fluctuations (Shestakova et al., 2016); however, temperature increases caused growth declines in some tropical forests (Feeley et al., 2007). This showed that trees in different climatic zones respond differently to climate change (Zhou et al., 2021). Altitudinal differences also increase the complexity of the response of tree growth to climate change (Zhang et al., 2022).
The altitude gradients are equivalent to natural cooling and humidifying platforms. With increasing altitude, the relationship between tree growth processes and climate factors changes (Salzer et al., 2014). Trees at different altitudes can adopt different ecological strategies to respond wisely to vertical zonal changes in mountain forest climates (Zhang et al., 2022). It is generally believed that the radial growth of trees at high altitudes is mainly limited by temperature, and at low altitudes, it is mainly related to precipitation (Wu et al., 2017). However, with global warming, the limiting effect of temperature on tree growth at high latitudes is gradually diminishing. The effect of precipitation on tree growth is becoming stronger (Büntgen et al., 2008; Loehle, 2009). This may be because there are different response mechanisms of tree radial growth to temperature and precipitation along altitudinal gradients. Therefore, it is important to explore how tree redial growth responds to climate change at different altitudes in the same region.
Currently, under the context of a warming and drying climate, the Qinling Mountains, an important geographic and climate dividing line in China, have attracted the interest of many scholars for their climate change and possible impacts (Li et al., 2022). Similarly, the Zhongtiao Mountains, which belong to the same warm-temperate to subtropical climate transition zone as the Qinling Mountains, are also a sensitive and vulnerable area for tree growth in response to climate change. In addition, P. tabuliformis is an important timber-tree species in North, Northwest and Northeast China, with a wide natural distribution area, which has a prominent strategic position in maintaining regional ecological security and significant economic advantage. Moreover, it is one of the dominant tree species in the Zhongtiao Mountains and the main forestation species in this area. Considering that tree rings have the characteristics of climatic sensitivity, accurate dating, strong continuity and high resolution, they have been widely used in dendroclimatological studies. For example, the tree rings of P. tabuliformis can be used to predict climate trends based on data about January-June precipitation changes over the last 296 years (1724-2019) (Cao et al., 2021) and the annual precipitation since 1686 in Ningwu County of China’s Shanxi Province (Li et al., 2006). However, the relationship between tree rings and micrometeorological factors at different altitudes is unclear. Therefore, it is necessary to explore the radial growth and dynamic response patterns of P. tabuliformis at different altitude gradients under climate change using the tree-ring method.
The aim of this study was to (1) investigate the pattern of radial growth of P. tabuliformis in response to climate factors at different altitude gradients; (2) quantify the contribution of climate factors to the radial growth of P. tabuliformis at different altitudes; and (3) reveal the radial growth trend of P. tabuliformis under the warming climate in the Zhongtiao Mountains. The study may provide sustainable management for forests with climate change in the Zhongtiao Mountains of China’s Shanxi Province.
2 Materials and methods
2.1 Study sites
Zhongtiao Mountain (34°36´-35°53´N and 115°15´-112°37´E) is located in southern Shanxi Province. It is a long, narrow, northeast−southwest trending mountain system, with the southern remnants of the Taihang Mountains on the east and the Yellow River cutting it off from the Huashan Mountains on the west. Lying in the transition zone from warm-temperate to subtropical climate zones, it belongs to a warm temperate continental monsoon climate. It is hot and rainy in summer and cold and dry in winter. According to the meteorological data of the nearest Qinshui County National Weather Station (35°41′ N, 112°11′ E), the mean annual temperature (MAT) was 3-11°C, with the coldest in January (-8.2°C) and the hottest in July (28.8°C). The mean annual precipitation (MAP) was 600-720 mm, mainly concentrated in June, July, August and September, accounting for 69.67% of the annual precipitation (Figure 1).
Figure 1 Climate records of monthly mean temperature (T), monthly mean minimum temperature (Tmmin), monthly mean maximum temperature (Tmmax), and monthly precipitation (P) during 1958-2017 at the Qinshui County Meteorological Station.
2.2 Tree-ring data and chronologies
In 2018, based on the distribution range of P. tabuliformis in the Zhongtiao Mountains, we selected 18 mature dominant trees as standard trees at low altitude (L, 1200-1300 m, n=7), medium altitude (M, 1300-1400 m, n=5) and high altitude (H, 1400-1500 m, n=6). The north−south direction and the location at breast height (1.3 m) were marked on the tree trunks. After felling, disks approximately 3-5 cm thick were intercepted at breast height (1.3 m) and marked with the standard tree number, north−south direction, etc. In the laboratory, the nonmarked surfaces were sanded and polished until the ring boundary of the disks was clearly visible. The disks were scanned by MICROTEK-3200A3L. Then, the ring width was measured using WinDENDEO (accuracy 0.001 mm) in eight directions. That is, we defined the due south direction as 0 degrees and recorded one core every 45 degrees in a clockwise direction, with a total of eight cores recorded for each disk. Meanwhile, we determined the accurate calendar age of each ring following a cross-dating procedure using the “dplR” package (Bunn, 2008) in R software (R Development Core Team, 2020). The age of each ring with dating errors was remeasured and cross-dated again. The unqualified cores that could not be dated (such as the severely damaged samples) were removed. Finally, 106 cores were used in this study.
The growth of trees is not only influenced by environmental factors but also controlled by their own physiology. To reduce the impact of individual differences between trees, the raw tree-ring width sequence was transformed into a standard and dimensionless ring-width index using the function “detrend” from the “dplR” package (Bunn, 2008) in R software. The standard chronology can accurately reflect the changes in tree growth patterns and the associated climatic factors (Wang and Yang, 2021). Additionally, the detrended tree-ring width indices were averaged using a robust biweight mean, with the function “chron” from the “dplR” package in R, to obtain the standard chronology. The quality of the chronologies was evaluated by several statistical parameters (Table 1). Specifically, mean sensitivity (MS) effectively reflects the interannual variability of the width between successive rings. For example, a high MS value indicates that the ring-width variation is more sensitive to climate than a low MS value. The signal-to-noise ratio (SNR) and expressed population signal (EPS) indicate the representativeness of the sample, with larger values indicating more environmental information for analysis in the chronology. For example, an EPS value above 0.85 generally indicates a qualified and reliable chronology (Wigley et al., 1984). The first-order autocorrelation (AC1) was generally used to identify whether the tree-ring width of the current year has a relationship with the climate factors of the previous year.
Table 1 Statistical parameters of the standardized chronology of P. tabuliformis at different altitude gradients.
2.3 Meteorological data
In this study, the daily meteorological data during 1958-2017 observed in the Qinshui County National Weather Station were used, as it is close to the study site and the data could represent the climate background in the study area. After calculation, we obtained climate indices, including monthly precipitation (P, mm), monthly temperature (i.e., mean temperature (T), maximum temperature (Tmax), minimum temperature (Tmin), mean maximum temperature (Tmmax), mean minimum temperature (Tmmin), °C), monthly relative humidity (RH, %), monthly sunshine hours (S, h), monthly average wind speed (WS, m/s), and standardized precipitation-evapotranspiration index (SPEI). Specifically, SPEI was calculated by the difference between standardized precipitation and potential evapotranspiration using the “SPEI” package (Beguería and Vicente-Serrano, 2009) in R software. For example, a negative SPEI value indicates a lack of water, which may depress tree growth. Conversely, a positive SPEI value indicates that tree growth can be promoted (Vicente-Serrano et al., 2010). Considering the “lag effect” of climate factors (Fritts, 1976), the climate data from the previous June to the current September during the whole period of 1958–2017 were used to investigate the effects of intra-annual climate factors on tree radial growth.
2.4 Data analyses
The Mann-Kendall test (Gocic and Trajkovic, 2013) was used to analyze the trend of meteorological data during 1958-2017. The correlation coefficient (r) of standard chronologies at different altitudes with monthly climate variables was calculated using the “treeclim” package (Bunn, 2008) in R software. Moreover, according to previous studies (Dormann et al., 2013; Harrison et al., 2018), |r| < 0.7 is an appropriate indicator that the tree’s radial growth at different altitudes is independent. To quantify the effects of climatic factors on tree radial growth, multivariate stepwise regression analysis was performed by the “step” function from the R package “stats” (Zhang et al., 2020) using the tree-ring width index as the predictor variable and climatic factors as explanatory variables to obtain the optimal regression model. Before regression, the values of climatic factors were standardized. To quantify the interpretation rate of tree radial growth by climatic factors, the R2 of each explanatory variable in the regression model was used as the explanatory rate (the sum of the explanatory rates of each factor was the total R2 of the regression model, which could actually reflect the explanatory amount of each factor for tree radial growth), and the percentage of the absolute value of the standardized regression coefficient of each explanatory variable to the sum of the absolute values of the regression coefficients of all explanatory variables was used as the contribution rate of climatic factors to tree radial growth (Le Bagousse-Pinguet et al., 2019). Data analysis was performed in R language, and the significance level was α=0.05.
3 Results
3.1 Climatic trends
According to the results of the Mann-Kendall test analysis of the meteorological data during 1958-2017 in the study area, the mean annual temperature (MAT) increased at a rate of approximately 0.15°C/decade, whereas the mean annual precipitation (MAP) decreased at a rate of 24.27 mm/decade (Figure 2). Additionally, the UF and UB curves of MAT and MAP had crossing points in 1996 and 1969, respectively (Figure 2). We considered 1996 to be an abrupt point in climate change because the UF value of MAT before 1996 was almost less than zero, and 1969 was the initial stage of the study.
Figure 2 The variation and abrupt change in (A) mean annual temperature and (B) mean annual precipitation during 1958–1996 (red circle) and 1997–2017 (blue square) based on the Mann-Kendall test. UF: sequence of forward test values, UB: sequence of reverse test values.
There was a decreasing trend for MAP in both periods (i.e., 1958-1996 and 1997-2017), which was similar to the MAP of the whole period (Figure 2B). Conversely, there were increasing trends for MAT in the period of 1997-2017 and the whole time, except for 1958-1996 (Figure 2A). Additionally, the MAT of 1997-2017 had a larger coefficient of variation (CV=10.70%) than that of 1958-1996 (CV=3.72%). Collectively, we divided the entire study period into relatively stable (1958-1996) and rapid change periods (1997-2017).
3.2 Chronology characteristics and tree radial growth along altitudinal gradients
In this study, the statistical parameters of the standard chronology had generally similar variation trends at three different altitude gradients (Table 1). For example, the values of first-order autocorrelation (0.826), signal-to-noise ratio (7.499) and expressed population signal (0.882) were the highest at medium altitude, followed by low altitude (0.559, 6.068 and 0.859, respectively), and finally high altitude (0.452, 5.972 and 0.857, respectively).
Comparing the characteristic curves of the ring width index (RWI) changes at different altitude gradients of the two periods, the RWI of low and medium altitudes increased slowly before 1996 (the slopes of the regression line were 9.34×10-4 and 2.64×10-3, respectively) and increased rapidly after 1996 (8.1×10-3 and 8.91×10-3, respectively), while the RWI of high altitude decreased at both periods (Figure 3).
Figure 3 Tree radial growth of P. tabuliformis in natural forests during 1958–1996 (red circle) and 1997–2017 (blue square) at (A) low altitude, (B) medium altitude and (C) high altitude.
3.3 Response of tree radial growth to climate at different altitude gradients
For the period 1958-1996, the tree radial growth of P. tabuliformis at the low altitude had a significantly negative correlation with monthly mean temperature in May (TC5, p < 0.05), maximum temperature in June (Tmax-C6, p < 0.05) and mean maximum temperature in May (Tmmax-C5, p < 0.01) but had a significantly positive correlation with monthly relative humidity in May (RHC5, p < 0.05) and June (RHC6, p < 0.05), as well as the SPEI in May (SPEIC5, p < 0.01) (Figure 4A). Consistent with low altitude, the tree radial growth at medium altitude had a significantly negative correlation with Tmax-C6 (p < 0.05) and a positive correlation with RHC6 (p < 0.05) (Figure 4B). However, the radial growth at high altitude had a significantly negative correlation with TC5 (p < 0.05), Tmmax-C5 (p < 0.05) and monthly minimum temperature in August of the previous year (Tmin-P8, p < 0.05) but had a significant positive correlation with TC8 (p < 0.01), Tmmax-C8 (p < 0.01), monthly maximum temperature in September of the previous year (Tmax-P9, p < 0.01), monthly mean precipitation in May (PC5, p < 0.05), as well as SPEI in May of the current year (SPEIC5) and September of the previous year (SPEIP9, p < 0.05) (Figure 4C).
Figure 4 Correlations between the ring width index (RWI) and monthly climatic factors at (A, D) low (L), (B, E) medium (M) and (C, F) high (H) altitudes during (A–C) 1958–1996 and (D–F) 1997–2017. *, p < 0.05; **, p < 0.01. RH, relative humidity; T, monthly mean temperature; P, monthly precipitation; SPEI, standardized precipitation-evapotranspiration index; WS, monthly mean wind speed; S, monthly mean sunshine hours; Tmax, monthly maximum temperature; Tmmax, monthly mean maximum temperature; Tmin, monthly minimum temperature; Tmmin, monthly mean minimum temperature. For the labels on the right, P indicates the previous year; C indicates the current year, and the numbers indicate the month.
Compared with the period of 1958-1996, the radial growth of P. tabuliformis at different altitudes showed more complex relationships with climatic factors during 1997-2017 (Figure 4). In the low-altitude zone, radial growth had a significantly negative correlation with monthly mean sunshine hours (SC5, p < 0.05) and TC5 (p < 0.01) but a positive correlation with RHC5 (p < 0.05), P-5 (p < 0.05), SPEIC5 (p < 0.05), Tmax-C5 (p < 0.05), Tmax-P7 (p < 0.01), Tmmin-C1 (p < 0.05) and wind speed in February of the current year (WDC2, p < 0.01).
Unlike at the low altitude, the radial growth at medium and high altitudes was significantly correlated with only a few meteorological factors (Figure 4). For example, the radial growth at medium altitudes had a significantly positive correlation with P-6, SPEIC6 and Tmmin-C6 but a negative correlation with Tmmin-P8 (p < 0.05, Figure 4E). Additionally, the radial growth at the high altitude only showed a significantly positive correlation with WSC2 (p < 0.01, Figure 4F).
3.4 Contribution of climatic factors to the radial growth of trees
The results of multivariate stepwise regression analysis showed that the influence of temperature and precipitation on the radial growth of P. tabuliformis varied with altitude (Table 2). Moreover, the contributions of temperature and precipitation to the radial growth of P. tabuliformis were different (Figure 5).
Table 2 Estimates of the multivariate stepwise regression model of the effect of climate factors on the radial growth of P. tabuliformis..
Figure 5 The contribution of different climatic factors in explaining the tree radial growth of P. tabuliformis at (A, D) low (L), (B, E) medium (M) and (C, F) high (H) altitudes during (A–C) 1958–1996 and (D–F) 1997–2017, respectively, based on the optimized regression model (see Table 2). See Figure 4 for abbreviations of meteorological factors.
During the period of 1958-1996, the current year’s temperature explained 67% (37% for Tmmax-C5 and 30% for Tmmin-C2), and precipitation explained 33% (RHC6) of the variation in radial growth of P. tabuliformis at low altitude, respectively (Figure 5A). For radial growth at the medium altitude, the contributions of temperature (Tmmin-C2) and precipitation (RHC6) in the current year were 49% and 51%, respectively (Figure 5B). At the high altitude, the temperature of the previous and current years explained 79% (34% for TC5, 21% for Tmax-C8 and 24% for Tmin-P8), and the precipitation of the previous year explained 21% (SPEIP9) of the radial growth variation in P. tabuliformis (Figure 5C).
During the period of 1997-2017, the contribution rate of temperature to the radial growth variation of P. tabuliformis at the low altitude was 66% (25% for Tmax-P7, 22% for Tmmin-C1 and 19% for Tmax-C4), and that of wind speed (WSC2) and sunshine hours (SC5) was 13% and 21%, respectively (Figure 5D). At the medium altitude, the contribution rates of temperature and precipitation to the variation in radial growth of P. tabuliformis were 53% (29% for Tmmin-P8 and 24% for Tmmin-C1) and 47% (SPEIC7), respectively (Figure 5E). At the high altitude, the radial growth variation of P. tabuliformis was only significantly affected by the wind speed in February of the current year (Figure 5F).
4 Discussion
4.1 Chronological characteristics and change patterns at different altitudes
In this study, most of the statistical parameters (e.g., SNR, EPS and AC1) of the standard chronology at the medium altitude exhibited higher values than those at low and high altitudes, but the value of mean sensitivity had the opposite trend (Table 1), which was consistent with the findings of previous studies (Fritts, 1976; Shrestha et al., 2017; Panthi et al., 2020). These results indicated that compared with that at low and high altitudes, the tree-ring width at the medium altitude exhibited lower sensitivity to climate change. This may be because the climatic conditions at medium altitudes were suitable, and other topographic or microenvironmental factors had an impact on tree growth (Weemstra et al., 2021).
Previous studies were conducted at sites in the arid zone, where spring precipitation may be the limiting factor for tree growth (Gou et al., 2004). However, our study area is in the transition zone from a warm temperate to subtropical climate zone, with relatively abundant precipitation, indicating that the growth of trees may be less affected by precipitation than other disturbance events. For example, our study area belongs to a warm temperature continental monsoonal climate zone, with cold and dry winters and windy springs. A previous study showed that the average daily wind speed in winter (pre-November to January) was approximately 3 m/s, which can directly cause the trunk and branches of the tree to break; furthermore, the wind also caused snow to accumulate, indirectly causing the breakage of the trunks and branches (Firm et al., 2009). Additionally, physiological drought caused by strong winds may also be a main cause affecting tree growth (Yu et al., 2016). In the present study, the radial growth of P. tabuliformis at high altitudes was more vulnerable to climate change than that at low and medium altitudes (Table 1), which was inconsistent with the previous findings by Zhang et al. (2015). This may be because the lowest elevation in Zhang et al.’s study was 2330 m, while the highest elevation in this study was 1500 m, indicating that tree radial growth at different altitudes was affected by microenvironmental factors.
4.2 Tree growth-climate response models at different altitudes
In this study, the radial growth of P. tabuliformis among the three altitudes had different degrees of variation and inconsistent relationships with climatic factors in the relatively stable (1958-1996) and rapid change periods (1997-2017) (Figure 3).
At the low altitude, the radial growth of P. tabuliformis was inhibited by the maximum temperature in May and June during the stable period (Table 2, Figure 4A), which was consistent with previous studies conducted in the eastern and northern Qinling Mountains (Shi et al., 2012; Sun and Liu, 2016). A study conducted in a lower montane site of the Swiss Plateau (805 m a.s.l.) showed that tree-ring growth had a significantly negative correlation with June temperature (Meyer and Braker, 2001). This may be because the higher temperature of the early growing season accelerated soil moisture loss, and tree roots could not obtain the water required for photosynthesis, leading to tree radial growth restriction (Peng et al., 2011). Furthermore, the radial growth of P. tabuliformis was enhanced by the relative humidity in May and June during the stable period (Table 2, Figure 4A). This is presumably because May and June are early in the growing season, when the water requirement for rapid shoot growth and dry matter accumulation increases (Yang et al., 2004). Meanwhile, the increase in relative humidity of the early growing season helps to replenish the water required for the physiological activities inside the trees and accelerates cell division, which leads to an increase in photosynthetic products and promotes the growth of the tree (Cui et al., 2021). Additionally, we found that tree radial growth at low altitudes was enhanced by the maximum temperature in April during the rapid change period (Table 2, Figure 4D). This phenomenon was also previously found in other regions of China (Yang et al., 2022). This is likely because the moderate warming in April was useful to break dormancy, promote snow melting and increase soil moisture, indirectly enhancing tree growth (Dang et al., 2007).
For the medium altitude, the radial growth of P. tabuliformis showed a significantly positive correlation with minimum temperature in February in the stable period (Table 2, Figure 4B). This was consistent with previous findings on the southern slope of the Qinling Mountains (Liu and Shao, 2003). The reason may be that February is the dormant period of plants, when the minimum temperature can affect the mechanical damage of trees and the reserve of tree water requirements (Liu et al., 2009). Conversely, the relationship between radial growth and the maximum temperature in June during the stable period was significantly negative (Figure 4B). For example, the increase in temperature of the growing season will result in the closure of stomata to reduce transpiration, thereby decreasing the photosynthetic rate and influencing organic matter synthesis and storage, which is detrimental to the growth of trees in the current year and even affects the growth of trees in the following year (Chen et al., 2011). In addition, the elevated summer temperature tended to induce drought events, which caused the embolism of xylem conduits and inhibited tree growth (Silva and Anand, 2013; Subedi and Sharma, 2013; Deslauriers et al., 2014). However, in the rapid change period, tree radial growth had a significant and negative correlation with the minimum temperature in August of the previous year (Figure 4E). A previous study showed that higher temperatures in the growing season caused stronger evapotranspiration, leading to soil water loss and stomatal closure, which in turn affected plant photosynthesis and inhibited the development of cells in the cambium (Tian et al., 2017).
Relative to the low and medium altitudes, tree radial growth at the high altitude during the stable period was affected by the temperature and precipitation of the current year as well as in the previous year (Figure 4C). For example, there was a significantly positive correlation between tree radial growth and the mean and maximum temperatures in August of the current year (Figure 4C). This may be related to the fact that higher temperatures can effectively prolong the length of the growing season in the high-altitude zone (Vitasse et al., 2018). Similarly, a previous study conducted at high altitudes in the Changbai Mountains showed that the higher temperature of the growing season was beneficial to the radial growth of trees (Wang et al., 2013). Moreover, there was a significantly positive correlation between tree radial growth at high altitudes and precipitation and SPEI in May during the stable period (Figure 4C). This may be because more precipitation in the early growing season helps to alleviate soil moisture deficits, increase the photosynthetic rate of trees, and facilitate carbohydrate synthesis and accumulation, thus promoting tree growth (Deslauriers et al., 2016; Chen et al., 2017). However, tree radial growth at the high altitude during the rapid change period was only influenced by the wind speed in February of the current year (Table 2, Figure 4F). One possible explanation is that winter precipitation in the study area mostly occurs in the form of snowfall, and low-temperature winds from the Siberian continent in February help the snowfall serve as an insulation layer for tree roots against frostbite (Reinmann and Templer, 2016).
4.3 Effect of warm and dry climate on the growth of P. tabuliformis
In our study, meteorological data recording showed an increase in temperature (from 10.25°C to 10.91°C) and a decrease in precipitation (from 614.39 mm to 564.70 mm) from 1958-1996 to 1997-2017 (Figure 2), leading to a trend toward a warmer and drier climate.
At the low altitude, the radial growth of P. tabuliformis in warm and dry climates (i.e., 1997-2017) had slightly lower mean values (0.97 vs. 1.00) and a higher coefficient of variation (CV, 27.01% vs. 12.40%) compared to the preclimate change period (i.e., 1958-1996) (Figure 3A). A previous study showed that the 600 mm precipitation isoline distinguished the response of tree ring width to climate, and the temperature south of this isoline significantly affected tree growth (Liu et al., 2019). This may be one reason for explaining our result, considering the area of our study located south of the 600 mm precipitation isoline.
Similar to the low altitude, the radial growth of P. tabuliformis at medium altitude in the period of 1997-2017 had higher CV (24.69% vs. 18.42%) and slightly higher mean values (0.99 vs. 0.97) than that of 1958-2016 (Figure 3B). Most likely, appropriate increases in temperature, within the threshold range of the optimal growth temperature of trees, can promote the absorption of water and nutrients, which in turn stimulates tree growth (Körner, 2007). As our study showed, the contribution of temperature for low (66% vs. 67%) and medium (49% vs. 53%) altitudes was almost equal over the two periods (Figure 5).
In contrast to the findings at medium and low altitudes, the radial growth of P. tabuliformis at the high altitude tended to decrease in both periods (Figure 3). This may be because high-altitude P. tabuliformis grows mainly in sandy and rocky areas with poor water storage capacity, where trees are susceptible to drought stress due to soil moisture being depleted, thus inhibiting tree radial growth (Girardin et al., 2016). Therefore, the negative effects of drought stress caused by rising temperatures on tree growth at the high altitude may be greater than the positive effects of increased precipitation with increasing altitude (Babst et al., 2019).
5 Conclusions
In this study, according to the results of the Mann-Kendall test, the climate background was divided into two periods: a relatively stable period (1958-1996) and a rapid change period (1997-2017). Changes in temperature (from 10.25°C to 10.91°C) and precipitation (from 614.39 mm to 564.70 mm) over the two periods indicated that the climate became warmer and drier in this region.
The radial growth of P. tabuliformis along the altitudinal gradients responded differently to the changing climate. Compared with the tree radial growth during the relatively stable period, the tree radial growth at low and medium altitudes increased faster with increasing temperature during the rapid change period. Conversely, a decreasing trend in radial growth at high altitudes was found throughout the whole period.
Additionally, the relationship between tree radial growth and climatic factors varied with altitude and varied in different climatic periods. For the trees at the low altitude, tree radial growth in the rapid change period was affected by more meteorological factors than that in the relatively stable period, although the temperature was the main influencing factor in both periods. Similar phenomena could be found at medium altitudes. However, tree radial growth at high altitudes during the rapid change period was only affected by wind speed in February, despite the temperature explaining 79% of the variation in radial growth in the relatively stable period.
Therefore, this study indicated that climate response models for tree radial growth along altitudinal gradients will be useful for providing more accurate suggestions for the future implementation of forest management policies in P. tabuliformis plantations under the background of global warming. Furthermore, it is helpful to better understand the dynamic responses of forest ecosystems to climate change, as well as to predict the impact of future global change on terrestrial ecosystems.
Data availability statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Author contributions
PN (First Author):Methodology, Software, Investigation, Formal Analysis, Writing-Original Draft; MZ: Investigation; TB: Investigation; BZ: Data Curation; LY: Investigation; SD: Writing, Editing; XY: Software; RG (Corresponding Author):Conceptualization, Supervision, Editing. All authors contributed to the article and approved the submitted version.
Funding
This project was supported by the National Natural Science Foundation of China (31400536).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
Babst, F., Bodesheim, P., Charney, N., Friend, A. D., Girardin, M. P., Klesse, S., et al. (2018). When tree rings go global: Challenges and opportunities for retro- and prospective insight. Quat. Sci. Rev. 197, 1–20. doi: 10.1016/j.quascirev.2018.07.009
Babst, F., Bouriaud, O., Poulter, B., Trouet, V., Girardin, M. P., Frank, D. C. (2019). Twentieth century redistribution in climatic drivers of global tree growth. Sci. Adv. 5, 1–10. doi: 10.1126/sciadv.aat4313
Beguería, S., Vicente-Serrano, S. M. (2009) SPEI: Calculation of the standardized precipitation evapotranspiration index. Available at: https://cran.r-project.org/package=SPEI.
Bonan, G. B. (2008). Forests and climate change: Forcings, feedbacks, and the climate benefits of forests. Science 320, 1444–1449. doi: 10.1126/SCIENCE.1155121
Bunn, A. G. (2008). A dendrochronology program library in r (dplR). Dendrochronologia. 26, 115–124. doi: 10.1016/j.dendro.2008.01.002
Büntgen, U., Frank, D., Wilson, R., Carrer, M., Urbinati, C., Esper, J. (2008). Testing for tree-ring divergence in the European Alps. Glob. Change Biol. 14, 2443–2453. doi: 10.1111/j.1365-2486.2008.01640.x
Cao, H. H., Zhao, X. E., Chen, F., Wang, S. J., Liu, X. H. (2021). Reconstructing January-June precipitation in southeastern shanxi over the past 296 years inferred from tree-ring records of pinus tabuliformis Carr. Chin. J. Appl. Ecol. 32, 3618–3626. doi: 10.13287/j.1001-9332.202110.018
Chen, L., Huang, J. G., Stadt, K. J., Comeau, P. G., Zhai, L., Dawson, A., et al. (2017). Drought explains variation in the radial growth of white spruce in western Canada. Agric. For. Meteorol. 233, 133–142. doi: 10.1016/j.agrformet.2016.11.012
Chen, L., Wu, S. H., Dai, E. F. (2011). Analysis of the tree-ring width chronologies of pinus koraiensis and larix olgensis on changbai mountains, northeast China. Geogr. Res. 30, 1147–1155.
Cui, J. Y., Peng, J. F., Li, J. R., Li, X., Peng, M., Yang, L. (2021). Responses of tree-ring width of pinus tabuliformis plantation to climatic factors in songshan mountains, central China. Chiness J. Appl. Ecol. 32, 3497–3504. doi: 10.13287/j.1001-9332.202110.002
Dang, H., Jiang, M., Zhang, Q., Zhang, Y. (2007). Growth responses of subalpine fir (Abies fargesii) to climate variability in the qinling mountain, China. For. Ecol. Manage. 240, 143–150. doi: 10.1016/j.foreco.2006.12.021
Deslauriers, A., Beaulieu, M., Balducci, L., Giovannelli, A., Gagnon, M. J., Rossi, S. (2014). Impact of warming and drought on carbon balance related to wood formation in black spruce. Ann. Bot. 114, 335–345. doi: 10.1093/aob/mcu111
Deslauriers, A., Huang, J. G., Balducci, L., Beaulieu, M., Rossi, S. (2016). The contribution of carbon and water in modulating wood formation in black spruce saplings. Plant Physiol. 170, 2072–2084. doi: 10.1104/pp.15.01525
Dormann, C. F., Elith, J., Bacher, S., Buchmann, C., Carl, G., Carré, G., et al. (2013). Collinearity: A review of methods to deal with it and a simulation study evaluating their performance. Ecography (Cop.). 36, 27–46. doi: 10.1111/j.1600-0587.2012.07348.x
Fan, Z. X., Bräuning, A., Tian, Q. H., Yang, B., Cao, K. F. (2010). Tree rings recorded may-august temperature variations since A.D. 1585 in the gaoligong mountains, southeastern Tibetan plateau. Paleogeogr. Paleoclimatol. Paleoecol. 296, 94–102. doi: 10.1016/j.paleo.2010.06.017
Feeley, K. J., Joseph Wright, S., Nur Supardi, M. N., Kassim, A. R., Davies, S. J. (2007). Decelerating growth in tropical forest trees. Ecol. Lett. 10, 461–469. doi: 10.1111/j.1461-0248.2007.01033.x
Firm, D., Nagel, T. A., Diaci, J. (2009). Disturbance history and dynamics of an old-growth mixed species mountain forest in the Slovenian Alps. For. Ecol. Manage. 257, 1893–1901. doi: 10.1016/j.foreco.2008.09.034
Girardin, M. P., Bouriaud, O., Hogg, E. H., Kurz, W., Zimmermann, N. E., Metsaranta, J. M., et al. (2016). No growth stimulation of canada’s boreal forest under half-century of combined warming and CO2 fertilization. Proc. Natl. Acad. Sci. U. S. A. 113, E8406–E8414. doi: 10.1073/pnas.1610156113
Gocic, M., Trajkovic, S. (2013). Analysis of changes in meteorological variables using Mann-Kendall and sen’s slope estimator statistical tests in Serbia. Glob. Planet. Change 100, 172–182. doi: 10.1016/j.gloplacha.2012.10.014
Gou, X. H., CHen, F. H., Yang, M. X., Peng, J. F., Qiang, W. Y., Chen, T. (2004). Analysis of the tree-ring width chronology of the qilian mountains at different elevations. Acta Ecol. Sin. 24, 172–176.
Harrison, X. A., Donaldson, L., Correa-Cano, M. E., Evans, J., Fisher, D. N., Goodwin, C. E. D., et al. (2018). A brief introduction to mixed effects modeling and multimodel inference in ecology. PeerJ 6, e4794. doi: 10.7717/peerj.4794
Jacoby, G. C., D’Arrigo, R. D. (1995). Tree ring width and density evidence of climatic and potential forest change in Alaska. Glob. Biogeochem. Cycles. 9, 227–234. doi: 10.1029/95GB00321
Keenan, T. F., Hollinger, D. Y., Bohrer, G., Dragoni, D., Munger, J. W., Schmid, H. P., et al. (2013). Increase in forest water-use efficiency as atmospheric carbon dioxide concentrations rise. Nature. 499, 324–327. doi: 10.1038/nature12291
Körner, C. (2006). “Significance of Temperature in Plant Life,” in Plant Growth and Climate Change, eds. Morison, J. I. L., Morecroft, M. D. (Blackwell Publishing Ltd), 48–69. doi: 10.1002/9780470988695.ch3
Le Bagousse-Pinguet, Y., Soliveres, S., Gross, N., Torices, R., Berdugo, M., Maestre, F. T. (2019). Phylogenetic, functional, and taxonomic richness have both positive and negative effects on ecosystem multifunctionality. Proc. Natl. Acad. Sci. U. S. A. 116, 8419–8424. doi: 10.1073/pnas.1815727116
Li, Q., Liu, Y., Cai, Q. F., Sun, J. Y., Yi, L. (2006). Reconstruction of annual precipitation since 1686 A.D. from NingWu region, shanxi province. Quat. Sci. 26, 999–1006. doi: 10.13287/j.1001-9332.202208.040
Li, Y., Qi, G., Feng, R., Liu, K. (2022). Responses of radial growth of pinus tabuliformis to climate at the northern of qinling mountains, China. Chines J. Appl. Ecol. 33, 2043–2050.
Liang, E. Y., Shao, X. M., Xu, Y. (2009). Tree-ring evidence of recent abnormal warming on the southeast Tibetan plateau. Theor. Appl. Climatol. 98, 9–18. doi: 10.1007/s00704-008-0085-6
Liu, Y., Linderholm, H. W., Song, H. M., Cai, Q. F., Tian, Q. H., Sun, J. Y., et al. (2009). Temperature variations recorded in pinus tabulaeformis tree rings from the southern and northern slopes of the central qinling mountains, central China. Boreas. 38, 285–291. doi: 10.1111/j.1502-3885.2008.00065.x
Liu, H. B., Shao, X. M. (2003). Reconstruction of January to April mean temperature at qingling mts from 1789 to 1992 using tree-ring chronologies. J. Appl. Meteorol. Sci. 14, 188–196.
Liu, Y., Song, H., Sun, C., Song, Y., Cai, Q., Liu, R., et al. (2019). The 600-mm precipitation isoline distinguishes tree-ring-width responses to climate in China. Natl. Sci. Rev. 6, 359–368. doi: 10.1093/nsr/nwy101
Loehle, C. (2009). A mathematical analysis of divergence problems in dendroclimatology. Clim. Change 94, 233–245. doi: 10.1007/s10584-008-9488-8
Matías, L., Linares, J. C., Sánchez-Miranda, Á., Jump, A. S. (2017). Contrasting growth forecasts across the geographical range of scots pine due to altitudinal and latitudinal differences in climatic sensitivity. Glob. Change Biol. 23, 4106–4116. doi: 10.1111/gcb.13627
Meyer, F. D., Braker, O. U. (2001). Climate response in dominant and suppressed spruce trees, picea abies (L.) karst., on a subalpine and lower montane site in Switzerland. Ecoscience. 8, 105–114. doi: 10.1080/11956860.2001.11682636
Panthi, S., Fan, Z. X., van der Sleen, P., Zuidema, P. A. (2020). Long-term physiological and growth responses of Himalayan fir to environmental change are mediated by mean climate. Glob. Change Biol. 26, 1778–1794. doi: 10.1111/gcb.14910
Peng, J. F., Yang, A. R., Tian, Q. H. (2011). Response of radial growth Chinese pine (Pinus tabulaeformis) to climate factors in wanxian mountain of he’nan province. Acta Ecol. Sin. 31, 5977–5983.
Rai, S., Dawadi, B., Wang, Y., Lu, X., Ru, H., Sigdel, S. R. (2020). Growth response of abies spectabilis to climate along an elevation gradient of the manang valley in the central Himalayas. J. For. Res. 31, 2245–2254. doi: 10.1007/s11676-019-01011-x
R Development Core Team (2020). R: A language and environment for statistical computing. Vienna, Austria: R Foundation for Statistical Computing, URL http://cran.r-project.org
Reed, C. C., Ballantyne, A. P., Cooper, L. A., Sala, A. (2018). Limited evidence for CO2-related growth enhancement in northern rocky mountain lodgepole pine populations across climate gradients. Glob. Change Biol. 24, 3922–3937. doi: 10.1111/gcb.14165
Reinmann, A. B., Templer, P. H. (2016). Reduced winter snowpack and greater soil frost reduce live root biomass and stimulate radial growth and stem respiration of red maple (Acer rubrum) trees in a mixed-hardwood forest. Ecosystems. 19, 129–141. doi: 10.1007/s10021-015-9923-4
Salzer, M. W., Larson, E. R., Bunn, A. G., Hughes, M. K. (2014). Changing climate response in near-treeline bristlecone pine with elevation and aspect. Environ. Res. Lett. 9, 114007. doi: 10.1088/1748-9326/9/11/114007
Shestakova, T. A., Gutiérrez, E., Kirdyanov, A. V., Camarero, J. J., Génova, M., Knorre, A. A., et al. (2016). Forests synchronize their growth in contrasting Eurasian regions in response to climate warming. Proc. Natl. Acad. Sci. U. S. A. 113, 662–667. doi: 10.1073/pnas.1514717113
Shi, J., Li, J., Cook, E. R., Zhang, X., Lu, H. (2012). Growth response of pinus tabulaeformis to climate along an elevation gradient in the eastern qinling mountains, central China. Clim. Res. 53, 157–167. doi: 10.3354/cr01098
Shrestha, K. B., Chhetri, P. K., Bista, R. (2017). Growth responses of abies spectabilis to climate variations along an elevational gradient in langtang national park in the central himalaya, Nepal. J. For. Res. 22, 274–281. doi: 10.1080/13416979.2017.1351508
Silva, L. C. R., Anand, M. (2013). Probing for the influence of atmospheric CO2 and climate change on forest ecosystems across biomes. Glob. Ecol. Biogeogr. 22, 83–92. doi: 10.1111/j.1466-8238.2012.00783.x
Subedi, N., Sharma, M. (2013). Climate-diameter growth relationships of black spruce and jack pine trees in boreal Ontario, Canada. Glob. Change Biol. 19, 505–516. doi: 10.1111/gcb.12033
Sun, C., Liu, Y. (2016). Climate response of tree radial growth at different timescales in the qinling mountains. PloS One 11, 1–15. doi: 10.1371/journal.pone.0160938
Tian, Q. Y., He, Z., Xiao, S., Peng, X. M., Ding, A. J., Lin, P. F. (2017). Response of stem radial growth of qinghai spruce (Picea crassifolia) to environmental factors in the qilian mountains of China. Dendrochronologia. 44, 76–83. doi: 10.1016/j.dendro.2017.04.001
Vicente-Serrano, S. M., Beguería, S., López-Moreno, J. I. (2010). A multiscalar drought index sensitive to global warming: The standardized precipitation evapotranspiration index. J. Clim. 23, 1696–1718. doi: 10.1175/2009JCLI2909.1
Vitasse, Y., Signarbieux, C., Fu, Y. H. (2018). Global warming leads to more uniform spring phenology across elevations. Proc. Natl. Acad. Sci. U. S. A. 115, 1004–1008. doi: 10.1073/pnas.1717342115
Wang, X., Yang, B. (2021). Divergent tree radial growth at alpine coniferous forest ecotone and corresponding responses to climate change in northwestern China. Ecol. Indic. 121, 107052. doi: 10.1016/j.ecolind.2020.107052
Wang, X., Zhang, M., Ji, Y., Li, Z., Li, M., Zhang, Y. (2017). Temperature signals in tree-ring width and divergent growth of Korean pine response to recent climate warming in northeast Asia. Trees. 31, 415–427. doi: 10.1007/s00468-015-1341-x
Wang, X., Zhao, X., Gao, L., Jiang, Q. (2013). Climatic response of betula ermanii tree-ring growth along an altitudinal gradient on the northern slope of the changbai mountains. Chin. J. Appl. Environ. Biol. 19, 929–934. doi: 10.3724/SP.J.1145.2013.00929
Weemstra, M., Freschet, G. T., Stokes, A., Roumet, C. (2021). Patterns in intraspecific variation in root traits are species-specific along an elevation gradient. Funct. Ecol. 35, 342–356. doi: 10.1111/1365-2435.13723
Wigley, T. M. L., Briffa, K. R., Jones, P. D. (1984). On the average value of correlated time series, with applications in dendroclimatology and hydrometeorology. J. Clim. Appl. Meteorol. 23, 201–213. doi: 10.1175/1520-0450(1984)023<0201:OTAVOC>2.0.CO;2
WMO (2022). Provisional state of the global climate in 2022. International Organization, Switzerland. https://public.wmo.int/en
Wu, Y., Wang, X., Ouyang, S., Xu, K., Hawkins, B. A., Sun, O. J. (2017). A test of BIOME-BGC with dendrochronology for forests along the altitudinal gradient of mt. changbai in northeast China. J. Plant Ecol. 10, 415–425. doi: 10.1093/jpe/rtw076
Yang, J. W., Liang, Z. S., Han, R. L., Wang, P. Z. (2004). Growth and water comsuption characteristics of Chinese pine under soil drought stress. J. Northwest A F Univ. (Natural Sci. Ed). 32, 88–92. doi: 10.3321/j.issn:1671-9387.2004.04.022
Yang, Y. R., Zhang, M. S., Zhang, L. N., Lu, Q. Q., Hong, Y. X., Liu, X. H. (2022). Different responses of radial growth of pinus tabuliformis to climate in the middle and western qinling mountains. Acta Ecol. Sin. 42, 1474–1486. doi: 10.5846/stxb202102060392
Yu, J., Xu, Q. Q., Liu, W. H., Luo, C. W., Yang, J. L., Li, J. Q., et al. (2016). Response of radial growth to climate change for larix olgensis along an altitudinal gradient on the eastern slope of changbai mountain, northeast China. Chin. J. Plant Ecol. 40, 24–35. doi: 10.17521/cjpe.2015.0216
Zhang, H., Cai, Q., Liu, Y. (2022). Altitudinal difference of growth-climate response models in the north subtropical forests of China. Dendrochronologia. 72, 125935. doi: 10.1016/j.dendro.2022.125935
Zhang, W. T., Jiang, Y., Wang, M. C., Zhang, L. N., Dong, M. Y. (2015). Responses of radial growth in larix principis-rupprechtii to climate change along an elevation gradient on the southern slope of luya mountain. Acta Ecol. Sin. 35, 6481–6488. doi: 10.5846/stxb201403140452
Zhang, X., Liu, S., Wang, J., Huang, Y., Freedman, Z., Fu, S., et al. (2020). Local community assembly mechanisms shape soil bacterial β diversity patterns along a latitudinal gradient. Nat. Commun. 11, 5428. doi: 10.1038/s41467-020-19228-4
Zhang, Q., Qiu, H. Y. (2007). A millennium-long tree-ring chronology of Sabina przewalskii on northeastern qinghai-Tibetan plateau. Dendrochronologia. 24, 91–95. doi: 10.1016/j.dendro.2006.10.009
Keywords: climate change, radial growth, altitude gradient, Pinus tabuliformis Carr., Zhongtiao Mountains
Citation: Ning P, Zhang M, Bai T, Zhang B, Yang L, Dang S, Yang X and Gao R (2023) Dendroclimatic response of Pinus tabuliformis Carr. along an altitudinal gradient in the warm temperate region of China. Front. Plant Sci. 14:1147229. doi: 10.3389/fpls.2023.1147229
Received: 18 January 2023; Accepted: 21 March 2023;
Published: 30 March 2023.
Edited by:
Hongbo Li, Chinese Academy of Agricultural Sciences (CAAS), ChinaReviewed by:
Mingjie Xu, Shenyang Agricultural University, ChinaYe Yuan, Anhui Normal University, China
Copyright © 2023 Ning, Zhang, Bai, Zhang, Yang, Dang, Yang and Gao. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Runmei Gao, sxndgrm@163.com