Alpine Karst Research Articles

Travertine deposition rather than tourism activity is the primary contributor to the microplastic risks in alpine karst lakes

Microplastics (MPs) are emerging as anthropogenic vectors to form plastisphere, facilitating microbiome colonization and pathogenic dissemination, thus contributing to environmental and health crises across various ecosystems. However, a knowledge gap persists regarding MPs risks and their driving factors in certain unique and vulnerable ecosystems, such as Karst travertine lakes, some of which are renowned World Natural Heritage Sites under ever-increasing tourism pressure. We hypothesized that tourism activities serve as the most important factor of MPs pollution, whereas intrinsic features, including travertine deposition can exacerbate potential environmental risks. Thus, metagenomic approaches were employed to investigate the geographical distribution of the microbiome, antibiotic resistance genes (ARGs), virulence factor genes (VFGs), and their combined environmental risks in Jiuzhaigou and Huanglong, two famous tourism destinations in Southwest China. The plastisphere risks were higher in Huanglong, contradicting our hypothesis that Jiuzhaigou would face more crucial antibiotic risks due to its higher tourist activities. Specifically, the levels of Lipopolysaccharide Lewis and fosD increased by sevenfold and 20-fold, respectively, from upstream to downstream in Huanglong, whereas in Jiuzhaigou, no significant accrual was observed. Structural equation modeling results showed that travertine deposition was the primary contributor to MPs risks in alpine karstic lakes. Our findings suggest that tourism has low impact on MPs risks, possibly because of proper management, and that travertine deposition might act as an MPs hotspot, emphasizing the importance of considering the unique aspects of travertine lakes in mitigating MPs pollution and promoting the sustainable development of World Natural Heritage Sites.

Erratum to: Recent podocopid ostracods from an alpine karst lake in Turkey

Erratum to: Recent podocopid ostracods from an alpine karst lake in Turkey

Unraveling the functioning of the vadose zone in alpine karst aquifers: New insights from a tracer test in the Migovec cave system (Julian alps, NW Slovenia)

The aquifers of alpine karst and high karst plateaus are abundant water resources. They are difficult to characterise due to their complex, partly glaciokarstic, evolution in active tectonic environments, and an unsaturated zone up to two kilometres thick. We present and discuss the results of a tracing test in the alpine karst of the Julian Alps (Slovenia), more precisely in the Migovec System, the longest cave system in Slovenia (length = 43 km, depth = 972 m). The cave extends below a mountain ridge that separates the Soča and Sava Valleys, thus forming a topographic divide between the Adriatic and Black Sea basins, which gives the test greater regional significance. In early September 2019, three kilograms of uranine were injected into a perched lake in a remote part of the system, approximately 900 metres below the plateau and 100 metres above the low water table. All known springs in the valleys on either side of the mountain were monitored by manual or instrumental sampling and a field fluorometer. Due to the unexpectedly dry season, no tracer was detected at any site for two months until a heavy rainfall event in early November. Subsequently, about 60-65 % of the tracer mass appeared within 60 hours in the Tolminka River. No tracer was detected at other sites, either becauseit was not present or because it was highly diluted. The study suggests that the lake containing the tracer is bypassed by the vadose flow and that the tracer was only mobilised during large events when the lake became part of the epihreatic flow. The linear peak flow velocity from the injection site to the Tolminka Spring was only about 1.7 m/h. However, assuming that the tracer was only mobilised by the large rain event, the velocity would be 70 m/h. The study highlights the challenges and pitfalls of water tracing in alpine karst systems and suggests ways to avoid them.

Improving understanding of groundwater flow in an alpine karst system by reconstructing its geologic history using conduit network model ensembles

Abstract. Reconstructing the geologic history of a karst area can advance understanding of the system's present-day hydrogeologic functioning and help predict the location of unexplored conduits. This study tests competing hypotheses describing past conditions controlling cave formation in an alpine karst catchment, by comparing an ensemble of modeled networks to the observed network map. The catchment, the Gottesacker karst system (Germany and Austria), is drained by three major springs and a paleo-spring and includes the partially explored Hölloch cave, which consists of an active section whose formation is well-understood and an inactive section whose formation is the subject of debate. Two hypotheses for the formation of the inactive section are the following: (1) glaciation obscured the three present-day springs, leaving only the paleo-spring, or (2) the lowest of the three major springs (Sägebach) is comparatively young, so its subcatchment previously drained to the paleo-spring. These hypotheses were tested using the pyKasso Python library (built on anisotropic fast-marching methods) to generate two ensembles of networks, one representing each scenario. Each ensemble was then compared to the known cave map. The simulated networks generated under hypothesis 2 match the observed cave map more closely than those generated under hypothesis 1. This supports the conclusion that the Sägebach spring is young, and it suggests that the cave likely continues southwards. Finally, this study extends the applicability of model ensemble methods from situations where the geologic setting is known but the network is unknown to situations where the network is known but the geologic evolution is not.

Aeolian inputs and dolostone dissolution involved in soil formation in Alpine karst landscapes (Corna Bianca, Italian Alps)

Very different soil types occur across a few tens of meters on dolostone, in a never glaciated karst landscape in the Lombard Pre-Alps (Selvino, Italy). This substrate is locally enriched in well-crystallized sand-grained quartz. The quartz content is responsible for the localized genesis of Podzols. Other soil types observed in the area include strongly rubified (Terra Rossa) horizons (paleosols), preserved in the most protected dolostone cracks. Non-rubified Luvisols and Cambisols were observed in karst dry valleys and dolines while Rendzic Leptosols/Phaeozems were common on the steepest slopes. Such a variety of soils was explained assuming different parent materials (dolostone, silica-rich dolostone, with different amounts of aeolian inputs, ascertained using textural properties, mineralogy, micromorphology, total element composition, mass balance calculations, rare earth elements, and stable elements. Many soils were highly polycyclic, with different layers associated to different parent materials and characterized by different pedogenic processes evidencing different ages. We were thus able to distinguish the horizons mainly developed from the dissolution of the dolostone from those formed in Pleistocene loess. The geochemistry of all surface soil horizons, including Podzols and Rendzic Leptosols/Phaeozems, apparently formed from pure or quartz-rich dolostone dissolution, has been influenced by recent aeolian additions (likely Saharan dust), with deeply modified effects according to different pedogenetic processes acting locally. Saharan dust, in fact, significantly increased metal and rare earth elements contents compared to the substrate, also in the youngest and least weathered soil types.

Hydrochemical Characteristics of the Large Karst Springs in the Catchment of Mesta River (Bulgaria)

Mesta river is a cross-border river with Greece, whose catchment area includes parts of Rila, Pirin and the Rhodope Mountains. The karst is developed in pre-Paleozoic marble that forms isolated outcrops with different basins. The largest of them forms two alpine karst basins located in Northern Pirin and Slavyanka Mountains. In the Rhodopes region, the marble outcrops are fragmented and embedded in non-karstic rocks. They are drained from springs with relatively low flow rates. The largest karst springs drain different hydrodynamic zones of the karst massif of Northern Pirin. From the karst basin of Slavyanka Mountain, only the large subthermal spring near the village of Musomishta falls into the catchment area of the Mesta river. The data on the chemical composition of some of the karst springs are summarized, paying attention to the largest of them. An analysis of the changes in the chemical composition and comparison is made. The changes of the hydrochemical parameters of the springs draining different vertical zones in Northern Pirin are compared. The obtained results prove that the main role in the formation of the chemical composition of karst waters has natural factors - interaction between water and marbles. It was found that the waters are unsaturated to carbonate minerals, which proves that even now there are active processes of karst formation.

High-resolution multi-parameter monitoring of microbial water quality and particles at two alpine karst springs as a basis for an early-warning system

Karst aquifers are important resources for drinking water supply and are very vulnerable to contamination. Microbial concentrations at karst springs, in particular, often vary quickly over a short period of time. In this study, the response of microbial water quality and particle-size distribution of two alpine karst springs to rainfall events was investigated to test and validate parameters that can be used as early-warning systems for fecal contamination. At both investigated karst springs, total organic carbon, particle-size distribution (especially small particle fractions), and particle load show a good correlation to the fecal indicator bacteria E. coli and can therefore be used as a real-time indicator of fecal contamination at the investigated springs. In addition to conventional bacterial determination methods, the β-D-glucuronidase activity, which can be measured in near real-time, was used as a novel indicator parameter for fecal contamination. At the event scale, the β-D-glucuronidase (GLUC) activity shows a good correlation to E. coli and can be used as an additional real-time indicator of fecal contamination. For the studied springs, when they show two peaks in turbidity and small particles, these two parameters are suitable for an early warning system because the bacterial contamination occurs during the secondary peak of these parameters. These results highlight the vulnerability of karst aquifers and demonstrate the applicability of advanced measurement techniques in detecting fecal contamination in real-time, which is especially important given the time-consuming nature of conventional bacterial detection methods.

Multi-year evolution of 75 snow and ice deposits in Schachtdolines and Shafts of recently deglaciated karst terrain: Observations from Mount Canin-Kanin, Julian Alps, Europe

Dolines and shafts are common geomorphologies in karst landscapes. When affected also by glacial processes, they are considered an example of glaciokarst landform. In high alpine karst, their typical vertical shape lends itself particularly well to the formation of snow and ice accumulations that can survive the ablation season, potentially becoming perennial. This study analyses the 12-years changes of several permanent snow-ice deposits in Schachtdolines and Shafts (SIDS) and compares them with the mass balance of surrounding external ice bodies. According to recent studies and in contrast with the rest of the Alpine glaciers, ice masses of the area are resilient to climate change and show a slightly positive mass balance in the last decade. The presence of a long term mass balance monitoring program provides an excellent basis for a robust comparison with 75 selected SIDS, all located in the Mount Canin-Kanin massif, in the Julian Alps. Six airborne LiDAR surveys performed from 2006 to 2018 at the end of the ablation season allowed accurate calculation of thickness and volume changes. The measured differences show a general increase in thickness of the SIDS with more than one meter gain, although multi-annual observed variability between the sites is high. The observed 2006–2018 positive mass balance is in agreement with the external cryosphere although SIDS appear even more resilient than the external ice patches to both interannual weather variability and climate forcing.

Impact of Cold Stress on Leaf Structure, Photosynthesis, and Metabolites in Camellia weiningensis and C. oleifera Seedlings

Camellia weiningensis Y. K. Li. sp. nov. (CW) is an endemic oil-tea species in Guizhou province, distributed in the alpine karst area, which exhibits cold resistance and better economic characters than C. oleifera (CO). The mechanism of cold response in CW seedlings has not been studied in depth. Herein, we performed anatomical, physiological, and metabolic analyses to assess the impact of cold stress on leaf structure, photosynthesis, and metabolites in CW and CO seedlings. Anatomical analysis of leaves showed CW seedlings had greater leaf and palisade thicknesses, tissue structure tightness, and palisade-spongy tissue ratio to enhance chilling stress (4 °C) tolerance, but freezing stress (−4 °C) caused loosening of the leaf tissue structure in both CW and CO seedlings. Photosynthetic analysis showed a reduction in the chlorophyll (Chl) fluorescence (Fv/Fm) and photosynthetic parameters under freezing stress in both CW and CO seedlings. Cold stress increased the abscisic acid (ABA) contents in both the Camellia species, and CW exhibited the highest ABA content under −4 °C treatment. Additionally, the indole-3-acetic acid (IAA) content was also increased in CW in response to cold stress. An obviously distinct metabolite composition was observed for CW and CO under different temperatures, and significantly changed metabolites (SCMs) were enriched under freezing stress. Prenol lipids, organooxygen compounds, and fatty acyls were the main metabolites in the two Camellia species in response to cold stress. The top key SCMs, such as medicoside G, cynarasaponin F, yuccoside C, and methionyl-proline were downregulated under freezing stress in both CW and CO. The contents of some key metabolites associated with sugar metabolism, such as UDP-glucose, UDP-D-apiose, and fructose 6-phosphate, were higher in CW than in CO, which may contribute to enhancing the cold resistance in CW. Our findings are helpful in explaining how CW adapt to alpine karst cold environments, and will provide a reference for cold tolerance improvement and application of stress-resistant breeding of Camellia in alpine and cold areas.

The stochastic simulation of karst conduit network structure using anisotropic fast marching, and its application to a geologically complex alpine karst system

Anisotropic fast-marching algorithms are computationally efficient tools for generating realistic maps of karst conduit networks, constrained by both the spatial extent and the orientation of karstifiable geologic units. Existing models to generate conduit network maps are limited either by high computational requirements (for chemistry-based models) or by their inability to incorporate the effects of elevation and orientation gradients (for isotropic fast-marching models). The new anisotropic fast-marching approach described here provides a significant improvement, though it imitates rather than reproduces actual speleogenetic processes. It can rapidly generate a stochastic ensemble of plausible networks from basic geologic information, which can also be used as input to karst-appropriate flow models. This paper introduces an open-source, easy-to-use implementation through the Python package pyKasso, then describes its application to a well-mapped geologically complex long-term study site: the Gottesacker alpine karst system (Germany/Austria). Groundwater flow in this system is exceptionally well understood from speleological investigations and tracer tests. Conduit formation primarily occurs at the base of the karst aquifer, following plunging synclines. Although previous attempts to reproduce the conduit network at this site yielded implausible network maps, pyKasso quickly generated networks faithful to the known conduit system. However, the model was only able to generate these realistic networks when the inlet-outlet connections of the system were correctly assigned, highlighting the importance of pairing modeling efforts with field tracer tests. Therefore, a model ensemble method is also presented, to optimize field efforts by identifying the most informative tracer tests to perform.

Relations between Benthic Diatom Community and Characteristics of Karst Ponds in the Alpine Region of Slovenia

The aim of this research was to investigate the structure of the benthic diatom community and its relations to selected environmental parameters. We collected samples in 16 karst ponds in the alpine region of Slovenia, where the Alpine karst is found. Since the predominating substrate in these ponds was clay, the epipelic community was analyzed. Hydromorphological characteristics, and physical and chemical conditions were also measured at each site. We found 105 species of diatoms, which belonged to 32 genera. The most frequent taxa were Gomphonema parvulum (Kützing) Kützing, Navicula cryptocephala Kützing, Sellaphora pupula (Kützing) Mereschkowsky (species group) and Achnanthidium pyrenaicum (Hustedt) Kobayasi. The pond with the lowest diversity was found at the highest altitude, while, on the other hand, the most species-rich pond was found at the lowest altitude. Regarding the ecological types, the most common were motile species. We confirmed a positive correlation between the number of diatom species and the saturation of water with oxygen, while correlation between species richness and NH4-N was negative. The content of NO3-N and NH4-N explained almost 20% of the total variability of diatom community. Unlike our expectations, we calculated a negative correlation between the diversity of macroinvertebrates and diatoms, which is probably a consequence of different responses to environmental conditions.

Unusual catchment runoff in a high alpine karst environment influenced by a complex geological setting (Northern Calcareous Alps, Tyrol, Austria)

Garber Schlag (Q-GS) is one of the major springs of the Karwendel Mountains, Tyrol, Austria. This spring has a unique runoff pattern that is mainly controlled by the tectonic setting. The main aquifer is a moderately karstified and jointed limestone of the Wetterstein Formation that is underlain by nonkarstified limestone of the Reifling Formation, which acts as an aquitard. The aquifer and aquitard of the catchment of spring Q-GS form a large anticline that is bound by a major fault (aquitard) to the north. Discharge of this spring shows strong seasonal variations with three recharge origins, based on δ18O and electrical conductivity values. A clear seasonal trend is observed, caused by the continuously changing portions of water derived from snowmelt, rainfall and groundwater. At the onset of the snowmelt period in May, the discharge is composed mainly of groundwater. During the maximum snowmelt period, the water is dominantly composed of water derived from snowmelt and subordinately from rainfall. During July and August, water derived from snowmelt continuously decreases and water derived from rainfall increases. During September and October, the water released at the spring is mainly derived from groundwater and subordinately from rainfall. The distinct discharge plateau from August to December and the following recession until March is likely related to the large regional groundwater body in the fissured and moderately karstified aquifer of the Wetterstein Formation and the tectonic structures (anticline, major fault). Only a small portion of the water released at spring Q-GS is derived from permafrost.

Applicability of Susceptibility Model for Rock and Loess Earthquake Landslides in the Eastern Tibetan Plateau

It is crucial to explore a suitable landslide susceptibility model with an excellent prediction capability for rapid evaluation and disaster relief in seismic regions with different lithological features. In this study, we selected two typical seismic events, the Jiuzhaigou and Minxian earthquakes, which occurred in the Alpine karst and loess regions, respectively. Eight influencing factors and five models were chosen to calculate the susceptibility of landslide, including the information (I) model, certainty factor (CF) model, logistic regression (LR) model, I + LR coupling model, and CF + LR coupling model. Then, the accuracy and the landslide susceptibility distribution of these models were assessed by the area under curve (AUC) and distribution criteria. Finally, the model with high accuracy and good applicability for the rock landslide or loess landslide regions was optimized. Our results showed that the accuracy of the coupling model is higher than that of the single models. Except for the LR model, the landslide susceptibility distribution for the above-mentioned models is consistent with universal cognition. The coupling models are generally better than their single models. Among them, the I + LR model can obtain the best comprehensive results for assessing the distribution and accuracy of both rock and loess landslide susceptibility, which is helpful for disaster relief and policy-making, and it can also provide useful scientific data for post-seismic reconstruction and restoration.

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

AbstractLarge 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 CO2 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.

Geochemical insights into groundwater movement in alpine karst, Bear River Range, Utah, USA

Alpine karst aquifers control the availability and longevity of some water resources, but are not well understood. A conceptual model of the alpine karst aquifer system in the Bear River Range of northern Utah (USA) has been developed by geochemical analysis (major ions, δ18O, δ2H and δ13C values) of seasonal water samples from seven perennial springs, and residence-time assessment (3H and CFCs) of two low- and two high-discharge springs. All spring data can be explained by reaction paths dominated by the dissolution of calcian dolomite. The δ13C values align well with reaction paths for open-system dissolution. Saturation indices and low Ca:Mg molar ratios indicate that incongruent dissolution exerts a strong control on water–rock interactions, complicating interpretation of natural solute tracers. Values of δ18O and δ2H in springs follow the Utah meteoric water line. Snow δ18O values correlate with elevation, but not with increasing rainout distance, providing qualitative estimates of recharge elevation that generally align with previous dye-traces to five of the seven springs. Concentrations of 3H and CFCs likely are best described by binary mixing of subannual recharge with 60–65-year-old groundwater, suggesting that the alpine karst aquifer system in the Bear River Range is best represented by a double-porosity model. Subannual recharge documented by dye traces implies that caverns are the primary flowpaths to the springs, but the presence of decadal-age water may indicate that lower permeability flowpaths dominate during baseflow. No evidence was found for a longer-residing flow component, suggesting high sensitivity to future climate variability.

A model ensemble generator to explore structural uncertainty in karst systems with unmapped conduits

Karst aquifers are characterized by high-conductivity conduits embedded in a low-conductivity fractured matrix, resulting in extreme heterogeneity and variable groundwater flow behavior. The conduit network controls groundwater flow, but is often unmapped, making it difficult to apply numerical models to predict system behavior. This paper presents a multi-model ensemble method to represent structural and conceptual uncertainty inherent in simulation of systems with limited spatial information, and to guide data collection. The study tests the new method by applying it to a well-mapped, geologically complex long-term study site: the Gottesacker alpine karst system (Austria/Germany). The ensemble generation process, linking existing tools, consists of three steps: creating 3D geologic models using GemPy (a Python package), generating multiple conduit networks constrained by the geology using the Stochastic Karst Simulator (a MATLAB script), and, finally, running multiple flow simulations through each network using the Storm Water Management Model (C-based software) to reject nonbehavioral models based on the fit of the simulated spring discharge to the observed discharge. This approach captures a diversity of plausible system configurations and behaviors using minimal initial data. The ensemble can then be used to explore the importance of hydraulic flow parameters, and to guide additional data collection. For the ensemble generated in this study, the network structure was more determinant of flow behavior than the hydraulic parameters, but multiple different structures yielded similar fits to the observed flow behavior. This suggests that while modeling multiple network structures is important, additional types of data are needed to discriminate between networks.

Geological and geomorphic exploration and formation mechanism of the Shilin Geopark in Zecha, Gansu Province

Zecha Shilin Geopark is located in Luqu County, Gansu Province, where the landform completely preserves the characteristics of alpine karst, with typicality and rarity. In order to reveal its natural scientific value, promote the development of natural resources and landscape protection, through field geological survey, taking Luqu County regional geology as the background, the karst landform types of Zecha Shilin Geopark are introduced, which are mainly plateau stone forest, a line of sky as the highlight, and small karst cave as the auxiliary. Then, the formation mechanism of karst landform in the park is analyzed from four aspects of geological structure, lithology, hydrological climate, and biological factors. It is of certain significance to improve the scientific research of karst series geological relics in high-altitude areas. At the same time, the ecological transformation of alpine karst has a good demonstration effect on the control of karst desertification in South China.

Effects of an extreme flood event on an alpine karst system

The effects of an extreme storm (501 mm of rainfall in less than five days) were monitored in an alpine show cave to assess the characteristics of its aquifer unsaturated and saturated zones. The selected cave, Bossea (SW Piedmont, Italy), is an excellent test site for hydrological investigation of karst systems, because it hosts an underground karst laboratory since the late 70 s. The investigation was carried out by means of an integrated approach that combined measurements of flow discharge, physico-chemical parameters, major chemical components, and trace elements (metals, Rb, Ba, Sr and REE). The hydrology and hydrochemistry of the main underground river and two of its secondary tributaries (a drip site and a small secondary inflow) draining the unsaturated zone were monitored during the November 2016 flood, an exceptional hydrological event (estimated recurrence time of 200 years) that caused severe damages in the whole southwestern Piedmont region. The results of the 2016 monitoring were compared with those of another less extreme flood occurred in 2011 (recurrence time of 20 years). The karst system showed an impulsive response to flooding at the catchment and at the individual inflow scales, but each site exhibited a characteristic response due to the complex geological and structural setting of the cave system. The development of this karst system is, in fact, related to the lateral and vertical juxtaposition of rocks with different lithologies and mechanical properties, which form different hydrogeologic compartments. Coupling hydrograph, chemograph analysis and REE normalized patterns permitted to assess which compartments were activated during each phase of the floods. In particular, it was possible to recognize a progressive increasing contribution of the non-carbonate lithologies to flow during the peak of the flood for both the two unsaturated inflows and the main underground river, suggesting the activation of larger portions of the aquifer. However, the response of each individual unsaturated inflow is more influenced by its recharge system architecture rather than the magnitude of the meteorological event. Similar complex geological and lithological karst systems are typical of many areas characterized by orogenetic processes (Alps, Rocky Mountains, Caledonides, etc.), but despite them being relatively widespread, the hydrodynamics of these aquifers is still poorly understood. This study affirms that only long-term and well-integrated monitoring and sampling can help unravel the behavior of such complex karst systems.

Snowmelt as a determinant factor in the hydrogeological behaviour of high mountain karst aquifers: The Garcés karst system, Central Pyrenees (Spain)

Time series of environmental tracers (groundwater stable isotope composition, electrical conductivity and temperature) and concentration breakthrough curves of artificial tracers (uranine, eosine, amino-G and naphtionate) have been analyzed to characterize fast preferential and slow matrix in-transit recharge flows in the Paleocene-Eocene limestone aquifer of the Ordesa and Monte Perdido National Park, an alpine karst system drained by a water table cave, a rare hydrological feature in high mountain karst systems with similar characteristics. Snowmelt favors the areal recharge of the system. This process is reflected in the large proportion of groundwater flowing through the connected porosity structure of the karst aquifer, which amounts the 75% of the total system water discharge. From the perspective of water resources recovery, the water capacity of the fissured-porous zone (matrix) represents 99% of the total karst system storage. The volume associated to the karst conduits is very small. The estimated mean travel times are 9 days for conduits and 475 days for connected porosity. These short travel times reveal high vulnerability of the karst system to pollutants in broad sense and a great impact of climate change on the associated water resources.

Complementary use of dating and hydrochemical tools to assess mixing processes involving centenarian groundwater in a geologically complex alpine karst aquifer

AbstractEnvironmental dating tracers (3H, 3He, 4He, CFC‐12, CFC‐11, and SF6) and the natural spring response (hydrochemistry, water temperature, and hydrodynamics) were jointly used to assess mixing processes and to characterize groundwater flow in a relatively small carbonate aquifer with complex geology in southern Spain. Results evidence a marked karst behaviour of some temporary outlets, with sharp and rapid responses to precipitation events, while some perennial springs show buffer and delayed variations with respect to recharge periods. The general geochemical evolution shows a pattern, from higher to lower altitudes, in which mineralization and the Mg/Ca ratio rise, evidencing longer water–rock interaction. The large SF6 concentrations in groundwater suggest terrigenic production, whereas CFC‐11 values are affected by sorption or degradation. The groundwater age in the perennial springs—as deduced from CFC‐12 and 3H/3He—points to mean residence times of several decades, although the large amount of radiogenic 4He in samples indicate a contribution of old groundwater (free of 3H and CFC‐12). Lumped parameter models and shape‐free models were created based on 3H, tritiogenic 3He, CFC‐12, and radiogenic 4He data in order to interpret the age distribution of the samples. Results evidence the existence of two mixing components, with an old fraction ranging between 160 and 220 years in age. The correlation of physicochemical parameters with some dating parameters, derived from the mixing models, serves to explain the hydrogeochemical processes occurring within the system. Altogether, long residence times are shown to be possible in small alpine systems with a clearly karst behaviour if the geological setting features highly tectonized media including units with diverse hydrogeological characteristics. These findings highlight the importance of applying different approaches, including groundwater dating techniques, when studying such groundwater flow regimes.