Soil Water Movement Changes Associated with Revegetation on the Loess Plateau of China

Soil water is the limitation factors in the semiarid region for vegetation growth. With the large scale “Grain for Green” implementation on the Loess Plateau of China, an amount of sloping cropland was converted to forestland, shrubland, and grassland. The spatial and temporal distribution of soil water was changed. However, the effect of revegetation on soil water movement is still unclear. In this study, we analyze the stable isotopes changes in precipitation and soil water in sloping cropland, forestland, shrubland, and grassland to trace the movement of moisture in soil. The results showed that δ18O in shallow layers (<20 cm depth) of sloping cropland, forestland, shrubland, and grassland were −3.54‰, −2.68‰, −4.00‰, and −3.16‰, respectively. The δ18O in these layers were higher than that in the lower layers, indicating that evaporation was mainly from the shallow layers. The δ18O for the soil water in the unsaturated zone in the grassland, shrubland, and forestland of the temporal variability decreases with depth and approaches a minimum value at 160 cm, 180 cm, and 200 cm, respectively, suggesting that the soil water is relatively stable many months or even longer. Precipitation was infiltrated with piston and preferential modes, and infiltration demonstrated obvious mixing. Present study demonstrated the δ18O was more sensitive than the soil water content for tracing the maximum infiltration depth of event water and recharge mechanisms. Consequently, we suggested that the land user management such as type, plant density should be considered in the revegetation.


Introduction
Acting as a limitation factors for the vegetation growth, soil water is crucial to the hydrological cycles in semiarid environments on the Loess Plateau of China (LPC).The Loess Plateau has the greatest erosion rate in the world.One of the causes of soil erosion is excessive use of land [1,2].With the large-scale "Grain for Green" implementation on the Loess Plateau, many of the sloping croplands were converted to forestland, shrubland, and grassland.The world is trying to achieve the sustainable goals for development by the United Nations and achieve the land degradation neutrality via the use of restoration [3,4].Land use type is one of the main driving forces affecting the evolution of soil and water resources in unsaturated zones [5][6][7][8].Understanding the movement of soil water in unsaturated zones in different land use types is important for understanding the distribution and cycling of water, salt, and other nutrients, and for improving land management for controlling erosion and runoff in revegetation and ecological construction [9,10].However, the influence of land use for the soil water movement through unsaturated zones on the Loess Plateau has been reported in a few papers.Accordingly, research on Water 2019, 11, 731 3 of 14 typical loessial hilly and gully landform.The difference in altitude from top to bottom is typically from 100 to 200 m.The soil of the study area has developed on wind-deposited loessial parental material [38].The most common soil is loess with silt contents range from 62% to 73% and clay contents range from 17% to 20% (Table 1).The Wangmaogou watershed responds to the country policy of Grain to Green Project (GTGP) since 1953.With the developing of GTGP on the Loess Plateau, many of the sloping croplands were converted to forestland, shrubland, and grassland.
Water 2019, 11, x FOR PEER REVIEW 3 of 14 bottom is typically from 100 to 200 m.The soil of the study area has developed on wind-deposited loessial parental material [38].The most common soil is loess with silt contents range from 62% to 73% and clay contents range from 17% to 20% (Table 1).The Wangmaogou watershed responds to the country policy of Grain to Green Project (GTGP) since 1953.With the developing of GTGP on the Loess Plateau, many of the sloping croplands were converted to forestland, shrubland, and grassland.
The results shown in Figure 3 indicate that the SWC located in the shallow layers presents a big difference, but the SWC in the deeper layers was most stable for grassland, forestland, and

Sampling Methods
Under comprehensive consideration of land use types of revegetation, soil samples were collected from the grassland, shrubland, forestland, and sloping cropland in May, August, and November 2016 by drilling with a hand auger of 9 cm in diameter, at 10 cm intervals to a depth of 100 cm and at 20 cm intervals from 100 to 400 cm.Three replicates were collected for each soil sample.The parts of soil samples were packed into an aluminum box sealed with plastic wrap and measured for gravimetric soil water content.The parts of soil samples were then packed into glass bottles sealed with Parafilm to prevent evaporation and promptly brought to State Key Laboratory of Eco-hydraulics in Northwest Arid Region of China, Xi'an University of Technology, stored at −4 • C until laboratory analysis.A Mastersizer 2000 laser particle size analyzer (Malvern Instruments, Malvern, UK) was used to measure the percentages by volume of the particle sizes.Four undisturbed soil samples per replicate treatment were randomly selected from the soil layer (0-10 cm) for laboratory experiments using 50 mm × 54 mm cylindrical cores.Soil bulk density was determined from oven-dried undisturbed cores as mass per volume of oven-dried soil.Samples of precipitation were collected from the six sampling sites from January to December 2016.This work collected the precipitation using the ''Rain Gauge'' available in an all weather, into which a small amount (≈ 3 mL) of mineral oil had been added to prevent evaporation.During the winter months, snow was collected in the same collector and melted to determine the water equivalent.A total of 50 precipitation samples were collected from each sampling site.Three replicates were collected for each precipitation sample.Samples of precipitation were collected in glass bottles sealed with Parafilm to prevent evaporation and promptly brought to State Key Laboratory of Eco-hydraulics in Northwest Arid Region of China, Xi'an University of Technology, stored at −4 • C until laboratory analysis.

Gravimetric Soil Water Content
The content of soil water was measured by oven-drying.The gravimetric soil water content (SWC) of each sample was measured by weighing the samples in the field, heating the sample for 24 h at 100 • C, and weighing the dry soil.
where SWC is the gravimetric soil water content (%), M w is the weight of the samples in the field (g), and M D is the weight of the dry soil (g).

Measurements of Hydrogen and Oxygen Isotopes
A LI-2000 (LICA, China) liquid water vacuum extraction method was used to extract the soil water.The stable hydrogen and oxygen isotope composition (δ 18 O and δ 2 H) of the soil water was determined using an LGR liquid water isotope analyzer (Los Gatos Research Inc., USA) at the State Key Laboratory of Eco-hydraulics in Northwest Arid Region of China, Xi'an University of Technology.
Each sample was injected six times, but the first two injections were discarded to eliminate cross contamination.The accuracy of the measurements was ±0.3‰ for δ 2 H and ±0.1‰ for δ 18 O.The sample-injection volume (2.50 × 10 16 to 4.50 × 10 16 water molecules, ±3.00%), temperature (±1 • C), and the accuracy of the δ 2 H and δ 18 O measurements were checked.The test results that did not satisfy these conditions were eliminated from the analysis.The results are expressed as δ values relative to Vienna standard mean ocean water (VSMOW) in per mil (‰): Water 2019, 11, 731 6 of 14 where δ sample is the deviation of the isotope ratio of a sample relative from that of VSMOW; R sample is the ratio of 2 H to 1 H (or 18 O to 16 O) in the sample; and R VSMOW is the ratio of 2 H to 1 H (or 18 O to 16 O) in VSMOW.

Source Partitioning Using Stable Isotopes
The isotopic mass balance was calculated using Equations ( 3) and (4) [39]: where f new and f old are the proportions of new and old water (f new + f old = 1), respectively; δ new and δ old are the isotopic signatures of the new and old water, respectively, and δ mix is the isotopic signature of the mixed water.

Calculation of Biomass
Tree biomass was calculated by: where H is the tree height, D denotes the diameter at breast height, and W 1 and W 2 are the above-and belowground biomasses, respectively [40,41].Samples of rootlet biomass from herbaceous plants were collected by drilling (inner diameter of 100 mm) at 10-cm intervals to a depth of 100 cm, and every horizontal was mixed in a bag.The soil was then soaked and then washed under the tap water with a 0.25 mm sieve to isolate plant roots.The roots were air-dried, weighed, and stored at 65 • C to calculate the biomass.

Statistical Analysis
All statistical analyses were performed using SPSS 19.0 and Origin 8.5.One-way ANOVAs followed by Tukey's HSD tests (p < 0.05) were used to compare the effect of seasonal changes and land use types on SWC.

Isotopic Composition of the Precipitation
The statistical characteristics of the precipitation isotopes are shown in Figure 2. The δ 18 O of the precipitation varied greatly during the year.The precipitation isotopes from June to November were more depleted and were enriched in the study area from December to May (Figure 2).

Biomass of Four Land Use Types
Table 3 shows the vegetation cover's root distributions in the four land uses.Mean maize rootlet biomass in the upper 30 cm of the sloping cropland accounted for 84% of the total biomass.The root distribution of Bothriochloa ischaemum was well-proportioned, except for the lowest values in the 80-100 cm layer (Figure 3).The roots of Pinus tabuliformis were concentrated in the 40-60 cm layers, corresponding to the peaks of water loss in the profiles.The source of the soil water became depleted in 18-oxygen from May to November.The peak negative value was recorded at depths of 40 cm in the grassland and forestland, 50 cm in the shrubland and sloping cropland; rain fell on 8 and 12 November 2016 (30 mm).The δ 18 O of the soil The results shown in Figure 3 indicate that the SWC located in the shallow layers presents a big difference, but the SWC in the deeper layers was most stable for grassland, forestland, and shrubland.Sloping cropland retained water the best, with the highest SWCs throughout the study period.SWC in the 40-100 cm layers show a different value with the change of growing season.

Biomass of Four Land Use Types
Table 3 shows the vegetation cover's root distributions in the four land uses.Mean maize rootlet biomass in the upper 30 cm of the sloping cropland accounted for 84% of the total biomass.The root distribution of Bothriochloa ischaemum was well-proportioned, except for the lowest values in the 80-100 cm layer (Figure 3).The roots of Pinus tabuliformis were concentrated in the 40-60 cm layers, corresponding to the peaks of water loss in the profiles.

Composition of δ 18 O and δ 2 H in the Four Land Use Types
The Local Meteoric Water Line (LMWL) is given in Figure 6.The Local Meteoric Water Line (LMWL, Figure 6) was deviation to the Global Meteoric Water Line (GMWL) (δ 2 H = 8δ 18   water hydrogen and oxygen stable isotopic composition in the unsaturated zone in the four land use types in May, August, and November.Cross-plotting the 2-hydrogen and 18-oxygen isotopic records with different slopes suggests seasonally variable evaporative enrichment of 18-oxygen (the red, olive, and blue lines represent the isotopic profiles for May, August, and November, respectively).The slope of the regression line for soil water evaporation differed among the four land use types.The effect of soil evaporation in May distribution from big to small was as follows: Sloping cropland, shrubland, grassland, and forestland.For August, the list was sloping cropland, shrubland, grassland, and forestland.Meanwhile, the evaporation effect of soil in November obey the following distribution: Sloping cropland, grassland, shrubland, forestland.

Composition of δ 18 O and δ 2 H in the Four Land Use Types
The Local Meteoric Water Line (LMWL) is given in Figure 6.The Local Meteoric Water Line (LMWL, Figure 6) was deviation to the Global Meteoric Water Line (GMWL) (δ 2 H = 8δ 18   water hydrogen and oxygen stable isotopic composition in the unsaturated zone in the four land use types in May, August, and November.Cross-plotting the 2-hydrogen and 18-oxygen isotopic records with different slopes suggests seasonally variable evaporative enrichment of 18-oxygen (the red, olive, and blue lines represent the isotopic profiles for May, August, and November, respectively).The slope of the regression line for soil water evaporation differed among the four land use types.The effect of soil evaporation in May distribution from big to small was as follows: Sloping cropland, shrubland, grassland, and forestland.For August, the list was sloping cropland, shrubland, grassland, and forestland.Meanwhile, the evaporation effect of soil in November obey the following distribution: Sloping cropland, grassland, shrubland, forestland.

Composition of δ 18 O and δ 2 H in the Four Land Use Types
The Local Meteoric Water Line (LMWL) is given in Figure 6.The Local Meteoric Water Line (LMWL, Figure 6) was deviation to the Global Meteoric Water Line (GMWL) (δ 2 H = 8δ 18 O + 10).The values of δ 18 O varied between −16.74 and −2.87‰ with a mean of −8.54‰.The δ 18 O for a storm (60 mm) on 12 September 2016 was −15.19‰.The δ 18 O for a larger storm (80 mm) on 25 July was −8.31‰.Figure 6 show the soil water hydrogen and oxygen stable isotopic composition in the unsaturated zone in the four land use types in May, August, and November.Cross-plotting the 2-hydrogen and 18-oxygen isotopic records with different slopes suggests seasonally variable evaporative enrichment of 18-oxygen (the red, olive, and blue lines represent the isotopic profiles for May, August, and November, respectively).The slope of the regression line for soil water evaporation differed among the four land use types.The effect of soil evaporation in May distribution from big to small was as follows: Sloping cropland, shrubland, grassland, and forestland.For August, the list was sloping cropland, shrubland, grassland, and forestland.Meanwhile, the evaporation effect of soil in November obey the following distribution: Sloping cropland, grassland, shrubland, forestland.

Features of the Precipitation Isotopes
Combining the isotopic values for precipitation from January to December 2016 (Figures 2c and  6), the local meteoric water Lline (LMWL) were given in Figure 6.The precipitation isotopes from June to November were more depleted and were enriched in the study area from December to May, maybe because the precipitation from June to November came from maritime air masses, and the precipitation from December to May was predominantly from continental air masses [42,43].The LMWL from the study site deviated from to the GMWL, which is attributed to enrichment due to evaporation (Figure 6).The precipitation isotopes in our study varied seasonally and were strongly affected by rainfall (Figure 2).Heavy rainfall generally depleted the isotopic compositions.For example, δ 18 O for a storm (60 mm) on 12 September 2016 was −15.19‰.The δ 18 O for a larger storm (80 mm) on 25 July, however, was −8.31‰.The δ 18 O for this sample therefore represented δ 18 O for the initial rainfall carrying abundant heavy isotopes, other than those of the whole rainfall.

Effect of Evaporation on Soil Water
In our study area, evaporation is a major factor for the enrichment of 18-oxygen of soil water in shallow layers.A comparison of the regression lines for soil water with GMWL and LMWL, specifically δ 2 H = 8δ 18 O + 10 [44] and δ 2 H = 6.99δ 18 O − 3.25 (R 2 = 0.96, n = 50), indicated that the effect of evaporation varied among the land use types (Figure 6).The regression lines from four land use types deviated from the LMWL, which arises from enrichment of 18-oxygen of soil water because of evaporation.The δ 18 O differentially varied with each season for the four land use types (Figure 5).The non-equilibrium fractionation of evaporation made 18-oxygen of soil water enrichment around the soil surface.The isotopic profiles for May, August, and November indicated that most of the variation in δ 18 O was in the upper 200 cm for the four land use types, due to evaporation, infiltration, as well as mixing.Dry conditions caused rainwater to strongly evaporate, and confined the SWC recharge, as suggested by the negative peaks in May and August.The δ 18 O of soil water in the unsaturated zone is also affected by precipitation.So, the decline in the soil water isotopic profile is not as gradual as Zimmerman predicted [14].The δ 18 O varied slightly seasonally under the 200 cm soil layers (Figures 4 and 5), indicating that evaporation had a weak effect on deep soil water.Evaporation caused 18-oxygen enrichment in soil water around the surface, which created δ 18 O profiles that decreased in the 0-100 cm layers.The land use types affected δ 18 O in the soil profile significantly; this phenomenon was caused by the differences of the soil physical and hydrological properties.Correspondingly, the mechanism of infiltration of soil water in the unsaturated zone may be affected by the land use types.An ANOVA analysis of δ 18 O in the upper 100 cm indicated that soil water δ 18 O values were significantly in the grassland, shrubland, forestland, and sloping cropland (P = 0.015, 0.01, 0.01, 0.04, respectively).The spatial and temporal differences in soil water isotopes in the unsaturated zone (Figures 4 and 5) were attributed to evaporation.

Features of the Precipitation Isotopes
Combining the isotopic values for precipitation from January to December 2016 (Figures 2c and 6), the local meteoric water Lline (LMWL) were given in Figure 6.The precipitation isotopes from June to November were more depleted and were enriched in the study area from December to May, maybe because the precipitation from June to November came from maritime air masses, and the precipitation from December to May was predominantly from continental air masses [42,43].The LMWL from the study site deviated from to the GMWL, which is attributed to enrichment due to evaporation (Figure 6).The precipitation isotopes in our study varied seasonally and were strongly affected by rainfall (Figure 2).Heavy rainfall generally depleted the isotopic compositions.For example, δ 18 O for a storm (60 mm) on 12 September 2016 was −15.19‰.The δ 18 O for a larger storm (80 mm) on 25 July, however, was −8.31‰.The δ 18 O for this sample therefore represented δ 18 O for the initial rainfall carrying abundant heavy isotopes, other than those of the whole rainfall.

Effect of Evaporation on Soil Water
In our study area, evaporation is a major factor for the enrichment of 18-oxygen of soil water in shallow layers.A comparison of the regression lines for soil water with GMWL and LMWL, specifically δ 2 H = 8δ 18 O + 10 [44] and δ 2 H = 6.99δ 18 O − 3.25 (R 2 = 0.96, n = 50), indicated that the effect of evaporation varied among the land use types (Figure 6).The regression lines from four land use types deviated from the LMWL, which arises from enrichment of 18-oxygen of soil water because of evaporation.The δ 18 O differentially varied with each season for the four land use types (Figure 5).The non-equilibrium fractionation of evaporation made 18-oxygen of soil water enrichment around the soil surface.The isotopic profiles for May, August, and November indicated that most of the variation in δ 18 O was in the upper 200 cm for the four land use types, due to evaporation, infiltration, as well as mixing.Dry conditions caused rainwater to strongly evaporate, and confined the SWC recharge, as suggested by the negative peaks in May and August.The δ 18 O of soil water in the unsaturated zone is also affected by precipitation.So, the decline in the soil water isotopic profile is not as gradual as Zimmerman predicted [14].The δ 18 O varied slightly seasonally under the 200 cm soil layers (Figures 4 and 5), indicating that evaporation had a weak effect on deep soil water.Evaporation caused 18-oxygen enrichment in soil water around the surface, which created δ 18 O profiles that decreased in the 0-100 cm layers.The land use types affected δ 18 O in the soil profile significantly; this phenomenon was caused by the differences of the soil physical and hydrological properties.Correspondingly, the mechanism of infiltration of soil water in the unsaturated zone may be affected by the land use types.An ANOVA analysis of δ 18 O in the upper 100 cm indicated that soil water δ 18 O values were significantly in the grassland, shrubland, forestland, and sloping cropland (p = 0.015, 0.01, 0.01, 0.04, respectively).The spatial and temporal differences in soil water isotopes in the unsaturated zone (Figures 4 and 5) were attributed to evaporation.

Infiltration Mechanisms based on δ 18 O Isotopic Characteristics
The main scope of the δ 18 O of soil water changed by the different land use types is located in the 0-100 cm layers.The δ 18 O of soil water in the shrubland was −13.22‰ at a depth of 30-40 cm in November, and the precipitation δ 18 O was −13.65‰, respectively.Obviously, the results shown in Figure 5 indicated that the δ 18 O of soil water is consistent with the precipitation δ 18 O.The δ 18 O of soil water in the shrubland and precipitation isotopic values in November suggested that piston flow occurred.The δ 18 O of soil water in the 90-100 cm layer in the shrubland in November corresponded with δ 18 O in the 20-30 cm layer in August, indicating that soil water migration along the profile existed with a lag effect.The isotopic profiles of the forestland and sloping cropland in August contained abruptly lower values at depths of 160 and 50 cm, respectively, indicating that preferential flow occurred in these land use types, perhaps due to the paths for preferential flow provided by root holes, wormholes, cracks, and fissures in the loess [19].Due to this, macropores in the soil provided pathways for the rainwater during the percolation of soil water.Low SWC, especially during the rainy season (Figure 3), would facilitate the preferential flow [45].In the process of rainwater infiltration, the piston flow and the preferential flow coexist, which was largely impacted by the land use type.Briefly, other than the sloping cropland, precipitation on the Loess Plateau rarely recharged the deep soil (>200 cm), and the types of infiltration and land use primarily impacted the isotope characteristics of the soil profile.
Generally, the δ 18 O of soil water at the four land use types were more depleted below the 200 cm soil layer and suggest greater contribution of rainy season precipitation to replenishment.According to the Nielson and Bouma (1985) classification system: CV ≤ 10%, 10% < CV <100%, and CV ≥ 100% indicate low, middle, and strong variability, respectively [46].The CV values of soil water δ 18 O were −4.42% in the grassland, −5.60% in the shrubland, and −4.91% in the forestland.The CV value of soil water δ 18 O belongs to weak variability in the grassland and shrubland forestland below 160 cm, 180 cm, and 200 cm, respectively, indicating that normal rains and even continuous rains rarely affected the soil water in the deeper layers, except for the sloping cropland.The CV value of soil water δ 18 O belongs to strong variability in the sloping cropland throughout the year, indicating that it may have a large impact on groundwater recharge.The isotope depth profiles could be used to calibrate a numerical model and derive this way time-variant travel times, while the soil water isotope sampling is limited to snapshots into the processes (e.g., travel time) [47].

Tracing the Mixing of Soil Water Based on δ 18 O Isotopic Characteristics
The hydrogeochemical characteristics of the soil water provided evidence of mixing.The amplitude of the seasonal cyclic variation in soil water δ 18 O for each layer gradually decreased with depth (Figures 4 and 5), and the difference in δ 18 O between precipitation and soil water could not account for evaporation under the climatic conditions in our study area, where temperatures were lower (13 • C) and relative humidity was higher (95%) in autumn.The δ 18 O was most negative in the shrubland at −13.22‰, similar to −13.65‰ for the precipitation, further illustrating that the evaporation of precipitation may have been negligible in the enrichment during November.The δ 18 O was most negative in the forestland at −12.67‰, in the grassland at −12.91‰, and more enriched to −13.65‰ for the precipitation, indicating that the difference may be because antecedent soil water was combined by rainwater.The δ 18 O profile in November (Figure 5) probably denotes antecedent soil water in summer ("old water") and a mixture of precipitation in autumn ("new water").The mixing of soil water probably increased slightly with depth.The δ 18 O for the soil water in the grassland, shrubland, and forestland of the temporal variability decreases with depth and approaches a minimum value at 160 cm, 180 cm, and 200 cm, respectively, suggesting that the residence times of this water could be several months or more.Mixing can be caused by the infiltration or percolation of newer water into older soil water, hydrodynamic dispersion, or the interaction of preferential flow, among other causes.Mixing with "old water" could explain the isotopic enrichment of the precipitation in the profile and the effect of evaporation on mixing, suggesting that δ 18 O is useful for tracing the mixing of soil water in the unsaturated zone.Identifying the dominant factor(s) for the mixing based on the amount of water and the δ 18 O data for the precipitation and soil water, however, is difficult.

Tracing Soil Water Movement by Cross Plotting δ 2 H and δ 18 O
The results [24] indicated that the water absorbed by plant root would not fractionate the soil water isotopes.The absorption of water by the vegetation, however, would affect the water distribution in the soil profile and lead to changes in the patterns of soil water movement, which would affect the isotopic profile.It is suggested here that evaporation impact was changed with various seasons and diverse four land use types through the comparison of the regression line for soil water with the LMWL, in particular, δ 2 H = 6.99(δ 18 O) − 3.25 (R 2 = 0.96, n = 50).In May, there was a strong evaporation phenomenon in the sloping cropland.The lower precipitation (<50 mm) and the unexcavated slope value were lower, which could be expressed as 5.69.Furthermore, forestland with dense canopies and thick litter had a higher slope value of 6.56.Evapotranspiration became stronger in August in the four land use types.The SWC of forestland and grassland decreased with the decrease of root distribution.However, this phenomenon has not occurred in the sloping cropland (Table 3 and Figure 3).This finding requires further research.The root distribution of plants is closely related to SWC, and the water absorption layer is often related to the root distribution area [48].The forestland, which contain Pinus tabuliformis as the dominant plants, have a high root biomass at 40-60 cm depth, which agree with the water loss in the profiles.Stolte et al. [10] reported a similar result.In November, the regression line for the soil water in the grassland, shrubland, forestland, and sloping cropland were expressed as y = 5.44x − 19.96 (R 2 = 0.94), y = 5.45x − 20.03 (R 2 = 0.94), y = 6.33x − 7.73 (R 2 = 0.96), and y = 4.75x − 24.88 (R 2 = 0.86), respectively, suggesting weak evaporation.The higher slope values and SWC indicated that the compensation effect occurs during the late growing season of the four land use types (Figures 3 and 6).Yet similar phenomena have not been observed in forestland.The reason may be that the dense of foliage and canopy reduce the infiltration rate of precipitation moisture retained on the surface.This assumption is supported by the SWC in November reaching 22.48% (Figure 3).And reference [49] indicated that the connectivity can provide a basis for the development of better measurement and modeling approach.This research can be conducive to have a further understanding about water quantifying and sediment transferring in catchment systems.

Conclusions
Oxygen isotopic compositions and mechanisms of soil water movement in the unsaturated zone were studied in four land use types in the Wangmaogou watershed on the Loess Plateau.The critical factors that controlled how water moved through, out of, and into the soil profile were identified.The strong evaporation on the Loess Plateau was likely the dominant factor accounting for 18-oxygen enrichment in the shallow layer, whereas the characteristics of precipitation infiltration in the various seasons were the critical factors controlling the profile dynamics of δ 18 O in deeper soil.The δ 18 O of the soil water were in a steady state in the grassland, shrubland, and forestland were 160 cm, 180 cm, and 200 cm, respectively.The composition of the lighter isotopes during infiltration suggested that large mixing with "old water" followed the infiltration.The mixing may have been due to the infiltration or percolation of "new water" into the "old water", hydrodynamic dispersion, or the interaction of preferential flow, among other causes.
The characteristics of the δ 18 O isotopic composition profile indicated that piston flow and preferential flow coexisted in the infiltration of rainwater.The comparison of four land use types shows that the isotopic composition of soil water may vary due to local drainage and soil characteristics.The higher δ 18 O in the drier land uses may partly attribute to the high evaporation 18-oxygen enrichment and different proportions of "new water" and "old water".This study did not address the time of "new water" and "old water" mixing residence and nature.Future studies using combinations of analytical techniques would help address these problems.To sum up, the δ 18 O is suggested to be more sensitive than SWC for tracing the depth of precipitation infiltration and recharge mechanisms.

Figure 1 .
Figure 1.Location of the study area in Shaanxi Province, China (a,b); sampling sites (c).

Figure 1 .
Figure 1.Location of the study area in Shaanxi Province, China (a,b); sampling sites (c).

Figure 2 .
Figure 2. Daily variation of (a) temperature and precipitation and monthly variation of precipitation amount (b) and precipitation δ 18 O (c) at the study site.

Figure 2 .
Figure 2. Daily variation of (a) temperature and precipitation and monthly variation of precipitation amount (b) and precipitation δ 18 O (c) at the study site.
cropland retained water the best, with the highest SWCs throughout the study period.SWC in the 40-100 cm layers show a different value with the change of growing season.

Figure 3 .
Figure 3. Distribution of soil water content (SWC) in the four land use types.The standard deviation represents the dispersion degree of the data.

Figures 4 and 5
Figures 4 and 5 show the characteristics of δ 18 O for the soil water in the unsaturated zone in the four land uses in May, August, and November.The δ 18 O in May ranged from −10.77‰ to 4.85‰, with a mean of −7.61‰, and averaged −7.67‰ in the grassland, −7.76‰ in the shrubland, −7.83‰ in the forestland, and −7.19‰ in the sloping cropland.The δ 18 O in August ranged from −10.65‰ to 0.10‰, with a mean of −8.24‰, and averaged −8.33‰ in the grassland, −8.12‰ in the shrubland, −8.14‰ in the forestland, and −8.35‰ in the sloping cropland.The δ 18 O in November ranged from −13.23 to −2.85‰, with a mean of −9.14‰, and averaged −9.02‰ in the grassland, −9.33‰ in the shrubland, −8.55‰ in the forestland, and −9.65‰ in the sloping cropland.The decreasing trend of δ 18 O was interrupted by precipitation early in August and November and peaked again in November in the 80-90 cm layer.The δ 18 O was stable in deep soil (>200 cm) and varied seasonally the most in the 0-10 cm layer.The δ 18 O decreased significantly with depth and stabilized in the deep layers.The source of the soil water became depleted in 18-oxygen from May to November.The peak negative value was recorded at depths of 40 cm in the grassland and forestland, 50 cm in the shrubland and sloping cropland; rain fell on 8 and 12 November 2016 (30 mm).The δ 18 O of the soil

Figure 3 .
Figure 3. Distribution of soil water content (SWC) in the four land use types.The standard deviation represents the dispersion degree of the data.

Figures 4 and 5
Figures4 and 5show the characteristics of δ18 O for the soil water in the unsaturated zone in the four land uses in May, August, and November.The δ 18 O in May ranged from −10.77‰ to 4.85‰, with a mean of −7.61‰, and averaged −7.67‰ in the grassland, −7.76‰ in the shrubland, −7.83‰ in the forestland, and −7.19‰ in the sloping cropland.The δ 18 O in August ranged from −10.65‰ to 0.10‰, with a mean of −8.24‰, and averaged −8.33‰ in the grassland, −8.12‰ in the shrubland, −8.14‰ in the forestland, and −8.35‰ in the sloping cropland.The δ 18 O in November ranged from −13.23 to −2.85‰, with a mean of −9.14‰, and averaged −9.02‰ in the grassland, −9.33‰ in the shrubland, −8.55‰ in the forestland, and −9.65‰ in the sloping cropland.The decreasing trend of δ18 O was interrupted by precipitation early in August and November and peaked again in November in the 80-90 cm layer.The δ 18 O was stable in deep soil (>200 cm) and varied seasonally the most in the 0-10 cm layer.The δ 18 O decreased significantly with depth and stabilized in the deep layers.
water in the 0-10 cm layer in the four land use types indicated 18-oxygen enrichment, and 18-oxygen was less abundant in the deepest layers.The δ 18 O was relatively stable in the layers below 200 cm and varied little seasonally.Soil water δ 18 O ranged from −9.93‰ to −8.28‰ in the grassland, ranged from −10.11‰ to −8.07‰ in the shrubland, and ranged from −9.57‰ to −7.92‰ in the forestland under 160 cm, 180 cm, and 200 cm, respectively.

Figure 4 .
Figure 4. Variation of soil-water δ 18 O in the profiles of the four land use types.

Figure 5 .
Figure 5. Seasonal variation of soil-water δ 18 O in the profiles of the four land use types.

Figure 4 .
Figure 4. Variation of soil-water δ 18 O in the profiles of the four land use types.

Figure 4 .
Figure 4. Variation of soil-water δ 18 O in the profiles of the four land use types.

Figure 5 .
Figure 5. Seasonal variation of soil-water δ 18 O in the profiles of the four land use types.

Figure 5 .
Figure 5. Seasonal variation of soil-water δ 18 O in the profiles of the four land use types.The source of the soil water became depleted in 18-oxygen from May to November.The peak negative value was recorded at depths of 40 cm in the grassland and forestland, 50 cm in the shrubland and sloping cropland; rain fell on 8 and 12 November 2016 (30 mm).The δ 18 O of the soil water in the 0-10 cm layer in the four land use types indicated 18-oxygen enrichment, and 18-oxygen was less abundant in the deepest layers.The δ 18 O was relatively stable in the layers below 200 cm and little seasonally.Soil water δ 18 O ranged from −9.93‰ to −8.28‰ in the grassland, ranged from −10.11‰ to −8.07‰ in the shrubland, and ranged from −9.57‰ to −7.92‰ in the forestland under 160 cm, 180 cm, and 200 cm, respectively.

Figure 6 .
Figure 6.The soil water isotopic composition of the four land use types.The black lines represent the global meteoric water line (GMWL) of Carig (1961) shown for reference.The violet lines represent the local meteoric water line (LMWL).

Figure 6 .
Figure 6.The soil water isotopic composition of the four land use types.The black lines represent the global meteoric water line (GMWL) of Carig (1961) shown for reference.The violet lines represent the local meteoric water line (LMWL).

Table 1 .
Soil particle-size distribution and bulk density in the six land uses.

Table 1 .
Soil particle-size distribution and bulk density in the six land uses.

Table 3 .
Biomasses in the six land uses.
3.4.Profile Variations of δ 18 O in the Four Land Use Types

Table 3 .
Biomasses in the six land uses.
3.4.Profile Variations of δ 18 O in the Four Land Use Types