Much attention has been paid on the soil-water characteristic curve (SWCC) during decades because it plays great roles in unsaturated soil mechanics. However, it is time-consuming and costly to obtain a series of entire saturation-suction data by experiments. The curves acquired by directly fitting empirical equations to limited experimental data are greatly different from the actual SWCC, and the relevant soil parameters obtained by inaccurate curve are also incorrect. Thus, an improved prediction method for more accurate entire SWCC was established. This novel method was based on the analysis of shape characteristics of SWCC with three critical points , , and under the hypothesis of geometrical symmetric relation. The theoretical computation was specifically deduced under conventional Gardner, VG, and FX models, respectively, and then inferred on different soil types of 45 collected SWCC datasets. This geometrical symmetric relation exhibited well in all these three conventional empirical equations, especially in Gardner equation. Finally, a series of filer paper tests on sand, silt, and clay were also carried out to acquire entire SWCC curve for the verification and evaluation of the proposed geometrical method. Results show that this improved prediction method effectively decreases deviation resulting from directly fitting empirical equations to limited data of wide types of soils. The averaged improvement was larger under VG equation than under Gardner and FX equation. It proved that the accuracy of predicting greatly depends on the shape characteristic point of maximum curve curvature (point ), other than the number of points. This research provides a novel computation method to improve prediction accuracy even under relative less experimental data.

1. Introduction

Most of the sedimentary soils widely distributed on the surface belong to unsaturated soils. With the acceleration of human development, more and more engineering practices are related to unsaturated soils, such as the rainfall-induced slope failures, construction about the road and railway subgrades [1], and even the realization of “sponge cities” [2]. Since the soil-water characteristic curve (SWCC) is used as the basis for the prediction of other unsaturated soil parameters, such as the permeability [3, 4], shear-strength functions [5], and diffusion coefficient [6], consequently, it is important to have a reasonably accurate characterization of soil-water characteristic curve. And a great number of research studies have focused on the problem about the acquisition of SWCC.

There are two common methods to obtain an entire SWCC experimental dataset. The first method is to directly measure suction through complex and time-consuming experiments. Actually, the pressure plate method is commonly used in the suction measurement of unsaturated soils [7], but the maximum measurement range of this method only can reach 2000 kPa [8], which is not enough for the soil that has a great quantity of fine particles, especially for clay. Other conventional experimental methods such as Tempe cell, solution method [9], filter paper method, and axis translation method are tedious and have narrow measurement range. Some monitors could also be used to simultaneously measure soil suction in the model test. Recently, some of newer equipment for determination of SWCC comprise HYPROP, small-scale centrifuge, and dew-point chilled mirror [10, 11]. However, each technique has its own limitation in terms of suction. Sometimes, multiple measurement methods are used in combination to acquire an entire SWCC. Consequently, the excessive costs associated with direct measurement of SWCC have encouraged the pursuit of a great quantities of indirect methods to acquiring SWCC based on some empirical relationship and statistical regularity of known soil parameters or limited experimental data, including pedotransfer functions [12], physical-empirical equation [13, 14], fractal equation [15, 16], and empirical formula method [17]. In order to get more accurate properties of SWCC, these methods are usually combined with each other [18, 19]. Additionally, the machine learning or intelligent algorithms are also applied with these methods [20, 21].

Besides, considering the relation between the grain-size distribution and SWCC [2224], some equations can be used to estimate unimodal SWCC with unimodal grain-size distribution [11, 25]. Pore size distribution function (PSDF) and volumetric shrinkage curve (VSC) also play an important role in estimating SWCC. Fredlund and Xing [26] proposed the SWCC equation by integrating the assumed pore size distribution function. The PSDF also had a significant effect on the saturated coefficient of permeability for sandy soils. The volumetric shrinkage curve (VSC) is commonly used to convert the SWCC in the form of gravimetric water content (w-SWCC) into a curve that is in the form of degree of saturation (S-SWCC) [27]. Lin and Cerato [28] showed that the VSC bears a relationship to the soil-water characteristic curve during drying from a near-saturated condition. Zhai et al. [29] put forward a framework used to measured w-SWCCs, and the VSC of soil in a loose condition is used as the input information for the estimation of the w-SWCCs of soil under different void ratios.

Because of its wide application and the simplicity of the calculation process, the empirical equation method has become the most common method for estimating SWCC. The common empirical equations that have been proposed include Gardner equation [30], van Genuchten equation (VG equation) [17], Fredlund and Xing equation (FX equation) [31], and Brooks–Corey equation [32]. Recently, some new equations were put forward. Considering the existence of bimodal SWCC and the limitation of conventional equations, Wijaya and Leong [33] proposed a new equation for both unimodal and bimodal SWCCs. However, there is an evident problem that final fitting curves would greatly deviate from the actual SWCC when directly fitting empirical equation to limited data acquired by limitation of suction measurement range of the frequently used suction measure method (Figure 1). Siller and Fredlund [35] stated that the lowest measured water contents and soil suctions appear to influence the magnitude of the residual water content. The fitting routine tends to overestimate the residual water content if the experimental data do not extend well into the zone of desaturation. Consequently, the important values of unsaturated soil predicted by the incorrect SWCC also greatly differ from actual values. Moreover, the related design and evaluation results based on the inaccurate values will be unsafe.

However, there are few scholars currently studying the above problems. Wang et al. [36] combined statistical algorithms to predict entire SWCC from the Bayesian perspective, nevertheless, that approach is relatively complicated. Ren et al. [37] proposed a method for estimating SWCC with limited data, but they did not quantitatively testify the feasibility of that method.

This paper tried to acquiring more accurate SWCC from the perspective of geometry. Observing Figure 1, there seems to be a symmetric relation between the point which is the largest absolute value of slope and two points and which are two extreme points of curvature on SWCC obtained by Gardner equation in semilogarithmic coordinate. Therefore, in the subsequent part of this paper, the geometrical relation between three points would be discussed, and a geometrical method based on this critical relation for preprocessing limited data before directly fitting with empirical equation was put forward.

2. Methodology

2.1. The Geometrical Relationship of , , and of SWCC under Gardner Equation and Validation

The special geometrical relation is found on the SWCC fitting by Gardner equation; thus, the characteristics of Gardner equation need to be studied. Gardner used a function of three parameters (equation (1)) to describe the soil-water characteristic curve, and it is one of the simplest empirical equations.where is the volumetric water content; is the saturated water content; is the residual water content as the ratio of the volumetric water gradient to suction approaches zero; is the fitting parameter related to the air-entry value of value (bubbling pressure); is the parameter that has relation with desorption rate when the matric suction is bigger than the air-entry value of value; and is the suction of soil, kPa.

Considering the possible geometrical symmetry of three points , , and which have specific mathematical significance on SWCC in semilogarithmic coordinate, theoretical derivation and checking calculation of substantial datasets were done to prove whether this relation exists.

Firstly, a great number of SWCC datasets of different type of soils (sand, silt, and clay) have been collected to do numerical calculation in order to qualitatively prove that assumption is supported by a quantity of practical experimental data. Every parameter in equation (1) and the suction value of three points are displayed in Table 1, and the value of and (equations (2) and (3)) is used to judge the geometrical relation between three points (Table 2). The calculation results of 45 datasets show that is just equal with , and the values of and in all datasets are 1, which means the symmetric relation between the point and two points and may exist.

Although the possible relation between three points have been testified by actual data, a series of mathematical derivations also are needed in order to completely and truly prove that geometrical symmetric relation. Thus, the curves of absolute value of slope in semilogarithmic coordinate can be acquired by the following equation:and the point is the extreme point of the curves of absolute value of slope in semilogarithmic that can be derived when the second derivation of (equation (5)) is zero.

It is obvious that achieves zero when , so the location of point can be determined as . Then, the expression of curvature of SWCC in semilogarithmic coordinate can also be defined as the following equation:

The points and are two extreme points of curvature on SWCC in semilogarithmic coordinate; consequently, the points meet the conditions that or . Although it is difficult to find all the analytic solution of , two solutions and can be found out, where is a complicate function about , , and . Thus, the points and can be determined as and . Then, the symmetric relation between the point and points and in semilogarithmic coordinate can be verified by equations (7) and (8). Plugged the specific coordinates of point and points and acquired above into equations (7) and (8), and it was found that those equalities held.

of FX equation can be acquired by the equation [56].

Through the above deduction, the symmetric relation between three points can be theoretically proved. In conclusion, the symmetric relation between point and points and in semilogarithmic coordinate can be quantitatively and qualitatively proved based on numerical checking of 45 SWCC experimental datasets and mathematical derivations. Consequently, the specific position of point can be acquired when the points and are known by making use of this relation.

2.2. The Geometrical Relationship of Points , , and of SWCC under VG Equation and FX Equation

Besides Gardner equation, VG equation put forward by van Genuchten (equation (9)) and FX equation (equation (10)) raised by Fredlund are also the popular empirical equations used to acquire SWCC. Considering the VG equation is one of the most frequently used empirical equations and the FX equation can give the best fit among the equations [56], they are selected to testify whether this graphic symmetric relationship exists in other equations.

The same calculations were done to 45 selected experimental datasets in order to find three special points , , and whether have the same relation on the curve fitting to the VG equation and FX equation. of points , , and is shown in Table 2, and the value of and (equations (2) and (3)) is also used to judge the geometrical symmetric relation of three points. It is found that and of Gardner equation and FX equation are not exactly equal to 1, but they are still very close to 1. In consideration of the fact that these empirical equations are all used to fit with the same dataset of discrete points, despite their slightly different forms, the curves acquired are roughly same. Thus, the positions of points , , and on the curves obtained by VG and FX equations are similar to those on the curve obtained by Gardner equation. In conclusion, the graphic symmetric relation between three special points , , and can also be considered generally valid on the curves fitting by VG equation and FX equation.

2.3. The Physical Meaning of Points , , and

Although the mathematical significance and geometrical symmetric relation of three points , , and have been proved by theoretical derivation and checking calculation of substantial datasets, the specific physical meaning of these points is also needed to be discussed. Some researchers have already studied SWCC from the geometrical perspective. Siller and Fredlund [35] used two points of extreme values of curvature to divide SWCC into three zones (capillary saturation zone, funicular zone, and zone of residual saturation). Fredlund and Xing [26] put forward a geometric solution to determine parameters , , and in equation (12) by drawing a tangent line through point and, respectively, intersecting the approximate straight lines of low-suction part (capillary saturation zone) and high suction part (zone of residual saturation) of SWCC. They considered that two points and related with the air-entry value (the matric suction where air starts to enter the largest pores in the soil) and the residual water content are the intersection points. However, the points and put forward in this paper are on the SWCC which are the extreme values of curvature. The points and are the intersection points of tangent line of point and lines of low-suction part which are out of SWCC. In view of graph, they are all around the curvy part of SWCC and close to each other (Figure 2). Additionally, the points and are on the SWCC, so they may have stronger physical connection with SWCC than points and , especially the change of curvature. Therefore, points and also have relationship with the air-entry value and the residual water content . The air-entry value and the residual water content are closely associated with the mineral composition, pore structure, and dense state; among them, the mineral composition and pore structure are the basic factors. Moreover, these two points are crucial for determining the position and shape of SWCC because of their important mathematical meaning. In other words, SWCC is a typical sigmoid curve; once the locations of mutation point of curvature are confirmed, the rough shape and position of SWCC can be obtained.

The point is the inflection point of SWCC (the point is largest absolute value of slope), from the perspective of pore size distribution; it corresponds to the point with the maximum value of frequency in the pore size distribution curve [57]. Because of the log-normal distribution of pore size of most soils, the maximum value of frequency in the pore size distribution curve can be considered as the mean pore size of the soil sample [37]. Besides this, Ren and Santamaina [58] and Ren et al. [59] mentioned that the mean pore size can be expressed by void ratio and specific surface area. Therefore, the position of point , void ratio, and specific surface area may have some inner link.

2.4. The Significance of Point

When the empirical equations are used to fitting available datasets, it is significant to have an effective and precise standard to judge whether datasets can generate accurate outcomes. Leong and Rahardjo [56] mentioned that the data used to obtain SWCC should include points after which is not directly available value before fitting. Fredlund and Xing [26] considered that the magnitude of will generally be in the range of 1500 to 3000 kPa but not an exact value of .

According to the above discussion about the physical and geometrical relation between the point and value of , this paper desired to quantitatively check whether point is also important to acquiring accurate SWCC. The SWCC datasets of three types of soils were chosen. Then, the half and quarter points of datasets were used to fit to equation (12), which is the empirical equation of low sensitivity according to the analyses carried out by Leong [56]. Moreover, the points before the point and the dataset including the point just larger than point were also been calculated. Because all the empirical equations are the overparameterized functions (the parameters of the functions are not independent), some different fitting parameter combinations result in the same curves for the same dataset. Thus, equation (13) was used to evaluate the difference between the curves obtained by different datasets, where and are two different SWCC suction-saturation equation needed to be compared and is the proximity of the two curves in some extent; the smaller is, the closer the two curves are. The calculating results of different types of soils are shown in Figure 3 and Table 3.

From Figure 3 and Table 3, it was easy to drive a conclusion that curves acquired by half points and a quarter points were almost coincided with the SWCC obtained by original dataset. However, the difference between the curves acquired by the points before point and the original SWCC was relatively large. And the gap reduced when the fitting dataset just including the point . In some occasions, the SWCC obtained by the dataset including point is almost similar to the original curve. Therefore, a conclusion could be drawn that the accuracy of empirical equations fitting depends on the dataset whether it includes point other than the number of points. Then, in this paper, the dataset which is of sigmoid shape and at least reaches the maximum curvature (point ) can be regarded as a “entire” SWCC dataset.

2.5. A Geometrical Method for Entire SWCC

As what are discussed above, the experimental datasets sometimes miss the points of part of funicular zone and zone of residual saturation because of the limitation of the range of the measuring instruments. Thus, this paper puts forward a method to predict the position of point and point based on the graphic symmetric relationship of points , , and to acquire more accurate SWCC.

Firstly, the first critical point is the first extreme value of curvature on the SWCC that is obtained by directly fitting empirical equations to limited data. because the upper part of SWCC acquired by fitting to limited dataset is almost coincident with that of actual SWCC. In most cases, the first point of extreme value of curvature () on these two curves is almost the same. The inflection point on the SWCC of Gardner equation is , and the inflection point acquired by other empirical equations is very close. Therefore, point can be seen as the inflection point on the SWCC, and can be easily obtained from the saturated samples. The value of is closely associated with the soil characteristics, like mineral composition, pore structure, and dense state. In consideration of that, of different datasets under different empirical equations was calculated (Table 1), and the average values of same dataset under three different empirical equations were gained. Afterwards, these 15 values of different types of soils were also averaged; thus, the statistical values of of three types of soils were obtained. It was found that these three statistical values of were in line with the characteristics of different soils that the finer particle size is, the higher values of is. Consequently, the empirical values of of sand, silt, and clay were, respectively, chosen as 0.054082, 0.062445, and 0.079245. The points , , and are almost in a line, considering the coincidence of the upper part of SWCC acquired from limited dataset and that of actual SWCC; therefore, the line of the SWCC acquired by limited data almost coincides with the line of the actual SWCC. In other words, predicting point is on the line . Consequently, could be derived by the empirical and known (equation (14)), and the value of could be gained by the line and the obtained value of (equation (15)). Eventually, the point could be acquired by predicting point and known point (equations (16) and (17)). Moreover, point was added to limited dataset, and the new dataset was fitted to empirical equations to get a more accurate SWCC.

However, in some cases, the lowest value of of limited dataset may be lower than the empirical value of this paper adopted. Then, could be chosen as . An example shows procedures and processing result (Figure 4), and diagrammatic sketch about specific procedures of this geometrical predicting method is shown in Figure 5.

3. Materials and Experiment

3.1. Materials

The soil sample used in this test was remolded soil. Sand, silt, and clay are all configured by forging kaolin and standard sand according to the predetermined mass ratio. Two kinds of particle sizes were used (0.045 mm and 0.001 mm). Standard sand is made according to the standard named GB/T17671-1999, and the grain diameter ranges from 0.08 mm to 2 mm. The specific ratio is shown in Table 4, and the grain-size distribution curves of three kinds of sample are shown in Figure 6.

3.2. Filter Paper Method

In the measurement of suction of unsaturated soils, compared with the pressure plate method or other methods, the filter paper method is more time-consuming and cost-saving, and it also has larger measuring range [60]. Therefore, this paper uses the filter paper method to obtain the entire suction-saturation datasets of sand, silt, and clay.

This method was first proposed by Gardner [61], and since then, many researchers have been conducting further research. Feuerharmel et al. [62] and Agus and Schanz [63] compared and analyzed the filter paper method and other methods and found out that the value obtained by the filter paper method has higher accuracy. Marinho and Oliveira[64] also made a deeper discussion on the operation procedures and control conditions of the filter paper method.

In filter paper experiment, the most important procedure is to determine the calibration curve. Fawcett and Collis [65], Leong and Rahardjo [56], and other scholars, respectively, used the solution method or pressure plate instrument and other methods to calibrate the Whatman No. 42 filter paper calibration curve and the applicable range. The specific calibration equation is given. Scholars in China such as Tang et al. [66], Bai et al. [67], and Wang et al. [68] have measured the calibration curves of domestic double-loop filter papers. The ASTM (American Society for Testing and Materials) has formulated the relevant test equipment and operating procedures of the filter paper method into a relevant standard-ASTM D5298-10, which also specifies the calibration curve of Whatman No. 42 filter paper. This paper used the calibration curve equation (equations (18) and (19)) put forward by Leong to calculate suction, where is the mass of wet filter paper; is the mass of dry filter paper; and is the water content of processed filter paper:

In addition, considering the hysteresis characteristics of SWCC, Houston et al. [69] and Leong et al. [70] have demonstrated through a large number of experiments that there is no significant difference in matrix suction during the hygroscopic process and dehumidifying process through a large number of experiments. Muñoz et al. [71] proposed that using wet filter paper can eliminate the hysteresis characteristic of SWCC, but the test results show that the results obtained with dry filter paper and wet filter paper are not much different. In summary, this paper uses the dry Whatman No. 42 filter to measure the suction of soil samples under different water content. The sketch of the ring-knife layout and the photos of experimental process are shown in Figure 7. And the results are shown in Table 5.

4. Results and Discussion

4.1. Verification Using the Experimental Results of Filter Paper Experiments

The datasets acquired by filter paper experiment covered entire suction range; thus, the curves they formed could be considered as the actual SWCC. Firstly, a part of the data of filter paper experiment were selected in order to make the curve they formed did not have two points of extreme value of curvature. Therefore, this processed dataset could be regarded as limited experimental data. The proposed method was applied to these limited data and compared quantitatively performance of this method applied to the experimental data by equations (13) and (20). In equation (20), is the index parameter that quantifies improvement degree of the proposed method; is the deviation between the curve obtained by directly fitting limited experimental datasets to empirical equation and actual SWCC; and is the deviation between the curve obtained by the proposed method and the actual SWCC. The larger value is, the greater the improvement the method makes. The specific comparison charts and calculation results charts are shown in Figure 8 and Table 6.

From Table 6 and Figure 7, it is found out that the improvement this method made on three empirical equations ranges from 0.948 to 0.459; in other words, the curves gained from the method are much closer to the actual SWCC than the curves obtained by directly fitting limited data to empirical equations. Thus, this method performs well when it is applied to the data from filter paper experiments on three different types of soils. However, the performance of this method on these three datasets is not stable, namely, it does not show any preference of applicability to types of soil. In the perspective of empirical equations, the value of and of FX equation is lower than that of VG equation and Gardner equation, which means the curves obtained by FX equation were closer to the actual SWCC. That is because the FX equation has a restriction point (1,000,000, 0) and the deviation resulting from directly fitting FX equations to limited data relatively decreases. When this method is used on FX equation, there two points are added to the limited data, besides point , so that the final obtained curve is closer to actual SWCC.

4.2. Verification Using the Previous Experimental Data

Considering the limitedness and occasionality of only three experimental datasets, this paper selected and collected the entire previous experimental data according to the standard mentioned above to verify the effectiveness of the proposed method. Similarly, the entire data were selected as above processes, and in the dataset of silt and clay, the points whose suction were lower than 2,000 kPa and were selected in order to simulate limitation of the measurement range. The calculation results and comparison charts are shown in Table 6 and Figure 9.

Comparing with the performance of this method on filter paper experimental data, its application of silt on FX equation was not well. Generally, some improvement of this method has been made at different degree ranging from 0.815 to 0.097, expect the result of silt on FX equation. However, the value of and under FX equation was also lower than that under VG equation and Gardner equation. From Figure 10, the performance of this method did not greatly vary with soil types, and the averaged improvement was the largest on VG equation among three empirical equations.

4.3. Discussion

Overall, the data processed by this method could acquire more accurate SWCC to some extent, comparing with the limited data. From a series of calculation, it is found that the values of and under FX equation are lower than those under VG equation and Gardner equation. That is because the restriction point (1,000,000, 0) makes the fitting results of FX equation more approximate than those of other empirical equations. For sake of this, the improvement of adding point on FX equation is not stable. Among three empirical equations, the averaged improvement of this method is the largest on VG equation. Therefore, it is suggested this method would achieve the better effectiveness when it is combined with VG equation.

In addition, this study was focussed on the prediction about typical unimodal SWCC. The empirical equations this method based on and discussed in this paper were not applicable to both unimodal SWCC and bimodal SWCC. Therefore, the proposed method is not suitable for the prediction of bimodal SWCC. SWCCs of residual and colluvial soils commonly show bimodal behavior; thus, the applicability of the proposed method on those soils is limited.

What this method needs to be improved in the future is how to make the value of or residual water content closer to reality because the value of is the averaged statistics of a quantity of entire datasets in this paper, other than some regularity of the soil characteristics. According to the analysis about the physical meaning of three points , , and and the regularity presented by calculation above, the value of may have some qualitative relations with soil properties (mineral composition, pore structure, and dense state). Thus, the future work can be focused on finding some specific relation between some soil index and the value of residual water content , aiming at obtaining more accurate whose physical meaning are also more specific.

5. Conclusions

This paper proposed a geometrical method based on a symmetric relation of three critical points (point which is the largest absolute value of slope; two points and which are two extreme points of curvature) on SWCC in semilogarithmic coordinate to preprocess limited data before directly fitting with conventional empirical equations. Meantime, a series of theoretical derivation and confirmation computation proved this symmetric relation of three critical points exhibited well in all three conventional empirical equations, especially in Gardner equation. The accuracy of SWCC obtained by empirical equations depends on the dataset whether the formed curve includes the critical shape characteristic point of maximum curvature (point ). If this point has not been reached in experimental data, even many numbers of points obtained, the fitting results directly by the empirical model could be large error. Three datasets of sand, silt, and clay acquired by the filter paper method were used to testify applicability of this method. The performance of this method did not greatly vary with soil types, but it was floating a lot on different empirical equations. Through quantitative analysis, the method generally performed well on sand, silt, and clay. In perspective of empirical equations, the evident improvement can be made when the proposed method is combined with VG equation, Gardner equation, and FX equation, but averaged improvement was the largest on the VG equation among three empirical equations. Generally, this method provides good evidence of wide application potential. And future work could be concentrated on finding out the relation between residual water content and soil characteristics instead of using averaged statistical value.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.


The research work herein was supported by the National Natural Science Foundation of China (No. 41702299) and Foundation of State Key Laboratory of Frozen Soil Engineering (No. SKLFSE201916).