# Segmental Modification of the Mualem Model by Remolded Loess

Authors: Le-fan Wang; Xing-zhong Weng; Zhi-hua Yao; Ren-yi Zhang; Wan Li; Ye Li
Journal: Mathematical Problems in Engineering
 (2017)

---

## Abstract

The measured diffusion coefficient and soil-water characteristic curve (SWCC) of remolded loess were used to modify the Mualem model for increasing its accuracy. The obtained results show that the goodness of fit between the Mualem model and the variable parameter-modified Mualem method comparing with the test results was not high. The saturation of 0.65 was introduced as the boundary to divide the curve of the measured diffusion coefficient into two segments. When the segmentation method combined with the variable parameter method was used to modify the Mualem model, the fitting correlation coefficient was increased to 0.921–0.998. The modified parametersKo and L corresponding to remolded loess were calculated for different dry densities. Based on the exponential function between Ko and dry density and the linear relation between L and dry density, the segmentally modified Mualem model was established for remolded loess by considering variation in dry density. The results of the study can be used for directly determining the unsaturated infiltration coefficient and for indirectly determining the SWCC through diffusion coefficient.


---

## Body

## 1. Introduction

The unsaturated infiltration coefficient is an important content in the study of unsaturated soil and is commonly predicted by using the soil-water characteristic curve (SWCC). Childs and Collis-George [1], Marshall [2], Xu [3], and Mei et al. [4] have achieved lots of beneficial results in the field of unsaturated infiltration coefficient prediction. Mualem [5] proposed a model, widely known as the Mualem model, for unsaturated infiltration coefficient prediction by studying 45 groups of undisturbed soil and remolded soil in 1976, which has been widely applied.The saturated infiltration coefficientKs is an important parameter of the Mualem model. Van Genuchten and Nielsen [6] believed that Ks mainly reflects the characteristic of moisture flow in soil pores incurred by external force, while unsaturated seepage is mainly caused by the matrix suction inside the soil mass. Schaap and Leij [7], by analyzing 235 groups of soil sample, stated that the saturated infiltration coefficient Ks in the Mualem model should be replaced by the modified infiltration coefficient Ko and that Ko is one order of magnitude smaller than Ks.Another important parameterL of the Mualem model was originally fixed at 0.5 [5]. However, Yates et al. [8] found the L value can be a variable greater than −3.31 by reanalyzing the data of Mualem in 1976. Likewise, Schuh and Cline [9] obtained the L value between −8.73 and 14.80 through analyzing 75 groups of data. It is found in these papers [6–9] that setting all the parameters of the Mualem model as fixed values reduces the model’s accuracy in predicting the unsaturated infiltration coefficient. Valiantzas [10] combined Brooks-Corey and Van Genuchten closed-form models improving prediction of conductivity values near the residual water content and saturation. Hu et al. [11] developed a new soil-water characteristic curve model for deformable soils. Zhang et al. [12] improved the accuracy of saturated hydraulic conductivity in the Mualem model.The unsaturated infiltration coefficient can also be of either the SWCC method or the diffusion coefficient method [13]. Wang et al. [14] and Yao et al. [15] used this method to analyze the infiltration characteristic of loess. Guided by the idea of integrating the above two methods of determining the unsaturated infiltration coefficient, this study conducted a horizontal one-dimensional (1D) seepage test to measure the diffusion coefficient and used SWCCs of remolded loess of different dry densities to modify the Mualem model. Subsequently, this paper proposes a dry density-based, remolded loess-modified Mualem model.


## 2. Theoretical Basis

### 2.1. Mualem Model and Proposed Modification of Its Variable Parameters

The Mualem model for predicting the unsaturated infiltration coefficient is presented as follows:(1)Kθ=KsSe0.5∫0SeψSe-1dSe∫01ψSe-1dSe2,where Ks is the saturated infiltration coefficient and ψ(Se) is the matrix suction corresponding to the effective saturation Se.On the other hand, Hoffmann-Riem proposed a universal model for predicting the unsaturated infiltration coefficient in 1999 as follows:(2)Kθ=KoSeL∫0Seψ-βdSe∫01ψ-βdSeγ.The present paper setβ=1 and γ=2 in (2) in correspondence with the Mualem model. This is equivalent to introducing two variable parameters, Ko and L, in the original Mualem model to obtain a diffusion coefficient-modified Mualem model.


### 2.2. Determination of the Diffusion Coefficient by Horizontal 1D Seepage Test [13]

Without considering the action of gravity or the deformation of soil mass, the differential equation and boundary conditions for moisture movement in the horizontal 1D seepage test without free head are shown as follows:(3)∂θ∂t=∂∂Dθ∂θ∂x,(4)θ=θax>0,t=0,(5)θ=θox=0,t>0,where θ is the volumetric moisture content (cm3·cm−3), Dθ is the diffusion coefficient (cm2·min−1), x is the distance from the water inlet (cm), θa and θo are, respectively, the initial volumetric moisture content of soil sample and the volumetric moisture content at water inlet (cm3·cm−3), and t is the infiltration time (min).After measuring the distribution of the volumetric moisture contentθ of soil columns at time t, λθ~θ can be calculated according to the Boltzmann parameter as λθ=xt-0.5. The diffusion coefficient Dθ can be calculated as follows: (6)Dθ=-12ΔλΔθ∑θaθλdθ.The relationship among the diffusion coefficientDθ, the SWCC ψθ, and the infiltration coefficient Kθ is given as(7)Kθ=-Dθ·dθdψθ.After measuring the diffusion coefficient and the SWCC, the unsaturated infiltration coefficient can be calculated from (7).


### 2.3. SWCC Model

In order to modify the Mualem model, the first step is to determine the SWCC model. For this purpose, the present paper adopted the Brook and Corey model, as shown below:(8)θ-θrθs-θr=Se=1,ψ≤ψa,(9)θ-θrθs-θr=Se=ψaψσ,ψ>ψa,where σ is the shape parameter of the SWCC, θs is the saturated volumetric moisture content, θr is the residual volumetric moisture content, and ψa is the air entry value.


## 2.1. Mualem Model and Proposed Modification of Its Variable Parameters

The Mualem model for predicting the unsaturated infiltration coefficient is presented as follows:(1)Kθ=KsSe0.5∫0SeψSe-1dSe∫01ψSe-1dSe2,where Ks is the saturated infiltration coefficient and ψ(Se) is the matrix suction corresponding to the effective saturation Se.On the other hand, Hoffmann-Riem proposed a universal model for predicting the unsaturated infiltration coefficient in 1999 as follows:(2)Kθ=KoSeL∫0Seψ-βdSe∫01ψ-βdSeγ.The present paper setβ=1 and γ=2 in (2) in correspondence with the Mualem model. This is equivalent to introducing two variable parameters, Ko and L, in the original Mualem model to obtain a diffusion coefficient-modified Mualem model.


## 2.2. Determination of the Diffusion Coefficient by Horizontal 1D Seepage Test [13]

Without considering the action of gravity or the deformation of soil mass, the differential equation and boundary conditions for moisture movement in the horizontal 1D seepage test without free head are shown as follows:(3)∂θ∂t=∂∂Dθ∂θ∂x,(4)θ=θax>0,t=0,(5)θ=θox=0,t>0,where θ is the volumetric moisture content (cm3·cm−3), Dθ is the diffusion coefficient (cm2·min−1), x is the distance from the water inlet (cm), θa and θo are, respectively, the initial volumetric moisture content of soil sample and the volumetric moisture content at water inlet (cm3·cm−3), and t is the infiltration time (min).After measuring the distribution of the volumetric moisture contentθ of soil columns at time t, λθ~θ can be calculated according to the Boltzmann parameter as λθ=xt-0.5. The diffusion coefficient Dθ can be calculated as follows: (6)Dθ=-12ΔλΔθ∑θaθλdθ.The relationship among the diffusion coefficientDθ, the SWCC ψθ, and the infiltration coefficient Kθ is given as(7)Kθ=-Dθ·dθdψθ.After measuring the diffusion coefficient and the SWCC, the unsaturated infiltration coefficient can be calculated from (7).


## 2.3. SWCC Model

In order to modify the Mualem model, the first step is to determine the SWCC model. For this purpose, the present paper adopted the Brook and Corey model, as shown below:(8)θ-θrθs-θr=Se=1,ψ≤ψa,(9)θ-θrθs-θr=Se=ψaψσ,ψ>ψa,where σ is the shape parameter of the SWCC, θs is the saturated volumetric moisture content, θr is the residual volumetric moisture content, and ψa is the air entry value.


## 3. Issues Related to the Modification of Mualem Model

### 3.1. Calculation Method of Modified Parameters

The unsaturated infiltration coefficient can be calculated after measuring the SWCC and the diffusion coefficient in the test. The Mualem model was modified directly based on the measured diffusion coefficient in this paper to reduce the steps of calculating the unsaturated infiltration coefficient and to avoid the accumulation of errors in the calculation process.A relational expression for the diffusion coefficient and the volumetric moisture content can be obtained, through transformation, by substituting (7) and (9) into (5):(10)Dθ=ψaσθs-θrKoθ-θrθs-θr1+L+1/σ.After obtaining the parameters of (10) through curve fitting the diffusion coefficient and volumetric moisture content measured in the horizontal 1D seepage test, we can substitute the SWCC parameters to calculate the modified parameters Ko and L of the diffusion coefficient-modified Mualem model.


### 3.2. Requirements on Horizontal 1D Seepage Test

The diffusion coefficient-modified Mualem model, as discussed in Section3.1, actually combines the two methods of indirectly determining the unsaturated infiltration coefficient. As indicated by (10), the diffusion coefficient of soil sample was zero when the volumetric moisture content was equal to or lower than the residual moisture content. In the horizontal 1D seepage test, the diffusion coefficient corresponding to the initial moisture content of soil sample was zero, so the combination of these two methods must satisfy(11)θr=θo.The residual moisture contentθr is an important parameter of the SWCC, so, in the horizontal 1D seepage test, the initial moisture content of soil sample should be precontrolled to be approximately equal to the residual moisture content.


## 3.1. Calculation Method of Modified Parameters

The unsaturated infiltration coefficient can be calculated after measuring the SWCC and the diffusion coefficient in the test. The Mualem model was modified directly based on the measured diffusion coefficient in this paper to reduce the steps of calculating the unsaturated infiltration coefficient and to avoid the accumulation of errors in the calculation process.A relational expression for the diffusion coefficient and the volumetric moisture content can be obtained, through transformation, by substituting (7) and (9) into (5):(10)Dθ=ψaσθs-θrKoθ-θrθs-θr1+L+1/σ.After obtaining the parameters of (10) through curve fitting the diffusion coefficient and volumetric moisture content measured in the horizontal 1D seepage test, we can substitute the SWCC parameters to calculate the modified parameters Ko and L of the diffusion coefficient-modified Mualem model.


## 3.2. Requirements on Horizontal 1D Seepage Test

The diffusion coefficient-modified Mualem model, as discussed in Section3.1, actually combines the two methods of indirectly determining the unsaturated infiltration coefficient. As indicated by (10), the diffusion coefficient of soil sample was zero when the volumetric moisture content was equal to or lower than the residual moisture content. In the horizontal 1D seepage test, the diffusion coefficient corresponding to the initial moisture content of soil sample was zero, so the combination of these two methods must satisfy(11)θr=θo.The residual moisture contentθr is an important parameter of the SWCC, so, in the horizontal 1D seepage test, the initial moisture content of soil sample should be precontrolled to be approximately equal to the residual moisture content.


## 4. Modification of the Mualem Model

### 4.1. Test Introduction

The loess used in the test was collected from Heping Town in Ningxia Province at a depth of 2-3 m. The horizontal 1D seepage test was carried out in an organic glass tube 186 mm in inner diameter and 1,000 mm in length, and the soil sample was placed in the organic glass tube in layers of 5 cm thickness and was compacted until the set dry density was achieved. A Mariotte bottle was used to supply water for the seepage test, and the free head was controlled at 10 cm. According to the literature [16], the SWCC was determined in the horizontal 1D seepage process with the preembedded TDR moisture meter and Fredlund thermal conduction suction probe jointly. After 10–12 hours of water supply and seepage, the soil sample was sliced and dried to calculate the diffusion coefficient.The loess used in the test was naturally air-dried loess, and its physical parameters are shown in Table1. According to the residual moisture content of remolded loess in [9] and the method of measuring the residual moisture content in [17], it was confirmed that the initial moisture content of soil sample had satisfied relevant requirements. Figure 1 shows the relationship between the diffusion coefficient and saturation S and Figure 2 shows the measured SWCC.Table 1
Physical properties of remolded loess.
Particulate density/g⋅cm−3 Mass moisture content/% W P/% W L/% I P 2.71 6.3 17.6 28.7 10.1Figure 1
Diffusion coefficient under different dry densities.Figure 2
SWCC under different dry densities.The measured SWCC was fitted according to the Brook and Corey model. The resulting parameters are shown in Table2.Table 2
SWCC fitting parameters.
Dry density/g⋅cm−3 θ s/cm3⋅cm−3 θ r/cm3⋅cm−3 ψ a/cm 1/σ R 2 1.25 0.539 0.079 169 0.90 0.956 1.35 0.502 0.083 227.91 0.94 0.978 1.45 0.465 0.091 290 1.01 0.972 1.55 0.428 0.098 373.79 0.98 0.972 1.65 0.391 0.104 487.7 0.91 0.981According to Table2, the correlation coefficient R2 between the Brook and Corey model data and the measured SWCC data exceeded 0.95, confirming that the Brook and Corey model applies to the remolded loess used in the test and suggesting that the initial moisture content had been appropriately controlled.


### 4.2. Applicability Test of Mualem Model and Variable Parameter Modification Method

Before modification, the applicability of the Mualem model and the variable parameter modification method was tested, thus to demonstrate the necessity of modification. The measured diffusion coefficient was fitted by (10), as shown in Figure 3.Figure 3
Curve fitting for the diffusion coefficient.The modified parametersKo and L of (7) were then calculated on the basis of Table 3.Table 3
Modified parametersKo and L.
Dry density/g⋅cm−3 K o/cm2⋅min−1 L R 2 1.25 1.95 E - 02 1.39 0.726 1.35 1.02 E - 02 1.63 0.790 1.45 9.76 E - 03 2.46 0.852 1.55 2.83 E - 03 2.66 0.902 1.65 1.82 E - 03 3.53 0.908Figure3 and Table 3 indicated that (i) the goodness of fit between (10) and the measured diffusion coefficient increased with an increase in dry density, (ii) the modified parameter Ko diminished with an increase in dry density, and (iii) the modification coefficient L increased with an increase in dry density. Two conclusions can be drawn. First, the Mualem model could not be directly applied to the remolded loess used in the test, as (2) had set β and γ in correspondence with the Mualem model and the L value obtained through curve fitting ranged between 1.39 and 3.53, which is significantly different from the L value of 0.5 in the Mualem model. Second, the method of purely introducing the variable parameters Ko and L to modify the Mualem model also needed to be improved, as the correlation between (10) and the measured diffusion coefficient was not high (the correlation coefficient ranging only between 0.726 and 0.908).


### 4.3. Segmental Modification of the Mualem Model

The low goodness of fit between the method of modifying the Mualem model with variable parameters and the measured value shown in Figure1 is due to the abrupt change in the diffusion coefficient of soil around a degree of saturation of 0.65, which has also been demonstrated by Chen et al. [18] through experimental tests.As can be found from the analysis of Figure1, the diffusion coefficient of the two individual segments before and after the degree of saturation of 0.65 and (10) had similar curve shapes. Thus, the measured diffusion coefficient was split into two segments separated at the degree of saturation of 0.65, and the Mualem model was segmentally modified by (10). The resulting curve fitting parameters are shown in Figure 4.Figure 4
Piecewise fitting of the diffusion coefficient.As shown in Figures3 and 4, the segmental modification method combined with the variable parameter method significantly increased the goodness of fit between the Mualem model data and the measured values. The modified parameters Ko and L of the modified Mualem models corresponding to various dry densities of the two segments were calculated based on Table 3 and are shown in Table 4.Table 4
Modified parameters of piecewise fitting.
Dry density/g⋅cm−3 Segment less than 0.65 in saturation Segment above 0.65 in saturation K o/cm2⋅min−1 L R 2 K o/cm2⋅min−1 L R 2 1.25 4.34 E - 03 0.39 0.921 1.31 10.83 0.998 1.35 2.51 E - 03 0.72 0.940 0.41 10.51 0.995 1.45 1.32 E - 03 0.96 0.972 7.27 E - 02 8.22 0.990 1.55 8.65 E - 04 1.64 0.973 1.93 E - 02 7.29 0.994 1.65 3.21 E - 04 2.00 0.962 7.68 E - 03 7.62 0.995As shown in Table4, Ko always diminished with an increase in dry density. However, L showed opposing trends in the two saturation segments for variation with an increase in dry density. In the curve fitting, the correlation coefficient ranged from 0.921 to 0.998. The results obtained by introducing the variable parameters Ko and L and adopting the saturation of 0.65 as the boundary for the segmental modification of the Mualem model fitted well with the test results.There are three types of unsaturated soil corresponding to three types of water connection states and transmission pressure capacities. Thus, the diffusion coefficient may also present a three-segment variation law. In the present paper, the initial moisture content of the horizontal 1D seepage test sample was controlled to be approximately equal to the residual moisture content. Moreover, the diffusion coefficients measured were the diffusion coefficients of the double-open unsaturated soil and gas-sealed unsaturated soil. Thus, good results were achieved by dividing the diffusion coefficient into two segments.


## 4.1. Test Introduction

The loess used in the test was collected from Heping Town in Ningxia Province at a depth of 2-3 m. The horizontal 1D seepage test was carried out in an organic glass tube 186 mm in inner diameter and 1,000 mm in length, and the soil sample was placed in the organic glass tube in layers of 5 cm thickness and was compacted until the set dry density was achieved. A Mariotte bottle was used to supply water for the seepage test, and the free head was controlled at 10 cm. According to the literature [16], the SWCC was determined in the horizontal 1D seepage process with the preembedded TDR moisture meter and Fredlund thermal conduction suction probe jointly. After 10–12 hours of water supply and seepage, the soil sample was sliced and dried to calculate the diffusion coefficient.The loess used in the test was naturally air-dried loess, and its physical parameters are shown in Table1. According to the residual moisture content of remolded loess in [9] and the method of measuring the residual moisture content in [17], it was confirmed that the initial moisture content of soil sample had satisfied relevant requirements. Figure 1 shows the relationship between the diffusion coefficient and saturation S and Figure 2 shows the measured SWCC.Table 1
Physical properties of remolded loess.
Particulate density/g⋅cm−3 Mass moisture content/% W P/% W L/% I P 2.71 6.3 17.6 28.7 10.1Figure 1
Diffusion coefficient under different dry densities.Figure 2
SWCC under different dry densities.The measured SWCC was fitted according to the Brook and Corey model. The resulting parameters are shown in Table2.Table 2
SWCC fitting parameters.
Dry density/g⋅cm−3 θ s/cm3⋅cm−3 θ r/cm3⋅cm−3 ψ a/cm 1/σ R 2 1.25 0.539 0.079 169 0.90 0.956 1.35 0.502 0.083 227.91 0.94 0.978 1.45 0.465 0.091 290 1.01 0.972 1.55 0.428 0.098 373.79 0.98 0.972 1.65 0.391 0.104 487.7 0.91 0.981According to Table2, the correlation coefficient R2 between the Brook and Corey model data and the measured SWCC data exceeded 0.95, confirming that the Brook and Corey model applies to the remolded loess used in the test and suggesting that the initial moisture content had been appropriately controlled.


## 4.2. Applicability Test of Mualem Model and Variable Parameter Modification Method

Before modification, the applicability of the Mualem model and the variable parameter modification method was tested, thus to demonstrate the necessity of modification. The measured diffusion coefficient was fitted by (10), as shown in Figure 3.Figure 3
Curve fitting for the diffusion coefficient.The modified parametersKo and L of (7) were then calculated on the basis of Table 3.Table 3
Modified parametersKo and L.
Dry density/g⋅cm−3 K o/cm2⋅min−1 L R 2 1.25 1.95 E - 02 1.39 0.726 1.35 1.02 E - 02 1.63 0.790 1.45 9.76 E - 03 2.46 0.852 1.55 2.83 E - 03 2.66 0.902 1.65 1.82 E - 03 3.53 0.908Figure3 and Table 3 indicated that (i) the goodness of fit between (10) and the measured diffusion coefficient increased with an increase in dry density, (ii) the modified parameter Ko diminished with an increase in dry density, and (iii) the modification coefficient L increased with an increase in dry density. Two conclusions can be drawn. First, the Mualem model could not be directly applied to the remolded loess used in the test, as (2) had set β and γ in correspondence with the Mualem model and the L value obtained through curve fitting ranged between 1.39 and 3.53, which is significantly different from the L value of 0.5 in the Mualem model. Second, the method of purely introducing the variable parameters Ko and L to modify the Mualem model also needed to be improved, as the correlation between (10) and the measured diffusion coefficient was not high (the correlation coefficient ranging only between 0.726 and 0.908).


## 4.3. Segmental Modification of the Mualem Model

The low goodness of fit between the method of modifying the Mualem model with variable parameters and the measured value shown in Figure1 is due to the abrupt change in the diffusion coefficient of soil around a degree of saturation of 0.65, which has also been demonstrated by Chen et al. [18] through experimental tests.As can be found from the analysis of Figure1, the diffusion coefficient of the two individual segments before and after the degree of saturation of 0.65 and (10) had similar curve shapes. Thus, the measured diffusion coefficient was split into two segments separated at the degree of saturation of 0.65, and the Mualem model was segmentally modified by (10). The resulting curve fitting parameters are shown in Figure 4.Figure 4
Piecewise fitting of the diffusion coefficient.As shown in Figures3 and 4, the segmental modification method combined with the variable parameter method significantly increased the goodness of fit between the Mualem model data and the measured values. The modified parameters Ko and L of the modified Mualem models corresponding to various dry densities of the two segments were calculated based on Table 3 and are shown in Table 4.Table 4
Modified parameters of piecewise fitting.
Dry density/g⋅cm−3 Segment less than 0.65 in saturation Segment above 0.65 in saturation K o/cm2⋅min−1 L R 2 K o/cm2⋅min−1 L R 2 1.25 4.34 E - 03 0.39 0.921 1.31 10.83 0.998 1.35 2.51 E - 03 0.72 0.940 0.41 10.51 0.995 1.45 1.32 E - 03 0.96 0.972 7.27 E - 02 8.22 0.990 1.55 8.65 E - 04 1.64 0.973 1.93 E - 02 7.29 0.994 1.65 3.21 E - 04 2.00 0.962 7.68 E - 03 7.62 0.995As shown in Table4, Ko always diminished with an increase in dry density. However, L showed opposing trends in the two saturation segments for variation with an increase in dry density. In the curve fitting, the correlation coefficient ranged from 0.921 to 0.998. The results obtained by introducing the variable parameters Ko and L and adopting the saturation of 0.65 as the boundary for the segmental modification of the Mualem model fitted well with the test results.There are three types of unsaturated soil corresponding to three types of water connection states and transmission pressure capacities. Thus, the diffusion coefficient may also present a three-segment variation law. In the present paper, the initial moisture content of the horizontal 1D seepage test sample was controlled to be approximately equal to the residual moisture content. Moreover, the diffusion coefficients measured were the diffusion coefficients of the double-open unsaturated soil and gas-sealed unsaturated soil. Thus, good results were achieved by dividing the diffusion coefficient into two segments.


## 5. The Mualem Model Modified by the Unsaturated Infiltration Coefficient of Loess

### 5.1. Combination of the Prediction Method and the Calculation Method for Determining the Unsaturated Infiltration Coefficient

Among the unsaturated infiltration coefficient, the SWCC, and the diffusion coefficient, there are two independent parameters. The first method, the prediction method, based on the assumption of the continuous distribution of soil pores, directly establishes the relationship between the SWCC and the unsaturated infiltration coefficient. The second method, the calculation method, solves the problem on the basis of the relationships among the three. The above two methods of determining the unsaturated infiltration coefficient both require the measurement of the SWCC. When the prediction method is combined with the calculation method and the parametersKo and L in (10) are given, the expressions of the unsaturated infiltration coefficient and the SWCC can be simultaneously obtained from the measured diffusion coefficient. When compared to the SWCC, the diffusion coefficient can be more easily determined.The combination of the two methods simplifies the determination of the unsaturated infiltration coefficient and the SWCC. The key lies in selecting reasonable prediction method for unsaturated infiltration coefficient. Modifying a prediction model for the unsaturated infiltration coefficient can not only increase the accuracy in predicting the unsaturated infiltration coefficient, but also lower the requirements imposed by the unsaturated soil seepage parameter test on test instruments and equipment.


### 5.2. Dry Density-Based Modified Mualem Model

It is found in the analysis thatKo and dry density show a power exponential relationship, while L and dry density show a linear relationship. For reference, the fitting charts are shown in Figures 5 and 6, and the fitting equations are shown in Table 5.Table 5
Expressions between the modified parameters and dry density.
Modified parameter Relational expression with dry density R 2 Segment when saturation is less than 0.65 K o = 0.04 × ρ d - 9 0.974 L = 5.31 × ρ d - 5.37 0.959 Segment when saturation is greater than 0.65 K o = 103.87 × ρ d - 19.24 0.994 L = - 9.64 × ρ d + 22.87 0.844Figure 5
Relationship between the modified parameterKo and dry density.Figure 6
Relationship between the modified parameterL and dry density.By substituting the fitting equations in Table5 into (1), the segmentally modified Mualem model can be obtained for the unsaturated infiltration coefficient of remolded loess. The model considers the variation in the unsaturated infiltration coefficient with the dry density of remolded loess and can be expressed as(12)Kθ=aρdbSecρd-d∫0Seψ-1dSe∫01ψ-1dSe2,where parameters a, b, c, and d can be found in Table 5.


## 5.1. Combination of the Prediction Method and the Calculation Method for Determining the Unsaturated Infiltration Coefficient

Among the unsaturated infiltration coefficient, the SWCC, and the diffusion coefficient, there are two independent parameters. The first method, the prediction method, based on the assumption of the continuous distribution of soil pores, directly establishes the relationship between the SWCC and the unsaturated infiltration coefficient. The second method, the calculation method, solves the problem on the basis of the relationships among the three. The above two methods of determining the unsaturated infiltration coefficient both require the measurement of the SWCC. When the prediction method is combined with the calculation method and the parametersKo and L in (10) are given, the expressions of the unsaturated infiltration coefficient and the SWCC can be simultaneously obtained from the measured diffusion coefficient. When compared to the SWCC, the diffusion coefficient can be more easily determined.The combination of the two methods simplifies the determination of the unsaturated infiltration coefficient and the SWCC. The key lies in selecting reasonable prediction method for unsaturated infiltration coefficient. Modifying a prediction model for the unsaturated infiltration coefficient can not only increase the accuracy in predicting the unsaturated infiltration coefficient, but also lower the requirements imposed by the unsaturated soil seepage parameter test on test instruments and equipment.


## 5.2. Dry Density-Based Modified Mualem Model

It is found in the analysis thatKo and dry density show a power exponential relationship, while L and dry density show a linear relationship. For reference, the fitting charts are shown in Figures 5 and 6, and the fitting equations are shown in Table 5.Table 5
Expressions between the modified parameters and dry density.
Modified parameter Relational expression with dry density R 2 Segment when saturation is less than 0.65 K o = 0.04 × ρ d - 9 0.974 L = 5.31 × ρ d - 5.37 0.959 Segment when saturation is greater than 0.65 K o = 103.87 × ρ d - 19.24 0.994 L = - 9.64 × ρ d + 22.87 0.844Figure 5
Relationship between the modified parameterKo and dry density.Figure 6
Relationship between the modified parameterL and dry density.By substituting the fitting equations in Table5 into (1), the segmentally modified Mualem model can be obtained for the unsaturated infiltration coefficient of remolded loess. The model considers the variation in the unsaturated infiltration coefficient with the dry density of remolded loess and can be expressed as(12)Kθ=aρdbSecρd-d∫0Seψ-1dSe∫01ψ-1dSe2,where parameters a, b, c, and d can be found in Table 5.


## 6. Conclusions

Taking advantage of its simple form, the Mualem model has been widely applied to the prediction of the unsaturated infiltration coefficient. In this paper, an applicability test and a modification in the Mualem model were performed through an experiment. The main conclusions are as follows:(1)
The results of the applicability analysis of the Mualem model indicated that the Mualem model could not be satisfactorily applied to the remolded loess tested in the present paper. The method of introducing only two variable parametersKo and L for modifying the Mualem model had a low correlation coefficient with the measured value.(2)
An abrupt change in the diffusion coefficient was observed corresponding to a saturation of 0.65. This was attributed to the change in water connection state in unsaturated soil.(3)
When the saturation of 0.65 was adopted as the boundary to divide the curve of measured diffusion coefficient into two segments, the fitting correlation coefficient ranged from 0.921 to 0.998. The approach of combining the segmental modification method and the variable parameter method fitted well with the test results. This paper obtained the values of modified parametersKo and L corresponding to remolded loess of five dry densities.(4)
Based on the power exponential function betweenKo and dry density and the linear relationship between L and dry density, a segmentally modified Mualem model was established by considering dry density changes of remolded loess.



---

*Source: 2768952-2017-01-15.xml*