Modeling and Evaluating Soil Salt and Water Transport in a Cultivated Land–Wasteland–Lake System of Hetao, Yellow River Basin’s Upper Reaches

: With the implementation of water-conservation projects in the Hetao Irrigation District (HID), great changes have taken place in the ecohydrological process. A cultivated land–wasteland– lake system in the upper Yellow River Basin (YRB) was chosen to study the soil salt and water transport process with the Hydrus-1D model. The model parameters were calibrated and validated by measuring the soil salt and water data. Measured values were in good agreement with the simulated values. The results showed, in the whole growth period, the deep percolation of cultivated land was 34–40% of the total applied water (rainfall and irrigation). The capillary rise in the cultivated land, wasteland, and lake boundary was 24%, 29–35%, and 62–68% of their own evapotranspiration, respectively. The capillary rise in the lake boundary was about 2 times that of the wasteland and 2.6 times that of the cultivated land. The salt storage in the 1 m soil zone of the lake boundary was more than 10% and 18% greater than that of the wasteland and cultivated land, respectively. The salt of the capillary rise in the lake boundary exceeded that of the wasteland by a factor of three. The salt accumulation in the upper soil zone of the cultivated land, wasteland, and lake boundary was 13%, 37%, and 48%. Soil salinization in the upper soil zone of the wasteland and lake boundary was serious, and some measures should be taken to reduce the salt content to prevent soil salinization. The results act as a theoretical basis for the ecohydrological control of the HID.


Introduction
Water deficits and soil salinization are the two major limitations on agricultural production and the ecological environment [1]. Irrigation water is seriously scarce in the semi-arid and arid areas of the world, and approximately one-fifth of irrigated agriculture is adversely affected by soil salinity [2]. Agricultural areas that are located in arid and semi-arid irrigation districts occupy the typical characteristics of shallow groundwater depth and fragmented land cover, for example, the Yellow River Basin (YRB) irrigation area, China [3], Central Asia's Fergana Valley [4], and India's Indo-Gangetic Plain [5], where the agro-hydrological cycles are extremely complex processes associated with soil water, salt migration, atmosphere, groundwater dynamics, and vegetation growth. Understanding these agro-hydrological processes is extremely significant for the efficient use of water resources and the control of salinization.
The Hetao Irrigation District (HID), in the upper reaches of the YRB, China, is a large salinized irrigation district. Due to poor irrigation and drainage practices, a large Hydrus-1D model in the semi-arid coastal plain of Korba in Tunisia [22]. Liu et al. (2019) applied the Hydrus-1D model to evaluate soil salt transport in wheat-maize cultivated lands with different texture layers in the lowland of the North China Plain [23]. However, the Hydrus-1D model has rarely been used to evaluate and analyze the internal relationship between salt and water dynamics in a cultivated land-wasteland-lake system in the HID.
Therefore, the primary aims of this study were to: (1) understand the mechanism of soil water consumption and salt accumulation; (2) reveal the relationship between soil salt and water in a cultivated land-wasteland-lake system; and (3) use the Hydrus-1D model, calibrated by field observation data, to evaluate the soil salt and water balance and quantitatively estimate the soil salt and water in the typical period. The results of this study provide a theoretical basis for the effective utilization of water resources and soil salinization control in the HID.

The Study Area
The HID is approximately 50 km long from north to south and 250 km wide from east to west ( Figure 1). The HID extends over 1.12 Mha; the cultivated land area is 5.74 × 10 5 hm 2 , the wasteland area is 2.10 × 10 5 hm 2 , and the water area is 1.30 × 10 4 hm 2 [5,24]. Wasteland and lakes are distributed in patches, and the formation of lakes mainly depends on the irrigation water. A large amount of salt has accumulated in the low-lying wetland [9]. Ground slope is from 1/5000 to 1/8000 from the west to the east and from 1/4000 to 1/8000 from the south to the north. The depth of groundwater varies between 0.5 and 3 m over the whole year [25]. Hetao has the characteristics of a continental arid and semi-arid climate. Annual precipitation is about 50-144.2 mm. The rainy season extends from June to August, accounting for 70% of annual precipitation. The study area has a mean annual temperature of 7.5 • C and the average annual evaporation is 2237 mm. Maize and sunflower are the main crops in Hetao.  (c) study area in red [26]. The study area is located at the Zhangliansheng site (40 • 54 36 N, 107 • 15 59 E), which is east of the Jiefangzha Irrigation District (JID) in Hetao ( Figure 1). The areas of the cultivated land, wasteland, and lake are 3.7 × 10 4 m 2 , 4.4 × 10 4 m 2 , and 5.12 × 10 5 m 2 , respectively ( Figure 2). The results of a real-time kinematic (RTK) showed that the maximum ground height difference in the farmland was 15 cm, and the average height of the farmland was 45 cm higher than that of the wasteland. The changes in the groundwater depth during the crop growth period were 50-220 cm for the cultivated land, 70-160 cm for the wasteland, and 30-140 cm for the lake boundary. The main crop in the study area is sunflower. The soil texture of the 300 cm soil zone in the study area is mainly sand and sandy loam (Table 1).

Meteorological Data
The daily meteorological data, including temperature, relative humidity, sunshine duration, wind speed, and rainfall, were collected. Based on the daily meteorological data and the Penman-Monteith equation, the daily reference evapotranspiration (ET0) was calculated [27]. The ET0 and rainfall in 2018 and 2019 are shown in Figure 3. The precipitation in 2018 and 2019 was 124 mm and 78 mm, respectively. The precipitation in the crop  The daily meteorological data, including temperature, relative humidity, sunshine duration, wind speed, and rainfall, were collected. Based on the daily meteorological data and the Penman-Monteith equation, the daily reference evapotranspiration (ET 0 ) was calculated [27]. The ET 0 and rainfall in 2018 and 2019 are shown in Figure 3. The precipitation in 2018 and 2019 was 124 mm and 78 mm, respectively. The precipitation in the crop growth period (1 June to 1 October) was 109.6 mm and 64.88 mm ( Figure 3). The average annual daily evapotranspiration was 2.7 mm/d. In the crop growth period, the average daily evapotranspiration was from 4.3 mm/d to 4.6 mm/d. growth period (1 June to 1 October) was 109.6 mm and 64.88 mm ( Figure 3). The average annual daily evapotranspiration was 2.7 mm/d. In the crop growth period, the average daily evapotranspiration was from 4.3 mm/d to 4.6 mm/d.

Groundwater Data
Seventeen groundwater observation wells were established in the study area, including seven key observation wells (in red circle) and ten regularly observed wells (in blue circle). Automatic groundwater sensors (CTD-10, Meter Company, San Francisco, CA, USA) were installed in key observation wells, and an EM50 logger (Meter Company) was used to record groundwater table, temperature, and electrical conductivity (EC) data in one-hour intervals. The groundwater depth of the regularly observed wells was measured every 7 days. Three key observation wells, A (cultivated land), B (wasteland), and C (lake boundary), were selected to reveal groundwater dynamics ( Figure 4). Figure 4 shows the dynamic process of groundwater depth during the growth period. The dynamics of the groundwater table were greatly affected by irrigation and rainfall.

Groundwater Data
Seventeen groundwater observation wells were established in the study area, including seven key observation wells (in red circle) and ten regularly observed wells (in blue circle). Automatic groundwater sensors (CTD-10, Meter Company, San Francisco, CA, USA) were installed in key observation wells, and an EM50 logger (Meter Company) was used to record groundwater table, temperature, and electrical conductivity (EC) data in one-hour intervals. The groundwater depth of the regularly observed wells was measured every 7 days. Three key observation wells, A (cultivated land), B (wasteland), and C (lake boundary), were selected to reveal groundwater dynamics ( Figure 4). Figure 4 shows the dynamic process of groundwater depth during the growth period. The dynamics of the groundwater table were greatly affected by irrigation and rainfall. growth period (1 June to 1 October) was 109.6 mm and 64.88 mm (Figure 3). The average annual daily evapotranspiration was 2.7 mm/d. In the crop growth period, the average daily evapotranspiration was from 4.3 mm/d to 4.6 mm/d.

Groundwater Data
Seventeen groundwater observation wells were established in the study area, including seven key observation wells (in red circle) and ten regularly observed wells (in blue circle). Automatic groundwater sensors (CTD-10, Meter Company, San Francisco, CA, USA) were installed in key observation wells, and an EM50 logger (Meter Company) was used to record groundwater table, temperature, and electrical conductivity (EC) data in one-hour intervals. The groundwater depth of the regularly observed wells was measured every 7 days. Three key observation wells, A (cultivated land), B (wasteland), and C (lake boundary), were selected to reveal groundwater dynamics ( Figure 4). Figure 4 shows the dynamic process of groundwater depth during the growth period. The dynamics of the groundwater table were greatly affected by irrigation and rainfall.

Soil Data
There were 63 soil observation points in the study area. The distance between each soil observation point was 50 m. Soil samples in each soil zone (0 to 10, 10 to 20, 20 to 40, 40 to 60, 60 to 80, and 80 to 100 cm) were collected in 10-day intervals over the period of crop growing. The oven-drying method was used to measure soil moisture, and the EC values of the soil extract solution with a soil-water ratio of 1:5 was determined using a conductivity meter (DDS-307A type, Shanghai Youke Instrument Company, Shanghai, China). In addition, there were seven key soil observation points. The soil samples at the key soil observation points were collected in each soil zone (0-20, 20-40, 40-60, 60-80, 80-100, 100-120, 120-160, 60-200, 200-250, and 250-300 cm). Soil automatic sensors (5TE, Meter Company) were installed at the key observation points to monitor soil moisture, EC, and temperature, which were recorded by the EM50 logger.
The conversion formula of the EC measured by the 5TE in the EC values of the soil extract solution with a soil-water ratio of 1:5 is [28]: where SSC is soil salt (g kg −1 ). Salt in the soil solution (CSW) can be calculated as [30]: where θ is the moisture of the soil (cm 3 cm −3 ), γ indicates dry soil bulk density (g cm −3 ), and Csw is salt in the soil solution (g L −1 ).
Three key soil points, A (cultivated land), B (wasteland), and C (lake boundary), were

Soil Data
There were 63 soil observation points in the study area. The distance between each soil observation point was 50 m. Soil samples in each soil zone (0 to 10, 10 to 20, 20 to 40, 40 to 60, 60 to 80, and 80 to 100 cm) were collected in 10-day intervals over the period of crop growing. The oven-drying method was used to measure soil moisture, and the EC values of the soil extract solution with a soil-water ratio of 1:5 was determined using a conductivity meter (DDS-307A type, Shanghai Youke Instrument Company, Shanghai, China). In addition, there were seven key soil observation points. The soil samples at the key soil observation points were collected in each soil zone (0-20, 20-40, 40-60, 60-80, 80-100, 100-120, 120-160, 60-200, 200-250, and 250-300 cm). Soil automatic sensors (5TE, Meter Company) were installed at the key observation points to monitor soil moisture, EC, and temperature, which were recorded by the EM50 logger.
The conversion formula of the EC measured by the 5TE in the EC values of the soil extract solution with a soil-water ratio of 1:5 is [28]: where EC 1:5 represents the EC values of the soil extract solution with a soil-water ratio of 1:5 (dS m −1 ) and EC 5TE is the EC values measured by 5TE (dS m −1 ).
The EC e of all the samples was estimated using Equation (2) [29]: The conversion formula of EC 1:5 in soil salt (SSC) is [29]: where SSC is soil salt (g kg −1 ). Salt in the soil solution (C SW ) can be calculated as [30]: where θ is the moisture of the soil (cm 3 cm −3 ), γ indicates dry soil bulk density (g cm −3 ), and C sw is salt in the soil solution (g L −1 ). Three key soil points, A (cultivated land), B (wasteland), and C (lake boundary), were selected to unveil soil textural properties. The soil particle sizes were measured utilizing a dry particle size analyzer (Helos & Rodos, Germany New Partec, Dresden, Germany). The saturated hydraulic conductivity, bulk density, and saturated moisture of undisturbed soil samples were measured (100 m 3 ). The results showed that the soil texture in the 0-300 cm soil zone of cultivated land was sandy loam. The soil profile of the wasteland and lake boundary could be described by two layers. The soil texture in the 0-80 cm soil zone of the wasteland was sandy loam, and it was sand in the 80-300 cm soil zone. The soil texture in the 0-20 cm soil zone of the lake boundary was sandy loam, and in the 20-300 cm soil zone, it was sand. Table 1 summarizes the important physical properties of soil at a depth of 0-300 cm in the cultivated land (A), wasteland (B), and lake boundary (C). Table 2 presents the irrigation schedule for 2018 and 2019. The scheduling of irrigation depends on both the experience of farmers and the conditions of canal water. A trapezoidal thin-wall weir was used to measure the irrigation volumes. Three irrigation events were conducted per year. Rainfall and irrigation water were collected on rainy days and during the irrigation period. The lake and groundwater samples were collected every 10 days. Each water sample was collected in triplicate. The EC (Electrical Conductivity) was measured using a conductivity meter (DDS-307A). The conversion formula of the EC and the total dissolved solids (TDS) is [31]: The irrigation schedule of the study area is shown in Table 2. The irrigation quota was 386 mm and 380 mm in 2018 and 2019. The TDS of irrigation water varied between 0.54 g/L and 0.65 g/L.

Crop Growth Period
The sunflower development stages were recorded ( Table 3). The planting time in 2018 and 2019 was 31 May and 2 June, respectively. It took 113 days and 111 days from sowing to maturity in 2018 and 2019.

Soil Salt and Water Dynamic Model
Hydrus-1D is a vertical, one-dimensional model that is used to simulate soil water flow, heat and solute transport, and root water uptake in a saturated-unsaturated soil profile [32]. With the improvement of the Hydrus-1D model, it has been widely applied to simulate carbon dioxide and major ion solution movement.

The Equation for Soil Water Flow
The governing equation for soil water flow is a modified form of the Richards equation: where θ is the water volumetric content (cm 3 cm −3 ), z represents the vertical space coordinate (cm) (positive upward), t is time (d), h is the water pressure head (cm), K is the function of hydraulic conductivity (cm d −1 ), and S represents a sink term (cm 3 cm −3 d −1 ) representing the rate of water uptake by roots. The hydraulic properties of soil are represented by the van Genuchten and Mualem [33,34] model: where S e represents effective saturation (−), θ r is the residual moisture (cm 3 cm −3 ), θ s is the saturated moisture (cm 3 cm −3 ), K s is hydraulic conductivity under a saturated condition (cm day −1 ), α (cm −1 ), n (−) and m (−) are statistical parameters, and l is the pore connectivity parameter (−). The impacts of osmotic and water stress on the uptake of water by roots are taken to be multiplicative: where T p is the plausible rate of the uptake of water (cm d −1 ), h φ represents the osmotic head (cm), b(z) is the distribution of the normalized uptake of water in the root zone (cm −1 ), and α rs (h φ ) and α rw (h) are osmotic response and water stress response functions, respectively. The present study selected the function form of Feddes to estimate α rw (h) [35] and the modified S-shape function [

The Equation for Soil Solute Transport
Solute transport was described using the convective-dispersive equation: where ρ b is bulk soil density (g cm −3 ), c is the solute concentration of the liquid phase (mg cm −3 or g L −1 ), s is the concentration of the adsorbent (g g −1 ), q represents the water vertical flux (cm d −1 ), D dif represents the coefficient of diffusion (cm 2 d −1 ), and D dis represents the coefficient of dispersion (cm 2 d −1 ).

Evaporation and Transpiration
The potential evapotranspiration of crops (ET p ) was estimated as ET p = K c × ET 0 , where K c is the sunflower crop coefficient. The K c values of sunflower were 0.35, 1.15, and 0.35 for the late, middle, and early growth stages, respectively, according to FAO56 recommendations [27]. In the Hydrus-1D model, ET P is categorized into potential transpiration (T P ) and potential soil evaporation (E P ) according to the leaf area index (LAI) [39,40]. where K represents the coefficient of the extinction of the crop canopy-K was set to 0.83 in this study-and LAI represents the leaf area index.

Profile Discretization
Three key observation points, A (cultivated land), B (wasteland), and C (lake boundary), were selected to simulate the salt and water transport of soil profiles. The soil profiles of the cultivated land, wasteland, and lake boundary were each set as 300 cm below the surface. The soil profiles of the cultivated land, wasteland, and lake boundary were divided into one horizon, two horizons, and two horizons, respectively, as per the texture of soil shown in Table 1. The 1D, vertical profile of the soil was subdivided into 301 nodes with 1 cm uniform spacing. The period of 1 June to 30 September was chosen as the period of simulation. The initial timestep, minimum timestep, and maximum timestep were 0.1 d, 0.001 d, and 5 d, respectively. The simulation timesteps and simulation domain were the same for the cultivated land, wasteland, and lake boundary.

Initial Condition and Boundary Condition
The atmosphere boundary condition containing the surface layer was set to the upper water flow boundary condition, and the variable pressure head was assigned as the bottom boundary condition. The atmospheric boundary conditions need the daily rainfall and irrigation flux (cultivated land), potential evaporation, and transpiration data. The groundwater depth (GWD) needed to be measured, and it was then input into the time-variable settings of the pressure head of the bottom boundary condition.
The upper and bottom boundary conditions were set as the concentration flux boundary. The upper boundary condition of the cultivated land was the irrigation water concentration, and for the wasteland and lake boundary, it was the rainfall concentration (the default value of 0). The bottom boundary condition for the cultivated land, wasteland, and lake boundary was the groundwater concentration. The soil salt and water initial conditions were established according to field observation values.
A classical logistic growth function was used to model root growth. The present study assumed that a 50% rooting depth would be reached by halfway through the growing season [32]. According to field observations, the maximum depth of the roots of sunflower was set at 60 cm. The piecewise linear water stress function proposed by Feddes et al. (1978) [35] and the S-shaped salt stress function proposed by van Genuchten et al. (1984) [36] were used to analyze the effect of soil salt and water dynamics on crop root water uptake.

Model Parameters
The initial soil hydraulic parameters (θ s , θ r , α, n, K s , and l) were determined according to the measured dry soil bulk density and the proportions of silt, clay, and sand, as well as soil water-retention data using Rosetta pedotransfer functions. The parameter for the transport of solutes (dispersivity length, L) was obtained according to prior studies [5,30,41].

Calibration and Verification Procedures
The measured soil water and EC values in different layers of soil (i.e., 0 to 10, 10 to 20, 20 to 40, 40 to 60, 60 to 80, and 80 to 100 cm) relative to 2018 and 2019 were utilized to calibrate and verify the model, respectively. The regression coefficient (b), determination coefficient (R 2 ), mean relative error (MRE), and root mean square error (RMSE) were utilized to assess the performance of the model [30]:

Research Process
Based on the collected soil, groundwater, meteorological, crop growth, and irrigation data, this research studied the water and salt transport of the cultivated land, wasteland, and lake boundary with the Hydrus-1D model. The study flowchart was shown in Figure 5. where

Research Process
Based on the collected soil, groundwater, meteorological, crop growth, and irrigation data, this research studied the water and salt transport of the cultivated land, wasteland, and lake boundary with the Hydrus-1D model. The study flowchart was shown in Figure  5.

Model Calibration
The model was calibrated using the experimental data from 2018. The simulated and measured values of the soil water content and EC in the different soil layers of the cultivated land, wasteland, and lake boundary are shown in Figures 6-8. The simulated values of the soil water content were consistent with the measured values. The MRE values of the cultivated land, wasteland, and lake boundary varied from −0.43% to 4.89%, the RMSE values varied from 0.01 to 0.06 cm 3 cm −3 , the R 2 values varied from 0.75 to 0.90, and the b values were about 1 ( Table 4). The salt transport processes were more complex than that of water; therefore, the simulated salt transport discrepancy was larger, especially in the 0-40 cm soil zone. Owing to the effects of rainfall, irrigation, and evaporation, the convective-dispersive process of salt transport was more complex and difficult to capture. The simulated salt transport discrepancies of the wasteland and the lake boundary, which were only affected by rainfall and evaporation, were smaller than those of the cultivated land. The simulated soil salt dynamics were also well-matched with observations. The MRE values of the cultivated land, wasteland, and lake boundary varied from −2.23% to −0.91%, the RMSE values were about 0.09 dS/m, the R 2 values varied from 0.89 to 0.91, and the b values varied from 0.95 to 0.97 ( Table 4). The performance of the model in representing the soil salt and water transport processes was acceptable. Table 5 summarizes the values of the calibrated parameters representing soil hydraulic and solute transport processes, and Table 6 shows the root water uptake parameters.

Model Verification
The model was verified against experimental data from 2019. All model statistical coefficients and input parameters were carried over from the calibrated model. The verification results are shown in Figure 9.

Quantitative Estimation of the Moisture of Soil in a Typical Period
Since the climate did not significantly change from 2018 to 2019, the results for 2018 will be discussed below. The study area was irrigated on 23 June and 5 July. This study selected 22 June, 28 June, 4 July, 10 July, and 30 September as the typical dates to study characteristics of the salt and water transport in the 1 m soil zone of the cultivated land, wasteland, and lake boundary. The dynamics of the soil water are shown in Figure 10.
For the cultivated land, before and after the first irrigation (22 June-28 June), the soil water in the 0-20 cm soil zone changed greatly; i.e., it increased by 41%, and the water of the 20-60 cm soil zone increased by 18-20%, while changes in water in the 60-100 cm soil zone were not obvious. The first irrigation was 104 mm, and the water storage in the 1 m soil zone increased by 81 mm on the fifth day after irrigation; therefore, the utilization efficiency of irrigation was 78%. On the 11th day after irrigation (11 July), the soil water storage decreased by 66 mm compared with that on 28 June, when water in the 0-20 cm Lake boundary Lake boundary Figure 9. A comparison between observed and model-simulated soil water content (left) and soil electrical conductivity (right) for model verification for different land use types.

Quantitative Estimation of the Moisture of Soil in a Typical Period
Since the climate did not significantly change from 2018 to 2019, the results for 2018 will be discussed below. The study area was irrigated on 23 June and 5 July. This study selected 22 June, 28 June, 4 July, 10 July, and 30 September as the typical dates to study characteristics of the salt and water transport in the 1 m soil zone of the cultivated land, wasteland, and lake boundary. The dynamics of the soil water are shown in Figure 10. soil layer declined by 27%; for the 20-60 cm layer, it declined by 14-16%, and the change in the 80-100 cm layer was not obvious. Before and after the second irrigation (4 July-10 July), the water storage in the 1 m soil zone increased by 69 mm because the water in the second irrigation was 96 mm, and the utilization coefficient of the irrigation water was 72%, which was 6% less than that of the first irrigation. The interval between the first and the second irrigation was short. After the second irrigation, the changes in the soil water storage were smaller. The water in the 0-20 cm soil zone increased by 27%; in the 20-60 cm soil zone, it increased by 14-16%; and the change was not visible in the 80-100 cm layer. On the 87th day (30 September) after irrigation, soil water storage decreased by 186 mm. Soil water in the 0-20 cm soil zone declined by 56%; in the 20-40 cm soil zone, it declined by 57%; in the 40-60 cm soil zone, it declined by 34%; and in the 60-100 cm soil zone, it declined by 13-18%.  For the cultivated land, before and after the first irrigation (22 June-28 June), the soil water in the 0-20 cm soil zone changed greatly; i.e., it increased by 41%, and the water of the 20-60 cm soil zone increased by 18-20%, while changes in water in the 60-100 cm soil zone were not obvious. The first irrigation was 104 mm, and the water storage in the 1 m soil zone increased by 81 mm on the fifth day after irrigation; therefore, the utilization efficiency of irrigation was 78%. On the 11th day after irrigation (11 July), the soil water storage decreased by 66 mm compared with that on 28 June, when water in the 0-20 cm soil layer declined by 27%; for the 20-60 cm layer, it declined by 14-16%, and the change in the 80-100 cm layer was not obvious. Before and after the second irrigation (4 July-10 July), the water storage in the 1 m soil zone increased by 69 mm because the water in the second irrigation was 96 mm, and the utilization coefficient of the irrigation water was 72%, which was 6% less than that of the first irrigation. The interval between the first and the second irrigation was short. After the second irrigation, the changes in the soil water storage were smaller. The water in the 0-20 cm soil zone increased by 27%; in the 20-60 cm soil zone, it increased by 14-16%; and the change was not visible in the 80-100 cm layer. On the 87th day (30 September) after irrigation, soil water storage decreased by 186 mm. Soil water in the 0-20 cm soil zone declined by 56%; in the 20-40 cm soil zone, it declined by 57%; in the 40-60 cm soil zone, it declined by 34%; and in the 60-100 cm soil zone, it declined by 13-18%.
For the wasteland, the storage of water in the 1 m layer of soil increased by 36 mm (22 June-28 June), while soil water in the 60-80 cm soil zone increased by 15%; in the 80-100 cm layer, it increased by 12%, while soil water changes in the 0-40 cm layer were not obvious. After irrigation, the groundwater level of the wasteland rose, which caused the water in the deep soil zone to increase. Water storage in the 1 m soil zone increased by 51 mm from 4 July to 10 July, which was 15 mm higher than that of the first irrigation. After the second irrigation, the groundwater level of the wasteland rose more than that of the first irrigation ( Figure 4). Water in the 40-100 cm layer increased by 12-14%, and in the upper soil zone, it showed a slightly increasing trend. Soil water storage decreased by 150 mm from 10 July to 30 September, when water in the 0-20, 20-60, and 60-100 cm soil zones decreased by 45%, 20%, and 37%, respectively.
For the lake boundary, the storage of water in the 1 m soil zone increased by 27 mm (22 June-28 June). Owing to shallow groundwater depth, water in the 60-100 cm layer was close to saturation, and water in the 0-60 cm soil zone increased by 6-9%. Water storage in the 1 m soil zone increased by 9 mm (4 July-10 July). Since the interval between the first and the second irrigation was short, soil water in the deep soil zone was saturated, and water in the 0-40 cm soil zone increased by 3-6%. The water storage in the 1 m soil zone decreased by 60 mm (10 July-30 September).

Quantitative Estimation of Soil Salt in a Typical Period
The dynamics of soil EC are shown in Figure 11. For cultivated land, soil salt in the 0-20 cm soil zone decreased by 26% (22 June-28 June), while in the 80-100 cm soil zone, it increased by 24%. Soil salt in the 0-20 cm soil zone decreased by 20% (4 July-11 July), while in the 80-100 cm soil zone, it increased by 9%. Since the interval between the first and the second irrigation was only 12 days, the amount of accumulated salt in the soil profile was very small. The leached salt from the upper soil zone after the second irrigation was 6% less than that after the first irrigation, while the salt in the deep soil zone increased by 15%. The salt in the 0-20 cm soil zone increased by 49% (10 July-30 September). The salt in the 0-20 cm soil zone increased by 13% (22 June-30 September), which shows a slightly increasing trend, and in the 20-40 cm and 40-60 cm soil zones, it increased by 5% and 26%, respectively.
For the wasteland, soil salt in the 0-10 cm soil zone on 22 June-June 28 and 4 July-10 July increased by 20% and 6%, respectively. After the first irrigation, the salt accumulation in the upper soil zone was 14% larger than that after the second irrigation. After the first irrigation, a large amount of groundwater salt in the cultivated land was transported to the wasteland groundwater, which caused the salt in the deep soil zone to increase (80 cm, 1.09-1.20 dS/m; 100 cm, 1.09-1.20 dS m −1 ). The salt in the upper soil zone increased by 15%, and in the deep soil zone, it increased by 11% from 10 July to 30 September. The salt in the upper soil zone increased by 37%, and in the deep soil zone, it increased by 15% from 22 June to 30 September. Owing to strong evaporation, the accumulated salt content in the upper soil zone was high.

Soil Salt and Water Balance Analysis
As shown in Table 7, there were significant differences in the water-salt balance in the 0-100 cm soil zone of the cultivated land (A), wasteland (B), and lake boundary (C). During the growth period in 2018 and 2019, the irrigation water was 200 mm and 210 mm, respectively. The cultivated land was mainly affected by irrigation, while the wasteland and the lake boundary were not. Therefore, the deep percolation and capillary rise in the cultivated land, wasteland, and lake boundary in the study area were quite different. During the growth period in 2018 and 2019, the deep percolation of the cultivated land was 105 mm and 110 mm, which was 34% and 40% of the applied water (rainfall and irrigation), respectively. The capillary rise in the cultivated land was 118 mm and 140 mm, which was 24% of the evapotranspiration. The capillary rises in the wasteland was 29% and 35% of their own evaporation, and those of the lake boundary were 62% and 68% of their own evaporation in 2018-2019, respectively. The salt storage in the 1 m soil zone of the cultivated land, wasteland, and lake boundary had an increasing trend, which was an average increase of 19%, 27%, and 37%, respectively. Under the effects of irrigation and rainfall on cultivated land, the amount of salt for percolation and capillary rise varied greatly from 2012 to 1568 g m −2 and from 2291 to 1904 g m −2 from 2018 to 2019, respectively. The capillary rise salt in the wasteland and the lake boundary was an average of 2519 and 7331 g m −2 from 2018 to 2019, respectively. Owing to the shallow groundwater depth and the high groundwater TDS in the lake boundary, the capillary rise salt in the lake bound- For the lake boundary, the salt in the upper soil zone increased by 48%, and in the deep soil zone, it increased by 13%. The amount of accumulated salt in the upper soil zone was 35% greater than in the deep soil zone from 22 June to 30 September. During the irrigation period, the groundwater salt in the wasteland was transported to the lakes, which caused salt accumulation in the deep soil zone (80-100 cm). Those factors, including shallow groundwater depth, sandy loam in the 0-20 cm soil zone, and sand in the 20-100 cm soil zone, are beneficial to salt accumulated in the upper soil zone.

Soil Salt and Water Balance Analysis
As shown in Table 7, there were significant differences in the water-salt balance in the 0-100 cm soil zone of the cultivated land (A), wasteland (B), and lake boundary (C). During the growth period in 2018 and 2019, the irrigation water was 200 mm and 210 mm, respectively. The cultivated land was mainly affected by irrigation, while the wasteland and the lake boundary were not. Therefore, the deep percolation and capillary rise in the cultivated land, wasteland, and lake boundary in the study area were quite differ-ent. During the growth period in 2018 and 2019, the deep percolation of the cultivated land was 105 mm and 110 mm, which was 34% and 40% of the applied water (rainfall and irrigation), respectively. The capillary rise in the cultivated land was 118 mm and 140 mm, which was 24% of the evapotranspiration. The capillary rises in the wasteland was 29% and 35% of their own evaporation, and those of the lake boundary were 62% and 68% of their own evaporation in 2018-2019, respectively. The salt storage in the 1 m soil zone of the cultivated land, wasteland, and lake boundary had an increasing trend, which was an average increase of 19%, 27%, and 37%, respectively. Under the effects of irrigation and rainfall on cultivated land, the amount of salt for percolation and capillary rise varied greatly from 2012 to 1568 g m −2 and from 2291 to 1904 g m −2 from 2018 to 2019, respectively. The capillary rise salt in the wasteland and the lake boundary was an average of 2519 and 7331 g m −2 from 2018 to 2019, respectively. Owing to the shallow groundwater depth and the high groundwater TDS in the lake boundary, the capillary rise salt in the lake boundary was three times that of the wasteland. In recent years, with the implementation of water-conservation projects, the salt and water balance system that had been formed over many years has been disrupted. On the basis of the Hydrus-1D model, a soil water-salt transport analysis and a balance analysis in a cultivated land-wasteland-lake system were performed.
This study found that the groundwater table and soil water content in the deep soil zone of the wasteland and lake boundary rose during the irrigation period, and soil water content in the 60-100 cm soil zone of the wasteland increased by 12-15%. Li et al. (2010) showed that groundwater lateral flow between the cultivated land and wasteland was engendered after irrigation, and they also found that soil water content in the 20-70 cm soil zone of the wasteland increased by 4-11% [17]. Our result was larger than that of Li because it is difficult to increase water in the 20-40 cm soil zone with capillary action. This study found that the utilization efficiency of irrigation was 78%. Wang et al. (2020), based on the δD and δ 18 O isotopes, revealed that the utilization efficiency of irrigation was 82% [41]. The results showed there was a sharp decline in soil EC with irrigation, after which it gradually increased prior to the subsequent irrigation. Hao et al. (2015) showed that irrigation water diluted the soil solution, and the salt was leached to the deep soil zone [42]. He et al. (2017) used the Hydrus-1D model to study the irrigated cropping system, and they showed that freshwater irrigation can leach out salts from the 100 cm root-profile depth [43]. Tong et al. (2018) found that the mean soil salt reduction rate under drip irrigation and sub-surface pipe drainage was 19.30% and 58.12%, respectively [44]. This study found that soil salt in the 0-20 cm soil zone of cultivated land decreased by 20~26% during the irrigation period. The salt reduction rate of the border irrigation was higher than that of drip irrigation but lower than that of the sub-surface pipe drainage. This study found that the salt storage in the 1 m soil zone of the cultivated land, wasteland, and lake boundary had an increasing trend, which was an average increase of 19%, 27%, and 37%, respectively. Wang et al. (2021) studied the soil water and salt transport process in a sand dune-wasteland-lake system, and they found that the salt accumulation rate of the wasteland in the 1 m soil zone was 25% during the growth period [45]. Their study results are similar to ours, which indicates that the wasteland was in a state of salt accumulation during the growth period. During the growth period, the lake boundary was in a state of salt accumulation. Those factors, including shallow groundwater depth, sandy loam in the 0-20 cm soil zone, and sand in the 20-100 cm soil zone, are beneficial to the salt accumulated in the upper soil zone; conversely, those factors are not good for the upward movement of salt and water due to weak capillary action. Chen et al. (2010) found that it is significant and necessary for salinity control to conduct large flood irrigation to decrease salt accumulation after harvest [46].

Soil Water and Salt Balance Analysis
We revealed that the deep percolation of cultivated land was 34-40% of the applied water (rainfall and irrigation). Ren et al. (2018) pointed out that about 36% of the total applied water was stored in shallow groundwater [31], which indicated that the results of this study are acceptable. Ren et al. (2016) found that capillary rise in watermelon, sunflower, and maize fields accounted for over 45%, 24%, and 33% of their respective total evapotranspiration [5]. Xu et al. (2015) found that cumulative capillary rise could reach 40% of evapotranspiration when the groundwater depth was set to 100 cm [47]. We found that capillary rise in the cultivated land, wasteland, and lake boundary was 24%, 29-35%, and 62-68% of their own evapotranspiration, respectively. These results are similar to prior studies, which were credible. The capillary rise in the lake boundary was larger than that of the wasteland since the groundwater table of the lake boundary was 30 cm~140 cm during the growth period (Figures 3 and 4).  pointed out that 40% of the total salt was retained in a wasteland, and 39% was retained in farmland [48]. Li et al. (2010) showed that the wasteland was in a state of salt accumulation, and salt accumulation in the upper soil zone was the most serious [17]. Based on the previous research, we further found that the salt accumulation in the upper soil zone of the cultivated land, wasteland, and lake boundary was 13%, 37%, and 48%, respectively, and in the deep soil zone, it was 34%, 15%, and 13%, respectively, which is significant in preventing soil salinization in the Hetao Irrigation District. This study suggests some measures should be taken to reduce the salt content in the upper soil zone of the wasteland and lake boundary and in the deep soil zone of the cultivated land to prevent soil salinization.

Conclusions
The Hydrus-1D model was successfully used to evaluate the soil salt and water transport characteristics of the cultivated land, wasteland, and lake boundary in the HID. The simulations of the model provide a good representation of measurements. The conclusions are as follows: (1) On the fifth day after irrigation, the soil water storage increments of the cultivated land, wasteland, and lake boundary were 69~81 mm, 36~51 mm, and 9~22 mm, respectively; in particular, the soil water content in the 60-100 cm soil zone of the wasteland increased by 12-15%. On the 87th day after irrigation, the water storage in the cultivated land, wasteland, and lake boundary decreased by 186 mm, 150 mm, and 60 mm; the soil water content in the 0-60 cm soil zone of the cultivated land and wasteland decreased by 34-56% and 20-45%.
(2) The deep percolation of cultivated land was 34-40% of the applied water (rainfall and irrigation). The capillary rise in the cultivated land, wasteland, and lake boundary was 24%, 29-35%, and 62-68% of their own evapotranspiration, respectively.
(3) In the whole growth period, the salt storage in the 1 m soil zone of the cultivated land, wasteland, and lake boundary, on average, increased by 19%, 27%, and 37%, respec-tively. Salt accumulation in the upper soil zone of the cultivated land, wasteland, and lake boundary was 13%, 37%, and 48%. Soil salinization in the upper soil zone of the wasteland and lake boundary was serious, so some measures should be taken to reduce the salt content to prevent soil salinization.
(4) Since the study scale of the cultivated land, wasteland, and lake systems was small, the simulation accuracy of the regional hydrological model is limited. It is necessary to expand the study area and establish a large hydrological monitoring network to study the hydrological process in the future.

Informed Consent Statement:
Written informed consent has been obtained from the patient(s) to publish this paper.