Spatial and Temporal Variation Characteristics of Stable Isotopes in Precipitation and Their Relationships with Meteorological Factors in the Shiyang River Basin in China

: Stable isotopes of precipitation play an important role in understanding hydrological and climatic processes of arid inland river basins. In order to better understand the difference of regional water cycle and precipitation patterns, precipitation samples were collected in Shiyang River Basin from June 2018 to May 2020, and spatial and temporal variation characteristics of stable isotopes in precipitation and their relationships with meteorological factors were discussed. The results showed that stable isotopes in precipitation showed seasonal variation characteristics, that δ 2 H and δ 18 O values were higher in summer and autumn and lower in winter and spring, and d-excess values were higher in spring and autumn and lower in winter and summer. Slopes and intercepts of the local meteoric water lines gradually increased with elevation increasing. δ 18 O signiﬁcantly showed a positive correlation with temperature but a negative correlation with precipitation in summer. Contrary to below 15 hPa, δ 18 O signiﬁcantly showed a negative correlation with water vapor pressure above 15 hPa. Due to the inﬂuence of meteorological factors, there are signiﬁcant differences in stable isotopes of precipitation in arid inland river basins, which were more inﬂuenced by local moisture recycling in upstream area but by below-cloud secondary evaporation in mid-downstream area.


Introduction
Atmospheric precipitation is an important part of the water cycle process and a major source of recharge for terrestrial water resources [1,2].Stable hydrogen and oxygen isotopes ( 2 H and 18 O) are the important constituents of water molecules.During isotopic fractionation, the lighter isotopes ( 1 H and 16 O) evaporate into the gas phase first, while the heavier isotopes ( 2 H and 18 O) condense into the liquid phase first. 18O and 2 H are constantly condensed preferentially from moist air when water vapor in clouds condenses to form raindrops.Hydrogen and oxygen stable isotopes ( 2 H and 18 O) are natural tracers of the water cycle, which are susceptible to various meteorological parameters, such as atmospheric temperature, relative humidity, precipitation and water vapor pressure, and decrease or increase with environmental factor changes [2][3][4], and precipitation isotope changes have continental effects, latitude effects, temperature effects, rainfall amount effects, seasonal effects, and elevation effects [5][6][7][8][9][10][11][12][13].Temperature effects were found mainly in mid-high-latitude continents [14], and rainfall amount effects were exhibited in the low-latitude coastal, island and monsoonal humid zones [14], and elevation effects existed in mountainous areas [6].Craig [15] determined the global meteoric water line (GMWL) equation as δ 2 H = 8δ 18 O + 10, but there is a significant deviation between the Water 2023, 15, 3836 2 of 21 local meteoric water line (LMWL) and the GMWL.To quantify the extent to which the LMWL deviates from the GMWL, Dansgaard [14] proposed a deuterium excess (defined as d-excess = δ 2 H − 8δ 18 O) to quantify the extent to which the LMWL deviates from the GMWL [1].The mean value of d-excess in global atmospheric precipitation is about 10‰, and its variation depends mainly on the relative humidity and temperature of the moisture source, with lower d-excess values under the conditions of higher moisture source humidity, lower temperature (with weaker evaporation from ocean to atmosphere), and vice versa [1,16].In addition, d-excess values are affected by moisture recirculation and secondary evaporation [17,18].Therefore, the study on isotope variation in precipitation can contribute to understanding the precipitation process and water cycle mechanism.
The arid area of the northwest China is distant from the sea and is difficult for water vapor to reach, thus the climate is arid and the precipitation is scarce.Therefore, the limited water resource is crucial for regional socio-economic development.Using stable isotope techniques, many scholars have conducted numerous studies on atmospheric precipitation and its hydrological processes in the area.For example, Juan et al. [2] studied the patterns of precipitation isotopes in the Qilian Mountains and their controlling factors; Wu et al. [19] investigated the effects of different condensation mechanisms or different moisture sources on precipitation isotopes in the Heihe River Basin; Pang et al. [9] showed that the isotopic composition in precipitation under arid climate conditions was affected by sub-cloud secondary evaporation and moisture recycling; Li et al. [20] found that oases and mountainous areas formed a complete local water cycle.The above studies investigated the variation of precipitation isotopes or water cycle processes in different locations in the arid area of the northwest China, which can be a reference for other basins in this arid area.
The Shiyang River Basin is located in the northwest China and the transition zone between the monsoon area and the arid area, including three major ecosystems of mountain, oasis and desert.Influenced by global warming, the Shiyang River Basin will suffer more evaporation during precipitation, thus making the water cycle process more complex [21].It has been shown that the changes of δ 18 O and δ 2 H values in precipitation in the Shiyang River Basin were significantly influenced by temperature, elevation and sub-cloud evaporation [22].Ma et al. [23] found that the temperature was the dominant factor controlling δ 18 O values in precipitation in the Shiyang River Basin.Yuan et al. [24] pointed out that there were regional differences of stable isotopes in precipitation in the Shiyang River Basin and decreased with the increase of elevation.Gui et al. [25] found that the seasonal variation of precipitation isotopes in the upstream area of Shiyang River was significant.Zhao et al. [26] showed that the temperature effect of precipitation isotopes in the downstream area of Shiyang River was significant.The above studies have researched preliminarily the characteristics and environmental effects of precipitation isotopes in the Shiyang River Basin in short time series and small space scales, but the studies on a large spatial scale are lacking.
The previous studies were relatively short in time scale and influenced by special precipitation events, which could not comprehensively reflect the objective pattern of stable isotopes of precipitation in the study area.Spatially, they focused on certain or a few sampling points and could not fully reflect the spatial differences in stable isotope changes in precipitation in the study area.Using stable isotope and meteorological data from 11 stations in the Shiyang River Basin from June 2018 to May 2020, the spatial and temporal variation characteristics of stable isotopes in precipitation and their relationships with meteorological factors were analyzed in order to understand the process and influencing factors of stable isotope changes in precipitation in arid inland river basins.This helps to realize the difference of regional water cycle and precipitation patterns in the upstream and mid-downstream of the Shiyang River Basin, and further to provide a theoretical basis for scientific management and utilization of water resources in arid inland river basins.

Study Area
The Shiyang River Basin is located in the inland of the northwest China and on the northern slope of the Qilian Mountains (Figure 1), with a total area of about 4.16 × 10 4 km 2 .The terrain sloping from southwest to northeast is high in the south and low in the north, with the highest elevation of 5125 m and the lowest elevation of 1254 m.The Shiyang River Basin has a temperate continental arid climate with strong solar radiation and evaporation as well as a little precipitation mainly from May to October.The vapor source of precipitation is mainly controlled by the westerly and influenced by the local recirculated moisture during its transportation [24].From south to north, the Shiyang River Basin can be roughly divided into three climatic and geomorphological units, which are the alpine semi-arid and semi-humid area in the Qilian Mountains, the warm and cool arid area in the central corridor plain, and the warm and arid area in the northern desert [27].In the southern Qilian Mountains, the central corridor plain and the northern desert, the average annual temperatures of them are 0.2 • C, 6.9 • C and 8.6 • C, respectively, the annual precipitation of them are 300-600 mm, 150-300 mm and 0-150 mm, respectively, the average annual evaporation of them are 700-1200 mm, 1300-2000 mm and 2000-2600 mm, respectively [28,29].The Shiyang River originates from the northern slope of Lenglongling in the eastern Qilian Mountains and is mainly recharged by precipitation and meltwater of glaciers and snow in the mountains.The upstream of the Shiyang River flows through the Qilian Mountains and its rich precipitation becomes the recharge source of the river, the midstream and downstream areas of the region are relatively flat, with more arable land and relatively intensive agricultural activity.Because midstream and downstream areas are surrounded by the Badanjilin Desert and the Tengger Desert, they become runoff dissipation areas for arid climates [30][31][32].
Water 2023, 15, x FOR PEER REVIEW 3 of 22

Study Area
The Shiyang River Basin is located in the inland of the northwest China and on the northern slope of the Qilian Mountains (Figure 1), with a total area of about 4.16 × 10 4 km 2 .The terrain sloping from southwest to northeast is high in the south and low in the north, with the highest elevation of 5125 m and the lowest elevation of 1254 m.The Shiyang River Basin has a temperate continental arid climate with strong solar radiation and evaporation as well as a little precipitation mainly from May to October.The vapor source of precipitation is mainly controlled by the westerly and influenced by the local recirculated moisture during its transportation [24].From south to north, the Shiyang River Basin can be roughly divided into three climatic and geomorphological units, which are the alpine semi-arid and semi-humid area in the Qilian Mountains, the warm and cool arid area in the central corridor plain, and the warm and arid area in the northern desert [27].In the southern Qilian Mountains, the central corridor plain and the northern desert, the average annual temperatures of them are 0.2 °C, 6.9 °C and 8.6 °C, respectively, the annual precipitation of them are 300-600 mm, 150-300 mm and 0-150 mm, respectively, the average annual evaporation of them are 700-1200 mm, 1300-2000 mm and 2000-2600 mm, respectively [28,29].The Shiyang River originates from the northern slope of Lenglongling in the eastern Qilian Mountains and is mainly recharged by precipitation and meltwater of glaciers and snow in the mountains.The upstream of the Shiyang River flows through the Qilian Mountains and its rich precipitation becomes the recharge source of the river, the midstream and downstream areas of the region are relatively flat, with more arable land and relatively intensive agricultural activity.Because midstream and downstream areas are surrounded by the Badanjilin Desert and the Tengger Desert, they become runoff dissipation areas for arid climates [30][31][32].

Sampling Procedure
In this study, 11 sampling points of precipitation collection were laid out along the Shiyang River (Figure 1), among which four sampling points (Lenglongling, Hulinzhan, Huajianxiang and Xiyingwugou) were located in the upstream area, four points (Xiyingzhen, Yangxiaba, Jiuduntan and Hongqigu) were located in the midstream area, and three points (Xuebaizhen, Datanxiang and Qingtuhu) were located in the downstream area.The terrain of the midstream and downstream areas is relatively flat, but due to being surrounded by the Badanjilin Desert and the Tengger Desert, they become runoff dissipation areas for arid climates, as well as irrigation, cropland and wind-sand areas

Sampling Procedure
In this study, 11 sampling points of precipitation collection were laid out along the Shiyang River (Figure 1), among which four sampling points (Lenglongling, Hulinzhan, Huajianxiang and Xiyingwugou) were located in the upstream area, four points (Xiyingzhen, Yangxiaba, Jiuduntan and Hongqigu) were located in the midstream area, and three points (Xuebaizhen, Datanxiang and Qingtuhu) were located in the downstream area.The terrain of the midstream and downstream areas is relatively flat, but due to being surrounded by the Badanjilin Desert and the Tengger Desert, they become runoff dissipation areas for arid climates, as well as irrigation, cropland and wind-sand areas [30,31].Because of their similar climatic and geomorphologic features [32], they are subject to similar local geographic environments.Therefore, they are referred to as the mid-downstream area and discussed as a whole.A total of 688 samples were collected during the two hydrological years from June 2018 to May 2020 (Table 1).Precipitation samples were collected using the standard rainfall barrel of the weather station, which was placed on 1.5 m high above the ground and without obstacles around.A precipitation event was considered from 20:00 on the day to 20:00 on the next day, and the water samples were filled into high-density polyethylene (HDPE) bottles and sealed with waterproof tape.Solid precipitation samples were sealed in low-density polyethylene (LDPE) bags and placed in a room to melt completely at room temperature before they were filled and sealed.During precipitation sample collection, automatic meteorological observation stations record meteorological parameters simultaneously and include temperature, dew point temperature, precipitation, relative humidity and barometric pressure.

Methods
The precipitation samples were tested in the stable isotope laboratory of Northwest Normal University with the LGR-DLT-100 liquid water isotope analyzer (DLT-100, Los Gatos Research Inc., Mountain View, CA, USA).The analytical precision of δ 2 H and δ 18 O are ±0.6‰ and ±0.2‰, respectively.The results of the sample analysis are expressed as the thousandths of the Vienna Standard Mean Ocean Water (V-SMOW) equivalent: where R s is the ratio of 18 where δ i is the isotope value and p i is the respective precipitation amount.
Craig [15] determined the global meteoric water line (GMWL) equation as: Dansgaard [14] defined excess deuterium (d-excess) based on isotope data of global precipitation.It can be calculated as: An ordinary least squares regression (OLSR) is widely used to determine the LMWL, and this regression method logically gives all data points equal weighting, which may be sensitive for some heavy and small precipitation events [33].In order to reduce the bias of small rainfall events, Hughes and Crawford [33] introduced the Precipitation Weighted Least Squares Regression (PWLSR).In our study, the PWLSR method was used to calculate LWML, and data from daily precipitation events were used to calculate LMWL (Including the LMWL of the upstream areas, Upstream, Mid-downstream and each station).The event-based datasets provide a wider range of values to better constrain the regression than that achieved using monthly data over a short period.
The Lifting Condensation Level (LCL) can be calculated as [34]: where is the dry adiabatic lapse rate, and Γ dew ( • C/m) is the wet adiabatic lapse rate.

Daily Variations
For all precipitation events, the daily variation of δ 18 O values ranged from −31.5 to 5.8‰, with a weighted-average value of −8.0‰, and the variation of δ 2 H values ranged from −238.6 to 38.6‰, with a weighted-average value of −51.5‰ (Figure 2). Figure 2 shows the variation of δ 18 O values over time at each sampling point, among sample data of three sampling points (Xiyingzhen, Yangxiaba and Xuebaizhen) in the mid-downstream area mainly in the first hydrologic year.The daily variation ranges of δ 18 O values were larger in the upstream area, among which it was the largest at Lenglongling with the largest difference of 34.6‰.While they were smaller in the mid-downstream area, among which it was the smallest at Datanxiang with a difference of 9.8‰.The highest and lowest values of δ 18 O were 3.1‰ (Lenglongling, 12 July 2018) and −31.5‰ (Lenglongling, 19 January 2019) respectively in the upstream area, and they were 5.8‰ (Xuebaizhen, 6 June 2018) and −25.7‰ (Jiuduntan, 10 January 2020) respectively in the mid-downstream area.In addition, the δ 18 O values were also lower at Xiyingwugou, which may be related to the local microclimate influenced by the Xiying reservoir nearby [35].The daily variations of δ 18 O and δ 2 H values had similarities, and all sampling points showed the characteristics of higher in the warm season (May-October) and lower in the cold season (November-April).For example, the δ 18  .For all precipitation events, 62.7% of d-values were greater than 10‰.In addition, 11.3% of the precipitation events had negative d-values, most of which occurred in summer.This may be related to the strong sub-cloud secondary evaporation effect caused by high temperatures.
occurred in summer.This may be related to the strong sub-cloud secondary evaporation effect caused by high temperatures.
The daily variation trend of δ 18 O values and d-values in the mid-downstream area was the opposite, while it was unobvious in the upstream area.This is because the belowcloud secondary evaporation is strong in the mid-downstream area, which makes δ 18 O values rise and d-values decrease, and further makes δ 18 O values and d-values show a reverse trend.However, as can be seen in Figure 3, the lifting condensation level (LCL) in the upstream region is lower than that in the middle-downstream region, so the distance of raindrop landing is shorter and the sub-cloud secondary evaporation is lower.The daily variation trend of δ 18 O values and d-values in the mid-downstream area was the opposite, while it was unobvious in the upstream area.This is because the below-cloud secondary evaporation is strong in the mid-downstream area, which makes δ 18 O values rise and d-values decrease, and further makes δ 18 O values and d-values show a reverse trend.However, as can be seen in Figure 3, the lifting condensation level (LCL) in the upstream region is lower than that in the middle-downstream region, so the distance of raindrop landing is shorter and the sub-cloud secondary evaporation is lower.

Monthly and Seasonal Variations
On the monthly scale (Figure 4), the maximum values of δ 2 H and δ 18 O values occurred in June and the minimum values occurred in January.δ 2 H and δ 18 O values showed consistent trends, maintaining high values from June to September, a decreasing trend from September to January and an increasing trend from January to June.On the seasonal scale, the weighted-average values of δ 2 H were −64.0‰, −39.3‰, −51.8‰ and −164.9‰ in spring, summer, autumn and winter respectively, and those of δ 18 O were −9.9‰, −6.3‰, −8.5‰ and −21.4‰ respectively.The seasonal variations of δ 2 H and δ 18 O values showed that summer > autumn > spring > winter.It can be seen that the seasonal variation characteristics of stable isotopes in precipitation were higher in summer and autumn but lower in winter and spring.The low values of δ 2 H and δ 18 O were mostly found in winter because of the low temperature and the depletion of moisture after long-distance transport [16].In the study area, moisture in winter is mainly transported by westerly, which originated from the Atlantic passes through Central Asia and then encounters mountains to form precipitation, leading to stable isotope depletion through Rayleigh fractionation during longer distance transportation [36].As the temperature increases in spring, the values of δ 2 H and δ 18 O gradually increase.This is because the higher temperature during condensation influences equilibrium fractionation between water vapor and raindrops, and raindrops are influenced by sub-cloud secondary evaporation during their fall to the Earth's surface, making the values of δ 2 H and δ 18 O gradually increase and appear high values in summer.The opposite was true after autumn when the decreasing temperature led to the values of δ 2 H and δ 18 O gradually decreasing.In summer, local moisture is continuously infiltrated along the way during transportation, which compensates for the depletion of the isotopes.
The maximum value of d-excess occurred in October and the minimum value occurred in January.The values of d-excess increased from June to October, decreased from October to January and increased again from January to May.The weighted-average values of d-excess were 14.8‰, 11.1‰, 16.2‰ and 6.5‰ in spring, summer, autumn and winter respectively.The seasonal variation of d-values showed autumn > spring > summer > winter and showed the characteristics of higher in spring and autumn and lower in winter and summer.Overall, spring and summer are primarily influenced by the moisture source area but autumn and winter are primarily influenced by secondary processes.Zhang et al. [36] pointed out that as the lower relative humidity of the westerly moisture source area causes higher d-values in spring, but influenced by the monsoon moisture, the

Monthly and Seasonal Variations
On the monthly scale (Figure 4), the maximum values of δ 2 H and δ 18 O values occurred in June and the minimum values occurred in January.δ 2 H and δ 18 O values showed consistent trends, maintaining high values from June to September, a decreasing trend from September to January and an increasing trend from January to June.On the seasonal scale, the weighted-average values of δ 2 H were −64.0‰, −39.3‰, −51.8‰ and −164.9‰ in spring, summer, autumn and winter respectively, and those of δ 18 O were −9.9‰, −6.3‰, −8.5‰ and −21.4‰ respectively.The seasonal variations of δ 2 H and δ 18 O values showed that summer > autumn > spring > winter.It can be seen that the seasonal variation characteristics of stable isotopes in precipitation were higher in summer and autumn but lower in winter and spring.The low values of δ 2 H and δ 18 O were mostly found in winter because of the low temperature and the depletion of moisture after long-distance transport [16].In the study area, moisture in winter is mainly transported by westerly, which originated from the Atlantic passes through Central Asia and then encounters mountains to form precipitation, leading to stable isotope depletion through Rayleigh fractionation during longer distance transportation [36].As the temperature increases in spring, the values of δ 2 H and δ 18 O gradually increase.This is because the higher temperature during condensation influences equilibrium fractionation between water vapor and raindrops, and raindrops are influenced by sub-cloud secondary evaporation during their fall to the Earth's surface, making the values of δ 2 H and δ 18 O gradually increase and appear high values in summer.The opposite was true after autumn when the decreasing temperature led to the values of δ 2 H and δ 18 O gradually decreasing.In summer, local moisture is continuously infiltrated along the way during transportation, which compensates for the depletion of the isotopes.
The maximum value of d-excess occurred in October and the minimum value occurred in January.The values of d-excess increased from June to October, decreased from October to January and increased again from January to May.The weighted-average values of d-excess were 14.8‰, 11.1‰, 16.2‰ and 6.5‰ in spring, summer, autumn and winter respectively.The seasonal variation of d-values showed autumn > spring > summer > winter and showed the characteristics of higher in spring and autumn and lower in winter and summer.Overall, spring and summer are primarily influenced by the moisture source area but autumn and winter are primarily influenced by secondary processes.Zhang et al. [36] pointed out that as the lower relative humidity of the westerly moisture source area causes higher d-values in spring, but influenced by the monsoon moisture, the d-values are lower in summer, however, the higher contribution of local recycling moisture causes

Interannual Variations
As stable isotope data of four sampling stations at Lenglongling, Hulinzhan, Xiyingwugou and Jiuduntan from June 2018 to May 2020 were continuous, we used data from these four stations to compare the characterization of precipitation isotopes for the two hydrologic years in June 2018/May 2019 and June 2019/May 2020 (Three of these sites, Lenglongling, Hulinzhan and Xiyingwugou, are located upstream, therefore yield results are more consistent with the characteristics of the upstream area).From the inter-annual trends (Figures 2l and 5), δ 2 H, δ 18 O and d-values all showed regular and periodic seasonal variations, i.e., δ 2 H and δ 18 O values were higher in summer and autumn and lower in winter and spring, while d-values were higher in spring and autumn and lower in winter and summer.However, inter-annual variations of δ 2 H and δ 18 O values were different because the seasonal fluctuations were stronger in the hydrological year of June 2018/May 2019 than in the hydrological year of June 2019/May 2020.The δ 2 H and δ 18 O values appeared a rapidly decreasing trend for the former but a slowly and steadily decreasing trend for the latter from June to August.As seen in Figure 5, the precipitation was more in 2018 than in 2019 and it was month-by-month increasing from June to August in 2018, which showed that it had a more significant precipitation effect on δ 2 H and δ 18

Interannual Variations
As stable isotope data of four sampling stations at Lenglongling, Hulinzhan, Xiyingwugou and Jiuduntan from June 2018 to May 2020 were continuous, we used data from these four stations to compare the characterization of precipitation isotopes for the two hydrologic years in June 2018/May 2019 and June 2019/May 2020 (Three of these sites, Lenglongling, Hulinzhan and Xiyingwugou, are located upstream, therefore yield results are more consistent with the characteristics of the upstream area).From the inter-annual trends (Figures 2l and 5), δ 2 H, δ 18 O and d-values all showed regular and periodic seasonal variations, i.e., δ 2 H and δ 18 O values were higher in summer and autumn and lower in winter and spring, while d-values were higher in spring and autumn and lower in winter and summer.However, inter-annual variations of δ 2 H and δ 18 O values were different because the seasonal fluctuations were stronger in the hydrological year of June 2018/May 2019 than in the hydrological year of June 2019/May 2020.The δ 2 H and δ 18 O values appeared a rapidly decreasing trend for the former but a slowly and steadily decreasing trend for the latter from June to August.As seen in Figure 5, the precipitation was more in 2018 than in 2019 and it was month-by-month increasing from June to August in 2018, which showed that it had a more significant precipitation effect on δ 2 H and δ 18

Spatial Variations of Stable Isotopes
Divided into all year, warm season (May-October), and cold season (November-April), the spatial variations of stable isotopes in the Shiyang River Basin were drawn (Figure 6).The high values of δ 18 O and δ 2 H in all year were mainly distributed in the middownstream area, and the low values of them were distributed in the upstream area, reflecting that δ 2 H and δ 18 O values decreased with the elevation increasing (Figure 6a,d 6‰ respectively.Compared to the hydrological year of June 2019/May 2020, though there were higher mean temperatures and fewer precipitation events, the more precipitation in the hydrological year of June 2018/May 2019 caused the isotope values to have smaller variation ranges and lower mean values, which indicated the precipitation had a significant effect on the variation of stable isotopes.

Spatial Variations of Stable Isotopes
Divided into all year, warm season (May-October), and cold season (November-April), the spatial variations of stable isotopes in the Shiyang River Basin were drawn (Figure 6).The high values of δ 18 O and δ 2 H in all year were mainly distributed in the mid-downstream area, and the low values of them were distributed in the upstream area, reflecting that δ 2 H and δ 18 O values decreased with the elevation increasing (Figure 6a,d).The spatial variation of d-values was opposite to δ 18 O and δ 2 H values (Figure 6g).During the warm season (Figure 6b,e,h) and cold season (Figure 6c,f,i), the spatial variations of δ 2 H, δ 18 O and dvalues were consistent with all year.It has been established that the stable isotope decreases with the elevation increasing when water vapor lifts along a certain slope [37], which exhibits the elevation effect.In all year, warm season and cold season, the elevation effects of δ 2 H in the Shiyang River Basin were −1.19‰/100 m, −1.55‰/100 m and −0.6‰/100 m respectively, and those of δ 18 O were −0.20‰/100 m, −0.23‰/100 m and −0.12‰/100 m respectively, and those of d-excess were 0.38‰/100 m, 0.28‰/100 m and 0.33‰/100 m respectively.The elevation effect of stable isotope in precipitation is actually a reflection of the temperature effect, where the temperature gradually decreases with elevation increasing and the stable isotope in precipitation also gradually decreases [6,22].Guo et al. [38] concluded that the mean δ 18 O value in precipitation in Guilin City was −5.6‰.Zhang et al. [39] concluded that the mean values of δ 18 O in precipitation in the Tsaidam-Qinghai Lake region and Junggar-Tuha region were −8.7‰ and −9.1‰, respectively.In this study, the mean value of δ 18 O in precipitation in the Shiyang River basin was −8.0‰, which was lower than that of the coastal area and higher than that of the inland areas, reflecting the gradual depletion of stable isotopes of precipitation from the coast to the inland.In the distribution map of precipitation stable isotope contours in China drawn by Kong et al. [40], δ 18 O and δ 2 H show a general trend of depletion from south to north and from east to west, of which the contour lines in the Shiyang River basin area show a decreasing trend from northeast to southwest.This is consistent with the decreasing trend of δ 18 O and δ 2 H with increasing altitude in this study.
Water 2023, 15, x FOR PEER REVIEW 10 of 22 −0.12‰/100 m respectively, and those of d-excess were 0.38‰/100 m, 0.28‰/100 m and 0.33‰/100 m respectively.The elevation effect of stable isotope in precipitation is actually a reflection of the temperature effect, where the temperature gradually decreases with elevation increasing and the stable isotope in precipitation also gradually decreases [6,22].Guo et al. [38] concluded that the mean δ 18 O value in precipitation in Guilin City was −5.6‰.Zhang et al. [39] concluded that the mean values of δ 18 O in precipitation in the Tsaidam-Qinghai Lake region and Junggar-Tuha region were −8.7‰ and −9.1‰, respectively.In this study, the mean value of δ 18 O in precipitation in the Shiyang River basin was −8.0‰, which was lower than that of the coastal area and higher than that of the inland areas, reflecting the gradual depletion of stable isotopes of precipitation from the coast to the inland.In the distribution map of precipitation stable isotope contours in China drawn by Kong et al. [40], δ 18 O and δ 2 H show a general trend of depletion from south to north and from east to west, of which the contour lines in the Shiyang River basin area show a decreasing trend from northeast to southwest.This is consistent with the decreasing trend of δ 18 O and δ 2 H with increasing altitude in this study.
As shown in Table 2, the slope of LMWL at Lenglongling was the highest in the whole basin, and that was the lowest at Xuebaizhen.It can be seen that the slope of LMWL increased with the elevation increasing, which was related to the decreasing of temperature and the weakening of the evaporation effect.The intercept of LMWL also showed an increasing trend with the elevation increasing.The intercept of LMWL was the highest at Lenglongling and the lowest at Qingtuhu.The intercept of LMWL was related to subcloud secondary evaporation and local moisture circulation.As the elevation increased, the temperature and evaporation gradually decreased while the precipitation, relative humidity and local moisture circulation gradually increased, which led to the slope and intercept of LMWL gradually increasing correspondingly.Overall, the slope and intercept of LMWL were higher in the upstream area than in the mid-downstream area.In the arid region of the northwest China, the slope and intercept of LMWL are generally low at the river downstream and generally high at the river upstream [43,44], which is consistent with the results of this study.It has been shown that local recycling moisture can lead to them increase, but the effect of sub-cloud secondary evaporation can cause them to decrease [45][46][47][48].The more rainfall and low temperature in the upstream area make δ 2 H values more enriched than δ 18 O values in the precipitation process [43], and the d-excess value increases significantly, further leading to the higher slope and intercept of LMWL [49], however, little precipitation, high temperature, strong evaporation and low relative humidity in the mid-downstream area led to them lower [43], furthermore, topography and subsurface were also important factors affecting them in the arid area of central Asia [50].In the upstream area, the higher elevation results in a lower temperature, weaker evaporation and more rainfall because of topographic uplift, and the higher vegetation cover causes a higher slope and intercept of the LMWL [43,49].In the mid-downstream area, the lower elevation, resulting in a higher temperature and stronger sub-cloud secondary evaporation, and the lower vegetation cover cause the lower slope and intercept of LMWL.Because of the good correlation between δ 18 O and δ 2 H, δ 18 O was used mainly for the analysis.The precipitation was divided into three groups according to the temperature: T < 0 °C, 0-10 °C and T > 10 °C (Table 3, Figure 8).From Figure 8a and Table 3, it can be seen that δ 18 O showed an increasing trend with the temperature increasing.Especially, when the temperature increased from 0 °C to 10 °C, the increasing extent of δ 18 O with temperature and the correlation coefficient r was the largest.The d-excess showed a decreasing trend when T > 0 °C, and the decreasing extent of d-excess was the largest when T > 10 °C.
At the whole basin scale, when T < 0 °C, the temperature effect was 0.52‰/°C, and the d-excess was positively correlated with T. At this temperature, precipitation was usually snowfall, and stable isotope changes were influenced by equilibrium fractionation [2,27,52].When 0 °C < T < 10 °C, the temperature effect was 0.77‰/°C, but the d-excess showed a weak negative correlation with T. At this temperature, the temperature effect was the most significant and the stable isotope concentration of precipitation was influenced by below-cloud secondary evaporation.When the temperature was greater than 10 °C, the temperature effect was not obvious, but the d-excess showed an obvious negative correlation with T and its temperature effect was −0.39‰/°C.The temperature in the Shiyang River Basin is relatively high in the summer half-year when the precipitation is relatively concentrated.The effect of sub-cloud secondary evaporation will be weakened with the increase in precipitation [53].At the same time, the higher temperature-induced isotope enrichment and the circulating water vapor-induced isotope depletion make the Because of the good correlation between δ 18 O and δ 2 H, δ 18 O was used mainly for the analysis.The precipitation was divided into three groups according to the temperature: T < 0 • C, 0-10 • C and T > 10 • C (Table 3, Figure 8).From Figure 8a and Table 3, it can be seen that δ 18 O showed an increasing trend with the temperature increasing.Especially, when the temperature increased from 0 • C to 10 • C, the increasing extent of δ 18 O with temperature and the correlation coefficient r was the largest.The d-excess showed a decreasing trend when T > 0 • C, and the decreasing extent of d-excess was the largest when T > 10 • C.  downstream area than in the upstream area.When T > 10 °C, the temperature effect weakened in the mid-downstream area, while δ 18 O and temperature showed a weak negative correlation in the upstream area.At this temperature, the negative correlation between dexcess and T was obvious in the mid-downstream area.This was because the stable isotopes in precipitation were more affected by sub-cloud secondary evaporation in the middownstream area, while they were more influenced by local recirculated moisture in the upstream area [43].At the whole basin scale, when T < 0 • C, the temperature effect was 0.52‰/ • C, and the d-excess was positively correlated with T. At this temperature, precipitation was usually snowfall, and stable isotope changes were influenced by equilibrium fractionation [2,27,52].When 0 • C < T < 10 • C, the temperature effect was 0.77‰/ • C, but the d-excess showed a weak negative correlation with T. At this temperature, the temperature effect was the most significant and the stable isotope concentration of precipitation was influenced by below-cloud secondary evaporation.When the temperature was greater than 10 • C, the temperature effect was not obvious, but the d-excess showed an obvious negative correlation with T and its temperature effect was −0.39‰/ • C. The temperature in the Shiyang River Basin is relatively high in the summer half-year when the precipitation is relatively concentrated.The effect of sub-cloud secondary evaporation will be weakened with the increase in precipitation [53].At the same time, the higher temperature-induced isotope enrichment and the circulating water vapor-induced isotope depletion make the variation of δ 18 O with T in a relatively stable state [54].Thus, the r between δ 18 O and T is small at this temperature.Influenced by monsoon moisture in the summer half-year [55], there was an obvious negative correlation between d-excess and T.
The temperature effects showed regional differences in different temperature conditions (Table 3, Figure 8).When T < 0 • C, the temperature effect was obvious in the upstream area.When 0 • C < T < 10 • C, the temperature effect was more significant in the mid-downstream area than in the upstream area.When T > 10 • C, the temperature effect weakened in the mid-downstream area, while δ 18 O and temperature showed a weak negative correlation in the upstream area.At this temperature, the negative correlation between d-excess and T was obvious in the mid-downstream area.This was because the stable isotopes in precipitation were more affected by sub-cloud secondary evaporation in the mid-downstream area, while they were more influenced by local recirculated moisture in the upstream area [43].

Relationship between Stable Isotopes in Precipitation and Precipitation Amount
At the daily scale, the linear equations of δ 2 H, δ 18 O and d-excess versus precipitation P (mm) were obtained as follows: δ 2 H = 1.06P − 62.31 (R 2 = 0.01, p < 0.01), δ 18 O = 0.08P − 8.85 (R 2 = 0.003), d-excess = 0.42P + 8.47 (R 2 = 0.04, p < 0.01).For all precipitation events, δ 2 H and δ 18 O showed positive correlations with precipitation, which showed that the stable isotopes of precipitation in the Shiyang River Basin didn't show the rainfall amount effect on the daily scale.Under the seasonal scale (Table 4), δ 2 H and δ 18 O showed a negative correlation with precipitation in summer (passed the 0.01 significance test), indicating that the variation of δ 2 H and δ 18 O showed the rainfall amount effect in summer.This was because the precipitation in the study area was mainly concentrated in summer when it was influenced by monsoon moisture [55], thus the rainfall amount effect was stronger.At the daily scale and the seasonal scale (except in winter), the rainfall amount effect between the d-excess with precipitation was stronger, which indicated the variation of d-excess was influenced by that of the precipitation.

Number of Samples
The Amount Effect of The Amount Effect of The Amount Effect of d-Excess (‰/mm) (R The rainfall amount effect varied in different regions.In summer 2019, the rainfall amount effect of δ 18 O was −0.27‰/mm (r = −0.38) in the upstream area (Figure 9a), and it was −0.25‰/mm (r = −0.26) in the mid-downstream area (Figure 9b).The negative correlation between δ 18 O and P was greater in the upstream area than in the mid-downstream area, indicating that the rainfall amount effect of δ 18 O was stronger in the former than in the latter.During the same period, the temperature effect of δ 18 O was 0.26‰/ • C (r = 0.24) in the upstream area, and it was 0.65‰/ • C (r = 0.45) in the mid-downstream area, indicating that the temperature effect of δ 18 O was greater in the latter than in the former.In the mid-downstream area, the raindrops were influenced by stronger sub-cloud secondary evaporation during the descent of the raindrops, which caused the temperature effect to mask the rainfall amount effect.It also has occurred in Southeast Asia [13].The d-excess showed a significant positive correlation with precipitation, which was consistent with the results in the Tolai River Basin [56].Under arid climate conditions, local moisture recirculation formed by evaporation and higher relative humidity brought by precipitation amount increases making the d-excess higher [20], which leads to the d-excess showing a positive correlation with precipitation.Additionally, the increased extent of d-excess with precipitation was smaller in summer than in spring and autumn, which was related to the decreasing of d-excess influenced by monsoon moisture [55].According to precipitation class (Precipitation classes are classified by the amount of precipitation accumulated in 24 h), precipitation can be categorized into three types: light rain (0.1-9.9 mm), moderate rain (10.0-24.9mm) and heavy rain (25.0-49.9mm).Since there were only two samples of heavy rain, the discussion was based on light and moderate rain (Table 5).When the rain was 0.1-9.9mm, the δ 18 O of the whole basin was positively correlated with P (r = 0.05), indicating there wasn't a rainfall amount effect between δ 18 O with P. When the rain increased to 10.0-24.9mm, the δ 18 O of the whole basin was negatively correlated with P (r = −0.11),indicating δ 18 O decreased with precipitation increasing.Changing from light rain to moderate rain, the rainfall amount effect of δ 18 O became stronger.As the rainfall amount was below 10 mm, the d-excess increased with the rainfall amount increasing.Most of the precipitation less than 10 mm appeared in spring or autumn, when the temperature gradually rose (spring) or was still high (autumn), but the monsoon had not arrived (spring) or had withdrawn (autumn), which caused stable isotopes were more influenced by temperature.Some precipitation of less than 10 mm appeared in summer, the high temperature and strong evaporation made isotope enrichment for below-cloud secondary evaporation [46,47,[57][58][59], which caused stable isotopes to be more affected by temperature.When the precipitation was greater than 10 mm, the relative humidity rose with the increase of precipitation, and stable isotopes were caused of depletion influenced by precipitation increasing.A similar phenomenon has been found that stable isotopes in smaller precipitation (1-10 mm) had a stronger temperature effect [60], and the higher relative humidity was prone to isotope dilution effects in precipitation greater than 10 mm [1].
Table 5. Relationship between stable isotopes of precipitation and precipitation amount under different rainfall conditions in the Shiyang River basin.

Area
The Amount Effect of The Amount Effect of The  According to precipitation class (Precipitation classes are classified by the amount of precipitation accumulated in 24 h), precipitation can be categorized into three types: light rain (0.1-9.9 mm), moderate rain (10.0-24.9mm) and heavy rain (25.0-49.9mm).Since there were only two samples of heavy rain, the discussion was based on light and moderate rain (Table 5).When the rain was 0.1-9.9mm, the δ 18 O of the whole basin was positively correlated with P (r = 0.05), indicating there wasn't a rainfall amount effect between δ 18 O with P. When the rain increased to 10.0-24.9mm, the δ 18 O of the whole basin was negatively correlated with P (r = −0.11),indicating δ 18 O decreased with precipitation increasing.Changing from light rain to moderate rain, the rainfall amount effect of δ 18 O became stronger.As the rainfall amount was below 10 mm, the d-excess increased with the rainfall amount increasing.Most of the precipitation less than 10 mm appeared in spring or autumn, when the temperature gradually rose (spring) or was still high (autumn), but the monsoon had not arrived (spring) or had withdrawn (autumn), which caused stable isotopes were more influenced by temperature.Some precipitation of less than 10 mm appeared in summer, the high temperature and strong evaporation made isotope enrichment for below-cloud secondary evaporation [46,47,[57][58][59], which caused stable isotopes to be more affected by temperature.When the precipitation was greater than 10 mm, the relative humidity rose with the increase of precipitation, and stable isotopes were caused of depletion influenced by precipitation increasing.A similar phenomenon has been found that stable isotopes in smaller precipitation (1-10 mm) had a stronger temperature effect [60], and the higher relative humidity was prone to isotope dilution effects in precipitation greater than 10 mm [1].

Relationship between Stable Isotopes in Precipitation and Relative Humidity
The linear equations of δ 18 O, δ 2 H and d-excess with the relative humidity RH in the Shiyang River Basin were calculated as: δ 18 O = 0.03H − 10.39 (R 2 = 0.004, n = 673); δ 2 H = 0.42H − 84.75 (R 2 = 0.01, p < 0.01); d-excess = 0.19H − 1.64 (R 2 = 0.08, p < 0.01).δ 18 O, δ 2 H and d-excess all increased with the RH increasing, but there was a weak positive Water 2023, 15, 3836 16 of 21 correlation between δ 18 O with the RH.To further analyze the effect of RH on stable isotopes, the RH was classified as 0-50%, 50-60%, 60-70%, 70-80%, and 80-100% (Table 6, Figure 10).correlation coefficient r decreased from 0.17 to 0.09).When the RH was 80-100%, δ 18 O decreased with the increase of RH and showed a weak negative correlation.Chen et al. [61] found a negative correlation between δ 18 O and the RH in the northwest China, which was also found in the Tibetan Plateau [62] and Bangladesh [63].In the range of 0-80%, the RH shows a positive correlation with δ 18 O, which was mainly influenced by the temperature effect.When the RH changed from 50% to 100%, the slope of δ 18 O versus RH decreased continuously, indicating that the increase of the RH weakened the temperature effect.This was because the below-cloud secondary evaporation effect gradually decreased with the precipitation amount and the RH increasing [27].In addition, the isotope dilution effect occurred when the RH was high and the precipitation was formed [1].As can be seen in Figure 10, the d-excess varied little with the RH below 50%, while the dexcess decreased with the RH increasing at 50-60%.Zhao et al. [58] found that the RH less than 60% could lead to a decrease of d-excess but lead to an increase of d-excess when the RH was greater than 60%.Crawford et al. [64] found a positive correlation between the RH and the d-excess in the western side of the Great Dividing Range, which was also found in the Shule River Basin [45] and the Binggou River Basin [44] in the Qilian Mountains.For the RH in the ranges of 0-50%, 50-60%, 60-70%, 70-80% and 80-100%, the weighted-average values of d-excess were 8.9‰, 8.2‰, 11.9‰, 13.8‰ and 15.4‰ respectively.This indicated that the d-values gradually increased with the RH increasing.The RH is an important indicator to control evaporation.Under low RH conditions, the subcloud secondary evaporation would be significantly enhanced [12], which caused the values of d-excess to decrease and those of δ 18 O to increase, it is vice versa.Combined with Table 6 and Figure 10, when the RH was below 50%, δ 18 O increased with the increase of RH, but the correlation coefficient was very small.When the RH changed from 50% to 80%, δ 18 O increased smaller and smaller with the increase of RH (the correlation coefficient r decreased from 0.17 to 0.09).When the RH was 80-100%, δ 18 O decreased with the increase of RH and showed a weak negative correlation.Chen et al. [61] found a negative correlation between δ 18 O and the RH in the northwest China, which was also found in the Tibetan Plateau [62] and Bangladesh [63].In the range of 0-80%, the RH shows a positive correlation with δ 18 O, which was mainly influenced by the temperature effect.When the RH changed from 50% to 100%, the slope of δ 18 O versus RH decreased continuously, indicating that the increase of the RH weakened the temperature effect.This was because the below-cloud secondary evaporation effect gradually decreased with the precipitation amount and the RH increasing [27].In addition, the isotope dilution effect occurred when the RH was high and the precipitation was formed [1].As can be seen in Figure 10, the d-excess varied little with the RH below 50%, while the d-excess decreased with the RH increasing at 50-60%.Zhao et al. [58] found that the RH less than 60% could lead to a decrease of d-excess but lead to an increase of d-excess when the RH was greater than 60%.Crawford et al. [64] found a positive correlation between the RH and the d-excess in the western side of the Great Dividing Range, which was also found in the Shule River Basin [45] and the Binggou River Basin [44] in the Qilian Mountains.For the RH in the ranges of 0-50%, 50-60%, 60-70%, 70-80% and 80-100%, the weighted-average values of d-excess were 8.9‰, 8.2‰, 11.9‰, 13.8‰ and 15.4‰ respectively.This indicated that the d-values gradually increased with the RH increasing.The RH is an important indicator to control evaporation.Under low RH conditions, the sub-cloud secondary evaporation would be significantly enhanced [12], which caused the values of d-excess to decrease and those of δ 18 O to increase, it is vice versa.

Relationship between Stable Isotopes in Precipitation and Water Vapor Pressure
The linear equations of δ 18 O, δ 2 H and d-excess with water vapor pressure e in Shiyang River Basin were calculated as: δ 18 O = 0.71e − 15.16 (R 2 = 0.33, n = 673, p < 0.01); δ 2 H = 5.48e − 107.58 (R 2 = 0.32, p < 0.01); d-excess = −0.22e+ 13.72 (R 2 = 0.02, p < 0.01).The δ 18 O and δ 2 H with the water vapor pressure showed a significant positive correlation, and this was opposite to d-excess.In the arid region of the northwest China, stable isotopes in precipitation are influenced by the water vapor pressure.For example, there was a positive correlation between δ 18 O and water vapor pressure in the central part of the Hexi Corridor [65] and in the Tolai River Basin [66], and a negative correlation between d-excess and water vapor pressure, which was consistent with the results of this study.
According to Table 7 and Figure 11, when the water vapor pressure was in the range of 0-5 hpa, δ 18 O and δ 2 H had positive correlations with it and the increasing extent of them was the largest, and the d-excess increased with its increasing (r = 0.36).At this time, the effects of sub-cloud secondary evaporation and local moisture recycling were weak because the average T was below 0 • C and the initial d-values were high in dry air masses influenced by the moisture source [66].When the water vapor pressure changed to 5-10 hPa, the increase extents of δ 18 O and δ 2 H decreased rapidly and the d-excess decreased with its increasing (r = −0.03).At this moment, the average T was 8.4 • C, and the RH was 70.31%, which made the d-excess decrease for evaporation.When the water vapor pressure changed to 10-15 hPa, the increased extents of δ 18 O and δ 2 H continued to decrease and the d-excess decreased significantly with its increasing (r = −0.11).At this time, the average T was 15.9 • C, and the RH was 70.29%, which caused the evaporation to enhance and the d-excess further to decrease.When the water vapor pressure was greater than 15 hPa, δ 18 O and δ 2 H with it changed to negative correlations (r = −0.23),while the d-excess increased with its increasing (r = 0.05).At this moment, the average T and the RH were the highest, and the water vapor pressure tended to saturate and the below-cloud evaporation effect was weak, which made the d-excess increase.It has been shown that the increases of temperature, the water vapor content and the water vapor pressure led to isotope enrichment in arid climates [65,67].As can be seen in Table 7, the time stage with higher temperature also had higher relative humidity and water vapor pressure as well as more precipitation.Ren et al. [53] found that water vapor gradually saturated with precipitation proceeding, which resulting in the δ 18 O values decreasing and d-values increasing.

Conclusions
Using stable isotope and meteorological record data, the spatial and temporal variation characteristics of stable isotopes in precipitation and their relationships with meteorological factors in the Shiyang River Basin were studied and the main conclusions are as follows.
On the temporal scale, the stable isotopes of precipitation in the Shiyang River Basin showed obvious and periodic seasonal variations, with δ 2 H and δ 18 O higher values in summer and autumn and lower values in winter and spring, and d-excess higher values in spring and autumn and lower values in winter and summer.On the spatial scale, δ 2 H and δ 18 O values in precipitation decreased with the elevation increasing, and the opposite was true for d-excess values.The slope and intercept of local meteoric water lines gradually increased with the elevation increasing.
δ 18 O in precipitation shows a significant positive correlation with temperature.From 0 °C to 10 °C, the positive correlation was the largest, and the temperature effect of δ 18 O

Conclusions
Using stable isotope and meteorological record data, the spatial and temporal variation characteristics of stable isotopes in precipitation and their relationships with meteorological factors in the Shiyang River Basin were studied and the main conclusions are as follows.
On the temporal scale, the stable isotopes of precipitation in the Shiyang River Basin showed obvious and periodic seasonal variations, with δ 2 H and δ 18 O higher values in summer and autumn and lower values in winter and spring, and d-excess higher values in spring and autumn and lower values in winter and summer.On the spatial scale, δ 2 H and δ 18 O values in precipitation decreased with the elevation increasing, and the opposite was true for d-excess values.The slope and intercept of local meteoric water lines gradually increased with the elevation increasing.
δ 18 O in precipitation shows a significant positive correlation with temperature.From 0 • C to 10 • C, the positive correlation was the largest, and the temperature effect of δ 18 O was 0.77‰/ • C. When the temperature was greater than 0 • C, the temperature effect was greater in the mid-downstream area than in the upstream area.δ 18 O in precipitation showed a weak negative correlation with precipitation when the precipitation was greater than 10 mm.On the seasonal scale, the rainfall amount effect to δ 18 O exhibited in summer is influenced by monsoonal moisture.
There are significant differences in stable isotopes of precipitation in the entire watershed of arid inland river basins, the upstream area was more influenced by local moisture recycling but the mid-downstream area was more influenced by sub-cloud secondary evaporation.Our research reveals the process and influencing factors of stable isotope changes in precipitation in the Shiyang River Basin, providing a basis for further research on the evolution mechanism of stable isotope in precipitation in arid inland river areas.At the same time, it helps to understand the hydrological and meteorological processes and provides a certain theoretical basis for the rational development and utilization of water resources in arid areas.

Figure 1 .
Figure 1.Location of Shiyang River Basin and distribution of sampling points.The background (left) is based on Natural Earth (https://www.naturalearthdata.com/accessed on 8 June 2023).

Figure 1 .
Figure 1.Location of Shiyang River Basin and distribution of sampling points.The background (left) is based on Natural Earth (https://www.naturalearthdata.com/accessed on 8 June 2023).
higher d-values in autumn.In winter, the raindrops that fall in the dry atmosphere are influenced by sub-cloud secondary evaporation, which reduces lower d-values.Water 2023, 15, x FOR PEER REVIEW 8 of 22 d-values are lower in summer, however, the higher contribution of local recycling moisture causes higher d-values in autumn.In winter, the raindrops that fall in the dry atmosphere are influenced by sub-cloud secondary evaporation, which reduces lower d-values.

Figure 4 .
Figure 4. Monthly variation of δ 18 O, δ 2 H and d-excess values in Shiyang River Basin (the average values of the box and whisker plots were based on values from June 2018 to May 2020).
O values during the summer in 2018.Meanwhile, the inter-annual variation of d-values was also different in that its value in June was lower in the hydrological year of June 2018/May 2019 than in the hydrological year of June 2019/May 2020 significantly.A total of two persistent precipitation events of more than 3 days occurred in June 2018, which all occurred at Lenglongling (16 June 2018 to 18 June 2018 and 23 June 2018 to 27 June 2018).A total of six persistent precipitation events of more than 3 days occurred in June 2019, of which four occurred at Lenglongling (3 June 2019 to 6 June 2019, 9 June 2019 to 11 June 2019, 13 June 2019 to 17 June 2019 and 20 June 2019 to 27 June 2019), once occurred at Xiyingwugou

Figure 4 .
Figure 4. Monthly variation of δ 18 O, δ 2 H and d-excess values in Shiyang River Basin (the average values of the box and whisker plots were based on values from June 2018 to May 2020).
O values during the summer in 2018.Meanwhile, the inter-annual variation of d-values was also different in that its value in June was lower in the hydrological year of June 2018/May 2019 than in the hydrological year of June 2019/May 2020 significantly.A total of two persistent precipitation events of more than 3 days occurred in June 2018, which all occurred at Lenglongling (16 June 2018 to 18 June 2018 and 23 June 2018 to 27 June 2018).A total of six persistent precipitation events of more than 3 days occurred in June 2019, of which four occurred at Lenglongling (3 June 2019 to 6 June 2019, 9 June 2019 to 11 June 2019, 13 June 2019 to 17 June 2019 and 20 June 2019 to 27 June 2019), once occurred at Xiyingwugou (21 June 2019 to 23 June 2019) and once occurred at Jiuduntan (25 June 2019 to 27 June 2019).
precipitation, the relative humidity increased and the local moisture recirculation was enhanced, which made the d-values higher.respectively,and the weighted-average values were −60.3‰, −9.41‰ and 15.03‰ respectively.For the hydrological year of June 2019/May 2020, the precipitation samples were 221, the annual precipitation was 362.5 mm, mean temperature was 1.5 °C, and the variations of δ 2 H, δ18 O and d-values ranged from −203.4 to 35.9‰, −25.7 to 3.6‰ and −13.2 to 30.6‰ respectively, and the weighted-average values were −55.4‰, −8.8‰ and 14.6‰ respectively.Compared to the hydrological year of June 2019/May 2020, though there were higher mean temperatures and fewer precipitation events, the more precipitation in the hydrological year of June 2018/May 2019 caused the isotope values to have smaller variation ranges and lower mean values, which indicated the precipitation had a significant effect on the variation of stable isotopes.

Figure 5 .
Figure 5. Monthly variations of δ 2 H values, δ 18 O values, d-excess values, temperature and total precipitation in the Shiyang River Basin at Lenglongling, Hulinzhan, Xiyingwugou and Jiuduntan from June 2018 to May 2020.

Figure 5 .
Figure 5. Monthly variations of δ 2 H values, δ 18 O values, d-excess values, temperature and total precipitation in the Shiyang River Basin at Lenglongling, Hulinzhan, Xiyingwugou and Jiuduntan from June 2018 to May 2020.For the hydrological year of June 2018/May 2019, the precipitation samples were 209, the annual precipitation was 386.4 mm, mean temperature was 1.6 • C, and the variations of δ 2 H, δ 18 O and d-values ranged from −238.6 to 22.3‰, −31.5 to 3.2‰ and −22.1 to 31.0‰ respectively, and the weighted-average values were −60.3‰, −9.41‰ and 15.03‰ respectively.For the hydrological year of June 2019/May 2020, the precipitation samples were 221, the annual precipitation was 362.5 mm, mean temperature was 1.5 • C, and the variations of δ 2 H, δ 18 O and d-values ranged from −203.4 to 35.9‰, −25.7 to 3.6‰ and −13.2 to 30.6‰ respectively, and the weighted-average values were −55.4‰, −8.8‰ and 14.6‰ respectively.Compared to the hydrological year of June 2019/May 2020, though there were higher mean temperatures and fewer precipitation events, the more precipitation in the hydrological year of June 2018/May 2019 caused the isotope values to have smaller variation ranges and lower mean values, which indicated the precipitation had a significant effect on the variation of stable isotopes.

Figure 7 .
Figure 7. Relationship between δ 18 O in precipitation and temperature in the Shiyang River Basin (a).monthly scale, (b).daily scale.

Figure 7 .
Figure 7. Relationship between δ 18 O in precipitation and temperature in the Shiyang River Basin (a).monthly scale, (b).daily scale.

Figure 8 .
Figure 8. Variation of δ 18 O and d-excess of precipitation with temperature in the Shiyang River Basin.

Figure 10 .
Figure 10.Variation of δ 18 O and d-excess in precipitation with the relative humidity in the Shiyang River Basin.

Figure 10 .
Figure 10.Variation of δ 18 O and d-excess in precipitation with the relative humidity in the Shiyang River Basin.

Figure 11 .
Figure 11.Variation of precipitation δ 18 O and d-excess with water vapor pressure in Shiyang River Basin.

Figure 11 .
Figure 11.Variation of precipitation δ 18 O and d-excess with water vapor pressure in Shiyang River Basin.

Table 1 .
Information of sampling points in the Shiyang River Basin.

Table 2 .
Isotope values and local meteoric water line in each sampling point of the Shiyang River Basin.

Table 3 .
Relationship between δ 18 O and d-excess in precipitation with temperature in different ranges in the Shiyang River Basin.

Table 4 .
Relationship between stable isotopes of precipitation and precipitation amount in the Shiyang River Basin during different seasons.

Table 5 .
Relationship between stable isotopes of precipitation and precipitation amount under different rainfall conditions in the Shiyang River basin.

Table 6 .
Relationship between stable isotopes of precipitation with the relative humidity and mean values of meteorological elements in the Shiyang River Basin.

Table 7 .
Relationship between stable isotopes of precipitation with water vapor pressure and mean values of meteorological elements in the Shiyang River Basin.

Table 7 .
Relationship between stable isotopes of precipitation with water vapor pressure and mean values of meteorological elements in the Shiyang River Basin.