Next Article in Journal
Estimating Real-Time Water Area of Dongting Lake Using Water Level Information
Next Article in Special Issue
Identification of Sulfate Sources and Biogeochemical Processes in an Aquifer Affected by Peatland: Insights from Monitoring the Isotopic Composition of Groundwater Sulfate in Kampinos National Park, Poland
Previous Article in Journal
The Costs of Living with Floods in the Jamuna Floodplain in Bangladesh
Previous Article in Special Issue
Local Meteoric Water Line of Northern Chile (18° S–30° S): An Application of Error-in-Variables Regression to the Oxygen and Hydrogen Stable Isotope Ratio of Precipitation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Stable Isotopes of Precipitation in China: A Consideration of Moisture Sources

1
Key Laboratory of Shale Gas and Geoengineering, Institute of Geology and Geophysics, Chinese Academy of Sciences, Beijing 100029, China
2
Institutions of Earth Science, Chinese Academy of Sciences, Beijing 100029, China
3
University of Chinese Academy of Sciences, Beijing 100049, China
4
College of Water Sciences, Beijing Normal University, Beijing 100875, China
*
Authors to whom correspondence should be addressed.
Water 2019, 11(6), 1239; https://doi.org/10.3390/w11061239
Submission received: 17 April 2019 / Revised: 2 June 2019 / Accepted: 5 June 2019 / Published: 13 June 2019
(This article belongs to the Special Issue Use of Water Stable Isotopes in Hydrological Process)

Abstract

:
An accurate representation of the spatial distribution of stable isotopes in modern precipitation is vital for interpreting hydrological and climatic processes. Considering the dominant impact of moisture sources in controlling water isotopes and deuterium excess, we conducted a meta-analysis of precipitation isotopes using instrumental data from 68 stations around China. The entire country is divided into five regions according to the major moisture sources: Region I (the westerlies domain), Region II (the arctic domain), Region III (the northeast domain), Region IV (the Pacific domain), and Region V (the Tibetan Plateau). Each region has unique features of spatial distribution and seasonal variation for stable precipitation isotopes and deuterium excess. In particular, seasonal variation in Region IV tracks the onset of Asian summer monsoons well. The regional meteoric water lines are presented for each region. A significant temperature effect is found in Regions I and III, with δ18O-temperature gradients of 0.13–0.68‰/°C and 0.13–0.4‰/°C, respectively. However, the reasons for the temperature effects are quite different. In Region I, this effect is caused by the seasonal shift of the westerlies, whereas in Region III, it is caused by the seasonal difference in moisture sources. The precipitation amount effect is most significant in the region along the southeast coast in China, where the δ18O-precipitation amount the gradient is −0.24 to −0.13‰/mm. The findings in our paper could serve as a reference for isotopic application in hydrological and paleo-climatic research.

Graphical Abstract

1. Introduction

The stable oxygen and hydrogen isotope compositions of precipitation are powerful tools for studying hydrological and climatic processes. Spatial variations in stable precipitation isotopes are related to various isotope effects, including the continental effect, the precipitation amount effect, and the altitude effect [1]. Generally, atmospheric moisture is derived primarily from low-latitude oceanic regions, where the initial composition is fixed by evaporation and boundary layer diffusion between the ocean and atmosphere. High δ2H or δ18O values occur near to source regions where vapor is sourced to the atmosphere, and low values occur far from source regions as a result of a progressive rainout history.
Access to data on spatial precipitation isotopes serves as a basis for the widespread application of water isotopes. Since 1961, the International Atomic Energy Agency (IAEA), in cooperation with the World Meteorological Organization (WMO), has been conducting a worldwide survey of the isotopic composition of monthly precipitation [2]. Based on the IAEA and WMO observations, several analyses have been conducted in Southeast Asia, Amazon catchments, Australia, and South Africa [3,4,5,6,7,8]. In China, systemic work has yielded the following findings (including but not limited to): (1) The moisture source is very important to the spatial variations in precipitation isotopes; (2) the temperature effect is significant in North China; and (3) the precipitation amount effect is significant in South China [6,9,10,11].
Despite large-scale precipitation isotope monitoring throughout the country, in West China, observations on precipitation isotopes are quite scarce relative to the huge area and its various terrains. Given the vital role of the Tibetan Plateau and Tianshan Mountains in affecting water resources and the economic issues of the billions of people in Asia, it is necessary to conduct long-term observations of water isotopes in West China. Fortunately, considerably more work has been done recently, and increasingly more isotopic data have been obtained in this region [12,13,14,15,16,17,18]. The addition of new data from West China lays the foundation for a more precise spatial isotopic pattern for China as a whole.
The most common applications of water isotope data are related to the global meteoric water line (GMWL) for hydrological research and the temperature effect for climatic research [1,19]. The use of the GMWL is based on the assumption that precipitation in the studied region has a meteorological pattern similar to the global average. However, this assumption is not valid on a local scale or even on a regional scale, owing to differences in moisture sources and geographical features. The temperature effect may be significantly altered as a result of the addition of another moisture source, especially recycled continental moisture. Therefore, it is necessary to obtain the regional meteoric water line (RMWL) within a region with similar meteorological conditions.
The derivative parameter “deuterium excess” (which was defined as d excess = δ2H–8*δ18O by Dansgaard [1]) is evidently inherited from the initial isotopic composition of an air mass and determined by the air–sea interaction regime [1]; it remains relatively invariant during Rayleigh rainout [20]. D excess in precipitation might help to address the moisture sources and reveal processes, including the condensation of atmospheric vapor, evaporation of water, sublimation of ice, and recycling of continental moisture. At least to the authors’ best knowledge, there has been no spatial analysis of d excess at the regional scale in China, and such an analysis would be undeniably helpful for improving our understanding of precipitation isotopes and hydro-meteorological processes.
In this paper, we divide China into 5 regions based primarily on the moisture sources. Then, based on the new regionalization, we provide spatial and seasonal distributions of water isotopes and d excess. Explanations of the isotopic patterns are related to moisture sources. The RMWLs are presented for each region, which can be used in comparing with the isotopic composition of surface water and groundwater. The relationship between isotopes and temperature (the temperature effect) and that between isotopes and precipitation amount (the precipitation amount effect) are discussed to shed light on their roles in interpreting paleoclimate proxies.

2. Materials and Methods

2.1. Isotopic Data and Analysis

The precipitation isotopic data available for stations in China from the IAEA-WMO Global Network for Isotopes in Precipitation (GNIP) are taken into account in our analysis. To eliminate the edge effect of interpolation data, isotopic data from the GNIP stations around China are adopted. Although this database contains the most comprehensive observation set for water isotopes in precipitation, data in West China are quite scarce relative to the huge area of this region. Recently, new data have been independently generated and published by Tian et al. [12] and Pang et al. [13] in West China, which improved the spatial integrity of precipitation isotopes in China. Although the data are obtained in different years, they are still comparable in analyzing the impact of moisture sources on it and the isotope–meteorological relationship [21,22]. In total, data from a total of 68 stations are used to characterize the isotopic distribution in China (see Table 1). All the δ18O and δ2H values of the precipitation samples are measured by traditional isotope-ratio mass spectrometry or by laser absorption spectroscopy. It should be noted that the data above are from different labs, but all the data were calibrated with respect to Vienna Standard Mean Ocean Water (VSMOW) standards using the same normalization method. Both Tian et al. [12] and Pang et al. [13] tested the same sample in different labs to ensure that the measurement results from different labs were consistent. In total, all the data have a precision within 0.2‰ for δ18O and 1‰ for δ2H.
The multi-year annual average isotopic data (δ18Oave) is calculated using precipitation amount as a weighting factor. Taking δ18O at a station for example,
δ18Oave = sum(δ18Oi, j*Pi, j)/sum(Pi, j)
where δ18Oi, j and Pi, j are the isotopic value and the precipitation amount at month i, year j, respectively.
The regional meteoric water lines (RMWLs) are determined by the least squares regression based on the monthly isotopic data.

2.2. Interpolation

The spatial distribution of precipitation δ18O, δ2H, and d excess is obtained by adopting the Kriging interpolation method. It takes both the distance and degree of variation among known points into consideration and then estimate the unknown points [23]. Among all the Kriging interpolating methods, the ordinary Kriging is selected due to its merits in providing the best linear unbiased predictions and the best sense of minimum variance [24], which helps us to get a regional overview of isotopes. The spherical model, which is one of the most commonly used models in characterizing precipitation isotopes distribution [25], was chosen to describe the experimental semi-variance in this work. Because we are interested in the impact of moisture sources on the precipitation isotopes (large scale atmospheric circulation, i.e., > 100 km), we chose to do the interpolation without the use of ancillary variables, such as elevation [26], which would introduce high amplitude variability into the interpolated surface over short scale lengths.
In order to make the interpolation more reasonable, all the data at the 68 stations were used to generate a spatial distribution of isotopes to the extent of covering all the stations. Then, we added the extent of China to it. All the interpolation and mapping processes were carried out using the ARCGIS software.

2.3. Correlation

The coefficient of determination (r2) between isotopic values, temperature, and precipitation amount is calculated as
r 2 = ( i = 1 n ( x i x ) ( y i y ) ) 2 i = 1 n ( x i x ) 2 i = 1 n ( y i y ) 2
where n is the number of monthly samples, x is the isotopic value, and y is the temperature or precipitation. Linear regression analysis is employed to compute the gradient between isotopic values, temperature, and precipitation. The t test is employed to test the significance of the correlation. Among all the stations in China, more than 90% stations have at least 19 samples (n = 19). According to the t test (α = 0.05), the correlation is significant r2 > 0.3 (n = 19).

3. Results and Discussion

3.1. Regionalization of China

Based on the climate, vegetation, and topography, Luo [27] divided China into 7 regions (Figure 1a). Following this work, Huang [28] re-divided China into 3 regions by merging Regions 1, 2, 3, 4, and part of 5 as the monsoon-affected region, leaving the rest as the Tibetan Plateau (Region 6 and part of 5) and the arid region (Region 7). Generally, these two regionalization schemes are used for the discussion of climatic and geographic features in China. Araguas-Araguas et al. [7] developed a new regionalization scheme based on the moisture sources dominating the precipitation in China (Figure 1b). There are five major air masses that dominate the pluviometric regime of China [7,9]: (1) the polar air mass originating in the Arctic; (2) the westerlies with recycled continental air mass over central Asia; (3) the tropical-maritime air mass originating in the northern Pacific; (4) the equatorial-marine air mass originating in the western equatorial Pacific; and (5) the equatorial-marine air mass originating in the Indian Ocean. Terzer et al. [29] developed a new model, named the regionalized cluster-based water isotope prediction approach (RCWIP), to predict point and large scale spatial-temporal patterns of the stable precipitation isotopes. In their work, they defined 36 climatic clusters, with 4 of them being related to China, whose representative stations are Chiang Rai, Erenhot, Lhasa, and Shanghai, respectively. Their criterion for the clusters were the differences of the climatic variables of temperature, precipitation amount, and vapor pressure between different stations.
The regionalization by Luo [27] incorporates geographical features in China, but it disregards the factor of moisture sources. As a result, one area that has a similar moisture source throughout is divided into two regions. For example, Regions 3 and 4 have the same moisture source from the South Pacific. However, in some cases, regions that have different moisture sources are treated as a single region. For instance, Region 7 is affected by both the westerlies and the arctic air masses. The regionalization by Huang [28] merges all the regions affected by monsoons but overlooks the different geographic features. Although the regionalization by Araguas-Araguas et al. [7] is based on moisture sources, some regions overlap (Figure 1). Furthermore, the extent of some moisture sources is still controversial. Tian et al. [12] found that the northern limit of summer monsoons from the Indian Ocean is in the middle of the Tibetan Plateau, at approximately 34°–35° N, which is farther north than the extent adopted by Araguas-Araguas et al. [7] and Johnson and Ingram [10]. Li et al. [30], Zhou et al. [31], and Xu et al. [32] found that the East Asian summer monsoons (EASM) can reach the Qaidam Basin (farther west than 100° E), whereas previous studies indicated that the direct influence of EASMs reach as far west as 100° E in China [7,33]. Pang [34] summarized the moisture sources over Northwest China on the basis of a meta-analysis of water isotopes in the region, demonstrating that the Arctic polar air mass could reach as far as the Junggar Basin and that summer monsoons from both the Indian and Pacific Oceans has little impact.
When comparing and summarizing from the regionalization schemes of Luo [27], Huang [28], and Araguas-Araguas et al. [7], we find that each region could have its own isotopic patterns affected by moisture sources and other factors. This has also been verified by the recent isotopic studies in China that showed the dominating factor affecting isotopes should be moisture sources and then the local geographical meteorological factors [11,35,36]. To understand the temporal and spatial variation of precipitation isotopes, we devise a new regionalization scheme in this work. The primary factor we consider for this scheme is the moisture source. Thus, we divide the northern continental portion of China into Northwest China (Region I: the westerlies domain) and North China (Region II: the arctic domain). Because the South and North Pacific-dominated regions partially overlap owing to their similar moisture sources, we treat Northeast China as a single region (Region III: the northeast domain) due to its climatological and geographical differences from the rest of the region (Region IV: the Pacific domain). As in the previous regionalization, the Tibetan Plateau is treated as a single region (Region V) because the climatology and tectonic features of the Tibetan Plateau are quite different from those of other regions in China. The Tibetan Plateau plays an important role in affecting regional climatology [7,12,37]. This region impacts both the westerlies and the summer monsoons. The high mountains block the westerlies and split the jet stream, which moves to the south and north of the plateau [7]. The elevated heating by the Tibetan Plateau of the atmosphere plays a fundamental role in the formation and maintenance of the summer circulation, at least over Asia. The onset of the Bay of Bengal monsoon (BOBM) and the EASM can be linked to the thermal and mechanical forces of the Tibetan Plateau [37].
Compared with the regionalization by Araguas-Araguas et al. [7], our regionalization keeps the primary criterion as the moisture sources. However, the improvements can be found in our regionalization by taking into account the geographical and climatic factors. The overlapping regions do not exist as a result of combining the regionalization schemes of Luo [27], Huang [28], and Terzer et al. [29]. The extent of Regions IV and V is ascertained by considering the recent findings of the Indian Ocean and Pacific monsoons mentioned above. The region of the Tibetan Plateau is highlighted in light of its effect on the regional climate.
Collectively, China is divided into 5 regions (Figure 2). All of the results and discussion regarding precipitation isotopes in China will be based on this regionalization.

3.2. Spatial Distribution of Stable Isotopes and D Excess

A clear pattern of the annual average isotopic distribution can be seen in Figure 3a,b. The minimum δ18O with a value of −15.9‰ is at Lhasa station, and the maximum value of −5.7‰ is at Haikou station. The minimum δ2H is −117.2‰ at Lhasa station, and the maximum is −34.7‰ at Guilin station. Isotopes become depleted from south to north overall, which is consistent with the latitude effect reported globally [1]. Though a similar depleted trend can be found from east to west, we could not attribute this trend to the continental or altitude effect because precipitation is formed by different moisture sources from East to West China.
The d excess ranges from 4.7‰ at Chengdu station to 14.8‰ at Guilin station (Figure 3c). Most of the values are close to the global average of 10‰ [1]. Comparing the spatial distribution of d excess with that of δ18O (Figure 2 and Figure 3), we find that in Regions I, III, and V, the contours of δ18O are intensive, whereas the contours of d excess are scarce; in contrast, in Region IV, the δ18O values change slightly, whereas the d excess varies as a complex feature. This pattern occurs because in Regions I, III, and V, δ18O declines gradually along the track of the air masses from the westerlies, the North Pacific, and the Indian Ocean, respectively, whereas in Region IV (the monsoon region), precipitation mainly occurs in summer caused by the monsoons, maintaining consistent δ18O values [38,39]. This pattern fits well with the dominant role of moisture sources in controlling precipitation isotopes.

3.3. Seasonal Variations of Precipitation Isotopes and D Excess

Stations with a full-year record of stable isotopes are selected to show the seasonal variations and their relationship with moisture sources in each region.
(1) Region I
The Wulumuqi and Zhangye stations are used to show the seasonal variations in precipitation isotopes (Figure 4). The δ18O values range from −21.0‰ to −3.7‰, increase gradually from January to July, and then decrease from August to December. The d excess in this region has the largest amplitude (32.0‰) among all the regions. In summer, the d excess is lower than the global average of 10, whereas in winter, it is greater than 10. As noted above, precipitation in this region is mainly formed by moisture from the westerlies. Measuring the isotopic composition of ice core samples taken from a high-altitude glacier in the Tianshan Mountains and the isotopic composition of snow and firn taken from a high mountain of the Siberian Altay, Kreutz et al. [40] and Aizen et al. [41] suggested that seasonal changes in moisture sources and recycling in the Caspian Sea region during transport from the Atlantic to the sampling site in the Tianshan Mountains are responsible for the variability in d excess. Pang et al. [13] and Kong et al. [15] further noted that moisture recycling during summer and autumn amplified the seasonal change in d excess in this region.
(2) Region II
The Altay and Baotou stations are selected to show the seasonal variations in precipitation isotopes (Figure 4). The isotopic change in Region II is very similar to that in Region I. Nevertheless, the d excess is very different—the d excess has a clear seasonal pattern in Region I, which does not exist in Region II. Comparing the d excess at Wulumuqi and Altay stations, we find they are similar in summer (8.4‰ at Wulumuqi and 7.8‰ at Altay) and different in winter (20.3‰ at Wulumuqi and 9.8‰ at Altay). However, the d excess at Baotou station is closer to that of Wulumuqi in winter, which are both larger than 14.0‰. This finding may show that in winter the only moisture source for Altay is the Arctic Ocean, whereas Baotou has an additional moisture source from the westerlies.
(3) Region III
The Qiqihar station is selected to show the seasonal variations in precipitation isotopes (Figure 4). The δ18O values in this region are lower than in other regions. The δ18O values (−9.6‰) are higher in summer and lower (−25.1‰) in winter. In contrast, the d excess values are lower (3.9) in summer and higher (9.3) in winter. The temporal variations of δ18O and d excess values in this region are similar to that in Region I. However, Figure 2 shows that moisture sources in Region III are quite different from those in Region I. In Region III, the moisture source in summer is the Pacific Ocean, which has lower d excess and high δ18O values; however, the source in winter is the Arctic Ocean, which has higher d excess and lower δ18O values [35,39]. This difference in moisture sources leads to the temporal variation of precipitation isotopes and d excess.
(4) Region IV
To track the onset of the monsoon season, Region IV is further divided into the sub-regions of South China (SC), Central China (CC), and North China (NC) (Figure 4).
The Haikou, Hong Kong, Guilin, and Liuzhou stations are selected to show the seasonal variations in precipitation isotopes in Region SC. Figure 4 shows that the δ18O values begin to decrease in April, which is the onset of the monsoons [36,42]. The observed decreasing δ18O values could be explained by the precipitation amount effect. The d excess values change very little throughout the year, exhibiting features of a marine moisture source.
The Nanjing, Fuzhou, Guiyang, and Wuhan stations are selected to show the seasonal variations in precipitation isotopes in Region CC. The δ18O values in CC begin to decrease in May, which is a month later than in Region SC. The d excess in winter (15.0‰) is higher than in summer (8.0‰), which is different from Region SC. This finding implies that precipitation isotopes in winter have been affected by moisture from the north.
The Shijiazhuang, Tianjin, Xian, and Zhengzhou stations are selected to show the seasonal variations in precipitation isotopes in Region NC. Both the δ18O and d excess values in this region vary more significantly than in Regions SC and CC. Obviously, the values are affected by moisture from the north and northwest. However, in summer, the moisture is mainly derived from the Pacific Ocean [38]. The δ18O values begin to decrease in May. Compared with the δ18O values in Regions SC and CC, the decreasing δ18O values reveal the onset of the monsoon season from April to June in different parts of Region IV.
(5) Region V
The Lhasa, Nyalam, and Yushu stations are selected to show the seasonal variations in precipitation isotopes (Figure 4). The seasonal variation in δ18O values in this region is similar to that in Region SC. However, the δ18O values are lower than in Region SC because of the higher altitude in Region V. The d excess values in this region are high in winter and low in summer. It should be noted that the d excess values at Yushu station are different from those at Lhasa and Nyalam stations. In summer, the d excess at Yushu station (14‰) is much larger than that at Lhasa and Nyalam stations (10‰ and 5‰) because the moisture is mainly from the Indian Ocean at Lhasa and Nyalam stations in summer, whereas Yushu station is less influenced by Indian Ocean moisture [12].
In summary, it is found that in each region, the moisture source plays a primary role in controlling the precipitation isotopes. What’s more, with the long-term observation at the GNIP stations (i.e. Hongkong, Wulumuqi), the seasonality remains robust. Compared with the existing regionalization scheme, such a new regionalization scheme helps to identify the factors affecting isotopic variability. However, when using the regionalization scheme, it should be noted that the similarity of isotopic seasonality at each station is incorporated as a factor for the regionalization scheme, but is not the only factor. For instance, the δ18O values are higher in summer and lower in winter at both Wulumuqi and Qiqihar stations, but we divide them into different regions due to different moisture sources and geographical factors. Even in one region, the same moisture source might play different roles at different stations, especially for the stations located at the boundary of two regions, such as Baotou in Region II, Shijiazhuang in Region IV, and Yushu in Region V. For example, both stations of Altay and Baotou are located in Region II of the Artic Region, but in winter the Artic moisture source affects the precipitation isotopes at Altay station more significantly than that at the Baotou station, while in summer the effects at both stations are similar. Therefore, the factors influencing the precipitation isotopes are complex. The regionalization scheme helps to illustrate the impact of moisture sources on the precipitation isotopes, but one has to keep in mind that other factors, including local metrological and geographical conditions, could also modify them.

3.4. Regional Meteoric Water Line (RMWL)

The regional δ2H-δ18O relationship, which is well known as the RMWL, provides a reference for interpreting the provenance of surface water and groundwater. The East Asia meteoric water line (EAMWL) of δ2H = 7.92δ18O + 9.2 was obtained by Araguas-Araguas et al. [7]. This line is indistinguishable from the GMWL of δ2H = 8δ18O + 10. For regions where isotopic data are unavailable, the EAMWL is very useful as a reference for groundwater studies; however, it cannot represent the features of local precipitation for a given region because precipitation isotopes at each region could be modified by local climatic factors, including the origin of the vapor mass, sub-cloud evaporation during rainfall, and the seasonality of precipitation [43]. Thus, in this case, the EAMWL cannot define the groundwater input. Taking these local factors into consideration, we present the meteoric water lines based on the monthly isotopic data for each region (Figure 2) with a dominant single moisture source and similar climatological characteristics (Table 2).
The slopes of RMWLs for Regions I–IV, which have an average of 7.19 with a deviation of 0.52, are slightly less than the slope (7.92) of EAMWL and the slope (8) of GMWL, while the slope of RMWL for Region V (the Tibetan Plateau) is 8.41, which is larger than all the other slopes (Table 2). It is known that the sub-cloud evaporation decreases the slope of the meteoric water line, whereas moisture recycling increases it [13,15,43,44]. Although the seasonality of precipitation could also modify the slope of RMWL, there are no significant differences of precipitation patterns between Region V and the other regions. Thus, we deduce that the Tibetan Plateau region (Region V) is affected by moisture recycling, whereas other regions are affected by sub-cloud evaporation. This conclusion maintains consistency with previous studies in the Tibetan Plateau, Northwest, North, South, and Eastern China [15,25,45,46,47].

3.5. Spatial Extent of the Isotope–Climate Relationship

Given the parallel behavior of δ18O and δ2H, we use the δ18O–temperature and δ18O–precipitation amount relationships to reveal the effect of climate on precipitation isotopes. Figure 5 shows the regional variation in correlation coefficients between δ18O, temperature, and precipitation amount. According to the available data, the correlation is significant when the correlation coefficient is larger than 0.3 (pale yellow region in Figure 5a and dark blue region in Figure 5b); otherwise, it is not significant.
Figure 5a demonstrates that the temperature effect is significant in Regions I, II, and III, whereas in the other regions, it is not significant. Rozanski et al. [21] claimed that the temperature effect was mainly observed in regions of middle and high latitudes, which is consistent with our findings. The gradient of the temperature effect in Region I is identified as 0.13–0.68‰/°C (Figure 6a). The westerlies are the dominant moisture source in Region I [13,15]. In the context of an individual moisture source in the region, a shift in moisture source to the North Atlantic Ocean and the alternation of trajectory in summer and winter are considered as the factors leading to enriched isotopes in summer and depleted isotopes in winter [12,13,40,47]. The temperature effect in Region III is reported as increasing in significance from south to north [39,48]. The moisture in Region III is controlled by the westerlies and polar air masses in winter and by the westerlies and the Pacific Ocean in summer (Figure 1 and Figure 2). In winter, when temperatures are low, isotopes incorporated in both westerly and polar air masses are depleted; in contrast, in summer, when temperatures are high, precipitation isotopes from the westerlies and Pacific Ocean are enriched. Therefore, the temperature effect in Region III is significant. From south to north in Region III, the difference in moisture sources becomes increasingly prominent, which leads to a more significant temperature effect. Thus, we can draw the conclusion that the temperature effect in Region I is caused by the seasonal shift of the westerlies, whereas in Region III, the temperature effect is mainly attributed to the seasonal differences in moisture sources.
The correlation coefficients between stable isotopes and precipitation amount demonstrate the precipitation amount effect (Figure 5b). In total, the correlation between isotopes and precipitation amount is not as significant as that between isotopes and temperature. The most significant parts are located in Southwest and Northeast China, but Figure 6b illustrates that in Southwest China, the slope is negative (anti-correlation), whereas in Northeast China, it is positive. Precipitation events in Northeast China occur mainly in summer, when temperatures are high, whereas in winter, both precipitation and temperature are low. Given the significant temperature effect in Northeast China, it is reasonable to consider that precipitation amount has little impact on precipitation isotopes in this region. A significant precipitation amount effect was observed only in Southwest China and in some regions along the southeast coast. This finding implies that precipitation amount has a stronger effect in monsoon-controlled South China. Johnson and Ingram [10] noted that in the monsoon-affected regions, multi-regression analysis may be particularly useful for capturing the effect of precipitation amount on isotopes. Considering the temporal and spatial variation of monsoon intensity, moisture source analysis is indispensable in addressing monsoon isotopic data.

4. Conclusions

The dominant effect of moisture sources on precipitation isotopes is revealed based on long-term observations at 68 stations in and around China. At a regional scale, China is divided into five moisture source regions: Region I (the westerlies domain), Region II (the arctic domain), Region III (the northeast domain), Region IV (the Pacific domain), and Region V (the Tibetan Plateau domain). In Regions I, III, and V, precipitation isotopes exhibit a clear variation trend and d excess almost remains constant; however, in Region IV, both isotopes and d excess show no clear spatial patterns. The seasonal δ18O and d excess changes in Region IV could be used to track the onset of the EASM.
Regional meteoric water lines (RMWLs) are more useful in comparing with the isotopes of surface water and groundwater. According to the region regionalization above, the RMWLs are δ2H = 7.24δ18O + 1.96, δ2H = 7.55δ18O + 3.91, δ2H = 6.32δ18O − 4.07, δ2H = 7.63δ18O + 8.03, and δ2H = 8.41δ18O + 16.72 in the westerlies domain, the arctic domain, the northeast domain, the Pacific domain, and the Tibetan Plateau domain, respectively.
The temperature effect in Regions I and III is significant for different reasons: in Region I, the δ18O-temperature gradient ranges from 0.13–0.68‰/°C and is larger in Xinjiang, with a value of 0.40-0.68‰/°C, which is caused by the seasonal shift of the westerlies. In Region III, the temperature effect is mainly attributed to the seasonal differences in moisture sources. The precipitation amount effect is most significant in the region along the southeast coast, where the δ18O-precipitation amount gradient is −0.24 to −0.13‰/mm and is affected by the intensity and extent of the East Asia Summer Monsoons (EASM). Each region has its own features of seasonal variation in δ18O and d excess. Our study highlights the importance of moisture sources in affecting the spatial and temporal distribution of precipitation isotopes in China. The results show that an analysis with the new regionalization scheme could capture the features of precipitation isotopes, and thus could help in interpreting the hydrological and climatic processes.

Author Contributions

Y.K. and Z.P. designed the framework. Y.K. wrote the paper. Y.K., K.W., and J.L. plotted the figures. All authors revised and approved the manuscript.

Funding

This study is supported by the National Natural Science Foundation of China (Grant 91647101 and U1703122).

Acknowledgments

This paper is dedicated to the memory of our colleague Prof. Klaus Froehlich from the IAEA, who was involved in many discussions and provided several suggestions for improving this manuscript. The authors would like to acknowledge Prof. Gabriel J. Bowen from the University of Utah, Dr. Tao Pu from the Cold and Arid Regions Environmental and Engineering Research Institute, Chinese Academy of Sciences, the academic editor and four anonymous reviewers for their comments, which greatly helped to increase the quality of the paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Dansgaard, W. Stable isotopes in precipitation. Tellus 1964, 16, 436–468. [Google Scholar] [CrossRef]
  2. IAEA. Isotope Hydrology Information System, the ISOHIS Database. 2006. Available online: http://www.iaea.org/water (accessed on 1 March 2013).
  3. Salati, E.; Dall’Olio, A.; Matsui, E.; Gat, J.R. Recycling of water in the Amazon Basin: An isotopic study. Water Resour. Res. 1979, 15, 1250–1258. [Google Scholar] [CrossRef]
  4. Gonfiantini, R. On the isotopic composition of precipitation in tropical stations. Acta Amaz. 1985, 5, 121–139. [Google Scholar] [CrossRef]
  5. Joseph, A.; Frangi, J.P.; Aranyossy, J.F. Isotope characteristics of meteoric water and groundwater in the Sahelo-Sudanese Zone. J. Geophys. Res. 1992, 97, 7543–7551. [Google Scholar] [CrossRef]
  6. Rozanski, K.; Araguas-Araguas, L. Spatial and temporal variability of stable isotope composition of precipitation over the South American continent. Bulletin De L’Institut Français D’Études Andines 1995, 24, 379–390. [Google Scholar]
  7. Araguas-Araguas, L.; Froehlich, K.; Rozanski, K. Stable isotope composition of precipitation over Southeast Asia. J. Geophys. Res. 1998, 103, 721–742. [Google Scholar] [CrossRef]
  8. Bowen, G.J.; Revenaugh, J. Interpolating the isotopic composition of modern meteoric precipitation. Water Resour. Res. 2003, 39, 1299. [Google Scholar] [CrossRef]
  9. Wang, D.; Wang, K. Isotopes in precipitation in China (1986–1999). Sci. China Ser. E 2001, 44, 48–51. [Google Scholar] [CrossRef]
  10. Johnson, K.; Ingram, B.L. Spatial and temporal variability in the stable isotope systematics of modern precipitation in China: Implications for paleoclimate reconstructions. Earth Planet. Sci. Lett. 2004, 220, 365–377. [Google Scholar] [CrossRef]
  11. Liu, J.; Song, X.; Yuan, G.; Sun, X.; Yang, L. Stable isotopic compositions of precipitation in China. Tellus 2014, 66, 22567. [Google Scholar] [CrossRef]
  12. Tian, L.; Yao, T.; MacClune, K.; White, J.W.C.; Schilla, A.; Vaughn, B.; Vachon, R.; Ichiyanagi, K. Stable isotopic variations in west China: A consideration of moisture sources. J. Geophys. Res. 2007, 112, D10112. [Google Scholar] [CrossRef]
  13. Pang, Z.; Kong, Y.; Froehlich, K.; Huang, T.; Yuan, L.; Li, Z.; Wang, F. Processes affecting isotopes in precipitation of an arid region. Tellus 2011, 63, 352–359. [Google Scholar] [CrossRef] [Green Version]
  14. Kong, Y.; Pang, Z. Evaluating the Sensitivity of Glacier Rivers to Climate Change based on Hydrograph Separation of Discharge. J. Hydrol. 2012, 434–435, 121–129. [Google Scholar] [CrossRef]
  15. Kong, Y.; Pang, Z.; Froehlich, K. Quantifying recycled moisture fraction in precipitation of an arid region using deuterium excess. Tellus 2013, 65, 19251. [Google Scholar] [CrossRef]
  16. Kong, Y.; Pang, Z. Statistical analysis of stream discharge in response to climate change for Urumqi river catchment, Tianshan Mountains, Central Asia. Quat. Int. 2014, 336, 44–51. [Google Scholar] [CrossRef]
  17. Wang, S.; Zhang, M.; Hughes, C.; Zhu, X.; Dong, L.; Ren, Z.; Chen, F. Factors controlling stable isotope composition of precipitation in arid conditions: An observation network in the Tianshan Mountains, central Asia. Tellus B 2016, 68, 26206. [Google Scholar] [CrossRef]
  18. Li, Z.X.; Feng, Q.; Li, Z. Climate Background, Facts, and Hydrological Effects of Multiphase Water Transformation in Cold Regions OF the Western China: A review. Earth-Sci. Rev. 2018. [Google Scholar] [CrossRef]
  19. Craig, H. Isotopic variations in meteoric waters. Science 1961, 133, 1702–1703. [Google Scholar] [CrossRef]
  20. Gat, J. Oxygen and hydrogen isotopes in the hydrological cycle. Annu. Rev. Earth Planet. Sci. 1996, 24, 225–262. [Google Scholar] [CrossRef]
  21. Rozanski, K.; Araguas-Araguas, L.; Gonfiantini, R. Relation between long-term trends of oxygen-18 isotope composition of precipitation and climate. Science 1992, 258, 981–985. [Google Scholar] [CrossRef]
  22. Stumpp, C.; Klaus, J.; Stichler, W. Analysis of long-term stable isotopic composition in German precipitation. J. Hydrol. 2014, 517, 351–361. [Google Scholar] [CrossRef]
  23. Hatvani, I.G.; Leuenberger, M.; Kohan, B.; Kern, Z. Geostatistical analysis and isoscape of ice core derived water stable isotope records in an Antarctic macro region. Polar Sci. 2017, 13, 23–32. [Google Scholar] [CrossRef]
  24. Oliver, M.A.; Webster, R. A tutorial guide to geostatistics: Computing and modelling variograms and kriging. Catena 2014, 113, 56–69. [Google Scholar] [CrossRef]
  25. Liu, J.R.; Song, X.F.; Yuan, G.F.; Sun, X.; Liu, X.; Wang, S. Characteristics of δ18O in precipitation over Eastern Monsoon China and the water vapor sources. Chin. Sci. Bull. 2010, 55, 200–211. [Google Scholar] [CrossRef]
  26. Bowen, G.J.; Wilkinson, B. Spatial distribution of δ18O in meteoric precipitation. Geology 2002, 30, 315–318. [Google Scholar] [CrossRef]
  27. Luo, K. Draft of natural geography regionalization of China. Acta Geogr. Sin. 1954, 20, 379–394. [Google Scholar]
  28. Huang, B. Draft of the complex physical geographical division of China. Chin. Sci. Bull. 1959, 18, 594–602. [Google Scholar]
  29. Terzer, S.; Wassenaar, L.I.; Araguás-Araguás, L.J.; Aggarwal, P.K. Global isoscapes for δ18O and δ2H in precipitation: Improved prediction using regionalized climatic regression models. Hydrol. Earth Syst. Sci. 2013, 17, 4713–4728. [Google Scholar] [CrossRef]
  30. Li, Z.X.; Yao, T.; Tian, L. Variation of δ18O in precipitation in annual timescale with moisture transport in Delingha region. Earth Sci. Front. 2006, 13, 330–334. [Google Scholar]
  31. Zhou, S.Q.; Nakawo, M.; Sakai, A.; Matsuda, Y.; Duan, K.Q.; Pu, J.C. Water isotope variations in the snow pack and summer precipitation at July 1 Glacier, Qilian Mountains in northwest China. Chin. Sci. Bull. 2007, 52, 2963–2972. [Google Scholar] [CrossRef]
  32. Xu, G.; Chen, T.; Liu, X.; An, W.; Wang, W.; Yun, H. Potential linkages between the moisture variability in the northeastern Qaidam Basin, China, since 1800 and the East Asian summer monsoon as reflected by tree ring δ18O. J. Geophys. Res. 2011, 116, D09111. [Google Scholar] [CrossRef]
  33. Winkler, M.G.; Wang, P.K. The Late-Quaternary vegetation climate of China. In Global Climates since the Last Glacial Maximum; University of Minnesota Press: St. Paul, MN, USA, 1993; pp. 221–261. [Google Scholar]
  34. Pang, Z. Mechanism of water cycle changes and implications on water resources regulation in Xinjiang Uygur Autonomous Region. Quat. Sci. 2014, 34, 907–917. [Google Scholar]
  35. Li, J.; Pang, Z.; Kong, Y.; Zhou, M.; Huang, T. Spatial distributions of stable isotopic composition and deuterium excess in precipitation during the summer and winter seasons in China, Fresenius Environmental Bulletin. Fresenius Environ. Bull. 2014, 23, 2074–2085. [Google Scholar]
  36. Kong, Y.; Pang, Z. A positive altitude gradient of isotopes in the precipitation over the Tianshan Mountains: Effects of moisture recycling and sub-cloud evaporation. J. Hydrol. 2016, 542, 222–230. [Google Scholar] [CrossRef]
  37. Wu, G.; Zhang, Y. Tibetan Plateau Forcing and the Timing of the Monsoon Onset over South Asia and the South China Sea. Mon. Weather Rev. 1998, 126, 913–927. [Google Scholar] [CrossRef]
  38. Yamanaka, T.; Shimada, J.; Hamada, Y.; Tanaka, T.; Yang, Y.; Zhang, W.; Hu, C. Hydrogen and oxygen isotopes in precipitation in the northern part of the North China Plain: Climatology and inter-storm variability. Hydrol. Process. 2004, 18, 2211–2222. [Google Scholar] [CrossRef]
  39. Liu, Z.; Bowen, G.J.; Welker, J.M. Atmospheric circulation is reflected in precipitation isotope gradients over the conterminous United States. J. Geophys. Res. 2010, 115, D22120. [Google Scholar] [CrossRef]
  40. Kreutz, K.J.; Wake, C.P.; Aizen, V.B.; Cecil, L.D.; Synal, H.A. Seasonal deuterium excess in a Tien Shan ice core: Influence of moisture transport and recycling in Central Asia. Geophys. Res. Lett. 2003, 30, 1922. [Google Scholar] [CrossRef]
  41. Aizen, V.; Aizen, E.; Fujita, K.; Nikitin, S.; Kreutz, K. Stable-Isotope Time Series and Precipitation Origin from Firn-Core and Snow Samples, Altai Glaciers, Siberia. J. Glaciol. 2005, 51, 637–654. [Google Scholar] [CrossRef]
  42. Ding, Y.; Chan, J. The East Asian summer monsoon: An overview. Meteorol. Atmos. Phys. 2005, 89, 117–142. [Google Scholar]
  43. Clark, I.; Fritz, P. Environmental Isotopes in Hydrogeology; Clark, I., Fritz, P., Eds.; Lewis Publishers: New York, NY, USA, 1997. [Google Scholar]
  44. Froehlich, K.; Kralik, M.; Papesch, W.; Rank, D.; Scheifinger, H.; Stichler, W. Deuterium excess in precipitation of Alpine regions-moisture recycling. Isot. Environ. Health Stud. 2008, 44, 61–70. [Google Scholar] [CrossRef] [PubMed]
  45. Chen, Z.; Nie, Z.; Zhang, G.; Wan, L.; Shen, J. Environmental isotopic study on the recharge and residence time of groundwater in the Heihe River Basin, Northwestern China. Hydrogeol. J. 2006, 14, 1635–1651. [Google Scholar] [CrossRef]
  46. Tian, L.; Yao, T.; Schuster, P.F.; White, J.W.C.; Ichiyanagi, K.; Pendall, E.; Pu, J.; Yu, W. Oxygen-18 concentrations in recent precipitation and ice cores on the Tibetan Plateau. J. Geophys. Res. 2003, 108, 4293. [Google Scholar] [CrossRef]
  47. Kong, Y.; Pang, Z.; Li, J.; Huang, T. Seasonal variations of water isotopes in the Kumalak River catchments, Western Tianshan Mountains. Central Asia. Fresenius Environ. Bull. 2014, 23, 169–174. [Google Scholar]
  48. Li, X.; Zhang, M.; Ma, Q.; Li, Y.; Wang, S.; Wang, B. Characteristics of Stable Isotopes in Precipitation over Northeast China and Its Water Vapor Sources. Chin. J. Environ. Sci. 2012, 33, 2924–2931. [Google Scholar]
Figure 1. Regionalization of China by Luo [27]. (a) Numbers 1, 2, 3, 4, 5, 6, and 7 indicate the Northeast region, the North region, the Central region, the South region, the Sichuan-Yunnan region, the Tibetan region, and the Mongolia-Xinjiang region, as defined by Araguas-Araguas et al. [7]. (b) Numbers 1, 2, 3, 4, and 5 indicate the Arctic region, the Westerly region, the North Pacific region, the South Pacific region, and the Indian Ocean region.
Figure 1. Regionalization of China by Luo [27]. (a) Numbers 1, 2, 3, 4, 5, 6, and 7 indicate the Northeast region, the North region, the Central region, the South region, the Sichuan-Yunnan region, the Tibetan region, and the Mongolia-Xinjiang region, as defined by Araguas-Araguas et al. [7]. (b) Numbers 1, 2, 3, 4, and 5 indicate the Arctic region, the Westerly region, the North Pacific region, the South Pacific region, and the Indian Ocean region.
Water 11 01239 g001
Figure 2. Regionalization of China based on moisture sources and meteorological conditions. The dashed line shows the extent of the sub-regions in Region IV.
Figure 2. Regionalization of China based on moisture sources and meteorological conditions. The dashed line shows the extent of the sub-regions in Region IV.
Water 11 01239 g002
Figure 3. Contours of stable isotopes of precipitation in China: (a) long-term annual average δ18O values; (b) long term annual average δ2H values; (c) d excess calculated based on δ18O and δ2H values.
Figure 3. Contours of stable isotopes of precipitation in China: (a) long-term annual average δ18O values; (b) long term annual average δ2H values; (c) d excess calculated based on δ18O and δ2H values.
Water 11 01239 g003
Figure 4. Seasonal variations in precipitation (pillars), δ18O (solid circles), and d excess (hollow circles) in each region of China.
Figure 4. Seasonal variations in precipitation (pillars), δ18O (solid circles), and d excess (hollow circles) in each region of China.
Water 11 01239 g004
Figure 5. Correlation coefficients between δ18O and temperature (a) and between δ18O and precipitation amount (b). The dotted lines indicate the range of each region.
Figure 5. Correlation coefficients between δ18O and temperature (a) and between δ18O and precipitation amount (b). The dotted lines indicate the range of each region.
Water 11 01239 g005
Figure 6. The gradients between (a) δ18O and temperature (unit: ‰/°C) and (b) between δ18O and precipitation amount (unit: ‰/mm). The dashed lines circle the region where the correlations between δ18O and temperature (a) and between δ18O and precipitation amount (b) are significant (larger than 0.3). The dotted lines indicate the range of each region.
Figure 6. The gradients between (a) δ18O and temperature (unit: ‰/°C) and (b) between δ18O and precipitation amount (unit: ‰/mm). The dashed lines circle the region where the correlations between δ18O and temperature (a) and between δ18O and precipitation amount (b) are significant (larger than 0.3). The dotted lines indicate the range of each region.
Water 11 01239 g006
Table 1. Summary of precipitation amount of weighted annual isotopes (δ18O and δ2H) and d excess with their standard deviations std1, std2, and std3, multi-year average annual precipitation (P), multi-year average annual temperature (T) at the selected stations, the number of months at the GNIP stations, and daily number of samples at stations from Tian et al. [12] and Pang et al. [13] (n).
Table 1. Summary of precipitation amount of weighted annual isotopes (δ18O and δ2H) and d excess with their standard deviations std1, std2, and std3, multi-year average annual precipitation (P), multi-year average annual temperature (T) at the selected stations, the number of months at the GNIP stations, and daily number of samples at stations from Tian et al. [12] and Pang et al. [13] (n).
Sample SiteLongitude (E)Latitude (N)Altitude (m)δ2H (‰)Std1 (‰)δ18O (‰)Std2 (‰)d excess (‰)Std3 (‰)P (mm)T (°C)nObservation PeriodSource
Omsk73.3855.0194−98.8/−13.5/8.9/361.02.781990GNIP
Enisejsk92.1558.4598−98.4/−13.3/7.9/609.00.5121990GNIP
Novosibirsk82.9055.03162−104.3/−14.6/12.8/479.02.5121990GNIP
Irkutsk104.3552.27485−94.513.1−12.20.93.25.7482.02.4141990GNIP
Habarovsk135.1748.5272−102.9/−14.3/11.4/911.02.2111971GNIP
Tashkent69.2741.27428−42.9/−7.0/13.1/316.08.671971GNIP
Kabul69.2134.571860−34.825.6−6.63.518.36.1348.911.8861967–1970, 1973–1974, 1982–1989GNIP
Karachi67.1324.9023−24.315.0−4.12.08.74.4438.824.2391961–1968, 1970, 1973GNIP
New Delihi77.2028.58212−34.618.2−5.42.48.34.0899.124.83041961–1969, 1973–1996, 2000–2008GNIP
Allahabad81.7325.4598−67.8/−9.5/8.1/1531.017.851980GNIP
Shillong91.8825.571598−34.430.2−5.73.811.12.81950.017.0301969–1970, 1973–1976, 1978GNIP
Salagiri79.4418.19259−33.4/−5.0/7.0/640.0/51977GNIP
Colombo79.876.917−23.51.0−4.10.28.90.52179.325.4541983−1987, 1992−1994GNIP
Hongkong114.1622.3166−42.19.8−6.61.210.31.92342.823.04081961−1965, 1973−2007GNIP
Pohang129.3836.036−51.817.4−7.82.510.34.51141.013.31101961−1966, 1973−1976GNIP
Ryori141.8139.03260−55.98.1−8.51.212.04.01400.610.21891979−1986, 1998−2006GNIP
Tokyo139.7735.684−46.57.0−7.20.811.42.81400.215.61831961−1969, 1973−1979GNIP
Yangoon96.1716.7720−29.47.8−4.50.56.33.52434.727.3181961−1963GNIP
Bangkok100.5013.732−42.97.7−6.63.99.56.61638.828.53471968−1970, 1973−2008GNIP
Ko Sichang100.8013.1726−38.24.9−6.10.610.21.71273.328.3401985, 1987−1991GNIP
Ko Samui100.059.477−28.310.4−4.91.310.70.11444.227.7271979−1980, 1982−1983GNIP
Singapore103.901.3532−46.52.5−7.50.213.11.42345.926.2521968−1969, 1973−1975GNIP
Luang-Prabang102.1319.88305−53.210.5−7.81.38.92.21407.926.2261961−1964GNIP
Qiqihar123.9147.38147−80.96.2−10.80.65.81.8570.23.9501988−1992GNIP
Haerbin126.6245.68172−78.04.3−10.50.95.84.7623.79.1321986−1990, 1996−1997GNIP
Hetian79.9337.131375−34.94.8−5.70.611.11.9217.89.2461988−1992GNIP
Wulumuqi87.6243.78918−74.015.7−10.92.112.92.9336.37.31231986−1992, 1996−1998, 2001−2003GNIP
Zhangye100.4338.931483−47.68.5−6.61.35.47.7150.37.8791986−1992, 1995−1996, 2001−2003GNIP
Lanzhou103.8836.051517−41.612.0−6.21.57.96.9345.810.9391985−1987, 1996−1990GNIP
Yinchuan106.2138.481112−43.512.3−6.82.211.26.6308.49.2301988−1992GNIP
Shijiazhuang114.4138.0380−56.87.7−8.10.88.33.6572.413.71291985−2003GNIP
Yantai121.4037.5347−50.64.8−7.30.97.78.8575.714.1441986−1991GNIP
Taiyuan112.5537.78778−62.115.1−8.71.87.21.3444.214.8201986−1988GNIP
Changchun125.2143.90237−63.89.8−9.50.812.511.8448.37.4221999−2001GNIP
Jinzhou121.1041.1366−54.65.2−7.41.14.84.6541.017.0121986−1989GNIP
Tianjin117.1639.103−50.03.6−7.70.511.71.4536.513.8641988−1992, 2000−2001GNIP
Baotou109.8540.671067−55.924.7−7.91.07.43.4300.67.0521986−1993GNIP
Lhasa91.1329.703649−117.235.7−15.94.610.31.8427.08.4411986−1992GNIP
Chengdu104.0230.67506-51.211.1-7.01.74.74.0854.216.3651986−1992, 1996−1999GNIP
Kunming102.6825.011892-69.212.1−10.11.711.25.51005.515.51481986−1991, 2000−2003GNIP
Xian108.9334.30397−50.718.9−7.72.711.23.4585.513.4481985−1993GNIP
Zhengzhou113.6534.72110−54.015.4−7.52.25.84.0642.614.1551985−1992GNIP
Wuhan114.1330.6223−48.019.2−6.72.35.45.01248.717.9491986−1998GNIP
Changqing106.6029.62192−69.6 −10.4 13.7 493.0/51992−1993GNIP
Changsha113.0628.2037−32.20.1−5.70.013.00.11280.717.1571988−1992GNIP
Zunyi106.8827.70844−55.67.0−8.41.011.71.9988.115.4701986−1992GNIP
Guiyang106.7126.581071−52.010.5−8.31.414.51.5970.015.3581988−1992GNIP
Guilin110.0825.07170−34.76.1−6.20.914.81.21523.619.0911983−1990GNIP
Nanjing118.1832.1826−53.06.7−8.20.812.61.01212.014.4581987−1992GNIP
Fuzhou119.2826.0816−43.014.0−6.61.59.83.31471.820.4711985−1992GNIP
Liuzhou109.4024.3597−43.910.8−6.51.07.93.01082.520.9451988−1992GNIP
Guangzhou113.3223.137−42.88.8−6.11.16.02.41830.823.0301986−1989GNIP
Haikou110.3520.0315−37.612.9−5.71.68.34.71923.324.9581988−1991, 1996−2000GNIP
Weathership V164.0031.00/−24.04.0−6.10.624.83.4943.419.1811962−1969GNIP
Taguac Guam IS.144.8313.55110−29.912.3−5.11.510.72.92686.126.01031962−1966, 1973−1977GNIP
Wake Island166.6519.283−9.88.7−2.21.27.62.6904.926.61221962−1969, 1973−1976GNIP
Truk151.857.472−32.510.8−5.41.410.31.93606.227.4721968−1969, 1972−1977GNIP
Yap138.099.4923−34.516.2−5.52.59.64.52913.327.0651968−1969, 1972−1976GNIP
Tarawa172.921.354−37.5/−6.1/11.3/3187.025.9171991GNIP
Manila121.0014.5214−44.212.7−6.81.59.83.81975.327.0461961−1965GNIP
Gaoshan86.8343.13545−53.6/−7.7/8.4/390-4.3592003−2004Pang et al. [13]
Houxia87.1843.282100−55.4/−7.9/7.7/4241.5882003−2004 Pang et al. [13]
Nyalam85.9728.183811−85.9/−12.1/10.6/650/2951998−2001 Tian et al [12]
Lhasa91.0329.723650−129.7/−17.0/5.8/411/1521998−2001Tian et al. [12]
Gaize84.0532.34416−92.0/−12.7/9.6/170/1401998−2001Tian et al. [12]
Shiquanhe80.132.54279−76.3/−11.2/12.9/71/621998−2001Tian et al. [12]
Yushu97.0233.023682−82.0/−12.2/15.8/48233391998−2001Tian et al. [12]
Altay88.0847.73737−100.4/−13.8/9.9/194/2261998−2001Tian et al. [12]
Table 2. Regional meteoric water lines for different regions of the westerlies domain (I), the Artic domain (II), the Northeast China domain (III), the Pacific domain (IV), and the Tibetan Plateau domain (V), respectively.
Table 2. Regional meteoric water lines for different regions of the westerlies domain (I), the Artic domain (II), the Northeast China domain (III), the Pacific domain (IV), and the Tibetan Plateau domain (V), respectively.
RegionsRMWLR2 (n)Level of Significance
Iδ2H = (7.24 ± 0.11)δ18O + 1.96 ± 1.380.95 (264)0.01
IIδ2H = ( 7.55 ± 0.10) δ18O + 3.91 ± 1.220.97 (182)0.01
IIIδ2H = (6.32 ± 0.12) δ18O − 4.07 ± 0.970.87 (460)0.01
IVδ2H = (7.63 ± 0.06) δ18O + 8.03 ± 0.410.93 (1246)0.01
Vδ2H = (8.41 ± 0.20)δ18O + 16.72 ± 2.980.97 (108)0.01

Share and Cite

MDPI and ACS Style

Kong, Y.; Wang, K.; Li, J.; Pang, Z. Stable Isotopes of Precipitation in China: A Consideration of Moisture Sources. Water 2019, 11, 1239. https://doi.org/10.3390/w11061239

AMA Style

Kong Y, Wang K, Li J, Pang Z. Stable Isotopes of Precipitation in China: A Consideration of Moisture Sources. Water. 2019; 11(6):1239. https://doi.org/10.3390/w11061239

Chicago/Turabian Style

Kong, Yanlong, Ke Wang, Jie Li, and Zhonghe Pang. 2019. "Stable Isotopes of Precipitation in China: A Consideration of Moisture Sources" Water 11, no. 6: 1239. https://doi.org/10.3390/w11061239

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop