Feasibility Investigation of Improving the Modified Green – Ampt Model for Treatment of Horizontal Infiltration in Soil

Water infiltration in soil is a complex process that still requires appreciation of interactions among three phases (soil particles, water and air) to enable accurate estimation of water transport rates. To simulate this process, the Green–Ampt (GA) model and the Modified Green-Ampt (MGA) model introduced in the paper “A new method to estimate soil water infiltration based on a modified Green–Ampt model” have been widely used. The GA model is based on the hypothesis that the advance of the wetting front in soil under matric suction can be treated as a rectangular piston flow that is instantaneously transformed after passage of the infiltration front, and the MGA model does not contain the influence of pore size change. This cannot accurately reflect the soil moisture change process from unsaturation to saturation. Due to soil stratification and other inhomogeneity, predictions produced with these models often differ widely from observations. To quickly obtain the soil moisture distribution after passage of the wetting front for horizontal infiltration, an improved modified Green–Ampt (IMGA) model is presented, which estimates the soil moisture profile along a horizontal column in a piecewise manner with three functions. A logarithmic function is used to describe the gradual soil saturation process in the transmission zone, and two linear functions are used to represent the wetting zone. The algorithm of the IMGA model for estimating the water infiltration rate and cumulative infiltration is configured. To verify the effectiveness of IMGA model, a lab model test was performed, and a numerical model was built to solve the horizontal one-dimensional Richards equation using the finite–element method. The results show that the IMGA model is more accurate than the GA and MGA models. The horizontal soil moisture profiles obtained by the IMGA model are closer to the measured data than the numerical simulation results. The relative errors of the MGA and IMGA models decrease with an increase in infiltration time, whereas that of the GA model first decreases and then increases with infiltration time. The primary novelty of this study is nonlinear description of soil moisture content distribution, and derivation of unit transfer coefficient.


Introduction
The unsaturated zone of the earth's surface is an important link between the atmosphere and deeper soil profiles.Water movement and its distribution within the unsaturated zone determine the physical and hydrological properties of the soil medium and its utility in many socio-economic activities such as environmental sustainability, geo-hazards management, stability of marine dike, and beach erosion [1][2][3][4].The unsaturated ground surface is also the most important homeland for insects, and the activities of these insects are changing water infiltration model [5].Weller et al. [6] proposed that knowledge about infiltration fronts patterns is helpful in understanding nutrients and pollutants translocation.Toro-Guerrero et al. [7] pointed out that the evaluation of the soil moisture distribution and its variation are crucial for understanding water movement in landslide dams.Yeh and Tsai [8] found that when the soil moisture content increases in a slope, the effective stress and shear strength distinctly decrease.The reduction of effective stress can lead to the collapse of slopes.Therefore, evaluation of the stability of landslide at different infiltration rates and the prediction of the failure time are important for the determination of reasonable evacuation plans.
Many models have been proposed to describe water transport in soil.On the basis of their axioms and assumptions, these models have been classified into three types: empirical, semi-empirical, and physical [9].Considering that the empirical and semi-empirical models cover a limited range of scenarios without deep rational treatment, the preferential use of physical models was suggested [9].The Richards and Green-Ampt (GA) models are two most popular physical models.
With the advantages of having a reasonable physical mechanism and easy-to-solve quantitative expressions, the GA model has been widely used in water infiltration calculations for a number of hydrologic scenarios, including homogeneous soil [10], inhomogeneous soil [11][12][13][14][15], ponding sites [16], non-ponding sites [17][18][19][20], steady rainfall, and unsteady rainfall [21,22].In addition, several methods have been proposed for determining key parameters of the GA model [23][24][25][26].With respect to the direction of water movement, there are three types of infiltration: vertical downward flow, vertical upward flow, and horizontal infiltration.In engineering practice, horizontal infiltration is frequently encountered.It underlies many phenomena, such as the leakage of earth-and rock-filled dams, sand beach erosion, piping in dikes, and landslides induced by the fluctuation of reservoir water levels [27].In spite of the high frequency of these phenomena, limited research has been performed on horizontal flow models [28].
Prevedello et al. [11,29,30] proposed an analytical solution of the GA model for horizontal infiltration, which has a wider range of applications than the traditional Philip's solution [31].This analytical solution was derived from the relationship between the Boltzmann transform and the matric potential.The soil moisture profile (horizontal distribution of soil moisture) at any moment can be estimated by combining Prevedello's solution with Philip's theories [31,32].However, Prevedello et al. [11] also emphasized that this analytical solution was restricted to sands with very abrupt soil water retention curves, which exist in coarse sand but rarely found in fine sand, silt, or clay.Barry [33] suggested a power law relationship between the hydraulic head and Boltzmann variable, and extended the theories of Prevedello et al. [11,29,30].However, neither Prevedello nor Barry considered the effect of gravity on horizontal infiltration.In their studies, the only force that drives water infiltration is matric suction.They characterized the moisture profile as being composed of three regions: (i) the back transmission zone (TZ), (ii) the middle wetting zone (WZ), and (iii) the wetting front (WF) [11,29,30].It was assumed that the TZ is fully saturated throughout the infiltration process, which is very similar to the "piston" hypothesis adopted by Green and Ampt [34].However, these assumptions cannot quantitatively describe cumulative moisture variation in the TZ.The results of laboratory experiments performed by Mao et al. [35] show that the soil moisture profiles in the TZ do not satisfy Prevedello's hypothesis.They discovered that the moisture in the TZ changes throughout the entire infiltration process, gradually moving from unsaturated to a saturated state.To reflect this gradual change, Mao et al. [35] proposed a modified GA model (MGA) that uses two linear functions to piecewise fit the soil moisture profiles during the horizontal infiltration process.As a trend descriptor of the primary moisture profile, the MGA model has two limitations: (i) the MGA model does not account for the influence of adsorbed water on soil water suction; and (ii) it can be clearly seen that the moisture profiles measured by Prevedello et al. [29,30] cannot be described by the straight lines of the MGA model.Hence, some improvements are needed to fill these gaps and make more accurate estimates of soil moisture changes in soil.
The purpose of this paper is to evaluate and propose an improved MGA model (IMGA) based on the work of Mao et al. [35] using the model test results.The mathematical algorithms presented herein, were derived to estimate the water infiltration rate and cumulative infiltration.By comparing the IMGA, MGA, and GA models, the validity of the IMGA model was assessed.

Theoretical Model Development
The GA model is based on the assumption that a rectangular piston moves into homogeneous isotropic soil with a uniform initial moisture content.It is assumed that there is a sharp WF in the soil that separates the soil moisture profile into a completely saturated zone and an unsaturated zone with initial moisture θ in .Mao et al. [35] modified these assumptions and used two linear functions, θ 1 (x) and θ 2 (x), to describe the soil moisture profiles in their MGA model, as seen in Figure 1. it can be clearly seen that the moisture profiles measured by Prevedello et al. [29,30] cannot be described by the straight lines of the MGA model.Hence, some improvements are needed to fill these gaps and make more accurate estimates of soil moisture changes in soil.The purpose of this paper is to evaluate and propose an improved MGA model (IMGA) based on the work of Mao et al. [35] using the model test results.The mathematical algorithms presented herein, were derived to estimate the water infiltration rate and cumulative infiltration.By comparing the IMGA, MGA, and GA models, the validity of the IMGA model was assessed.

Theoretical Model Development
The GA model is based on the assumption that a rectangular piston moves into homogeneous isotropic soil with a uniform initial moisture content.It is assumed that there is a sharp WF in the soil that separates the soil moisture profile into a completely saturated zone and an unsaturated zone with initial moisture in θ .Mao et al. [35] modified these assumptions and used two linear functions, 1 ( ) x θ and 2 ( ) x θ , to describe the soil moisture profiles in their MGA model, as seen in Figure 1.x is advancement distance of the intersection point of two linear functions in the wetting zone and in the transmission zone (m), i x is advancement distance of the WF (m), TZ is transmission zone, WZ is wetting zone, WF is wetting front, and i means a moment during the infiltration process.
As shown in Figure 1, the soil moisture profile function in the TZ, 1 ( ) x θ , is used to estimate the soil moisture profile shape of this zone which was defined by Prevedello et al. [29,30].Soil moisture profile function, 2 ( ) x θ , is used to reflect the shape of the WZ.In the MGA model, it is assumed that after the WF, the soil moisture increase rate decreases linearly with distance of advance.This means that the matric suction after the WF, is determined by the soil particle surfaces and hypothetical constant capillary pores among the particles, which does not consider size changes of those particles.
According to the theories introduced by Hillel [36] and Lu [37], water that exists in an unsaturated soil is subject to capillarity and adsorption, which combine to produce a matric suction ψ (kPa), which can be described by the Laplace equation: As shown in Figure 1, the soil moisture profile function in the TZ, θ 1 (x), is used to estimate the soil moisture profile shape of this zone which was defined by Prevedello et al. [29,30].Soil moisture profile function, θ 2 (x), is used to reflect the shape of the WZ.In the MGA model, it is assumed that after the WF, the soil moisture increase rate decreases linearly with distance of advance.This means that the matric suction after the WF, is determined by the soil particle surfaces and hypothetical constant capillary pores among the particles, which does not consider size changes of those particles.According to the theories introduced by Hillel [36] and Lu [37], water that exists in an unsaturated soil is subject to capillarity and adsorption, which combine to produce a matric suction ψ (kPa), which can be described by the Laplace equation: where T s is surface tension of water (N/m 2 ), r soil is the radius of an idealized soil pore for fine soil (mm), and γ w is the unit weight of water.Figure 2a shows that r soil decreases with the growth of adsorbed water around soil particles, which causes ψ and the rate of the saturation process to increase.Soils at different distances from the water supply, have different values of r soil -the closer to the water supply, the smaller the value of r soil .Therefore, to improve the MGA model, θ 1 (x) should be redefined to take this effect into account.In addition, the data obtained by Prevedello et al. [29,30] indicate that the soil moisture profile in the WZ is not linear; it can be approximately fitted with two straight lines with different slopes.The ordinates of the intersections of these two lines are basically the same and are located in the moisture transition range.T is surface tension of water (N/m 2 ), soil r is the radius of an idealized soil pore for fine soil (mm), and w γ is the unit weight of water.Figure 2a shows that soil r decreases with the growth of adsorbed water around soil particles, which causes ψ and the rate of the saturation process to increase.Soils at different distances from the water supply, have different values of soil r -the closer to the water supply, the smaller the value of soil r .Therefore, to improve the MGA model, 1 ( ) x θ should be redefined to take this effect into account.In addition, the data obtained by Prevedello et al.
[ 29,30] indicate that the soil moisture profile in the WZ is not linear; it can be approximately fitted with two straight lines with different slopes.The ordinates of the intersections of these two lines are basically the same and are located in the moisture transition range.A, B, C, and D are the four points of inflection of the soil moisture profile.A logarithmic function is adopted to reflect the soil moisture change in the TZ (section AB in Figure 3): where S θ is saturated soil moisture (m 3 /m 3 ) that is a constant for a specific soil, and ( ) b t is a variable related to the elapsed time, which should be calibrated by lab tests.The cumulative infiltration can be estimated as [35]: where I is cumulative infiltration (mm). 1 S , 2 S , and 3 S can be further expressed as: where tr C is the unit transfer coefficient, which equals 1.509 here (the calculation process is introduced in Section 3.4).c θ is critical soil moisture between TZ and WZ (m 3 /m 3 ), b θ is critical soil moisture between in the middle of WZ (m 3 /m 3 ), in θ is initial soil moisture (m 3 /m 3 ), 1 ( ) In view of these observations, some improvements of the MGA model were made to transform the MGA model into the new IMGA model.Three functions, namely θ 1 (x), θ 2 (x), and θ 3 (x), are used to describe the soil moisture profiles throughout the zones, as shown in Figure 3.In this figure , A, B, C, and D are the four points of inflection of the soil moisture profile.A logarithmic function is adopted to reflect the soil moisture change in the TZ (section AB in Figure 3): where θ S is saturated soil moisture (m 3 /m 3 ) that is a constant for a specific soil, and b(t) is a variable related to the elapsed time, which should be calibrated by lab tests.The cumulative infiltration can be estimated as [35]: where I is cumulative infiltration (mm).S 1 , S 2 , and S 3 can be further expressed as: where C tr is the unit transfer coefficient, which equals 1.509 here (the calculation process is introduced in Section 3.4).θ c is critical soil moisture between TZ and WZ (m 3 /m 3 ), θ b is critical soil moisture x is advancement distance of the point of intersection of two linear functions in the wetting zone (m).

Setup of Calibration and Validation Tests
The horizontal infiltration tests were performed in a polyvinyl chloride (PVC) pipe.The pipe has a length of 8.2 m and an inner diameter of 15 cm, as shown in Figure 4.The PVC tube was divided into two sections by two filter papers that allow water to flow freely and prevent soil particles from passing through.The left space (with a length of 0.1 m) depicted in Figure 4 was a water tank for supplying water to the soil column during the following tests.The middle space with a length of 8 m was used for packing for soil samples.The right space with a length of 0.1 m was empty for draining.A piece of surface of the PVC pipe was removed, as shown in Figure 4b.Prior to testing, the soil samples were stirred by a shovel to make them homogeneous enough, and then backfilled into the tube.A carbon fiber heated cable (CFHC) was installed in the middle of the tube (Figure 4a).After all the soil samples were filled into the pipe, a transparent membrane was pasted on the ground surface to prevent water evaporation from soil and runoff on the soil surface.On the membrane, the removed PVC plate was covered by tape during the infiltration process.Prevedello et al. [29,30] proposed that the area of the WZ is constant during the infiltration process when elapsed time is large enough.Thus, the terms (x bi − x ci ), θ b +θ c 2 − θ in , (θ b − θ in ), and (x i − x bi ) are constants.Equation ( 7) is further simplified as: According to the Newton-Leibniz rule, Equation ( 8) can be rewritten as: To compare the IMGA model with the GA and MGA models, the relationship between the advancement distance x ci and the elapsed time t is expressed in the functional form that was used by Mao et al. [35]: where k 1 , k 2 , k 3 , and k where a is a constant for a specific soil (a > 0).Substituting Equations ( 11) and (10) into Equation ( 9) yields: The second term describes the influence of θ 1 (x) change rate on the water infiltration rate, which is one magnitude less than the first term of Equation ( 12), so Equation ( 12) can be further simplified as: Equations ( 13)-( 15) are the relationships between the infiltration rate and the elapsed times estimated by the IMGA, MGA, and GA models, respectively.

Setup of Calibration and Validation Tests
The horizontal infiltration tests were performed in a polyvinyl chloride (PVC) pipe.The pipe has a length of 8.2 m and an inner diameter of 15 cm, as shown in Figure 4.The PVC tube was divided into two sections by two filter papers that allow water to flow freely and prevent soil particles from passing through.The left space (with a length of 0.1 m) depicted in Figure 4 was a water tank for supplying water to the soil column during the following tests.The middle space with a length of 8 m was used for packing for soil samples.The right space with a length of 0.1 m was empty for draining.A piece of surface of the PVC pipe was removed, as shown in Figure 4b.Prior to testing, the soil samples were stirred by a shovel to make them homogeneous enough, and then backfilled into the tube.A carbon fiber heated cable (CFHC) was installed in the middle of the tube (Figure 4a).After all the soil samples were filled into the pipe, a transparent membrane was pasted on the ground surface to prevent water evaporation from soil and runoff on the soil surface.On the membrane, the removed PVC plate was covered by tape during the infiltration process.x is advancement distance of the point of intersection of two linear functions in the wetting zone (m).

Setup of Calibration and Validation Tests
The horizontal infiltration tests were performed in a polyvinyl chloride (PVC) pipe.The pipe has a length of 8.2 m and an inner diameter of 15 cm, as shown in Figure 4.The PVC tube was divided into two sections by two filter papers that allow water to flow freely and prevent soil particles from passing through.The left space (with a length of 0.1 m) depicted in Figure 4 was a water tank for supplying water to the soil column during the following tests.The middle space with a length of 8 m was used for packing for soil samples.The right space with a length of 0.1 m was empty for draining.A piece of surface of the PVC pipe was removed, as shown in Figure 4b.Prior to testing, the soil samples were stirred by a shovel to make them homogeneous enough, and then backfilled into the tube.A carbon fiber heated cable (CFHC) was installed in the middle of the tube (Figure 4a).After all the soil samples were filled into the pipe, a transparent membrane was pasted on the ground surface to prevent water evaporation from soil and runoff on the soil surface.On the membrane, the removed PVC plate was covered by tape during the infiltration process.During the tests, a distributed moisture sensing (DMS) interrogator (Model No. NZS-DMS-02, Suzhou NanZee Sensing Technology Co., Ltd, Suzhou, China) was used to automatically collect the soil moisture data.The primary performance parameters are shown in Table 1.This interrogator has the integrated functions of heating the CFHC, measuring the distribution of temperature, deducing the soil moisture, and demonstrating the results.The working principle of this interrogator is based During the tests, a distributed moisture sensing (DMS) interrogator (Model No. NZS-DMS-02, Suzhou NanZee Sensing Technology Co., Ltd, Suzhou, China) was used to automatically collect the soil moisture data.The primary performance parameters are shown in Table 1.This interrogator has the integrated functions of heating the CFHC, measuring the distribution of temperature, deducing the soil moisture, and demonstrating the results.The working principle of this interrogator is based on the actively heated fiber optic (AHFO) method [38], which is suitable not only for laboratory experiments [38,39], but also for in situ tests [40][41][42][43].Hence, the results measured with the AHFO method in these laboratory experiments can be compared with those in the field tests.The soil used in the experiments was a mixture of quartz sand and kaolin at a mass ratio of 9:1.The initial soil moisture was 0.111 m 3 /m 3 and the dry bulk density was 1.46 g/cm 3 .The initial soil moisture was the natural moisture content, which was determined by the air humidity and measured by the oven-drying method.The grain size distribution and soil water characteristic curve (SWCC) are shown in Figures 5 and 6, respectively.The soil matric suctions under different soil moistures were measured using the filter paper method [44].The SWCC was fitted using the Van Genuchten function [45].
During the test, under the pressure of the capillary force, the water in the left water supply tank entered the soil.The water supply tank was fully filled with water during the whole experiment, and the amount of water added to the water tank was recorded in real time.In general, the influence of water level on the horizontal infiltration is not considered.The horizontal infiltration experiments are often conducted in a liner tube that the influence of tube diameter can be neglected.Thus, the smaller the linear tube, the less the influence of the water level.In this study, the diameter (15 cm) is much smaller than the length of the soil column (8 m), thus, the impact of water tank size is not discussed.The soil moisture distribution was measured once every hour.The total time for every measurement is 20 min.The position of the WF was measured using a ruler with an accuracy of ±1 mm.The total period of each experiment was 24 h, and the soil moisture data were collected once per hour.on the actively heated fiber optic (AHFO) method [38], which is suitable not only for laboratory experiments [38,39], but also for in situ tests [40][41][42][43].Hence, the results measured with the AHFO method in these laboratory experiments can be compared with those in the field tests.The soil used in the experiments was a mixture of quartz sand and kaolin at a mass ratio of 9:1.The initial soil moisture was 0.111 m 3 /m 3 and the dry bulk density was 1.46 g/cm 3 .The initial soil moisture was the natural moisture content, which was determined by the air humidity and measured by the oven-drying method.The grain size distribution and soil water characteristic curve (SWCC) are shown in Figures 5 and 6, respectively.The soil matric suctions under different soil moistures were measured using the filter paper method [44].The SWCC was fitted using the Van Genuchten function [45].During the test, under the pressure of the capillary force, the water in the left water supply tank entered the soil.The water supply tank was fully filled with water during the whole experiment, and the amount of water added to the water tank was recorded in real time.In general, the influence of water level on the horizontal infiltration is not considered.The horizontal infiltration experiments are often conducted in a liner tube that the influence of tube diameter can be neglected.Thus, the smaller the linear tube, the less the influence of the water level.In this study, the diameter (15 cm) is much smaller than the length of the soil column (8 m), thus, the impact of water tank size is not discussed.The soil moisture distribution was measured once every hour.The total time for every measurement is 20 min.The position of the WF was measured using a ruler with an accuracy of ±1 mm.The total period of each experiment was 24 h, and the soil moisture data were collected once per hour.

Parameter Determination Method
When determining soil moisture using the DMS, it is needed to calibrate the parameters.Soil moisture ( θ ) was calculated according to the temperature characteristic value ( t T ), which was defined as the average temperature during a time period in the heating process [46].Sayde et al. [38,47] suggested adopting a piecewise function to describe the relationship between t T and θ .In this study, if the soil moisture was larger than the threshold moisture ( 0 θ ), a logarithmic function was used; otherwise, a linear function was adopted, as expressed in Equation 16: where 1 y , 2 y , 3 y 4 y , and 5 y are parameters determined by soil grain sizes, pore structure, composition, etc.In order to obtain these parameters, the actually measured soil moisture was measured by the oven-drying method via collected soil samples.The soil samples at different positions were collected every 30 min in the first four hours of the experiment.The relationship between t T directly recorded by interrogator and θ was fitted using the least square method.Then, the soil moisture during the whole experiment was inferred from the calibrated relationship.The constants in the functions 1 ( ) x θ , 2 ( ) x θ , and 3 ( ) x θ were determined by the method of undetermined coefficients.After obtaining the actual measured soil moisture distribution by the DMS, the inflection points can be located through observation.b θ and c θ were the soil moisture of the inflection points.Thus the coordinates of the points A, B, C and D in Figure 3 could be determined.The slopes of lines BC and CD in Figure 3 could be calculated through the coordinates of the points Figure 6.Soil water characteristic curves for the soil samples.θ S (drying) and θ S (wetting) are the saturated soil moistures under drying and wetting paths (m 3 /m 3 ), respectively; θ r is the residual moisture (m 3 /m 3 ); ψ b is the air entry pressure (kPa); and ψ w is the water entry pressure (kPa).

Parameter Determination Method
When determining soil moisture using the DMS, it is needed to calibrate the parameters.Soil moisture (θ) was calculated according to the temperature characteristic value (T t ), which was defined as the average temperature during a time period in the heating process [46].Sayde et al. [38,47] suggested adopting a piecewise function to describe the relationship between T t and θ.In this study, if the soil moisture was larger than the threshold moisture (θ 0 ), a logarithmic function was used; otherwise, a linear function was adopted, as expressed in Equation ( 16): where y 1 , y 2 , y 3 , y 4 , and y 5 are parameters determined by soil grain sizes, pore structure, composition, etc.In order to obtain these parameters, the actually measured soil moisture was measured by the oven-drying method via collected soil samples.The soil samples at different positions were collected every 30 min in the first four hours of the experiment.The relationship between T t directly recorded by interrogator and θ was fitted using the least square method.Then, the soil moisture during the whole experiment was inferred from the calibrated relationship.The constants in the functions θ 1 (x), θ 2 (x), and θ 3 (x) were determined by the method of undetermined coefficients.After obtaining the actual measured soil moisture distribution by the DMS, the inflection points can be located through observation.θ b and θ c were the soil moisture of the inflection points.Thus the coordinates of the points A, B, C and D in Figure 3 could be determined.The slopes of lines BC and CD in Figure 3 could be calculated through the coordinates of the points A, B, C and D. The functions θ 1 (x), θ 2 (x), and θ 3 (x) could be further determined using the method of undetermined coefficients.The calibration operation can be referred to [39].

Numerical Model
The theoretical soil moisture profile at any moment can be obtained by solving the Richards' equation.In this study, only the horizontal mode is considered, which can be expressed as [37]: Water 2019, 11, 645 where D(θ) is the soil water diffusivity (m 2 /s), t is elapsed time (h), and x is horizontal distance (m).The solution of Equation ( 17) is obtained by use of a finite element-based numerical method.The numerical model, which has the same size as the physical model in Figure 4, was built using the finite element method.The numerical model uses Richards' equation as its governing equation.The unit grid is an equilateral triangle with a length of 0.015 m.The Brooks & Corey function and Van Genuchten function has been the most common models used in solving the Richards' equation [45,48].In this paper, the Van Genuchten function was selected.The parameters of the Van Genuchten function are listed in Table 2, which are determined using the method that was introduced by Carsel & Parrish [49].

Parameters Calibration Results
The specific expression of Equation ( 16) were fitted as Equation ( 18) with an R 2 of 0.987: The horizontal soil moisture profiles at different infiltration moment were shown in Figure 7.The WZ during the first 2 h is not a constant that is not discussed in this study.Averaging the values of θ b and θ c that were measured during the period of 2-24 h in Figure 7

Horizontal distance (m) 0 h 1 h 2 h 4 h 6 h 9 h 12 h 18 h 22 h 24 h θc zone
According to Equation ( 15), 1 ( ) x θ can be calculated at any moment.This means that the soil moisture profiles are only determined by ( ) b t during this period.

Infiltration Rate
By substituting the parameters that were calibrated in Section 2.3 into Equations ( 13)-( 15), the following expressions are obtained: The parameters k 1 , k 2 , k 3 , and k 4 in Equation (10) were fit using the least-squares method, as shown in Figure 8.The fitted result, with R 2 = 0.981, is expressed as: It should be pointed out that only for the soil used in this study, k 2 equals 0. The unit of hour is used for t.For the soil used in this study, by substituting θ b = 0.25 m 3 /m 3 and θ c = 0.33 m 3 /m 3 into Equation ( 2), the following expression is obtained: According to Equation ( 15), θ 1 (x) can be calculated at any moment.This means that the soil moisture profiles are only determined by b(t) during this period.

Infiltration Rate
By substituting the parameters that were calibrated in Section 2.3 into Equations ( 13)-( 15), the following expressions are obtained: Equations ( 21)-( 23) indicate that, at any time, the infiltration rate estimated by the GA model is the largest, followed by those of the IMGA and MGA models.This is due to the difference in the assumed soil moisture profiles, as shown in Figure 9. ( Equations ( 21)- (23) indicate that, at any time, the infiltration rate estimated by the GA model is the largest, followed by those of the IMGA and MGA models.This is due to the difference in the assumed soil moisture profiles, as shown in Figure 9.

Cumulative Infiltration
Since the analytical integral process of Equation ( 8) is complicated and time consuming in applications, Equations ( 21)-( 23) were directly integrated to obtain the cumulative infiltration in this study: Mao et al. [35] suggested that the error of the MGA model could be evaluated by comparing the predicted and actual values, but they did not specify which moment should be selected.To fill this gap, the predicted values of the IMGA, MGA, and GA models and the actual values that were measured using the graduated cylinder throughout the entire experimental process were selected, as shown in Figure 10.
It can be seen from Figure 10 that the deviations between the estimated cumulative infiltrations of the various models and the measured ones are different at different times.To quantify the deviation, the relative error ( δ ), which was configured by Mao et al. [35], was used.It is defined as: where J is the cross-sectional area of the soil column (m 2 ) packed by the PVC tube, V is the total supplied water volume (m 3 /m 3 ), and δ is the relative error.The comparison of δ among the IMGA, MGA, and GA models is shown in Figure 10.

Cumulative Infiltration
Since the analytical integral process of Equation ( 8) is complicated and time consuming in applications, Equations ( 21)-( 23) were directly integrated to obtain the cumulative infiltration in this study: Mao et al. [35] suggested that the error of the MGA model could be evaluated by comparing the predicted and actual values, but they did not specify which moment should be selected.To fill this gap, the predicted values of the IMGA, MGA, and GA models and the actual values that were measured using the graduated cylinder throughout the entire experimental process were selected, as shown in Figure 10. Figure 10 shows that at the initial time, t = 0, the initial δ values for these models are very large.This is attributed to the weakness of this evaluation method.When the infiltration time is 0, the estimated cumulative infiltration is equal to 0 because the water has not flowed.However, it is necessary to add a certain amount of water to the tank to keep the water level at the required height.Therefore, there is a difference between the estimated and measured water volumes.When the elapsed time is larger than 0, the test periods can be divided into two stages: I and II, according to the relative error trend.The separation between these stages is at hour three.During the first stage, the relative errors of the three models decrease rapidly with time.During the second stage, for the GA model, δ increases with time, but for the MGA and IMGA models, it decreases with time.It can It can be seen from Figure 10 that the deviations between the estimated cumulative infiltrations of the various models and the measured ones are different at different times.To quantify the deviation, the relative error (δ), which was configured by Mao et al. [35], was used.It is defined as: where J is the cross-sectional area of the soil column (m 2 ) packed by the PVC tube, V is the total supplied water volume (m 3 /m 3 ), and δ is the relative error.The comparison of δ among the IMGA, MGA, and GA models is shown in Figure 10.
Figure 10 shows that at the initial time, t = 0, the initial δ values for these models are very large.This is attributed to the weakness of this evaluation method.When the infiltration time is 0, the estimated cumulative infiltration is equal to 0 because the water has not flowed.However, it is necessary to add a certain amount of water to the tank to keep the water level at the required height.Therefore, there is a difference between the estimated and measured water volumes.When the elapsed time is larger than 0, the test periods can be divided into two stages: I and II, according to the relative error trend.The separation between these stages is at hour three.During the first stage, the relative errors of the three models decrease rapidly with time.During the second stage, for the GA model, δ increases with time, but for the MGA and IMGA models, it decreases with time.It can also be found that at any moment, δ of the IMGA model is less than that of the MGA model.Therefore, the elapsed time and experimental conditions should be considered when using the method proposed by Mao et al. [35] to evaluate the models' accuracies.
Prevedello et al. [29,30] found that the TZ zone dominates the cumulative infiltration.The GA model assumes that the TZ is saturated at all times, which contradicts the actual measured values in Figure 7. Hence, with increased time, the differences between the GA model and the measured values accumulate, leading to increasingly large values of δ.Both the MGA and IMGA models can precisely account for the true soil moisture change in the TZ.So, with increased time, the influence of the initial quantity of water that was added decreases.When calculating the infiltration rate and cumulative infiltration, the unit transfer coefficient is an essential parameter.Therefore, it is important to accurately determine this parameter.

Unit Transfer Coefficient C tr
Mao et al. [35] took C tr = 10 in their study; however, they did not describe the specific calculation process of this parameter.C tr is determined by many factors, such as particle size, gradation, mineral composition, dry density, and model size.In this study, C tr was calibrated by laboratory experiments using the GA model.In the GA model, it is assumed that the TZ and WZ are saturated, so C tr can be calculated based on the mass conservation law.At moment t, the following holds: The general expressions of Equations ( 23) and ( 26) are: where E and F are constants.Here, t = 24 h was taken as an example to illustrate the derivation process.It can be inferred from Equation ( 18) that x ci (24 h) = 7.87 m.By substituting x ci (24 h) = 7.87 m into Equations (29), it was found that E = 3.813, as shown in Equation (26).By calculating the different coefficients of Equation ( 30), it was found that F = 0.321, as shown in Equation (23).Equation ( 13) can be further rewritten as: Water 2019, 11, 645 13 of 19 when t = 24 h, it can be determined from Equation (31) that i GA (24 h) is 0.04275 mm.Considering that the term θ s − θ in is a constant equal to 0.31 m 3 /m 3 , we obtained C tr = 1.509.C tr is an indicator of water driving force in TZ.According to the theories introduced by Hillel [36] and Lu [37], in the unsaturated soil without consideration of gravity, the driving force is soil suction.Thus, C tr is determined by many factors, such as absolute temperature, specific volume of water, water density, molecular mass of water vapor, partial pressure of pore-water vapor, saturation pressure of water vapor over a flat surface of pure water, pore structure, and grain size.In general, the smaller the soil particle is, the larger the C tr [36].

Moisture Profiles of Different Soils during the Infiltration Process
After investigating the infiltration characteristics of three soils, Mao et al. [35] found that different soils have different parameters in the MGA model.To verify the validity and reliability of the IMGA model in different soils, the data collected by Mao et al. [35] was fitted with the parameters of the IMGA model, as shown in Figure 11.
where E and F are constants.Here, 24 h t = was taken as an example to illustrate the derivation process.It can be inferred from Equation ( 18) that  29), it was found that E = 3.813, as shown in Equation (26).By calculating the different coefficients of Equation (30), it was found that F = 0.321, as shown in Equation (23).Equation (13)  of water driving force in TZ.According to the theories introduced by Hillel [36] and Lu [37], in the unsaturated soil without consideration of gravity, the driving force is soil suction.Thus, tr C is determined by many factors, such as absolute temperature, specific volume of water, water density, molecular mass of water vapor, partial pressure of pore-water vapor, saturation pressure of water vapor over a flat surface of pure water, pore structure, and grain size.In general, the smaller the soil particle is, the larger the tr C [36].

Moisture Profiles of Different Soils during the Infiltration Process
After investigating the infiltration characteristics of three soils, Mao et al. [35] found that different soils have different parameters in the MGA model.To verify the validity and reliability of the IMGA model in different soils, the data collected by Mao et al. [35] was fitted with the parameters of the IMGA model, as shown in Figure 11.As shown in Figure 11, the results of the IMGA model (red curves) fit well with the measured data.Different soils have different line shapes with different parameters, as listed in Table 3.It can be seen that for each soil, the R 2 of the IMGA model is larger than that of the MGA model, indicating that the former has a higher accuracy.The differences of the parameters among different soils suggest that these parameters should be calibrated for each soil.As shown in Figure 11, the results of the IMGA model (red curves) fit well with the measured data.Different soils have different line shapes with different parameters, as listed in Table 3.It can be seen that for each soil, the R 2 of the IMGA model is larger than that of the MGA model, indicating that the former has a higher accuracy.The differences of the parameters among different soils suggest that these parameters should be calibrated for each soil.Both the IMGA model and Richards' equation are physical models that describe horizontal infiltration of soil.To analyze and investigate their differences, the soil moisture profiles calculated by these two methods are shown in Figure 12.
Both the IMGA model and Richards' equation are physical models that describe horizontal infiltration of soil.To analyze and investigate their differences, the soil moisture profiles calculated by these two methods are shown in Figure 12.As shown in Figure 12a, the numerical solutions of Richards' equation do not conform to the measured values in two respects.First, the shapes of horizontal soil moisture profiles obtained by Richards' equation and measured data are different.It can be distinctly pick out the TZ, WZ, and WF on the measured curves; however, the smooth simulation curves have no obvious demarcations (critical soil moisture), which is agrees well with the results of [50][51][52].Mao et al. [35] and Herrada et al. [53] proposed that although the simplified numerical simulations can produce analytical solutions, they may lead to significant deviations from actually measured data.Second, the positions of the WF that were estimated by the numerical solutions of Richards' equation differ from the observed positions.For the sake of analyzing the difference between these WF positions, a new variable x Δ (m) was used, which was defined as: where N x is the advancement distance is estimated by Richards' equation (m) and M x is the measured distance (m).Substituting N x and M x , at the moments 2 h, 4 h, 6 h, 12 h, and 24 h into equation (32) yields: Figure 12b shows that the soil profiles obtained by the IMGA model is closer than that of the simulated results of Richards equation.This is because when solving Richards' equation, the Van Genuchten [45] function was used to describe the relationship between matric suction and moisture content.This function ascribes water infiltration to a matric suction imbalance without any consideration of the influence of gravity, pore size distribution, or interface.As a matter of fact, soils with different compositions, structures, and bulk densities have different SWCCs.In addition, it is difficult to ensure that the soil in the infiltration test is absolutely uniform and that the conditions in the infiltration test are the same as those under which the SWCC was measured.Mao et al. [35] also found that in porous soil, water tends to move in the direction of decreasing potential energy, and the gradient of potential energy with distance is the primary force of water flow.Therefore, various uncontrollable factors can cause the soil matric suction to change, and then lead to errors in the values predicted using Richards' equation [54,55].As shown in Figure 12a, the numerical solutions of Richards' equation do not conform to the measured values in two respects.First, the shapes of horizontal soil moisture profiles obtained by Richards' equation and measured data are different.It can be distinctly pick out the TZ, WZ, and WF on the measured curves; however, the smooth simulation curves have no obvious demarcations (critical soil moisture), which is agrees well with the results of [50][51][52].Mao et al. [35] and Herrada et al. [53] proposed that although the simplified numerical simulations can produce analytical solutions, they may lead to significant deviations from actually measured data.Second, the positions of the WF that were estimated by the numerical solutions of Richards' equation differ from the observed positions.For the sake of analyzing the difference between these WF positions, a new variable ∆x (m) was used, which was defined as: Figure 12b shows that the soil profiles obtained by the IMGA model is closer than that of the simulated results of Richards equation.This is because when solving Richards' equation, the Van Genuchten [45] function was used to describe the relationship between matric suction and moisture content.This function ascribes water infiltration to a matric suction imbalance without any consideration of the influence of gravity, pore size distribution, or interface.As a matter of fact, soils with different compositions, structures, and bulk densities have different SWCCs.In addition, it is difficult to ensure that the soil in the infiltration test is absolutely uniform and that the conditions in the infiltration test are the same as those under which the SWCC was measured.Mao et al. [35] also found that in porous soil, water tends to move in the direction of decreasing potential energy, and the gradient of potential energy with distance is the primary force of water flow.Therefore, various uncontrollable factors can cause the soil matric suction to change, and then lead to errors in the values predicted using Richards' equation [54,55].

Conclusions
In this paper, the soil moisture profiles were measured by using actively heated fiber optic technology in laboratory, and some improvements have been made to the Green-Ampt (GA) and modified Green-Ampt (MGA) models.A new model, i.e., the improved modified Green-Ampt (IMGA) model, has been proposed to describe water infiltration in landslide dams.This model estimates the horizontal soil moisture profiles using three functions.A logarithmic function was used to reflect the soil moisture profile in the transition zone, and two linear functions were proposed for the wetting zone.The relationship between the advancement distance of the wetting front and time was fitted with an exponential function.After comparing and analyzing the IMGA, MGA, GA models, and the solutions of the Richards Equation, the following conclusions are drawn:

•
The IMGA model enables accurate estimation of the soil moisture profile at any moment.The logarithmic function adopted by the IMGA model can capture changes in moisture content in the transition zone with a reflection of its decrease as the wetting front advances.

•
By comparing the estimated cumulative infiltration data obtained using the models with those obtained by measurements, it was found that the IMGA model has the lowest relative error, followed by the MGA and GA models.It was also discovered that the relative errors of the models change with the infiltration time.The relative errors of the IMGA and MGA models monotonically decrease with time; however, for the GA model, the relative error first increases and then decreases.Therefore, the infiltration time should be specified when using this method to evaluate the errors of the models.

•
It was demonstrated that the unit transfer coefficient C tr can be inferred from the GA model using the mass conservation law.Soils with different compositions, grain size and structures also have different C tr values, which can be determined by the method suggested in this study.The smaller the grain size, the larger C tr is.

•
Due to the unavoidable influence of interfaces and inhomogeneity in soil, the advancement distance of water infiltration calculated using Richards' equation is shorter than that estimated by the IMGA model and the measured values.Another feature of the numerical solutions of Richards' equation is that it is unable to reflect transmission zone, wetting zone, and demarcation point between these two zones.The horizontal soil moisture profiles obtained by the IMGA model is closer to the measured data than that of the numerical simulation results.
Although the IMGA model enables an accurate description of the horizontal soil moisture profiles during an infiltration process, similar to the MGA and GA models, there are still some deficiencies that limit its further application.For example, it cannot predict the infiltration in heterogeneous soil, which exists extensively in nature, or estimate two-or three-dimensional soil moisture profiles.

Nomenclature and Abbreviations
The following notations are used in this paper: water advancement distance at the bottom of the polyvinyl chloride (PVC) column (m)

Nomenclature
x bi advancement distance of the intersection point between the transmission zone and the wetting zone (m) x ci advancement distance of the intersection point of two linear functions in the wetting zone (m) x i advancement distance of the WF (m) x T water advancement distance at the top of the PVC column (m) y [1][2][3][4] parameters that described relationship between T t and θ Greek letters

Figure 1 .
Figure1.Assumed soil moisture profiles in the Green-Ampt (GA) and modified Green-Ampt (MGA) models.θ S is saturated soil moisture (m 3 /m 3 ), θ c is critical soil moisture between transmission zone (TZ) and wetting zone (WZ) (m 3 /m 3 ), θ in is initial soil moisture (m 3 /m 3 ), θ 1 (x) is soil moisture profile function in the transmission zone (m 3 /m 3 ), θ 2 (x) soil moisture profile function connected with TZ in the wetting zone (m 3 /m 3 ), x ci is advancement distance of the intersection point of two linear functions in the wetting zone and in the transmission zone (m), x i is advancement distance of the WF (m), TZ is transmission zone, WZ is wetting zone, WF is wetting front, and i means a moment during the infiltration process.

Figure 2 .
Figure 2. Water in an unsaturated soil: (a) actual soil moisture distribution introduced by Hillel [36] after a wetting front; (b) soil moisture distribution described by the GA model; and (c) soil moisture distribution described by the MGA model.

Figure 2 .
Figure 2. Water in an unsaturated soil: (a) actual soil moisture distribution introduced by Hillel [36] after a wetting front; (b) soil moisture distribution described by the GA model; and (c) soil moisture distribution described by the MGA model.

7 ) 19 Figure 3 .
Figure 3. Soil moisture profiles simulated with the improved modified Green-Ampt (IMGA) model over time (the dashed line is the soil moisture profile at time i, and the solid line is the moisture profile at a later time f).bθ is critical soil moisture in the middle of WZ (m 3 /m 3 ), 3( )  x θ

Figure 4 .
Figure 4. Configuration of the experimental apparatus: (a) soil column and cross section of the carbon

Figure 3 .
Figure 3. Soil moisture profiles simulated with the improved modified Green-Ampt (IMGA) model over time (the dashed line is the soil moisture profile at time i, and the solid line is the moisture profile at a later time f ).θ b is critical soil moisture in the middle of WZ (m 3 /m 3 ), θ 3 (x) is soil moisture profile function connected with WF in the wetting zone (m 3 /m 3 ), x bi is advancement distance of the point of intersection of two linear functions in the wetting zone (m).
4 are constants related to soil properties and the water head difference, which should be calibrated by experiments.b(t) denotes the soil moisture change rate in the TZ, which decrease with time.Here, a quadratic polynomial is suggested to be used to describe the relationship between b(t) and t: b(t) = −at 2(11)

Figure 3 .
Figure 3. Soil moisture profiles simulated with the improved modified Green-Ampt (IMGA) model over time (the dashed line is the soil moisture profile at time i, and the solid line is the moisture profile at a later time f).bθ is critical soil moisture in the middle of WZ (m 3 /m 3 ), 3( )  x θ

Figure 4 .
Figure 4. Configuration of the experimental apparatus: (a) soil column and cross section of the carbon fiber heated cable (CFHC); (b) cut region of the polyvinyl chloride (PVC) tube.

Figure 4 .
Figure 4. Configuration of the experimental apparatus: (a) soil column and cross section of the carbon fiber heated cable (CFHC); (b) cut region of the polyvinyl chloride (PVC) tube.

Water 2019 ,
11, x FOR PEER REVIEW 7 of 19
soil moistures under drying and wetting paths (m 3 /m 3 ), respectively; r θ is the residual moisture (m 3 /m 3 ); b ψ is the air entry pressure (kPa); and w ψ is the water entry pressure (kPa).
produced θ b = 0.25 m 3 /m 3 and θ c = 0.33 m 3 /m 3 .The slopes of lines BC and CD in Figure 3 are k BC =0.09 and k CD =0.35 with R 2 of 0.987 and 0.975.θ s and θ in were obtained by the drying method; here, they are 0.42 m 3 /m 3 and 0.11 m 3 /m 3 , respectively.It was further determined that (x bi − x ci ) =1.5 m, (x i − x bi ) = 0.4 m , θ b +θ c 2 − θ i = 0.19 m 3 /m 3 , and (θ b − θ in ) = 0.14 m 3 /m 3 .The relationship between advancement distance of the wetting front and elapsed time is shown in Figure 8. Water 2019, 11, x FOR PEER REVIEW 10

Figure 7 .
Figure 7. Soil moisture profiles measured by the DMS interrogator.c θ zone is the soil moisture range c θ located in, and b θ zone is the soil moisture range b θ located in.

Figure 7 .Figure 7 .
Figure 7. Soil moisture profiles measured by the DMS interrogator.θ c zone is the soil moisture range θ c located in, and θ b zone is the soil moisture range θ b located in.

Figure 8 .
Figure 8. Advancement distance of the wetting front with elapsed time.

Figure 8 .
Figure 8. Advancement distance of the wetting front with elapsed time.

Figure 9 .
Figure 9. Water infiltration rates estimated by the IMGA, MGA, and GA models.

Figure 9 .
Figure 9. Water infiltration rates estimated by the IMGA, MGA, and GA models.

Figure 10 .
Figure10.Relative errors of the IMGA, MGA, and GA models as a function of time.I is the period in which the error of GA model degreases.II is the period in which the error of GA model increases.

Figure 10 .
Figure10.Relative errors of the IMGA, MGA, and GA models as a function of time.I is the period in which the error of GA model degreases.II is the period in which the error of GA model increases.

Figure 11 .
Figure 11.Soil moisture profiles measured by Mao et al. [35] and the fit results of the IMGA model.BJ, AS and YL represent soil samples collected in Beijing, Ansai, and Yangling (China), respectively.

Figure 11 .
Figure 11.Soil moisture profiles measured by Mao et al. [35] and the fit results of the IMGA model.BJ, AS and YL represent soil samples collected in Beijing, Ansai, and Yangling (China), respectively.

Figure 12 .
Figure 12.Soil moisture profile comparison: (a) wetting front comparison between actually measured values (M) and numerical solutions (N) of Richards' equation, and (b) soil moisture profile comparison among the IMGA model, solutions of Richards' equation and actually measured values.

Figure 12 .
Figure 12.Soil moisture profile comparison: (a) wetting front comparison between actually measured values (M) and numerical solutions (N) of Richards' equation, and (b) soil moisture profile comparison among the IMGA model, solutions of Richards' equation and actually measured values.

θ 1
(x) soil moisture profile function in the transmission zone (m 3 /m 3 ) θ 2 (x) soil moisture profile function collected with TZ in the wetting zone (m 3 /m 3 ) θ 3 (x) soil moisture profile function collected with WF in the wetting zone (m 3 /m 3 ) θ b critical soil moisture in the middle of WZ (m 3 /m 3 ) θ c critical soil moisture between TZ and WZ (m 3 /m 3 ) θ in initial soil moisture (m 3 /m 3 ) θ r residual moisture (m 3 /m 3 ) θ S (drying) saturated soil moisture under a drying path (m 3 /m 3 ) θ S (wetting) saturated soil moisture under a wetting path (m 3 /m

Table 2 .
Parameter values for the van Genuchten model function.
can be further rewritten as:

Table 3 .
Parameter values for different soils.BJ, AS and YL represent soil samples collected in Beijing, Ansai, and Yangling (China).

Table 3 .
Parameter values for different soils.BJ, AS and YL represent soil samples collected in Beijing, Ansai, and Yangling (China).