Characterization of effective saturated hydraulic conductivity in an agricultural field using Karhunen‐Loève expansion with the Markov chain Monte Carlo technique

[1] Process-based soil hydrologic models require input of saturated hydraulic conductivity (Ksat). However, model users often have limited access to measured data and thus use published or estimated values for many site-specific hydrologic and environmental applications. We proposed an algorithm that uses the Karhunen-Loeve expansion (KLE) in conjunction with the Markov chain Monte Carlo (MCMC) technique, which employs measured soil moisture values to characterize the saturated hydraulic conductivity of an agricultural field at a 30 m resolution. The study domain is situated in the Walnut Creek watershed, Iowa, with soybean crop (in 2005) and well-defined top (atmospheric) and bottom (groundwater) boundary conditions. The KLE algorithm parameterizes and generates Ksat fields with random correlation lengths that are used in the SWMS_3D model for predicting the soil moisture dynamics for two different scenarios: (1) the van Genuchten soil hydraulic parameters (except Ksat) are constant and are based on the soil type of the grid block within the domain, and (2) Ksat is correlated with the van Genuchten parameter α as Ksat ∝ α2. The predicted soil moisture fields for both the scenarios are evaluated with the measured soil moisture in the MCMC algorithm for acceptance (or rejection) of the Ksat fields. The accepted Ksat fields are evaluated against the laboratory-measured Ksat at specific locations as well as with a large Ksat data set measured in situ in a nearby field with similar hydrologic conditions, and the comparisons show reasonably good agreement. The KLE-MCMC algorithm was further tested in the same study domain for another year (2002) having different vegetation (corn) and local forcings. The algorithm shows potential to characterize the effective saturated hydraulic conductivity fields at 30 m resolution using inexpensive and more regularly measured soil moisture measurements. Further studies are required to incorporate variability in different hydroclimatic regions and diverse topography to extend the application of this algorithm.

[1]  Y. C. Oq.  Variation of hydraulic conductivity in a tilled soil , 2002 .

[2]  B. Mohanty.  Scaling hydraulic properties of a macroporous soil , 1999 .

[3]  G. J. Wall,et al.  Stochastic Analysis of Groundwater Levels in a Temperate Climate , 1992 .

[4]  Michel Loève,et al.  Probability Theory I , 1977 .

[5]  M. Vanclooster,et al.  Three‐Dimensional Modeling of the Scale‐ and Flow Rate‐Dependency of Dispersion in a Heterogeneous Unsaturated Sandy Monolith , 2006 .

[6]  B. Mohanty,et al.  Scaling of near‐saturated hydraulic conductivity measured using disc infiltrometers , 1998 .

[7]  D. Rolph,et al.  Infiltration Rate as Affected by an Alfalfa and No-till Cotton Cropping System , 1990 .

[8]  Dongxiao Zhang,et al.  Stochastic analysis of saturated–unsaturated flow in heterogeneous media by combining Karhunen-Loeve expansion and perturbation method , 2004 .

[9]  G. J. Wall,et al.  COMPARISON OF SATURATED HYDRAULIC CONDUCTIVITY MEASURED BY VARIOUS FIELD METHODS , 1993 .

[10]  Keith Loague,et al.  R‐5 revisited: 1. Spatial variability of infiltration on a small rangeland catchment , 1990 .

[11]  D. Hillel,et al.  COMPARISON OF THREE METHODS FOR ASSESSING SOIL HYDRAULIC PROPERTIES , 1993 .

[12]  John Geweke,et al.  Evaluating the accuracy of sampling-based approaches to the calculation of posterior moments , 1991 .

[13]  A. W. Warrick,et al.  13 – Spatial Variability of Soil Physical Properties in the Field , 1980 .

[14]  B. Mohanty,et al.  Saturated hydraulic conductivity and soil water retention properties across a soil‐slope transition , 2000 .

[15]  M. Celia,et al.  A General Mass-Conservative Numerical Solution for the Unsaturated Flow Equation , 1990 .

[16]  B. Mohanty,et al.  A Robust-Resistant Approach to Interpret Spatial Behavior of Saturated Hydraulic Conductivity of a Glacial Till Soil Under No-Tillage System , 1991 .

[17]  Alex. B. McBratney,et al.  Some considerations on methods for spatially aggregating and disaggregating soil information , 1998, Nutrient Cycling in Agroecosystems.

[18]  Stephan T. Grilli,et al.  Probabilistic analysis of flow in random porous media by stochastic boundary elements , 1997 .

[19]  N. Metropolis,et al.  The Monte Carlo method. , 1949 .

[20]  B. Mohanty,et al.  Spatial variability of residual nitrate-nitrogen under two tillage systems in central Iowa: A composite three-dimensional resistant and exploratory approach , 1994 .

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

[22]  Three-dimensional structure characterisation and transient flow modelling of a variably saturated heterogeneous monolith , 2006 .

[23]  D. R. Nielsen,et al.  Spatial variability of field-measured soil-water properties , 1973 .

[24]  D. Jaynes,et al.  Spatial variability of hydraulic conductivity in a cultivated field at different times , 1996 .

[25]  Marcel G. Schaap,et al.  Description of the unsaturated soil hydraulic database UNSODA version 2.0 , 2001 .

[26]  M. Vauclin,et al.  Modeling of unsaturated water flow in double-porosity soils by the homogenization approach , 2004 .

[27]  C. Dirksen,et al.  Hydraulic Conductivity and Diffusivity: Laboratory Methods , 2018, SSSA Book Series.

[28]  M. Sharma,et al.  Areal distribution of infiltration parameters and some soil physical properties in lateritic catchments , 1987 .

[29]  B. Mohanty,et al.  Soil Hydraulic Conductivities and their Spatial and Temporal Variations in a Vertisol , 2006 .

[30]  Warren J. Busscher,et al.  Simulation of Field Water Use and Crop Yield , 1980 .

[31]  Bernd Diekkrüger,et al.  Effective Soil Water Characteristics and Ensemble Soil Water Profiles in Heterogeneous Soils , 1996 .

[32]  R. Carsel,et al.  Developing joint probability distributions of soil water retention characteristics , 1988 .

[33]  Dongxiao Zhang,et al.  Conditional simulations of flow in randomly heterogeneous porous media using a KL-based moment-equation approach , 2004 .

[34]  Ramesh S. Kanwar,et al.  Comparison of Saturated Hydraulic Conductivity Measurement Methods for a Glacial-Till Soil , 1994 .

[35]  B. Mohanty,et al.  Spatial analysis of hydraulic conductivity measured using disc infiltrometers , 1994 .

[36]  D. McLaughlin,et al.  A Reassessment of the Groundwater Inverse Problem , 1996 .

[37]  R. Horton,et al.  EFFECT OF CULTIVATION ON HYDRAULIC PROPERTIES OF AN IOWA SOIL USING TENSION INFILTROMETERS1 , 1998 .

[38]  Yoram Rubin,et al.  HYDRO_GEN: A spatially distributed random field generator for correlated properties , 1996 .

[39]  H. P. Denton,et al.  Influence of Cover Crop and Wheel Traffic on Soil Physical Properties in Continuous No‐Till Corn , 1989 .

[40]  J. Bouma.  Field measurement of soil hydraulic properties characterizing water movement through swelling clay soils , 1980 .