Delayed nitrate dispersion within a coastal aquifer provides constraints on land-use evolution and nitrate contamination in the past

Identifying sources of anthropogenic pollution, and assessing the fate and residence time of pollutants in aquifers is important for the management of groundwater resources, and the ecological health of groundwater dependent ecosystems. This study investigates anthropogenic contamination in the shallow alluvial aquifer of the Marana-Casinca, hydraulically connected to the Biguglia lagoon (Corsica, France). A multi-tracer approach, combining geochemical and environmental isotopic data (δ¹⁸O-H₂O, δ²H-H₂O, ³H, δ¹⁵N-NO₃^-, δ¹⁸O-NO₃^-, δ¹¹B), and groundwater residence-time tracers (³H and CFCs) was carried out in 2016, and integrated with a study of land use evolution in the catchment during the last century. Groundwater NO₃^- concentrations, ranged between 2 mg/L and up to 30 mg/L, displaying the degradation of groundwater quality induced by anthropogenic activities (agricultural activities). Comparatively high δ¹⁵N-NO₃^- values (up to 19.7‰) in combination with δ¹¹B values that were significantly lower (between 23‰ and 26‰) than the seawater background are indicative of sewage contamination. The ongoing deterioration of groundwater quality can be attributed to the uncontrolled urbanization development all over the alluvial plain, with numerous sewage leakages from the sanitation network and private sewage systems. Integration of contaminant and water-residence time data revealed a progressive accumulation of pollutants with time in the groundwater, particularly in areas with major anthropogenic pressure and slow dynamic groundwater flow. Our approach provides time-dependent insight into nitrogen pollution in the studied aquifer over the past decades, revealing a systematic change in the dominant NO₃^- source, from agricultural to sewage contamination. Yet, today's low groundwater quality is to large parts due to legacy pollution from land-use practices several decades ago, underlining the poor self-remediating capacity of this hydrosystem. Our results can be taken as warning that groundwater pollution that happened in the recent past, or today, may have dire impacts on the quality of groundwater-dependent ecosystems in the future.

Impact of irrigated agriculture on groundwater resources in a temperate humid region

The groundwater irrigation expansion, and its multiple potential impacts on the quantity and quality of water resources, is not just restricted to areas that are water limited. In this study we present the seasonal impacts irrigation practices can have on groundwater resources in a temperate humid region, where the average annual rain/PET ratio is 1.0. In this system the irrigation expansion is solely supported by groundwater pumping, but despite this only 5 boreholes are monitored for hydraulic head data. In this study, we compensate the scarce hydrophysical dataset by incorporating environmental tracers (major ions, δ¹⁸O, δ²H and δ¹³C) and dating tracers (³H, CFC, SF₆ and ¹⁴C). Results indicate that at 9 of the 15 irrigation sites investigated, groundwater pumping for irrigation has induced the mixing of recent groundwater (up to <1year) with older waters. The origin of the older waters was from either the deeper marl aquifer, or the shallow sand-clay aquifer (SCB) that has a ¹⁴C mean residence time (MRT) of up to 9700years. Secondly, although high nitrate loads in infiltrating waters were being diverted via the artificial subsurface drainage system, increases in fertiliser loads have resulted in higher NO₃ concentrations in younger groundwater (NO₃: 9-45mg/L, MRT <20years), compared with older groundwater (NO₃≤9mg/L, MRT>20years). The changes in flow pathways, induced by irrigation, also results in seasonal declines in groundwater NO₃ concentrations due to mixing with older waters. In temperate humid areas, such evaluations of the seasonal evolution of water residence time, mixing process, and agrochemical contaminants are an important contribution to real water resources management in irrigated catchments.

Strontium isotopes as tracers of water-rocks interactions, mixing processes and residence time indicator of groundwater within the granite-carbonate coastal aquifer of Bonifacio (Corsica, France)

This study aims at identifying the water-rock interactions and mixing rates within a complex granite-carbonate coastal aquifer under high touristic pressure. Investigations have been carried out within the coastal aquifer of Bonifacio (southern Corsica, France) mainly composed of continental granitic weathering products and marine calcarenite sediments filling a granitic depression. A multi-tracer approach combining physico-chemical parameters, major ions, selected trace elements, stable isotopes of the water molecule and ⁸⁷Sr/⁸⁶Sr ratios measurements is undertaken for 20 groundwater samples during the low water period in November 2014. 5 rock samples of the sedimentary deposits and surrounding granites are also analysed. First, the water-rock interactions processes governing the groundwater mineralization are described in order to fix the hydrogeochemical background. Secondly, the flow conditions are refined through the quantification of inter aquifer levels mixing, and thirdly, the kinetics of water-rock interaction based on groundwater residence time from a previous study using CFCs and SF₆ are quantified for the two main flow lines. A regional contrast in the groundwater recharge altitude allowed the oxygene-18 to be useful combined with the ⁸⁷Sr/⁸⁶Sr ratios to differentiate the groundwater origins and to compute the mixing rates, revealing the real extension of the watershed and the availability of the resource. The results also highlight a very good correlation between the groundwater residence time and the spatial evolution of ⁸⁷Sr/⁸⁶Sr ratios, allowing water-rock interaction kinetics to be defined empirically for the two main flow lines through the calcarenites. These results demonstrate the efficiency of strontium isotopes as tracers of water-rock interaction kinetics and by extension their relevance as a proxy of groundwater residence time, fundamental parameter documenting the long term sustainability of the hydrosystem.

Glacial recharge, salinisation and anthropogenic contamination in the coastal aquifers of Recife (Brazil)

Implying large residence times and complex water origins deep coastal aquifers are of particular interest as they are remarkable markers of climate, water use and land use changes. Over the last decades, the Metropolitan Region of Recife (Brazil) went through extensive environmental changes increasing the pressure on water resources and giving rise to numerous environmental consequences on the coastal groundwater systems. We analysed the groundwater of the deep aquifers Cabo and Beberibe that are increasingly exploited. The processes potentially affecting groundwater residence times and flow paths have been studied using a multi-tracer approach (CFCs, SF6, noble gases, 14C, 2H and 18O). The main findings of these investigations show that: (1) Groundwaters of the Cabo and Beberibe aquifers have long residence times and were recharged about 20,000years ago. (2) Within these old groundwaters we can find palaeo-climate evidences from the last glacial period at the tropics with lower temperatures and dryer conditions than the present climate. (3) Recently, the natural slow dynamic of these groundwater systems was significantly affected by mixing processes with contaminated modern groundwater coming from the shallow unconfined Boa Viagem aquifer. (4) The large exploitation of these aquifers leads to a modification of the flow directions and causes the intrusion through palaeo-channels of saline water probably coming from the Capibaribe River and from the last transgression episodes. These observations indicate that the current exploitation of the Cabo and Beberibe aquifers is unsustainable regarding the long renewal times of these groundwater systems as well as their ongoing contamination and salinisation. The groundwater cycle being much slower than the human development rhythm, it is essential to integrate the magnitude and rapidity of anthropogenic impacts on this extremely slow cycle to the water management concepts.

Groundwater contamination in coastal urban areas: Anthropogenic pressure and natural attenuation processes. Example of Recife (PE State, NE Brazil)

In a context of increasing land use pressure (over-exploitation, surface-water contamination) and repeated droughts, identifying the processes affecting groundwater quality in coastal megacities of the tropical and arid countries will condition their long-term social and environmental sustainability. The present study focuses on the Brazilian Recife Metropolitan Region (RMR), which is a highly urbanized area (3,743,854 inhabitants in 2010) on the Atlantic coast located next to an estuarial zone and overlying a multi-layered sedimentary system featured by a variable sediment texture and organic content. It investigates the contamination and redox status patterns conditioning potential attenuation within the shallow aquifers that constitute the interface between the city and the strategic deeper semi-confined aquifers. These latter are increasingly exploited, leading to high drawdown in potenciometric levels of 20-30m and up to 70m in some high well density places, and potentially connected to the surface through leakage. From a multi-tracer approach (major ions, major gases, δ(11)B, δ(18)O-SO4, δ(34)S-SO4) carried out during two field campaigns in September 2012 and March 2013 (sampling of 19 wells and 3 surface waters), it has been possible to assess the contamination sources and the redox processes. The increasing trend for mineralization from inland to coastal and estuarial wells (from 119 to around 10,000μS/cm) is at first attributed to water-rock interactions combined with natural and human-induced potentiometric gradients. Secondly, along with this trend, one finds an environmental pressure gradient related to sewage and/or surface-channel network impacts (typically depleted δ(11)B within the range of 10-15‰) that are purveyors of chloride, nitrate, ammonium and sulfate. Nitrate, ammonium and sulfate (ranging from 0 to 1.70mmol/L, from 0 to 0,65mmol/L, from 0.03 to 3.91mmol/L respectively are also potentially produced or consumed through various redox processes (pyrite oxidation, denitrification, dissimilatory nitrate reduction to ammonium) within the system, as is apparent within a patchwork of biogeochemical reactors. Furthermore, intensive pumping in the coastal area with its high well density punctually leads to temporary well salinization ([Cl] reaching temporarily 79mmol/L). Our results, summarized as a conceptual scheme based on environmental conditions, is a suitable basis for implementing sustainable management in coastal sedimentary hydrosystems influenced by highly urbanized conditions.

Reaction chain modeling of denitrification reactions during a push-pull test

Field quantitative estimation of reaction kinetics is required to enhance our understanding of biogeochemical reactions in aquifers. We extended the analytical solution developed by Haggerty et al. (1998) to model an entire 1st order reaction chain and estimate the kinetic parameters for each reaction step of the denitrification process. We then assessed the ability of this reaction chain to model biogeochemical reactions by comparing it with experimental results from a push-pull test in a fractured crystalline aquifer (Ploemeur, French Brittany). Nitrates were used as the reactive tracer, since denitrification involves the sequential reduction of nitrates to nitrogen gas through a chain reaction (NO3(-)→NO2(-)→NO→N2O→N2) under anaerobic conditions. The kinetics of nitrate consumption and by-product formation (NO2(-), N2O) during autotrophic denitrification were quantified by using a reactive tracer (NO3(-)) and a non-reactive tracer (Br(-)). The formation of reaction by-products (NO2(-), N2O, N2) has not been previously considered using a reaction chain approach. Comparison of Br(-) and NO3(-) breakthrough curves showed that 10% of the injected NO3(-) molar mass was transformed during the 12 h experiment (2% into NO2(-), 1% into N2O and the rest into N2 and NO). Similar results, but with slower kinetics, were obtained from laboratory experiments in reactors. The good agreement between the model and the field data shows that the complete denitrification process can be efficiently modeled as a sequence of first order reactions. The 1st order kinetics coefficients obtained through modeling were as follows: k1=0.023 h(-1), k2=0.59 h(-1), k3=16 h(-1), and k4=5.5 h(-1). A next step will be to assess the variability of field reactivity using the methodology developed for modeling push-pull tracer tests.

Nitrate dynamics in agricultural catchments deduced from groundwater dating and long-term nitrate monitoring in surface- and groundwaters

Although nitrate export in agricultural catchments has been simulated using various types of models, the role of groundwater in nitrate dynamics has rarely been fully taken into account. We used groundwater dating methods (CFC analyses) to reconstruct the original nitrate concentrations in the groundwater recharge in Brittany (Western France) from 1950 to 2009. This revealed a sharp increase in nitrate concentrations from 1977 to 1990 followed by a slight decrease. The recharge concentration curve was then compared with past chronicles of groundwater concentration. Groundwater can be interpreted as resulting from the annual dilution of recharge water in an uncontaminated aquifer. Two aquifers were considered: the weathered aquifer and the deeper fractured aquifer. The nitrate concentrations observed in the upper part of the weathered aquifer implied an annual renewal rate of 27 to 33% of the reservoir volume while those in the lower part indicated an annual renewal rate of 2-3%. The concentrations in the deep fractured aquifer showed an annual renewal rate of 0.1%. The river concentration can be simulated by combining these various groundwater reservoirs with the recharge. Winter and summer waters contain i) recharge water, or water from the variably saturated zone with rapid transfer and high nitrate concentrations, and ii) a large contribution (from 35 to 80% in winter and summer, respectively) from the lower part of the aquifer (lower weathered aquifer and deep fractured aquifer). This induces not only a relatively rapid response of the catchment to variations in agricultural pressure, but also a potential inertia which has to be taken into account.

A modified Lagrangian-volumes method to simulate nonlinearly and kinetically sorbing solute transport in heterogeneous porous media

Transport in subsurface environments is conditioned by physical and chemical processes in interaction, with advection and dispersion being the most common physical processes and sorption the most common chemical reaction. Existing numerical approaches become time-consuming in highly-heterogeneous porous media. In this paper, we discuss a new efficient Lagrangian method for advection-dominated transport conditions. Modified from the active-walker approach, this method comprises dividing the aqueous phase into elementary volumes moving with the flow and interacting with the solid phase. Avoiding numerical diffusion, the method remains efficient whatever the velocity field by adapting the elementary volume transit times to the local velocity so that mesh cells are crossed in a single numerical time step. The method is flexible since a decoupling of the physical and chemical processes at the elementary volume scale, i.e. at the lowest scale considered, is achieved. We implement and validate the approach to the specific case of the nonlinear Freundlich kinetic sorption. The method is relevant as long as the kinetic sorption-induced spreading remains much larger than the dispersion-induced spreading. The variability of the surface-to-volume ratio, a key parameter in sorption reactions, is explicitly accounted for by deforming the shape of the elementary volumes.

Deep water circulation, residence time, and chemistry in a karst complex

We investigated the hydrochemistry of a complex karst hydrosystem made of two carbonate units along a coastal lagoon. Ground water emerges on the lagoon floor from a submarine spring. In addition, thermal waters circulate through the limestone and mix with karst water near the lagoon shore. A distinction between the water from the two carbonate units is related to marine influences and human activities. In one of the massifs, the data show an incongruent dissolution of dolomite with time. In the other system, a slight contamination by saline fluids from the thermal reservoir has led to high calcium and magnesium concentrations. 36Cl, 14C, and 3H data constrain the residence time of the water, and allow for the distinguishing of four circulation types: (1) shallow surface circulation (primarily above sea level) in the karstic units with short residence times (<20 years); (2) shallow subsurface circulation (approximately 0 to -50 m) below the karstic units with residence time in the order of 50 years; (3) deep circulation at depth of 700 to 1500 m in the Jurassic limestones below thick sedimentary cover, with residence time of several thousand years for a part of the water; and (4) deep circulation at a depth of approximately 2500 m, which represents the thermal reservoir in the Jurassic units with residence time of approximately 100,000 years. An interpretative hydrogeological framework is based on the constraints of the geochemical analyses of the deep thermal system, and by water flow from the surface to the deep parts of the carbonate formations.