Next Article in Journal
A Coupled Model for Simulating Water and Heat Transfer in Soil-Plant-Atmosphere Continuum with Crop Growth
Next Article in Special Issue
Isotopes in Hydrology and Hydrogeology
Previous Article in Journal
Recent Advances in Water and Wastewater Treatment with Emphasis in Membrane Treatment Operations
Previous Article in Special Issue
Waters from the Djiboutian Afar: A Review of Strontium Isotopic Composition and a Comparison with Ethiopian Waters and Red Sea Brines
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Classification of Stable Isotopes and Identification of Water Replenishment in the Naqu River Basin, Qinghai-Tibet Plateau

1
School of Hydraulic Engineering, Dalian University of Technology, Dalian 116024, China
2
Department of Water Conservancy Engineering, North China University of Water Resources and Electric Power, Zhengzhou 450046, China
3
Collaborative Innovation Center of Water Resources Efficient Utilization and Support Engineering, Zhengzhou 450046, China
*
Author to whom correspondence should be addressed.
Water 2019, 11(1), 46; https://doi.org/10.3390/w11010046
Submission received: 29 August 2018 / Revised: 28 November 2018 / Accepted: 24 December 2018 / Published: 28 December 2018
(This article belongs to the Special Issue Isotopes in Hydrology and Hydrogeology)

Abstract

:
The stable isotopic study of the mechanism of runoff replenishment in the Qinghai-Tibet Plateau is a time-consuming and complicated process requiring complex monitoring data and scientific evaluation methods. Based on the data of water stable isotopes (18O and 2H) in the Naqu River basin, the present paper developed a framework of the variable fuzzy evaluation model (VFEM) to provide a method to classify stable isotopes and generalize the source identification of water replenishment by rainfall or snowmelt in the Naqu River basin. The grade eigenvalues of tributaries were ranked from low to high as follows: 1, 1.005, 1.089, 1.151, 1.264, 1.455 and 2.624. Three sets of tributaries were distinguished. The grade eigenvalues of the Najinqu, Bazongqu, Mumuqu, Chengqu and Gongqu Rivers were small, indicating that these tributaries were strongly supplemented by precipitation and snowmelt; the grade eigenvalue of the Zongqingqu River was in the medium range (1.455); the third group included the Mugequ River with a high status value (2.624). This study mainly highlighted the combination of the classification of stable isotopes and plots of δ2H vs. δ18O in the source identification of water replenishment, which will be helpful for studying runoff replenishment and the evolution mechanism in the Qinghai-Tibet Plateau.

1. Introduction

The study of water cycle evolution to solve water resource and related environmental issues has become a significant task for hydrological observation processes. However, traditional hydrological observation is inconvenient in ungauged catchments, especially in the arid region with extreme climate conditions [1].
Variations in the stable isotopes of water, such as 18O and 2H, are widely used as environmental tracers in studying hydrological processes [2,3,4,5,6]. Short and long-term isotope signals in precipitation are thus transmitted through the whole catchment [7]. Isotopes of particular interest for hydrological studies include the stable isotopes of water (18O, 2H), which are incorporated within the water molecule (H218O, 1H2H16O), and exhibit systematic spatial and temporal variations as a result of isotope fractionations that accompany water-cycle phase changes and diffusion [8]. Stable- and radio-isotope tracers have been widely applied in Earth systems studies, including hydrological and climatological research, for their ability to provide a sharper focus on some of the underlying processes that control the chemical and physical behaviour of elements and compounds in nature. Spatial and temporal variations in hydrogen and oxygen isotopic compositions (δ2H and δ18O) could be attributed to equilibrium exchanges occurring at all phase transitions and kinetic isotopic fractionation during transformation from the ocean surface to the atmosphere, small rainfall amounts when rain drops fall below the cloud [9]. For example, ratios of stable isotopes of hydrogen and oxygen (2H/1H and 18O/16O) in river waters were measured to investigate the hydrological pathway of the Xijiang River, in southwest China [10]; groundwater was found to be caused principally by the mixing of rainfall and water stored before the event in the Meilin watershed [11]; δ2H and δ18O were used for evaluating the source of water [12,13,14,15] and exploring the interaction between surface water and groundwater [16,17,18,19] by using the change in isotopic composition due to evaporation or mixing by different waters.
The Qinghai-Tibet Plateau has long been the focus of international academic concern because it possesses a unique high plateau climate and is very sensitive to global climate change and human activities. Isotope studies on the processes of rivers and lakes in the plateau area have achieved significant progress since the 1960s [20]. Zhang et al. [21] investigated the variation of δ2H and δ18O in glacial meltwater at Mount Everest and found that δ2H in ice and snow is larger than that in river and lake water compared to the change of δ18O. Regarding precipitation, the variation of δ18O is subject to the effects of both temperature and rainfall quantity [22]. By analyzing the variation of δ18O in the Naqu River basin in the plateau, Tian et al. [23] proved that the effects of lake water and river runoff on stable isotopes differ considerably. Pang et al. [24] undertook a detailed examination of the isotope fractionation process occurring at the ice–snow phase transition zones at the Baishui Glacier. While isotopic characteristics of river water have been described in general, few studies have investigated the classification of stable isotopes for the source identification of water replenishment.
Due to the particularity of its geographical location and the lack of the data from hydrological stations, it is very difficult to explore the hydrological cycle in the Naqu River basin. The variable fuzzy evaluation model (VFEM) can reasonably identify the relationship degree and relationship function between the sample index and standard interval of each level index. The assessment level of each sample by varying the model and its parameters are properly determined, so that the reliability of assessment can be improved. The proposed method has been applied to assess the water resources carrying capacity of the Huaihe River basin [25].
This paper presents a set of criteria using VFEM to determine the classification of stable isotopes for the identification of water replenishment in the Naqu River basin. It can be concluded that VFEM is useful for the classification of water replenishment. The research result may provide support to runoff replenishment by rainfall or snowmelt in the Naqu River basin.

2. Description of the Study Area

The Naqu River basin is located in the North Tibet Plateau in the Tibet Autonomous Region of China with its domain between 30°54′–32°43′ N and 91°12′–92°54′ E (Figure 1), and with a catchment area of 16,350 km2. The average elevation of the study area is higher than 4600 m. The length of the river is 460 km with a 920-m drop. The mean annual rainfall is 477.8 mm, and the mean annual temperature is −0.6 °C. The annual rainfall is concentrated in the summer (from June to October) and accounts for about 81.9% of the total annual rainfall. Thus, the period from June to October has been defined as the wet season, and the remaining months as the dry season. Although snowfall mostly occurs in the non-flood season and total snowfall accounts for less than 20% of the precipitation amount, snowmelt flow contributes significantly to the streamflow in May and June. The Mugequ, Chengqu and Gongqu Rivers are the main sub-basins of the Naqu River basin and have catchment areas of 2103, 1090 and 1232 km2, respectively. The average annual runoff of the Naqu River basin is 3.1 billion cubic meters.

3. Materials and Methods

3.1. Sampling and Analytical Methods

River and soil waters were sampled on 20–23 August 2016 in the Naqu River basin, and further sample collection from upstream to downstream in the Naqu River basin was conducted on 4–10 January 2017. In total, we collected 27 samples of surface water (14 samples in January, 13 samples in August), 13 samples of soil water in August, and 5 samples of snowmelt in January. Thirteen sample sites were chosen around the mainstream and tributary of the Naqu River basin. The sampling sites are shown in Figure 1. Twelve water samples were collected from the main stream and three tributaries, and then stored in glass sampling bottles for 2H and 18O analysis.
Twelve soil samples were also collected near the water sampling sites for water extraction. The soil was drilled at the top layer of soil (0–10 cm) and then immediately stored in polytetrafluoroethylene bottles to prevent evaporation. Soil samples were frozen and then thawed overnight before water was extracted using the cryogenic vacuum distillation method.
Water samples were analyzed for 2H and 18O through wavelength-scanned cavity ring down spectroscopy (WS-CRDS) (Picarro L1115-I, Picarro, Santa Clara, CA, USA). The δ2H and δ18O were corrected using VSMOW (Vienna Standard Mean Ocean Water) (δ2H = 0‰, δ18O = 0‰) and Standard Light Antarctic Precipitation (δ2H = −428‰, δ18O = −55.5‰). The analytical precision was generally 0.5‰ for δ2H and 0.1‰ for δ18O. (2H/1H)sample or (18O/16O)sample is the stable isotope ratio in the sampled water, and (2H/1H)standard or (18O/16O)standard is the stable isotope ratio in standard mean ocean water. The final results are presented as
δ 2 H V S M O W = ( 2 H / 1 H s a m p l e 2 H / 1 H s t a n d a r d 1 ) × 1000 ( )
δ 18 O V S M O W = ( 18 O / 16 O s a m p l e 18 O / 16 O s t a n d a r d 1 ) × 1000 ( )

3.2. Methods

The concept of fuzzy sets depicting imprecision or vagueness was introduced by Zadeh, which gives fuzziness a scientific description and has generated great interest in academic circles. An extended Zadeh’s fuzzy set theory and the engineering fuzzy set theory were established according to relative membership function. Afterwards, the theory of variable fuzzy sets (VFS) was presented by Chen [26], which can scientifically and reasonably determine membership degrees and relative membership functions of disquisitive objectives (or indices) at level interval relating to the engineering system. This can fully use one’s experience and knowledge as well as qualitative and quantitative information of the index system to obtain weights of objectives (or indices) [27]. In this paper, a VFEM based on variable fuzzy sets was proposed to identify the classification of stable isotopes in the Naqu River basin, and the theory and methodology of the proposed model is presented here.

3.2.1. The Definition of Variable Fuzzy Sets

In order to define the concept of VFS, we suppose that U is a fuzzy concept, and A and A C represent attractability and repellency, respectively. To any elements u ( u U ), μ A ( u ) and μ A C ( u ) are the relative membership function of elements u to A and A C that express degrees of attractability and repellency, respectively, where, μ A ( u ) + μ A C ( u ) = 1 , and 0 μ A ( u ) 1 , 0 μ A C ( u ) 1 .
Let
D A ( u ) = μ A ( u ) μ A C ( u )
where D A ( u ) is defined as the relative difference degree of u to A .
Then, mapping,
{ D A ( u ) : D [ 1 , 1 ] u | D A ( u ) [ 1 , 1 ]
is defined as the relative difference function of u to A . Since μ A ( u ) + μ A C ( u ) = 1 , then:
D A ( u ) = 2 μ A ( u ) 1
Let,
V ~ = { ( u , D ) | u U , D ( u ) = μ A ( u ) μ A C ( u ) , D [ 1 , 1 ] }
A + = { u | u U , 0 < D ( u ) < 1 }
A = { u | u U , 1 < D ( u ) < 0 }
A 0 = { u | u U , D ( u ) = 0 }
where V ~ denotes VFS, and A + , A , and A 0 indicate attracting sets, repelling sets and the balance and qualitative change boundary of VFS, respectively [28]. The details of the definition of VFS are also available in these references.

3.2.2. The Relative Membership Degree of Indices

Suppose the object evaluated has m indices eigenvalue as follows:
X = { x 1 , x 2 , , x m }
The indices can be evaluated by c levels, and the criteria interval matrices of m indices and c levels can be expressed as follows:
I a b = ( [ a i h , b i h ] )
where i = 1 , 2 , , m ;   h = 1 , 2 , , c . Level 1 means the superior level and level c means the inferior level. If the index is the larger the better, then a > b . If the index is the smaller the better, then a < b . For every [ a i h , b i h ] , the range of interval [ c i h , d i h ] can be determined according to the upper and lower bound of its adjacent intervals.
I c d = ( [ c i h , d i h ] )
(1) When h = 1 , then its criteria interval of the index I a b = [ a i 1 , b i 1 ] . Obviously, for A , the relative membership of the upper bound a i 1 , μ A ( a i 1 ) equals 1, and the relative membership degree of the lower bound b i 1 , μ A ( b i 1 ) equals 0.5.
Suppose that M i 1 is a point value belonging to the interval [ a i 1 , b i 1 ] , and its relative membership degree to A is 1. Since level 1 is the superior level, and the relative membership degree of the left bound point a i 1 to A is 1, then M i 1 = a i 1 ; that is, M i 1 is located at the left bound point of [ a i 1 , b i 1 ] .
When the index eigenvalue x i of the object evaluated belongs to the interval [ a i 1 , b i 1 ] , x i must be located at the right of M i 1 . The relative membership degree of x i to level 1 can be calculated as follows:
μ A ( x i ) 1 = 0.5 [ 1 + ( x i b i 1 M i 1 b i 1 ) β ] ; x i [ M i 1 , b i 1 ]
The above equation satisfies the conditions that when x i = a i 1 , μ A ( a i 1 ) = 1 , and when x i = b i 1 , μ A ( b i 1 ) = 0.5 .
When x i belongs to adjacent level interval [ a i 2 , b i 2 ] , x i must be located at the right of M i 1 . The relative membership degree of x i to level 1 can be calculated as follows:
μ A ( x i ) 1 = 0.5 [ 1 + ( x i b i 1 d i 1 b i 1 ) β ] ; x i [ b i 1 , b i 1 ]
The above equation satisfies the conditions that when x i = b i 1 , μ A ( b i 1 ) = 0.5 , and when x i = d i 1 , μ A ( d i 1 ) = 0 where β is a non-negative number. When β = 1 , Equations (11) and (12) become linear functions.
(2) When h = c , then its criteria interval of the index I a b = [ a i c , b i c ] . Obviously, for A   , the relative membership of the upper bound b i c , μ A   ( b i c ) equals 1, and the relative membership degree of the upper bound a i 1 , μ A   ( a i c ) equals 0.5.
Suppose that M i c is a point value belonging to the interval [ a i c , b i c ] , and its relative membership degree to A   is 1. Since level c is the superior level, and the relative membership degree of the left bound point b i c to A   is 1, then M i c = b i c ; that is, M i c is located at the right bound point of [ a i c , b i c ] .
When x i belongs to the interval [ a i c , b i c ] , x i must be located at the left of M i c . The relative membership degree of x i to level c can be calculated as follows:
μ A C ( x i ) c = 0.5 [ 1 + ( x i b i c M i c b i c ) β ] ; x i [ M i c , b i c ]
The above equation satisfies the conditions that when x i = b i c , μ A C ( b i c ) = 1 , and when x i = a i c , μ A C ( a i c ) = 0.5 .
When x i belongs to adjacent level interval [ a i ( c 1 ) , b i ( c 1 ) ] , x i must be located at the left of M i c . The relative membership degree of x i to level c can be calculated as follows:
μ A C ( x i ) c = 0.5 [ 1 + ( x i a i 1 c i c a i c ) β ] ; x i [ c i c , a i c ]
(3) When c is the middle level, suppose that c is an odd number, and the middle level l = ( c + 1 ) / 2 means not good and not bad. For the criteria interval of the index [ a i l , b i l ] , its relative membership degree of the upper bound a i c and lower bound b i c , μ A ( a i c ) equals 0.5 and μ A C ( b i c ) equals 0.5, respectively. However, they are not the same as the relationship with the adjacent level, viz., a i l and b i ( l 1 ) are coincident, and b i l and a i ( l + 1 ) are coincident. The relative membership degree of the former is less than 0.5 to level ( l 1 ) , and the relative membership degree of the latter is less than 0.5 to level ( l + 1 ) .
Suppose that M i l is a point value of μ A ( u ) = 1 in the interval [ a i l , b i l ] , and M i l can be selected as the midpoint value of the interval [ a i l , b i l ] because the change from the superior to inferior levels is gradual; that is, M i l = ( a i l , b i l ) / 2 .
If x i belongs to the interval [ a i l , b i l ] , and x i is located at the left of M i l , the relative membership degree of x i to level l can be calculated as follows:
μ A ( x i ) l = 0.5 [ 1 + ( x i a i l M i l a i l ) β ] ; x i [ a i l , M i l ]
If x i belongs to the left adjacent level interval [ a i ( l 1 ) , b i ( l 1 ) ] , x i must also be located at the left of M i c . The relative membership degree of x i to level l can be calculated as follows:
μ A ( x i ) l = 0.5 [ 1 ( x i a i l c i l a i l ) β ] ; x i [ c i l , a i l ]
If x i belongs to the interval [ a i l , b i l ] , and x i is located at the right of M i l , the relative membership degree of x i to level l can be calculated as follows:
μ A ( x i ) l = 0.5 [ 1 + ( x i b i l M i l b i l ) β ] ; x i [ m i l , b i l ]
If x i belongs to the right adjacent level interval [ a i ( l + 1 ) , b i ( l + 1 ) ] , x i must also be located at the right of M i c . The relative membership degree of x i to level l can be calculated as follows:
μ A ( x i ) l = 0.5 [ 1 ( x i b i l d i l b i l ) β ] ; x i [ b i l , d i l ]
The M i h is an important parameter and can be obtained according to the standard eigenvalue a i h and b i h as follows:
M i h = c h c 1 a i h + h 1 c 1 b i h
Equation (19) satisfies the above three suppositions: (1) when h = 1 , then M i 1 = a i 1 ; (2)when h = c , then M i c = b i c ; (3) when h = l = + 1 2 , then M i l = a i l + b i l 2 .

3.2.3. Variable Fuzzy Evaluation Model (VFEM) for Classification of Stable Isotopes

Assume that there are n samples of stable isotopes which need to be assessed as follows:
X = { x 1 , x 2 , , x n }
The eigenvector of sample j is denoted by means of m indices eigenvalue as follows:
x j = ( x 1 j , x 2 j , , x m j ) T
Then, the matrix of eigenvalues to all samples is obtained as follows:
X = { x i j }
where x i j is the eigenvalue of sample j of the j th index, i = 1 , 2 , , m and j = 1 , 2 , , n . The evaluated samples are recognized according to c grades, so the standard eigenvalue matrix can be written as
Y = { y i h }
where y i h is the standard eigenvalue of index i about the level h , and h = 1 , 2 , , c .
According to the standard eigenvalue matrix of indices and practical situation, the matrix of attraction range I a b , the matrix of range I c d and the point value matrix M can be computed as follows:
I a b = ( [ a i h , b i h ] ) ,   I c d = ( [ c i h , d i h ] )
M = ( M i h )
where [ a i h , b i h ] and [ c i h , d i h ] are the attraction scope and the range of M i h is a point belonging to [ a i h , b i h ] , for which its relative membership degree equals 1.
The matrix of the relative membership degrees can be computed by Equations (11)–(18).
[ U h ] = ( μ A ( x i j ) h )
where μ A ( x i j ) h is the relative membership degree of the sample j to the index i about the grade h .
Assume that the weight vector of the indices is:
ω = ( ω 1 , ω 2 , , ω m )
and satisfies i = 1 m ω i = 1 , which can be obtained according to the order consistency theorem of the importance of deterministic index.
A synthetic relative membership degree vector can be obtained by the variable fuzzy evaluation model.
i u = 1 / [ 1 + j = 1 m [ w j ( 1 μ A ( x i j ) h ) ] p j = 1 m [ w j μ A ( x i j ) h ] p ] α / p
where α is the rule parameter of model optimization ( α = 1 for the least-single method and α = 2 for least squares method), and p is distance parameter ( p = 1 denotes hamming distance and p = 2 denotes Euclidean distance).
After normalizing i u , the normalized synthetic relative membership degree of each index about the grade h is obtained:
U = ( i u )
i μ h = i μ h / h = 1 c i μ h
The grade eigenvalues can be computed using Equation (27), and the evaluation level can be determined according to the value H .
H = ( 1 , 2 , , c ) U

3.2.4. Determination of the Weight of Index

In order to determine the weight vector of the two indices (δ2H and δ18O), the consistency scale matrix of the importance of two indices is obtained according to the order consistency theorem of the importance of deterministic index:
F = | 0.5 1 1 0.5 |
According to the ranking of the importance of matrix F and empirical knowledge, the index δ18O is in a “slightly” important position compared with δ2H. The non-normalized weight vector of the two evaluation indices can be obtained based on the relationship table between the mood operator and the relative membership degree [25]:
w = | 1 0.55 |
Further, the weights of δ2H and δ18O can be obtained by normalization:
w = | 0.65 0.35 |

4. Results

4.1. Classification of Stable Isotopes in the Naqu River Basin Based on the Variable Fuzzy Evaluation Method

Based on the comprehensive analysis of the influence factors of stable isotopic concentrations in water in the Naqu River basin, referring to the range of the water isotopes in Qinghai-Tibet plateau, the selected index system included δ18O and δ2H. Figure 2a,b show the variation range of δ2H and δ18O in water. The variation range of δ18O was −15.60‰ to 10.52‰, and the variation range of δ2H was −117.7‰ to −88.6‰. Figure 2a,b refer to the sampled waters. The stable isotopes values are shown in Table 1, which were from the sampled waters.
Two indices were divided into four levels according to the principle of proportion (Table 2). The higher the level was, the greater the concentration of heavy isotopes in the river after evaporation. This means that the river was under a non-equilibrium process. If the level was lower, it meant that the river was characterized by isotope dilution and replenished by precipitation or snowmelt.
According to the above analysis, matrices I a b , M , and I c d can be determined as follows
I a b = | [ 15.6 , 14.6 ] [ 14.6 , 13.6 ] [ 13.6 , 12.6 ] [ 12.6 , 10.52 ] [ 117.73 , 110.73 ] [ 110.73 , 103.73 ] [ 103.73 , 96.73 ] [ 96.73 , 88.62 ] |
I c d = | [ 15.6 , 13.6 ] [ 15.6 , 12.6 ] [ 14.6 , 10.52 ] [ 13.6 , 10.52 ] [ 117.73 , 103.73 ] [ 117.73 , 96.73 ] [ 110.73 , 88.62 ] [ 103.73 , 88.62 ] |
M = | 15.6 14.6 13.6 10.52 117.73 110.73 103.73 88.62 |
In the first place, it should be judged whether the evaluating index x i j is located at the left or right of point M i h . Then, the relative membership degree can be computed by using Equations (11)–(18).
U 1 = | 0.84 0.63 0.92 0.72 0 0 0.59 0 1 0.69 0.87 0.63 0.98 0.84 0 0.23 0.93 0 1 0.48 |
U 2 = | 0.16 0.38 0.09 0.28 0 0.30 0.42 0 0 0.32 0.13 0.37 0.02 0.16 0 0.73 0.07 0.21 0 0.98 |
U 3 = | 0 0 0 0 0 0.71 0 0.43 0 0 0 0 0 0 0 0.27 0 0.71 0 0 |
U 4 = | 0 0 0 0 0 0 0 0 1 1 0.21 0 0 0 0.35 0.29 0 0 0 0 |
When α = 1 and p = 2, according to Equations (28) and (29), the synthetically relative membership degree matrix U can be determined as follows:
U = | 0.970 0.736 0.995 0.911 0 0.007 0.849 0 1 0.545 0.030 0.264 0.005 0.089 0 0.385 0.151 0.008 0 0.455 0 0 0 0 0 0.585 0 0.735 0 0 0 0 0 0 1 0.023 0 0.256 0 0 |
The eigenvector of level H was computed according to Equation (30)
H = ( 1 2 3 4 ) × | 0.970 0.736 0.995 0.911 0 0.007 0.849 0 1 0.545 0.030 0.264 0.005 0.089 0 0.385 0.151 0.008 0 0.455 0 0 0 0 0 0.585 0 0.735 0 0 0 0 0 0 1 0.023 0 0.256 0 0 |
= ( 1.089 1.030 1.264 1.005 4 2.624 1.151 3.248 1 1.455 )
The grade eigenvalues of the stable isotopes in the river were obtained according to VFEM, as shown in Table 3. The grade eigenvalues in Table 3 based on Equation (30) were from sampled river waters in the summer. According to VFEM, this paper realized the quantitative classification of stable isotopes in different river sections. The status value in different river sections was ranked from low to high as follows: 1, 1.005, 1.089, 1.151, 1.264, 1.455 and 2.624.

4.2. Local Meteoric Water Line (LMWL) in the Naqu River Basin

The δ2H and δ18O in precipitation are currently measured at over 300 stations across the globe as part of the Global Net-work for Isotopes in Precipitation (GNIP), which cooperates with the International Atomic Energy Agency (IAEA) and the World Meteorological Organization (WMO) [29]. Despite the complexity of global hydrological processes, Craig [30] found that stable isotope ratios of hydrogen and oxygen in monthly precipitation correlate on a global scale in a linear relationship known as the global meteoric water line (GMWL):
δ2H = 8δ18O + 10
The best-fit line of this regression is the local meteoric water line (LMWL), which can be compared to the GMWL [31]. The isotopic characteristics of precipitation can be shown by the regression of the LMWL [32]. However, there is no suitable station having at least one full-year record in the Naqu River basin. So we analyzed the data for Lhasa to replace the LMWL because Lhasa and the Naqu River basin are located near each other in the middle of the Qinghai–Tibet Plateau, with a distance of approximately 300 km. Furthermore, the water vapor source of these two places is similar. The Lhasa LMWL can be described by the following equation [33]:
δ2H = 7.9δ18O + 6.29‰
The stable isotope values of precipitation are mainly determined by air temperature and moisture. An LMWL can reveal the evaporation capacity in the localities. The slope and intercept of LMWL (7.9,6.29) are both smaller than the GMWL (8,10) which shows the stable isotopic characteristics of precipitation in an arid region. The heavy isotope is enriched because of evaporation during rainfall in an arid climate. The slope of the regression line corresponds to the GMWL, indicating that condensation and precipitation occur at full equilibrium between the vapor and precipitation phases [34]. In contrast to the GMWL, the slope is slightly lower, showing drier and stronger evaporation conditions. The intercept is higher, likely due to a higher deuterium d-excess [35,36].
Deuterium excess is defined as:
d = δ2H − 8δ18O
It is used to measure the deviation of a set of data points from a line with a slope of 8 in δ2H vs. 8δ18O plots [37]. The d-excess value of GMWL is 10. The deuterium excess varies with humidity, wind speed and sea surface temperature during primary evaporation and averages to 10‰ on a global scale [38].

4.3. Stable Isotopic Composition of Surface Water, Soil Water and Snowfall

The stable isotopic composition of the the Naqu River is highly variable and is determined by the flow of water generated by heavy rainfall, soil water, and snowmelt events. These numbers came from the average of the sampled waters, which were from sampled river waters, soil water and snowfall. As shown in Table 4, in August, the δ18O values of the river waters varied from −15.9‰ to −10.5‰, and the δ2H values of the river waters varied from −119‰ to −89‰. The composition of heaver isotopes did not follow the LMWL.
The river water line is defined as:
δ2H = 5.5964δ18O − 30.942‰ (R2 = 0.94)
The regression line for the soil water line is:
δ2H = 7.6452 δ18O − 8.8688‰ (R2 = 0.91)
The regression line for the snow line is:
δ2H = 4.9117 δ18O − 24.915‰ (R2 = 0.91)
Figure 3 is related to the sampled waters. As shown in Figure 3, the slope and intercept of both the river water and soil water lines were lower than those of GMWL (8,10) and the LMWL (7.9,6.29) in August. The average deuterium excess of surface water was 3.17‰, which was smaller than the GMWL (10‰), implying a slower evaporation process in August. The slope of soil water regression was 7.65, which was larger than the Naqu River water line (5.60). The evaporation of soil water was not obvious. The slope (4.9) and intercept (−24.9) of the snow line were lower than the slope of GMWL (8,10) and the LMWL (7.9,6.29) in the winter.
Compared to the LMWL, three sets of isotopic data were distinguished. The first group included data points with low δ18O values and plots lay above the LMWL, which was from snowmelt. The second group, with medium δ18O values, fit the LMWL, which was from river water, which means that the water originated from local rainfall. The third group with high δ18O values, was below the LMWL, which was from soil water, signifying the effect of intensive evaporation processes.

5. Discussion

5.1. Analysis of Stable Isotopic Composition of Surface Water, Soil Water and Snowfall

The stable isotopic composition of the Naqu River is highly variable. There is no significant difference in δ18O and δ2H values between river waters and soil waters. However, there is a significant difference in δ18O and δ2H values between river waters and snowfall during the cold months, of which the mean value is −11.4‰ for δ18O and −81‰ for δ2H. The characteristic analysis of the snow samples showed the enrichment of stable isotopes. We assumed that the snow samples collected were the condensation of water vapor generated by continental water vapor. In the winter, it has less precipitation and is relatively depleted because of the lower temperature and water vapor concentration. During periods of low flow in an arid environment, evaporation produces characteristic heavy-isotope enrichment in residual surface waters.
In arid or semi-arid climates, the stable isotopic composition of surface waters can be modified by non-equilibrium evaporation, leading the more highly-evaporated river waters to deviate substantially from the LMWL. In this non-equilibrium process, the boundary layer effects induce the kinetic fraction of hydrogen and oxygen isotopes, which acts unequally on both isotopes. After the evaporation of waters, the heavier isotopes (18O and 2H) in the remaining water are enriched. As shown in Figure 3, the slope and intercept of both the river water and soil water lines were lower than those of the GMWL (8,10) and LMWL (7.9,6.29) in August, indicating that precipitation after evaporation was the main source of river water in the study area. The average deuterium excess of surface water was 3.17‰, which was smaller than the GMWL (10‰) implying a slower evaporation process in August. The intersection of the Naqu River water line and soil water line indicates that the surface water of the Mugequ River is replenished by the soil water. The slope of soil water regression was 7.65, which was larger than the Naqu River water line (5.60), indicating that surface water evaporated after soil water infiltration. The evaporation of soil water was not obvious. We assumed that soil water was replenished by rainfall. The slope (4.9) and intercept (−24.9) of the snow line were lower than the slope of the GMWL (8,10) and LMWL (7.9,6.29) in the winter. We assumed that the snow samples collected were the condensation of water vapor generated by continental water vapor, and that the river was only replenished by soil water in the winter. After the monsoon retreat, the Naqu River basin is mainly affected by continental water vapor, the atmospheric humidity is reduced, and the heavy isotope values of atmospheric water vapor are relatively high.
The study in the Naqu River basin showed the significant difference in mean values of both δ2H and δ18O along the river [38]. Therefore, the slope and intercept of the LMWL were lower than those of the refined GMWL. This deviation from the GMWL was mainly due to the different source of water replenishment with evaporation-induced kinetic fractionation effects [29]. Evaporation significantly affected both δ2H and δ18O values of flowing water and soil water, which were recharged by early snowmelt in the summer.
Compared to the LMWL, three sets of isotopic data are distinguished. Data points above the LMWL represent the surface water that mainly originated from snowmelt in mountains with low temperatures and absolute moisture content of air. Because the water exiting from the snowpack first during snowmelt was highly depleted isotopically, most of the river water plots above the LMWL were from the early melting water of snowpack and had depleted isotope values owing to the fractionation effect during the melting process [39]. The data fitting the LMWL means that the water was supplemented by local rainfall. The last group, which was from Tsonag Lake water and the main stream, was below the LMWL, signifying the effect of intensive evaporation processes (Figure 3).
Figure 3 shows the relationship between δ2H and δ18O values of precipitation in the Naqu River basin. Referring to the LMWL, there was a higher evaporation process in the main stream and the Mugequ River. The intersection of the Naqu River water line and soil water line indicateds that the surface water of Mugequ River is replenished by the soil water. δ2H and δ18O values between the Mugequ River and soil water were at the same level, which indicates the resupply between river and soil water. The Najinqu, Bazongqu, Mumuqu, Gongqu and Zongqingqu Rivers showed a lower level of heavy isotopes, which was the opposite of soil water, signifying the recharge of precipitation or snowmelt.

5.2. Estimation in Classification of Stable Isotopes in the Naqu River Basin Based on the Variable Fuzzy Evaluation Method

As shown in the Figure 4, it could be seen that the status values of the Najinqu, Bazongqu, Mumuqu, Chengqu and Gongqu Rivers were small, indicating that these tributaries were strongly supplemented by precipitation and snowmelt. The level of status value in the Mugequ River was 2.62 with a high level. In the lower reaches of the lake, the stable isotopes in the river were more abundant due to the evaporation of the lake and the role of stable isotopic fractionation.
The present paper developed a framework for analyzing stable isotopes according to VFEM, and generalized the source identification of water replenishment by rainfall or snowmelt in the Naqu River basin. This approach combined δ2H and δ18O together for the comprehensive evaluation of stable isotopes. Thus, the runoff replenishment by rainfall or snowmelt in the river can be visually determined by the status value of stable isotopes.
The contribution of this work is to transform double indices into the status value for the classification of stable isotopes in the source identification of water replenishment. For example, δ2H and δ18O were used for evaluating the source of water [12,13,14] by using the change in stable isotopic composition due to evaporation or mixing by the different waters, which mainly focused on the exploration of the characteristics of water isotopes; Zhang et al. [21] investigated the variation of δ2H and δ18O in glacial meltwater at Mount Everest and found that the δ2H in ice and snow was larger than in river and lake water relative to the change of δ18O, and δ2H and δ18O studied independently. This study mainly highlighted the discernible speed of runoff replenishment by rainfall or snowmelt based on the combination of the classification of stable isotopes and plot of δ2H vs. δ18O in the Naqu River basin.
However, according to VFEM and the results adopted in this paper, the rationality of the demarcation of the index standard value needs to be further improved. Furthermore, the data of stable isotope samples of river water, soil water and precipitation might need to be enlarged.

6. Conclusions

The classification of stable isotopes was obtained based on VFEM in the Naqu River basin. The grade eigenvalues in different river sections were ranked from low to high as follows: 1, 1.005, 1.089, 1.151, 1.264, 1.455 and 2.624. Compared to the LMWL, three sets of isotopic data were distinguished based on the plots of δ2H vs. δ18O because of different sources of water replenishment in the Naqu River basin. In the winter, snow is produced by the condensation of water vapor generated by continental water vapor, and the river is only replenished by soil water in the winter. In the summer, the status values of the Najinqu, Bazongqu, Mumuqu, Chengqu and Gongqu Rivers were small, indicating that these tributaries were strongly supplemented by precipitation and snowmelt. The level of status value in the Mugequ River was 2.62 with a high level. In the lower reaches of the lake, the stable isotopes in the river were more abundant due to the evaporation of the lake and the role of stable isotopic fractionation.
The present paper developed a framework of VFEM to provide a new method for the classification of stable isotopes. This research will be helpful for the study of runoff replenishment and the evolution mechanism in the Qinghai-Tibet Plateau. Furthermore, long-term and synchronous analysis of precipitations and runoff should be evaluated to enable a better understanding of runoff replenishment in the Naqu River basin.

Author Contributions

X.C. and F.W. conceived, designed, and drafted the manuscript; X.C. and G.W. planned and designed the methodology; X.C. revised the manuscript; F.W. guided and supervised the whole process; and all authors read and approved the final manuscript.

Funding

This work was supported by the Major Research Plan of the National Natural Science Foundation of China (91547209), National Key Research and Development Program of China (2016YFC0401401), the National Natural Science Foundation of People’s Republic of China (51579101,51709111), and the Distinguished Young Scholar of Science and Technology Innovation (184100510014).

Acknowledgments

The authors would like to express their sincere gratitude to the anonymous reviewers for their constructive comments and useful suggestions that helped us improve our paper.

Conflicts of Interest

No conflict of interest exists in the submission of this manuscript, and the manuscript is approved by all authors for publication.

References

  1. Edwards, T.W.D.; Bursey, G.G.; Prowse, T.D. Estimating Evaporation Using Stable Isotopes: Quantitative Results and Sensitivity Analysis for Two Catchments in Northern Canada. Hydrol. Res. 1993, 24, 79–94. [Google Scholar]
  2. Buttle, J.M. Fundamentals of Small Catchment Hydrology. In Isotope Tracers in Catchment Hydrology; Elsevier Science: Amsterdam, The Netherlands, 1998; Chapter 1; pp. 1–49. [Google Scholar]
  3. Barbieri, M.; Boschetti, T.; Petitta, M.; Tallini, M. Stable isotope (2H, 18O and 87Sr/86Sr) and hydrochemistry monitoring for groundwater hydrodynamics analysis in a karst aquifer (Gran Sasso, Central Italy). Appl. Geochem. 2005, 20, 2063–2081. [Google Scholar] [CrossRef]
  4. Clark, I.D.; Fritz, P. Environmental Isotopes in Hydrogeology; Lewis Publishers: Boca Raton, FL, USA, 1997; pp. 148–154, 172–186. [Google Scholar]
  5. Antúnez Contreras, S.P. Hydrogeological conceptual model of groundwater from carbonate aquifers using environmental isotopes (18O, 2H) and chemical tracers: A case study in southern Latium region, central Italy. J. Water Resour. Prot. 2012, 4, 695–716. [Google Scholar]
  6. Calligaris, C.; Mezga, K.; Slejko, F.; Urbanc, J.; Zini, L. Groundwater Characterization by Means of Conservative (δ18O and δ2H) and Non-Conservative (87Sr/86Sr) Isotopic Values: The Classical Karst Region Aquifer Case (Italy–Slovenia). Geosciences 2018, 8, 321. [Google Scholar] [CrossRef]
  7. Rank, D.; Wyhlidal, S.; Schott, K.; Weigand, S.; Oblin, A. Temporal and spatial distribution of isotopes in river water in Central Europe: 50 years experience with the Austrian network of isotopes in rivers. Isot. Environ. Health Stud. 2018, 54, 115–136. [Google Scholar] [CrossRef] [PubMed]
  8. Gibson, J.J.; Edwards, T.W.D.; Birks, S.J.; St Amour, N.A.; Buhay, W.M.; McEachern, P.; Wolfe, B.B.; Peters, D.L. Progress in isotope tracer hydrology in Canada. Hydrol. Process. 2005, 19, 303–327. [Google Scholar] [CrossRef]
  9. Wu, H.; Zhang, X.; Li, X.; Li, G.; Huang, Y. Seasonal variations of deuterium and oxygen-18 isotopes and their response to moisture source for precipitation events in the subtropical monsoon region. Hydrol. Process. 2015, 29, 90–102. [Google Scholar] [CrossRef]
  10. Han, G.; Lv, P.; Tang, Y.; Song, Z. Spatial and temporal variation of H and O isotopic compositions of the Xijiang River system, Southwest China. Isot. Environ. Health Stud. 2018, 54, 137–146. [Google Scholar] [CrossRef]
  11. Qu, S.; Wang, Y.; Zhou, M.; Liu, H.; Shi, P.; Yu, Z.; Xiang, L. Temporal 18O and deuterium variations in hydrologic components of a small watershed during a typhoon event. Isot. Environ. Health Stud. 2017, 53, 172–183. [Google Scholar] [CrossRef]
  12. Meredith, K.T.; Hollins, S.E.; Hughes, C.E.; Cendón, D.I.; Chisari, R.; Griffiths, A.; Crawford, J. Evaporation and concentration gradients created by episodic river recharge in a semi-arid zone aquifer: Insights from Cl, δ18O, δ2H, and 3H. J. Hydrol. 2015, 529, 1070–1078. [Google Scholar] [CrossRef]
  13. Biggs, T.W.; Lai, C.; Chandan, P.; Lee, R.M.; Messina, A.; Lesher, R.S.; Khatoon, N. Evaporative fractions and elevation effects on stable isotopes of high elevation lakes and streams in arid western Himalaya. J. Hydrol. 2015, 522, 239–249. [Google Scholar] [CrossRef]
  14. Yang, Y.; Xiao, H.; Qin, Z.; Zou, S. Hydrogen and oxygen isotopic records in monthly scales variations of hydrological characteristics in the different landscape zones of alpine cold regions. J. Hydrol. 2013, 499, 124–131. [Google Scholar] [CrossRef]
  15. Edirisinghe, E.A.; Pitawala, H.M.; Dharmagunawardhane, H.A.; Wijayawardane, R.L. Spatial and temporal variation in the stable isotope composition (δ18O and δ2H) of rain across the tropical island of Sri Lanka. Isot. Environ. Health Stud. 2017, 53, 1–18. [Google Scholar] [CrossRef] [PubMed]
  16. Sanchez-Murillo, R.; Esquivel-Hernandez, G.; Saenz-Rosales, O.; Piedra-Marin, G.; Fonseca-Sanchez, A.; Madrigal-Solis, H.; Ulloa-Chaverri, F.; Rojas-Jimenez, L.D.; Vargas-Viquez, J.A. Isotopic composition in precipitation and groundwater in the northern mountainous region of the Central Valley of Costa Rica. Isot. Environ. Health Stud. 2017, 53, 1–17. [Google Scholar] [CrossRef]
  17. Jin, L.; Siegel, D.I.; Lautz, L.K.; Lu, Z. Identifying streamflow sources during spring snowmelt using water chemistry and isotopic composition in semi-arid mountain streams. J. Hydrol. 2012, 470–471, 289–301. [Google Scholar] [CrossRef]
  18. Wassenaar, L.I.; Athanasopoulos, P.; Hendry, M.J. Isotope hydrology of precipitation, surface and ground waters in the Okanagan Valley, British Columbia, Canada. J. Hydrol. 2011, 411, 37–48. [Google Scholar] [CrossRef]
  19. Rock, L.; Mayer, B. Isotope hydrology of the Oldman River basin, southern Alberta, Canada. Hydrol. Process. 2007, 21, 3301–3315. [Google Scholar] [CrossRef] [Green Version]
  20. Yao, Z.; Liu, J.; Huang, H.; Song, X.; Dong, X.; Liu, X. Characteristics of isotope in precipitation, river water and lake water in the Manasarovar basin of Qinghai–Tibet Plateau. Environ. Geol. 2009, 57, 551–556. [Google Scholar] [CrossRef]
  21. Zhang, S.; Yu, W.X.; Zhang, Q.L. The distribution of deuterium and heavy oxygen in snow and ice, in the Jolmo Lungma regions of the southern Tibet. Sci. China B 1973, 4, 430–433. [Google Scholar]
  22. Zhang, X.P.; Zhong, W.Z.Y.; Yao, T.D.; Han, J.K.; Xie, Z.C. The stable isotope transformation of space-time in meteoric water in Qinghai–Tibet Plateau and its adjacent region. Sci. China D 2001, 31, 353–361. [Google Scholar]
  23. Tian, L.D.; Yao, T.D.; Shen, Y.P.; Yang, M.X.; Bai-Sheng, Y.E.; Tsujimura, M. Study on stable isotope in river water and precipitation in Naqu River basin, Tibetan Plateau. Adv. Water Sci. 2002, 13, 206–210. [Google Scholar]
  24. Pang, H.; He, Y.; Lu, A.; Zhao, J.; Ning, B.; Yuan, L.; Song, B.; Zhang, N. Comparisons of stable isotopic fractionation in winter and summer at Baishui Glacier No. 1, Mt. Yulong, China. J. Geogr. Sci. 2006, 16, 306–314. [Google Scholar]
  25. Chen, S.Y.; Hu, J.M. Variable fuzzy assessment method and its application in assessing water resources carrying capacity. J. Hydraul. Eng. 2006, 37, 264–277. [Google Scholar]
  26. Chen, S.Y. Theory and model of engineering variable fuzzy sets—Mathematical basis for fuzzy hydrology and water resources (in Chinese). J. Dalian Univ. Technol. 2005, 45, 308–312. [Google Scholar]
  27. Wang, W.; Xu, D.; Chau, K.; Lei, G. Assessment of river water quality based on theory of variable fuzzy sets and fuzzy binary comparison method. Water Resour. Manag. 2014, 28, 4183–4200. [Google Scholar] [CrossRef]
  28. Zou, Z.; Zhou, C.; Song, L.X.; Guo, J. Comprehensive flood risk assessment based on set pair analysis-variable fuzzy sets model and fuzzy AHP. Stoch. Environ. Res. Risk Assess. 2013, 27, 525–546. [Google Scholar] [CrossRef]
  29. Price, R.M.; Swart, P.K.; Willoughby, H.E. Seasonal and spatial variation in the stable isotopic composition (δ18O and δ2H) of precipitation in south Florida. J. Hydrol. 2008, 358, 193–205. [Google Scholar] [CrossRef]
  30. Craig, H. Isotopic variations in meteoric waters. Science 1961, 133, 1702–1703. [Google Scholar]
  31. Li, S.G.; Romero-Saltos, H.; Tsujimura, M.; Sugimoto, A.; Sasaki, L.; Davaa, G.; Oyunbaatar, D. Plant water sources in the cold semiarid ecosystem of the upper Kherlen River catchment in Mongolia: A stable isotope approach. J. Hydrol. 2007, 333, 109–117. [Google Scholar] [CrossRef]
  32. Šturm, M.B.; Ganbaatar, O.; Voigt, C.C.; Kaczensky, P. First field-based observations of δ2H and δ18O values of event-based precipitation, rivers and other water bodies in the Dzungarian Gobi, SW Mongolia. Isot. Environ. Health Stud. 2017, 53, 157–171. [Google Scholar]
  33. Tian, L.D.; Yao, T.D.; Sun, W.Z.; Stievenard, M.; Jouzel, J. Relationship between δ2H and δ18O and water vapour circulation in the precipitation of the Qinghai-Tibet plateau. Sci. Sin. Terrae 2001, 31, 214–220. [Google Scholar]
  34. Dansgaard, W. Stable isotopes in precipitation. Tellus 1964, 16, 436–468. [Google Scholar] [CrossRef]
  35. Clark, I.; Fritz, P. Environmental Isotopes in Hydrogeology; Lewis Publishers: Boca Raton, FL, USA, 1964; 328p. [Google Scholar]
  36. Wu, J.K.; Ding, Y.J.; Yang, J.H.; Liu, S.W.; Chen, J.Z.; Zhou, J.X.; Qin, X. Spatial variation of stable isotopes in different waters during melt season in the Laohugou Glacial Catchment, Shule River basin. J. Mt. Sci. 2016, 13, 1453–1463. [Google Scholar] [CrossRef]
  37. Jeelani, G.; Saravana Kumar, U.; Kumar, B. Variation of δ18O and δ2H in precipitation and stream waters across the Kashmir Himalaya (India) to distinguish and estimate the seasonal sources of stream flow. J. Hydrol. 2013, 481, 157–165. [Google Scholar] [CrossRef]
  38. Mayr, C.; Lücke, A.; Stichler, W.; Trimborn, P.; Ercolano, B.; Oliva, G.; Ohlendorf, C.; Soto, J.; Fey, M.; Haberzettl, T.; et al. Precipitation origin and evaporation of lakes in semi-arid Patagonia (Argentina) inferred from stable isotopes (δ18O, δ2H). J. Hydrol. 2007, 334, 53–63. [Google Scholar] [CrossRef]
  39. Wang, Y.; Chen, Y.; Li, W. Temporal and spatial variation of water stable isotopes (18O and 2H) in the Kaidu River basin, Northwestern China. Hydrol. Process. 2014, 28, 653–661. [Google Scholar] [CrossRef]
Figure 1. Research area and the sampling positions distribution.
Figure 1. Research area and the sampling positions distribution.
Water 11 00046 g001
Figure 2. Spatial variation of stable isotopic enrichment in the Naqu River basin.
Figure 2. Spatial variation of stable isotopic enrichment in the Naqu River basin.
Water 11 00046 g002
Figure 3. Plot of δ2H vs. δ18O for all water bodies in the Naqu River basin.
Figure 3. Plot of δ2H vs. δ18O for all water bodies in the Naqu River basin.
Water 11 00046 g003
Figure 4. Spatial variation of the status values in the Naqu River basin.
Figure 4. Spatial variation of the status values in the Naqu River basin.
Water 11 00046 g004
Table 1. The characteristic values of δ18O and δ2H in the rivers.
Table 1. The characteristic values of δ18O and δ2H in the rivers.
The Evaluation IndexEvaluation of Regional (‰)
Najinqu RiverSangqu RiverBazongqu RiverMumuqu RiverMiddle StreamMugequ RiverChengqu RiverThe DownstreamGongqu RiverZongqingqu River
δ18O−15.3−14.9−15.4−15.0−10.5−13.2−14.8−12.3−15.6−15.0
δ2H−116−113−117−115−89−107−117−100−118−110
Table 2. The standard value of stable isotopes in the river.
Table 2. The standard value of stable isotopes in the river.
The Evaluation IndexIndex Standard Value (‰)
Level 1Level 2Level 3Level 4
δ18O
δ2H
≤−14.6
≤−111
(−14.6) to (−13.6)
(−111) to (−104)
(−13.6) to (−12.6)
(−104) to (−97)
0 ≥ −12.6
≥−97
Table 3. The grade eigenvalues of the stable isotopes in the river.
Table 3. The grade eigenvalues of the stable isotopes in the river.
Evaluation of RegionalStatus ValueRating
Najinqu River1.0891
Sangqu River1.031
Bazongqu River1.2641
Mumuqu River1.0051
Middle stream44
Mugequ River2.6242
Chengqu River1.1511
The downstream3.2483
Gongqu River11
Zongqingqu River1.4551
Table 4. The statistical table of stable isotopic composition in the Naqu River basin.
Table 4. The statistical table of stable isotopic composition in the Naqu River basin.
TimeSamplesRange of Variation (‰)
δ18Oδ2H
AugustSurface water−15.6 to 10.5−118 to −89
JanuarySurface water−16.6 to −9.1−122 to −81
AugustSoil water−17.3 to −11.2−151 to −96
JanuarySnowfall−15.0 to −7.6−101 to −63

Share and Cite

MDPI and ACS Style

Chen, X.; Wang, G.; Wang, F. Classification of Stable Isotopes and Identification of Water Replenishment in the Naqu River Basin, Qinghai-Tibet Plateau. Water 2019, 11, 46. https://doi.org/10.3390/w11010046

AMA Style

Chen X, Wang G, Wang F. Classification of Stable Isotopes and Identification of Water Replenishment in the Naqu River Basin, Qinghai-Tibet Plateau. Water. 2019; 11(1):46. https://doi.org/10.3390/w11010046

Chicago/Turabian Style

Chen, Xi, Guoli Wang, and Fuqiang Wang. 2019. "Classification of Stable Isotopes and Identification of Water Replenishment in the Naqu River Basin, Qinghai-Tibet Plateau" Water 11, no. 1: 46. https://doi.org/10.3390/w11010046

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