Improved Prediction of Unsaturated Hydraulic Conductivity with the Mualem‐van Genuchten Model

In many vadose zone hydrological studies, it is imperative that the soil's unsaturated hydraulic conductivity is known. Frequently, the Mualem-van Genuchten model (MVG) is used for this purpose because it allows prediction of unsaturated hydraulic conductivity from water retention parameters. For this and similar equations, it is often assumed that a measured saturated hydraulic conductivity (K s ) can be used as a matching point (K o ) while a factor S L e is used to account for pore connectivity and tortuosity (where S e is the relative saturation and L = 0.5). We used a data set of 235 soil samples with retention and unsaturated hydraulic conductivity data to test and improve predictions with the MVG equation. The standard practice of using K o = K, and L = 0.5 resulted in a root mean square error for log(K) (RMSE K ) of 1.31. Optimization of the matching point (K o ) and L to the hydraulic conductivity data yielded a RMSE K of 0.41. The fitted K 0 were, on average, about one order of magnitude smaller than measured K s . Furthermore, L was predominantly negative, casting doubt that the MVG can be interpreted in a physical way, Spearman rank correlations showed that both K 0 and L were related to van Genuchten water retention parameters and neural network analyses confirmed that K 0 and L could indeed be predicted in this way. The corresponding RMSE K was 0.84, which was half an order of magnitude better than the traditional MVG model. Bulk density and textural parameters were poor predictors while addition of K s improved the RMSE K only marginally. Bootstrap analysis showed that the uncertainty in predicted unsaturated hydraulic conductivity was about one order of magnitude near saturation and larger at lower water contents.

[1]  G. G. Stokes.  "J." , 1890, The New Yale Book of Quotations.

[2]  Feike J. Leij,et al.  Analysis of Measured, Predicted, and Estimated Hydraulic Conductivity Using the RETC Computer Program , 1992 .

[3]  E. Bruce Jones,et al.  Watershed Management in the Eighties , 1985 .

[4]  I. Fatt,et al.  Relative Permeability Studies , 1951 .

[5]  S. Powers,et al.  Non–aqueous phase liquid dissolution in heterogeneous systems: Mechanisms and a local equilibrium modeling approach , 1998 .

[6]  John A. Nelder,et al.  A Simplex Method for Function Minimization , 1965, Comput. J..

[7]  J. Feyen,et al.  Solute transport in a heterogeneous soil for boundary and initial conditions: Evaluation of first‐order approximations , 1998 .

[8]  Peter Finke,et al.  Comparison of class and continuous pedotransfer functions to generate soil hydraulic characteristics , 1995 .

[9]  S. Hyakin,et al.  Neural Networks: A Comprehensive Foundation , 1994 .

[10]  H. Vereecken.  Estimating the unsaturated hydraulic conductivity from theoretical models using simple soil properties , 1995 .

[11]  J. C. van Dam,et al.  Comparison of six methods to determine unsaturated soil hydraulic conductivity. , 1994 .

[12]  N. T. Burdine.  Relative Permeability Calculations From Pore Size Distribution Data , 1953 .

[13]  Harry Vereecken,et al.  ESTIMATING THE SOIL MOISTURE RETENTION CHARACTERISTIC FROM TEXTURE, BULK DENSITY, AND CARBON CONTENT , 1989 .

[14]  T. Vogel,et al.  On the reliability of unsaturated hydraulic conductivity calculated from the moisture retention curve , 1988 .

[15]  W. Durner.  Hydraulic conductivity estimation for soils with heterogeneous pore structure , 1994 .

[16]  D. R. Nielsen,et al.  On describing and predicting the hydraulic properties of unsaturated soils , 1985 .

[17]  P. McCullagh.  Tensor Methods in Statistics , 1987 .

[18]  G. Dagan,et al.  Hydraulic Conductivity of Soils: Unified Approach to the Statistical Models1 , 1978 .

[19]  D. R. Nielsen,et al.  A consistent set of parametric models for the two‐phase flow of immiscible fluids in the subsurface , 1989 .

[20]  J. Feyen,et al.  ESTIMATING UNSATURATED HYDRAULIC CONDUCTIVITY FROM EASILY MEASURED SOIL PROPERTIES , 1990 .

[21]  D. K. Cassel,et al.  EVALUATION OF SPATIAL DISTRIBUTION OF HYDRAULIC CONDUCTIVITY USING EFFECTIVE POROSITY DATA , 1989 .

[22]  Y. Pachepsky,et al.  Artificial Neural Networks to Estimate Soil Water Retention from Easily Measurable Data , 1996 .

[23]  M. Schaap,et al.  Modeling water retention curves of sandy soils using neural networks , 1996 .

[24]  Van Genuchten,et al.  A closed-form equation for predicting the hydraulic conductivity of unsaturated soils , 1980 .

[25]  Wolfgang Durner,et al.  State-of-the-Art in Inverse Modeling of Inflow/Outflow Experiments , 1999 .

[26]  Marcel G. Schaap,et al.  Database-related accuracy and uncertainty of pedotransfer functions , 1998 .

[27]  M. Schaap,et al.  Neural network analysis for hierarchical prediction of soil hydraulic properties , 1998 .

[28]  Ken'ichirou Kosugi,et al.  General model for unsaturated hydraulic conductivity for soils with lognormal pore-size distribution. , 1999 .

[29]  Y. Mualem.  A New Model for Predicting the Hydraulic Conductivity , 1976 .

[30]  W. Rawls,et al.  Estimating generalized soil-water characteristics from texture , 1986 .

[31]  D Or,et al.  Microgravity effects on water flow and distribution in unsaturated porous media: Analyses of flight experiments , 1999, Water resources research.

[32]  Karsten H. Jensen,et al.  Laboratory investigations of effective flow behavior in unsaturated heterogeneous sands , 1999 .

[33]  W. M. Schuh,et al.  Effect of Soil Properties on Unsaturated Hydraulic Conductivity Pore-Interaction Factors , 1990 .

[34]  Ken'ichirou Kosugi,et al.  Three‐parameter lognormal distribution model for soil water retention , 1994 .

[35]  W. M. Schuh,et al.  Effect of Soil Properties on Hydraulic Conductivity‐Moisture Relationships , 1986 .