Probabilistic sinkhole modelling for hazard assessment

A method for quantitatively assessing sinkhole susceptibility (spatial probability) and hazard (spatio‐temporal probability) has been developed and independently tested in a 50 km² sector of the Ebro Valley evaporite karst. Three genetic types of sinkholes have been mapped in the floodplain and a terrace surface: 947 small cover‐collapse sinkholes (type 1, terrace), large collapse sinkholes (type 2, floodplain) and large subsidence depressions (type 3, floodplain). The type 1 sinkhole inventory includes two temporal populations: 447 sinkholes formed before 24 November 2005, and 500 between that date and 2 November 2006. Sinkhole susceptibility models have been elaborated analysing the statistical relationships between the sinkholes of the 2005 inventory and a set of potential conditioning factors. The independent evaluation (validation) of the susceptibility models by means of several strategies (random, sequentially excluded, and temporal) has allowed us to select the most significant variables for each sinkhole type and assess quantitatively the quality of models; which are reasonable for the three sinkhole types. Validation has also provided information on the contribution of specific variables and the effect of changing their accuracy to the prediction capability of models. Susceptibility models for type 3 sinkholes have been validated satisfactorily with the 2006 sinkhole inventory (temporal validation). The best susceptibility model has been transformed into a hazard map considering the frequency of sinkholes that occurred in each susceptibility class between 2005 and 2006, as well as their average size. The susceptibility and hazard models obtained could be used as an objective basis for the application of mitigation measures, either of preventive or corrective nature. Copyright © 2008 John Wiley & Sons, Ltd.     

Deep seated fault-related volcanogenic H2S as the key agent of high sinkhole concentration areas

Karst collapse sinkholes (KCS) are a peculiar karst morphology resulting from the collapse of the topographic surface caused by subsurface karstification. In the Southern Apennines these phenomena are not randomly distributed but concentrate in several zones, named High Sinkhole Concentration Areas (HSCA), showing peculiar geological, structural and hydrogeological conditions. Gas vents and mineral springs occur widely in these HSCA, and are often stricty related to the KCS. Starting from four representative areas, the aim of this study is to ascertain if there is a genetic link between peculiar mineralogical vs geochemical features of densely fractured/altered rocky masses in KCS, gas vents/springs occurrences and diffuse hypogenic karstification. By means of a multidisciplinary approach (geological/geomorphological, mineralogical and geochemical) we selected and analysed four different kinds of samples related to KCS, gas vent and spring occurrences: (i) altered to deeply altered limestones sampled in the sinkholes; (ii) unaltered limestones sampled in close proximity to sinkhole areas; (iii) gypsum-rich crusts/patinae, precipitated near to gas vents; (iv) gypsum- and calcite-rich precipitates occurring at springs or inside caves. Among neoformed non-carbonate minerals, gypsum is virtually ubiquitous, halides (fluorite, halite, sylvite) also occur in small to trace amounts. The mineralogical assemblages of the different samples show similarities and may also be compatible with hypogenic speleogenesis and with a process of alteration of the carbonate bedrock by means of uprising mineralizing fluids along structural discontinuities. Stable isotopic compositions (S, O) display strong variability in δ³⁴S and δ¹⁸O for sulfate in the different areas, but a deep-seated sulfur source can be hypothesized for many of the studied KCS-related samples. This study has important implications for the relationships between areas of high concentration of sinkholes, regional fault systems, mechanical characteristics of rocks and the high seismicity typical of these areas of the Southern Apennines. © 2018 John Wiley & Sons, Ltd.     

Geoelectric structure of the Ein Gedi sinkhole occurrence site at the Dead Sea shore in Israel

Hazardous sinkholes started to appear in alluvial fans and unconsolidated sediments along the western Dead Sea coast in 1990. Since then hundreds of sinkholes have appeared from north to south along the shoreline. The Electrical Resistivity Tomography (ERT) method was used to achieve a better understanding of the subsurface geoelectric structure at the sinkhole development sites, taking into account that electric parameters (such as resistivity or conductivity) are very sensitive to formation properties and their variations in time. Fifteen image lines were surveyed at the Ein Gedi area during a period of active sinkhole development (in 2001–2002) over an area of 300 × 550 m². Resistivity cross-sections and maps were constructed from 2-D linear surveys. The process of sinkhole formation in the surveyed area is located in a strip 50–70 m wide and 300–500 m long, extending approximately in a north–south direction. The sinkholes are arranged along a tortuous line within this strip. On resistivity maps and sections this U-shaped zone appears as an alternation of high resistivity anomalies of 350–1000 Ωm (at sinkhole group locations) with narrow background resistivity zones of 50–100 Ωm. The large size of resistivity anomalies (250 × 300 m²), which are considerably greater than those of the sinkholes, form one of the features of the sinkhole sites in the Ein Gedi area. The anomalies continue down to the water table or even deeper (maximum of 25–35 m depth). A low resistivity layer of 1–8 Ωm underlies them. The combined analysis of the image results and other geophysical data shows that high resistivity anomalies are associated with the decompaction of the soil mass at the sinkhole development sites and surrounding areas. Recent studies have shown that sinkholes in the Ein Gedi area are developing along the salt western edge located at a depth of 50 m. The subsurface high resistivity anomaly conforms to the sinkhole line (and salt boundary). They are presumably located above the great dissolution caverns at the salt edge. The heterogeneity of the resistivity structure within the high resistivity anomaly (seen in both lateral and vertical planes) confirms that a disintegration of internal formation structure takes place. Away from the sinkhole sites the subsurface resistivity distribution is homogeneous.     

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.     

Evaporite dissolution and development of karst features on the rolling plains of the texas panhandle

More than 400 dolines, consisting of collapse sinkholes and closed depressions, occur in Hall and Briscoe Counties, Texas. Of these, at least 36 sinkholes and two depressions formed between 1940 and 1972. Collapse sinks are typically circular and range up to 100 m in diameter and 15 m in depth. Closed depressions are irregularly-shaped, internally-drained depressions up to 2.4 km in length. Surface fractures up to 1.5 km in length are associated with certain closed depressions. Karst features are forming as a result of intrastratal dissolution of Upper Permian evaporites, primarily halite and secondarily gypsum. Salt dissolution and overburden collapse are evident from geophysical log cross-sections and occur at depths ranging from 195 to 275 m. Dissolution is progressing down dip and mean dissolution rates for salt beds in the Hall County area range from 0.29 to 7.73 cm/yr. The valley of the Prairie Dog Town Fork of the Red River is an area of ground-water discharge for flow systems moving eastward from topographically high areas of the Rolling Plains and from the eastern rim of the High Plains. Halite beds in contact with these waters undergo dissolution. Salt springs and brine seeps, some associated with sinkholes and depressions, occur in the discharge area and indicate that dissolution is an ongoing process.     

TEM study of the geoelectrical structure and groundwater salinity of the Nahal Hever sinkhole site, Dead Sea shore, Israel

Since the 1990s a large number of sinkholes have appeared in the Dead Sea (DS) coastal area. Sinkhole development was triggered by the lowering of the DS level. In the literature the relationship between the sinkholes and the DS level is explained by intrusion of relatively fresh water into the aquifer thereby dramatically accelerating the salt dissolution with creation of subsurface caverns, which in turn cause sinkholes. The main goal of our project was detection and localization of relatively fresh groundwater. During our study we used the transient electromagnetic method (TEM) to measure the electrical resistivity of the subsurface. As a test site we selected Nahal Hever South which is typical for the DS coast. Our results show that resistivity of the shallow subsurface reflects its vertical and lateral structure, e.g., its main hydrogeological elements explain the inter-relations between geology, hydrogeology, and sinkholes. The TEM method has allowed detailed differentiation of layers (clay, salt, etc.) in the subsurface based on their bulk resistivity. The 10 m-thick salt layer composed of idiomorphic crystals of halite deposited during the earlier Holocene period was extrapolated from borehole HS-2 through the study area. It was found that in Nahal Hever the typical value of the bulk resistivity of clay saturated with the DS brine varies between 0.2 and 0.3 Ωm, whereas saturated gravel and sandy sediments are characterized by resistivity between 0.4 and 0.6 Ωm. The high water salinity of the aquifer (enveloping the salt layer) expressed in terms of resistivity is also an important characterization of the sinkhole development mechanism. The electrical resistivity of the aquifer in the vicinity of the salt unit and its western border did not exceed 1 Ωm (in most cases aquifer resistivity was 0.2-0.6 Ωm) proving that, in accordance with existing criteria, the pores of the alluvial sediments are filled with highly mineralized DS brine. However, we suggest that the criterion of the aquifer resistivity responsible for the salt dissolution should be decreased from 1 Ωm to 0.6 Ωm corresponding to the chloride concentration of approximately 100 g/l (the chloride saturation condition reaches 224 g/l in the northern DS basin and 280 g/l in the southern one).Based on TEM results we can reliably conclude that in 2005, when most of sinkholes had appeared at the surface, salt was located within a very low resistivity environment inside sediments saturated with DS brine. Intrusion of relatively fresh groundwater into the aquifer through the 600 × 600 m²area affected by sinkholes has not been observed.     

Geological controls on the formation of alluvial meanders and floodplain wetlands: the example of the Klip River,eastern Free State,South Africa

Floodplain wetlands are common features of rivers in southern Africa, but they have been little studied from a geological or geomorphological perspective. Study of the upper Klip River, eastern Free State, South Africa, indicates strong geological controls on the formation of alluvial meanders and associated floodplain wetlands. Along this river, pronounced and abrupt changes in valley width are strongly linked to lithological variations. Where weakly cemented sandstone crops out, the Klip has laterally eroded bedrock and carved valleys up to 1500 m wide. In these valleys, the river meanders (sinuosity up to ～1·75) on moderate gradients (<0·001) within extensive floodplains marked by numerous oxbow lakes, backswamps and abandoned channels, many of which host substantial wetlands. In contrast, where highly resistant dolerite crops out, lateral erosion of bedrock is restricted, with the Klip tending instead to erode vertically along joints or fractures. Here, valleys are narrower (<200 m), channel‐bed gradients are steeper (>0·003), the river follows a much straighter course (sinuosity ～1·10–1·34), and floodplains are restricted in width. Long‐term landscape development in the Klip and numerous similar catchments depends on the interaction between fluvial processes in the sandstone and dolerite valleys. In the sandstone valleys, vertical erosion rates are controlled by erosion rates of the more resistant dolerites downstream. Hence, in the short‐ to medium‐term (decades to tens of thousands of years), lateral erosion dominates over vertical erosion, with the river concomitantly planing sandstone in the channel floor and reworking floodplain sediments. The thickness of alluvial fill in the sandstone valleys is limited (<4 m), but the resultant meanders are naturally dynamic, with processes such as point bar deposition, cutoff formation and channel avulsion resulting in an assemblage of fluvial landforms. In the longer term (greater than tens of thousands of years), however, vertical erosion will occur in the sandstone valleys as the downstream dolerites are lowered by erosion, resulting in channel incision, floodplain abandonment, and desiccation of the wetlands. Identification of the geological controls on meander and wetland formation provides information vital for the design of effective management guidelines for these ecologically rich habitats, and also contributes to a better understanding of rivers that are intermediate between fully alluvial and fully bedrock. Copyright © 2002 John Wiley & Sons, Ltd.     

Geomorphology,hydrology, and aquatic vegetation drive seasonal hyporheic flow patterns across a gravel‐dominated floodplain

Across 1·7 km² of the Umatilla River floodplain (Oregon, USA), we investigated the influences of an ephemeral tributary and perennial ‘spring channel’ (fed only by upwelling groundwater) on hyporheic hydrology. We derived maps of winter and summer water‐table elevations from data collected at 46 monitoring wells and 19 stage gauges and used resulting maps to infer groundwater flow direction. Groundwater flow direction varied seasonally across the floodplain and was influenced by main channel stage, flooding, the tributary creek, and the location and direction of hyporheic exchange in the spring channel. Hyporheic exchange in the spring channel was evaluated with a geochemical mixing model, which confirmed patterns of floodplain groundwater movement inferred from water‐table maps and showed that the spring channel was fed predominantly by hyporheic water from the floodplain aquifer (87% during winter, 80% during summer), with its remaining flow supplied by upslope groundwater from the adjacent catchment aquifer. Summertime growth of aquatic macrophytes in the spring channel also influenced patterns of hyporheic exchange and groundwater flow direction in the alluvial aquifer by increasing flow resistance in the spring channel, locally raising surface water stage and adjacent water‐table elevation, and thereby altering the slope of the water‐table in the hyporheic zone. The Umatilla River floodplain is larger than most sites where hyporheic hydrology has been investigated in detail. Yet, our results corroborate other research that has identified off‐channel geomorphic features as important drivers of hyporheic hydrology, including previously published modeling efforts from a similar river and field observations from smaller streams. Copyright © 2007 John Wiley & Sons, Ltd.     

Evaluation of trenching,ground penetrating radar (GPR) and electrical resistivity tomography (ERT) for sinkhole characterization

This paper explores the suitability and advantages of combining the trenching technique with geophysical surveys [ground penetrating radar (GPR) and electrical resistivity tomography (ERT)] for sinkhole characterization in a mantled karst area. The approach is applied to two active sinkholes concealed by anthropogenic deposits and formed by contrasting subsidence mechanisms; collapse and sagging. The ERT section acquired across the collapse sinkhole images the clayey fill of the depression as an obvious low resistivity area, showing the approximate location of the sinkhole edges. Spatially dense GPR surveys provide information on the position of the boundaries of the concealed subsidence structures and their three‐dimensional (3D) internal geometry, revealing the dominant subsidence mechanism. We illustrate the impact of several factors on the quality of the GPR data such as sinkhole size, nominal frequency of the antennas, antenna shielding, and the presence of backfilled excavations and above‐surface objects. Trenches provided detailed information on the subsurface structure of the sinkhole, subsidence magnitude, partitioning of the strain, and the position of the sinkhole edges, especially when they are deep enough and excavated across the central sector and perpendicular to the boundaries. The stratigraphic and structural relationships observed in the trench were then used to infer the spatial evolution of the sinkholes (e.g. enlargement), their kinematic behavior (episodic versus progressive), and to differentiate discrete subsidence events and their associated magnitude. Numerical dates were used to estimate average subsidence rates and the recurrence of subsidence events. Such integrated data sets may be used as an objective basis to forecast the future behavior of potentially damaging sinkholes and to assess the associated hazard and risk. Copyright © 2013 John Wiley & Sons, Ltd.     

Detection of preferential infiltration pathways in sinkholes using joint inversion of self-potential and EM-34 conductivity data

The percolation of water in the ground is responsible for measurable electric potentials called self‐potentials. These potentials are influenced by the distribution of the electrical conductivity of the ground. Because sinkholes are associated both with self‐potential and electrical conductivity anomalies, a joint inversion of EM‐34 conductivity and self‐potential data is proposed as a way of delineating the location of these features. Self‐potential and EM conductivity data were obtained at a test site in Normandy (France) where sinkholes and crypto‐sinkholes are present over a karstic area in a chalk substratum overlain by clay‐with‐flint and loess covers. The presence of sinkholes and crypto‐sinkholes is associated with negative self‐potential anomalies with respect to a reference electrode located outside the area where the sinkholes are clustered. The sinkholes also have a conductivity signature identified by the EM‐34 conductivity data. We used the simulated‐annealing method, which is a global optimization technique, to invert jointly EM‐34 conductivity and self‐potential data. Self‐potential and electrical conductivity provide clear complementary information to determine the interface between the loess and clay‐with‐flint formations. The sinkholes and crypto‐sinkholes are marked by depressions in this interface, focusing the groundwater flow towards the aquifer contained in the chalk substratum.     

Influence of bedrock groundwater on streamflow characteristics in a volcanic catchment

Groundwater movements in volcanic mountains and their effects on streamflow discharge and representative elementary area (REA) have remained largely unclear. We surveyed the discharge and chemical composition of spring and stream water in two catchments: the Hontani river (NR) catchment (6.6 km²) and the Hosotani river (SR) catchment (4.0 km²) at the southern part of Daisen volcano, Japan. Daisen volcano is a young volcano (17 × 10³ years) at an early stage of erosion. Our study indicated that deep groundwater that moved through thick lava and pyroclastic flows and that could not be explained by shallow movements controlled by surface topography contributed dominantly to streamflow at larger catchment areas. At the NR catchment, the deep groundwater contribution clearly increased at a catchment boundary defined by an area of 3.0 km² and an elevation of 800 m. At the SR catchment, the contribution deep groundwater to the stream also increased suddenly at a boundary threshold of 2.0 and 700 m. Beyond these thresholds, the contributions of deep bedrock groundwater remained constant, indicating that the REA is between 2 and 3 km² at the observed area. These results indicate that the hydrological conditions of base flow were controlled mainly by the deep bedrock groundwater that moved through thick lava and pyroclastic flows in the undissected volcanic body of the upper part of the catchment. Our study demonstrates that deep and long groundwater movements via a deep bedrock layer including thick deposits of volcanic materials at the two catchments on Daisen volcano strongly determined streamflow discharge instead of the mixing of small‐scale hydrological conditions. Copyright © 2015 John Wiley & Sons, Ltd.     

Morphology and distribution of dolines on ultramafic rocks from airborne LiDAR data: the case of southern Grande Terre in New Caledonia (SW Pacific)

Dolines are closed geomorphological depressions which are surface manifestations of karstic systems. Usually developed on limestones, they also typify the morphology of the New Caledonian landscape, particularly on the southern massif of the main island (known as Massif du Sud). The specificity of dolines here lies in their development on ultramafic rocks. They are evidences of subsidence, suffosion and collapse phenomena resulting from dissolution weathering of peridotites. However, extensive underground drainage systems are still not yet recognized. Semi‐automatic mapping of dolines is carried out on a 148 km² area of the Massif du Sud from a high accuracy LiDAR digital elevation model. In this area 8601 dolines ranging from 1 m² to 2 km² are identified and morphologically characterized with precision. Most are small, shallow and round‐shaped, yet more complex shapes are locally observed. Size distribution analysis allows the setting of a threshold of 20 000 m² above which surface processes rather than chemical weathering control doline evolution. Doline density analysis reveals high concentrations on flat areas where ferricrete overlies the complete weathering profile, especially in the case of elevated rainy watersheds. Dolines are aligned and elongated along a north 135° ± 5° major fracture direction, which is inherited from the obduction of the Pacific Plate upper mantle in the Late Eocene. Finally, we propose a pioneering morphometric typology of dolines that provides important clues as to pseudokarstic activity. We define collapse, bowl‐shaped and flat bottom dolines. Collapse and bowl‐shaped dolines are assumed to denote active pseudokarst. They may widen and deepen, or eventually be filled by sediments. They are distinguished from flat bottom dolines that are partially to completely filled, which suggests that they are associated with paleo‐pseudokarsts. However the groundwater flow paths associated with the genesis and evolution of dolines must be clarified, thus collapse and bowl‐shaped dolines should be hydrologically monitored. Copyright © 2016 John Wiley & Sons, Ltd.     

Effects of bedrock groundwater discharge on spatial variability of dissolved carbon,nitrogen, and phosphorous concentrations in stream water within a forest headwater catchment

The quantitative evaluation of the effects of bedrock groundwater discharge on spatial variability of stream dissolved organic carbon (DOC), dissolved inorganic nitrogen (DIN) and dissolved inorganic phosphorous (DIP) concentrations has still been insufficient. We examined the relationships between stream DOC, DIN and DIP concentrations and bedrock groundwater contribution to stream water in forest headwater catchments in warm-humid climate zones. We sampled stream water and bedrock springs at multiple points in September and December 2013 in a 5 km² forest headwater catchment in Japan and sampled groundwater in soil layer in small hillslopes. We assumed that stream water consisted of four end members, groundwater in soil layer and three types of bedrock groundwater, and calculated the contributions of each end member to stream water from mineral-derived solute concentrations. DOC, DIN and DIP concentrations in stream water were compared with the calculated bedrock groundwater contribution. The bedrock groundwater contribution had significant negative linear correlation with stream DOC concentration, no significant correlation with stream DIN concentration, and significant positive linear correlation with stream DIP concentration. These results highlighted the importance of bedrock groundwater discharge in establishing stream DOC and DIP concentrations. In addition, stream DOC and DIP concentrations were higher and lower, respectively, than those expected from end member mixing of groundwater in soil layer and bedrock springs. Spatial heterogeneity of DOC and DIP concentrations in groundwater and/or in-stream DOC production and DIP uptake were the probable reasons for these discrepancies. Our results indicate that the relationships between spatial variability of stream DOC, DIN and DIP concentrations and bedrock groundwater contribution are useful for comparing the processes that affect stream DOC, DIN and DIP concentrations among catchments beyond the spatial heterogeneity of hydrological and biogeochemical processes within a catchment.     

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.     

Electrical resistivity imaging and hydrochemical analysis for groundwater investigation in Kuala Langat,Malaysia

ABSTRACT

Integrated two-dimensional electrical resistivity imaging (ERI) and hydrochemical surveys were used to investigate the groundwater alluvial aquifer in Kuala Langat, Malaysia. The study in the Langat basin considered the thickness of the aquifer, the depth of the bedrock, the regions influenced by seawater intrusion, and the monitoring of water levels. The resistivity imaging results show that the upper layer consists of clay, while the second layer is an aquifer whose thickness varies mostly in the range of 10–30 m, and in some cases extends to 40 m. The bedrock depth varies from 30 to 65 m. The chemical analyses were carried out on groundwater samples from nine boreholes collected between 2008 and 2012. The analyses indicate that the total dissolved solids (TDS) exceed 1000 mg L^-1 near the coastal area and are often less than 500 mg L^-1 further inland. The ERI and hydrochemical analyses reveal that groundwater in the study area, especially towards the coast, is a mixture of brackish and fresh waters.

EDITOR D. Koutsoyiannis; ASSOCIATE EDITOR M.D. Fidelibus     

Isotopic and hydrochemistry spatial variation of sulfate for groundwater characterization in karstic aquifers

The Sierra Gorda aquifer is one of the most extensive of southern Spain. The main groundwater discharge is produced at its northern boundary through several high‐flow springs. In this study, stable isotopes of dissolved sulfate (δ³⁴S and δ¹⁸O) and groundwater chemistry were used to determine the origin of the sulfate and to characterize the groundwater flow. We sampled the main springs, as well as other minor outlets related to perched water tables, in order to determine the different sources of SO₄^2? (e.g., dissolution of evaporites and atmospheric deposition). The substantial difference in the amount of dissolved SO₄^2? between the springs located in its northwestern part (≈25 mg/L) and those elsewhere in the northern part (≈60 mg/L) suggests zones with separate groundwater flow systems. A third group of springs, far from the northeastern boundary of the permeable outcrops, shows higher SO₄^2? content than the rest (≈125 mg/L). The isotopic range of sulfate (?0.3‰ to 14.82‰ V‐CTD) points to several sources, including dissolution of Triassic or Miocene evaporites, atmospheric deposition, and decomposition of organic material in the soil. Among these, the dissolution of Triassic gypsum—which overlies the saturated zone as a consequence of the folds and faults that deform the aquifer—is the main source of SO₄^2? (range from 12.79‰ to 14.82‰ V‐CTD). This range is typical for Triassic gypsum. The higher karstification in the western sector, together with important differences in the saturated thickness between the western and eastern sectors, would also be due to the tectonic structure and could explain the difference in SO₄^2? contents in the water. This singular arrangement may cause a higher residence time of groundwater in the eastern sector; thus, a higher contact time with Triassic evaporitic rocks is inferred. Accordingly, the stable isotopes of SO₄^2? are found to be a valuable tool for identifying areas with different flow systems in the saturated zone of karstic aquifers, as well as for evaluating aspects such as the degree of karstification .     

Using δ18O and δ2H to Detect Hydraulic Connection Between a Sinkhole Lake and a First-Magnitude Spring

Oxygen and hydrogen isotopes were used in this study to detect a hydraulic connection between a sinkhole lake and a karst spring. In karst areas, surface water that flows to a lake can drain through sinkholes in the lakebed to the underlying aquifer, and then flows in karst conduits and through aquifer matrix. At the study site located in northwest Florida, USA, Lake Miccosukee immediately drains into two sinkholes. Results from a dye tracing experiment indicate that lake water discharges at Natural Bridge Spring, a first-magnitude spring 32 km downgradient from the lake. By collecting weekly water samples from the lake, the spring, and a groundwater well 10 m away from the lake during the dry period between October 2019 and January 2020, it was found that, when rainfall effects on isotopic signature in spring water are removed, increased isotope ratios of spring water can be explained by mixing of heavy-isotope-enriched lake water into groundwater, indicating hydraulic connection between the lake and the spring. Such a detection of hydraulic connection at the scale of tens of kilometers and for a first-magnitude spring has not been previously reported in the literature. Based on the isotope ratio data, it was estimated that, during the study period, about 8.5% the spring discharge was the lake water that drained into the lake sinkholes.     

A coupled vertically integrated model to describe lateral exchanges between surface and subsurface in large alluvial floodplains with a fully penetrating river

This paper presents a vertically averaged model for studying water and solute exchanges between a large river and its adjacent alluvial aquifer. The hydraulic model couples horizontal 2D Saint Venant equations for river flow and a 2D Dupuit equation for aquifer flow. The dynamic coupling between river and aquifer is provided by continuity of fluxes and water level elevation between the two domains. Equations are solved simultaneously by linking the two hydrological system matrices in a single global matrix in order to ensure the continuity conditions between river and aquifer and to accurately model two‐way coupling between these two domains. The model is applied to a large reach (about 36 km²) of the Garonne River (south‐western France) and its floodplain, including an instrumented site in a meander. Simulated hydraulic heads are compared with experimental measurements on the Garonne River and aquifer in the floodplain. Model verification includes comparisons for one point sampling date (27 piezometers, 30 March 2000) and for hydraulic heads variations measured continuously over 5 months (5 piezometers, 1 January to 1 June 2000). The model accurately reproduces the strong hydraulic connections between the Garonne River and its aquifer, which are confirmed by the simultaneous variation of the water level in the river and in piezometers located near the river bank. The simulations also confirmed that the model is able to reproduce groundwater flow dynamics during flood events. Given these results, the hydraulic model was coupled with a solute‐transport component, based on advection‐dispersion equations, to investigate the theoretical dynamics of a conservative tracer over 5 years throughout the 36 km² reach studied. Meanders were shown to favour exchanges between river and aquifer, and although the tracer was diluted in the river, the contamination moved downstream from the injection plots and affected both river banks. Copyright © 2008 John Wiley & Sons, Ltd.     

The impact of groundwater–surface water interactions on the water balance of a mesoscale lowland river catchment in northeastern Germany

The glacially formed northeastern German lowlands are characterized by extensive floodplains, often interrupted by relatively steep moraine hills. The hydrological cycle of this area is governed by the tight interaction of surface water dynamics and the corresponding directly connected shallow groundwater aquifer. Runoff generation processes, as well as the extent and spatial distribution of the interaction between surface water and groundwater, are controlled by floodplain topography and by surface water dynamics. A modelling approach based on extensive experimental analyses is presented that describes the specific water balance of lowland areas, including the interactions of groundwater and surface water, as well as reflecting the important role of time‐variable shallow groundwater stages for runoff generation in floodplains. In the first part, experimental investigations of floodplain hydrological characteristics lead to a qualitative understanding of the water balance processes and to the development of a conceptual model of the water balance and groundwater dynamics of the study area. Thereby model requirements which allow for an adequate simulation of the floodplain hydrology, considering also interactions between groundwater and surface water have been characterized. Based on these analyses, the Integrated Modelling of Water Balance and Nutrient Dynamics (IWAN) approach has been developed. This consists of coupling the surface runoff generation and soil water routines of the deterministic, spatially distributed hydrological model WASIM‐ETH‐I with the three‐dimensional finite‐difference‐based numerical groundwater model MODFLOW and Processing MODFLOW. The model was applied successfully to a mesoscale subcatchment of the Havel River in northeast Germany. It was calibrated for two small catchments (1·4 and 25 km²), where the importance of the interaction processes between groundwater and surface waters and the sensitivity of several controlling parameters could be quantified. Validation results are satisfying for different years for the entire 198 km² catchment. The model approach was further successfully tested for specific events. The experimental area is a typical example of a floodplain‐dominated landscape. It was demonstrated that the lateral flow processes and the interactions between groundwater and surface water have a major importance for the water balance and periodically superimposed on the vertical runoff generation. Copyright © 2006 John Wiley & Sons, Ltd.