Fully-coupled physically-based approach for modeling conventional and managed subsurface drainage

Subsurface drainage is an essential water management practice that is partly responsible for the legendary agricultural productivity of the “corn-belt” region in the United States. Recent researches indicate that the conventional drainage practice has contributed a large percentage of the surface water pollution by agricultural chemicals while the managed drainage practice can address the environmental concerns as well as the needs for crop yields. The fully-coupled physically-based hydrology models are important tools for evaluating these environmental issues and agroeconomic benefits by considering the subsurface drainage practices into the integrated surface and subsurface flow environment at the agricultural fields. Our research on the applications of the integrated hydrology model HydroGeoSphere to the real tile-drained farmlands in Indiana and Illinois yields a new tile drain module that can simulate both the conventional and managed drainage practices. The new tile drain module features the embedded node method for integrating tile drains into the subsurface, the resistance adjustment method for computing lateral flow into tile drains, the Preissmann slot method for combining free surface and pressurized flow routings inside tile drains, and a ghost node flow boundary condition for the free fall and submergence conditions at the drainage outlet. The new tile drain module is verified with both hypothetical examples and laboratory experiments. Applications of HydroGeoSphere with the new tile drain module to 2D and 3D agricultural drainage fields are demonstrated and discussed.

[1]  Derek Winstanley,et al.  Hypoxia in the Gulf of Mexico , 1999 .

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

[3]  Rao S. Govindaraju,et al.  Dynamics of Moving Boundary Overland Flows Over Infiltrating Surfaces at Hillslopes , 1991 .

[4]  B. Mohanty,et al.  New piecewise‐continuous hydraulic functions for modeling preferential flow in an intermittent‐flood‐irrigated field , 1997 .

[5]  John L. Nieber,et al.  Drains as a Boundary Condition in Finite Elements , 1986 .

[6]  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 .

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

[8]  G. Ampt,et al.  Studies on Soil Physics: Part II — The Permeability of an Ideal Soil to Air and Water , 1912, The Journal of Agricultural Science.

[9]  R. Wayne Skaggs Effect of Drain Tube Openings on Water-Table Drawdown , 1978 .

[10]  A comparison of observed and simulated hydrograph separations for a field-scale rainfall-runoff experiment. , 2000 .

[11]  A. M. Garcia,et al.  Drainage equations for random and irregular tile drainage systems , 2001 .

[12]  David A. Woolhiser,et al.  Overland Flow on an Infiltrating Surface , 1971 .

[13]  Hongbin Zhan,et al.  Theoretical and experimental studies of coupled seepage-pipe flow to a horizontal well , 2003 .

[14]  R. Maxwell,et al.  The groundwater land-surface atmosphere connection: Soil moisture effects on the atmospheric boundary layer in fully-coupled simulations , 2007 .

[15]  T. Steenhuis,et al.  A simple mixing layer model predicting solute flow to drainage lines under preferential flow , 1996 .

[16]  Charles S. Sawyer,et al.  Productivity Comparison of Horizontal and Vertical Ground Water Remediation Well Scenarios , 1998 .

[17]  R. D. Black,et al.  An Experimental Investigation of Runoff Production in Permeable Soils , 1970 .

[18]  D. Kirkham Flow of ponded water into drain tubes in soil overlying an impervious layer , 1949 .

[19]  Peter A. Forsyth,et al.  Incomplete Factorization Methods for Fully Implicit Simulation of Enhanced Oil Recovery , 1984 .

[20]  E. Youngs The contribution of physics to land drainage , 1983 .

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

[22]  Arlen W. Harbaugh,et al.  A modular three-dimensional finite-difference ground-water flow model , 1984 .

[23]  G. Pinder,et al.  Computational Methods in Subsurface Flow , 1983 .

[24]  J. A. Cunge,et al.  Intégration numérique des équations d'écoulement de barré de Saint-Venant par un schéma implicite de différences finies , 1964 .

[25]  P. A. Forsyth Comparison of the single-phase and two-phase numerical model formulation for saturated-unsaturated groundwater flow , 1988 .

[26]  Peter A. Forsyth,et al.  Flow and transport in fractured tuff at Yucca Mountain: numerical experiments on fast preferential flow mechanisms , 2000 .

[27]  Horst H. Gerke,et al.  Evaluation of a first-order water transfer term for variably saturated dual-porosity flow models , 1993 .

[28]  Hongbin Zhan,et al.  Groundwater flow to a horizontal or slanted well in an unconfined aquifer , 2002 .

[29]  E. Sudicky,et al.  Three-dimensional analysis of variably-saturated flow and solute transport in discretely-fractured porous media , 1996 .

[30]  G. Chavent,et al.  Application of the mixed hybrid finite element approximation in a groundwater flow model: Luxury or necessity? , 1994 .

[31]  Integrated flow and transport processes in subsurface -drained agricultural fields , 2003 .

[32]  C. W. Thornthwaite An Approach Toward a Rational Classification of Climate , 1948 .

[33]  David A. Kovacic,et al.  Nitrogen Balance in and Export from an Agricultural Watershed , 1997 .

[34]  Larry S.K. Fung,et al.  Reservoir Simulation With a Control-Volume Finite-Element Method , 1992 .

[35]  K. Loague,et al.  Hydrologic‐Response simulations for the R‐5 catchment with a comprehensive physics‐based model , 2001 .

[36]  J. Stillman,et al.  A semi-analytical model for transient flow to a subsurface tile drain , 2006 .

[37]  Arlen W. Harbaugh,et al.  MODFLOW-2000, The U.S. Geological Survey Modular Ground-Water Model - User Guide to Modularization Concepts and the Ground-Water Flow Process , 2000 .

[38]  Peter A. Forsyth,et al.  A two‐phase, two‐component model for natural convection in a porous medium , 1991 .

[39]  Ben Chie Yen,et al.  Modeling of conjunctive two-dimensional surface-three-dimensional subsurface flows , 2002 .

[40]  R. Horton The Rôle of infiltration in the hydrologic cycle , 1933 .

[41]  L. S. Pereira,et al.  Crop evapotranspiration : guidelines for computing crop water requirements , 1998 .

[43]  J. Gallichand Representing subsurface drains in a finite difference model , 1993 .

[44]  E. O'Loughlin,et al.  Saturation regions in catchments and their relations to soil and topographic properties , 1981 .

[45]  René Therrien,et al.  Improved three-dimensional finite-element techniques for field simulation of variably saturated flow and transport , 1993 .

[46]  E. Sudicky,et al.  On the Incorporation of Drains into Three‐Dimensional Variably Saturated Groundwater Flow Models , 1996 .

[47]  P. Kalita,et al.  Simulation of base-flow and tile-flow for storm events in a subsurface drained watershed , 2009 .

[48]  Markus Flury,et al.  Experimental evidence of transport of pesticides through field soils - a review , 1996 .

[49]  R. Allan Freeze,et al.  Role of subsurface flow in generating surface runoff: 1. Base flow contributions to channel flow , 1972 .

[50]  Peter A. Forsyth,et al.  Practical considerations for adaptive implicit methods in reservoir simulation , 1986 .

[51]  F. Sartoretto,et al.  Linear Galerkin vs mixed finite element 2D flow fields , 2009 .

[52]  D. Rose,et al.  Some errors estimates for the box method , 1987 .

[53]  L. Durlofsky Accuracy of mixed and control volume finite element approximations to Darcy velocity and related quantities , 1994 .

[54]  J. Vanderkwaak Numerical simulation of flow and chemical transport in integrated surface-subsurface hydrologic systems , 1999 .

[55]  T. Narasimhan,et al.  AN INTEGRATED FINITE DIFFERENCE METHOD FOR ANALYZING FLUID FLOW IN POROUS MEDIA , 1976 .

[56]  G. O. Schwab,et al.  Effect of Openings on Inflow Into Corrugated Drains , 1977 .

[57]  W. Dierickx Electrolytic analogue study of the effect of openings and surrounds of various permeabilities on the performance of field drainage pipes , 1980 .

[58]  J. Hewlett Factors affecting the response of small watersheds to precipitation in humid areas , 1967 .

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

[60]  William John Northcott Modeling Water Quantity and Quality on Tile Drained Watersheds With a GIS Coupled DRAINMOD , 1999 .

[61]  R. McCuen,et al.  Comparison of Overland Flow Hydrograph Models , 1987 .

[62]  Jirka Simunek,et al.  Single-porosity and dual-porosity modeling of water flow and solute transport in subsurface-drained fields using effective field-scale parameters , 2005 .

[63]  Yousef Saad,et al.  Iterative methods for sparse linear systems , 2003 .

[64]  Eungyu Park,et al.  Hydraulics of a finite-diameter horizontal well with wellbore storage and skin effect , 2002 .

[65]  W. Kinzelbach,et al.  Continuous Groundwater Velocity Fields and Path Lines in Linear, Bilinear, and Trilinear Finite Elements , 1992 .

[66]  P. Huyakorn,et al.  Techniques for Making Finite Elements Competitve in Modeling Flow in Variably Saturated Porous Media , 1984 .

[67]  Richard P. Hooper,et al.  Flux and Sources of Nutrients in the Mississippi-Atchafalaya River Basin , 1999 .

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

[69]  M. B. David,et al.  Anthropogenic inputs of nitrogen and phosphorus and riverine export for Illinois, USA. , 2000 .

[70]  J. Nash,et al.  River flow forecasting through conceptual models part I — A discussion of principles☆ , 1970 .

[71]  Jim E. Jones,et al.  Newton–Krylov-multigrid solvers for large-scale, highly heterogeneous, variably saturated flow problems , 2001 .

[72]  L. B. Leopold,et al.  Water In Environmental Planning , 1978 .

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

[74]  Ben Chic Yen,et al.  Open-Channel Flow Equations Revisited , 1973 .

[75]  G. M. Pohll,et al.  Modeling Regional Flow and Flow to Drains , 1994 .

[76]  H. Gerke,et al.  Spatial and Temporal Dynamics of Preferential Bromide Movement towards a Tile Drain , 2005, Vadose Zone Journal.

[77]  R. W. Skaggs,et al.  Drain tube opening effects on drain inflow , 1983 .

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

[79]  P. Huyakorn,et al.  A fully coupled physically-based spatially-distributed model for evaluating surface/subsurface flow , 2004 .

[80]  G. Gertner,et al.  Eutrophication: Nitrate flux in the Mississippi River , 2001, Nature.

[81]  R. Maxwell,et al.  Integrated surface-groundwater flow modeling: A free-surface overland flow boundary condition in a parallel groundwater flow model , 2006 .

[82]  M. W. Kawecki Transient Flow to a Horizontal Water Well , 2000 .

[83]  C. Zheng,et al.  MODFLOW 2001 and Other Modeling Odysseys , 2003 .

[84]  Fred J. Molz,et al.  A physically based, two-dimensional, finite-difference algorithm for modeling variably saturated flow , 1994 .

[85]  Reed M. Maxwell,et al.  Quantifying the effects of three-dimensional subsurface heterogeneity on Hortonian runoff processes using a coupled numerical, stochastic approach , 2008 .

[86]  Peter A. Forsyth,et al.  A Control Volume Finite Element Approach to NAPL Groundwater Contamination , 1991, SIAM J. Sci. Comput..

[87]  Y. K. Tang,et al.  Effect of Drain Diameter, Openings and Envelopes on Water Table Drawdown , 1979 .