PaPRIKa: a method for estimating karst resource and source vulnerability��application to the Ouysse karst system (southwest France)

The intrinsic vulnerability mapping method, PaPRIKa, is proposed as a common basis for karst groundwater protection in France. PaPRIKa is a specialized method for studying karst aquifers, derived from updating the RISKE and EPIK methods. Both the structure and functioning of karst aquifers are considered in order to develop a resource and source-vulnerability mapping method. PaPRIKa means Protection of aquifers from the assessment of four criteria: P for protection (considering the most protective aspects among parameters related to soil cover, unsaturated zone and epikarst behavior), R for rock type, I for infiltration and Ka for karstification degree. The Ouysse karst system, located in the Causses area in southwest France, is one of the nine pilot sites where this method was tested and standardized. The specificities of the Ouysse system such as the size of the catchment area, the spatial variability of the karst network development, the thick infiltration zone and the system??s dual character (both karst and non-karst areas), have provided a valuable field of application. The vulnerability of the resource was assessed for the entire catchment area, while source-orientated cartography was attempted for the catchment areas of the three different capture works used for drinking water.     

Regional-scale analysis of karst underground flow deduced from tracing experiments: examples from carbonate aquifers in Malaga province,southern Spain

Tracer concentration data from field experiments conducted in several carbonate aquifers (Malaga province, southern Spain) were analyzed following a dual approach based on the graphical evaluation method (GEM) and solute transport modeling to decipher flow mechanisms in karst systems at regional scale. The results show that conduit system geometry and flow conditions are the principal factors influencing tracer migration through the examined karst flow routes. Solute transport is mainly controlled by longitudinal advection and dispersion throughout the conduit length, but also by flow partitioning between mobile and immobile fluid phases, while the matrix diffusion process appears to be less relevant. The simulation of tracer breakthrough curves (BTCs) suggests that diffuse and concentrated flow through the unsaturated zone can have equivalent transport properties under extreme recharge, with high flow velocities and efficient mixing due to the high hydraulic gradients generated. Tracer mobilization within the saturated zone under low flow conditions mainly depends on the hydrodynamics (rather than on the karst conduit development), which promote a lower longitudinal advection and retardation in the tracer migration, resulting in a marked tailing effect of BTCs. The analytical advection-dispersion equation better approximates the effective flow velocity and longitudinal dispersion estimations provided by the GEM, while the non-equilibrium transport model achieves a better adjustment of most asymmetric and long-tailed BTCs. The assessment of karst underground flow properties from tracing tests at regional scale can aid design of groundwater management and protection strategies, particularly in large hydrogeological systems (i.e. transboundary carbonate aquifers) and/or in poorly investigated ones.     

Carbone organique total (COT) et magnésium (Mg2+) : deux traceurs complémentaires du temps de séjour dans l'aquifère karstique

The Total Organic Carbon (TOC) is an interesting tracer of fast infiltration within karstic systems [3,7]. Regular sampling on several aquifers, from the experimental site of Vaucluse, made it possible to demonstrate the high sensitivity of this tracer compared with other commonly used chemical and isotopic tracers in karstic hydrogeology. The complementarity of magnesium, indicator of the residence time of water within the system, and TOC appears as a relevant tool in order to characterize the behaviour of the aquifer, to differentiate the water types which participate to the karstic flow (fast infiltration, unsaturated zone, saturated zone) and then, to evaluate their vulnerability.     

Development of a coupled surface-groundwater-pipe network model for the sustainable management of karstic groundwater

This paper considers the hydrogeological simulation of groundwater movement in karstic regions using a hydrological modelling system (SHETRAN) which has been adapted for modelling flow in karstic aquifers. Flow and transport through karstic aquifers remains poorly understood, yet quantitative hydrogeological models are essential for developing and implementing groundwater protection policies. The new model has been developed and used within the STALAGMITE (Sustainable Management of Groundwater in Karstic Environments) project, funded by the European Commission. The SHETRAN model is physically based insofar as most of the parameters have some physical meaning. The SHETRAN model represents all of the key processes in the hydrological cycle, including subsurface flow in the saturated and unsaturated zones, surface flow over the ground surface and in channels, rainfall interception by vegetation canopies, evapotranspiration, snow-pack development and snowmelt. The modifications made to SHETRAN to simulate karstic aquifers are (1) the coupling of a pipe network model to a variably saturated, three-dimensional groundwater component (the VSS-NET component), to simulate flow under pressure in saturated conduits; (2) the coupling of surface water features (e.g. sinking streams or "ponors", and spring discharges) to the conduit system; (3) the addition of a preferential "bypass" flow mechanism to represent vertical infiltration through a high-conductivity epikarst zone. Lastly, a forward particle tracking routine has been developed to trace the path of hypothetical particles with matrix and pipe flow to springs or other discharge points. This component allows the definition of groundwater protection zones around a source for areas of the catchment (watershed) which are vulnerable to pollution from non-point sources (agriculture and forestry).     

Temporal variation of chemical and isotopic signals in major discharges of an alpine karst aquifer in Turkey: implications with respect to response of karst aquifers to recharge

Proper management of karst aquifers requires a better understanding of flow and transport mechanisms in these systems. Flow in karst aquifers is inherently very complex due to the non-linear and non-stationary relationship between recharge and discharge. Information on this relationship has been acquired for a large (1,000 km²), mountainous (>3,500 m asl) karst aquifer with a deep unsaturated zone (>2,000 m) in the Aladaglar mountain range of south-central Turkey. All major discharges from the aquifer, which drain almost all the recharge, have been observed periodically for specific electrical conductivity, tritium and oxygen-18 variations during a period of 12 months. Observations reveal that the system’s response to recharge depends strongly on the competition between the infiltration and drainage velocities. These velocities, which are controlled by variables such as the time of precipitation, time of infiltration, intensity, and continuity of recharge, determine the degree of dominance of different types of flow mechanisms in the aquifer. Bypass, well-mixed and piston flow mechanisms are used to explain the response of the aquifer to the spatio-temporal variations in recharge. It appears that the aquifer switches among these flow mechanisms depending on the prevailing recharge mode and the competition between infiltration and drainage velocities.     

Comparative application of four methods of groundwater vulnerability mapping in a Slovene karst catchment

Four methods of groundwater vulnerability mapping have been applied in a Slovene karst catchment and validated by tracer tests. The test site is characterised by high water table fluctuations, manifested in intermittent lakes and variable drainage divides. A first multi-tracer test (two injections) allowed subdivision of the catchment into zones of different degrees of contribution (‘inner zone’ and ‘outer zone’). For vulnerability mapping, only methods that consider the specific nature of karst aquifers such as heterogeneity and duality of infiltration processes, were selected: EPIK, PI, the ‘Simplified Method’ and the ‘Slovene Approach’. For validation, a second multi-tracer test (four injections) was carried out. The time of first detection and the normalised recovery were used as validation criteria. The results suggest that EPIK and the Simplified Method sometimes overestimate vulnerability, while PI and the Slovene Approach tend to deliver more realistic results, at least during low-flow conditions. The Slovene Approach gives the clearest guidance on how to deal with hydrologic variability, for example by assigning lower vulnerability to occasionally active sinking surface waters than to permanent ones. As a conclusion, commonly accepted validation techniques are needed and should be applied by default to evaluate different vulnerability mapping methods and the resulting maps.     

Source vulnerability mapping in carbonate (karst) aquifers by extension of the COP method: application to pilot sites

A step from resource to source vulnerability mapping is presented, based on the European COST Action 620 approach for karst groundwater protection. Guidelines on vulnerability assessment of the horizontal groundwater flow path within the karst saturated zone (K factor) are proposed. By integrating this into the previously existing COP method for intrinsic resource vulnerability mapping, adequate source protection can be assessed. The proposed “Karst saturated zone (K) factor” assessment considers groundwater travel time (t subfactor), connection and contribution to the source (r subfactor) and active conduit or fissured network (n subfactor). The extended COP method was applied in two carbonate aquifers in southern Spain with different geological, hydrogeological and climate settings. The results are coherent with previous research results of the studied areas. Moreover, they are consistent with the occasional groundwater contamination detected in one of the springs. On the other hand, an absence of contamination, despite high risk, justifies the lower degrees of vulnerability assigned to the sources surveyed. The source vulnerability maps obtained can thus be used as a basis for the delineation of protection zones.     

Role of karstification and rainfall in the behavior of a heterogeneous karst system

 The complex functioning of karst systems depends on several factors (e.g., geology, karstification status, climate) which influence flow conditions. Understanding a system requires monitoring that includes a sampling interval well adapted to the system's size, and the use of hydrodynamic and geochemical approaches. From our own observations, the general physical and geochemical characteristics of the aquifers are perennial and representative of karst evolution at a given moment. The Notre-Dame-des-Anges karst system, which is situated in France on the rim of the Vaucluse karst region, was studied during 2 water years. This study demonstrates that several parameters are not definitive, and that they do not provide information about the actual status of the karst evolution of the system. We were able to observe that hydrological parameters in particular are dependent on the configuration of the outlets, which can dramatically change. For example, a heavy storm (300 mm of rain in 4 h), which occurred in the intake area on 22 September 1992, produced a piston effect in the aquifer, which was expressed by extensive unclogging of the whole aquifer. The appearance of a new gryphon at the outlet induced a change in the flow rate and added to the unclogging. The other major result of this study is the demonstration of intensive, even preponderant, participation of water originating in the infiltration zone in the behavior of the system. This contribution, often minimized in previous studies, is closely dependent upon flow conditions and especially upon the recharge status of the infiltration zone. Physical and chemical parameters allowed us to trace this water, and to demonstrate the existence of an actual hanging reserve which plays not a negligible part in flow, during both high and low water levels. This paper demonstrates the necessity of coupling a hydrodynamic study with a hydrochemical approach, and of questioning the representativeness of parameters during a given period, for a better understanding of the functioning of the karst system. This is especially true if the infiltration zone of the karst system is heavily developed, and if external factors (e.g., river, landslip, clogging) influence the outlet Received: 10 June 1996 / Accepted: 11 November 1996     

Hydrological modeling in the karst area,Rižana spring catchment,Slovenia

Karst aquifers are known for their heterogeneity and irregular complex flow patterns which make them more difficult to model and demand specific modeling approaches. This paper presents one such approach which is based on a conceptual model. The model was applied in a karst area of the catchment of Rižana spring (200 km2). It is based on the MIKE SHE code and incorporates the main hydrological processes and geological features of the karst aquifer (diffuse and concentrated infiltration, allogenic recharge, quick and slow groundwater flow, shifting groundwater divides and groundwater outflow from the catchment area). Modeling of evapotranspiration and flow in the upper part of the unsaturated zone is more detailed. For the modeling of groundwater flow in the karst aquifer, a conceptual model was applied which uses drainage function for the simulation of groundwater flow through large conduits (karst channels and large fissures). The model was calibrated and validated against the observed Rižana spring discharge which represents a measured response of the aquifer. The results of validation show that the model is able to adequately simulate temporal evolution of the spring discharge, measured by Nash–Sutcliffe coefficient (0.82) as well as overall water balance.     

Simulation of the development of karst aquifers: role of the epikarst

 The evolution of a karst aquifer is modelled taking into account the karst groundwater flow as well as the dissolution kinetics of calcite. In particular, infiltration of water from the epikarst is simulated which controls the temporal and spatial distribution of recharge to the phreatic zone. The results show that the evolution of karst conduits is initiated in the spring. The existence of preferential flow paths leads to the evolution of highly conductive so-called dendritic cave systems, i.e., single passages which concentrate the flow and drain the catchment. With time, the amount of undersaturated water flowing directly into the conduit system is increased leading to an acceleration of the conduits enlargement. Three phases are identified for the evolution of karst aquifers: (a) an initiation stage; (b) an enlargement stage; and (c) a stagnation phase. Received: 4 August 1997 / Accepted: 19 January 1998     

Comparative application of two methods (COP and PaPRIKa) for groundwater vulnerability mapping in Mediterranean karst aquifers (France and Spain)

A comparative test of two vulnerability mapping methods (COP and PaPRIKa) specifically dedicated to for karst aquifers was carried out on two Mediterranean carbonate aquifers. The vulnerability maps obtained for each aquifer present important differences. To identify and determine the origin of these differences, the results were statistically analyzed using sensitivity analysis, coefficients of determination and scatter graphs. In addition, the global vulnerability (Gv) parameter was used to measure the general vulnerability of the aquifer and to compare the results obtained. This statistical analysis led us to conclude that the main cause of differences between these two methods used to assess aquifer vulnerability lie in the relative importance of the parameters employed in calculating the vulnerability index. For the PaPRIKa method, the variable related to infiltration (slope and karst features) has the most influence, with less weight being assigned to the protective capacity of layers overlying the aquifer. For the COP method, the most influent variable is defined by the layers overlying the aquifer, together with infiltration characteristics, determined by the relative importance of different forms of infiltration in each aquifer. The vulnerability mappings performed using the COP method present greater coherence with the known hydrogeological behavior of the study areas, especially the Spanish aquifers. Nevertheless, further hydrogeological investigations are needed, such as ones to validate the obtained vulnerability maps.     

Hydrochemical dynamics of TOC and NO3 − contents as natural tracers of infiltration in karst aquifers

Analysis of Total Organic Carbon (TOC) and NO₃ ^? contents in the water at six springs in the province of Malaga (Southern Spain), sampled under different hydrologic conditions, revealed two different hydrochemical behaviour patterns for these natural tracers of infiltration. TOC content increased during every recharge period, following the rapid arrival of water infiltrating through the soil. On the other hand, NO₃ ^? content only rose during the first flood episodes (normally in autumn), and fell during the winter and spring ones. This difference is consequence of the distinctive biogeochemical kinetics of nitrogen with respect to organic carbon, both in the soil and within the aquifer. Unlike the mineralisation undergone by TOC from the surface to the spring, the NO₃ ^? ion remains in the aquifer almost unaffected, due to the oxidizing conditions prevailing within the karst medium, which do not allow its denitrification and favour its preservation within the saturated zone. In non-polluted aquifers, TOC and NO₃ ^? have a common origin in the soil and can be used to determine infiltration processes and the hydrogeological functioning of karst aquifers. Their different hydrochemical evolution provides information about mineralization and degradation processes of organic matter within karst aquifers, which can be used to validate the vulnerability to contamination in this type of medium.     

岩溶地下水位对降雨响应的时空变异特征及成因探讨——以广西桂林甑皮岩为例

岩溶地下水位对降雨响应具有时空变异性,甑皮岩遗址地下水动力系统结构的认识存在分歧。利用高分辨率降雨水位数据,将研究区分割为不同含水体,通过水位动态、相关分析、滑动窗口采样相关分析等方法,探讨岩溶地下水对降雨响应时空变异特征及成因。结果表明,岩溶强发育、扩散流导水的含水体水位对降雨的响应表现为缓升缓降,水位自相关性强;发育岩溶管道的含水体水位表现为陡升陡降,水位自相关系数衰减速率快,对降雨响应的滞后时间短,互相关函数图呈多峰型;岩溶发育的极不均匀性是造成空间响应差异的主要原因。雨季地下水位对降雨响应的滞后时间远小于枯季;雨季累积降雨量大、水位埋深浅、包气带长期处于饱和或者近饱和状态,降雨垂直入渗补给历时短;雨季暴雨频繁导致含水体地下水短期内形成较大水力梯度,径流补给速度加快。综合分析认为,甑皮岩遗址地下水动力系统由NE向岩溶管道、NS向管道-裂隙以及NE向强径流带3个子径流系统组成。     

Water vulnerability assessment in karst environments: a new method of defining protection areas using a multi-attribute approach and GIS tools (EPIK method)

 Groundwater resources from karst aquifers play a major role in the water supply in karst areas in the world, such as in Switzerland. Defining groundwater protection zones in karst environment is frequently not founded on a solid hydrogeological basis. Protection zones are often inadequate and as a result they may be ineffective. In order to improve this situation, the Federal Office for Environment, Forests and Landscape with the Swiss National Hydrological and Geological Survey contracted the Centre of Hydrogeology of the Neuchatel University to develop a new groundwater protection-zones strategy in karst environment. This approach is based on the vulnerability mapping of the catchment areas of water supplies provided by springs or boreholes. Vulnerability is here defined as the intrinsic geological and hydrogeological characteristics which determine the sensitivity of groundwater to contamination by human activities. The EPIK method is a multi-attribute method for vulnerability mapping which takes into consideration the specific hydrogeological behaviour of karst aquifers. EPIK is based on a conceptual model of karst hydrological systems, which suggests considering four karst aquifer attributes: (1) Epikarst, (2) Protective cover, (3) Infiltration conditions and (4) Karst network development. Each of these four attributes is subdivided into classes which are mapped over the whole water catchment. The attributes and their classes are then weighted. Attribute maps are overlain in order to obtain a final vulnerability map. From the vulnerability map, the groundwater protection zones are defined precisely. This method was applied at several sites in Switzerland where agriculture contamination problems have frequently occurred. These applications resulted in recommend new boundaries for the karst water supplies protection-zones. Received: 27 October 1997 · Accepted: 4 July 1998     

Regional-scale airborne electromagnetic surveying of the Yucatan karst aquifer (Mexico): geological and hydrogeological interpretation

Geometry and connectivity of high-permeability zones determine groundwater flow in karst aquifers. Efficient management of karst aquifers requires regional mapping of preferential flow paths. Remote-sensing technology provides tools to efficiently map the subsurface at such scales. Multi-spectral remote sensing imagery, shuttle radar topography data and frequency-domain airborne electromagnetic (AEM) survey data were used to map karst-aquifer structure on the Yucatan Peninsula, Mexico. Anomalous AEM responses correlated with topographic features and anomalous spectral reflectance of the terrain. One known preferential flow path, the Holbox fracture zone, showed lower bulk electrical resistivity than its surroundings in the AEM surveys. Anomalous structures delineated inland were sealed above by a low-resistivity layer (resistivity: 1–5 Ωm, thickness: 5–6?m). This layer was interpreted as ejecta from the Chicxulub impact (Cretaceous/Paleogene boundary), based on similar resistivity signatures found in borehole logs. Due to limited sensitivity of the AEM survey, the subsurface configuration beneath the low-resistivity layer could not be unambiguously determined. AEM measurements combined with remote-sensing data analysis provide a potentially powerful multi-scale methodology for structural mapping in karst aquifers on the Yucatan Peninsula and beyond.     

Preferential flow,diffuse flow,and perching in an interbedded fractured-rock unsaturated zone

Layers of strong geologic contrast within the unsaturated zone can control recharge and contaminant transport to underlying aquifers. Slow diffuse flow in certain geologic layers, and rapid preferential flow in others, complicates the prediction of vertical and lateral fluxes. A simple model is presented, designed to use limited geological site information to predict these critical subsurface processes in response to a sustained infiltration source. The model is developed and tested using site-specific information from the Idaho National Laboratory in the Eastern Snake River Plain (ESRP), USA, where there are natural and anthropogenic sources of high-volume infiltration from floods, spills, leaks, wastewater disposal, retention ponds, and hydrologic field experiments. The thick unsaturated zone overlying the ESRP aquifer is a good example of a sharply stratified unsaturated zone. Sedimentary interbeds are interspersed between massive and fractured basalt units. The combination of surficial sediments, basalts, and interbeds determines the water fluxes through the variably saturated subsurface. Interbeds are generally less conductive, sometimes causing perched water to collect above them. The model successfully predicts the volume and extent of perching and approximates vertical travel times during events that generate high fluxes from the land surface. These developments are applicable to sites having a thick, geologically complex unsaturated zone of substantial thickness in which preferential and diffuse flow, and perching of percolated water, are important to contaminant transport or aquifer recharge.     

DRASTIC-Fm: a modified vulnerability mapping method for structurally controlled aquifers in the southern Gulf Islands,British Columbia,Canada

DRASTIC, the methodology for mapping the intrinsic vulnerability of aquifers, is modified to incorporate the structural characteristics of fractured bedrock aquifers. In these aquifers, groundwater flow is predominantly through fractures, with large-scale fracture zones and faults acting as primary conduits for flow at the regional scale. The methodology is applied to the southern Gulf Islands region of southwestern British Columbia, Canada. Bedrock geology maps, soil maps, structural measurements, mapped lineaments, water-well information and topographic data, assembled within a comprehensive GIS database, form the basis for assigning traditional DRASTIC indices, while adding the structural indices necessary for capturing the importance of regional structural elements in recharge and well capture zone determinations.     

Use of artificial and natural tracers to assess groundwater transit-time distribution and flow systems in a high-alpine karst system (Wetterstein Mountains,Germany)

Groundwater in mountainous karst regions is vital for regional water budgets and freshwater supply. Owing to increasing water demand and climate change, detailed knowledge of the highly heterogeneous alpine aquifer systems is required. Multi-tracer analyses have been conducted in the steep karstic Wetterstein Mountains, which includes Germany’s highest summit, Zugspitze (2,962 m asl). Results of artificial tracer tests demonstrate well-developed flow paths through the unsaturated zone (up to 1,000 m thickness). Flow paths cross topographic divides and contribute to deep drainage systems underneath alpine valleys. Cross-formational flow has been identified. Quantitative analysis of tailing-dominated breakthrough curves and stable isotopes (¹⁸O) has enabled determination of the mean transit-time distribution. A fast-flow component with transit times between 3 and 13 days was found in karst conduits and open fissures, dependent on flow conditions. An intermediate-flow component, showing mean transit times of about 2.9–4.9 months, was found in well-drained fissures and fractures. A slow-flow component with mean transit times greater than 1 year is attributable to slow flow and low storage in the poorly drained fissures and rock matrix. The conceptual model enables a better understanding of drainage, water resources and vulnerability of the high-alpine karst system.     

The importance of protection and management of Karst water as drinking water resources in Iran

Karst water is one of important potable water especially in arid and semiarid zones. Due to its vulnerability and high sensitivity to pollution, water resources protection is one of important measures to be observed in karst resources management. One case study in the The Maharlu Basin, Iran is discussed. Findings are used to postulate regional guidelines for establishing protection and management measures in important karst pumping fields. Suggested locations for new wells in valleys or depressions should be selected on the basis of specialized hydrogeological investigation. Recommendations for first, second and third degree levels of protection based on European standards are suggested and then defined by model. The first degree protection zone around the pumping field is recommended to be 10 × 10 m². The second degree protection zone should be established to protect primarily from biological pollution. The third degree zone is defined to protect from chemical parameters.     

塔北哈拉哈塘地区奥陶系层间改造区岩溶古水文条件分析

哈拉哈塘地区自北往南可分为潜山岩溶区、层间岩溶顺层改造区、层间岩溶台缘叠加区、层间岩溶斜坡区。层间岩溶-顺层改造区为潜山岩溶区和层间岩溶区过渡区域,文章对前志留纪和良里塔格组岩溶期沉积间断古岩溶水文条件特征进行了详细分析,认为地表河流特征总体表现为自北往南径流,深切河谷决定着局部排泄基准面的深度。在此基础上对地下暗河特征进行刻画,本区地下暗河类型多样,总体表现为流入型和流出型暗河管道系统发育,多为单流入、流出系统,地震反射暗河特征较明显,多为连续或断续串珠强反射,以良里塔格组岩溶期暗河最为发育,且前志留纪岩溶期“断头河”多为部分伏流出口。地下暗河管道附近往往发育大型缝洞体,是目前勘探开发重点研究攻坚对象,因此精确恢复和刻画古岩溶流域水文系统的分布特征,寻找油气藏有利储层是广大岩溶油藏专家所关注的重点。