Can we simulate regional groundwater flow in a karst system using equivalent porous media models? Case study, Barton Springs Edwards aquifer, USA

Various approaches can be used to simulate groundwater flow in karst systems, including equivalent porous media distributed parameter, lumped parameter, and dual porosity approaches, as well as discrete fracture or conduit approaches. The purpose of this study was to evaluate two different equivalent porous media approaches: lumped and distributed parameter, for simulating regional groundwater flow in a karst aquifer and to evaluate the adequacy of these approaches. The models were applied to the Barton Springs Edwards aquifer, Texas. Unique aspects of this study include availability of detailed information on recharge from stream-loss studies and on synoptic water levels, long-term continuous water level monitoring in wells throughout the aquifer, and spring discharge data to compare with simulation results. The MODFLOW code was used for the distributed parameter model. Estimation of hydraulic conductivity distribution was optimized by using a combination of trial and error and automated inverse methods. The lumped parameter model consists of five cells representing each of the watersheds contributing recharge to the aquifer. Transient simulations were conducted using both distributed and lumped parameter models for a 10-yr period (1989 ‐ 1998). Both distributed and lumped parameter models fairly accurately simulated the temporal variability in spring discharge; therefore, if the objective of the model is to simulate spring discharge, either distributed or lumped parameter approaches can be used. The distributed parameter model generally reproduced the potentiometric surface at different times. The impact of the amount of pumping on a regional scale on spring discharge can be evaluated using a lumped parameter model; however, more detailed evaluation of the effect of pumping on groundwater levels and spring discharge requires a distributed parameter modeling approach. Sensitivity analyses indicated that spring discharge was much more sensitive to variations in recharge than pumpage, indicating that aquifer management should consider enhanced recharge, in addition to conservation measures, to maintain spring flow. This study shows the ability of equivalent porous media models to simulate regional groundwater flow in a highly karstified aquifer, which is important for water resources and groundwater management. q 2003 Elsevier Science B.V. All rights reserved.

[1]  AN EXTENDED DOUBLE-POROSITY CONCEPT AS A PRACTICAL MODELLING APPROACH FOR A KARSTIFIED TERRAIN , 1993 .

[2]  C. Wicks,et al.  Modeling of storm responses in conduit flow aquifers with reservoirs , 1998 .

[3]  Michael E. Barrett,et al.  A parsimonious model for simulation of flow and transport in a karst aquifer , 1996 .

[4]  J. Pankow,et al.  An evaluation of contaminant migration patterns at two waste disposal sites on fractured porous media in terms of the equivalent porous medium (EPM) model , 1986 .

[5]  J. L. Robinson,et al.  Descriptions of anisotropy and heterogeneity and their effect on ground-water flow and areas of contribution to public supply wells in a karst carbonate aquifer system , 1996 .

[6]  Michael E. Barrett,et al.  A parsimonious model for simulating flow in a karst aquifer , 1997 .

[7]  R. M. Slade,et al.  Simulation of the flow system of Barton Springs and associated Edwards Aquifer in the Austin area, Texas , 1985 .

[8]  S. P. Neuman,et al.  Stochastic continuum representation of fractured rock permeability as an alternative to the REV and fracture network concepts , 1988 .

[9]  Arlen W. Harbaugh,et al.  User's documentation for MODFLOW-96, an update to the U.S. Geological Survey modular finite-difference ground-water flow model , 1996 .

[10]  A. Shapiro,et al.  Hydrogeologic characterization of the Minnelusa and Madison aquifers near Spearfish, South Dakota , 1999 .

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

[12]  Er-Wei Bai,et al.  Simulation Of Spring Discharge From A Limestone Aquifer In Iowa, USA , 1996 .

[13]  Marie Larocque,et al.  Determining Karst Transmissivities with Inverse Modeling and an Equivalent Porous Media , 1999 .

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

[15]  S. Hovorka,et al.  Permeability Structure of the Edwards Aquifer, South Texas: Implications for Aquifer Management , 1998 .

[16]  P. Ryder.  Hydrology of the Floridan Aquifer system in west-central Florida , 1985 .

[17]  P. A. Hsieh,et al.  Documentation of a computer program to simulate horizontal-flow barriers using the U.S. Geological Survey's modular three-dimensional finite-difference ground-water flow model , 1993 .

[18]  C. Kreitler,et al.  Hydrogeology of the Edwards Aquifer, Austin Area, Central Texas , 1984 .

[19]  Chin-Fu Tsang,et al.  Flow and Contaminant Transport in Fractured Rock , 1993 .

[20]  E. Greene.  Hydraulic properties of the Madison aquifer system in the western Rapid City area, South Dakota , 1993 .

[21]  J. E. Warren,et al.  The Behavior of Naturally Fractured Reservoirs , 1963 .

[22]  M. Sauter.  DOUBLE POROSITY MODELS IN KARSTIFIED LIMESTONE AQUIFERS: FIELD VALIDATION AND DATA PROVISION , 1993 .

[23]  Nico M. Hauwert,et al.  Geologic framework and hydrogeologic characteristics of the Edwards aquifer outcrop (Barton Springs segment), northeastern Hays and southwestern Travis Counties, Texas , 1996 .

[24]  Andrew V. Wolfsberg,et al.  Rock Fractures and Fluid Flow: Contemporary Understanding and Applications , 1997 .

[25]  W. Dragoni,et al.  The Problem of Modeling Limestone Springs: The Case of Bagnara (North Apennines, Italy) , 1997 .

[26]  Roberto Anaya,et al.  A Lumped Parameter Model for the Edwards Aquifer , 1993 .

[27]  W. Kinzelbach.  Applied groundwater modeling — Simulation of flow and advective transport , 1992 .

[28]  T. J. Larkin,et al.  Climatic atlas of Texas , 1983 .