Julie Lions
email: j.lions@brgm.fr

Nicol Devau
email: n.devau@brgm.fr

Daniel Elster
email: daniel.elster@geologie.ac.at

Denitza D Voutchkova

Birgitte Hansen

Jörg Schullehner

Tanja Petrov Pantić

Katarina Atanasković Samolov
email: k.samolov@gmail.com

Victor Camps
email: victor.camps@icgc.cat

Georgina Arnó
email: georgina.arno@icgc.cat

Ignasi Herms
email: ignasi.herms@icgc.cat

Nina Rman
email: nina.rman@geo-zs.si

Sonja Cerar
email: sonja.cerar@geo-zs.si

Jua Grim
email: j.grima@igme.es

Elena Giménez-Forcada

Juan Antonio Luque-Espinar

Eline Malcuit
email: e.malcuit@brgm.fr

Laurence Gourcy
email: l.gourcy@gmail.fr

B.; D ; Hansen

Marco Rotiroti

M Teresa

Condesso De Melo

Klaus Hinsby

A Broad-Scale Method for Estimating Natural Background Levels of Dissolved Components in Groundwater Based on Lithology and Anthropogenic Pressure

Keywords: natural background level, groundwater, water quality, anthropogenic pressure, trace element, groundwater monitoring

The Water Framework Directive (WFD) requires EU member states to assess the chemical status of groundwater bodies, a status defined according to threshold values for harmful elements and based on/the natural background level (NBL). The NBL is defined as the expected value of the concentration of elements naturally present in the environment. The aim of this study is to propose a methodology that will be broadly applicable to a wide range of conditions at the regional and national scale. Using a statistical approach, the methodology seeks to determine NBLs for SO 4 , As, Cd, Cr, Cu, Ni, Zn, and F based on the lithology of aquifers from which groundwater monitoring data were collected. The methodology was applied in six EU countries to demonstrate validity for a wide range of European regions. An average concentration was calculated for each parameter and chosen water point and linked to a lithology. Based on the dataset created, significant differences between lithologies and pressure categories (urban, agricultural, industrial, and mining) were tested using a nonparametric test. For each parameter, 90th percentiles were calculated to provide an estimation of the maximum natural concentrations possible for each lithology.

Introduction

The Water Framework Directive 2000/60/EC (WFD) [1] requires EU member states to assess the chemical status of groundwater bodies, a status defined by threshold values that should take into account concentrations of elements that may be both naturally present in aquatic environments and discharged by human activities. Member states must therefore characterize the natural background levels (NBLs) for each component and groundwater body while recognizing that a wide range of approaches have been developed. The NBL Water 2021, 13, 1531 2 of 21 is defined as the concentration of a component in a groundwater body that corresponds to zero, or only a very minor anthropogenic alteration or undisturbed conditions [2]. A natural concentration is defined as an element that results entirely from a natural source, whether geological, biological, or atmospheric, under conditions not disturbed by human activities. In addition to natural, a semi-natural origin is also possible, when elements derive naturally from the soil matrix or from the aquifer, but are mobilized by a change of conditions following human intervention, e.g., aquifer dewatering, an accident involving acid or a plume that contains salt and/or water rich in organic matter [1,[START_REF] Shand | The Baseline Inorganic Chemistry of European Groundwaters[END_REF].

Since 2000, the WFD has promoted the establishment of river basin monitoring networks, resulting in a large amount of available data. These data, if well-selected, can contribute to the identification of the expected trace element concentrations in groundwater and a determination of the NBLs. However, the methodology needs to be tailored and harmonized among river basins by taking into account aquifer conditions and anthropogenic inputs. Moreover, groundwater chemistry in a given aquifer provides an important basis for defining the nature and extent of current pollution, and aids in the identification of past contributions and necessary corrective measures.

A key objective of the WFD is the achievement of "good ecological and chemical status" for surface water bodies based on environmental quality standards (EQS) [4] and a "good quantitative and chemical status" for groundwater bodies based on quality standards or threshold values (mainly based on standards listed in the Drinking Water Directive 2020/2184) [START_REF]Directive 2020/2184 of the European Parliament and of the Council of 16 December 2020 on the Quality of Water Intended for Human Consumption[END_REF]. According to the Groundwater Directive 2006/118/EC (GWD) [2] requirements, (1) assessment of flow direction and (2) exchange rates between groundwater bodies and associated surface water are crucial aspects of river basin management plans [START_REF] Ballesteros-Navarro | Methodological proposal for conceptualization and classification of interactions between groundwater and surface water[END_REF]. Table 1 illustrates the standard groundwater and drinking water values. For several elements (i.e., Zn and Cu), a significant difference exists between groundwater threshold values that are usually applied and the EQS applied to surface water [4,7]. This fact underlines the importance of defining NBLs, and is possible when groundwater supplies are associated surface water (e.g., streams, lakes, wetlands, etc.), and may affect chemicallydependent surface water status. In addition, it is useful to know natural concentration ranges and average values to guide, support, and justify the chemical status assessment of surface water and groundwater bodies.

Table 1. Drinking water (DW) standards [START_REF]Directive 2020/2184 of the European Parliament and of the Council of 16 December 2020 on the Quality of Water Intended for Human Consumption[END_REF] and groundwater threshold value (GW-TV). Previous EU studies focused on WFD procedures for estimating the chemical status of groundwater and consideration of NBLs. The BRIDGE (Background Criteria for the Identification of Groundwater Thresholds) approach is based on the generalization that aquifers Water 2021, 13, 1531 3 of 21 with similar petrographic properties have similar water composition under analogous hydrodynamic and hydrologic conditions [START_REF] Appelo | Groundwater and Pollution[END_REF]. The BRIDGE approach uses the BRIDGE aquifer typology to distinguish between NBLs [START_REF] Blum | Background levels under the Water Framework Directive[END_REF][START_REF] Wendland | European aquifer typology: A practical framework for an overview of major groundwater composition at European scale[END_REF].

DW Standards GW-TV

BRIDGE methods consist of using only sampling points unaffected by anthropogenic factors (natural areas, groundwater with low NO 3 concentrations, etc.) [START_REF] Hinsby | European case studies supporting the derivation of natural background levels and groundwater threshold values for the protection of dependent ecosystems and human health[END_REF]. As a result, this approach limits NBL determination to specific areas. Lack of data and analytical method limitations (reporting limit being too high) made it difficult to determine NBLs for most trace elements. For this reason, first estimations were based primarily on bibliographic review [START_REF] Brenot | Identification des Zones à Risque de Fond Géochimique Élevé en Éléments Traces Dans les Cours D'eau et les Eaux Souterraines[END_REF][START_REF] Hobiger | Österreichweite Abschätzung von Regionalisierten[END_REF].

Groundwater chemistry variation depends on a number of factors, including rainfall composition, aquifer lithology, groundwater flow pathways, and residence time [START_REF] Edmunds | The natural (baseline) quality of groundwater: A UK pilot study[END_REF]. Accordingly, each European country developed its own studies to meet its obligation to determine NBLs; several key European studies have been published [START_REF] Shand | The Natural (Baseline) Quality of Groundwater in England and Wales[END_REF][START_REF] Wendland | Derivation of natural background levels and threshold values for groundwater bodies in the Upper Rhine Valley (France, Switzerland and Germany)[END_REF][START_REF] Molinari | Natural background levels and threshold values of chemical species in three large-scale groundwater bodies in Northern Italy[END_REF][START_REF] Coetsiers | Natural background levels and threshold values for groundwater in fluvial Pleistocene and Tertiary marine aquifers in Flanders, Belgium[END_REF][START_REF] Urresti-Estala | Determination of background levels on water quality of groundwater bodies: A methodological proposal applied to a Mediterranean River basin (Guadalhorce River, Málaga, southern Spain)[END_REF][START_REF] Sellerino | Identification of the natural background levels in the Phlaegrean fields groundwater body (Southern Italy)[END_REF][START_REF] Busico | Multivariate statistical analysis to characterize/discriminate between anthropogenic and geogenic trace elements occurrence in the Campania Plain, Southern Italy[END_REF][START_REF] Rotiroti | a new tool to derive natural background levels by the component separation approach: Application in two different hydrogeological contexts in northern Italy[END_REF][START_REF] Nisi | Analysis of complex regional databases and their support in the identification of background/baseline compositional facies in groundwater investigation: Developments and application examples[END_REF]. More broadly, NBLs have been determined for the purposes of groundwater and health risk assessment and drinking water quality management [START_REF] Morgenstern | Groundwater age for identification of baseline groundwater quality and impacts of land-use intensification-The National Groundwater Monitoring Programme of New Zealand[END_REF][START_REF] Missimer | Natural Background and Anthropogenic Arsenic Enrichment in Florida Soils, Surface Water, and Groundwater: A Review with a Discussion on Public Health Risk[END_REF][START_REF] Li | Identifying anthropogenic sources of groundwater contamination by natural background levels and stable isotope application in Pinggu basin, China[END_REF]. Based on the numerous data, statistical methods have been applied, e.g., frequency distributions, probability plots, box plots, and histograms. Concentrations that deviate from the basic distribution (between the 10th and 90th percentiles) are generally excluded from the derivation of background concentrations. Masciale et al. [START_REF] Masciale | Assessing Natural Background Levels in the Groundwater Bodies of the Apulia Region (Southern Italy)[END_REF] used Huber's non-parametric test to identify and eliminate anomalous data before NBL calculation, also introducing a confidence level for NBL values. Recent studies have also considered model-based statistical approaches, such as iterative 2-δ techniques, maximum normalized residual tests, and geostatistical multi-model approaches [START_REF] Li | Identifying anthropogenic sources of groundwater contamination by natural background levels and stable isotope application in Pinggu basin, China[END_REF][START_REF] Guadagnini | Probabilistic assessment of spatial heterogeneity of natural background concentrations in large-scale groundwater bodies through Functional Geostatistics[END_REF][START_REF] Kim | Determination of natural backgrounds and thresholds of nitrate in South Korean groundwater using model-based statistical approaches[END_REF][START_REF] Gao | Assessing natural background levels in shallow groundwater in a large semiarid drainage Basin[END_REF][START_REF] Molinari | Geostatistical multimodel approach for the assessment of the spatial distribution of natural background concentrations in large-scale groundwater bodies[END_REF]. Furthermore, an online tool evaluation of natural background levels (eNaBLe) was developed in Italy in 2020 [START_REF] Parrone | an On-Line Tool to Evaluate Natural Background Levels in Groundwater Bodies[END_REF].

The purpose of this study is to propose and test a methodology common to multiple countries, taking into account national specificities and available data related to groundwater, lithology, and anthropogenic pressure. Prior to this work, previous studies had been conducted at different scales and in specific areas (high mountain karst aquifers, specific geologic units) using data provided by groundwater monitoring networks [START_REF] Herms | Evaluation of natural background levels of high mountain karst aquifers in complex hydrogeological settings. A Gaussian mixture model approach in the Port del Comte (SE, Pyrenees) case study[END_REF][START_REF] Serianz | Hydrogeochemical characterization and determination of natural background levels (NBL) in groundwater within the main lithological units in Slovenia[END_REF][START_REF] Brielmann | Hydrochemie und Hydrogeologie der Österreichischen Grundwässer und Deren Natürliche Metall-und Nährstoffgehalte (Update Geohint 2018[END_REF][START_REF] Lions | Evaluation des Concentrations de Rérérence des Fonds Hydrogéochimiques des Eaux Souterraines par Lithologie D'aquifères[END_REF][START_REF] Voutchkova | Roadmap for Determining Natural Background Levels of Trace Metals in Groundwater[END_REF][START_REF] Cruz | Monitoring of the groundwater chemical status in the Azores archipelago (Portugal) in the context of the EU water framework directive[END_REF]. Each study has adapted the BRIDGE approach methodology in different ways. In some cases, data from wells with only a median nitrate concentration below 10 mg/L were used (e.g., Slovenia, France, Austria). In others, such as Slovenia, sites with up to 20% anthropogenic impact in recharge areas were considered [START_REF] Mezga | Natural Hydrogeochemical Background and Dynamics of Groundwater in Slovenia[END_REF]. In this study, a new methodology aims to work at a large scale (river basin or national scale) and focus on the components presented in Table 1.

However, this exercise faces two primary challenges in particular, because the data used are provided by monitoring under the WFD framework, which was not primarily designed to characterize NBLs. The challenges are:

•

Monitoring network sampling points that are selected to be representative for a groundwater body, but may not be representative for the NBL calculation.

•

Laboratory reporting limits (LRLs) that may not be in the most appropriate ranges for NBL calculation, because they were initially designed to verify compliance with drinking water standards.

Methodology

Background

Reference Concentrations

Several approaches can be used to estimate the ranges of natural concentrations for water management plans, either by determining the average value and standard deviation or by describing low and high values after eliminating outliers. Spatial heterogeneity and temporal variability of groundwater composition make it difficult to determine a single value. Heterogeneity is scale-dependent, which leads to the consideration of a concentration range or a confidence limit adapted to each scale of work [START_REF] Grima | Methodological approach for the analysis of groundwater quality in the framework of the Groundwater Directive[END_REF]. Indeed, upper limits for groundwater NBLs may exceed environmental or drinking water quality standards, particularly for elements such as As and F [START_REF] Preziosi | Natural background levels and threshold values derivation for naturally As, V and F rich groundwater bodies: A methodological case study in Central Italy[END_REF][START_REF] Alcaine | Hydrogeochemical controls on the mobility of arsenic, fluoride and other geogenic co-contaminants in the shallow aquifers of northeastern La Pampa Province in Argentina[END_REF].

Censored Data

A primary challenge in using statistical approaches to process geochemical datasets is the presence of values that are below the limit of quantification (LOQ)/detection limit (LOD), also known as left-censored values. In some datasets, more than 50% of the values may be left-censored. Sophisticated methods require, for example, work on the data distribution law (distributional methods) or an extrapolation from >LOQ data to <LOQ data (robust methods), and are therefore less direct than the substitution of the value by a constant [START_REF] Helsel | Techniques of Water-Resources Investigations of the United States Geological Survey. Hydrologic Analysis and Interpretation: Statistical Methods in Water Resources[END_REF]. The data substitution method also depends on the amount of data below the LOQ. For example, to estimate the median, Antweiler (2015) [START_REF] Antweiler | Evaluation of Statistical Treatments of Left-Censored Environmental Data Using Coincident Uncensored Data Sets. II. Group Comparisons[END_REF] confirmed that, for a low-censoring dataset, LOQ/2 can be used. In addition, when 50% of the data are censored, no technique gives a reliable distribution estimate. Other more robust approaches for the imputation of left-censored values take multivariate approaches into account [START_REF] Palarea-Albaladejo | zCompositions-R package for multivariate imputation of left-censored data under a compositional approach[END_REF].

Analytical methods have evolved significantly in recent years, and the quality of the results has improved significantly, reducing uncertainties, detection limits, and/or quantification limits. For these reasons, the temporal evolution of the methods and laboratory reporting limits needs to be reviewed before selecting the dataset (see Section 2.3).

Uncertainties

Before interpreting water quality data, it is important to consider the uncertainty associated with the measurement. Traditionally, only analytical uncertainties have been applied to water quality interpretation [START_REF] Guigues | Estimation of the measurement uncertainty, including the contribution arising from sampling, of water quality parameters in surface waters of the Loire-Bretagne river basin, France[END_REF]. However, these uncertainties represent only part of the uncertainty related to interpreting environmental measurements. Water sample collection is also part of the measurement acquisition process. Ghestem (2009) [START_REF] Ghestem | Incertitudes Liées à L'échantillonnage: Exemples D'estimation sur eau de Surface et eau Souterraine, Aquaref[END_REF] evaluated the overall uncertainty in groundwater measurement. Overall variations were estimated at 5-10% for most substances, including the major components (SO 4 , NO 3 , Cl) and a range of 10 to 35% for trace elements. For As, Cu, Cr, Ni, and Zn (>5 µg/L), the uncertainties are between 15 and 20%.

Proposed Methodology

Based on the WFD context and data limitations, the following methodology (Figure 1) shows the steps for determining NBLs based on groundwater network data and a consideration of the potential anthropogenic effects.

Groundwater Quality Dataset

The method involves processing data stored in national/regional groundwater monitoring databases (Table 2), which usually contain large datasets. The first step (step 1, Figure 1) is an inventory of the available and usable data; this step involves the selection of relevant parameters for NBL calculation and information on the lithological and hydrogeological context of each sampling point. Additional chemical parameters are needed to determine the hydrochemical groundwater conditions, such as pH, redox potential, NO 3 , and dissolved O 2 and Fe. Thermal and mineral water are excluded from this assessment. After considering data availability and laboratory reporting limits (LOD/LOQ), step 2 is an examination of the time period for which the data are considered (reference period) for each dataset (step 2).

Groundwater Quality Dataset

The method involves processing data stored in national/regional groundwater monitoring databases (Table 2), which usually contain large datasets. The first step (step 1, Figure 1) is an inventory of the available and usable data; this step involves the selection of relevant parameters for NBL calculation and information on the lithological and hydrogeological context of each sampling point. Additional chemical parameters are needed to determine the hydrochemical groundwater conditions, such as pH, redox potential, NO3, and dissolved O2 and Fe. Thermal and mineral water are excluded from this assessment. After considering data availability and laboratory reporting limits (LOD/LOQ), step 2 is an examination of the time period for which the data are considered (reference period) for each dataset (step 2). Dataset construction and formatting are based on the extraction of data and metadata (e.g., sampling date, unit, analytical method, LOD/LOQ) from national or regional databases (step 3) (Table 2), followed by the use of a set of controls and corrections prior to any treatment (step 4). Controls are applied to check systemic errors that may be observed in large databases (unit mg/L instead of µg/L), duplicated values (multiple uploads), and anomalous data (pH, temperature); when an error can be corrected (generic error), it is possible to substitute the value. Otherwise, the sample is deleted.

A large proportion of trace elements were censored at various LRLs. The censoring level was adjusted (Table 2), thus providing a lower LRL for each element. This is a simple way to deal with non-detected data, but it is a rather arbitrary approach depending on the analytical method and laboratory. In the case of the LOQ exceeding expected NBL ranges, it is more appropriate to remove the censored value with such a high LOQ from the dataset to avoid bias [START_REF] Lions | Evaluation des Concentrations de Rérérence des Fonds Hydrogéochimiques des Eaux Souterraines par Lithologie D'aquifères[END_REF]. This removal may require using an iterative approach.

The number of analyses and parameters varies for each dataset. The objective is to have one value for each parameter per sampling point (step 5), so that each sampling point has the same representativeness in the dataset. When several analyses (>3) were available for water points, we chose to calculate a median for these datasets so as to assign a single value to each sampling point. The median value is preferred because it is a more robust value than the mean as a result of being less affected by outliers. When several LOQ levels were reported for one sampling point, we applied the median LOQ.

Summary statistics, which include the minimum, maximum, and common percentiles, were computed for all parameters. Particular attention must be paid to the number of analyses available before conducting this first statistical analysis. In addition, redox and pH types were defined for each sampling point (step 6) according to the following classes (Tables 3 and4): pH and redox states that control mineral solubility and trace metal mobility, often through sorption/desorption processes for pH and the redox process for redox potential. If no direct measurements are available, dissolved O 2 , NO 3 , SO 4 , and Fe can be used to estimate the groundwater redox state based on the definition in Table 4.

Hydrogeological Characteristics of Sampling Points

For each sampling point, hydrogeological information (aquifer, lithology, depth) are linked to the data (step 7). In addition (step 8), lithological information related to water wells aids in the definition of local geology at the sampling point and in controlling consistency with the assigned groundwater body. For a common approach among all investigated aquifers, major categories were defined on the basis of eight major simplified lithologies (sedimentary: sand, gravel, carbonates, clay/marls, others; volcanic rocks; crystalline bedrock; and metamorphic rocks). This global approach makes it possible to define NBLs for simplified lithologies. Finally, each sampling point is linked to its groundwater catchment or associated surface watershed (step 9). This information is necessary for the assessment of anthropogenic pressure at each sampling point. For superficial aquifers, the link is quite important, whereas, for deep aquifers that are disconnected from the surface, the link is not critical. This link is also an illustration of groundwater vulnerability for the following steps.

Anthropogenic Pressure

Evaluating anthropogenic pressure and relating activities to specific dissolved components would make it possible to determine the expected NBLs in the area under agricultural, industrial (including mining), and urban conditions. The accuracy of this approach would depend on the information regarding anthropogenic activities and the relationship between the activities and dissolved components released. The pressure inventory (step 10) is provided for anthropogenic pressure, such as diffuse pollution (agriculture, urban), point source pollution (industrial), and mining areas. Anthropogenic pressure is determined on the basis of the Corine Land Cover (CLC) inventory and the following classification: urban pressure; industrial pressure; agricultural pressure; and mining activities. This inventory is complete when databases related to industrial sites and service activities, as well as a register of pollution discharges, are available (Table 2).

All information for each water point and anthropogenic pressure is summarized using geographical information systems (GISs) (step 11). This facilitates a mapping approach to check areas with abnormal concentrations or exceedances of DWS. Due to the heterogeneity of the available data and anthropogenic conditions among countries, different GIS approaches have been tailored to determine the prevailing pressure at each sampling point (step 12). This approach can be performed at the groundwater body scale, by lithology, or for the entire dataset. The approach for determining the pressure area around a water point will vary depending on hydrogeological and watershed characteristics. The goal of this step is to identify elements affected by anthropogenic activity (after a statistical test to evaluate the impact on concentration distribution) and to discard data for which an anthropogenic pressure is recognized.

Statistical Treatments

Using one average value for each water point, descriptive statistics (step 13) are calculated for each parameter:

•

Median-Quartile 1 and 3-10th percentile and 90th percentile • Percentage of quantified data by group in order to assess the weight of unquantified (censored) values in the calculation

The description of the distribution is illustrated by boxplots or a cumulative frequency plot. Cumulative probability curves may be preferred to facilitate data comparison, because they are useful for evaluating different data populations. Spatial distributions were highlighted by mapping with GIS tools. Outliers potentially illustrate water with anthropogenic contamination or a high geochemical background, which are to be distinguished from the NBLs defined in this study.

Discriminant function analysis (DFA) (step 14) is applied to find patterns within the dataset and classify them. DFA is useful for determining whether water points can be grouped into one water family for NBL determination or if subgroups are more appropriate.

Distribution comparison and non-parametric tests (step 15) may be applied to each parameter to determine if sampling points have the same concentrations under anthropogenic pressure or to compare different classes of sampling points using simplified lithology or geochemical type (acidic or basic water, redox conditions) so as to determine whether lithology or geochemical parameters may affect NBLs. This step makes it possible to identify water families for which NBLs need to be determined. Assumptions of data normality and homogeneity of variances between types within a domain are not always met. In this situation, a parametric test is less powerful than non-parametric alternatives. In a dataset containing groups with large sample numbers (>100), a parametric approach, such as ANOVA, may still be suitable. Nevertheless, groups with a limited amount of data (five to 20) are included in the datasets in this study. For this reason, non-parametric one-way tests, such as Mann-Whitney and Kruskal-Wallis tests, were applied to test the null hypothesis assuming identical distribution for different groups [START_REF] Helsel | Techniques of Water-Resources Investigations of the United States Geological Survey. Hydrologic Analysis and Interpretation: Statistical Methods in Water Resources[END_REF]. Results from the non-parametric tests were analyzed with the understanding that they were less powerful than parametric tests. Non-parametric statistics do not require distributional assumptions, just mutual independence of samples. The significance of differences is established from p < 0.05. In the event that the Kruskal-Wallis test was significant, that is, the null hypothesis could be clearly rejected, a post hoc analysis (Nemenyi test) was performed to determine which groups were identified in the dataset.

NBL Calculation

Based on the results of step 15, sampling points identified as being affected by anthropogenic pressure were removed from each parameter data subset (step 16).

NBL determination (step 17) consists of a calculation of natural maximum concentrations representative of a water family. A water family is a lithology type or groundwater body with a dominant lithology. NBLs are calculated for SO 4 , As, Cd, Cr, Cu, Ni, Zn, and F within each lithology. When geochemical parameters, such as pH and redox, are identified as affecting the data distribution (based on the Kruskal-Wallis test result, steps 14-15), the NBL is determined for subgroups. For each water family, the 90th percentile is determined. This percentile is usually applied as a cutoff to define natural concentrations by lithology without any anthropogenic effects or geogenic anomaly; the 90th percentile defines the concentration at which 90% of sampling points associated with a simplified lithology are of a lower concentration. The concentration range between the 10th and 90th percentiles represents the element's "expected" concentration [START_REF] Hinsby | European case studies supporting the derivation of natural background levels and groundwater threshold values for the protection of dependent ecosystems and human health[END_REF]. For this reason, the 90th percentile was selected for NBL determination in this study.

Results of Methodology Application to Regional Studies

The NBL study was conducted in six [START_REF] Ballesteros-Navarro | Methodological proposal for conceptualization and classification of interactions between groundwater and surface water[END_REF] countries; a summary of groundwater chemistry data is shown in the following tables and detailed in the Lions et al. report [START_REF] Lions | HOVER Deliverable 3.3-Data Set of the Results of the Statistical Data Treatment Allowing the Preparation of the Raw Elements for the Tasks 4 and 5 i.e., Concentrations of Elements of Natural Origin Per Typologies[END_REF]. Each dataset was based on the available analysis with consideration of the specificities and limits of each dataset (Tables S1 to S7). In several cases, the LRL varies significantly from year to year and between laboratories. For example, LRLs for As can range from 0.1 to 10 µg/L in France, whereas, in Denmark, they range from 0.01 to 1 µg/L. For the highest LOQs, which may be higher than NBLs, only data with a maximum LOQ were considered. LOQ cutoff values were optimized according to data distribution, because data with a large amount of high LOQs (i.e., higher than NBLs) disturb the natural concentration distribution [START_REF] Lions | Evaluation des Concentrations de Rérérence des Fonds Hydrogéochimiques des Eaux Souterraines par Lithologie D'aquifères[END_REF]. In the internal basins of Catalonia, assuming that a censored value may imply decreased data quality, only parameters (As, Cu, Zn, SO 4 ) with a LOD/LOQ below 33% were considered in the analysis. By considering a large amount of data < LOQs, 90th percentiles may be identical to the LOQ/2 or <LOQ; in this case, the NBL is indicated as <LOQ, and no concentration is calculated.

Anthropogenic Pressure (Steps 10-15)

Various GIS approaches have been used to determine the pressure at each sampling point.

Watershed (Loire-Bretagne)

In the Loire-Bretagne river basin, the methodology consists of linking the prevailing pressure to watershed units at each sampling point (Table S8). Watershed units are determined using the TauDEM (Terrain Analysis Using Digital Elevation Models) toolbox plugin combined with the Digital Terrain Model (DEM). This DEM (BD ALTI) has a spatial resolution of 25 m [START_REF] Devau | Etude par Approche Globale des Fonds Hydrogéochimiques des Eaux Souterraines Sur Le Bassin Loire-Bretagne[END_REF]. The territory was divided into 36,453 basins and sub-basins. Upstream and downstream basins were identified for each basin. Sampling points and pollution sources, such as discharge points, industrial sites, mining areas, non-point sources, or diffuse pollution per municipality and CLC [START_REF] Land | [END_REF], were linked to each watershed.

Individual Analysis

In Slovenia, for groundwater bodies that were declared to be of poor quality, all sampling points were listed as having anthropogenic impacts. Surrounding sites were investigated where potential pollution sources or urban areas were identified. This approach was used because the number of sampling points was limited. This data management led to the removal of 24 sampling points from the original dataset of 203 and one sampling point as an outlier for As, Cd, Cr, and Fe.

Buffer around the Sampling Point

For each sampling point, the areal proportion of CLC inventory types was calculated for the area surrounding sampling points (Table 2, Tables S9 andS10). Based on the large number of wells studied in Denmark (Figure 2) and Catalonia (>1000), it was decided that a one (1) km buffer would be an adequate proxy for well catchment zones, because actual well catchment areas are unknown at these national scales (Figure 3). In Austria, CLC [START_REF] Land | [END_REF] information was extracted as point data (10 m buffer), because most sampling points are found in gravel aquifers with complex catchment areas.

DRASTIC Method

The possibility of contaminants infiltrating from the surface in the Fruška Gora national park area was studied using the DRASTIC method [START_REF] Aller | A Standardized System for Evaluating Groundwater Pollution Potential Using Hydrogeologic Settings[END_REF], which was developed to evaluate groundwater vulnerability and make vulnerability maps [START_REF] Pantić | Hydrogeology and vulnerability of groundwater to polluting in the area of national park "Fruška Gora[END_REF][START_REF] Veljković | Vulnerability of groundwater pollutants in the area of national park "Fruška Gora[END_REF]. This study noted that, in the central part of Fruška Gora (ridge), the vulnerability level was low (75-100) to very low (<75). Down the Fruška Gora slope, the vulnerability level increased; thus, ridges in the mountains, which are densely populated, are the most vulnerable (150-185). By combining DRASTIC with CLC, anthropogenic influences (discontinuous urban fabric, industrial units, mineral extraction sites, agriculture areas) can be seen/are clearly observable in areas of high vulnerability. This approach can be useful for an improved determination of NBLs, because prevailing pressure depends on all DRASTIC factors and hydrogeological settings.

Buffer around the Sampling Point

For each sampling point, the areal proportion of CLC inventory types was calculated for the area surrounding sampling points (Table 2, Tables S9 andS10). Based on the large number of wells studied in Denmark (Figure 2) and Catalonia (>1000), it was decided that a one (1) km buffer would be an adequate proxy for well catchment zones, because actual well catchment areas are unknown at these national scales (Figure 3). In Austria, CLC [START_REF] Land | [END_REF] information was extracted as point data (10 m buffer), because most sampling points are found in gravel aquifers with complex catchment areas.

DRASTIC Method

Buffer around the Sampling Point

For each sampling point, the areal proportion of CLC inventory types was calculated for the area surrounding sampling points (Table 2, Tables S9 andS10). Based on the large number of wells studied in Denmark (Figure 2) and Catalonia (>1000), it was decided that a one (1) km buffer would be an adequate proxy for well catchment zones, because actual well catchment areas are unknown at these national scales (Figure 3). In Austria, CLC [START_REF] Land | [END_REF] information was extracted as point data (10 m buffer), because most sampling points are found in gravel aquifers with complex catchment areas.

DRASTIC Method

Geochemistry (pH, Redox, Lithology) (Steps 13-15)

Data analysis addresses the geological composition of aquifers and geochemical parameters (pH, redox), which are among the main factors that affect trace element occurrence. Some examples are presented here; for additional details, see Lions et al. [START_REF] Lions | HOVER Deliverable 3.3-Data Set of the Results of the Statistical Data Treatment Allowing the Preparation of the Raw Elements for the Tasks 4 and 5 i.e., Concentrations of Elements of Natural Origin Per Typologies[END_REF].

Lithology is the primary factor that controls water composition, including pH, conductivity, major ions, and trace elements. In Austria, major cations Na, K, Ca, and Mg were considered as quantitative variables to characterize aquifer lithology using a confusion matrix (Table S11). DFA (specifically, linear discriminant analysis) showed that, for Denmark, based on predictor variables (pH, SO 4 , As, Ni, F, Cl, O 2 , NO 3 , Fe), it was possible to discriminate easily between the two lithology classes sedimentary: carbonates and sedimentary: sand (accuracy and 95% confidence interval: 81.7% (79.3-83.9%)). Therefore, these lithology classes were used for determining Danish NBLs. For Slovenia, the Mann-Whitney U test results showed significant differences between NBLs of carbonates and sands at the 95% confidence level (α = 0.05). This was true for all elements except Cd, due to its low concentration (0.01 µg/L).

In Loire Bretagne, the Kruskal-Wallis test and Nemenyi post hoc test showed that trace elements belong to different groups that overlap when the lithology is not well-constrained, such as for clays and marls or gravels (Figure 4). Results show that volcanic rocks belong to an individual group for components such as Ni, Zn, F, or SO 4 . Crystalline rocks are in a group of their own for elements such as As (including the group sedimentary-others, which contains crystalline and sedimentary formations) and Cr, whereas Ni is individualized for metamorphic rocks and sand. The distribution of an element across lithology classes was tested to justify the NBL definition for each lithology. In addition, descriptive statistics showed that crystalline and metamorphic rocks were primarily represented by acidic water (>90% with a pH < 7). This fact is directly related to the water-rock interaction in these lithologies. Conversely, carbonates and clays/marls were represented by neutral water (7-7.5) (>83%), whereas basic water was present at 12% of the sampling points in carbonates. Water in gravels can be either acidic or neutral, whereas sand and volcanic aquifers are represented in each family (acidic, basic, neutral). Redox conditions will depend more on hydrogeological settings (depth, confined aquifer) than on lithology. The Kruskal-Wallis test was used to evaluate the effect of pH and redox on trace element concentrations for each element (Figure S1). Reduced water predominates in Denmark, whereas oxic water is predominant in Loire-Bretagne. However, As is affected by the redox class, as illustrated in Figure 5 for sedimentary aquifers, whereas pH is not discriminant. The Kruskal-Wallis test was used to evaluate the effect of pH and redox on trace element concentrations for each element (Figure S1). Reduced water predominates in Denmark, whereas oxic water is predominant in Loire-Bretagne. However, As is affected The Kruskal-Wallis test was used to evaluate the effect of pH and redox on trace element concentrations for each element (Figure S1). Reduced water predominates in Denmark, whereas oxic water is predominant in Loire-Bretagne. However, As is affected by the redox class, as illustrated in Figure 5 for sedimentary aquifers, whereas pH is not discriminant.

NBLs (Steps 16-17)

The following NBLs were determined for each dataset based on aquifer lithology. Concentrations that deviated due to anthropogenic factors were excluded from the dataset prior to NBL derivation according to the conclusions of each dataset. Only families with more than five (5) sampling points are presented in Tables 5 and6. For Denmark, the minimum number of sampling points for a lithological group was set at 20, but when the pH and redox were used, this threshold was increased to 50.

The resulting NBLs considering pH and redox conditions were determined for As, Ni, SO 4 , and F (Table S12). For the Loire-Bretagne area, the results showed that Ni and Zn were clearly affected by acidic conditions with higher concentrations. Reduced conditions led to high concentrations of Fe, Mn, and As. For this dataset, NBLs were determined for As based on redox only and for Ni and Cu based on redox and pH. For carbonates, only neutral and basic waters were considered. NBLs for Zn in sandy aquifers were considered based on pH (Table S13).

Discussion

Lithology

For the elements in this study, lithology is one of the primary factors that controls NBLs. The results show that using simplified aquifer lithologies is appropriate for comparing large datasets from different countries. Further subdivisions of lithology groups would lead in many cases to very low numbers of data per lithology. The NBL variations shown in Tables 5 and6 can be partly explained by the broad spectrum of elements found. NBLs can also show multiple and different geological settings (e.g., paleo-environment) that explain the possible variations of a parameter across several countries. However, the approach of this study highlights valuable primary characteristics of lithology groups, including NBL ranges for specific elements.

A discussion of lithology groups by country may be helpful in clarifying NBL variations that can be linked to lithology.

For Denmark, it is important to distinguish between different depositional ages for the sedimentary: sand class. For example, the DFA analysis showed a difference between quaternary sands (mostly of glacial origin) and pre-quaternary sands when pH, SO 4 , As, Ni, F, Cl, O 2 , NO 3 , and Fe are used as predictors [START_REF] Lions | HOVER Deliverable 3.3-Data Set of the Results of the Statistical Data Treatment Allowing the Preparation of the Raw Elements for the Tasks 4 and 5 i.e., Concentrations of Elements of Natural Origin Per Typologies[END_REF].

We correctly identified some national specificity. In Denmark and Serbia, the NBL of As in sandy aquifers exceeds threshold values. In Denmark, groundwater resources have often geogenic As with high concentrations [START_REF] Ramsay | Drinking Water Criteria for Arsenic in High-Income, Low-Dose Countries: The Effect of Legislation on Public Health[END_REF]. The Serbian dataset refers to a location north of Fruška Gora (Vojvodina), where groundwater usually contains elevated As [START_REF] Dalmacija | Environment in Autonomous Province of Vojvodina; Provincial Secretariat for Urbanism, Construction and Environmental Protection[END_REF][START_REF] Jovanović | Arsenic occurrence in drinking water supply systems in ten municipalities in Vojvodina Region, Serbia[END_REF][START_REF] Petrović | Macro-and microelements in bottled and tap waters of Serbia[END_REF].

The sedimentary: others class includes a wide range of area-specific lithologies, usually described as having distinct characteristics linked to local geology (specific environment, mineralization). Therefore, it is not a well-defined lithology, and NBL determination will reflect each regional lithology. In Loire-Bretagne, the sedimentary: others class includes aquifers composed of a mixture of sedimentary and other lithologies, such as volcanic strata or weathered crystalline or metamorphic rocks. In the Duero River Basin, there is an undifferentiated sequence of interbedded sedimentary sand, gravel, and clay. In such multilayered aquifer faces, NBLs are calculated for all sedimentary units. For this reason, all of them were aggregated into a single category (cenozoic). In Slovenia, the sedimentary: others class points that sample quaternary alluvial aquifers are mixtures of sand and gravel. Results for the clays/marls class are not presented here, because they were found only in France and parts of Austria [START_REF] Lions | HOVER Deliverable 3.3-Data Set of the Results of the Statistical Data Treatment Allowing the Preparation of the Raw Elements for the Tasks 4 and 5 i.e., Concentrations of Elements of Natural Origin Per Typologies[END_REF].

Water 2021, 13, 1531 14 of 21 Furthermore, regarding the crystalline or metamorphic classes, all formations cannot be merged into only one lithological class, because each geological domain (e.g., variscan, hercynian) has specific petrologies and mineralization that may affect the NBLs differently. This approach also excludes the contribution of local mineralization, faults, and mineral deposits. This study provides several standard NBLs for crystalline and metamorphic rocks, but further studies are needed to characterize local NBLs.

Effects of Physico-Chemical Parameters

Lithology and water chemistry control pH buffering in groundwater: acidic water in crystalline aquifers and neutral to basic water in carbonate aquifers. These factors will directly control the geochemical equilibrium for trace element release and resulting NBLs. In addition, redox potential may be directly controlled by depth or confined/unconfined conditions. In this case, redox potential is an indicator of groundwater conditions, which is correctly demonstrated by reduced water speciation (low NO 3 , low Fe, low dissolved O 2 ) in deep or confined aquifers. In the Loire-Bretagne area and Denmark, F concentrations differ according to redox classes. For this element, redox is not directly involved in F mobilization, because it is not sensitive to redox processes; however, F is higher in confined aquifers and deep groundwater.

Anthropogenic Effects

This land use-based method for assessing anthropogenic effects aids in the identification of statistical differences between sampling points subject to prevailing agricultural pressure compared to those affected by industrial and/or urban pressures. Sampling point distribution according to prevailing pressures is presented in [START_REF] Lions | HOVER Deliverable 3.3-Data Set of the Results of the Statistical Data Treatment Allowing the Preparation of the Raw Elements for the Tasks 4 and 5 i.e., Concentrations of Elements of Natural Origin Per Typologies[END_REF] and Table S8-S10. Regarding statistical significance among various lithologies and anthropogenic impacts, it is possible to conclude that prevailing pressures may affect trace element concentrations in several cases.

In the Duero River Basin (Spain), most sampling points are mainly affected by agricultural pressure. In other words, within a 1 km buffer around each sampling point, most of the land is used for agricultural purposes. The pre-selection of sampling locations to include only drinking water supply wells may explain the absence of points affected by industry, mining, or urban activities.

In the Danish case, less than 1% of the sampling points were devoid of any anthropogenic pressure (i.e., are natural). All are located near the coast, and are not considered to be representative of undisturbed conditions. Depending on the components considered, between 74.8 and 85.8% of the sampling points show effects of agricultural pressure. This finding limits the application of the method, because it means that a statistical test for comparing elemental distributions cannot be used against the natural group. As a result, the normal state is assumed to be dominant agricultural pressure. In this case, only sampling point groups for which industrial or urban pressure predominates can be compared to the agricultural group. The shortcoming of this method is that it fails to distinguish between polluted and potentially polluted sampling points, but excludes the entire group with the specific dominating pressure. As a consequence, an entire data group is excluded, and some NBLs increase because sampling points with low concentrations (potentially the NBL) were excluded.

In Slovenia, a previous study [START_REF] Cruz | Monitoring of the groundwater chemical status in the Azores archipelago (Portugal) in the context of the EU water framework directive[END_REF] used only sampling points with a minimum probability of being affected by anthropogenic activities, and therefore determined NBLs for NO 3 (3.8 mg/L) and SO 4 (6.9 mg/L). [START_REF] Serianz | Hydrogeochemical characterization and determination of natural background levels (NBL) in groundwater within the main lithological units in Slovenia[END_REF] [START_REF] Serianz | Hydrogeochemical characterization and determination of natural background levels (NBL) in groundwater within the main lithological units in Slovenia[END_REF] assigned 6.0 mg/L for NO 3 after removing all analyses exceeding 10 mg/L of nitrate; SO 4 NBLs were 17.0 mg/L (by the probability method) and 26.0 mg/L (by the pre-selection method). In this study, the natural group NBLs are 15.5 mg/L for NO 3 and 14.0 mg/L for SO 4 , but as much as 41.5 mg/L and 25.0 mg/L for the anthropogenic group, respectively. This result indicates that many points in the Slovenian national monitoring network used for groundwater body assessment are noticeably affected by anthropogenic pressures, so trace element concentrations should be evaluated with care. As a result, only 87 sites without obvious anthropogenic impact were used for the NBL calculation. Thus, it is important to carefully assess the statistical test results (step 15, Figure 1), which can serve as a first screening followed by expert assessment before data removal. In addition, a more detailed assessment of anthropogenic pressure and exclusion of polluted sites can be done for a target area within a specific groundwater body or aquifer.

The results of each prevailing pressure assessment (Table 7) could not be generalized to other groundwater bodies, because they are directly linked to the distribution and nature of anthropic activities. However, anthropogenic impacts have been statically observed for each trace components (As, Cd, Cr, Cu, Ni, Zn, F, SO 4 ) under the following pressures: agriculture (As, Cd, Cu, Ni, Zn), urban and industrial (all), and mining (As, Cu). Other trends may exist, but they have not been observed in the investigated datasets.

NBLs

Discriminant analyses of lithology classes (qualitative variables) and selected major and minor elements (quantitative variables) were performed to characterize the most prevalent water families in each dataset. While the median is representative of the average concentration, the 90th percentile gives a high range for natural concentrations [START_REF] Griffioen | Comparison of Selection Methods to Deduce Natural Background Levels for Groundwater Units[END_REF]. The 90th percentile makes it possible to estimate the possible concentrations of natural origin in each group. Values defined by higher percentiles, for example, 95 or 97.7%, would be also appropriate as a reference for maximum natural concentrations, but higher percentiles require larger datasets and well-selected sampling points for higher confidence intervals [START_REF] Edmunds | Natural Groundwater Quality[END_REF]. When a large variability is observed, NBLs could be much higher than the average concentration in groundwater. Therefore, NBLs are not average concentrations, but concentrations in a higher range. These values are helpful to emphasize anomalous concentrations in groundwater. The NBLs obtained from these seven case studies are well within the range of NBLs determined by international studies [START_REF] Tedd | Research 183: Assessing and Developing Natural Background Levels for Chemical Parameters in Irish Groundwater[END_REF].

In this study, we were able to apply the proposed method and establish NBL values that are within the expected range compared to previous studies [START_REF] Serianz | Hydrogeochemical characterization and determination of natural background levels (NBL) in groundwater within the main lithological units in Slovenia[END_REF][START_REF] Brielmann | Hydrochemie und Hydrogeologie der Österreichischen Grundwässer und Deren Natürliche Metall-und Nährstoffgehalte (Update Geohint 2018[END_REF][START_REF] Lions | Evaluation des Concentrations de Rérérence des Fonds Hydrogéochimiques des Eaux Souterraines par Lithologie D'aquifères[END_REF][START_REF] Voutchkova | Roadmap for Determining Natural Background Levels of Trace Metals in Groundwater[END_REF][START_REF] Cruz | Monitoring of the groundwater chemical status in the Azores archipelago (Portugal) in the context of the EU water framework directive[END_REF][START_REF] Grima | Methodological approach for the analysis of groundwater quality in the framework of the Groundwater Directive[END_REF][START_REF] Catalana De L'aigua | (ACA). Pla de Gestió del Districte de la Conca Fluvial de Catalunya: Període 2016-2021[END_REF]. In the case of the Duero River Basin, NBLs are compared to threshold values provided by local authorities. In the case of SO 4 , the NBL obtained when considering the cenozoic aquifer is about 241 mg/L (anthropogenically affected value), while the value without the anthropogenically affected data is about 100 mg/L.

For groundwater management, these NBLs help to highlight groundwater bodies where high concentrations need further investigation, so as to identify the origin of the anomaly. An anomalous concentration in a groundwater body may be a result of groundwater contamination, but may also come from a natural occurrence in a specific aquifer or in a particular geological context, such as mineralization or thermal waters [START_REF] Binda | Anomalously high concentration of Ni as sulphide phase in sediment and in water of a mountain catchment with serpentinite bedrock[END_REF][START_REF] Hossain | Hydrogeochemical characterisation and health hazards of fluoride enriched groundwater in diverse aquifer types[END_REF][START_REF] Cinti | Spatial distribution of arsenic, uranium and vanadium in the volcanic-sedimentary aquifers of the Vicano-Cimino Volcanic District (Central Italy)[END_REF]. Even so, the NBLs established in this study should not be applied directly to management and planning purposes at the groundwater body scale, because the data processing was conducted at a different scale (regional, national). In addition, in the case of high geological heterogeneity, such as in the Catalonian aquifers, NBL estimation needs more criteria to define water families for groundwater management implementation.

Conclusions

The methodology developed in this study was demonstrated in multiple countries and groundwater monitoring contexts. The use of the methodology made it possible to distinguish water families on the basis of lithology and geochemical conditions (pH, redox), for which different NBLs are expected. The proposed method can be used for the broad contexts investigated in this study; it provides coherent results. Therefore, it can be used as a step to explore a large dataset from groundwater monitoring networks. The methodology is an overall lithological approach; its results can be extrapolated to aquifers of comparable lithology by excluding local specificity links, for example, to the presence of mineralization, mining (metals and metalloids), or evaporites (gypsum). For each aquifer, this first estimation is open to the addition of new qualitative and hydrogeological data to refine the preliminary NBLs and the area of application.

The statistical tests applied depend on the amount of available data and the number of hydrogeological and lithological families. This approach may be limited for some hydrogeological entities (e.g., heterogeneous aquifers) depending on the number of available sampling points and characterization of hydrogeological settings (depth, confined/unconfined, pumped layer). NBLs determined for families characterized by only a few water points or that are not representative of the groundwater body should be viewed with caution. NBLs may be refined on the basis of additional data, the definition of lithology or groundwater entities, and the lithologic and petrographic description of the sampled aquifer. Hydrogeological settings for each groundwater need to be carefully addressed.

Determination of NBLs in alluvial aquifers is challenging because of obvious contamination and presence of anthropic activities. This caution is also relevant for unconsolidated and unprotected shallow aquifers, such as the aeolian or glacio-fluvial aquifers in Denmark [START_REF] Hansen | Nitrate vulnerability assessment of aquifers[END_REF]. In general, these vulnerable aquifers are usually impacted by nitrates or/and organic pollutants according to prevailing pressures [START_REF] Lopez | Screening of French groundwater for regulated and emerging contaminants[END_REF]. Conversely, confined aquifers are naturally protected from surface infiltration and, as a result, anthropogenic contamination is less common.

Using monitoring network data and comparing measured concentrations to calculated NBLs, it is possible to identify groundwater bodies for additional investigation to identify the origin of anomalous concentrations or anthropogenic inputs. Delineation of these entities and the interpretation of exceedances are beyond the scope of this study, because they require local investigation. Specific geogenic features (mineralization, evaporites, thermal waters) require a different methodology from the lithological NBL as defined in this study. The results of each prevailing pressure assessment cannot be generalized to other groundwater bodies, because trace element contaminations are directly linked to the type of anthropic activity; only general trends for each element can be highlighted.

This work constitutes a first step in the definition of NBLs. At this scale, data analytical methods might also be accompanied by process-oriented analyses and supported by expert background knowledge at the local scale. This methodology will depend strongly on the selection of relevant water points and the use of sampling and analytical methods suitable for reducing uncertainty in NBL characterization. To provide better estimates of trace element NBLs in the future, monitoring networks should be improved or additional sites besides the national monitoring network should be used with the selection of water points less subject to anthropogenic impacts to provide additional information regarding trace element concentrations. In addition, targeted analytical methods should be used to obtain lower LOQs and thereby provide more useful datasets.

Supplementary Materials:

The following are available online at https://www.mdpi.com/article/ 10.3390/w13111531/s1, Table S1. Denmark: The number of chemical analyses for each element (n), account of LODs in the dataset, range of LODs, the number and percentage of analyses < LOD, the max substitute value, and the number of analyses with the max substituted. Table S2. Serbia: The number of chemical analyses for each element (n), account of LODs in the dataset, the percentage of analyses < LOD, the max substitute value, and max LOD substituted. Table S3. Slovenia: The number of chemical analyses for each element, account of LOQs in the dataset (two were reported in 2016), range of LOQs, the number and percentage of analyses < LOQ, the max substitute value, and the number of analyses with the max substitute. Table S4. Loire-Bretagne, France: The number of chemical analyses for each element (n), account of value > LOQs in the dataset, range of LOQs (low, high), percentage of analyses < LOQ, percentage with high LOQ, and the maximum value substituted. Table S5. Internal Basins, Catalonia: The number of analyses for each element (n) and number and percentage of analyses < LOD, range of LOD, and maximum value imputed with each method. Table S6. Austria: Dataset overview, ranked by percentage of values below LOD or LOQ. Table S7. Spain, DRB: The number and percentage of chemical analyses below the limit of quantification (LOQ),

Figure 1 .

 1 Figure 1. Flowchart for NBL calculation applied in this study.

  Groundwater Quality Dataset Construction (Steps 1-9)

Figure 2 .

 2 Figure 2. Sampling points in Denmark classified on the basis of (a) aquifer type and (b) prevailing pressures.

Figure 3 .

 3 Figure 3. Example of 1 km buffer around sampling points. Abbreviations in DK model geology: kalk-carbonate aquifers; ks-quaternary sand, ps-pre-quaternary sand (figures from [56]).

Figure 2 .

 2 Figure 2. Sampling points in Denmark classified on the basis of (a) aquifer type and (b) prevailing pressures.

Figure 2 .

 2 Figure 2. Sampling points in Denmark classified on the basis of (a) aquifer type and (b) prevailing pressures.

Figure 3 .

 3 Figure 3. Example of 1 km buffer around sampling points. Abbreviations in DK model geology: kalk-carbonate aquifers; ks-quaternary sand, ps-pre-quaternary sand (figures from [56]).

Figure 3 .

 3 Figure 3. Example of 1 km buffer around sampling points. Abbreviations in DK model geology: kalk-carbonate aquifers; ks-quaternary sand, ps-pre-quaternary sand (figures from [56]).

  Water 2021, 13, x FOR PEER REVIEW 12 of 22

Figure 4 .

 4 Figure 4. Results for the Kruskal-Wallis test and the Nemenyi post hoc test for lithology families in the Loire-Bretagne area (number of samples in brackets). The categories sharing letters (above the boxplot) are not significantly different at the 95% confidence level based on the post hoc Nemenyi test [56].

Figure 4 .

 4 Figure 4. Results for the Kruskal-Wallis test and the Nemenyi post hoc test for lithology families in the Loire-Bretagne area (number of samples in brackets). The categories sharing letters (above the boxplot) are not significantly different at the 95% confidence level based on the post hoc Nemenyi test [56].

  class, as illustrated in Figure5for sedimentary aquifers, whereas pH is not discriminant.

Figure 4 .

 4 Figure 4. Results for the Kruskal-Wallis test and the Nemenyi post hoc test for lithology families in the Loire-Bretagne area (number of samples in brackets). The categories sharing letters (above the boxplot) are not significantly different at the 95% confidence level based on the post hoc Nemenyi test [56].

Figure 5 .

 5 Figure 5. Results for Kruskal-Wallis test and Nemenyi post hoc test for lithology families according to pH and redox classes in the Loire-Bretagne area (number of samples in brackets). The categories sharing letters (above the boxplot) are not significantly different at the 95% confidence level based on the post hoc Nemenyi test [56].

Figure 5 .

 5 Figure 5. Results for Kruskal-Wallis test and Nemenyi post hoc test for lithology families according to pH and redox classes in the Loire-Bretagne area (number of samples in brackets). The categories sharing letters (above the are not significantly different at the 95% confidence level based on the post hoc Nemenyi test [56].

(µg/L) Austria, France, Spain, Slovenia Denmark Serbia c

	As	10 (5) a	10	5	10
	Cd	5	5	0.5	3
	Cr	50	50	25	50
	Cu	2000	2000	100	2000
	Ni	20	20	10	20
	Zn	5000 b	not defined (5000-France)	100	3000
	F	1500	1500	1500	1200
	SO 4	250,000	250,000	250,000
	a Danish legislation for drinking water supply, LBK nr 118, 22 February 2018, and Ministerial Order BEK nr
	1070, 28 October 2019;

b

French decree of 11 January 2007; c Official Gazette of FRY volume 42/98 and 44/99 and Official Gazette of RS. Regulation on the hygienic acceptability of potable water.

Table 2 . Dataset descriptions and methodologies applied by area. Simplified lithology refers to the lithology classes de- fined in this study 1 . Area Reference Period Sampling Points/A nalyses Land Use and Prevailing Pressure (Methodology) Lithologies (Methodology and Classes) LOQ Treatment (Methodology) Data Source

 2 Flowchart for NBL calculation applied in this study.

				Watershed approach
	Loire-Bretagne basin, France (155,000 km 2 )	2009-2020	4200/7855 9	based on "TAUDEM" module, CLC 2012 [48] completed by inventories (agricultural uses: Agreste; industrial loca-SOL; mining inventory: tion: ICPE, BASIAS, BA-	Sampling point link to morphic, and volcanic BDLisa (1:50 000) simplified lithology: sedi-sand), crystalline, meta-mentary:(carbonate, gravel,	Substituted value LOQ removed. with LOQ/2-high	ADES
				SIG Mines)
	Nationwide Denmark (42,933 km 2 )	2009-2018	6388/1251 06	Buffer of 1 km, pressures extracted from CLC 2012 [48]	Sampling point link to the groundwater bodies delin-model [49]. eated based on the DK	Values < LOD sub-stituted with LOQ/2, removed 3*LOD. High LOD where LOQ =	JUPITER

Table 2 .

 2 Dataset descriptions and methodologies applied by area. Simplified lithology refers to the lithology classes defined in this study 1 .

	Area

Reference Period Sampling Points/Analyses Land Use and Prevailing Pressure (Methodology) Lithologies (Methodology and Classes) LOQ Treatment (Methodology) Data Source

	Water 2021, 13, 1531						6 of 21
				Watershed approach based on "TAUDEM" module, CLC 2012 [48] completed by	Sampling point link to BDLisa (1:50 000) simplified lithology:
	Loire-Bretagne basin, France (155,000 km 2 )	2009-2020	4200/78559	Agreste; industrial inventories (agricultural uses:	gravel, sand), sedimen-tary:(carbonate,	LOQ removed. Substituted value with LOQ/2-high	ADES
				location: ICPE,	crystalline,
				BASIAS, BASOL; mining inventory: SIG Mines)	metamorphic, and volcanic
					Sampling point link
					to the groundwater
	Nationwide Denmark (42,933 km 2 )	2009-2018	6388/125106	Buffer of 1 km, pressures extracted from CLC 2012 [48]	bodies delineated based on the DK model [49]. simplified lithology:	Values < LOD substituted with LOQ/2, where LOQ = 3*LOD. High	JUPITER
					sedimen-	LOD removed
					tary:(carbonate,
					sand)
					Derived from national
	Nationwide, Austria (83,879 km 2 )	2010-2020	2024/604353	Point data (10 m buffer) pressures extracted from CLC 2018 [50]	hydrogeological map (1:500.000) [51] sedimentary: gravel, simplified lithology:	Values < LOD or LOQ LOD/2 or LOQ/2 substituted with	WISA
					metamorphic
						"Log-Ratio EM
	Internal basins of Catalonia (Spain) (32,108 km 2 )	1994-2018 (long period due to low number of data available)	1336/8316	Buffer of 1 km, pressures extracted from CLC 2018 [50]	Aquifer map produced by the ACA (1: 50.000); simplified lithology: carbonate. sedimentary: crystalline rocks and	Algorithm" method with the lrEM function (zCompositions method applied by transformation Wilcoxon v1.3.4) [45] or the	SDIM-ACA
						calculating < LOD/2.

Table 2 .

 2 Cont.

	Area	Reference Period	Sampling Points/Analyses	Land Use and Prevailing Pressure (Methodology)	Lithologies (Methodology and Classes)	LOQ Treatment (Methodology)	Data Source
	Mountain Fruška Gora (Serbia)	2006-2017	33/132	Drastic vulnerability map, pressures 2006 [52,53] extracted from CLC	Simplified lithology: (quaternary deposits). sedimentary: sand	Values < LOD LOD/2 substituted with	Studies
					Lithostratigraphy
					according to Basic
					Geological Map
				Water body pressure	(1:100.000 and
	Nationwide, Slovenia (20,271 km 2 )	2016	203/2848	and land use from the Management Plan National Water	groundwater bodies. 1:250.000) [54] and points linked to 21	LOQ/2 Values < LOQ substituted with	ARSO
				2016-2021	Simplified lithology:
					sedimen-
					tary:(carbonate,
					sands)
	Spain, Duero River Basin (78,859 km 2 )	2010-2020	465/6457	Point data (1 km buffer) pressures extracted from CLC 2018 [50]	Simplified lithology: metamorphic and (Cenozoic) sedimentary: others	Values < LOQ LOQ/2 substituted with	IGME

1 Abbreviations: CLC-Corine land cover map (https://land.copernicus.eu/pan-european/corine-land-cover (accessed on 20 March 2021); LOQ-limit of quantification, LOD-limit of detection.

Table 3 .

 3 Redox and pH types applied.

	pH Type (n = 3)	Redox Types (n = 3)
	Acidic (pH < 7)	Oxic or anoxic (A, B redox types)
	Neutral (pH ∈ [7, 7.5])	Weakly or strongly reduced (C, D redox types)
	Basic (pH > 7.5)	Mixed (X redox type)

Table 4 .

 4 Definition of redox types (modified from[START_REF] Hansen | Kemisk Grundvandskortlaegning[END_REF]).

	REDOX TYPE	Redox Condition	NO 3 -mg/L	Fe mg/L	O 2 mg/L	SO 4 mg/L
	A	Oxic water	>1	<0.2	≥1	-
	B	Nitrate-reducing anoxic water	>1	<0.2	<1	-
	C	Weakly reduced water	≤1	≥0.2	-	≥20
	D	Strongly reduced	≤1	≥0.2	-	<20

Table 5 .

 5 NBLs for sedimentary aquifers. DK-Denmark, SR-Serbia, SI-Slovenia, FR-Loire-Bretagne, Cat-Internal Basins Catalonia, A-Austria, (n, number of sampling points). Values higher than the respective threshold value are in bold.

		DK			SR	SI		FR		DK		Cat		SI		FR		A		FR
		n	NBL	n	NBL	n	NBL	n	NBL	n	NBL	n	NBL	n	NBL	n	NBL	n	NBL	n	NBL
	As (µg/L)	3639	6.8	24	29		<0.5	343	1.46	1830	3.9	124	3.5		<0.5	554	3	1168	1.1	33	5
	Cd (µg/L)	259	0.01	24	<1	9	0.01	221	<0.5	77	0.03	_	_	51	<0.5	723	<0.5	135	0.1	723	<0.5
	Cr (µg/L)	186	0.2	13	20.3	9	2.4	5	<0.45	55	0.18	_	_	51	0.5	17	<0.4	673	0.6	<5	_
	Cu (µg/L)	218	1	15	4.6	9	1.3	39	4.4	67	1.3	229	8	51	0.6	131	3.1	126	2.2	9	1.4
	Ni (µg/L)	3661	1.8	14	6.3	9	3.3	231	2.3	1856	5.3	_	_	51	0.5	375	2.8	135	0.5	27	7.8
	Zn (µg/L)	271	8.1	19	41	9	380	38	12.3	85	20	235	99	51	4.5	133	24	_	_	10	65
	F (mg/L)	4046	0.4	24	0.3		<0.1	388	0.9	2213	1.2	_	_		<0.1	673	0.3	_	_	51	0.2
	SO 4 (mg/L)	3642	87	24	35.8	9	16	394	69.9	1834	89	8	60	47	8.2	686	44	_	_	52	56

Table 6 .

 6 NBLs for sedimentary: others (SI-quaternary sand; DRB-cenozoic), crystalline, metamorphic, and volcanic aquifers: SI-Slovenia, DRB-Duero River Basin (Spain), FR-Loire-Bretagne, Cat-Internal Basins Catalonia, A-Austria (n, number of sampling points).

		Sedimentary: Others		Crystalline Rocks				Metamorphic Rocks			Volcanic Rocks
		SI		DRB	Cat		FR			DRB	A		FR		FR
		n	NBL	n	NBL	n	NBL	n	NBL	n	NBL	n	NBL	n	NBL	n	NBL
	As (µg/L)	_	_	202	8.1	90	4.2	336	3.6	16	3.3	53	3.4	1321	6.7	373	2.4
	Cd (µg/L)	22	0.02	203	0.025	_	_	380	<0.5	16	<0.05	35	0.04	1437	<0.5	378	<0.5
	Cr (µg/L)	22	2	_	_	_	_	31	<0.5	_	_	50	0.5	123	<0.5	21	1
	Cu (µg/L)	22	0.7	203	2.7	237	10.8	25	9.6	16	< 5	34	0.5	95	11	13	<1
	Ni (µg/L)	22	15	_	_	_	_	211	5.1	_	_	35	0.5	741	5	139	0.3
	Zn (µg/L)	22	27	_	_	234	44.8	40	23.8	_	_	_	_	139	40.4	19	<5
	F (mg/L)	_	_	197	0.95	_	_	318	0.1	16	0.27	_	_	1146	0.1	310	0.2
	SO 4 (mg/L)	22	15.6	203	100.3	43	33.8	383	24	16	30.9	_	_	1469	18	392	7.3

Table 7 .

 7 Differences statistically observed for anthropogenically influenced values.

		SI	FR	Cat	DRB	A
	As		Significant differences between natural state and anthropogenic influences (agriculture, urban, industrial, mining)	Significant differences between natural state and anthropogenic influences (mining)		No significant differences between anthropogenic influences and natural state
	Cd	No significant differences between anthropogenic influences and natural state				Significant differences between natural state and anthropogenic influences (agricultural, industrial, and urban)
		Significant				Significant differences
		differences between				between natural state
	Cr	anthropogenic				and anthropogenic
		influences and				influences (industrial
		natural state				and urban)
	Cu	Significant differences between anthropogenic influences and natural state	Significant differences between natural state and anthropogenic influences (mining)	Significant differences between natural state and anthropogenic influences (mining)		Significant differences between natural state and anthropogenic influences (agricultural, industrial, and urban)
	Ni	Significant differences between anthropogenic influences and natural state	Significant differences between natural state and anthropogenic influences (agriculture)			Significant differences between natural state and anthropogenic influences (agricultural, industrial, and urban)
		Significant	Significant differences	No significant
		differences between	between natural state	differences between
	Zn	anthropogenic	and anthropogenic	anthropogenic
		influences and	influences (agriculture,	influences and
		natural state	urban, industrial)	natural state

Table 7 .

 7 Cont.

		SI	FR	Cat	DRB	A
	F		Significant differences between natural state and anthropogenic influences (urban)		Significant differences between anthropogenic influences and natural state
				Significant
	SO 4	Significant differences between anthropogenic influences and natural state	Significant differences between natural state and anthropogenic influences (urban)	differences between natural state and anthropogenic influences (agriculture, urban,	Significant differences between anthropogenic influences and natural state	Significant differences between anthropogenic influences and natural state
				industrial, mining)

Data Availability Statement: Final Report for HOVER WP3.3 is available here: https://repository. europe-geology.eu/egdidocs/hover/hover_33_statistical_data_treatment_backgroundlev.pdf. For France, monitoring data (ADES) are available https://ades.eaufrance.fr/ (accessed on 12 February 2021) and information for aquifer are available on https://bdlisa.eaufrance.fr/ (accessed on 15 September 2020). For Denmark, database JUPITER is available at https://eng.geus.dk/productsservices-facilities/data-and-maps/national-well-database-jupiter (accessed on 3 July 2019). Austrian data are available on the Austrian groundwater monitoring network (WISA): Wasserinformationssystem; Gewässerzustandsüberwachungsverordnung, BGBl. II Nr. 479/06 i.d.g.F. Catalonian data from SDIM-ACA are available on http://aca-web.gencat.cat/sdim21/ (accessed on 25 September 2020). Data from the Slovenian Environmental Agency are available on https://www.arso.gov.si/vode/ podatki/ (accessed on 15 July 2020) and land use on www.mkgp.gov and https://rkg.gov.si/vstop/ (accessed on access on 21 September 2020). Data from Duero River Basin (Spain) are available on https://www.chduero.es/web/guest/red-control-estado-quimico (accessed on 30 October 2020) and land use on http://centrodedescargas.cnig.es/CentroDescargas/catalogo.do?Serie=SIOSE (accessed on 29 December 2020).

range of LOQs in the dataset, and min and max substituted value. Table S8. Loire-Bretagne, France: The number of water sampling points with data for each element and each group of prevailing pressure. Table S9. Denmark: The number of water sampling points with data for each element and each group of prevailing pressure and results from the Kruskal-Wallis rank sum test. Table S10. Spain, DRB: The number of water sampling points with data for each element and each group of prevailing pressure and results from the Kruskal-Wallis rank sum test. Table S11. Austria: The confusion matrix for lithology. Table S12. Denmark: Natural background levels (NBLs) for sedimentary aquifers combined with pH and redox conditions. Table S13. Loire-Bretagne, France: Natural background levels (NBLs) for the sedimentary aquifers combined with pH and redox conditions.

Conflicts of Interest:

The authors declare no conflict of interest.