Next Article in Journal
Conservation and Restoration of Mangroves in Response to Invasion of Spartina alterniflora Based on the MaxEnt Model: A Case Study in China
Previous Article in Journal
Effect of Samples Length on the Characteristics of Moisture Transfer and Shrinkage of Eucalyptus urophylla Wood during Conventional Drying
Previous Article in Special Issue
NDVI-Based Vegetation Dynamics and Response to Climate Changes and Human Activities in Guizhou Province, China
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Radial Growth of Juniperus squamata Showed Sharp Increase in Response to Climate Warming on the Three-River Headwaters Region of Tibetan Plateau since the Early 21st Century

1
College of Water and Soil Conservation, Beijing Forestry University, Beijing 100083, China
2
Ji County Station, Chinese National Ecosystem Research Network (CNERN), Beijing Forestry University, Beijing 100083, China
3
State Key Laboratory of Urban and Regional Ecology, Research Center for Eco-Environmental Sciences, Chinese Academy of Sciences, Beijing 100085, China
4
State Key Laboratory of Desert and Oasis Ecology, Xinjiang Institute of Ecology and Geography, Chinese Academy of Sciences, Urumqi 830011, China
*
Author to whom correspondence should be addressed.
Forests 2023, 14(6), 1219; https://doi.org/10.3390/f14061219
Submission received: 3 April 2023 / Revised: 12 May 2023 / Accepted: 2 June 2023 / Published: 12 June 2023
(This article belongs to the Special Issue Responses of Forest Ecosystems to Climate Change)

Abstract

:
In order to explore the impact of climate change on the ecosystem at high altitudes, dendroclimatology was used to study the response of radial growth of Juniperus squamata Buch.-Ham. ex D.Don to the rapid warming in Nangqian County over the past 60 years, and a tree-ring width chronology for 115 years was established. (1) Meteorological data showed that the temperature in Nangqian County of the Tibetan Plateau has increased continuously during the past 60 years, and the minimum temperature has had the most significant change (0.63 °C/10a), especially between 2000–2019 (0.12 °C/a). Over the same time period precipitation has not changed significantly (0.94 mm/a, p > 0.10). The standard chronology was used to reconstruct the mean temperature series from July to September in Nangqian meteorological station during the past 115 years (1905–2019). The explained variance of the reconstructed equation was 42.8% (40.8%, after adjusting for degrees of freedom). The reconstructed temperature series can be roughly divided into two stages: from 1905 to 1999, the temperature fluctuated around the average value, 12.10 °C, and from 2000 to 2019, the temperature showed a significant upward trend. (2) The analysis of the climate-tree growth relationship indicated that the response of radial growth of Juniperus squamata to temperature was significantly stronger than the response to precipitation; especially in the last 20 years, when the radial growth of Juniperus squamata was positively correlated with temperature (p < 0.01). Compared to the maximum temperature and mean temperature, the correlation between radial growth of Juniperus squamata and minimum temperature was more significant. (3) Under the background of climate warming, the radial growth trend of Juniperus squamata in Nangqian county was consistent with temperature changes. Particularly in the past 20 years, the radial growth of Juniperus squamata showed a significantly increased trend and entered a rapid growth period.

1. Introduction

As the highest plateau in the world, the Tibetan Plateau (TP) is an important barrier for China’s ecological and national security, with an area of about 2.5 million km2 and an average altitude of more than 4000 m [1,2,3,4]. It has a profound influence on China’s climate background. The TP significantly affects the Asian monsoon system and northern hemisphere atmospheric circulation due to its unique topography and dynamic and thermal effects [5,6,7,8]. It is also the source of major rivers in China, such as the Yangtze River, Yellow River, and Lancang River, and is a key area for the global energy and water cycle [9]. Furthermore, the TP is also the initiator, amplifier, and regulator of global climate change [8,10]. Therefore, it is of great significance to study the climate change of TP in the historical period. However, due to the scarcity of meteorological stations and observation data, the study of TP climate change and its cause based on instrumental data is limited, and it is insufficient to understand the whole picture of TP climate change.
Tree rings have the advantages of accurate dating, good continuity, and high resolution when they are used to reconstruct past climatic factors in a certain area according to climate-tree growth relationship analysis [11,12,13,14,15,16]. Numerous studies have obtained complete historical temperature and precipitation data based on tree ring width (TRW) data, which lays a foundation for future predictions [17]. For example, Cufar [18] and Büntgen [19] used oak trees to reconstruct summer meteorological data; D’arrigo [20], Buckley [21], and Youngblut [22] reconstructed the temperature data of North America for more than 300 years with TRW data. Linderholm et al. [23] reconstructed temperature and precipitation data from 1786 to 2000 using Pinus sylvestris L. Buckley et al. [24] established a TRW chronology in northern Thailand and studied regional dry-wet changes. Wilson et al. and Anchukaitis et al. reconstructed the last millennium’s summer temperature data with tree rings in the Northern Hemisphere [25,26]. Long time-scale climate change research is crucial to understand the climate characteristics and change patterns of the TP [27]. The reconstruction of historical average temperatures could enrich the understanding of regional thermal variation from a long-term perspective and from current warming trends, in addition to forest growth responses to such warming trends [15,16,28].
According to the fifth report of IPCC, global average surface temperatures have risen by 0.85 °C since 1880 [29]. Studies show that the temperature rise in the TP is about 30 years ahead of global warming; the average temperature increase in each decade is more than 0.3 °C, which is approximately twice the global average temperature increase over the same period, especially in winter [30,31,32]. Vegetation is the most active part of terrestrial ecosystem; it connects the atmosphere, hydrosphere, pedosphere, and biosphere. Tree growth is sensitive to climate change. Therefore, clarifying tree-climate dynamic response is helpful to evaluations of the impact of climate change on the forest ecosystem [33,34,35,36], to predictions of the change in the forest ecosystem carbon sink, and to sustainable forest management [37]. It is generally thought that warming will enhance growth of trees in cold regions, unless there is temperature-induced drought stress [38,39,40,41]. Numerous studies show that tree growth above the tree line is mainly limited by temperature, and the annual minimum temperature is the most significant indicator [42]. In alpine mountains, above the forest line, shrubs such as Juniperus pingii W. C. Cheng, Hippophae rhamnoides L., and Rhododendron simsii Planch., are controlled by the low temperature in the growing seasons [43]. Zhang et al. [44] used tree rings from Qamdo, eastern Tibet andfound that the mean minimum temperature was the most significant factor, and used tree rings to reconstruct the mean minimum temperature for the past 413 years. Gou et al. [45] found that in the source of the Yellow River, standard TRW chronology is significantly correlated to the minimum temperature. Zhu et al. [46] analyzed climate-Juniperus przewalskii Komarov. growth in Wulan area of Qinghai Province and found that the growth was limited by temperature. In the past decade, many ecologists have used changes in shrub growth and range as sensitive ecological indicators to reflect global climate change [47]. Compared to trees, shrubs distributed in alpine regions have better adaptability to alpine climates [43,48], and in most high-altitude regions, shrub species are more abundant than trees. The shrub distribution area is the concentrated population distribution area of the TP. In this paper, the radial growth process of Juniperus squamata and its correlation with climate factors, such as temperature and precipitation, were studied. Furthermore, attention was paid to the long-term change of the annual mean temperature in this area, to provide a scientific basis for the assessment of the suitability of shrub growth and human settlement on the TP, considering climate change effects, under the background of a warmer and more humid climate.
Juniperus squamata Buch.-Ham. ex D.Don is a shrub in the genus Cypress in the family Cypressidae, and is the representative cold-temperate needle-leaved forest in the eastern TP [49]. It has the characteristics of a deep root system, with strong stress resistance and long life span. Juniperus squamata are widely distributed in the source region of the Lancang River, China, and they play an important role in soil and water conservation, biodiversity conservation, and climate control. The objectives of the present study were to: (1) develop a TRW chronology of Juniperus squamata in the Three-River headwaters of the TP, and analyze its growth status, (2) determine the key climatic factors that limit radial growth of Juniperus squamata in the Three-River headwaters of the TP, and (3) reconstruct the mean temperature variability of the TP from July to September in the past 115 years in the context of climate warming.

2. Study Area

The sampling site is located at Nangqian County, Yushu Tibetan Autonomous Prefecture, Qinghai Province (Figure 1), with an average elevation of 4309 m. Located in the southeastern TP, Nangqian county is a transition zone from canyon area to the plateau. The landform is mountainous, high in the northwest and low in the southeast. The five major rivers, Baqu River, Ziqu River, Jiqu River, Requ River, and Zaqu River in Nangqian belong to the Lancang River system.
The climate of the area is under the influence of the westerly and monsoon circulations of the Indian and Pacific Oceans, with pronounced dry and wet seasonality [14,50,51,52]. According to the climate data (time span 1957–2019) from meteorological stations (96°16′48″ E, 32°7′12″ N, 3643.7 m) close to tree-ring sampling sites, from 1957 to 2019, the extreme maximum and minimum temperatures were 28.1 °C and −24.1 °C, respectively (Figure 2). The mean annual temperature was 4.4 °C, the total multi-years’ annual precipitation was 535.8 mm, and the mean annual relative humidity was 53% (Figure 3). The major shrub species in this area is Juniperus squamata; other shrub species mainly include Dasiphora fruticose L., Hippophae rhamnoides Linn., etc. [53]. The relatively steep topography at high altitudes in southeastern TP is not favorable for soil development; therefore, a thin soil layer consists mainly of Alpine desert soil, Alpine meadow soil, Chernozem, Mountain meadow soil, Grey cinnamon soil, Alpine steppe soil, Chestnut soil [54].

3. Materials and Methods

3.1. Tree-Ring Data

In July 2019, we collected 20 Juniperus squamata discs in an area with little human interference, in Nangqian County (96°25′30″ E, 31°56′51″ N), on the southeastern TP. The samples were mainly selected from the main branches or side branches with healthy growth and strong stems, and the discs were intercepted along the base of the stem. The collected discs were preprocessed according to the international tree ring sample processing standard. In the laboratory, after the shrub discs were naturally dried, we polished them with P400, P800, and P1200 sandpapers, in turn, until the tracheary element, parenchyma, and ring boundaries of Juniperus squamata could be clearly seen under the microscope [51]. We took a radius in the east-west and north-south directions of the discs’ facet, and the width of the rings was measured using WinDENDRO system with a resolution a 0.01 mm. The COFECHA software was utilized to check the results of cross-dating [52]. Finally, 40 radii of 20 discs were used for further study. A negative exponential curve or a linear regression with a negative slope was applied to remove the age-related trend for each series, using the ARSTAN program [50,51,52]. Tukey’s biweight robust mean method was applied to average the dimensionless TRW indices, to achieve a standard chronology. The reliability and signal strength of standard chronology was assessed using a 50-year moving expressed population signal (EPS) with a 25-year overlap, and the mean series inter-correlations.

3.2. Climate Data

Monthly maximum, minimum, mean temperatures, and precipitation climate data were obtained from the Nangqian meteorological station (96°16′48″ E, 32°7′12″ N; 3643.7 m altitude), which was the station nearest to the tree ring sampling site (straight-line distance, 28.94 km). We obtained interval data from 1957 to 2019 (National Meteorological Science Data Center: http://data.cma.cn/, accessed on 12 December 2021).

3.3. Data Analysis

To determine the key climatic factors influencing radial growth of tree, we conducted growth-climate correlation analyses between TRW chronology and climate variables (maximum temperature, minimum temperature, mean temperature, and precipitation) using the correlation and response function analyses provided in the Dendroclim 2002 software. We used a bootstrap response function analysis with 1000 resampling, at the statistically significant 95% confidence level, for the 1957–2019 period [51]. Considering that tree growth is not only influenced by climatic conditions during the current year, but also from the previous year (legacy effect), our climate-tree growth relationship analyses included monthly windows from the previous year’s June to the current year’s October. A moving correlation analysis between TRW chronology and annual mean temperature was conducted at 34-year intervals to investigate the temporal variability of the correlation between radial growth of Juniperus squamata and each limiting factor. A moving correlation analysis was also conducted between TRW chronology and precipitation to investigate the dynamic influences of precipitation on alpine shrub growth, and to test whether water stress influenced shrub growth.

3.4. Climate Reconstruction and Testing Methods

The precondition of tree-ring climate reconstruction is to determine the suitable meteorological factors, to establish the functional equations of the meteorological factor series and the chronological series through the correlation screening between radial growth and climate response analysis, and then, to extrapolate the region’s historical climate. In tree-ring climatology, most of the models used in climate reconstruction are linear regression models, such as:
C l i m = a × c r n + b
In this model, Clim is the meteorological factor series, crn is the chronological series, and a and b are the constants. The constants are derived from the common time period between the least squares chronology and the climate factors. The corresponding chronology data representing the blank climate factor series with no data in the same time period are substituted into a function to obtain the reconstructed climate factor series. Adjusting the equation for the size of the equation R2, allows the addition of independent variables to the right of the equation, such as chronological data from the previous year. The advantage of a multivariate linear equation is the ability to identify variables in a stepwise regression manner.
Due to the short recording period (1960–2019) at the Nangqian meteorological station, we used the “leave-one-out” calibration to test the stability of the reconstructed model. The specific implementation is described below.
In an overlapping period, eliminate a certain year, establish a regression equation, and then substitute the chronological value of the eliminated year into this equation to obtain estimated values of the year. Repeat the above steps until the estimated sequence of climatic elements is obtained, and compare against the measured values to test the stability [55]. In this study, the specific methods used, and parameters obtained for testing were the sign test, average value of product, reduced error, and correlation coefficient. The test results show if the equation is stable at the 0.05 significance level. Finally, compare the measured value with the reconstructed value and, accordingly, determine if the result is reliable.

4. Results

4.1. Climate Change Trend

The multi-year mean temperature of Nangqian from 1957 to 2019 was 4.43 ± 1.61 °C, of which, the lowest temperature was 2.82 °C, recorded in 1963 (Figure 3a). The temperature showed an obvious increasing trend from 1957 to 2019 (p < 0.001). Minimum temperature, mean temperature, and maximum temperature increased at 0.63 °C/10a, 0.36 °C/10a, and 0.34 °C/10a (Figure 3d), respectively. Increases have been more obvious in the minimum temperature (Tmin) than in the other indicators, especially during 2000–2019 (Figure 3b).
The annual precipitation showed an increasing trend (0.94 mm/a, p > 0.10, Figure 3e). The multi-year mean precipitation from 1957 to 2019 was 514.2 ± 169.3 mm. Meanwhile, the relative humidity (−0.03 mm/a) showed a decreasing trend (Figure 3f).

4.2. Climate-Tree Growth Response Analysis

The correlation analysis between the TRW chronology and the monthly meteorological factors showed that the radial growth of Juniperus squamata was more sensitive to temperature than to precipitation, in the Three-River Headwaters Region of the TP.
According to the results, the correlation between TRW chronology and Tmin data was significant at the 99% confidence level throughout all months, and was stronger than the correlations between TRW chronology and Tmax or Tmean. The correlation between TRW chronology and Tmean was significant at the 95% confidence level in previous October and current March, significant at the 99% confidence level from July to September (previous and current year), from the previous November to current February, and in April and May of that current year. The correlation between TRW chronology and Tmax data was also significant at the 95% confidence level in previous November and current February, at the 99% confidence level in the previous December and January and July to September (previous and current year) (Figure 4a). In addition, the result of the climate response function analysis between chronology and temperature from July to September (previous and current year) was statistically significant.
In general, the correlation between TRW chronology and precipitation was weak. The TRW chronology was positively correlated with precipitation in spring months (March, April, May) and in October (previous and current year) at the 99% confidence level; the correlation between TRW chronology and precipitation was strongest in the previous year’s October (Figure 4b).

4.3. Annual Tmean Reconstruction

Based on the ring width of Juniperus squamata, the TRW chronology was synthesized. According to the relevant results in the previous section, the standard chronology was significantly correlated with the mean temperature from July to September of the same year in the Three-River Headwaters Region of TP. It provided a basis to reconstruct the mean temperature of the growing season (July–September) in this region. Therefore, the current year chronology and the previous year chronology were used as independent variables in this study, and the mean temperature from July to September was reconstructed as the dependent variable to establish the regression equation, and the stability of the equation was tested using “leave-one-out” calibration.
We reconstructed historical (1905–2019) July to September Tmean data based on the relationship between TRW chronology and Tmean (1960–2019) in the Three-River Headwaters Region. The transfer function is as follows:
T7~9 = 10.638 + 0.536CSTDx + 1.129 CSTDx−1.
In the equation, T7–9 is the mean temperature series from July to September in the Three-River Headwaters Region, CSTDx is the standard chronological series of the current year, and CSTDx−1 is the standard chronological series of the previous year
The climate reconstruction data’s statistical parameters are listed in Table 1, and the results indicated that the reconstructions were reliable. For the reconstructed sequence (1960–2019), the correlation coefficient of the equation is 0.654, the explained variance 42.8%, and the corrected explained variance 40.8%, which indicates a good fit of the tree ring chronology to the mean temperature from July to September, which can explain the main variance of the measured sequence. F is the linear regression’s parameter that is tested for statistical significance; the value of F in this study was 21.33, indicating that the data fit the regression model very well, and the confidence in the regression equation is very high. The results of the sign test and product average test for the reconstructed equation were statistically significant at the 1% significance level, indicating that the reconstructed sequence is consistent with the measured sequence, regardless of high or low frequency. The value of the reduced error was 0.246. It is generally believed that the detected data are reliable when the reduction error is positive. The product average T value is reasonable, indicating that the consistency between the two sequences is reliable. Therefore, the reconstruction dates back to 1905 and tracked well to both high- and low- frequency variations of the mean temperature from July to September (Figure 5a).
The comparison (1960–2019) between the observed temperature series (black line) and the reconstructed series (red line) showed that the TRW chronology better reflects the upward trend of temperature over the past 60 years (Figure 5a). However, the inter-annual variation between the reconstructed series and the observed series differed greatly in individual years (>1.5 °C in 1965, 1976, and 1994). The largest difference occurred in 1976 (the reconstructed value is 1.7 °C higher than the measured value). Further comparison of first-differences of observed and reconstructed temperatures revealed high-frequency consistencies, except for individual years (1965, 1975–1976, 1980, 1993–1994) (Figure 5b).
From the July–September Tmean reconstruction series (1905 to 2019), the temperature change during this period can be roughly divided into two stages: from 1905 to 1999, when the temperature was within one standard deviation, and fluctuated around the average value,12.10 °C, and from 2000 to 2019, when the temperature showed a rapid rising trend and the growth rate was unprecedented over the past 95 years (Figure 5a).
In the 21st century, the radial growth of Juniperus squamata has accelerated dramatically as temperatures have risen, and the levels of increase were unprecedented in the last 95 years (Figure 6). The correlation results showed that the TRW chronology in the study area was closely related to the Tmean variation, reaching a 95% confidence level in 2008. In addition, the moving correlation results between the chronology and Tmean data gradually increased during the investigated period (1960–2019) (Figure 7). The results showed that the radial growth of Juniperus squamata was correlated with mean temperature.

5. Discussion

5.1. Chronology Statistics Characteristics

We established a standard TRW chronology of Juniperus squamata, using the sample size with a subsample signal strength (SSS > 0.85), as the starting point of the chronology (Figure 6). The chronology can reflect the inter-annual fluctuation characteristics of tree radial growth. The standard deviation (SD) value was 0.31, which contains more climate information [56]. The first-order autocorrelation (AR1) value was 0.69, indicating that the growth was strongly affected by the previous year, which meant a relatively strong “legacy effect” [57,58]. The mean inter-series correlation (Rbt) value is 0.24, indicating that the change of TRW between different discs is consistent. The expressed population signal (EPS) value was 0.89, higher than the “0.85” threshold value, indicating that the samples could represent the local environmental characteristics [59,60]. The signal-to-noise ratio (SNR) value was 7.96, higher than the “4” threshold value, indicating that the samples contain abundant climate signals [61]. The statistical characteristics indicated that the TRW data contained reliable climate information (Table 2).

5.2. Response to Temperature

The correlation analysis between the TRW chronology and the monthly meteorological factors showed that the radial growth of Juniperus squamata was more sensitive to temperature than to precipitation in the Three-River Headwaters Region of the TP. Moreover, the radial growth of Juniperus squamata was positively correlated with the temperature in most months (June of last year–October of current year), especially the monthly mean minimum temperature (Figure 4a).
The chronology was significantly positively correlated with the mean temperature from July to September. Research shows in some alpine and humid regions, temperature rise is conducive to the acceleration of the radial growth and the extension of the growing season, which increases forest productivity [35]. This is because temperature rise during the growing season (July to September) can promote the synthesis and transport of auxin, and the division of cambium cells, as well as the accumulation of carbohydrates. Some studies have pointed out the mean temperature threshold for the start time of tree cambium cells activity is 6–8 °C, In cold and humid regions, the end time of radial growth of trees is mainly controlled by temperature [35,61,62,63]. Therefore, the critical period of growth and development of Juniperus squamata in Nangqian County is mainly concentrated in May to September. The radial growth of Juniperus squamata was positively correlated to the mean temperature in March. A study of Smith fir trees in Sejila Mountain in southeastern Tibet showed that the lowest temperature threshold limiting the onset of wood partialization is 0.7 ± 0.4 °C [63,64,65], which is far lower than temperatures previously reported. Therefore, under the ideal state, the beginning of the growth period of the Juniperus squamata is around March (1.07 °C). There is a significant positive correlation between the radial growth of Juniperus squamata and temperature in December. Juniperus squamata is deep-rooted plant, driven by photosynthesis; its main nutrients come from the soil. In Nangqian County, the ground temperature in winter is higher than that in spring [66]. In spring, the ground temperature in most parts of the TP is still at a subzero state. Permafrost requires so much heat to melt, that its temperature cannot rise quickly. In winter, the ground freezes and releases heat, which makes its temperature decrease slowly and maintains a relatively high temperature [51,52,67,68]. This is a unique phenomenon in seasonal frozen soil or permafrost area of the TP. The surface temperature is mainly affected by atmospheric temperature. The radial growth of Juniperus squamata has a significant positive correlation with the maximum temperature from November of the previous year to February of the current year. Generally, the xylem is affected by the accumulated temperature in winter and spring, and the warm winter prevents the leaf tissue from freezing and maintains the normal metabolism of trees, which is beneficial to the next year tree growth [9,46,69,70]. Therefore, the higher the maximum temperature is, the more favorable to the xylem growth of Juniperus squamata, and the easier it is to form wide rings.
Based on the correlation between TRW chronology and climate variables, the annual Tmin was the key climatic factor influencing radial growth of tree in conifer species at higher elevations on the southeastern TP. The conclusion is consistent with other studies [15,16,37,71,72,73]. Critical low temperature limits xylem growth on the forest line in the southeastern TP, especially Tmin in July [3,65,74,75]. Temperature is the main factor restricting Sabina przewalskii on the northeastern margin of the Tibetan Plateau. The temperature of July and August is significantly positively correlated with its annual ring index [76,77,78]. The thermal index required for alpine shrubs at an altitude of 4000–4500 m is the annual average temperature of −1~2 °C, the average temperature of the hottest month is 7~10 °C, the average temperature of the coldest month is −11~−8 °C, and the extreme minimum temperature is −30~−8 °C [9,35,56,65]. In Nangqian County, from 1957 to 2019, the mean temperature was 4.43 °C, the hottest mean monthly temperature was 25.16 °C, and the coldest mean monthly temperature was −19.80 °C. Among the three indicators, both the annual mean temperature and the hottest mean monthly temperature met the thermal index. Clearly, the minimum temperature was the main factor restricting growth, which is consistent with other results.

5.3. Response to Precipitation

The correlation analysis between the TRW chronology and the monthly meteorological factors showed that the radial growth of Juniperus squamata had a weak correlation with precipitation in the Three-River Headwaters Region of the TP (Figure 4b).
The radial growth of Juniperus squamata is positively correlated with the precipitation from March to May, indicating that its increases make the division of cambium cells more vigorous. Moreover, the precipitation and humid air environment in early spring are conducive to the budding of leaves [79,80,81]. The radial growth of Juniperus squamata is positively correlated with the precipitation in October because the cambium cells is still active at the end of the growing season, and rainwater can promote photosynthesis and accumulate nutrients. Moreover, precipitation makes the soil hydration sufficient, and the Juniperus squamata does not need to use the water in the body, so fallen leaves are reduced [82].

5.4. Fluctuation Characteristics of the Reconstruction Series

The reconstructed series reflects the temperature fluctuations from July to September in past 115 years in the Three-River headwaters region, which can be roughly divided into two stages (Figure 5a). From 1905 to 1999, the mean temperature was within one standard deviation, and fluctuated around the average value; this is similar to the trend of mean maximum temperature from May to June between 1360 and 2005 with Sabina Tibetica tree rings reconstructed by Shi et al. in Zaduo, Qinghai Province [83]. From 2000 to 2019, the temperature showed a rapid rising trend, with a growth rate that was unprecedented in the past 95 years. Kyimu et al. found that the minimum temperature has also been increasing in the past two decades when the minimum temperature was reconstructed with tree rings in the southeast of the Tibetan Plateau [15]. The explained variance of our July–September mean temperature reconstruction (42.8% of actual temperature) is comparable to, or even higher than, previous climate reconstructions of the TP and nearby regions [46,51,52]. In the 1980s, 1987 was a turning point, as the plateau monsoon had reached its peak in 1984. With stronger monsoons, temperature increases, and there is a three-year lag in its impact on summer and autumn temperatures. There are significant differences in the temperature distribution in the region around 1987 [51,79,84,85,86], and this point in time is consistent with our study. From 1987 to 1997 in the Three-River Headwaters Region, there were 7 years with long durations and wide ranges of snow cover, specifically in 1987, 1989, 1990, 1993, 1994, 1995 and 1997 [28,74,77,87]. Therefore, after 1987, there was a 7-year cold period. After 1998, the increasing rate of summer temperatures has accelerated significantly [88,89]. This year was in the late-1990s, which was an important turning point (Figure 5a). The average temperature from July to September could not represent the complete summer or autumn. The growth of grassland is significantly correlated with temperature in the same period, and the maximum forage yield in the region is in the months of July and August [75]. Since 2006, the forage yield of different grassland types has shown an increasing trend, and vegetation is getting better, which has a good correspondence with the sequence of this study. The temperature series shows that after 2006, the climate has entered to a warm period.
Further comparison of first-differences of observed and reconstructed temperatures revealed high-frequency consistencies, except for individual years (1965, 1975–1976, 1980, 1993–1994). Although there was some difference in the first-differences between the reconstruction value and the observed value in these years (>1.5 °C in 1976, 1980, and 1993), the largest difference occurred in 1975 (the reconstruction value was 1.8 °C higher than the observed value), but the difference was within the acceptable range (Figure 5b). This is similar to the first-differences analysis performed by Panthi et al. who described retrospective winter temperatures based on the ring width of Rhododendron campanulatum D.Don shrubs in the central Himalaya [90]. Considering that tree growth is affected by a variety of external factors, and temperature from July to September is but one of them, it cannot be ruled out that trees are impacted by other factors in individual years. Generally, however, the reconstructed series change trend fit with the observed temperature (from July to September). Furthermore, it was shown that the shrub Juniperus squamata was closely correlated with temperature and had good growth potential.

6. Conclusions

In order to investigate the effects of climate change on alpine shrubs, we have established a standard chronology of Juniperus squamata tree-ring width in the Three-River Headwaters Region of the Tibetan Plateau by using a dendroclimatology method, and its response to climate change was analyzed. Based on the relationship between temperature and tree growth, the mean temperature from July to September in the past 115 years was reconstructed. There were three primary findings.
(1) Meteorological data show that the temperature in Nangqian County of the Tibetan Plateau has increased continuously in the past 60 years, and the minimum temperature has the most significant change (0.63 °C/10a), especially between 2000–2019 (0.12 °C/a), while the precipitation has not changed significantly (0.94 mm/a, p > 0.10). The changes in our reconstructed temperature series can be roughly divided into two stages: from 1905 to 1999, when the temperature fluctuated around the average value, 12.10 °C, and from 2000 to 2019, when the temperature showed a significant upward trend (0.08 °C/a). The similarity and verification results between the current reconstructed series and the previously reported historical temperature records suggests reliability of the new reconstructed series.
(2) During 1957–2019, the response of radial growth of Juniperus squamata to temperature was significantly stronger than that of precipitation, especially in the last 20 years; the radial growth of Juniperus squamata was positively correlated with temperature (p < 0.01). Compared to the maximum temperature and mean temperature, the correlation between radial growth of Juniperus squamata and minimum temperature was more significant.
(3) According to our results, radial growth of high altitude shrubs in the study area kept pace with rising temperature, which could be divided into two stages: a steady increase trend from 1905–1999; a significant increase from 2000–2019, where temperature entered a rapid growth period.

Author Contributions

Conceptualization, G.Z. and Z.X.; methodology, Z.L.; software, M.K.; validation, Z.X., Z.L. and M.K.; formal analysis, Z.X.; investigation, J.L. and Y.H.; resources, J.L. and Y.H.; data curation, G.Z.; writing—original draft preparation, G.Z.; writing—review and editing, Z.X., Z.L. and M.K.; visualization, G.Z.; supervision, G.Z.; project administration, G.Z.; funding acquisition, Z.X. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Second Tibetan Plateau Scientific Expedition and Research Program (STEP), (2019QZKK0608), the National Natural Science Foundation of China (42177319), and the Beijing Municipal Education Commission for Inter-disciplinary Program “Ecological Restoration Engineering”.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Sun, H.; Zheng, D.; Yao, T.; Zhang, Y. Protection and Construction of the National Ecological Security Shelter Zone on Tibetan Plateau. Acta Geogr. Sin. 2012, 67, 3–12. [Google Scholar]
  2. Yan, S.; Wang, L.; Zou, L.; Pan, W. Relationship between soil organic carbon and total nitrogen and soil properties under different use patterns of grassland in the Qinghai-Tibet Plateau. Ecol. Sci. 2019, 38, 105–111. [Google Scholar]
  3. Gao, Z.; Xie, Y.; Wang, N.; Jiang, G.; Zhou, P. Response of Three Global DEM Data Accuracy to Different Terrain Factors in Qinghai-Tibet Plateau. Bull. Soil Water Conserv. 2019, 39, 190–197. [Google Scholar]
  4. Feng, Z.; Li, W.; Li, P.; Xiao, C. Relief degree of land surface and its geographical meanings in the Qinghai Tibetan Plateau, China. J. Geogr. 2020, 75, 1359–1372. [Google Scholar]
  5. Wang, L.; Zheng, Q.; Song, Q. Study on the influence of the lower mat of the Qinghai-Tibet Plateau on summer circulation in China. J. Nanjing Inst. Meteorol. 2002, 25, 186–191. [Google Scholar]
  6. Xu, Y.; Zhao, Z.; Li, D. Simulation analysis of climate change in the next 50 years on the Qinghai-Tibet Plateau and along the railway. Plateau Meteorol. 2005, 24, 700–707. [Google Scholar]
  7. Duan, A.; Wu, G.; Zhang, Q.; Liu, Y. New evidence that climate warming on the Qinghai-Tibet Plateau is the result of intensifying greenhouse gas emissions. Chin. Sci. Bull. 2006, 51, 989–992. [Google Scholar] [CrossRef]
  8. Ma, Y.; Hu, Z.; Tian, L.; Zhang, F.; Duan, A.; Yang, K.; Zhang, Y.; Yang, Y. Study Progresses of the Tibet Plateau Climate System Change and Mechanism of Its Impact on East Asia. Adv. Earth Sci. 2014, 29, 207–215. [Google Scholar]
  9. Yang, Y.; Hu, Z.; Lu, F.; Cai, Y.; Yu, H.; Guo, R.; Fu, C.; Fan, W.; Wu, D. Progress of Recent 60 Years’ Climate Change and Its Environmental Impacts on the Qinghai-Xizang Plateau. Plateau Meteorol. Syst. 2022, 41, 1–10. [Google Scholar]
  10. Pan, B.; Li, J. Qinghai-Tibetan Plateau: A Driver and Amplifier of the Global Climate Change—III. The effects of the uplift of Qinghai-Tibetan Plateau on Climate Changes. J. Lanzhou Univ. 1996, 32, 108–115. [Google Scholar]
  11. Kang, X.; Zhang, Q.; Graumlich, L.J.; Sheppard, P. Reconstruction of A 1835 a Past Climate for Dulan, Qinghai Province, Using Tree-Ring. J. Glacial Permafr. 2000, 22, 65–72. [Google Scholar]
  12. Smith, K.T. An organismal view of dendrochronology. Dendrochronologia 2008, 26, 185–193. [Google Scholar] [CrossRef]
  13. Balouet, J.C.; Smith, K.T.; Vroblesky, D.; Oudijk, G. Use of Dendrochronology and Dendrochemistry in Environmental Forensics: Does It Meet the Daubert Criteria? Environ. Forensics 2009, 10, 268–276. [Google Scholar] [CrossRef]
  14. Keyimu, M.; Li, Z.S.; Zhang, G.S.; Fan, Z.X.; Wang, X.C.; Fu, B.J. Tree ring–based minimum temperature reconstruction in the central Hengduan Mountains, China. Theor. Appl. Climatol. 2020, 141, 359–370. [Google Scholar] [CrossRef]
  15. Keyimu, M.; Li, Z.; Liu, G.; Fu, B.; Fan, Z.; Wang, X.; Wu, X.; Zhang, Y.; Halik, U. Tree-ring based minimum temperature reconstruction on the southeastern Tibetan Plateau. Quaternary Sci. Rev. 2021, 251, 106712. [Google Scholar] [CrossRef]
  16. Keyimu, M.; Li, Z.S.; Fu, B.J.; Liu, G.H.; Zeng, F.J.; Chen, W.L.; Fan, Z.X.; Fang, K.Y.; Wu, X.C.; Wang, X.C. A 406-year non-growing season precipitation reconstruction in the southeastern Tibetan Plateau. Clim. Past 2021, 17, 2381–2392. [Google Scholar] [CrossRef]
  17. Creasman, P.P.; Bannister, B.; Towner, R.H.; Dean, J.S.; Leavitt, S.W. Reflections On the Foundation, Persistence, and Growth of the Laboratory of Tree-Ring Research, Circa 1930–1960. Tree-Ring Res. 2012, 68, 81–89. [Google Scholar] [CrossRef] [Green Version]
  18. Cufar, K.; De Luis, M.; Berdajs, E. Main patterns of variability in beech tree-ring chronologies from different sites in Slovenia and their relation to climate. Zb. Gozdarstva Lesar. 2008, 87, 123–134. [Google Scholar]
  19. Büntgen, U.; Tegel, W.; Nicolussi, K.; McCormick, M.; Frank, D.; Trouet, V.; Kaplan, J.O.; Herzig, F.; Heussner, K.U.; Wanner, H.; et al. 2500 Years of European Climate Variability and Human Susceptibility. Science 2011, 331, 578–582. [Google Scholar] [CrossRef] [Green Version]
  20. D’Arrigo, R.D.; Jacoby, G.C. Northern North American tree-ring evidence for regional temperature changes after major volcanic events. Clim. Chang. 1999, 41, 1–15. [Google Scholar] [CrossRef]
  21. Buckley, B.M.; Wilson, R.J.; Kelly, P.E.; Larson, D.W.; Cook, E.R. Inferred summer precipitation for southern Ontario back to AD 610, as reconstructed from ring widths of Thuja occidentalis. Can. J. For. Res. 2004, 34, 2541–2553. [Google Scholar] [CrossRef]
  22. Youngblut, D.; Luckman, B. Maximum June–July temperatures in the southwest Yukon over the last 300 years reconstructed from tree rings. Dendrochronologia 2008, 25, 153–166. [Google Scholar] [CrossRef]
  23. Linderholm, H.; Niklasson, M.; Molin, T. Summer Moisture Variability in East Central Sweden Since the Mid-Eighteenth Century Recorded in Tree Rings. Geogr. Ann. Ser. A Phys. Geogr. 2004, 86, 277–287. [Google Scholar] [CrossRef]
  24. Buckley, B.M.; Palakit, K.; Duangsathaporn, K.; Sanguantham, P.; Prasomsin, P. Decadal scale droughts over northwestern Thailand over the past 448 years: Links to the tropical Pacific and Indian Ocean sectors. Clim. Dyn. 2007, 29, 63–71. [Google Scholar] [CrossRef]
  25. Wilson, R.; Anchukaitis, K.J.; Briffa, K.R.; Buntgen, U.; Cook, E.; D’Arrigo, R.; Davi, N.; Esper, J.; Frank, D.; Gunnarson, B. Last millennium northern hemisphere summer temperatures from tree rings: Part Ⅰ: The long term context. Quat. Sci. Rev. 2016, 134, 1–18. [Google Scholar] [CrossRef] [Green Version]
  26. Anchukaitis, K.J.; Wilson, R.; Briffa, K.R.; Buntgen, U.; Cook, E.; D’Arrigo, R.; Davi, N.; Esper, J.; Frank, D.; Gunnarson, B.; et al. Last millennium Northern Hemisphere summer temperatures from tree rings: Part Ⅱ: Spatially resolved reconstructions. Quat. Sci. Rev. 2017, 163, 1–22. [Google Scholar] [CrossRef] [Green Version]
  27. Marcott, S.; Shakun, J.; Clark, P.; Mix, A. A Reconstruction of Regional and Global Temperature for the Past 11,300 Years. Science 2013, 339, 1198–1201. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Yeshey, K.; Anan, P.; Supat, I.N.A. Dendroclimatic Reconstruction of Mean Annual Temperatures over Treeline Regions of Northern Bhutan Himalayas. Forests 2022, 13, 1794. [Google Scholar]
  29. Dahe, Q.; Stocker, T. Highlights of the IPCC Working Group 1 Fifth Assessment Report. Prog. Clim. Chang. Res. 2014, 10, 1–6. [Google Scholar]
  30. Pauli, H.; Gottfried, M.; Reiter, K.; Klettner, C.; Grabherr, G. Signals of range expansions and contractions of vascular plants in the high Alps: Observations (1994–2004) at the GLORIA master site Schrankogel, Tyrol, Austria. Glob. Chang. Biol. 2010, 13, 147–156. [Google Scholar] [CrossRef]
  31. Duan, A.; Xiao, Z. Does the climate warming hiatus exist over the Tibetan Plateau? Sci. Rep. 2015, 5, 13711. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Chen, D.; Xu, B.; Yao, T.; Guo, Z.; Cui, P.; Chen, F.; Zhang, R.; Zhang, X.; Zhang, Y.; Fan, J.; et al. Assessment of past, Present and Future environmental changes on the Tibetan Plateau. Chin. Sci. Bull. 2015, 60, 3025–3035 + 1–2. [Google Scholar]
  33. Mann, M.E.; Bradley, R.S.; Hughes, M.K. Northern hemisphere temperatures during the past millennium: Inferences, uncertainties, and limitations. Geophys. Res. Lett. 1999, 26, 759–762. [Google Scholar] [CrossRef] [Green Version]
  34. Fan, G.; Hua, W.; Huang, X.; Wang, Y.; Liu, Y.; Li, X.; Li, H. Advances of the Study on Influence of Vegetation Change over Tibetan Plateau on Regional Climate. Plateau Mt. Meteorol. Res. 2008, 28, 72–80. [Google Scholar]
  35. Jin, M.; Li, J.; Che, Z.; Wang, F.; Zhang, J.; Gou, X. Intra-annual radial growth responses of Qilian juniper (Juniperus przewalskii) to climate factors in the central Qilian Mountains. Acta Ecol. Sin. 2020, 40, 7699–7708. [Google Scholar]
  36. Yu, J.; Chen, J.; Zhou, G.; Liu, G.; Wang, Y.; Li, J.; Liu, Q. Response of Radial Growth of Abies forrestii and Picea Likiangensis to Climate Factors in the Central Hengduan Mountains, Southwest China. For. Sci. 2020, 56, 28–38. [Google Scholar]
  37. Guo, X.; Yu, B.; Zhang, S.; Li, J.; Wang, J.; Huang, J. Research Progress on Xylem Formation Dynamics and Its Regulation Mechanism. J. Trop. Subtrop. Bot. 2019, 27, 541–547. [Google Scholar]
  38. Barber, V.A.; Juday, G.P.; Finney, B.P. Reduced growth of Alaskan white spruce in the twentieth century from temperature-induced drought stress. Nature 2000, 405, 668–673. [Google Scholar] [CrossRef]
  39. Rehfeldt, G.E.; Tchebakova, N.M.; Parfenova, Y.I. Genetic responses to climate and climate-change in conifers of the temperate and boreal forests. Recent Res. Dev. Hum. Genet. 2004, 1, 113–130. [Google Scholar]
  40. Wilmking, M.; Juday, G.P.; Barber, V.A.; Zald, H.S.J. Recent climate warming forces contrasting growth responses of white spruce at treeline in Alaska through temperature thresholds. Glob. Chang. Biol. 2004, 10, 1724–1736. [Google Scholar] [CrossRef]
  41. Reich, P.B.; Oleksyn, J. Climate warming will reduce growth and survival of Scots pine except in the far north. Ecol. Lett. 2008, 11, 588–597. [Google Scholar] [CrossRef] [PubMed]
  42. Zhang, Z.; Kang, D.; Guo, W.; Zhao, L.; Cui, L.; Li, J. Climate sensitivity of purple cone spruce (Picea purpurea) across an altitudinal gradient on the eastern Tibetan Plateau. Dendrochronologia 2019, 56, 125586. [Google Scholar]
  43. Lu, X.; Liang, E. Progresses in dendrochronology of shrubs. Acta Ecol. Sin. 2013, 33, 1367–1374. [Google Scholar]
  44. Zhang, R.; Yuan, Y.; Wei, W.; Yu, S.; Shang, H.; Wang, L.; Chen, F. Analysis on Mean Minimum Temperature in the East Tibet from Autumn to Winter in Tree Ring of the past 400 Years. Plateau Meteorol. 2010, 29, 359–365. [Google Scholar]
  45. Gou, X.; Chen, F.; Yang, M.; Peng, J.; Zhang, Y.; Zhang, Y. The temperature rise in the source region of the Yellow River recorded by tree rings in the last hundred years. Collect. Abstr. 2006 Annu. Conf. Geogr. Soc. China 2006, 2, 60–61. [Google Scholar]
  46. Zhu, H.; Zheng, Y.; Shao, X.; Liu, X.; Xu, Y.; Liang, E. Millennial temperature reconstruction based on tree-ring widths of Qilian juniper from Wulan, Qinghai Province, China. Chin. Sci. Bull. 2008, 53, 3914–3920. [Google Scholar] [CrossRef] [Green Version]
  47. Xiao, S.; Peng, X.; Ding, A.; Tian, Q.; Han, C. Research Advance of Shrub Dendrochronology in the Cold and Arid Regions of China. Prog. Earth Sci. 2020, 35, 561–567. [Google Scholar]
  48. Myers-Smith, I.H.; Hallinger, M.; Blok, D.; Sass-Klaassen, U.; Rayback, S.A.; Weijers, S.; Trant, A.J.; Tape, K.D.; Naito, A.T.; Wipf, S. Methods for measuring arctic and alpine shrub growth: A review. Earth-Sci. Rev. 2015, 140, 1–13. [Google Scholar] [CrossRef]
  49. Song, Y.C.; Xu, G.S. A Scheme of Vegetation Classification of Taiwan, China. Acta Bot. Sin. 2003, 45, 883–895. [Google Scholar]
  50. Li, Z.; Zhang, Q.; Ma, K. Tree-ring reconstruction of summer temperature for A.D. 1475–2003 in the central Hengduan Mountains, Northwestern Yunnan, China. Clim. Chang. 2012, 110, 455–467. [Google Scholar] [CrossRef]
  51. Liang, H.; Lyu, L.; Wahab, M. A 382-year reconstruction of August mean minimum temperature from tree-ring maximum latewood density on the southeastern Tibetan Plateau, China. Dendrochronologia 2016, 37, 1–8. [Google Scholar] [CrossRef]
  52. Huang, R.; Zhu, H.; Liang, E.; Liu, B.; Shi, J.; Zhang, R.; Yuan, Y.; Griessinger, J. A tree ring-based winter temperature reconstruction for the southeastern Tibetan Plateau since 1340 CE. Clim. Dyn. 2019, 53, 3221–3233. [Google Scholar] [CrossRef]
  53. Cui, X.; Ma, H.; Huang, G.; Hou, M.; Xu, M.; Zheng, G.; Cui, B.; Zhuo, L.; Liao, C. Research on the Land Suitable for Planting 6 Major Tree Species in Qinghai Province. For. Resour. Manag. 2016, 61, 74–78. [Google Scholar]
  54. Zhou, Y.; Wang, B.; Liu, Y.; Wang, S.; Zhou, J. Geochemical evaluation of soil quality and land-use regionalization of selenium-rich soils in cultivated area around Nangqian County, Qinghai. Resour. Environ. Arid. Areas 2020, 34, 93–101. [Google Scholar]
  55. Liu, H.; Shao, X. Reconstruction of early-spring temperature at Zhen’an from 1755 using tree ring chronology. Acta Meteorol. Sin. 2000, 58, 223–233. [Google Scholar]
  56. Liu, M.; Ma, Z. Responses of biomass allocation to simulated warming in an alpine scrubland of eastern Qinghai-Tibetan Plateau. Acta Ecol. Sin. 2020, 41, 1421–1430. [Google Scholar]
  57. Zhang, T.; Yuan, Y.; He, Q.; Wei, W.; Diushen, M.; Shang, H.; Zhang, R. Development of tree-ring width chronologies and tree-growth response to climate in the mountains surrounding the Issyk-Kul Lake, Central Asia. Dendrochronologia 2014, 32, 230–236. [Google Scholar] [CrossRef]
  58. Zhao, W.; Luo, T.; Zhang, L. Relative impact of climate change and grazing on NDVI variations in typical alpine desert grassland in Tibet. Acta Ecol. Sin. 2019, 39, 8494–8503. [Google Scholar]
  59. Wigley, T.; Briffa, K.R.; Jones, P.D. On the average value of correlated time series with applications in dendroclimatology and hydrometeorology. J. Appl. Meteorol. 1984, 23, 201–213. [Google Scholar] [CrossRef]
  60. Marcon, A.K.; Longhi-Santos, T.; Galvão, F.; Martins, K.G.; Botosso, P.C.; Blum, C.T. Climatic Response of Cedrela fissilis Radial Growth in the Ombrophilous Mixed Forest, Paraná, Brazil. Floresta E Ambiente 2019, 26, e20180361. [Google Scholar] [CrossRef] [Green Version]
  61. Gao, L.; Gou, X.; Deng, Y.; Yang, M.; Zhang, F. Assessing the influences of tree species, elevation and climate on tree-ring growth in the Qilian Mountains of northwest China. Trees 2017, 31, 393–404. [Google Scholar] [CrossRef]
  62. Rossi, S.; Deslauriers, A.; Anfodillo, T.; Carraro, V. Evidence of threshold temperatures for xylogenesis in conifers at high altitudes. Oecologia 2007, 152, 1–12. [Google Scholar] [CrossRef] [PubMed]
  63. Li, X.; Liang, E.; Gričar, J.; Prislan, P.; Rossi, S.; Cufar, K. Age dependence of xylogenesis and its climatic sensitivity in Smith fir on the south-eastern Tibetan Plateau. Tree Physiol. 2013, 33, 48–56. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  64. Peng, J.; Gou, X.; Chen, F.; Li, J.; Liu, P.; Zhang, Y. Altitudinal variability of climate-tree growth relationships along a consistent slope of Anyemaqen Mountains, northeastern Tibetan Plateau. Dendrochronologia 2008, 26, 87–96. [Google Scholar] [CrossRef]
  65. Li, X.; Liang, E.; Gricar, J.; Rossi, S.; Cufar, K.; Ellison, A.M. Critical minimum temperature limits xylogenesis and maintains treelines on the southeastern Tibetan Plateau. Sci. Bull. 2017, 62, 804–812. [Google Scholar] [CrossRef] [Green Version]
  66. Kang, X.; Zhang, Q. The Growing Period of Cypress in High Elevation Area of the Tibetan Plateau. J. Glacial Permafr. 2001, 23, 149–155. [Google Scholar]
  67. Gou, X.; Peng, J.; Chen, F.; Yang, M.; Levia, D.F.; Li, J. A dendrochronological analysis of maximum summer half-year temperature variations over the past 700 years on the northeastern Tibetan Plateau. Theor. Appl. Climatol. 2008, 93, 195–206. [Google Scholar] [CrossRef]
  68. Linderholm, H.W. Annual temperature during the past 2485 years in the mid-eastern Tibetan Plateau inferred from tree rings. Sci. China Ser. D: Earth Sci. 2009, 52, 348–359. [Google Scholar]
  69. Bräuning, A. Summer temperature and summer monsoon history on the Tibetan Plateau during the last 400 years recorded by tree rings. Geophys. Res. Lett. 2004, 31, L24205. [Google Scholar] [CrossRef]
  70. Shao, X.; Huang, L.; Liu, H.; Liang, E.; Fang, X.; Wang, L. Reconstruction of precipitation variation from tree rings in recent 1000 years in Delingha, Qinghai. Sci. China (Ser. D) 2005, 48, 939–949. [Google Scholar] [CrossRef] [Green Version]
  71. Zheng, S. Study on Relationship between Tree-Ring and Climate Change. Sci. Silvae Sin. 2006, 42, 100–107. [Google Scholar]
  72. Wang, X.; Zhang, Q.; Ma, K.; Xiao, S. A tree-ring record of 500-year dry-wet changes in northern Tibet, China. Holocene 2008, 18, 579–588. [Google Scholar] [CrossRef] [Green Version]
  73. Tian, Q.; Gou, X.; Zhang, Q.; Wang, Y.; Fan, Z. May–June mean temperature reconstruction over the past 300 years based on tree rings in the Qilian Mountains of the northeastern Tibetan Plateau. IAWA J. 2009, 30, 421–434. [Google Scholar] [CrossRef]
  74. Wang, X.; Zhang, L.; Fang, J. Geographical Differences in Alpine Timberline and Its Climatic Interpretation in China. Acta Geogr. Sin. 2004, 59, 871–879. [Google Scholar]
  75. Han, B.; Zhou, B.; Zhao, H.; Shi, M.; Zhou, H.; Niu, D.; Fu, H. Analysis of Spatiotemporal Variation of Grassland Vegetation and its Impact Factors over Three-river Resource Region. Chin. J. Grassl. 2020, 42, 77–85. [Google Scholar]
  76. Zhang, R.; Yuan, Y.; Wei, W.; Gou, X.; Yu, S.; Shang, H.; Chen, F.; Zhang, T.; Qin, L. Dendroclimatic reconstruction of autumn–winter mean minimum temperature in the eastern Tibetan Plateau since 1600 AD. Dendrochronologia 2015, 33, 1–7. [Google Scholar] [CrossRef]
  77. Bi, Y.; Xu, J.; Yang, J.; Li, Z.; Gebrekirstos, A.; Liang, E.; Zhang, S.; Yang, Y.; Yang, Y.; Yang, X. Ring-widths of the above tree-line shrub Rhododendron reveal the change of minimum winter temperature over the past 211 years in Southwestern China. Clim. Dyn. 2017, 48, 3919–3933. [Google Scholar] [CrossRef] [Green Version]
  78. Cao, J.; Qi, R.; Liu, T.; Li, B.; Gao, B.; Chen, X.; Zhao, Y.; Zhao, Z. Patterns of species and phylogenetic diversity in Picea purpurea forests under different levels of disturbance on the northeastern Qinghai-Tibetan Plateau. Glob. Ecol. Conserv. 2021, 30, e01779. [Google Scholar] [CrossRef]
  79. Liu, Y.; Xin, Z.; Li, Z.; Keyimu, M. Climate effect on radial growth of Populus simonii in Northwest of Hebei for last four decades. Acta Ecol. Sin. 2020, 40, 9108–9119. [Google Scholar]
  80. Buchwal, A.; Rachlewicz, G.; Fonti, P.; Cherubini, P.; Gärtner, H. Temperature modulates intra-plant growth of Salix polaris from a High Arctic site (Svalbard). Polar Biol. 2013, 36, 1305–1318. [Google Scholar] [CrossRef] [Green Version]
  81. Stine, A.R.; Huybers, P. Arctic tree rings as recorders of variations in light availability. Nat. Commun. 2015, 5, 3836. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  82. Zhang, Y.; Yin, D.; Tian, K.; Zhang, W.; He, R.; He, W.; Sun, J.; Liu, Z. Radial growth responses of Picea Iikiangensis to climate variabilities at different altitudes in Yulong Snow Mountain, southwest China. Chin. J. Plant Ecol. 2018, 42, 629–639. [Google Scholar]
  83. Shi, X.H.; Qin, N.S.; Zhu, H.F.; Shao, X.M.; Wang, Q.C.; Zhu, X.D. May-June mean maximum temperature change during 1360-2005 as reconstructed by tree rings of Sabina Tibetica in Zaduo, Qinghai Province. Chin. Sci. Bull. 2010, 55, 3023–3029. [Google Scholar] [CrossRef]
  84. Hou, X.; Shi, J.; Li, L.; Lu, H. Growth response of Abies fargesii to climate in Shennongjia Mount of Hubei Province, Southeastern China. Chin. J. Appl. Ecol. 2015, 26, 689–696. [Google Scholar]
  85. Shao, X.; Fan, J. Past climate on the west Sichuan Plateau as reconstructed from ring-widths of dragon spruce. Quat. Sci. 1999, 19, 81–89. [Google Scholar]
  86. Li, Z.; Liu, G.; Zhang, Q.; Hu, C.; Luo, S.; Liu, X.; He, F. Tree ring reconstruction of summer temperature variations over the past 159 years in Wolong National Natural Reserve, western Sichuan, China. J. Plant Ecol. 2010, 34, 628–641. [Google Scholar]
  87. Gou, X.; Chen, F.; Jacoby, G.; Cook, E.; Yang, M.; Peng, H.; Zhang, Y. Rapid tree growth with respect to the last 400 years in response to climate warming, northeastern Tibetan Plateau. Int. J. Climatol. 2007, 27, 1497–1503. [Google Scholar] [CrossRef]
  88. Dai, S.; Li, L.; Liu, C.; Shi, X.; Yang, Y. Characteristics and Prediction Model of Summer Drought in Qinghai Province. J. Glacial Permafr. 2012, 34, 1433–1440. [Google Scholar]
  89. Yang, B. Spatial and temporal patterns of climate variations over the Tibetan Plateau during the period 1300–2010. Quat. Sci. 2012, 32, 81–94. [Google Scholar]
  90. Panthi, S.; Fan, Z.X.; Bräuning, A. Ring widths of Rhododendron shrubs reveal a persistent winter warming in the central Himalaya. Dendrochronologia 2021, 65, 125799. [Google Scholar] [CrossRef]
Figure 1. Location map of the tree-ring sample site and meteorological station in the Three-River Headwaters Region of the Tibetan Plateau.
Figure 1. Location map of the tree-ring sample site and meteorological station in the Three-River Headwaters Region of the Tibetan Plateau.
Forests 14 01219 g001
Figure 2. Seasonal distribution of temperature and precipitation in the Three-River Headwaters Region of the Tibetan Plateau during 1957 to 2019. Total monthly precipitation (blue bars), mean minimum temperature (yellow squares), mean maximum temperature (red circles) and mean temperature (green triangles).
Figure 2. Seasonal distribution of temperature and precipitation in the Three-River Headwaters Region of the Tibetan Plateau during 1957 to 2019. Total monthly precipitation (blue bars), mean minimum temperature (yellow squares), mean maximum temperature (red circles) and mean temperature (green triangles).
Forests 14 01219 g002
Figure 3. Climatic change trend in the Three-River Headwaters Region of the Tibetan Plateau from 1957 to 2019. (a) Mean temperature, Monthly mean temperature; (b) Min. temperature, Monthly mean minimum temperature; (c) Max. temperature, Monthly mean maximum temperature; (d) Temperature, Monthly mean maximum, mean, and minimum temperature; (e) Precipitation, Total annual precipitation; (f) Relative humidity, Monthly mean relative humidity.
Figure 3. Climatic change trend in the Three-River Headwaters Region of the Tibetan Plateau from 1957 to 2019. (a) Mean temperature, Monthly mean temperature; (b) Min. temperature, Monthly mean minimum temperature; (c) Max. temperature, Monthly mean maximum temperature; (d) Temperature, Monthly mean maximum, mean, and minimum temperature; (e) Precipitation, Total annual precipitation; (f) Relative humidity, Monthly mean relative humidity.
Forests 14 01219 g003aForests 14 01219 g003b
Figure 4. Results of correlation and response function analyses between TRW chronology of Juniperus squamata and (a) monthly temperatures (Tmin, minimum; Tmean, mean; and Tmax, maximum) and (b) precipitation (PRE) from prior June to current October for the period 1957–2019. The bar graphs denote the results of correlation analysis. The line graphs denote the results of response function analysis. The horizontal dashed and solid lines denote the 0.05 and 0.01 significance levels of the correlations, respectively. The “−” before the month denotes the previous year, For example, −6 denote June of the previous year.
Figure 4. Results of correlation and response function analyses between TRW chronology of Juniperus squamata and (a) monthly temperatures (Tmin, minimum; Tmean, mean; and Tmax, maximum) and (b) precipitation (PRE) from prior June to current October for the period 1957–2019. The bar graphs denote the results of correlation analysis. The line graphs denote the results of response function analysis. The horizontal dashed and solid lines denote the 0.05 and 0.01 significance levels of the correlations, respectively. The “−” before the month denotes the previous year, For example, −6 denote June of the previous year.
Forests 14 01219 g004
Figure 5. July–September mean temperature reconstruction in the Three-River Headwaters Region of the Tibetan Plateau. (a) Comparison between observed and reconstructed values of July–September mean temperature, (b) comparison between the first-difference data of observed and reconstructed July–September mean temperature.
Figure 5. July–September mean temperature reconstruction in the Three-River Headwaters Region of the Tibetan Plateau. (a) Comparison between observed and reconstructed values of July–September mean temperature, (b) comparison between the first-difference data of observed and reconstructed July–September mean temperature.
Forests 14 01219 g005
Figure 6. Standard TRW chronology and number of the discs sampled from Juniperus squamata in the Three-River Headwaters Region of the Tibetan Plateau.
Figure 6. Standard TRW chronology and number of the discs sampled from Juniperus squamata in the Three-River Headwaters Region of the Tibetan Plateau.
Forests 14 01219 g006
Figure 7. Results of moving correlation analysis between TRW chronology and annual mean temperature (Tmean) and precipitation (PRE) (using 34 years moving interval). The horizontal dashed and solid lines denote 95% and 99% confidence interval, respectively.
Figure 7. Results of moving correlation analysis between TRW chronology and annual mean temperature (Tmean) and precipitation (PRE) (using 34 years moving interval). The horizontal dashed and solid lines denote 95% and 99% confidence interval, respectively.
Forests 14 01219 g007
Table 1. Statistical parameters of mean temperature reconstruction process from July to September in the Three-River Headwaters Region of the Tibetan Plateau (1960–2019). r, correlation coefficient; R2, explained variance; R2adj, explained variance after adjustment; F, regression test statistic; ST, sign test; Er, reduced error; T, the product average; **, p < 0.01.
Table 1. Statistical parameters of mean temperature reconstruction process from July to September in the Three-River Headwaters Region of the Tibetan Plateau (1960–2019). r, correlation coefficient; R2, explained variance; R2adj, explained variance after adjustment; F, regression test statistic; ST, sign test; Er, reduced error; T, the product average; **, p < 0.01.
rR2R2adjFSTErT
0.6540.4280.40821.3330+/300.2466.88 **
Table 2. Statistics of the standard TRW chronology of Juniperus squamata in the Three-River Headwaters Region of the Tibetan Plateau. SD, standard deviation; AR1, first-order autocorrelation; Rbt, mean inter-series correlation; SNR, signal-to-noise ratio; EPS, express population signal.
Table 2. Statistics of the standard TRW chronology of Juniperus squamata in the Three-River Headwaters Region of the Tibetan Plateau. SD, standard deviation; AR1, first-order autocorrelation; Rbt, mean inter-series correlation; SNR, signal-to-noise ratio; EPS, express population signal.
LocationElevation(m)Time LengthSamples(n)Mean Value
Type96.25° E 31.56° N4309.01957–2019400.93
SDAR1RbtSNREPS
Juniperus squamata0.310.690.247.960.89
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Zhao, G.; Xin, Z.; Liu, J.; Huang, Y.; Keyimu, M.; Li, Z. The Radial Growth of Juniperus squamata Showed Sharp Increase in Response to Climate Warming on the Three-River Headwaters Region of Tibetan Plateau since the Early 21st Century. Forests 2023, 14, 1219. https://doi.org/10.3390/f14061219

AMA Style

Zhao G, Xin Z, Liu J, Huang Y, Keyimu M, Li Z. The Radial Growth of Juniperus squamata Showed Sharp Increase in Response to Climate Warming on the Three-River Headwaters Region of Tibetan Plateau since the Early 21st Century. Forests. 2023; 14(6):1219. https://doi.org/10.3390/f14061219

Chicago/Turabian Style

Zhao, Guoqing, Zhongbao Xin, Jinhao Liu, Yanzhang Huang, Maierdang Keyimu, and Zongshan Li. 2023. "The Radial Growth of Juniperus squamata Showed Sharp Increase in Response to Climate Warming on the Three-River Headwaters Region of Tibetan Plateau since the Early 21st Century" Forests 14, no. 6: 1219. https://doi.org/10.3390/f14061219

APA Style

Zhao, G., Xin, Z., Liu, J., Huang, Y., Keyimu, M., & Li, Z. (2023). The Radial Growth of Juniperus squamata Showed Sharp Increase in Response to Climate Warming on the Three-River Headwaters Region of Tibetan Plateau since the Early 21st Century. Forests, 14(6), 1219. https://doi.org/10.3390/f14061219

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop