Inverse estimation of parameters for multidomain flow models in soil columns with different macropore densities

Soil and crop management practices have been found to modify soil structure and alter macropore densities. An ability to accurately determine soil hydraulic parameters and their variation with changes in macropore density is crucial for assessing potential contamination from agricultural chemicals. This study investigates the consequences of using consistent matrix and macropore parameters in simulating preferential flow and bromide transport in soil columns with different macropore densities (no macropore, single macropore, and multiple macropores). As used herein, the term"macropore density" is intended to refer to the number of macropores per unit area. A comparison between continuum-scale models including single-porosity model (SPM), mobile-immobile model (MIM), and dual-permeability model (DPM) that employed these parameters is also conducted. Domain-specific parameters are obtained from inverse modeling of homogeneous (no macropore) and central macropore columns in a deterministic framework and are validated using forward modeling of both low-density (3 macropores) and high-density (19 macropores) multiple-macropore columns. Results indicate that these inversely modeled parameters are successful in describing preferential flow but not tracer transport in both multiple-macropore columns. We believe that lateral exchange between matrix and macropore domains needs better accounting to efficiently simulate preferential transport in the case of dense, closely spaced macropores. Increasing model complexity from SPM to MIM to DPM also improved predictions of preferential flow in the multiple-macropore columns but not in the single-macropore column. This suggests that the use of a more complex model with resolved domain-specific parameters is recommended with an increase in macropore density to generate forecasts with higher accuracy.

[1]  M. Ye,et al.  Estimation of effective unsaturated hydraulic conductivity tensor using spatial moments of observed moisture plume , 2004 .

[2]  J. Vrugt,et al.  A formal likelihood function for parameter and predictive inference of hydrologic models with correlated, heteroscedastic, and non‐Gaussian errors , 2010 .

[3]  Christer Jansson,et al.  Preferential water flow in a glacial till soil , 2005 .

[4]  J. Böhlke,et al.  Groundwater recharge and agricultural contamination , 2002 .

[5]  M. Brusseau,et al.  Modeling solute transport in structured soils: a review , 1990 .

[6]  Viliam Novák,et al.  Infiltration of water into soil with cracks. , 2000 .

[7]  R. Jamieson,et al.  Movement and persistence of fecal bacteria in agricultural soils and subsurface drainage water: A review , 2002 .

[8]  N. Jarvis.  A review of non‐equilibrium water flow and solute transport in soil macropores: principles, controlling factors and consequences for water quality , 2007 .

[9]  Evaluation of Sorption Models in the Simulation of Naphthalene Transport Through Saturated Soils , 1991 .

[10]  R. Freeze.  A stochastic‐conceptual analysis of one‐dimensional groundwater flow in nonuniform homogeneous media , 1975 .

[11]  J. W. Biggar,et al.  Miscible Displacement: II. Behavior of Tracers1 , 1962 .

[12]  Yu-Shu Wu,et al.  Sensitivity analysis of hydrological parameters in modeling flow and transport in the unsaturated zone of Yucca Mountain, Nevada, USA , 2006 .

[13]  J. Šimůnek,et al.  Bromide transport at a tile-drained field site: experiment, and one- and two-dimensional equilibrium and non-equilibrium numerical modeling , 2006 .

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

[15]  N. J. Jarvis,et al.  The MACRO Model (version 3.1). Technical description and sample simulations. , 1994 .

[16]  M. V. Genuchten,et al.  Review and comparison of models for describing non-equilibrium and preferential flow and transport in the vadose zone , 2003 .

[17]  B. Mohanty,et al.  Water flow processes in a soil column with a cylindrical macropore: Experiment and hierarchical modeling , 2005 .

[18]  K. Beven,et al.  Macropores and water flow in soils , 1982 .

[19]  K. E. Johnsen,et al.  Macropore Transport of a Surface-Applied Bromide Tracer: Model Evaluation and Refinement , 1995 .

[20]  M. V. Genuchten,et al.  Estimating unsaturated soil hydraulic properties from laboratory tension disc infiltrometer experiments , 1999 .

[21]  Peter J. Shouse,et al.  Lateral Water Diffusion in an Artificial Macroporous System: Modeling and Experimental Evidence , 2003 .

[22]  M. Muir.  Physical Chemistry , 1888, Nature.

[23]  M. Thévenot,et al.  Evaluating equilibrium and non-equilibrium transport of bromide and isoproturon in disturbed and undisturbed soil columns. , 2007, Journal of contaminant hydrology.

[24]  John Maximilian Köhne,et al.  A review of model applications for structured soils: a) Water flow and tracer transport. , 2009, Journal of contaminant hydrology.

[25]  J. Šimůnek,et al.  Multi-process herbicide transport in structured soil columns: experiments and model analysis. , 2006, Journal of contaminant hydrology.

[26]  Markus Weiler,et al.  An infiltration model based on flow variability in macropores: development, sensitivity analysis and applications , 2005 .

[27]  Wolfgang Durner,et al.  Determination of parameters for bimodal hydraulic functions by inverse modeling , 1998 .

[28]  Marco Bittelli,et al.  Conceptual Models of Flow and Transport in the Fractured Vadose Zone. , 2002 .

[29]  Jirka Simunek,et al.  Numerical evaluation of a second‐order water transfer term for variably saturated dual‐permeability models , 2004 .

[30]  Peter F. Germann,et al.  Kinematic Wave Approach to Infiltration and Drainage into and from Soil Macropores , 1985 .

[31]  F. Stagnitti,et al.  A new modelling approach to simulate preferential flow and transport in water repellent porous media: Parameter sensitivity, and effects on crop growth and solute leaching , 2005 .

[32]  Y. Mualem,et al.  A New Model for Predicting the Hydraulic Conductivity , 1976 .

[33]  H. Gerke,et al.  Estimating the hydraulic functions of dual‐permeability models from bulk soil data , 2002 .

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

[35]  J. Feyen,et al.  Modelling Water Flow and Solute Transport in Heterogeneous Soils: A Review of Recent Approaches☆ , 1998 .

[36]  Horst H. Gerke,et al.  Preferential flow descriptions for structured soils , 2006 .

[37]  M. V. Genuchten,et al.  A dual-porosity model for simulating the preferential movement of water and solutes in structured porous media , 1993 .

[38]  N. J. Jarvis,et al.  Quantifying the impact of preferential flow on solute transport to tile drains in a sandy field soil , 1999 .

[39]  M. V. Genuchten,et al.  Modeling flow and transport in a two-dimensional dual-permeability system with spatially variable hydraulic properties , 2000 .

[40]  Ole Wendroth,et al.  Non‐equilibrium water flow characterized by means of upward infiltration experiments , 2001 .

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

[42]  K. Jensen,et al.  Flow and transport processes in a macroporous subsurface-drained glacial till soil II. Model analysis , 1998 .

[43]  J. J. McKeown.  Non-Linear Parameter Estimation , 1972 .

[44]  Hannes Flühler,et al.  Transport of Chemicals Through Soil: Mechanisms, Models, and Field Applications , 1992 .

[45]  Willem Bouten,et al.  Toward Improved Identifiability of Soil Hydraulic Parameters: On the Selection of a Suitable Parametric Model , 2003 .

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

[47]  H. Gerke,et al.  Dual-permeability modeling of preferential bromide leaching from a tile-drained glacial till agricultural field , 2004 .

[48]  D. R. Nielsen,et al.  MISCIBLE DISPLACEMENT IN SOILS: I. EXPERIMENTAL INFORMATION , 1961 .

[49]  Jan Feyen,et al.  Simultaneous inverse estimation of soil hydraulic and solute transport parameters from transient field experiments: Homogeneous soil , 2003 .

[50]  C. S. Simmons,et al.  Recent Studies of Flow and Transport in the Vadose Zone , 1991 .

[51]  M. V. Genuchten,et al.  Macroscopic representation of structural geometry for simulating water and solute movement in dual-porosity media , 1996 .

[52]  Junfeng Zhu,et al.  Traditional analysis of aquifer tests: Comparing apples to oranges? , 2005 .

[53]  R. Müller,et al.  Soil and macro-pores under uniaxial compression. I. Mechanical stability of repacked soil and deformation of different types of macro-pores , 2008 .

[54]  References , 1971 .

[55]  L. C. Brown,et al.  Pesticide Transport to Subsurface Tile Drains in Humid Regions of North America , 2001 .

[56]  D. Legates,et al.  Evaluating the use of “goodness‐of‐fit” Measures in hydrologic and hydroclimatic model validation , 1999 .

[57]  A. Franzluebbers,et al.  Tillage-induced seasonal changes in soil physical properties affecting soil CO2 evolution under intensive cropping , 1995 .

[58]  M. V. Genuchten,et al.  Mass transfer studies in sorbing porous media. I. Analytical solutions , 1976 .

[59]  J. Šimůnek,et al.  Impact of Soil Micromorphological Features on Water Flow and Herbicide Transport in Soils , 2008 .

[60]  R. Wagenet,et al.  Simulation of water and chemicals in macropore soils Part 1. Representation of the equivalent macropore influence and its effect on soilwater flow , 1992 .

[61]  F. Larney,et al.  PREFERENTIAL LEACHING IN LARGE UNDISTURBED SOIL' BLOCKS FROM CONVENTIONAL TILLAGE AND NO-TILL FIELDS IN SOUTHERN ALBERTA , 1999 .