Suitability of Soil Water Retention Characteristic Models (SWRC) in Regions and Soil Depth

DOI: 10.4236/jwarp.2019.116044   PDF   HTML   XML   336 Downloads   675 Views  


Soil Water Retention Characteristics (SWRC) models have been widely used in many applications. Presently, there are many models in the literature and many more still being developed so much so that it is confusing which model to prefer. The current choice of the appropriate model to use has not been well guided by any incisive research on the predictive performance of these models. Consequently, SWRC model applications have been largely moved by convenience. This study used a global dataset to evaluate 12 commonly used SWRC models. The measured data onto which the models were evaluated was grouped into different soil depths and different regions of the world. The evaluation used correlation, Nash-Sutcliffe efficiency, and residual standard error statistics to choose the best overall performing model and models for each category. It gives an indication of the type of SWRC models to use in different regions of the world and depths of sampling. The suitability of the models to regions showed that the Fredlund and Xing model had the best performance in subsoils in Africa; Omuto in Southern Asia; and van Genuchten in subsoils of the other regions. It is recommended that many more models be tested using the procedures in this study so that benchmarks can be established on SWRC model selection suitable for various regions.

Share and Cite:

Too, V. , Omuto, C. , Biamah, E. and Obiero, J. (2019) Suitability of Soil Water Retention Characteristic Models (SWRC) in Regions and Soil Depth. Journal of Water Resource and Protection, 11, 740-747. doi: 10.4236/jwarp.2019.116044.

1. Introduction

Water is held in the soil by gravity, matric suction, and osmotic forces. The models describing different levels of these forces and the corresponding amount of water held in the soil are known as soil water retention characteristics (SWRC) models. There are two major categories of these models in the literature: pedotransfer functions (PTFs) which relate easily measurable soil properties to SWRC [1] and mathematical functions that can be fitted to experimental SWRC [2] . Predictive models under these two categories have been widely used in the literature in areas such as irrigation water management, crop-yield estimation, global circulation modeling, land degradation assessment, geotechnical engineering, etc. [3] [4] [5] [6] [7] . The major difference between these two model categories is that the mathematical functions contain fitting parameters which can be related to soil physical properties while PTFs do not necessarily have this relationship. Consequently, the mathematical functions have been used as a fundamental input in the various applications using SWRC. PTFs have been used largely as an alternate method to derive the parameters of the mathematical functions. This explains why there is still an active research to produce better user-friendly and accurate mathematical models for SWRC than the existing models [8] [9] [10] .

There are many documented mathematical models for SWRC. Reviews of literature on these models show that the popular ones are those that contain three-, four- and five-parameters [11] [12] . For example, a review by Leong and Rahardjo [11] showed that the popular models are Brooks-Corey, Farrel-Larson, Fredlund-Xing, Gardner, McKee-Bumb, van Genuchten, and Williams SWRC models. Another review by Bullied et al. [13] portrayed Campbell, Russo and Tani models as popular in addition to some of the models contained in the review by Leong and Rhardjo [11] . Another aspect of these reviews shows that the popular models mostly have high success rates in areas where they were developed. It is important to note that they were developed and tested in different regions of the world. Consequently, some of them have been successful in those areas, others have been modified to extend their applications in independent areas, while some have been used in different places without regard to whether they are successful or not. Still, many more are being developed in an attempt to improve on the existing models. However, there is no exhaustive study to compare their performance in different regions of the world. There is a need to determine the relative performance of these models in order to guide their application or encourage further development of new models.

Since the advent of digital soil mapping (DSM) and digital soil assessment (DSA) paradigms, there has been an upsurge of development of soil inference models and the need to increase accurate application of soil mapping products [14] . Soil hydraulic parameters and hydrologic functions cannot be left behind in this regard. There is a need to improve accuracy of their assessment and application within the realms of DSM and DSA. This present study contributed towards this goal by testing the performance of commonly used SWRC models on a global dataset. The focus on the use of a global dataset was to establish regions where these models can register the best performance and consequently inform users of these models where they can guarantee a high level of certainty.

2. Materials and Methods

2.1. SWRC Models Analyzed

The SWRC models analyzed in this study and their curve-fitting parameters are given in Table 1. The models were also placed into three groups: five-parameter, four-parameter, and three-parameter SWRC models. This grouping was necessary because models with many parameters may tend to perform better than those with few parameters and lead to biased performance comparison.

The five-parameter models were those proposed by Fredlund and Xing [2] , Omuto [15] and van Genutchen [16] . The van Genuchten [16] model is that in which n and m parameters are independent. The four-parameter models were those proposed by Kosugi [8] , Dexter et al. [9] , van Genuchten [16] (in which the parameter m = 1 − 1/n), Gardner [17] , Brooks and Corey [18] , and Russo [19] . The three-parameter models were those proposed by Campbell [20] , Tani [21] , and McKee and Bumb [22] .

2.2. Input Data

The global data used for analysis of the SWRC models was downloaded on 6th May 2013. The data contained soil water content (θ in cm3∙cm−3) at various levels

Table 1. SWRC models tested.

of suction potential (h in cm). Only samples with at least eight (8) levels of suction potential and water contents were considered. The samples were grouped into topsoil (at depths between 0 cm and 50 cm from the soil surface) and subsoil (at depths greater than 50 cm from the soil surface). This grouping of samples into various depths was mainly to provide suitability of the models to various uses like irrigation modeling, watershed management and geotechnical investigations since soil depths are considered. These two groups were further split into samples from Central America, South America, Africa, Europe, Middle East, North Asia, and South Asia (Figure 1).

Data description and associated measurement methods have been given in Batjes [23] . There were at least nine (9) samples with complete SWRC data (i.e. with eight suction potential levels) for topsoil and at least six (6) samples with complete SWRC data for subsoils in each region. Preliminary analysis of the data showed that the samples had at least one measurement around saturation (i.e. when the suction potential is very low), at least one inflection point, and at least one measurement around the dry end (i.e. when the suction potential is very high). In addition, the majority of the samples had an inflection point around −1.0 m suction potential.

2.3. Data Analysis

The models in Table 1 were fitted to measured SWRC data using nonlinear curve fitting method. The analysis of these models in predicting measured SWRC for various regions and depths were analyzed using the following popular statistics:

1) Residual standard error, RSE

i = 1 n ( y i y ^ i ) 2 / ( n 2 ) (1)

2) Nash-Sutcliffe model efficiency, EF

1 [ i = 1 n ( y i y ^ i ) 2 / i = 1 n ( y i y ¯ i ) 2 ] (2)

3) Correlation, r2

Figure 1. Spatial distribution of samples for tested SWRC data.

( n i = 1 n y i y ^ i ( y i ) ( y ^ i ) ) 2 / [ ( n ( y i 2 ) ( y ) 2 ) ( n ( y ^ i 2 ) ( y ^ i ) 2 ) ] (3)

where yi is the measured soil water contents, y ^ i are the fitted soil water content, y ¯ i is the mean measured soil water contents, and n are the number of data points (i.e. soil moisture levels for each soil sample). In analysis, values of EF and r2 close to unity and low RSE are preferred for a good model.

3. SWRC Models Suitability in Various Regions and Depths

The suitability of SWRC models in various regions was done by comparing their performance with data from different regions of the world and depths of soil sample below the soil surface. Using correlation statistics, the performance of the five-parameter models using the subsoils data showed that the Fredlund and Xing model had the best performance (R2 = 0.9964) in Africa; Omuto in Southern Asia (R2 = 0.9987) and van Genuchten in Central America, South America, Europe and Northern Asiaregions (Table 2). In the four-parameter category, The van Genuchten, four-parameter model, did well among the other four-parameter models in subsoil Africa, America, Europe, Asia, M. East regions with R2 of 0.9931, 0.9973, 0.9922, 0.9935, 0.9899, 0.9903 and 0.9977 respectively.

Using topsoil samples, in the five-parameter category, Omutomodel performed well in North and South Asiaregions with R2 of 0.9991 and 0.9986 respectively, Fredlund and Xing model did well in Africa and south America regions with R2 of 0.9960 and 0.9976 respectively. The five-parameter van Genuchtenmodel performed well in Central America and Europe with R2 of 0.9948 and 0.9971 while the four-parameter van Genuchtenmodel was found to perform well in all the other regions except in Central America and Europe on the topsoil samples. The three-parameter exponential model proposed by McKee and Bumb outperformed the other models in the three-parameter category in all regions and soil depths. This performance pattern was also portrayed by the other statistical indices (i.e. RSE and EF).

The group-results shown in Table 2 may be used to guide models selection in case available SWRC data has limited suction potential levels. For example in Africa, if SWRC data has four measured suction potential levels then the exponential model by McBee and Bumb would give reliable results; in the case of SWRC data having five measured suction levels then the four-parameter van Genuchten model should be preferred; and in the case of SWRC data having at least six measured suction potential levels then Fredlund and Xing model should be preferred. It should be noted that, in order to guarantee best predictive results with these models, the minimum data-points for the SWRC data should contain at least three important points: at least a point around the saturation, at least a point at any of the inflection points (e.g. at air-entry potential), and a point around the dry end. These are anchoring points which are important in SWRC model fitting.

Table 2. Summary statistics of SWRC models performance.

Best performing model in the category.

4. Conclusion

This study used a global measured dataset of SWRC to evaluate popular models for fitting SWRC. The analysis serves as a guide for selecting the models to be preferred for fitting SWRC models to use in different regions of the world and depths of sampling. Van Genuchten model (both the four and five parameter models) had the best performance in most of the regions and soil depth. This was notable in Africa, Central America, South America, Europe and Northern Asia regions. The three-parameter McKee and Bumb model also performed well in all regions and soil depths. Owing to the good performance of exponential-based models, the analysis also floated a suggestion for future models to develop models based on exponential pore-size distribution. It is recommended that many more models be tested using the procedures used in this study so that benchmarks can be established on SWRC model selection for different applications.


The soil data used in this study were obtained from The International Soil Reference and Information Centre (ISRIC). The organization and countries which contributed to this dataset are highly acknowledged.

Conflicts of Interest

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


[1] Minasny, B. and McBratney, A.B. (2007) Estimating the Water Retention Shape Parameter from Sand and Clay Content. Soil Science Society of America Journal, 71, 1105-1110.
[2] Fredlund, D.G. and Xing, A. (1994) Estimation of Soil Water Characteristic Curve. Canadian Geotechnical Journal, 31, 521-532.
[3] Vanapalli, S.K., Fredlund, D.G., Pufahl, D.E. and Clifton, A.W. (1996) Model for Prediction of Shear Strength with Respect to Matric Suction. Canadian Geotechnical Journal, 33, 379-392.
[4] Liu, J., Goering, C.E. and Tian, L. (2001) Neural Network for Setting Target Corn Yields. Transactions of American Society of Agricultural Engineers, 44, 705-713.
[5] Dexter, A.R. (2004) Soil Physical Quality. Part I. Theory, Effects of Soil Texture, Density, and Organic Matter, and Effects on Root Growth. Geoderma, 120, 201-214.
[6] Chotpantarat, S., Limpakanwech, C., Siriwong, W., Siripattankul, S. and Sutthirat, C. (2011) Effects of Soil Water Characteristic Curves on Simulation of Nitrate Vertical Transport in a Thai Agricultural Soil. Sustainable Environmental Research, 21, 187-193.
[7] Patil, N.G., Pal, K.D., Mandal, C. and Mandal, D.K. (2012) Soil Water Retention Characteristics of Vertisols and Pedotransfer Functions Based on Nearest Neighbor and Neural Networks Approaches to Estimate AWC. Journal of Irrigation and Drainage Engineering, 138, 177-184.
[8] Kosugi, K. (1996) Lognormal Distribution Model for Unsaturated Soil Hydraulic Properties. Water Resources Research, 32, 2697-2703.
[9] Dexter, A.R., Czyz, E.A., Richard, G. and Reszkwoska, A. (2008) A User-Friendly Water Retention Function That Takes Account of the Textural and Structural Pore Spaces in Soil. Geoderma, 143, 243-253.
[10] Omuto, C.T. and Gumbe, L.O. (2009) Estimating Water Infiltration and Retention Characteristics Using a Computer Program in R. Computers and Geosciences, 35, 579-585.
[11] Leong, E.C. and Rahardjo, H. (1997) review of Soil Water Characteristic Curve Equations. Journal of Geotechnical and Geo-Environmental Engineering, 123, 1106-1117.
[12] Malaya, C. and Sreedeep, S. (2012) Critical Review on the Parameters Influencing Soil-Water Characteristic Curve. Journal of Irrigation and Drainage Engineering, 138, 55-62.
[13] Bullied, W.J., Bullock, R.P. and van Acker, R.C. (2012) Modelling Soil Water Retention for Weed Seed Germination Sensitivity to Water Potential. Applied Environmental and Soil Science, 2012, Article ID: 812561.
[14] Minasny, B., Malone, B.P. and McBratney, A.B. (2012) Digital Soil Assessments and Beyond: Proceedings of the 5th Global Workshop on Digital Soil Mapping 2012, Sydney, Australia. CRC Press, Boca Raton.
[15] Omuto, C.T. (2009) Biexponential Model for Water Retention Characteristics. Geoderma, 149, 235-242.
[16] van Genuchten, M.T. (1980) A Closed-Form Equation for Predicting the Hydraulic Conductivity of Unsaturated Soils. Soil Science Society of America Journal, 44, 892-898.
[17] Gardner, W.R. (1958) Some Steady State Solutions of the Unsaturated Moisture Flow Equation with Application to Evaporation from a Water Table. Soil Science, 85, 228-232.
[18] Brooks, R.H. and Corey, A.T. (1964) Hydraulic Properties of Porous Medium. Hydrology Paper Number 3. Colorado State University, Fort Collins.
[19] Russo, D. (1988) Determining Soil Hydraulic Properties by Parameter Estimation: On the Selection of a Model for the Hydraulic Properties. Water Resources Research, 24, 453-459.
[20] Campbell, G.S. (1974) A Simple Method for Determining Unsaturated Conductivity from Moisture Retention Data. Soil Science, 117, 311-314.
[21] Tani, M. (1982) The Properties of Water-Table Rise Produced by a One-Dimensional, Vertical, Unsaturated Flow. Journal of Japan Forestry Society, 64, 409-418.
[22] McKee, C.R. and Bumb, A.C. (1984) The Importance of Unsaturated Flow Parameters in Designing a Monitoring System for Hazardous Wastes and Environmental Emergencies. Proceedings of Hazardous Materials and Control Research Institute of Nature Conference, Houston, Tex, 50-58.
[23] Batjes, N.H. (1995) Modelling Soil Water Retention for Weed Seed Germination Sensitivity to Water Potential Version 1. ISRIC, Wageningen, 47.

comments powered by Disqus

Copyright © 2020 by authors and Scientific Research Publishing Inc.

Creative Commons License

This work and the related PDF file are licensed under a Creative Commons Attribution 4.0 International License.