USGS Science⌕ Search

USGS · 70012577

Effects of karst and geologic structure on the circulation of water and permeability in carbonate aquifers

Abstract

The results of the natural processes caused by solution and leaching of limestone, dolomite, gypsum, salt and other soluble rocks, is known as karst. Development of karst is commonly known as karstification, which may have a pronounced effect on the topography, hydrology and environment, especially where such karst features as sinkholes and vertical solution shafts extend below the land surface and intersect lateral solution passages, cavities, caverns and other karst features in carbonate rocks. Karst features may be divided into two groups: (1) surficial features that do not extend far below the surface; and (2) karst features such as sinkholes that extend below the surface and affect the circulation of water below. The permeability of the most productive carbonate aquifers is due chiefly to enlargement of fractures and other openings by circulation of water. Important controlling factors responsible for the development of karst and permeability in carbonate aquifers include: (1) climate, topography, and presence of soluble rocks; (2) geologic structure; (3) nature of underground circulation; and (4) base level. Another important factor is the condition of the surface of the carbonate rocks at the time they are exposed to meteoric water. A carbonate rock surface, with soil or relatively permeable, less soluble cover, is more favorable for initiation of karstification and solution than bare rocks. Water percolates downward through the cover to the underlying carbonate rocks instead of running off on the surface. Also, the water becomes more corrosive as it percolates through the permeable cover to the underlying carbonate rocks. Where there is no cover or the cover has been removed, the carbonate rocks become case hardened and resistant to erosion. However, in regions underlain not only by carbonate rocks but also by beds of anhydrite, gypsum and salt, such as the Hueco Plateau in southeastern New Mexico, subsurface solution may occur where water without natural acids moves down from bare rock surfaces through cracks to the beds that are more soluble than carbonate rocks. For example, in the area of Carlsbad Caverns in southeastern New Mexico, much of the water responsible for solution that formed the caverns apparently entered the groundwater system through large open fractures and did not form sinkhole topography. East of the Carlsbad Caverns, however, in the Pecos River Valley where the carbonate rocks are overlain by the less soluble Ogallala Formation of Late Tertiary age, solution began along escarpments as the Pecos River and its tributaries cut through the less soluble cover. As these escarpments retreated, sinkholes and other karst features developed. Joints or fractures are essential for initiation of downward percolation of water in compact carbonate rocks such as some Paleozoic limestone in which there is no intergranular permeability. Also joints or fractures and bedding planes may be essential in the initiation of lateral movement of water in the zone of saturation. Where conditions of recharge and discharge are favorable, groundwater may move parallel to the dip. However, the direction of movement of water in most carbonate rocks is not necessarily down dip or parallel to the dip. The general direction of movement of both surface and groundwater may be parallel to the strike in a breached anticline. Faults may restrict the lateral movement of water, especially if water-bearing beds are faulted against relatively impervious beds. Conversely, some fault may serve as avenues through which water may move as, for example, in the Cretaceous Edwards aquifer in the San Antonio area, Texas. Karst aquifers, chiefly carbonate rocks, may be placed in three groups according to water-bearing capacity. Water in aquifers of group 1 occurs chiefly in joints, fractures, and other openings that have not been enlarged by solution. The yield of wells is small. Aquifers in group 2 , with low to intermediate yields, are those in which water occurs in joints and fractures with some cavities and channels enlarged by solution. Aquifers in group 3 are those in which the yield of wells and springs range from intermediate to very large. This group includes five of the most productive aquifers in the U.S.A. The water-bearing beds of all of these productive aquifers, except the Biscayne aquifer in southeastern Florida, contain buried paleokarst in which the permeability has been reactivated and enlarged by the present circulation system.

Explore related subjects

Keep this discovery

Explore connections, maps & timelines

BibTeXRIS

V. T. Stringfield, J. R. Rapp, R.B. Anders. 1979. Effects of karst and geologic structure on the circulation of water and permeability in carbonate aquifers. https://doi.org/10.1016/0022-1694(79)90178-1

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

KEEP EXPLORING

Related USGS reports

Inverse modeling for post-wildfire hydrologic process understanding over snow-dominated watersheds

Wildfire-induced changes to vegetation and soil properties can alter snowpack dynamics and hydrologic processes in ways that complicate predictions of post-wildfire water availability. To better understand these effects, we coupled the Parameter ESTimation in C++ (PEST++) software with the Weather Research & Forecasting Hydrologic model (WRF-Hydro) across three burned, snow-dominated watersheds in the western United States. We incorporated three burn-severity-based land cover categories into WRF-Hydro and used an iterative Ensemble Smoother (iES) within PEST++ to estimate model parameters controlling snow accumulation and melt, infiltration, and streamflow generation. This framework supported three complementary objectives: (1) assessing parameter sensitivities under pre- and post-wildfire conditions, (2) comparing pre- and post-wildfire parameter ensembles to identify wildfire-driven shifts in hydrologic processes, and (3) using counterfactual simulations to isolate the hydrologic response to wildfire from interannual climate variability. Our findings indicate that model parameter sensitivity varied by domain and shifted following wildfire. Maximum snow albedo and snow albedo decay were most influential in the higher-elevation East Troublesome domain, while the rain-snow partitioning temperature threshold dominated in the lower-elevation Chippy Creek and Las Conchas domains. Following wildfire, streamflow sensitivity shifted away from snowpack parameters toward soil hydraulic parameters, consistent with wildfire-driven reductions in evapotranspiration and increases in soil moisture. Incorporating wildfire effects led to complex changes in simulated snowpack accumulation and streamflow timing, generally leading to earlier snowmelt freshet initiation and cessation and reduced late season flows. Simulations were more sensitive to interannual climate variations than wildfire land cover changes. These results highlight both the potential and challenges of ensemble-based inverse modeling for understanding post-wildfire hydrologic change in snow-dominated watersheds.

Colorado, Montana, New Mexico↗

Reservoir releases and land cover interact to drive event-scale nitrate export in a large agricultural basin

Understanding the drivers of nitrate export in rivers is critical for developing effective nutrient management strategies. However, few studies have explored event-scale drivers of export in large river basins with human modifications like reservoirs. Here, we analyzed nitrate concentration-discharge (C-Q) relationships from 215 events at the outlet of the Kansas River Basin, USA (155,690 km 2 ) from 2014 to 2022 to (1) characterize event-scale nitrate export behaviors in a large agricultural basin, and (2) determine how different event characteristics are linked to these behaviors. We found that C-Q behaviors varied greatly, with 60% of events producing nitrate enrichment (n = 130) and 40% of events producing nitrate dilution (n = 85). These behaviors were correlated with complex spatial interlinkages between climate and land cover: across most of the basin, nitrate enrichment was correlated with drier antecedent conditions, but in wetter areas with higher proportions of urban and forested land cover, enrichment was more strongly correlated with precipitation magnitude/intensity. This difference in hydroclimatic controls on nitrate export might be related to differential distributions of nitrate sources within these land covers. Upstream of major reservoirs, however, neither variable was strongly correlated with C-Q behavior, suggesting that nitrate attenuation within reservoirs decouples event-scale concentration signals in upstream waters from those downstream. Reservoir outflows had variable impacts on C-Q behavior, reflecting reservoir-specific variations in nitrate attenuation efficiency. Together, these results identify specific complex interactions between hydroclimate, land cover, reservoir positioning, and individual reservoir properties that control event-scale nitrate export from large basins.

Colorado, Kansas, Nebraska↗

Deep learning error post-processing improves stochastic watershed modeling

Hydrologic extremes, including floods and droughts, pose substantial societal risks that are expected to intensify with climate change. Deterministic watershed models (DWMs) remain a mainstay for modeling these extremes, but lack explicit representation of uncertainty, limiting their utility for risk-informed planning. Stochastic watershed models (SWMs) address this limitation by generating ensembles of streamflow via models of observed DWM residuals. However, most SWMs struggle with the complex dependence between DWM residuals and the underlying hydrologic state, which can complicate stochastic simulations under nonstationary climates. Deep learning (DL) models, whether used as standalone models or post-processors for process-based DWMs, offer a pathway to address this challenge by reducing conditional dependence. In this study, we evaluate SWMs applied to seven models: three process-based models (PRMS, Hymod, and HBV), their hybrid process-DL counterparts, and a pure DL DWM, focusing on daily simulations and extremes under both historical conditions and synthetic climate change scenarios. Results for a case study watershed in Massachusetts show that SWMs applied to hybrid or pure DL DWMs consistently outperform those applied to process-based DWMs. However, an SWM applied to the pure DL model exhibits weaknesses at low flows for this study basin, underscoring the value of hybrid approaches. Extending the analysis across 73 additional basins demonstrates that these improvements are robust and generalizable statewide. This work highlights the potential of a DL-enhanced stochastic watershed modeling framework to advance hydrologic risk prediction under changing climate conditions, offering a scalable methodology for integrating uncertainty into watershed modeling for long-term planning.

Journal of Hydrology↗