The O and H isotope characteristics of water from major rivers in China

•••••

The O and H isotope characteristics of water from major rivers in China

Si-Liang Li•Fu-Jun Yue•Cong-Qiang Liu•Hu Ding•Zhi-Qi Zhao•Xiaodong Li

Isotopic compositions of the water sources and evapotranspiration in large rivers of China are not well explored.In this investigation,δ18O and δD values of water from major large rivers of China are integrated to understand thespatialandtemporalvariationsintheisotopesofwaterand d-excess.River water samples were collected from the channelsandtributariesofthe Changjiang,Yellow,Songhua, andLiaoRiversinChina.Waterfromupstreamwasgenerally enriched with light isotopes relative to that from downstream areas.Water from the Songhua River and upper reach of the Changjiang River has relatively negative isotopic values, impacted by latitude and altitude effects,respectively.The spatial pattern showed that δ18O-H2O values in rivers mainly ranged from-9 to-7‰in most of China based on roughly mapped contours,and rivers enriched in light isotopes were primarily located in the Tibetan Plateau and northeast China. The slopes of linear regression from the Yellow and Liao Rivers suggested that the loss of water by direct evaporation shifted the isotopes of water,especially in some small tributaries.Additionally,there was a narrow range of low δ18OH2O values from the main channel of rivers during the lowfowseasonrelativetothehigh-fowseasonduetoinfuenceof themonsoonclimate.Nodistinctspatialvariationofd-excess values was found in the rivers of China.However,water samples in some rivers with low d-excess values had high δ18O values,indicating that evaporation in the tributaries has led to signifcant water loss,especially in tributaries of the Yellow and Liao Rivers.

Stable isotopes·Spatial variation·d-Excess· Major rivers·China

1 Introduction

River water originates from alpine snow/ice melting,recent precipitation,and groundwater in the basin,and serves as the major water source for cities and ecological systems located alongside the rivers.Therefore,it is important to understand the sources and evapotranspiration of water to improve the effects of best management practices on such resources.Isotopes are ideal tracers for investigation of hydrological processes such as transport and mixing because they provide a conservativefngerprintthatisintrinsictothewatermolecules (Gat 1996;Clark and Fritz 1997;Criss 1999).The δ18O and δDvaluesofwaterinriversandstreamsarecommonlyusedto estimate water sources,through examination of proportions among different reservoirs and evaporation processes within catchments(Dalaietal.2002;Lambsetal.2009;Bowenetal. 2011).Accordingly,spatial isotopic data can be used to reconstruct continental palaeoclimate and palaeohydrology since climatic signals controlling the isotopic composition of precipitation samples have been preserved in river samples (Hren et al.2009;Bershaw et al.2012).Additionally,the isotopicfngerprintofwatersamplescanhelpevaluateaspects of biogeochemical cycling of nutrient and microbial production that are dependent on the water cycle,such as the formation of nitrate,sulfate,and phosphate via oxidation processes(Michener and Lajtha 2007).

The isotopic compositions of water from different areas are affected by meteorological processes,which provide a characteristic fngerprint of their origin(Clark and Fritz 1997).The isotopic values of river water decrease northward roughly parallel to latitude in the eastern UnitedStates,while they decrease eastward from the Pacifc Ocean to the top of the Rocky Mountains and then increase eastward along the eastern fanks of the Rockies in the western United States(Kendall and Coplen 2001).The waters of major Indian rivers are enriched with light isotopes in Himalayan rivers and with heavy isotopes in peninsular rivers in a northwest-southeast gradient(Lambs et al.2005).The isotope signature of precipitation was found to be reasonably well preserved in both monsoon and non-monsoon seasons according to the regression analysis of δD-δ18O data describing the Yamuna River(Dalai et al. 2002).Moreover,isotope balance calculations in the Indus River basin indicated that up to 72%of the Indus discharge close to its mouth might be derived from the Mediterranean end-member(Karim and Veizer 2002).

Many studies have investigated the isotopic composition of precipitation and controlling factors in China,especially in the Himalayas and western China(Ding et al.2009;Tian et al.2001,2007;Yao et al.2009;Wen et al.2012;Zhao et al.2011,2012).However,a few data are available for individual rivers such as the Yellow River(Su et al.2003), Changjiang(Li et al.2010;Ding et al.2013),and other small rivers in western China(Liu et al.2009;Yu et al. 2010).Moreover,the spatial isotopic compositions of rivers have not been explored well in China.Therefore,in this study,data from our studies and other studies of major rivers of China are integrated to enable a better understanding of spatial and temporal variations in isotopes of water as well as d-excess.

2 Materials and methods

2.1 Characteristics of the major rivers in China

The major rivers in China are shown in Fig.1.The two largest rivers(the Changjiang River and the Yellow River) originate from the Tibetan Plateau,and the basin area drained by these two rivers is approximately 25%of the total area of mainland China.The Changjiang and Yellow Rivers are the two longest rivers in China.Both rivers originate in the eastern Tibetan Plateau at an elevation higher than 4,000 m.The Changjiang River has a basin area of 180×104km2with a mainstream river length of over 6,300 km,while the Yellow River has a basin area of 75.2×104km2with a mainstream river length of over 5,460 km(Table 1;Zhu 2007).The Songhua River Basin and Liao River Basin are located in northeast China.TheSonghua River Basin is bounded by the Changbai Mountains in the east,Daxingan Mountains in the west,Xiaoxingan Mountains in the north,and Songliao watershed divide in the south;the Songhua drains an area of 54.6×104km2and has a length of 2,030 km.The Liao River has an area of 2.3×105km2and a length of 1,430 km.The Pearl River originates in Maxiong Mountain of the Yunnan-Guizhou Plateau at above 2,000 m and has a basin area 44.2×104km2and a mainstream river length of 2,214 km that runs through southern China;it has second greatest discharge of all the rivers in China.

Fig.1 Major rivers,spatial variation of annual rainfall (dashed line)and moisture masses in China

Table 1 The characteristics of major rivers in China(Zhu 2007)

Southern and eastern China are impacted by a southeast monsoon from the Pacifc Ocean and a southwest monsoon from the Indian Ocean,which bring high rainfall in summer (Fig.1).Additionally,the westerly and minor polar air mass in winter provide moisture for northern and western China (Tian et al.2007).The majority of the Changjiang and Pearl River Basins are located in a humid subtropical zone subject to a typical monsoon climate,with precipitation that generallyrangesfrom800 to 2,000 mm/a.TheYellowRiver fows across the dry Loess Plateau of northern China,where precipitation is relatively low(approximately 500 mm/a).The Songhua and the Liao Rivers belong to the temperate continental monsoon climate zone,with some areas belonging to the cold temperate climate zone,which receives precipitation of 400 mm/a in the west and 800 mm/a in the east.

2.2 Sampling and analysis method

Water samples were collected by us from the main channel and the major tributaries of the Changjiang(Li et al.2010), Liao(Yue et al.2013),Yellow,and Songhua Rivers(Yue et al.2014).One sampling campaign was carried out for the Changjiang during August,2006.For other rivers, water samples were collected during the high fow season of 2010 and the low fow season of 2010,respectively.The δ18O and δD of water were measured using a Micromass IsoPrime mass spectrometer coupled to an automated line after equilibration between O-water and CO2gas(δ18O), and between H-water and H2gas with a Pt catalyst(δD. The analytical precision of the δ18O-H2O was generally better than 0.2‰.δD-H2O has a precision of 3‰ for water samples in the Changjiang River(Li et al.2010)and 2‰for other rivers.

3 Results

There was wide variation in isotopic values from-14.7 to -3.0‰forδ18Oand-102.4to-32.2‰forδDintherivers. The δ18O-H2O ranged from-14.7 to-6.0‰with a mean value of-9.7‰,while the δD-H2O ranged from-101.2 to -32.2‰with a mean value of-72.4‰in the Changjiang River.Waterfromupstreamwasgenerallyenrichedwithlight isotopesrelativetowaterfromdownstreamintheChangjiang. Theδ18Oofwaterincreasedfrom-13.8to-8.3‰alongthe mainchanneloftheChangjiangRiverduetothehighisotopic values of water from tributaries(Fig.2).The δ18O and δDH2Owerefoundtorangefrom-10.7to-3.0‰and-74.7to -35.8‰in the YellowRiver,with meanvaluesof-9.3 and -64.9‰,respectively.An increasing trend in light δ18OH2O was observed in the main channel of the Yellow River during the high fow season(Fig.3),similar to that of the Changjiang River.During the low fow season,the water hadslight negative values of δ18O-H2O relative to the high fow season in the Yellow River.

Fig.2 Oxygen isotopic compositions and altitude as functions of the distance from the source of the Changjiang River.Circlesrepresents data from Li et al.(2010),squaresandtrianglesrepresent data between 2003 and 2007 from Ding et al.(2013)

Fig.3 Oxygen isotopic compositions and altitude as functions of the distance from the source of the Yellow River

In the Songhua River,the river water contained stable hydrogen and oxygen isotope ratios ranging from-102.4 to -49.5‰with a mean value of-83.8‰and from-14.5 to -5.8‰with a mean value of-11.5‰,respectively.The waterwithlowδ18O-H2Ovalueswaslocatedattributariesand headwaters originating from the mountain.There was a narrow range of δ18O-H2O values in water collected during low fowseasonfromtheSonghuaRivermainchannel,withwater samples having δ18O-H2O values between-11.6 and -9.9‰.The water samples were generally enriched with heavy isotopes of water from the Liao River when compared withwaterfromSonghuaRiver,whichhadδ18O-andδD-H2O values ranging from-10.7 to-3.7‰with a mean value of -8.0‰and from-77.6 to-34.6‰with a mean value of -60.4‰,respectively.

4 Discussion

4.1 Comparison of isotopes of river water and precipitation in different basins

Stable isotope ratios in precipitation found to form the Global Meteoric Water Line(GMWL)were defned as δD=8δ18-O+10 on a global scale(Craig 1961).The more accurate GMWL of δD=8.17δ18O+10.35 was then given based on the large database from Global Network of Isotope in Precipitation(GNIP)sites located worldwide(Rozanskiet al.1993). Localprecipitationshouldbeamajorsourceofwaterforrivers, while evaporation would increase the oxygen and hydrogen isotope ratios of the residual fraction and the residual surface water to produce systematic linear deviations from the local meteoric water line(LMWL),which can be used to track the source of water and as tracers of hydrogeological processes (Clark and Fritz 1997;Wassenaar et al.2011).

LMWLs in different basins were generated from the database of GNIP sites(IAEA 2006)based on locations in China(Fig.4).A signifcant difference between LMWL and GMWLindifferentbasinswasobservedexceptLMWLofthe Changjiang Basin,which showed local water sources were impacted by different air masses and evaporation processes. TheδD-δ18Oplotsofthedatafromdifferentriversare shown in Fig.4.The Changjiang River water line produced by data from previous studies(Li et al.2010;Ding et al.2013)is δD=6.62δ18O-2.71(R2=0.88,p＜0.005),which differed from that of the Chanjiang Basin somewhat,with an LMWLofδD=7.83δ18O+9.65(R2=0.92)beingpresent in the Changjiang Basin.These fndings suggest that multi water sources in the upper reach and evaporation in some tributariesorlakesinthemiddleandlowerreachesresultedin the difference.Linear regression of the Yellow River water line results in the equation of δD=4.91δ18O-19.55 (R2=0.83,p＜0.005),which has a slope lower than that from the LMWL of δD=6.86δ18O-0.95(R2=0.93)in the Yellow Basin.The river water lines of the Songhua and Liao Rivers generated from water isotopic analyses were δD=5.62δ18O-19.00 (R2=0.92, p＜0.001) and δD=4.91δ18O-20.93(R2=0.90,p＜0.001),respectively.ThelocalmeteoricwaterlinesoftheSonghuaandLiao Basinswere δD=7.48δ18O-0.98 (R2=0.97)and δD=4.28δ18O-22.57(R2=0.84),respectively.In western China,evaporation was the major reason for the shift in isotopes of water due to limited precipitation and high evapotranspiration rates,which lead to residual surface water beingenrichedwithheavyisotopes(Suetal.2003;Zhaoetal. 2011).Indeed,LMWL slopes of less than fve have been reportedinthe Great Basin andother aridpartsofthe western UnitedStates(KendallandCoplen2001).Theslopesoflinear regression from the Yellow River and the Songhua River are lower than those of the LMWL,which suggests that signifcant loss of water by direct evaporation shifts the isotopes of water,especially in some individual tributaries.Another possibleexplanationforthedifferencebetweentheriverwater line and LMWL in the Songhua Basin might be the limited data describing precipitation from just two cities in GNIP sites.Similar slopes of regression analysis of δD-δ18O data between river water and precipitation were observed at New Delhiduringthesameperiod(Dalaietal.2002).Itshowedthat there was little evaporation from river surfaces or soils in the hydrological continuum with surface water based on similar slopesofthelinearregressionbetweentheriverwaterlineand LMWL in the Indus River(Karim and Veizer 2002).The slopes of the Liao River Water Line and the LMWL were generallylessthanfve,andtherewasnosignifcantdifference inslopesbetweenthe LiaoRiver Water Line andthe LMWL. Thepatternofslopemightindicatethatevaporationprocesses from surface water did not play a major role in formation ofisotopic compositions in the Liao River during the sampling periods.

Fig.4 Scatter plots of δ18O and δD in different rivers of China.The slopes of regression in the Changjiang River and Yellow River are generated by all plots in the fgure

4.2 Spatial variation of δ18O-H2O in major rivers of China and controlling factors

The isotopic compositions of precipitation are infuenced by temperature,altitude,continental effects,latitude,and amount and source of vapor(Dansgaard 1964;Rozanski et al.1993;Clark and Fritz 1997),which should be responsible for distribution of isotopes in river water being mainly derived from modern precipitation.The moisture for precipitation in China mainly included the southwest monsoon from the Indian Ocean,the southeast monsoon from the Pacifc Ocean,and moisture from the westerlies (Numaguti 1999).The moisture from the westerlies would produce precipitation with relatively heavy oxygen in northwest China due to signifcant local moisture recycling (Tian et al.2001)when compared with precipitation from monsoons,especially the southwest monsoon.As indicated by Fig.5,the δ18O value of weighted precipitation is generally between-9 and-6‰in most areas located in the south,east,southwest,and northwest of China.Meteoric water in the Tibetan Plateau and mountain areas of northeast and western China generally has δ18O values lower than-10‰(Luo et al.2008;Zhang and Yao 1998; Zhao et al.2012).Specifcally,the lowest δ18O values were observed in rivers in the southern Tibetan Plateau due to movement of water vapor from the southeast monsoon along the southern boundary of the Asian Continent and the southern slope of the Himalayas(Bershaw et al.2012;Yao et al.2009).

Water enriched in light isotopes in the upper reaches of the Changjiang and Yellow Rivers originated from the Tibetan Plateau.The δ18O value of the frst sample from the main channel of the Changjiang River was-13.8‰, which was lower than that of other water samples from the main channel due to the river originating from a mountain area.During most months,the δ18O value of precipitation weighted monthly for Yushu station located on the upper tributary of Changjiang River had values lower than -12‰ (Tian et al.2007).The isotopic pattern of these water samples suggested that alpine precipitation and snowmelt were the major source for the upper reach of the Changjiang River during summer,mainly refecting the altitude effect.The increasing trend was observed from different sampling periods(Fig.2)due to water contribution from tributaries with relatively high δ18O values,which is similar to the results of other studies(Dalai et al. 2002;Karim and Veizer 2002;Lambs et al.2009).There were two distinct increases in δ18O values in the main channel of the Changjiang River.The frst was due to water from other tributaries in the Sichuan Basin,such as Minjiang River(-9.9‰),Jialinjiang River(-8.8‰),and Wujiang River(-7.1‰).The latter increase was due to water from Dongting Lake and nearby tributaries,which have a large discharge in the middle reach and generally show δ18O values near-6.5‰.In the lower reach,the δ18O values of water in the main channel of the Changjiang River were stable due to there being few large tributaries.

Fig.5 The contour of δ18O-H2O values(‰)in mainland of China.a precipitation,modifed from studies(Zhang and Yao 1998;Luo et al. 2008);b river waters based on data during summer mentioned in the text

The δ18O value of the frst sample from the main channel of the Yellow River was near-11‰,which was lower than those of water samples from the main channel (Fig.3),and similar to the trend observed in the Changjiang River.The headwater in the Yellow River has higher δ18O values than that from the Changjing River,and there is a signifcantly negative relationship between precipitation δ18O and altitude in the northern Tibetan Plateau(Luo et al.2008;Zhao et al.2012).During summer,the moisture at the south of the Tibetan Plateau is affected by the southwest monsoon,while the mid-latitude westerlies provide moisture for the northern part of west China.The pattern of precipitation infuenced by altitude and mixing of airmasses should be responsible for the distribution of isotopic compositions in the Tibetan Plateau(Tian et al. 2007;Bershaw et al.2012),which lead to water in the northern Tibetan Plateau being enriched with heavy isotopes.Collectively,these effects of different δ18O of precipitation at different places might be responsible for the higher δ18O values in headwater of the Yellow River than that in the Changjiang River.The δ18O values have a relatively narrow range between-10.7 and-9.0‰ in the main channel of the Yellow River when compared with the Changjiang River.However,there was a wide variance of δ18O values in water samples collected from tributaries ranging from-10.3 to-3.0‰,which indicated that evaporation leads to large water loss in some tributaries due to the arid climate in northwest China.

There is an increasing trend of δ18O values with a range from-9.5 to-6.2‰in the main channel of the Liao River. The river water that originated in the mountains in the east of the Liao Basin had lower δ18O values than that from the west of the Liao Basin owing to the impact of the westerlies. Evaporation was likely another reason for the high δ18O values above-6‰observed in several tributaries that have low discharge of water in the river channel when sampling. However,only a slight increasing trend was observed in δ18O values in the main channel of the Songhua River during the sampling period.The water from tributaries originating in the Daxingan Mountains in the west had low δ18O values that primarily ranged from-14 to-12‰, while water from the tributaries originating from the Changbai Mountains in the east had δ18O values primarily ranging from-13 to-11‰ in the Songhua Basin.The distribution of δ18O values in these two rivers is in accordance with the difference in weighted δ18O values in precipitation based on the database for these basins(IAEA 2006),which showed that the weighted δ18O values of precipitation were-8.3‰for Liao Basin and-10.2‰for the Songhua Basin.The relatively high annual precipitation contents in the Songhua Basin were likely responsible for this pattern due to the basin being more strongly infuenced by monsoons than that in the Liao Basin.

There are few δ18O-H2O data available for the Pearl River.The water from several small streams fowing into Pearl River has δ18O values ranging from-9 to-6‰in the southern Guizhou Plateau(Li et al.2013),while the majority ranged from-7 to-6‰for river water fowing into Pearl River.The δ18O values from underground river water in Guilin in the Pearl River basin ranged from-6.6to-5.8‰ (Huang et al.2013),while the weighted precipitation δ18O values were close to-6‰ at Guilin (-6.3‰),Liuzhou(-6.4‰),and Guangzhou(-5.9‰) in the Pearl River basin(IAEA 2006).These data suggested that the δ18O values in the Pearl River might be close to -6‰,especially for water in the middle and lower reaches.Many studies based on dual isotopes of water have been conducted in rivers of southern Tibet for paleoelevation reconstruction and to determine methods of controlling factors for patterns of distribution of δ18O-H2O (Ding et al.2009;Hren et al.2009).The results of these studies have suggested that the unique pattern of δ18O-H2O values in Tibetan rivers was mainly infuenced by elevation and mixing between the southwest monsoon from the Indian Ocean and westerly moisture sources from central Asia.The Brahmaputra River originates in the Himalayas and fows across southern Tibet into India.The river has a basin area of 24.1×104km2with a mainstream river length of over 2,210 km in China.The isotopic compositions of water in the upper Brahmaputra have depleted δ18O values ranging from-20 to-15‰in southern Tibet (Ding et al.2009;Hren et al.2009).An increasing trend of δ18O values from south to north were observed with increasing altitude and δ18O values of the river water in the north of the Tibetan Plateau primarily ranged from-12 to -8‰(Ding et al.2009)due to isotopic characteristics of different air masses in Tibet.In Xinjiang Province of west China,the δ18O value was found to be close to-9‰in the Kaidu River in the south of Xinjiang Province(Liu et al.2009),while the headwater of the Urumqi River had a δ18O value of about-10‰ in the middle of Xinjiang Province(Zhang et al.2003).It has been reported that the mean value of δ18O-H2O was-6.3‰in surface water of Beijing in northern China(Zhai et al.2011),while water from the mouth of Minjiang River in southern China has a δ18O value of-6.5‰ (Zhang et al.2012).Based on the above discussion,the contour of δ18O-H2O for river water in summer could be roughly mapped(Fig.5).This map was generated using data collected from rivers across in China,excluding data from the main channels of large rivers.The spatial variation in δ18O values for river water ft the distribution pattern of δ18O of the meteoric water well regardless of limited data sources and lack of a monitoring network.In most areas of China,δ18O-H2O values in rivers mainly ranged from-9 to-7‰.The river waters enriched with light isotopes were mainly located in the Tibetan Plateau and northeast China.Water with more positive values above-7‰ in rivers was found in southern China and in the Bohai Sea surrounding northern China.However,the headwaters of the large river might be infuenced by snowmelt and lake recharge water with long term evaporation(Ding et al.2013),which would produce a slight deviation for δ18O-H2O,especially for streams derived from the Tibetan Plateau.Another possible deviation was caused by over evaporation in small tributaries in northern and western China due to the arid and semiarid climate.Of course,contouring the spatial distribution of the isotopic composition of river water would present a number of challenges(Kendall and Coplen 2001)due to possible temporal variation and different size basins,different extent of evaporation,as well as unweighted data in the large area.More studies and data sources are needed to understand spatial variation of water isotopes and sources.

4.3 Seasonal variation of isotopic composition of river water

In southern and eastern China,the precipitation generally has more negative δ18O values in summer than winter due to the summer monsoon(Zhao et al.2012).There were no signifcant differences in the isotopic values found in river water from Changjiang rivers between spring and summer throughout the fve-year sampling period(Ding et al. 2013).The lowest δ18O-H2O values were found in August in Chongqing station of the upper reaches based on semimonthly sampling in the Changjiang River(Lu et al.2012), which might have been due to the contribution of water from snowmelt and precipitation in Tibet during summer. Poor seasonal correspondence of δ18O values has been reported for many sites in the United States(Kendall and Coplen 2001).The average δ18O values of-11.4‰and δD-H2O values of-82.2‰in the dry season were similar to those for δ18O of-11.5‰and δD-H2O of-82.6‰during the wet season based on a comparison of data from the same sampling sites in the Songhua River.A possible reason for this seasonal pattern is that monsoon precipitation is the major river water supply in the monsoon domain,which could account for up to 80%of the total river water supplied annually(Yao et al.2009).In general, the low amount of precipitation should contribute less water to rivers during the dry season.This pattern differs from that observed in other studies showing that heavy isotopes in the Tien and Hau rivers were enriched during the dry season(November to April),indicating seasonal variation in the Mekong River Basin of Southeast Asia (Nguyen et al.2012).However,slightly depleted δ18OH2O values were observed in the middle reaches of the Yellow River(Fig.3),which might have been due to evaporation in summer and to that portion of the river having a different water source.In the Yellow River,the mean values of δ18O and δD-H2O in the dry season ranged from-9.4 to-66.1‰,respectively,slightly lower than the values of δ18O and δD-H2O of-9.2 and-63.7‰, respectively,during the wet season based on data from the same sampling sites.These results are similar to those of previous study(Su et al.2003)in which no distinctseasonal variation of isotopic compositions in water from the Yellow River were observed.

The δ18O values of precipitation in winter were lower than those in summer throughout the west and most of the North China due to precipitation being infuenced by the polar air mass in winter(Tian et al.2001;Zhao et al.2012). It has been reported that the values of δ18O in April to June are lower due to snowmelt,while they are higher at the end of summer due to local rainfall in the Garonne River (Lambs et al.2009).The alpine river water in Lulang in southeast Tibet was found to have more negative δ18O values between August and November,indicating that monsoon precipitation was the major source of water to this stream(Yang et al.2009).These fndings revealed that moisture entering the headwaters of the Heihe River Basin was derived predominantly from the westerly and polar air masses,leading to the more depleted δ18O values during the dry season(Zhao et al.2011).Additionally,evaporation will lead to residual waters being enriched with heavy isotopes under high temperature conditions during the wet season(Clark and Fritz 1997).During the dry season,the water had mean δ18O and δD-H2O values of-8.7 and -61.6‰,respectively,which were lower than the mean values of δ18O and δD-H2O of-7.7 and-58.9‰, respectively,that were observed during the wet season at the same sampling sites.The seasonal pattern of dual isotopes of river water samples might suggest that evaporation and different water mass infuenced the isotopic composition in the Liao River,with the former being responsible for the small difference in δD values between seasons.Moreover,long term monitoring should be conducted to enable a better understanding of temporal variations in water sources based on isotopic analysis of rivers in China.

4.4 The variation of d-excess in river water

The deuterium excess(d-excess)in precipitation is defned as the y-intercept when the slope=8,expressed by d=δD-8×δ18O(Dansgaard 1964).This equation could refect differences in humidity between vapor source and evaporation(Clark and Fritz 1997;Rozanski et al.1993) although the analyzed precision of δ18O and δD-H2O would introduce uncertainties of the calculation of d-excess values.The d-excess values in water from Changjiang River ranged from-2.2 to 10‰ with a mean value of 4.5‰,which is different from the mean of 10.6‰determined based on the isotopic data obtained from multiple sampling campaigns(Ding et al.2013),suggesting the infuence by hydrological condition under different sampling campaigns.The water in the Yellow River was found to have d-excess values ranging from-11.8 to 17.3‰with a mean value of 9.6‰,which is similar to the range of-2.9 to 19.7‰with a mean value of 8.5‰ observed for the Songhua River.The lowest mean value of 3.8‰for d-excess values was found in Liao River waters,where they ranged from-11.9 to 12.1‰.More negative dexcess values might suggest that evapotranspiration is impacting the water isotopic compositions.The high dexcess and low slope of the δ18O-δD regression line(＜8)in precipitation suggest that kinetic evapotranspiration infuenced the isotopic shift in the source region under low humidity conditions(Clark and Fritz 1997).The d-excess values were higher than the long-term average value of the Indian monsoon(8‰),which refected contributions from water vapor originating in the Mediterranean(22‰)or other inland seas(Karim and Veizer 2002).Several water samples from the tributaries of the Yellow River,the Songhua River,and the Liao River had d-excess values higher than 15,which was far above the mean value of the rivers.These fndings may indicate that the water in these regions was partly derived from vapor recycling due to a portion of these basins being impacted by semi-arid and arid climates.

As indicated in Fig.6,there were signifcant negative relationships between d-excess and δ18O values in the Yellow River(d=-3.1δ18O-19.4,R2=0.73,p＜0.001), the Songhua River (d=-2.3δ18O-19.0, R2=0.68,p＜0.001),and the Liao River (d= -3.1δ18O-20.9,R2=0.78,p＜0.001).The signifcant inverse trend between δ18O values and d-excess in the Yamuna River system suggests different mixing proportions of snow/glacial melt and precipitation(Dalai et al. 2002).In this study,water samples with d-excess values below 0 generally had high δ18O values,and samples with these values were only collected from small tributaries or tributaries discharged from lake systems.This pattern suggests that the evaporation of surface water leads to water loss and shifting isotopic compositions in these rivers.However,the relationship between d-excess and δ18O values in the Changjiang River did not show a similarpattern(d=-0.4δ18O+0.87,R2=0.10).A slight negative trend was found based on data from another study (Ding etal.2013),in which d=-1.8δ18O-5.3 (R2=0.39).These results suggest that multiple factors infuenced isotopic compositions of river water in the Changjiang Basin,such as mixing between glacial melt and different water masses from tributaries and local precipitation,as well as evaporation effects.Meanwhile,the high precision of δD analysis by laser spectroscopy analyzer (Wassenaar et al.2011)would be better to understand the variations of d-excess values and the controlling factors. Moreover,there was no clear relationship between dexcess and δ18O values for the entire dataset,which is similar to earlier results(Kendall and Coplen 2001).

Fig.6 δ18O andd-excess values in different rivers of China

There was no signifcant spatial pattern of d-excess values for precipitation based on data at GNIP sites in China(IAEA 2006).The monsoonal precipitation in summer at the Tibetan Plateau has low d-excess and dual isotopic values when compared with precipitation in winter due to westerly effects(Yao et al.2009).In stream water in the southern Tibetan Plateau,the deuterium excess values are characterized by d values above+10‰in water from below 4 km and a wide range of values from-7 to+15‰(Hren et al. 2009).The d-excess values ranged from 6.9 to 14.5‰with a mean value of 10.2‰in the Lasha River in the southern Tibetan Plateau(Yu et al.2010).Most stream water was found to have d-excess values close to 8‰based on the data from southern Guizhou province in the Pearl River Basin(Li et al.2013).Based on this rough assessment,river water generally has d-excess values below 10‰in southern and eastern China and a wide range of values in western and northern China due to evapotranspiration of air masses and evaporation from the water surface under the arid/semi-arid conditions of the region.

5 Conclusions

Water was sampled from the Changjiang,Yellow,Songhua,and Liao Rivers to evaluate spatial and temporal variation of isotopes of water as well as d-excess in rivers of China.There is a wide range of isotopic values in water in major rivers of China.Within a river,the altitude effect infuences the isotopic characteristics signifcantly,showing water from upstream was generally enriched with light isotopes.The slopes of equations describing river water values were generally slightly lower than the LMWL generated by data from GNIP sites.In the four investigated rivers,the evaporation process remains limited in the Changjiang and Songhua Rivers,and the Yellow and Liao Rivers present relatively high evaporation as deduced from the slope close to 5.The rough contour of δ18O values in summer for river water fts the distribution pattern of δ18O of the meteoric water in China regardless of limited data sources and lack of a monitoring network.The spatial pattern has showed that δ18O-H2O values in rivers mainly range from-9 to-7‰for most of China,but samples from the Tibetan Plateau and northeast China have lower values.The δ18O values in the middle and lower reaches of the Pearl River and other rivers in southern China might be enriched in heavy isotopes due to relatively high isotopic values in the rain water.Meanwhile,some rivers in northern and northwest China have relative high δ18O values due to the evaporation effect.In general,there is a wide variation of isotopic values of river waters in the summer compared with those in the winter.However,the slight distinct seasonal variation of isotopic compositions found in these rivers(aside from the Songhua)might be caused by river water mainly derived from summer monsoons across the majority of China.River water with relatively low d-excess values has more positive δ18O values in rivers of northern and western China,indicating that the evaporation in the tributaries of large rivers leads to water loss.Furthermore,long-term monitoring networks should be established and more samples should be evaluated to enable a better understanding of the water source and circulation at a regional scale,which could improve the management of water resources in major rivers of China.

AcknowledgmentsWe thank Drs.Zhang Z,Deng Y,Liu B,Hu J, and Li J for their help during sample collection.This study is fnancially supported by National Natural Science Foundation of China (Grant Nos.41210004 and 41130536)and Chinese Academy of Sciences through grants KZCX2-EW-102.

Bershaw J,Penny SM,Garzione CN(2012)Stable isotopes of modern water across the Himalaya and eastern Tibetan Plateau:implications for estimates of paleoelevation and paleoclimate.J Geophys Res.doi:10.1029/2011JD016132

Bowen GJ,Kennedy CD,Liu Z,Stalker J(2011)Water balance model for mean annual hydrogen and oxygen isotope distributions in surface waters of the contiguous United States. J Geophys Res 116:G04011.doi:10.1029/2010JG001581

Clark ID,Fritz P(1997)Environmental Isotopes in Hydrogeology. Lewis Publishers,New York

Craig H(1961)Isotopic variation in meteoric waters.Science 133:1702–1703

Criss RE(1999)Principles of stable isotope distribution.Oxford University Press,New York

Dalai TK,Bhattacharya SK,Krishnaswami S(2002)Stable isotopes in the source waters of the Yamuna and its tributaries:seasonal and altitudinal variations and relation to major cations.Hydrol Process 16:3345–3364

Dansgaard W (1964)Stable isotopesin precipitation.Tellus 16:436–468

Ding L,Xu Q,Zhang L,Yang D,Lai Q,Huang F,Shi R(2009) Regional variation of river water oxygen isotope and empirical elevation prediction models in Tibetan Plateau.Quat Sci 29:1–12 (in Chinese)

Ding T,Gao J,Shi G,Chen F,Luo X,Wang C,Hang D(2013) Spatial and temporal variations of H and O isotope compositions of the Yangtze River water and their environmental implications. Acta Geol Sinica 87:661–675(in Chinese)

Gat JR(1996)Oxygen and hydrogen isotopes in the hydrologic cycle. Annu Rev Earth Planet Sci 24:225–262

Hren MT,Bookhagen B,Blisniuk PM,Booth AL,Chamberlain CP (2009)δ18O and δD of streamwaters across the Himalaya and Tibetan Plateau:implications for moisture sources and paleoelevation reconstructions.Earth Planet.Sci.Lett.288:20–32

Huang Q,Qin X,Tang P,Liu P(2013)The characteristic and signifcance of carbon isotope(δ13CDIC)and oxygen isotope (δ18O)value in different type of karst water in Guilin. Geochimica 42:64–72(in Chinese)

IAEA(2006)Global Network of Isotopes in Precipitation.The GNIP Database.http://isohis.iaea.org

Karim A,Veizer J(2002)Water balance of the Indus River Basin and moisture source in the Karakoram and western Himalayas: implications from hydrogen and oxygen isotopes in river water. J Geophys Res.doi:10.1029/2000JD000253

Kendall C,Coplen TB (2001)Distribution of oxygen-18 and deuterium in river waters across the United States.Hydrol Process 15:1363–1393

Lambs L,Balakrishna K,Brunet F,Probst JL(2005)Oxygen and Hydrogen isotopic composition of major Indian rivers:a frst global assessment.Hydrol Process 19:3345–3355

Lambs L,Brunet F,Probst JL(2009)Isotopic characteristics of the Garonne River and its tributaries.Rapid Commun Mass Spectrom 23:2543–2550

Li SL,Liu CQ,Li J,Liu XL,Chetelat B,Wang BL,Wang FS(2010) Assessment of the sources of nitrate in the Changjiang River, China using a nitrogen and oxygen isotopic approach.Environ Sci Technol 44:1573–1578

Li SL,Liu CQ,Li J,Xue Z,Guan J,Lang Y,Ding H,Li LB(2013) Evaluation of nitrate source in surface water of southwestern China based on stable isotopes.Environ Earth Sci 68:219–228

Liu Y,Jiang G,Jin M,Li X(2009)Features of environmental isotopes of hydrogen and oxygen of water and their indication in Yanqi Basin,China.Geol Sci Technol Inf 28:89–93(in Chinese)

Lu B,Sun T,Wang C,Dai S,Kuang J,Wang J(2012)Temporal and spatial variations of δ18O along the main stem of Yangtze River, China.In:IAEA(eds)Monitoring isotopes in rivers:creation of the Global Network of Isotopes in Rivers(GNIR).Vienna, pp 211–220

Luo W,Wang S,Liu X(2008)Regional characteristics of modern precipitation δ18O values and implications for paleoclimate research in China.Earth Environ 36:47–55(in Chinese)

Michener R,Lajtha K(2007)Stable isotopes in ecology and environmental science,2nd edn.Blackwell Publishing Ltd,Oxford

Nguyen KC,Huynh L,Le DC,Nguyen VN,Tran BL(2012)Isotope composition of Mekong River fow water in south of Vietnam.In: IAEA(eds)Monitoring isotopes in rivers:creation of the Global Network of Isotopes in Rivers(GNIR).Vienna,pp 197–210

Numaguti A(1999)Origin and recycling processes of precipitating water over the Eurasian continent:experiments using an atmospheric general circulation model.J Geophys Res 104:1957–1972

Rozanski K,Araguas-Araguas L,Gonfantini R(1993)Isotope patterns in modern global precipitation.In:Swart PK et al (eds)Climate change in continental isotope records,geophysical monograph series,vol 78.AGU,Washington,pp 1–36

Su X,Lin X,Liao Z,Wang J(2003)Variation of isotopes in the Yellow River along the fow path and its affecting factors. Goechimica 32:349–357(in Chinese)

Tian L,Masson-Delmotte V,Stievenard M,Yao T,Jouzel J(2001) Tibetan Plateau summer monsoon northward extent revealed by measurements ofwaterstable isotopes.J Geophys Res 106:28,081–28088

Tian L,Yao T,MacClune K,White JWC,Schilla A,Vaughn B, Vachon R,Ichiyanagi K(2007)Stable isotopic variations in west China:a consideration of moisture sources.J Geophys Res 112:D10112.doi:10.1029/2006JD007718

Wassenaar LI,Athanasopoulos P,Hendry MJ(2011)Isotope hydrology of precipitation,surface and ground waters in the Okanagan Valley,British Columbia,Canada.J Hydrol 411:37–48

Wen R,Tian L,Weng Y,Liu Z,Zhao Z(2012)The altitude effect of δ18O in precipitation and river water in the Southern Himalayas. Chin Sci Bull 57:1693–1698

Yang X,Xu B,Yang W,Qu D,Lin P(2009)Study of altitudinal lapse rates of δ18O in precipitation/river water with seasons on the southeast Tibetan Plateau.Chin Sci Bull 54:2742–2750

Yao T,Zhou H,Yang X(2009)Indian monsoon infuences altitude effect of δ18O in precipitation/river water on the Tibetan Plateau. Chin Sci Bull 54:2724–2731

Yu T,Gan Y,Zhou A,Liu C,Liu Y,Li X,Cai H(2010) Characteristics of oxygen and hydrogen isotope distribution of surface runoff in the Lhasa River Basin.J China Univ Geosci 35:873–878(in Chinese)

Yue F-J,Li S-L,Liu C-Q,Zhao Z-Q,Hu J(2013)Using dual isotopes to evaluate sources and transformation of nitrogen in the Liao River,northeast China.Appl Geochem 36:1–9

Yue F-J,Liu C-Q,Li S-L,Zhao Z-Q,Liu X-L,Ding H,Liu B-J, Zhong J(2014)Analysis of δ15N and δ18O to identify nitrate sources and transformations in Songhua River,Northeast China. J Hydrol 519:329–339

Zhai Y,Wang J,Teng Y,Zuo R(2011)Variations of δD and δ18O in water in Beijing and their implications for the local water cycle. Resour Sci 33:92–97(in Chinese)

Zhang X,Yao T(1998)Distributional features of δ18O in precipitation in China.Acta Geograph Sin 53:356–364(in Chinese)

Zhang X,Yao T,Tian L,Liu J(2003)Stable oxygen isotope in water mediums in Urumqi River basin.Adv Water Sci 14:50–56(in Chinese)

Zhang B,Guo Z,Gao A,Yuan X,Li K,Ma Z(2012)Estimating groundwater discharge into Minjiang River estuary based on stable isotopes deuterium and oxygen-18.Adv Water Sci 23:539–548(in Chinese)

Zhao L,Yin L,Xiao H,Cheng G,Zhou M,Yang Y,Li C,Zhou J (2011)Isotopic evidence for the moisture origin and composition of surface runoff in the headwaters of the Heihe River basin. Chin Sci Bull 56:406–416

Zhao L,Xiao H,Zhou M,Cheng G,Wang L,Yin L,Ren J(2012) Factors controlling spatial and seasonal distributions of precipitation δ18O in China.Hydrol Process 26:143–152

Zhu DQ(2007)Dictionary of the Chinese River.Qingdao Press, Qingdao(In Chinese)

Received:4 April 2014/Revised:30 June 2014/Accepted:22 July 2014/Published online:23 December 2014 ©Science Press,Institute of Geochemistry,CAS and Springer-Verlag Berlin Heidelberg 2014

S.-L.Li(✉)·F.-J.Yue·C.-Q.Liu·H.Ding·Z.-Q.Zhao· X.Li

The State Key Laboratory of Environmental Geochemistry, Institute of Geochemistry,Chinese Academy of Sciences, Guiyang 550002,China

e-mail:lisiliang@vip.skleg.cn