Multiple Linear Regression (MLR), Artificial Neural Network (ANN) and Rosetta model were employed to develop pedotransfers functions (PTFs) for soil moisture prediction using available soil properties for northern soils of Iran. The Rosetta model is based on ANN works in a hierarchical approach to predict water retention curves. For this purpose, 240 soil samples were selected from the south of Guilan province, Gilevan region, northern Iran. The data set was divided into two subsets for calibration and testing of the models. The general performance of PTFs was evaluated using coefficient of determination (R2), root mean square error (RMSE) and mean biased error between the observed and predicted values. Results showed that ANN with two hidden layers, Tan-sigmoid and linear functions for hidden and output layers respectively, performed better than the others in predicting soil moisture. In the other hand, ANN can model non-linear functions and showed to perform better than MLR. After ANN, MLR had better accuracy than Rosetta. The developed PTFs resulted in more accurate estimation at matric potentials of 100, 300, 500, 1000, 1500 kPa. Whereas, Rosetta model resulted in slightly better estimation than derived PTFs at matric potentials of 33 kPa. This research can provide the scientific basis for the study of soil hydraulic properties and be helpful for the estimation of soil water retention in other places with similar conditions, too.
Additional key wordsmultiple linear regressionneural networkspedotransfer functionRosettasoil moisture curveAbbreviations usedANN (artificial neural network)CEC (cation exchange capacity)FWC (field water capacity)MBE (mean bias error)MLP (multi-layer perceptron)MLR (multiple linear regressions)PTFs (pedotransfer functions)PWP (permanent wilting point)RMSE (root mean square error)SWR (soil water retention)This research was supported by the University of Guilan, Iran.
Competing interests: The authors have declared that no competing interests exist.
Introduction
Soil hydrodynamic properties drive the flow of water in the soil-plant-atmosphere system, and hence control processes such as aquifer recharge or nutrient fluxes between soil and vegetation. Knowledge of soil hydrodynamics is important for modeling physical processes related to soil water content. Despite great advances in measurement methods, it is still difficult to determine soil hydraulic properties accurately, especially for undisturbed soils and in the dry range. However, the measurement of the soil hydraulic properties is time-consuming, labor-intensive and expensive (Merdun, 2010). If the area being evaluated is large enough to exhibit substantial spatial variability of soil water availability, it is costly and time consuming to perform adequate measurements to provide good estimates within the temporal and financial constraints of the project (Givi et al., 2004).
Recently, an alternative, indirect estimation of soil hydraulic properties from widely available or more easily measured basic soil properties using pedotransfer functions (PTFs) has attracted considerable attention of researchers in a variety of fields such as soil scientists, hydrologists, and agricultural and environmental engineers (Minasny et al., 2004; Huang et al., 2010).PTFs are based on physical approaches or on empirical regression equations that link soil physical and/or chemical characteristics, easily measured, to hydrodynamic properties (Bouma, 1989; Hamblin, 1991) that just exist as a narrow relationship among these soil properties.
PTFs are classified in three main groups, i.e., class PTFs (Baker & Ellison, 2008); point PTFs (Ghanbarian-Alavijeh & Millán, 2010); and parametric or function PTFs (Santra & Das, 2008). The classes PTFs are used to estimate an average value of soil hydraulic parameters for each class of soil texture. These PTFs are inexpensively developed; however, their accuracy is less than point and parametric PTFs (Abbasi et al., 2011). Point PTFs functions estimate the water contents at fixed pressure heads, often including the water content at the field water capacity (FWC) and the permanent wilting point (PWP) or the water content at a given matric potential (Givi et al., 2004). Function PTFs predict the parameters of a closed-form analytical equation or empirical parameters of water retention curve models, such as the model of Brooks & Corey (1964) or the van Genuchten equation (Van Genuchten, 1980).
At present, there are two common methods to develop PTFs for point and/or function estimations, which are the Multiple Linear Regression (MLR) method (Merdun et al., 2006) and the Artificial Neural Network (ANN) method (Minasny et al., 2004; Sarmadian & Mehrjardi, 2008). In recent years, the ANN method was used to improve the prediction accuracy of PTFs (Kumar et al., 2010; Rezaei-Arshad et al., 2013; Moghimi et al., 2014;Mukhlisin & Abd Rahman, 2014; Rezaeianzadeh et al., 2014).
An advantage of ANN, as compared to other methods, is that it does not require a priori models. The optimal, possibly nonlinear, relations which link the input data (bulk density, particle-size data, etc.) to output data (soil water retention, FWC, etc.) are obtained and implemented in an iterative calibration procedure (Schaap et al., 2001). As a result, ANN method can typically extract the maximum amount of information from the data. For these reasons, the ANN method had a wide application (Merdun et al., 2006; Yi et al., 2013), and many scholars implicate that the general prediction of ANN is better than the MLR method (Moghimi et al., 2014).
However, this method also has some significant disadvantages which must be taken into consideration (Minasny et al., 2004). First of all, the interpretation of ANN is often difficult and subjective, because the fitting with the transfer function is a black-box approach. In addition, as it is usually the case in optimization, the sets of optimized weighting factors are not mathematically unique because the likelihood of convergence is at the local minimum. Consequently, different initial weight values may yield various results that deviate from the global minimum. Nevertheless, the ANN method showed good prediction function, and it would play a greater role in the prediction of soil hydraulic properties with the improvement of scientific technology.
Najafi & Givi (2006) used the ANN and PTF methods for prediction of soil bulk density. They pointed out that the ANNs are able to predict the soil bulk density better than the PTFs. Amini et al. (2005) estimated the cation exchange capacity (CEC) in the central part of Iran using soil organic matter and clay contents. They used the ANN and five experimental models that were on the basis of regression methods for their predictions. They showed that a neural network PTF with eight hidden neurons was able to predict CEC better than the regression PTFs. Also the ANN model significantly improved the accuracy of the prediction by up to 25%. They concluded that network models are in general more suitable for capturing the non-linearity of the relationship between variables. Jain &Kumar (2006) indicated that the ANN technique can be successfully employed for the purpose of calibration of infiltration equations. They had also found that the ANNs are capable of performing very well in situations of limited data availability. In contrast, Merdun (2010) pointed out that although the differences between regression and ANN models were not statistically significant, regression predicted point and parametric variables of soil hydraulic parameters better than ANN.
In order to make the PTFs as widely applicable as possible, Schaap et al. (2001) developed the Rosetta software by using a large number of soil hydraulic data and corresponding predictive soil properties. The Rosetta model based on ANNs works in a hierarchical approach which employs five different PTFs to predict water retention curves (Schaap et al., 2001).Rosetta software is used widely by many scientists for the estimation of soil water retention (SWR), saturated hydraulic conductivity and unsaturated hydraulic conductivity parameters (Schaap et al., 2001; Givi et al., 2004; Minasny et al., 2004; Stumpp et al., 2009).
Although there are many studies on developing and using PTFs as listed above, there is no universal method for the prediction of soil hydraulic parameters. Moreover, the existing PTFs for the estimation of soil hydraulic properties in the literature were not always applicable in other regions with acceptable accuracy (Cornelis et al., 2001; Wagner et al., 2001; Nemes et al., 2003). So, it is necessary to develop and evaluate the PTFs in different regions. Moreover, there are few studies comparing the performance of different methods simultaneously in the development of PTFs in the study area and its surrounding areas. Therefore, the objectives of this study were (1) to develop and validate the point PTFs using Multi-Layer Perceptron (MLP), MLR and function PTFs using the Rosetta software for the estimation of soil moisture in various matric potentials, and (2) to compare the predictive capabilities of the three methods using selected evaluation criteria in the Gilevan Region of Guilan Province, northern Iran.
Material and methodsStudy area and data collection
Study area is located in south of Guilan province, Gilevan region, northern Iran (36° 54´ 10” to 36° 50´ 00” N, 49° 02´ 30” to 49° 16´ 08” E) (Fig. 1). The climate is aridic. The annual precipitation is 245.1 ± 3 mm, and the average temperature is17.4 ± 2°C.
Location of study area in northern Iran.
Totally 240 samples were taken from 0–30 cm depth and air dried. Some clods were used to measure soil bulk density using clod method. Samples were passed through a 2mm sieve to determine particle-size distribution by the pipette method in combination with sieving method. The organic matter content was analyzed with the Walkly-Black method. Calcium carbonate was determined based on calcimetery method (Burt, 2004). Disturbed soil samples were used to measure soil water content at −10, −33, −100, −300 kPa matric potentials using pressure plate and at −500, −1000, and −1500 kPa using pressure membrane. The resulting curves were fitted to the Van Genuchten model using the RETC software package (Van Genuchten, 1980) that provides four parameters (θr, θs, α and n). The first step for using statistical methods is to study the normality of data. In order to know whether the data were normal, Kolmogorov-Smirnov test was used. The data sets were divided into two subsets, randomly. One subset that includes 80% of the total data was used for calibrating selected PTFs, and the other subset was used for testing the calibrated models includes 20% of total data.
Multiple linear regression method (MLR)
The most common method used in point PTFs is to employ MLR. Gupta & Larson (1979) used the MLR method to estimate the SWR characteristics according to the data, such as bulk density, organic matter and so on. This method also became one of the most popular methods. Therefore, based on this model, the following function was used to develop the PTFs in this study (Eq. [1]):
where, θp is the soil water content at specific matric potentials; a0 is the regression constant; a1, a2, a3, a4, a5 are the regression coefficients; OC and BD represent the organic carbon and bulk density, respectively.
Artificial neural networks (ANN) model - MLP
Neural networks consist of a large class of different architectures. In many cases, the issue is approximating a static nonlinear, mapping fx with a neural network fNNx, where x∈RK. The most useful neural network in function approximation is Multi-Layer Perceptron (MLP) network. A MLP consists of an input layer, several hidden layers, and an output layer. Node i, also called a neuron, in a MLP network is shown in Suppl. Fig. S1 [pdf online].
It includes a summer and a nonlinear activation function g. The inputs, xk, k = 1,…,K to the neuron are multiplied by weights wki and summed up together with the constant bias term θi. The resulting ni is the input to the activation function g. The activation function was originally chosen to be a relay function, but for mathematical convenience a hyberbolic tangent (tan h) or a sigmoid function are most commonly used. Hyberbolic tangent is defined as:
The output of node i become:
Connecting several nodes in parallel and series, a MLP network is formed (Heikki, 2008). In MLP, the weighted sum of the inputs and bias term are passed to activation level through a transfer function to produce the output, and the units are arranged in a layered feed-forward topology called Feed Forward Neural Network (Venkatesan & Anitha, 2006).The three-layer neural network was used in this study: input, hidden and output layer. This type of network generally provides better performances in comparison to other types. The total number of nodes in the input and output layers coincides with the number of input and output variables in the data set. The ideal number of nodes in the hidden layer has to be found through trial-and-error method (Chang et al., 2004). The number of neurons in the hidden layer is of great importance, as too many neurons may cause over-fitting problems (Huang & Foo, 2002). For this purpose, variables that were been finally selected in MLR for developing PTFs, were used as inputs in ANN. Also, we used 1-10 nodes, and the ideal node number selected based on the least RMSE value. Each node carries a weight. By changing these weights, the input–output relation can be simulated. Also, we used Tan-sigmoid (non-linear) and linear activation function for the hidden and output layers, respectively. For ANN modeling, the computer software MatLab and the Neural Network Toolbox were used (Demuth & Beale, 2004).
Rosetta method
To facilitate application of the PTFs, Schaap et al. (2001) developed “Rosetta”, a computer program that implements some of the models published by Schaap & Leij (2000). Rosetta is able to estimate the Van Genuchten water retention parameters (Van Genuchten, 1980) and saturated hydraulic conductivity (Ks) (Mualem, 1976) based on pore-size model. The retention function is given by Eq. [4]:
where θ(h) is the measured volumetric water content (cm3/cm3) at the suction h (cm, taken positive for increasing suctions). The parameters θs (cm3/cm3) and θr (cm3/cm3) are saturated and residual water contents respectively; α>0 (in cm−1) is related to the inverse of the air entry suction; and n>1 is a measure of the pore-size distribution (Van Genuchten, 1980). A hierarchical approach with limited or more extended sets of predictors was used to estimate the water retention parameters (θs, θr, α, n) and Ks (Schaap et al., 2001).
Evaluation criteria
Accuracy of the MLR and Rosetta methods for derivation of PTFs was evaluated by using the coefficient of determination (R2), root mean square error (RMSE) and mean biased error (MBE) between the measured and predicted values of a given hydraulic parameter. The R2, RMSE and MBE are expressed as:
where yi denotes the measured value, is the predicted value, is the average of the measured value, and n is the total number of observations. The MBE characterizes the mean difference between the calculated and measured data; hence, it is a criterion of systematic error in the model fitting. Negative and positive values of MBE indicate under and over estimation of PTFs for a given parameter respectively.
Results
In this paper, use of MLR, ANN and Rosetta models, for the prediction of soils water content, was described and compared. Descriptive statistics for soil properties are summarized in Table 1. The studied soils can be classified as clayey, clay loam, silty clay loam and sandy clay loam. In general, clay, silt and sand content were between 18% to 48%, 15% to 45.5%, and 14% to 57% respectively. Average bulk density was 1.42 g/cm3 with values varying between 1.28 to 1.6 g/cm3.
Descriptive statistics of variables used for the development of PTFs
Correlation coefficients among the water content in each potential and soil physical and chemical properties were calculated and are reported in Table 2. From the values of correlation coefficients in Table 2, the soil bulk density (BD), clay and sand content have close relationships with the water content at 10 and 33 kPa at the significance level consistent with other scholars (Nemes & Rawls, 2006) and clay and sand content are only suitable for PTFs developing to predict the water content of 100, 300 and 500 kPa. The water content at 1000 and 1500 kPa were only controlled by clay content. The absolute effect of BD decreased with hm in agreement with Nemes & Rawls (2006). The correlation coefficient of BD was negative for low hm and positive for high hm (i.e. 1000 and 1500 kPa). Dexter et al. (2008) found that BD decreased the SWR at low hm and increased it at high hm. Aina & Periaswamy (1985) also found a negative correlation between θ33 kPa and BD for western Nigerian soils. However, Hutson & Cass (1987) reported that BD had a non-significant effect on the SWR. Clay and silt increased the SWR for the major of hm values due to their effects on water retention surfaces. But the effect of silt was not significant at 0.01. This is in agreement with Hutson & Cass (1987). Aina & Periaswamy (1985) found positive relations of θ33 kPa with silt and clay contents and of θ1500 kPa with clay content. Kern (1995) found that effects of silt content and BD on θ33 kPa and θ1500 kPa and effect of sand content on θ1500 kPa were not significant but that sand content significantly affected θ33 kPa. Nemes & Rawls (2006) reported that sand and silt contents decreased and clay content increased θFC and θPWP. Dexter et al. (2008) found a positive relation between θPWP and clay content. Decreasing effect of sand content on the SWR was obvious and had a descending trend fromθ10 to θ1500 kPa. Its effect on θ1000 and θ1500 kPa was not significant. The influence of CaCO3 was minor and had no significant effect on SWR. It increased the SWR due to its impacts on aggregation and flocculation (by soluble Ca2+). Rajkai & Varallyay (1992) found that CaCO3 content was the second important independent variable entering PTFs for θ1500 kPa. Khodaverdiloo & Homaee (2004), however, observed that carbonates decrease the SWR at high hm because the carbonates with clay size behave like silt in water retention. Effect of OC on the SWR was not significant. Increasing effect of organic matter on the SWR is dominant at low hm. Gupta & Larson (1979) also reported that organic matter highly increases the SWR at low hm. Bell & van Keulen (1995) found a positive and significant relation between θPWP and organic matter in four Mexican soils.
Correlation coefficients among the water content in each potential and soil properties from total measured data.
Based on the collected training data, the following PTFs using MLR method were developed and are listed in Table 3.The independent variables were selected at a 0.01 significance level for inclusion in the regression model. The values of adjusted coefficients of determination (R2adj) given in Table 3 indicate the range of the dependent variable that is explained by the independent variables. These results show that at low pressure heads, the clay, sand and bulk density are appropriate predictors of water content. At high pressure heads (1000 and 1500 kPa), the clay content was a better predictor. These regression equations and their statistics are shown in Table 3. Table 3 showed that clay content is an important factor for predicting water content at all potential points. Also, the studied soils showed a high water content (saturation ≅ 0.60 cm3/cm3 and wilting point ≅ 0.20 cm3/cm3).This feature can be related to high silt and clay particle size content (e.g.,Rawls et al., 2003), which was always greater than 65% for all samples.
The point PTFs derived from multiple linear regression (MLR) method (training data).
Clay type plays an important role in the retention properties of soils. So, soils in the humid tropics can have a much lower capacity to retain water than soils in the temperate regions with the same clay content but with a different type of clays (Christopher et al., 2013). Soils in southern Guilan province are principally dominated by montmorillonite clay (Fig. 2) which is expanding type clay minerals and swell by absorbing substantial amounts of water in its interlayer space. Presence of 2:1 smectitic clays, especially, montmorillonite in these soils was also an important specific soil characteristic that contributed to increase the water retention and hygroscopic water content, i.e., the water content in the driest part of the curve. Fig. 2 illustrates the X-ray difractogram of the clay fraction and shows presence of smectite clay minerals in studied soils.
X-ray diffractogram of the clay fraction (“d” means interlayer spacing).
Discussion
Comparison among different methods for prediction of soil water content from testing data is presented in Suppl. Figs. S2, S3 and S4 [pdfs online] and summarized in Table 4.
Comparison among different methods for prediction of soil water content from testing data.
These figures and Table 4 imply that the general prediction was good for the three methods, but the ANN method had excellent performance and the prediction values were close to the measured values.
The RMSE values for the ANN were smaller than that for the derived point PTFs and Rosetta model in all matric potentials. MLR and Rosetta software hold the second and third places, respectively. The RMSE values for the derived point PTFs were smaller than that for the Rosetta model, except kPa that Rosetta model had better estimation than regression (Table 4). The R2 values of ANN for all potentials were greater than regressions and regressions were greater than Rosetta, except kPa that Rosetta model had greater than regression. So, the accuracy of the ANN is better than that of the derived point PTFs and Rosetta model. This result is in line with the work done by Minasny et al. (2004) and Amini et al. (2005). But, with regard to evaluation criteria, our study is more accurate than other mentioned researches, because, selection of inputs and model designing (such as type of learning algorithm, the number of layers and neurons, etc,) is carried out more accurately.
The reason of this superior efficiency of ANNs models compared with the basic regression equations is probably because the PTFs that have derived from various areas have different efficiencies. On the other hand, according to the hypothesis of Schaap et al. (2001), for designing of a neural network we do not need a special equation. They also believe that by creating a suitable equation between input and output data we are able to achieve the best results. Also, due to the occurring of nonlinear equations between dependent variables and predicting variables, the neural network have the better efficiency compared with the basic regression equations. Yi et al. (2013) investigated the accuracy of ANN and analyzed the regression method using correlation coefficient and the RMSE. They reported that the neural network is able to predict the easily measurable soil parameters with more accuracy and less error. Similar results have been reported by Mukhlisin & Abd Rahman (2014) as well. They found that using ANN leads to less RMSE values than the MLR. They also reported that the neural network has not better efficiency than linear regression models in occasion of high stability of data. However, the high accuracy of data leads to more efficiency of neural network and also, shows the proper selection of testing and training data. Analysis of the ANN parameters suggested that more input variables were necessary to improve the prediction of soil parameters (Mermoud & Xu, 2006). The RMSE of the different neurons in hidden layer is plotted in Fig. 3. This figure illustrated that the best model obtained with 10, 5, 8, 4, 10, 9 and 7 neurons for 10, 33, 100, 300, 500, 1000 and 1500 kPa respectively.
Relationship between root mean square error (RMSE) and number of hidden layer neurons in different matric potentials in artificial neural network (ANN) method. Red bars are number of neurons in hidden layer with the lowest RMSE.
The RMSE values of different ANN-PTFs and regression-PTFs were lower in the prediction of volumetric water content at PWP than the others. Likewise, based on MBE values, all these PTFs especially at PWP, showed slight underestimation of volumetric water content. But this underestimation of ANN model was very low and could be ignored. Also, in the evaluation study of Tomasella & Hodnett (2004), the general trend of lower RMSE values at permanent wilting point compared to field capacity conditions was observed for all the PTFs without any exception. Recently, Vereecken et al. (2010), also found that RMSE values at -1500 kPa are generally lower than at -33 kPa. The water content in the wet range of the retention curve is related more to the soil structural properties whereas in the dry range it depends more on the particle size distribution. Thus, the point PTFs were able to describe the water retention at different soil water pressure heads by incorporating different independent variables. The relationships between water retention parameters and basic soil properties are complex; consequently, the performance of point PTFs is better than that of Rosetta model. A similar comparison was made by Tomasella et al. (2003) both of whom reported similar differences between these two types of soil water prediction. The results of our investigation showed that in the reliability test, all the derived PTFs were associated with negative values of the MBE. Therefore, for the independent data set, the derived PTFs tend to underestimate the water retention curve (Table 4). As mentioned earlier, the range in clay content was not the major cause of its low predictive capability. Many authors, as Medina et al. (2002), indicated that clay type plays a vital role in the retention and transmission properties of a given soil. This is why soils in the humid tropics can have a much lower capacity to retain water than soils in the temperate regions with the same clay content but with a different type of clays (Christopher et al., 2013). Soils in southern Guilan province are principally dominated by montmorillonite clay (Fig. 3) which is expanding type clay minerals and swell by absorbing substantial amounts of water in its interlayer space. The prediction errors of these methods increased with increasing water content (high matric potentials) probably because at high matric potentials the water content is more related to the soil structure attributes. Because there was no predictor to adequately explain the soil structure effects on soil water content, and soil bulk density alone seems to be insufficient to represent these effects, the error prediction of the PTFs increased at high matric potentials. Comparison of these values with those obtained for the derived PTFs (Table 4) shows that all of the derived PTFs performed better than the Rosetta package. This may be related to the pedogenic differences between the data sets used to derive the Rosetta package and those for the reliability test. These low adjustments of Rosetta model, especially at–1500 kPa, could be related to the textural properties of these soils. Using a different particle size distribution, as does the Canadian Society of Soils Science (McKeague, 1978), a detailed silt particle-size distribution is obtained, where one can observe that exactly the fine silt content (5–2 μm) was between 57%and 84%. Thus, the hydrodynamic behavior of fine silt content, that could be similar to a clay particle size, would increase the water retention content near to PWP, and this would therefore justify the poor fits of Rosetta, because Rosetta predicts according to USDA textural distribution using three well-differentiated particle sizes (sand, silt and clay).The MBE values show that the Rosetta package overestimated water retention, especially at high pressure heads, except at 100 kPa pressure. This was also reported by Schaap et al. (2001).
With decreasing of matric potential (from 10 kPa to 1500 kPa), correlation coefficient increased in all three methods. In other words, the accuracy of prediction increased with increasing in soil suction. Because soils of study area had low content of organic matter and so, soil structure was very weak. On the other hand, the hydrodynamic behavior of fine silt content, that could be similar to a clay particle size, would increase the water retention content near to permanent wilting point. The content of fine silt of soils was between 57% and 84%.
In this study, multivariate linear regression and neural network model were employed to develop a pedotransfer function for predicting soil moisture using available soil properties. This neural network was consisted of three hidden layers, a sigmoid activation function in hidden and linear function in output layer. Results showed that artificial neural network gave the best model with Levenberg–Marquardt learning algorithm and tangent sigmoid (tansig) transfer function. Multi-Layer Perceptron architecture of different matric potentials according to number of inputs, number of neurons in the hidden layer and output parameter were as: 3-10-1, 3-5-1, 2-8-1, 2-4-1, 2-10-1, 1-9-1 and 1-7-1 for 10, 33, 100, 300, 500, 1000 and 1500 kPa, respectively. The ANN model was more suitable for capturing the non-linearity of the relationship between variables than multivariate regression and Rosetta, and can model non-linear functions and has been shown to perform better than linear regression. With regarding to the evaluation criteria, the results of this study revealed that the ANNs had superiority to the basic regression equations for prediction of mentioned soil parameter. This is a crucial result because, since ANN– PTFs formed from local data produce more accurate predictions than those built from data spread from a wider area, the concept of data conservation becomes a critical factor in ANN–PTFs construction (Baker & Ellison, 2008). However, due to difficulties of direct measurement of soil parameters, we recommend using of neuro-fuzzy models in the future studies for obtaining the logical equations of other soil parameters, especially soil hydraulic properties, in each area.
Acknowledgements
This paper is a part of M.Sc thesis of corresponding author. So, the authors thank for laboratory assistance of Soil Science Department technicians in University of Guilan. The authors are grateful to anonymous reviewers who considerably improved the quality of the manuscript.
ReferencesAbbasiYGhanbarian-AlavijehBLiaghatAMShorafaMEvaluation of pedotransfer functions for estimating soil water retention curve of saline and saline-alkali soils of Iran 2011212230237http://dx.doi.org/10.1016/S1002-0160(11)60122-7AinaPOPeriaswamySPEstimating available water-holding capacity of western Nigerian soils from soil texture and bulk density, using core and sieved samples198514015558http://dx.doi.org/10.1097/00010694-198507000-00007AminiMAbbaspourKCKhademiHFathianpourNAfyuniMSchulinRNeural network models to predict cation exchange capacity in arid regions of Iran200556551559http://dx.doi.org/10.1111/j.1365-2389.2005.0698.xBakerLEllisonDThe wisdom of crowds. Ensembles and modules in environmental modelling200814717http://dx.doi.org/10.1016/j.geoderma.2008.07.003BellMAvanKeulenHSoil pedotransfer functions for four Mexican soils199559865871http://dx.doi.org/10.2136/sssaj1995.03615995005900030034xBoumaJUsing soil survey data for quantitative land evaluation19899177213http://dx.doi.org/10.1007/978-1-4612-3532-3_4BrooksRHCoreyATHydrology Paper No. 31964Colorado Sta. UnivFort Collins, CO, USABurtRSoil survey laboratory methods manual2004No. 42, V. 4US Dept. Agr. Nat. Resour. Conserv. Serv. Nat. Soil Surv. CentChangYMChangLCChangFJComparison of static feed forward and dynamic-feedback neural networks for rainfall runoff modeling2004290297311http://dx.doi.org/10.1016/j.jhydrol.2003.12.033ChristopherTBeckettSCharlesEAugardeEPrediction of soil water retention properties using pore-size distribution and porosity2013504435450http://dx.doi.org/10.1139/cgj-2012-0320CornelisWMRonsynJMeirvenneMVEvaluation of pedotransfer functions for predicting the soil moisture retention curve200165638648http://dx.doi.org/10.2136/sssaj2001.653638xDemuthHBealeM2004The Math Works Inc.Natick, MA, USADexterARCzyEARichardGA user-friendly water retention function that takes account of the textural and structural pore spaces in soil2008143243253http://dx.doi.org/10.1016/j.geoderma.2007.11.010Ghanbarian-AlavijehBMillánHPoint pedotransfer functions for estimating soil water retention curve201024243251GiviJPrasherSOPatelRMEvaluation of pedotransfer functions in predicting the soil water contents at field capacity and wilting point2004708396http://dx.doi.org/10.1016/j.agwat.2004.06.009GuptaSCLarsonWEEstimating soil water characteristic from particle size distribution, organic matter percent, and bulk density19791516331635http://dx.doi.org/10.1029/WR015i006p01633HamblinASustainable agricultural systems, what are the appropriate measures for soil structure?199129709715http://dx.doi.org/10.1071/SR9910709HeikkiNK2008HuangMBFredlundDGFredlundMDComparison of measured and PTF predictions of SWCCs for Loess soils in China201028105117http://dx.doi.org/10.1007/s10706-009-9284-xHuangWFooSNeural network modelling of salinity variation in Apalachicola River2002361356362http://dx.doi.org/10.1016/S0043-1354(01)00195-6HutsonJLCassAA retentivity function for use in soil-water simulation models198738105113http://dx.doi.org/10.1111/j.1365-2389.1987.tb02128.xJainAKumarAAn evaluation of artificial neural network technique for the determination of infiltration model parameters200663272282http://dx.doi.org/10.1016/j.asoc.2004.12.007KernJSEvaluation of soil water retention models based on basic soil physical properties19955911341141http://dx.doi.org/10.2136/sssaj1995. 03615995005900040027xKhodaverdilooHHomaeeMWhrleNSheurerMPedotransfer functions of some calcareous soils2004Int ConfSeptember 4-12Freiburg, Germany1027111KumarKKumarYVenkateshVEstimation of liquid and plastic limit using artificial neural network modelsIndian Geotechnical ConferenceDecember 16–182010873876McKeagueJAManual on soil sampling and methods of analysis1978Ottawa, CanadaMedinaHTarawallyMdel ValleARuizMEEstimating soil water retention curve in Rhodic Ferralsols from basic soil data2002108277285http://dx.doi.org/10.1016/S0016-7061(02)00135-0MerdunHAlternative methods in the development of pedotransfer functions for soil hydraulic characteristics2010436271http://dx.doi.org/10.1134/S1064229310010084MerdunHCinarOMeralRApanMComparison of artificial neural network and regression pedotransfer functions for prediction of soil water retention and saturated hydraulic conductivity200690108116http://dx.doi.org/10.1016/j.still.2005.08.011MermoudAXuDComparative analysis of three methods to generate soil hydraulic functions20068789100http://dx.doi.org/10.1016/j.still.2005.02. 034MinasnyBHopmansJWHarterTNeural networks prediction of soil hydraulic functions for alluvial soils using multistep outflow data200468417429http://dx.doi.org/10.2136/sssaj2004.4170MoghimiSMahdianMParviziYMasihabadiMEstimating effects of terrain attributes on local soil organic carbon content in a semi-arid pastureland20145297106MualemYA new model for predicting the hydraulic conductivity of unsaturated porous media197612513522http://dx.doi.org/10.1029/WR012i003p00513MukhlisinMAbd RahmanASBPrediction of Atterberg limits via ANN and ANFIS: a comparison2014Proc Int Conf on Environmental Science and Geoscience (ESG ‘14)Venice, Italy6974Available in http://www.europment.org/library/2014/venice/ENVIR.pdfNajafiMGiviJEvaluation of prediction of bulk density by artificial neural network and PTFs200610th Iranian Soil Sci ConfKaraj26-28 Aug680681In PersianNemesARawlsWJEvaluation of different representations of the particle-size distribution to predict soil water retention20061324758http://dx.doi.org/10.1016/j.geoderma.2005.04.018NemesASchaapMGWöstenJHMFunctional evaluation of pedotransfer functions derived from different scales of data collection20036710931102http://dx.doi.org/10.2136/sssaj2003.1093RajkaiKVarallyayGVan GenuchtenMThLeijFJLundLJEstimating soil water retention from simpler properties by regression techniques1992Riverside, CA, USA11-13101992417426RawlsWJPachepskyYARitchieJCEffect of soil organic carbon on soil water retention20031166176http://dx.doi.org/10.1016/S0016-7061(03)00094-6Rezaei-ArshadRSayyadGhMosaddeghiMGharabaghiBPredicting Saturated Hydraulic Conductivity by Artificial Intelligence and Regression Models2013Hindawi Publishing Corporationhttp://dx.doi.org/10.1155/2013/308159RezaeianzadehMTabariHArabiYazdiAIsikSKalinLFlood flow forecasting using ANN, ANFIS and regression models2014252537http://dx.doi.org/10.1007/s00521-013-1443-6SantraPDasBSPedotransfer functions for soil hydraulic properties developed from a hilly watershed of Eastern India2008146439448http://dx.doi.org/10.1016/j.geoderma.2008.06.019SarmadianFMehrjardiRTModeling of some soil properties using artificial neural network and multivariate regression in Gorgan province, north of Iran200823035SchaapMGLeijFJImproved prediction of unsaturated hydraulic conductivity with the Mualem-Van Genuchten model200064843851http://dx.doi.org/10.2136/sssaj2000.643843xSchaapMGLeijFJVan GenuchtenMThRosetta, a computer program for estimating soil hydraulic parameters with hierarchical pedotransfer functions2001251163176http://dx.doi.org/10.1016/S0022-1694(01)00466-8StumppCEngelhardtSHofmannMEvaluation of pedotransfer functions for estimating soil hydraulic properties of prevalent soils in a catchment of the Bavarian Alps2009128609620http://dx.doi.org/10.1007/s10342-008-0241-7TomasellaJHodnettMGPachepskyYARawlsWJPedotransfer functions for tropical soils2004ElsevierAmsterdam415429http://dx.doi.org/10.1016/S0166-2481(04)30021-8TomasellaJPachepskyYCrestanaSRawlsWJComparison of two techniques to develop pedotransfer functions for water retention20036710851092http://dx.doi.org/10.2136/sssaj2003.1085Van GenuchtenMThA closed-form equation for predicting the hydraulic conductivity of unsaturated soils198044892898http://dx.doi.org/10.2136/sssaj1980.03615995004400050002xVenkatesanPAnithaSApplication of a radial basis function neural network for diagnosis of diabetes mellitus200691911951199VereeckenHWeynantsMJavauxMPachepskyYSchaapMGVan GenuchtenMThUsing pedotransfer functions to estimate the van Genuchten-Mualem soil hydraulic properties: a review20109126http://dx.doi.org/10.2136/vzj2010.0045WagnerBTarnawskiVRHenningsVEvaluation of pedo-transfer functions for unsaturated soil hydraulic conductivity using an independent data set2001102275297http://dx.doi.org/10.1016/S0016-7061(01)00037-4YiXLiGYinYComparison of three methods to develop pedotransfer functions for the saturated water content and field water capacity in permafrost region2013881016http://dx.doi.org/10.1016/j.coldregions.2012.12.005