Identifying three-dimensional nested groundwater flow systems in a Tóthian basin

Abstract Nested groundwater flow systems have been revealed in Toth's theory as the structural property of basin-scale groundwater circulation but were only well known with two-dimensional (2D) profile models. The method of searching special streamlines across stagnation points for partitioning flow systems, which has been successfully applied in the 2D models, has never been implemented for three-dimensional (3D) Tothian basins because of the difficulty in solving the dual stream functions. Alternatively, a new method is developed to investigate 3D nested groundwater flow systems without determination of stagnation points. Connective indices are defined to quantify the connection between individual recharge and discharge zones along streamlines. Groundwater circulation cells (GWCCs) are identified according to the distribution of the connective indices and then grouped into local, intermediate and regional flow systems. This method requires existing solution of the flow velocity vector and is implemented via particle tracking technique. It is applied in a hypothetical 3D Tothian basin with an analytical solution of the flow field and in a real-world basin with a numerical modeling approach. Different spatial patterns of flow systems compared to 2D profile models are found. The outcrops boundaries of GWCCs on water table may significantly deviate from and are not parallel to the nearby water table divides. Topological network is proposed to represent the linked recharge–discharge zones through closed and open GWCCs. Sensitivity analysis indicates that the development of GWCCs depends on the basin geometry, hydraulic parameters and water table shape.

[1]  David W. Pollock,et al.  User's guide for MODPATH/MODPATH-PLOT, Version 3; a particle tracking post-processing package for MODFLOW, the U.S. Geological Survey finite-difference ground-water flow model , 1994 .

[2]  John David Wilson,et al.  MODFLOW-2000, The U.S. Geological Survey Modular Ground-Water Model -- GMG Linear Equation Solver Package Documentation , 2004 .

[3]  O. Batelaan,et al.  Regional groundwater discharge: phreatophyte mapping, groundwater modelling and impact analysis of land-use change , 2003 .

[4]  Xu-sheng Wang,et al.  On the use of late-time peaks of residence time distributions for the characterization of hierarchically nested groundwater flow systems , 2016 .

[5]  Hailong Li,et al.  An analytical study on groundwater flow in drainage basins with horizontal wells , 2014, Hydrogeology Journal.

[6]  A. Wörman,et al.  The use of spectral analysis-based exact solutions to characterize topography-controlled groundwater flow , 2011 .

[7]  D. Solomon,et al.  An integrated environmental tracer approach to characterizing groundwater circulation in a mountain block , 2005 .

[8]  Series solutions for seepage in three dimensional aquifers , 2005 .

[9]  P. Goderniaux,et al.  Groundwater flow systems theory: research challenges beyond the specified-head top boundary condition , 2016, Hydrogeology Journal.

[10]  S. Ge,et al.  An analytical study on stagnation points in nested flow systems in basins with depth‐decaying hydraulic conductivity , 2011 .

[11]  S. Ge,et al.  A quantitative study on accumulation of age mass around stagnation points in nested flow systems , 2012 .

[12]  W. Zijl Scale aspects of groundwater flow and transport systems , 1999 .

[13]  Hiroshi Kashiwagi,et al.  Visualization techniques for 3D vector fields: an application to electrostatic fields of molecules , 2001, Comput. Animat. Virtual Worlds.

[14]  M. Perk,et al.  The stability of groundwater flow systems in unconfined sandy aquifers in the Netherlands , 2008 .

[15]  J. Harvey,et al.  Exact three‐dimensional spectral solution to surface‐groundwater interactions with arbitrary surface topography , 2006 .

[16]  M. Jin,et al.  Numerical simulation of groundwater flow patterns using flux as upper boundary , 2013 .

[17]  Xu-sheng Wang,et al.  Semi-empirical equations for the systematic decrease in permeability with depth in porous and fractured media , 2010 .

[18]  T. Gleeson,et al.  Regional groundwater flow in mountainous terrain: Three‐dimensional simulations of topographic and hydrogeologic controls , 2008 .

[19]  F. W. Schwartz,et al.  Fundamentals of Ground Water , 2002 .

[20]  M. Anderson,et al.  Seasonal reversals of groundwater flow around lakes and the relevance to stagnation points and lake budgets , 1981 .

[21]  C. Duffy,et al.  Groundwater circulation in a closed desert basin: Topographic scaling and climatic forcing , 1988 .

[22]  Xiao-Wei Jiang,et al.  A new analytical solution of topography‐driven flow in a drainage basin with depth‐dependent anisotropy of permeability , 2011 .

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

[24]  P. Goderniaux,et al.  Partitioning a regional groundwater flow system into shallow local and deep regional flow compartments , 2013 .

[25]  R. Nativ,et al.  Evidence for Ground‐Water Circulation in the Brine‐Filled Aquitard, Oak Ridge, Tennessee , 1997 .

[26]  Numerical simulations based on stream functions and velocities in three—Dimensional groundwater flow , 1986 .

[27]  Xingcheng Lu,et al.  Direct observation of complex Tóthian groundwater flow systems in the laboratory , 2010 .

[28]  J. Tóth A Theoretical Analysis of Groundwater Flow in Small Drainage Basins , 1963 .

[29]  J. Tóth Groundwater as a geologic agent: An overview of the causes, processes, and manifestations , 1999 .

[30]  G. Matanga,et al.  Stream functions in three‐dimensional groundwater flow , 1993 .

[31]  A. Love,et al.  Hidden channels of groundwater flow in Tóthian drainage basins , 2013 .

[32]  M. Cardenas Potential contribution of topography‐driven regional groundwater flow to fractal stream chemistry: Residence time distribution analysis of Tóth flow , 2007 .

[33]  Xiao-Wei Jiang,et al.  Field identification of groundwater flow systems and hydraulic traps in drainage basins using a geophysical method , 2014 .