Quantifying Recycled Moisture in Precipitation in Qilian Mountains

: Studies about the hydrological cycle based on basin or regional scales often ignore the uniqueness of recycling moisture in mountain areas, and little effort has been made to understand the impact of the local recycled moisture on precipitation in mountain areas. We collected and analyzed a series of samples (stable isotope of precipitation, soil water, plant water, runoff, and groundwater) in the Qilian Mountains, northwest of China. Based on the isotopic mixing model, the characteristics of recycled moisture in the Qilian Mountains were evaluated. The results showed that lateral advection moisture is the primary source of precipitation (83.5~98.38%). The contribution rate of recycled moisture to precipitation was higher in the spring, summer, and autumn (2.05~16.5%), and lower in the winter (1.62~3.32%). The contribution of recycled moisture to precipitation in the high-elevation areas (>2400 m) was higher than that in the foothills area (2100~2400 m). The contribution of vegetation transpiration ( f Tr ) to precipitation in the east of Qilian Mountain was higher than that of the land surface evaporation ( f E v ). These proved that in the eastern part of Qilian Mountain, the arge-scale water cycle has a greater impact on precipitation in the area. The inﬂuence of local circulating water on precipitation dominated in the summer half of the year. Understanding the contribution of local circulating water to precipitation in the eastern part of Qilian Mountain will help us to understand the local hydrothermal conditions better and provide a basis for rationally arranging local agricultural production activities.


Introduction
Precipitation comes from the antecedent atmospheric moisture, lateral advection, and local evapotranspiration [1,2].Research has shown that lateral advection is the major moisture source of precipitation globally, with f Tr (proportion of plant transpiration water vapor in precipitation) and f Ev (proportion of surface evaporation water vapor in precipitation) being the second and third.However, there is considerable spatiotemporal variation among the three types of moisture in the world [3,4], which profoundly affects the global and local water cycle [5].Stable hydrogen and oxygen isotopes can play an important role in the quantitative research of the hydrologic cycle [6,7].The contribution of different sources of moisture to precipitation has been a hot topic in isotope hydrology.In general, it is feasible to observe and calculate the antecedent atmospheric water vapor and lateral advection directly.The difficulty is to determine the contribution of recycled moisture.Recycled moisture mainly comes from f Tr (proportion of plant transpiration water vapor in precipitation) and f Ev (proportion of surface evaporation water vapor in precipitation) [8][9][10][11].The linear mixing model for isotopes is an effective method to study the contribution of recycling moisture in different regions.Linear mixing model has been applied around the world [12], such as the Great Lakes region in North America [13], the Slave River Delta of Canada [14], the Amazon Basin of South America [15], the Nam Co Basin in the Qinghai-Tibetan Plateau [16], the Lake Shorty in Madagascar [17] and Lake Kasumiguara [18].
Most research on the moisture cycle focuses on the source of advected water vapor and its transport in different regions [8,19].Recently, more studies have considered the transport and conversion of basin recycled moisture [4,6].Research showed that the contribution of recycled moisture to precipitation varies greatly in time and space.In arid inland river basins, the contribution of recycled moisture is less than 30% [6,20], In some small oases, the recycled moisture contribution is less than 5% [21,22].However, the contribution of recycled moisture can reach 62% in the Tibetan Plateau, but is much lower in winter than that in summer [6,23].Some studies have focused on the causes of the spatial and temporal differences in recycled moisture, such as relative air humidity [19], soil water content [24], land use types and land cover changes [25].
Many studies have estimated the contribution rate of recycled moisture and its influencing factors in the watershed or regional scale [4,6,21,22,26].Due to the obvious difference between climate and environment in the mountainous-oasis-desert regions in arid areas, studies based on basins or regional scales tend to overlook the uniqueness of recycling moisture in mountain areas.The upper Shiyang River in the northeast of Qilian Mountain is the transition zone between the Qinghai-Tibet Plateau and the arid zone.Precipitation affects the development of the oasis and desert in the middle and lower reaches.Clarifying the characteristics and influencing factors of upstream water resources changes will contribute to a reasonable solution to the demands of the middle and lower reaches of the water resources.Therefore, this study used the stable isotope data of precipitation, soil water, plant water, surface water, and groundwater from 2016 to 2018 in the upper Shiyang River in the eastern part of the Qilian Mountains to calculate the proportion of plant transpiration water vapor (f Tr ), surface evaporation water vapor (f Ev ) and advection water vapor (f Adv ) in precipitation.We try to explore the source of moisture in precipitation in the mountain areas of the arid inland river basin, and to reveal the characteristics of the mountain water cycle, reasonable assessment of regional water resources.

Study Area
The study area is located in the Xiying River (XYR) basin in the Qilian Mountains (Figure 1).XYR is the main tributary of the Shiyang River.It originates below the glaciers in the northern slopes of the Qilian Mountains and eventually disappears into the desert.The elevation of the XYR basin ranges from 1510 to 4874 m above sea level.The average annual temperature is 6.3 • C and the average annual precipitation ranges from 200 to 700 mm [8,27].The upper reaches of the Shiyang River basin are located in the East Asian monsoon crisscross zone, controlled by the East Asian monsoon and plateau monsoon.Cold high pressure appears on the plateau in winter, and the air flows from the plateau to the surroundings.In summer, hot low air pressure appears on the plateau, and the air flows from all directions to the plateau [6], which is a typical continental alpine climate.Vegetation is mainly distributed in areas between 2000 m and 3600 m above sea level; the basin is affected by multiple sources of moisture [6,27].

Observation Network, Sampling, and Analysis
From April 2016 to October 2018, four sampling sites were established in the XY Basin (Figure 1).A total of 867 samples were collected from the upstream mountain are (Table 1).Two-hundred forty-five precipitation samples at sampling sites were collecte immediately after the end of each precipitation process using a rain gauge.The precipit tion samples were put into a 50 mL polyethylene sampling bottle.The bottle cap w tightened, and the bottle mouth was sealed with a sealing film and stored in cold storag until analysis.Surface water samples were sealed and stored in cold storage after eac collection.Meanwhile, automatic meteorological observation instruments recorded met orological elements such as temperature, precipitation, relative humidity and atmospher pressure [6].All samples were analyzed for δ 2 H and δ 18 O using liquid water.For runo the bottle was placed under the water surface until the container was filled.A total of runoff samples were collected at different elevations (Table 1).For soil water, soil sampl were collected from 10 cm to 100 cm, using soil spirals every 10 cm.All samples we placed in 100 mL polyethylene vials.A total of 450 soil samples were collected (Table 1 For plant water, typical vegetation at each sampling site was selected, and scissors we used to cut the xylem or branches of the same vegetation.A total of 90 plant samples we collected (Table 1).All samples were stored in a mobile freezer (−5 °C) and transferred a freezer laboratory (−15 °C) within a week after collection.

Observation Network, Sampling, and Analysis
From April 2016 to October 2018, four sampling sites were established in the XYR Basin (Figure 1).A total of 867 samples were collected from the upstream mountain areas (Table 1).Two-hundred forty-five precipitation samples at sampling sites were collected immediately after the end of each precipitation process using a rain gauge.The precipitation samples were put into a 50 mL polyethylene sampling bottle.The bottle cap was tightened, and the bottle mouth was sealed with a sealing film and stored in cold storage until analysis.Surface water samples were sealed and stored in cold storage after each collection.Meanwhile, automatic meteorological observation instruments recorded meteorological elements such as temperature, precipitation, relative humidity and atmospheric pressure [6].All samples were analyzed for δ 2 H and δ 18 O using liquid water.For runoff, the bottle was placed under the water surface until the container was filled.A total of 82 runoff samples were collected at different elevations (Table 1).For soil water, soil samples were collected from 10 cm to 100 cm, using soil spirals every 10 cm.All samples were placed in 100 mL polyethylene vials.A total of 450 soil samples were collected (Table 1).For plant water, typical vegetation at each sampling site was selected, and scissors were used to cut the xylem or branches of the same vegetation.A total of 90 plant samples were collected (Table 1).All samples were stored in a mobile freezer (−5 • C) and transferred to a freezer laboratory (−15 • C) within a week after collection.
The samples were melted at room temperature (20~25 • C) before further analyses.The plant and soil water were extracted using a cryogenic vacuum distillation apparatus (LI-2100, LICA United Technology Limited, China).All samples were analyzed using a liquid water isotope analyzer DLT-100 (Los Gatos Research, Inc.) in the Stable Isotope Laboratory, College of Geography and Environmental Science, Northwest Normal University.The isotope ratios of samples are expressed as parts per mil (‰) relative to the Vienna Standard Mean Ocean Water (V-SMOW) using δ notation: δ(‰) = (R S /R V-SMOW −1) × 1000, where δs is the isotope ratio of the samples relative to V-SMOW and Rs is the ratio of D/H or 18 O/ 16 O in the samples.The precision of the measurements was ±0.6% for δ 2 H and ±0.2% for δ 18 O, respectively.Meteorological data were obtained from four weather stations of the XYR basin.These weather stations record meteorological data such as temperature, precipitation, wind speed, evaporation, relative humidity, and soil moisture every 15 min.

Three-Component Mixing Model
Moisture from precipitation is derived from local evapotranspiration and advection [28,29].The linear mixing model can be used to calculate the contribution rate of each water source.The three-component mixing model based on δ 18 O and δ 2 H values can be expressed as follows.
In formula ( 1)-(3), the contribution rate of each moisture source is expressed in terms of f while the subscript indicates the source of moisture.Subscript P v indicates the local atmospheric moisture.Subscripts Tr, Ev, Adv, indicate moisture from transpiration, evaporation, and advection, respectively.
The contributions of f Tr , f Ev , and f Adv can be calculated as follows.
The δ 18 O pv and δD pv value can also be derived from the value of the stable isotope from local precipitation and the equilibrium fractionation factor: ) The specific formula is as follows [30,31]: The δ 18 O Ev and δD ev can be expressed by the value of the stable isotope from advection (δ 18 O Adv and δD Adv ) and local the surface water (δ 18 O s or δD s ), mean relative humidity (h), the sum of equilibrium (ε eq ) and kinetic (∆ε): (13) Based on previous studies on δ 18 O Adv and δD Adv , we used the following formula to investigate the characteristics in the XYR Basin.
In the above formulas, the F indicates the ratio between the initial and the final vapor, which is estimated by the precipitable water amount in the two sites.Previous research in the Tianshan Mountains concluded that precipitable water correlates with moisture pressure (c = 1.657e,where c indicates atmospheric moisture content in mm and e indicates surface moisture pressure in hPa, r 2 = 0.94) [32].Hence, the surface moisture pressure ratio between the two sampling sites is equal to the value of F. Since the isotope ratios in precipitating vapor at the sampling sites are much more depleted than the upwind station, the Rayleigh distillation equation was also applied.In this research, formula 6 was used to calculate the stable isotope ratios in advection vapor.If there is no significant depletion of isotopes between the sampling site and the upwind station, then the stable isotope ratios in advection vapor at the sampling site are considered to be the same as that in the precipitating vapor at the upwind station.
Because there is no fractionation during the transport of moisture from f Tr to the atmosphere [33], the value of δ 18 O Tr and δD Tr is the same as that of local water used by plants (δ 18 O w and δD w ).This research calculated the average value of δ 18 O (δ 2 H) from the soil surface to a depth of 40cm below and xylem, to obtain δ 18 O Tr and δD Tr .The results for each sampling site in this study area are shown in Table 2.

Hysplit Model and the Upper Wind Direction
We applied the HYSPLIT model to simulate the moisture sources in the Qilan Mountains [34][35][36][37].We found that westerly winds, southeast monsoons, and plateau monsoons all affect the Qilian Mountains in summer.In winter, westerly winds mainly affect the Qilian Mountains.
According to the clustering of air mass in different seasons, the air mass gathered at the northern foot of Qilian Mountain and then moved from a low elevation to a high elevation along the valley.Therefore, sampling site A was used as an upwind station for spring, summer, and autumn.In winter, the study area was dominated by a westerly wind, and Urumqi and Hotan (GNIP) in Central Asia were regarded as the upwind direction stations.

Method of Reducing Experimental Error
The isotope mixture model is an effective method for specific regional moisture recirculation research, but the inherent uncertainty of the model still needs to be considered.In this study, we focused on the correction of experimental errors in plant xylem moisture.
The water obtained in plant Xylem contained organic pollutants such as methanol and ethanol by the low-temperature vacuum distillation extraction technology, which caused deviations in the measurement of the laser isotope analyzer.This error has led to significant differences in the estimation of the amount of vegetation evapotranspiration.In different studies, some unreasonable calculation results of negative f Tr value will appear.
In this study, deionized water was mixed with methanol and ethanol (pure chromatographical) of different concentrations.The spectral software (LWIA-Spectral Contamination Identifier v1.0, Los Gatos company) was used to determine the spectral measurement of the pollution degree of methanol (NB), and ethanol (BB), the δD and δ 18 O spectral pollution correction methods were established [38][39][40].The configuration of the concentration of methanol and ethanol solutions during the calibration process is the same as the related experiments by Meng et al. (2012).For the calibration result of methanol, the logarithm of the broadband metric NB metric and ∆δ 2 H and ∆δ 18 O have significant quadratic curve relationships: ∆δ 2 H = 0.018 (ln.NB) 3 + 0.092 (ln.NB) 2 + 0.388ln.NB + 0.785 (R 2 = 0.991, p < 0.0001) (21) ∆δ 18 O = 0.017 (ln.NB) 3 + 0.017 (ln.NB) 2 + 0.545ln.NB + 1.356 (R 2 = 0.998, p < 0.0001) (22) For the calibration results, the broadband metric BB metric has a quadratic curve and linear relationship with ∆δ 2 H and ∆δ 18 O, respectively: After correction, the calculation result of negative f Tr value was eliminated.Of course, there are also other uncertainties in model research, such as driving data sam-pling and experimental errors, structural errors in the physical mechanism of the model, and parameter errors.

Isotopic Composition of Different Water Bodies
The differences in the Local Meteoric Water Line (LMWL) are explained by the differences in topography, underlying surface, climate, and other natural environmental factors.From sampling points A to D, both the temperature and the evaporation decrease because of the increasing elevation.The slope of the LMWL gradually increases, and the extent of soil water and plant water deviating from the LMWL also increases.These factors indicate that as the elevation increases, the evapotranspiration of soil and vegetation gradually decreases.
In sampling point D (River source), the δ 2 H and δ 18 O in different water bodies are similar.The δ 2 H and δ 18 O of soil moisture, vegetation moisture, and surface runoff of each sampling site fall to the lower right of the LMWL, indicating that different water bodies are all recharged by precipitation in the source areas, which is significantly different from sampling point A, B, and C.
The characteristics of δ 2 H and δ 18 O are similar in the different water bodies in the mountain area in sampling points A, B, and C. The δ 2 H and δ 18 O of the surface water samples plot near or above the LMWL (Figure 2), and the soil and vegetation water plot to the lower right of the LMWL (Figure 2).The distribution of samples indicates that the land surface water in the mountain area is mainly recharged by precipitation, while soil water and plant water experience different degrees of evaporation during water body migration and transformation.

Isotopic Composition of Different Water Bodies
The differences in the Local Meteoric Water Line (LMWL) are explained by the differences in topography, underlying surface, climate, and other natural environmental factors.From sampling points A to D, both the temperature and the evaporation decrease because of the increasing elevation.The slope of the LMWL gradually increases, and the extent of soil water and plant water deviating from the LMWL also increases.These factors indicate that as the elevation increases, the evapotranspiration of soil and vegetation gradually decreases.
In sampling point D (River source), the δ 2 H and δ 18 O in different water bodies are similar.The δ 2 H and δ 18 O of soil moisture, vegetation moisture, and surface runoff of each sampling site fall to the lower right of the LMWL, indicating that different water bodies are all recharged by precipitation in the source areas, which is significantly different from sampling point A, B, and C.
The characteristics of δ 2 H and δ 18 O are similar in the different water bodies in the mountain area in sampling points A, B, and C. The δ 2 H and δ 18 O of the surface water samples plot near or above the LMWL (Figure 2), and the soil and vegetation water plot to the lower right of the LMWL (Figure 2).The distribution of samples indicates that the land surface water in the mountain area is mainly recharged by precipitation, while soil water and plant water experience different degrees of evaporation during water body migration and transformation.

The Path of Moisture Transport
In spring, there are two main air mass movement paths moving from the western section of the study area: (1) The air mass originating from Central Asia moves along the edge of the Qinghai-Tibetan Plateau after entering the Tarim Basin; (2) The air mass originating from the West Siberian Plain, which travels the western arid region of China and then arrives at the XYR basin.In summer, the air mass mainly comes from the east of the study area, and the air mass moves from the east of the study area, similar to the pathway

The Path of Moisture Transport
In spring, there are two main air mass movement paths moving from the western section of the study area: (1) The air mass originating from Central Asia moves along the edge of the Qinghai-Tibetan Plateau after entering the Tarim Basin; (2) The air mass originating from the West Siberian Plain, which travels the western arid region of China and then arrives at the XYR basin.In summer, the air mass mainly comes from the east of the study area, and the air mass moves from the east of the study area, similar to the pathway in spring.However, as air mass from the west of the study area decreased in occurrence, the air mass from the southeast increased and the air mass from Central Asia affects the study area along the northwest edge of the Qinghai-Tibetan plateau.In autumn, there are two main air mass movement paths in the west of the study area: (1) Air mass originating in Central Asia and the air mass originating in the Xinjiang move along the edge of the Qinghai-Tibetan Plateau; (2) The other is the air mass arises in Central Asia, crosses the Kunlun Mountains and the Qinghai-Tibetan Plateau, moves through the Qaidam basin and then reaches the study area.In winter, the air mass in the study area is controlled by the westerly wind, while the east wind of the study area has little effect on the region (Figure 3). in spring.However, as air mass from the west of the study area decreased in occurren the air mass from the southeast increased and the air mass from Central Asia affects study area along the northwest edge of the Qinghai-Tibetan plateau.In autumn, there two main air mass movement paths in the west of the study area: (1) Air mass origina in Central Asia and the air mass originating in the Xinjiang move along the edge of Qinghai-Tibetan Plateau; (2) The other is the air mass arises in Central Asia, crosses Kunlun Mountains and the Qinghai-Tibetan Plateau, moves through the Qaidam ba and then reaches the study area.In winter, the air mass in the study area is controlled the westerly wind, while the east wind of the study area has little effect on the reg (Figure 3).

Spatial and Temporal Differences of Recycled Moisture
In spring (Figure 4a), the calculated contribution of recycled moisture in differ sampling points was 9.9% and 16.5%, higher than that in other seasons.The contribu was largest in sampling point B (Arbor belt), and lowest in sampling point D (river sou area).The melting of snow leads to higher soil moisture content in spring, which lead higher evaporation of the soil.Plants at low elevations region begin to grow first in spr which increases the vegetation evapotranspiration in the low elevation region.The p tion of recycled moisture decreased with elevation increased (Table 2).

Spatial and Temporal Differences of Recycled Moisture
In spring (Figure 4a), the calculated contribution of recycled moisture in different sampling points was 9.9% and 16.5%, higher than that in other seasons.The contribution was largest in sampling point B (Arbor belt), and lowest in sampling point D (river source area).The melting of snow leads to higher soil moisture content in spring, which leads to higher evaporation of the soil.Plants at low elevations region begin to grow first in spring, which increases the vegetation evapotranspiration in the low elevation region.The portion of recycled moisture decreased with elevation increased (Table 2).
In summer (Figure 4b), the contribution of recycled moisture for different sampling points varied between 3.52% and 12.87%.The contribution was highest for sampling point D (River source) and the lowest for sampling point B (Arbor belt).At the sampling point D (River source), where the plant begins to grow first in July, the contribution of local transpiration f Tr was 4.86%, and the contribution of evaporation f Ev was 7.99%.The soil moisture content is high in the river source area, and the frozen soil has thawed in summer, so that the evapotranspiration increases rapidly, which leads to a larger fraction of recycled moisture in precipitation in the river source area.Since there is less precipitation and soil water in the low-elevation regions, the contribution of recycled moisture is higher in high-elevation areas than in low-elevation areas in summer.
In autumn (Figure 4c), the contribution of recycled moisture in the study area varied between 2.05% and 9.63%, which is lower than in spring and summer.The contribution was highest at sampling point C (Shrub belt) and lowest at sampling point B (Arbor belt) the vegetation growth in the river source stagnated, and the soil began to freeze.The sampling point C (Shrub belt) has high evapotranspiration, and the evapotranspiration was reduced due to the lack of soil moisture in sampling points A and B. In summer (Figure 4b), the contribution of recycled moisture for different sampling points varied between 3.52% and 12.87%.The contribution was highest for sampling point D (River source) and the lowest for sampling point B (Arbor belt).At the sampling point D (River source), where the plant begins to grow first in July, the contribution of local transpiration fTr was 4.86%, and the contribution of evaporation fEv was 7.99%.The soil moisture content is high in the river source area, and the frozen soil has thawed in summer, so that the evapotranspiration increases rapidly, which leads to a larger fraction of recycled moisture in precipitation in the river source area.Since there is less precipitation and soil water in the low-elevation regions, the contribution of recycled moisture is higher in high-elevation areas than in low-elevation areas in summer.
In autumn (Figure 4c), the contribution of recycled moisture in the study area varied between 2.05% and 9.63%, which is lower than in spring and summer.The contribution was highest at sampling point C (Shrub belt) and lowest at sampling point B (Arbor belt) the vegetation growth in the river source stagnated, and the soil began to freeze.The sampling point C (Shrub belt) has high evapotranspiration, and the evapotranspiration was reduced due to the lack of soil moisture in sampling points A and B.
In winter, due to the influence of the cold westerly air mass, the rate of evapotranspiration weakened.The external air mass exerts a dominant influence on the moisture in winter.
Overall, the contribution rate of recycled moisture is subject to local vegetation cover, soil moisture content, and other climatic and hydrological conditions.At all sampling points, the contribution fTr moisture was higher than that of fEv moisture (Figure 4).In winter, due to the influence of the cold westerly air mass, the rate of evapotranspiration weakened.The external air mass exerts a dominant influence on the moisture in winter.
Overall, the contribution rate of recycled moisture is subject to local vegetation cover, soil moisture content, and other climatic and hydrological conditions.At all sampling points, the contribution f Tr moisture was higher than that of f Ev moisture (Figure 4).

Differences between Mountainous Areas in Inland River Basins
The XYR Basin is located in the east of the Qilian Mountains and is a marginal area of the summer monsoon.Its moisture source is significantly different from that in most arid regions.The portion of moisture from the southeast and southwest in the monsoon period (June to September) is higher than in most inland river basins [6].The sum of f Tr and f Ev recirculated moisture in the XYR Basin is about 10.34%, which is lower than that in Yeniugou, Hulugou, and Pailigou, which are located in the central Qilian Mountains [4], and higher than the Tianshan Mountains [21].The main reason for the difference in the contribution of recirculated moisture in different regions is f Tr .The f Tr values calculated in various studies vary greatly, while the proportion of f Ev is close to other stations.For XYR basin f Ev is 4.81%, for bison ditch 3.6%, for Hulu ditch 5.9%, and for the dew ditch 0.9%.Controlled by meteorological factors, there are also obvious differences in recycled moisture at different elevations in the same watershed.In high elevation areas, f Tr and f Ev will be higher due to summer plant growth and frozen soil melting.At sampling point A (Foothill), f Tr is relatively higher in spring and autumn, but soil moisture is lower due to less precipitation, resulting in a lower f Ev .
In terms of spatial distribution, the vegetation cover is higher in site B and C than D. However, the recycling ratio in D is the highest in all four seasons.A possible reason accounting for this spatial pattern is that recycling ratio is scale-dependence [41].We used site A as the unique upwind site for all sites B-D.As a result, the recycling ratio in D is the accumulated recycled moisture from site A to D and accumulated from site A to B for recycling ratio in site B. That is why recycling ratio increases with elevation rise, while vegetation and evapotranspiration are actually low in high-elevation regions due to the low air temperature.Scale-dependence is an important issue in recycling research but often is overlooked in isotope-based studies.This is also the reason why the isotope-based result is much smaller than model-based results.

Recycled Moisture Contribution in Mountainous, Oasis, and Desert Areas
The average contribution of f Tr , f Ev , and f Adv to precipitation was 5.53%, 4.8%, 90.89% in the mountains area, 21.9%, 7%, and 72% in the oasis region and 10%, 5% and 85% in the desert and gobi region [6].In Heihe River Basin, the average contributions of f Tr , f Ev , and f Adv vapor to precipitation were 24.15%, 26.9%, and 51.05% in the oasis region, 15.1%, 6.3%, and 21.4% in the desert region [4].In the Urumqi River Basin, the average contribution of oasis f Tr and f Ev to precipitation was 15.09% [21].Our study mainly studies the Xiying River Basin in the eastern part of the Qilian Mountains, this area is mountainous, the f Tr (proportion of plant transpiration water vapor in precipitation) and f Ev (proportion of surface evaporation water vapor in precipitation) is less than 10% throughout the year.It is in line with the local water cycle characteristics of the mountain areas in arid inland river basins.In general, the proportion of recycled water in the mountain areas of arid inland river basins are lower than those in the oasis and desert areas.
With the increase in population and utilization of the oasis area, the land use has changed.The cultivated land and ecological land must be maintained by artificial irrigation in the oasis area.Since canals and flooding can affect local evaporation, the growth of crops and forests can affect f Tr , which will change the contribution of recycled moisture rates in oasis areas, and this is different from the characteristics of moisture recirculation in the mountainous area.

Recycled Moisture in Precipitation in Different Regions
The contribution of recycled moisture is lowest in the North and South poles, while the highest values are found in the Tibetan Plateau, the Patagonia Plateau, and the Andes Mountains.The vast ocean provides a large amount of moisture for the moisture cycle in various regions in the world.At the same time, the latent heat from condensation absorbed and released by the vapor phase transition also promotes the flow of global energy.There is a continuous high-value area in the mid-latitudes regions [42].In tropical coastal areas, the contribution of recirculating moisture in different landscapes is only marginally different, with a contribution around 31~37% [19].In the marginal zone of the temperate monsoon, the source of water is complex, and the climate is changeable.The contribution of circulating water in mountainous areas, oasis and deserts, and its temporal and spatial changes are quite large [6].In addition, lakes have a significant impact on recycled moisture, contributing 5-16% in temperate continental monsoon climate zones [13], and 10-20% in temperate marine climate zones [43], 16-50% in tropical islands [17], and 3~37.9% are in the Qinghai-Tibet Plateau [16].Our study area is the eastern part of the Qilian Mountains.This area is on the edge of the East Asian monsoon.The temporal and spatial changes of the local water cycle are complex.The proportion of local circulating water in precipitation is smaller than that in mid-latitude regions and tropical coastal regions.However, for small areas, the circulating water in mountain areas is greater than in oasis and desert areas.
The δ 2 H was used in the ice core to estimate that the contribution of recycled moisture in the Qinghai-Tibet Plateau has increased in the past few decades [44,45].It is believed that the increase in the global temperature leads to strong local surface moisture evaporation and local moisture recycling [46].Secondly, any increase in the vegetation coverage on the land surface and the evapotranspiration associated with this increase also has a particularly strengthening effect on local moisture recycling [6].When we study the local water cycle, we should also consider the effects of long-term climate change and local human activities.

Conclusions
The contribution of recycled moisture to precipitation was estimated using the linear mixing model in different seasons for different elevations.The results showed that the contribution rate of recycled moisture to precipitation was higher in the spring, summer, and autumn (2.05~16.5%)and low in winter (1.62~3.32%).In the high-elevation areas (>2400 m), the contribution rate of recycled moisture to precipitation was higher than the foothills area (2100~2400 m).The contribution of the f Tr moisture was higher than that of the f Ev to precipitation in the mountain areas.The contribution of recycled moisture increased with the elevation in the Qilian mountains.The water obtained from the xylem of plants contains organic pollutants (methanol and ethanol) that may cause of greater uncertainty in f Tr in different studies.The quantification of the contribution rate of recycled moisture requires the cooperation of multi-source data, but the high-elevation regions often lack systematic mufti-element observation data.To better understand the contribution mechanism and influencing factors of recycled moisture, a long-term field monitoring system should be established to obtain comprehensive first-hand data.This study can provide a reference for the study of recycled moisture in other mountain areas.

Figure 2 .
Figure 2. The relationship between δ 2 H and δ 18 O in different water bodies in the Xiying River.Basin.

Figure 2 .
Figure 2. The relationship between δ 2 H and δ 18 O in different water bodies in the Xiying River.Basin.

Figure 3 .
Figure 3. Movement path and cluster of the Xiying River basin sampling station (A, B, C, D from April 2016 to October 2018).

Figure 3 .
Figure 3. Movement path and cluster of the Xiying River basin sampling station (A, B, C, D from April 2016 to October 2018).

Sustainability 2021 , 13 Figure 4 .
Figure 4. Schematic diagram of recycled water vapor contribution rates for each season.

5 . Discussion 5 . 1 .
Recycled Moisture in Inland River Basins 5.1.1.Differences between Mountainous Areas in Inland River BasinsThe XYR Basin is located in the east of the Qilian Mountains and is a marginal area of the summer monsoon.Its moisture source is significantly different from that in most

Figure 4 .
Figure 4. Schematic diagram of recycled water vapor contribution rates for each season.

Table 1 .
Basic information for each sampling site.

Table 1 .
Basic information for each sampling site.

Table 2 .
The data needed to calculate the recycled water vapor contribution rate and the calculation results.