Author Correction: Global karst springs hydrograph dataset for research and management of the world's fastest-flowing groundwater.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Global karst springs hydrograph dataset for research and management of the world's fastest-flowing groundwater.

Karst aquifers provide drinking water for 10% of the world's population, support agriculture, groundwater-dependent activities, and ecosystems. These aquifers are characterised by complex groundwater-flow systems, hence, they are extremely vulnerable and protecting them requires an in-depth understanding of the systems. Poor data accessibility has limited advances in karst research and realistic representation of karst processes in large-scale hydrological studies. In this study, we present World Karst Spring hydrograph (WoKaS) database, a community-wide effort to improve data accessibility. WoKaS is the first global karst springs discharge database with over 400 spring observations collected from articles, hydrological databases and researchers. The dataset's coverage compares to the global distribution of carbonate rocks with some bias towards the latitudes of more developed countries. WoKaS database will ensure easy access to a large-sample of good quality datasets suitable for a wide range of applications: comparative studies, trend analysis and model evaluation. This database will largely contribute to research advancement in karst hydrology, supports karst groundwater management, and promotes international and interdisciplinary collaborations.

Combining monitoring and modelling tools as a basis for city-scale concepts for a sustainable thermal management of urban groundwater resources.

Increasing anthropogenic impacts lead to elevated temperatures of the shallow subsurface, including the unsaturated and groundwater saturated zone, in many urban areas in comparison to unaffected natural thermal states. The "current thermal state" of four groundwater bodies in the urban area of Basel-City, Switzerland, was investigated by means of high-resolution multilevel temperature wells and numerical 3D groundwater flow and heat transport models. The calibrated and validated numerical groundwater flow and heat transport models allow evaluating and comparing groundwater and heat fluxes for the investigated groundwater bodies and defined cross-sections for differing urban districts, e.g. residential and industrial areas under development. We present the overall and the specific advective heat fluxes within two urban districts, which will be restructured in the near future. The management of groundwater resources in urban areas plays an important role not only for groundwater quantity but also for its quality, i.e. thermal subsurface and groundwater regimes. We demonstrate how monitoring and modelling tools can be the basis for a sustainable management of complex urban groundwater resources. Furthermore, we argue that such tools should be integrated in the thermal management of urban groundwater bodies. Such tools also allow integrating the potentially available energy of shallow subsurface resources into energetic management strategies on the urban scale.

Process-based monitoring and modeling of Karst springs - Linking intrinsic to specific vulnerability.

The presented work illustrates to what extent field investigations as well as monitoring and modeling approaches are necessary to understand the high discharge dynamics and vulnerability of Karst springs. In complex settings the application of 3D geological models is essential for evaluating the vulnerability of Karst systems. They allow deriving information on catchment characteristics, as the geometry of aquifers and aquitards as well as their displacements along faults. A series of Karst springs in northwestern Switzerland were compared and Karst system dynamics with respect to qualitative and quantitative issues were evaluated. The main objective of the studies was to combine information of catchment characteristics and data from novel monitoring systems (physicochemical and microbiological parameters) to assess the intrinsic vulnerability of Karst springs to microbiological contamination with simulated spring discharges derived from numerical modeling (linear storage models). The numerically derived relation of fast and slow groundwater flow components enabled us to relate different sources of groundwater recharge and to characterize the dynamics of the Karst springs. Our study illustrates that comparably simple model-setups were able to reproduce the overall dynamic intrinsic vulnerability of several Karst systems and that one of the most important processes involved was the temporal variation of groundwater recharge (precipitation, evapotranspiration and snow melt). Furthermore, we make a first attempt on how to link intrinsic to specific vulnerability of Karst springs, which involves activities within the catchment area as human impacts from agriculture and settlements. Likewise, by a more detailed representation of system dynamics the influence of surface water, which is impacted by release events from storm sewers, infiltrating into the Karst system, could be considered. Overall, we demonstrate that our approach can be the basis for a more flexible and differentiated management and monitoring of raw-water quality of Karst springs.

Spatiotemporal scales of river-groundwater interaction - The role of local interaction processes and regional groundwater regimes.

Drinking water production in the vicinity of rivers not only requires the consideration of different spatiotemporal scales and settings of river-groundwater interaction processes, but also of local and regional scale groundwater regimes. Selected case studies in combination with field-experiments and the setup of high-resolution groundwater flow models enabled the investigation of the spatiotemporal development of microbial (classical fecal indicator bacteria and total cell counts) and selected organic micropollutants in riverine and regional groundwater for different hydrological settings, including low and high flow conditions. Proxy indicators suitable as surrogates for the diverse contaminations in alluvial aquifers with different settings could be identified. Based on the study results, the basic elements for both groundwater management and river restoration concepts are derived, which include the: (1) compilation and evaluation of the "current state" concerning hydrogeology, microbiology and contamination by organic micropollutants, (2) definition of field-experiments to qualitatively assess variability related to the "current state", and (3) quantitative assessment of groundwater regimes, including variability of groundwater components and inflow areas, by application of high-resolution groundwater flow models. The validity and transferability of the concept and inferred controls (specifically drivers and controls of river-groundwater interaction) are tested by evaluations derived from hydraulic relationships to river sections with comparable settings and regional groundwater flow regimes in general. The results of our investigations illustrate the influence of dynamic hydrologic boundary conditions on river-groundwater interaction and of regional scale groundwater flow regimes on the water composition of riverine groundwater systems. It is demonstrated how to identify river sections and their variations with intensified river-groundwater exchange processes and how to quantify the transient character of the different groundwater components that constitute the raw water quality of drinking water wells near rivers.

Comprehensive micropollutant screening using LC-HRMS/MS at three riverbank filtration sites to assess natural attenuation and potential implications for human health.

Riverbank filtration (RBF) is used worldwide to produce high quality drinking water. With river water often contaminated by micropollutants (MPs) from various sources, this study addresses the occurrence and fate of such MPs at three different RBF sites with oxic alluvial sediments and short travel times to the drinking water well down to hours. A broad range of MPs with various physico-chemical properties were analysed with detection limits in the low ng L-1 range using solid phase extraction followed by liquid chromatography coupled to tandem high resolution mass spectrometry. Out of the 526 MPs targeted, a total of 123 different MPs were detected above the limit of quantification at the three different RBF sites. Of the 75-96 MPs detected in each river, 43-59% were attenuated during RBF. The remaining total concentrations of the MPs in the raw drinking water accounted to 0.6-1.6 µgL-1 with only a few compounds exceeding 0.1 µgL-1, an often used threshold value. The attenuation was most pronounced in the first meters of infiltration with a full elimination of 17 compounds at all three sites. However, a mixing with groundwater related to regional groundwater flow complicated the characterisation of natural attenuation potentials along the transects. Additional non-target screening at one site revealed similar trends for further non-target components. Overall, a risk assessment of the target and estimated non-target compound concentrations finally indicated during the sampling period no health risk of the drinking water according to current guidelines. Our results demonstrate that monitoring of contamination sources within a catchment and the affected water quality remains important in such vulnerable systems with partially short residence times.

Online analysis: Deeper insights into water quality dynamics in spring water.

We have studied the dynamics of water quality in three karst springs taking advantage of new technological developments that enable high-resolution measurements of bacterial load (total cell concentration: TCC) as well as online measurements of abiotic parameters. We developed a novel data analysis approach, using self-organizing maps and non-linear projection methods, to approximate the TCC dynamics using the multivariate data sets of abiotic parameter time-series, thus providing a method that could be implemented in an online water quality management system for water suppliers. The (TCC) data, obtained over several months, provided a good basis to study the microbiological dynamics in detail. Alongside the TCC measurements, online abiotic parameter time-series, including spring discharge, turbidity, spectral absorption coefficient at 254nm (SAC254) and electrical conductivity, were obtained. High-density sampling over an extended period of time, i.e. every 45min for 3months, allowed a detailed analysis of the dynamics in karst spring water quality. Substantial increases in both the TCC and the abiotic parameters followed precipitation events in the catchment area. Differences between the parameter fluctuations were only apparent when analyzed at a high temporal scale. Spring discharge was always the first to react to precipitation events in the catchment area. Lag times between the onset of precipitation and a change in discharge varied between 0.2 and 6.7h, depending on the spring and event. TCC mostly reacted second or approximately concurrent with turbidity and SAC254, whereby the fastest observed reaction in the TCC time series occurred after 2.3h. The methodological approach described here enables a better understanding of bacterial dynamics in karst springs, which can be used to estimate risks and management options to avoid contamination of the drinking water.

The thermal impact of subsurface building structures on urban groundwater resources - A paradigmatic example.

Shallow subsurface thermal regimes in urban areas are increasingly impacted by anthropogenic activities, which include infrastructure development like underground traffic lines as well as industrial and residential subsurface buildings. In combination with the progressive use of shallow geothermal energy systems, this results in the so-called subsurface urban heat island effect. This article emphasizes the importance of considering the thermal impact of subsurface structures, which commonly is underestimated due to missing information and of reliable subsurface temperature data. Based on synthetic heat-transport models different settings of the urban environment were investigated, including: (1) hydraulic gradients and conductivities, which result in different groundwater flow velocities; (2) aquifer properties like groundwater thickness to aquitard and depth to water table; and (3) constructional features, such as building depths and thermal properties of building structures. Our results demonstrate that with rising groundwater flow velocities, the heat-load from building structures increase, whereas down-gradient groundwater temperatures decrease. Thermal impacts on subsurface resources therefore have to be related to the permeability of aquifers and hydraulic boundary conditions. In regard to the urban settings of Basel, Switzerland, flow velocities of around 1 md-1 delineate a marker where either down-gradient temperature deviations or heat-loads into the subsurface are more relevant. Furthermore, no direct thermal influence on groundwater resources should be expected for aquifers with groundwater thicknesses larger 10m and when the distance of the building structure to the groundwater table is higher than around 10m. We demonstrate that measuring temperature changes down-gradient of subsurface structures is insufficient overall to assess thermal impacts, particularly in urban areas. Moreover, in areas which are densely urbanized, and where groundwater flow velocities are low, appropriate measures for assessing thermal impacts should specifically include a quantification of heat-loads into the subsurface which result in a more diffuse thermal contamination of urban groundwater resources.

Online flow cytometry reveals microbial dynamics influenced by concurrent natural and operational events in groundwater used for drinking water treatment.

Detailed measurements of physical, chemical and biological dynamics in groundwater are key to understanding the important processes in place and their influence on water quality - particularly when used for drinking water. Measuring temporal bacterial dynamics at high frequency is challenging due to the limitations in automation of sampling and detection of the conventional, cultivation-based microbial methods. In this study, fully automated online flow cytometry was applied in a groundwater system for the first time in order to monitor microbial dynamics in a groundwater extraction well. Measurements of bacterial concentrations every 15 minutes during 14 days revealed both aperiodic and periodic dynamics that could not be detected previously, resulting in total cell concentration (TCC) fluctuations between 120 and 280 cells µL-1. The aperiodic dynamic was linked to river water contamination following precipitation events, while the (diurnal) periodic dynamic was attributed to changes in hydrological conditions as a consequence of intermittent groundwater extraction. Based on the high number of measurements, the two patterns could be disentangled and quantified separately. This study i) increases the understanding of system performance, ii) helps to optimize monitoring strategies, and iii) opens the possibility for more sophisticated (quantitative) microbial risk assessment of drinking water treatment systems.

On the effects of subsurface parameters on evaporite dissolution (Switzerland).

Uncontrolled subsurface evaporite dissolution could lead to hazards such as land subsidence. Observed subsidences in a study area of Northwestern Switzerland were mainly due to subsurface dissolution (subrosion) of evaporites such as halite and gypsum. A set of 2D density driven flow simulations were evaluated along 1000 m long and 150 m deep 2D cross sections within the study area that is characterized by tectonic horst and graben structures. The simulations were conducted to study the effect of the different subsurface parameters that could affect the dissolution process. The heterogeneity of normal faults and its impact on the dissolution of evaporites is studied by considering several permeable faults that include non-permeable areas. The mixed finite element method (MFE) is used to solve the flow equation, coupled with the multipoint flux approximation (MPFA) and the discontinuous Galerkin method (DG) to solve the diffusion and the advection parts of the transport equation. Results show that the number of faults above the lower aquifer that contains the salt layer is considered as the most important factor that affects the dissolution compared to the other investigated parameters of thickness of the zone above the halite formation, a dynamic conductivity of the lower aquifer, and varying boundary conditions in the upper aquifer.

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.

Enhanced vulnerability assessment in karst areas by combining mapping with modeling approaches.

The objective of this work is to facilitate a sustainable regional planning of water resources in karst areas by providing a conceptual framework for an integrative vulnerability assessment. A combined mapping and modeling approach is proposed, taking into account both spatial and temporal aspects of karst groundwater vulnerability. The conceptual framework comprises the delineation of recharge areas, vulnerability mapping, numerical flow and transport modeling and the integration of information into a combined vulnerability map and time series. The approach is illustrated at a field site in northwest Switzerland (Gempen plateau). The results show that the combination of vulnerability mapping and numerical modeling allows the vulnerability distribution, both in the recharge and discharge areas, to be identified, and at the same time, the time dependence of karst groundwater vulnerability to be assessed. The combined vulnerability map and time series provide a quantitative basis for drinking water management and for regional planning.