Hostname: page-component-586b7cd67f-dsjbd Total loading time: 0 Render date: 2024-11-22T05:48:00.972Z Has data issue: false hasContentIssue false

Wavelet analysis reveals periodic oscillations in a 1700 year ice-core record from Guliya, China

Published online by Cambridge University Press:  14 September 2017

Meixue Yang
Affiliation:
Key Laboratory of Cryosphere and Environment, Cold and Arid Regions Environmental and Engineering Research Institute, Chinese Academy of Sciences, 260 Donggang West Road, Lanzhou 730000, China E-mail: [email protected] Institute of Tibetan Plateau Research, Chinese Academy of Sciences, Beijing 100085, China The Nansen-Zhu International Research Center, Institute of Atmospheric Physics, Chinese Academy of Sciences, Beijing 100101, China
Tandong Yao
Affiliation:
Institute of Tibetan Plateau Research, Chinese Academy of Sciences, Beijing 100085, China
Xiaohua Gou
Affiliation:
Center for Arid Environment and Paleoclimate Research, Lanzhou University, 298 Tianshui Road, Lanzhou 730000, China
Huijun Wang
Affiliation:
The Nansen-Zhu International Research Center, Institute of Atmospheric Physics, Chinese Academy of Sciences, Beijing 100101, China
Thomas Neumann
Affiliation:
Department of Geology, University of Vermont, Burlington, VT 05405-0122, USA
Rights & Permissions [Opens in a new window]

Abstract

Ice cores contribute important records of past climate changes. As one of the thickest ice caps in central Asia, the Guliya ice cap (35°17′ N, 81°29′ E) provides valuable information for this critical region about the past climate and environment changes. We used wavelet analysis to examine periodic temperature and precipitation oscillations over the past 1700 years recorded in the Guliya ice core. The results show non-linear oscillations in the ice-core records, with multiple timescales. Temperature records indicate persistent oscillations with periodicities of approximately 200, 150 and 70 years. Precipitation records show significant periodicities at 200, 100, 150 and 60 years. However, the amplitude modulation and frequency vary with time. Wavelet analysis can explore these time series in greater detail and furnish additional useful information.

Type
Research Article
Copyright
Copyright © The Author(s) [year] 2006

1. Introduction

The study of climate oscillations on decadal to century timescales supports important aspects of climate research. The Climate Variability and Predictability Program (CLIVAR) considers improved understanding of such oscillations as one of its three areas of special interest. According to the International Geosphere–Biosphere Program, climate-change studies can be grouped into three primary timescales: 20, 2 and 0.2 kyr. Of these, 2 kyr is a challenging timescale because it contains decadal- to century-scale climate changes and requires relatively high-resolution proxy data. Although China has abundant historical archives, reliable climate projections back 2000 years are difficult. This is especially true in western China, due to the shortage of historical records. Therefore, proxy data such as those from ice cores and tree rings are important for examining climatic and environmental changes in western China (Reference ThompsonThompson and others, 1989, Reference Thompson1995, Reference Thompson1997, Reference Thompson, Yao, Mosley-Thompson, Davis, Henderson and Lin2000; Reference Yao and ThompsonYao and Thompson, 1992; Reference Yao, Thompson, Jiao, Mosley-Thompson and YangYao and others, 1995, 1996, Reference Yao, Thompson, Shi, Qin and Jiao1997, Reference Yao, Xu and Pu2001, Reference Yao2002; Reference Shi, Yao and YangShi and others, 1999; Reference Yang, Braeuning, Johnson and YafengYang and others, 2002; Reference Gou, Chen, Yang, Li, Peng and JinGou and others, 2005, Reference Gou, Chen, Jacoby and Yang2006).

Power-spectrum and maximum-entropy-spectrum analysis only provide the average periods of the whole time series. They cannot distinguish between oscillations with different phases. However, the wavelet transform can provide more information about the oscillations within the time series (Reference HolschneiderHolschneider, 1995). The amplitude modulation, phase and frequency are important aspects of non-linear oscillations in a complex climate system with multiple timescales. Wavelet analysis is becoming a common tool for analyzing localized power variations within a time series. By transforming a time series into time–frequency space, one is able to determine both the dominant modes of variability and how those modes vary with time. The wavelet transform has been used for numerous studies in geophysics, including tropical convection (Reference Weng and LauWeng and Lau, 1994), the El Niño–Southern Oscillation (ENSO) (Reference Gu and PhilanderGu and Philander, 1995; Reference Wang and WangWang and Wang, 1996; Reference Jevrejeva, Moore and GrinstedJevrejeva and others, 2003), the ENSO– monsoon systems (Reference Torrence and WebsterTorrence and Webster, 1999), atmospheric cold fronts (Reference Gamage and BlumenGamage and Blumen, 1993), the dispersion of ocean waves (Reference Meyers, Kelly and O’BrienMeyers and others, 1993), wave growth and breaking (Reference Liu, Foufoula-Georgiou and KumarLiu, 1995) and coherent structures in turbulent flows (Reference FargeFarge, 1992). See Reference Foufoula-Georgiou and KumarFoufoula-Georgiou and Kumar (1995) for a complete description of geophysical applications, and Reference DaubechiesDaubechies (1992) for a theoretical treatment of wavelet analysis. In our analysis, we use the wavelet method introduced by Reference Torrence and CompoTorrence and Compo (1998) to examine the periodic climate oscillations over the past 1700 years recorded in the Guliya ice core.

2. Guliya Ice Core

Reference Hansen and LebedeffHansen and Lebedeff (1987) drew attention to the potential sensitivity of the Qinghai–Tibetan Plateau, China, with an average altitude of >4500 m, to the anticipated global warming of the next century. Ice cores drilled from several places on the Qinghai–Tibetan Plateau have been used to reconstruct climatic changes in the last few decades (Reference ThompsonThompson and others, 1989, Reference Thompson1995, Reference Thompson1997, Reference Thompson, Yao, Mosley-Thompson, Davis, Henderson and Lin2000; Reference Yao, Xie, Wu and ThompsonYao and others, 1990, Reference Yao, Yang and Lu1994, Reference Yao, Thompson, Jiao, Mosley-Thompson and Yang1995, 1996, Reference Yao, Thompson, Shi, Qin and Jiao1997, Reference Yao, Xu and Pu2001, Reference Yao2002; Reference Yao and ThompsonYao and Thompson, 1992).

Reference ThompsonThompson and others (1995) successfully retrieved three ice cores, with respective lengths of 308.7, 93.2 and 34.5 m, from 6200–6700m a.s.l. on the Guliya ice cap (35˚17 N, 81˚29 E), west Kunlun Shan, Qinghai–Tibetan Plateau (Fig. 1). The longest core reached bedrock at 308.7 m. The Guliya ice cap is the largest in central Asia, with an area of 376 km2 (Reference Yao, Thompson, Shi, Qin and JiaoYao and others, 1997). The 308.7m ice core at 6200ma.s.l. was drilled using a mechanical corer in a dry borehole to 200 m depth, and using a thermal core drill with an alcohol-water mixture from 200 m to bedrock. The ice temperatures in the core were −15.6˚C, −5.9˚C and −2.1 ˚C at 10 and 200 m depth and the bottom, respectively. Visible horizontal layers exist throughout the core, with no apparent breaks. Mass-balance measurements in 1990 and 1991 indicate that the ice cap receives ∼200mmw.e. of accumulation per year (Reference ThompsonThompson and others, 1997; Reference Yao, Thompson, Shi, Qin and JiaoYao and others, 1997).

Fig. 1. The location of the Guliya (G) ice cap and the site where the 308.7 m Guliya ice core was drilled, at 6200 ma.s.l. (after Reference ThompsonThompson and others, 1997).

Reference Yao, Xie, Wu and ThompsonYao and others (1990) dated the Guliya ice cores by counting the distinct annual dust layers and the seasonal variations in δ18O. Reference ThompsonThompson and others (1995) describe the details of the field and laboratory methods used. The temperature record in the Guliya ice core is derived from δ18O measurements. δ18O in precipitation is a reliable indicator of temperature on the northern Tibetan Plateau, as described by Yao and others (1996). First, according to measurements of δ18O and temperature change during precipitation events and the seasonal changes of δ18O and temperature during several years, δ18O in precipitation is positively correlated with temperature, i.e. low temperature in winter corresponds to low δ18O value in precipitation, and high temperature in summer corresponds to high δ18O value in precipitation. Second, the positive correlation between δ18O and temperature can be expressed quantitatively. Whenever δ18O in precipitation increases (or decreases) by 1%, temperature increases (or decreases) by about 1.6˚C (Reference Yao, Thompson, Jiao, Mosley-Thompson and YangYao and others, 1995). In this way, we use δ18O measurements from the Guliya ice core to derive a time series of temperature.

The precipitation time series is derived from three sources: (1) accumulation stakes, (2) visible stratigraphy in snow pits, and (3) insoluble particulates, and both p and tritium horizons in the Guliya ice core (Reference ThompsonThompson and others, 1995). We estimate precipitation on the Guliya ice cap from the net accumulation between annual horizons in the ice core. As studied by Yao and others (1996), net accumulation on this ice cap approximates the actual total precipitation, since sublimation is negligible. The Guliya ice-core records have supplied data for many climatic and environmental studies (Yao and others, 1992, Reference Yao, Yang and Lu1994, Reference Yao, Thompson, Jiao, Mosley-Thompson and Yang1995, 1996; Reference Li, Yao and HuangLi and others, 1993, Reference Li, Yao and Xie1994; Reference ThompsonThompson and others, 1995 Reference Thompson1997; Reference Yang, Yao, He and ThompsonYang and others 2000).

3. Periodic Climate Oscillations Over the Past 1700 Years

Figure 2 shows the variations of δ18O and net accumulation (10 year average) reconstructed from the Guliya ice core for the past 1700 years. Yao and others (1996, Reference Yao, Thompson, Shi, Qin and Jiao1997) examined the climatic and environmental characteristics of the Guliya ice cap from the climate variations on decadal and century timescales since the 3rd century. They concluded that the temperature varies more rapidly (shorter period cycles) than precipitation. They also reported that precipitation changes lagged behind temperature changes. Based on their work, Reference Shi, Yao and YangShi and others (1999) studied the decadal climatic characteristics recorded in the Guliya ice cap and compared the records from the Guliya ice core with historical records from eastern China. This comparison showed similar dating of major climatic events in the ice-core record and the historical documents. Their results suggest that consistencies and similarities in climatic variation extend to the lower and mid-latitudes of central Asia.

Fig. 2. Variations of the 10 year average net accumulation and δ18O in the Guliya ice core over the past 1700 years.

Yao and others (1996, Reference Yao, Thompson, Shi, Qin and Jiao1997) introduced maximum entropy spectrum analysis and highlighted the periodic oscillations recorded in the Guliya ice core over the past 2000 years. However, the maximum entropy spectrum analysis only shows the average periods over the length of the time series. By transforming the time series into time– frequency space, wavelet analysis can also determine the dominant modes of variability and how those modes vary over time.

Figure 3 shows the wavelet power spectrum for the 10 year average temperature time series using the Morlet wavelet for the past 1700 years. The cross-hatched region is the cone of influence, where zero padding has reduced the variance. The black contour lines represent the 10% significance level, using the global wavelet as the background spectrum. The global wavelet variance highlights that there is an approximately 300 year period during the past 1700 years. Around AD600, this period length was approximately 360 years; the period length decreases and the amplitude increases towards the present. At the end of the 20th century, this period length decreased to 300 years and the amplitude was the largest. A 200 year period also exists, but only lasts to AD 1100. After AD 1700, the 200 year period appears again. A period of 100 years also exists only until AD 1100. After examining the variations of temperature and precipitation recorded in the Guliya ice core, Yao and others (1996) pointed out that AD 1100 was a key transition era. Before AD 1100, cold and dry weather prevailed. After AD 1100, temperatures rose and precipitation increased.

Fig. 3. Portrayal of the wavelet analysis for 10 year average temperature over the last 1700 years. The periods shown on the left indicate the wavelet power spectrum. The power has been scaled by the global wavelet spectrum (shown at right). The cross-hatched region is the cone of influence, where zero padding has reduced the variance. The black contour lines represent the 10% significance level, using the global wavelet as the background spectrum. The right plot shows the global wavelet power spectrum.

From AD 1100 to 1350, a 70 year cycle becomes significant (above 10% significance level), while the 200 and 100 year periods disappeared. However, after AD 1400, the 70 year period essentially disappears, while the 150 year period is significant from AD 1400 to 1750. This period was especially strong around AD 1600 (above 10% significance level). Additionally, the wavelet power spectrum (Fig. 3) highlights a significant period of 40years around AD1500 (above 10% significance level).

Figure 4 shows the wavelet power spectrum for 10 year average precipitation for the past 1700 years using the Morlet wavelet. The cross-hatched region is the cone of influence, where zero padding has reduced the variance. The black contour lines represent the 10% significance level, using the global wavelet as the background spectrum. Around AD 400, the precipitation records show periods around 350 years. The period length shortened and the amplitude decreased from that time forward, disappearing around AD 1100. Beginning around AD 1300, periods around 200 years become significant. Periods around 100 years are strongest from AD 1400 to 1600 (above 10% significance level). The amplitude of the 150 year period is also large from AD 300 to 400, AD 900 to 1200, AD1400 to 1600, and AD 1900 to 2000. From AD 700 to 800, the amplitude of the periods around 60 years is large and these periods appear significant (above 1 0% significance level). Short periods between 20 and 30 years also strengthen from AD 300 to 500, and AD 1400 to 2000. Such periods also commonly exist in other proxy data (Reference Taylor, Rose and LamoreyTaylor and others, 1992; Reference Cook, Meko, Stahle, Cleaveland, Dean, Meko and SwetnamCook and others, 1996; Villalba and others, 1996; Reference Baliunas, Frick, Sokoloff and SoonBaliunas and others, 1997; Reference Biondi, Isaacs, Hughes, Cayan and BergerBiondi and others, 2000).

Fig. 4. Portrayal of the wavelet analysis for 10 year average precipitation for the last 1700 years. The periods shown in the left plot indicate the wavelet power spectrum. The power has been scaled by the global wavelet spectrum (shown at right). The cross-hatched region is the cone of influence, where zero padding has reduced the variance. The black contour lines represent the 1 0% significance level, using the global wavelet as the background spectrum. The right plot shows the global wavelet power spectrum.

The wavelet power spectra of the two time series of temperature and precipitation proxy data show that some periodicities such as 200–300 years existed over the entire length of record investigated. Around AD 1500, a 100– 150 year periodicity is significant in both records; but for other periods (e.g. the 60–70 year cycle), although they are significant, they do not appear at the same time.

4. Conclusions

The power-spectrum and maximum-entropy-spectrum analysis only reveal the average periodicities for the whole time series. They do not distinguish the different phases of the oscillations. Therefore, using the power-spectrum or maximum-entropy-spectrum analysis approach is not sufficient to characterize the oscillations within the time series. The wavelet transform can highlight localized intermittent oscillations and provide more information about the oscillations within the time series (Reference HolschneiderHolschneider, 1995). By transforming a time series into time–frequency space, one is able to determine both the dominant modes of variability and how those modes vary in time.

The wavelet power spectrum for temperature (derived from δ18O measurements) using the Morlet wavelet highlighted the intermittent existence of 100 and 200 year periods. From AD 1100 to 1350, periods around 70 years appear very significant, while the 200 and 100 year periods disappear. But a period around 150 years is also evident in the past 1700 years. However, the amplitudes and the period lengths vary with time. Precipitation also shows important periods around 200 and 100 years, as well as notable periods around 150 and 60 years. Again, the amplitudes and periods vary with time.

In this paper, the existence of the above periodic oscillations is highlighted and their variations in amplitude and periodicity are also addressed. However, it is more difficult to determine the external or internal forces that cause such oscillations. What is the relationship of these oscillations to meteorological observations and other climate proxy data? Are they also significant over the same time intervals in other climate proxy data? We will try to answer these questions in a future study.

Acknowledgements

We thank the referees for valuable and constructive criticism. This research was jointly supported by the National Key Basic Research Program (2005CB422004), the Chinese National Science Foundation (40201012, 40571036, 90102005), the Innovation Group Fund (40121101), the Key Projects of the Chinese Academy of Sciences (CAS) (KZCX2-SW-118), the Century Program of CAS (2004401), the Outstanding Youth Funds of the Chinese Natural Science Foundation (40125014), the Chinese Postdoctoral Science Foundation (20040350410) and the Talent Program of the Cold and Arid Regions Environmental and Engineering Research Institute, CAS. Wavelet software was provided by C. Torrence and G. Compo, and is available at http://paos.colorado.edu/research/wavelets/

References

Baliunas, S., Frick, P., Sokoloff, D. and Soon, W.. 1997. Time scales and trends in the Central England Temperature data (1659–1990): a wavelet analysis. Geophys. Res. Lett., 24(11), 13511352.Google Scholar
Biondi, F., Isaacs, C., Hughes, M.K., Cayan, D.R. and Berger, W.H.. 2000. The near-1600 dry/wet knock-out: linking terrestrial and near-shore ecosystems. In Proceedings of the 24th Annual Climate Diagnostics and Prediction Workshop. Springfield, VA, US Department of Commerce, 7679.Google Scholar
Cook, E.R., Meko, D.M., Stahle, D.W. and Cleaveland, M.K.. 1996. Tree-ring reconstructions of past drought across the coterminous United States: tests of a regression method and calibration/ verification results. In Dean, J.S., Meko, D.M. and Swetnam, T.W., eds. Tree rings, environment, and humanity. Tucson, AZ, Radiocarbon, 155169.Google Scholar
Daubechies, I. 1992. Ten lectures on wavelets. Second edition. Philadelphia, PA, Society for Industrial and Applied Mathematics.Google Scholar
Farge, M. 1992. Wavelet transforms and their applications to turbulence. Ann. Rev. Fluid Mech., 24, 395457.Google Scholar
Foufoula-Georgiou, E. and Kumar, P., eds. 1995. Wavelets in geophysics. Beijing, Academic Press.Google Scholar
Gamage, N. and Blumen, W.. 1993. Comparative analysis of low-level cold fronts: wavelet, Fourier, and empirical orthogonal function decompositions. Mon. Weather Rev., 121(10), 28672878.2.0.CO;2>CrossRefGoogle Scholar
Gou, X., Chen, F., Yang, M., Li, J., Peng, J. and Jin, L.. 2005. Climatic response of thick leaf spruce (Picea crassifolia) tree-ring width at different elevations over Qilian Mountains, northwestern China. J. Arid Environ., 61(4), 513524.Google Scholar
Gou, X., Chen, F., Jacoby, G. and Yang, M.. 2006. A comparison of tree-ring records and glacier variations over the past 700 years, northeastern Tibetan Plateau. Ann. Glaciol., 43 (see paper in this volume).Google Scholar
Gu, D. and Philander, S.G.H.. 1995. Secular changes of annual and interannual variability in the Tropics during the past century. J. Climate, 8(4), 864876.2.0.CO;2>CrossRefGoogle Scholar
Hansen, J. and Lebedeff, S.. 1987. Global trends of measured surface air temperature. J. Geophys. Res., 92(D11), 13,34513,372.Google Scholar
Holschneider, M. 1995. Wavelets: an analysis tool. New York, Oxford University Press.Google Scholar
Jevrejeva, S., Moore, J.C. and Grinsted, A.. 2003. Influence of the Arctic Oscillation and El Niño–Southern Oscillation (ENSO) on ice conditions in the Baltic Sea: the wavelet approach. J. Geophys. Res., 108(D21), 4677. (10.1029/2003JD003417.)Google Scholar
Li, Y., Yao, T. and Huang, C.. 1993. The chemical components spatial distribution on the Guliya ice cap. J. Glaciol. Geocryol., 15(3), 467473. [In Chinese with English abstract.]Google Scholar
Li, Z., Yao, T. and Xie, Z.. 1994. The modern atmospheric environmental records in Guliya ice core, Tibetan Plateau. Chinese Sci. Bull., 39(23), 22002201.Google Scholar
Liu, P.C. 1995. Wavelet spectrum analysis and ocean wind waves. In Foufoula-Georgiou, E. and Kumar, P., eds. Wavelets in geophysics. Beijing, Academic Press, 151166.Google Scholar
Meyers, S.D., Kelly, B.G. and O’Brien, J.J.. 1993. An introduction to wavelet analysis in oceanography and meteorology: with application to the dispersion of Yanai waves. Mon. Weather Rev., 121(10), 28582866.Google Scholar
Shi, Y., Yao, T. and Yang, B.. 1999. Decadal climatic variations recorded in Guliya ice core and comparison with the historical documentary data from East China during the last 2000 years. Sci. China D, 42(Supp.), 91100.Google Scholar
Taylor, K., Rose, M. and Lamorey, G.. 1992. Relationship of solar activity and climatic oscillations on the Colorado Plateau. J. Geophys. Res., 97(D14), 15,80315,811.Google Scholar
Thompson, L.G. and 9 others. 1989. Holocene–Late Wisconsin Pleistocene climatic ice core records from Qinghai–Tibetan Plateau. Science, 246(4929), 474–477.CrossRefGoogle ScholarPubMed
Thompson, L.G. and 6 others. 1995. A 1000 year climatic ice-core record from the Guliya ice cap, China: its relationship to global climate variability. Ann. Glaciol., 21, 175181.Google Scholar
Thompson, L.G. and 9 others. 1997. Tropical climate instability: the last glacial cycle from a Qinghai–Tibetan ice core. Science, 276(5320), 18211825.Google Scholar
Thompson, L.G., Yao, T., Mosley-Thompson, E., Davis, M.E., Henderson, K.A. and Lin, P.. 2000. A high-resolution millennial record of the south Asian monsoon from Himalayan ice cores. Science, 289(5486), 19161919.Google Scholar
Torrence, C. and Compo, G.P.. 1998. A practical guide to wavelet analysis. Bull. Am. Meteorol. Soc., 79, 6178.Google Scholar
Torrence, C. and Webster, P.J.. 1999. Interdecadal changes in the ENSO–monsoon system. J. Climate, 12, 26792690.Google Scholar
Villaba, R. and 6 others. 1996. Interdecadal climatic variations in millennial temperature reconstructions from southern South America. In Jones, P.D., Bradley, R.S. and Jouzel, J., eds. Climatic variations and forcing mechanisms of the last 2000 years. Berlin, Springer, 161189. (NATO ASI Series 141.)Google Scholar
Wang, B. and Wang, Y.. 1996. Temporal structure of the Southern Oscillation as revealed by waveform and wavelet analysis. J. Climate, 9(7), 15861598.Google Scholar
Weng, H. and Lau, K.-M.. 1994. Wavelets, period doubling, and time–frequency localization with application to organization of convection over the tropical western Pacific. J. Atmos. Sci., 51(17), 25232541.2.0.CO;2>CrossRefGoogle Scholar
Yang, B., Braeuning, A., Johnson, K.R. and Yafeng, S.. 2002. General characteristics of temperature variation in China during the last two millennia. Geophys. Res. Lett., 29(9), 1324. (10.1029/ 2001GL014485.)CrossRefGoogle Scholar
Yang, M., Yao, T., He, Y. and Thompson, L.G.. 2000. ENSO events recorded in the Guliya ice core. Climatic Change, 47(4), 401409.Google Scholar
Yao, T. and Thompson, L.G.. 1992. Trends and features of climatic changes in the past 5000 years recorded by the Dunde ice core. Ann. Glaciol., 16, 2124.Google Scholar
Yao, T., Xie, Z., Wu, X. and Thompson, L.G.. 1990. Climatic records since the Little Ice Age from the Dunde Ice Cap. Sci. China B, 20(11), 11961201.Google Scholar
Yao, T., Yang, Z. and Lu, J.. 1994. The warming in Tibetan Plateau revealed by ice core records. Chinese Sci. Bull., 39(5), 438441.Google Scholar
Yao, T., Thompson, L.G., Jiao, K., Mosley-Thompson, E. and Yang, Z.. 1995. Recent warming as recorded in the Qinghai–Tibetan cryosphere. Ann. Glaciol., 21, 196200.Google Scholar
Yao, T., Qin, D. and Tian, L.. 1997. The variations of the temperature and precipitation in the past 2 ka in Tibetan Plateau – Guliya ice core records. Sci. China D, 39(4), 425433.Google Scholar
Yao, T., Thompson, L.G., Shi, Y., Qin, D. and Jiao, K.. 1997. Climatic variation since the last interglacial recorded in the Guliya ice core. Sci. China D, 40(6), 662668.Google Scholar
Yao, T., Xu, B. and Pu, J.. 2001. Climatic changes on orbital and suborbital time scale recorded by the Guliya ice core in Tibetan Plateau. Sci. China D, 43(Suppl.), 360369.Google Scholar
Yao, T. and 7 others. 2002. Temperature and methane changes over the past 1000 years recorded in Dasuopu glacier (central Himalaya) ice core. Ann. Glaciol., 35, 379383.Google Scholar
Figure 0

Fig. 1. The location of the Guliya (G) ice cap and the site where the 308.7 m Guliya ice core was drilled, at 6200 ma.s.l. (after Thompson and others, 1997).

Figure 1

Fig. 2. Variations of the 10 year average net accumulation and δ18O in the Guliya ice core over the past 1700 years.

Figure 2

Fig. 3. Portrayal of the wavelet analysis for 10 year average temperature over the last 1700 years. The periods shown on the left indicate the wavelet power spectrum. The power has been scaled by the global wavelet spectrum (shown at right). The cross-hatched region is the cone of influence, where zero padding has reduced the variance. The black contour lines represent the 10% significance level, using the global wavelet as the background spectrum. The right plot shows the global wavelet power spectrum.

Figure 3

Fig. 4. Portrayal of the wavelet analysis for 10 year average precipitation for the last 1700 years. The periods shown in the left plot indicate the wavelet power spectrum. The power has been scaled by the global wavelet spectrum (shown at right). The cross-hatched region is the cone of influence, where zero padding has reduced the variance. The black contour lines represent the 1 0% significance level, using the global wavelet as the background spectrum. The right plot shows the global wavelet power spectrum.