Spatiotemporal Variations of Water Stable Isotope Compositions in Nujiang Headwaters, Qinghai-Tibetan Plateau

: The variations of the stable isotope compositions in water provide critical information on hydroclimatic mechanisms. The climatological and hydrological processes in the Nujiang headwaters in the central Qinghai–Tibetan Plateau are extremely complex and are controlled by alternating continental / local recycled and maritime moisture. However, previous studies have only derived limited data from di ﬀ erent types of water in the Nujiang headwaters. Therefore, aiming to understand the sources of stable oxygen ( δ 18 O) and hydrogen ( δ 2 H) isotopes’ compositional variability and how these are related to hydroclimatic processes, we measured δ 18 O and δ 2 H values from surface waters, snow and precipitation across the Nujiang headwaters from April to September 2018. We found higher δ 18 O ( − 13.7% (cid:24) ), δ 2 H ( − 101.8% (cid:24) ) and deuterium excess (d-excess; 7.6% (cid:24) ) values in the non-monsoon season and lower values in the summer monsoon season. Our ﬁndings indicated that the δ 18 O and δ 2 H compositions were signiﬁcantly a ﬀ ected by di ﬀ erent moisture sources in this region. The slope (6.66) and intercept ( − 14.90) of the surface water line (SWL: δ 2 H = 6.66 δ 18 O − 14.90, R 2 = 0.98) were lower than those of the local meteoric water line (LMWL: δ 2 H = 9.50 δ 18 O + 41.80, R 2 = 0.99) and global meteoric water line (GMWL), indicating that precipitation was the primary water vapor source for surface water, and evaporation was the dominant hydrological process for the Nujiang headwaters. In general, δ 18 O and δ 2 H tended to be negatively correlated with precipitation and air temperature. In addition, δ 18 O and δ 2 H values in the Nagqu River were inversely correlated with the intensity of discharge, highlighting a precipitation-driven isotope-discharge pattern. Our ﬁndings provide a theoretical basis for the hydroclimatic mechanisms occurring in the Nujiang headwaters and further augment our understanding of the southern–middle–northern hydroclimate in the Qinghai–Tibetan Plateau.


Introduction
The Qinghai-Tibetan Plateau, known as the "Asian water tower" in the Himalayas region, is the headwater source for the major Asian rivers (e.g., the Yangtze, Yellow, Yarlung Zangbo and Nujiang) [1][2][3]. The hydrology of the plateau has been adversely affected by rapid warming periods and other anthropogenic activities [1,4]. These alterations have induced local climatological, hydrological and ecosystem changes in the headwater region [1,[5][6][7][8]. For instance, the glacier coverage in this region has undergone a massive reduction, which caused increased precipitation and lake areas and higher river runoff [2,6,9]. Therefore, this region has been a research hot spot due to its unique geographic location and huge water circulation [1,10,11]. However, due to lack of interest and knowledge, the mechanisms of these climatological and hydrological changes remain unclear, particularly for the Nujiang catchment located in the core area of the plateau.
Hydrograph separation employing stable isotopes has been one of the effective methods used in hydrographic cycle studies. Precipitation, ground snow, glacial meltwater and groundwater can be differentiated because of their unique profiles of stable isotopes [12]. The relative contributions of different components of runoff during rainfall or snowmelt events can be identified to evaluate streamflow generation mechanisms [8,12]. Recently, stable oxygen (δ 18 O) and hydrogen (δ 2 H) isotope compositions have been widely used in watershed biogeochemical processes [12][13][14][15], particularly in terms of revealing moisture sources and tracing climatological and hydrological processes [8,[16][17][18]. In particular, deuterium excess (d-excess), as the second-order index of δ 18 O and δ 2 H, is a diagnostic tool that can trace water vapor sources in the atmosphere to a specific location [19,20]. Numerous isotopic investigations in the Qinghai-Tibetan Plateau have mainly focused on precipitation mechanisms; for example, isotopes have a significant "amount effect" and "altitude effect" [5,[21][22][23][24][25][26]. However, these stable isotope compositional differences are also related to temperature dependence and changes when the westerlies or Indian monsoon is the sole dominant atmospheric process or when the air temperature is low enough [21,22]. In addition, recent studies have indicated that the stable isotope compositions of streams, rivers and lakes in the Yarlung Zangbo headwaters and Yamdrok Tso Basin could be used to trace and display seasonal variations [8,12,27].
The above previous studies have documented systematic variation of stable isotope compositions in precipitation in the southern Qinghai-Tibetan Plateau; however, scant information was reported specifically on the components and the generation mechanism of precipitation in the central Qinghai-Tibetan Plateau, particularly in rivers and lakes [28]. The Nujiang River is one of the largest rivers in the Qinghai-Tibetan Plateau, and its headwaters are supplied by rain, snow and glacial meltwater [29]. The moisture transport over this region is primarily controlled by the westerlies and Indian monsoons and is strikingly different from that in the northern and southern regions [22]. Owing to the notable influence of both continental air masses and monsoon precipitation on isotopic composition, the behavior of the stable isotopes in precipitation is complex. This complexity is most likely due to the lack of stable isotope composition change tracing information in the Nujiang headwaters, even though this information is immensely important in regional geology, hydrology and climatology [12].
In this study, we analyzed the δ 18 O and δ 2 H values of precipitation, snow, glacier ice, lakes and rivers, which were sampled from April to September 2018. We compared these data with meteorological and hydrological measurements, including air temperature, precipitation and river discharge. Our goal was to generate a preliminary grasp of stable isotope compositions of different waters in the region to broaden the understanding of the stable isotope composition evolution in the water cycle of the Qinghai-Tibetan Plateau.

Study Area
The study was conducted in the Nujiang headwaters of the Qinghai-Tibetan Plateau (Figure 1), which is one of the largest rivers (1.65 × 10 3 km) on the plateau, with an average runoff of 4.96 × 10 10 m 3 , a drainage area of 1.36 × 10 5 km 2 and an average elevation of more than 4000 m [30]. Upstream, the Nujiang River originates from the Tanggula Mountain Glacier (6000 m) and flows through the alpine meadows of the central plateau. Further downstream, the river passes through major cities such as Nagqu, Baoshan and Tengchong and finally flows into South Asia. The Ciqu and Seni Rivers comprise the primary and secondary tributaries of the main flow (Nagqu River) of the Nujiang, with lengths of Sustainability 2020, 12, 6654 3 of 15 104 and 30 km, respectively. These two rivers run through Nagqu city, and water quality is strongly affected by anthropogenic activities. Cuona Lake absorbs the main river flows and is the highest (4650 m) freshwater lake in the world, covering an area of 300 km 2 .
Sustainability 2020, 12, x FOR PEER REVIEW 3 of 16 is strongly affected by anthropogenic activities. Cuona Lake absorbs the main river flows and is the highest (4650 m) freshwater lake in the world, covering an area of 300 km 2 .
Overall, most of the flow of these selected rivers and this lake originated from precipitation, glaciers, snowmelt and underground water, which were influenced by climate change [1,29]. As mentioned earlier, global warming and anthropogenic activities led to a reduction in the size of the plateau's glaciers, which significantly increased local precipitation and river runoff [2,9,31].

Sampling and Isotopic Analysis
Water samples were collected from the Nagqu (sites 1/4/5/11), Ciqu (6/7) and Seni (8/9/10) Rivers and Cuona Lake (2/3) into 60 mL polyethylene bottles with a 60 mL syringe at a water depth of 0-30 cm. Sampling at each site was repeated thrice, and then 20 mL water samples were collected and mixed together to make a composite sample. To avoid sample leakage and evaporation, the bottles were filled completely and immediately sealed. Ice samples were collected from vertical sections of an ice body at a depth of up to 50 cm, with each sample separated by 10 cm in depth, and collected in three replicates and mixed together to make a representative sample for analysis. Collected samples were first sealed in double-layer preservation plastic bags and melted at room temperature. Then, the melted water samples were transferred to polyethylene bottles as per above. The precipitation sample collection process was similar to that of rivers and lakes, while the process of snowfall sample collection was similar to that of ice. There were 53 samples in total, including 26 Overall, most of the flow of these selected rivers and this lake originated from precipitation, glaciers, snowmelt and underground water, which were influenced by climate change [1,29]. As mentioned earlier, global warming and anthropogenic activities led to a reduction in the size of the plateau's glaciers, which significantly increased local precipitation and river runoff [2,9,31].

Sampling and Isotopic Analysis
Water samples were collected from the Nagqu (sites 1/4/5/11), Ciqu (6/7) and Seni (8/9/10) Rivers and Cuona Lake (2/3) into 60 mL polyethylene bottles with a 60 mL syringe at a water depth of 0-30 cm. Sampling at each site was repeated thrice, and then 20 mL water samples were collected and mixed together to make a composite sample. To avoid sample leakage and evaporation, the bottles were filled completely and immediately sealed. Ice samples were collected from vertical sections of an ice body at a depth of up to 50 cm, with each sample separated by 10 cm in depth, and collected in three replicates and mixed together to make a representative sample for analysis. Collected samples were first sealed in double-layer preservation plastic bags and melted at room temperature. Then, the melted water samples were transferred to polyethylene bottles as per above. The precipitation sample collection process was similar to that of rivers and lakes, while the process of snowfall sample collection was similar to that of ice. There were 53 samples in total, including 26 samples of rivers, 8 samples of the lake, 8 samples of snow, 6 samples of precipitation and 5 samples of ice. All water samples were stored at 4 • C until analysis.
Two-milliliter water samples from the polyethylene bottles were extracted with a 2.5 mL syringe and filtered with a membrane filter with a 0.22 µm pore diameter. The filtrates were injected into a 2 mL sample bottle, and the cap was quickly screwed on. The δ 18 O and δ 2 H values of all water samples were analyzed by cavity ring-down laser spectroscopy using the Vario MICRO Cube elemental analyzer (Elementar Analysensysteme GmbH, Langenselbold). The analytical uncertainties were ±0.1% for δ 18 O and ±0.5% for δ 2 H. The reported results were relative to the standard Vienna Standard Mean Ocean Water (VSMOW), defined as

Environmental Variables
The meteorological and hydrological data (e.g., air temperature, precipitation and discharge) were obtained from weather and hydrometric stations in Nagqu city. At the same time and location, the water temperature and pH of the water samples were assistive measured with a calibrated portable handheld meter (Thermo Fisher Scientific, Singapore).

Statistical Analysis
All statistical analyses were performed on IBM SPSS Statistics version 19.0 (IBM, Armonk, NY, USA). Linear regression was used to determine whether there was a linear correlation between isotope composition and measured environmental variables. Figure 2 illustrates the daily variations in air temperature, relative humidity, discharge and precipitation in the Nujiang headwaters from April to September 2018. The average air temperature was the highest in July (11.3 • C) and the lowest in April (0.4 • C) (Figure 2a). The relative humidity ranged from 44.69% to 69.80%, with the highest occurring in September and lowest occurring in April ( Figure 2b). The discharge of the Nagqu River from July to September (average 61.94 m 3 s −1 ) was significantly greater than that from April to June (average 8.11 m 3 s −1 ) (Figure 2c). During the six-month study period, July (187.2 mm) and August (120.0 mm) had the most monthly precipitation, accounting for 36.7% and 23.5% of the total, respectively ( Figure 2d).

Meteorological and Hydrological Variations during the Study Period
The water temperatures of the Nagqu River, the Ciqu River and Cuona Lake were 5.44-15.86 • C, 8.75-12.33 • C and 10.18-20.62 • C, with averages of 11.63, 11.33 and 14.73 • C, respectively ( Table 1). The pH of Cuona Lake water ranged from 8.35 to 9.17 (average 8.84) and was the highest from June to August and the lowest from April to May and in September. The pH variation of the Nagqu River was similar to that of Cuona Lake, but the pH of the Ciqu River did not vary appreciably from month to month.

Variations of δ 18 O and δ 2 H Values
All measured δ 18 O and δ 2 H values in different waters are listed in Table 2. We examined the δ 18      The isotopic values from upstream to downstream ending at Cuona Lake in the Nujiang headwaters could be clearly distinguished. The δ 18 O and δ 2 H values of Cuona Lake were significantly higher than those of the Nagqu and Ciqu Rivers. The relative values of δ 18 O (−11.7‰) and δ 2 H (−95.9‰) at Nagqu River sampling site 4 were also higher than those of all other sites. For comparison, sampling sites downstream of monitoring site 4 displayed no significant variations in isotopic values of the Nagqu and Ciqu Rivers (Figure 4a,b).
Isotopic tracers of δ 18 O and δ 2 H indicated that these values of snow were directly related to elevation (Figure 4c). By contrast, δ 18   The isotopic values from upstream to downstream ending at Cuona Lake in the Nujiang headwaters could be clearly distinguished. The δ 18 O and δ 2 H values of Cuona Lake were significantly higher than those of the Nagqu and Ciqu Rivers. The relative values of δ 18 O (−11.7% ) and δ 2 H (−95.9% ) at Nagqu River sampling site 4 were also higher than those of all other sites. For comparison, sampling sites downstream of monitoring site 4 displayed no significant variations in isotopic values of the Nagqu and Ciqu Rivers (Figure 4a,b).
Isotopic tracers of δ 18 O and δ 2 H indicated that these values of snow were directly related to elevation (Figure 4c). By contrast, δ 18 O and δ 2 H values varied between −109.6% and −112.8% at −10 to −40 cm under the surface ice and increased slightly at −50 cm (Figure 4d).

The Relationship between δ 18 O and δ 2 H Values
Based on the stable isotope composition of precipitation, snow and surface water (ice, the lake and the rivers) in the Nujiang headwaters, we established the local meteoric water line (LMWL) and the surface water line (SWL) as δ 2 H = 9.

Correlations of Isotopic Composition and Environmental Factors
Previous studies have demonstrated that δ 18 O and δ 2 H values in water were linked to air temperature, relative humidity, precipitation and river discharge, which represented the

Correlations of Isotopic Composition and Environmental Factors
Previous studies have demonstrated that δ 18 O and δ 2 H values in water were linked to air temperature, relative humidity, precipitation and river discharge, which represented the "temperature effect", "amount effect" and "altitude effect" [8,33,34]. In general, we found that δ 18 O and δ 2 H in all waters were negatively correlated with precipitation and air temperature (Figure 6a,b,d). Meanwhile, we also found δ 18 O and δ 2 H values in the Nagqu River were inversely proportional to the intensity of discharge (Figure 6e,f).

Correlations of Isotopic Composition and Environmental Factors
Previous studies have demonstrated that δ 18 O and δ 2 H values in water were linked to air temperature, relative humidity, precipitation and river discharge, which represented the "temperature effect", "amount effect" and "altitude effect" [8,33,34]. In general, we found that δ 18 O and δ 2 H in all waters were negatively correlated with precipitation and air temperature ( Figure  6a,b,d). Meanwhile, we also found δ 18 O and δ 2 H values in the Nagqu River were inversely proportional to the intensity of discharge (Figure 6e,f).

Stable Isotopic Precipitation Patterns
Seasonal variations of isotopic composition for precipitation provide records of moisture sources and transport paths of air masses [5,8,19]. During the monitoring period, we found that precipitation

Stable Isotopic Precipitation Patterns
Seasonal variations of isotopic composition for precipitation provide records of moisture sources and transport paths of air masses [5,8,19]. During the monitoring period, we found that precipitation showed distinct changes in isotopic composition. δ 18 O and δ 2 H values were elevated in the pre-and post-monsoon seasons (April−May, September) and depressed in the monsoon season (June−August), which was directly related to precipitation (Figure 3). The high isotopic values were coincident with low precipitation, and the low values were coincident with high precipitation, suggesting that δ 18 O and δ 2 H were closely related to precipitation [24,25,35,36]. Similar results have been found for atmospheric relative humidity in the Qinghai−Tibetan Plateau [36]. These patterns revealed higher isotopic values of continental moisture in the non-monsoon season and lower values of maritime moisture in the monsoon season [22,35,37,38]. As the summer monsoon strengthened, precipitation and relative humidity increased correspondingly (111.5 and 398.1 mm and 55.12% and 65.94%, respectively) (Figure 2b,d). Long-distance precipitation transport results in relatively lower isotopic values in the Nujiang headwaters during the monsoon season, which progressively produces rainfall from the air mass during orographic lifting and adiabatic cooling [8,39,40]. From 2000 to 2003, the weighted average annual precipitation for δ 18 O was −16.2% in Nagqu [21]. In comparison, δ 18 O in the average precipitation in the study area was −19.3% (Figure 3a) and significantly reduced by 19.1%. As previously reported, the warming rate in this region was as high as 0.2 • C per decade from 1960 to 2010 [41], indicating that δ 18 O in precipitation was inversely correlated with air temperature. We also found a significant correlation between δ 18 O in precipitation and air temperature (R 2 = 0.80, P < 0.05) (Figure 6c).
Precipitation generally displayed high d-excess values (17.6% ) in the non-monsoon season and low (6.3% ) values in the monsoon season (Figure 3c). These results were similar to those of the southern plateau regions of Nyalam and Lhasa [25,26,42] and the central region of Nagqu [19]. The global average of d-excess in precipitation was 10% , and seasonal variations could provide information about moisture sources, the hydrological cycle and local climatic conditions. Lower d-excess values are indicative of marine moisture and higher values of continental and local recycled moisture [19,25,32,38,[42][43][44][45]. Our results indicated that summer moisture was derived from marine moisture in the Indian Ocean, while spring and autumn moisture was dominated by continental and local recycled moisture as well as westerly transport. There were two different paths for marine moisture to reach the Nujiang headwaters. One was to move along the Brahmaputra-Yarlung Zangbo River, and the other was over the Indian Peninsula crossing the Himalaya Mountains [19,46]. In the Nujiang headwaters, precipitation d-excess ranged from 5.5% to 20.2% (average 12.0% ), which was significantly higher than that of the Himalayas and the southern plateau (0.8% -6% ) and globally (10% ) [25,37,38,42]. Compared with other global precipitation investigations, precipitation in the Nujiang source region had a relatively higher d-excess, presumably associated with the lower local condensation temperature [47], dry atmospheric vapor and higher elevation [8]. Our investigation area was located in the core of the Qinghai−Tibetan Plateau, with an average elevation of more than 4500 m [48]. This unique geographical location let the region maintain low temperatures throughout the year and had a shorter distance between precipitation and clouds, which resulted in higher δ 18 O values in precipitation due to weaker evaporation [25,26].

Stable Isotopic Pattern in Snow and Ice
Although we used only one sampling site for five different depths of ice, a systematic isotopic pattern could still be observed. The stable isotope compositions of the ice varied between −112.8% and −109.6% for δ 2 H and −15.0% and −14.5% for δ 18 O from 10 to 40 cm under the surface ice and then up to 50 cm (Figure 4d). In general, the values of stable isotope compositions measured from ice were inversely correlated with the ice depth. By contrast, this pattern was not found in the northeast region (Qiyi Glacier) of the plateau, where δ 18 O values of ice increased with the increased depth [49]. This discrepancy was most likely associated with the new snow covering the ice surface and snowmelt flows due to rising temperatures migrating to the bottom of the snow pit along with the profile [49]. In our study, we found that the isotopic enrichment values (δ 18 O: −11.8% , δ 2 H: −70.9% ) of snow water were significantly higher than those of ice water (δ 18 O: −14.8% , δ 2 H: −111.4% ) (Figure 4a). Remarkably, the δ 18 O values of snow increased with the increasing altitude in the Nujiang headwaters (4560-4950 m; Figure 6a), which was contrary to other results (0.12% -0.21% /100 m) [24]. This significant difference was closely related to evaporation, melting, sublimation and recrystallization, leading to dynamic alterations in snow water isotopic compositions [12].
The average d-excess (δ 2 H − 8 × δ 18 O) in precipitation (12.0% ) and snowmelt (19.9% ) was significantly enriched over that of ice water (7.1% ) (Figure 3c), and the slope values of precipitation (9.43) and snowmelt (9.24) were significantly greater than that of ice water (5.58) ( Figure 5). This suggested that the intermediate values of d-excess in ice water were indicative of mixing between precipitation and snowmelt (Ren et al., 2016(Ren et al., , 2017. Moreover, the extremely high d-excess (22.24% ) in the snowmelt near Nujiang pointed to the predominance of continental recycled moisture as the source during the sampling period. Likewise, the high d-excess could evaluate the contribution of snowmelt to river discharge, which reflected the delayed recharge caused by spring precipitation in snow.

Stable Isotopic Patterns in Surface Water
Our sampling points only included two major tributaries (the Ciqu and Seni Rivers) and one lake (Cuona Lake), but the δ 18 O and δ 2 H values were still consistent with the stable isotopic trend of the main flow of the Nagqu River. From sampling sites 1 to 4, the δ 18 O and δ 2 H values of the main flow increased from −15.5% to −11.7% and −115.5% to −95.9% (Figure 4a,b), respectively, indicating a strong isotope enrichment, which was similar to other lakes and rivers of the Qinghai−Tibetan Plateau [12,27]. The phenomenon can be attributed to the strong evaporation processes over the open water, and the flow-through lake will preferentially leave heavy isotopes in the lake [12,27]. The isotope compositions of Cuona Lake were similar to the GMWL, and the lower slope indicated strong surface water evaporation (Figure 5b) [12]. In addition, δ 18 O in lake water was generally more than 10% higher than in other waters [27]. Accordingly, we also found that the main flow values for δ 18 O and δ 2 H downstream of site 4 declined sharply and then were similar to those of the tributaries. The Ciqu River and Seni River were the two major tributaries of the main flow, where the stable isotope composition was low and strongly affected by anthropogenic activities. Thus, the dilution caused by the input of these tributaries was clearly the leading cause for the progressive depletion of δ 18 O and δ 2 H in the main flow [12]. However, the δ 18 O and δ 2 H values of the upper Nujiang headwaters decreased and were characteristic of the downstream values, indicating the effect of elevation on precipitation isotopes and the cumulative effect of evaporation as water flows downstream [12].
We also collected discharge data from the Nagqu River throughout the monitoring period to examine the isotopic seasonal patterns in the main flow and their response to discharge variability. The discharge of the Nagqu River in the monsoon season was significantly higher than that of the non-monsoon season ( Figure 2). The δ 18 O and δ 2 H values of the main flow were enriched from May to June (Figure 3), which was generally consistent with the seasonal variations of precipitation isotopes in this region, indicating the δ 18 O and δ 2 H values were inversely correlated with the intensity of the discharge (Figure 6e,f). The monsoon precipitation-driven isotope-discharge pattern was similar to those previously reported [8,14,40]. The variations in isotopic values of the tributaries were almost identical to those in the main flow. In comparison, Cuona Lake showed no discernible variation (Figure 3). From a climatological perspective, the discharge of the main stream was primarily derived from continental or local recycled moisture and westerly transport in the non-monsoon season. During the monsoon season, the discharge was derived from the Indian Ocean's marine moisture. However, moisture originating from these two climatic systems exhibited distinct isotopic signatures, which were higher during the westerlies and lower during the summer monsoon [8,37]. This significant seasonal difference was due to the seasonality of the mainstream in the transition zone between the westerlies and the summer monsoon [21,22]. There was less flow in the channel course during the spring and early summer, while melting ice and snowmelt as well as monsoon moisture significantly increased the flow from July to September (Figure 2). In this study, δ 18 O and δ 2 H values of precipitation and snow were relatively lower in the monsoon season, leading to lower values in the main flow. Remarkably, the lowest discharge in June was not accompanied by the highest isotopic values in our study. Additionally, the highest discharge values in the Yarlung Zangbo River Basin did not result in its lowest isotopic values, which may be partly due to isotope depletion or enrichment of water sources, such as groundwater and delayed runoff from pre-monsoon precipitation that accumulated at higher elevations [8].

Surface Water Lines
Linear regression analysis was performed on the measured surface samples in the Nujiang headwaters to establish the SWLs as δ 2 H = 6.66 δ 18 O − 14.90 (R 2 = 0.91; Figure 5b). In general, δ 18 O and δ 2 H from ice, lake and river water varied in a systematic manner and fell close to the GMWL of δ 2 H = 8 δ 18 O + 10 [32]. The GMWL has proven to be a useful reference line to understand the spatial variability of water isotope content, especially the evaporation effect based on the deviation from water [25,32,50]. In our study, the slope (6.66) and intercept (−14.90) of the Nujiang SWL were lower than those of the GMWL. Evaporation could cause deviations in this arid climate, resulting in a reduced slope and increased isotopic enrichment [8,25,26]. Water lines varied from region to region and were closely related to geographical settings and individual climatic conditions (e.g., moisture source, condensation temperature and subcloud evaporation) [8]. For instance, the SWLs in the Yarlung Zangbo River Basin were defined as δ 2 H = 9.25 δ 18 O + 24.07 and δ 2 H = 8.16 δ 18 O + 12.77 [8,12]. Interestingly, as the lines moved northward, the SWL varied systematically and showed a progressively lower slope and intercept [12,25,26,51]. Our results also confirmed this, and the slope and intercept continued to decrease to 6.66 and −14.90, respectively. The Nujiang headwater region, located in the north of the Yarlung Zangbo River Basin, had higher evaporation than the Yarlung Zangbo River Basin. After the Indian summer monsoon crossed the Himalayas, oceanic moisture would have a longer transport distance from the south to the observation area and consequently experience stronger evaporation, especially under arid climate conditions [25,26]. Therefore, compared with the Yarlung Zangbo River Basin, summer moisture transported by the southwest monsoon would result in lower slope and intercept values in the Nujiang Basin. The δ 18 O−δ 2 H relationship of the Yarlung Zangbo River Basin deviated from the GMWL because of the rapid flow rates due to steep elevation gradients and cold temperatures [12]. However, the same elevation gradients and cold temperatures existed in the Nujiang Basin, indicating that the δ 18 O−δ 2 H slopes of the regression lines were less than eight due to evaporation. The atmospheric relative humidity in the Nujiang Basin was 60.56%, which was significantly lower than that in the Yarlung Zangbo River Basin and strongly supported our conclusion ( Figure 2).
A previous study had stated that the correlations between water δ 18 O and δ 2 H values could be effectively used to identify hydrological processes along river channels [8]. We found that precipitation and snow demonstrated strong linear correlations for δ 2 H = 9.50 δ 18 O + 41.80 (R 2 = 0.99), while surface water demonstrated weak linear correlations for δ 2 H = 6.66 δ 18 O − 14.90 (R 2 = 0.98). As a comparison, the slope and intercept of the SWL were much lower than those of the LMWL, indicating precipitation and snow were the primary water vapor source for surface water. Additionally, the SWLs of the Nagqu River, the Ciqu River, the Seni River, Cuona Lake and ice were δ 2 H = 7.01 δ 18 [8,12]. In addition, the slopes of the SWLs were lower than that of the GMWL (Figure 5), which also indicated that evaporation was the main hydrological process in the surface waters. Moreover, the spatiotemporal trends of the water parameters supported the above conclusion [8,52,53]. In our study, compared with Cuona Lake, the water temperature and pH in the main flow showed an increasing trend during the monitoring period. From April to May, the water temperature and pH in the main flow also exhibited an increasing trend, while in September, there was a decreasing trend compared to the tributaries (Table 1).

Conclusions
Our present work is the first systematic study of δ 18 O and δ 2 H values in various types of waters that supply the Nujiang headwaters in the central Qinghai−Tibetan Plateau. We found that the moisture in the Nujiang headwaters was controlled alternately by the westerlies and the summer monsoon, with higher isotopic values in the non-monsoon season. Across the Nujiang headwaters, δ 18 O and δ 2 H values in lake water were more enriched than local precipitation and other surface waters due to stronger evaporation. In addition, the isotopic values of the main flow were significantly affected by Cuona Lake and were higher upstream of the lake. We established an SWL of δ 2 H = 6.66 δ 18 O − 14.90 for the Nujiang headwaters, which was lower than the GMWL. Based on the above results, we concluded that water mixing and evaporation appeared to be the dominant hydrological processes in this region. In general, δ 18 O and δ 2 H values in the water tended to be negatively correlated with precipitation and air temperature. Moreover, the isotopic values of the Nagqu River were inversely proportional to the intensity of discharge, highlighting the precipitation-driven isotope-discharge pattern.