Simulating coupled surface and subsurface water flow in a tile-drained agricultural catchment

Summary The impact of shallow tile drainage networks on groundwater flow patterns, and associated transport of chemical species or sediment particles, needs to be quantified to evaluate the effects of agricultural management on water resources. A current challenge is to represent tile drainage networks in numerical models at the scale of a catchment for which proper simulation and reproduction of subsurface drainage water dynamics are essential. This study investigates the applicability of various tile drainage modeling concepts by comparing their results to a reference model. The models were setup for a 3.5 ha regularly monitored tile-drained area located in the clayey till Lillebaek catchment in Denmark, where the main input of water to the streams originates from subsurface drainage. Several simulations helped identify the most efficient way of modeling tile-drained areas using the integrated surface water and groundwater HydroGeoSphere code, which also simulates one-dimensional water flow in tile drains. The aim was also to provide insight on the design of larger scale models of complex tile drainage systems, for which computational time can become very large. A reference model that simulated coupled surface flow, groundwater flow and flow in tile drains was setup and showed a rapid response in drain outflow when precipitation is applied. A series of additional simulations were performed to test the influence of flow boundary conditions, temporal resolution of precipitation data, conceptual representations of clay tills and drains, as well as spatial discretization of the mesh. For larger scale models, the simulations suggest that a simplification of the geometry of the drainage network is a suitable option for avoiding highly discretized meshes. Representing the drainage network by a high-conductivity porous medium layer reproduces outflow simulated by the reference model. This approach greatly reduces the mesh complexity and simulation time and thus represents an alternative to a discrete representation of drains at the field scale, but specifying the hydraulic properties of the layer requires calibration against observed drain discharge.

[1]  P. Rasmussen Monitoring shallow groundwater quality in agricultural watersheds in Denmark , 1996 .

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

[3]  R. Sidle,et al.  Mass transport and scale-dependent hydraulic tests in a heterogeneous glacial till–sandy aquifer system , 2001 .

[4]  M. Vauclin,et al.  A simplified modelling approach for pesticide transport in a tile-drained field: The PESTDRAIN model , 2009 .

[5]  G. Sands,et al.  Effects of Agricultural Drainage on Aquatic Ecosystems: A Review , 2009 .

[6]  K. Jensen,et al.  Importance of including small‐scale tile drain discharge in the calibration of a coupled groundwater‐surface water catchment model , 2013 .

[7]  L. B. Hansen,et al.  Pesticides in Streams and Subsurface Drainage Water within Two Arable Catchments in Denmark: Pesticide Application, Concentration, Transport and Fate , 2003 .

[8]  Jeffrey G. Arnold,et al.  DEVELOPMENT AND APPLICATION OF SWAT TO LANDSCAPES WITH TILES AND POTHOLES , 2005 .

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

[10]  D. Durnford,et al.  Drainage Design for Salinity Control , 2015 .

[11]  M. Voltz,et al.  MHYDAS-DRAIN : A spatially distributed model for small, artificially drained lowland catchments , 2007 .

[12]  Nicola Fohrer,et al.  Incorporating landscape depressions and tile drainages of a northern German lowland catchment into a semi‐distributed model , 2010 .

[13]  N. Fausey,et al.  Effect of subsurface drainage on streamflow in an agricultural headwater watershed , 2014 .

[14]  M. V. Genuchten,et al.  Two-dimensional modelling of preferential water flow and pesticide transport from a tile-drained field , 2006 .

[15]  Antti Taskinen,et al.  Simulation of water balance in a clayey, subsurface drained agricultural field with three-dimensional FLUSH model , 2013 .

[16]  C. Simmons,et al.  HydroGeoSphere: A Fully Integrated, Physically Based Hydrological Model , 2012 .

[17]  S. Frey,et al.  Viscosity dependent dual‐permeability modeling of liquid manure movement in layered, macroporous, tile drained soil , 2012 .

[18]  L. McKay,et al.  Aquifer Vulnerability to Pesticide Migration Through Till Aquitards , 2004, Ground water.

[19]  Feng Yue,et al.  Fully-coupled physically-based approach for modeling conventional and managed subsurface drainage , 2010 .

[20]  R. W. Skaggs,et al.  A water management model for shallow water table soils , 1978 .

[21]  J. Rozemeijer,et al.  Integrated modeling of groundwater–surface water interactions in a tile‐drained agricultural field: The importance of directly measured flow route contributions , 2010 .

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

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

[24]  Peter A. Forsyth,et al.  A parallel computational framework to solve flow and transport in integrated surface-subsurface hydrologic systems , 2012, Environ. Model. Softw..

[25]  Philippe Renard,et al.  Modeling Fine‐Scale Geological Heterogeneity—Examples of Sand Lenses in Tills , 2013, Ground water.

[26]  Nadia Carluer,et al.  Assessment and modelling of the influence of man-made networks on the hydrology of a small watershed: implications for fast flow components, water quality and landscape management , 2004 .

[27]  Andrea Viezzoli,et al.  Nitrate reduction in geologically heterogeneous catchments--a framework for assessing the scale of predictive capability of hydrological models. , 2014, The Science of the total environment.

[28]  P. Gravesen,et al.  Fractures and biopores in weichselian clayey till aquitards at Flakkebjerg, Denmark , 1999 .

[29]  M. Mccabe,et al.  Assessing the impact of model spin‐up on surface water‐groundwater interactions using an integrated hydrologic model , 2012 .

[30]  M. Schaap,et al.  ROSETTA: a computer program for estimating soil hydraulic parameters with hierarchical pedotransfer functions , 2001 .

[31]  Garey A. Fox,et al.  Numerical Simulation of Flow Dynamics during Macropore–Subsurface Drain Interactions Using HYDRUS , 2008 .

[32]  K. Jensen,et al.  Numerical analysis of water and solute transport in variably-saturated fractured clayey till. , 2009, Journal of contaminant hydrology.

[33]  U. Sunday Tim,et al.  Coupling Vadose Zone Models with GIS: Emerging Trends and Potential Bottlenecks , 1996 .

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

[35]  M. Houmark‐Nielsen,et al.  Southwest Scandinavia, 40–15 kyr BP: palaeogeography and environmental change , 2003 .

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

[37]  I. Ginzburg,et al.  Field-scale modeling of subsurface tile-drained soils using an equivalent-medium approach , 2007 .

[38]  W. Wallender,et al.  Describing Near Surface, Transient Flow Processes in Unconfined Aquifers below Irrigated Lands: Model Application in the Western San Joaquin Valley, California , 2004 .

[39]  John R. Williams,et al.  LARGE AREA HYDROLOGIC MODELING AND ASSESSMENT PART I: MODEL DEVELOPMENT 1 , 1998 .

[40]  Erwin Zehe,et al.  A novel explicit approach to model bromide and pesticide transport in connected soil structures , 2011 .

[41]  J. Ayars,et al.  Controlled drainage for improved water management in arid regions irrigated agriculture , 2006 .

[42]  Andreas Bauwe,et al.  Application of the SWAT Model for a Tile-Drained Lowland Catchment in North-Eastern Germany on Subbasin Scale , 2013, Water Resources Management.

[43]  R. Therrien,et al.  Channel Representation in Physically Based Models Coupling Groundwater and Surface Water: Pitfalls and How to Avoid Them , 2014, Ground water.

[44]  B. Lesaffre,et al.  Subsurface Drainage Peak Flows in Shallow Soil , 1988 .

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