Using inverse methods for estimating soil hydraulic properties from field data as an alternative to direct methods

Water and solute transport in the vadose zone greatly depends on me physical and chemical properties of the soil, which generally exhibit high variability. Additionally, the experimental determination of those properties in the field or laboratory is tedious, time-consuming and involves considerable uncertainty for most practical applications. Recently, inverse modeling has been introduced to estimate effective properties in situ by deducing them from, e.g. a measured time series of soil water content. Inverse methods combine forward soil water flow models with appropriate optimization algorithms to find the best parameter set that minimizes an objective function. Global optimization methods are suitable for locating a global optimum for a given set of conditions (number of parameters, boundary conditions, etc.). In this raper we estimate the soil hydraulic properties of a sprinkler fertigated banana plot in the North of Tenerife (Canary Islands) in a direct and inverse way. For the inverse method, use was made of the measured time series of soil water content at three different depths. The forward model is the numerical solution of the Richards equation as implemented in the agro-environmental model WAVE. Two inverse methods are compared: the traditional "trial and error" method and an inverse method using a global search algorithm referred to as the global multilevel coordinate search combined sequentially with the Nelder-Mead simplex algorithm (GMCS-NMS). The global search is shown to be a relatively efficient procedure for estimating the soil hydraulic properties from measured soil water contents inthe banana plantation. However, some ill-posedness problems are identified which should be solved by upgrading the quality of the experimental programs in irrigated banana plantations. (C) 2002 Elsevier Science B.V. All rights reserved.

[1]  Marnik Vanclooster,et al.  A global multilevel coordinate search procedure for estimating the unsaturated soil hydraulic properties , 2002 .

[2]  S. P. Neuman,et al.  Estimation of aquifer parameters under transient and steady-state conditions: 2 , 1986 .

[3]  H. van Keulen,et al.  A summary model for crop growth , 1982 .

[4]  R. Rabbinge,et al.  Simulation and systems management in crop protection , 1989, Plant Growth Regulation.

[5]  M. J. Hall,et al.  Rainfall-Runoff Modelling , 2004 .

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

[7]  Reinder A. Feddes,et al.  Simulation model of the water balance of a cropped soil: SWATRE , 1983 .

[8]  J. Dam,et al.  From one-step to multi-step. Determination of soil hydraulic functions by outflow experiments. , 1990 .

[9]  M. Vanclooster,et al.  Simulating water and nitrogen behaviour in soils cropped with winter wheat , 2005, Fertilizer Research.

[10]  J. Šimůnek,et al.  In situ estimation of soil hydraulic functions using a multistep soil‐water extraction technique , 1998 .

[11]  Michael H. Young,et al.  Estimation of depth averaged unsaturated soil hydraulic properties from infiltration experiments , 2001 .

[12]  J. H. Dane,et al.  In-Situ Determination of Soil Hydraulic Properties during Drainage1 , 1983 .

[13]  Jan W. Hopmans,et al.  Optimization of Hydraulic Functions from Transient Outflow and Soil Water Pressure Data , 1993 .

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

[15]  D. L. Thomas,et al.  Agricultural Non-Point Source Water Quality Models: Their Use and Application , 1998 .

[16]  A. P. Annan,et al.  Electromagnetic determination of soil water content: Measurements in coaxial transmission lines , 1980 .

[17]  Wim G.M. Bastiaanssen,et al.  Inversely estimating soil hydraulic functions using evapotranspiration fluxes , 2002 .

[18]  Marnik Vanclooster,et al.  Evaluation of the WAVE model , 2001 .

[19]  Arnold Neumaier,et al.  Global Optimization by Multilevel Coordinate Search , 1999, J. Glob. Optim..

[20]  Uri Shani,et al.  Analyses of infiltration events in relation to determining soil hydraulic properties by inverse problem methodology , 1991 .

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

[22]  M. T. van Genuchten,et al.  Estimating Unsaturated Soil Hydraulic Properties from Tension Disc Infiltrometer Data by Numerical Inversion , 1996 .

[23]  M. Vanclooster,et al.  A deterministic evaluation analysis applied to an integrated soil-crop model , 1995 .

[24]  W. Yeh Review of Parameter Identification Procedures in Groundwater Hydrology: The Inverse Problem , 1986 .

[25]  W. Bouten,et al.  Potentials and limitations of modeling vertical distributions of root water uptake of an Austrian pine forest on a sandy soil , 2000 .

[26]  R. Muñoz‐Carpena,et al.  Nitrogen evolution and fate in a Canary Islands (Spain) sprinkler fertigated banana plot , 2002 .