Unfolding the spatial heterogeneity of the natural background level of arsenic in groundwater at the mesoscale using data from sites under remediation

The Natural Background Level (NBL) of contaminants in groundwater is typically determined using regional-scale monitoring networks or site-specific studies. However, regional scale values are limited in their ability to capture natural heterogeneities that affect contaminant mobility at smaller scales, potentially leading to local over- or underestimation of the natural contaminant concentration. Conversely, site-specific studies can be expensive and time-consuming, with limited use outside the specified case study. To overcome this issue, a study was conducted in a 2600 km2 area, analyzing arsenic concentration values from monitoring networks of sites under remediation as an alternative source of information. The main drawbacks of the alternative dataset were the lack of information on monitoring procedures at the remediation sites or potential anthropogenic influences on the concentration data. However, these limitations were adequately managed with a thorough data pre-treatment procedure informed by a conceptual model of the study area. The NBLs estimated with the alternative dataset were more reliable than that from the regional monitoring network, which, in the worst case (i.e., in the area with the highest geological and geochemical heterogeneity), the NBL of one order of magnitude was underestimated. As a future step, the project seeks to incorporate geological and geochemical heterogeneities as secondary variables in a geostatistical analysis to produce a continuous distribution of arsenic concentrations at the mesoscale. This would provide a useful tool for managing contaminated sites and a reproducible protocol for NBL derivation in different areas, overcoming the scale issue.