Using Stormwater Hysteresis to Characterize Karst Spring Discharge

Discharge from karst springs contains a mixture of conduit and matrix water, but the variations in groundwater mixing are poorly known. Storm events present an opportunity to try to map flow components because water entering during storms is more dilute and provides a tracer as it mixes with pre‐event water along the flowpath from the recharge area to discharge at a spring. We used hysteresis plots of Mg/Ca ratios in a spring in the Cumberland Valley of Pennsylvania to map conduit (higher Ca) vs. diffuse (higher Mg) sources of recharge. We observed two types of temporal heterogeneity: within a storm event and from storm to storm. The timing of the variation in Mg/Ca suggested sources of mixing waters. An increase in the Mg/Ca ratio at the beginning of some storms while conductivity declined suggested diffuse recharge through the epikarst. The rapid changes in Mg/Ca ratios for low‐intensity events probably occurred as the rainfall waxed and waned and illustrate that a variety of flowpaths are available at this spring because additional flushing of Mg occurred. In contrast, the conductivity hysteresis began with dilute water initially and rotation was similar from storm to storm. Hysteresis plots of the Mg/Ca ratio have the potential of revealing more of the complexity in discharge than conductivity alone. A better understanding of flow components in karst is needed to protect these aquifers as a groundwater resource.     

Upscaling Point Velocity Measurements to Characterize a Glacial Outwash Aquifer

Small‐scale point velocity probe (PVP)‐derived velocities were compared to conventional large‐scale velocity estimates from Darcy calculations and tracer tests, and the possibility of upscaling PVP data to match the other velocity estimates was evaluated. Hydraulic conductivity was estimated from grain‐size data derived from cores, and single‐well response testing or slug tests of onsite wells. Horizontal hydraulic gradients were calculated using 3‐point estimators from all of the wells within an extensive monitoring network, as well as by representing the water table as a single best fit plane through the entire network. Velocities determined from PVP testing were generally consistent in magnitude with those from depth specific data collected from multilevel monitoring locations in the tracer test, and similar in horizontal flow direction to the average hydraulic gradient. However, scaling up velocity estimates based on PVP measurements for comparison with site‐wide Darcy‐based velocities revealed issues that challenge the use of Darcy calculations as a generally applicable standard for comparison. The Darcy calculations were shown to underestimate the groundwater velocities determined both by the PVPs and large‐scale tracer testing, in a depth‐specific sense and as a site‐wide average. Some of this discrepancy is attributable to the selective placement of the PVPs in the aquifer. Nevertheless, this result has important implications for the design of in situ treatment systems. It is concluded that Darcy estimations of velocity should be supplemented with independent assessments for these kinds of applications.     

Multiple Wavelet Coherence to Evaluate Local Multivariate Relationships in a Groundwater System

Groundwater level fluctuations are affected by surface properties due to complex correlations of groundwater-surface water interaction and/or other surface processes, which are usually hard to be accurately quantified. Previous studies have assessed the relationship between groundwater level fluctuations and specific controlling factors. However, few studies have been conducted to explore the impact of the combination of multiple factors on the groundwater system. Hence, this paper tries to explore the localized and scale-specific multivariate relationships between the groundwater level and controlling factors (such as hydrologic and meteorological factors) using bivariate wavelet coherence and multiple wavelet coherence. The groundwater level fluctuations of two wells in areas covered by different plant densities (i.e., the riparian zone of the Colorado River, USA) are analyzed. Main findings include three parts. First, barometric pressure and river stage are the best factors to interpret the groundwater level fluctuations at small scales (<1 day) and large scales (>1 day) at the well of low-density plants stand, respectively. Second, at the well of high-density plants stand, the best predictors to control the groundwater level fluctuations include barometric pressure (<1 day), the combination of barometric pressure and temperature (1-7 days), temperature (7-30 days), and the combination of barometric pressure, temperature, and river stage (>30 days). The best predictor of groundwater head fluctuations depends on the variance of the vegetation coverage and hydrological processes. Third, these results provide a suite of factors to explain the groundwater level variations, which is an important topic in water-resource prediction and management.     

Simulating a Highly Parameterized Groundwater Model for a Tropical Eogenetic Karst Aquifer Using Physically-Based Numerical Modeling and Inverse Method

Water Resources - Formulation of effective water management based on the understanding of the interaction between climate and sub-surface components as well as its causative consequences to...     

Differential Transport of Escherichia coli Isolates Compared to Abiotic Tracers in a Karst Aquifer

Lack of filtration and rapid transport of groundwater and particulate matter make karst aquifers susceptible to bacterial contamination. This study utilized quantitative polymerase chain reaction (qPCR) to examine the transport and attenuation of two nonvirulent isolates of Escherichia coli (E. coli) in relation to traditional groundwater tracers (rhodamine WT dye and 1-µm diameter latex microspheres) in a karst-conduit aquifer in central Kentucky. Bacterial isolates were labeled with stable isotopes (¹⁵N and ¹³C). All tracers were detected more than 6 km downstream from the injection site and demonstrated overlapping breakthrough curves, with differential transport observed between the two bacterial strains. The E. coli isolate containing the kps gene (low attachment) arrived at sampling sites 1.25 to 36 h prior to the bacterial isolate containing the iha gene (high attachment) and was detected in samples collected following storm events in which the iha isolate was not detected. The storage potential of contaminants within karst systems was demonstrated by the remobilization of all tracers during storm events more than 1 month after injection. Bacteria-sized microspheres were more easily remobilized during periods of increased discharge compared to other tracers. The study demonstrated that molecular biology techniques such as qPCR can be utilized as a sensitive analysis of bacterial tracers in karst aquifers and may prove to be a more sensitive analytical technique than stable isotope analysis for field-scale traces.     

Comparison of Recharge Estimation Methods During a Wet Period in a Karst Aquifer

Management of water resources, implying their appropriate protection, calls for a sound evaluation of recharge. Such assessment is very complex in karst aquifers. Most methods are developed for application to detrital aquifers, without taking into account the extraordinary heterogeneity of porosity and permeability of karst systems. It is commonly recommended to estimate recharge using multiple methods; however, differences inherent to the diverse methods make it difficult to clarify the accuracy of each result. In this study, recharge was estimated in a karst aquifer working in a natural regime, in a Mediterranean‐type climate, in the western part of the Sierra de las Nieves (southern Spain). Mediterranean climate regions are characterized by high inter‐annual rainfall variability featuring long dry periods and short intense wet periods, the latter constituting the most important contribution to aquifer water input. This paper aims to identify the methods that provide the most plausible range of recharge rate during wet periods. Six methods were tested: the classical method of Thornthwaite‐Mather, the Visual Balan code, the chloride balance method, and spatially distributed methods such as APLIS, a novel spatiotemporal estimation of recharge, and ZOODRM. The results help determine valid methods for application in the rest of the unit of study and in similar karst aquifers.     

Percolation and Particle Transport in the Unsaturated Zone of a Karst Aquifer

Recharge and contamination of karst aquifers often occur via the unsaturated zone, but the functioning of this zone has not yet been fully understood. Therefore, irrigation and tracer experiments, along with monitoring of rainfall events, were used to examine water percolation and the transport of solutes, particles, and fecal bacteria between the land surface and a water outlet into a shallow cave. Monitored parameters included discharge, electrical conductivity, temperature, organic carbon, turbidity, particle-size distribution (PSD), fecal indicator bacteria, chloride, bromide, and uranine. Percolation following rainfall or irrigation can be subdivided into a lag phase (no response at the outlet), a piston-flow phase (release of epikarst storage water by pressure transfer), and a mixed-flow phase (increasing contribution of freshly infiltrated water), starting between 20 min and a few hours after the start of recharge event. Concerning particle and bacteria transport, results demonstrate that (1) a first turbidity signal occurs during increasing discharge due to remobilization of particles from fractures (pulse-through turbidity); (2) a second turbidity signal is caused by direct particle transfer from the soil (flow-through turbidity), often accompanied by high levels of fecal indicator bacteria, up to 17,000 Escherichia coli /100 mL; and (3) PSD allows differentiation between the two types of turbidity. A relative increase of fine particles (0.9 to 1.5 μm) coincides with microbial contamination. These findings help quantify water storage and percolation in the epikarst and better understand contaminant transport and attenuation. The use of PSD as "early-warning parameter" for microbial contamination in karst water is confirmed.     

Storage and Drainage Characteristics of a Highly Heterogeneous Karst Aquifer in Houzhai Basin

The characteristics of karst aquifers are difficult to be determined due to their heterogeneous physical properties and lack of hydrogeological information. In this case study, we applied two methods for a comparative analysis of storage and drainage characteristics in upstream, midstream, and downstream of Houzhai cave stream basin. In the first method, Minimum Smoothed Method (MSM) is used to determine the proportion of baseflow to the total flow (Baseflow Index, BFI). In the second method, a bicarbonate‐base two‐end member mixing model is used to quantify the slow flow component and fast flow component. For both methods, slow flow and quick flow are quantified at three sampling sites, which provide useful information for the analysis of storage and drainage characteristics. The results from flow separation method and hydrogeochemical analysis show a consistently increasing trend of the proportion of slow flow to total flow from the upstream to downstream which indicates that the voids of highly conductive conduits and well‐connected fissures decrease along the flow paths in the Houzhai cave stream basin in southwest China. The upstream areas have a low proportion of baseflow which indicates a high drainage capacity due to high permeable conduits and well‐connected fissures. The downstream areas, on the contrary, have a high proportion of baseflow which indicates a high storage capacity and slow infiltration due to the predominant presence of matrix and poorly‐connected fissures. These numerical methods provide alternative ways to investigate the storage and drainage characteristics of karst aquifers where direct measurement are not available.     

Estimating Aquifer Sensitivity to Nitrate Contamination Using Geochemical Information

Methods for predicting aquifer sensitivity to contamination typically ignore geochemical factors that affect the occurrence of contaminants such as nitrate. Use of geochemical information offers a simple and accurate method for estimating aquifer sensitivity to nitrate contamination. We developed a classification method in which nitrate-sensitive aquifers have dissolved oxygen concentrations > 1.0 mg/L, Eh values >250 mV, and either reduced iron concentrations < 0.1 mg/L or total iron concentrations < 0.7 mg/L. We tested the method in four Minnesota aquifer systems having different geochemical and hydrologic conditions. A surficial sand aquifer in central Minnesota exhibited geochemical zonation, with a rapid shift from aerobic to anaerobic conditions 5 m below the water table. A fractured bedrock aquifer in east-central Minnesota remained aerobic to depths of 50 m, except in areas where anaerobic ground water discharged upward from an underlying aquifer. A bedrock aquifer in southeast Minnesota exhibited aerobic conditions when overlain by surficial deposits lacking shale, whereas anaerobic conditions occurred under deposits that contained shale. Surficial sand aquifers in northwest Minnesota contained high concentrations of sulfate and were anaerobic throughout their extent. Nitrate-nitrogen was detected at concentrations exceeding 1 mg/L in 135 of 149 samples classified as sensitive. Nitrate was not detected in any of the 109 samples classified as not sensitive. We observed differences between our estimates of sensitivity and existing sensitivity maps, which are based on methods that do not consider aquifer geochemistry. Because dissolved oxygen, reduced iron, and Eh are readily measured in the field, use of geochemistry provides a quick and accurate way of assessing aquifer sensitivity to nitrate contamination.     

Groundwater Quality: Analysis of Its Temporal and Spatial Variability in a Karst Aquifer

This study develops an approach based on hierarchical cluster analysis for investigating the spatial and temporal variation of water quality governing processes. The water quality data used in this study were collected in the karst aquifer of Yucatan, Mexico, the only source of drinking water for a population of nearly two million people. Hierarchical cluster analysis was applied to the quality data of all the sampling periods lumped together. This was motivated by the observation that, if water quality does not vary significantly in time, two samples from the same sampling site will belong to the same cluster. The resulting distribution maps of clusters and box‐plots of the major chemical components reveal the spatial and temporal variability of groundwater quality. Principal component analysis was used to verify the results of cluster analysis and to derive the variables that explained most of the variation of the groundwater quality data. Results of this work increase the knowledge about how precipitation and human contamination impact groundwater quality in Yucatan. Spatial variability of groundwater quality in the study area is caused by: a) seawater intrusion and groundwater rich in sulfates at the west and in the coast, b) water rock interactions and the average annual precipitation at the middle and east zones respectively, and c) human contamination present in two localized zones. Changes in the amount and distribution of precipitation cause temporal variation by diluting groundwater in the aquifer. This approach allows to analyze the variation of groundwater quality controlling processes efficiently and simultaneously.