Search USGSSearch

USGS · 70019549

Use of chemical and isotopic tracers to characterize the interactions between ground water and surface water in mantled karst

Abstract

In the mantled karst terrane of northern Florida, the water quality of the Upper Floridan aquifer is influenced by the degree of connectivity between the aquifer and the surface. Chemical and isotopic analyses [ 18 O/ 16 O (δ 18 O), 2 H/ 1 H (δD), 13 C/ 12 C (δ 13 C), tritium ( 3 H), and strontium‐87/strontium‐86 ( 87 Sr/ 86 Sr)] along with geochemical mass‐balance modeling were used to identify the dominant hydrochemical processes that control the composition of ground water as it evolves downgradient in two systems. In one system, surface water enters the Upper Floridan aquifer through a sinkhole located in the Northern Highlands physiographic unit. In the other system, surface water enters the aquifer through a sinkhole lake (Lake Bradford) in the Woodville Karst Plain. Differences in the composition of water isotopes (δ 18 O and <δD) in rainfall, ground water, and surface water were used to develop mixing models of surface water (leakage of water to the Upper Floridan aquifer from a sinkhole lake and a sinkhole) and ground water. Using mass‐balance calculations, based on differences in δ 18 O and δD, the proportion of lake water that mixed with meteoric water ranged from 7 to 86% in water from wells located in close proximity to Lake Bradford. In deeper parts of the Upper Floridan aquifer, water enriched in 18 O and D from five of 12 sampled municipal wells indicated that recharge from a sinkhole (1 to 24%) and surface water with an evaporated isotopic signature (2 to 32%) was mixing with ground water. The solute isotopes, δ 13 C and 87 Sr/ 86 Sr, were used to test the sensitivity of binary and ternary mixing models, and to estimate the amount of mass transfer of carbon and other dissolved species in geochemical reactions. In ground water downgradient from Lake Bradford, the dominant processes controlling carbon cycling in ground water were dissolution of carbonate minerals, aerobic degradation of organic matter, and hydrolysis of silicate minerals. In the deeper parts of the Upper Floridan aquifer, the major processes controlling the concentrations of major dissolved species included dissolution of calcite and dolomite, and degradation of organic matter under oxic conditions. The Upper Floridan aquifer is highly susceptible to contamination from activities at the land surface in the Tallahassee area. The presence of post‐ 1950s concentrations of 3 H in ground water from depths greater than 100 m below land surface indicates that water throughout much of the Upper Floridan aquifer has been recharged during the last 40 years. Even though mixing is likely between ground water and surface water in many parts of the study area, the Upper Floridan aquifer produces good quality water, which due to dilution effects shows little if any impact from trace elements or nutrients that are present in surface waters.

Explore related subjects

Keep this discovery

Explore connections, maps & timelines

BibTeXRIS

B. G. Katz, T.B. Coplen, T.D. Bullen, J. Hal Davis. 2005-08-04. Use of chemical and isotopic tracers to characterize the interactions between ground water and surface water in mantled karst. https://doi.org/10.1111/j.1745-6584.1997.tb00174.x

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

KEEP EXPLORING

Related USGS reports

Evaluating micrometeorological estimates of groundwater discharge from Great Basin desert playas

Groundwater availability studies in the arid southwestern United States traditionally have assumed that groundwater discharge by evapotranspiration (ETg) from desert playas is a significant component of the groundwater budget. However, desert playa ETg rates are poorly constrained by Bowen Ratio energy budget (BREB) and eddy-covariance (EC) micrometeorological measurement approaches. Best attempts by previous studies to constrain ETg from desert playas have resulted in ETg rates that are within the measurement error of micrometeorological approaches. This study uses numerical models to further constrain desert playa ETg rates that are within the measurement error of BREB and EC approaches, and to evaluate the effect of hydraulic properties and salinity-based groundwater-density contrasts on desert playa ETg rates. Numerical models simulated ETg rates from desert playas in Death Valley, California and Dixie Valley, Nevada. Results indicate that actual ETg rates from desert playas are significantly below the uncertainty thresholds of BREB- and EC-based micrometeorological measurements. Discharge from desert playas likely contributes less than 2 percent of total groundwater discharge from Dixie and Death Valleys, which suggests discharge from desert playas also is negligible in other basins. Simulation results also show that ETg from desert playas primarily is limited by differences in hydraulic properties between alluvial fan and playa sediments and, to a lesser extent, by salinity-based groundwater density contrasts.

Ground Water

A semi-structured MODFLOW-USG model to evaluate local water sources to wells for decision support

In order to better represent the configuration of the stream network and simulate local groundwater-surface water interactions, a version of MODFLOW with refined spacing in the topmost layer was applied to a Lake Michigan Basin (LMB) regional groundwater-flow model developed by the U.S. Geological. Regional MODFLOW models commonly use coarse grids over large areas; this coarse spacing precludes model application to local management issues (e.g., surface-water depletion by wells) without recourse to labor-intensive inset models. Implementation of an unstructured formulation within the MODFLOW framework (MODFLOW-USG) allows application of regional models to address local problems. A &ldquo;semi-structured&rdquo; approach (uniform lateral spacing within layers, different lateral spacing among layers) was tested using the LMB regional model. The parent 20-layer model with uniform 5000-foot (1524-m) lateral spacing was converted to 4 layers with 500-foot (152-m) spacing in the top glacial (Quaternary) layer, where surface water features are located, overlying coarser resolution layers representing deeper deposits. This semi-structured version of the LMB model reproduces regional flow conditions, whereas the finer resolution in the top layer improves the accuracy of the simulated response of surface water to shallow wells. One application of the semi-structured LMB model is to provide statistical measures of the correlation between modeled inputs and the simulated amount of water that wells derive from local surface water. The relations identified in this paper serve as the basis for metamodels to predict (with uncertainty) surface-water depletion in response to shallow pumping within and potentially beyond the modeled area, see Fienen et al. (2015a).

Illinois, Indiana, Michigan, Wisconsin

Extending the MODPATH algorithm to rectangular unstructured grids

The recent release of MODFLOW-USG, which allows model grids to have irregular, unstructured connections, requires a modification of the particle-tracking algorithm used by MODPATH. This paper describes a modification of the semi-analytical particle-tracking algorithm used by MODPATH that allows it to be extended to rectangular-based unstructured grids by dividing grid cells with multi-cell face connections into sub-cells. The new method will be incorporated in the next version of MODPATH which is currently under development.

Ground Water