Introduction
Many researchers have used stable isotopes as a tool for reconstructing paleoclimate variations represented in ice-core records (e.g. Reference Cuffey, Alley, Grootes, Bolzan and AnandakrishnanCuffey and others, 1994; Reference Johnsen, Dahl-Jensen, Dansgaard and GundestrupJohnsen and others, 1995). However, the interpretation of stable-isotope records in ice cores must be based on knowledge of the variations of stable isotopes in precipitation with meteorological conditions (Reference DansgaardDansgaard, 1964; Reference Merlivat and JouzelMerlivat and Jouzel, 1979; Reference Jouzel and MerlivatJouzel and Merlivat, 1984; Reference Rozanski, Araguás-Araguás and GonfiantiniRozanski and others, 1992). The International Atomic Energy Agency (IAEA), in cooperation with the World Meteorological Organization (WMO), set up 550 meteorological stations that have measured the isotopic composition of monthly precipitation since 1961. However, in some regions, such as on the Tibetan Plateau (TP), few stations belong to the IAEA/WMO global network for monitoring isotopes in precipitation. IAEA/WMO set up one station at Lhasa on the TP in 1986, but this operated for only 7 years and was of insufficient duration or scope to study the isotopic composition of precipitation over such an extensive region as the TP. In 1991, the Chinese Academy of Sciences launched a program to establish a new observation network to continually survey isotopic compositions on the TP and adjacent regions. Due to the TP’s unique geographic position and topography, its moisture origins are complex. The different moisture sources interact to affect the spatial and temporal variations of stable isotopes in precipitation. Stable isotopes in precipitation show an apparent precipitation ‘amount effect’ in the south of the TP, related to the effect of monsoon activities (Reference Araguás-Araguás, Froelich and RozanskiAraguás-Araguás and others, 1998; Reference Tian, Yao, Numaguti and SunTian and others, 2001b, Reference Tian, Masson-Delmotte, Stievenard, Yao and Jouzelc, Reference Tian2003). On the northern TP, beyond the extent of the monsoon precipitation, temperature effects control δ18O in precipitation (Reference Zhang, Shi and YaoZhang and others, 1995; Reference Yao, Thompson, Mosley-Thompson, Yang, Zhang and LinYao and others, 1996, Reference Yao, Masson, Jouzel, Stiévenard, Sun and Jiao1999; Reference Tian, Yao, Numaguti and SunTian and others, 2001b, Reference Tian2003).
To observe δ18O in precipitation over a more extensive region on the TP, we established the Yushu station in the east, and the Gaize and Shiquanhe stations in the west of the TP. Based on precipitation sampling at Yushu, this paper presents the relationship between δ18O in precipitation and meteorological conditions, and discusses the variation of δ18O with the effect of moisture transport in this region. In addition, we compare the δ18O variation at Yushu with that at the Gaize and Shiquanhe stations, on the western TP, and at the Tuotuohe and Nagqu stations, in the central part of the TP. These stations all lie within a narrow band of latitude.
Methods and Materials
Study area
Yushu lies in the eastern part of the TP, between the Bayan Har Shan and the Tanggula Shan, with an average altitude of 4500 m. The sources of the Yangtze, Yellow and Lantsang rivers are traced to basins at the foot of these mountains (Fig. 1). This region has an annual average temperature of –8˚C and an average annual rainfall of 464 mm.
Precipitation sampling and measurement
The Yushu meteorological station (33˚01’ N, 97˚01’ E; 3682 m a.s.l.) launched a program to monitor oxygen isotopes in precipitation in November 2000 (data analyzed come from the sampling period November 2000–November 2002). We sampled each precipitation event, and 216 samples were collected. Rainfall samples were collected and immediately sealed in plastic bottles. Snow and other solid precipitation was collected on clean porcelain plates, put into clean plastic bags and sealed. After the samples melted at room temperature, they were processed the same as rainfall. All the samples were stored below freezing in a cold laboratory until analyzed. The Key Laboratory of Cryosphere and Environment, Lanzhou, measured oxygen isotope compositions of all the precipitation samples using a MAT-252 mass spectrometer with a precision of 0.2%. The measured ratio of the stable oxygen isotopes in precipitation samples is expressed as parts per mil of their deviation relative to the Vienna Standard Mean Ocean Water (VSMOW).
Meteorological data
The precipitation-sampling program recorded the duration of each precipitation event, the rainfall amount, wind direction and wind speed (before and after each precipitation event). Other relevant meteorological data, such as dew point and surface pressure, were obtained from the National Meteorological Information Center of China. In addition, we used re-analysis datasets from the US National Centers for Environmental Prediction (NCEP) and the US National Center for Atmospheric Research (NCAR) to calculate the wind field and humidity field for 500 hPa over the TP and adjacent regions, to determine the origin of moisture and to analyze the effect of moisture transport upon δ18O of precipitation in this region. The NCEP/NCAR datasets include monthly-averaged wind direction, monthly-averaged wind speed and monthly-averaged relative humidity re-analysis data, with a spatial precision of 2.5˚ longitude–latitude grids (available at http://www.cdc.noaa.gov/cgi-bin/db_search/SearchMenus.pl).
Results
Figure 2 shows the relationship between δ18O in precipitation and air temperature in recent precipitation events at the Yushu meteorological station. The δ18O values in precipitation and temperature (T) in this area correlate with a regression that can be expressed as: δ18O = 0.17T– 12.48 (r = 0.20). The correlation coefficient is 0.20, within a 0.01 confidence limit.
Figure 3 compares the variation of monthly-averaged δ18O and monthly-averaged temperature at the Yushu station. We find a better correlation between monthly-averaged δ18O and monthly-averaged temperature than that for individual precipitation events (above). Results from a linear regression least-squares fit using monthly-averaged data can be expressed as follows:
Clearly, monthly averages of δ18O in precipitation temperature on the eastern TP produce the best least-squares fit and a better δ18O–Tcorrelation. This shows a temperature effect on δ18O in precipitation in the period of sampling from the seasonal variation in this region.
The variations of the δ18O values of precipitation correlate not only with temperature, but also with dew point and surface pressure. The relationship between δ18O in precipitation and dew point (TDP) in recent precipitation events at Yushu is
Although the correlation coefficient is very small (r = 0.12), there is a positive correlation trend between δ18O in precipitation and dew point. Compared with the relationship for individual precipitation events, the correlation between monthly-averaged δ18O and monthly-averaged dew point is better. The correlation coefficient is 0.61, within a 0.01 confidence limit:
There is a positive correlation between δ18O in precipitation and surface pressure (P) at Yushu:
Moreover, there is a stronger correlation using monthly-averaged surface pressures, with a correlation coefficient of 0.45:
At the same time, we calculated the correlations of the precipitation-weighted monthly means of δ18O with air temperature, dew point and surface pressure, and found all the correlations of weighted means of δ18O and meteorological factors were poorer than those of monthly-averaged δ18O and meteorological factors (equations omitted). The results may be attributed to undue weight assigned to correlation coefficients with regard to the Yushu data. In brief, the variation of δ18O in precipitation at the Yushu station is closely related to the variations of air temperature, dew point and surface pressure, especially for the monthly averages. That is, air temperature, dew point and surface pressure are the main environmental variables correlated with δ18O in precipitation at the Yushu station. Accordingly, we attempt to deal with these data (monthly-averaged) for multiple regression analysis. Multiple regression analysis is a convenient method for viewing stable isotopes in precipitation in more than two dimensions, treating δ18O (Reference Johnson and IngramJohnson and Ingram, 2004). We model δ18O as a function of air temperature, dew point and surface pressure by using a simple multiple regression model:
In this study, we used the SPSS default method, Enter, which is a standard approach in regression models. Based on the output of this model, the standardized coefficients, or beta coefficients, are 0.188, 0.363 and 0.109 for the air-temperature, dew-point and surface-pressure variables, respectively. Hence, Equation (6) can be expressed as:
The value of the multiple correlation coefficient (R) between all of the independent variables (air temperature, dew point and surface pressure) and the dependent variable (δ18O) is 0.62. The F statistic is significant at less than the 5% level. Furthermore, the beta coefficient for dew point is higher than those for air temperature and surface pressure. Accordingly, the δ18O values of precipitation correlate more closely to dew point than to air temperature and surface pressure.
Discussion
In recent years, results of some studies confirmed that the region 32−33˚ N (the Tanggula Shan range, extending from west to east) forms an important climatic divide. These massive mountains block the flow of the southwest monsoon from the south. As a result, the monsoon prevails in the south of the Tanggula Shan, whereas the north is chiefly influenced by the continental air mass (Reference YaoYao, 1991; Reference Tian, Yao and YangTian and others, 1995). Hence, the temperature effect controls δ18O in precipitation in the north, whereas there is an ‘amount effect’ for the δ18O values related to the monsoon precipitation in the south. The Yushu station lies within this climatic divide, at the foot of the northeast Tanggula Shan. Hence, Yushu sits in a transition zone between the regions in the south dominated by the monsoon and those in the north dominated by the westerlies. In this zone, the westerlies and monsoon intermingle, with the monsoon dominating in summer. The rest of the year, the westerlies mainly control the moisture origins at Yushu. Consequently, δ18O in precipitation shows temperature dependence at the annual scale. The monsoon activities in summer weaken the correlation between δ18O and air temperature. Other stations lying in the climatic divide include Tuotuohe (34˚13’N, 92˚26’E; 4533ma.s.l.) and Nagqu (31˚29’N, 92˚04’E; 4508ma.s.l.) on the central TP, and Gaize (32˚18’N, 84˚04’E; 4430 m a.s.l.) and Shiquanhe (32˚30’N, 80˚05’E; 4278ma.s.l.) in the west. All these stations lie in the transition belt between ‘temperature effect’ regions in the north and ‘amount effect’ regions in the south, and the variation patterns of δ18O in precipitation at all these stations are similar. Except for Nagqu, δ18O variations in precipitation relate closely to the variations of temperature on a scale of years (see Table 1), and the δ18O values are low in summer, due to the monsoon precipitation (Reference TianTian and others, 2003; Reference Yu, Yao, Tian, Wang and SunYu and others, 2004). The closer proximity to the monsoon region results in the variations of δ18O in precipitation at Nagqu being more similar to the δ18O variations in the south of the TP (Reference Tian, Yao, Sun, Numaguti and ToshioTian and others, 2001a, Reference Tian2003). However, before the monsoon onset, air temperature remains the dominant factor that controls the variations of δ18O in this area (Yu and others, unpublished information; see Table 2). Table 2 also shows that the correlations of δ18O and air temperature before the monsoon onset are more significant than those on a scale of years at these stations; furthermore, the slopes (dd18O/δT) are higher than those for the entire year. Clearly, the monsoon activities in summer affect the relationships between δ18O in precipitation and air temperature in the transition zone.
Figure 4 shows the seasonal variation of δ18O in precipitation with air temperature, dew point and surface pressure at the Yushu station. The bar graph at the bottom of the figure shows daily rainfall. Note that air-temperature (Fig. 4a) and dew-point (Fig. 4b) trends follow similar tendencies, with high values in summer and low values in winter. Although the details of the curve of surface pressure are complex, there also exists a trend of high values in summer and low values in winter (Fig. 4c).
The comparison of surface-pressure and δ18O variations in precipitation is more complex (Fig. 4d). From November 2000 to early May 2001, before the southwest monsoon onset, the δ18O values increased gradually, consistent with the air-temperature, dew-point and surface-pressure trends. Apparently, temperature effects control δ18O in precipitation in this area, but dew point and surface pressure also cause variations in δ18O. Some of the trends of surface pressure and δ18O in precipitation are synchronous: both surface pressure and δ18O values increased gradually in February, then decreased in March and increased again in April. Before the monsoon onset, the moisture origins are controlled by the westerlies, the frequency of precipitation events is very low and the rainfall of each precipitation event is very small (Fig. 4e).
During the rainy season between late May and August 2001, the monsoon activities were very strong at Yushu, and δ18O values were low (they reach a minimum at the end of July or in early August), while summer air temperatures and dew points remained high. After the monsoon onset, however, the surface pressure decreased gradually, coincident with the δ18O values. Apparently, the monsoon precipitation has less effect on the δ18O–P than on the δ18O–T and δ18O–T DP correlations. The southwest monsoon activities in summer cause the frequency of precipitation events to increase rapidly, and rainfall amounts also increase during the monsoon (Fig. 4e).
With the attenuation of monsoon precipitation, the δ18O values increased gradually for a short time in early September 2001. After October 2001, because of the complete monsoon retreat, the temperature remained the dominant factor controlling the variation of δ18O. With air temperature descending gradually, the δ18O values decreased correspondingly, reaching a minimum in early December 2001, when air temperature was also at a low point in its curve. After the monsoon retreat, the westerlies mainly control the moisture origins in the area, and the content of atmospheric water vapor decreased immediately, resulting in sparse precipitation in this region. Similar to air temperature, dew point and surface pressure also descended gradually, contributing to the lowering of δ18O in precipitation.
In 2002, seasonal variations of δ18O with air temperature, dew point and surface pressure were very similar to those in 2001 (Fig. 4). Before the monsoon onset and after the monsoon retreat, the positive correlations between δ18O in precipitation and air temperature, dew point and surface pressure are all significant, and δ18O in precipitation shows very strong temperature dependence. During the monsoon precipitation, however, the δ18O values in precipitation were low, and they correlate much more poorly with air temperature, dew point and surface pressure. This relates to the moisture being transported from the Indian Ocean into the Yushu region. Similarly, before the monsoon onset and after the monsoon retreat, the moisture origins are controlled by the westerlies, resulting in reduced rainfall at the Yushu station in 2002. About 80% of the annual rainfall occurred during the monsoon, as documented in Figure 4e.
Table 3 lists the correlation coefficients and slopes (dδ18O/δT) in the equations of δ18O in precipitation with air temperature, dew point and surface pressure, before monsoon activities (January–May), during the rainy season (June–August), after monsoon retreat (September–December) and in the whole precipitation sampling period (January–December) for both years at the Yushu station. Before the monsoon onset and after the monsoon retreat, the δ18O–T and δ18O–T DP correlations are more significant than those for the whole sampling period, and the slopes (dδ18O/dT) are also higher than those for the full year. However, during the rainy season, the values of the correlation coefficients and slopes in the equations of δ18O and air temperature and dew point are negative (see Table 3). That is, the monsoon precipitation in summer weakens the correlations between δ18O and air temperature and dew point in this region. Compared with air temperature and dew point, the correlation with surface pressure in the whole precipitation sampling period (R = 0.20) appears similar in the period before the monsoon onset (R = 0.29) to that after the monsoon retreat (R = 0.36), even though during the rainy season the correlation coefficient and slope values in the equations of δ18O and surface pressure are positive (see Table 3). We conclude that the southwest monsoon activities have less effect upon the correlation between δ18O and surface pressure than upon the correlations between δ18O and temperature and dew point. As previously shown (Fig. 4c and d), during the monsoon precipitation in summer, the trend of surface pressure approximates that of δ18O in precipitation at the Yushu station.
In general, the variation of δ18O over the measurement period at the Yushu station was controlled by the temperature effect. However, the southwest monsoon transports oceanic moisture with low δ18O into the Yushu region in summer, especially from June to August, resulting in lower δ18O values than would be predicted from air temperature alone. Accordingly, the monsoon activities weaken the relationships between δ18O and air temperature and dew point. Before the monsoon onset and after the monsoon retreat, the periods when the westerlies dominate the moisture origins in this area, the correlations of δ18O with air temperature and dew point are significant.
Based on NCEP/NCAR re-analysis datasets, we calculated the wind field and humidity field at 500 hPa for January and July in 2001 and 2002, representing winter and summer conditions, respectively. Results (January and July in 2001) are shown in Figure 5; the arrow signals indicate wind direction, and the densely shaded areas indicate humidity. This figure shows that in winter 2001 the moisture in the TP and adjacent regions was controlled by the westerlies. This results in a very strong temperature dependence of δ18O in precipitation at all of the stations mentioned across the TP: Yushu in the east, Tuotuohe and Nagqu in the middle, and Gaize and Shiquanhe in the west. In summer 2001, the monsoon activities transported oceanic moisture into these study regions, and heavy isotopes in oceanic moisture were more depleted along the trajectories, resulting in low δ18O values in precipitation and weakening the correlation between δ18O and air temperature. The results in 2002 are very similar to those in 2001 (figures omitted).
Conclusions
Stable isotopes in precipitation depend on many factors such as air temperature, moisture origins and precipitation processes. Here, we discuss the variation of δ18O with air temperature, dew point and surface pressure, and discuss the effect of the southwest monsoon on the δ18O values in precipitation at the Yushu meteorological station.
At the Yushu station, δ18O in precipitation positively correlates with air temperature, especially for monthly averages. The δ18O values also positively correlate with dew point and surface pressure in this region. In summer, however, similar to the Nagqu, Tuotuohe, Gaize and Shiquanhe stations, all lying within a narrow band of latitude, the southwest monsoon brings high precipitation rates, which deplete the heavy isotopes, lowering δ18O values. Both the Yushu station and other stations are in the transition zone between the regions in the south dominated by the monsoon and the regions in the north dominated by the westerlies. Hence, before the monsoon onset and after the monsoon retreat, the correlations between δ18O in precipitation and air temperature and dew point are stronger than for the entire annual cycle at the Yushu station. Accordingly, the monsoon precipitation weakens the relationships between δ18O in precipitation and air temperature and dew point. Clearly, dew point, surface pressure and the monsoon activities all contribute to controlling the δ18O values in precipitation at Yushu.
Acknowledgements
This work was supported by the Innovation Program of the Chinese Academy of Sciences (KZCX3-SW-339), the Collective Innovation of the National Natural Science Foundation of China (40121101), the Ministry of Science and Technology of the People’s Republic of China (2005CB422004) and the National Natural Science Foundation of China (40571039, 40471025). We are greatly indebted to B. Currie and W. Theakstone for constructive comments during reviewing. Some meteorological data used in this paper are from the National Meteorological Information Center of China. The NCEP and NCAR provided the reanalysis datasets.