Diffraction imaging of sub‐vertical fractures and karst with full‐resolution 3D Ground‐Penetrating Radar

Vertical fractures with openings of less than one centimetre and irregular karst cause abundant diffractions in Ground‐Penetrating Radar (GPR) records. GPR data acquired with half‐wavelength trace spacing are uninterpretable as they are dominated by spatially undersampled scattered energy. To evaluate the potential of high‐density 3D GPR diffraction imaging a 200 MHz survey with less than a quarter wavelength grid spacing (0.05 m × 0.1 m) was acquired at a fractured and karstified limestone quarry near the village of Cassis in Southern France. After 3D migration processing, diffraction apices line up in sub‐vertical fracture planes and cluster in locations of karstic dissolution features. The majority of karst is developed at intersections of two or more fractures and is limited in depth by a stratigraphic boundary. Such high‐resolution 3D GPR imaging offers an unprecedented internal view of a complex fractured carbonate reservoir model analogue. As seismic and GPR wave kinematics are similar, improvements in the imaging of steep fractures and irregular voids at the resolution limit can also be expected from high‐density seismic diffraction imaging.     

An analytic benchmark test for karst-aquifer flow

A benchmark test for flow in karstic aquifers is presented in the form of an exact solution of the harmonic variations of water flux and head within a karst conduit that is imbedded within a three-dimensional porous matrix having a free surface. The variations are driven by a prescribed variation of head applied at one end of the conduit. The benchmark consists of expressions for the spring discharge as a function of time and the conduit head and flux as functions of distance along the conduit and time. These expressions contain three dimensionless parameters, permitting development of a wide range of specific benchmark tests. The expressions are particularly simple in the case of an infinitely deep aquifer. This limiting solution should provide the most severe test for two-dimensional models of karst aquifer flow. Another limiting case of interest is that in which the conduit diameter is equal to the water depth. This limiting solution should provide the easiest test for two-dimensional models.     

From rainfall to spring discharge: Coupling conduit flow,subsurface matrix flow and surface flow in karst systems using a discrete–continuum model

Physics-based distributed models for simulating flow in karst systems are generally based on the discrete–continuum approach in which the flow in the three-dimensional fractured limestone matrix continuum is coupled with the flow in discrete one-dimensional conduits. In this study we present a newly designed discrete–continuum model for simulating flow in karst systems. We use a flexible spatial discretization such that complicated conduit networks can be incorporated. Turbulent conduit flow and turbulent surface flow are described by the diffusion wave equation whereas laminar variably saturated flow in the matrix is described by the Richards equation. Transients between free-surface and pressurized conduit flow are handled by changing the capacity term of the conduit flow equation. This new approach has the advantage that the transients in mixed conduit flow regimes can be handled without the Preissmann slot approach. Conduit–matrix coupling is based on the Peaceman’s well-index such that simulated exchange fluxes across the conduit–matrix interface are less sensitive to the spatial discretization. Coupling with the surface flow domain is based on numerical techniques commonly used in surface–subsurface models and storm water drainage models. Robust algorithms are used to simulate the non-linear flow processes in a coupled fashion. The model is verified and illustrated with simulation examples.     

Deep water circulation, residence time, and chemistry in a karst complex

We investigated the hydrochemistry of a complex karst hydrosystem made of two carbonate units along a coastal lagoon. Ground water emerges on the lagoon floor from a submarine spring. In addition, thermal waters circulate through the limestone and mix with karst water near the lagoon shore. A distinction between the water from the two carbonate units is related to marine influences and human activities. In one of the massifs, the data show an incongruent dissolution of dolomite with time. In the other system, a slight contamination by saline fluids from the thermal reservoir has led to high calcium and magnesium concentrations. 36Cl, 14C, and 3H data constrain the residence time of the water, and allow for the distinguishing of four circulation types: (1) shallow surface circulation (primarily above sea level) in the karstic units with short residence times (<20 years); (2) shallow subsurface circulation (approximately 0 to -50 m) below the karstic units with residence time in the order of 50 years; (3) deep circulation at depth of 700 to 1500 m in the Jurassic limestones below thick sedimentary cover, with residence time of several thousand years for a part of the water; and (4) deep circulation at a depth of approximately 2500 m, which represents the thermal reservoir in the Jurassic units with residence time of approximately 100,000 years. An interpretative hydrogeological framework is based on the constraints of the geochemical analyses of the deep thermal system, and by water flow from the surface to the deep parts of the carbonate formations.     

Ground penetrating radar imaging of cap rock, caliche and carbonate strata

Field experiments show ground penetrating radar (GPR) can be used to image shallow carbonate stratigraphy effectively in a variety of settings. In south Florida, the position and structure of cap rock cover on limestone can be an important control on surface water flow and vegetation, but larger scale outcrops (tens of meters) of cap rock are sparse. GPR mapping through south Florida prairie, cypress swamp and hardwood hammock resolves variations in thickness and structure of cap rock to 3 m and holds the potential to test theories for cap rock–vegetation relationships. In other settings, carbonate strata are mapped to test models for the formation of local structural anomalies. A test of GPR imaging capabilities on an arid caliche (calcrete) horizon in southeastern Nevada shows depth penetration to 2 m with resolution of the base of caliche. GPR profiling also succeeds in resolving more deeply buried (5 m) limestone discontinuity surfaces that record subaerial exposure in south Florida.     

Transport of suspended solids from a karstic to an alluvial aquifer: the role of the karst/alluvium interface

This study focuses on the coupled transport of dissolved constituents and particulates, from their infiltration on a karst plateau to their discharge from a karst spring and their arrival at a well in an alluvial plain. Particulate markers were identified and the transport of solids was characterised in situ in porous and karstic media, based on particle size analyses, SEM, and traces. Transport from the sinkhole to the spring appeared to be dominated by flow through karst: particulate transport was apparently conservative between the two sites, and there was little difference in the overall character of the particle size distribution of the particulates infiltrating the sinkhole and of those discharging from the spring. Qualitatively, the mineralogy of the infiltrating and discharging material was similar, although at the spring an autochthonous contribution from the aquifer was noted (chalk particles eroded from the parent rock by weathering). In contrast, transport between the spring and the well appears to be affected by the overlying alluvium: particles in the water from the well, showed evidence of considerable size-sorting. Additionally, SEM images of the well samples showed the presence of particles originating from the overlying alluvial system; these particles were not found in samples from the sinkhole or the spring. The differences between the particulates discharging from the spring and the well indicate that the water pumped from the alluvial plain is coming from the karst aquifer via the very transmissive, complex geologic interface between the underlying chalk formation and the gravel at the base of the overlying alluvial system.     

Challenges of using electrical resistivity method to locate karst conduits—A field case in the Inner Bluegrass Region, Kentucky

Conduits serve as major pathways for groundwater flow in karst aquifers. Locating them from the surface, however, is one of the most challenging tasks in karst research. Geophysical methods are often deployed to help locate voids by mapping variations of physical properties of the subsurface. Conduits can cause significant contrasts of some physical properties that can be detected; other subsurface features such as water-bearing fractures often yield similar contrasts, which are difficult to distinguish from the effects of the conduits. This study used electrical resistivity method to search for an unmapped karst conduit that recharges Royal Spring in the Inner Bluegrass karst region, Kentucky, USA. Three types of resistivity techniques (surface 2D survey, quasi-3D survey, and time-lapse survey) were used to map and characterize resistivity anomalies. Some of the major anomalies were selected as drilling targets to verify the existence of the conduits. Drilling near an anomaly identified by an electrical resistivity profile resulted in successful penetration of a major water-filled conduit. The drilling results also suggest that, in this study area, low resistivity anomalies in general are associated with water-bearing features. However, differences in the anomaly signals between the water-filled conduit and other water-bearing features such as water-filled fracture zones were undistinguishable. The electrical resistivity method is useful in conduit detection by providing potential drilling targets. Knowledge of geology and hydrogeology about the site and professional judgment also played important roles in locating the major conduit.     

Estimating Preferential Flow in Karstic Aquifers Using Statistical Mixed Models

Karst aquifers are highly productive groundwater systems often associated with conduit flow. These systems can be highly vulnerable to contamination, resulting in a high potential for contaminant exposure to humans and ecosystems. This work develops statistical models to spatially characterize flow and transport patterns in karstified limestone and determines the effect of aquifer flow rates on these patterns. A laboratory‐scale Geo‐HydroBed model is used to simulate flow and transport processes in a karstic limestone unit. The model consists of stainless steel tanks containing a karstified limestone block collected from a karst aquifer formation in northern Puerto Rico. Experimental work involves making a series of flow and tracer injections, while monitoring hydraulic and tracer response spatially and temporally. Statistical mixed models (SMMs) are applied to hydraulic data to determine likely pathways of preferential flow in the limestone units. The models indicate a highly heterogeneous system with dominant, flow‐dependent preferential flow regions. Results indicate that regions of preferential flow tend to expand at higher groundwater flow rates, suggesting a greater volume of the system being flushed by flowing water at higher rates. Spatial and temporal distribution of tracer concentrations indicates the presence of conduit‐like and diffuse flow transport in the system, supporting the notion of both combined transport mechanisms in the limestone unit. The temporal response of tracer concentrations at different locations in the model coincide with, and confirms the preferential flow distribution generated with the SMMs used in the study.     

Using 2D inversion of magnetic resonance soundings to locate a water-filled karst conduit

A new methodology for magnetic resonance sounding (MRS) data acquisition and interpretation was developed for locating water-filled karst cavities. This methodology was used to investigate the Ouysse karst system in the Poumeyssens shaft in the Causse de Gramat (France). A new 2D numerical MRS response model was designed for improved accuracy over the previous 1D MRS approach. A special survey performed by cave divers confirmed the accuracy of the MRS results. Field results demonstrated that in favourable conditions (a low EM noise environment and a relatively shallow, large target) the MRS method, used with a coincident transmitter/receiver loop, can be an effective tool for locating a water-filled karst conduit. It was shown numerically that because an a priori orientation of the MRS profile with the karst conduit is used in the inversion scheme (perpendicular for instance), any error in this assumption introduces an additional error in locating the karst. However, the resulting error is within acceptable limits when the deviation is less than 30°. The MRS results were compared with an electrical resistivity tomography (ERT) survey. It was found that in Poumeyssens, ERT is not able to locate the water-filled karst. On the other hand, ERT provides additional information about heterogeneities in the limestone.     

Identifying and characterizing solution conduits in karst aquifers through geospatial (GIS) analysis of porosity from borehole imagery: An example from the Biscayne aquifer,South Florida (USA)

We apply geospatial analysis to borehole imagery in an effort to develop new techniques to evaluate the spatial distribution and internal structure of karst conduits. Remote sensing software is used to classify a high resolution, digital borehole image of limestone bedrock from the Biscayne aquifer (South Florida, USA) into a binary image divided into cells of rock matrix and pores. Within a GIS, 2D porosity is calculated for a series of rectangular sampling windows placed over the binary image and then plotted as a function of depth. Potential conduits that intersect the borehole are identified as peaks of high porosity. A second GIS technique identifies a conduit as a continuous object that spans the entire borehole width. According to these criteria, geospatial analysis reveals ∼10 discrete conduits along the ∼15 m borehole image. Continuous sampling of the geologic medium intersected by the borehole provides insight into the internal structure of karst aquifers and the evolution of karst features. Most importantly, this pilot study demonstrates that GIS-based techniques are capable of quantifying the depths, dimensions, shapes, apertures and connectivity of potential conduits, physical attributes that impact flow in karst aquifers.     

Combining Electrical Resistivity Tomography and Ground Penetrating Radar to study geological structuring of karst Unsaturated Zone

This paper highlights the efficiency and complementarity of a light package of geophysical techniques to study the structure of karst Unsaturated Zone (UZ) in typical Mediterranean environment where soil cover is thin or absent. Both selected techniques, 2D Ground Penetrating Radar (GPR) and Electrical Resistivity Tomography (ERT), are widely used in environmental studies and their application is accessible for a lot of scientists/engineers. However, GPR or ERT alone is not able to provide an enhanced characterization of geological features in karst media. In the present study, GPR results supply a near surface high resolution imaging and thus can provide relevant geological information such as stratifications and fractures. Despite the quality of the results GPR's investigation depth remains limited to around 12 m. Apparent and inverted resistivity provided by ERT surveys shows strong lateral and vertical variations. These variations can inform about general geological structuring and feature orientation. ERT is able to prospect down to 40 m but it's a low resolution integrative technique. In the study area the investigated limestone is a commonly electrical resistive formation (more than 2000 Ω.m). However deeper than 5–7 m, the ERT profiles reveal several zones of moderate resistivity (around 900 Ω.m). In these zones a stratification change corresponding to slanted bedding is clearly identified by GPR results. The combination of both GPR and ERT results can allow a well-established geological interpretation. These moderate resistivity zones with slanted beddings can explain the presence of a perennial water flow point 35 m below the surface of the studied site within the underground gallery of the Low-Noise Underground Laboratory (LSBB).     

A note on the role of the lichen collema auriforma in solution basin development on a carboniferous limestone substrate

A study carried out on Carboniferous limestone in the north and west of Ireland supports the idea that rock substrate is removed by the direct mechanical action of lichens. An experiment in which the lichen Collema auriforma was subjected to a number of wetting-drying cycles, showed, using scanning electron microscopy, that contraction of the lichen thallus during the drying phase plucked rock fragments from the substrate surface. This process could contribute to the formation of karstic features including solution basins.     

Isotope hydrological study of mean transit times in an alpine basin (Wimbachtal, Germany)

Measurements of tritium and ¹⁸O concentrations in precipitation and runoff were used to provide further insight into the groundwater storage properties of the Wimbachtal Valley, a catchment area of 33.4 km², extending between 636 and 2713 m a.s.l. in the Berchtesgaden Alps. The catchment includes three aquifer types: a dominant porous aquifer; a fractured dolomite; a karstic limestone aquifer. Employing a simple hydrological model, information about mean transit times of environmental tracers is derived for the groundwater runoff component and several karst springs from the application of the exponential and dispersion flow models to the isotopic input and output data. The mean transit times calculated from a dispersion model with transit times of 4.1 years for ¹⁸O and 4.2 years for tritium, which agree well, allow calculation of total (mobile + stagnant) groundwater storage volume, which is equivalent to 6.6 m of water depth. Direct runoff appears negligible as in many other cases.     

Spatially varying dispersion to model breakthrough curves

Often the water flowing in a karst conduit is a combination of contaminated water entering at a sinkhole and cleaner water released from the limestone matrix. Transport processes in the conduit are controlled by advection, mixing (dilution and dispersion), and retention-release. In this article, a karst transport model considering advection, spatially varying dispersion, and dilution (from matrix seepage) is developed. Two approximate Green's functions are obtained using transformation of variables, respectively, for the initial-value problem and for the boundary-value problem. A numerical example illustrates that mixing associated with strong spatially varying conduit dispersion can cause strong skewness and long tailing in spring breakthrough curves. Comparison of the predicted breakthrough curve against that measured from a dye-tracing experiment between Ames Sink and Indian Spring, Northwest Florida, shows that the conduit dispersivity can be as large as 400 m. Such a large number is believed to imply strong solute interaction between the conduit and the matrix and/or multiple flow paths in a conduit network. It is concluded that Taylor dispersion is not dominant in transport in a karst conduit, and the complicated retention-release process between mobile- and immobile waters may be described by strong spatially varying conduit dispersion.     

Detection of model voids by identifying reverberation phenomena in GPR records

In order to improve the reliability of the ground penetrating radar (GPR) method in identifying subsurface sinkholes and karst cavities, laboratory investigations have been performed. The main objective of this work was to examine the relationship between horizontal/vertical voids dimensions and wavelengths of various antennas, and the corresponding GPR responses. Emphasis was given to the investigation of the factors that cause the appearance of reverberation phenomena in the signal pattern.The tests were conducted in 5 m × 10 m area by 2-m-deep trench filled with homogenous, dry sand. The voids models (empty fiberglass cylinders in diameters of 0.6 m, 1.0 m, 1.5 m and 2.4 m, with various heights) were buried vertically with their tops at depths of between 0.7 and 1.5 m. Investigations were performed for the various model conditions by towing 500, 300 and 100 MHz antennas along a pre-established grid, for the various model conditions.The GPR data collected using the 500 MHz bistatic antenna above the 1.0-m- and the 1.5-m-diameter cylinders, and using the 300 MHz bistatic antenna above the 1.5-m-diameter cylinder, confirmed the presence of a reverberation phenomenon, i.e. a strong convex signal pattern, containing a series of high amplitude extending oscillations with reduced frequency.Based on past practical GPR experience of void detection and presently obtained experimental data, two rules of thumbs may be adopted for the prediction of the appearance of resonant radar pictures:

1. The void diameter larger than the wavelength in air of the antenna used.

2. The vertical size of the empty void not significantly smaller than its horizontal dimension.

The strong reverberations generated by the inner surface of the void targets were found to approximate standing waves generated in cylindrical waveguides and waveguide resonators. The theoretical, experimental and practical results obtained concur.     

The application of GPR and ERI in combination with exposure logging and retrodeformation analysis to characterize sinkholes and reconstruct their impact on fluvial sedimentation

This work illustrates the practicality of investigating sinkholes integrating data gathered by ground penetrating radar (GPR), electrical resistivity imaging (ERI) and trenching or direct logging of the subsidence‐affected sediments in combination with retrodeformation analysis. This mutidisciplinary approach has been tested in a large paleosinkhole developed during the deposition of a Quaternary terrace on salt‐bearing evaporites. The subsidence structure, exposed in an artificial excavation, is located next to Puilatos, a village that was abandoned in the 1970s due to severe subsidence damage. Detailed logging of the exposure revealed that the subsidence structure corresponds to an asymmetric sagging and collapse paleosinkhole with no clear evidence of recent activity. The sedimentological and structural relationships together with the retrodeformation analysis indicate that synsedimentary subsidence controlled channel location, the development of a palustrine environment and local changes in the channel pattern. GPR profiles were acquired using an array of systems with different antenna frequencies, including some recently developed shielded antennas with improved vertical resolution and penetration depth. Although radargrams imaged the faulted sagging structure and provided valuable data on fault throw, they did not satisfactorily image the complex architecture of the fluvial deposit. ERI showed lower resolution but higher penetration depth when compared to GPR, roughly capturing the subsidence structure and yielding information on the thickness of the high‐resistivity alluvium and the nature of the underlying low‐resistivity karstic residue developed on top of the halite‐bearing evaporitic bedrock. Data comparison allows the assessment of the advantages and limitations of these complementary techniques, highly useful for site‐specific sinkhole risk management. Copyright © 2016 John Wiley & Sons, Ltd.     

Shallow karst exploration using MT-VLF and DC resistivity methods1

A geophysical test was carried out over a well-located and fairly embedded karstic conduit. The MT-VLF method was selected because of its high resolution and its ability to provide a resistivity parameter sensitive to water and Clay. This method was used together with DC resistivity methods which allow vertical adjustment of the VLF data and show consistency between the investigation and target depths. After correcting the deformations due to the polarization of the primary field, the MT-VLF data show clearly, in the Central part of the site, that the conduit does not coincide with an anomaly axis but coincides with the boundary between a conductive area and a resistive area. 2D modelling confirms that direct detection of the conduit is not feasible and that the conduit is located close to a conductive zone corresponding to a completely Clay-filled fractured zone. This situation was observed on the whole site and the conduit seemed systematically to thread its way between the conductive zones to join the outlet. The distribution of the conductive fractured zones and the direction of the hydraulic gradient were two important elements in predicting the location of the conduits. A 3D approach would increase the probability of finding the conduits in such a case.     

Effects of glaciations and permafrost upon the development of karst in Canada

Canadian examples suggest that karst landforms may be divided into eight types in terms of their temporal relationships to the record of repeated Quaternary glaciations. Two types are postglacial, two are subglacial, one type occurs where glacial features are adapted to karstic drainage, and three types display sequences of karstic and glacial action. Glacier effects upon karst landforms and their underlying aquifers display the gamut of possibilities. They may destroy, inhibit, preserve, or stimulate karst development. Where continuous permafrost is maintained when covered by glacier ice, postglacial karst is limited to the active layer epikarst. Where permafrost is thawed beneath ice or during deglaciation there are a variety of postglacial karst developments, depending in part upon climate and in part upon local lithologic and relief conditions.     

Hydrodynamic behaviour of karstic aquifers in Boroujerd,western Iran

Abstract

Considering the geological conditions of the southwest of Boroujerd and northwest of Doroud, Iran, karst development is analysed with respect to the hydrodynamic behaviour of the main draining springs of the units and the karstic aquifers are classified as either those developed in Cretaceous limestone or those developed in Oligomiocene limestone. For this purpose, the yields of the main karstic springs of the region—Absardeh and Zoorabad (Cretaceous karstic limestone aquifer), Kalamsooz and Azizabad (Oligomiocene karstic limestone aquifer)—were measured and analysed. Analysis of the recession curve is used for hydrodynamical analysis and to construct the conceptual model for estimation of karst development in the selected aquifers. Based on the results, the dynamic storage capacity of the saturated zone in Cretaceous limestone is evaluated as low to medium and that in Oligomiocene limestone as medium to high. The dynamic storage capacity of the unsaturated zone in Cretaceous limestone is evaluated as high and that in Oligomiocene limestone as low to medium. Moreover, the contribution of quickflow in karstic aquifers developed in the Cretaceous limestone drained by the Absardeh and Zoorabad springs is 23.5 and 82.2%, respectively, and that for the Kalamsooz and Azizabad springs (Oligomiocene limestone) is 5.7 and 22.5%, respectively. Flow in the Cretaceous limestone aquifer drained by the Zoorabad Spring is of concentrated type and the main flow occurs in the well-developed karstic conduits. The main flow in the Oligomiocene limestone aquifer, drained by the Kalamsooz Spring, occurs in a network of joints and fractures and the contribution of concentrated flow is very low. The transmissivity of the saturated zone in the karstic aquifer drained by the Zoorabad and Absardeh springs is medium to high and that for the Kalamsooz and Azizabad springs is found to be low.