Mathematical Problems in Engineering

Mathematical Problems in Engineering / 2021 / Article

Research Article | Open Access

Volume 2021 |Article ID 9956824 | https://doi.org/10.1155/2021/9956824

Jie Zhou, Junjie Ren, Zeyao Li, "An Improved Prediction Method of Soil-Water Characteristic Curve by Geometrical Derivation and Empirical Equation", Mathematical Problems in Engineering, vol. 2021, Article ID 9956824, 22 pages, 2021. https://doi.org/10.1155/2021/9956824

An Improved Prediction Method of Soil-Water Characteristic Curve by Geometrical Derivation and Empirical Equation

Academic Editor: Paolo Crippa
Received09 Mar 2021
Revised15 May 2021
Accepted08 Jun 2021
Published17 Jun 2021

Abstract

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.


TypeData resourceGardner equationVG equationFX equation

Sand[38]0.21060.03397.503.140.03604.51100.560.010.034745.525.0912.820.400.0349
[38]0.42460.01531.116.790.01421.008.180.590.01101.211.1032.790.800.0135
[26]0.38840.073522.144.120.104458.693.3721.240.104330.7218.645.120.950.0941
[39]0.34640.06552.863.020.05591.64196.120.010.06076.62  1071.925.480.650.0607
[40]0.32630.080223.714.470.067820.265.200.560.034077.3418.965.910.740.0607
[40]0.26560.08621.714.350.05681.236.500.260.06566.67  10101.295.460.620.0695
[41]0.34180.08602.177.860.07411.68282.430.010.09017.31  10101.8020.320.440.0834
[42]0.29590.08540.510.430.07400.0317.090.010.07563.020.0732.040.120.07830.054082
[20]0.24990.11095.884.520.08293.7635.210.040.098723.374.3510.780.290.0975
[20]0.34480.10044.301.770.08421.7126.440.030.06237.68  10102.164.250.460.0823
[43]0.45270.02632.511.140.01931.4025.110.030.01237.27  1071.541.691.280.0193
[44]0.35230.0437180.152.170.040676.67115.400.010.0367104.89109.2910.400.330.0403
[44]0.35980.0447183.692.330.0377104.533.990.300.011088.64120.656.060.530.0311
[45]0.18380.01922.961.030.01460.8669.990.010.01726.421.304.310.550.0170
[45]0.24440.03563.641.550.02831.4837.790.020.021917.683.0715.090.410.0286
[34]0.40280.03771769.980.910.04042.75  1040.849.860.0362333.005152.190.904.270.0381
[34]0.40780.03911280.151.220.06941.32  1070.986488.410.0694165.45867.802.350.650.0593
[7]0.42260.0769742.251.940.0754541.702.210.640.0710259.60611.8299.270.140.0744
[7]0.38410.03401993.110.640.09282.32  1080.601034.880.04615.27  105250.670.860.970.0576
[38]0.33000.01563927.640.630.03829117.250.591.520.01961415.781329.200.571.100.0245
[38]0.30410.00176104.890.500.0127101.350.960.070.007647.281.13  1047.620.130.0073
[46]0.48040.11866.082.820.12187.782.501.610.11980.744.827.900.390.1201

Silt[47]0.45560.0341514.290.650.078853.341.510.380.05118532.0436.422.020.500.05470.062445
[47]0.44170.0141877.770.570.020733.601.580.200.01791823.6233.721.790.450.0176
[48]0.37940.036170.110.910.039761.190.950.940.0387213.4929.231.011.180.0382
[44]0.32200.0325513.710.750.0397847.530.721.350.025214.16317.821.990.390.0325
[45]s0.34560.045663.621.270.016521.192.050.280.0681191.3828.541.640.920.0434
[45]s0.43710.092323.361.220.096444.961.101.700.09491.6428.0962.510.160.0945
[49]0.39910.116021.973.110.199212.96161.640.010.11362146.4114.038.070.230.1429
[50]0.43680.127325.971.160.13395.621.830.150.13361.81  1087.451.530.350.1316
[42]0.49560.0498302.630.960.0532281.190.990.970.066537.96183.241.250.850.0565
[42]0.43640.0662945.310.520.04452741.660.522.220.07964.08641.511.730.970.0634
[51]0.69450.01022477.660.620.00394397.950.591.320.008033.373141.601.500.610.0074

Clay[7]0.60180.03432363.760.340.03734.80  1050.393.710.03381.72  105122.361.050.480.03510.079245
[7]0.44820.07377430.100.670.06472858.870.770.620.0835235.383.83  1052.02−2.620.0740
[14]0.52040.01082961.290.880.01323007.030.891.020.0143297.142599.031.121.070.0128
[43]0.49130.243466.270.620.27622.14  1060.53210.140.29262.681017.68135.80−0.030.2707
[44]0.44050.1567256.591.010.085434.7669.130.000.14804702.9353.192.830.300.1300
[44]0.55300.1594432.810.720.1179153.420.830.500.172516.30226.3717.840.050.1499
[52]0.46000.164711.421.560.09894.622.300.250.120996.735.191.950.560.1282
[53]0.37810.09152916.261.310.0432981.451.180.480.06165925.741024.262.710.400.0654
[53]0.48880.11962569.651.320.0841758.273.390.170.07283787.67866.213.430.330.0922
[54]0.90310.06203535.160.730.05371463.290.870.610.06137600.92492.901.100.650.0590
[54]0.64670.02948007.590.740.02884786.720.810.760.02491.91  1041268.510.940.760.0277
[55]0.47020.01591757.230.770.01711909.350.761.050.016160.421264.731.660.570.0164


TypeData resourceGardner equationVG equationFX equation

Sand[38]7.504.8111.691.001.004.704.504.951.341.175.624.866.441.010.94
[38]1.110.801.541.001.001.070.791.501.151.161.130.961.321.300.92
[26]22.1414.8233.091.001.0023.6915.9637.280.921.1521.6615.2630.501.020.98
[39]2.861.754.681.001.002.862.045.080.931.702.321.683.161.010.95
[40]23.7116.6733.741.001.0022.6616.0932.421.071.0522.2616.5029.691.010.96
[40]1.711.222.391.001.001.511.102.101.111.051.571.152.090.990.94
[41]2.171.692.781.001.001.911.672.310.981.471.911.682.201.231.15
[42]0.510.0210.921.001.000.490.076.820.831.390.080.070.081.040.94
[20]5.884.308.031.001.004.733.746.611.101.424.954.185.770.980.92
[20]4.301.999.281.001.003.951.698.890.520.963.871.939.140.711.23
[43]2.510.758.351.001.002.060.926.261.771.382.180.924.970.980.95
[44]180.1594.04345.121.001.00104.5680.12140.330.891.11124.03103.02147.961.020.95
[44]183.6998.81341.481.001.00140.3685.91234.851.111.05144.38107.77191.001.010.96
[45]2.960.8210.711.001.002.911.2310.960.821.541.681.063.061.291.30
[45]3.641.538.651.001.003.611.888.760.751.363.332.923.801.050.98
[34]1769.98400.897814.731.001.001786.54511.896112.250.930.981937.77419.467910.790.900.92
[34]1280.15416.533934.361.001.001707.71606.614685.020.910.981204.32624.032315.771.010.99
[7]742.25356.391545.851.001.00679.01333.141396.411.001.01669.57469.321209.230.741.66
[7]1993.11249.911.59  1041.001.002318.54458.221.16  1040.930.99632.86115.832923.600.940.90
[38]3927.64470.873.28  1041.001.004914.27619.073.99  1040.951.003532.32500.162.55  1041.031.01
[38]6104.89431.308.64  1041.001.006651.27482.948.31  1040.880.966488.40472.308.47  1041.310.98
[46]6.083.5210.521.001.006.443.6411.190.960.975.634.407.161.040.98

Silt[47]514.2965.784020.821.001.00101.6832.10325.641.051.0165.1329.16135.120.960.91
[47]877.7784.289141.811.001.0093.3724.84355.601.081.01149.0726.903686.971.021.87
[48]70.1116.15304.371.001.0065.5315.57275.971.001.0057.1013.74219.800.980.95
[44]513.7187.463017.351.001.00549.3496.483123.201.001.00501.14103.323107.871.131.16
[45]63.6221.83185.381.001.0039.5615.63101.321.071.0147.0918.84109.720.960.92
[45]23.367.6571.371.001.0027.719.0184.620.970.9928.7827.4430.671.441.33
[49]21.9719.5524.701.001.0013.3012.9413.761.411.2016.8213.3820.720.980.91
[50]25.978.2581.731.001.0015.834.5855.081.091.0117.746.0247.150.960.90
[42]302.6373.521245.691.001.00291.4472.391174.081.001.00279.6177.88911.760.950.92
[42]945.3273.251.22  1041.001.00596.6965.225426.520.971.00808.67370.401918.411.081.11
[51]2477.66270.502.27  1041.001.003137.80189.462.80  1040.790.782012.34135.662.16  1041.200.88

Clay[7]2363.7550.311.11  1051.001.001.61  104943.472.69  1050.960.992385.6673.031.10  1050.781.10
[7]7430.101014.255.44  1041.001.005347.42765.623.75  1041.021.001.40  104232.702.53  1050.670.71
[14]2961.28629.331.39  1041.001.002945.55634.5113669.101.001.003020.23662.401.15  1040.900.88
[43]66.277.92554.331.001.0091.2714.69563.920.941.0064.346.96527.201.070.95
[44]256.5967.82970.841.001.00237.6064.70940.800.721.06148.0047.51908.480.871.60
[44]432.8167.722766.341.001.00233.8525.602179.711.261.01248.25223.76272.490.940.90
[52]11.424.7427.531.001.0010.644.6425.770.941.079.264.0319.430.940.89
[53]2916.261035.338214.351.001.002694.63967.337548.710.881.012651.25978.067808.750.751.08
[53]2569.65903.527308.181.001.002281.71871.067791.750.841.282284.17865.657995.740.751.29
[54]3535.16513.192.44  1041.001.002560.38419.471.60  1041.031.011348.69297.344066.300.830.73
[54]8007.591254.265.11  1041.001.006715.741105.244.11  1041.011.004612.57733.924.04  1041.231.18
[55]1757.23304.321.01  1041.001.001877.85325.121.08  1040.981.001836.05349.761.18  1041.061.12

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.


Soil typePoints situation

Sand [38]24 points5.0377.4210.551
13 points4.8856.1460.5860.00508
6 points4.9826.6120.5900.01032
Before 5.02511.5520.4170.04922
Including 5.07210.7990.4520.02913

Silt [38]20 points10394.2030.4803.354
11 points8089.2300.5083.3480.03059
6 points9028.3410.4853.3930.02031
Before 2095.4110.5471.8800.02828
Including 7254.0440.5252.2230.01836

Clay [53]20 points1097.9021.7560.740
11 points1120.3151.7670.7500.00575
6 points1149.2661.7470.7750.02988
Before 878.8772.6530.4930.20576
Including 903.9822.0480.6290.16672

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.


TypeSpecific gravity of soil solidsDry densityIngredients
Sand (%)Silt (%)Clay (%)

Sand sample2.631.6577158
Silt sample2.661.30405010
Clay sample2.691.20232057

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.


SandSiltClay
Water content of filter paper (%)Water content of soil (%)Suction (kPa)Water content of filter paper (%)Water content of soil (%)Suction (kPa)Water content of filter paper (%)Water content of soil (%)Suction (kPa)

11.8510.00514042.38823.9170.0142164.4998.8690.00522288.953
23.6740.0092247.75324.9880.0141833.61012.6230.00912458.153
28.5230.0181060.18128.5580.0301054.53524.3010.0212039.631
38.2680.035234.18235.9680.058334.44830.2900.043806.218
46.5670.05264.72241.2790.089146.86535.4220.085363.979
58.2140.07037.66146.8600.11861.85140.8700.128156.466
74.6650.08815.81852.8070.14750.08748.5200.17262.789
87.5940.1068.00058.4040.17737.28755.9470.21542.443
101.2260.1153.89966.6980.20624.07863.0470.25829.190
116.1860.1281.77174.7220.23915.77172.7850.30017.467
126.3570.1421.03675.9930.25414.74976.1540.32214.624
142.3410.1500.44686.1140.2698.65099.9000.3434.181
155.3740.1600.224135.2650.2800.648187.9960.4270.040
145.3000.2950.382

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.


Data resourceVG equationGardner equationFX equation

Sand (filter paper)0.2470.0660.7350.1580.0080.9480.1240.0070.943
Silt (filer paper)0.5610.2100.6260.3550.1830.4840.1870.0370.800
Clay (filter paper)1.0120.4680.5380.9420.4070.5680.3540.1920.459
Sand [46]0.7580.3080.5930.7900.1460.8150.1290.1170.097
Silt [47]0.6390.1240.8060.4440.2970.3310.0150.189−11.694
Clay [51]0.7200.1720.7610.6120.1680.7250.3340.1570.529

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.

Acknowledgments

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).

References

  1. X. G. Zhong, X. Zeng, and J. G. Rose, “Shear modulus and damping ratio of rubber-modified asphalt mixes and unsaturated subgrade soils,” Journal of Materials in Civil Engineering, vol. 14, no. 6, pp. 496–502, 2002. View at: Publisher Site | Google Scholar
  2. Z. D. Sun, Z. X. Li, H. Q. Song et al., “Effect of soil-water characteristic curve on runoff in sponge city construction,” China Water and Wastewater, vol. 36, no. 9, pp. 117–122, 2020. View at: Google Scholar
  3. Y. Mualem and G. Dagan, “Hydraulic conductivity of soils: unified approach to the statistical models,” Soil Science Society of America Journal, vol. 42, no. 3, pp. 392–395, 1978. View at: Publisher Site | Google Scholar
  4. B. Ghanbarian, V. Taslimitehrani, G. Dong, and Y. A. Pachepsky, “Sample dimensions effect on prediction of soil water retention curve and saturated hydraulic conductivity,” Journal of Hydrology, vol. 528, pp. 127–137, 2015. View at: Publisher Site | Google Scholar
  5. D. G. Fredlund, N. R. Morgenstern, and R. A. Widger, “The shear strength of unsaturated soils,” Canadian Geotechnical Journal, vol. 15, no. 3, pp. 313–321, 1978. View at: Publisher Site | Google Scholar
  6. P. C. Lim, S. L. Barbour, and D. G. Fredlund, “The influence of degree of saturation on the coefficient of aqueous diffusion,” Canadian Geotechnical Journal, vol. 35, no. 5, pp. 811–827, 1998. View at: Publisher Site | Google Scholar
  7. R. Solone, M. Bittelli, F. Tomei, and F. Morari, “Errors in water retention curves determined with pressure plates: effects on the soil water balance,” Journal of Hydrology, vol. 470-471, pp. 65–74, 2012. View at: Publisher Site | Google Scholar
  8. B. Gong, H. Wu, and B. Wang, “Influence of stress states on soil-water characteristics of expansive soils,” Rock and Soil Mechanics, vol. 25, no. 12, pp. 1915–1918, 2004. View at: Google Scholar
  9. S. S. Agus, E. C. Leong, and H. Rahardjo, “Soil-water characteristic curves of Singapore residual soils,” Unsaturated Soil Concepts and Their Application in Geotechnical Practice, vol. 19, no. 3, pp. 285–309, 2001. View at: Publisher Site | Google Scholar
  10. P. K. Patrick, H. W. Olsen, and J. D. Higgins, “Comparison of chilled-mirror measurements and filter paper estimates of total soil suction,” Geotechnical Testing Journal, vol. 30, no. 5, pp. 360–367, 2007. View at: Publisher Site | Google Scholar
  11. D. G. Fredlund, “Unsaturated soil mechanics in engineering practice,” Journal of Geotechnical and Geoenvironmental Engineering, vol. 132, no. 3, pp. 286–321, 2006. View at: Publisher Site | Google Scholar
  12. J. H. M. Wösten, Y. A. Pachepsky, and W. J. Rawls, “Pedotransfer functions: bridging the gap between available basic soil data and missing soil hydraulic characteristics,” Journal of Hydrology, vol. 251, no. 3-4, pp. 123–150, 2001. View at: Publisher Site | Google Scholar
  13. L. M. Arya and J. F. Paris, “A physicoempirical model to predict the soil moisture characteristic from particle-size distribution and bulk density data,” Soil Science Society of America Journal, vol. 45, no. 6, pp. 1023–1030, 1981. View at: Publisher Site | Google Scholar
  14. S. Rafraf, L. Guellouz, H. Guiras, and R. Bouhlila, “A new model using dynamic contact angle to predict hysteretic soil water retention curve,” Soil Science Society of America Journal, vol. 80, no. 6, pp. 1433–1442, 2016. View at: Publisher Site | Google Scholar
  15. S. W. Tyler and S. W. Wheatcraft, “Application of fractal mathematics to soil water retention estimation,” Soil Science Society of America Journal, vol. 53, no. 4, pp. 987–996, 1989. View at: Publisher Site | Google Scholar
  16. B. Ghanbarian-alavijeh, A. Liaghat, G.-H. Huang, and M. T. Van Genuchten, “Estimation of the van Genuchten soil water retention properties from soil textural data,” Pedosphere, vol. 20, no. 4, pp. 456–465, 2010. View at: Publisher Site | Google Scholar
  17. M. T. van Genuchten, “A closed-form equation for predicting the hydraulic conductivity of unsaturated soils,” Soil Science Society of America Journal, vol. 44, no. 5, pp. 892–898, 1980. View at: Publisher Site | Google Scholar
  18. Y. W. Zhang, Z. P. Song, X. L. Weng, and Y. L. Xie, “A new soil-water characteristic curve model for unsaturated loess based on wetting-induced pore deformation,” Geofluids, vol. 2019, Article ID 1672418, 14 pages, 2019. View at: Publisher Site | Google Scholar
  19. G. L. Tao, Y. Y. Chen, H. L. Xiao, Y. Chen, and W. Peng, “Comparative analysis of soil-water characteristic curve in fractal and empirical models,” Advances in Materials Science and Engineering, vol. 2020, Article ID 1970314, 10 pages, 2020. View at: Publisher Site | Google Scholar
  20. K. Pham, D. Kim, Y. Yoon, and H. Choi, “Analysis of neural network based pedotransfer function for predicting soil water characteristic curve,” Geoderma, vol. 351, pp. 92–102, 2019. View at: Publisher Site | Google Scholar
  21. L. Wang, W. G. Zhang, and F. Y. Chen, “Bayesian approach for predicting soil-water characteristic curve from particle-size distribution data,” Energies, vol. 12, no. 15, 2019. View at: Publisher Site | Google Scholar
  22. C. P. K. Gallage and T. Uchimura, “Effects of dry density and grain size distribution on soil-water characteristic curves of sandy soils,” Soils and Foundations, vol. 50, no. 1, pp. 161–172, 2010. View at: Publisher Site | Google Scholar
  23. L. Zhang and Q. Chen, “Predicting bimodal soil-water characteristic curves,” Journal of Geotechnical and Geoenvironmental Engineering, vol. 131, 2005. View at: Publisher Site | Google Scholar
  24. Q. Zhai, H. Rahardjo, and A. Satyanaga, “Estimation of air permeability function from soil-water characteristic curve,” Canadian Geotechnical Journal, vol. 56, no. 4, pp. 505–513, 2019. View at: Publisher Site | Google Scholar
  25. Y. Perera, C. Zapata, W. Houston, and S. Houston, “Prediction of the soil-water characteristic curve based on grain-size-distribution and index properties,” in Proceedings of the Geo-Frontiers Congress 2005, Austin, TX, USA, October 2005. View at: Publisher Site | Google Scholar
  26. D. G. Fredlund and A. Xing, “Equations for the soil-water characteristic curve,” Canadian Geotechnical Journal, vol. 31, no. 4, pp. 521–532, 1994. View at: Publisher Site | Google Scholar
  27. N. Lu and Y. Dong, “Correlation between soil-shrinkage curve and water-retention characteristics,” Journal of Geotechnical and Geoenvironmental Engineering, vol. 143, no. 9, 2017. View at: Publisher Site | Google Scholar
  28. B. Lin and A. B. Cerato, “Hysteretic soil water characteristics and cyclic swell-shrink paths of compacted expansive soils,” Bulletin of Engineering Geology and the Environment, vol. 72, no. 1, pp. 61–70, 2013. View at: Publisher Site | Google Scholar
  29. Q. Zhai, H. Rahardjo, A. Satyanaga, G. Dai, and Y. Zhuang, “Framework to estimate the soil-water characteristic curve for soils with different void ratios,” Bulletin of Engineering Geology and the Environment, vol. 79, no. 8, pp. 4399–4409, 2020. View at: Publisher Site | Google Scholar
  30. W. R. Gardner, “Some steady-state solutions of the unsaturated moisture flow equation with application to evaporation from a water table,” Soil Science, vol. 85, no. 4, pp. 228–232, 1958. View at: Publisher Site | Google Scholar
  31. D. G. Fredlund, A. Xing, M. D. Fredlund, and S. L. Barbour, “The relationship of the unsaturated soil shear strength to the soil-water characteristic curve,” Canadian Geotechnical Journal, vol. 33, no. 3, pp. 440–448, 1996. View at: Publisher Site | Google Scholar
  32. R. H. Brooks and A. T. Corey, “Hydraulic properties of porous media,” Hydrology Paper, vol. 7, pp. 26–28, 1964. View at: Publisher Site | Google Scholar
  33. M. Wijaya and E. C. Leong, “Equation for unimodal and bimodal soil-water characteristic curves,” Soils and Foundations, vol. 56, no. 2, pp. 291–300, 2016. View at: Publisher Site | Google Scholar
  34. M. Aubertin, J.-F. Ricard, and R. P. Chapuis, “A predictive model for the water retention curve: application to tailings from hard-rock mines,” Canadian Geotechnical Journal, vol. 35, no. 1, pp. 55–69, 1998. View at: Publisher Site | Google Scholar
  35. W. S. Sillers and D. G. Fredlund, “Statistical assessment of soil-water characteristic curve models for geotechnical engineering,” Canadian Geotechnical Journal, vol. 38, no. 6, pp. 1297–1313, 2001. View at: Publisher Site | Google Scholar
  36. L. Wang, Z.-J. Cao, D.-Q. Li, K.-K. Phoon, and S.-K. Au, “Determination of site-specific soil-water characteristic curve from a limited number of test data - a bayesian perspective,” Geoscience Frontiers, vol. 9, no. 6, pp. 1665–1677, 2018. View at: Publisher Site | Google Scholar
  37. X. Ren, J. Kang, J. Ren, X. Chen, and M. Zhang, “A method for estimating soil water characteristic curve with limited experimental data,” Geoderma, vol. 360, Article ID 114013, 2020. View at: Publisher Site | Google Scholar
  38. M. Aubertin, M. Mbonimpa, B. Bussière, and R. P. Chapuis, “A model to predict the water retention curve from basic geotechnical properties,” Canadian Geotechnical Journal, vol. 40, no. 6, pp. 1104–1122, 2003. View at: Publisher Site | Google Scholar
  39. S. K. Vanapalli, D. G. Fredlund, D. E. Pufahl, and A. W. Clifton, “Model for the prediction of shear strength with respect to soil suction,” Canadian Geotechnical Journal, vol. 33, no. 3, pp. 379–392, 1996. View at: Publisher Site | Google Scholar
  40. X. Zhang, “Research on the principle problem relating to soil-water characteristic curve and its application in unsaturated soil mechanics,” Beijng Jiaotong University, Beijing, China, 2010, Doctoral Dissertation. View at: Google Scholar
  41. S. Assouline, D. Tessier, and A. Bruand, “A conceptual model of the soil water retention curve,” Water Resources Research, vol. 34, no. 2, pp. 223–231, 1998. View at: Publisher Site | Google Scholar
  42. X. Li, L. M. Zhang, and D. G. Fredlund, “Wetting front advancing column test for measuring unsaturated hydraulic conductivity,” Canadian Geotechnical Journal, vol. 46, no. 12, pp. 1431–1445, 2009. View at: Publisher Site | Google Scholar
  43. C. M. P. Vaz, M. de Freitas Iossi, J. de Mendonça Naime et al., “Validation of the arya and paris water retention model for Brazilian soils,” Soil Science Society of America Journal, vol. 69, no. 3, pp. 577–583, 2005. View at: Publisher Site | Google Scholar
  44. C. Chen, K. Hu, W. Li, G. Wang, and G. Liu, “Estimating the wet-end section of soil water retention curve by using the dry-end section,” Soil Science Society of America Journal, vol. 78, no. 6, pp. 1878–1883, 2014. View at: Publisher Site | Google Scholar
  45. A. G. Hunt and G. W. Gee, “Water-Retention of fractal soil models using continuum percolation theory: tests of hanford site soils,” Vadose Zone Journal, vol. 1, no. 2, pp. 252–260, 2002. View at: Publisher Site | Google Scholar
  46. M. H. Mohammadi and M. Vanclooster, “Predicting the soil moisture characteristic curve from particle size distribution with a simple conceptual model,” Vadose Zone Journal, vol. 10, no. 2, pp. 594–602, 2011. View at: Publisher Site | Google Scholar
  47. D. K. Jensen, M. Tuller, L. W. de Jonge, E. Arthur, and P. Moldrup, “A new two-stage approach to predicting the soil water characteristic from saturation to oven-dryness,” Journal of Hydrology, vol. 521, pp. 498–507, 2015. View at: Publisher Site | Google Scholar
  48. W. M. Cornelis, J. Ronsyn, M. Van Meirvenne, and R. Hartmann, “Evaluation of pedotransfer functions for predicting the soil moisture retention curve,” Soil Science Society of America Journal, vol. 65, no. 3, pp. 638–648, 2001. View at: Publisher Site | Google Scholar
  49. L. A. George, M. M. Dewoolkar, and D. Znidarcic, “Simultaneous laboratory measurement of acoustic and hydraulic properties of unsaturated soils,” Vadose Zone Journal, vol. 8, no. 3, pp. 633–642, 2009. View at: Publisher Site | Google Scholar
  50. Q. Y. Mu, H. Dong, H. J. Liao, Y. J. Dang, and C. Zhou, “Water-retention curves of loess under wetting−drying cycles,” Géotechnique Letters, vol. 10, no. 2, pp. 135–140, 2020. View at: Publisher Site | Google Scholar
  51. M. A. Dafalla, A. M. Al-Mahbashi, A. Almajed, and M. Al-Shamrani, “Predicting soil-water characteristic curves of clayey sand soils using area computation,” Mathematical Problems in Engineering, vol. 2020, Article ID 4548912, 9 pages, 2020. View at: Publisher Site | Google Scholar
  52. N. R. A. Bird, E. Perrier, and M. Rieu, “The water retention function for a model of soil structure with pore and solid fractal distributions,” European Journal of Soil Science, vol. 51, no. 1, pp. 55–63, 2008. View at: Publisher Site | Google Scholar
  53. K. M. A. Al Haj and J. R. Standing, “Soil water retention curves representing two tropical clay soils from Sudan,” Géotechnique, vol. 66, no. 1, pp. 71–84, 2016. View at: Publisher Site | Google Scholar
  54. E. Romero, G. Della Vecchia, and C. Jommi, “An insight into the water retention properties of compacted clayey soils,” Géotechnique, vol. 61, no. 4, pp. 313–328, 2011. View at: Publisher Site | Google Scholar
  55. D. Sun, Z. Zhang, Y. Gao, and B. Chen, “Determination of basic parameters of swcc for clayer soils,” Journal of Shanghai University (Natural Science), vol. 25, no. 6, pp. 957–964, 2019. View at: Google Scholar
  56. E. C. Leong and H. Rahardjo, “Review of soil-water characteristic curve equations,” Journal of Geotechnical and Geoenvironmental Engineering, vol. 123, no. 12, pp. 1106–1117, 1997. View at: Publisher Site | Google Scholar
  57. M. D. Fredlund, G. W. Wilson, and D. G. Fredlund, “Use of the grain-size distribution for estimation of the soil-water characteristic curve,” Canadian Geotechnical Journal, vol. 39, no. 5, pp. 1103–1117, 2002. View at: Publisher Site | Google Scholar
  58. X. W. Ren and J. C. Santamarina, “The hydraulic conductivity of sediments: a pore size perspective,” Engineering Geology, vol. 233, pp. 48–54, 2018. View at: Publisher Site | Google Scholar
  59. X. Ren, Y. Zhao, Q. Deng, J. Kang, D. Li, and D. Wang, “A relation of hydraulic conductivity - void ratio for soils based on kozeny-carman equation,” Engineering Geology, vol. 213, pp. 89–97, 2016. View at: Publisher Site | Google Scholar
  60. H. Rahardjo, A. Satyanaga, H. Mohamed, S. C. Yee Ip, and R. S. Shah, “Comparison of soil-water characteristic curves from conventional testing and combination of small-scale centrifuge and dew point methods,” Geotechnical and Geological Engineering, vol. 37, no. 2, pp. 659–672, 2019. View at: Publisher Site | Google Scholar
  61. R. Gardner, “A method of measuring the capillary tension of soil moisture over a wide moisture range,” Soil Science, vol. 43, no. 4, pp. 277–284, 1937. View at: Publisher Site | Google Scholar
  62. C. Feuerharmel, W. Y. Y. Gehling, and A. V. D. Bica, “The use of filter-paper and suction-plate methods for determining the soil-water characteristic curve of undisturbed colluvium soils,” Geotechnical Testing Journal, vol. 29, no. 5, pp. 419–425, 2006. View at: Publisher Site | Google Scholar
  63. S. S. Agus and T. Schanz, “Comparison of four methods for measuring total suction,” Vadose Zone Journal, vol. 4, no. 4, pp. 1087–1095, 2005. View at: Publisher Site | Google Scholar
  64. F. A. M. Marinho and O. M. Oliveira, “The filter paper method revisited,” Geotechnical Testing Journal, vol. 29, no. 3, pp. 250–258, 2006. View at: Publisher Site | Google Scholar
  65. R. Fawcett and N. Collis-George, “A filter-paper method for determining the moisture characteristics of soil,” Australian Journal of Experimental Agriculture, vol. 7, no. 25, pp. 162–167, 1967. View at: Publisher Site | Google Scholar
  66. D. Tang, D. Q. Li, H. B. Kim, and E. Ganju, “Research on calibration curves of home-made “double circle” filter papers,” Engineering Journal of Wuhan University, vol. 49, no. 1, pp. 1–8, 2016. View at: Google Scholar
  67. F. Bai, S. Liu, and J. Yuan, “Calibration test for total suction wetting curve of filter paper,” Rock and Soil Mechanics, vol. 32, no. 8, pp. 2336–2340, 2011. View at: Google Scholar
  68. Z. Wang, J. X. Yang, J. J. Kuang, J. Y. An, and Y. D. Luo, “Application of filter paper method in field measurement of matric suction,” Chinese Journal of Geotechnical Engineering, vol. 25, no. 4, pp. 405–408, 2003. View at: Google Scholar
  69. S. L. Houston, W. N. Houston, and A. M. Wagner, “Laboratory filter-paper suction measurements,” Geotechnical Testing Journal, vol. 17, no. 2, pp. 185–194, 1994. View at: Publisher Site | Google Scholar
  70. E. C. Leong, L. He, and H. Rahardjo, “Factors affecting the filter paper method for total and matric suction measurements,” Geotechnical Testing Journal, vol. 25, no. 3, pp. 322–333, 2002. View at: Publisher Site | Google Scholar
  71. J. A. Munoz-Castelblanco, J. M. Pereira, P. Delage, and Y. J. Cui, “The water retention properties of a natural unsaturated loess from northern France,” Géotechnique, vol. 62, no. 2, pp. 95–106, 2012. View at: Publisher Site | Google Scholar

Copyright © 2021 Jie Zhou et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Related articles

No related content is available yet for this article.
 PDF Download Citation Citation
 Download other formatsMore
 Order printed copiesOrder
Views378
Downloads354
Citations

Related articles

No related content is available yet for this article.

Article of the Year Award: Outstanding research contributions of 2021, as selected by our Chief Editors. Read the winning articles.