Variability of δ 2 H and δ 18 O in Soil Water and Its Linkage to Precipitation in an East Asian Monsoon Subtropical Forest Plantation

: The linkage between δ 2 H and δ 18 O of soil water and precipitation provides a way of understanding precipitation inﬁltration, residence time, and soil water source. Soil water at 0–5, 15–20, and 40–45 cm depths and event-based precipitation were collected in a subtropical forest plantation. Correlations between the δ 18 O of soil water and precipitation on the same day were used to determine the critical threshold of precipitation inﬁltration. Residence time of precipitation in soil was determined with correlations between the δ 18 O of soil water and cumulative precipitation before sampling. Soil water source was determined by the intersection points of Soil Water Evaporation Lines (SEL) and local meteoric water lines. The results showed that precipitation >5–6 mm could pass through canopy and litter, and inﬁltrate into soil. Residence times varied from a few days to several months, and increased with soil depth. The model-based approach for SEL estimation were more robust than the regression-based approach due to the inverse variability in the δ 2 H and δ 18 O of soil water source and soil evaporative fractionation. Soil water at a 0–5 cm depth originated mainly from precipitation in the current season, while those at 15–20 and 40–45 cm depths originated mainly from precipitation in the previous season. subtropical high and subtropical East Asian monsoon [23]. Annual precipitation and mean air temperature for the period of 1985–2020 were 1407.3 ± 300.2 mm and 18.1 ± 0.5 ◦ C, respectively, according to the meteorological records of CERN (Table 1). The whole year was divided into winter (December–February), spring (March–May), summer (June–August), and autumn (September–November) according to air temperatures [24,25]. More than half of the annual precipitation (68%) takes place in spring (32%) and summer (36%), while the mean air temperature peaked in summer (Table 1).


Introduction
Soil water, as a key link in the hydrological cycle, partitions precipitation into evaporation, transpiration, and runoff [1]. Precipitation is the main source of soil water replenishment, and it controls the mixing and redistribution of soil water [2]. The linkage of soil water with precipitation is increasingly relevant due to the changes in rain patterns caused by climate change [3]. δ 2 H and δ 18 O in water have been widely used as ideal tracers of water movement and mixing among different water pools [4]. δ 2 H and δ 18 O in soil water can largely retain the isotopic signal of precipitation due to the infiltration of precipitation larger than a certain amount, and by the lack of fractionation during root water uptake, except that by some halophytes or woody xerophytes [5,6]. Nevertheless, soil evaporation can result in the enrichment of soil water in heavy isotopes [7]. The linkage between δ 2 H and δ 18 O in soil water and precipitation is useful for the understanding of precipitation infiltration, residence time, and soil water source [8].
The replenishment of soil water occurs only via precipitation that exceeds the critical threshold-the amount that is lost to canopy interception or litter retention [9]. The determination of the critical threshold for the precipitation recharge of soil water is central to the accurate estimation of soil water replenishment and water balance [2,9]. Some authors have estimated this threshold with throughfall measurements using rain-gauges and determined that it corresponded to throughfall when precipitation was 0 based on linear regression between the two [10]. Others used the relationship between the variability in soil water content and precipitation amount to estimate this threshold [11]. However, δ 2 H and δ 18 O in precipitation directly influences those in soil water due to precipitation infiltration, and their linkage can be used for analyzing soil water replenishment such as the ecohydrological separation [12] or connectivity [13] between soil-bound and mobile water.
Knowledge of the residence time of precipitation in the soil profile helps to better understand hydrological processes and timescales of transport, and to improve hydrological models [14]. Tracers and hydrological models are usually used for estimating residence times of precipitation in soil [15]. Natural variation in δ 2 H and δ 18 O in soil water can be used to date water for up to 5 years, depending on the mixing and dispersion across the soil profile [4]. Hydrological models (e.g., HYDRUS-1D, which can simulate the onedimensional movement of water and multiple solutes in variably saturated porous media) can also be used to simulate residence times by tracking the travel time of water particles as they flow through a soil profile [16]. Although these models may allow researchers to link residence times with specific hydrological processes, they are based on the assumption that soil water is well mixed in pore space, which is in contrast to recent evidence that soil water can be distinguished into soil-bound and mobile water [1]. Existing studies showed that residence times of precipitation in soils typically span from hours to years due to the integrated influences of precipitation input, connectivity of soil pores, plants water uptake, and so on [15].
Soil water may be recharged by precipitation from the past rather than recent seasons due to the interactions between the mixtures of two soil water pools and variations in flow pathways, and the understanding of the seasonal origins of soil water remains limited [17]. Soil Water Evaporation Lines (SEL, the trend line of δ 2 H and δ 18 O in soil water in dualisotope space) deviated from Local Meteoric Water Lines (LMWL, a regression line between δ 2 H and δ 18 O values in precipitation) due to fractionation effects of soil evaporation [17,18]. To compensate for these fractionation effects, the intersection points of SEL and LMWL can be interpreted as the δ 2 H and δ 18 O in the soil water source [18,19]. SEL has been estimated by using theoretical model-based approaches and regression-based approaches with data from soil water samples [20]. Some authors reported that when isotopic variability in precipitation or the seasonal variability in evaporation fractionation of soil water were large, the SEL calculated with theoretical model-based approaches was more robust than that calculated with the regression-based approaches [21]. Using the Craig and Gordon model for the theoretical modeling of evaporative isotopic fractionation for diffusioncontrolled soil evaporation scenarios in estimating SEL, Xiang et al. [22] found that deep soil water in the Loess Plateau of China originated mainly from rainy season precipitation.
In this study, we analyzed δ 2 H and δ 18 O in soil water extracted from soil samples at 0-5, 15-20, and 40-45 cm depths, and event-based precipitation in an East Asian monsoon subtropical forest plantation. The critical threshold for the precipitation recharge of soil water was determined with correlations between δ 2 H and δ 18 O in soil water and precipitation on the same day. The residence time of precipitation in soil was determined with correlations between δ 2 H and δ 18 O in soil water and cumulative precipitation before soil sampling. Soil water source was determined by using the intersection points of LMWL and SEL calculated with regression-based and Craig and Gordon model-based approaches. The objectives of our study were: (1) to investigate how much precipitation could pass through the canopy and litter, and infiltrate into the soil; (2) to determine seasonal variability in the residence times of precipitation at different soil depths; (3) to compare the regression-based and Craig and Gordon model-based approaches for estimating SEL, and analyze the seasonal origins of soil water.

Study Site
This study was conducted at Qianyanzhou (QYZ) Ecological Experimental Station (26 •  subtropical high and subtropical East Asian monsoon [23]. Annual precipitation and mean air temperature for the period of 1985-2020 were 1407.3 ± 300.2 mm and 18.1 ± 0.5 • C, respectively, according to the meteorological records of CERN (Table 1). The whole year was divided into winter (December-February), spring (March-May), summer (June-August), and autumn (September-November) according to air temperatures [24,25]. More than half of the annual precipitation (68%) takes place in spring (32%) and summer (36%), while the mean air temperature peaked in summer (Table 1).

Study Site
This study was conducted at Qianyanzhou (QYZ) Ecological Experimental Station (26°44′52″ N, 115°39′47″ E, and elevation 102 m) of the Chinese Ecosystem Research Network (CERN), a member of ChinaFLUX, located in Taihe County, Jiangxi Province in southern China (Figure 1). The climate in the area is controlled by the Western Pacific subtropical high and subtropical East Asian monsoon [23]. Annual precipitation and mean air temperature for the period of 1985-2020 were 1407.3 ± 300.2 mm and 18.1 ± 0.5 °C, respectively, according to the meteorological records of CERN (Table 1). The whole year was divided into winter (December-February), spring (March-May), summer (June-August), and autumn (September-November) according to air temperatures [24,25]. More than half of the annual precipitation (68%) takes place in spring (32%) and summer (36%), while the mean air temperature peaked in summer (Table 1).  ipitation; Ta is air temperature; RH is relative humidity; Pa is atmospheric pressure; WS is wind speed; BD is soil sity; Ts is soil temperature; TD is tree density; TB is total biomass; LAImax is maximum of leaf area index.  P is precipitation; T a is air temperature; RH is relative humidity; P a is atmospheric pressure; WS is wind speed; BD is soil bulk density; T s is soil temperature; TD is tree density; TB is total biomass; LAI max is maximum of leaf area index.
The study plot (100 × 100 m 2 ) was at the top of a hill located within the Songtang catchment of QYZ (Figure 1), where soil depth is less than 100 cm, and average groundwater level is about 3.43 m. The soil at the site is red earth weathered from sandstone, sandy conglomerate, mudstone, and alluvium. Soil bulk density, sand, silt, and clay contents, and soil temperature are shown in Table 1 [26]. The subtropical forest plantation at the site was planted around 1985, and the dominant tree species were Masson pine (Pinus massoniana L.), slash pine (Pinus elliottii E.), and Chinese fir (Cunninghamia lanceolata L.) [27]. The growing season lasts from March to October. The tree density, the total biomass, and the maximum of leaf area index are shown in Table 1 [28].
The QYZ station (http://qya.cern.ac.cn, accessed on 26 August 2021) operates an above-canopy flux system for ecosystem evapotranspiration and CO 2 flux measurements. The flux system was mounted on a tower at 39.6 m, and consisted of an open-path CO 2 /H 2 O analyzer (LI-7500, Licor Inc., Lincoln, NE, USA) and a three-dimensional sonic anemometer (CSAT3, Campbell Scientific Inc., Logan, UT, USA). Flux variables were sampled at 10 Hz using a CR5000 datalogger (Model CR5000, Campbell Scientific Inc., Logan, UT, USA), and 30 min mean fluxes were calculated [28]. Air temperature and relative humidity were measured with a temperature and humidity sensor (HMP45C, Vaisala Inc., Helsinki, Finland). Precipitation was monitored with a rain gauge (RG13H, Vaisala Inc., Helsinki, Finland). A pressure sensor (DPA501, Delta Electronics Inc., Taibei, China) was used for measuring atmospheric pressure. Wind speed was monitored with an anemoscope (WAA151, Vaisala Inc., Helsinki, Finland). These meteorological variables were collected every 1 h. Precipitation was collected in a polyethylene bottle fitted with a funnel and topped with a ping-pong ball to prevent evaporation [26]. The end of a precipitation event was identified as the point after 6 h of no precipitation from the moment this precipitation had stopped. Precipitation samples were collected after each precipitation event from January 2011 to December 2017. All precipitation samples were refrigerated at 4 • C before isotope analysis.

Sample Collection and Measurements
In total, 3514 soil samples and 520 precipitation samples were collected during the study period. A week after sampling, soil water was extracted from soil samples with a cryogenic vacuum distillation system, with heating at >90 • C and an extraction time of 0.5-1.5 h, depending on soil water content [29]. In order to determine soil water extraction efficiency, pre-and post-water extraction soil sample weights as well as sample weights after additional oven-drying (105 • C, 48 h) were compared [29]. In this study, the extraction efficiency of water from soil samples was higher than 98.0%. Soil water and precipitation were filtered through a 0.45 µm mixed cellulose membrane (Jiuding Gaoke Co. Ltd., Beijing, China), and 2 mL of water samples was used for the analysis of δ 2 H and δ 18 O. δ 2 H and δ 18 O were analyzed using an Isotopic Ratio Infrared Spectroscopy (IRIS) system (Model DLT-100; Los Gatos Research Inc., San Jose, CA, USA) [30]. The number of injections into IRIS were 6, and the results of the last three injections were used for analysis. The results were normalized to Vienna Standard Mean Ocean Water (VSMOW), and expressed in the standard δ notation (in‰): where δ sample is the δ 2 H and δ 18 O of the sample, and R sample and R VSMOW are the ratio of 18

. Critical Thresholds for Precipitation Recharge of Soil Water
Same-day soil water and precipitation samples (n = 94) were selected in order to determine the critical thresholds for the precipitation recharge of soil water-when precipitation larger than this threshold can pass through the canopy and litter, and infiltrate into the soil. Correlation coefficients between δ 2 H and δ 18 O in selected soil water at each depth (0-5, 15-20, and 40-45 cm) and precipitation were calculated. Precipitation was classified into different precipitation amounts, including ≥0.2, 0.4, 0.6, 0.8, 1, 1.2, 1.6, 2, 2.4, 2.8, 4, 5, 6, 8, 11, and 15.5 mm. The threshold was determined to be the precipitation amount when correlation coefficients between δ 2 H and δ 18 O in soil water and precipitation that is larger than the precipitation amount were more than the coefficients between δ 2 H and δ 18 O in soil water and precipitation less than the precipitation amount, based on the graphical inference or statistical methods.
There were no significant differences (p > 0.05) among the triplicates of soil water samples according to multiple comparisons. Weighted mean δ 2 H and δ 18 O for soil water (δ s,mean ) was obtained by weighting the soil water content of each sample, and using the following equations: where θ i is the soil water content of the ith soil sample replicate; δ s,i is the δ 2 H and δ 18 O of the ith soil sample replicate; i is the number of the soil sample replicate.

Determining Residence Times of Precipitation in Soil
Correlation coefficients between δ 2 H and δ 18 O in soil water and the cumulative precipitation before soil sampling were calculated. Only precipitation larger than the critical thresholds for the precipitation recharge of soil water were analyzed. Cumulative precipitation was calculated for periods of 0, 2, 7, 15, 30, 45, 60, 75, 90, 105, 120, 135, and 150 days before sampling due to the collection frequency of soil at an average of 2 or 8-12 times per month. The residence time of precipitation after entering the soil profile was determined as that period of time (from 0 to 150 days) when correlation coefficients were the highest. Residence times at 0-5, 15-20, and 40-45 cm depths during winter, spring, summer, and autumn, respectively, were calculated.
The weighted means of δ 2 H and δ 18 O for precipitation (δ p,mean ) were obtained by weighting the amount of precipitation for each precipitation event, and calculated as: where δ p,i is the δ 2 H and δ 18 O of the ith precipitation event, PPT i is the amount of the ith precipitation event. The precipitation samples that significantly deviated from the LMWL were eventually excluded. Machine learning [31] was used for predicting the missing δ 2 H and δ 18 O of the precipitation event after data quality control. We used random forests to train data and learn relationships between predictors and outcomes. Precipitation, air temperature, vapor concentration, atmospheric pressure, relative humidity, atmospheric vapor pressure deficit (VPD), wind speed, evapotranspiration (ET), and outgoing longwave radiation value (OLR, downloaded from NOAA Climate Data Record) were selected as predictors for the training of machine learning models. The calendar seasons were added as fuzzy sets in the model. A total of 80% of the precipitation event samples was randomly selected as the training subset, and the remaining 20% of the samples was selected as the validation subset. The Root Mean Squared Errors (RMSE) of prediction using random forests were 16.90 for δ 2 H, and 2.04 for δ 18 O, respectively. There was no significant difference between the magnitude of results obtained by predicting the missing data and only the measured data. In order to determine if the regression-based approach was valid, SEL was calculated by the regression-based and Craig and Gordon model-based approaches. The regressionbased approach is defined as a linear regression between the δ 2 H and δ 18 O values in soil water taken from a depth profile [32]: where a SEL and b SEL represent the slope and intercept of SEL.
For the Craig and Gordon model-based approach, the slope of SEL (S SEL ) in a soil water pool fed by local precipitation is predicted by the Craig and Gordon model [33]: where h is the relative humidity, set to the mean relative humidity during the residence time; δ P is the δ 2 H and δ 18 O in precipitation, set to the weighted average of δ 2 H and δ 18 O in precipitation during the residence time; δ A is the δ 2 H and δ 18 O in ambient atmospheric vapor, and is determined using the precipitation-equilibrium assumption: α + is the liquid-vapor equilibrium isotopic fractionation, estimated from empirical relationships by Horita where T s is soil temperature at the evaporation front, set to the mean air temperature during the residence time. ε + is the equilibrium isotopic separation between liquid and vapor, calculated as ε + (‰) = (α + − 1) × 10 3 . ε K is the equivalent kinetic isotopic separation based on wind tunnel experiments [35], calculated as: where C 0 K is 25.0 and 28.6% for deuterium and oxygen-18, respectively; n = 1 for soil water; is an advection term to account for the potential influence of humidity, often set to θ ≈ 1 for small water bodies, and h' is the adjusted humidity of the downwind atmosphere following the admixture of evaporating moisture over the surface.
Intercept of SEL (I SEL ) is calculated with δ 2 H and δ 18 O in soil water (δ 2 H s and δ 18 O s ) as follows: (2) δ 2 H and δ 18 O in Soil Water Source δ 2 H and δ 18 O in soil water source were calculated using the intersection points of SEL with LMWL [19]: where a SEL and b SEL represent the slope and intercept of SEL, and a and b represent the slope and intercept of LMWL. δ 2 H and δ 18 O in precipitation >5 mm were used for fitting LMWL according to the following analysis of the critical thresholds for the precipitation recharge of soil water. The slopes or intercepts of LMWL in spring, summer, autumn, and winter were calculated using an ordinary least squares regression (OLSR). However, the precipitation amount weighted reduced major axis regression (PWRMA) may be more suitable for calculating LMWL, according to [36]. OLSR, reduced major axis regression (RMA), major axis regression (MA), and the corresponding precipitation weighted regressions (PWLSR, PWRMA, and PWMA) were compared by analyzing the calculated slopes and intercepts of LMWL, the average Root Mean Sum of Squared Error value of the fit (RMSE av ), and the significance in the difference of each regression to that of OLSR, according to [36] (Table S1). Although the slopes produced by OLSR were slightly less than those determined by PWRMA with the lowest RMSE av , there was no significant difference between those produced by OLSR and PWRMA (p > 0.05).

Statistical Analyses
Linear regression was used for analyzing correlations between δ 2 H and δ 18 O in soil water and precipitation with the SPSS 22.0 software. For the comparison of δ 2 H and δ 18 O in soil water or soil water source at different soil depths with δ 2 H and δ 18 O in precipitation, a multiple comparison Tukey's post hoc test was used with the SPSS 22.0 software. p < 0.05 was considered statistically significant at a 95% confidence interval. Pearson correlation coefficients were calculated using the R studio 3.5.0.  (Table S2). In addition, with an increase in the amount of precipitation, the weighted mean δ 18 O in precipitation gradually declined from −5.23‰, −5.05‰, −5.68‰, −6.50‰, −6.06‰ to −7.55‰ when the amount of precipitation was 0-5, 5-10, 10-20, 20-30, 30-40, and >40 mm/day, respectively (Table S2). The weighted mean δ 18 O in precipitation for large rainfall events (>40 mm/day) was significantly more depleted than that of small rainfall events (<10 mm/day).

Seasonal Variability of δ 2 H and δ 18 O in Soil Water and Precipitation
Seasonal variability of the weighted average δ 18 O in soil water at a depth of 0-5 cm, with −4.37‰, −7.94‰, −6.61‰, and −6.00‰ during spring, summer, autumn, and winter, respectively, was similar to that in precipitation ( Figure 2 and Table S2). The weighted average δ 18 O in soil water at a depth of 15-20 cm was most enriched in spring at −4.62‰, followed by −7.83‰, −8.13‰, and −8.38‰ in winter, summer, and autumn, respectively. They were similar to δ 18 O in precipitation during spring and summer, and more depleted than δ 18 O in precipitation during autumn and winter (Figure 3). At a depth of 40-45 cm, the weighted average δ 18 O in soil water was most enriched in spring at −5.40‰, followed by −7.20‰, −8.02‰, and −8.15‰ in summer, winter, and autumn, respectively. They were more enriched during summer, and more depleted during spring, autumn, and winter as compared with δ 18 O in precipitation ( Figure 3). However, in autumn and winter, the weighted average δ 18 O in soil water at a depth of 0-5 cm was significantly more enriched than that at 15-20 and 40-45 cm depths (Table S2). In spring, the weighted average δ 18 O in soil water at 0-5 and 15-20 cm depths were significantly more enriched than that at the 40-45 cm depth. Nevertheless, the weighted average δ 18 O in soil water at depths of 0-5 and 15-20 cm in summer was significantly more depleted than that at the 40-45 cm depth. more enriched during summer, and more depleted during spring, autumn, and winter as compared with δ 18 O in precipitation ( Figure 3). However, in autumn and winter, the weighted average δ 18 O in soil water at a depth of 0-5 cm was significantly more enriched than that at 15-20 and 40-45 cm depths (Table S2). In spring, the weighted average δ 18 O in soil water at 0-5 and 15-20 cm depths were significantly more enriched than that at the 40-45 cm depth. Nevertheless, the weighted average δ 18 O in soil water at depths of 0-5 and 15-20 cm in summer was significantly more depleted than that at the 40-45 cm depth.

Critical Thresholds of Precipitation for Infiltration into the Soil Profile
Correlation coefficients between δ 18 O in same-day soil water and precipitation larger than 6 mm were higher than those between δ 18 O in same-day soil water and precipitation less than 6 mm, at 0-5, 15-20, and 40-45 cm soil depths (Figure 4). This indicates that precipitation larger than 6 mm passed through the canopy and litter layer, infiltrated into the soil, and significantly changed δ 18 O in soil water. A critical threshold of precipitation for infiltration into the soil profile was determined to be 5-6 mm, and there were insufficient data to calculate correlation coefficients between δ 18 O in same-day soil water and

Critical Thresholds of Precipitation for Infiltration into the Soil Profile
Correlation coefficients between δ 18 O in same-day soil water and precipitation larger than 6 mm were higher than those between δ 18 O in same-day soil water and precipitation less than 6 mm, at 0-5, 15-20, and 40-45 cm soil depths (Figure 4). This indicates that precipitation larger than 6 mm passed through the canopy and litter layer, infiltrated into the soil, and significantly changed δ 18 O in soil water. A critical threshold of precipitation for infiltration into the soil profile was determined to be 5-6 mm, and there were insufficient data to calculate correlation coefficients between δ 18 O in same-day soil water and precipitation within the precipitation amount interval of 5-6 mm. Correlation coefficients between δ 18 O in same-day soil water and precipitation increased with an increase in precipitation amount, and reached their largest values when precipitation amounts were ≥11 mm. The peak values of correlation coefficients decreased with an increase in soil depth.

Critical Thresholds of Precipitation for Infiltration into the Soil Profile
Correlation coefficients between δ 18 O in same-day soil water and precipitation larger than 6 mm were higher than those between δ 18 O in same-day soil water and precipitation less than 6 mm, at 0-5, 15-20, and 40-45 cm soil depths (Figure 4). This indicates that precipitation larger than 6 mm passed through the canopy and litter layer, infiltrated into the soil, and significantly changed δ 18 O in soil water. A critical threshold of precipitation for infiltration into the soil profile was determined to be 5-6 mm, and there were insufficient data to calculate correlation coefficients between δ 18 O in same-day soil water and precipitation within the precipitation amount interval of 5-6 mm. Correlation coefficients between δ 18 O in same-day soil water and precipitation increased with an increase in precipitation amount, and reached their largest values when precipitation amounts were ≥ 11 mm. The peak values of correlation coefficients decreased with an increase in soil depth.

Seasonal Variability in Residence Times of Precipitation in Soil
In spring, correlation coefficient between δ 18 O in soil water at the 0-5 cm depth and δ 18 O in cumulative precipitation (R 0-5 ) was highest at 0.89 at 0 days, indicating that δ 18 O in soil water at the 0-5 cm depth was dominated by δ 18 O in precipitation on that sampling day ( Figure 5 0.88 at 30, 60, and 105 days, respectively. Residence times of precipitation in soil in autumn were more than those in winter at depths of 0-5, 15-20, and 40-45 cm. Nevertheless, peak R0-5, R15-20, and R40-45 at 0.67, 0.70, and 0.41 in 30, 60, and 105 days were low in summer as compared with those in other seasons, and it may be due to the combined influences of seasonal precipitation, soil evaporation, and plant transpiration. Residence times of precipitation in soil water increased with an increase in soil depth.

Seasonal Variability of δ 2 H and δ 18 O in Soil Water Source
In winter, slopes of SEL calculated with the Craig and Gordon model-based approach at the 0-5 cm depth (CG 0-5) were significantly lower than those at depths of [15][16][17][18][19][20] and [40][41][42][43][44][45] cm ( Figure 6). Slopes of CG 0-5, CG 15-20, and CG 40-45, which were 4.27, 5.27, and 5.65, respectively, were significantly lower than those calculated with the regression-based approach (7.05) in winter. Slopes of CG 0-5, CG 15-20, and CG 40-45, which were 4.13, 4.34, and 4.40, respectively, did not substantially differ, and were significantly lower than that calculated with the regression-based approach (7.00) in spring. Slopes of CG 0-5, CG 15-20, CG 40-45 and calculated with the regression-based approach in summer, which were 5.61, 5.94, 6.34, and 7.71, respectively, were similar to those in winter. In autumn, slopes of CG 0-5 (5.46) were significantly lower than those of CG 15-20 (5.90) and CG 40-45 (6.35) and calculated with the regression-based approach (6.47), and the latter three had no significant difference. A lower slope indicates a more enriched δ 2 H and δ 18 O in soil water source; δ 2 H and δ 18 O in the soil water source calculated with the two approaches for 0-5, 15-20, and 40-45 cm depths had opposite patterns (Table S3), except in autumn, when they were similar. Therefore, we conclude that both of the approaches can be used for calculating δ 2 H and δ 18 O in soil water source in autumn.
Based on the relationships between δ 18 O in soil water source and precipitation ( Figure 7 and Table S3) those of CG 15-20 (5.90) and CG 40-45 (6.35) and calculated with the regression-based approach (6.47), and the latter three had no significant difference. A lower slope indicates a more enriched δ 2 H and δ 18 O in soil water source; δ 2 H and δ 18 O in the soil water source calculated with the two approaches for 0-5, 15-20, and 40-45 cm depths had opposite patterns (Table S3), except in autumn, when they were similar. Therefore, we conclude that both of the approaches can be used for calculating δ 2 H and δ 18 O in soil water source in autumn.

Responses of Soil Water to Rainfall Events
Precipitation >5-6 mm passed through the canopy and litter layer, and infiltrated into soil in this subtropical plantation (Figure 4). This precipitation amount was slightly higher than that reported by Dai et al. [11], who showed that precipitation >3.7 mm led to a change in surface soil water content based on the regression equation of the variations in the soil water content from two consecutive days and daily precipitation. It may be due to more LAI (5.6 m 2 /m 2 ) in this subtropical plantation than that (3.6 m 2 /m 2 ) reported by Dai et al. [11]. However, Mello et al. [2] found that precipitation >1.58 mm could pass through the canopy layer (5.25 m 2 /m 2 ) based on the linear regression of throughfall and precipitation, without considering the litter layer. Existing studies have found that higher amounts of precipitation increased precipitation infiltration depths in grasses with a soil bulk density of 1.13-1.26 g/cm 3 [37], while the thresholds of precipitation for infiltration

Responses of Soil Water to Rainfall Events
Precipitation >5-6 mm passed through the canopy and litter layer, and infiltrated into soil in this subtropical plantation (Figure 4). This precipitation amount was slightly higher than that reported by Dai et al. [11], who showed that precipitation >3.7 mm led to a change in surface soil water content based on the regression equation of the variations in the soil water content from two consecutive days and daily precipitation. It may be due to more LAI (5.6 m 2 /m 2 ) in this subtropical plantation than that (3.6 m 2 /m 2 ) reported by Dai et al. [11]. However, Mello et al. [2] found that precipitation >1.58 mm could pass through the canopy layer (5.25 m 2 /m 2 ) based on the linear regression of throughfall and precipitation, without considering the litter layer. Existing studies have found that higher amounts of precipitation increased precipitation infiltration depths in grasses with a soil bulk density of 1.13-1.26 g/cm 3 [37], while the thresholds of precipitation for infiltration into soil at 0-5, 15-20, and 40-45 cm depths exhibited no significant differences in our study. However, the critical thresholds of precipitation for infiltration into the soil were closely related to vegetation and soil characteristics, and more LAI may cause a larger threshold.
Correlation coefficients between same-day δ 18 O in soil water and precipitation increased with the increase in the amount of precipitation, and reached peak values when the amount of precipitation was 11 mm (Figure 4). This indicates that precipitation may mainly recharge soil-bound water before the formation of mobile water, or that hydrologic connectivity occurred, which was defined as the mixture between soil-bound water and soil mobile water [38], and the connectivity may increase with an increase in precipitation [39,40] when precipitation amount was less than a certain value. In addition, we found that peak correlation coefficients decreased with an increase in soil depths (Figure 4). One reason for this may be that before the formation of mobile water, precipitation infiltrating into the soil profile decreases due to an increase in initial soil water content with soil depths [41,42]. Another reason may be that there is little connectivity between soil-bound and mobile water due to ecohydrologic separation after the formation of mobile water, or some preferential flows when the amount of precipitation is larger than a certain value [40,43].
Nevertheless, it should be noted that the processes of evaporation, mixture with soil water having different signals of δ 2 H and δ 18 O, and preferential flow may change the δ 2 H and δ 18 O of precipitation during the infiltration and influence the analysis of the precipitation threshold using the linkage between δ 2 H and δ 18 O in soil water and precipitation [8]. In addition, the incomplete extraction of soil water when soils contain a large percentage of silt and clay particles, or evaporation during the storage of soil samples, may cause certain uncertainties in the results [4].

Influencing Factors of Residence Times of Precipitation in Soil
In our study, residence times of precipitation in soil varied between a few days and several months, increased with the increase in soil depths, and were the shortest in spring and largest in autumn ( Figure 5). It was possible that the connectivity of soil pores for precipitation infiltration, soil evaporation, and plant transpiration decreased with the increase in soil depth [15]. For a depth of 0-5 cm, the integrative effects of evapotranspiration and precipitation infiltration may result in a short residence time (<30 days). The seasonal variation of residence times at 15-20 and 40-45 cm depths may be caused by the pattern of seasonal precipitation in the East Asian monsoon region [11]. Precipitation inputs are the fundamental supply of "new water" in soil, and thus its variability can drive variations in the residence times of old water [44,45]. High precipitation amounts in spring in our study area resulted in short residence times in the soil in spring due to a large influx of "new water" and high soil water recharge, which could be retained until autumn. This was not the case in winter.
Nevertheless, the δ 2 H and δ 18 O in soil water was not only determined by the isotopic signals of precipitation inputs, but also by the mixture between the two and soil evaporation [46]. Therefore, the best correlation coefficients may not provide correct information about the residence time. In our study, R0-5, R15-20, and R40-45 in summer were low, and this may have been due to the influence of high soil evaporation. In addition to this correlation coefficient method, other approaches such as the simple classical sine curve or the new ensemble hydrograph separation approaches can also be used to determine the mean transit times [15].

Seasonal Origins of Soil Water
Slopes of SEL calculated with a regression-based approach were always higher than those calculated with the Craig and Gordon model-based approach in spring, summer, and winter, while those slopes of SEL calculated with both approaches for autumn differed very little ( Figure 6). An inverse variability in δ 2 H and δ 18 O in soil water source and soil evaporation fractionation may lead to higher slopes of SEL calculated with the regression-based approach than that calculated with the Craig and Gordon model-based approach [11,21,22,47]. In winter and spring, δ 2 H and δ 18 O in soil water source at the 0-5 cm depth were more enriched than those at depths of 15-20 and 40-45 cm due to recharge by enriched winter and spring precipitation, while the evaporation fractionation of soil water at 15-20 and 40-45 cm depths was higher than that at a depth of 0-5 cm due to the retention of autumn precipitation, which experienced high evaporation fractionation. In summer, the decrease in the depleted summer precipitation recharge of soil water led to a gradual enrichment of soil water source with soil depth, while the evaporation fractionation of δ 2 H and δ 18 O in soil water gradually decreased. Nevertheless, in autumn, both δ 2 H and δ 18 O in soil water source and evaporation fractionation at the 0-5 cm depth were more enriched than those at 15-20 and 40-45 cm depths. Above all, the occurrence of soil samples with depleted soil water source and high evaporation fractionation, and those with enriched soil water source and low evaporation fractionation along a soil profile may cause the overestimation of SEL slopes with the regression-based approach. Therefore, the SEL slopes calculated with the Craig and Gordon model-based approach were more robust than those calculated with the regression-based approach.
In winter, soil water at the 0-5 cm depth originated mainly from autumn and winter precipitation, while that at 15-20 and 40-45 cm depths originated primarily from summer and autumn precipitation (Figure 7). In spring, soil water at 0-5 and 15-20 cm depths originated mainly from winter and spring precipitation, while soil water at a depth of 40-45 cm only originated from winter precipitation. In summer, soil water in all layers was derived mainly from current precipitation. In autumn, soil water at the 0-5 cm depth originated mainly from autumn precipitation, while soil water at 15-20 and 40-45 cm depths originated mainly from summer and autumn precipitation. We conclude that soil water at the 0-5 cm depth originated primarily from precipitation in the current season, while soil water at 15-20 and 40-45 cm depths originated mainly from precipitation in the previous season in all seasons, except summer. Seasonal origins of soil water were similar to seasonal patterns of residence times. The widespread presence of summer and autumn precipitation in deep soils in winter indicates that these waters often resided in soils for several months [20,47].

Conclusions
The results of this study indicate that precipitation in this area need to be larger than 5-6 mm to pass through the canopy and litter layer, and then infiltrate into the soil. When precipitation amount was less than 11 mm, precipitation may mainly recharge soil-bound water before the formation of mobile water, or hydrologic connectivity between soil-bound and mobile water occurred, and connectivity increased with increasing precipitation. Residence times of precipitation in soil varied between a few days and several months, increased with soil depth, and were the shortest in spring and largest in autumn. These results may be due to the connectivity of soil pores for precipitation infiltration, soil evaporation, plant transpiration, and the seasonal pattern of precipitation in the East Asian monsoon region. Slopes of SEL calculated with the Craig and Gordon model-based approach were more robust than those calculated with the regression-based approach due to the inverse variation in δ 2 H and δ 18 O in soil water source and soil evaporation fractionation. Nevertheless, the SEL slopes calculated with the two approaches differed little in autumn. δ 2 H and δ 18 O in soil water source and precipitation indicated that soil water at the 0-5 cm depth originated mainly from precipitation in the current season, while soil water in all seasons except summer at depths of 15-20 and 40-45 cm originated mainly from precipitation in the previous season. Our results highlight that precipitation in the previous season is important for alleviating the decreasing water availability during the autumn and winter seasons in the East Asian monsoon region.