Journal cover Journal topic
Hydrology and Earth System Sciences An interactive open-access journal of the European Geosciences Union
Journal topic
Hydrol. Earth Syst. Sci., 22, 4621-4632, 2018
https://doi.org/10.5194/hess-22-4621-2018
Hydrol. Earth Syst. Sci., 22, 4621-4632, 2018
https://doi.org/10.5194/hess-22-4621-2018

Research article 04 Sep 2018

Research article | 04 Sep 2018

# Predicting the soil water retention curve from the particle size distribution based on a pore space geometry containing slit-shaped spaces

Predicting the soil water retention curve
Chen-Chao Chang1,2 and Dong-Hui Cheng1,2 Chen-Chao Chang and Dong-Hui Cheng
• 1School of Environmental Sciences and Engineering, Chang'an University, Xi'an, 710054, China
• 2Key Laboratory of Subsurface Hydrology and Ecological Effects in Arid Region, Chang'an University, Ministry of Education, Xi'an, China
Abstract

Traditional models employed to predict the soil water retention curve (SWRC) from the particle size distribution (PSD) always underestimate the water content in the dry range of the SWRC. Using the measured physical parameters of 48 soil samples from the UNSODA unsaturated soil hydraulic property database, these errors were proven to originate from an inaccurate estimation of the pore size distribution. A method was therefore proposed to improve the estimation of the water content at high suction heads using a pore model comprising a circle-shaped central pore connected to slit-shaped spaces. In this model, the pore volume fraction of the minimum pore diameter range and the corresponding water content were accordingly increased. The predicted SWRCs using the improved method reasonably approximated the measured SWRCs, which were more accurate than those obtained using the traditional method and the scaling approach in the dry range of the SWRC.

1 Introduction

The soil water retention curve (SWRC), which represents the relationship between the water pressure and water content, is fundamental to researching water flow and chemical transport in unsaturated media (Pollacco et al., 2017). Direct measurements of the SWRC consume both time and money (Arya and Paris, 1981; Mohammadi and Vanclooster, 2011), while estimating the SWRC from the particle size distribution (PSD) is both rapid and economical. Therefore, a number of associated conceptual and physical models have been proposed.

The first attempt to directly translate a PSD into an SWRC was made by Arya and Paris (1981) (hereafter referred to as the AP model). In this model, the PSD is divided into multiple size fractions and the bulk and particle densities of the natural-structure samples are uniformly applied to each particle size fraction, from which it follows that the pore fraction and the corresponding solid fraction are equal. Thus, the degree of saturation can be set equal to the cumulative PSD function. The soil suction head can be obtained using the capillary equation based on a “bundle of cylindrical tubes” model, and the pore size in the equation is determined by scaling the pore length and pore volume (Arya et al., 2008). Based on the principle of the AP model, many researchers have focused on improving the suction head calculations, which are commonly based on the capillary equation; but methods that are used to translate the particle diameter into the pore diameter are different (Haverkamp and Parlange, 1986; Zhuang et al., 2001; Mohammadi and Vanclooster, 2011; Jensen et al., 2015). Some models estimate the pore diameter based on the particle packing patterns (e.g., the MV model) (Meskini-Vishkaee et al., 2014), while others utilize the proportionality factor between the pore size and the associated particle diameter (e.g., the HP model and the two-stage approach) (Haverkamp and Parlange, 1986; Jensen et al., 2015). However, the scheme employed to estimate the water content has not been modified and follows the approach of the AP model. The SWRC prediction models which use the same scheme to predict the water content and only improve the suction head calculation are referred to as the traditional models in the following text.

Many traditional models are based on a “bundle of cylindrical tubes” representation of the pore space geometry (Arya and Paris, 1981; Zhuang et al., 2001), which results in intrinsic errors when predicting the water flow in variably saturated soils. Consequently, some researchers have considered pore networks as bundles of triangular tubes, which could incorporate the contribution of water in pore corners to the water content (Helland and Skjæveland, 2007). A new pore geometry model comprised of a polygon-shaped central pore connected to slit-shaped spaces was proposed by Tuller et al. (1999) to provide a more realistic representation of natural pore spaces (Tuller et al., 1999; Or and Tuller, 1999; Tuller and Or, 2001). This pore model could represent a foundation for accurately describing the water status in natural soils, particularly in arid environments.

Therefore, the objectives of this study were to evaluate the leading factors that lead to an underestimation of the water content in the dry range of the predicted SWRC using traditional methods and to furthermore propose a method for accurately estimating the water content using a pore space geometry containing slit-shaped spaces to improve the prediction of the SWRC.

2 Basic descriptions

The relationship between the PSD and the PoSD is a fundamental element when predicting the SWRC from the PSD. Hwang and Powers (2003) found that the nonlinear relationship between the PSD and the PoSD is more appropriate than the linear relationship applied in the AP model and therefore described both the PSD and the PoSD as lognormal distributions. However, since the PSD and PoSD of soils do not strongly follow a lognormal distribution, this model performed very poorly for moderately fine-textured soils (Hwang and Choi, 2006). Obtaining an accurate PoSD from the PSD of a soil is highly difficult, and the errors that arise from this approach could cause inevitable errors in the predicted SWRC. However, the underestimation of the water content in the dry range of an SWRC has not been comprehensively evaluated from this perspective.

In this study, the measured PoSDs of 48 soil samples were compared with the PoSDs calculated using a traditional model (they were actually the corresponding PSDs) to identify the origins of the errors and their effects on the accuracy of the predicted SWRC. The provided 48 soil samples exhibited a wide range of physical properties (Table 1) and were selected from the UNSODA unsaturated soil hydraulic property database, which contains 790 soil samples with general unsaturated soil hydraulic properties and basic soil properties (e.g., water retention, hydraulic conductivity, soil water diffusivity, PSD, bulk density, and organic matter content) (Nemes et al., 2001). The maximum, minimum and mean values of the soil bulk density and the percentages of clay and sand of the used soil samples for the calibration stage are presented in Table 2.

Table 1Codes and texture classes of the 48 soils selected from UNSODA.

Table 2Basic soil properties of 48 samples for the model calibration.

Figure 1Measured vs. calculated pore volume fraction curves using the traditional method for (a) sand (code: 3172) and (b) clay (code: 2360). The measured and calculated PoSDs are in the insets at the top of the figures.

1. Calculating the PoSD using a traditional model

Traditional models commonly assume that the pore volume fraction of each size fraction can be set equal to the corresponding solid fraction (Arya and Paris, 1981). Thus, the cumulative pore volume fraction can take the following form:

$\begin{array}{}\text{(1)}& \sum _{j=\mathrm{1}}^{j=i}{\mathit{\nu }}_{j}=\sum _{j=\mathrm{1}}^{j=i}{\mathit{\omega }}_{j};\phantom{\rule{0.125em}{0ex}}\phantom{\rule{0.125em}{0ex}}\phantom{\rule{0.125em}{0ex}}i=\mathrm{1},\mathrm{2},\phantom{\rule{0.25em}{0ex}}\mathrm{\dots }\phantom{\rule{0.25em}{0ex}}n,\end{array}$

where ωj is the solid fraction of the jth particle fraction, νj is the pore volume fraction associated with the jth fraction, and n is the total number of size fractions in the PSD.

The routine procedures employed among the several traditional models to translate a particle diameter into a pore diameter are different. The equivalent pore diameter can be derived from physical properties, including the bulk density and the particle density, or from the proportionate relationship between the pore size and associated particle diameter. Although the former can logically characterize a pore, a complicated pattern can slightly reduce the model performance, while the latter approach is easy to use, and its rationality has been demonstrated by some researchers (Hamamoto et al., 2011; Sakaki et al., 2014). Here, the latter technique is applied, and it can be expressed as

$\begin{array}{}\text{(2)}& {d}_{i}=\mathrm{0.3}{D}_{i},\end{array}$

where Di is the mean particle diameter of the ith fraction (µm) and di is the corresponding equivalent pore diameter (µm). Inputting the PSD data, the calculated pore diameters are sequentially paired with corresponding pore volume fractions to obtain a calculated PoSD.

2. Estimating the PoSD from the measured SWRC

It is generally difficult to measure the PoSD of a soil; however, the PoSD can be indirectly obtained using the measured water content and suction head (Jayakody et al., 2014). The cumulative pore volume fraction of the ith fraction is equal to the ratio of the measured water content to the saturated water content (Eq. 3):

$\begin{array}{}\text{(3)}& {\mathit{\theta }}_{i}={\mathit{\theta }}_{s}\sum _{j=\mathrm{1}}^{j=i}{\mathit{\nu }}_{j};\phantom{\rule{0.125em}{0ex}}\phantom{\rule{0.125em}{0ex}}\phantom{\rule{0.125em}{0ex}}i=\mathrm{1},\mathrm{2},\mathrm{\dots }n,\end{array}$

where θs is the saturated water content (cm3 cm−3) and θi is the measured water content (cm3 cm−3).

Meanwhile, the corresponding pore diameters are derived on the basis of Laplace's equation and Eq. (4):

$\begin{array}{}\text{(4)}& {\mathit{\psi }}_{i}=\frac{\mathrm{2}\mathit{\sigma }\mathrm{cos}\mathit{\epsilon }}{{r}_{i}g{\mathit{\rho }}_{\mathrm{w}}},\end{array}$

where ψi is the suction head (mH2O), σis the surface tension (kg s−2), ε is the contact angle between the soil particle and water, ri is the pore radius (m), and ρw is the density of water (kg m−3). Assuming for water at 20 C $\mathit{\sigma }=\mathrm{7.275}×{\mathrm{10}}^{-\mathrm{2}}$ kg s−2, ρw=998.9 kg m−3, g=9.81 m s−2, and $\mathit{\epsilon }=\mathrm{0}{}^{\circ }$ (Mohammadi and Vanclooster, 2011), then transforming ri to di and substituting numerically the values of the constant yield a simplified expression as Eq. (5):

$\begin{array}{}\text{(5)}& {\mathit{\psi }}_{i}=\frac{\mathrm{3000}}{{d}_{i}},\end{array}$

where ψi is the suction head (cmH2O) and di is the pore diameter (µm). Then, the pore diameters calculated by Eq. (5) were sequentially paired with the cumulative pore volume fractions calculated by Eq. (3) to obtain a PoSD, which could be considered a measured PoSD.

The calculated and measured PoSD data were fitted using a modified logistic growth model (Eq. 6) (Liu et al., 2003):

$\begin{array}{}\text{(6)}& {w}_{i}=\frac{\mathrm{1}}{\mathrm{1}+a\mathrm{exp}\left(-b{d}_{i}^{c}\right)},\end{array}$

where wi is the cumulative pore volume fraction with diameters smaller than di (%) and a, b, and c are the fitting parameters (dimensionless). This model produced a good fit for the PoSD data employed in this study with a coefficient of determination (r2) that ranged from 0.972 to 0.999.

Figure 2Pore model containing two slit-shaped spaces (d denotes the diameter of the central pore, and αd and βd denote the widths and lengths of the slit-shaped spaces, respectively).

The measured pore volume fraction curves for the typical samples, namely, sand (code: 3172) and clay (code: 2360), and their calculated curves using the traditional model are presented in Fig. 1. The small maps embedded in Fig. 1 exhibit the measured and calculated PoSD curves, which show that the calculated PoSD curves approximately coincide with the measured curves in the larger pore diameter range, while in the smaller range, which corresponds to the higher suction head range on the SWRC, the calculated values are obviously smaller than the measured values. The underestimation of the pore volume fraction in the smaller pore diameter range can consequently lead to an underestimation of the water content at high suction heads. In particular, the calculated pore volume fraction associated with the smallest pore diameter (d≤0.6µm) was far less than the measured value. These results illustrated that the underestimation of the pore volume fraction with respect to the smallest pore diameter (d≤0.6µm) was a key factor with regard to the underestimation of the water content in the dry range of the SWRC. In addition, the underestimation of the pore volume fraction is associated with an oversimplified pore space geometry, which traditional models have generally characterized as a bundle of cylindrical capillaries. The measured and calculated pore curves of the other 46 soil samples behaved in the same fashion, and those curves are provided in the Supplement (Fig. S1).

3 Improved method

## 3.1 Estimating the pore volume fraction

In this study, the soil pore structure was conceptualized within a pore model in which the elementary unit cell is composed of a relatively larger circle-shaped central pore connected to two slit-shaped spaces (see Fig. 2). Relative to the polygonal central pore connected to the slit-shaped spaces as described by Or and Tuller (1999), both the slit width and the slit length are proportional to the diameter of the associated central pore d and are therefore expressed as αd and βd, respectively.

When estimating the pore volume fraction using the pore model described above, the volume fractions of the central pore and slit-shaped spaces are distinguished. Considering that the sizes of the slit-shaped spaces are smaller than that of the minimum central pore, the slit-shaped spaces are accordingly classified into it. Therefore, the pore volume fractions of the soil samples were simplified into those of the central pores, but the volume fraction of the minimum central pores included that of all slit-shaped spaces. Using the geometric relationship described in Fig. 2 and the traditional assumption that the volume fraction of each unit cell (i.e., the central pore connecting to two slit-shaped spaces) is equal to the corresponding particle mass fraction, the pore volume fractions with respect to different sizes can be readily obtained.

The procedure utilized to calculate the pore volume fractions is shown in Fig. 3. Assuming that the soil pores are composed of numerous unit cells with various sizes, the fraction of the ith unit cell is equal to the relative particle mass fraction ωi. The addition of the volume fraction of the smallest unit cell (ω1) and the sum of the slit volume fractions of various sizes (${\mathit{\zeta }}_{\mathrm{2}}+{\mathit{\zeta }}_{\mathrm{3}}+\mathrm{\dots }+{\mathit{\zeta }}_{i}$) result in the volume fraction of the smallest pore (ν1). Successively accumulating it with the volume fractions of other central pore (i.e., ν2, ν3, ν4 …) provides the PoSD of a sample. The volume fraction of the slit-shaped spaces, ζi, the volume fraction of the smallest pore ν1 and the volume fractions of the other pores νi were calculated using Eqs. (7), (8) and (9), respectively:

$\begin{array}{}\text{(7)}& & {\mathit{\zeta }}_{i}={\mathit{\omega }}_{i}\frac{\mathrm{2}\mathit{\alpha }\mathit{\beta }{d}_{i}^{\mathrm{2}}}{\mathrm{2}\mathit{\alpha }\mathit{\beta }{d}_{i}^{\mathrm{2}}+\frac{\mathit{\pi }}{\mathrm{4}}{d}_{i}^{\mathrm{2}}},\text{(8)}& & {\mathit{\nu }}_{\mathrm{1}}={\mathit{\omega }}_{\mathrm{1}}+\sum _{\mathrm{2}}^{n}{\mathit{\zeta }}_{i},\text{(9)}& & {\mathit{\nu }}_{i}={\mathit{\omega }}_{i}-{\mathit{\zeta }}_{i},\end{array}$

where ζi is the volume fraction of the slit-shaped spaces, νi is the volume fraction of the ith pore fraction, and α and β are the scaling parameters of the slit width and the slit length, respectively.

Figure 3Schematic of the procedure used to calculate the pore volume fraction.

## 3.2 Values of α and β

To obtain the values of α and β, an expression containing both of these parameters with respect to the specific surface area (SSA) was applied here. The SSA of the pore as shown in Fig. 2 can be described using a geometrical relationship as follows:

$\begin{array}{}\text{(10)}& {S}_{\mathrm{SA}}=\frac{\mathit{\varphi }}{\mathrm{1000}{\mathit{\rho }}_{b}}\sum _{i=\mathrm{1}}^{n}{\mathit{\omega }}_{i}\frac{\mathrm{4}\mathit{\beta }{d}_{i}+\mathit{\pi }{d}_{i}}{\mathrm{2}\mathit{\alpha }\mathit{\beta }{d}_{i}^{\mathrm{2}}+\frac{\mathit{\pi }}{\mathrm{4}}{d}_{i}^{\mathrm{2}}};\phantom{\rule{0.125em}{0ex}}\phantom{\rule{0.125em}{0ex}}\phantom{\rule{0.125em}{0ex}}i=\mathrm{1},\mathrm{2},\mathrm{\dots }n,\end{array}$

where SSA is the specific surface area (m2 g−1), di is the pore diameter (m), ρb is the bulk density (kg m−3) and Φ is the measured porosity. Therefore, an important requirement for the calculation of the α and β values is an estimation of the SSA at sample scale. Here, a power equation was applied as follows (Sepaskhah et al., 2010):

$\begin{array}{}\text{(11)}& {S}_{\mathrm{SA}}=\mathrm{3.89}{d}_{\mathrm{g}}^{-\mathrm{0.905}},\end{array}$

where SSA is the estimated specific surface area (m2 g−1), and dg is the geometric mean particle size diameter (mm) obtained using Eq. (12) (Shirazi and Boersma, 1984):

$\begin{array}{}\text{(12)}& {d}_{\mathrm{g}}=\mathrm{exp}\left({f}_{\mathrm{c}}\mathrm{ln}{M}_{\mathrm{c}}+{f}_{\mathrm{si}}\mathrm{ln}{M}_{\mathrm{si}}+{f}_{\mathrm{sa}}\mathrm{ln}{M}_{\mathrm{sa}}\right),\end{array}$

where fc, fsi and fsa are the clay, silt and sand fractions (%) of the soil sample, respectively; Mc, Msi and Msa are the mean diameters of clay, silt and sand that are empirically taken as 0.001, 0.026 and 1.025 mm, respectively.

Consequently, the quantitative relationship between the parameters α and β can be obtained using Eq. (10), associated with the additional constraint of Eq. (11), and the values of α and β can be theoretically solved if the measured volume fraction of the slit-shaped pore or the measured SWRC is known. However, an analytical solution is difficult to derive due to the high nonlinearity of both equations. Here, a trial and error approach was adopted that was much easier than the analytical method. Conveniently, the UNSODA database provided a great deal of soil information, including the measured SWRCs and diverse soil physical properties.

The routine procedure for handling a soil sample involved the following steps. First, given the initial value of α, the value of β was calculated using Eqs. (10)–(12), after which the PoSD was predicted using Eqs. (7)–(9). Subsequently, the SWRC was estimated using the method described in Sect. 3.3. Finally, the values of α were changed repeatedly until the newer predicted SWRC was in good agreement with the measured SWRC and the water content corresponding to a suction head of 5000 cmH2O was within 90 % of the measured data (see Fig. S2). The results for the 48 soil samples indicated that the β values exhibited a broad variation range for all samples, while the α values showed regular changes with the soil texture. The relationship between the sand contents and the α values for 48 samples is shown in Fig. 4, which clearly demonstrates that the values of α are similar to samples with specific sand contents.

Figure 4The α values for 48 soil samples with different sand content. The α values for clay, silt, and fine sand particles of specific samples are listed in Fig. 4, except those of the coarse sand particle, which are the same value of 0.0004 for all of the samples. For the samples with sand content ranging from 10 % to 40 %, two sets of α value are observed. The α values for the silt content less than and more than 50 % are highlighted in red and blue, respectively, thereby reflecting the dominant functions of the silt or clay particle on the hydraulic properties of typical samples.

Therefore, the approach was simplified by setting α as a constant for similar soil textures. The corresponding detailed descriptions are summarized in Table 3. The values of α were inside the range of $\mathrm{3.34}×{\mathrm{10}}^{-\mathrm{05}}$ to $\mathrm{2.12}×{\mathrm{10}}^{-\mathrm{02}}$, which were estimated by Or and Tuller (1999) using a pore model comprising a polygon-shaped central pore connected to the slit-shaped spaces. According to the sand contents of the samples, Table 3 is a reference for determining the α values that serve as input parameters in predicting the SWRC from the PSD hereafter.

Table 3The estimated values of α for various soil textures.

## 3.3 Estimating the SWRC

The values of α and β for the various soil samples facilitated the acquisition of the volume fractions of the slit pores using Eq. (7) and the PoSD using Eqs. (8) and (9). The water contents associated with different pore filling stages could be estimated by substituting the PoSD into Eq. (3), and the pore size and the corresponding suction head could be calculated using Eqs. (2) and (5). The SWRC could be ultimately obtained using the calculated suction heads and water contents.

4 Model validation

## 4.1 Data sources

Twenty-nine soil samples with a wide range of physical properties were also selected from the UNSODA database to validate the model; the codes of the samples are summarized in Table 4 and their detailed information is presented in Table 5. For the soil samples that were not provided with a saturated water content θs, the first data point of the measured SWRC corresponding to the lowest suction head was regarded as θs.

Table 4Codes of the 29 soil samples selected from UNSODA for the model validation.

To generate a detailed PSD, a modified logistic growth model (Eq. 6) was used to fit the measured PSD data. Here, the detailed PSD was generated at diameter classes of 2, 5, 10, 15, 20, 30, 40, 50, 60, 70, 80, 90, 100, 200, 500, 1000 and 2000 µm. The values of α were chosen from Table 3 according to the sand contents of the soil samples. The values of β were obtained by substituting the SSA values predicted using Eq. (11) into Eq. (10). Then, the PoSD was predicted using Eqs. (7)–(9). Finally, the SWRC was estimated using the methods, as described in Sect. 3.3.

The SWRC was also predicted using the traditional method presented in Sect. 2. In the traditional method, the predicted PoSD was equivalent to the PSD (Eq. 1) and was substituted into Eq. (3) to obtain the water content. The corresponding suction heads were predicted using Eqs. (2) and (5).

A scaling approach proposed by Meskini-Vishkaee et al. (2014) was used to compare with the proposed method to demonstrate its prediction performance. The detailed calculation procedures were described by Meskini-Vishkaee et al. (2014).

The van Genuchten equation (Eq. 13) was used to fit the predicted SWRC calculated via the three models (van Genuchten, 1980):

$\begin{array}{}\text{(13)}& \frac{\mathit{\theta }-{\mathit{\theta }}_{\mathrm{r}}}{{\mathit{\theta }}_{\mathrm{s}}-{\mathit{\theta }}_{\mathrm{r}}}={\left[\frac{\mathrm{1}}{\mathrm{1}+{\left(a\mathit{\psi }\right)}^{n}}\right]}^{m},\end{array}$

where θ is the water content (cm3 cm−3), θr is the residual water content (cm3 cm−3), and a, n, m, and θr are fitting parameters. The predicted SWRCs of 29 samples exhibited good fits with an average r2 value of greater than 0.999.

For each set of predictions, the agreement between the predicted and measured water contents was expressed in terms of the root mean square error (ERMS), which is given by

$\begin{array}{}\text{(14)}& {E}_{\mathrm{RMS}}=\sqrt{\frac{\mathrm{1}}{N}\sum _{i=\mathrm{1}}^{N}\left({\mathit{\theta }}_{\mathrm{pi}}-{\mathit{\theta }}_{\mathrm{mi}}{\right)}^{\mathrm{2}}},\end{array}$

where N is the number of measured data points, θpi is the predicted water content and θmi is the measured water content.

## 4.2 Results

The predicted and measured SWRCs in Fig. 5 show that the improved method exhibited good fits with the measured data in the entire range of the SWRC; moreover, the improved method was clearly better than the traditional method and the scaling approach, especially in the dry range (the other 25 samples are listed in Fig. S3). In this study, the scaling approach, which improved the performance of the original MV–VG model via scaling of the parameter n in the van Genuchten equation, performed better than the traditional method for clay (code: 1360), loamy (code: 3190) and loamy sand (code: 3160). However, it performed worse for coarse-textured soil (e.g., sand (code: 3144)), which may result from the relatively small scaling degree of the parameter n and the poor fit of the fitting equation to the measured PSD data in their study. In general, the improved method here applied well to a wide range of soils, while the scaling approach performed better for fine- and medium-textured soils.

Figure 5Measured and predicted SWRCs for clay (code: 1360), loam (code: 3190), loamy sand (code: 3160) and sand (code: 3144).

Table 6 shows the ERMS of the improved method, the scaling approach and the traditional method for samples used in model validation. The ERMS values range from 0.017 to 0.054 for the improved method (with an average of 0.028), from 0.026 to 0.060 for the scaling approach (with an average of 0.037) and from 0.040 to 0.106 for the traditional method (with an average of 0.061). In terms of the ERMS, the improved method provided the best predictions and the traditional method performed the worst. The results also showed that there is a significant difference between the performance of the improved method and the traditional method (p=0.001). Only for the sand samples does the performance of the improved method and scaling approach exhibit a significant statistical difference (p=0.01).

Table 5Basic soil properties of 29 samples for the model validation.

Figure 6The measured and predicted pore volume fraction curves using the improved method and traditional method for clay (code: 1360), loam (code: 3190), loamy sand (code: 3160) and sand (code: 3144).

Table 6The root mean square errors (ERMS) of the predicted SWRC using the improved method, the scaling approach and the traditional method.

The accuracy of the predicted SWRC using the improved method depends on the accuracy of the corresponding predicted pore volume fractions. The calculated and measured pore volume fraction curves in Fig. 6 indicate that the predicted curves using the improved method are more similar to the measured data than those predicted using the traditional method, thereby showing that the proposed method performed better. The errors in the predicted pore fractions using the traditional method mainly occur at the minimum pore size (d≤0.6µm), which proves that the errors of the predicted SWRC using the traditional method originate from the neglect of small pores, such as slit-shaped spaces in a natural sample. The proposed method used the pore model containing slit-shaped spaces to represent the pore space geometry and consequently improved the prediction of the SWRC. However, the uncertainties are unavoidable when choosing the parameters α for unknown media, which is the main factor affecting the accuracy of the predicted SWRCs.

Figure 7The effects of a change in the estimated SSA on the SWRC for loamy sand (code: 3170) and clay (code: 4680). The SSA denotes the accurate value of the specific surface area.

Figure 8The calculated slit width αd, slit length βd and SSAi for loamy sand (code: 3170).

## 4.3 Discussion

### 4.3.1 The suction head calculation in the slit-shaped spaces

When capillary water coexists with adsorptive water in the narrow pores, the capillary and surface forces, including ionic–electrostatic, molecular, structural, and adsorption forces, contribute to the potential energy of water in the slit-shaped pores (Tuller et al., 1999; Iwamatsu and Horii, 1996). When considering only the capillary forces, the drainage potential in slit-shaped pores is given as Eq. (15) (Derjaguin and Churaev, 1992):

$\begin{array}{}\text{(15)}& \mathit{\mu }=\frac{-\mathrm{2}\mathit{\sigma }}{\mathit{\rho }\mathit{\alpha }d},\end{array}$

where μ is the critical potential (J kg−1).

However, the applicability of this formula is limited by the width of the slit. Tuller and Or (2001) defined a critical slit spacing (αd*) by Eq. (16) that classifies the slit sizes responding to capillary drainage and adsorption-dominated drainage. In the case of slit-shaped spaces greater than αd*, the capillary-based slit drainage is applied.

$\begin{array}{}\text{(16)}& \mathit{\alpha }{d}^{*}=\sqrt{-\frac{\mathrm{9}{A}_{\mathrm{svl}}}{\mathrm{4}\mathit{\pi }\mathit{\sigma }}},\end{array}$

where Asvl is the Hamaker constant for solid–vapor interactions through the intervening liquid, usually set as $-\mathrm{6.0}\phantom{\rule{0.125em}{0ex}}×\phantom{\rule{0.125em}{0ex}}{\mathrm{10}}^{-\mathrm{20}}$ J (Tuller and Or, 2001). The value of αd* is 0.591 nm, which means that Eq. (15) could be applied to calculate the drainage potential in the slit-shaped spaces in this study.

In our study, the critical drainage suction head for the minimum central pore calculated using Eq. (5) is 5000 cmH2O, while that of the widest slit-shaped spaces calculated using Eq. (15) is 6202 cmH2O (the potential is converted to the suction head). This result illustrates that all slit-shaped spaces are still filled with water when the suction head is up to the critical drainage suction head for the minimum central pores. On the other hand, the largest slit width calculated from the parameters in Table 3 is 0.24 µm, which is smaller than the minimum pore diameter of 0.6 µm. According to the above analysis, it is reasonable that the volume fractions of the minimum pores include the volume fractions of the minimum central pore and all slit-shaped spaces.

### 4.3.2 The effects of the estimated SSA

The SSA values estimated using Eq. (11) could affect the accuracy of the predicted SWRCs. Figure 7 shows that an overestimation of the SSA prompts the dry range of the SWRC to move in the direction of larger water content, and vice versa. When the estimated SSA value was altered by 10 % and −10 % of its accurate value for loamy sand (code: 3170), the water contents with respect to the highest suction head were higher and lower, respectively, by approximately 0.007 cm3 cm−3 than those of the original SWRC. For clay (code: 4680), the water contents were higher and lower by approximately 0.009 cm3 cm−3 at the same 10 % and −10 % alterations, respectively. Consequently, for coarse-textured soil, the water content and prediction error of the SWRC changed relatively little for the same degree of change in the SSA. Figure 7 also showed that a relatively small error appeared between the calculated and measured SWRCs when the error of the estimated SSA was within 20 %.

Previous work has shown that the SSA of soil is closely dependent upon the soil texture and could be estimated from the soil properties and PSD (Sepaskhah and Tafteh, 2013; Resurreccion et al., 2015). The method used to estimate the SSA in Sect. 3.2 was presented by Sepaskhah et al. (2010), who estimated the SSA based on the geometric mean particle size diameter as shown in Eq. (12) with an r2 value of 0.88. Moreover, the appropriateness of this equation was validated using 64 soil samples by Fooladmand (2011). Sepaskhah et al. (2010) pointed out that the estimation deviations increased distinctly for measured SSA greater than 200 m2 g−1. In the proposed method, the estimated SSA is mainly used to gain the parameter α and β and to estimate the volume fractions of the slit-shaped spaces; thus, the estimation accuracy of SSA influences the dry range of the SWRC (Fig. 7) and equivalently influences the degree of improvement in the predicted SWRC. Therefore, more effort should be made toward developing a more accurate transformation from the soil physical properties to SSA to further improve the prediction of the SWRC.

### 4.3.3 The slit-shaped spaces and the SSA at the sample scale

Since the central pore diameter d is proportional to the corresponding particle diameter D, the slit width αd, the slit length βd and the specific surface area SSAi of each unit cell are associated with the particle size. The calculated values of αd, βd and SSAi of clay, silt, fine sand and coarse sand particles for the loamy sand (code: 3170) are listed in Fig. 8. The results confirm that the pores formed by bigger soil particles are large with a correspondingly large slit width αd; this is similar to the results of Or and Tuller (1999), and the values are on the same order of magnitude. It is common knowledge that larger soil particles tend to have large surface areas, and therefore, the slit length formed by the contact of soil particle edges should be relatively long, leading to the positive relationship between the slit length βd and the particle diameter as shown in Fig. 8. This result is different from that in Or and Tuller (1999), where the slit length βd was inversely proportional to the particle diameter. In addition, the SSAi of the ith particle fractions decreased with an increase in the particle diameter, which is consistent with the findings of Or and Tuller (1999) and is in accordance with the general understanding of the SSA.

5 Conclusions

The traditional models employed to translate the PSD into the SWRC underestimate the water content in the dry range of the SWRC. The errors originate from a setting that the cumulative PoSD is equal to the corresponding PSD, which resulted in an underestimate of the pore volume fraction of the minimum pore diameter range and consequently the water content in the dry range of the SWRC. If slit-shaped pore spaces are taken into consideration when estimating the PoSD with a pore model comprising a circle-shaped central pore connected to slit-shaped spaces, the pore volume fraction of the minimum pore diameter range will be accordingly increased; therefore, the SWRC can be more accurately predicted. The estimation of the α and β values is a key step to predicting the SWRC in the proposed method. The α values were obtained using 48 measured soil samples, and those values served as input parameters for predicting the SWRC; then, the β values were readily calculated using a constraint on the estimated SSA. The validation results illustrate that the SWRCs predicted using the proposed method provided the best predictions of the SWRCs, closely followed by the scaling approach, and the traditional method performed worst.

Data availability
Data availability.

The unsaturated soil hydraulic properties and basic soil properties of samples are available from the UNSODA database (https://www.ars.usda.gov/research/datasets/, last access: 19 August 2018).

Supplement
Supplement.

Author contributions
Author contributions.

CCC developed the method and wrote the paper. DHC performed the analysis and contributed ideas and comments on the method and writing.

Competing interests
Competing interests.

The authors declare that they have no conflict of interest.

Acknowledgements
Acknowledgements.

This research was partially supported by the Special Fund for Basic Scientific Research of Central Colleges (310829162015) and by the National Natural Science Foundation of China (41472220). The authors thank Kang Qian for providing the UNSODA unsaturated soil hydraulic property database.

Edited by: Roberto Greco
Reviewed by: Fatemeh Meskini-Vishkaee and
one anonymous referee

References

Arya, L. M. and Paris, J. F.: A physicoempirical model to predict the soil moisture characteristic from particle-size distribution and bulk density, Soil Sci. Soc. Am. J., 45, 1023–1030, https://doi.org/10.2136/sssaj1981.03615995004500060004x, 1981.

Arya, L. M., Bowman, D. C., Thapa, B. B., and Cassel, D. K.: Scaling soil water characteristics of golf course and athletic field sands from particle-size distribution, Soil Sci. Soc. Am. J., 72, 25–32, https://doi.org/10.2136/sssaj2006.0232, 2008.

Derjaguin, B. V. and Churaev, N. V.: Polymolecular adsorption and capillary condensation in narrow slit pores, Prog. Surf. Sci., 40, 173–191, https://doi.org/10.1016/0079-6816(92)90045-J, 1992.

Fooladmand, H. R.: Estimating soil specific surface area using the summation of the number of spherical particles and geometric mean particle-size diameter, Afr. J. Agr. Res., 6, 1758–1762, 2011.

Hamamoto, S., Moldrup, P., Kawamoto, K., Jonge, L. W. D., Schjønning, P., and Komatsu, T.: Two-region extended archie's law model for soil air permeability and gas diffusivity, Soil Sci. Soc. Am. J., 75, 795–806, https://doi.org/10.2136/sssaj2010.0207, 2011.

Haverkamp, R. and Parlange, J.-Y.: Predicting the water-retention curve from particle-size distribution: 1. sandy soils without organic matter1, Soil Sci., 142, 325–339, https://doi.org/10.1097/00010694-198612000-00001, 1986.

Helland, J. O. and Skjæveland, S. M.: Relationship between capillary pressure, saturation, and interfacial area from a model of mixed-wet triangular tubes, Water Resour. Res., 43, 398–408, https://doi.org/10.1029/2006WR005698, 2007.

Hwang, S. I. and Choi, S. I.: Use of a lognormal distribution model for estimating soil water retention curves from particle-size distribution data, J. Hydrol., 323, 325–334, https://doi.org/10.1016/j.jhydrol.2005.09.005, 2006.

Hwang, S. I. and Powers, S. E.: Using particle-size distribution models to estimate soil hydraulic properties. Soil Sci. Soc. Am. J., 67, 1103–1112, https://doi.org/10.2136/sssaj2003.1103, 2003.

Iwamatsu, M. and Horii, K.: Capillary condensation and adhesion of two wetter surfaces, J. Colloid Interf. Sci., 182, 400–406, https://doi.org/10.1006/jcis.1996.0480,1996..

Jayakody, K. P. K., Shimaoka, T., Komiya, T., and Ehler, P.: Laboratory determination of water retention characteristics and pore size distribution in simulated MSW landfill under settlement, Int. J. Environ. Res., 8, 79–84, https://doi.org/10.22059/IJER.2014.696, 2014.

Jensen, D. K., Tuller, M., Jonge, L. W. D., Arthur, E., and Moldrup, P.: A new Two-Stage Approach to predicting the soil water characteristic from saturation to oven-dryness, J. Hydrol., 521, 498–507, https://doi.org/10.1016/j.jhydrol.2014.12.018, 2015.

Lebeau, M. and Konrad, J. M.: A new capillary and thin film flow model for predicting the hydraulic conductivity of unsaturated porous media, Water Resour. Res., 46, W12554, https://doi.org/10.1029/2010WR009092, 2010.

Liu, J. L., Xu, S. H., and Liu, H.: Investigation of different models to describe soil particle-size distribution data, Advances in Water Science, 35, 68–76, https://doi.org/10.3321/j.issn:1001-6791.2003.05.010, 2003.

Meskini-Vishkaee, F., Mohammadi, M. H., and Vanclooster, M.: Predicting the soil moisture retention curve, from soil particle size distribution and bulk density data using a packing density scaling factor, Hydrol. Earth Syst. Sci., 18, 4053–4063, https://doi.org/10.5194/hess-18-4053-2014, 2014.

Mohammadi, M. H. and Meskini-Vishkaee, F.: Predicting the film and lens water volume between soil particles using particle size distribution data, J. Hydrol., 475, 403–414, https://doi.org/10.1016/j.jhydrol.2012.10.024, 2012.

Mohammadi, M. H. and Meskini-Vishkaee, F.: Predicting soil moisture characteristic curves from continuous particle-size distribution data, Pedosphere, 23, 70–80, https://doi.org/10.1016/S1002-0160(12)60081-2, 2013.

Mohammadi, M. H. and Vanclooster, M.: Predicting the soil moisture characteristic curve from particle size distribution with a simple conceptual model, Vadose Zone J., 10, 594–602, https://doi.org/10.2136/vzj2010.0080, 2011.

Nemes, A., Schaap, M. G., Leij, F. J., and Wösten, J. H. M.: Description of the unsaturated soil hydraulic database UNSODA version 2.0, J. Hydrol., 251, 151–162, https://doi.org/10.1016/S0022-1694(01)00465-6, 2001.

Or, D. and Tuller, M.: Liquid retention and interfacial area in variably saturated porous media: Upscaling from single-pore to sample-scale model, Water Resour. Res., 35, 3591–3605, https://doi.org/10.1029/1999WR900262, 1999.

Pollacco, J. A. P., Webb, T., McNeill, S., Hu, W., Carrick, S., Hewitt, A., and Lilburne, L.: Saturated hydraulic conductivity model computed from bimodal water retention curves for a range of New Zealand soils, Hydrol. Earth Syst. Sci., 21, 2725–2737, https://doi.org/10.5194/hess-21-2725-2017, 2017.

Resurreccion, A. C., Moldrup, P., Tuller, M., Ferré, T. P. A., Kawamoto, K., Komatsu, T., and Jonge, L. W. D.: Relationship between specific surface area and the dry end of the water retention curve for soils with varying clay and organic carbon contents, Water Resour. Res., 47, 240–250, https://doi.org/10.1029/2010WR010229, 2015.

Sakaki, T., Komatsu, M., and Takahashi, M.: Rules-of-Thumb for predicting air-entry value of disturbed sands from particle size, Soil Sci. Soc. Am. J., 78, 454–464, https://doi.org/10.2136/sssaj2013.06.0237n, 2014.

Sepaskhah, A. R. and Tafteh, A.: Pedotransfer function for estimation of soil-specific surface area using soil fractal dimension of improved particle-size distribution, Arch. Agron. Soil Sci., 59, 1–11, https://doi.org/10.1080/03650340.2011.602632, 2013.

Sepaskhah, A. R., Tabarzad, A., and Fooladmand, H. R.: Physical and empirical models for estimation of specific surface area of soils, Arch. Agron. Soil Sci., 56, 325–335, https://doi.org/10.1080/03650340903099676, 2010.

Shahraeeni, E. and Or, D.: Pore-scale analysis of evaporation and condensation dynamics in porous media, Langmuir the Acs Journal of Surfaces & Colloids, 26, 13924–13936, https://doi.org/10.1021/la101596y, 2010.

Shirazi, M. A. and Boersma, L.: A unifying quantitative analysis of soil texture, Soil Sci. Soc. Am. J., 48, 142–147, https://doi.org/10.2136/sssaj1984.03615995004800010026x, 1984.

Tuller, M. and Or, D.: Hydraulic conductivity of variably saturated porous media: Film and corner flow in angular pore space, Water Resour. Res., 37, 1257–1276, https://doi.org/10.1029/2000WR900328, 2001.

Tuller, M. and Or, D.: Water films and scaling of soil characteristic curves at low water contents, Water Resour. Res., 41, 319–335, https://doi.org/10.1029/2005WR004142, 2005.

Tuller, M., Or, D., and Dudley, L. M.: Adsorption and capillary condensation in porous media: Liquid retention and interfacial configurations in angular pores, Water Resour. Res., 35, 1949–1964, https://doi.org/10.1029/1999WR900098, 1999.

van Genuchten, M. T.: A closed-form equation for predicting the hydraulic conductivity of unsaturated soils, Soil Sci. Soc. Am. J., 44, 892–898, https://doi.org/10.2136/sssaj1980.03615995004400050002x, 1980.

Zhuang, J., Jin, Y., and Miyazaki, T.: Estimating water retention characteristic from soil particle-size distribution using a non-similar media concept, Soil Sci., 166, 308–321, https://doi.org/10.1097/00010694-200105000-00002, 2001.