Advances in Meteorology
Volume 2015 (2015), Article ID 830306, 10 pages
http://dx.doi.org/10.1155/2015/830306

Research Article

Hydrograph Separation in the Headwaters of the Shule River Basin: Combining Water Chemistry and Stable Isotopes

Jiaxin Zhou,1,2 Jinkui Wu,1,2 Shiwei Liu,1,2 Guoxiong Zeng,³ Jia Qin,1,2 Xiuna Wang,1,2 and Qiudong Zhao1,2

¹Laboratory of Watershed Hydrology and Ecology, Cold and Arid Regions Environmental and Engineering Research Institute, Chinese Academy of Sciences, Lanzhou 730000, China
²State Key Laboratory of Cryospheric Science, Cold and Arid Regions Environmental and Engineering Research Institute, Chinese Academy of Sciences, Lanzhou 730000, China
³Shule River Basin Water Resources Administration of Gansu Province, Yumen 735200, China

Received 25 December 2014; Revised 20 April 2015; Accepted 28 April 2015

Academic Editor: Fengjing Liu

Copyright © 2015 Jiaxin Zhou et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

The runoff components were identified in the headwater area of Shule River Basin, using isotopic and chemical tracing with particular focus on the temporal variations of catchment sources. A total of 95 samples, including precipitation, groundwater, and glacial meltwater, were collected and analyzed for stable water isotopes (¹⁸O and ²H) and major chemical ion parameters (potassium, sodium, calcium, magnesium, sulfate, chloride, and bicarbonate). Based on the isotope and water chemistry data, we applied end member mixing analysis (EMMA) to identify and quantify the major runoff generating sources and their contributions. The contributions of groundwater, precipitation, and glacial meltwater were 66.7%, 19.9%, and 13.4%, respectively. The study indicated that groundwater dominated runoff in the headwater area of Shule River Basin. The roles of glacier meltwater should be remarkable in water resource management in this basin. The uncertainties of the EMMA method were summarized and estimated via a classical Gaussian error propagation technique. Analyses suggested that the uncertainty in the measurement method was less important than that in the temporal and spatial variations of tracer concentrations. The uncertainty was sensitive when the difference between mixing components was small. Therefore, the variation of tracers and the difference of mixing components should be considered when hydrograph separation was applied in the basin.

1. Introduction

The quantification of catchment response to rainfall or snowmelt events in terms of water fluxes and chemical composition is an important issue in catchment hydrology. In particular, during flooding periods, different interacting processes occur that are spatially distributed within the catchment [1, 2]. These processes are defined by physiographic characteristics. In addition, runoff generation depends on the initial state of the various hydrological reservoirs and on the characteristics of the hydrological input (precipitation or snowmelt). Due to these factors, it is difficult to identify the dominant runoff generation processes [3].

The water shortage and low use efficiency make china thirsty, and the loss of glacier and wetland in the western plateau will exaggerate this thirst in the future [4], while the same situation happens in other places of the world [5]. Therefore, it is important to understand the runoff generation mechanism. In arid or semiarid regions, water is a key factor affecting the biomass production. A better understanding of runoff generation processes as well as catchment function is important for improved water resources management [6]. The hydrograph separation technique using natural tracers, in which different runoff components are quantified according to their chemical signature, is a widely used method for investigating runoff generation processes at the catchment scale [7]. Isotopes were used to quantify the interaction of different end member in glaciated catchments for longer time periods (monthly) and/or larger catchment areas (>1000 km²) [8–10]. Isotope techniques can be easily and successfully used to study the origin and dynamics of surface water and groundwater, evaporation of water bodies, and mixing processes between various water sources [11–13]. To obtain both temporal and spatial origins, some investigations using stable isotopes associated with chemical tracers have been undertaken in several different basins [14, 15]. Hydrochemical tracers, such as PH, electrical conductivity, or the concentration of different anions and cations [16], have been used to determine the origin of runoff components. In recent years, geochemical methods and environmental isotope techniques have been used increasingly to determine runoff components in various catchments under different environmental conditions [17, 18]. In the arid and semiarid areas, a combination of hydrologic and environmental isotope methods (¹⁸O, D) has been proved to be a valuable tool for studying processes within the water cycle and in isotope hydrology [19–21].

One common tool to identify runoff sources and flux components and calculate their contributions to the stream discharge is end member mixing analysis (EMMA) [22]. EMMA techniques have been applied in varieties of studies, at both the plot and the catchment scales. Only few studies have applied EMMA in the arid or semiarid regions [23]. It has been applied in many studies to identify end members at small catchment scales that describe the vertical sequence of water storages to flow contribution. These vertical end members are, for example, rain, soil water, and groundwater [24] or overland flow, soil water, and hillslope water. Fewer studies have applied EMMA at larger catchment scales of hundreds or even thousands of square kilometers [25]. This approach is based on three mass conservations: one for water, one for isotopic tracer, and one for geochemical tracer. It allows separating the relative contribution of the different components, which correspond to different reservoirs or different contributive areas. The use of isotopic tracers allows separating the runoff hydrograph into preevent water and event water while the use of geochemical tracers allows identifying the three origins of the runoff components [26, 27].

There is a clear need to develop predictive capabilities related to the identification of runoff generating sources in large ungauged basins, particularly in emergent countries such as China [28]. Runoff generation and dynamics is an important issue in watershed and water resource management. On the one hand, knowledge about runoff generation processes and flow pathways is crucial for evaluating the vulnerability of surface and groundwater system [29]. On the other hand, such knowledge helps to develop and validate hydrological models. Since arid and semiarid basins usually have more severe natural conditions and scarcities of observation data, the application of isotope techniques in catchment hydrology study seems to be a more economic and helpful tool [30]. Understanding hydrological processes will significantly add our ability to evaluate potential tradeoffs between social development and water availability. We expect that scientific results will provide an insight for water resource and watershed management in a large-area.

In this study, we applied the EMMA method to identify and quantify the major runoff generating sources in a three end members’ system. The objectives of this study are (1) to identify runoff producing sources using ¹⁸O and chloride ion as tracers in the headwater area of the Shule River Basin, (2) to investigate the applicability of the EMMA method in semiarid catchments, (3) to calculate the contributions of the three components of runoff. We developed a conceptual hydrograph separation technique, namely, three components’ mixing model. It is based on the steady-state mass balance equations of water and tracer fluxes in a catchment. In addition, uncertainties analyses were performed for the hydrograph separation.

2. Study Areas

The Shule River, the third largest inland rivers basin with the whole catchment area of approximately 14.21 × 10⁴ km², is located in the Western Qilian Mountains. The upstream of the Shule River (Figure 1), at 96.6°E~99.0°E and 38.2°N~40.0°N, with the area of 1.14 × 10⁴ km² and the mean elevation of 3885 m, is located in the Tianjun, Qinghai Province, Northeastern margin of the Tibetan Plateau. The mainstream annual runoff varying significantly during different years is 10.83 × 10⁸ m³. The 53% of the total annual runoff concentrated between July and September. However, the runoff distributed unevenly over time. The runoff in spring and winter only makes up for 8.5% and 10%, respectively. According to nearly 40 years’ hydrological data of Changmabao gauge station, the annual runoff of dry years and wet years is 5.36 × 10⁸ m³ and 15.07 × 10⁸ m³, respectively.

Figure 1: Site map showing the upstream of Shule River catchment and sampling sites.

In the headwater area of Shule River Basin, our study area, Gahe, at 96.49°E~98.58°E and 38.02°N~39.12°N, with the area 4 096 km², there are 347 glaciers and the area of glaciers is 29.45 km², which accounts for 0.72% of the headwater area (Figure 2). Glaciers are mainly distributed above elevation 4 500 m, which is located in Shule Nanshan and Tuolai Nanshan [31]. Annual sunshine time is 3 033–3 246 hours. The mean annual elevation is 4 000~4 500 m. The mean annual air temperature is approximately −5°C, the annual precipitation is 100–300 mm and mainly falling between May and September, and annual evaporation is about 1 200 mm [32]. The temperatures of the hottest month (July) and the coldest month (January) are 7.5°C and −17.5°C, respectively. The study area belongs to the continental arid desert climate region which is characterized by cold, dry winters and relative warm, wet summers [33]. In the growing season of May through September, the plentiful sunshine and rainfall (80% of annual total precipitation) allows plants to grow efficiently.

Figure 2: Glacier distribution of Gahe region.

The Quaternary sediments, comprising diluvial-alluvial, aeolian, and lacustrine deposits, form the main aquifers in the basin. These sediments are enriched in calcite, gypsum, and mirabilite in parts of the middle reaches, and soil salinization occurs widely in the middle and lower reaches. The depth to the water table is 5–10 m. There formed a large area of swamp in the source places. The aquifer system includes a thick unconfined zone consisting of coarse-grained gravel, sand and a confined part consisting of medium to fine and silty sand. The landscape is characterized by large mountain ranges with steep valleys and gorges interspersed with relatively level and wide intermountain grassland basins [34].

3. Material and Methods

3.1. Field Sampling

Intensive synoptic sampling was carried out between April and September 2009 in Gahe, the headwater area of the Shule River, consisting of 95 samples. Precipitation, glacial meltwater, groundwater, and river water were sampled. The number of four kinds of samples is precipitation 15, river water 30, groundwater 31, and glacial meltwater 19, respectively. Precipitation, glacial meltwater, groundwater, and river water were sampled and analyzed for stable water isotopes (¹⁸O and ²H), major ion chemistry parameters as well. Samples were collected in polyethylene bottles and filtered through 0.45 mm Millipore membrane for major element analyses. Meteorological parameters and hydrology data were measured continuously by means of an automatic weather station and gauge station.

Precipitation samples were collected immediately after each precipitation event in order to minimize the alteration of heavy isotopes by evaporation with plastic basin sets. River water and groundwater samples were collected once a week. Due to the limitations of some nature conditions, we cannot get to glaciers distributed around Gahe. Hence, we collected the glacial meltwater samples in Laohugou Glacier number 12. Due to the background of the same atmospheric circulation, the moisture of the two sites comes from the same source. We considered the substitute is feasible.

3.2. Laboratory Analyses

All samples were kept in near-frozen condition and transported to the State Key Laboratory of Cryospheric Science, Cold and Arid Regions Environmental and Engineering Research Institute, Chinese Academy of Sciences, for test.

Concentration of anions (Cl⁻, ) was analyzed by Ion Chromatography (IC, DX-120, Dionex, Germany) while and were analyzed by the titration method. Cations, K⁺, Na⁺, Ca²⁺, and Mg²⁺ were analyzed by using Atomic Absorption Spectroscopy (AAS) method. Every sample value represents the mean of two consecutive measurements. Measurement errors were less than 1%. The detection limits of all ions were lower than 0.1 mg/L. Chloride ion and ¹⁸O were finally selected to assess the different contributing sources using mass balance equations and end member mixing diagrams.

The δD and δ ¹⁸O composition of all water samples were analyzed by Liquid-Water Isotope Analyzer (DLT 100, Los Gatos, USA) based on off-axis integrated cavity output spectroscopy (OA-ICOS). Each sample is injected six times to avoid memory effect between samples. The isotopic ratios were expressed in per mil () units relative to Vienna Standard Mean Ocean Water (V-SMOW):

where is the ration ¹⁸O/¹⁶O or ²H/¹H. Precision of δD and δ ¹⁸O was ±0.6 and ±0.2, respectively.

3.3. Hydrograph Separation Method and Uncertainties Analysis

In general, hydrograph separations are based on the steady-state mass balance equations of water and tracer fluxes in a catchment [34]. Following are the equations:

We converted the equations to vector form:

Then the right side equations are divided by :

where is the discharge and and are the concentration of tracer chloride ion and ¹⁸O, respectively. Subscripts , and refer to river water, groundwater, precipitation, and glacial meltwater, respectively. The application of these equations is based on certain assumptions, which are discussed, for instance, by Hinton et al. [35], Buttle [1], or Rodhe [36]:

(1)there is a significant difference between the tracer concentrations of the different components;(2)the tracer concentrations are constant in space and time, or any variations can be accounted for;(3)contributions of an additional component must be negligible, or the tracer concentrations must be similar to that of another component;(4)the tracers must mix conservatively;(5)the tracer concentrations of the components are not collinear.

Recent focus of hydrograph separation has been on uncertainty analysis. Several approaches are available for calculating uncertainty. Genereux (1998) suggested a general uncertainty propagation technique using Gaussian error estimators for two- and three-component separations [37]. However, an extensive overview of all possible causes of hydrograph separation uncertainties during different periods of a given event is still lacking. A classical Gaussian error propagation technique was applied to quantify the uncertainty of tracer-based hydrograph separations. This technique is generally used in other scientific and engineering problems. Errors of all separation equation variables are considered. Assuming that the uncertainty in each variable is independent of the uncertainty in the others, the relative error of the contribution of a specific runoff component is related to the uncertainty in each of the variables by the following [37]:

where represents the uncertainty in the variable specified in the subscript, is the concentration of corresponding tracer, represents the event water, and represents the preevent water. In the results the relative error is given as percentage value.

It is demonstrated that large relative uncertainties must be considered for the quantification of runoff components. Uncertainties are caused by (1) tracer analysis and discharge measurement; (2) intrastorm variability of ¹⁸O; (3) elevation effect of ¹⁸O and chloride; (4) solution of minerals during runoff formation; and (5) general spatial heterogeneity of tracer concentrations. The last source of error was the most significant. An investigation on the dominating runoff generation processes in the catchment, before a model is set up, would reduce such uncertainties.

4. Results

4.1. Temporal Variance of Runoff

The temporal variance of runoff is showed in Figure 3. The runoff showed a significant seasonal variation. The runoff varied in the range of 3.82~338.09 m³/s with an average of 71.57 m³/s. There is a minor peak in April since the snowmelt peak usually occurs in spring. We could see from the figure that most of peak flows were corresponding with the big rainfall events between June and September. It means the significant increase of runoff is the results of precipitation event. The average runoff is dominated by a snowmelt peak in spring followed by a decline in discharge over the growing season. From June to September, when most rainstorms occur, there is a considerable increase in discharge, followed by an again declining hydrograph until October, when the river itself freezes.

Figure 3: Runoff and rainfall from April to October 2009.

4.2. Isotopic Composition

4.2.1. Isotopic Composition of River Water

The isotopic composition of river water during April and September in the Gahe station shows a steady variability, ranging from −9.9 to −8.5 in δ ¹⁸O and from −68.9 to −58.2 in δD, respectively. The possibility of differential isotopic evaporation of samples can be analyzed by comparing the samples to the meteoric water line which is formed by plotting δD against δ ¹⁸O [38]. The local meteoric water line (LMWL) is commonly used as indicators of water vapor source, source of the humidity, and kinetic conditions in a number of fields including isotope hydrology [39]. The relationship between river water and the local meteoric water line (LMWL) was displayed by bivariate plot of δ ¹⁸O versus δD (Figure 4). According to the distribution of river water in the space of δ ¹⁸O versus δD, most of the river water sample points were located approaching the local meteoric water line (LMWL) (δD = 8.11δ ¹⁸O + 11.40, , ). Also, the slope of the regression line was fairly close to the multiple-year observed values in Northwest China (7.88) and in Heihe River Basin, an inland river basin neighboring the study area (7.82) [40].

Figure 4: Stable isotope (δD and δ ¹⁸O) compositions of river water.

4.2.2. Isotopic Signature of Precipitation, Groundwater, and Glacial Meltwater

The isotopic composition of precipitation shows a relative significant variance. The values of δ ¹⁸O fluctuate in −13.0~−8.3 and −95.1~−54.2 in δD. The equation between δ ¹⁸O and δD (δD = 7.68δ ¹⁸O + 9.29, , ) (Figure 5(a)). The temporal and spatial variability of ¹⁸O in precipitation are relatively high. This is caused by fractionation during evapotranspiration and condensation due to lower saturated vapor pressure of water molecules containing the heavier ¹⁸O isotope than that of water molecules containing the lighter ¹⁶O isotope. As a result, the δ ¹⁸O in precipitation decreases with decreasing air temperature, increasing elevation, increasing latitude, increasing distance of vapor transport through the atmosphere, and increasing precipitation amounts.

(a)

(b)

(c)

Figure 5: Stable isotope (δD and δ ¹⁸O) compositions of precipitation (a), groundwater (b), and glacial meltwater (c).

The isotopic composition of groundwater ranges in −9.5~−6.7 (in δ ¹⁸O) and −68.2~−45.3 (in δD). The equation between δ ¹⁸O and δD (δD = 8.28δ ¹⁸O + 10.74, , ) (Figure 5(b)). It is fairly close to that of river water. The value of δ ¹⁸O ranges in −14.7~−12.3 and δD in −105.0~−85.1 in glacial meltwater. The equation between δ ¹⁸O and δD (δD = 7.55δ ¹⁸O + 7.34, , ) (Figure 5(c)). The stable isotope ratios of hydrogen and oxygen of water samples can provide essential information about water dynamics within a given watershed. In general, this is from isotope fractionation by evaporation, altitude effects, and different water sources they received [23]. The slope and the intercept of LMWL were slightly lower, showing drier and stronger local evaporation conditions. Evaporation caused a differential increase in the δD and δ ¹⁸O values of the remaining water, resulting in a lower slope for the linear relationship between δD and δ ¹⁸O values [41].

4.3. Temporal Variance of Cl⁻

It can be assumed that mixing processes in the catchment determine the isotopic concentration of total runoff. However, the hydrochemical composition of water is essentially changed as a result of interactions with organic and inorganic material during its passage through the unsaturated and saturated zones. The concentration of Cl⁻ in river water fluctuates in 9.4~13.4 mg·L⁻¹, with an average of 11.2 mg·L⁻¹. The variance of Cl⁻ concentration has much relationship with runoff (Figure 6). In spring, the spring flood caused by snowmelt water makes the soil chemical ions into the river, so the concentration of Cl⁻ is relatively high. Afterwards, with the increase of snowmelt water, glacier meltwater, and precipitation, the runoff has been showing a different amplitudes increase. With the increase of runoff, the dilution effect of ions has also increased, so the concentration of Cl⁻ decreased. Although groundwater recharged, the dilution effect outweighs the supply effect. During August and September, the runoff has a considerable decrease; the concentration of Cl⁻ decreased because of the weakening of dilution effect.

Figure 6: Monthly runoff and concentration of Cl⁻ (with error bars).

4.4. Identification of End Members

Applying the method of geochemical and isotopic tracing, in this paper the runoff characteristics and hydrological law at the Gahe station are investigated at different periods in 2009. We conducted a performance of principal component analysis (PCA) on the concentration data of chloride ion. The results showed that the chemical tracer exhibits conservative behavior, whereas isotopes are geographical source tracers and only change composition due to slow fractionation processes [42]. ¹⁸O belongs to the group of stable environmental isotopes occurring naturally in water. It has been widely used to separate storm flow into proportions of event and preevent water [1]. As part of the water molecule, ¹⁸O behaves conservatively; that is, the combination of chemical and isotopic tracers allows identifying the origin of water pathways.

For the three-component hydrograph separation, the choice of a suitable tracer constellation to explain the chemical changes in discharge during a storm as well as to determine and to identify dominant sources, flow paths, and residence times in the catchment becomes increasingly important. The study suggests that hydrological tracer chlorine concentration and ¹⁸O can be used under certain hydrological and lithological conditions. A system of algebraic equations is introduced that enables a three-component hydrograph separation by using ¹⁸O and chlorine. These alternative tracers should, however, be verified against more conventional tracers before use, as the behavior depends on specific characteristics of solutes.

The basic assumption in EMMA is that the stream water is a discrete mixture of its sources. The sources must therefore be of sufficiently different concentrations compared with the stream water. We projected the average values of tracer chloride ion and ¹⁸O of three end members in triangle to test and verify the independence (Figure 7). It shows that most of stream water observations fall into the triangle that is spanned by three end members (precipitation, groundwater, and glacial meltwater). However, there exist some stream water observations that lie outside of the triangle. In many other studies that apply EMMA such a situation has been described [14, 22–24, 35]. These outliers result from a number of factors including (1) uncertainty in field sampling and laboratory analyses, (2) lack of temporal invariance of end members, or (3) the expression of different end member in the mixture as water source areas change temporally. Overall, the result can lead to over- or underprediction of the contributions of each end member to the stream water and should be understood as a source of uncertainty.

Figure 7: Stream water observations and the average values of tracers Cl⁻ and ¹⁸O of three end members that spanned the triangle.

4.5. Contribution of End Members to Runoff

Owing to the geological and geomorphological genesis of the study site there are at least three runoff sources having distinct hydrological characteristics. Results obtained by the use of the three-component mixing model are shown in Figure 8. Based on the concentration data of isotope and chloride ion (Table 1), the contributions of each end member to river water were calculated according to the steady-state mass balance equations of water and tracer fluxes (equations (2)). Isotopic hydrograph separation shows that the contribution of groundwater, precipitation, and glacial meltwater is 66.7%, 19.9%, and 13.4%, respectively. The study indicated that groundwater dominated runoff in the headwater area of Shule River Basin. And the roles of glacier meltwater should be significantly noticed in water resource management in this catchment. The glaciers are the headwaters of many rivers and they affect the water discharge of large rivers [43].

Table 1: Mean, maximum, minimum, and standard deviation values for the concentration of isotope and chloride ion in river water, precipitation, groundwater, and glacial meltwater.

δ ¹⁸O ()
δD ()
Cl⁻ (mg·L⁻¹)

Mean
Max
Min
SD
Mean
Max
Min
SD
Mean
Max
Min
SD

River water
−9.2
−8.5
−9.9
0.36
−62.7
−58.2
−68.9
3.49
11.2
13.4
9.4
1.17

Precipitation
−10.2
−8.3
−13.0
1.56
−69.7
−54.1
−95.1
12.5
0.5
0.8
0.2
0.22

Ground water
−8.1
−6.7
−9.5
0.76
−56.4
−45.3
−68.2
6.44
16.3
19.7
11.0
2.20

Glacial meltwater
−13.2
−12.3
−14.7
0.73
−92.1
−85.1
−105.0
5.63
2.2
3.2
1.2
0.56

Figure 8: The conclusion of hydrograph separation.

Despite the reasonable illustration of the qualitative behavior of runoff components, an exact quantification of runoff components contributing remained difficult and is strongly related to the determination of tracer concentrations in both runoff sources. Owing to the presence of various uncertainties, only qualitative results were achieved. Further experimental investigations are needed to define the tracer concentrations and their variability with greater accuracy. For hydrograph separations in larger scale basins, an extensive consideration of the spatial variability along with the superposition of spatially distributed runoff components is a challenging task for future research. Our results prompt us to focus future work on understanding interannual changes in end member contribution especially in semiarid regions.

5. Discussion and Conclusions

The isotopic and chemical values originate from measurement methods, field data, or the expert knowledge of the investigators. This is reasonable even if the implications of the effects are the same, since all of them cause an uncertain estimation of the end member concentrations and thus of the contribution of different runoff components. So-called end member concentrations need to be defined for every tracer of a specific runoff component for each separation time step in order to calculate the component proportions using mass balance equations for the tracers and the water. However, the determination of these concentrations is often problematic, as it has been shown that they may exhibit high temporal and spatial variability and always include errors caused by the analysis [18, 21]. Therefore, the uncertainty of hydrograph separation results must be addressed. In general, large relative uncertainties must be considered while performing hydrograph separations. Predictive uncertainty is the primary impediment to progress in this area, but continued progress is being made to more fully quantify uncertainty and more fully explore its implications. The importance of reducing errors that have the largest impact on uncertainty is clearly demonstrated. Therefore, future investigations are needed to define with greater accuracy end member tracer concentrations and their spatial and temporal variability. Moving towards an understanding of uncertainties complexity is a challenging and important task for future research in catchment hydrology.

Assumed values of the uncertainty in isotopic composition were , [43]. We calculated the uncertainty of tracers itself to be 9%. Analyses suggested that the uncertainty in the measurement method was less important than that in the temporal and spatial variations of tracer concentrations. The uncertainty terms for precipitation were generally higher than 80% of the total uncertainty, indicating that the δ ¹⁸O values of precipitation account for the majority of uncertainty. The uncertainty was sensitive when the difference between mixing components was small. Therefore, the variation of tracers and the difference of mixing components should be considered when hydrograph separation was applied in the basin.

Hydrograph separation shows that the contribution of groundwater, precipitation, and glacial meltwater is 66.7 ± 6.02%, 19.9 ± 1.79%, and 13.4 ± 1.20%, respectively. There are 347 glaciers and the area of glaciers is 29.45 km², which accounts for 0.72% of the headwater area. Under the background of global warming, rising temperature leads to the increase of snowmelt and accelerating the retreat of glaciers, which will have a significant impact on regional runoff. The roles of glacier meltwater should be significantly noticed in water resource management in this catchment. In addition to temporal variability, a superposition of spatial and temporal distributed runoff components needs to be considered. Moving towards an understanding of this complexity is a challenging and important task for future research in catchment hydrology.

Several studies compared the results of two- and three-component separation. Wenninger et al. (2004) showed a difference of 10% in preevent water contributions between the two methods, because the three-component separation accounted for snow and rain inputs together while the two-component separation accounted for rain inputs only [44]. Dense temporal sampling of hydrographs is often challenging especially at remote locations. Therefore studies that investigate runoff generation in alpine catchment are rare. Pionke et al. showed that for a 7.4 km² watershed three of four monitored storms were dominated by preevent water (55–94% in total) [45]. DeWalle showed that in a smaller catchment (0.198 km²) storm runoff was also dominated by preevent water contributions, 90% over the course of the hydrograph [46]. For a 45 km² catchment, Buda found more than 80% preevent water contributions to the channel storm flow with 67% during the peak flow in a 6 ha catchment [47]. A similar range of preevent water contributions (80%, 60%) was reported by Munyanzea et al. (2012) for two mesoscale catchments (129.3 km², 257.4 km²) in Rwanda [48]. Dong et al. have used the recursive digital filter method and smoothed minimum method to separate base flow based on daily runoff data from 1954 to 2009 in the upper reaches of the Shule River Basin, the recursive digital filter and smoothed minimum method were used for base flow separation. The base flow index is different between the calculation results from the two methods (0.77 and 0.66) [49].

One shortcoming of this study is that the identification of end members is limited to data collected during the vegetation period which comprises only 6 months of the year. Moreover, we only have the isotope data and chemical parameters of one hydrologic section, so we cannot analyze the spatial variability. The importance of reducing errors that have the largest impact is clearly demonstrated; therefore, a targeted sampling strategy is required. In order to fully characterize the range of climatic variability, our results emphasize the need of continued development of the long-term measurement. Our results prompt us to focus future work on understanding interannual changes in end member contribution especially in semiarid regions.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

This research was financially supported by the National Natural Science Foundation of China (nos. 41130638; 41271085; 41401039). The authors are also grateful to all participants in the field (J. H. Yang, H. Wu) and in the laboratory (R. Xu) for their contribution which allows this study to progress in good conditions.

References

J. M. Buttle, “Isotope hydrograph separations and rapid delivery of pre-event water from drainage basins,” Progress in Physical Geography, vol. 18, no. 1, pp. 16–41, 1994.

M. Bonell, “Selected challenges in runoff generation research in forests from the hillslope to headwater drainage basin scale,” Journal of the American Water Resources Association, vol. 34, no. 4, pp. 765–785, 1998.

S. Hoeg, S. Uhlenbrook, and Ch. Leibundgut, “Hydrograph separation in a mountainous catchment—combining hydrochemical and isotopic tracers,” Hydrological Processes, vol. 14, no. 7, pp. 1199–1216, 2000.

Z. S. Wang, C. F. Zhou, B. H. Guan et al., “The headwater loss of the western plateau exacerbates China's long thirst,” AMBIO, vol. 35, no. 5, pp. 271–272, 2006.

B. G. Mark, “Hot ice: glaciers in the tropics are making the press,” Hydrological Processes, vol. 16, no. 16, pp. 3297–3302, 2002.

F. K. Barthold, J. Wu, K. B. Vaché, K. Schneider, H.-G. Frede, and L. Breuer, “Identification of geographic runoff sources in a data sparse region: hydrological processes and the limitations of tracer-based approaches,” Hydrological Processes, vol. 24, no. 16, pp. 2313–2327, 2010.

S. Uhlenbrook and S. Hoeg, “Quantifying uncertainties in tracer-based hydrograph separations: a case study for two-, three- and five-component hydrograph separations in a mountainous catchment,” Hydrological Processes, vol. 17, no. 2, pp. 431–453, 2003.

J. Cable, K. Ogle, and D. Williams, “Contribution of glacier meltwater to streamflow in the Wind River Range, Wyoming, inferred via a Bayesian mixing model applied to isotopic measurements,” Hydrological Processes, vol. 25, no. 14, pp. 2228–2236, 2011.

Y. Kong and Z. Pang, “Evaluating the sensitivity of glacier rivers to climate change based on hydrograph separation of discharge,” Journal of Hydrology, vol. 434-435, pp. 121–129, 2012.

T. Pu, Y. He, G. Zhu, N. Zhang, J. Du, and C. Wang, “Characteristics of water stable isotopes and hydrograph separation in Baishui catchment during the wet season in Mt.Yulong region, south western China,” Hydrological Processes, vol. 27, no. 25, pp. 3641–3648, 2013.

C. Kendall and T. B. Coplen, “Distribution of oxygen-18 and deuteriun in river waters across the United States,” Hydrological Processes, vol. 15, no. 7, pp. 1363–1393, 2001.

F. J. Liu, M. Williams, and N. Caine, “Sources waters and flow paths in an alpine catchment, Colorado Front Range, USA,” Water Resources Research, vol. 40, Article ID W09401, 2004.

J. J. Gibson, T. W. D. Edwards, S. J. Birks et al., “Progress in isotope tracer hydrology in Canada,” Hydrological Processes, vol. 19, no. 1, pp. 303–327, 2005.

C. Wels, R. J. Cornett, and B. D. Lazerte, “Hydrograph separation: a comparison of geochemical and isotopic tracers,” Journal of Hydrology, vol. 122, no. 1–4, pp. 253–274, 1991.

P. Durand, M. Neal, and C. Neal, “Variations in stable oxygen isotope and solute concentrations in small sub-mediterranean montane streams,” Journal of Hydrology, vol. 144, no. 1–4, pp. 283–290, 1993.

N. C. Christophersen, C. Neal, R. P. Hooper, R. D. Vogt, and S. Andersen, “Modelling streamwater chemistry as a mixture of soilwater end- members—a step towards second-generation acidification models,” Journal of Hydrology, vol. 116, no. 1–4, pp. 307–320, 1990.

G. F. Pinder and J. F. Jones, “Determination of the ground-water component of peak discharge from the chemistry of total runoff,” Water Resources Research, vol. 5, no. 2, pp. 438–445, 1969.

J. J. McDonnell, M. Bonell, M. K. Stewart, and A. J. Pearce, “Deuterium variations in storm rainfall: implications for stream hydrograph separation,” Water Resources Research, vol. 26, no. 3, pp. 455–458, 1990.

W. G. Mook, “Environmental isotopes in the hydrological cycle, principles and applications, Volume III: surface water. International Hydrological Programme,” Technical Documents in Hydrology 39, IAEA, Vienna, Austria, 2001.

P. K. Aggarwal, “Isotope hydrology at the International Atomic Energy Agency,” Hydrological Processes, vol. 16, no. 11, pp. 2257–2259, 2002.

N. Christophersen and R. P. Hooper, “Multivariate analysis of stream water chemical data: the use of principal components analysis for the end-member mixing problem,” Water Resources Research, vol. 28, no. 1, pp. 99–107, 1992.

F. J. Liu, R. Parmenter, P. D. Brooks, M. H. Conklin, and R. C. Bales, “Seasonal and interannual variation of streamflow pathways and biogeochemical implications in semi-arid, forested catchments in Valles Caldera, New Mexico,” Ecohydrology, vol. 1, no. 3, pp. 239–252, 2008.

J. Chaves, C. Neill, S. Germer, S. G. Neto, A. Krusche, and H. Elsenbeer, “Land management impacts on runoff sources in small Amazon watersheds,” Hydrological Processes, vol. 22, no. 12, pp. 1766–1775, 2008.

V. Acuña and C. N. Dahm, “Impact of monsoonal rains on spatial scaling patterns in water chemistry of a semiarid river network,” Journal of Geophysical Research G: Biogeosciences, vol. 112, no. 4, Article ID G04009, 2007.

B. Ladouche, A. Probst, D. Viville et al., “Hydrograph separation using isotopic, chemical and hydrological approaches (Strengbach catchment, France),” Journal of Hydrology, vol. 242, no. 3-4, pp. 255–274, 2001.

D. Yang, C. Li, B. Ye et al., Chinese Perspectives on PUB and the Working Group Initiative. Predictions in Ungauged Basins: International Perspectives on the State of the Art and Pathways Forward, vol. 301, IAHS Publication, 2005.

J. Klaus and J. J. McDonnell, “Hydrograph separation using stable isotopes: review and evaluation,” Journal of Hydrology, vol. 505, pp. 47–64, 2013.

C. H. Leibundgut, “Tracer-based assessment of vulnerability in mountainous headwaters,” in Hydrology, Water Resources and Ecology in Headwaters, vol. 248, pp. 317–326, 1998.

Z. Kattan, “Characterization of surface water and groundwater in the Damascus Ghotta basin: hydrochemical and environmental isotopes approaches,” Environmental Geology, vol. 51, no. 2, pp. 173–201, 2006.

W. Liu, S. Chen, X. Qin et al., “Storage, patterns, and control of soil organic carbon and nitrogen in the northeastern margin of the Qinghai–Tibetan Plateau,” Environmental Research Letters, vol. 7, no. 3, Article ID 035401, 2012.

M. J. Liu, T. D. Han, J. Wang et al., “Variations of the components of radiation in permafrost region of the upstream of Shule River,” Plateau Meteorology, vol. 32, no. 2, pp. 411–422, 2013.

Y. Sheng, J. Li, J. C. Wu, B. S. Ye, and J. Wang, “Distribution patterns of permafrost in the upper area of Shule River with the application of GIS technique,” Journal of China University of Mining and Technology, vol. 39, no. 1, pp. 32–40, 2010.

X. Xie, G. J. Yang, Z.-R. Wang, and J. Wang, “Landscape pattern change in mountainous areas along an altitude gradient in the upper reaches of Shule River,” Chinese Journal of Ecology, vol. 29, no. 7, pp. 1420–1426, 2010.

M. G. Sklash and R. N. Farvolden, “The role of groundwater in storm runoff,” Journal of Hydrology, vol. 43, no. 1–4, pp. 45–65, 1979.

M. J. Hinton, S. L. Schiff, and M. C. English, “Examining the contributions of glacial till water to storm runoff using two- and three-component hydrograph separations,” Water Resources Research, vol. 30, no. 4, pp. 983–993, 1994.

A. Rodhe, The origin of streamwater traced by oxygen-18 [Ph.D. thesis], Uppsala University, 1987, UNGI Report Series A No. 41.

D. Genereux, “Quantifying uncertainty in tracer-based hydrograph separations,” Water Resources Research, vol. 34, no. 4, pp. 915–919, 1998.

A. L. James and N. T. Roulet, “Investigating the applicability of end-member mixing analysis (EMMA) across scale: a study of eight small, nested catchments in a temperate forested watershed,” Water Resources Research, vol. 42, no. 8, Article ID W08434, 2006.

W. Dansgaard, “Stable isotopes in precipitation,” Tellus, vol. 16, no. 4, pp. 436–468, 1964.

J. Jouzel, K. Froehlich, and U. Schotterer, “Deuterium and oxygen-18 in present-day precipitation: data and modelling,” Hydrological Sciences Journal, vol. 42, no. 5, pp. 747–763, 1997.

J. Wu, Y. Ding, B. Ye, Q. Yang, X. Zhang, and J. Wang, “Spatio-temporal variation of stable isotopes in precipitation in the Heihe River Basin, Northwestern China,” Environmental Earth Sciences, vol. 61, no. 6, pp. 1123–1134, 2010.

J. K. Wu, Y. J. Ding, J. H. Yang et al., “Stable isotopes in different waters during melt season in Laohugou Glacial Catchment, Shule River basin, Northwestern China,” Journal of Mountainous Science. In press

T. D. Yao, L. Thompson, W. Yang et al., “Different glacier status with atmospheric circulations in Tibetan Plateau and surroundings,” Nature Climate Change, vol. 2, pp. 663–667, 2012.

J. Wenninger, S. Uhlenbrook, N. Tilch, and C. Leibundgut, “Experimental evidence of fast groundwater responses in a hillslope/floodplain area in the Black Forest Mountains, Germany,” Hydrological Processes, vol. 18, no. 17, pp. 3305–3322, 2004.

H. B. Pionke, W. J. Gburek, and G. J. Folmar, “Quantifying stormflow components in a Pennsylvania watershed when ¹⁸O input and storm conditions vary,” Journal of Hydrology, vol. 148, no. 1–4, pp. 169–187, 1993.

D. R. Dewalle and B. R. Swistock, “Differences in oxygen-18 content of throughfall and rainfall in hardwood and coniferous forests,” Hydrological Processes, vol. 8, no. 1, pp. 75–82, 1994.

A. R. Buda and D. R. DeWalle, “Dynamics of stream nitrate sources and flow pathways during stormflows on urban, forest and agricultural watersheds in central Pennsylvania, USA,” Hydrological Processes, vol. 23, no. 23, pp. 3292–3305, 2009.

O. Munyaneza, J. Wenninger, and S. Uhlenbrook, “Identification of runoff generation processes using hydrometric and tracer methods in a meso-scale catchment in Rwanda,” Hydrology and Earth System Sciences, vol. 16, no. 7, pp. 1991–2004, 2012.

W. W. Dong, Y. J. Ding, and X. Wei, “Variation of the base flow and its causes in the upper reaches of the Shule River in the Qilian Mountains,” Journal of Glaciology and Geocryology, vol. 36, no. 3, pp. 661–669, 2014.