Studying the Flow Dynamics of a Karst Aquifer System with an Equivalent Porous Medium Model

The modeling of groundwater flow in karst aquifers is a challenge due to the extreme heterogeneity of its hydraulic parameters and the duality in their discharge behavior, that is, rapid response of highly conductive karst conduits and delayed drainage of the low‐permeability fractured matrix after recharge events. There are a number of different modeling approaches for the simulation of the karst groundwater dynamics, applicable to different aquifer as well as modeling problem types, ranging from continuum models to double continuum models to discrete and hybrid models. This study presents the application of an equivalent porous model approach (EPM, single continuum model) to construct a steady‐state numerical flow model for an important karst aquifer, that is, the Western Mountain Aquifer Basin (WMAB), shared by Israel and the West‐Bank, using MODFLOW2000. The WMAB was used as a catchment since it is a well‐constrained catchment with well‐defined recharge and discharge components and therefore allows a control on the modeling approach, a very rare opportunity for karst aquifer modeling. The model demonstrates the applicability of equivalent porous medium models for the simulation of karst systems, despite their large contrast in hydraulic conductivities. As long as the simulated saturated volume is large enough to average out the local influence of karst conduits and as long as transport velocities are not an issue, EPM models excellently simulate the observed head distribution. The model serves as a starting basis that will be used as a reference for developing a long‐term dynamic model for the WMAB, starting from the pre‐development period (i.e., 1940s) up to date.     

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.     

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.     

Unsaturated Zone Flow Processes and Aquifer Response Time in the Chalk Aquifer,Brighton, South East England

The Chalk aquifer is one of the main sources of water in South East England. The unsaturated zone in the aquifer plays an important role controlling the time and magnitude of recharge and is major pathway for contaminant transport to the water table. A range of previous work has addressed flow processes in the Chalk unsaturated zone, but physical understanding is still incomplete. Here we present the results of a study on flow mechanism in the Chalk unsaturated zone using a combination of statistical analysis and novel laboratory methods. The study was undertaken at three sites (North Heath Barn [NHB], Pyecombe East [PE], and Preston Park [PP]) on the Chalk of the Brighton block, South East England. Daily and hourly time series data of groundwater level and rainfall were correlated. The results show that a slower groundwater level response to rainfall occurs during dry seasons (summer and autumn) when the amount of effective rainfall is less than 4 mm/day, with a thicker and drier unsaturated zone. A faster response occurs during wet seasons (winter and spring) when the daily effective rainfall exceeds 4 mm/day with a thinner and wetter unsaturated zone. Periods of very rapid response (within 15 h) were observed during wet seasons at NHB and PE sites, with unsaturated hydraulic conductivity (K_u) inferred to reach 839 mm/day. A slower response was observed at an urbanized site (PP) as a result of reduction in direct recharge due to reduced infiltration, due to presences of impermeable infrastructure covering the area around PP borehole. Laboratory measurements of K_u of the Chalk matrix using a geotechnical centrifuge show variation from 4.27 to 0.07 mm/day, according to the level of saturation. Thus, the rapid responses cannot be linked to matrix flow only but indicate the contribution of fracture and karstic flow processes in conducting water.     

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.     

Field Testing the Hypothesis of Darcian Flow Through a Carbonate Aquifer

The acceptability of the hypothesis of Darcian flow through a semiconfined carbonate aquifer was tested prior to running a multiple-day aquifer test in Pinellas County, Florida. The approach used to test the hypothesis was to run a number of hour-long aquifer tests at different discharges with drawdown measured at the same time during each test in two observation wells, one at 35 feet and the other at 733 feet from the pumped well. If the hypothesis were acceptable, a plot of drawdown versus discharge should describe a straight line. A linear relation was described in 15-minute and 1-hour data plots of drawdown versus discharge for both observation wells. Thus, the hypothesis of Darcian flow through the semiconfined carbonate aquifer was deemed acceptable. Drawdown data from the distant observation well collected during the following multiple-day aquifer test were then analyzed for aquifer-property values of transmissivity, storage coefficient, and leakance coefficient using standard methods. Discharge for the multiple-day test was within the discharge range of the hour-long tests.     

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.     

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.     

Characterization of Flow Dynamics and Vulnerability in a Coastal Aquifer System

Traditional aquifer vulnerability techniques primarily rely on spatial property data for a region and are limited by their ability to directly or indirectly assess flow and transport processes occurring from the surface to depth within an aquifer system. The main objective of this study was to investigate groundwater vulnerability in terms of aquifer interconnectivity and flow dynamics. A combination of stable isotopes, groundwater age‐dating (radiocarbon), and geomorphic/geogenic spatial analyses was applied to a regional, highly developed coastal aquifer to explain the presence of nitrate at depth. The average δ¹³C value (?17.3 ± 2‰ VPDB, n = 27) is characteristic of groundwater originating from locally infiltrated precipitation through extensively cultivated soils. The average δ¹⁸O and δD values (?4.0 ± 0.1‰ VSMOW, n = 27; δD: ?19.3 ± 1‰ VSMOW, n = 27, respectively) are similar to precipitation water derived from maritime sources feeding the region's surface water and groundwater. Stable and radioactive isotopes reveal significant mixing between shallow and deep aquifers due to high velocities, hydraulic connection, and input of local recharge water to depths. Groundwater overdevelopment enhances deeper and faster modern water downward flux, amplifying aquifer vulnerability. Therefore, aquifer vulnerability is a variable, dependent on the type and degree of stress conditions experienced by a groundwater system as well as the geospatial properties at the near surface.     

Nonlinear Flow in Karst Formations

The variation of effective hydraulic conductivity as a function of specific discharge in several 0.2-m and 0.3-m cubes of Key Largo Limestone was investigated. The experimental results closely match the Forchheimer equation. Defining the pore-size length scale in terms of Forchheimer parameters, it is demonstrated that significant deviations from Darcian flow will occur when the Reynolds number exceeds 0.11. A particular threshold model previously proposed for use in karstic formations does not show strong agreement with the data near the onset of nonlinear flow.     

Using Wavelet Coherence to Characterize Surface Water Infiltration into a Low-Lying Karst Aquifer

Karstified carbonate aquifers may receive significant recharge contributions from losing streams, hence, the knowledge about surface water-groundwater (SW-GW) interactions is crucial with regard to water management (e.g., source protection zone delineation). The dynamics of SW-GW interactions may depend on factors such as the relative water levels between streams and aquifers, resulting in a temporal variation of exchange, which imposes complexity to the understanding of such dynamics. This study highlights the use of high-resolution time series and multiresolution analysis to help to gain insights into such complex dynamics. Wavelet coherence is applied on hourly time series of rainfall, stream, and spring discharges of a low-lying karstified spring catchment to yield a correlation in the time-frequency domain. This analysis provides comprehensive information on the overall impact of the river on the spring, which is supported by the cross-correlation function, as well as by more detailed information, including time-variant influences such as a threshold level of influence. Field observations of turbidity sampling at the spring appear to support this interpretation. This innovative approach relies on basic hydrological parameters, water level, or discharge, and is therefore applicable to many other systems with such existing time series.