Translator Disclaimer
1 November 2007 Dendroclimatic Temperature Record Derived from Tree-Ring Width and Stable Carbon Isotope Chronologies in the Middle Qilian Mountains, China
Author Affiliations +

Using long-lived Qilian juniper (Sabina przewalskii Kom.) in the middle Qilian Mountains, the temperature variations in the last 1000 yr were reconstructed. We find that the annual growth ring width and δ13C series mainly reflect variations in regional temperature. Except in May, warmer temperatures indicate greater growth over the period from December to April, and δ13C values in tree-rings are higher for years with higher temperature. The notable features in the temperature reconstruction are the occurrence of the Little Ice Age from A.D. 1600 to 1880 and the abrupt warming over the end of past millennia. The comparison of our chronology to a Northern Hemispheric temperature proxy shows that our tree-ring data will facilitate intercontinental differentiation of large-scale synoptic climate variability.


Tree-ring chronologies provide a unique environmental archive because of their high resolution and the possibility of constructing precisely dated records covering thousands of years. Width and stable carbon isotope ratios of tree rings have been used to reconstruct past climates by correlating with temperature (Wilson and Grinsted, 1977; Sheu et al., 1996; Aucour et al., 2002), relative humidity (Stuiver and Braziunas, 1987; Loader et al., 1995; Liu, et al., 2003) and human activities (Juknys et al., 2002; Jedrysek et al., 2003).

Over the last two decades, the discovery of long-lived juniper (Wang et al., 1982) has demonstrated significant potential for the paleoclimatic reconstruction from Qilian juniper, a tree species found in China (Zhang and Wu., 1997; Kang et al., 2002, Zhang, et al., 2003). Tree-ring widths of Qilian juniper have proven particularly useful as a proxy for water status, such as moisture index and river run-off (Zhang and Wu., 1997; Kang et al., 2002). So far, however, research has been mainly confined to the width index. In the Qilian Mountains, there is a need for complementary investigations at this key region for climatic reconstructions using tree-ring width and isotope techniques in combination. This paper presents 1000-yr temperature-sensitive tree-ring width and stable carbon isotope chronologies derived from living trees of Qilian juniper (Sabina przewalskii Kom.) at timberline. The results differ from Kang et al. (2002) in middle Qilian Mountains, China. Climate reconstructions in the region will provide insights into the spatiotemporal patterns of past climate changes.

Materials and Methods

Study Area and Tree-ring Samples

The Qilian Mountains (37–39°N, ∼99–103°E) are located at the northeastern edge of the Tibetan Plateau, at the convergence of the Qinghai-Xizang (Tibet) Plateau, the Inner Mongolia–Xinjiang Plateau, and the Loess Plateau (Fig. 1). High temperature extremes range from 28.5 to 32.4°C and the low temperature extremes from −27.8 to −29.0°C at Sidalong meteorological station (2600 m elevation) near the sampling site. The annual precipitation ranges from 401.9 to 632.3 mm, the annual evaporation varies between 1041.2 and 1234.2 mm, and the relative humidity is about 57% (Fu and Che, 1990).

Figure 1.

Location of tree-ring sampling site and meteorological stations used in this study.


Following the International Tree-Ring Data Bank (ITRDB) standards (Grissino-Mayer and Fritts, 1997), Qilian juniper samples were collected at Sidalong Forest Center of Sunan County (38°26′N, 99°56′E) in October 2000. Our sampling site was about 15 km away from and on average 300 m higher elevation than the study site of Kang et al. (2002). The sampled trees are generally growing on steep (42°), south-facing slopes near timberline where there is little disturbance due to human activities. Fifty cores were taken from 25 trees at elevations from 3400 to 3550 m in a relatively open stand. An increment borer was used to enable cross dating and to establish tree-ring width chronology. Samples included four discs for isotope analysis.

Tree-ring Width Chronology

Using basic procedures of tree-ring analysis (Fritts, 1976) and of tree-ring research in arid zone (Shao et al., 2003), the skeleton-plot technique was used for the primary dating, and the cross-dating was further checked with the computer program COFECHA (Holmes, 1983).

In order to retain low-frequency information, we used conservative detrending methods to remove age-related growth trends (Fritts, 1976). Each tree-ring series was fitted with a negative exponential curve, while preserving variations of all frequencies that may be related to climate (Cook and Kairiukstis, 1990). Using the ARSTAN function (Cook and Holmes, 1986), we synthesized a standard tree-ring width chronology using double-weighted average methods for detrending series as shown in Figure 2.

Figure 2.

Ring-width chronology developed for Qilian juniper in middle Qilian Mountains, along with the sample depth plot.


The mean sensitivity of width chronology is 0.18 with mean variance 0.21, and first order autocorrelation coefficient is 0.49, which suggests that there is a 24.0% lag between previous and current years. The S/N ratio of standard chronology is 8.9 and the explained variance of first principal component is 40.84% (Liu, et al., 2005).

Tree-ring Carbon Stable Isotope Chronology

Leavitt and Long (1984) suggested that four trees are adequate to assess δ13C within a site. The rings in this study were very narrow, and even under a binocular lens, it was not feasible to attempt to separate annual rings without risking contamination from the previous and following years in the wood sample. So increment cores from four discs and cores from five other different trees were cut into 3-yr ring groups and the samples in the same period were combined covering the period from A.D. 1060 to 2000. All samples were ground in a mill to 60-mesh size. Then the α-cellulose was extracted using the method modified by Mullane et al. (1988) from Green (1963) with three steps: (1) Removal of the extractives with a 1:1 toluene-ethanol mixture, (2) removal of the lignin with a NaClO2-acetic acid solution, and (3) removal of the hemicellulose with NaOH. The cellulose samples were dried and about 6 to 7 mg were weighed and vacuum-sealed in Pyrex tubes with an excess (about 1.5 g) copper oxide (CuO) wire as oxide and some platinum wire as catalyst. Each of the cellulose samples was then combusted at 550°C for 4 h to produced CO2 (Li et al., 1994). Other studies have used a shorter combustion at higher temperature, but the procedure described here obviates the need for quartz-glass tube, which are more expensive and difficult to prepare. Purified CO2 was analyzed on a Delta-plus mass spectrometer (Finnigan MAT) in the State Key Laboratory of Cryospheric Sciences, Chinese Academy of Sciences. The standard deviation during combustion and analysis is below ±0.1‰. The results were corrected and expressed as deviations from the PeeDee Belemnite (PDB) standard: δ13C = (R(sample)/R(standard) − 1) × 103 per mil, where R = 13C/12C.

Figure 3 illustrates a 940-yr δ13C record from AD 1060 to 2000 with 3-yr resolution. The δ13C series varies by about 3.03‰ in the 940 yr. Linear regression of δ13C with time yields a slope of 0.01‰ yr−1 from 1800 to 2000. The long-term decrease from 1750 to 2000, however, is smaller than the results from Mt. Helan and Huangling (Liu et al., 2004). The first-order correlation coefficient of autoregression is 0.65, which implies that there is a significant 42.2% isotope lag effect between the each 3-yr sample. It is noted that this results not removed the effects of averaging 3 yr to one value.

Figure 3.

Raw averaged 3-yr groups δ13C-records of tree-ring cellulose from Qilian juniper and δ13C in atmospheric CO2 covering the period from 1060 to 2000. Open points show the δ13C data in atmospheric CO2 from ice-core (Francey et al., 1999; Trudinger et al., 1999).


Since the beginning of the 19th century, the atmospheric carbon dioxide concentration has increased and its carbon stable isotope ratio (δ13Catm) has decreased (Stuiver and Braziunas, 1987). Both affect the δ13C ratio of plants because the carbon isotope composition of photosynthetically produced organic matter is determined by the δ13C value of the atmospheric CO2 and the leaf intercellular to atmospheric CO2 concentration (Farquhar et al., 1982). These effects will mask the natural climatic signal to a greater or lesser degree. When considering the relationships between cellulose isotopic composition and meteorological parameters, this long-term effect has to be corrected. In this study, values of δ13Catm were taken from ice-core (Trudinger et al., 1999; Francey et al., 1999), as shown in Figure 3, and a correction can be applied for the past 1000 yr as follows (Farquhar et al., 1982):


where δ13Cplant and δ13Catm are stable carbon stable isotope ratios in tree-rings and in atmospheric CO2, respectively. Then using equation (1), we obtain the isotope shift Δ13C series shown in Figure 4, in which the effects of atmospheric CO2 have been removed, for further climatic analyses.

Figure 4.

Carbon isotope discrimination (Δ13C) series using equation (1) to eliminate the effects of declining atmospheric δ13CO2, along with the number of tree cores used for isotope analysis (right axis).


Tree-ring Width and Stable Carbon Isotope/climate Relationships

Climate data from the Zhangye observation station (38°56′N, 100°26′E, 1483 m elevation, period of record 1951–2000), and the Minle observation station (38°27′N, 100°49′E, 2271 m elevation, period of record 1958–2000) are used for climate response analysis. Before the data were used, tests of homogeneity and randomness against the trend were performed used the Mann-Kendall statistic and double-mass analysis. It should be noted that the record from the Zhangye station has missing values in 1951–1953. The reference stations used in these tests (see Fig. 1 for locations) were Wuwei (38°56′N, 100°2′E, 1532 m elevation, period of record of 1951–2000), Sunan (38°50′N, 99°37′E, 2312 m, period of record of 1957–2000), and Qilian (38°11′N, 100°15′E, 2787 m elevation, period of record of 1957–2000).

The statistical tests show that temperature data among the six stations has no significant inhomogeneity and also precipitation trends are concurrent (Table 1), except data from Wuwei station due to its greater distance from Qilian Mountains. We noted that the distance between Sidalong station and our sampling site is small, but the climatic data of this station were not chosen for response analysis because they are incomplete and the record period is short.

Table 1

Correlation matrix of temperature and precipitation among neighboring stations covering the common period.


We calculated simple linear correlation coefficients between width and stable carbon isotope ratios, and monthly temperature and precipitation of the current and previous biological years of all nearby stations; the results show the closer relationships between the two proxies (width and stable carbon isotope ratios) and climate data from Zhangye (1954–2000) and Minle (1958–2000). So we combined records of the two stations from 1958 to 2000 to represent regional climatic variability.

Relationships between tree-ring width index and stable carbon isotope to climatic variables are identified using the correlation analysis by calculating the simple correlation coefficients (Blasing et al., 1984). Since the climatic conditions in the previous year usually have an effect on tree-ring growth of the current year (Fritts, 1976), the monthly climatic variables over a 12-mo period, from October of the prior growth year to September of the current growth year, are used as predictor variables to determine their significance in affecting the concurrent ring width and carbon stable isotope ratio.

The results of correlation analysis show that radial growth is positively correlated with temperature from previous October to September of the current year, except in May (Fig. 5a), with higher growth indicating warmer temperatures and lower growth cooler conditions. The most important periods of temperature influence on tree growth can be divided into two stages: December of previous year to April, and July to September. In addition, during the December of prior year to January and from March to July, tree growth is correlated positively with precipitation. The relationships between carbon stable isotopes and temperature shows similar relationships to width and temperature: in May the Δ13C is related negatively to temperature, but in other periods, the relationship is positive, namely high temperature corresponds to higher δ13C values. The climate variables and stable carbon isotope relationships are relatively complex at various periods, and in April, May, and June, the effects of temperature and precipitation on Δ13C are reversed, but in period from July to September, there were concurrent effects on δ13C variables. In addition, the simple linear correlation coefficient between ring width and Δ13C spanning the period from 1065 to 2000 at a 3-yr time resolution shows a significant positive trend (r = 0.419, P < 0.001), suggesting that the two series contain similar information. As shown in Figure 5, the relationship of ring width (yearly values) and carbon stable isotope (3-yr mean) to temperature from December of the previous year to April were significant with correlation coefficients of 0.599 (P < 0.001) and 0.577 (P < 0.05), respectively. The results indicate that temperature is a major factor limiting tree growth and carbon metabolism; we interpret the tree-growth rate and stable carbon isotope variations in tree-rings at this tree-line site to be primarily controlled by temperature.

Figure 5.

Monthly correlation coefficients of tree-ring width chronology and stable carbon isotope series to combined instrumental data from Zhangye and Minle meteorological station from 1958 to 2000. (a) Correlation coefficients for ring width (1 yr) to climate data. (b) Correlation coefficients for stable carbon isotope (3 yr) to climate data. Dashed lines show the 95% confidence limits. p: the climate data of previous growth year.


Regional Variability In Reconstructed Temperature

Having established good relationships between ring width, δ13C and December–April temperature, we now reconstruct the past 1000 yr temperature changes with 3-yr resolution. The reconstruction equation is


where Y is average temperature, WD is tree-ring width index, and SCIR is stable carbon isotope ratio value, Δ13C. Figure 6 shows the comparison of the observed mean temperature from previous December to current April with the reconstructed temperature in the calibration period. It can be seen that there is a very good association between them (Fig. 6a). The reconstruction captures 75.6% (adjusted for loss of degrees of freedom, R2adj = 71.2%, F = 17.039, P < 0.0001) of temperature variance over the calibration period 1960–2000 (n = 14) with 3-yr resolution. These statistics suggested that the reconstruction formula was reliable, and it could be used to reconstruct past temperature variations.

Figure 6.

Comparison of observed and reconstructed temperature from previous December (12p) to current April (4). (a) Time series plot. Solid line indicates observed value and dash line indicates reconstructed value. (b) Scatter plot.


The temperature standard deviation anomalies are shown in Figure 7 after reconstruction obtained from equation (2). The reconstruction of temperature in middle Qilian Mountains (Fig. 7a) shows that the temperature was relatively high above the mean during the first century of the record, corresponding to the Medieval Warm Period (Zhang, 1994; Shi and Zhang, 1996; Hughes and Diaz, 1994). From 1200 to 1600, the temperature shows greater amplitude fluctuations. After 1600, the record turns to a state of lower temperatures until the later part of 19th century, namely the Little Ice Age (Yao et al., 1991). This is followed by a period of temperature increase to the present time. The δ18O data (indicating temperature fluctuation) in ice-core drilled from Mt. Qilian (Yao et al., 1991), the nearby long- and high-resolution proxy data currently available in this region, show good agreement with our reconstruction except for the period of around 1690. This agreement between such widely separated sites further suggests a common climatic signature.

Figure 7.

Comparison of (a) temperature variations (standard deviations) around mean obtained from tree-ring width and carbon stable isotope in tree-rings in middle Qilian Mountains, and (b) temperature proxy “North Hemispheric regional curve standardization chronology” from Esper et al. (2002). The dashed lines show the trends over the entire periods of the series.


Climatic change studies (Mann et al., 1999; Esper et al., 2002) have documented the occurrence of the Medieval Warm Period and Little Ice Age in the Northern Hemisphere over the last 1000 yr. Figure 7 compares the detailed variations in temperature series from Esper et al. (2002) and our reconstruction. The two chronologies are very similar over the past ∼1000 yr, although some differences in the timing of peak conditions are evident. Esper et al.'s (2002) chronology shows evidence for above average temperatures during the Medieval Warm Period (900–1300), while our reconstruction indicated a warmer period covering the interval 1065–1150. After the year 1200, the two chronologies lock together remarkably well at multidecadal and centennial time scales. There is strong evidence for inferred below-average temperature over the much of the 1580–1880 interval, which may be regarded as an expression of the Little Ice Age. Temperature fluctuation in our series shows lower values with episodes of 1160s–1330s, 1440s–1500s, and 1580s–1880s. These three episodes correspond well with that of by Esper et al. (2002) (Fig. 7b), suggesting that temperature change was synchronous to some extent over a large scale. Our analysis also indicated that two warmest periods cover the intervals 1060–1150 and 1900–2000, with the peak occurring around the 1100 and 1999. The comparison suggests our reconstructed temperature has not included all of the Medieval Warm Period and, perhaps, not even its warmest period.


Additional paleorecords are essential for evaluation of climate change, including detection of possible recent anthropogenic influences. Our 1000-yr tree-ring width and carbon stable isotope chronologies are currently the longest temperature-sensitive tree-ring record on Mt. Qilian. This record shows that climate in this region has undergone oscillations, sometimes of large amplitude during the last millennia. The greatest fluctuation in the temperature during the past 1000 yr occurred in 1390 to 1600. The reconstructed climate data reveal that the Medieval Warm Period and Little Ice Age were synchronous in China and the Northern Hemisphere, and they help us better understand the longer climate change history over space and time.

Tree-ring width and stable carbon isotope studies for Mt. Qilian along east to west transects are in progress, and we believe that they will provide longer detailed climatic information.


This research was support by the Knowledge Innovation Project of Chinese Academy of Sciences (KZCX2-XB2-04-03), the National Natural Science Foundations of China (40501076), and the Foundation of the State Key Laboratory of Cryosphere and Environment (058143).

References Cited


A. M. Aucour, F. X. Tao, S. M. F. Sheppard, N. W. Huang, and C. Q. Liu . 2002. Climatic and monsoon isotopic signals (δD, δ13C) of northeastern China tree rings. Journal of Geophysical Research 107.D9doi:10.1029/2001JD000464. Google Scholar


T. J. Blasing, A. M. Solomon, and D. N. Duvick . 1984. Response function revisited. Tree-Ring Bulletin 44:1–15. Google Scholar


E. R. Cook and R. L. Holmes . 1986. User manual for ARSTAN Tucson Laboratory of Tree Ring Research, University of Arizona. Google Scholar


E. R. Cook and L. A. Kairiukstis , editors. 1990. Methods of Dendrochronology Dordrecht Kluwer Academic Publishers. 394. Google Scholar


J. Esper, E. R. Cook, and F. H. Schweingruber . 2002. Low-frequency signals in long tree-ring chronologies for reconstructing past temperature variability. Science 295:2250–2253. Google Scholar


G. D. Farquhar, M. H. O'Leary, and J. A. Berry . 1982. On the relationship between carbon isotope discrimination and intercellular carbon dioxide concentration in leaves. Australian Journal of Plant Physiology 9:121–137. Google Scholar


R. J. Francey, C. E. Allison, D. M. Etheridge, C. M. Trudinger, I. G. Enting, M. Leuenberger, R. L. Langenfelds, E. Michel, and L. P. Steele . 1999. A 1000-year high precision record of δ13C in atmospheric CO2. Tellus Series B-Chemical and Physical Meteorology 51:2170–193. Google Scholar


H. C. Fritts 1976. Tree Rings and Climate London Academic Press. 567. Google Scholar


E. H. Fu and K. J. Che . 1990. Study on frost-hydrological effect in Qilianshan (north slope). Journal of Lanzhou University (Special issue) 20:33–39. Google Scholar


J. W. Green 1963. Wood cellulose. 9–21. In R. L. Whistler , editor. Methods of Carbohydrate Chemistry III Academic Press. London. Google Scholar


H. D. Grissino-Mayer and H. C. Fritts . 1997. The international tree-ring data bank: an enhanced global database serving the global scientific community. The Holocene 7:235–238. Google Scholar


R. L. Holmes 1983. Computer assisted quality control in tree ring dating and measurement. Tree-Ring Bulletin 43:69–75. Google Scholar


M. K. Hughes and H. F. Diaz . 1994. Was there a Medieval Warm Period, and if so, where and when?. Climatic Change 26:109–142. Google Scholar


M. O. Jedrysek, M. Krapiec, G. Skrzypek, and A. Kaluzny . 2003. Air-pollution effect and paleotemperature scale versus δ13C records in tree-rings and in a peat core (Southern Poland). Water, Air, and Soil Pollution 145:359–375. Google Scholar


R. Juknys, V. Stravinskiene, and J. Vencloviene . 2002. Tree-ring analysis for the assessment of anthropogenic changes and trends. Environmental Monitoring and Assessment 77:81–97. Google Scholar


X. C. Kang, G. D. Cheng, E. S. Kang, and Q. H. Zhang . 2002. Mountain pass run off reconstruction of Heihe river during past 1000 year using tree-ring. Science in China (Series D) 32:8675–685. (In Chinese.). Google Scholar


S. W. Leavitt and A. Long . 1984. Sampling strategy for stable carbon isotope analysis of tree rings in pine. Nature 311:5982145–147. Google Scholar


Z. H. Li, R. M. Liu, Z. S. An, X. D. Wu, and Y. Liu . 1994. Evidence of increasing concentration of atmospheric CO2 infer from stable carbon isotope in tree-rings since Industrial revolution. Chinese Science Bulletin 39:23485–493. (In Chinese.). Google Scholar


Y. Liu, L. M. Ma, S. W. Leavitt, Q. F. Cai, and W. G. Liu . 2004. A preliminary seasonal precipitation reconstruction from tree-ring stable carbon isotopes at Mt. Helan, China, since AD 1804. Global and Planetary Change 41:229–239. Google Scholar


X. H. Liu, D. H. Qin, X. M. Shao, T. Chen, and J. Ren . 2003. Climatic significance of stable carbon isotope in tree-rings of Abies spectabilis in southeastern Tibet. Chinese Science Bulletin 48:182000–2004. Google Scholar


X. H. Liu, D. H. Qin, X. M. Shao, T. Chen, and J. W. Ren . 2005. Temperature variations recovered from tree-rings in the middle Qilian Mountains over the last millennium. Science in China Series D–Earth Sciences 48:4521–529. Google Scholar


N. J. Loader, V. R. Switsur, and E. M. Field . 1995. High-resolution stable isotope analysis of tree-rings: implications of ‘microdendroclimatology’ for palaeoenvironmental research. Holocene 5:457–460. Google Scholar


M. E. Mann, R. S. Bradley, and M. K. Hughes . 1999. Northern hemisphere temperatures during the past millennium: inferences, uncertainties, and limitations. Geophysical Research Letters 26:759–762. Google Scholar


M. V. Mullane, J. S. Waterhouse, and V. R. Switsur . 1988. On the development of a novel method for the determination of stable oxygen isotope ratios in cellulose. Applied Radiation and Isotopes 39:1029–1035. Google Scholar


X. M. Shao, X. Q. Fang, H. B. Liu, and L. Huang . 2003. Dating the 1000-year-old Qilian Juniper in mountains along the eastern margin of the Qaidam Basin. Acta Geographica Sinica 58:190–100. (In Chinese.). Google Scholar


D. D. Sheu, P. Kou, C. H. Chiu, and M. J. Chen . 1996. Variability of tree-ring δ13C in Taiwan fir: Growth effect and response to May–October temperatures. Geochimica et Cosmochimica Acta 60:1171–177. Google Scholar


Y. F. Shi and P. Y. Zhang . 1996. Climate Change History in China Jinan Shandong Sci-tech Press. 1533. (In Chinese.). Google Scholar


M. Stuiver and T. F. Braziunas . 1987. Tree cellulose 13C/12C isotope ratios and climate change. Nature 328:258–60. Google Scholar


A. M. Trudinger, I. G. Enting, R. J. Francey, D. M. Etheridge, and P. J. Rayner . 1999. Long-term variability in the global carbon cycle inferred from a high-precision CO2 and δ13C ice-core record. Tellus Series B–Chemical and Physical Meteorology 51:233–248. Google Scholar


Y. X. Wang, G. Y. Liu, X. G. Zhang, and C. F. Li . 1982. The relationship between the tree ring of Sabina przewalskii and climate changes of thousand years, glacier activity in Qilian Mountains. Chinese Science Bulletin 27:211316–1319. (In Chinese.). Google Scholar


A. T. Wilson and M. J. Grinsted . 1977. 12C/13C in cellulose and lignin as paleothermometers. Nature 265:133–135. Google Scholar


T. D. Yao, Z. C. Xie, X. L. Wu, and L. G. Thompson . 1991. Climatic change since Little Ice Age recorded by Dunde ice cap. Science in China Series B–Chemistry Life Sciences & Earth Sciences 34:6760–767. Google Scholar


D. E. Zhang 1994. Evidence for the existence for the Medieval Warm Period in China. Climate Change 26:289–298. Google Scholar


Q. B. Zhang, G. Cheng, T. Yao, X. Kang, and A. Huang . 2003. A 2326-year tree-ring record of climate variability on the northeastern Qinghai-Tibetan Plateau. Geophysical Research Letters 30:141739. doi: 10.1029/2003GL017425. Google Scholar


Z. H. Zhang and X. D. Wu . 1997. Climatic reconstruction in Qilian Mountains during past 700 years using tree-ring data. Chinese Science Bulletin 42:8849–851. (In Chinese.). Google Scholar
Xiaohong Liu, Xuemei Shao, Liangju Zhao, Dahe Qin, Tuo Chen, and Jiawen Ren "Dendroclimatic Temperature Record Derived from Tree-Ring Width and Stable Carbon Isotope Chronologies in the Middle Qilian Mountains, China," Arctic, Antarctic, and Alpine Research 39(4), 651-657, (1 November 2007).[LIU-X]2.0.CO;2
Accepted: 1 June 2007; Published: 1 November 2007

Back to Top