Development of Pedotransfer Functions for Saturated Hydraulic Conductivity

摘要

The purpose of the study was to develop pedotransfer functions for determining saturated hydraulic conductivity (Ks). Pedotransfer functions (PTFs) for predicting soil physical properties used in determining saturated hydraulic conductivity, based on moisture retention characteristics, were developed. The van Genuchten moisture retention equation was fitted to measured moisture retention properties obtained from International Soil Reference and Information Centre (ISRIC) soils data base in order to determine parameters in the equation i.e. saturated soil moisture content (θs), residual soil moisture (θr), air entry parameter (α) and the pore size distribution parameter (n). 457 samples drawn from the data base were used to be the maximum possible sample size that contained the measured soils characteristics data required. Using statistical regression, mathematical relationships were developed between moisture retention parameters (response variables) and appropriately selected transformed basic soil properties (predictor variables). The developed PTFs were evaluated for accuracy and reliability. It was found that pedotransfer functions developed for θs produced the best performance in reliability compared to the remaining parameters yielding a correlation coefficient value of coefficient of determination (R2 = 0.76), RMSE = 2.09, NSE = 0.75 and RSR = 0.5 indicating good performance. Relatively poorest performance was obtained from the pedotransfer function developed for α which yielded a correlation coefficient, R2= 0.06, RMSE = 0.85 and a NSE of 0.02 reflecting the best possible equation derived for the parameter for use in predicting hydraulic conductivity. Out of the pedotransfer functions developed for each of the moisture retention parameters, the best performing PTF was identified for each parameter. The accuracy of the pedotransfer functions assessed based on R2 were for θs (R2 = 0.80), θr (R2 = 0.42), α (R2 = 0.04) and for n (R2 = 0.30), when the variables were expressed directly in terms of the selected transformations of the basic soil properties.