Response of Liquid Water and Vapor Flow to Rainfall Events in Sandy Soil of Arid and Semi-Arid Regions

: In arid and semi-arid regions, rainfall takes on a critical signiﬁcance to both agricultural and engineering construction activities, and the transport process and driving mechanism of soil water under rainfall conditions are in need of further investigation. To clarify the variations in soil moisture, temperature, and liquid and vapor ﬂux under various rainfall scenarios, the Mu Us Sandy Land was selected as the study region, and a water–vapor–heat transport model was established using the Hydrus-1D software with in situ observed soil and meteorological data. The simulated results were in good agreement with the measured data during both the calibration and validation periods, suggesting that the model was accurate and applicable to the study region. The variations in the selected dry and rainy periods proved the signiﬁcant effect of rainfall events on soil matric potential, temperature, and driving forces. When rainfall occurred, the hydraulic conductivity for liquid water rose by three to ﬁve orders of magnitude, driving the liquid water ﬂow downward. In contrast, the vapor ﬂux played a vital role in soil water movement, accounting for about 15% of the total water ﬂux in the shallow layer when the soil was dry, while it became non-signiﬁcant during rainy periods due to the reduction in hydraulic conductivity for vapor and the temperature gradient. These results clariﬁed the mechanisms of soil liquid water and vapor movement in arid areas, which could provide scientiﬁc support for future studies on vegetation restoration and ecosystem sustainability in ecologically fragile areas.


Introduction
Arid and semi-arid regions are extensively distributed around the world, covering over 30% of the Earth's surface, and their proportion continues to expand [1,2].Affected by strong evaporation and scarce rainfall, the soil is relatively dry and vegetation is sparse, leading to gradual land desertification [3][4][5].In these regions, rainfall infiltration has long been regarded as the greatest challenge to the hydrological cycle, which plays a critical role in agricultural activities, as well as affecting soil water distribution [6][7][8].Meanwhile, soil water in these dry regions is a limited resource, which is thought to be the primary regulator for land desertification, and is a paramount water source for the growth and transpiration of many plants, especially in areas with a deep groundwater depth [9,10].Since rainfall-derived soil water is more readily available, studies have proved that plants preferentially utilize soil water over groundwater [11].Therefore, clarifying the mechanism of soil water movement and its response to rainfall events can provide more insights Agronomy 2023, 13, 2424 2 of 21 into the hydrological cycle in arid and semi-arid regions, which is critical to vegetation restoration and ecosystem sustainability in ecologically fragile areas [12][13][14].
As the key zone connecting various strata, the vadose zone provides a place for the transmission of soil water and heat, the exchange of water and energy between the environment and soil, and the root water uptake [15][16][17].A considerable amount of research has been conducted on the vadose zone, especially regarding the mechanism of rainfall infiltration and soil water movement [8,18].For instance, Cheng et al. [19] estimated the soil water response to precipitation based on the water balance method in the Mu Us Sandy Land.They proposed that the wetting front for 53.8 and 24.2 mm of accumulative rainfall could penetrate to 200 cm depth when the antecedent soil water contents were approximately 0.065 and 0.102 cm 3 cm −3 , respectively.In dry sand dunes, it was calculated that the wetting front could soon recharge to 100 cm depth after 100 mm of rainfall and reached 140 cm several days later [20].With the rapid advances in computation technology, numerical simulation has been increasingly used in relevant research.Simulation research using numerical software has become more prevalent because fieldwork is complicated and unpredictable [21,22].The emergence of relevant software, such as SHAW, Hydrus-1D, CoupModel, and STEMMUS, has significantly improved the convenience of research related to soil water movement (e.g., improved dependability and increased accuracy of predictions) [23][24][25].
Based on Richard's equation, Philip and de Vries [26] initially developed the theory of coupled liquid and vapor transfer under isothermal conditions.Cass et al. [27] proposed the enhancement factor of vapor flow and estimated it using different calculation methods.Cahill and Parlange [28] quantitatively estimated the effect of vapor flow, suggesting that vapor flux can account for an appreciable 25% of total water flux at the surface.Meanwhile, vapor flux in the upper layer could transport a significant amount (up to 50%) of the overall heat flux.Since then, vapor flow driven by the temperature gradient has been considered an important component of soil water flow, especially in drylands with low soil moisture [18,[29][30][31][32]. Nevertheless, research on the coupled liquid and vapor flow during various dry and wet conditions is still limited, especially for studies based on in situ experiments rather than laboratory tests [33].The soil water transport process and its driving factors, including liquid water and vapor water, have not been fully explored under drought and various rainfall conditions.
To gain more insights into these issues, the Mu Us Sandy Land in northwest China was selected as the research object.Through field investigation, in situ and laboratory tests, and numerical simulation, this study attempted to clarify the response of both liquid and vapor flow to rainfall events.The specific research goals of this work were (i) to establish a coupled water, vapor, and heat transport model using Hydrus-1D and verify its accuracy and applicability; (ii) to clarify the distribution and variations in soil moisture, temperature, and water flux before and after rainfall; and (iii) to discuss the impact of various rainfall events on soil water movement and the significance of vapor flux in soil water transfer.

Field Tests 2.1.1. Description of the Study Site
The Mu Us Sandy Land, one of the four largest sandy lands in China, covering an area of 42,200 km 2 , is located between the Ordos Plateau and the Loess Plateau of northwestern China, as shown in Figure 1.It is an ecologically fragile zone in the Yellow River basin and one of the focal points of the country's sandstorm defense.Sand is the dominant soil in areas of the Mu Us Sandy Land, and most of the sandy land belongs to the temperate semi-arid region [11,14,25,34].In this study, an in situ observation site was selected in the southeast of the Mu Us Sandy Land (Figure 1, adapted from Chen et al. [14]), located in Yulin Desert Ecological Station of the State Forestry Administration (38°23′ N, 109°42′ E).As indicated by the collected meteorological data since 1961 (from the local weather station), the average yearly temperature is 6.4 °C, and the annual evaporation and precipitation are 2343 and 360 mm, respectively.Note the precipitation distribution is uneven, mainly concentrated from July to September, accounting for more than 70% of the total precipitation.The groundwater depth usually fluctuates between 8 and 10 m around the study site.

In Situ Experiments and Observation
At the study site, the soil surface is flat with vegetation covering less than 5%, meaning that the impact of vegetation on soil water transport can be ignored.To obtain soil physical properties, samples at different depths (with a 20 cm interval from the surface to 300 cm depth and a 50 cm interval from 300 cm to 630 cm depth) were collected using the soil drilling method and then investigated in the laboratory to determine their physical properties.According to the soil classification standard of the United States Agriculture Department, sand is the dominating soil [35].The soil profile was divided into five layers, as listed in Table 1.For determining the soil hydraulic properties, both experimental (including oven-drying and double-ring infiltration tests) and predictive methods were applied.The specific parameter determination processes and adopted values are introduced in the following section.In this study, an in situ observation site was selected in the southeast of the Mu Us Sandy Land (Figure 1, adapted from Chen et al. [14]), located in Yulin Desert Ecological Station of the State Forestry Administration (38 • 23 N, 109 • 42 E).As indicated by the collected meteorological data since 1961 (from the local weather station), the average yearly temperature is 6.4 • C, and the annual evaporation and precipitation are 2343 and 360 mm, respectively.Note the precipitation distribution is uneven, mainly concentrated from July to September, accounting for more than 70% of the total precipitation.The groundwater depth usually fluctuates between 8 and 10 m around the study site.

In Situ Experiments and Observation
At the study site, the soil surface is flat with vegetation covering less than 5%, meaning that the impact of vegetation on soil water transport can be ignored.To obtain soil physical properties, samples at different depths (with a 20 cm interval from the surface to 300 cm depth and a 50 cm interval from 300 cm to 630 cm depth) were collected using the soil drilling method and then investigated in the laboratory to determine their physical properties.According to the soil classification standard of the United States Agriculture Department, sand is the dominating soil [35].The soil profile was divided into five layers, as listed in Table 1.For determining the soil hydraulic properties, both experimental (including oven-drying and double-ring infiltration tests) and predictive methods were applied.The specific parameter determination processes and adopted values are introduced in the following section.
Subsequently, an observation well (with a diameter of 150 cm) was dug in which to place measuring instruments (as shown in Figure 1).The Hydra Probe II sensors were installed horizontally at 10, 50, 100, 200, 400, and 630 cm depths through the wall of the monitoring well to observe changes in soil temperature and soil moisture.To avoid the influence of the well itself on the observation data, the distance of the sensors from the wellbore exceeded 50 cm, and the observation started 6 months after the installation to allow the soil to establish better contact with the sensors.Before installation, these sensors were calibrated using gravimetric measurements (using the oven-drying method) taken from the same soil horizons.For soil moisture measurements, the principle of the Hydra Probe II sensors is based on the frequency-domain (FD) method since the permittivity of liquid water (≈78) is significantly higher than the other components (e.g., air, soil matrix).Due to their relatively high accuracy, FD sensors are widely applied for field soil water content observation [36,37].In addition, a micro-meteorology station was also installed at the study site to monitor meteorological variables (i.e., rainfall, air temperature, relative humidity, and wind speed at 2 m above the surface).After installation, the in situ measuring sensors for soil and meteorology were linked to the solar-powered automatic data loggers CR1000 and Em50, respectively, and the observation data were recorded at a 10 min interval.In this study, the period of 1 April-31 October 2018, was selected for detailed analysis.The reasons for this selection are mainly that approximately 90% of the total rainfall is concentrated from April to October and vegetation growth occurs during this period of the year.Figure 2 presents the variations in the observed meteorological variables during the analyzed period.It can be seen that the precipitation at the test site reached 640.2 mm, far exceeding the local average annual precipitation.The average air temperature, relative humidity, and wind speed during the analyzed period were 17.1 • C, 62.6%, and 1.5 m s −1 , respectively.Subsequently, an observation well (with a diameter of 150 cm) was dug in which to place measuring instruments (as shown in Figure 1).The Hydra Probe II sensors were installed horizontally at 10, 50, 100, 200, 400, and 630 cm depths through the wall of the monitoring well to observe changes in soil temperature and soil moisture.To avoid the influence of the well itself on the observation data, the distance of the sensors from the wellbore exceeded 50 cm, and the observation started 6 months after the installation to allow the soil to establish better contact with the sensors.Before installation, these sensors were calibrated using gravimetric measurements (using the oven-drying method) taken from the same soil horizons.For soil moisture measurements, the principle of the Hydra Probe II sensors is based on the frequency-domain (FD) method since the permittivity of liquid water (≈78) is significantly higher than the other components (e.g., air, soil matrix).Due to their relatively high accuracy, FD sensors are widely applied for field soil water content observation [36,37].In addition, a micro-meteorology station was also installed at the study site to monitor meteorological variables (i.e., rainfall, air temperature, relative humidity, and wind speed at 2 m above the surface).After installation, the in situ measuring sensors for soil and meteorology were linked to the solar-powered automatic data loggers CR1000 and Em50, respectively, and the observation data were recorded at a 10 min interval.
In this study, the period of 1 April-31 October 2018, was selected for detailed analysis.The reasons for this selection are mainly that approximately 90% of the total rainfall is concentrated from April to October and vegetation growth occurs during this period of the year.Figure 2 presents the variations in the observed meteorological variables during the analyzed period.It can be seen that the precipitation at the test site reached 640.2 mm, far exceeding the local average annual precipitation.The average air temperature, relative humidity, and wind speed during the analyzed period were 17.1 °C, 62.6%, and 1.5 m s −1 , respectively.

Model Establishment 2.2.1. Mathematical Equations
For one-dimensional vertical soil water movement, considering the effect of temperature and vapor transport, the governing equation could be expressed using the revised Richard's equation, as indicated below [38,39]: where θ l and θ v represent the liquid water and vapor content (cm 3 cm −3 ), ρ l and ρ v represent the density of liquid water and water vapor (g cm −3 ), t is time (h), z is the spatial coordinate (cm), h is the pressure head (cm), T is temperature (K), K lh and K vh are the isothermal hydraulic conductivities for liquid water and vapor (cm h −1 ), K lT and K vT are the thermal hydraulic conductivities for liquid water and vapor (cm 2 K −1 h −1 ), and S is the sink term accounting for root water uptake (h −1 ).The conductivities for liquid water and vapor transfer can be expressed as follows: where K s is the saturated hydraulic conductivity (cm h −1 ), S e is the effective liquid saturation, l and m are empirical parameters as defined by Mualem [40], G wt is the gain factor, γ is the surface tension of soil water (J cm −2 ), γ 0 is the surface tension at 25 • C (=71.89 g s −2 ), D is the vapor diffusivity in soil (cm 2 h −1 ), ρ sv is the saturated vapor density (g cm −3 ), H r is the relative humidity, M is the molecular weight of water (=0.018015kg mol −1 ), g is the gravitational acceleration, R is the universal gas constant (=8.315J mol −1 K −1 ), and η is the enhancement factor.The K lh is calculated using the van Genuchten model as follows [41]: where θ r and θ s are the residual and saturated water contents (cm 3 cm −3 ), respectively, α (cm −1 ) and n are empirical parameters.
To determine the soil hydraulic properties at the study site, the oven-drying and the double-ring infiltration methods were utilized to measure saturated water content (θ s ) and the saturated hydraulic conductivity (K s ), while the values of other parameters (including θ r , a, and n) of the van Genuchten model were obtained using the neural network method with the Rosetta module.This module is directly implemented in the Hydrus-1D program, and estimates hydraulic parameters according to the measured data of soil particle composition and bulk density (as listed in Table 1) [35].Subsequently, to obtain simulation results with better fitting accuracy, the Levenberg-Marquardt algorithm, which is also directly implemented into the inverse module of Hydrus-1D, is used to optimize these parameters using the observed soil water content data [42].The estimation process is as follows: where f s is the cumulative error between measured and simulated values, p m (z, t) and p s (z, t) are the measured and simulated soil water content at depth z at time t, respectively, and m q represents the number of data used for the inverse.As indicated by Equation ( 5), this method is based upon the minimization of the discrepancy between the observed values and simulated water content at selected depths.According to a parameter sensitivity analysis, the empirical parameters α and n have a significant impact on the simulation results.Therefore, these two parameters were mainly identified during the optimization process.The optimized parameters for numerical simulation at different soil layers are listed in Table 2. On the other hand, considering the effect of temperature and vapor flow, the governing equation for vertical soil heat flow in the vadose zone could be expressed using the energy conservation equation, as indicated below [38,39]: where C p , C w , and C v are the volumetric heat capacities of moist soil, liquid water and vapor phase (J cm −3 K −1 ), respectively.L 0 is the latent heat of vaporization of liquid water (J cm −3 ), q l and q v are the flux density of liquid water and vapor (cm h −1 ), and λ is the thermal conductivity of soil (J cm −1 h −1 K −1 ), which can be estimated using the Chung and Horton equation [43]: Since the vegetation is quite sparse at the study site, covering less than 5% of the soil surface, the effect of root water uptake for water flow and heat transport is not considered in the simulation process.Detailed calculations for the above-mentioned soil hydraulic and thermal parameters can be found in Zheng et al. [25] and Saito et al. [38].

Establishment of Hydrus-1D Model
Due to its convenient operation and comprehensive functions, research on soil water movement has made substantial use of the Hydrus-1D program [44,45].On the basis of soil water transport research, Scanlon et al. [17] firstly considered vapor flow and its phase-transition effects in Hydrus-1D modeling research.Then, Saito et al. [38] took surface water and thermal equilibrium procedures into account in the Hydrus-1D code, and improved the flexibility of setting meteorological conditions in the model.Since then, the standard Hydrus-1D model could be utilized to investigate the coupled water-vapor-heat transport process.
By collecting the in situ monitoring soil and meteorological data, a model for the combined transfer of liquid water, water vapor, and heat established using the Hydrus-1D software is used to obtain the liquid water and vapor flux.To gain detailed diurnal variations, the unit of the simulation period was hours (h) with a total time of 5136 h and 2184 h for the calibration (1 April-31 October 2018) and validation (1 April-30 June 2019) periods, respectively.Moreover, spatial discretization was carried out at an interval of 2 cm from the surface to 630 cm below, with a total of 316 nodes.For the initial water and heat conditions of the model (Equation ( 8)), the water content and temperature data observed on 1 April at depths of 10, 50, 100, 200, 400, and 630 were used for linear interpolation to obtain the initial conditions of the 0-630 cm soil profile.As indicated in Equations ( 9) and ( 10), the upper boundary for water flow is controlled by meteorological factors and determined by rainfall and evaporation conditions, and free drainage is the lower boundary condition because it is far from the groundwater level (nearly 10 m).For soil heat transport, the calculated surface temperature (using soil temperature data observed in the shallow layer [46]) and measured bottom temperature (from the Hydra Probe II sensor located at 630 cm) values are employed as the circumstances at the upper and lower boundaries, respectively (Equation ( 11)).
where E m is the maximum infiltration or evaporation rate at the current moment (cm h −1 ).
For the atmospheric boundary condition, the variation in rainfall was automatically monitored, while hourly potential evaporation (E p ) was estimated as follows: where ρ vs and ρ va are vapor densities (g cm −3 ) at the soil surface and the 20 cm height (where the air temperature is measured), respectively, and r v and r s are the aerodynamic and soil surface resistances (h cm −1 ) for vapor movement, respectively.These parameters can be estimated using the observed soil water content and meteorological data, as introduced by Bittelli et al. [31] and Kroener et al. [47].The calculated values of E p were in close agreement with the observed data at the Yulin meteorological station.
The main program unit of the Hydrus-1D graphical user interface defines the overall computational environment of the system, containing a project manager and both the preprocessing and post-processing units.The pre-processing unit includes specification of all necessary parameters (as mentioned above) to successfully run the Hydrus-1D FORTRAN codes.The established coupled water-vapor-heat transport model could be utilized to investigate the distribution and variations in soil moisture, temperature, and water flux before and after rainfall.

Model Evaluation
To evaluate the model's accuracy, the measured soil water content and temperature data at depths of 10, 50, 100, 200, 400, and 630 cm during 1 April-31 October 2018 and 1 April-30 June 2019 were used for model calibration and validation, respectively.In addition, mean absolute error (MAE) and root mean square error (RMSE) were the two statistical indicators employed to assess the model's accuracy for describing soil water content and temperature variations.
The observed and simulated soil water contents and temperatures at the research site over the calibration period are compared in Figure 3.The simulated results accurately captured the ups and downs of real soil water content change.Since the middle and bottom soil layers were less affected by external environmental factors, the fitting results of the soil water content were relatively good (e.g., RMSE = 0.005 and 0.002 cm 3 cm −3 at 200 and 400 cm depths, respectively).Considering the major impact of outside environmental influences, such as strong precipitation and evaporation, the fitting results for the shallow soil water content were relatively poor (e.g., MAE = 0.011 cm 3 cm −3 and RMSE = 0.017 cm 3 cm −3 at 10 cm depth), whereas the error was within an acceptable range.On the other hand, it could be found that soil temperature during the analyzed period first increased and then decreased as a result of the impact of changing air temperature.Compared with its soil water content simulation results, the fitted soil temperature matched very well with the measured data, with calculated MAE and RMSE values less than 1 • C for most depths.
captured the ups and downs of real soil water content change.Since the middle and bottom soil layers were less affected by external environmental factors, the fitting results of the soil water content were relatively good (e.g., RMSE = 0.005 and 0.002 cm 3 cm −3 at 200 and 400 cm depths, respectively).Considering the major impact of outside environmental influences, such as strong precipitation and evaporation, the fitting results for the shallow soil water content were relatively poor (e.g., MAE = 0.011 cm 3 cm −3 and RMSE = 0.017 cm 3 cm −3 at 10 cm depth), whereas the error was within an acceptable range.On the other hand, it could be found that soil temperature during the analyzed period first increased and then decreased as a result of the impact of changing air temperature.Compared with its soil water content simulation results, the fitted soil temperature matched very well with the measured data, with calculated MAE and RMSE values less than 1 °C for most depths.In addition, the simulated variations were in good agreement with the observed data during the validation period as well (as shown in Figure 4), with the average RMSE values of 0.017 cm 3 cm −3 and 1 °C for soil water content and temperature, respectively.Taken together, the established model is proved to have high accuracy overall, which suggests that it can be utilized to analyze the coupled liquid water, water vapor, and heat transport inside the arid areas of Mu Us Sandy Land.In addition, the simulated variations were in good agreement with the observed data during the validation period as well (as shown in Figure 4), with the average RMSE values of 0.017 cm 3 cm −3 and 1 • C for soil water content and temperature, respectively.Taken together, the established model is proved to have high accuracy overall, which suggests that it can be utilized to analyze the coupled liquid water, water vapor, and heat transport inside the arid areas of Mu Us Sandy Land.

Selection of Simulation Periods
amounts of rainfall in detail, 5 day short periods with the above-mentioned rainfall scenarios were selected.The selected periods should be as unaffected as possible by previous rainfall, with the soil in a relatively dry state (with the soil water content in the shallow layer less than 0.05 cm 3 cm −3 ) at the beginning of each period.Meanwhile, a drought period (no rainfall for more than 10 days) was selected for comparison.Detailed information for the selected periods is listed in Table 4.

Selection of Simulation Periods
By collecting the measured precipitation data during the calibration period, the statistical information of the specific daily rainfall amount was sorted, as indicated in Table 3.It is evident that the rainfall amount of 25-50 mm occupies the greatest proportion of the total rainfall amount, followed by the amounts of 10-25 mm and less than 10 mm.To investigate the response of soil hydrothermal and water flux to scenarios with different amounts of rainfall in detail, 5 day short periods with the above-mentioned rainfall scenarios were selected.The selected periods should be as unaffected as possible by previous rainfall, with the soil in a relatively dry state (with the soil water content in the shallow layer less than 0.05 cm 3 cm −3 ) at the beginning of each period.Meanwhile, a drought period (no rainfall for more than 10 days) was selected for comparison.Detailed information for the selected periods is listed in Table 4.
In the following sections, the characteristics of the distribution and variation in soil water content, temperature, driving force for water movement, and liquid and vapor flux under non-isothermal conditions on an hourly scale were specifically analyzed.Note that the influence of rainfall was mainly concentrated in the upper section of the unsaturated zone, but not notable for the lower part.Meanwhile, Figure 2 indicates that during the simulation period, meteorological factors mainly affect the soil moisture and temperature at the soil profile of 0-200 cm, with little impact on soil moisture and temperature below Agronomy 2023, 13, 2424 10 of 21 200 cm.Therefore, the following detailed analysis focused on the layer from the surface to 200 cm soil depth.In the selected four periods, the spatial and temporal distributions of the soil matric potential in the profile are depicted in Figure 5. Influenced by evaporation throughout the dry period, the shallow soil's matric potential (mainly for the top 10 cm layer) was obviously lower than that of other soil layers, with a value of −10,000 cm (which was set as the minimum value in the Hydrus-1D model).To clarify the distribution characteristics at various levels of the profile, the minimum value in the figure was set to −1000 cm.As depicted in Figure 5, there were two significant characteristics of the matric potential distribution.In the following sections, the characteristics of the distribution and variation in soil water content, temperature, driving force for water movement, and liquid and vapor flux under non-isothermal conditions on an hourly scale were specifically analyzed.Note that the influence of rainfall was mainly concentrated in the upper section of the unsaturated zone, but not notable for the lower part.Meanwhile, Figure 2 indicates that during the simulation period, meteorological factors mainly affect the soil moisture and temperature at the soil profile of 0-200 cm, with little impact on soil moisture and temperature below 200 cm.Therefore, the following detailed analysis focused on the layer from the surface to 200 cm soil depth.

Soil Matric Potential
In the selected four periods, the spatial and temporal distributions of the soil matric potential in the profile are depicted in Figure 5. Influenced by evaporation throughout the dry period, the shallow soil's matric potential (mainly for the top 10 cm layer) was obviously lower than that of other soil layers, with a value of −10,000 cm (which was set as the minimum value in the Hydrus-1D model).To clarify the distribution characteristics at various levels of the profile, the minimum value in the figure was set to −1000 cm.As depicted in Figure 5, there were two significant characteristics of the matric potential distribution.First, the distribution of the matric potential in the profile during dry and wet days was completely different.In P1, the soil maintained a dry state, and the surface's matric potential decreased to its lowest level.Below the top layers, the matric potential tended to decrease due to the increasing dryness intensity (e.g., decreasing from −1058 cm to −1180 cm at 10 cm).When rainfall occurred, the matric potential distribution in the profile changed significantly.The matric potential of shallow soil increased instantaneously, and the isolines at the depths of 0-50 cm became very dense, with the maximum values reaching −250, −135, and −79 cm, respectively, for P2, P3, and P4.Secondly, different rainfall intensities have different effects on matric potential changes (such as infiltration duration and depth).For instance, in P3, since it only rained on the second and third days, the soil on the surface gradually dried after the rainfall, and the matric potential decreased significantly due to the influence of water infiltration and evaporation.As for P4, the rainfall intensity was relatively high on the second and third days (Table 4), resulting in the matric potential in the profile changing over a wide range.Comparing the results at 100 cm, the matric potential in P4 increased from −370 cm before rainfall to −267 cm at the end of this period, while the matric potential slightly decreased in P2 and P3, reflecting the influence of different rainfall intensities on infiltration depth.First, the distribution of the matric potential in the profile during dry and wet days was completely different.In P1, the soil maintained a dry state, and the surface's matric potential decreased to its lowest level.Below the top layers, the matric potential tended to decrease due to the increasing dryness intensity (e.g., decreasing from −1058 cm to −1180 cm at 10 cm).When rainfall occurred, the matric potential distribution in the profile changed significantly.The matric potential of shallow soil increased instantaneously, and the isolines at the depths of 0-50 cm became very dense, with the maximum values reaching −250, −135, and −79 cm, respectively, for P2, P3, and P4.Secondly, different rainfall intensities have different effects on matric potential changes (such as infiltration duration and depth).For instance, in P3, since it only rained on the second and third days, the soil on the surface gradually dried out after the rainfall, and the matric potential decreased significantly due to the influence of water infiltration and evaporation.As for P4, the rainfall intensity was relatively high on the second and third days (Table 4), resulting in the matric potential in the profile changing over a wide range.Comparing the results at 100 cm, the matric potential in P4 increased from −370 cm before rainfall to −267 cm at the end of this period, while the matric potential slightly decreased in P2 and P3, reflecting the influence of different rainfall intensities on infiltration depth.With increasing depth, the gradient of the matric potential decreases significantly.For example, at the beginning of P1, the matric potential gradients at 10, 20, 50, and 100 cm are 48.3,16.5, 3.9, and 1.2 cm cm −1 , respectively.The decrease in the gradient reflects the reduction in the driving force.There exists a divergent zero-gradient plane around 140~180 cm depth, suggesting that the layer has a relatively high matric potential.Since With increasing depth, the gradient of the matric potential decreases significantly.For example, at the beginning of P1, the matric potential gradients at 10, 20, 50, and 100 cm are 48.3,16.5, 3.9, and 1.2 cm cm −1 , respectively.The decrease in the gradient reflects the reduction in the driving force.There exists a divergent zero-gradient plane around 140~180 cm depth, suggesting that the layer has a relatively high matric potential.Since there was no rainfall in P1, the soil gradually dried, leading to a gradual decline in the zero-gradient plane.Similar to Figure 5, the effect of different rainfall intensities on the matric potential distribution is notable, e.g., the maximum downward gradient value reached −43.8 cm cm −1 in P4.During the rainy period, a convergent zero-gradient plane was also identified, which continuously migrated downward from the surface.When the rainfall stopped, a divergent zero-gradient plane reappeared in the shallow layer due to soil evaporation.there was no rainfall in P1, the soil gradually dried, leading to a gradual decline in the zero-gradient plane.Similar to Figure 5, the effect of different rainfall intensities on the matric potential distribution is notable, e.g., the maximum downward gradient value reached −43.8 cm cm −1 in P4.During the rainy period, a convergent zero-gradient plane was also identified, which continuously migrated downward from the surface.When the rainfall stopped, a divergent zero-gradient plane reappeared in the shallow layer due to soil evaporation.Influenced by the air temperature, there was a significant diurnal temperature change in the soil above the shallow 50 cm soil layer, while the changes in soil temperature were relatively smaller in deep layers.Taking period 1 as an example, the surface temperature changed periodically, ranging from 23.3 °C to 41.7 °C, while the soil temperature at 100 cm increased from 25.2 °C to 25.9 °C with a variation of only 0.7 °C.When it started to rain, the amplitude of the soil temperature change significantly decreased.For instance, in P4, the variation in the surface temperature reached 15.7 °C when it did not rain on the first day, while during the following two rainy days, the surface temperature only showed a gradually decreasing trend.

Soil Matric Potential Gradient
The soil temperature distribution in different periods was also significantly influenced by seasonal changes in air temperature.Deep soil layers' temperatures for P2 were higher than those of shallow soil, which was not consistent with the other three periods.The major reason for the above result was that the vadose zone began to release heat with the decrease in air temperature, causing a gradual decline in temperature in the shallow soil layer.Influenced by the air temperature, there was a significant diurnal temperature change in the soil above the shallow 50 cm soil layer, while the changes in soil temperature were relatively smaller in deep layers.Taking period 1 as an example, the surface temperature changed periodically, ranging from 23.3 • C to 41.7 • C, while the soil temperature at 100 cm increased from 25.2 • C to 25.9 • C with a variation of only 0.7 • C. When it started to rain, the amplitude of the soil temperature change significantly decreased.For instance, in P4, the variation in the surface temperature reached 15.7 • C when it did not rain on the first day, while during the following two rainy days, the surface temperature only showed a gradually decreasing trend.
The soil temperature distribution in different periods was also significantly influenced by seasonal changes in air temperature.Deep soil layers' temperatures for P2 were higher than those of shallow soil, which was not consistent with the other three periods.The major reason for the above result was that the vadose zone began to release heat with the decrease in air temperature, causing a gradual decline in temperature in the shallow soil layer.

Soil Temperature Gradient
The soil temperature gradients across time and space in the profile are presented in Figure 8. Similar to the temperature distribution, the gradient of temperature within the shallow layer also indicated obvious diurnal changes, and the changes were significantly weakened when rainfall occurred.

Soil Temperature Gradient
The soil temperature gradients across time and space in the profile are presented in Figure 8. Similar to the temperature distribution, the gradient of temperature within the shallow layer also indicated obvious diurnal changes, and the changes were significantly weakened when rainfall occurred.There were two main zero-gradient planes in the profile.The first one was the divergent plane, which mainly appeared at night.Since the soil temperature of the shallow layer was lower at night, the gradient beyond this plane was positive.The second one was the convergent plane, which mainly exists in the daytime.On that basis, in response to the increase in surface temperature, there was a negative temperature gradient in the shallow layer.Driven by the temperature gradient, both the flux density and transfer direction of liquid water and vapor flow showed significant diurnal variations, which are elucidated in the following section.With the increase in soil depth, the temperature gradient decreased rapidly.It was usually smaller than 0.05 °C cm −1 below 100 cm depth, resulting in a decrease in the driving force for liquid and vapor flow.

Isothermal Flux
The variations in isothermal liquid (qlh) and vapor (qvh) fluxes driven by the soil water potential gradient during the selected periods are depicted in Figures 9 and 10, respectively, and 2 cm (shallow part), 10 cm (upper part), 50 cm (middle part), and 150 cm (lower part) depths in the profile were employed as representative depths for specific analysis.
It is evident that as the depth increased, the flux drastically dropped, which was largely attributed to the reduction in the soil matric potential and temperature gradients.For qlh, there were obvious daily variations in the shallow layer stratum on dry days, and its values in the upper 100 cm soil layer were all positive, indicating that liquid water There were two main zero-gradient planes in the profile.The first one was the divergent plane, which mainly appeared at night.Since the soil temperature of the shallow layer was lower at night, the gradient beyond this plane was positive.The second one was the convergent plane, which mainly exists in the daytime.On that basis, in response to the increase in surface temperature, there was a negative temperature gradient in the shallow layer.Driven by the temperature gradient, both the flux density and transfer direction of liquid water and vapor flow showed significant diurnal variations, which are elucidated in the following section.With the increase in soil depth, the temperature gradient decreased rapidly.It was usually smaller than 0.05 • C cm −1 below 100 cm depth, resulting in a decrease in the driving force for liquid and vapor flow.

Isothermal Flux
The variations in isothermal liquid (q lh ) and vapor (q vh ) fluxes driven by the soil water potential gradient during the selected periods are depicted in Figures 9 and 10, respectively, and 2 cm (shallow part), 10 cm (upper part), 50 cm (middle part), and 150 cm (lower part) depths in the profile were employed as representative depths for specific analysis.
It is evident that as the depth increased, the flux drastically dropped, which was largely attributed to the reduction in the soil matric potential and temperature gradients.For q lh , there were obvious daily variations in the shallow layer stratum on dry days, and its values in the upper 100 cm soil layer were all positive, indicating that liquid water flowed upward affected by evaporation.When rainfall occurred, the flux became downward, and the flux showed a downward trend density that varied due to different rainfall events.For instance, the maximum values at 2 cm reached −0.087, −0.319 and −0.897 cm h −1 for P2, P3, and P4, respectively, revealing the instantaneous rainfall intensity.
The qvh flux was significantly less than qlh, e.g., usually four and seven orders of magnitude smaller in the profile's top and lower portions, respectively.The maximum value of qvh was only 1.5 × 10 −5 cm h −1 in the whole selected period, exerting a very limited effect on the total water flux.This suggests that the vapor flux driven by the matric potential gradient can be ignored.

Thermal Flux
Under the influence of the temperature gradient, there were variations in thermal liquid (qlT) as well as vapor (qvT) fluxes at different depths, indicated in Figures 11 and 12, respectively.
During the dry period, the qvT flux in the top soil layer was greater than qlT, and both fluxes showed regular daily variations.With the increase in solar radiation in the daytime, the shallow temperature reduced significantly, and the downward flux value was rela- h −1 for P2, P3, and P4, respectively, revealing the instantaneous rainfall intensity.
The qvh flux was significantly less than qlh, e.g., usually four and seven orders of magnitude smaller in the profile's top and lower portions, respectively.The maximum value of qvh was only 1.5 × 10 −5 cm h −1 in the whole selected period, exerting a very limited effect on the total water flux.This suggests that the vapor flux driven by the matric potential gradient can be ignored.

Thermal Flux
Under the influence of the temperature gradient, there were variations in thermal liquid (qlT) as well as vapor (qvT) fluxes at different depths, indicated in Figures 11 and 12, respectively.
During the dry period, the qvT flux in the top soil layer was greater than qlT, and both fluxes showed regular daily variations.With the increase in solar radiation in the daytime, the shallow temperature reduced significantly, and the downward flux value was rela- The q vh flux was significantly less than q lh , e.g., usually four and seven orders of magnitude smaller in the profile's top and lower portions, respectively.The maximum value of q vh was only 1.5 × 10 −5 cm h −1 in the whole selected period, exerting a very limited effect on the total water flux.This suggests that the vapor flux driven by the matric potential gradient can be ignored.

Thermal Flux
Under the influence of the temperature gradient, there were variations in thermal liquid (q lT ) as well as vapor (q vT ) fluxes at different depths, indicated in Figures 11 and 12, respectively.
pletely different.The qvT value in the shallow layer decreased significantly, which was mainly because the increase in water content reduced the hydraulic conductivity of the vapor.At this time, the influence of the vapor flux can be ignored.On the contrary, qlT significantly increased, and its value became much higher than the vapor flux.For example, the maximum value in P4 at the 2 cm depth reached 0.043 cm h −1 , while it was only 8.7 × 10 −4 cm h −1 in dry conditions.

Migration of Wetting Front under Different Rainfall Events
Comparing the simulation results of different periods, the impact of rainfall on soil water transfer differed significantly, and the migration of the wetting front could be used to describe the infiltration process.As depicted in Figures 5 and 6, the migration rate of the wetting front significantly changed under different rainfall intensities.For example, in P2, due to the relatively low rainfall intensity, the wetting front only moved to 28 cm at 46 pletely different.The qvT value in the shallow layer decreased significantly, which was mainly because the increase in water content reduced the hydraulic conductivity of the vapor.At this time, the influence of the vapor flux can be ignored.On the contrary, qlT significantly increased, and its value became much higher than the vapor flux.For example, the maximum value in P4 at the 2 cm depth reached 0.043 cm h −1 , while it was only 8.7 × 10 −4 cm h −1 in dry conditions.

Migration of Wetting Front under Different Rainfall Events
Comparing the simulation results of different periods, the impact of rainfall on soil water transfer differed significantly, and the migration of the wetting front could be used to describe the infiltration process.As depicted in Figures 5 and 6, the migration rate of the wetting front significantly changed under different rainfall intensities.For example, in P2, due to the relatively low rainfall intensity, the wetting front only moved to 28 cm at 46 During the dry period, the q vT flux in the top soil layer was greater than q lT , and both fluxes showed regular daily variations.With the increase in solar radiation in the daytime, the shallow temperature reduced significantly, and the downward flux value was relatively larger in comparison.The maximum q vT flux at 2 cm appeared at 14:00 daily, reaching −0.04 cm h −1 .In this condition, the q vT flux maintained was of the same order of magnitude as the q lh flux, reflecting that q vT was a significant part of the total water flux during dry periods.
After being affected by rainfall events, the changing trends in q vT and q lT were completely different.The q vT value in the shallow layer decreased significantly, which was mainly because the increase in water content reduced the hydraulic conductivity of the vapor.At this time, the influence of the vapor flux can be ignored.On the contrary, q lT significantly increased, and its value became much higher than the vapor flux.For example, the maximum value in P4 at the 2 cm depth reached 0.043 cm h −1 , while it was only 8.7 × 10 −4 cm h −1 in dry conditions.

Migration of Wetting Front under Different Rainfall Events
Comparing the simulation results of different periods, the impact of rainfall on soil water transfer differed significantly, and the migration of the wetting front could be used to describe the infiltration process.As depicted in Figures 5 and 6, the migration rate of the wetting front significantly changed under different rainfall intensities.For example, in P2, due to the relatively low rainfall intensity, the wetting front only moved to 28 cm at 46 h (10 h after the beginning of rainfall).The infiltration rate tended to decrease, and the wetting front stopped and began to fall at 100 h (located at 74 cm).On the contrary, the wetting front of P4 migrated to 68 cm at 50 h (10 h after the beginning of rainfall), and it eventually exceeded 200 cm.
To quantitatively describe the characteristics of the wetting front, Table 5 lists the calculation results for the infiltration rate and cumulative infiltration amount at different depths, in which the cumulative infiltration amount only considered the downward flux at this depth.With the increase in depth, there was a discernible declining trend in both the maximum infiltration rate and the cumulative amount.The average correlation coefficients of the maximum infiltration rate and the cumulative amount calculated during these three periods were 0.92 and −0.97, respectively, suggesting that these two variables exhibited significant variation with increasing depth.The reason for the decrease could mainly be attributed to the following two reasons.On one hand, some infiltration water was retained in the upper soil layer after redistribution, increasing the soil moisture there.Then, the evaporation process would dry the wetted soil and form a drying front when rainfall stopped, causing soil water to migrate upward again [8].On the other hand, there was a significant contrast in the calculation results of various periods due to the influence of different rainfall intensities.For example, the maximum instantaneous rate of the three periods at 10 cm were −0.037, −0.222, and −0.748 cm h −1 , respectively.Although rainfall was larger for P3 than for P2, the cumulative infiltration amount at the 50 cm depth in P2 was greater, which was mainly due to the higher soil water content in September of P2 than in April of P3.These quantitative results imply that the established model offers a suitable way to estimate the infiltration rate and amount, which are crucial for clarifying the soil hydrological cycle.When rainfall penetrated the vadose zone, soil liquid water and vapor flow were mainly impacted by two factors.As indicated by the simulation results, soil water potential and temperature gradients changed significantly when rainfall occurred, and the variations were closely related to the rainfall intensity.The driving force not only dominated the movement direction of soil liquid and vapor, but also affected the flux density.In addition, the hydraulic conductivity was another factor that affected soil water movement.Since the vapor flux has been confirmed to be mainly driven by temperature gradients, the effect of isothermal hydraulic conductivity for water vapor (K vh ) can be ignored, while the isothermal hydraulic conductivity for liquid water (K lh ) and thermal hydraulic conductivities for liquid water (K lT ) and water vapor (K vT ) play critical roles in soil liquid and vapor transfer [48,49].Because the calculations for K lh and K lT were closely related (Equation ( 2)), the variations in these two variables were similar, which were significantly increased (usually three to five orders of magnitude) when rainfall occurred.This significant increase indicated that K lh should be the dominant factor affecting the flux density during the rainfall periods, while the variation in the soil matric potential gradient played a greater role in liquid water flux on dry days.In contrast, K vT decreased in response to rainfall events, which can be attributed to the reduction in pore air content and vapor diffusivity because of the increase in soil water content [31].Unlike K lh and K lT , it was the reduction in both K vT and the soil temperature gradient that resulted in the decrease in vapor flux [48,50].
During the analyzed period, the total rainfall amount reached 640.2 mm, which far exceeds the local average rainfall.Affected by frequent rainfall events, the isothermal liquid water flux in deep layer varied significantly as well.The calculated cumulative amount at the 200 cm depth was approximately 90 mm, and most of the leakage amount would eventually recharge the groundwater.With the established model, the liquid and vapor transfer in different years with various rainfall amounts can be investigated; these are crucial for guiding agricultural irrigation and rational utilization of water resources in arid areas.

The Role of Vapor Flux in the Soil Hydrological Cycle
Based on the simulation results, the characteristics of vapor flow during both dry and rainy conditions are described in detail.Since the vapor density in the vadose zone is almost saturated most of the time, water vapor will evaporate or condense accordingly as long as the soil temperature changes [51].Driven by the temperature gradient, vapor mainly flows downward from the surface towards the deep soil layer in the daytime, while it moves upward at night.The transfer and condensation of water vapor has a significant impact on soil water flow in the upper and lower layers of the vadose zone.For the shallow layer, the vapor flux dominated by the temperature gradient is relatively larger.For instance, the average proportion of the thermal vapor flux in the total water flux accounted for 15.3% at the 2 cm depth of P1, which corresponds with the calculation results of Zhang et al. [33] and Deb et al. [51].With the increase in dry periods as well as the temperature gradient, the instantaneous proportion could reach nearly 40%, and the flux density would continue to increase.Taking the Gobi Desert of northwestern China as an example, where the annual precipitation is only 35.2 mm, Du et al. [52] found that the quantity of vapor flux exceeds the liquid flux, controlling soil water transfer.The above finding confirms that the impact of vapor flow should not be ignored in arid and semi-arid regions, especially when the soil has a long-term dry condition [32].Although the vapor flux is relatively lower in deep layers due to the decreased driving force, vapor driven by the temperature gradient continuously migrates and accumulates towards the deep layer, which should not be ignored [53].
In addition to affecting soil water transfer during the unfrozen period, the effect of vapor flow should be significant when soil is frozen [54].For the Mu Us Sandy Land and similar arid areas in northwestern China, it is noted that these regions belong to seasonally frozen areas.During the freezing period, the soil liquid flux is significantly reduced due to the presence of ice in soil pores, leading to the impact of vapor flow being more significant.The maximum frozen depth usually reached 150 cm in the study area, and vapor mainly flow downward towards the frozen layer with a lower temperature [25].Therefore, the seasonal vapor transfer during the freezing-thawing period has a significant impact on soil hydrological cycle processes, such as influencing soil moisture distribution and promoting the formation of a higher-water-content layer.

Study Limitations and Future Research Prospects
Based on field observations and experiments, a coupled water-vapor-heat transport model was established using the Hydrus-1D program in this study.For the shallow soil layer, the fitting accuracy was affected by external environmental factors.While at 100 cm depth, the relatively large discrepancy between the observed and simulated soil water content was likely caused by a sudden change in soil hydraulic properties above and below 80 cm depth due to different soil textures (Tables 1 and 2).By installing more monitoring points in the soil profile and determining soil hydraulic parameters (e.g., the vapor diffusivity in soil and the enhancement factor) using experimental rather than empirical methods, the model's accuracy is expected to be further improved [39,55].In addition, it is noted that the simulated results were obtained based on in situ measurements at a site with bare land, thus similar results as those obtained could be expected for other areas with similar soil conditions.For areas with other soil conditions or shallow groundwater levels, the established model may be inappropriate and may need to be modified.Nevertheless, the simulation and analyzed results confirmed that the model considering coupled liquid water-vapor-heat transport should be appropriate for evaluating soil water flow in arid areas.
In arid and semi-arid regions, the capacity of soil water storage and rainfall infiltration are of great significance in determining vegetation states [56].The soil layers at shallow depths are the main layers receiving precipitation and providing a regular water source for vegetation growth.In past studies, research has mainly focused on the impact of liquid water on vegetation growth, while ignoring the effect of water vapor.The condensation and aggregation of vapor in Mu Us Sandy Land has been proven to be crucial for preserving the desert ecosystem, which could provide an important water source for the growth of vegetation such as Salix psammophila and Artemisia ordosica.The impact of vegetation on soil water is mainly because of its influence on infiltration rates and transpiration, and investigating the response of vegetation to precipitation is a key challenge [57,58].A root water uptake model can clarify the influence of water uptake on the soil water flow dynamics, and a precise and proper model would be significant for investigating the vegetation water mechanism in arid regions.By coupling the root water uptake process into the coupled water, vapor, and heat transport model, the specific effect of vapor flux on plant growth and the ecological significance of soil water could be clarified in detail in future research, which would be crucial for vegetation restoration in arid areas.

Conclusions
To better understand soil liquid and vapor movement and their response to rainfall events in arid and semi-arid areas, the typical Mu Us Sandy Land in northwest China was selected as the study site, and a combination of in situ observations, field and laboratory tests, and numerical simulation methods were adopted.Based on the observed data, a combined water, vapor, and heat transport model was established using the Hydrus-1D software, and the model's accuracy was proved using the in situ observation data.The comparison of the simulation results in different periods indicated that rainfall completely changed the distribution of soil matric potential, temperature, and the driving forces, and the variations were closely correlated with the rainfall intensity.The isothermal liquid flux was the most significant component of the total water flux in the vadose zone, followed by the thermally driven liquid and vapor fluxes.The significance of the vapor flux was demonstrated by the fact that it accounted for roughly 15% of the shallow layer's total water flux in dry conditions, while the vapor flux became negligible during rainy periods.In addition to the driving forces, the significant changes in hydraulic conductivity also contributed to the variation in liquid water and vapor flux.
The findings of this research may shed further light on how the soil water cycle works in the Mu Us Sandy Land, which could provide a theoretical basis to further explore the specific effect of soil water on vegetation growth.By coupling the vegetation module into the coupled water, vapor, and heat transport process, the effect of vapor flux on plant growth and the ecological significance of soil water could be clarified in detail, which will be crucial for vegetation restoration and ecological environment protection in arid areas.

Figure 1 .
Figure 1.Location and the experiment layout of the study site.

Figure 1 .
Figure 1.Location and the experiment layout of the study site.

Figure 2 .
Figure 2. Variations in observed meteorological variables during 1 April-31 October 2018.Note P1, P2, P3, and P4 in the precipitation figure indicate the selected dry, light rainfall, moderate rainfall, and heavy rainfall periods, respectively.

where b 1 ,
b 2 , and b 3 represent empirical factors that reflect soil texture, which are defined in the Hydrus-1D program and used in the established model as 0.228, −2.406, and 4.909 W m −1 K −1 for sands, respectively.

Figure 3 .
Figure 3.Comparison of simulated and measured soil water content and temperature during the calibration period (1 April-31 October 2018).The units of MAE and RMSE for soil water content and temperature are cm 3 cm −3 and °C, respectively.

Figure 3 .
Figure 3.Comparison of simulated and measured soil water content and temperature during the calibration period (1 April-31 October 2018).The units of MAE and RMSE for soil water content and temperature are cm 3 cm −3 and • C, respectively.

Figure 4 .
Figure 4. Comparison of simulated and measured soil water content and temperature during the validation period (1 April-30 June 2019).The units of MAE and RMSE for soil water content and temperature are cm 3 cm −3 and °C, respectively.

Figure 4 .
Figure 4. Comparison of simulated and measured soil water content and temperature during the validation period (1 April-30 June 2019).The units of MAE and RMSE for soil water content and temperature are cm 3 cm −3 and • C, respectively.

Figure 6
Figure 6 presents the spatial and temporal distribution of the soil matric potential gradient in the profile.To be specific, the positive values in the figure indicate that the gradient direction ascended while the negative values indicate the opposite result.Likewise, in order to show the distribution characteristics more clearly, the maximum gradient value in the dry period was located at 40 cm cm −1 in the figure.

Figure 6
Figure 6 presents the spatial and temporal distribution of the soil matric potential gradient in the profile.To be specific, the positive values in the figure indicate that the gradient direction ascended while the negative values indicate the opposite result.Likewise, in order to show the distribution characteristics more clearly, the maximum gradient value in the dry period was located at 40 cm cm −1 in the figure.

3. 3 .
Spatial and Temporal Distribution of Soil Temperature and Its Gradients 3.3.1.Soil Temperature For the selected four typical periods, the profile's spatial and temporal distributions of soil temperature are shown in Figure 7. Agronomy 2023, 13, x FOR PEER REVIEW 12 of 22

3. 3 .
Spatial and Temporal Distribution of Soil Temperature and its Gradients 3.3.1.Soil Temperature For the selected four typical periods, the profile's spatial and temporal distributions of soil temperature are shown in Figure 7.

Figure 9 .
Figure 9.Diurnal variations in the isothermal liquid flux simulated by the Hydrus-1D model at different depths.

Figure 10 .
Figure 10.Diurnal variations in the isothermal vapor flux simulated by the Hydrus-1D model at different depths.

Figure 9 .
Figure 9.Diurnal variations in the isothermal liquid flux simulated by the Hydrus-1D model at different depths.

Figure 9 .
Figure 9.Diurnal variations in the isothermal liquid flux simulated by the Hydrus-1D model at different depths.

Figure 10 .
Figure 10.Diurnal variations in the isothermal vapor flux simulated by the Hydrus-1D model at different depths.

Figure 10 .
Figure 10.Diurnal variations in the isothermal vapor flux simulated by the Hydrus-1D model at different depths.

Figure 11 .
Figure 11.Diurnal variations in the thermal liquid flux simulated by the Hydrus-1D model at different depths.

Figure 12 .
Figure 12.Diurnal variations in the thermal vapor flux simulated by the Hydrus-1D model at different depths.

Figure 11 .
Figure 11.Diurnal variations in the thermal liquid flux simulated by the Hydrus-1D model at different depths.

Figure 11 .
Figure 11.Diurnal variations in the thermal liquid flux simulated by the Hydrus-1D model at different depths.

Figure 12 .
Figure 12.Diurnal variations in the thermal vapor flux simulated by the Hydrus-1D model at different depths.

Figure 12 .
Figure 12.Diurnal variations in the thermal vapor flux simulated by the Hydrus-1D model at different depths.

Table 1 .
Soil physical analysis in different layers of the study area.

Table 1 .
Soil physical analysis in different layers of the study area.

Table 2 .
Optimized soil hydraulic parameters of different layers.

Table 3 .
Characteristics of rainfall distribution during the simulation periods.

Table 3 .
Characteristics of rainfall distribution during the simulation periods.

Table 4 .
Detailed information for selected periods.

Table 4 .
Detailed information for selected periods.

Table 5 .
Infiltration process of different rainfall periods.