Characterizing moldic and vuggy pore space in karst aquifers using borehole-wall,slabbed-core and thin-section images

Carbonate aquifers are prolific and important sources of potable water in many parts of the world owing to enlarged dissolution features that enhance porosity and interconnectivity. To better understand the variations of pore space in different karst aquifers, image and geospatial analyses are used to analyze pore attributes (i.e., pore area and perimeter) in images of vuggy aquifers. Pore geometry and 2D porosity derived from images of the moldic Castle Hayne and vuggy Biscayne aquifers are analyzed at three scales of observation: borehole televiewer, core and thin-section. The Castle Hayne and Biscayne aquifers are the foci of this study because the pore spaces that control the hydrologic properties in each of these aquifers are markedly different even though both of these carbonate reservoirs are prolific aquifers. Assessments of pore area, perimeter and shape index (a measure of shape complexity) indicate that pore geometries and pore complexities vary as a function pore type and scale of observation. For each aquifer type, the areas, perimeters and complexities of pores are higher at the larger scale of observation (e.g., borehole) than the smaller scale of observation (e.g., thin section). When the complexity of the moldic pores is compared to the complexity of vuggy pores, the results indicate that moldic pores are generally more complex than vuggy pores at the same scale of observation. Whereas estimates of 2D porosity from the borehole televiewer image of the vuggy aquifer are higher than those derived from the moldic aquifer, the range of 2D porosities is larger in core and thin section images for the vuggy aquifer than the moldic aquifer. A model for the development of pores is presented that suggests that the coalescence of small pores with simple shapes leads to the growth of larger pores with more complex shapes. The model suggests that the younger Biscayne aquifer is a more mature karst than the Castle Hayne aquifer.     

Geostatistical Borehole Image‐Based Mapping of Karst‐Carbonate Aquifer Pores

Quantification of the character and spatial distribution of porosity in carbonate aquifers is important as input into computer models used in the calculation of intrinsic permeability and for next‐generation, high‐resolution groundwater flow simulations. Digital, optical, borehole‐wall image data from three closely spaced boreholes in the karst‐carbonate Biscayne aquifer in southeastern Florida are used in geostatistical experiments to assess the capabilities of various methods to create realistic two‐dimensional models of vuggy megaporosity and matrix‐porosity distribution in the limestone that composes the aquifer. When the borehole image data alone were used as the model training image, multiple‐point geostatistics failed to detect the known spatial autocorrelation of vuggy megaporosity and matrix porosity among the three boreholes, which were only 10 m apart. Variogram analysis and subsequent Gaussian simulation produced results that showed a realistic conceptualization of horizontal continuity of strata dominated by vuggy megaporosity and matrix porosity among the three boreholes.     

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.     

Review of Laboratory Scale Models of Karst Aquifers: Approaches,Similitude, and Requirements

This review focuses on investigations of groundwater flow and solute transport in karst aquifers through laboratory scale models (LSMs). In particular, LSMs have been used to generate new data under different hydraulic and contaminant transport conditions, testing of new approaches for site characterization, and providing new insights into flow and transport processes through complex karst aquifers. Due to the increasing need for LSMs to investigate a wide range of issues, associated with flow and solute migration karst aquifers this review attempts to classify, and introduce a framework for constructing a karst aquifer physical model that is more representative of field conditions. The LSMs are categorized into four groups: sand box, rock block, pipe/fracture network, and pipe-matrix coupling. These groups are compared and their advantages and disadvantages highlighted. The capabilities of such models have been extensively improved by new developments in experimental methods and measurement devices. Newer technologies such as 3D printing, computed tomography scanning, X-rays, nuclear magnetic resonance, novel geophysical techniques, and use of nanomaterials allow for greater flexibilities in conducting experiments. In order for LSMs to be representative of karst aquifers, a few requirements are introduced: (1) the ability to simulate heterogeneous distributions of karst hydraulic parameters, (2) establish Darcian and non-Darcian flow regimes and exchange between the matrix and conduits, (3) placement of adequate sampling points and intervals, and (4) achieving some degree of geometric, kinematic, and dynamic similitude to represent field conditions.     

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.     

Comparison of flowpaths to a well and spring in a karst aquifer

The permeability of some karst aquifers consists of networks of poorly integrated conduits and dissolution-widened fractures. The flow includes conduit flow, especially during storm recharge, but lacks the focused recharge into single master conduits that occurs in more highly developed karst systems. The proportions of conduit and dispersed flow are difficult to quantify in such systems. This study examines the flowpaths in a small karst watershed, based on comparing the physical and chemical response to storm flow at both a spring and a well. By conducting continuous monitoring at both locations, a better understanding of the flowpaths in a poorly integrated network was obtained. A more permeable flowpath to the spring leads to faster storm response and lower ion concentrations. The flowpath to and from the well is more complicated. The higher ion content and slower storm response suggest slower, more dispersed flowpaths. However, the well has greater variation in ion chemistry. Periodic recharge may dilute well concentrations due to faster (conduit or fracture) flowpaths. Although karst systems such as this are difficult to characterize, applying a variety of geochemical and physical monitoring techniques at multiple locations illustrates that the flowpaths can vary in both space and time.     

A geographic data model for representing ground water systems

The Arc Hydro ground water data model is a geographic data model for representing spatial and temporal ground water information within a geographic information system (GIS). The data model is a standardized representation of ground water systems within a spatial database that provides a public domain template for GIS users to store, document, and analyze commonly used spatial and temporal ground water data sets. This paper describes the data model framework, a simplified version of the complete ground water data model that includes two-dimensional and three-dimensional (3D) object classes for representing aquifers, wells, and borehole data, and the 3D geospatial context in which these data exist. The framework data model also includes tabular objects for representing temporal information such as water levels and water quality samples that are related with spatial features.     

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.     

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.     

Insights into agricultural influences and weathering processes from major ion patterns

Karst areas and their catchments pose a great challenge for protection because fast conduit flow results in low natural attenuation of anthropogenic contaminants. Studies of the hydrochemistry of karst sources and river solutes are an important tool for securing and managing water resources. A study of the geochemical downriver evolution of the Wiesent River and its tributaries, located in a typical karst terrain, revealed unexpected downstream decreases of nitrate with maximum mean values of 30 mg/L at the source to minimum values of 18 mg/L near the river mouth. This trend persisted over the length of the river even though increased agricultural activities are evident in the downstream section of the catchment. This pattern is caused by fertilizer inputs via diffusive and fast conduits flow from karst lithology in the upstream area that may have reached the river's source even from beyond the hydrological catchment boundaries. Further downstream, these influences became diluted by tributary inputs that drain subcatchments dominated by claystone and sandstone lithologies that increased potassium and sulphate concentrations. Our findings indicate that bedrock geology remains the dominant control on the major ion chemistry of the Wiesent River and that agricultural influences are strongest near the headwaters despite increased land use further downstream, due to long‐term storage and accumulation in karst aquifers. This feature may not be unique to the Wiesent River system, as carbonates cover significant portions of the Earth's surface and subsequent work in other river systems could establish whether such patterns are ubiquitous worldwide.     

Validating the KAGIS black‐box GIS‐based model in a Mediterranean karst aquifer: Case of study of Mela aquifer (SE Spain)

Karst Aquifer GIS‐based model (KAGIS model) is developed and applied to Mela aquifer, a small karst aquifer located in a Mediterranean region (SE Spain). This model considers different variables, such as precipitation, land use, surface slope and lithology, and their geographical heterogeneity to calculate both, the run‐off coefficients and the fraction of precipitation which contributes to fill the soil water reservoir existing above the aquifer. Evapotranspiration uptakes deplete water, exclusively, from this soil water reservoir and aquifer recharge occurs when water in the soil reservoir exceeds the soil field capacity. The proposed model also obtains variations of the effective porosity in a vertical profile, an intrinsic consequence of the karstification processes. A new proposal from the Nash–Sutcliffe efficiency index, adapted to arid environments, is presented and employed to evaluate the model's ability to predict the water table oscillations. The uncertainty in the model parameters is determined by the Generalized Likelihood Uncertainty Estimation method. Afterwards, when KAGIS is calibrated, wavelet analysis is applied to the resulting data in order to evaluate the variability in the aquifer behaviour. Wavelet analysis reveals that the rapid hydrogeological response, typical of a wide variety of karst systems, is the prevailing feature of Mela aquifer. This study proves that KAGIS is a useful tool to quantify recharge and discharge rates of karst aquifers and can be effectively applied to develop a proper management of water resources in Mediterranean areas.     

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

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

Storage and Drainage Characteristics of a Highly Heterogeneous Karst Aquifer in Houzhai Basin

The characteristics of karst aquifers are difficult to be determined due to their heterogeneous physical properties and lack of hydrogeological information. In this case study, we applied two methods for a comparative analysis of storage and drainage characteristics in upstream, midstream, and downstream of Houzhai cave stream basin. In the first method, Minimum Smoothed Method (MSM) is used to determine the proportion of baseflow to the total flow (Baseflow Index, BFI). In the second method, a bicarbonate‐base two‐end member mixing model is used to quantify the slow flow component and fast flow component. For both methods, slow flow and quick flow are quantified at three sampling sites, which provide useful information for the analysis of storage and drainage characteristics. The results from flow separation method and hydrogeochemical analysis show a consistently increasing trend of the proportion of slow flow to total flow from the upstream to downstream which indicates that the voids of highly conductive conduits and well‐connected fissures decrease along the flow paths in the Houzhai cave stream basin in southwest China. The upstream areas have a low proportion of baseflow which indicates a high drainage capacity due to high permeable conduits and well‐connected fissures. The downstream areas, on the contrary, have a high proportion of baseflow which indicates a high storage capacity and slow infiltration due to the predominant presence of matrix and poorly‐connected fissures. These numerical methods provide alternative ways to investigate the storage and drainage characteristics of karst aquifers where direct measurement are not available.     

Application of a Discrete-Continuum Model to Karst Aquifers in North China

A generalized discrete-continuum model is developed to simulate ground water flow in the karst aquifers of North China. The model is a hybrid numerical flow model, which takes into account both quick conduit flow and diffusive fissure flow. The conduit flow is represented by a discrete network model, and the fissure flow is modeled by a continuum approach. The developed model strongly emphasizes the function of the conduits in the flow fields. They control the general drainage pattern, as demonstrated in the simulation of a complex karst aquifer in North China. The model reproduces reasonably well the flow field in response to an unanticipated discharge of ground water from the karst aquifer into an underground mine based on the aquifer parameters that are manually calibrated from a multiple-well pumping test. Sensitivity of the model to the aquifer parameters was evaluated in the context of the case study.     

Use of shallow groundwater temperature profiles to infer climate and land use change: interpretation and measurement challenges

Temperature is often used to infer the effect of land use and climate conditions on aquifers. Reliable data are needed to examine the temperature behaviour in the subsurface; thus, the use of robust acquisition techniques is unavoidable. Three temperature measurement techniques were applied to assess the sources of bias that could occur during temperature logging in a shallow Quaternary coastal aquifer in Ferrara (Northern Italy). Open borehole temperature logging, multilevel sampling straddle packers isolated temperature measurements within a flow cell above ground and multilevel sampling straddle packers isolated temperature measurements via an in‐well level logger (MLS‐IW) were compared for several coastal monitoring wells to gain insights on the limitations of each technique. Results show that the source of bias between the three applied techniques are different: (i) the open borehole temperature logging method tends to record heat convection through the open borehole and is not representative of the aquifer temperature distribution; (ii) the multilevel sampling straddle packers isolated temperature measurements within a flow cell above ground method is swayed by the air temperature and the heating of the submersible pump used to lift groundwater above ground; and (iii) the MLS‐IW provides the most reliable vertical thermal profiling both in summer and winter, because groundwater temperature is directly measured at the selected monitoring depth. The implementation of a 1D flow model demonstrates that if precise temperature profiles are needed to infer the influence that land use and climate changes have on groundwater, the MLS‐IW method is a reliable method that could be applied to existing monitoring wells. Copyright © 2016 John Wiley & Sons, Ltd.     

Geostatistics applied to cross-well reflection seismic for imaging carbonate aquifers

Cross-well seismic reflection data, acquired from a carbonate aquifer at Port Mayaca test site near the eastern boundary of Lake Okeechobee in Martin County, Florida, are used to delineate flow units in the region intercepted by two wells. The interwell impedance determined by inversion from the seismic reflection data allows us to visualize the major boundaries between the hydraulic units. The hydraulic (flow) unit properties are based on the integration of well logs and the carbonate structure, which consists of isolated vuggy carbonate units and interconnected vug systems within the carbonate matrix. The vuggy and matrix porosity logs based on Formation Micro-Imager (FMI) data provide information about highly permeable conduits at well locations. The integration of the inverted impedance and well logs using geostatistics helps us to assess the resolution of the cross-well seismic method for detecting conduits and to determine whether these conduits are continuous or discontinuous between wells. A productive water zone of the aquifer outlined by the well logs was selected for analysis and interpretation. The ELAN (Elemental Log Analysis) porosity from two wells was selected as primary data and the reflection seismic-based impedance as secondary data. The direct and cross variograms along the vertical wells capture nested structures associated with periodic carbonate units, which correspond to connected flow units between the wells. Alternatively, the horizontal variogram of impedance (secondary data) provides scale lengths that correspond to irregular boundary shapes of flow units. The ELAN porosity image obtained by cokriging exhibits three similar flow units at different depths. These units are thin conduits developed in the first well and, at about the middle of the interwell separation region, these conduits connect to thicker flow units that are intercepted by the second well. In addition, a high impedance zone (low porosity) at a depth of about 275 m, after being converted to ELAN porosity, is characterized as a more confined low porosity structure. This continuous zone corresponds to a permeability barrier in the carbonate aquifer that separates the three connected conduits observed in the cokriging image. In the zones above and below this permeability barrier, the water production is very high, which agrees with water well observations at the Port Mayaca aquifer.     

Hydraulically Controlled Discrete Sampling from Open Boreholes

Groundwater sampling from open boreholes in fractured‐rock aquifers is particularly challenging because of mixing and dilution of fluid within the borehole from multiple fractures. This note presents an alternative to traditional sampling in open boreholes with packer assemblies. The alternative system called ZONFLO (zonal flow) is based on hydraulic control of borehole flow conditions. Fluid from discrete fractures zones are hydraulically isolated allowing for the collection of representative samples. In rough‐faced open boreholes and formations with less competent rock, hydraulic containment may offer an attractive alternative to physical containment with packers. Preliminary test results indicate a discrete zone can be effectively hydraulically isolated from other zones within a borehole for the purpose of groundwater sampling using this new method.     

Soil Gas Screening for Chlorinated Solvents at Three Contaminated Karst Sites in Tennessee

Soil gas was sampled using active sampling techniques and passive collectors at three sites in Tennessee to evaluate the effectiveness of these techniques for locating chlorinated solvent sources and flowpaths in karst aquifers. Actively collected soil gas samples were analyzed in the field with a portable gas chromatograph, and the passive soil gas collectors were analyzed in the lab with gas chromatography/mass spectrometry. Results of the sampling indicate that the effectiveness of both techniques is highly dependent on the distribution of the contaminants in the subsurface, the geomorphic and hydrogeologic characteristics of the site, and, in one case, on seasonal conditions. Both active and passive techniques identified areas of elevated subsurface chlorinated solvent concentrations at a landfill site where contamination remains concentrated in the regolith. Neither technique detected chlorinated solvents known to be moving in the bedrock at a manufacturing site characterized by thick regolith and an absence of surficial karst features. Passive soil gas sampling had varied success detecting flowpaths for chloroform in the bedrock at a train derailment site characterized by shallow regolith and abundant surficial karst features. At the train derailment site, delineation of the contaminant flowpath through passive soil gas sampling was stronger and more detailed under winter conditions than summer.     

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.     

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.