1. bookVolume 45 (2018): Issue 1 (January 2018)
Journal Details
License
Format
Journal
eISSN
1897-1695
First Published
04 Jul 2007
Publication timeframe
1 time per year
Languages
English
access type Open Access

Picea schrenkiana tree-ring chronologies development and vegetation index reconstruction for the Alatau Mountains, Central Asia

Published Online: 02 Jun 2018
Volume & Issue: Volume 45 (2018) - Issue 1 (January 2018)
Page range: 107 - 118
Received: 08 Jan 2018
Accepted: 16 Apr 2018
Journal Details
License
Format
Journal
eISSN
1897-1695
First Published
04 Jul 2007
Publication timeframe
1 time per year
Languages
English
Abstract

In this study, a total of 176 tree cores from Schrenk spruce (Picea schrenkiana) were used to establish a tree-ring chronology and a 167-year July–October normalized differential vegetation index (NDVI) for the Alatau Mountains in Central Asia was reconstructed using this newly developed chronology. The tree-ring based NDVI reconstruction tracks the observed data well (r=0.577, p<0.01, n=25) and precisely captures the drought events recorded in historical documents that occurred over a large area in 1917 and 1938. After applying a 21-year moving average, three dense (1860–1870, 1891–1907, and 1950–1974) and three sparse (1871–1890, 1908–1949, and 1975–2006) vegetation coverage periods were found in this reconstruction. Spatial correlation proves that the reconstructed NDVI series contains climatic signals representative for a large area including southern Kazakhstan. Although a comparison between this reconstruction and four climatic reconstructions for southeastern Kazakhstan, Nilka (in the Ili region), the Issyk Lake, and the Aksu region reveals similar variations, the coherence between these reconstructions become weak with the increase in spatial distance from north to south. In addition to the local representation, it was also demonstrated that the newly developed NDVI index can indicate the large-scale circulations over Eurasia, with the higher NDVI associated with stronger westerly winds from the Atlantic to the Alatau Mountains, and the lower NDVI associated with the weaker winds.

Keywords

Introduction

Terrestrial vegetation is regarded as a sensitive indicator of climate change (Pauli et al., 2002; Nagy, 2006; Jiapaer et al., 2015). The insights from vegetation phenology may help to increase the understanding of the interaction between climate and vegetation and the exchanges of energy and matter in a given area (Hmimina et al., 2013). The normalized difference vegetation index (NDVI) was developed as a measure of vegetation greenness and is typically used to reflect ecological environment changes and evaluate land surface phenology at local, regional, and global scales (Gamon et al., 2015). The near-infrared and red-light reflectance ratio is based on satellite sensor data reflected by the vegetation status. In this ratio, NDVI=(NIR–RED)/(NIR+RED), where RED represents the amount of red light, and NIR represents the amounts of near-infrared light (Pettorelli et al., 2005). However, these data, provided since the 1980s, do not adequately describe the dynamic variations of vegetation coverage over a long time scale. To address this information gap, proxy data, such as tree rings, are frequently selected. These data can be used to evaluate the relationship between the radial growth of coniferous species and vegetation coverage and reflect the seasonal variability of vegetation coverage for hundreds of years (Malmström et al., 1997; D’Arrigo et al., 2000; Leavitt et al., 2008). The above studies mostly focused on high-latitude (Kaufmann et al., 2004; Lopatin et al., 2006; Beck et al., 2013), high-altitude (He and Shao, 2006; Shang et al., 2016), and arid areas (Liang et al., 2009; Wang et al., 2014), and relied on samples of tree rings from coniferous species.

The widespread coniferous forests in Central Asia, comprising mostly local tree species (Juniperus turkestanica Kom., Larix sibirica Ledeb., Pinus wallichiana and Picea schrenkiana) provide a good opportunity for dendrochronological studies (Esper, 2000; Chen et al., 2012; Zhang et al., 2014; Seim et al., 2016; Opała et al., 2017). Wooden cores derived from Central Asia have been systematically collected since the 1990s (Bräuning, 1994; Esper et al., 2007). Many dendrochronological studies related to ring-width variations (Esper, 2000), tree growth climate response (Esper et al., 2003; Winter et al., 2009) and hydro-climate reconstruction (Zhang et al., 2015; Chen et al., 2017a) have been carried out. However, when compared with the previous dendrochronological studies addressing hydro-climatological reconstructions, studies of tree-ring based vegetation coverage variations for Central Asia are relatively rare. In addition, the knowledge of the corresponding relationships between the radial growth of conifers and the vegetation coverage in Central Asia is still limited.

Therefore, the aims of this study are to: (1) develop a tree-ring width-based regional chronology for the Alatau Mountains; (2) reconstruct historical NDVI series for the study area and explore its variations, and (3) assess climatic signals inhering in the newly developed NDVI reconstruction.

Materials And Methods
Sample collection

The Alatau Mountains form a boundary between China and Kazakhstan and are part of the Tian Shan Mountains of Central Asia. The Alatau Mountains extend more than 450 km from east to west, and the maximum elevation is approximately 4464 m a.s.l. The Alatau Mountains experience a temperate continental climate. The northern slopes of the Alatau Mountains have more annual precipitation (900–1000 mm) than the southern slopes and a lower snowline (3600 m a.s.l.). Approximately two-thirds of the mountains’ glaciers are distributed in the northern areas (He et al., 2014). The widespread, long-living, coniferous forest growing on the Alatau Mountains provides a good opportunity for carrying out dendrochronological studies to evaluate vegetation characteristics over a long-term period.

A long-lived Schrenk spruce (Picea schrenkiana) was selected as our study target. These spruces often grow up to 40 m in height and frequently live more than 200 years. More than 90% of the forest from 1200 to 2600 m a.s.l. in the Tien Shan Mountains contains Schrenk spruce, and it is the dominant tree species in the forests of the mountain zone. Tree-ring cores from three sites were sampled in July 2016. The study area and location of the sampling sites are shown in Fig. 1. Forest stands are moderately open and canopy densities are low. We selected healthy spruces with little evidence of damage by bushfire, landslides, earthquakes, human destruction, or animal invasion to avoid sampling non-climatic effects on radial growth. For most sampled trees, we parallelly extracted two cores from one tree at chest height. The two cores were collected from different directions at two flanks of a given tree. To collect tree cores that contained consistent climate signals, the altitude difference between the highest and lowest locations of one sampling site was less than 100 m. In total, we obtained 182 cores from 96 trees using increment borers at sites Amanboktek (code: AMA; location: 45.26°N–80.08°E; elevation: ~2050 m a.s.l.), Kikbay (KBZ; 45.21°N–79.97°E; ~1750 m a.s.l.), and Basika (BSK; 45.25°N–80.15°E; ~1450 m a.s.l.).

Fig. 1

(a) Location of the study area; (b) map of the study area; and (c) sampling sites: Amanboktek (AMA), Kikbay (KBZ), and Basika (BSK). The selected NDVI grid (44.932°–46.023°N, 78.986°–80.950°E) is represented by the dash-block diagram.

Ring-width measurements and chronology development

Following normal dendrochronological techniques (Speer, 2010), the sampled tree-ring cores were dried naturally and mounted on a wooden plank with grooves. Then, each core was sanded with abrasive paper and marked with needles under a microscope. Every ring on the sanded cores was measured under a binocular microscope using a Velmex Measuring System at a resolution 0.001 mm. This measuring system, the standard of North America’s Dendrological Research Community, is specifically designed for the researcher to perform non-contact measurement analysis. The system is composed of the following three main components: a UniSlide® dovetail rapid advance screw motion assembly, a linear encoder with one micron (0.001 mm) resolution, and a Velmex VRO™ Encoder Readout. Two programs, COFECHA (Grissino-Mayer, 2001) and ARSTAN (Cook and Krusic, 2005), were utilized for cross-dating quality control and to develop chronologies. To eliminate age trends that were affected by factors other than climate, the tree-ring width series were detrended using a negative exponential method, i.e. a conservative function. We used an expressed population signal (EPS) to evaluate the reliability of the chronology of the tree rings. The statistical analyses were done at 20-year intervals with an overlap of 10 years across the chronology. An EPS of ≥ 0.85 was used to ensure a reliable chronology length (Wigley et al., 1984).

Climatic data and vegetation index

For further analyses, we selected the precipitation, mean temperature, mean maximum temperature and mean minimum temperature per month of the Climatic Research Unit Time-Series (CRU TS 4.00; New et al., 2000), and the Palmer Drought Severity Index per month (scPDSI, CRU TS 3.25; Schrier et al., 2013) of the gridded 0.5° × 0.5° dataset (45°–46°N, 79°–81°E), available from 1901 to 2015, to describe the climatic conditions in the study area. Monthly normalized differential vegetation index (NDVI) data (44.932°–46.023°N, 78.986°–80.950°E) available from 1982 to 2006 were used to describe the vegetation conditions. The above climatic data and vegetation index were acquired from the Royal Netherlands Meteorological Institute (KNMI) Climate Explorer (http://climexp.knmi.nl).

Statistical analysis

We used 13-year reciprocal filters to decompose the developed chronologies of the tree rings into high- and low-pass components to evaluate the patterns of variation in various frequency ranges (Yuan et al., 2013). Variation patterns in these original and decomposed chronologies were assessed using Pearson correlations in the original, high-frequency, and low-frequency domains. Correlation analyses were also used to quantify the strengths of climatic signals and the vegetation conditions inherent in chronologies from spruces at the study site. We identified the strongest seasonal relationship between tree-ring width and vegetation data, and then a linear regression model was employed for reconstruction. We used Bootstrap (Young, 1994), Leave-one-out cross-validation (Michaelsen, 1987), and split-sample calibration-verification test (Meko and Graybill, 1995) methods to evaluate the statistical reliability of the reconstruction model. During the split-sample calibration-verification tests, the period of climatic data was split into two parts for calibration and verification. Several statistical parameters, including a reduction in error, coefficient of efficiency, and product mean test were calculated to evaluate the observed and estimated data (Cook et al., 1999). Spatial correlation was used to identify coherence between our reconstruction and the gridded 0.5° × 0.5° CRU self-calibrating PDSI 3.21 dataset (Wells et al., 2004) over a large region. In addition, statistical data from the Royal Netherlands Meteorological Institute (KNMI) were also used. Some tree-ring based climatic reconstructions and the newly developed vegetation coverage series were standardized using zero-mean normalization and decomposed using 13-year reciprocal filters to assess their coherence in a larger spatial context. Reanalysis 1 from the National Centers for Environmental Prediction/Department of Energy (NCEP/DOE) (Kalnay et al., 1996) was utilized to derive the large-scale circulation associated with the variation in the NDVI index. The Reanalysis is a meteorological data product that aims to represent the historical state of the atmosphere. Reanalysis 1 from the NCEP/DOE is one example of these Reanalysis data, and has a horizontal resolution of 2.5° × 2.5°. The climatology data from 1951 to 2000 has been removed to obtain the anomalies.

Results
Meteorological data and NDVI analysis

Fig. 2a shows that the highest temperature periods in the study area were in summertime (from June to August), with peaks in July. Most of the annual precipitation fell in May, June, and July (120.1 mm). The two peak values of precipitation in May (42.2 mm) and November (30.1 mm), which account for 12.6% and 8.9%, respectively, of the total annual precipitation. Wintertime has less precipitation than summer. August is typically the wettest month of the year (PDSI=–0.46) for the study area (Fig. 2b). The climate data for 1901–present reveal remarkable rising tendencies in the study area for the following: annual total precipitation (Y=0.5466X–734.15, R2=0.066, p<0.01) (Fig. 2d), annual mean temperature (Y=0.0168X–29.944, R2=0.357, p<0.001) (Fig. 2e), annual mean minimum temperature (Y=0.0222X–47.377, R2=0.500, p<0.001) (Fig. 2f), annual mean maximum temperature (Y=0.0113X–12.533, R2=0.163, p<0.001) (Fig. 2g), and annual mean PDSI (Y=0.0136X–27.148, R2=0.069, p<0.01) (Fig. 2h).

Fig. 2

Gridded climatic and NDVI data for the Alatau Mountain: (a) Precipitation, mean temperature, mean maximum temperature and mean minimum temperature per month. Bars indicate precipitation (in mm), and curves with different colors represent mean temperature, mean minimum temperature and mean maximum temperature (in °C); (b) Monthly PDSI; (c) Monthly NDVI; (d) Annual precipitation during 1901–2015; (e–g) Annual mean temperature, annual mean minimum temperature and annual mean maximum temperature during 1901–2015; (h) Annual mean PDSI during 1901–2015; and (i) Annual mean NDVI during 1982–2006. Dashed lines represent trends in climate and NDVI data.

Fig. 2c indicates that the monthly NDVI values from January to March are relatively low (not exceeding 0.140) and without obvious fluctuations, but have an evident increase starting at April (0.248), and peaking in July (0.536), when the study area reaches the heathiest vegetation period of the year. Thereafter, the monthly NDVI values maintain a long-term decrease until December. An increasing trend of annual mean NDVI since 1982 is less significant (Y=0.0002X–0.0673, R2=0.013, p<0.05) (Fig. 2i). The variations in NDVI in the study area are in accordance with findings for Central Asia reported by Yin et al. (2017). The correlation coefficients between NDVI and annual precipitation and annual mean PDSI are 0.352 (p<0.1, n=25) and 0.440 (p<0.05, n=25) respectively, while the correlation coefficients between NDVI and annual mean temperature, annual mean minimum temperature and annual mean maximum temperature are not significant.

Statistical characteristics of chronologies and development of a regional chronology

After cross-dating quality control, two cores from the AMA site and four cores from the KBZ site were rejected because of the low correlations between the subseries and the master series. Ultimately, 50 cores from 28 spruces at the AMA site, 81 cores from 42 spruces at the KBZ site, and 45 cores from 23 spruces at the BSK site were used to develop ring-width chronologies. The oldest tree at the KBZ site was nearly 195 years old. In order to carry out a further analysis, we utilized the standard (STD) chronologies that contained the common variations of respective series of tree samples and retained a low-through high-frequency common variance. This variance is hypothesized as being dependent on climate (Cook, 1985). Three chronologies, the depths of their samples, EPS, and Rbar are demonstrated in Fig. 3. General statistics of these chronologies for a common period of analysis (from 1900 to 2010) are listed in Table 1. The reliable lengths of AMA, KBZ, and BSK chronologies were 107 (1910–2016), 167 (1850–2016), and 137 (1880–2016) years, respectively.

Fig. 3

Chronologies (thin lines) with their sample depth (dashed lines). Red and blue lines represent EPS and Rbar data: (a) AMA; (b) KBZ; (c) BSK; and (d) ARC.

Statistical characteristics of chronologies over the common period of 1900–2010

StatisticAMAKBZBSKARC
Standard deviation (SD)0.1980.2240.4140.208
Mean sensitivity (MS)0.1800.1860.3790.180
First-order autocorrelation (AC1)0.3960.4960.5770.458
Interseries correlation (trees)0.4170.4840.7380.418
Interseries correlation (all series)0.4310.4880.7530.422
Mean within-tree correlation0.7860.8050.9610.808
Signal-to-noise ratio (SNR)9.85862.01218.29362.025
Expressed population signal (EPS)0.9080.9840.9480.984
The first principal component (PC#1)0.4920.5180.7990.461
First year EPS >0.851910185018801850

In addition to the reciprocal filters, we also used Pearson correlation coefficients to analyze the three sets of data, which included original unfiltered data, high-pass filtered data, and low-pass filtered data. As shown in Table 2, the correlations for the three chronologies employing the low-pass filter are relatively low, but the confidence levels for the correlation coefficients were all over 99.9% in the case of the original, high-frequency, and low-frequency domains over the common period of 1910–2016. In the common period, the years of 1924 (value: 1.452), 1970 (1.421), 1993 (1.319), 1973 (1.307), 1994 (1.292), 1950 (1.282), 1951 (1.275), 1923 (1.27), 1948 (1.264), and 1990 (1.257) are regarded as 10 highest-value years of AMA chronology, and the 10 lowest-value years appear in 2008 (0.511), 1945 (0.607), 1913 (0.614), 1917 (0.625), 1914 (0.632), 1911 (0.645), 2000 (0.675), 1927 (0.694), 1995 (0.705), and 2014 (0.708). The 10 highest-value years of KBZ chronology are 1993 (1.499), 1994 (1.448), 1968 (1.419), 1962 (1.395), 1973 (1.392), 1969 (1.369), 1924 (1.349), 2005 (1.341), 1989 (1.336), and 1970 (1.304), and the 10 lowest-value years appear in 1945 (0.457), 1927 (0.496), 1917 (0.509), 1944 (0.554), 1943 (0.596), 2000 (0.602), 1928 (0.626), 1938 (0.648), 1946 (0.66), and 2015 (0.686). The 10 highest-value years of BSK chronology are 1929 (1.678), 1935 (1.573), 1970 (1.556), 1952 (1.517), 1964 (1.502), 1960 (1.473), 1962 (1.473), 1993 (1.467), 2010 (1.456), and 1954 (1.44), and the 10 lowest-value years appear in 2008 (0.001), 1985 (0.137), 2000 (0.211), 1997 (0.226), 1927 (0.249), 1991 (0.249), 1917 (0.286), 1977 (0.304), 1974 (0.322), and 2014 (0.334). After comparing these extreme-value years, two highest-value years (1970 and 1993) and three lowest-value years (1917, 1927 and 2000) were observed in the three chronologies. The results indicated a good coherence of extreme-values among these chronologies. As a result, we combined all the ring-width data from the AMA, KBZ, and BSK sites in order to establish a regional chronology (ARC) that was longer and more replicated.

Correlation coefficients for three chronologies over the common period of 1910–2016.

Original (n=107)High-frequency (n=95)Low frequency (n=95)
AMAKBZBSKAMAKBZBSKAMAKBZBSK
AMA10.783

Significant at p<0.001.

0.660

Significant at p<0.001.

10.842

Significant at p<0.001.

0.682

Significant at p<0.001.

10.641

Significant at p<0.001.

0.683

Significant at p<0.001.

KBZ/10.624

Significant at p<0.001.

/10.759

Significant at p<0.001.

/10.388

Significant at p<0.001.

BSK//1//1//1

Note: The table shows correlation coefficients. Results for the original, high-, and low-pass filtered chronologies are shown.

The ARC chronology and the depths of the samples are shown in Fig. 3d, and the general descriptive statistics are listed in Table 1. Values of standard deviation (SD) and mean sensitivity (MS) for the chronologies obtained using the ARC chronology are slightly lower than those of the chronologies for single sampling sites because of the combination of samples. The first-order autocorrelation (AC1) assesses relationships between tree growth in a current year and previous growth. These values ranged from 0.396 to 0.577, and revealed that chronologies possessed low-frequency variance, which was affected by the lag effects of climate and tree physiology. The interseries correlations of the three chronologies obtained from sites of single sampling and the regional chronology all exceed 0.4. These values indicated a good coherence in the tree-growth. The ARC chronology with higher signal-to-noise ratios (SNR) and EPS demonstrated more climatic signals in the regional chronology. The reliable length of the ARC chronology was 167 (1850–2016) years based on the initial year of EPS>0.85.

Correlation analysis and physiological mechanisms

A strong biological lag effect was indicated by the high values of AC1 from 0.396 to 0.577 for the chronologies from single and composite sites (Table 1). Therefore, the climatic data per month (1901–2015) from July in the previous year to October in the current year were applied to assess the influence of climatic factors on radial growth of spruces in our study area. The results of the correlation analysis revealed that the relationship between the radial growth of spruces and precipitation was positive in general, and that significant correlations appeared in July of the previous year and from March to July of the current year (Fig. 4). There was a negative correlation between the tree-ring width of spruces and the temperature. Generally, these chronologies had very negative correlations with the temperature at the end of the growing season of the previous year and in the middle of the growing season of the current year. Furthermore, almost all of monthly the positive monthly correlation coefficients between chronologies and PDSI from July of the previous year to October of the current exceeded the 0.05 significance level. The above results demonstrated that moisture was the main climatic limitation on development of spruces’ tree rings in the study area. More precipitation at the end of the growing season of the previous year and the fast growing season of the current year may enhance the potential for accumulating water reserves in the soil, thus resulting in larger leaves (Liu et al., 2011), larger roots (Lebourgeois et al., 2004), and more nutrient in buds (Barbaroux et al., 2003). Wider tree rings may form in the current year if the above assumptions are correct. However, higher temperatures in these periods may increase evapotranspiration rates and enhance moisture stress leading to narrow rings. The tree growth-climate relationship of spruce trees established in this research is similar to that obtained in previous research in Central Asia (Zhang et al., 2015; Chen et al., 2017a). It is noteworthy that, in general, the influence of precipitation decreases with increased elevation, while the influence of temperature increases.

Fig. 4

Pearson correlations between tree-ring chronologies and climatic data (precipitation, mean temperature, mean minimum temperature, mean maximum temperature, and PDSI). Lowercase and capitals represent months of the previous and current years, respectively.

The optimal temperature for evergreen conifers to conduct photosynthesis ranges from 10°C to 25°C. If the temperature falls below –3°C to –5°C or raises above 35°C to 42°C, photosynthesis by conifers may cease (Wang, 2000). Fig. 2a reveals that the mean temperature ranged from 4.3°C to 18.3°C in April to October, and from –5.4°C to –4.5°C in March to November, respectively. Therefore, April–October was taken as the growing season for spruces in our study area. The monthly NDVI data (1982–2006) from April to October were utilized to assess the connection between vegetation coverage and tree growth (Fig. 5). Significant correlations were found from July to October based on a 0.05 significance level and peaked in September. By testing various combinations of months, we found that the maximal correlation coefficient was between the ARC chronology and the July–October NDVI (r=0.577, p<0.01, n=25). The relationship between ring width and vegetation coverage in the growing season for our study is similar to that found for other species in previous studies in northern China: the species include Sabina przewalskii Kom. and Sabina tibetica Kom. in the Tibetan Plateau (He and Shao, 2006; Shang et al., 2016), Pinus tabulaeformis Carr. in the Ortindag Sand Land (Liang et al., 2009), and Pinus sylvestris var. mongolica in the Greater Higgnan Mountains (He et al., 2005a).

Fig. 5

Pearson correlation coefficients for the ARC chronology and monthly NDVI from April to October. J–O specifies the period from July to October.

NDVI reconstruction and stability tests

We computed the correlation coefficients between spruces’ growth in a radial pattern within our study area and various assemblages of months for the NDVI for the period 1982–2006, so that the most suitable season for reconstruction could be selected. To do this, we reconstructed the NDVI from July to October by using of the regional chronology. We also used a linear regression model to portray the relationship between the ARC chronology and the NDVI. The model is as follows:

NJuly.Oct.=0.359+0.061×ARC,$$\begin{array}{} \displaystyle {\it N}\rm _{July.-Oct.}=0.359+0.061\times ARC, \end{array}$$

(n=25, r=0.577, R2=33.3%, R2adj=30.4%, SE=0.017, F=11.489, and Durbin-Watson=1.818)

where NJuly.–Oct. is the July–October NDVI for the study area and ARC refers to the regional chronology. The model accounts for 33.3% of the NDVI variance during the calibration period 1982–2006 for the ARC chronology. Fig. 6a shows that the reconstructed NDVI tracks the observations well. The results of the Leave-one-out and Bootstrap tests (100 iterations in the recomputation process) revealed that values of r, R2, R2adj, standard error (SE), F-value (F), and Durbin-Watson (D/W) are nearly equal to those of the original regression model (1) (Table 3). Table 4 shows statistics from the split-sample calibration-verification tests for the reconstructed NDVI series. The values for reduction in error (RE) and coefficient of efficiency (CE) exceed or equal 0. Values of the product mean test (t) are positive, which indicating significant accuracy in the tree-ring estimates.

Fig. 6

(a) Comparison between the reconstructed (red line) and observed (blue line) July–October NDVI. (b) July–October NDVI reconstruction for the Alatau Mountains since 1850 (solid line). The thick line represents the 21-year moving average. Ten red and ten blue dots represent the lowest and highest value years, respectively.

Verification results from the Leave-one-out and Bootstrap tests for the NDVI reconstruction.

StatisticCalibrationLeave-one-out Mean (Range)Bootstrap (100 Iterations) Mean (Range)
r0.580.58 (0.52–0.63)0.57 (0.34–0.77)
Squared multiple correlation (R2)0.330.33 (0.27–0.40)0.34 (0.12–0.60)
Adjusted squared multiple correlation (Radj2$\begin{array}{} \displaystyle {\it R}\rm ^{2}_{adj} \end{array}$)0.300.30 (0.24–0.37)0.31 (0.08–0.58)
Standard error (SE)0.0170.017 (0.015–0.018)0.16 (0.11–0.19)
F-value (F)11.4511.06 (8.10–14.56)13.21 (3.17–32.39)
Durbin-Watson (D/W)1.821.84 (1.60–2.05)1.22 (0.49–1.95)

Statistics of split-sample calibration-verification tests for the NDVI reconstruction.

StatisticCalibration (1982–1994)Verification (1995–2006)Calibration (1994–2006)Verification (1982–1993)Full calibration (1982–2006)
r0.670.410.440.680.58
R20.430.170.190.460.33
Radj2$\begin{array}{} \displaystyle {\it R}\rm ^{2}_{adj} \end{array}$0.38/0.12/0.30
RE/0.10/0.35/
CE/0.09/0.34/
t/3.01/3.49/

The above results, which indicates significant skill in the tree-ring estimates, all demonstrate that the model (1) is stable and reliable. Therefore, we used model 1 to reconstruct the NDVI around the Alatau Mountains for July to October during the period of 1850–2016; an average value of 0.418 and a standard deviation of σ=0.012 (Fig. 6b) were obtained.

Characteristics of the reconstruction

The values for extreme years of the NDVI reconstruction are now listed. The 10 highest value years were 1897 (value: 0.447), 1993 (0.442), 1924 (0.441), 1970 (0.440), 1962 (0.440), 1994 (0.440), 1973 (0.438), 1952 (0.437), 1852 (0.437), and 1922 (0.435), and the 10 lowest-value years were 1945 (0.389), 1927 (0.389), 1917 (0.390), 1879 (0.390), 1866 (0.391), 2008 (0.393), 1944 (0.396), 2000 (0.396), 1938 (0.398), and 1867 (0.399). It can be seen that the difference between the highest-value (1897) and lowest-value (1945) years is 0.058.

The interdecadal variability was conspicuous. When a 21-year moving average was applied to the reconstruction (Fig. 6b), three periods of dense and sparse vegetation coverage could be differentiated. The periods of dense vegetation coverage (above the mean value of the reconstruction) were 1860–1870 (average value: 0.419), 1891–1907 (0.421), and 1950–1974 (0.423). The sparse vegetation coverage periods (below the mean) were 1871–1890 (0.417), 1908–1949 (0.416), and 1975–2006 (0.416).

Discussion

The lowest-value years of 1917 and 1938 in the newly reconstructed NDVI series coincide with the historical documents that record a lack of rainfall in summer and autumn leading to a great drought occurring in the Ili region (Wen et al., 2006), a region near the present study area. In particular, a period of extreme drought occurred during the 1920s in a wide area of northern China, which is mentioned in historical records (Shi et al., 1991). Hydro-climatic reconstructions based on tree rings (Liang et al., 2006; Gou et al., 2007) are also corroborated by the third lowest-value year (1917) in the NDVI reconstruction. This shows that our reconstructed NDVI series can capture the signals of drought disasters for the Alatau Mountains.

Spatial correlation was carried out to evaluate the regional significance of the ring-width based NDVI reconstruction. The NDVI series has a correlations (>0.3) with the gridded PDSI data for July–October for 1901–2015 in a large area between approximately 40°–48°N and 68°–82°E. The highest correlations (>0.4) appear in southern Kazakhstan (Fig. 7). The newly reconstructed NDVI series for the Alatau Mountains (NAM) were compared with four climatic reconstructions based on tree-ring data: (1) August–January standardized precipitation evapotranspiration index for the southeastern Kazakhstan (SKZ, 1785–2014; Chen et al., 2017b); (2) July–August precipitation for Nilka in the Ili region (PNL, 1671–2006; Zhang et al., 2010); (3) July–June precipitation for the Issyk Lake (PIL, 1756–2012; Zhang et al., 2015); and (4) August–April precipitation for the Aksu region (PAK, 1396–2005; Zhang et al., 2009). The locations of these reconstructions are indicated in Fig. 7. There is a good coherence between the standardized NDVI reconstruction and the climatic reconstructions in the original, high-frequency, and low-frequency domains (Fig. 8). Table 5 shows that, although the coherence between the NAM and the PAK reconstructions are relatively weak, all of the correlation coefficients between the RIL and the other climatic reconstructions in the original, high-frequency and low-frequency domains exceed the 0.05 significance level in the common period of 1850–2005. The correlations between the NAM and the other climatic reconstructions weakened with the increasing north to south spatial distance. The best correlation relationship is between the NAM and the SKZ reconstructions and is a result of their proximity. The weakest correlation is between the NAM and the PAK reconstructions, because these areas are the furthest apart. Fig. 8a shows that the NDVI reconstruction and the climatic reconstructions from previous studies corresponded well in the common period. Furthermore, most of the high and low values in these reconstructions appeared in similar years in the high-frequency domain (Fig. 8b). In addition, the variation modes and long-term trends of these reconstructions were roughly synchronous in the low-frequency domain (Fig. 8c).

Fig. 7

Spatial correlation between the reconstructed NDVI series ana gridded July–October PDSI dataset. The solid circle represents the July–October NDVI reconstruction, and the dashed circles represent the four climatic reconstructions (SKZ, PNL, PIL, and PAK) used in the comparison analyses.

Fig. 8

Graphical comparison of the NDVI reconstruction for the Alatau Mountains (black line) and four climatic reconstructions for Central Asia (SKZ: red line; PNL: blue line; PIL: gray line; PAK: green line).

Correlation coefficients between the NDVI and climatic reconstructions over the common period of 1850–2005.

Original(n = 156)High-frequency (n = 144)Low-frequency (n = 144)
SKZPNLPILPAKSKZPNLPILPAKSKZPNLPILPAK
NAM0.5070.460

Significant at p < 0.01

0.371

Significant at p < 0.01

0.197

Significant at p < 0.05

0.580

Significant at p < 0.01

0.508

Significant at p < 0.01

0.332

Significant at p < 0.01

0.189

Significant at p < 0.05

0.498

Significant at p < 0.01

0.326

Significant at p < 0.01

0.452

Significant at p < 0.01

0.190

Significant at p < 0.05

To further analyze the large-scale climate anomalies associated with the newly developed NDVI series, two typical years (1993 and 2008) were investigated. As shown in Fig. 9a, above normal rainfall was evident around the Alatau Mountains, which accounts for the highest NDVI in 1993. An anomalous westerly wind was pronounced over Eurasia between 40°–60°N, which brought air rich in water vapor from the Atlantic Ocean, and the Mediterranean and Caspian Seas. As a result, a rainy belt was observed from Western Europe to the Alatau Mountains. In contrast, below normal rainfall was evident around the Alatau Mountains in 2008 (Fig. 9b), which is well represented by the associated lowest NDVI since the 1950s. Unlike the effects of the remarkable zonal water vapor transport in 1993, precipitation in the Alatau Mountains is normally controlled by dry air from high latitudes (He et al., 2005b).

Fig. 9

Patterns of precipitation anomalies (shading; mm/month) and horizontal wind anomalies at 850hPa (vector; m/s) in (a) the high NDVI year of 1993, and (b) the low NDVI year of 2008. The red block indicates the location of the present study.

Conclusions

A 167-year, regional, ring-width chronology was developed from 93 living and healthy spruces on the northern slopes of the Alatau Mountains. The results of correlation analyses revealed that more precipitation and lower temperatures in the growth season of the previous and current years might help to form a wider tree ring in the target coniferous species. Because of the coherence between ring width and vegetation coverage under the same water limitations, a July–October NDVI reconstruction for the Alatau Mountains was developed using the regional chronology. The reconstruction matched the observed data well, and its lowest-value years precisely captured drought events mentioned in historical documents. A spatial analysis and correlations with other tree-ring based climatic reconstructions indicated that the newly reconstructed NDVI series was affected by large-scale climatic oscillations.

The preliminary study in this paper provides new knowledge about the variations in vegetation coverage for the Alatau Mountains and will be helpful in improving the tree-ring dataset for Central Asia. However, this study is limited because it is local research based on only three sampling areas. Thus, further studies should be carried out at a larger spatial scale to verify the results presented in this paper.

Fig. 1

(a) Location of the study area; (b) map of the study area; and (c) sampling sites: Amanboktek (AMA), Kikbay (KBZ), and Basika (BSK). The selected NDVI grid (44.932°–46.023°N, 78.986°–80.950°E) is represented by the dash-block diagram.
(a) Location of the study area; (b) map of the study area; and (c) sampling sites: Amanboktek (AMA), Kikbay (KBZ), and Basika (BSK). The selected NDVI grid (44.932°–46.023°N, 78.986°–80.950°E) is represented by the dash-block diagram.

Fig. 2

Gridded climatic and NDVI data for the Alatau Mountain: (a) Precipitation, mean temperature, mean maximum temperature and mean minimum temperature per month. Bars indicate precipitation (in mm), and curves with different colors represent mean temperature, mean minimum temperature and mean maximum temperature (in °C); (b) Monthly PDSI; (c) Monthly NDVI; (d) Annual precipitation during 1901–2015; (e–g) Annual mean temperature, annual mean minimum temperature and annual mean maximum temperature during 1901–2015; (h) Annual mean PDSI during 1901–2015; and (i) Annual mean NDVI during 1982–2006. Dashed lines represent trends in climate and NDVI data.
Gridded climatic and NDVI data for the Alatau Mountain: (a) Precipitation, mean temperature, mean maximum temperature and mean minimum temperature per month. Bars indicate precipitation (in mm), and curves with different colors represent mean temperature, mean minimum temperature and mean maximum temperature (in °C); (b) Monthly PDSI; (c) Monthly NDVI; (d) Annual precipitation during 1901–2015; (e–g) Annual mean temperature, annual mean minimum temperature and annual mean maximum temperature during 1901–2015; (h) Annual mean PDSI during 1901–2015; and (i) Annual mean NDVI during 1982–2006. Dashed lines represent trends in climate and NDVI data.

Fig. 3

Chronologies (thin lines) with their sample depth (dashed lines). Red and blue lines represent EPS and Rbar data: (a) AMA; (b) KBZ; (c) BSK; and (d) ARC.
Chronologies (thin lines) with their sample depth (dashed lines). Red and blue lines represent EPS and Rbar data: (a) AMA; (b) KBZ; (c) BSK; and (d) ARC.

Fig. 4

Pearson correlations between tree-ring chronologies and climatic data (precipitation, mean temperature, mean minimum temperature, mean maximum temperature, and PDSI). Lowercase and capitals represent months of the previous and current years, respectively.
Pearson correlations between tree-ring chronologies and climatic data (precipitation, mean temperature, mean minimum temperature, mean maximum temperature, and PDSI). Lowercase and capitals represent months of the previous and current years, respectively.

Fig. 5

Pearson correlation coefficients for the ARC chronology and monthly NDVI from April to October. J–O specifies the period from July to October.
Pearson correlation coefficients for the ARC chronology and monthly NDVI from April to October. J–O specifies the period from July to October.

Fig. 6

(a) Comparison between the reconstructed (red line) and observed (blue line) July–October NDVI. (b) July–October NDVI reconstruction for the Alatau Mountains since 1850 (solid line). The thick line represents the 21-year moving average. Ten red and ten blue dots represent the lowest and highest value years, respectively.
(a) Comparison between the reconstructed (red line) and observed (blue line) July–October NDVI. (b) July–October NDVI reconstruction for the Alatau Mountains since 1850 (solid line). The thick line represents the 21-year moving average. Ten red and ten blue dots represent the lowest and highest value years, respectively.

Fig. 7

Spatial correlation between the reconstructed NDVI series ana gridded July–October PDSI dataset. The solid circle represents the July–October NDVI reconstruction, and the dashed circles represent the four climatic reconstructions (SKZ, PNL, PIL, and PAK) used in the comparison analyses.
Spatial correlation between the reconstructed NDVI series ana gridded July–October PDSI dataset. The solid circle represents the July–October NDVI reconstruction, and the dashed circles represent the four climatic reconstructions (SKZ, PNL, PIL, and PAK) used in the comparison analyses.

Fig. 8

Graphical comparison of the NDVI reconstruction for the Alatau Mountains (black line) and four climatic reconstructions for Central Asia (SKZ: red line; PNL: blue line; PIL: gray line; PAK: green line).
Graphical comparison of the NDVI reconstruction for the Alatau Mountains (black line) and four climatic reconstructions for Central Asia (SKZ: red line; PNL: blue line; PIL: gray line; PAK: green line).

Fig. 9

Patterns of precipitation anomalies (shading; mm/month) and horizontal wind anomalies at 850hPa (vector; m/s) in (a) the high NDVI year of 1993, and (b) the low NDVI year of 2008. The red block indicates the location of the present study.
Patterns of precipitation anomalies (shading; mm/month) and horizontal wind anomalies at 850hPa (vector; m/s) in (a) the high NDVI year of 1993, and (b) the low NDVI year of 2008. The red block indicates the location of the present study.

Correlation coefficients for three chronologies over the common period of 1910–2016.

Original (n=107)High-frequency (n=95)Low frequency (n=95)
AMAKBZBSKAMAKBZBSKAMAKBZBSK
AMA10.783

Significant at p<0.001.

0.660

Significant at p<0.001.

10.842

Significant at p<0.001.

0.682

Significant at p<0.001.

10.641

Significant at p<0.001.

0.683

Significant at p<0.001.

KBZ/10.624

Significant at p<0.001.

/10.759

Significant at p<0.001.

/10.388

Significant at p<0.001.

BSK//1//1//1

Statistics of split-sample calibration-verification tests for the NDVI reconstruction.

StatisticCalibration (1982–1994)Verification (1995–2006)Calibration (1994–2006)Verification (1982–1993)Full calibration (1982–2006)
r0.670.410.440.680.58
R20.430.170.190.460.33
Radj2$\begin{array}{} \displaystyle {\it R}\rm ^{2}_{adj} \end{array}$0.38/0.12/0.30
RE/0.10/0.35/
CE/0.09/0.34/
t/3.01/3.49/

Statistical characteristics of chronologies over the common period of 1900–2010

StatisticAMAKBZBSKARC
Standard deviation (SD)0.1980.2240.4140.208
Mean sensitivity (MS)0.1800.1860.3790.180
First-order autocorrelation (AC1)0.3960.4960.5770.458
Interseries correlation (trees)0.4170.4840.7380.418
Interseries correlation (all series)0.4310.4880.7530.422
Mean within-tree correlation0.7860.8050.9610.808
Signal-to-noise ratio (SNR)9.85862.01218.29362.025
Expressed population signal (EPS)0.9080.9840.9480.984
The first principal component (PC#1)0.4920.5180.7990.461
First year EPS >0.851910185018801850

Verification results from the Leave-one-out and Bootstrap tests for the NDVI reconstruction.

StatisticCalibrationLeave-one-out Mean (Range)Bootstrap (100 Iterations) Mean (Range)
r0.580.58 (0.52–0.63)0.57 (0.34–0.77)
Squared multiple correlation (R2)0.330.33 (0.27–0.40)0.34 (0.12–0.60)
Adjusted squared multiple correlation (Radj2$\begin{array}{} \displaystyle {\it R}\rm ^{2}_{adj} \end{array}$)0.300.30 (0.24–0.37)0.31 (0.08–0.58)
Standard error (SE)0.0170.017 (0.015–0.018)0.16 (0.11–0.19)
F-value (F)11.4511.06 (8.10–14.56)13.21 (3.17–32.39)
Durbin-Watson (D/W)1.821.84 (1.60–2.05)1.22 (0.49–1.95)

Correlation coefficients between the NDVI and climatic reconstructions over the common period of 1850–2005.

Original(n = 156)High-frequency (n = 144)Low-frequency (n = 144)
SKZPNLPILPAKSKZPNLPILPAKSKZPNLPILPAK
NAM0.5070.460

Significant at p < 0.01

0.371

Significant at p < 0.01

0.197

Significant at p < 0.05

0.580

Significant at p < 0.01

0.508

Significant at p < 0.01

0.332

Significant at p < 0.01

0.189

Significant at p < 0.05

0.498

Significant at p < 0.01

0.326

Significant at p < 0.01

0.452

Significant at p < 0.01

0.190

Significant at p < 0.05

Barbaroux C, Breda N and Dufrene E, 2003. Distribution of above-ground and below-ground carbohydrate reserves in adult trees of two contrasting broad-leaved species (Quercus petraea and Fagus sylvatica). New Phytologist 157(3): 605–615, 10.1046/j.1469-8137.2003.00681.x.BarbarouxCBredaNDufreneE2003Distribution of above-ground and below-ground carbohydrate reserves in adult trees of two contrasting broad-leaved species (Quercus petraea and Fagus sylvatica)New Phytologist157360561510.1046/j.1469-8137.2003.00681.xOpen DOISearch in Google Scholar

Beck PSA, Andreu-Hayles L, D’Arrigo R, Anchukaitis KJ, Tucker CJ, Pinzón JE and Goetz SJ, 2013. A large-scale coherent signal of canopy status in maximum latewood density of tree rings at arctic treeline in North America. Global and Planetary Change 100: 109–118, 10.1016/j.gloplacha.2012.10.005.BeckPSAAndreu-HaylesLD’ArrigoRAnchukaitisKJTuckerCJPinzónJEGoetzSJ2013A large-scale coherent signal of canopy status in maximum latewood density of tree rings at arctic treeline in North AmericaGlobal and Planetary Change10010911810.1016/j.gloplacha.2012.10.005Open DOISearch in Google Scholar

Bräuning A, 1994. Dendrochronology for the last 1400 years in Eastern Tibet. GeoJournal 34(1): 75–95, 10.1007/BF00813972.BräuningA1994Dendrochronology for the last 1400 years in Eastern TibetGeoJournal341759510.1007/BF00813972Open DOISearch in Google Scholar

Chen F, Yuan YJ, Wei WS, Wang L, Yu SL, Zhang RB, Fan ZA, Shang HM, Zhang TW and Li Y, 2012. Tree ring density-based summer temperature reconstruction for Zajsan Lake area, East Kazakhstan. International Journal of Climatology 32: 1089–1097, 10.1002/joc.2327.ChenFYuanYJWeiWSWangLYuSLZhangRBFanZAShangHMZhangTWLiY2012Tree ring density-based summer temperature reconstruction for Zajsan Lake area, East KazakhstanInternational Journal of Climatology321089109710.1002/joc.2327Open DOISearch in Google Scholar

Chen F, He Q, Bakytbek E, Yu SL and Zhang RB, 2017a. Reconstruction of a long streamflow record using tree rings in the upper Kurshab River (Pamir-Alai Mountains) and its application to water resources management. International Journal of Water Resources Development 33(6): 976–986, 10.1080/07900627.2016.1238347.ChenFHeQBakytbekEYuSLZhangRB2017aReconstruction of a long streamflow record using tree rings in the upper Kurshab River (Pamir-Alai Mountains) and its application to water resources managementInternational Journal of Water Resources Development33697698610.1080/07900627.2016.1238347Open DOISearch in Google Scholar

Chen F, Mambetov B, Maisupova B and Kelgenbayev N, 2017b. Drought variations in Almaty (Kazakhstan) since AD 1785 based on spruce tree rings. Stochastic Environmental Research and Risk Assessment 31(8): 2097–2105, 10.1007/s00477-016-1290-y.ChenFMambetovBMaisupovaBKelgenbayevN2017bDrought variations in Almaty (Kazakhstan) sincen AD, 1785 based on spruce tree ringsStochastic Environmental Research and Risk Assessment3182097210510.1007/s00477-016-1290-yOpen DOISearch in Google Scholar

Cook ER, 1985. A time-series analysis approach to tree-ring standardization. PhD dissertation, The University of Arizona.CookER1985A time-series analysis approach to tree-ring standardizationPhD dissertationThe University of ArizonaSearch in Google Scholar

Cook ER, Meko DM, Stahle DW and Cleaveland MK, 1999. Drought reconstructions for the continental United States. Journal of Climate 12(4): 1145–1162, 10.1175/1520-0442(1999)012<1145:DRFTCU>2.0.CO;2.CookERMekoDMStahleDWCleavelandMK1999Drought reconstructions for the continental United StatesJournal of Climate1241145116210.1175/1520-0442(1999)012<1145:DRFTCU>2.0.CO;2Open DOISearch in Google Scholar

Cook ER and Krusic PJ, 2005. Program ARSTAN: A Tree-Ring Standardization Program Based on Detrending and Autoregressive Time Series Modeling, with Interactive Graphics. Lamont-Doherty Earth Observatory, Columbia University, Palisades, New York.CookERKrusicPJ2005Program ARSTAN: A Tree-Ring Standardization Program Based on Detrending and Autoregressive Time Series Modeling, with Interactive GraphicsLamont-Doherty Earth ObservatoryColumbia UniversityPalisades, New YorkSearch in Google Scholar

D’Arrigo RD, Malmstrom CM, Jacoby GC, Los SO and Bunker DE, 2000. Correlation between maximum latewood density of annual tree rings and NDVI based estimates of forest productivity. International Journal of Remote Sensing 21: 2329–2336, 10.1080/01431160050029611.D’ArrigoRDMalmstromCMJacobyGCLosSOBunkerDE2000Correlation between maximum latewood density of annual tree rings and NDVI based estimates of forest productivityInternational Journal of Remote Sensing212329233610.1080/01431160050029611Open DOISearch in Google Scholar

Esper J, 2000. Long-term tree-ring variations in Juniperus at the upper timber-line in the Karakorum (Pakistan). Holocene 10: 253–260, 10.1191/095968300670152685.EsperJ2000Long-term tree-ring variations in Juniperus at the upper timber-line in the Karakorum (Pakistan)Holocene1025326010.1191/095968300670152685Open DOISearch in Google Scholar

Esper J, Shiyatov SG, Mazepa VS, Wilson RJS, Graybill DA and Funkhouser G, 2003. Temperature-sensitive Tien Shan tree ring chronologies show multi-centennial growth trends. Climate Dynamics 21: 699–706, 10.1007/s00382-003-0356-y.EsperJShiyatovSGMazepaVSWilsonRJSGraybillDAFunkhouserG2003Temperature-sensitive Tien Shan tree ring chronologies show multi-centennial growth trendsClimate Dynamics2169970610.1007/s00382-003-0356-yOpen DOISearch in Google Scholar

Esper J, Frank DC, Wilson RJS, Büntgen U and Treydte K, 2007. Uniform growth trends among central Asian low- and high-elevation juniper tree sites. Trees-Structure and Function 21: 141–150, 10.1007/s00468-006-0104-0.EsperJFrankDCWilsonRJSBüntgenUTreydteK2007Uniform growth trends among central Asian low- and high-elevation juniper tree sitesTrees-Structure and Function2114115010.1007/s00468-006-0104-0Open DOISearch in Google Scholar

Gamon JA, Kovalchuck O, Wong CYS, Harris A and Garrity SR, 2015. Monitoring seasonal and diurnal changes in photosynthetic pigments with automated PRI and NDVI sensors. Biogeosciences 12: 4149–4159, 10.5194/bgd-12-2947-2015.GamonJAKovalchuckOWongCYSHarrisAGarritySR2015Monitoring seasonal and diurnal changes in photosynthetic pigments with automated PRI, NDVI sensorsBiogeosciences124149415910.5194/bgd-12-2947-2015Open DOISearch in Google Scholar

Gou XH, Chen FH, Cook ER, Gordon J, Yang MX and Li JB, 2007. Streamflow variations of the Yellow River over the past 593 years in western China reconstructed from tree rings. Water Resources Research 43(6): 113–119, 10.1029/2006WR005705.GouXHChenFHCookERGordonJYangMXLiJB2007Streamflow variations of the Yellow River over the past 593 years in western China reconstructed from tree ringsWater Resources Research43611311910.1029/2006WR005705Open DOISearch in Google Scholar

Grissino-Mayer HD, 2001. Evaluating crossdating accuracy: a manual and tutorial for the computer program COFECHA. Tree-Ring Research 57: 205–221.Grissino-MayerHD2001Evaluating crossdating accuracy: a manual and tutorial for the computer program COFECHATree-Ring Research57205221Search in Google Scholar

He JC, Wang LL and Shao XM, 2005a. The relationships between mongolian scotch pine tree ring indices and normalized difference vegetation index in Mohe, China. Quaternary Sciences 25: 252–257 (in Chinese, with English abstract).HeJCWangLLShaoXM2005aThe relationships between mongolian scotch pine tree ring indices and normalized difference vegetation index in Mohe, ChinaQuaternary Sciences25252257(in Chinese, with English abstract)Search in Google Scholar

He JH, Liu YY and Chang Y, 2005b. Analysis of summer precipitation anomaly and the feature of water vapor transport and circulation in Northwest China. Arid Meteorology 23(1): 10–16 (in Chinese, with English abstract).HeJHLiuYYChangY2005bAnalysis of summer precipitation anomaly and the feature of water vapor transport and circulation in Northwest ChinaArid Meteorology2311016(in Chinese, with English abstract)Search in Google Scholar

He JC and Shao XM, 2006. Relationships between tree-ring width index and NDVI of grassland in Delingha. Chinese Science Bulletin 51(9): 1106–1114, 10.1007/s11434-006-1106-4.HeJCShaoXM2006Relationships between tree-ring width index and NDVI of grassland in DelinghaChinese Science Bulletin5191106111410.1007/s11434-006-1106-4Open DOISearch in Google Scholar

He Y, Yang TB, Ji Q and Chen J, 2014. Glacier shrinkage and its dependence on climate in the Central Asia Alatau region. Arid Land Geography 37(5): 908–915 (in Chinese, with English abstract).HeYYangTBJiQChenJ2014Glacier shrinkage and its dependence on climate in the Central Asia Alatau regionArid Land Geography375908915(in Chinese, with English abstract)Search in Google Scholar

Hmimina G, Dufrêne E, Pontailler JY, Delpierre N, Aubinet M, Caquet B, Grandcourt A, Burban B, Flechard C, Granier A, Gross P, Heinesch B, Longdoz B, Moureaux C, Ourcival JM, Rambal S, Saint André L and Soudani K, 2013. Evaluation of the potential of MODIS satellite data to predict vegetation phenology in different biomes: An investigation using ground-based NDVI measurements. Remote Sensing of Environment 132: 145–158, 10.1016/j.rse.2013.01.010.HmiminaGDufrêneEPontaillerJYDelpierreNAubinetMCaquetBGrandcourtABurbanBFlechardCGranierAGrossPHeineschBLongdozBMoureauxCOurcivalJMRambalSSaint AndréLSoudaniK2013Evaluation of the potential of MODIS satellite data to predict vegetation phenology in different biomes: An investigation using ground-based NDVI measurementsRemote Sensing of Environment13214515810.1016/j.rse.2013.01.010Open DOISearch in Google Scholar

Jiapaer G, Liang S, Yi Q and Liu J, 2015. Vegetation dynamics and responses to recent climate change in Xinjiang using leaf area index as an indicator. Ecological Indicators 58: 64–76, 10.1016/j.ecolind.2015.05.036.JiapaerGLiangSYiQLiuJ2015Vegetation dynamics and responses to recent climate change in Xinjiang using leaf area index as an indicatorEcological Indicators58647610.1016/j.ecolind.2015.05.036Open DOISearch in Google Scholar

Kalnay E, Kanamitsu M, Kistler R, Collins W, Deaven D, Gandin L, Iredell M, Saha S, White G, Woollen J, Zhu Y, Leetmaa A, Reynolds B, Chelliah M, Ebisuzaki W, Higgins W, Janowiak J, Mo KC, Ropelewski C, Wang J, Jenne R and Joseph D, 1996. The NCEP/NCAR 40-year reanalysis project. Bulletin of the American Meteorological Society 77(3): 437–472, 10.1175/1520-0477(1996)077<0437:TNYRP>2.0.CO;2.KalnayEKanamitsuMKistlerRCollinsWDeavenDGandinLIredellMSahaSWhiteGWoollenJZhuYLeetmaaAReynoldsBChelliahMEbisuzakiWHigginsWJanowiakJMoKCRopelewskiCWangJJenneRJosephD1996The NCEP/NCAR 40-year reanalysis projectBulletin of the American Meteorological Society77343747210.1175/1520-0477(1996)077<0437:TNYRP>2.0.CO;2Open DOISearch in Google Scholar

Kaufmann RK, D’Arrigo RD, Laskowski C, Myneni RB, Zhou L and Davi NK, 2004. The effect of growing season and summer greenness on northern forests. Geophysical Research Letters 31: L09205, 10.1029/2004GL019608.KaufmannRKD’ArrigoRDLaskowskiCMyneniRBZhouLDaviNK2004The effect of growing season and summer greenness on northern forestsGeophysical Research Letters31L0920510.1029/2004GL019608Open DOISearch in Google Scholar

Leavitt SW, Chase TN, Rajagopalan B, Lee E and Lawrence PJ, 2008. Southwestern U.S. tree-ring carbon isotope indices as a possible proxy for reconstruction of greenness of vegetation. Geophysical Research Letters 35: L12704, 10.1029/2008GL033894.LeavittSWChaseTNRajagopalanBLeeELawrencePJ2008Southwestern U.S. tree-ring carbon isotope indices as a possible proxy for reconstruction of greenness of vegetationGeophysical Research Letters35L1270410.1029/2008GL033894Open DOISearch in Google Scholar

Lebourgeois F, Cousseau G and Ducos Y, 2004. Climate-tree growth relationships of Quercus petraea Mill. stand in the Forest of Bercé (“Futaie des Clos”, Sarthe, France). Annals of Forest Science 61(4): 361–372, 10.1051/forest:2004029.LebourgeoisFCousseauGDucosY2004Climate-tree growth relationships of Quercus petraea Mill. stand in the Forest of Bercé (“Futaie des Clos”, Sarthe, France)Annals of Forest Science61436137210.1051/forest:2004029Open DOISearch in Google Scholar

Liang E, Liu X, Yuan Y, Qin N, Fang X, Huang L, Zhu H, Wang L and Shao X, 2006. The 1920s drought recorded by tree rings and historical documents in the semi-arid and arid areas of northern China. Climatic Change 79(3): 403–432, 10.1007/s10584-006-9082-x.LiangELiuXYuanYQinNFangXHuangLZhuHWangLShaoX2006The 1920s drought recorded by tree rings and historical documents in the semi-arid and arid areas of northern ChinaClimatic Change79340343210.1007/s10584-006-9082-xOpen DOISearch in Google Scholar

Liang EY, Eckstein D and Liu HY, 2009. Assessing the recent grassland greening trend in a long-term context based on tree-ring analysis: A case study in North China. Ecological Indicators 9: 1280–1283, 10.1016/j.ecolind.2009.02.007.LiangEYEcksteinDLiuHY2009Assessing the recent grassland greening trend in a long-term context based on tree-ring analysis: A case study in North ChinaEcological Indicators91280128310.1016/j.ecolind.2009.02.007Open DOISearch in Google Scholar

Liu Y, Wang CY, Hao WJ, Song HM, Cai QF, Tian H, Sun B and Linderholm HW, 2011. Tree-ring-based annual precipitation reconstruction in Kalaqin, Inner Mongolia for the last 238 years. Chinese Science Bulletin 56: 2995–3002, 10.1007/s11434-011-4706-6.LiuYWangCYHaoWJSongHMCaiQFTianHSunBLinderholmHW2011Tree-ring-based annual precipitation reconstruction in Kalaqin, Inner Mongolia for the last 238 yearsChinese Science Bulletin562995300210.1007/s11434-011-4706-6Open DOISearch in Google Scholar

Lopatin E, Kolström T and Spiecker H, 2006. Determination of forest growth trends in Komi Republic (northwestern Russia): combination of tree-ring analysis and remote sensing data. Boreal Environment Research 11: 341–353.LopatinEKolströmTSpieckerH2006Determination of forest growth trends in Komi Republic (northwestern Russia): combination of tree-ring analysis and remote sensing dataBoreal Environment Research11341353Search in Google Scholar

Malmström CM, Thompson MV, Juday G, Los SO, Randerson JT and Field CB, 1997. Interannual variation in global-scale net primary production: testing model estimates. Global Biogeochemical Cycles 11: 367–392, 10.1029/97GB01419.MalmströmCMThompsonMVJudayGLosSORandersonJTFieldCB1997Interannual variation in global-scale net primary production: testing model estimatesGlobal Biogeochemical Cycles1136739210.1029/97GB01419Open DOISearch in Google Scholar

Meko DM and Graybill DA, 1995. Tree-ring reconstruction of Upper Gila River discharge. Journal of the American Water Resources Association 31(4): 605–616, 10.1111/j.1752-1688.1995.tb03388.x.MekoDMGraybillDA1995Tree-ring reconstruction of Upper Gila River dischargeJournal of the American Water Resources Association31460561610.1111/j.1752-1688.1995.tb03388.xOpen DOISearch in Google Scholar

Michaelsen J, 1987. Cross-validation in statistical climate forecast models. Journal of Applied Meteorology 26(11): 1589–1600, 10.1175/1520-0450(1987)026<1589:CVISCF>2.0.CO;2.MichaelsenJ1987Cross-validation in statistical climate forecast modelsJournal of Applied Meteorology26111589160010.1175/1520-0450(1987)026<1589:CVISCF>2.0.CO;2Open DOISearch in Google Scholar

Nagy L, 2006. European high mountain (alpine) vegetation and its suitability for indicating climate change impacts. Biology and Environment 106B(3): 335–341, 10.3318/BIOE.2006.106.3.335.NagyL2006European high mountain (alpine) vegetation and its suitability for indicating climate change impactsBiology and Environment106B333534110.3318/BIOE.2006.106.3.335Open DOISearch in Google Scholar

New M, Hulme M and Jones P, 2000. Representing twentieth-century space–time climate variability. Part II: Development of 1901–96 monthly grids of terrestrial surface climate. Journal of Climate 13: 2217–2238, 10.1175/1520-0442(2000)013<2217:RTCSTC>2.0.CO;2.NewMHulmeMJonesP2000Representing twentieth-century space–time climate variability. Part II: Development of 1901–96 monthly grids of terrestrial surface climateJournal of Climate132217223810.1175/1520-0442(2000)013<2217:RTCSTC>2.0.CO;2Open DOISearch in Google Scholar

Opała M, Niedźwiedź T, Rahmonov O, Owczarek P and Małarzewski Ł, 2017. Towards improving the Central Asian dendrochronological network—New data from Tajikistan, Pamir-Alay. Dendrochronologia 41: 10–23, 10.1016/j.dendro.2016.03.006.OpałaMNiedźwiedźTRahmonovOOwczarekPMałarzewskiŁ2017Towards improving the Central Asian dendrochronological network—New data from Tajikistan, Pamir-AlayDendrochronologia41102310.1016/j.dendro.2016.03.006Open DOISearch in Google Scholar

Pauli H, Gottfried M, Reiter K and Grabherr G, 2002. High Mountain Summits as Sensitive Indicators of Climate Change Effects on Vegetation Patterns: The “Multi Summit-Approach” of GLORIA (Global Observation Research Initiative in Alpine Environments). Advances in Global Change Research 9: 45–51, 10.1007/0-306-48051-4_6.PauliHGottfriedMReiterKGrabherrG2002High Mountain Summits as Sensitive Indicators of Climate Change Effects on Vegetation Patterns: The “Multi Summit-Approach” of GLORIA (Global Observation Research Initiative in Alpine Environments)Advances in Global Change Research9455110.1007/0-306-48051-4_6Open DOISearch in Google Scholar

Pettorelli N, Vik JO, Mysterud A, Gaillard J, Tucker CJ and Stenseth NC, 2005. Using the satellite-derived NDVI to assess ecological responses to environmental change. Trends in Ecology & Evolution 20(9): 503–510, 10.1016/j.tree.2005.05.011.PettorelliNVikJOMysterudAGaillardJTuckerCJStensethNC2005Using the satellite-derived NDVI to assess ecological responses to environmental changeTrends in Ecology & Evolution20950351010.1016/j.tree.2005.05.011Open DOISearch in Google Scholar

Schrier GVD, Barichivich J, Briffa KR and Jones PD, 2013. A scPDSI-based global data set of dry and wet spells for 1901–2009. Journal of Geophysical Research Atmospheres 118(10): 4025–4048, 10.1002/jgrd.50355.SchrierGVDBarichivichJBriffaKRJonesPD2013A scPDSI-based global data set of dry and wet spells for 1901–2009Journal of Geophysical Research Atmospheres118104025404810.1002/jgrd.50355Open DOISearch in Google Scholar

Seim A, Tulyaganov T, Omurova G, Nikolyai L, Botman E and Linderholm HW, 2016. Dendroclimatological potential of three juniper species from the Turkestan range, northwestern Pamir-Alay Mountains, Uzbekistan. Trees 30(3): 733–748, 10.1007/s00468-015-1316-y.SeimATulyaganovTOmurovaGNikolyaiLBotmanELinderholmHW2016Dendroclimatological potential of three juniper species from the Turkestan range, northwestern Pamir-Alay Mountains, UzbekistanTrees30373374810.1007/s00468-015-1316-yOpen DOISearch in Google Scholar

Shang HM, Wei WS, Yuan YJ, Yu SL, Zhang RB, Hong JC, Chen F, Zhang TW and Fan ZA, 2016. Normalized vegetation variation index reconstruction based on the tree-ring width in central Tibet. Journal of Lanzhou University: Natural Sciences 52(1): 18–24 (in Chinese, with English abstract).ShangHMWeiWSYuanYJYuSLZhangRBHongJCChenFZhangTWFanZA2016Normalized vegetation variation index reconstruction based on the tree-ring width in central TibetJournal of Lanzhou University: Natural Sciences5211824(in Chinese, with English abstract)Search in Google Scholar

Shi FC, Wang GA, Gao ZD, Mu P and Ma GA, 1991. Recurrence probability of 11-year continuous low water period (1922–1932 A.D.) in the Yellow River. Advances in Water Science 2(4): 258–263 (in Chinese, with English abstract).ShiFCWangGAGaoZDMuPMaGA1991Recurrence probability of 11-year continuous low water period (1922–1932 A.D.) in the Yellow RiverAdvances in Water Science24258263(in Chinese, with English abstract)Search in Google Scholar

Speer JH, 2010. Fundamentals of Tree-ring Research. The University of Arizona Press, Tucson.SpeerJH2010Fundamentals of Tree-ring ResearchThe University of Arizona PressTucsonSearch in Google Scholar

Wang Z, 2000. Phytophysiology. China Agriculture Press, Beijing (in Chinese).WangZ2000PhytophysiologyChina Agriculture PressBeijing(in Chinese)Search in Google Scholar

Wang YJ, Lu RJ, Ma YZ, Meng HW and Gao SY, 2014. Response to climate change of different tree species and NDVI variation since 1923 in the middle arid region of Ningxia, China. Sciences in Cold and Arid Regions 6(1): 30–36, 10.3724/SP.J.1226.2014.00030.WangYJLuRJMaYZMengHWGaoSY2014Response to climate change of different tree species and NDVI variation since 1923 in the middle arid region of Ningxia, ChinaSciences in Cold and Arid Regions61303610.3724/SP.J.1226.2014.00030Open DOISearch in Google Scholar

Wells N, Goddard S and Hayes MJ, 2004. A self-calibrating Palmer Drought Severity Index. Journal of Climate 17: 2335–2351, 10.1175/1520-0442(2004)017<2335:ASPDSI>2.0.CO;2.WellsNGoddardSHayesMJ2004A self-calibrating Palmer Drought Severity IndexJournal of Climate172335235110.1175/1520-0442(2004)017<2335:ASPDSI>2.0.CO;2Open DOISearch in Google Scholar

Wen KG, Shi YG and Ren YY, 2006. Documents of Chinese Meteorological Disaster: Volume of Xinjiang. Meteorological Publishers, Beijing (in Chinese).WenKGShiYGRenYY2006Documents of Chinese Meteorological Disaster: Volume of XinjiangMeteorological PublishersBeijing(in Chinese)Search in Google Scholar

Wigley TML, Briffa KR and Jones PD, 1984. On the average value of correlated time series, with applications in dendroclimatology and hydrometeorology. Journal of Climate and Applied Meteorology 23(2): 201–213, 10.1175/1520-0450(1984)023<0201:OTAVOC>2.0.CO;2.WigleyTMLBriffaKRJonesPD1984On the average value of correlated time series, with applications in dendroclimatology and hydrometeorologyJournal of Climate and Applied Meteorology23220121310.1175/1520-0450(1984)023<0201:OTAVOC>2.0.CO;2Open DOISearch in Google Scholar

Winter MB, WolV B, Gottschling H and Cherubini P, 2009. The impact of climate on radial growth and nut production of Persian walnut (Juglans regia L.) in Southern Kyrgyzstan. European Journal of Forest Research 128: 531–542, 10.1007/s10342-009-0295-1.WinterMBWolVBGottschlingHCherubiniP2009The impact of climate on radial growth and nut production of Persian walnut (Juglans regia L.) in Southern KyrgyzstanEuropean Journal of Forest Research12853154210.1007/s10342-009-0295-1Open DOISearch in Google Scholar

Yin G, Meng XY, Wang H, Hu ZY and Sun ZQ, 2017. Spatial-temporal variation of vegetation and its correlation with climate change in Central Asia during the period of 1982–2012. Acta Ecologica Sinica 37(9): 3149–3163 (in Chinese, with English abstract).YinGMengXYWangHHuZYSunZQ2017Spatial-temporal variation of vegetation and its correlation with climate change in Central Asia during the period of 1982–2012Acta Ecologica Sinica37931493163(in Chinese, with English abstract)10.5846/stxb201601240164Search in Google Scholar

Young GA, 1994. Bootstrap: more than a stab in the dark. Statistical Science 9: 382–415, 10.1214/ss/1177010383.YoungGA1994Bootstrap: more than a stab in the darkStatistical Science938241510.1214/ss/1177010383Open DOISearch in Google Scholar

Yuan YJ, Zhang TW, Wei WS, Nievergelt D, Verstege A, Yu SL, Zhang RB and Esper J, 2013. Development of tree-ring maximum latewood density chronologies for the western Tien Shan Mountains, China: Influence of detrending method and climate response. Dendrochronologia 31: 192–197, 10.1016/j.dendro.2013.05.004.YuanYJZhangTWWeiWSNievergeltDVerstegeAYuSLZhangRBEsperJ2013Development of tree-ring maximum latewood density chronologies for the western Tien Shan Mountains, China: Influence of detrending method and climate responseDendrochronologia3119219710.1016/j.dendro.2013.05.004Open DOISearch in Google Scholar

Zhang RB, Wei WS, Yuan YJ, Yu SL and Yang Q, 2009. A precipitation series of A.D. 1396–2005 in Aksu River Basin on the southern slopes of Tianshan Mountains: reconstruction and analysis. Journal of Glaciology and Geocryology 31: 27–33 (in Chinese, with English abstract).ZhangRBWeiWSYuanYJYuSLYangQ2009A precipitation series of A.D. 1396–2005 in Aksu River Basin on the southern slopes of Tianshan Mountains: reconstruction and analysisJournal of Glaciology and Geocryology312733(in Chinese, with English abstract)Search in Google Scholar

Zhang L, Yuan YJ, Wei WS and Qin L, 2010. Reconstruction and analysis of the 336a July and August precipitation series in Nilka County, Xinjing. Journal of Glaciology and Geocryology 32(5): 914–920 (in Chinese, with English abstract).ZhangLYuanYJWeiWSQinL2010Reconstruction and analysis of the 336a July and August precipitation series in Nilka County, XinjingJournal of Glaciology and Geocryology325914920(in Chinese, with English abstract)Search in Google Scholar

Zhang TW, Yuan YJ, He Q, Wei WS, Diushen M, Shang HM and Zhang RB, 2014. Development of tree-ring width chronologies and tree-growth response to climate in the mountains surrounding the Issyk-Kul Lake, Central Asia. Dendrochronologia 32: 230–236, 10.1016/j.dendro.2014.03.002.ZhangTWYuanYJHeQWeiWSDiushenMShangHMZhangRB2014Development of tree-ring width chronologies and tree-growth response to climate in the mountains surrounding the Issyk-Kul Lake, Central AsiaDendrochronologia3223023610.1016/j.dendro.2014.03.002Open DOISearch in Google Scholar

Zhang TW, Zhang RB, Yuan YJ, Gao YQ, Wei WS, Diushen M, He Q, Shang HM and Wang J, 2015. Reconstructed precipitation on a centennial timescale from tree rings in the western Tien Shan Mountains, Central Asia. Quaternary International 358: 58–67, 10.1016/j.quaint.2014.10.054.ZhangTWZhangRBYuanYJGaoYQWeiWSDiushenMHeQShangHMWangJ2015Reconstructed precipitation on a centennial timescale from tree rings in the western Tien Shan Mountains, Central AsiaQuaternary International358586710.1016/j.quaint.2014.10.054Open DOISearch in Google Scholar

Recommended articles from Trend MD

Plan your remote conference with Sciendo