Monitoring groundwater nitrate attenuation in a regional system coupling hydrogeology with multi-isotopic methods: The case of Plana de Vic (Osona, Spain)

This paper describes an integrated application of classical hydrogeological methods and multi-isotopic methods (δ¹⁵N, δ O N O 3 18 , δ³⁴S, δ O S O 4 18 , δ¹³C) to assess the fate of groundwater nitrate in the Osona area, declared vulnerable to nitrate pollution by the Catalan Government in 1998, where nitrate is derived from intensive pig farming activities. Previous studies, involving a small area, indicated the occurrence of denitrification processes and their relationship with pyrite oxidation [Vitòria, L., Soler, A., Canals, A., Otero, N., 2008. Environmental isotopes (N, S, C, O, D) to determine natural attenuation processes in nitrate contaminated waters: example of Osona (NE Spain). Appl. Geochem. 23, 3597-3611]. For the present study, groundwater samples were collected at 60 production wells at three different periods between April 2005 and May 2006 to confirm that denitrification takes place in a larger area than that studied by Vitòria et al. [Vitòria, L., Soler, A., Canals, A., Otero, N., 2008. Environmental isotopes (N, S, C, O, D) to determine natural attenuation processes in nitrate contaminated waters: example of Osona (NE Spain). Appl. Geochem. 23, 3597-3611]. The aim of the study was to characterize the denitrification processes that control natural attenuation and to study their spatial and temporal variations. Nitrate concentration ranged from 10 to 529mg/l, with 82% of the wells above the drinking water threshold of 50mg NO₃/l. Nitrate isotopic composition ranged from +5.3[per thousand] to +35.3[per thousand] for δ¹⁵N and from +0.4[per thousand] to +17.6[per thousand] for δ O N O 3 18 , and the samples showed a positive correlation between δ¹⁵N and δ O N O 3 18 , with a εN/εO ratio of 1.8, consistent with denitrification processes. The link between denitrification and pyrite oxidation is demonstrated by coupling chemical data with nitrate and sulfate isotopes. Furthermore, a spatial distribution of samples with significant denitrification was observed, allowing us to determine two main hydrogeological zones where natural attenuation was most effective. In several of the studied points, denitrification processes related to pyrite oxidation predominated and an estimation of the isotopic enrichment factors was performed using the temporal variations of nitrate concentration and the isotopic composition of dissolved nitrate ( δ N N O 3 15 and δ O N O 3 18 ). Finally, using estimated isotopic enrichment factors, an approximation of the degree of natural attenuation of nitrate was performed on those samples showing clear denitrification, and a median value of 30% of contaminant diminution was obtained.