Combining natural background levels (NBLs) assessment with indicator kriging analysis to improve groundwater quality data interpretation and management

Combination of hydrochemical graphical methods and multivariate statistical analysis to delineate groundwater bodies

The identification of hydrogeological boundaries and the assessment of the quantitative and qualitative state of groundwater are necessary for the delineation of groundwater bodies, according to the European Guidelines. In this context, this study attempts to verify the current delineation of groundwater bodies (GWBs) focusing on their hydrochemical features through the application of hydrogeochemical methods and multivariate statistical analyses. The areas of interest are three adjacent GWBs located in the north-western sector of Campania Region in Southern Italy: the Volturno-Regi Lagni Plain, a coastal plain constituted of fluvial, pyroclastic and marine sediments; the Plain of Naples, an innermost plain of fluvial and pyroclastic sediments and the Phlegrean Fields, an active volcanic area with a series of monogenic volcanic edifices. Physicochemical data measured in groundwater samples collected at more than 200 sampling points were considered. Results reveal five different hydrogeochemical processes variably influencing the chemical features of the three GWBs: dissolution of carbonate rocks, influence of volcanic deposits and/or gas/hydrothermal fluid upwelling, salinisation, reducing conditions, and anthropogenic contamination. Combining hydrochemical diagrams and multivariate statistical analysis (i.e., factor analysis) allows depicting areas characterised by one or more hydrogeochemical processes, mostly reflecting known processes, but also highlighting the influence of groundwater flowpaths on water chemistry. The current delineation of the three GWBs should be revised considering that some hydrochemical features and processes are peculiar to each GWB, but others are in common between two or more GWBs. Uncertainties in the delineation of GWBs based on their hydrochemical features could be reduced by improving the distribution of the monitoring well network to ensure a more homogeneous coverage.

A regional investigation on natural background levels and contamination assessment of shallow groundwater contaminants in a coastal urbanized area

Knowledge on natural background levels (NBLs) is indispensable for evaluating groundwater contamination at regional scales. This study used a combination of modified oxidation capacity method and the box plot method to assess NBLs of common contaminants including iron, nitrite, sodium, total hardness, fluoride, sulfate, and selenium in groundwater of a coastal urbanized area of China (Pearl River Delta [PRD]), discussed geochemical factors controlling NBLs of these related contaminants in various groundwater units via multivariate statistical techniques, and evaluated groundwater comprehensive contamination in the PRD on the basis of groundwater NBLs. Results showed that NBLs of iron, nitrite, sodium, total hardness, fluoride, and sulfate in Groundwater Unit I of the PRD were higher than those in other groundwater units. Conversely, NBLs of selenium in Groundwater Units II and III of the PRD were higher than that in other groundwater units. Both of NBLs of iron and total hardness in Groundwater Unit I were higher than their allowable limits in groundwater quality standards recommended by China, while others were not. The high NBL of iron in Groundwater Unit I was probably attributed to the reductive dissolution of iron minerals in Quaternary sediments induced by the mineralization of organic matters in overlaid marine sediments. By contrast, the irrigation of river water and the leaching of trapped seawater in overlying marine sediment were likely responsible for the high NBL of total hardness in Groundwater Unit I. Using a groundwater contamination evaluation method on the basis of NBLs shows that uncontaminated groundwater accounted for 17.6%, 32.4%, 45.9%, and 53.6% in urban lands, peri-urban lands, agricultural lands, and other lands of the PRD, respectively. Groundwater with moderate to extremely high contamination was mainly distributed in Foshan, Dongguan, southern part of Guangzhou, western part of Shenzhen, northern part of Zhongshan, southern part of Zhuhai, and downstream of these cities. The leakage of industrial wastewater and domestic sewage and the irrigation of contaminated river water were likely main factors controlling groundwater comprehensive contamination in the PRD. PRACTITIONER POINTS: Groundwater natural background levels (NBLs) of seven contaminants in a coastal urbanized area (Pearl River Delta) were assessed. Geochemical factors controlling groundwater NBLs of seven contaminants in the Pearl River Delta were discussed. Groundwater comprehensive contamination in the Pearl River Delta was assessed on the basis of groundwater NBLs.

Natural background levels and driving factors of aluminum in shallow groundwater of an urbanized delta: Insight from eliminating anthropogenic-impacted groundwaters

Knowledge on natural background levels (NBLs) of aluminum (Al) in groundwater can accurately assess groundwater Al contamination at a regional scale. However, it has received little attention. This study used a combination of preselection and statistic methods consisting of the oxidation capacity and the boxplot iteration methods to evaluate the NBL of shallow groundwater Al in four groundwater units of the Pearl River Delta (PRD) via eliminating anthropogenic-impacted groundwaters and to discuss driving factors controlling high NBLs of Al in groundwater in this area. A total of 280 water samples were collected, and 18 physico-chemical parameters including Redox potential, dissolved oxygen, pH, total dissolved solids, HCO3 -, NH4 +, NO3 -, SO4 2-, Cl-, NO2 -, F-, K+, Na+, Ca2+, Mg2+, Fe, Mn, and Al were analyzed. Results showed that groundwater Al NBLs in groundwater units A-D were 0.11, 0.16, 0.15, and 0.08 mg/L, respectively. The used method in this study is acceptable for the assessment of groundwater Al NBLs in the PRD, because groundwater Al concentrations in various groundwater units in residual datasets were independent of land-use types, but they were opposite in the original datasets. The dissolution of Al-rich minerals in sediments/rocks was the major source for groundwater Al NBLs in the PRD, and the interaction with Al-rich river water was secondary one. The high groundwater Al NBL in groundwater unit B was mainly attributed to the acid precipitation and the organic matter mineralization inducing the release of Al in Quaternary sediments. By contrast, the high groundwater Al NBL in groundwater unit C mainly was ascribed to the release of Al complexes such as fluoroaluminate from rocks/soils into groundwater induced by acid precipitation, but it was limited by the dissolution of Mg minerals (e.g., dolomite) in aquifers. This study provides not only useful groundwater Al NBLs for the evaluation of groundwater Al contamination but also a reference for understanding the natural geochemical factors controlling groundwater Al in urbanized deltas such as the PRD. PRACTITIONER POINTS: The natural background level (NBL) of groundwater aluminum in the Pearl River Delta (PRD) was evaluated. The dissolution of aluminum-rich minerals in sediments/rocks was the major source for groundwater aluminum NBLs in the PRD. The acid precipitation and organic matter mineralization contribute to high groundwater Al NBL in the groundwater unit B. The acid precipitation contributes to high groundwater Al NBL in the groundwater unit C, while dissolution of magnesium minerals limits it.

Evaluating natural background levels of heavy metals in shallow groundwater of the Pearl River Delta via removal of contaminated groundwaters: Comparison of three preselection related methods

Assessing natural background levels (NBLs) in groundwater is a global concern. Knowledge on groundwater NBLs in urbanized areas is challenging due to the impact of complex human activities. Preselection related methods are common ones for assessing groundwater NBLs. The present study used three preselection related methods to assess groundwater heavy metals (lead, zinc, barium) NBLs in four groundwater units of the Pearl River Delta (PRD) where urbanization continues, and to identify the best one for assessing groundwater NBLs in urbanized areas. Here, methods include a preselection method (method-P), a preselection dominated method (method-PD), and a statistic dominated method (method-SD). Results showed that the method-PD was better than other two methods for assessing groundwater NBLs of heavy metals in the PRD. This is supported by the evidence that differences among heavy metals concentrations in various land-use types in residual datasets formed by the method-PD were insignificant. NBLs of lead in groundwater units I to IV assessed by the method-PD were 2.8 μg/L, 5.9 μg/L, 5.8 μg/L, and 2.6 μg/L, respectively. NBLs of zinc in groundwater units I to IV assessed by the method-PD were 30 μg/L, 180 μg/L, 160 μg/L, and 100 μg/L, respectively. NBLs of barium in groundwater units I to IV assessed by the method-PD were 120 μg/L, 120 μg/L, 90 μg/L, and 50 μg/L, respectively. Compared to the method-PD, the method-SD often underestimates groundwater NBLs of heavy metals because of using the experiential evaluation for residual datasets. The method-P also has an inaccurate evaluation of groundwater NBLs of heavy metals in comparison with the method-PD, owing to both of using the experiential evaluation and the absence of a function for outliers test. The method-P combining with an outliers test would be better than itself for assessing groundwater NBLs. Therefore, the method-PD is the first choice to be recommended for assessing groundwater NBLs in urbanized areas such the PRD. However, this method should not be taken into account for assessing groundwater NBLs in areas where groundwater Cl/Br mass ratios are invalid. Instead, the method-SD and the method-P combining with one outliers test may be choices, because no constraint for these two methods.

Assessing natural background levels of geogenic contaminants in groundwater of an urbanized delta through removal of groundwaters impacted by anthropogenic inputs: New insights into driving factors

Knowledge on driving forces controlling natural background levels (NBLs) of geogenic contaminants (GCs) in groundwater of coastal urbanized areas are still limited because of complex hydrogeological conditions and anthropogenic activities. This study assesses NBLs of two GCs including arsenic (As) and manganese (Mn) in four groundwater units of the Pearl River Delta (PRD) with large scale urbanization by using a preselection method composed of the chloride/bromide mass ratio versus chloride concentration and the oxidation capacity with the combination of Grubbs' test. More importantly, driving factors controlling NBLs of As/Mn in groundwater of the PRD are discussed. Results showed that groundwater As/Mn concentrations in residual datasets were independent of land-use types, while those in original datasets in different land-use types were distinct because of various human activities, indicating that the used preselection method in this study is valid for NBLs-As/Mn assessment in groundwater of the PRD. NBL-As in coastal-alluvial aquifers was >6 times that in other groundwater units. NBL-Mn in coastal-alluvial aquifers was 1.4 times that in alluvial-proluvial aquifers, and both were >4 times that in other two groundwater units. High NBLs-As/Mn in coastal-alluvial aquifers is mainly attributed to reduction of FeMn oxyhydr(oxides) induced by mineralization of organic matter in Quaternary sediments. Elevated pH also contributes higher NBL-As in coastal-alluvial aquifers. By contrast, higher NBL-Mn in alluvial-proluvial aquifers than in other two groundwater units mainly ascribes to reduction of FeMn oxyhydr(oxides) in Quaternary sediments triggered by irrigation of reducing river waters. In addition, more occurrence of As/Mn-rich sediments and the infiltration of As/Mn-rich river water are also important factors for high NBLs-As/Mn in coastal-alluvial aquifers. This study shows that revealing natural driving factors of GCs-rich groundwater in coastal urbanized areas on the basis of identification of contaminated groundwaters via the used preselection methods is acceptable.

BRIDGE methodology-based quality standards to assess aquifer chemical status in the southwest Bengal Basin, Bangladesh

Assessment of natural background levels (NBLs) of compositional groundwater parameters helps to identify the potential threats to groundwater resources. This study is the first attempt to apply the pre-selection-based BRIDGE (Background cRiteria for the IDentification of Groundwater thrEshold) methodology to calculate the NBLs and threshold values (TVs) of major groundwater constituents in the southwest Bengal Basin, Bangladesh. A database consisting of 78 groundwater samples was used to assess the NBLs and associated TVs of the major groundwater parameters (EC, Ca²⁺, Mg²⁺, Na⁺, K⁺, Cl^-, NO₃^-, SO₄^2-, PO₄^3-, Mn²⁺, and Fe²⁺). NBLs were derived based on 90th and 97.7th percentiles. The status of regional groundwater resources was assessed by applying 90th percentile NBL on a regional dataset (n = 196). Results revealed the "poor" chemical status of shallow aquifers denoting heavy deterioration of the groundwater quality due to anthropogenic interventions. Nitrate contamination and salinization were identified as the major threats to the deep groundwater of the southwest Bengal Basin. Finally, to verify the chemical status of groundwater in a heavily urbanized area, derived TVs were applied throughout the experimental site Khulna. Twenty-five deep groundwater samples were collected for this purpose. Though most of the parameters exhibited "good" chemical status, nitrate demonstrated anthropogenic groundwater contamination in Khulna City. Thus, the developed TVs would provide an early warning system of pollution. On a national scale, it is expected to facilitate the sustainable groundwater management of the country and contribute to achieving the Sustainable Development Goals (SDG) of the United Nations (UN) in Bangladesh.

A double pre-selection method for natural background levels assessment in coastal groundwater bodies

To evaluate the chemical status of groundwater bodies (GWB) according to the European Groundwater Directive, EU Member States are required to take into account natural background levels (NBLs) where needed. Assessing the NBLs in coastal GWBs is complicated by seawater intrusion which can be amplified by groundwater withdrawals increasing the salinization of such groundwater systems. This paper proposes a new method for the NBLs assessment in coastal areas based on a double pre-selection (PS) with fixed/dynamic limits. A case study in the Apulia region, located in southeastern Italy, is proposed, where we investigated four adjacent GWBs which form the complex karst, fractured Murgia aquifer, hosted in the Jurassic-Cretaceous carbonate platform, bounded by two seas and sustained by saltwater of marine intrusion in the coastal areas. Data related to 139 monitoring stations (MSs) of the regional groundwater monitoring network were used. The first PS, "static", based on a fixed limit of anthropogenic contamination markers (NO₃ and NH₄), allows for the elimination of MSs impacted by human activities. On these, the second PS, "dynamic", based on the identification of Cl anomalous values, allows for the identification of additional MSs affected by saline contamination. The residual dataset of MSs was used for the definition of NBLs of Cl, SO₄, F and B. A statistical comparison with historical Cl observations finally allowed us to verify if the salinity of current groundwater is representative of pristine conditions. The calculated NBLs of salinity parameters are higher for the two coastal GWBs, with chloride values between 0.8 and 2 mg/L. Conversely, fluorides always show very low NBLs. The double PS approach seems more effective for NBLs calculation in coastal aquifers affected by saline contamination, where the use of a fixed Cl limit fails. It may respond to the international needs for a standardized procedure for NBL assessment.

Synoptic Risk Assessment of Groundwater Contamination from Landfills

Waste management in Europe has improved in recent years, reducing the amount of waste disposed at landfills. However, there are still many landfills in the countries. It is well known that landfills that do not have measures in place to control leachate entering groundwater can contaminate groundwater long after the landfill is closed. Collecting monitoring results from all landfills allows permitting and management agencies to improve action plans. This relies on a synoptic risk assessment that allows prioritization and milestones to be set for required actions. The developed method of synoptic risk assessment is based on a conceptual model of the landfill and the results of chemical groundwater monitoring tested at 69 landfills in Slovenia. The study confirms that most landfills have a direct or indirect impact on groundwater quality. All landfills were classified into three priority classes on the basis of the synoptic risk assessment. The results show that a total of 24 landfills have a clearly pronounced impact on groundwater. A total of 31 landfills have a less pronounced impact due to the favorable natural attenuation capacity of the soil or the technically appropriate design of the landfill itself. A total of 14 landfills have a less pronounced or negligible impact on groundwater.

Prediction of tempo-spatial patterns and exceedance probabilities of atmospheric corrosion of Q235 carbon steel across China

To reduce the losses caused by the atmospheric corrosion of carbon steels, it is important to establish a prediction model to determine the corrosion rate of carbon steels in natural environments. In this study, a prediction model of atmospheric corrosion of Q235 carbon steel (PMACC-Q235) in China was established by coupling the mean impact value algorithm and back propagation artificial neural network. Tempo-spatial patterns of corrosion rates in five long-exposure time categories across China were analyzed. Ten main factors affecting the atmospheric corrosion of Q235 were identified. The corrosion rates in a single year were similar (approximately 30 μm/a) and larger than those for 2 (25.30 μm/a) and 3 years (21.66 μm/a). The spatial corrosion rates in the northwestern areas were primarily lower than those in southeastern coastal areas. This could be influenced by climatic factors, such as temperature, humidity, and precipitation. All corrosion rates reached the C2 level (>1.3 μm/a), and there was some possibility that they reached higher corrosion levels. The largest probability for the C3 level in all periods was an average of 0.91, and that for the C4 level was 0.83. Spatially, higher probabilities were mainly located in the southern area, especially in Hainan, located in the south and surrounded by sea. Corrosion rates largely varied among climatic zones, and mean corrosion rates in the tropical monsoon climate zone were the largest (average of three periods 33.39 μm/a). SO₂ and soluble-dust fall had the largest impact on the variations in the corrosion rates among different climatic zones.

Natural background levels in groundwater in the Pearl River Delta after the rapid expansion of urbanization: A new pre-selection method

Establishment of natural background levels (NBL) of groundwater in urbanized areas such as the Pearl River Delta (PRD) is challenging. Pre-selection methods are the most common approaches for NBL assessment, but it will overestimate (or underestimate) contaminated groundwater in urbanized areas by using present pre-selection methods with empirical definite values because of complicated human activities. Unlike present pre-selection methods, this study aims to establish a new pre-selection method with the indicative of Cl/Br ratios to identify contaminated groundwaters with convincing evidences. Specifically, this new method consists of indicatives of the oxidation capacity and the Cl/Br ratio combining with contaminated-markers. In addition, factors controlling NBL of Cl and NO₃ in groundwater in various hydrogeological units in the PRD were also discussed. Main procedures of this new method: contaminated-markers in various hydrogeological units are extracted by a hierarchical cluster analysis, thereby determining threshold values of Cl/Br ratios and Cl concentration in various hydrogeological units for identifying contaminated groundwater; After that, groundwater chemical datasets was selected by the oxidation capacity, and then tested by Grubbs' test until normal distributions. Groundwater Cl and NO₃ concentrations in datasets before and after this new method are dependent and independent of urbanization levels, respectively, indicating that the new method is useful for groundwater NBL assessment in urbanized areas such as the PRD. Both the seawater intrusion and the diffusion of Cl from marine deposits are likely to be responsible for the much higher NBL-Cl in coastal-alluvial and marine aquifers than in other hydrogeological units. Groundwater Cl enrichment resulted from groundwater recharge and evaporation is mainly responsible for the higher NBL-Cl in alluvial-proluvial aquifers than in lacustrine aquifers, fissured aquifers, as well as karst aquifers. More than double times NBL-NO₃ in alluvial-proluvial and fissured aquifers than in other hydrogeological units is probably attributed to more oxidizing conditions of their vadose zones and groundwaters.

A multivariate non-parametric approach for estimating probability of exceeding the local natural background level of arsenic in the aquifers of Calabria region (Southern Italy)

The concept of natural background level (NBL) aims at distinguishing the natural and anthropogenic contributions to concentrations of specific contaminants, as groundwater management and protection tools. This is usually defined as a unique value at a regional scale, even when the hydrogeological and geochemical features of a certain territory are far from homogeneous. The concentration of target contaminants is affected by multiple hydrogeochemical processes. This is the case of arsenic in the Calabria region, where concentrations are definitely variable in groundwater. To overcome the limitation of a traditional approach and to include the intrinsic hydrogeological and geochemical heterogeneity into the definition of the natural contribution to As content in groundwater, an integrated probabilistic approach to the NBL assessment combining aquifer-based preselection criteria and multivariate non-parametric geostatistics was proposed. In detail, different NBL values were selected, based on the aquifer type and/or hydrogeochemical features. Then, these aquifer-based NBL values of arsenic were used in the Probability Kriging method to map the probability of exceedance and to provide contamination risk management tools. This multivariate geostatistical approach that takes advantage of the physico-chemical variables used in the aquifer-based NBL values definition allowed mapping the probability of exceedance of As in a physically-based way. The hydrogeochemical diversity of the study area and all the processes affecting As concentrations in the aquifers have been considered too. As a result, the obtained map was characterized by a short-range and long-range variability due to local hydrogeochemical anomalies and water-rock interaction and/or atmospheric precipitation. By this approach, the NBL exceedance probability maps proved to be less "noisy", because the local hydrogeochemical conditions were filtered, and more capable of pointing out anthropogenic inputs or very anomalous natural contributions, which need to be investigated more in detail and properly managed.

Linking local natural background levels in groundwater to their generating hydrogeochemical processes in Quaternary alluvial aquifers

Calculating natural background levels (NBLs) in groundwater is vital for supporting a sustainable use of groundwater resources. Although NBLs are often assessed through a unique concentration value per groundwater body, where hydrogeochemical features are highly variable, spatial heterogeneity needs to be accounted for, leading to the calculation of so-called "local" NBLs. Despite much research devoted to the identification of the best performing techniques for local NBLs spatialization, a deep understanding of the link between local NBL values and their generating hydrogeochemical processes is often lacking and so is addressed here for the redox-sensitive species As, NH₄, Fe and Mn in the groundwater bodies of Lombardy region, N Italy. Local NBLs were calculated by a tired approach involving the hybridization of preselection and probability plot methods. Since the spatial variability of the target species depends mainly on redox conditions, a redox zonation was performed using multivariate statistical analysis. A conceptual model was developed and improved combing factor and cluster analysis. Results showed that NBLs for arsenic were up to 291 μg/L, reached in groundwaters under methanogenesis, a condition related to the prolonged degradation of peat buried in aquifer sediments. Ammonium NBLs up to 6.62 mg/L were generated by the upwelling of fluids from deep sediments hosting petroleum systems; ammonium NBLs up to 4.48 mg/L were generated as the accumulation of by-products of peat degradation. Iron and manganese NBLs up to, respectively, 6.0 and 1.51 mg/L were generated by the oxidation of younger and less stable Mn and Fe oxides within river valleys, mostly the Po River valley. The evaluation of local NBLs, and their association to generating natural hydrogeochemical processes/conditions, achieves a step forward from the commonly used approach of a single NBL per groundwater body, improving decision-support tools for sustainable groundwater management and protection.

Assessment of the Impact of Climate Extremes on the Groundwater of Eastern Croatia

This paper analyzes the groundwater in the deep Quaternary aquifer of Eastern Croatia. These waters are collected at the Vinogradi Pumping Station (Osijek, Croatia) for the needs of public water supply. This research aimed to assess the impact of climate extremes, namely, high air temperatures and low rainfall, on the quantity and quality of groundwater. On the basis of data from the Vinogradi Pumping Station in the period 1987–2015, three extremely warm and low-water years were singled out. For these three years, the following were analyzed: climate diagrams, groundwater levels (in the piezometers closest to and farthest from the pumping station), and the quality of the affected groundwater. The results of this research indicate that the reaction of aquifers to the analyzed extreme climatic conditions for the observed period was manifested in the variation of the amplitude of groundwater levels by a maximum of 4–5 m. Considering the total thickness of the affected layers (60–80 m), this variation is not a concern from the point of view of water supply. As for the quality of groundwater, it was found to be of constant quality in its composition and was not affected by climatic extremes.

Source apportionment and natural background levels of major ions in shallow groundwater using multivariate statistical method: A case study in Huaibei Plain, China

Understanding the sources, natural background levels (NBLs), and threshold values (TVs) of the major ions in groundwater is essential for the effective protection of water resources. In this study, a total of 70 shallow groundwater samples were collected in Suzhou, Huaibei Plain, China. A variety of statistical methods and cumulative probability distribution techniques were performed to identify the sources, NBLs, and TVs of the major ions. The major ion concentrations found in decreasing order as follows: HCO₃^- > SO₄^2- > NO₃^- > Cl^- and Na⁺ > Ca²⁺ > Mg²⁺. Piper diagram for hydrochemical types shows that groundwater types were Mg-HCO₃ (36%), Ca-HCO₃ (34%), and Na-HCO₃ (30%). According to the factor and the Unmix model analysis, anthropogenic (agriculture-related) and geogenic source (water-rock interactions-related) were identified to be responsible for the chemical composition of the groundwater in the study area, and their mean contributions for the major ion concentrations are 47.9% and 52.1%, respectively. The NBLs for Na⁺, Ca²⁺, Mg²⁺, Cl^-, SO₄^2-, and NO₃^- were determined to be 29.5-44.2, 26.2-38.9, 18.9-39.5, 1.0-9.9, 12.9-19.4, and 2.1-16.5 mg/L, respectively, and the TVs were calculated as 122.1, 169.5, 39.5, 129.6, 134.7, and 18.3 mg/L, respectively. Moreover, this study shows the feasibility and reliability of using these multivariate statistical methods and natural background levels to evaluate the status of groundwater quality.

Identification of Groundwater Contamination in a Rapidly Urbanized Area on a Regional Scale: A New Approach of Multi-Hydrochemical Evidences

Efficient identification of groundwater contamination is a major issue in the context of groundwater use and protection. This study used a new approach of multi-hydrochemical indicators, including the Cl-Br mass ratio, the hydrochemical facies, and the concentrations of nitrate, phosphate, organic contaminants, and Pb in groundwater to identify groundwater contamination in the Pearl River Delta (PRD) where there is large scale urbanization. In addition, the main factors resulting in groundwater contamination in the PRD were also discussed by using socioeconomic data and principal component analysis. Approximately 60% of groundwater sites in the PRD were identified to be contaminated according to the above six indicators. Contaminated groundwaters commonly occur in porous and fissured aquifers but rarely in karst aquifers. Groundwater contamination in porous aquifers is positively correlated with the urbanization level. Similarly, in fissured aquifers, the proportions of contaminated groundwater in urbanized and peri-urban areas were approximately two times that in non-urbanized areas. Groundwater contamination in the PRD was mainly attributed to the infiltration of wastewater from township-village enterprises on a regional scale. In addition, livestock waste was also an important source of groundwater contamination in the PRD. Therefore, in the future, the supervision of the wastewater discharge of township-village enterprises and the waste discharge of livestock should be strengthened to protect against groundwater contamination in the PRD.

Monitoring spatiotemporal variation of groundwater level and salinity under land use change using integrated field measurements, GIS, geostatistical, and remote-sensing approach: case study of the Feija aquifer, Middle Draa watershed, Moroccan Sahara

The cultivation of watermelons has been a fast growing agriculture industry in the arid, desert regions of Morocco, relying on groundwater pumping and transformation of rangelands to farms due to growing demand for the fruit in national and international markets. This study aims to measure the impact of watermelon expansion on groundwater resources in the Feija Basin, which is one of the largest watermelon cultivation areas in Southern Morocco. Field measurements, statistics, Kriging interpolation, and regression methods were used to measure the temporal variations in the groundwater level (GL) and salinity between 2013 and 2018 to determine the correlation between different parameters. Remote sensing data was also used to monitor the watermelon cultivation expansion. Results show a rapid expansion of agricultural areas from just 185.11 ha in 2007 to 2560.1 ha in 2018. The groundwater level declined rapidly by about 10 m below ground level during the 5 years of the study period. Additionally, the decline was accompanied by a significant increase in electrical conductivity (salinity) values over the same time interval from 1077.55 to 1211.9 µS/cm. As a consequence of the continuous overexploitation and unsustainable management, a lot of wells have run dry and there have been drinking water shortages in the city of Zagora, the closest city nearby. Results can help target efforts to improve the implementation of conservation strategies to ensure the sustainability of water use and food production in this region of Morocco.

Hydrogeochemical Processes and Natural Background Levels of Chromium in an Ultramafic Environment. The Case Study of Vermio Mountain, Western Macedonia, Greece

The hydrogeochemical processes and natural background levels (NBLs) of chromium in the ultramafic environment of Vermio Mountain, Western Macedonia, Greece, were studied. Seventy groundwater samples were collected from 15 natural springs between 2014–2020, and an extensive set of physical and chemical parameters were determined. The ultramafic-dominated environment of western Vermio Mt. favors elevated groundwater concentrations of dissolved magnesium (Mg2+), silicon (Si), nickel (Ni), and Cr in natural spring waters. Chromium was the principal environmental parameter that exhibited a wide range of concentrations, from 0.5 to 131.5 μg/L, systematically exceeding the permissible limit of 50 μg/L for drinking water. Statistical evaluation of hydrogeological, hydrochemical, and hydrological data highlighted the water-ultramafic rock process as the predominant contributor of Cr in groundwater. The NBL assessment for Cr and Cr(VI) was successfully applied to the typical ultramafic-dominated spring “Potistis” that satisfied all the methodology criteria. The NBLs of Cr and Cr(VI) were defined at 130 μg/L and 100 μg/L, respectively, revealing that a natural ultramafic-dominated environment exhibits the geochemical potential to contribute very high concentrations of geogenic Cr to groundwater. The holistic methodology, proposed herein, could be implemented in any catchment scale to assess geogenic and anthropogenic Cr-sources that degrade groundwater quality.

Natural Background Level and Contamination of Shallow Groundwater Salinity in Various Aquifers in a Coastal Urbanized Area, South China

Assessing natural background levels (NBLs) of chemical components in groundwater is useful for the evaluation of groundwater contamination in urbanized areas. The present study assessed the NBL of total dissolved solids (TDS) in various groundwater units in the Pearl River Delta (PRD) where urbanization is a large scale and discussed factors controlling groundwater salinity contamination in the PRD. Results showed that the NBL of TDS in groundwater in the coastal-alluvial plain was more than 1.5 times that in other groundwater units because of the seawater intrusion in this groundwater unit. By contrast, interactions of water and soils/rocks were the main factors controlling the NBLs of TDS in other groundwater units. Groundwater salinity contamination in the PRD was positively correlated with the urbanization level. Wastewater from township-village enterprises and industrial wastewater were likely to be the main sources for groundwater salinity contamination in the PRD. Moreover, the wastewater leakage from sewer systems was one of the main pathways for groundwater salinity contamination in urbanized areas, because the proportion of groundwater salinity contamination in urbanized areas formed in 1988–1998 was more than 1.5 times that in urbanized areas formed in 1998–2006 regardless of groundwater units. Besides, sewage irrigation and leakage of landfill leachate were also important sources for groundwater salinity contamination in the PRD.

Assessing Natural Background Levels in the Groundwater Bodies of the Apulia Region (Southern Italy)

Defining natural background levels (NBL) of geochemical parameters in groundwater is a key element for establishing threshold values and assessing the environmental state of groundwater bodies (GWBs). In the Apulia region (Italy), carbonate sequences and clastic sediments host the 29 regional GWBs. In this study, we applied the Italian guidelines for the assessment of the NBLs, implementing the EU Water Framework Directive, in a south-European region characterized by the typical Mediterranean climatic and hydrologic features. Inorganic compounds were analyzed at GWB scale using groundwater quality data measured half-yearly from 1995 to 2018 in the regional groundwater monitoring network (341 wells and 20 springs). Nitrates, chloride, sulfate, boron, iron, manganese and sporadically fluorides, boron, selenium, arsenic, exceed the national standards, likely due to salt contamination along the coast, agricultural practices or natural reasons. Monitoring sites impacted by evident anthropic activities were excluded from the dataset prior to NBL calculation using a web-based software tool implemented to automate the procedure. The NBLs resulted larger than the law limits for iron, manganese, chlorides, and sulfates. This methodology is suitable to be applied in Mediterranean coastal areas with high anthropic impact and overexploitation of groundwater for agricultural needs. The NBL definition can be considered one of the pillars for sustainable and long-term groundwater management by tracing a clear boundary between natural and anthropic impacts.

Deriving Natural Background Levels of Arsenic at the Meso-Scale Using Site-Specific Datasets: An Unorthodox Method

Arsenic is found in groundwater above regulatory limits in many countries and its origin is often from natural sources, making the definition of Natural Background Levels (NBLs) crucial. NBL is commonly assessed based on either dedicated small-scale monitoring campaigns or large-scale national/regional groundwater monitoring networks that may not grab local-scale heterogeneities. An alternative method is represented by site-specific monitoring networks in contaminated/polluted sites under remediation. As a main drawback, groundwater quality at these sites is affected by human activities. This paper explores the potential for groundwater data from an assemblage of site-specific datasets of contaminated/polluted sites to define NBLs of arsenic (As) at the meso-scale (order of 1000 km2). Common procedures for the assessment of human influence cannot be applied to this type of dataset due to limited data homogeneity. Thus, an “unorthodox” method is applied involving the definition of a consistent working dataset followed by a statistical identification and critical analysis of the outliers. The study was conducted in a highly anthropized area (Ferrara, N Italy), where As concentrations often exceed national threshold limits in a shallow aquifer. The results show that site-specific datasets, if properly pre-treated, are an effective alternative for the derivation of NBLs when regional monitoring networks fail to catch local-scale variability.

eNaBLe, an On-Line Tool to Evaluate Natural Background Levels in Groundwater Bodies

Inorganic compounds in groundwater may derive from both natural processes and anthropogenic activities. The assessment of natural background levels (NBLs) is often useful to distinguish these sources. The approaches for the NBLs assessment can be classified as geochemical (e.g., the well-known pre-selection method) or statistical, the latter involving the application of statistical procedures to separate natural and anthropogenic populations. National Guidelines for the NBLs assessment in groundwater have been published in Italy (ISPRA 155/2017), based mainly on the pre-selection method. The Guidelines propose different assessment paths according to the sample size in spatial/temporal dimension and the type of the distribution of the pre-selected dataset, taking also into account the redox conditions of the groundwater body. The obtained NBLs are labelled with a different confidence level in function of number of total observations/monitoring sites, extension of groundwater body and aquifer type (confined or unconfined). To support the implementation of the Guidelines, the on-line tool evaluation of natural background levels (eNaBLe), written in PHP and using MySQL as DBMS (DataBase Management System), has been developed. The main goal of this paper is to describe the functioning of eNaBLe and test the tool on a case study in central Italy. We calculated the NBLs of As, F, Fe and Mn in the southern portion of the Mounts Vulsini groundwater body, within the volcanic province of Latium (Central Italy), also separating the reducing and oxidizing facies. Specific results aside, this study allowed to verify the functioning and possible improvements of the online tool and to identify some criticalities in the procedure NBLs assessment at the groundwater body scale

Probabilistic assessment of spatial heterogeneity of natural background concentrations in large-scale groundwater bodies through Functional Geostatistics

We propose and exemplify a framework to assess Natural Background Levels (NBLs) of target chemical species in large-scale groundwater bodies based on the context of Object Oriented Spatial Statistics. The approach enables one to fully exploit the richness of the information content embedded in the probability density function (PDF) of the variables of interest, as estimated from historical records of chemical observations. As such, the population of the entire distribution functions of NBL concentrations monitored across a network of monitoring boreholes across a given aquifer is considered as the object of the spatial analysis. Our approach starkly differs from previous studies which are mainly focused on the estimation of NBLs on the basis of the median or selected quantiles of chemical concentrations, thus resulting in information loss and limitations related to the need to invoke parametric assumptions to obtain further summary statistics in addition to those considered for the spatial analysis. Our work enables one to (i) assess spatial dependencies among observed PDFs of natural background concentrations, (ii) provide spatially distributed kriging predictions of NBLs, as well as (iii) yield a robust quantification of the ensuing uncertainty and probability of exceeding given threshold concentration values via stochastic simulation. We illustrate the approach by considering the (probabilistic) characterization of spatially variable NBLs of ammonium and arsenic detected at a monitoring network across a large scale confined groundwater body in Northern Italy.

An investigation into seasonal variations of groundwater nitrate by spatial modelling strategies at two levels by kriging and co-kriging models

Nitrate pollution of groundwater through spatial models is investigated in this paper by using a sample of nitrate values at monitoring wells using the data from four seasons of a year, in which data are sparse. Two spatial modelling strategies are formulated at two levels, in which Strategy 1 comprises: three variations of kriging-based models (ordinary kriging, simple kriging and universal kriging), which are constructed at Level 1 to predict nitrate concentrations; and a Multiple Co-Kriging (MCoK) model is used at Level 2 to enhance the accuracy of the predictions. Strategy 2 is also at two levels but employs Indicator Kriging (IK) at Level 1 as a probabilistic spatial model to predict areas at risk of exceeding two thresholds of 37.5 mg/L and 50 mg/L of nitrate concentration, and Multiple Co-Indicator Kriging (MCoIK) at Level 2 for a better accuracy. The improvements at Level 2 for both strategies are remarkable and hence they are used to gain an insight into inherent problems. The results of a study delineate areas with excessive nitrate concentrations, which are in the vicinity of urban areas and hence reflect poor planning practices since the 1990s. The results further reveal the patterns on sensitivities to seasonal variations driven by aquifer recharge and strong dilution processes in spring times; and on the role of pumpage impacting aquifers giving rise to possible hotspots of nitrate concentrations.

Combining geospatial analyses to optimize quality reference values of rare earth elements in soils

High pedological and geological variability can trigger the formation of REE hotspots, causing a need to optimize the establishment of quality reference values (QRVs). Thus, we determined the background concentrations of REEs in the soils of an emerging Brazilian state and used a combination of Moran's I and indicator kriging to identify REE hotspots and determine QRVs. A total of 100 composite soil samples was collected at a 0.20 m depth to establish background concentrations, QRVs, and spatial distribution and to elaborate probability maps for REEs. The QRVs established for soils were the following (mg kg^-1): La (27.21), Ce (57.26), Pr (10.49), Nd (24.29), Sm (4.75), Eu (0.90), Gd (4.22), Tb (0.82), Dy (1.54), Ho (0.38), Er (1.23), Yb (1.07), Lu (0.24), Y (10.65), and Sc (3.70). It was possible to draw attention to the Northwest and Southwest regions of the Rio Grande do Norte (RN) state, due to the formation of REE hotspots, indicated by Moran's I, and a high tendency to exceed the QRVs, confirmed by the indicator kriging. The high background concentrations and geochemical patterns for REEs showed that a single QRV for each REE and the entire state can neglect specific environmental characteristics and misrepresent the natural geochemistry of the soil. Thus, specific QRVs were established to optimize the monitoring of natural REE values by identifying hotspot areas. The criteria established here may be useful for other groups of potentially toxic elements, provided that observations meet the requirements of the spatial autocorrelation and kriging analyses. Graphical abstract.

Meteorological Variability and Groundwater Quality: Examples in Different Hydrogeological Settings

Rainfall and temperature variability causes changes in groundwater recharge that can also influence groundwater quality by different processes. The aim of this study is the analysis of the hydrogeochemical variations over time due to meteorological variability in two different study areas in Italy: an alluvial aquifer in the Piedmont Po plain and an alluvial-pyroclastic aquifer in the Campanian plain. The examined plains show groundwater with natural quality not satisfying the European drinking water standards, or anthropogenic contamination. The peculiar natural quality is due, in the Campanian plain, to the closeness of volcanic areas, and to the presence of reducing conditions. In Piedmont plain a test site is characterized by a point-source contamination by heavy metals, due to the presence of past industrial activities. In all the examined areas there is a diffuse nitrate contamination. The fluctuations of the ions As, F, Fe, Mn, Cr VI, NO3, and Cl were analyzed and compared, using statistical methods, with the variations over time in precipitation, temperature, and piezometric levels, sometimes significant. Results highlight the importance of the groundwater and meteorological monitoring and the key role of the recharge variation in the hydrogeochemical processes. The linking degree between rainfall/temperature variability and hydrogeochemistry is variable, in function of the typology of chemical species, their origin, and of the aquifer characteristics. The fluctuation of climate variables determines sudden changes in the geochemistry of shallow unconfined aquifers (e.g., in the Piedmont plain), while semiconfined or confined aquifers (e.g., in the Volturno-Regi Lagni plain) react with a greater delay to these variations. Moreover, natural quality is more affected by climatic variations than anthropogenic contamination, which is the result of multiple environmental and anthropic factors.

Advanced redox zonation of the San Pedro Sula alluvial aquifer (Honduras) using data fusion and multivariate geostatistics

The incorrect wastewater management and the land use distribution lead to severe environmental problems, creating heavy eutrophication condition in surface-water; when surface-water/groundwater relationships exist, the organic matter transferred to the aquifer oxidizes and triggers redox processes (i.e. Terminal Electron Accepting Processes, TEAPs), that provoke severe groundwater quality modifications and complicate its exploitation and management. For this reason, the definition of the redox zonation within an aquifer is an effective tool for the identification of the contamination sources and for the conceptual model refinement, when remediation strategies need to be planned. Although the redox processes are dynamic reactions, the aquifer redox zonation is generally aimed to identify homogenous zones, characterized by a predominant TEAP. To overcome this methodological approach, the Multi-Collocated Factorial Kriging (MCFK) was applied to redox-related physico-chemical parameters, which allowed identifying their spatial relationships at different scales, transferring this method from precision agriculture and soil science to hydrogeochemistry. The selected study area is the San Pedro Sula aquifer (Honduras), a multi-layer alluvial aquifer characterized by well-known surface-water/groundwater interactions and heavy eutrophicated streams. Here, high concentrations of Mn and Fe were found in the aquifer. The MCFK results identified a short-range (2300 m) factor, highlighting a strong relation between Mn concentrations and anoxic conditions, due to the organic matter transfer from eutrophicated surface-water into the aquifer. Simultaneously, the relationship between Fe and turbidity is related to a fine Fe(III) oxi-hydroxide colloidal phase, developed when different redox conditions of groundwater mix up in the wells. The long-range (6000 m) factor points out that Fe is related to redox processes at a wider scale, especially in the northern San Pedro Sula alluvial plain. These results are supported by both the Principal Component Analysis and the hydrogeochemical numerical modeling. As a result, different TEAPs occur simultaneously in contaminated areas, acting at multiple scales.

Groundwater Autochthonous Microbial Communities as Tracers of Anthropogenic Pressure Impacts: Example from a Municipal Waste Treatment Plant (Latium, Italy)

The groundwater behavior at a municipal solid waste disposal dump, located in Central Italy, was studied using a multi-parameter monitoring over 1 year consisting of 4 seasonal samples. The hydrological and hydrogeological dynamics of water circulation, microbiological parameters (microbial abundance and cell viability of the autochthonous microbial community), dissolved organic carbon, and several contaminants were evaluated and related to the geological structures in both two and three dimensions and used for geostatistical analysis in order to obtain 3D maps. Close relationships between geological heterogeneity, water circulation, pollutant diffusion, dissolved organic carbon, and cell viability were revealed. The highest cell viability values were found with dissolved organic carbon (DOC) values ≤0.5 mg/L; above this value, DOC negatively affected the microbial community. The highest DOC values were detected in groundwater at some sampling points within the site indicating its probable origin from the waste disposal dump. Although legislation limits for the parameters measured were not exceeded (except for a contaminant in one piezometer), the 1-year multi-parameter monitoring approach made it possible to depict both the dynamics and the complexity of the groundwater flux and, with “non-legislative parameters” such as microbial cell viability and DOC, identify the points with the highest vulnerability and their origin. This approach is useful for identifying the most vulnerable sites in a groundwater body.

A multi-method approach for the assessment of natural background levels in groundwater

The assessment of geochemical Natural Background Levels (NBLs) in groundwater, aims at distinguishing the naturally high levels of geogenic compounds from anthropogenic pollution. This is a fundamental issue in groundwater management, in particular when the concentration of inorganic compounds exceeds the threshold values set for the evaluation of the groundwater chemical status, as requested by environmental regulations. In this paper, we describe a new procedure that integrates the pre-selection method and statistical techniques, using the example of two case studies. The pre-selection aims to identify suitable groundwater samples for the NBLs assessment. The nitrate concentration threshold, for the removal of the groundwater samples affected by human activities, is established locally through different graphical and statistical approaches. Then, the statistical distribution of each compound is analyzed and the outliers are identified. Normality tests on the datasets allow one to select the most appropriate value, e.g. one percentile, to be adopted as NBL within the data distribution. In the selected case studies, we have defined the NBLs for As, F, Mn, Fe and SO₄. The two sites are part of a volcanic-sedimentary aquifer in central Italy, where the geochemical background is frequently well above the standards for human consumption. The results of the simple and easily reproducible pre-selection method are strengthened by integration with statistical techniques, notably in selecting the appropriate percentile. New criteria are suggested for the choice of the nitrate threshold to be used for the pre-selection of uncontaminated samples.

Spatial analysis of total dissolved solids in Dezful Aquifer: Comparison between universal and fixed rank kriging

The spatial structure and auto-correlation of total dissolved solids (TDS) in an aquifer located in southwestern part of Iran were investigated by both Moran's index and variography. Since the feature of interest was non-stationary so, conventional methods of spatial analysis were not applicable and Universal kriging (UK) as a common method for spatial prediction of features with a spatial trend along with a novel geostatistical method known as fixed rank kriging (FRK) were utilized in this respect. The results of Moran's index were consistent with that of spatial analysis by geostatistical methods indicating the dominance of spatial clusters within the extent of study area. The spatial analysis by FRK was more efficient than that of its UK counterpart however the performance of UK was reasonable enough, as well. A variable selection by random forest (RF) was applied on eleven other water quality parameters that were the main constituents of TDS to identify the main parameters influencing the observed variability of TDS. It was turned out that RF is a viable method for variable selection in the realm of environmental sciences.

Geostatistical multimodel approach for the assessment of the spatial distribution of natural background concentrations in large-scale groundwater bodies

Quantification of the (spatially distributed) natural contributions to the chemical signature of groundwater resources is an emerging issue in the context of competitive groundwater uses as well as water regulation and management frameworks. Here, we illustrate a geostatistically-based approach for the characterization of spatially variable Natural Background Levels (NBLs) of target chemical species in large-scale groundwater bodies yielding evaluations of local probabilities of exceedance of a given threshold concentration. The approach is exemplified by considering three selected groundwater bodies and focusing on the evaluation of NBLs of ammonium and arsenic, as detected from extensive time series of concentrations collected at monitoring boreholes. Our study is motivated by the observation that reliance on a unique NBL value as representative of the natural geochemical signature of a reservoir can mask the occurrence of localized areas linked to diverse strengths of geogenic contributions to the groundwater status. We start from the application of the typical Pre-Selection (PS) methodology to the scale of each observation borehole to identify local estimates of NBL values. The latter are subsequently subject to geostatistical analysis to obtain estimates of their spatial distribution and the associated uncertainty. A multimodel framework is employed to interpret available data. The impact of alternative variogram models on the resulting spatial distributions of NBLs is assessed through probabilistic weights based on model identification criteria. Our findings highlight that assessing possible impacts of anthropogenic activities on groundwater environments with the aim of designing targeted solutions to restore a good groundwater quality status should consider a probabilistic description of the spatial distribution of NBLs. The latter is useful to provide enhanced information upon which one can then build decision-making protocols embedding the quantification of the associated uncertainty.

Disentangling natural and anthropogenic impacts on groundwater by hydrogeochemical, isotopic and microbiological data: Hints from a municipal solid waste landfill

Within human-impacted areas, high levels of inorganic compounds in groundwater are broadly and preventively attributed to local anthropogenic pollution, thoroughly disregarding geogenic natural background levels. Particularly in landfills, a proper evaluation of the significant adverse environmental effects should be completed through a detailed groundwater characterization, and appropriate reference values established prior to landfill onset. However, the monitoring network may lack a full hydrogeological representativeness of the site and of the background conditions of groundwater. This study aimed at disentangling natural and anthropogenic impacts through a synoptic analysis of hydrogeochemical, isotopic and microbiological characteristics of groundwaters from a municipal solid waste landfill area in Central Italy. Samples were collected during four seasonal monitoring surveys from the mostly anoxic aquifer underlying the target area. Field parameters, inorganic and organic compounds, environmental isotopes, faecal contamination, and microbial community characteristics were determined, along with a detailed hydrogeological conceptual model. Key inorganic contaminants (As, Fe and Mn) exceeded the local threshold values in most of the sampling points, while organic contamination was generally very low. Stable isotopes suggested that groundwater originated mainly from local rainfall, except at one monitoring points where tritium levels might indicate moderate impact. Microbiological data and the microbial community characterization, assessed by flow cytometry and BIOLOG assays, provided further supportive information, also highlighting fundamental effects of groundwater quality alterations. Overall, an integrated multi-parametric approach proved suitable to distinguish geogenic and anthropogenic impacts, thus improving strategies and schemes for protection and management of groundwaters in landfills and waste related industrial areas.

Nitrate trends in groundwater of the Campania region (southern Italy)

The Environmental Protection Agency of the Campania region in Italy (ARPAC) manages a groundwater quality monitoring network. For almost all the polluted waters, the key parameter driving the classification is the concentration of nitrate; hence, the Campania region, in coherence with the EU regulations, outlined the vulnerable areas and undertook remediation policies. The best groundwater quality is recorded for carbonate aquifers of the Apennine chain; on the contrary, the Tyrrhenian coastal plains are affected by severe contamination, with a locally very contaminated groundwater of the shallow and also the deeper aquifers. The study is especially focused on a large coastal plain of Campania region, where nitrate concentration sometimes exceeds 200 mg/L. The study, based on almost 200 sampling points for the whole region during the period 2003-2015 (approx two samples per year), verified the effectiveness of the groundwater monitoring network, the present distribution of nitrate in groundwater, and the evolution of nitrate trends at different scales: regional, groundwater body, and single well, using spatial and time series statistical approaches. Significant variations in contamination evolution within the study area have been observed and the correlation with land use has been highlighted.

Spatio-Temporal Analysis of Natural and Anthropogenic Arsenic Sources in Groundwater Flow Systems

The presence of arsenic in groundwater constitutes a hazard for the environment and human health, and the determination of its source has become a global challenge, which can be approached by defining the natural background levels (NBL) in conjunction with the indicator kriging method, with the aim of delineating anthropogenically contaminated areas. However, having a unique value of NBL for large areas can generate interpretation errors. This research integrates the determination of the flow systems present in the Calera Aquifer, and the definition of the natural background levels in each flow system by making estimation maps in ArcGIS using two databases, 10 years apart, to evaluate the spatio-temporal variation of arsenic in groundwater. The results indicate a notable increase in the probability of exceeding the arsenic NBL, mainly in the intermediate flow, which may be due to movement resulting from mining activities as well as a mixture of regional and intermediate flows caused by the extraction of water for agriculture and drinking water supplies. The presented values exceed the maximum limits allowed for human consumption, as stated by the World Health Organization.

Determination of trigger levels for groundwater quality in landfills located in historically human-impacted areas

Landfills are one of the most recurrent sources of groundwater contamination worldwide. In order to limit their impacts on groundwater resources, current environmental regulations impose the adoption of proper measures for the protection of groundwater quality. For instance, in the EU member countries, the calculation of trigger levels for identifying significant adverse environmental effects on groundwater generated by landfills is required by the Landfill Directive 99/31/EC. Although the derivation of trigger levels could be relatively easy when groundwater quality data prior to the construction of a landfill are available, it becomes challenging when these data are missing and landfills are located in areas that are already impacted by historical contamination. This work presents a methodology for calculating trigger levels for groundwater quality in landfills located in areas where historical contaminations have deteriorated groundwater quality prior to their construction. This method is based on multivariate statistical analysis and involves 4 steps: (a) implementation of the conceptual model, (b) landfill monitoring data collection, (c) hydrochemical data clustering and (d) calculation of the trigger levels. The proposed methodology was applied on a case study in northern Italy, where a currently used lined landfill is located downstream of an old unlined landfill and others old unmapped waste deposits. The developed conceptual model stated that groundwater quality deterioration observed downstream of the lined landfill is due to a degrading leachate plume fed by the upgradient unlined landfill. The methodology led to the determination of two trigger levels for COD and NH₄-N, the former for a zone representing the background hydrochemistry (28 and 9 mg/L for COD and NH₄-N, respectively), the latter for the zone impacted by the degrading leachate plume from the upgradient unlined landfill (89 and 83 mg/L for COD and NH₄-N, respectively).

Geostatistics as a tool to improve the natural background level definition: An application in groundwater

The Natural Background Level (NBL), suggested by UE BRIDGE project, is suited for spatially distributed datasets providing a regional value that could be higher than the Threshold Value (TV) set by every country. In hydro-geochemically dis-homogeneous areas, the use of a unique regional NBL, higher than TV, could arise problems to distinguish between natural occurrences and anthropogenic contaminant sources. Hence, the goal of this study is to improve the NBL definition employing a geostatistical approach, which reconstructs the contaminant spatial structure accounting geochemical and hydrogeological relationships. This integrated mapping is fundamental to evaluate the contaminant's distribution impact on the NBL, giving indications to improve it. We decided to test this method on the Drainage Basin of Venice Lagoon (DBVL, NE Italy), where the existing NBL is seven times higher than the TV. This area is notoriously affected by naturally occurring arsenic contamination. An available geochemical dataset collected by 50 piezometers was used to reconstruct the spatial distribution of arsenic in the densely populated area of the DBVL. A cokriging approach was applied exploiting the geochemical relationships among As, Fe and NH4+. The obtained spatial predictions of arsenic concentrations were divided into three different zones: i) areas with an As concentration lower than the TV, ii) areas with an As concentration between the TV and the median of the values higher than the TV, and iii) areas with an As concentration higher than the median. Following the BRIDGE suggestions, where enough samples were available, the 90th percentile for each zone was calculated to obtain a local NBL (LNBL). Differently from the original NBL, this local value gives more detailed water quality information accounting the hydrogeological and geochemical setting, and contaminant spatial variation. Hence, the LNBL could give more indications about the distinction between natural occurrence and anthropogenic contamination.

An Integrated Approach for the Environmental Characterization of a Wide Potentially Contaminated Area in Southern Italy

This paper deals with the environmental characterization of a large and densely populated area, with a poor reputation for contamination, considering the contribution of environmental features (air, soil, soil hydraulic and groundwater) and the potential effects on human health. The use of Geographic Information System (GIS) has made possible a georeferenced inventory and, by overlaying environmental information, an operational synthesis of comprehensive environmental conditions. The cumulative effects on environmental features were evaluated, taking into account superposition effects, by means of the Spatial MultiCriteria Decision Analysis (S-MCDA). The application of the S-MCDA for converging the combination of heterogeneous factors, related to soil, land and water, deeply studied by heterogeneous groups of experts, constitutes the novelty of the paper. The results confirmed an overall higher potential of exposure to contaminants in the environment and higher mortality rates in the study area for some tumours, but hospital admissions for tumours were generally similar to the regional trend. Besides, mortality data may be strictly dependent on the poor socioeconomic conditions, quality of therapy and a lack of welfare in the area relative to the rest of Italy. Finally, as regards the possible relationship between presence of contaminants in the environment and health conditions of the population no definite conclusions can be drawn, although the present study encourages the use of the new proposed methods, that increase the possibilities for studying the combined effect of more environmental factors.