Next Article in Journal
Sensitivity Analysis of Extreme Daily Rainfall Depth in Summer Season on Surface Air Temperature and Dew-Point Temperature
Previous Article in Journal
RZ-TRADEOFF: A New Model to Estimate Riparian Water and Air Quality Functions
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Water Stable Isotopes in an Alpine Setting of the Northeastern Tibetan Plateau

1
College of Geography and Environmental Science, Northwest Normal University, Lanzhou 730070, China
2
Laboratory of Atmospheric Physics, University of Patras, GR 265 00 Patras, Greece
3
State Key Laboratory of Cryospheric Sciences, Northwest Institute of Eco-Environment and Resources, Chinese Academy of Sciences, Lanzhou 730000, China
*
Author to whom correspondence should be addressed.
Water 2019, 11(4), 770; https://doi.org/10.3390/w11040770
Submission received: 13 March 2019 / Revised: 9 April 2019 / Accepted: 10 April 2019 / Published: 13 April 2019
(This article belongs to the Section Hydrology)

Abstract

:
Hydrological processes produce effects on water resources in inland mountainous regions. To perform a comprehensive investigation of the important segments of the water cycle, using the Qilian Mountains as a case study, precipitation, soil, plant, river, and groundwater were collected during the plant growing season of 2016. All samples were collected on a monthly basis, except precipitation, which was collected on a per event basis. The results showed that: the “temperature effect” was apparent, which suggested a drier climate background; there were differences in the slope and intercept of the local meteoric water line, using different regression methods; and the δ18O of soil water varied greatly in the topsoil, tended to be similar in the deep soil, and became increasingly depleted as the soil depth increased. The responses of the soil water isotopes to precipitation pulses had different boundaries. The major water source for Caragana Fabr. in no-precipitation month was located in the 0–30 cm soil layer, but was different in months when precipitation occurred. Overall, the findings from the stable isotopes provide insights into hydrological processes and offer a platform to understand mountainous water cycle in arid areas.

1. Introduction

As componets of the water molecule, the stable isotopes 18O and 2H (D) are useful tracers to describe regional water cycle because they reveal information about hydrological processes including precipitation, infiltration, evaporation, and transpiration [1]. In spite of their low abundance in nature, these heavy isotopes are sensitive to changes in climate and hydrology. The stable isotopes of different water bodies can be used to investgate atmospheric moisture transport [2,3], plant water sources and water use patterns [4,5,6], recharge mechanisms and transport of soil water [7,8], runoff generation and confluence [9,10], and recharge source and groundwater mechanisms [11,12], etc. Thus, stable isotopes are popular in hydrological and climatic study.
The Qilian Mountains are located in the northeastern edge of the Tibetan Plateau. Lofty mountains and deep valleys occur in a staggered formation, which creates a special geographic position. The mountains are the headwaters of the Shiyang, Heihe, Shule, and many other inland rivers. Therefore, the Qilian Mountains may be considered a “wet island” in the western Gobi desert, and are the water source of the Hexi Corridor oasis. They also play a crucial role in the ecological balance of the Badain Jaran and Tengger deserts. Due to the development of stable isotope research techniques in recent decades, stable isotopes in water bodies have been studied for different areas of the Qilian Mountains [5,13,14,15,16,17]. However, a comprehensive investigation of the different water bodies, from precipitation to surface, soil, and plant waters (sometimes known as soil-plant-atmosphere continuum, SPAC), is still limited, since most isotope sampling studies have only been conducted once or over a very short period of time. In this study three sites were selected: one located on a north-facing slope, the second on a south-facing slope, and the third on the mountain ridge. Soils, plants, river water, and groundwater were sampled during the growing season (i.e., from May to September), while precipitation was sampled throughout the complete year. δ18O and δ2H were measured in order to trace a series of key parameters of the water cycle process in this arid mountainous area, extract baseline information, and study the variations in precipitation and other water isotope characteristics.
The objective of this work is to analyze: (1) the applicability of different local meteoric water line (LMWL) regressions compared to the most commonly used method (ordinary least squares regression, OLSR); (2) the variation in isotopic profiles for soil water at different sites and its influencing factors; (3) how plants use the soil water from different depths in an arid mountainous area, and whether soil water use is the same for one type of plant grown at different sites; (4) the water recharge connection between precipitation, river water, and groundwater; and (5) the relationship of the above-mentioned water bodies from an isotopic perspective. From this information, we can describe the characteristics and interactions of stable isotopes of different water bodies, and provide a reference for the research of water cycle processes in similar mountainous regions outside China.

2. Materials and Methods

2.1. Study Area

The Qilian Mountains (93°30′–103° E, 36°30′–39°30′ N) are located in the northeastern edge of the Tibetan Plateau. They are a geosyncline originating in the Paleozoic era that later formed a fold belt due to Caledonian and Variscan orogeny. They consist of several parallel mountains and valleys, about 850 km long and 200–300 km wide, resulting from the movement of fault blocks beginning in the Cretaceous period. The altitude of the ridges ranges between 4000 and 6000 m. Mountains higher than 4000 m are covered perennially with snow and develop glaciers. Many rivers originate from these mountains and, hence, they are known as the solid water reservoirs of the high mountains. The climate of the Qilian Mountains is either temperate or semi-arid, since it is affected by the continental and the plateau climates. In the cold climate region, winter is long and summer is short, where the average annual temperature is below 4 °C (below −11 °C in January and below 15 °C in July). The annual mean precipitation amounts decrease from 700 mm at the south slope to 400 mm at the north slope. The Qilian Mountains are the headwaters of many rivers of the Hexi Corridor, hence, they are extremely important for this region. As a result of the complex natural conditions and large differences in water-heat conditions of the region, plenty of vegetation species can be found here, and many soil types have been formed, with obvious vertical gradients and horizontal differences. The main soil types found here are alpine cold desert soil, alpine meadow soil, mountain meadow soil, gray cinnamon soil, chernozem, chestnut soil, etc. The vertical vegetation belt, starting at a height of 2000 m, consists mainly of desert steppe belt, steppe belt, forest steppe belt, shrub steppe belt, meadow steppe belt, snow belt, etc. During our spatially-intensive sampling, water samples for isotopic analysis were collected at three locations covering the north slope (Gulang), the south slope (Tianzhu), and the mountain ridge (Wushaoling) (Figure 1 and Table 1). During our temporally-intensive sampling, the one-year precipitation samples were collected on an event basis. All other samples were collected monthly throughout the growing season.

2.2. Sampling and Data

Precipitation samples were collected by the meteorological stations’ personnel in 2016. Precipitation was collected in a clean, large caliber container, positioned on flat ground before the start of each precipitation event, and stored in either plastic bags and plastic bottles. Solid precipitation samples were sealed in plastic bags, left to melt at room temperature, and transferred into bottles and sealed. Liquid precipitation samples were directly stored in bottles and sealed. At the same time, relevant meteorological parameters (air temperature, relative humidity, precipitation amount, etc.) were recorded.
During the growing season, from May to September 2016, we collected soil, plant, river water, and groundwater samples at the three sites in the Qilian Mountains. Samples were collected on a monthly basis. Soil samples were collected by digging a soil pit (1 m × 1 m × 1 m), and collecting samples at every 10 cm vertical interval. The topsoil sampling began at 2 cm below the surface, in order to avoid samples biased from continuous exposure to the free atmosphere.
A representative plant (Caragana Fabr.) was selected for plant xylem sampling. We collected branches of about 0.3–0.5 cm in diameter and 3–5 cm long, put the xylem into sample bottles after removing the bark and phloem, and sealed the bottles.
When collecting river water and groundwater (well or spring water), we first rinsed the sample bottle 3 times with the river water or groundwater, collected water from a few centimeters below the surface (in order to ensure that the water was fully mixed and to avoid evaporation of the surface water, causing isotope fractionation) into the sample bottle, and sealed it.
In total, 896 samples were collected (excluding repeat samples, Table 2). All samples were numbered and labeled according to its type when collected; the date and meteorological and environmental parameters were also recorded. To avoid isotope fractionation caused by evaporation, all samples were sealed immediately using insulation tape. Soil and plant samples were frozen, and the water samples were refrigerated.
Soil and plant water were extracted using the automatic water extraction system LI-2100 (LICA United Technology Limited, Beijing, China). The δ18O and δ2H of soil water, stem water, precipitation, river water, and groundwater were measured using the liquid-water isotope analyzer DLT-100 (Los Gatos Research, CA, USA) at the Stable Isotope Laboratory, College of geography and environmental science, Northwest Normal University. The ratio of 18O/16O and 2H/H isotopes in water samples were measured accurately and expressed as δ18O and δ2H, relative to the Vienna Standard Mean Ocean Water (V-SMOW) using,
δ 18 O = [ R sample R standard 1 ] × 1000
δ 2 H = [ R sample R standard 1 ] × 1000
where Rsample is the ratio of 18O/16O (or 2H /H) in the water sample, and Rstandard is the ratio of 18O/16O (2H/H) in the V-SMOW. The measurement accuracy for δ18O and δ2H is 0.3‰ and 1‰, respectively.
Each sample and standard sample was injected and sequentially analyzed 6 times using a microliter syringe. The results of the first two analyses were discarded since these were performed in order to eliminate any residual water and instrument memory. The final result for each sample was the mean of the measurements based on the last four injections. Quality control was performed via the LWIA Post Analysis (version 3.1.0.9) software, developed by Los Gatos Research (LGR), Inc. (San Jose, CA, USA). The isotope standards used were provided by Los Gatos Research, Inc. According to the isotope ratio range of our samples, three standards were selected from the five LGR standards (expressed relative to the V-SMOV), which include for δ2H and δ18O, respectively: lgr1a: −154.30‰ and −19.50‰, lgr2a: −123.60‰ and −16.14‰, lgr3a: −96.40‰ and −13.10‰, lgr4a: −51.00‰ and −7.69‰, and lgr5a: −9.50‰ and −2.80‰. The complete protocol for the analysis was provided by Los Gatos Research, Inc.
Soil water content (SWC) was calculated by measuring each sample’s wet weight immediately after unsealing it, drying the samples at 105 °C until the weight stabilized, and then measuring the soil’s dry weight.

2.3. LWML Regressions

The relationship between δ2H and δ18O in precipitation at a certain site, known as the local meteoric water line (LMWL), is commonly determined using the ordinary least squares regression (OLSR), which gives equal weighting to each sample [19]. Alternatively, reduced major axis (RMA) regression and major axis (MA) regression can be used to determine LMWL without precipitation weighting [19,20]. More recently, a new precipitation-weighted least squares regression method (PWLSR) was proposed by Hughes and Crawford [21], with corresponding amount-weighted regressions defined by Crawford et al. [20]: precipitation-weighted reduced major axis (PWRMA) regression and precipitation-weighted major axis (PWMA) regression. All these six regressions can be calculated using the Local Meteoric Water Line Freeware software [20].

2.4. IsoSource Model

IsoSource is a procedure developed by Phillips, to calculate the range of all possible source contributions for a mixture when the number of potential sources is greater than n + 1. It has to satisfy the isotope quality conservation [22]:
{ δ 2 H = f 1 δ 2 H 1 + f 2 δ 2 H 2 + f 3 δ 2 H 3 + δ 18 O = f 1 δ 18 O 1 + f 2 δ 18 O 2 + f 3 δ 18 O 3 + f 1 + f 2 + f 3 + = 1
where δ2H and δ18O are the ratios of the stable hydrogen and oxygen isotope in stem water, δ2H1 (δ18O1), δ2H2 (δ18O2), δ2H3 (δ18O3), δ2Hn (δ18On) are the ratios of the stable hydrogen and oxygen isotope in n different water sources, and f1, f2, f3, ……, fn are the proportions of n water sources in the total amount of stem water.

3. Results

3.1. Stable Isotopes in Different Water Bodies

3.1.1. Precipitation

The basic information about δ18O and δ2H in precipitation at the three sites are presented (Table 3). The differences in the mean δ18O, δ2H, and D-excess (i.e., δ2H–8δ18O) [23] among these sites were small, with the greatest δ18O and δ2H values found in Tianzhu, suggesting a more intense under-cloud evaporation during precipitation. The mean D-excess values at the three sites were all higher than that of the Global Meteoric Water Line (GMWL, i.e., 10). This finding was in accordance with the meteorological conditions of the area, which are characterized by low relative humidity and strong evaporation under kinetic fractionation [14]. The correlation and regression relationship between δ18O and precipitation, and between δ18O and temperature are showed in Table 4. The “amount effect” was not significant in the Qilian Mountains, but the “temperature effect” was apparent, which suggested evaporation enrichment.
Craig [24], using 400 precipitation samples from around the globe, found that there is linear relationship between δ2H and δ18O, called the Global Meteoric Water Line (GMWL), where δ2H = 8δ18O + 10. The line calculated using samples collected at a specific region, called the Local Meteorological Water Line (LMWL), often deviates from the GMWL due to different climatic and geographical factors (water vapor sources, air humidity, water vapor path, temperature, altitude, etc.).
The LMWL of the three sampling sites, using six regressions (OLSR, MA, RMA, PWLSR, PWMA, PWRMA), are presented in Figure 2. The slope differences calculated among the various regression methods were the least for Wushaoling, followed by Tianzhu, with the largest differences in slope at Gulang. Crawford et al. [20] suggested that the differences in slope, when calculated by different regressions, were smaller for sites where a positive correlation existed between the precipitation amount and δ18O. However, the conditions in the Qilian Mountains were not like this: the precipitation “amount effect” was not significant here, but the “temperature effect” was apparent (Table 4). The correlation coefficient between δ18O and temperature was the smallest in Gulang, corresponding to the largest slope differences among various regressions. For unweighted regressions, the slope and intercept produced by OLSR were the smallest, and the greatest was produced by MA for all sites, so did the precipitation-weighted methods (Table 5).

3.1.2. Soil Water

Figure 3 presents the isotopic profile and soil water content (SWC) as a function of depth (0–100 cm) during the growing season (May to September). δ18O in soil water was mainly influenced by mixing with precipitation and evaporation [25]. Evaporation from the topsoil during dry intervals between rain events, resulted in heavy isotope enrichment of the residual water. δ18O in soil water changed significantly in the topsoil and tended to remain constant in the deep soil. From the vertical profile of δ18O in soil water (Figure 3), it can be seen that isotope values became increasingly depleted as soil depth increased. This is explained by the fact that topsoil is directly influenced by, and responds quickly to, precipitation, thus modifying the δ18O values. The isotope values in deep soil, which mix with topsoil water during the infiltration of precipitation, represent a hysteresis response. The intra-seasonal variation in isotope values was more pronounced at Gulang and Tianzhu, rather than at Wushaoling, depending on the month and location, likely reflecting the variation in the isotopic signatures of precipitation, as indicated by topsoil water content. These variations were reduced as δ18O became increasingly depleted with soil depth.
The δ18O values of soil water in Gulang and Tianzhu were higher during the summer (June–August) when compared to the rest of the year. This was due to the higher temperature that increased evaporation of soil moisture after precipitation, thus contributing to isotope enrichment. It should be noted that the δ18O in soil water in Wushaoling varied significantly in the 0–70 cm depth, which indicated that precipitation infiltrated the soil up to this depth. Wushaoling is located in a mountain area at the boundary between the monsoon and the non-monsoon zones (between the arid and semi-arid areas), and consequently has an abundant rainfall.
The monthly mean soil water isotope values, weighted by SWC, were calculated, the measured near-surface temperature values were recorded, and any precipitation events one week before sampling were also marked (Figure 4). The changes in monthly δ18O and δ2H were minor, but soil temperature obviously changed. These changes showed that temperature was not the sole factor that impacted soil water isotopes, but that some other factors, like vegetation coverage, also played a part.
The correlation coefficients between δ18O and soil water content are provided in Table 6. Gulang and Tianzhu generally presented positive in May and negative correlation coefficients from June to September. This corresponded to the seasonal changes in precipitation amount. Frequent summer precipitation led to isotopic depleted rain infiltration into deeper soil. However, the infiltration process would decrease the extent of, or slow down the time of, evaporation. However, precipitation in May was less frequent and its amount was lower, therefore soil water could not infiltrate to deep soil layers because of intense evaporation. Consequently, with the increase in SWC, stable isotopes were enriched. On the contrary, the correlation coefficients for Wushaoling were positive in rainy July and August, and negative during the rest of the year. However, in months with frequent precipitation, the relative high temperature, single meadow vegetation, and loose soil texture still allowed for rapid evaporation.
The sites that had multiple precipitation events in the one week period before sampling are used to discuss the response of soil water isotopes to precipitation pulses of different amounts (Wushaoling and Tianzhu) (Figure 5). The response of each layer’s SWC to precipitation was positive: the SWC value increased after precipitation, and the higher the precipitation amount, the bigger the SWC value. While the response was different for the δ18O and δ2H of soil water, apparent boundaries existed in the different soil layers’ response to precipitation. For Wushaoling, the variation could be classified into three boundaries: the 0–20 cm layer varied with SWC, the 30–60 cm layer had a reverse response, and the 70–100cm layer had few changes. For Wushaoling, the variations in both δ18O and δ2H were classified into two boundaries: the 0–40 cm layer varied with precipitation amount, whereas the 50–100 cm layer changes were small. Hence, different sites had different responses to precipitation.

3.1.3. Groundwater and River Water

Rivers in the Qilian Mountains are mainly derived from ice and snow melt, precipitation, and groundwater recharge, so their stable isotopes can reflect the characteristics of water sources and evaporation. The spatial and temporal differences in stable isotopes of river water are influenced by recharge source, season, stream segment, and evaporation. For groundwater, besides its recharge source, the recharge mechanism and the water table can also cause differences in stable isotopes.
In terms of each site, the variation trend in δ18O and δ2H for river water and groundwater were generally consistent, and indicated a similar water source for the two water bodies (Figure 6). In most months, the isotope values of river water were higher than groundwater, which showed that river water experienced more evaporation due to its exposure to air. There was an intersection between the river water and groundwater sampled, which showed a water exchange and mutual-supply relationship (the intersection occurred in June at Gulang and Tianzhu, and in August at Wushaoling). The correlation analysis between amount-weighted precipitation, river water, and groundwater suggested that the proportion of river water sourced from precipitation was almost as much as that exchanged with groundwater (Table 7). However, the correlation coefficients between river water and groundwater at Wushaoling and Tianzhu were both high, which demonstrated a very strong hydraulic connection between river water and groundwater. Moreover, the means of δ18O and δ2H for the river water and groundwater in Tianzhu were greater than in Gulang and Wushaoling, which reflected a drier climate (Table 8).

3.1.4. Relationships between the Isotope Signatures of Various Water Bodies

Figure 7 shows the linear relationship between δ18O and δ2H in precipitation, soil water, stem water (the plant species collected are presented in Table 9), river water, and groundwater. The isotope values of the different water bodies at the three sites were generally all below the LMWL, i.e., more enriched when compared with the isotope values of precipitation, which indicated that all of the water sources sampled, especially soil and xylem water, had undergone evaporation.
The isotope values of the river water in Gulang (Figure 7a) were generally located along the LMWL, which indicated that river water came from local precipitation [5]. The groundwater and soil water isotope values were quite similar, which suggested a recharge connection. Few of the stem water values are close to soil water values, which meant that the direct contribution of soil water to plant water was low. Since plant sampling sites were far from the rivers and groundwater, we could not trace the other plant water sources. The isotope values for soil water and river water in Wushaoling (Figure 7b) were situated between those for precipitation and groundwater, which suggested that the soil and river water originate from both precipitation and groundwater recharge. In Tianzhu (Figure 7c), the isotope values for river water lay between those of precipitation and groundwater, which indicated that river water resulted from precipitation. The specific recharge relationship between river water and groundwater requires further research. Stem water isotope values coincided with those of soil water, and were rather distinct from those of precipitation, which indicated that stem water came mainly from soil water.

3.2. The Fraction of Plant Water Uptake from Different Soil Depths

3.2.1. Isotopic Characteristics of Stem Water

The variations in δ18O and δ2H for the stem water of a plant genus (Caragana Fabr.) common to all three sites were analyzed as a case study (Figure 8). The isotope values for both δ18O and δ2H were higher at south slope than the north slope and mountain locations, which indicated higher evaporation on the south slope. The variation in δ18O ranged between −7.66‰ and −1.89‰, and was lower than that of δ2H, which ranged between −64.17‰ and −45.43‰, which suggested that the variation in δ18O was smaller than in δ2H. The maximum value of δ18O was observed in Gulang and Tianzhu in September, and in Wushaoling in July. The maximum δ2H value was observed in Gulang and Wushaoling in July, and in Tianzhu in August, which indicated that, although water was extracted from the same plant species, their distinct isotopic signatures reflect the differences in location and climate.
The major meteorological parameters that were analyzed in order to explain the relationship between stem water stable isotopes and meteorological conditions (Figure 9). δ18O was generally negatively correlated with temperature, water vapor, and precipitation, whereas δ2H responded to the variation in relative humidity.

3.2.2. Plant Water Use Pattern

As all plant sampling sites were located far from the rivers and groundwater, our measurements could only be used to calculate the fraction of water uptake from different soil depths, and the relationship between plant water and precipitation (if precipitation had occurred one week before sampling). Calculations were performed using the IsoSource software (Figure 10)—a stable isotope mixing model [22]. The original set of analytical solution was too vast to display, so the distribution of the default characteristic values (the mean, the minimum, the maximum, and the quartile values) are shown. Some factors, such as soil texture, vegetation coverage, and meteorological parameters, would result in different water use patterns by plants.
In different months, the water use pattern of the same plant were distinct. According to the mean value and the 50% quartile, the major water source for Caragana Fabr. at Gulang was from soil water, by month, at the 30–70 cm, 0–30 cm, 70–100 cm, 0–30 cm, and 0–30 cm layer from May to September, respectively (Figure 10). From these results, we found that in months when no precipitation occurred before sampling, the major water source was the 0–30 cm soil layer (in June, August and September). For the months when precipitation occurred before sampling, the more often that precipitation occurred, the deeper the soil water layer that the plant used (30–70 cm after three precipitation events in May, and 70–100 cm after four precipitation events in July). Overall, plants used direct precipitation less.
In August, when no precipitation occurred at Wushaoling (Figure 11), the major water source was also the 0–30 cm layer, the same as Gulang. In the remaining months, although there were several precipitation events, the water source was distinct for each month: the 0–30 cm layer in May, similar proportions from each soil layer in June, the precipitation event one day before sampling occurred in July, and the 70–100 cm layer in September.
For Tianzhu, in the no-precipitation month the result was the same as the above for the other two locations (Figure 12). In the other months, the water use pattern varied. The contribution of each water source to the plant water was similar in May with one precipitation event, but, in June, the proportion of water uptake from the 0–30 cm layer increased. The contribution from each water source to plant water was small in July, even with six precipitation events, but in September, with five precipitation events, the major source of plant water was the 0–30 cm layer.

4. Discussion

With regards to the water extraction system for soil and plant, Orlowski et al. [26] conducted an inter-laboratory comparison of cryogenic water extraction systems for stable isotope analysis of soil water. They showed large differences in isotopic signatures among participating laboratories, with mean differences ranging from +18.1‰ to −108.4‰ for δ2H and +11.8‰ to −14.9 ‰ for δ18O across all laboratories when compared to the reference water. The observed differences in isotope data were influenced by interactions between multiple factors (soil type and properties, soil water content, system setup, extraction efficiency, extraction system leaks, and the lab’s internal accuracy). They suggested that laboratories might have to establish calibration functions for their specific extraction system for each natural, individual soil type.
Here, the results produced by our extraction system (LI-2100, LICA United Technology Limited, China) are evaluated, including the extraction efficiency and isotopic variation, after extraction compared to the reference water. According to its user manual, the extraction system considers temperature, vacuum, extraction time, and set build-in parameters, which users choose according to sample types.
Our test samples, including the reference water (tap water) and the soil samples with different SWC (20%, 50% and 80%) filled by reference water, all had three repeats. The information on extraction efficiency is presented in Table 10, from which the extraction rate was an average of 98%, so the water loss was very small. The change in isotope values after extraction are shown in Figure 13. The greatest difference in isotope value, before and after extraction, was 1.05 ‰ for δ18O and 6.66 ‰ for δ2H, which were both lower than the mean differences reported by Orlowski et al. [26] when conducted by the 16 laboratories. Thus, the extraction results of our system were credible. In addition, the higher the SWC, the closer the isotope values were when compared to the original reference water after extraction.
It should be noted that the soil water and the plant stem water measured by the Liquid Water Isotope Analyzer need to be detected for water contamination. The LWIA Spectral Contamination Identifier (LWIA-SCI) software is a simple tool, provided by Los Gatos Research (LGR), to identify features that are consistent with water contamination in the spectra recorded by LGR’s liquid water isotope analyzer. The LWIA-SCI analyzes the data files and the recorded spectra from the LWIA to produce a metric that indicates the likelihood of contamination from either narrow-band or broad-band absorbers. This metric is compared to the known clean samples (general standards) in the data set and produces a good or bad flag to indicate the reliability of the associated isotopic ratio measurements.
Finally, we should note that, since only a one-year dataset was used and the number of sampling sites was limited, longer datasets and more sampling sites are required to confirm these findings. Also, the use of water by plants from different soil depths was studied, without considering other water sources. In future sampling campaigns, additional sampling points near rivers and wells should be considered, in order to take into account all of the potential water sources and the proportion of their contribution to plant water. In order to better understand the changes in the isotopic signatures of different water bodies before and after precipitation, especially the change in water use strategy by plants, several typical precipitation events should be selected and samples collected before and after the precipitation event.

5. Conclusions

In this study we selected three sampling sites, specifically Gulang, Wushaoling and Tianzhu, located in the northeastern Tibetan Plateau. Event-based precipitation samples were collected for a complete year. Monthly samples of plant water, soil water, river water, and groundwater were collected for a complete growing season. The isotopic characteristics of the various water bodies, the recharge relationships among them, and the fraction of water uptake from different soil layers were studied.
The “temperature effect” was apparent at the three sites, suggesting an evaporation enrichment. Using the six regressions (OLSR, MA, RMA, PWLSR, PWMA, PWRMA), the slope differences among the various methods were tested, and were the smallest for Wushaoling, followed by Tianzhu, with the largest difference in Gulang. The δ18O of soil water varied significantly in the topsoil layers, while they tended to be similar in deep soil. The isotopes became increasingly depleted along the vertical soil profile, from top to bottom. The changes in monthly δ18O and δ2H were minor, but soil temperature visibly changed, which showed that temperature was not the sole factor that impacted soil water isotopes. The soil water isotopes at the different sites had different responses to precipitation, and they had variable soil layer boundaries, where Wushaoling’s soil boundary layers were located at 20 and 60 cm below the surface, while for Tianzhu this was located at 40 cm. The consistent variation in δ18O and δ2H for river water and groundwater indicated a similar source water. The similarities between these two water bodies showed a water exchange process. The relationship among precipitation, river water, and groundwater was distinct at all three sites. The hydraulic connection of river water and groundwater was stronger at Wushaoling and Tianzhu than at Gulang. Using the isotope data from the xylem of the same plant genus (Caragana Fabr.) sampled at three sites and applying the IsoSource model, it was found that the water use pattern of plants varied with both month and site location, but was the same between sites in no-precipitation months, i.e., plant water relied on soil water at 0–30 cm.

Author Contributions

Conceptualization, M.Z.; methodology, X.Q.; investigation, R.G.; resources, R.C. and H.M.; writing—original draft preparation, X.Q.; writing—review and editing, S.W. and A.A.

Funding

This research was funded by National Natural Science Foundation of China, grant number 41461003; 41771035 and Scientific Research Program of Higher Education Institutions of Gansu Province, grant number 2018C-02.

Acknowledgments

The authors greatly thank all the meteorological stations for their warm help in field sampling. We also thank the colleagues in the Northwest Normal University for in laboratory analyzing and field work, especially thank Xiuxiu Yu, Sue Zhou and Yaning Zhang.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Tan, H.; Wen, X.; Rao, W.; Bradd, J.; Huang, J. Temporal variation of stable isotopes in a precipitation-groundwater system: Implications for determining the mechanism of groundwater recharge in high mountain-hills of the Loess Plateau, China. Hydrol. Process. 2016, 30, 1491–1505. [Google Scholar] [CrossRef]
  2. He, S.; Richards, K. Stable isotopes in monsoon precipitation and water vapour in Nagqu, Tibet, and their implications for monsoon moisture. J. Hydrol. 2016, 540, 615–622. [Google Scholar] [CrossRef]
  3. Hu, H.; Dominguez, F. Evaluation of oceanic and terrestrial sources of moisture for the north American monsoon using numerical models and precipitation stable isotopes. J. Hydrometeorol. 2015, 16, 19–35. [Google Scholar] [CrossRef]
  4. Evaristo, J.; McDonnell, J.J.; Scholl, M.A.; Bruijnzeel, L.A.; Chun, K.P. Insights into plant water uptake from xylem-water isotope measurements in two tropical catchments with contrasting moisture conditions. Hydrol. Process. 2016, 30, 3210–3227. [Google Scholar] [CrossRef]
  5. Yang, B.; Wen, X.; Sun, X. Irrigation depth far exceeds water uptake depth in an oasis cropland in the middle reaches of Heihe River Basin. Sci. Rep. 2015, 5, 15206. [Google Scholar] [CrossRef] [Green Version]
  6. Yang, B.; Wen, X.; Sun, X. Seasonal variations in depth of water uptake for a subtropical coniferous plantation subjected to drought in an East Asian monsoon region. Agric. For. Meteorol. 2015, 201, 218–228. [Google Scholar] [CrossRef]
  7. Wan, H.; Liu, W. An isotope study (18O and D) of water movements on the Loess Plateau of China in arid and semiarid climates. Ecol. Eng. 2016, 93, 226–233. [Google Scholar] [CrossRef]
  8. Lee, K.S.; Kim, J.M.; Lee, D.R.; Kim, Y.; Lee, D. Analysis of water movement through an unsaturated soil zone in Jeju Island, Korea using stable oxygen and hydrogen isotopes. J. Hydrol. 2007, 345, 199–211. [Google Scholar] [CrossRef]
  9. Marchina, C.; Bianchini, G.; Natali, C.; Pennisi, M.; Colombani, N.; Tassinari, R.; Knoeller, K. The Po river water from the Alps to the Adriatic Sea (Italy): New insights from geochemical and isotopic (δ18O-δD) data. Environ. Sci. Pollut. Res. 2015, 25, 5184–5203. [Google Scholar] [CrossRef] [PubMed]
  10. Chiogna, G.; Santoni, E.; Camin, F.; Tonon, A.; Majone, B.; Trenti, A.; Bellin, A. Stable isotope characterization of the Vermigliana catchment. J. Hydrol. 2014, 509, 295–305. [Google Scholar] [CrossRef]
  11. Liu, J.; Chen, Z.; Zhang, Y.; Liu, F. Stable isotope evidences on sources and mechanisms of groundwater recharge in Hohhot basin, China. Environ. Earth Sci. 2016, 75, 410–420. [Google Scholar] [CrossRef]
  12. Jasechko, S.; Birks, S.J.; Gleeson, T.; Wada, Y.; Fawcett, P.J.; Sharp, Z.D.; McDonnell, J.J.; Welker, J.M. The pronounced seasonality of global groundwater recharge. Water Resour. Res. 2014, 50, 8845–8867. [Google Scholar] [CrossRef] [Green Version]
  13. Li, Z.; Gui, J.; Wang, X.; Feng, Q.; Zhao, T.; Ouyang, C.; Guo, X.; Zhang, B.; Shi, Y. Water resources in inland regions of central Asia: Evidence from stable isotope tracing. J. Hydrol. 2019, 70, 1–16. [Google Scholar] [CrossRef]
  14. Li, Z.; Feng, Q.; Yong, S.; Wang, Q.J.; Jiao, Y.; Li, Y.; Li, J.; Guo, X. Stable isotope composition of precipitation in the south and north slopes of Wushaoling Mountain, northwestern China. Atmos. Res. 2016, 182, 87–101. [Google Scholar] [CrossRef]
  15. Ma, J.; Zhang, P.; Zhu, G.; Wang, Y.; Edmunds, W.M.; Ding, Z.; He, J. The composition and distribution of chemicals and isotopes in precipitation in the Shiyang River system, northwestern China. J. Hydrol. 2012, 436–437, 92–101. [Google Scholar] [CrossRef]
  16. Zhao, L.; Yin, L.; Xiao, H.; Cheng, G.; Zhou, M.; Yang, Y.; Li, C.; Zhou, J. Isotopic evidence for the moisture origin and composition of surface runoff in the headwaters of the Heihe River basin. Chin. Sci. Bull. 2011, 56, 406–416. [Google Scholar] [CrossRef]
  17. Wu, J.; Ding, Y.; Ye, B.; Yang, Q.; Zhang, X.; Wang, J. Spatio-temporal variation of stable isotopes in precipitation in the Heihe River Basin, Northwestern China. Environ. Earth Sci. 2010, 61, 1123–1134. [Google Scholar] [CrossRef]
  18. Wang, Y.; Wang, J.; Qi, Y.; Yan, C. Dataset of Desert Distribution in China (1:100,000). Cold and Arid Regions Science Data Center at Lanzhou (CARD), 2005. Available online: http://doi.org/10.3972/westdc.006.2013.db (accessed on 30 December 2018).
  19. IAEA. Statistical Treatment of Data on Environmental Isotopes in Precipitation; Technical Report Series; IAEA: Vienna, Austria, 1992; Volume 331, p. 781. [Google Scholar]
  20. Crawford, J.; Hughes, C.E.; Lykoudis, S. Alternative least squares methods for determining the meteoric water line, demonstrated using GNIP data. J. Hydrol. 2014, 519, 2331–2340. [Google Scholar] [CrossRef]
  21. Hughes, C.E.; Crawford, J. Anew precipitation weighted method for determining the meteoric water line for hydrological applications demonstrated using Australian and global GNIP data. J. Hydrol. 2012, 464–465, 344–351. [Google Scholar] [CrossRef]
  22. Phillips, D.L.; Gregg, J.W. Source partitioning using stable isotopes: Coping with too many sources. Oecologia 2003, 136, 261–269. [Google Scholar] [CrossRef]
  23. Dansgaard, W. Stable isotopes in precipitation. Tellus 1964, 16, 436–468. [Google Scholar] [CrossRef]
  24. Craig, H. Isotopic variations in meteoric waters. Science 1961, 133, 1702–1703. [Google Scholar] [CrossRef]
  25. Gat, J.R. Oxygen and hydrogen isotopes in the hydrologic cycle. Annu. Rev. Earth Planet. Sci. 1996, 24, 225–262. [Google Scholar] [CrossRef]
  26. Orlowski, N.; Breuer, L.; Angeli, N.; Boeckx, P.; Brumbt, C.; Cook, C.S.; Dubbert, M.; Dyckmans, J.; Gallagher, B.; McDonnell, J.J.; et al. Inter-laboratory comparison of cryogenic water extraction systems for stable isotope analysis of soil water. Hydrol. Earth Syst. Sci. 2018, 22, 3619–3637. [Google Scholar] [CrossRef]
Figure 1. A map showing the locations of the sampling sites in the Qilian Mountains. The boundary of the mountainous area is marked by the black line. The satellite-derived background was acquired from Natural Earth (http://www.naturalearthdata.com). The desert areas are based on Wang et al. [18]. The inset map identifies the location of the Qilian Mountains relative to the Tibetan Plateau.
Figure 1. A map showing the locations of the sampling sites in the Qilian Mountains. The boundary of the mountainous area is marked by the black line. The satellite-derived background was acquired from Natural Earth (http://www.naturalearthdata.com). The desert areas are based on Wang et al. [18]. The inset map identifies the location of the Qilian Mountains relative to the Tibetan Plateau.
Water 11 00770 g001
Figure 2. The Local Meteorological Water Line (LMWL) at the three sites; the Global Meteoric Water Line (GMWL) was also shown for reference. (ac) are the LMWL at Gulang, Wushaoling, and Tianzhu, respectively.
Figure 2. The Local Meteorological Water Line (LMWL) at the three sites; the Global Meteoric Water Line (GMWL) was also shown for reference. (ac) are the LMWL at Gulang, Wushaoling, and Tianzhu, respectively.
Water 11 00770 g002
Figure 3. Mean δ18O, δ2H, and soil water content (SWC) for soil water at a depth of 0–100 cm in Gulang (ac), Wushaoling (df) and Tianzhu (gi).
Figure 3. Mean δ18O, δ2H, and soil water content (SWC) for soil water at a depth of 0–100 cm in Gulang (ac), Wushaoling (df) and Tianzhu (gi).
Water 11 00770 g003
Figure 4. Monthly variation of SWC-weighted soil water δ18O and δ2H, and near-surface temperature. (ac) are the variations at Gulang, Wushaoling and Tianzhu, respectively.
Figure 4. Monthly variation of SWC-weighted soil water δ18O and δ2H, and near-surface temperature. (ac) are the variations at Gulang, Wushaoling and Tianzhu, respectively.
Water 11 00770 g004
Figure 5. The response of soil water, δ18O, and δ2H to different amounts. (ac) are the response of SWC, δ18O, and δ2H to precipitation pulsse at Wushaoling, and (df) are for the responses at Tianzhu.
Figure 5. The response of soil water, δ18O, and δ2H to different amounts. (ac) are the response of SWC, δ18O, and δ2H to precipitation pulsse at Wushaoling, and (df) are for the responses at Tianzhu.
Water 11 00770 g005
Figure 6. δ18O and δ2H in river water and groundwater where (a,b) are the δ18O and δ2H measurements, respectively, in river water and groundwater at Gulang; (c,d) are the measurements for Wushaoling; and (e,f) are the measurements for Tianzhu.
Figure 6. δ18O and δ2H in river water and groundwater where (a,b) are the δ18O and δ2H measurements, respectively, in river water and groundwater at Gulang; (c,d) are the measurements for Wushaoling; and (e,f) are the measurements for Tianzhu.
Water 11 00770 g006
Figure 7. The linear relationship between δ18O and δ2H of the different water bodies sampled in (a) Gulang, (b) Wushaoling, and (c) Tianzhu.
Figure 7. The linear relationship between δ18O and δ2H of the different water bodies sampled in (a) Gulang, (b) Wushaoling, and (c) Tianzhu.
Water 11 00770 g007
Figure 8. The monthly variations in stem water isotopes. (a) is δ18O variation and (b) is δ2H variation.
Figure 8. The monthly variations in stem water isotopes. (a) is δ18O variation and (b) is δ2H variation.
Water 11 00770 g008
Figure 9. The monthly variations in major meteorological parameters: (a) is air temperature, (b) is relative humidity, (c) is vapor pressure, and (d) is precipitation amount.
Figure 9. The monthly variations in major meteorological parameters: (a) is air temperature, (b) is relative humidity, (c) is vapor pressure, and (d) is precipitation amount.
Water 11 00770 g009
Figure 10. The fraction of plant water use from the different soil layers (0–30 cm, 30–70 cm, 70–100 cm) and precipitation at Gulang. (ae) represent the months of May, June, July, August, and September, respectively.
Figure 10. The fraction of plant water use from the different soil layers (0–30 cm, 30–70 cm, 70–100 cm) and precipitation at Gulang. (ae) represent the months of May, June, July, August, and September, respectively.
Water 11 00770 g010
Figure 11. The fraction of plant water use from the different soil layers (0–30 cm, 30–70 cm, 70–100 cm) and precipitation at Wushaoling. (ae) represent the months of May, June, July, August, and September, respectively.
Figure 11. The fraction of plant water use from the different soil layers (0–30 cm, 30–70 cm, 70–100 cm) and precipitation at Wushaoling. (ae) represent the months of May, June, July, August, and September, respectively.
Water 11 00770 g011
Figure 12. The fraction of plant water use from the different soil layers (0–30 cm, 30–70 cm, 70–100 cm) and precipitation at Tianzhu. (ae) represents May, Jun, July, August and September, respectively.
Figure 12. The fraction of plant water use from the different soil layers (0–30 cm, 30–70 cm, 70–100 cm) and precipitation at Tianzhu. (ae) represents May, Jun, July, August and September, respectively.
Water 11 00770 g012
Figure 13. The isotopic differences for the soil samples with different SWC and the reference water. (a,b) are for δ18O and δ2H, respectively.
Figure 13. The isotopic differences for the soil samples with different SWC and the reference water. (a,b) are for δ18O and δ2H, respectively.
Water 11 00770 g013
Table 1. The locations and long term climatology (T = air temperature and P = precipitation amount) during 1981–2010 at the sampling sites in the Qilian Mountains.
Table 1. The locations and long term climatology (T = air temperature and P = precipitation amount) during 1981–2010 at the sampling sites in the Qilian Mountains.
Sampling SitesLatitudeLongitudeAltitude (m)T (°C)P (mm)
Gulang37°29′102°54′2072.45.7352.3
Wushaoling37°12′102°52′3045.10.3407.4
Tianzhu36°59′103°11′2484.83.5450.0
Table 2. Number of samples collected in the Qilian Mountains from 2015 to 2016.
Table 2. Number of samples collected in the Qilian Mountains from 2015 to 2016.
Sampling SitesPrecipitationPlantSoilRiverGroundwater
Gulang601516555
Wushaoling1562016555
Tianzhu1071316555
Table 3. Isotopic signatures of precipitation in the Qilian Mountains from 2015 to 2016.
Table 3. Isotopic signatures of precipitation in the Qilian Mountains from 2015 to 2016.
δ18O (‰)δ2H (‰)D-Excess (‰)
Gulang−19.91−2.40−135.86−14.28−10.96−29.41
(mean −8.58)(mean −54.21)(mean 14.46)
Wushaoling−31.11−2.58−232.15−28.64−8.77−31.43
(mean −8.61)(mean −54.78)(mean 14.10)
Tianzhu−19.67−7.89−158.44−60.81−10.01−31.46
(mean −7.53)(mean −48.48)(mean 11.79)
Table 4. The relationship between δ18O and precipitation amount (P), and between δ18O and temperature (T).
Table 4. The relationship between δ18O and precipitation amount (P), and between δ18O and temperature (T).
Siteδ18O−Pδ18O−T
Correlation CoefficientRegression Coefficient (‰·mm−1)Correlation CoefficientRegression Coefficient (‰·oC−1)
Gulang−0.15−0.110.402 *0.124
P = 0.397 > 0.05P = 0.017 < 0.05
Wushaoling0.0950.1130.483 **0.226
P = 0.289 > 0.05P = 0.000
Tianzhu−0.096−0.0660.562 **0.362
P = 0.367 > 0.05P = 0.000
* means statistically significant at the 0.05 level. ** means statistically significant at the 0.01 level.
Table 5. Parameters of the fit for OLSR (ordinary least squares regression), RMA (reduced major axis regression), MA (major axis regression), PWLSR (precipitation weighted least squares regression), PWRMA (precipitation weighted reduced major axis regression), PWMA (precipitation weighted major axis regression) regressions at each site in the Qilian Mountains.
Table 5. Parameters of the fit for OLSR (ordinary least squares regression), RMA (reduced major axis regression), MA (major axis regression), PWLSR (precipitation weighted least squares regression), PWRMA (precipitation weighted reduced major axis regression), PWMA (precipitation weighted major axis regression) regressions at each site in the Qilian Mountains.
SitesAlgorithmSlopeInterceptrmSSEavt-Value
Gulang N = 43OLSR8.31 ± 0.3218.86 ± 3.151.01-
RMA8.56 ± 0.3121.17 ± 3.081.000.55
MA8.81 ± 0.3323.48 ± 3.241.011.09
PWLSR9.01 ± 0.3624.24 ± 3.551.001.50
PWRMA9.31 ± 0.3727.00 ± 3.581.032.07
PWMA9.60 ± 0.3729.78 ± 3.661.082.59
Wushaoling N = 127OLSR8.09 ± 0.1014.99 ± 1.021.00-
RMA8.17 ± 0.1015.64 ± 1.021.000.54
MA8.24 ± 0.1016.27 ± 1.031.001.06
PWLSR8.09 ± 0.1215.21 ± 1.111.000.01
PWRMA8.21 ± 0.1216.15 ± 1.121.000.84
PWMA8.32 ± 0.1317.07 ± 1.131.011.65
Tianzhu N = 90OLSR7.83 ± 0.1611.18 ± 1.451.01-
RMA7.96 ± 0.1512.22 ± 1.431.000.61
MA8.10 ± 0.1613.25 ± 1.471.011.21
PWLSR8.05 ± 0.1712.60 ± 1.561.001.01
PWRMA8.20 ± 0.1713.89 ± 1.571.011.68
PWMA8.35 ± 0.1715.17 ± 1.591.032.32
N is the number of data points used in the regression, rmSSEav is the average value of root mean sum of squared error and t-value is a value used in a t-test of the significance in the difference of each regression to that of OLSR.
Table 6. Correlation coefficients between soil water δ18O and soil water content (SWC).
Table 6. Correlation coefficients between soil water δ18O and soil water content (SWC).
MonthMayJuneJulyAugustSeptember
Gulang0.781 **−0.733 *−0.291−0.655 *−0.057
Wushaoling−0.125−0.5240.2270.907 **−0.709 *
Tianzhu0.589−0.640 *0.040−0.571−0.024
* means statistically significant at the 0.05 level. **means statistically significant at the 0.01 level.
Table 7. The correlation relationships between precipitation, river water, and groundwater.
Table 7. The correlation relationships between precipitation, river water, and groundwater.
SitesPrecipitation-River WaterPrecipitation-GroundwaterRiver Water-Groundwater
Gulang0.3860.1180.358
Wushaoling0.4130.1420.893 *
Tianzhu−0.455−0.3980.798
* means statistically significant at the 0.05 level.
Table 8. The mean δ18O and δ2H for river water and groundwater.
Table 8. The mean δ18O and δ2H for river water and groundwater.
SiteRiver WaterGroundwater
δ18O (‰)δ2H (‰)δ18O (‰)δ2H (‰)
Gulang−8.19144−50.7811−8.36746−52.4299
Wushaoling−8.05623−49.566−8.34914−53.3473
Tianzhu−7.62581−47.0895−7.65622−47.5145
Table 9. Identification of the plant species collected in the Qilian Mountains in 2016.
Table 9. Identification of the plant species collected in the Qilian Mountains in 2016.
Sampling SitesLatin Name of Plant Species
GulangArmeniaca vulgaris Lam.
Caragana korshinskii Kom.
Populus davidiana
WushaolingCaragana stenophylla Pojark.
Iris lactea Pall. Var. chinensis (Fisch.) Koidz.
Achnatherum inebrians (Hance) Keng
TianzhuPopulus davidiana
Caragana stenophylla Pojark.
Ceratoides arborescens (Losinsk.) Tsien et C.G.Ma
Picea crassifolia
Table 10. Information about the extraction efficiency of LI-2100.
Table 10. Information about the extraction efficiency of LI-2100.
NumberSoil Weight (g)Water Weight Filled by Reference Water (g)Water Weight after Extraction (g)Extraction Rate (%)
20-115.743.153.1299
20-214.983.002.9498
20-314.522.902.7695
50-18.884.444.42100
50-210.245.125.10100
50-37.973.983.9298
80-17.455.965.8097
80-26.975.585.5099
80-37.546.035.9398

Share and Cite

MDPI and ACS Style

Qiu, X.; Zhang, M.; Wang, S.; Argiriou, A.A.; Chen, R.; Meng, H.; Guo, R. Water Stable Isotopes in an Alpine Setting of the Northeastern Tibetan Plateau. Water 2019, 11, 770. https://doi.org/10.3390/w11040770

AMA Style

Qiu X, Zhang M, Wang S, Argiriou AA, Chen R, Meng H, Guo R. Water Stable Isotopes in an Alpine Setting of the Northeastern Tibetan Plateau. Water. 2019; 11(4):770. https://doi.org/10.3390/w11040770

Chicago/Turabian Style

Qiu, Xue, Mingjun Zhang, Shengjie Wang, Athanassios A. Argiriou, Rong Chen, Hongfei Meng, and Rong Guo. 2019. "Water Stable Isotopes in an Alpine Setting of the Northeastern Tibetan Plateau" Water 11, no. 4: 770. https://doi.org/10.3390/w11040770

APA Style

Qiu, X., Zhang, M., Wang, S., Argiriou, A. A., Chen, R., Meng, H., & Guo, R. (2019). Water Stable Isotopes in an Alpine Setting of the Northeastern Tibetan Plateau. Water, 11(4), 770. https://doi.org/10.3390/w11040770

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