Self‐Organized Permeability in Carbonate Aquifers

Advances over the past 40 years have resulted in a clear understanding of how dissolution processes in carbonate rocks enhance aquifer permeability. Laboratory experiments on dissolution rates of calcite and dolomite have established that there is a precipitous drop in dissolution rates as chemical equilibrium is approached. These results have been incorporated into numerical models, simulating the effects of dissolution over time and showing that it occurs along the entire length of pathways through carbonate aquifers. The pathways become enlarged and integrated over time, forming self‐organized networks of channels that typically have apertures in the millimeter to centimeter range. The networks discharge at point‐located springs. Recharge type is an important factor in determining channel size and distribution, resulting in a range of aquifer types, and this is well demonstrated by examples from England. Most carbonate aquifers have a large number of small channels, but in some cases large channels (i.e., enterable caves) can also develop. Rapid velocities found in ground water tracer tests, the high incidence of large‐magnitude springs, and frequent microbial contamination of wells all support the model of self‐organized channel development. A large majority of carbonate aquifers have such channel networks, where ground water velocities often exceed 100 m/d

Conceptualization and simulation of the Edwards aquifer, San Antonio region, Texas

A new numerical ground-water-flow model (Edwards aquifer model) that incorporates important components of the latest information and plausible conceptualization of the Edwards aquifer was developed. The model includes both the San Antonio and Barton Springs segments of the Edwards aquifer in the San Antonio region, Texas, and was calibrated for steady-state (1939–46) and transient (1947–2000) conditions, excluding Travis County. Transient simulations were conducted using monthly recharge and pumpage (withdrawal) data. The model incorporates conduits simulated as continuously connected (other than being separated in eastern Uvalde and southwestern Medina Counties), one-cell-wide (1,320 feet) zones with very large hydraulic-conductivity values (as much as 300,000 feet per day). The locations of the conduits were based on a number of factors, including major potentiometricsurface troughs in the aquifer, the presence of sinking streams, geochemical information, and geologic structures (for example, faults and grabens). The simulated directions of flow in the Edwards aquifer model are most strongly influenced by the presence of simulated conduits and barrier faults. The simulated flow in the Edwards aquifer is influenced by the locations of the simulated conduits, which tend to facilitate flow. The simulated subregional flow directions generally are toward the nearest conduit and subsequently along the conduits from the recharge zone into the confined zone and toward the major springs. Structures simulated in the Edwards aquifer model influencing groundwater flow that tend to restrict flow are barrier faults. The influence of simulated barrier faults on flow directions is most evident in northern Medina County. A water budget is an accounting of inflow to, outflow from, and storage change in the aquifer. For the Edwards aquifer model steady-state simulation, recharge (from seepage losses from streams and infiltration of rainfall) accounts for 93.5 percent of the sources of water to the Edwards aquifer, and inflow thr

Conceptualization and simulation of the Edwards aquifer, San Antonio region, Texas

A new numerical ground-water-flow model (Edwards aquifer model) that incorporates important components of the latest information and plausible conceptualization of the Edwards aquifer was developed. The model includes both the San Antonio and Barton Springs segments of the Edwards aquifer in the San Antonio region, Texas, and was calibrated for steady-state (1939–46) and transient (1947–2000) conditions, excluding Travis County. Transient simulations were conducted using monthly recharge and pumpage (withdrawal) data. The model incorporates conduits simulated as continuously connected (other than being separated in eastern Uvalde and southwestern Medina Counties), one-cell-wide (1,320 feet) zones with very large hydraulic-conductivity values (as much as 300,000 feet per day). The locations of the conduits were based on a number of factors, including major potentiometricsurface troughs in the aquifer, the presence of sinking streams, geochemical information, and geologic structures (for example, faults and grabens). The simulated directions of flow in the Edwards aquifer model are most strongly influenced by the presence of simulated conduits and barrier faults. The simulated flow in the Edwards aquifer is influenced by the locations of the simulated conduits, which tend to facilitate flow. The simulated subregional flow directions generally are toward the nearest conduit and subsequently along the conduits from the recharge zone into the confined zone and toward the major springs. Structures simulated in the Edwards aquifer model influencing groundwater flow that tend to restrict flow are barrier faults. The influence of simulated barrier faults on flow directions is most evident in northern Medina County. A water budget is an accounting of inflow to, outflow from, and storage change in the aquifer. For the Edwards aquifer model steady-state simulation, recharge (from seepage losses from streams and infiltration of rainfall) accounts for 93.5 percent of the sources of water to the Edwards aquifer, and inflow thr

Eogenetic Karst from the Perspective of an Equivalent Porous Medium

The porosity of young limestones experiencing meteoric diagenesis in the vicinity of their deposition (eogenetic karst) is mainly a double porosity consisting of touching-vug channels and preferred passageways lacing through a matrix of interparticle porosity. In contrast, the porosity of limestones experiencing subaerial erosion following burial diagenesis and uplift (telogenetic karst) is mainly a double porosity consisting of conduits within a network of fractures. The stark contrast between these two kinds of karst is illustrated by their position on a graph showing the hydraulic characteristics of an equivalent porous medium consisting of straight, cylindrical tubes (n-D space, where n is porosity,D is the diameter of the tubes, and logn is plotted against logD). Studies of the hydrology of small carbonate islands show that large-scale, horizontal hydraulic conductivity (K) increases by orders of magnitude during the evolution of eogenetic karst. Earlier petrologic studies have shown there is little if any change in the total porosity of the limestone during eogenetic diagenesis. The limestone of eogenetic karst, therefore, tracks horizontally inn-D space. In contrast, the path from initial sedimentary material to telogenetic karst comprises a descent on the graph with reduction ofn during burial diagenesis, then a sideways shift with increasingD due to opening of fractures during uplift and exposure, and finally an increase inD andn during development of the conduits along the fractures. Eogenetic caves are mainly limited to boundaries between geologic units and hydrologic zones: stream caves at the contact between carbonates and underlying impermeable rocks (and collapse-origin caves derived therefrom); vertical caves along platform-margin fractures; epikarst; phreatic pockets (banana holes) along the water table; and flank margin caves that form as mixing chambers at the coastal freshwater-saltwater “interface”. In contrast, the caverns of telogenetic karst are part of a system of interconnected conduits that drain an entire region. The eogenetic caves of small carbonate islands are, for the most part, not significantly involved in the drainage of the island