Description and evaluation of numerical groundwater flow models for the Edwards Aquifer, south-central Texas

Richard J. Lindgren; Charles J. Taylor; Natalie A. Houston

doi:10.3133/sir20095183

Home/Archive/USGS

USGS · sir20095183

Description and evaluation of numerical groundwater flow models for the Edwards Aquifer, south-central Texas

Richard J. Lindgren·Charles J. Taylor·Natalie A. Houston

View original source ↗Export metadata ↓

Abstract

A substantial number of public water system wells in south-central Texas withdraw groundwater from the karstic, highly productive Edwards aquifer. However, the use of numerical groundwater flow models to aid in the delineation of contributing areas for public water system wells in the Edwards aquifer is problematic because of the complex hydrogeologic framework and the presence of conduit-dominated flow paths in the aquifer. The U.S. Geological Survey, in cooperation with the Texas Commission on Environmental Quality, evaluated six published numerical groundwater flow models (all deterministic) that have been developed for the Edwards aquifer San Antonio segment or Barton Springs segment, or both. This report describes the models developed and evaluates each with respect to accessibility and ease of use, range of conditions simulated, accuracy of simulations, agreement with dye-tracer tests, and limitations of the models. These models are (1) GWSIM model of the San Antonio segment, a FORTRAN computer-model code that pre-dates the development of MODFLOW; (2) MODFLOW conduit-flow model of San Antonio and Barton Springs segments; (3) MODFLOW diffuse-flow model of San Antonio and Barton Springs segments; (4) MODFLOW Groundwater Availability Modeling [GAM] model of the Barton Springs segment; (5) MODFLOW recalibrated GAM model of the Barton Springs segment; and (6) MODFLOW–DCM (dual conductivity model) conduit model of the Barton Springs segment. The GWSIM model code is not commercially available, is limited in its application to the San Antonio segment of the Edwards aquifer, and lacks the ability of MODFLOW to easily incorporate newly developed processes and packages to better simulate hydrologic processes. MODFLOW is a widely used and tested code for numerical modeling of groundwater flow, is well documented, and is in the public domain. These attributes make MODFLOW a preferred code with regard to accessibility and ease of use. The MODFLOW conduit-flow model incorporates improvements over previous models by using (1) a user-friendly interface, (2) updated computer codes (MODFLOW–96 and MODFLOW–2000), (3) a finer grid resolution, (4) less-restrictive boundary conditions, (5) an improved discretization of hydraulic conductivity, (6) more accurate estimates of pumping stresses, (7) a long transient simulation period (54 years, 1947–2000), and (8) a refined representation of high-permeability zones or conduits. All of the models except the MODFLOW–DCM conduit model have limitations resulting from the use of Darcy’s law to simulate groundwater flow in a karst aquifer system where non-Darcian, turbulent flow might actually dominate. The MODFLOW–DCM conduit model is an improvement in the ability to simulate karst-like flow conditions in conjunction with porous-media-type matrix flow. However, the MODFLOW–DCM conduit model has had limited application and testing and currently (2008) lacks commercially available pre- and post-processors. The MODFLOW conduit-flow and diffuse-flow Edwards aquifer models are limited by the lack of calibration for the northern part of the Barton Springs segment (Travis County) and their reliance on the use of the calibrated hydraulic conductivity and storativity values from the calibrated Barton Springs segment GAM model. The major limitation of the Barton Springs segment GAM and recalibrated GAM models is that they were calibrated to match measured water levels and springflows for a restrictive range of hydrologic conditions, with each model having different hydraulic conductivity and storativity values appropriate to the hydrologic conditions that were simulated. The need for two different sets of hydraulic conductivity and storativity values increases the uncertainty associated with the accuracy of either set of values, illustrates the non-uniqueness of the model solution, and probably most importantly demonstrates the limitations of using a one-layer model to represent the heterogeneous hydrostratigraphic units composing the Edwards aquifer. In general, the best matches or agreement between groundwater flow directions inferred by numerical model simulation, and by dye-tracer tests, are observed where model outputs accurately reproduce the configuration of the potentiometric surface with regard to the positions of major and minor groundwater troughs and divides. None of the models, with the possible exception of the MODFLOW–DCM conduit model, has a documented capability to accurately simulate travel times for conduit-dominated velocities in the Edwards aquifer. Public water system assessments of wells in the Barton Springs segment of the Edwards aquifer, and elsewhere where conduit-flow conditions are thought to dominate aquifer hydraulic behavior, might be enhanced by use of either the MODFLOW–DCM model or the newly developed U.S. Geological Survey MODFLOW Conduit-Flow Process module for MODFLOW–2005, because each incorporates a type of dual or triple hydraulic conductivity approach and has the capability to explicitly simulate turbulent flow and conduit hydraulic characteristics.

Explore related subjects

Scientific Investigations Report Report

Source-reported bounding extent: 28° to 31° latitude; -101° to -97° longitude. This indicates report coverage, not an exact sampling location. View area on OpenStreetMap.

Keep this discovery

Explore connections, maps & timelines

BibTeX RIS

Richard J. Lindgren, Charles J. Taylor, Natalie A. Houston. 2009. Description and evaluation of numerical groundwater flow models for the Edwards Aquifer, south-central Texas. https://doi.org/10.3133/sir20095183

Cite the original work for its findings. Save a collection to share your selection of sources.

KEEP EXPLORING

Related USGS reports

USGS2026

Monitoring recreation on federally managed lands and waters—Aspects of visitor use

Federally managed public lands and waters receive around 1 billion recreational visits each year. Data on these visitors can aid in guiding policy decisions, managing resources effectively, and communicating the economic contributions of lands and waters. This report explores the methods used by agencies to collect data on aspects of recreational visitor use to Federal lands and waters (apart from visitation numbers, which are the focus of a companion publication). Aspects of recreational visitor use include visitor demographics, recreational activity participation, visitor satisfaction, visitor attitudes and experiences, trip characteristics, and economic contributions. We review practices used to understand aspects of visitor use across seven Federal agencies, revealing similarities such as the use of visitor intercept surveys and coverage of similar topic areas, and differences in how survey programs are operationalized and how specific questions on visitor surveys are worded. We also evaluate emerging technologies, such as geolocated social media and mobile device location data, for their potential to aid in understanding aspects of visitor use. This report concludes with potential opportunities to enhance data collection and coordination, ensuring cost-effective data collection and informed decision making.

Scientific Investigations Report↗

USGS2026

Evaluation of best management practices at an edge-of-field site in the Eagle Creek watershed, Ohio, 2012–20

In 2010, the U.S. Geological Survey worked in partnership with the Great Lakes Restoration Initiative and the Natural Resources Conservation Service to identify farm fields in priority watersheds—watersheds critically important to Great Lake health. Intensive best management practices were implemented in five States to test the efficacy of best management practices, also referred to as agricultural conservation practices, to reduce sediment and nutrient runoff on agricultural fields. These fields were chosen as representatives of the priority watersheds because their farming practices and geographical conditions were common among farms in those watersheds. A farm site located near Findlay, Ohio, in the Eagle Creek watershed used two conservation practices, cover crop (Natural Resources Conservation Service practice standard 340) and variable rate technology (VRT; Natural Resources Conservation Service practice standard 590), which represent common practices in the Maumee priority watershed. This study monitored surface runoff and subsurface tile runoff at the Eagle Creek watershed site. The effects of cover crop and VRT nutrient application on sediment and nutrient runoff were assessed from October 2012 to September 2020. Cover crops were applied in fall of 2016, 2017, and 2018, coinciding with VRT applications. The following parameters were analyzed as part of the study: runoff discharge; peak discharge; and concentrations and loads of total phosphorus, particulate phosphorus, dissolved reactive phosphorus, and suspended sediment. High peak runoff events disproportionally affected nutrient concentrations and loads in runoff. A threshold of the 85th percentile of the peak discharge was chosen to assess the effect of cover crop and VRT nutrient application on nutrient mitigation at each gage. Runoff events below the 85th percentile were considered small, and runoff events over the 85th percentile were considered large. Concentrations from collected water-quality samples were used for the analysis; estimated concentrations were excluded from the statistical evaluation. All surface runoff parameters during small runoff events, except surface runoff volume, were significantly lower in the period with cover crops and VRT implementation. Contrarily, subsurface runoff parameters during small events were not statistically different between the periods with and without cover crops and VRT implementation. Surface and subsurface parameters during large events were also not statistically different between the periods with and without cover crops and VRT implementation. These results indicate cover crops and VRT may improve water quality during small runoff events, but additional best management practices that mitigate large runoff events may lead to greater water-quality improvement given the contribution of large events to overall losses.

Ohio↗

USGS2026

Continuous monitoring and temporal variability of fluorescence of dissolved organic matter in karst groundwater of the Edwards aquifer, south-central Texas, 2019–24

Widespread urbanization on the Edwards aquifer recharge zone in south-central Texas has prompted the development of new approaches and tools for evaluating the current and future status of water quality in the San Antonio and Barton Springs segments of the Edwards aquifer. The U.S. Geological Survey, in cooperation with the San Antonio Water System and the City of Austin, applied continuous monitoring of fluorescence of dissolved organic matter (fDOM) to characterize the sources and transport of dissolved organic matter and associated constituents in the Edwards aquifer. Continuous fDOM time-series data were adjusted using empirically derived temperature, turbidity, and inner-filter-effect correction factors prior to analysis and interpretation. Fully corrected fDOM time-series data from the study sites were evaluated for short- and long-term temporal variability in fDOM in the context of the hydrologic and climatic conditions that occurred during the 5-year study period to understand vulnerability of the aquifer to potential contaminants in recharge. Study results indicate that fDOM is an effective proxy for dissolved organic carbon that, in turn, is indicative of the influx of recent surface water and associated contaminants, and thus, aquifer vulnerability to those contaminants. Although other continuously monitored water-quality parameters provide insights into changes in water quality in response to varying hydrologic conditions, the quantitative relation between fDOM and dissolved organic carbon directly indicates the timing and magnitude of pulses of organic constituents derived from the land surface, whereas other water-quality parameters do not.

Texas↗