1. Introduction
Risk based approach to drinking water safety is taken up world-wide to complement endpoint monitoring (World Health Organization 2022). Monitoring at the tap is still seen as an indispensable element of drinking water control but shifts from being the primary means of safeguarding drinking water quality towards becoming a tool for validating water safety plans (or similar frameworks). The recent recast of the European Union (EU) 2020/2184 Drinking Water Directive (DWD) introduced for the first time the obligation of risk assessment for the entire water supply chain [
1]. As drinking water quality is strongly influenced by the water source, better understanding of source water quality and the associated hazards is a cornerstone of the risk-based approach [
2]. Hazard assessment of the catchments of drinking water abstraction points is also a requirement under the DWD [
1].
Drinking water quality regulations addressing parameters derived from the source water usually focus on potential anthropogenic pollutants and a few geological compounds of known health impact, such as arsenic. Monitoring requirements under the former EU regulation in the European Union followed the same approach, focusing mainly on heavy metals associated primarily with human pollution (e.g. mercury, lead, cadmium, chromium, nickel) [
3]. The number of geogenic pollutants for which substantial monitoring data is available in the EU is limited, comprising generally of arsenic, selenium, boron, and fluoride, though national legislations may include additional parameters. The recast DWD extended monitoring requirements (starting from 2026) to new components, including metals of health relevance. Uranium was added with a parametric value of 30 µg/L following the recommendation of the World Health Organisation (WHO) [
4]. Calcium, potassium and magnesium were also included in the list of parameters to be monitored, but only for consumer information, and no parametric values have been assigned [
1].
The newly introduced natural elements have a confirmed or presumptive impact on human health. Uranium is a nephrotoxic pollutant of primarily geological origin [
4,
5]. Its concentration in natural waters is usually low (<1.0 µg/L), but depending on the geological framework (e.g. granite rich in uranium) and the mobilizing capacity of groundwater flow systems high concentrations can also occur in certain areas (up to 57 µg/L in Mongolia, 85 µg/L in Italy, 288 µg/L in China and 750 µg/L in Norway) [
6,
7,
8,
9,
10].
Potassium, calcium and magnesium are essential elements. Potassium concentrations also vary widely in natural waters, but the concentrations present in drinking water are generally not a risk to human health and may be of significance only for sensitive individuals (e.g. those suffering from renal failure) [
4].
The inverse association between magnesium and calcium in drinking water and cardiovascular mortality is well established. The protective effect is clearer for magnesium, while further studies are needed for calcium [
11,
12,
13,
14]. Low magnesium and calcium intake via drinking water have also been associated with high blood pressure, metabolic syndrome, preeclampsia and amyotrophic sclerosis [
15]. For these elements, prior data is scarce, though many jurisdictions, including the Hungarian, require monitoring of total hardness, which primarily consists of magnesium and calcium salts.
Vanadium is not regulated by the DWD, but some EU Member States include it in national legislation, e.g. Italy, with a parametric value of 140 µg/L [
10]. The health effects of vanadium are unclear, but therapeutic and toxic effects on human health are possible (e.g. bone, cardiovascular, neurological, immune system, and body weight effects) [
16,
17]. Its concentration varies in a wide range in surface waters (0.010 µg/L to 68 µg/L), reaching even higher levels in groundwaters, especially in volcanic areas (up to 350 µg/l) [
10,
18].
Lithium is also an unregulated element of health relevance. The concentration of lithium in drinking water in European countries ranged from <1 µg/L to 20–60 µg/L in northern England, Lithuania and Italy [
19,
20,
21], over 100 µg/L in Greece [
22], and more than 1000 µg/L in Austria [
23]. Several studies have found a protective association between lithium in drinking water and suicide mortality [
20,
23,
24,
25,
26,
27].
Molybdenum and selenium are essential elements. The concentration of Mo in drinking water is typically low (<10 µg/L), although near mining areas, this value can be much higher (>200 µg/L). The recommended health-based guideline value of WHO is 70 µg/L. The optimal range of Se intake is narrow. Selenium deficiency has been associated with cardiomyopathy and joint disease, while symptoms of selenosis include brittle hair and hair loss, nail malformation and diarrhea [
28]. The parametric value of Se has been increased from 10 to 20 µg/L during the recast, based on recent scientific evidence [
1,
29]. WHO guideline value is 40 µg/L. Environmental concentrations are generally well below 10 µg/L, except in a few seleniferous areas.
The health impact of barium is unclear. There is no evidence of its carcinogenic and mutagenic properties, but there are contradictory results for cardiovascular toxicity (hypertension) and it has been shown to cause nephropathy in animal studies. The WHO recommends a guideline value of 1300 µg/L in drinking water [
4,
17].
In Hungary, the primary source of drinking water is groundwater (93% of total production volume), including 31% bank filtered water and a smaller proportion of karst water (12%). Surface water accounts for only 7% of the water supply [
30]. There are 1403 public utility water supply systems registered in the National Drinking Water Database, operated by 38 water suppliers. Of these systems, 36 supplies use mainly bank filtration, 116 karst water and 1201 other groundwater sources, while 22 rely on surface water abstraction.
Drinking water quality is regularly monitored by the water suppliers and the local public health authorities. Monitoring does not cover elements of natural origin beyond the EU DWD requirements, although total hardness – comprising primarily of calcium and magnesium – is regulated as an indicator parameter (acceptable range: 50–350 mg/L CaO equivalent). Drinking water in Hungary is typically moderately hard (between 100 and 200 mg CaO/L) or hard (>200 mg CaO/L), with total hardness ranging from <1.0 to 393 mg CaO/L [
30]. For other unregulated elements, only local studies are available. Previous targeted investigations detected uranium in concentrations up to 25 µg/L in some drinking water sources in Hungary, though not exceeding the new parametric value [
31,
32,
33,
34].
The aim of the present study was to carry out a nationwide survey mapping the presence and distribution of (previously) unregulated elements of natural origin in Hungarian drinking water sources. We aimed to identify if the observed distribution of the elements was associated with the water source type or the geographical location of the water supply, reflecting differences in (hydro)geological factors. It was evaluated whether any of these elements poses a potential hazard for drinking water production and requires further attention in risk assessment (water safety plans) on a national or local scale. The survey outcomes also provide representative data for epidemiological studies assessing potential health risks or gains associated with these compounds in drinking water.
2. Materials and Methods
2.1. Study Design
The survey was initiated by the National Center for Public and Pharmacy (NCPHP) as the national regulatory authority. All public utility drinking water suppliers were invited to participate. Participation was voluntary. Water suppliers responding positively received sample containers, instructions for sampling, and a data sheet to be completed during sampling. The sampling data sheet recorded the date and exact location of sampling, the type of raw water (surface, deep groundwater (second or deeper aquifer or depth >50 m), shallow groundwater (from the first aquifer or depth <50 m), karst, bank filtered, other/mixed) and information on water treatment and disinfection.
2.2. Sampling
Water samples were collected between March 2018 and June 2022 by the drinking water supply operators. At least one source water sample was taken for each water supply. Water suppliers were instructed to take a sample representative of the source water composition before water treatment. From water supplies with multiple abstraction points, either a representative subset of the wells or mixed well water representing normal operation were sampled. Samples from water supplies using surface water abstraction were collected after abstraction and mechanical filtration. Samples from water supplies using bank filtration were taken from selected bank filtration wells or their mixed water. Most water supplies were sampled twice in the sampling period. Samples were collected in 50 ml PP centrifuge tubes, and the pH of the samples was reduced below 2 with high-purity nitric acid for preservation (ISO 5667-11:2009). Samples were transported to the laboratory of the NCPHP and stored at room temperature until analysis within 30 days.
2.3. Chemical Analysis
The analytical tests were performed at the NCPHP’s laboratory. Elements were analyzed by ICP-MS ((iCAP RQ, Thermo Fisher Scientific, Germany) (ISO 17294-2:2016). The concentrations of the following elements were measured: boron (B), barium (Ba), beryllium (Be), calcium (Ca), cobalt (Co), potassium (K), lithium (Li), magnesium (Mg), molybdenum (Mo), sodium (Na), selenium (Se), strontium (Sr), titanium (Ti), uranium (U) and vanadium (V).
2.4. Statistical and Geospatial Methods
Statistical evaluation of the results was carried out using the Statistica program (TIBCO Software Inc. (2020). Data Science Workbench, version 14.
http://tibco.com.). Values under the lower limit of quantification (LOQ) were replaced by LOQ/√2 in the statistical calculations [
35].
Hierarchical cluster analysis (HCA) was used to group the samples based on their similar characteristics [
36]. In HCA, each sampling point is initially a separate cluster, and then in iterative steps, the two closest clusters are merged into a single cluster. The number of clusters is thus reduced in each step, eventually leaving a single cluster. Ward’s method and quadratic Euclidean distance were used for merging the clusters. Discriminant analysis (DA) was used to test the appropriateness of the clusters formed by HCA.
Principal component analysis (PCA) was used to reveal the basic characteristics and hidden relationships of the large amount of data and for data reduction [
37].
For multivariate statistical analyses (both PCA and HCA), the results were standardised by z-score transformation to reduce the bias resulting from the often orders of magnitude differences in the concentration of the elements.
Elemental concentrations of the different clusters were compared using the non-parametric Kruskal-Wallis H-test (p-value: 0.05), and the correlation between elements was tested using Pearson's correlation test (p-value: 0.05).
4. Discussion
Drinking water supply in Hungary relies almost entirely on groundwater sources. Surface water abstraction is limited and, in some areas, seasonal. For example, around Lake Balaton, an important holiday destination, it is used in the summer periods to compensate for the temporary increase of the population. Although by volume, bank filtration contributes significantly to drinking water production, by the number of supplies, deep groundwaters comprise an overwhelming majority. We have hypothesised that the natural element composition of water types liable to surface impact (surface water, bank filtered water, karst water and shallow groundwater) would differ from deep groundwater, but only karst waters exhibited a unique profile, dominated by Ca and Mg due to the predominantly dolomite karst water reservoirs. The similarity in the elemental compositions of surface waters (specifically Lake Balaton) and groundwaters in might be an indication of their interconnection. The significance of groundwater discharge into lakes in Hungary was already demonstrated for Lake Balaton [
38] and also for Lake Velence [
33].
In groundwater, the prevalence of natural elements is determined not only by the geological composition of the aquifer but also by the groundwater flow systems that are characterised by different physicochemical conditions (e.g. pH, temperature, redox potential, dissolved gas content), changing from recharge to discharge areas and from local to regional flow systems. These systematic changes control the dissolution of elements and the resulting differences in the composition of the groundwater even within the same aquifer [
39].
Of the geological features of the Pannonian Basin, the sediments originating from the ancient Lake Pannon are among the most important water reservoirs of drinking water and thermal water. As the Lake basin was filled, the deep basin and tubidite formations, followed by the various delta facies change into sediments of fluvial origin at the end of the sequence. These regions, to a large extent, overlap with water sources of high (up to 100 µg/L) arsenic concentration [
40]. For decades, arsenic has been the primary focus of investigations and drinking water quality improvement programmes in Hungary. Since ample historical data was already available for arsenic, it was not part of the current study, but potential overlaps were investigated. High arsenic water sources were evenly distributed between Clusters 3, 4 and 5, and no clear correlation was observed at this level of analysis with the prevalence of other elements. In the ranges of Hungary, carbonates (dolomite, limestone) and volcanics, lesser extent metamorphic and magmatic (granite) rocks are most prevalent. The latter ones can be the source of various elements, among others uranium.
According to the current scientific consensus, U has the highest health relevance of the investigated parameters. Observed uranium concentrations in Hungary are low compared to other European countries [
6,
7,
8,
9,
10]. The highest uranium levels (in Cluster 6) corresponded spatially to locations of elevated gross alpha activity used for the determination of the indicative dose in drinking water [
41]. Some of these locations were investigated previously in detail, approaching the groundwater quality problem from the perspective of groundwater flow systems, taking into account (besides geological factors) the subsurface residence time and flow path of groundwater by a basin-scale survey [
33,
38]. These local studies confirmed that uranium contributed most to the measured elevated gross alpha activity (up to 0.8 Bq/L), which originated from the mobilizing effect of local flow systems due to oxidising conditions. Neither the calculated indicative dose nor U as a nephrotoxic element were identified as a source of significant health hazard.
Se, V and Mo are typically present at concentrations below the LOQ or in low concentrations in Hungarian water sources. Higher values (>5 µg/L) were measured in a limited number of water supplies (8 and 3 for V and Se, respectively), stipulating a local source. In the case of Mo and Se, similar results were obtained in Croatia by Ćurković et al. (2016) in drinking water (median concentration Se <1.0 µg/L; Mo <1.0 µg/L) [
42]. Smedley et al. reported low Mo concentrations in drinking water in England and Wales (<1.0-1.5 µg/L) [
43]. Studies from other European countries identified considerably higher levels (e.g. up to 350 µg/L V in Italy and 40 µg/L Se in France). However, these problems in other countries are also occur in relatively confined geographic areas (e.g. V together with U in volcanic areas in central Italy and Se in the Chalk Aquifer of Northern France). [
10,
29]. In the present study, V occurred in areas localised in the Northern Range volcanic complexes. Moderately strong correlation was seen between U and Se, which occur in the southern foreland of the Velence Hills built-up of granitic rocks.
Lithium was also measured in significant concentrations, even exceeding 250 µg/L in specific locations, which is relatively high compared to most other European studies [
19,
20,
22,
23]. Li is characteristic for the deep groundwater sources, suggesting longer groundwater residence time (i.e. more time for rock-water interactions). Higher concentrations than the ones measured in Hungary were only reported from Austria in Europe [
23]. The protective effect of higher lithium intake via drinking water was also confirmed in a previous Hungarian study using a subset of data from the present study [
24]. The presence of Be, Co and Mo is negligible in Hungarian drinking water sources. Neither Be nor Ba exceeded in any of the samples and Mo only four times (0.3%) the guideline value recommended by WHO (12, 1300 and 70 µg/L, respectively) [
4].
Ca and Mg were dominant in karst waters. However, other aquifers hosted by sedimentary rocks, such as sandstone with carbonate cement may exhibit higher abundance of Ca beside Na and K as dominant ions.
For both Ca and Mg, median concentrations and lower quartile values exceed the minimum concentrations recommended by several EU Member States (30 mg/L and 10 mg/L, respectively) [
15]. However, in-country variability was high, and some regions, especially the Great Plains in southeastern Hungary, population is at risk of insufficient intake of Mg and Ca from drinking water.
5. Conclusions
The study confirmed that none of the investigated parameters pose a risk to drinking water quality and consequently to the health of the consumers on a national scale. However, local hydrogeological conditions may lead to emergence of natural elements in higher concentration in some water supplies. Such local exceedances e.g. of U or V require further investigation and if needed, risk mitigation measures.
In hydrogeological studies, usually the major elements (Ca, Mg, Na, K, Cl, SO4, HCO3) are used for characterization of composition of waters. The presented study highlighted the importance of trace elements and the multivariate data analysis tools, which can reveal links and explain geographical distributions. Nevertheless, involving more data into the analysis (well depth, lithology, flow system evaluation results) would result in more detailed picture about the origin of these elements.
Author Contributions
Conceptualization, B. I. and M.V.; methodology, B.I..; formal analysis, B.I..; investigation, B.I..; data curation, I.B., N.E., K.H-CS..; writing—original draft preparation, B.I..; writing—review and editing, P.B., V.M., A.E..; visualization, I.B., K.H-CS.; supervision, M.V. All authors have read and agreed to the published version of the manuscript.