Effects of turbulence on hydraulic heads and parameter sensitivities in preferential groundwater flow layers

[1] A conduit flow process (CFP) for the Modular Finite Difference Ground-Water Flow model, MODFLOW-2005, has been created by the U.S. Geological Survey. An application of the CFP on a carbonate aquifer in southern Florida is described; this application examines (1) the potential for turbulent groundwater flow and (2) the effects of turbulent flow on hydraulic heads and parameter sensitivities. Turbulent flow components were spatially extensive in preferential groundwater flow layers, with horizontal hydraulic conductivities of about 5,000,000 m d−1, mean void diameters equal to about 3.5 cm, groundwater temperature equal to about 25°C, and critical Reynolds numbers less than or equal to 400. Turbulence either increased or decreased simulated heads from their laminar elevations. Specifically, head differences from laminar elevations ranged from about −18 to +27 cm and were explained by the magnitude of net flow to the finite difference model cell. Turbulence also affected the sensitivities of model parameters. Specifically, the composite-scaled sensitivities of horizontal hydraulic conductivities decreased by as much as 70% when turbulence was essentially removed. These hydraulic head and sensitivity differences due to turbulent groundwater flow highlight potential errors in models based on the equivalent porous media assumption, which assumes laminar flow in uniformly distributed void spaces.

[1]  Edward R. German,et al.  Regional evaluation of evapotranspiration in the Everglades , 1996 .

[2]  Mary C. Hill,et al.  PRECONDITIONED CONJUGATE-GRADIENT 2 (PCG2), a computer program for solving ground-water flow equations , 1990 .

[3]  Georg Teutsch,et al.  Hydraulic boundary conditions as a controlling factor in karst genesis: A numerical modeling study on artesian conduit development in gypsum , 2003 .

[4]  J. Bear.  Dynamics of Fluids in Porous Media , 1975 .

[5]  M. L. Merritt.  Simulation of the water-table altitude in the Biscayne aquifer, southern Dade County, Florida, water years 1945-89 , 1995 .

[6]  W. Dreybrodt.  Processes in Karst Systems , 1988 .

[7]  Georg Teutsch,et al.  Simulation of the development of karst aquifers using a coupled continuum pipe flow model , 2003 .

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

[9]  Sebastian Bauer,et al.  Documentation of a Conduit Flow Process (CFP) for MODFLOW-2005 , 2007 .

[10]  Michael E. Barrett,et al.  Can we simulate regional groundwater flow in a karst system using equivalent porous media models? Case study, Barton Springs Edwards aquifer, USA , 2003 .

[11]  Mary C. Hill,et al.  Methods and Guidelines for Effective Model Calibration , 2000 .

[12]  R. Price,et al.  Use of tritium and helium to define groundwater flow conditions in Everglades National Park , 2003 .

[13]  S. Birk.  Characterisation of Karst Systems by Simulating Aquifer Genesis and Spring Responses : Model Development and Application to Gypsum Karst , 2002 .

[14]  E. Poeter,et al.  Documentation of UCODE; a computer code for universal inverse modeling , 1998 .

[15]  H. Vacher,et al.  Eogenetic karst from the perspective of an equivalent porous medium , 2002, Carbonates and Evaporites.

[16]  A. W. Harbaugh.  MODFLOW-2005 : the U.S. Geological Survey modular ground-water model--the ground-water flow process , 2005 .

[17]  C. Voss,et al.  Inverse modeling for seawater intrusion in coastal aquifers: Insights about parameter sensitivities, variances, correlations and estimation procedures derived from the Henry problem , 2006 .

[18]  Paul Williams,et al.  Karst Geomorphology and Hydrology , 1989 .

[19]  Michael A. Wacker,et al.  Assessing the Vulnerability of a Municipal Well Field to Contamination in a Karst Aquifer , 2005 .

[20]  R. Price,et al.  Identification of a hydrodynamic threshold in karst rocks from the Biscayne Aquifer, south Florida, USA , 2008 .

[21]  Sebastian Bauer,et al.  Modeling of karst aquifer genesis: Influence of exchange flow , 2003 .

[22]  P. Jeannin.  Modeling flow in phreatic and epiphreatic Karst conduits in the Hölloch Cave (Muotatal, Switzerland) , 2001 .

[23]  G. Schneebeli.  EXPÉRIENCES SUR LA LIMITE DE VALIDITÉ DE LA LOI DE DARCY ET L'APPARITION DE LA TURBULENCE DANS UN ÉCOULEMENT DE FILTRATION , 1955 .

[24]  W. Dreybrodt.  Processes in karst systems : physics, chemistry, and geology , 1988 .

[25]  Michael A. Wacker,et al.  Application of carbonate cyclostratigraphy and borehole geophysics to delineate porosity and preferential flow in the karst limestone of the Biscayne aquifer, SE Florida , 2006 .

[26]  J. Bear.  Hydraulics of Groundwater , 1979 .

[27]  John King Vennard,et al.  Elementary Fluid Mechanics , 1940 .

[28]  Michael A. Wacker,et al.  A cyclostratigraphic and borehole-geophysical approach to development of a three-dimensional conceptual hydrogeologic model of the karstic Biscayne aquifer, southeastern Florida , 2006 .