A 406-year non-growing-season precipitation reconstruction in the southeastern Tibetan Plateau

Trees record climatic conditions during their growth, and tree rings serve as proxy to reveal the features of the historical climate of a region. In this study, we collected tree-ring cores of hemlock forest (Tsuga forrestii) from the northwestern Yunnan area of the southeastern Tibetan Plateau (SETP) and created a residual tree-ring width (TRW) chronology. An analysis of the relationship between tree growth and climate revealed that precipitation during the non-growing season (NGS) (from November of the previous year to February of the current year) was the most important constraining factor on the radial tree growth of hemlock forests in this region. In addition, the influence of NGS precipitation on radial tree growth was relatively uniform over time (1956–2005). Accordingly, we reconstructed the NGS precipitation over the period spanning from 1600– 2005. The reconstruction accounted for 28.5 % of the actual variance during the common period of 1956–2005. Based on the reconstruction, NGS was extremely dry during the years 1656, 1694, 1703, 1736, 1897, 1907, 1943, 1982 and 1999. In contrast, the NGS was extremely wet during the years 1627, 1638, 1654, 1832, 1834–1835 and 1992. Similar variations of the NGS precipitation reconstruction series and Palmer Drought Severity Index (PDSI) reconstructions of early growing season from surrounding regions indicated the reliability of the present reconstruction. A comparison of the reconstruction with Climate Research Unit (CRU) gridded data revealed that our reconstruction was representative of the NGS precipitation variability of a large region in the SETP. Our study provides the first historical NGS precipitation reconstruction in the SETP which enriches the understanding of the long-term climate variability of this region. The NGS precipitation showed slightly increasing trend during the last decade which might accelerate regional hemlock forest growth.

. These long-term records enable us to identify the interannual, decadal and multi-decal variability of historical climatic conditions. They also provide a reference to better understand the nature of current climatic conditions (warming/cooling, drying/wetting) and to project the future regional climate, as well as the dynamic response of Earth processes (e.g., forest growth, glacier retreat/advance, stream flow, drought frequency and forest fires) to climate change.
Being the "third pole" of the Earth, the Tibetan Plateau (TP) (average of 4000 m a.s.l.) is particularly sensitive to climate change and is one of the fastest warming places in the world (Chen et al., 2020). The average decadal temperature increase on the TP is 0.33 • C, which is higher than the world's average decadal temperature increase of 0.20 • C (Yan and Liu, 2014). Because of its geographical extent and position within the global circulation system, the TP plays a key role in regional and global atmospheric circulation patterns (Griessinger et al., 2017), not only affecting the midlatitude westerlies but also influencing the Asian monsoon circulation through its thermodynamical feedbacks (Duan et al., 2006;Rangwala, 2009;Wu et al., 2015).
There are large areas of coniferous forest distributed at high altitudes in the southeastern Tibetan Plateau (SETP). Due to their age and relative lack of disturbance, they are a source of proxy material (tree rings) that can be used to reveal the past climatic conditions in this region (Bräuning and Mantwill, 2004;Fan et al., 2009;Fang et al., 2010;Li et al., 2011;Wang et al., 2015;Shi et al., 2017;Huang et al., 2019;Shi et al., 2019;Keyimu et al., 2021). Many dendroclimatological reconstructions of hydroclimatic variables have also been conducted in the SETP (Fan et al., 2008;Zhang et al., 2015;Wernicke et al., 2015;Griessinger et al., 2017;He et al., 2018). However, few studies have focused on the reconstruction of precipitation history (He et al., 2012;Wernicke et al., 2015). The nongrowing season (NGS) of vegetation (from November of the previous year to February of the current year) includes the winter monsoon and pre-summer monsoon seasons in the SETP, and water availability during the NGS might therefore have a constraining effect on radial tree growth (Linderholm and Chen, 2005). It is important to understand the long-term precipitation variations during the NGS to evaluate the current trend of precipitation variation and estimate its future patterns and to determine the future responses of the forest ecosystem under the changing precipitation trend. To our knowledge, however, there have been no reports of the reconstruction of NGS precipitation in this area. This hinders our understanding of NGS variability from a long-term perspective.
In this study, we collected tree-ring cores of hemlock forest from the Xinzhu village of northwestern Yunnan in the SETP. The main objectives of the present study were to (1) develop a new tree-ring chronology and identify the re-sponses of hemlock forest radial growth to climate in the investigation area, (2) reconstruct the historical NGS precipitation and evaluate the recent NGS precipitation change in the long-term context and (3) validate the reliability of the reconstruction. Our results not only enrich the historical hydroclimatic information available in the SETP but also provide a basis to understand the current trend of regional NGS precipitation variation, which is relevant for evaluating the future development of regional forest ecosystem.

Study area and sampling sites
Tree-ring core samples were collected from Xinzhu Village in Lijiang County in northwestern Yunnan. The sample site was in the Hengduan Mountains in the SETP (Fig. 1). The climate of the study area is regulated by a westerly circulation and the monsoon circulations of the Indian and Pacific oceans. "Hengduan" means "transverse" in the Chinese language, which implies that the mountains in this region lie in the transverse direction from south to north, and the area is a passageway for the Indian monsoon to flow in and climb up to the TP and other parts of the mainland. The SETP is susceptible to monsoon flow and atmospheric circulations (Bräuning and Mantwill, 2004). According to the Weixi meteorological station of the China Meteorological Administration, which was the closest station to our sampling site, the mean annual precipitation was 953 mm from 1955 to 2016. Most of the annual precipitation (Nearly 70 %) concentrated in the monsoon season from May to October in this region (Fig. 2), and thus tree growth is usually constrained by water availability during the non-growing season. The coldest temperature was −2.9 • C in January and the warmest temperature was 18.6 • C in July. The topography of the sampling area is relatively steep, and it is not in favor of soil development; hence, a thin soil layer of alpine meadow soil (Chinese soil taxonomy) covers the bedrock. Hemlock is the dominant tree species of the sampling site, and its tree-ring cores were collected from trees which are healthy and relatively isolated, an optimal condition for maximizing climate signals in tree rings . The elevation of the sampling site was 2966 m a.s.l. A total of 48 tree-ring cores were extracted from 48 trees using a 5.1 mm diameter increment borer. We have used one sampling method per tree to improve the spatial representativity of radial tree growth. Sampling was conducted along an axis perpendicular to the slope inclination to avoid the impact of tension wood (Keyimu et al., 2020).

Establishment of the tree-ring chronology
The tree-ring samples were treated with standard dendrochronological procedures. They were first glued onto wooden holders and air dried, and then polished to a flat surface with sand paper until the tree rings were clearly visible.  Zhang et al., 2015). The blue dot is the meteorological station in Weixi County. On the right is the landscape image of tree-ring sampling site. The LINTAB 6.0 tree-ring measurement system (Rinntech Inc., Heidelberg, Germany) was used to measure the treering width (TRW). We have marked the tree rings of each sample at each 10-year interval and visually checked the treering pattern matching among samples, then confirmed the crossdating quality using the COFFECHA program (Holmes, 1983) (Arizona, USA). In total, 38 of the tree-ring cores were adopted for a further analysis after excluding the bad quality samples and the uncrossdated samples. The tree-ring series was detrended with a negative exponential model to remove the age dependency of tree growth (Cook et al., 1995). We have used the residual chronology since it removes the autocorrelation in tree-ring growth and captures a high-frequency climate signal. The "dplR" software toolkit (Bunn, 2018) within the R software environment (Bates et al., 2020) was used for detrending and chronology establishment. The reliable period of the chronology was determined based on the criterion of expressed population signal (EPS) > 0.85 (Wigley, 1984).

Climate data
Temperature and precipitation records were obtained from the Weixi meteorological station (27.17 • N, 99.28 • E; 2326 m a.s.l.) operated by the China Meteorological Administration. Data were available for the period of 1955-2005. Climate data (including the maximum, minimum and average temperatures, and precipitation) were provided by the China Meteorological Data Sharing Service Platform. A self-calibrated Palmer Drought Severity Index (scPDSI) was downloaded from the 3.26e gridded dataset of the Climate Research Unit (CRU) via the Royal Netherlands Meteorological Institute (KNMI) climate explorer (data accessed on 23 December 2020, data re-accessed for the updated version (CRU scPDSI 4.05 early) of PDSI data on 20 April 2021) using the coordinates of the tree-ring sampling site. The range of CRU grid box is 27.0-27.5 • N, 99.0-99.5 • E.

Tree growth and climate relationship analysis
We analyzed the relationship between climate and tree growth using DENDROCLIM2002 software (Biondi and Waikul, 2004). Pearson correlation values and response function values were calculated for the relationships between TRW indices and climate variables for the period of 1955-2005. Due to the carryover effect of the climatic conditions of the previous year on the current-year tree growth (Fritts, 1976), the tree growth-climate relationship analysis spanned a 16-month period from June of the previous year to September of the current year. We also used the seasonalized climate variables because it made more ecophysiological sense for growth than single months. To observe the temporal stability of the climate influence on radial tree growth, we conducted a moving correlation analysis at a moving interval of 32 years. All the correlation results were considered significant at the 95 % confidence level.

Statistics of chronology and climate reconstruction
We have used the expressed population signal (EPS) to determine the reliable period of the chronology; mean interseries correlation (Rbar), signal-to-noise ratio (SNR) and variance of first eigenvector (VFE) to evaluate the common signal among measurement series; standard deviation (SD) and mean sensitivity (MS) to show the degree of interannual variability of the chronology. According to the analysis of the relationship between the TRW indices and constraining climatic factors, we developed a linear regression model (Cook and Kairiukstis, 1990) for the climate reconstruction. As in many other tree-ring-based climate reconstructions, we tested the goodness of fit of the model using the leave-oneout cross-validation method (Michaelsen, 1987). We used the Pearson's correlation coefficient (r), explained variance (R 2 ), adjusted explained variance (R adj 2 ), reduction of error (RE), sign test (ST), coefficient of efficiency (CE) and product mean test (Pmt) to evaluate the fidelity of the reconstruction model (Fritts et al., 1990).

Characteristics of the TRW chronology
Residual TRW chronology of hemlock forests from the investigation area was established (Fig. 3). The descriptive statistics of the chronology were presented in Table 1. According to the criteria of EPS > 0.85, the most reliable length of the TRW chronology was 406 years (1600-2005). The mean correlation among tree-ring series (Rbar) was 0.48, and the variance in the first eigenvector (VFE) was 27 %, which implied a relatively strong common signal among individual trees constituting the chronology. The relatively low interannual variability of the chronology was expressed by the small mean sensitivity value (0.23). The EPS and SNR values (average EPS and SNR were 0.89 and 6.87 for the total length chronology, respectively) further implied the existence of the common signal among each individual measurement series. In general, all the statistical parameters indicated the potential climate signal imprinted in our TRW chronology.

Tree growth and climate relationship analysis
According to the results of the tree growth and climate relationship analyses (Fig. 4), the precipitation during the NGS was the most important constraining factor (R = 0.56, p < 0.001) on the radial growth of hemlock forests in the study area. The results of a response function analysis further confirmed the strong correlation between NGS precipitation and hemlock forest radial growth. The results of a moving correlation analysis between TRW chronology and instrumental NGS precipitation record (Fig. 5) were positively significant (at 99 %) during the investigated period , indicating that the NGS precipitation influence was stationary over time.

Non-growing-season precipitation reconstruction
According to the relationship between the TRW chronology and NGS precipitation, we developed a linear regression model (y = 229.94x−109.45 mm) and reconstructed the historical NGS precipitation series, which extended back to 1600 (Fig. 6b). In the model, y is the NGS precipitation, and x is the TRW index. The reconstruction accounted for 28.5 % of the instrumental NGS precipitation variability during the common time span . Figure 6a shows the similarities between the instrumental and reconstructed NGS precipitation series. We used a leave-one-out crossverification method to evaluate the legitimacy of the reconstruction model ( Table 2). The positive RE and CE values (0.18 and 0.15, respectively) were indicative of the legitimacy of the reconstruction. The significant value (at 95 %) of the sign test implied that the model-predicted values were generally in line with the variation trend of instrumental values. In addition, the significant values of the F test (at 99 %) and PM test (at 95 %) further confirmed the validity of the reconstruction. Overall, the statistics indicated that the reconstruction model possessed good predictive skills. Figure 6b shows the reconstructed NGS precipitation over the past 406 years . The mean of the reconstructed NGS precipitation series was 118 mm, and the standard deviation (SD) was 26 mm. We pre-defined the years that had NGS precipitation above 144 mm (mean + SD) as wet NGS years, and above 170 mm (mean + 2 SD) as extremely wet years (Table 3), whereas we defined years that had precipitation below 92 mm (mean − SD) as dry NGS years, and below 66 mm (mean − 2 SD) as extremely dry NGS years (Table 3). The dry/wet periods and some of the extreme dry/wet NGS periods in the present reconstruction were synchronized with dry/wet periods and extreme dry/wet periods in previously reported PDSI reconstruction from the surrounding region (Fig. 7, Tables S2 and S3), though some  dissimilarities also existed. As shown in Fig. 8, the instrumental (a, c) and reconstructed (b, d) NGS precipitation series could represent the climatic conditions over a similar area in the SETP.

Tree growth and climate relationship
The results of the tree growth and climate relationship analyses suggested that the hemlock forest radial growth in the northwestern Yunnan region of the SETP was strongly constrained by hydroclimatic factors. According to the Pearson correlation analysis, the influence of precipitation during the NGS on radial tree growth was greater than that of any other investigated climate variables and any correlation window. The response function analysis further confirmed the strong impact of NGS precipitation. In addition, the results of the 32-year interval of moving correlation analysis (Fig. 5) suggested the temporally consistent influence of NGS precipitation on hemlock forest radial growth in this region. The importance of NGS precipitation on the radial tree growth could be attributed to the fact that precipitation during the NGS compensated for the soil moisture, which was crucially important for supporting tree growth in the following season (Linderholm and Chen, 2005;Treydte et al., 2006;Wu et al., 2019;Li et al., 2021). This is because tree growth is often water stressed in the early stages of its growth in each year on the SETP when the monsoon precipitation does not arrive (Bräuning and Mantwill, 2004;Zhang et al., 2015), and the earlywood of tree rings mainly use spring meltwater (Zhu et al., 2021). The ecophysiological importance of NGS precipitation on tree growth and tree water usage was also revealed by isotope-ratio method-based investigations. Brinkmann et al.'s (2018) study showed that nearly 40 % of the water taken up by Fagus sylvatica and Picea abies trees in a temperate forest of middle Europe is sourced from NGS precipitation. Tree-ring oxygen isotope ratios (δ 18 O) are demonstrated to contain NGS precipitation signals in the Himalayan region Zhu et al., 2021). Huang et al.'s (2019) study revealed that NGS precipitation (snowfall) increased the snow depth and the later snowmelt compensated soil moisture in the spring and early summer, which was a crucially important water source for the Juniper growth in the southwestern Tibetan Plateau. Zhu et al.'s (2021) investigation in the western Himalaya revealed that formation of earlywood in tree rings of Pinus wallichiana depended on the snowmelt that originated from NGS precipitation. The weak influence of precipitation on regional hemlock forest growth during March and April and strong influence during May were connected with the saddle-shaped monthly rainfall pattern of this area (Fig. 2). The highest correlation between precipitation and TRW chronology was observed in May of the current year. This is because the xylogenous activity to form earlywood coincided with the low precipitation in this month (Fig. 2). In addition, the meltwater was probably used Note: R is the correlation coefficient, R 2 is the explained variance, R adj 2 is the adjusted explained variance, F is the F test, sign test is the sign of paired observed and estimated departures from their mean on the basis of the number of agreements/disagreements, Pmt is the product mean test, RE is the reduction of error, and CE is the coefficient of efficiency; * p < 0.05, * * p < 0.01. up (tree uptake combined with evaporation) during the early spring. Therefore, water stressed was increased during the late spring (May). The correlations between precipitation and the TRW chronology were not significant during the growing season (June-September) because an adequate water supply was available in the summer monsoon season. Precipitation during the NGS over the SETP falls as snow. According to Sommerfeld et al. (1993) and Stadler et al. (1996), the development of a snowpack insulates the underlying soil from freezing temperatures, which creates unfrozen soil conditions and most of the soil processes that are active during warmer conditions also persist under snow cover, albeit at a reduced rate (Edwards, 2007). Unfrozen soil can reduce the cold and frost damage to the shallow root systems of conifer trees in this region (Schenk and Jackson, 2002). A reduction in the cold damage to roots decreases the energy required to form new roots in the following growth year (Pederson et al., 2004), with the saved energy potentially used to initiate xylogenesis and form earlywood cells. Evergreen tree species are known to carry out year-round photosynthetic activity (Oquist and Huner, 2003;Prats and Brodersen, 2020), albeit at a slower rate during the NGS, and therefore the higher moisture availability contributes to the carbohydrate and energy accumulation process of hemlock forests in the investigation area. In contrast, the radial tree growth was negatively correlated to temperature in most correlation windows (Fig. 4). This can be explained by the fact that higher temperature enhances evapotranspiration, and thus decreases water availability, which eventually constrains tree growth. The negative impact of NGS temperature on radial tree growth was obvious because the strengthened evaporation due to higher temperatures might reduce the moisture compensation to the soil layer and cause water stress during the early stage of the following growth season.

Validity of the reconstructed precipitation series
We have tried to validate the fidelity of the newly reconstructed series from different aspects. Although we used the residual TRW chronology in the present study, which removes autocorrelation (Cook and Kairiukstis, 1990) to capture the high-frequency climate signals as in Fan et al. (2008) and Chen et al. (2016), the variability of dry and wet NGS at different scales was still retained in the reconstructed series. The reconstructed series in the present study demonstrated the variation in dry and wet NGS years (Fig. 6b). As in many other proxy-based historical climate reconstruction studies, we compared our NGS precipitation series with other hydroclimatic reconstructions from the surrounding areas to investigate the reliability of our reconstruction. However, there was no reported historical NGS precipitation record in the SETP, and we had to compare the present reconstruction series with available hydroclimatic reconstructions, e.g., PDSI. There are only countable numbers of PDSI recon-structions in the nearby region. Hence, we could only compare the present NGS precipitation reconstruction with existing PDSI reconstructions (Fig. 7) which are of spring or early summer. Dry/wet climate during these seasons is usually associated with the winter precipitation; hence, it makes certain sense to carry out the comparison. The correlation coefficients between our NGS precipitation reconstruction and the PDSI reconstructions of Fan et al. (2008), Fang et al. (2010), Zhang et al. (2015) and  were 0.51 (n = 702), 0.35 (n = 1062), 0.25 (n = 1062) and 0.22 (n = 1016) (p < 0.001). We have extracted the drought series of Asian Monsoon Atlas  from the nearest point to our investigation site and compared it with the NGS precipitation reconstruction in present study (R = 0.35, n = 1062, p < 0.001). As can be observed from Fig. 7, there were dry and wet periods in compared reconstruction series which were consistent with the NGS precipitation variabilities. These similarities indicated the reliability of our NGS precipitation reconstruction to some extent. The correlation  Zhang et al. (2015). (e) The current-year April-June average PDSI reconstruction in . (f) Drought series extracted from Asian Monsoon Atlas from the nearest point . The blue and purple bars show the common wet and dry periods of the different reconstructions, respectively. coefficients for the present reconstruction with those of Fan et al. (2008) and Fang et al. (2010) were greater than those with  and Zhang et al. (2015). These differences were probably due to the different distances among the study sites. Although the major dry and wet periods were similar in the hydroclimatic reconstructions referenced above, there were still certain discrepancies in duration and the strength of the dry/wet climatic conditions. This is probably because of the differences in the types of hydroclimatic variables (precipitation, PDSI), specific seasons reconstructed (annual, seasonal), tree species (species with different drought tolerances), chronology recording methods (standard chronology, residual chronology), length of calibration period, sample replication and the geomorphic differences of the tree-ring sampling sites (altitude, slope) (Table S1 in the Supplement).
In addition, we uploaded both of the instrumental and reconstructed NGS precipitation data for the same period of 1956-2005 on the KNMI website and conducted a spatial correlation analysis with the CRU gridded climate dataset. The similar patterns of spatial correlation between the instrumental and reconstructed data and their first-differenced data (Fig. 8) indicated that the present reconstruction was reliable and could represent the NGS precipitation over a large area in the SETP. Besides, the occurrence of some great historical drought events in the Asian monsoon area (Cook et al.,Figure 8. Spatial correlations of actual (a raw data; c first-differenced data) and reconstructed (b raw data; d first-differenced data) NGS precipitation with a gridded dataset of the NGS precipitation (average from November of the previous year to February of the current year) during their overlapping periods . The black square indicates the location of the study site. 2010; Kang et al., 2013), i.e., those of 1756-1768 (strange parallels drought), 1790, 1792-1796 (east India drought) and the 1920s (China mega-drought), matched the dry NGS periods in our reconstruction, which also further confirmed the reliability of our reconstruction.

Conclusions
In this study, we investigated 406 years of residual TRW chronology of hemlock forests in the SETP, China. The climate and tree growth relationship analyses showed that the TRW chronology was mostly negatively correlated with the thermal variable (temperature), whereas it was positively correlated with hydroclimatic variables (precipitation) and PDSI, indicating that hydroclimatic conditions determined the radial growth of hemlock forests in this region. Accordingly, we derived a linear model of the relationship between climate and tree growth, which accounted for 28.5 % of the actual NGS precipitation variance , and we used the model to reconstruct the historical (1600-2005) NGS precipitation. The reconstructed series showed that the NGS was extremely dry during the years 1656, 1694, 1703, 1736, 1897, 1907, 1943, 1982 and 1999. In contrast, the NGS was extremely wet during the years 1627, 1638, 1654, 1832, 1834-1835 and 1992. A comparison between the NGS precipitation reconstruction in this study and PDSI reconstructions from nearby regions revealed a coherency in the timing of dry and wet episodes, suggesting the reliability of our reconstruction. Our results showed that the NGS precipitation has demonstrated a slightly increasing trend since the 1980s, which is in favor of the future forest ecosystem development. In the future, more efforts should be made to collect the wide area of tree-ring data and develop more proxy chronologies that will enable us to reveal historical precipitation variability on a longer and wider scale in the SETP. Data availability. The climate reconstruction series in this study is provided in the Supplement.
Author contributions. ZL and MK conceived the study; ZL, ZF and XWa collected the tree-ring data; MK, ZL, ZF, KF and XWu elaborated on the methodology; MK, LZ and WC analyzed the data; MK and ZL led the writing of the manuscript; ZL and ZF revised the manuscript; BF and GL and FZ validated the final manuscript.
Competing interests. The authors declare that they have no conflict of interest.
Disclaimer. Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.