An isotope mixing-based clustering approach for an improved apportionment of nitrate sources in groundwater systems

Several groundwater nitrate (NO₃^-) studies have focused on apportioning NO₃^- sources for various sets of groundwater samples statistically clustered (i.e., grouped) based on hydrochemical data. This can be useful when NO₃^- concentrations in groundwater systems are influenced by the aquifer aqueous geochemistry. Nonetheless, several studies observed independent behavior of NO₃^- reflected by its insignificant correlation with most major ions in groundwater. Therefore, this study introduces a novel clustering approach using NO₃^- concentrations and δ¹⁵N_NO3 values for classifying samples controlled by common mixing of NO₃^--containing groundwaters. This clustering, constituting a first step of an approach allowing an improved apportionment of NO₃^- sources by MixSIAR modeling, is followed by (ii) NO₃^- source identification based not only on δ¹⁵N_NO3 and δ1⁸O_NO3 tracers, but also on δ¹¹B, which is particularly able to distinguish manure from sewage; (iii) constraining the model by selected NO₃^- sources and site-specific isotope fractionation effects; and (vi) a separate NO₃^- sources apportionment processing step for different clustered groundwater samples. The approach is applied to a dataset from a coastal Mediterranean agricultural area affected by several sources of NO₃^-. Two mixing scenarios with independent clusters were established, and the most realistic one that integrates four different clusters was retained for the subsequent steps. Mixing in three clusters is controlled by two NO₃^- sources (i.e., manure/sewage for two clusters, and manure/synthetic fertilizers for one cluster), whereas the fourth cluster is governed by mixing of three NO₃^- sources (i.e., manure, sewage, and synthetic fertilizers). MixSIAR modeling reveals manure as the major source contributing NO₃^- to groundwater according to the four clusters, but with different contributions varying from 61 to 71 %; consistent with differences in land use of the study area. Ultimately, the approach presented in this study provides a process to implement selective NO₃^- mitigation strategies, adapted to discriminated sites, rather than a general NO₃^- mitigation strategy applied to the whole study area.