Improved understanding of dynamic water and mass budgets of high-alpine karst systems obtained from studying a well-defined catchment area

Large areas of Europe, especially in the Alps, are covered by carbonate rocks and in many alpine regions, karst springs are important sources for drinking water supply. Because of their high variability and heterogeneity, the understanding of the hydrogeological functioning of karst aquifers is of particular importance for their protection and utilisation. Climate change and heavy rainfall events are major challenges in managing alpine karst aquifers which possess an enormous potential for future drinking water supply. In this study, we present research from a high-alpine karst system in the UNESCO Biosphere Reserve Großes Walsertal in Austria, which has a clearly defined catchment and is drained by only one spring system. Results show that (a) the investigated system is a highly dynamic karst aquifer with distinct reactions to rainfall events in discharge and electrical conductivity; (b) the estimated transient atmospheric CO₂ sink is about 270 t/a; (c) the calculated carbonate rock denudation rate is between 23 and 47 mm/1000a and (d) the rainfall-discharge behaviour and the internal flow dynamics can be successfully simulated using the modelling package KarstMod. The modelling results indicate the relevance of matrix storage in determining the discharge behaviour of the spring, particularly during low-flow periods. This research and the consequent results can contribute and initiate a better understanding and management of alpine karst aquifers considering climate change with more heavy rainfall events and also longer dry periods.     

Validation of a Numerical Indicator of Microbial Contamination for Karst Springs

Rapid changes in spring water quality in karst areas due to rapid recharge of bacterially contaminated water are a major concern for drinking water suppliers and users. The main objective of this study was to use field experiments with fecal indicators to verify the vulnerability of a karst spring to pathogens, as determined by using a numerical modeling approach. The groundwater modeling was based on linear storage models that can be used to simulate karst water flow. The vulnerability of the karst groundwater is estimated using such models to calculate criteria that influence the likelihood of spring water being affected by microbial contamination. Specifically, the temporal variation in the vulnerability, depending on rainfall events and overall recharge conditions, can be assessed and quantified using the dynamic vulnerability index (DVI). DVI corresponds to the ratio of conduit to diffuse flow contributions to spring discharge. To evaluate model performance with respect to predicted vulnerability, samples from a spring were analyzed for Escherichia coli, enterococci, Clostridium perfringens, and heterotrophic plate count bacteria during and after several rainfall events. DVI was shown to be an indication of the risk of fecal contamination of spring water with sufficient accuracy to be used in drinking water management. We conclude that numerical models are a useful tool for evaluating the vulnerability of karst systems to pathogens under varying recharge conditions     

A high‐altitude temporary spring in a compartmentalized carbonate aquifer: the role of low‐permeability faults and karst conduits

The aim of this research was to refine the actual conceptual model related to the activation of high‐altitude temporary springs within the carbonate Apennines in southern Italy. The research was carried out through geophysical, hydrogeological, hydrochemical and isotopic investigations at the Acqua dei Faggi experimental site during five hydrologic years. The research demonstrated that, in carbonate aquifers where low‐permeability faults cause the aquifer system to be compartmentalized, high‐altitude temporary springs may be recharged by groundwater. In such settings, neither surface water infiltration in karst systems nor perched temporary aquifers play a role of utmost importance. The rare (once or a few time a year) activation of such springs is due to the fact that groundwater unusually reach the threshold head that allows the spring to flow. The activation of the studied high‐altitude temporary spring also depended on relationships between a low‐permeability fault core and a karst system that locally interrupts the low‐permeability barrier. In fact, when the hydraulic head did not reach the karst system, the concentrated head loss within the fault core did not allow the spring to flow, because the groundwater entirely flowed through the fault towards the downgradient compartment. Copyright © 2009 John Wiley & Sons, Ltd.     

Self‐Organized Permeability in Carbonate Aquifers

Advances over the past 40 years have resulted in a clear understanding of how dissolution processes in carbonate rocks enhance aquifer permeability. Laboratory experiments on dissolution rates of calcite and dolomite have established that there is a precipitous drop in dissolution rates as chemical equilibrium is approached. These results have been incorporated into numerical models, simulating the effects of dissolution over time and showing that it occurs along the entire length of pathways through carbonate aquifers. The pathways become enlarged and integrated over time, forming self‐organized networks of channels that typically have apertures in the millimeter to centimeter range. The networks discharge at point‐located springs. Recharge type is an important factor in determining channel size and distribution, resulting in a range of aquifer types, and this is well demonstrated by examples from England. Most carbonate aquifers have a large number of small channels, but in some cases large channels (i.e., enterable caves) can also develop. Rapid velocities found in ground water tracer tests, the high incidence of large‐magnitude springs, and frequent microbial contamination of wells all support the model of self‐organized channel development. A large majority of carbonate aquifers have such channel networks, where ground water velocities often exceed 100 m/d.     

A Waste Treatment/Disposal Site Evaluation Process for Areas Underlain by Carbonate Aquifers

Fractured carbonate or karst aquifers are highly susceptible to contamination. In Minnesota, the areas underlain by karst aquifers are often characterized by thin soil cover, poorly cased wells and/or collapse features. When coupled with improper land use decisions these aquifers have become extensively contaminated.
Evaluation of existing or potentially contaminated sites has required development of unique site evaluation techniques or combinations of techniques. Shallow geophysics have been used successfully to determine subsurface karst features, pinpoint locations for additional investigations and evaluate pump-test data. Soil texture, density and relative moisture content have been closely examined and criteria established to evaluate subsurface erosion susceptibility. Percent rock core recovery, mechanical response of drilling and drilling fluid loss have been used to characterize the secondary porosity and potential contamination susceptibility of the carbonate rock. Dye tracing via surface sinks and well injection have determined primary flow paths and have been useful in estimating flow times.
These techniques used solely or in conjunction have been very useful in evaluating contamination existing in the aquifers and as a predictive tool to evaluate the proposed land uses impacts.     

A Tool for Forecasting the Arrival Time of a Tracer or a Pollutant at a Karst Spring

Karst aquifers are a groundwater resource of global importance. Protection and contamination of karst water resources are often tackled by means of land-use planning and monitoring of physico-chemical and microbial water properties coupled to early-warning systems. However, these solutions of management do not integrate observational data on karst hydrodynamic, for example, tracer tests. Thus, this technical note aims at introducing a straightforward method for forecasting the arrival time of a tracer or a pollutant at a karst spring. Deliverables from the proposed method can assist professionals in the case of a pollution involving a swallow hole of the karst network. The approach is based on repeated tracer testing, gauging of spring flow rate and mathematical curve fitting. Tracer tests and gauging of spring flow rates are repeated under low, middle, and high flow periods. Arrival times of tracer from a given swallow hole of the network are plotted as a function of measured flow rates at the karst spring. The observational data are finally fitted by a mathematical function. Flow rate monitoring at high-temporal resolution at the karst spring and the mathematical functions allow forecasting the arrival of a tracer or a pollutant at the spring in the case of an accident involving a swallow hole of the network or estimating the duration of a pollution if hydrodynamic conditions do not change abruptly. The applicability of the method is documented by three study sites in Switzerland.     

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

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

Simulation of spring flows from a karst aquifer with an artificial neural network

In China, 9·5% of the landmass is karst terrain and of that 47,000 km² is located in semiarid regions. In these regions the karst aquifers feed many large karst springs within basins of thousands of square kilometres. Spring discharges reflect the fluctuation of ground water level and variability of ground water storage in the basins. However, karst aquifers are highly heterogeneous and monitoring data are sparse in these regions. Therefore, for sustainable utilization and conservation of karst ground water it is necessary to simulate the spring flows to acquire better understanding of karst hydrological processes. The purpose of this study is to develop a parsimonious model that accurately simulates spring discharges using an artificial neural network (ANN) model. The karst spring aquifer was treated as a non‐linear input/output system to simulate the response of karst spring flow to precipitation and applied the model to the Niangziguan Springs, located in the east of Shanxi Province, China and a representative of karst springs in a semiarid area. Moreover, the ANN model was compared with a previous time‐lag linear model and it was found that the ANN model performed better. Copyright © 2007 John Wiley & Sons, Ltd.     

Colonization by aerobic bacteria in karst: laboratory and in situ experiments

Experiments were carried out to investigate the potential for bacterial colonization of different substrates in karst aquifers and the nature of the colonizing bacteria. Laboratory batch experiments were performed using limestone and PVC as substrates, a natural bacterial isolate and a known laboratory strain (Escherichia coli [E. coli]) as inocula, and karst ground water and a synthetic formula as growth media. In parallel, fragments of limestone and granite were submerged in boreholes penetrating two karst aquifers for more than one year; the boreholes are periodically contaminated by enteric bacteria from waste water. Once a month, rock samples were removed and the colonizing bacteria quantified and identified. The batch experiments demonstrated that the natural isolate and E. coli both readily colonized limestone surfaces using karst ground water as the growth medium. In contrast, bacterial colonization of both the limestone and granite substrates, when submerged in the karst, was less intense. More than 300 bacterial strains were isolated over the period sampled, but no temporal pattern in colonization was seen as far as strain, and colonization by E. coli was notably absent, although strains of Salmonella and Citrobacter were each observed once. Samples suspended in boreholes penetrating highly fractured zones were less densely colonized than those in the borehole penetrating a less fractured zone. The results suggest that contamination of karst aquifers by enteric bacteria is unlikely to be persistent. We hypothesize that this may be a result of the high flow velocities found in karst conduits, and of predation of colonizing bacteria by autochthonous zooplankton.     

Dissolution in a variably confined carbonate platform: effects of allogenic runoff,hydraulic damming of groundwater inputs,and surface–groundwater exchange at the basin scale

In variably confined carbonate platforms, impermeable confining units collect rainfall over large areas and deliver runoff to rivers or conduits in unconfined portions of platforms. Runoff can increase river stage or conduit heads in unconfined portions of platforms faster than local infiltration of rainfall can increase groundwater heads, causing hydraulic gradients between rivers, conduits and the aquifer to reverse. Gradient reversals cause flood waters to flow from rivers and conduits into the aquifer where they can dissolve limestone. Previous work on impacts of gradient reversals on dissolution has primarily emphasized individual caves and little research has been conducted at basin scales. To address this gap in knowledge, we used legacy data to assess how a gradient of aquifer confinement across the Suwannee River Basin, north‐central Florida affected locations, magnitudes and processes of dissolution during 2005–2007, a period with extreme ranges of discharge. During intense rain events, runoff from the confining unit increased river stage above groundwater heads in unconfined portions of the platform, hydraulically damming inputs of groundwater along a 200 km reach of river. Hydraulic damming allowed allogenic runoff with SI_CAL < ?4 to fill the entire river channel and flow into the aquifer via reversing springs. Storage of runoff in the aquifer decreased peak river discharges downstream and contributed to dissolution within the aquifer. Temporary storage of allogenic runoff in karst aquifers represents hyporheic exchange at a scale that is larger than found in streams flowing over non‐karst aquifers because conduits in karst aquifers extend the area available for exchange beyond river beds deep into aquifers. Post‐depositional porosity in variably confined carbonate platforms should thus be enhanced along rivers that originate on confining units. This distribution should be considered in models of porosity distribution used to manage water and hydrocarbon resources in carbonate rocks. Copyright © 2013 John Wiley & Sons, Ltd.     

Estimation of groundwater vulnerability to pollution based on DRASTIC in the Niipele sub-basin of the Cuvelai Etosha Basin,Namibia

Surface water is a scarce resource in Namibia with about sixty percent of Namibia's population dependent on groundwater for drinking purposes. With increasing population, the country faces water challenges and thus groundwater resources need to be managed properly. One important aspect of Integrated Water Resources Management is the protection of water resources, including protection of groundwater from contamination and over-exploitation. This study explores vulnerability mapping as a basic tool for protecting groundwater resources from pollution. It estimates groundwater vulnerability to pollution in the upper Niipele sub-basin of the Cuvelai-Etosha in Northern Namibia using the DRASTIC index. The DRASTIC index uses GIS to estimate groundwater vulnerability by overlaying different spatially referenced hydrogeological parameters that affect groundwater contamination. The study assesses the discontinuous perched aquifer (KDP) and the Ohangwena multi-layered aquifer 1 (KOH-1). For perched aquifers, point data was regionalized by a hydrotope approach whereas for KOH-1 aquifer, inverse distance weighting was used. The hydrotope approach categorized different parts of the hydrogeological system with similar properties into five hydrotopes. The result suggests that the discontinuous perched aquifers are more vulnerable than Ohangwena multi-layered aquifer 1. This implies that vulnerability increases with decreasing depth to water table because contaminants have short travel time to reach the aquifer when they are introduced on land surface. The nitrate concentration ranges between 2 and 288 mg/l in perched aquifers while in Ohangwena multi-layered aquifer 1, it ranges between 1 and 133 mg/l. It was observed that perched aquifers have high nitrate concentrations than Ohangwena 1 aquifer, which correlates well with the vulnerability results.     

Formation of phreatic caves in an eogenetic karst aquifer by CO2 enrichment at lower water tables and subsequent flooding by sea level rise

Formation of extensive phreatic caves in eogenetic karst aquifers is widely believed to require mixing of fresh and saltwater. Extensive phreatic caves also occur, however, in eogenetic karst aquifers where fresh and saltwater do not mix, for example in the upper Floridan aquifer. These caves are thought to have formed in their modern settings by dissolution from sinking streams or by convergence of groundwater flow paths on springs. Alternatively, these caves have been hypothesized to have formed at lower water tables during sea level low‐stands. These hypotheses have not previously been tested against one another. Analyzing morphological data and water chemistry from caves in the Suwannee River Basin in north‐central Florida and water chemistry from wells in the central Florida carbonate platform indicates that phreatic caves within the Suwannee River Basin most likely formed at lower water tables during lower sea levels. Consideration of the hydrological and geochemical constraints posed by the upper Floridan aquifer leads to the conclusion that cave formation was most likely driven by dissolution of vadose CO₂ gas into the groundwater. Sea level rise and a wetter climate during the mid‐Holocene lifted the water table above the elevation of the caves and placed the caves tens of meters below the modern water table. When rising water tables reached the land surface, surface streams formed. Incision of surface streams breached the pre‐existing caves to form modern springs, which provide access to the phreatic caves. Phreatic caves in the Suwannee River Basin are thus relict and have no causal relationship with modern surficial drainage systems. Neither mixing dissolution nor sinking streams are necessary to form laterally extensive phreatic caves in eogenetic karst aquifers. Dissolution at water tables, potentially driven by vadose CO₂ gas, offers an underappreciated mechanism to form cavernous porosity in eogenetic carbonate rocks. Copyright © 2012 John Wiley & Sons, Ltd.     

Investigation of karstic hydrological processes of Niangziguan Springs (North China) using wavelet analysis

Niangziguan Spring complex is the largest karst spring in North China. We investigate the karst hydrological processes by using Morlet wavelet transform analysis and cross wavelet analysis based on monthly precipitation from 1958 to 2010 and spring discharge from 1958 to 2009. From Morlet wavelet transform coefficients of precipitation and the spring discharge in Niangziguan Springs Basin, we find that the precipitation and discharge are characterized by the multi‐scale features in the time domain, and the energy distribution of the signal is highly irregular across scales. Although precipitation eventually becomes spring discharge by infiltrating and propagating through karst formations, the signals are attenuated. The results also show that the precipitation of Niangziguan Springs Basin has the main periodic components of 1‐, 5‐, 12‐, and 17‐year periods with alternating wet–drought cycle. Similarly, the spring discharge of Niangziguan Springs has the main components of 17‐year periods, but the 1‐, 5‐, and 12‐year periodicity of precipitation are not reflected in spring discharge, which is filtered by the aquifers. The results of cross wavelet analysis reveal that the precipitation and spring discharge share the common periodicity of 17 years. This means that those signals with high energy and long timescales can penetrate through the aquifer and be reflected in spring discharge, whereas other signals are filtered and modified. Copyright © 2011 John Wiley & Sons, Ltd.     

Natural responses of Neoproterozoic dynamic karst springs to rainfall events,São Miguel Watershed,Minas Gerais,Brazil

Karst aquifers consist of complex networks of conduits in which groundwater flows and recharge/discharge processes are generally more dynamic than in other types of aquifers. Due to their intrinsic heterogeneity and anisotropy, monitoring, quantifying, and analysing natural responses of karst springs is an efficient tool. Unlike Cenozoic and Mesozoic rocks, in Neoproterozoic karst systems, groundwater circulates and stores generally in dissolution features known as tertiary porosity, as the rock's primary porosity is recrystallized, considered negligible. This article studies the hydrodynamics of a karst portion of the São Miguel River basin, southwest of the state of Minas Gerais, Brazil. The region is predominantly composed of Neoproterozoic carbonate rocks, dating from about 570 to 540 million years ago. During a hydrological year (2019–2020), three karst springs (S1, S2, and S3) were daily monitored through their natural responses (variations of electrical conductivity, EC, temperature, T, and discharge, Q) to rainfall episodes. The data were interpreted based on the analysis of spring hydrographs, time series, recession curves (seasonal and intra-annual), and statistics of EC, T, and Q variations. The results show the three springs generally exhibit quick flow, typically karstic, in the case of hydrosystems with a well-structured and functional underground drainage network. The time series indicate the hydrosystem drained by S1 presents slower circulation and a lower degree of linearity, resulting from the higher sinuosity of the system, while the hydrosystems of S2 and S3 have similar behaviours, of quick water circulations immediately after a rainy episode. The degrees of karstification classify S1 and S2 as complex and extensive karst systems consisting of several subsystems, and S3 as a system in which the conduit network is more developed at the upper epiphreatic zone than near the outlet.     

Interpretation of groundwater level monitoring results in karst aquifers: examples from the Dinaric karst

The paper presents an attempt to determine the characteristics of karst aquifers using information on groundwater level (GWL) in natural holes and boreholes with different data quantity and time resolution of GWL measurements. In this paper the particulars of karst aquifers were analysed for four examples from the Dinaric karst. In all four study areas, aquifers are formed in bare, deep and well‐developed Dinaric karst consisting of Cretaceous limestones. The first example represents a wide area of Imotsko polje in the karst. The aquifer was analysed on the basis of infrequent water level monitoring in natural karst water features (jamas, lakes, wells) and discharges of springs and rivers. The karst aquifer in this example is complex, non‐homogenous and variable in space and time, which is frequent in the Dinaric karst. Regardless of the aforementioned it was possible to determine its elementary characteristics. The second example represents 10 wells used for the water supply for the city of Pula. The GWL and salinity were measured once a week in the period between 1981 and 1996. Even though these measurements were relatively infrequent in space and time, they served as bases for assessment of average and maximum aquifer conditions as well as boundaries of saltwater intrusion. In the third example only a portion of aquifer of the karst spring Blaz, which is in the contact with the Adriatic Seas, has been analyzed. It is a spring with an intrusion of salt water. For purposes of study of saltwater intrusion, 26 piezometers were drilled in its vicinity in which GWL, salinity and temperature were measured once a day during 168 days, a period comprising one complete cycle of seawater intrusion and retreat. These measurements proved the existence of dispersed discharge from the aquifer into the sea and its non‐homogeneity in space. In the fourth example GWL was measured continuously in 10 deep (up to 300 m) piezometers in the hinterland of the Ombla Spring catchment. The measurement period lasted 2 years (January 1988 to December 1989). The analyses are made with hourly data. The results made it possible to determine numerous characteristics of the karst aquifer and a significant non‐homogeneity of groundwater distribution in karst aquifers, depending more on the underground karst phenomena than the surface karst forms. Copyright © 2000 John Wiley & Sons, Ltd.     

Characterizing the heterogeneity of karst critical zone and its hydrological function: An integrated approach

Spatial heterogeneity in the subsurface of karst environments is high, as evidenced by the multiphase porosity of carbonate rocks and complex landform features that result in marked variability of hydrological processes in space and time. This includes complex exchange of various flows (e.g., fast conduit flows and slow fracture flows) in different locations. Here, we integrate various “state‐of‐the‐art” methods to understand the structure and function of this poorly constrained critical zone environment. Geophysical, hydrometric, and tracer tools are used to characterize the hydrological functions of the cockpit karst critical zone in the small catchment of Chenqi, Guizhou Province, China. Geophysical surveys, using electrical resistivity tomography (ERT), inferred the spatial heterogeneity of permeability in the epikarst and underlying aquifer. Water tables in depression wells in valley bottom areas, as well as discharge from springs on steeper hillslopes and at the catchment outlet, showed different hydrodynamic responses to storm event rainwater recharge and hillslope flows. Tracer studies using water temperatures and stable water isotopes (δD and δ¹⁸O) could be used alongside insights into aquifer permeability from ERT surveys to explain site‐ and depth‐dependent variability in the groundwater response in terms of the degree to which “new” water from storm rainfall recharges and mixes with “old” pre‐event water in karst aquifers. This integrated approach reveals spatial structure in the karst critical zone and provides a conceptual framework of hydrological functions across spatial and temporal scales.     

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

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

Heterogeneous aquitard properties in sedimentary successions in the Apennine chain: case studies in southern Italy

The aim of this study was to analyse the hydrogeologic role of siliciclastic deposits and marly‐clayey‐calcareous successions within the carbonate Apennine chain (southern Italy). The study was carried out along the northern part of the Matese carbonate massif through (1) the hydraulic characterization of siliciclastic rocks in a test site, by means of Lugeon tests, and (2) the identification of the groundwater flow system discharging at an important spring located within a marly‐clayey‐calcareous succession in a second test site, by means of isotopic investigations. The results showed that the investigated siliciclastic deposits and marly‐clayey‐calcareous successions may allow significant groundwater discharge from carbonate aquifers. Thus, they do not everywhere behave as aquitard, contrary to the previous model. Instead, groundwater flows through the upper part of these successions, where stress release fracturing enhanced rock permeability in the near‐surface bedrock. Thus, these successions may locally be a new groundwater source within the southern Apennine chain. Copyright © 2009 John Wiley & Sons, Ltd.     

Coastal and inland karst morphologies driven by sea level stands: a GIS based method for their evaluation

Sea level is the base level for groundwater circulation in coastal aquifers. The evolution of karst surface landforms and subsurface drainage systems in these aquifers has been conditioned in geological time by tectonics and glacio‐eustatic sea‐level changes. Present morpho‐structural settings and the type/distribution of karst surface and subsurface forms have developed in different carbonate formations according to differences in lithology, climate and exposure time, all driving the intensity of morphologic and karst processes. The repeated and significant changes of groundwater level linked to ‘sea‐level changes’ have had the most important role in driving the continuous evolution of karstic drainage systems, and has resulted in most cases in a multiphase karst. This study aims at defining a general method for identifying, in karst coastal settings, the elevations of flat or low topographic gradient surfaces (using morphometric analysis of Digital Elevation Models (DEMs) and geographical information systems (GISs), and their comparison with elevations of distinctive karstic levels (passages, lateral solution cavities) observed in vertical shafts and horizontal caves. Of the elevations of flat or low topographic gradient surfaces only those agreeing, within ±10 m or ±20 m, with elevation ranges marked by the high frequency of distinctive karst levels were considered as representative of the more probable past sea‐level stands. The method is applied to a regional coastal carbonate formation in southern Italy, by using a 10 m DEM and information on 140 complex caves and 85 shafts. Of the 15 elevations indicated by DEM analysis [620, 600, 470, 450, 425, 385, 355, 315, 270, 250, 205, 180, 150, 110, and 70 m above sea level (a.s.l.)], 13 match clearly those highlighted by significant frequencies of distinctive karstic levels. These elevations are validated by comparison to the elevation of terraces and karst plains indicated in the literature. Copyright © 2012 John Wiley & Sons, Ltd.