Deep denitrification: Stream and groundwater biogeochemistry reveal contrasted but connected worlds above and below.

Excess nutrients from agricultural and urban development have created a cascade of ecological crises around the globe. Nutrient pollution has triggered eutrophication in most freshwater and coastal ecosystems, contributing to a loss in biodiversity, harm to human health, and trillions in economic damage every year. Much of the research conducted on nutrient transport and retention has focused on surface environments, which are both easy to access and biologically active. However, surface characteristics of watersheds, such as land use and network configuration, often do not explain the variation in nutrient retention observed in rivers, lakes, and estuaries. Recent research suggests subsurface processes and characteristics may be more important than previously thought in determining watershed-level nutrient fluxes and removal. In a small watershed in western France, we used a multi-tracer approach to compare surface and subsurface nitrate dynamics at commensurate spatiotemporal scales. We combined 3-D hydrological modeling with a rich biogeochemical dataset from 20 wells and 15 stream locations. Water chemistry in the surface and subsurface showed high temporal variability, but groundwater was substantially more spatially variable, attributable to long transport times (10-60 years) and patchy distribution of the iron and sulfur electron donors fueling autotrophic denitrification. Isotopes of nitrate and sulfate revealed fundamentally different processes dominating the surface (heterotrophic denitrification and sulfate reduction) and subsurface (autotrophic denitrification and sulfate production). Agricultural land use was associated with elevated nitrate in surface water, but subsurface nitrate concentration was decoupled from land use. Dissolved silica and sulfate are affordable tracers of residence time and nitrogen removal that are relatively stable in surface and subsurface environments. Together, these findings reveal distinct but adjacent and connected biogeochemical worlds in the surface and subsurface. Characterizing how these worlds are linked and decoupled is critical to meeting water quality targets and addressing water issues in the Anthropocene.

A hillslope-scale aquifer-model to determine past agricultural legacy and future nitrate concentrations in rivers.

The long-term fate of agricultural nitrate depends on rapid subsurface transfer, denitrification and storage in aquifers. Quantifying these processes remains an issue due to time varying subsurface contribution, unknown aquifer storage and heterogeneous denitrification potential. Here, we develop a parsimonious modelling approach that uses long-term discharge and river nitrate concentration time-series combined with groundwater age data determined from chlorofluorocarbons in springs and boreholes. To leverage their informational content, we use a Boussinesq-type equivalent hillslope model to capture the dynamics of aquifer flows and evolving surface and subsurface contribution to rivers. Nitrate transport was modelled with a depth-resolved high-order finite-difference method and denitrification by a first-order law. We applied the method to three heavily nitrate loaded catchments of a crystalline temperate region of France (Brittany). We found that mean water transit time ranged 10-32 years and Damköhler ratio (transit time/denitrification time) ranged 0.12-0.55, leading to limited denitrification in the aquifer (10-20%). The long-term trajectory of nitrate concentration in rivers appears determined by flows stratification in the aquifer. The results suggest that autotrophic denitrification is controlled by the accessibility of reduced minerals which occurs at the base of the aquifer where flows decrease. One interpretation is that denitrification might be an interfacial process in zones that are weathered enough to transmit flows and not too weathered to have remaining accessible reduced minerals. Consequently, denitrification would not be controlled by the total aquifer volume and related mean transit time but by the proximity of the active weathered interface with the water table. This should be confirmed by complementary studies to which the developed methodology might be further deployed.

What do we need to predict groundwater nitrate recovery trajectories?

Nitrate contamination affects many of the Earth's aquifers and surface waters. Large-scale predictions of groundwater nitrate trends normally require the characterization of multiple anthropic and natural factors. To assess different approaches for upscaling estimates of nitrate recovery, we tested the influence of hydrological, historical, and biological factors on predictions of future nitrate concentration in aquifers. We tested the factors with a rich hydrogeological dataset from a fractured bedrock catchment in western France (Brittany). A sensitivity analysis performed on a calibrated model of groundwater flow, denitrification, and nitrogen inputs revealed that trends in nitrate concentration can effectively be approximated with a limited number of key parameters. The total mass of nitrate that entered the aquifer since the beginning of the industrial period needs to be characterized, but the shape of the historical nitrogen input time series can be largely simplified without substantially altering the predictions. Aquifer flow and transport processes can be represented by the mean and standard deviation of the residence time distribution, offering a tractable tool to make reasonable predictions at watershed to regional scales. Apparent sensitivity to denitrification rate was primarily attributable to time lags in oxygen depletion, meaning that denitrification can be simplified to an on/off process, defined only by the time needed for nitrate to reach the hypoxic reactive layer. Obtaining these key parameters at large scales is still challenging with currently available information, but the results are promising regarding our future ability to predict nitrate concentration with integrated monitoring and modeling approaches.

Genome reconstruction reveals distinct assemblages of Gallionellaceae in surface and subsurface redox transition zones.

Fe-oxidizing bacteria of the family Gallionellaceae are major players in the Fe biogeochemical cycle in freshwater. These bacteria thrive in redox transition zones where they benefit from both high Fe concentrations and microaerobic conditions. We analysed the Gallionellaceae genomic diversity in an artesian hard-rock aquifer where redox transition zones develop (i) in the subsurface, where ancient, reduced groundwater mixes with recent oxygenated groundwater, and (ii) at the surface, where groundwater reaches the open air. A total of 15 new draft genomes of Gallionellaceae representing to 11 candidate genera were recovered from the two redox transition zones. Sulfur oxidation genes were encoded in most genomes while denitrification genes were much less represented. One genus dominated microbial communities belowground and we propose to name it 'Candidatus Houarnoksidobacter'. The two transition zones were populated by completely different assemblages of Gallionellaceae despite the almost constant upward circulation of groundwater between the two zones. The processes leading to redox transition zones, oxygen diffusion at the surface or groundwater mixing in subsurface, appear to be a major driver of the Gallionellaceae diversity.

Kinetic Study on Clogging of a Geothermal Pumping Well Triggered by Mixing-Induced Biogeochemical Reactions.

The sustainability of ground-source geothermal systems can be severely impacted by microbially mediated clogging processes. Biofouling of water wells by hydrous ferric oxide is a widespread problem. Although the mechanisms and critical environmental factors associated with clogging development are widely recognized, effects of mixing processes within the wells and time scales for clogging processes are not well characterized. Here we report insights from a joint hydrological, geochemical, and metagenomics characterization of a geothermal doublet in which hydrous ferric oxide and hydrous manganese oxide deposits had formed as a consequence of mixing shallow groundwater containing dissolved oxygen and nitrate with deeper, anoxic groundwater containing dissolved iron (FeII) and manganese (MnII). Metagenomics identify distinct bacteria consortia in the pumping well oxic and anoxic zones, including autotrophic iron-oxidizing bacteria. Batch mixing experiments and geochemical kinetics modeling of the associated reactions indicate that FeII and MnII oxidation are slow compared to the residence time of water in the pumping well; however, adsorption of FeII and MnII by accumulated hydrous ferric oxide and hydrous manganese oxide in the well bore and pump riser provides "infinite" time for surface-catalyzed oxidation and a convenient source of energy for iron-oxidizing bacteria, which colonize the surfaces and also catalyze oxidation. Thus, rapid clogging is caused by mixing-induced redox reactions and is exacerbated by microbial activity on accumulated hydrous oxide surfaces.

Stratification of reactivity determines nitrate removal in groundwater.

Biogeochemical reactions occur unevenly in space and time, but this heterogeneity is often simplified as a linear average due to sparse data, especially in subsurface environments where access is limited. For example, little is known about the spatial variability of groundwater denitrification, an important process in removing nitrate originating from agriculture and land use conversion. Information about the rate, arrangement, and extent of denitrification is needed to determine sustainable limits of human activity and to predict recovery time frames. Here, we developed and validated a method for inferring the spatial organization of sequential biogeochemical reactions in an aquifer in France. We applied it to five other aquifers in different geological settings located in the United States and compared results among 44 locations across the six aquifers to assess the generality of reactivity trends. Of the sampling locations, 79% showed pronounced increases of reactivity with depth. This suggests that previous estimates of denitrification have underestimated the capacity of deep aquifers to remove nitrate, while overestimating nitrate removal in shallow flow paths. Oxygen and nitrate reduction likely increases with depth because there is relatively little organic carbon in agricultural soils and because excess nitrate input has depleted solid phase electron donors near the surface. Our findings explain the long-standing conundrum of why apparent reaction rates of oxygen in aquifers are typically smaller than those of nitrate, which is energetically less favorable. This stratified reactivity framework is promising for mapping vertical reactivity trends in aquifers, generating new understanding of subsurface ecosystems and their capacity to remove contaminants.

Autotrophic denitrification supported by biotite dissolution in crystalline aquifers (1): New insights from short-term batch experiments.

We investigate denitrification mechanisms through batch experiments using crushed rock and groundwater from a granitic aquifer subject to long term pumping (Ploemeur, France). Except for sterilized experiments, extensive denitrification reaction induces NO3 decreases ranging from 0.3 to 0.6mmol/L. Carbon concentrations, either organic or inorganic, remain relatively stable and do not document potential heterotrophic denitrification. Batch experiments show a clear effect of mineral dissolution which is documented through cation (K, Na, Ca) and Fluoride production. These productions are tightly related to denitrification progress during the experiment. Conversely, limited amounts of SO4, systematically lower than autotrophic denitrification coupled to sulfur oxidation stoichiometry, are produced during the experiments which indicates that sulfur oxidation is not likely even when pyrite is added to the experiments. Analysis of cation ratios, both in isolated minerals of the granite and within water of the batch, allow the mineral dissolution during the experiments to be quantified. Using cation ratios, we show that batch experiments are characterized mainly by biotite dissolution. As biotite contains 21 to 30% of Fe and 0.3 to 1.7% of F, it constitutes a potential source for these two elements. Denitrification could be attributed to the oxidation of Fe(II) contained in biotite. We computed the amount of K and F produced through biotite dissolution when entirely attributing denitrification to biotite dissolution. Computed amounts show that this process may account for the observed K and F produced. We interpret these results as the development of microbial activity which induces mineral dissolution in order to uptake Fe(II) which is used for denitrification. Although pyrite is probably available, SO4 and cation measurements favor a large biotite dissolution reaction which could account for all the observed Fe production. Chemical composition of groundwater produced from the Ploemeur site indicates similar denitrification processes although original composition shows mainly plagioclase dissolution.

Dating groundwater with dissolved silica and CFC concentrations in crystalline aquifers.

Estimating intermediate water residence times (a few years to a century) in shallow aquifers is critical to quantifying groundwater vulnerability to nutrient loading and estimating realistic recovery timelines. While intermediate groundwater residence times are currently determined with atmospheric tracers such as chlorofluorocarbons (CFCs), these analyses are costly and would benefit from other tracer approaches to compensate for the decreasing resolution of CFC methods in the 5-20 years range. In this context, we developed a framework to assess the capacity of dissolved silica (DSi) to inform residence times in shallow aquifers. We calibrated silicate weathering rates with CFCs from multiple wells in five crystalline aquifers in Brittany and in the Vosges Mountains (France). DSi and CFCs were complementary in determining apparent weathering reactions and residence time distributions (RTDs) in shallow aquifers. Silicate weathering rates were surprisingly similar among Brittany aquifers, varying from 0.20 to 0.23 mg L-1 yr-1 with a coefficient of variation of 7%, except for the aquifer where significant groundwater abstraction occurred, where we observed a weathering rate of 0.31 mg L-1 yr-1. The silicate weathering rate was lower for the aquifer in the Vosges Mountains (0.12 mg L-1 yr-1), potentially due to differences in climate and anthropogenic solute loading. Overall, these optimized silicate weathering rates are consistent with previously published studies with similar apparent ages range. The consistency in silicate weathering rates suggests that DSi could be a robust and cheap proxy of mean residence times for recent groundwater (5-100 years) at the regional scale. This methodology could allow quantification of seasonal groundwater contributions to streams, estimation of residence times in the unsaturated zone and improve assessment of aquifer vulnerability to anthropogenic pollution.

Multi-layered water resources, management, and uses under the impacts of global changes in a southern coastal metropolis: When will it be already too late? Crossed analysis in Recife, NE Brazil.

Coastal water resources are a worldwide key socio-environmental issue considering the increasing concentration of population in these areas. Here, we propose an integrative transdisciplinary approach of water resource, water management and water access in Recife (NE Brazil). The present-day water situation is conceptualized as an imbricated multi-layered system: a multi-layered water resource, managed by a multi-layered governance system and used by a multi-layered social population. This allows identifying processes of quantitative, qualitative, and sanitary conflicts between governance and population strategies regarding water supply, as well as the institutional and individual denials of these conflicts. Based on this model, we anticipate future water-related problematic fates. Concerning the water resource system, the rapid groundwater level decrease due to unsustainable water predatory strategies, and the very low recharge rate have drastically modified the aquifer system functioning, inducing hydraulic connection between shallow groundwater (contaminated and locally salty) and deep ones (mostly fresh, with local inherited salinity), threatening the deep strategic water resource. Concerning the water governance system, the investments to increase the capacity storage of surface water, the water regulation agencies and the public/private partnership should shortly improve the water supply and wastewater issue. Nevertheless, the water situation will remain highly fragile due to the expected water demand increase, the precipitation decrease and the sea-level increase. Concerning the water access system, the population variably perceives these current and further effects and the possible mitigation policies, and develops alternative individual strategies. Authorities, policymakers and water managers will have to implement a well-balanced water governance, taking into account the specificities of the PPP, public and private groundwater users, and with a strong political willingness for a sustainable water management to ensure water supply for all the population. In other words, an anticipatory and integrated vision is necessary to reduce the discrepancies in this complex system.

Autotrophic denitrification supported by biotite dissolution in crystalline aquifers: (2) transient mixing and denitrification dynamic during long-term pumping.

We investigated the mixing and dynamic of denitrification processes induced by long-term pumping in the crystalline aquifer of Ploemeur (Brittany, France). Hydrological and geochemical parameters have been continuously recorded over 15 boreholes in 5km2 on a 25-year period. This extensive spatial and temporal monitoring of conservative as well as reactive compounds is a key opportunity to identify aquifer-scale transport and reactive processes in crystalline aquifers. Time series analysis of the conservative elements recorded at the pumped well were used to determine mixing fractions from different compartments of the aquifer on the basis of a Principal Component Analysis approach coupled with an end-member mixing analysis. We could reveal that pumping thus induces a thorough reorganization of fluxes known as capture, favoring infiltration and vertical fluxes in the recharge zone, and upwelling of deep and distant water at long-term time scales. These mixing fractions were then used to quantify the extent of denitrification linked to pumping. Based on the results from batch experiments described in a companion paper, our computations revealed that i) autotrophic denitrification processes are dominant in this context where carbon sources are limited, that ii) nitrate reduction does not only come from the oxidation of pyrite as classically described in previous studies analyzing denitrification processes in similar contexts, and that iii) biotite plays a critical role in sustaining the nitrate reduction process. Both nitrate reduction, sulfate production as well as fluor release ratios support the hypothesis that biotite plays a key role of electron donor in this context. The batch-to-site similarities support biotite availability and the role by bacterial communities as key controls of nitrate removal in such crystalline aquifers. However, the long term data monitoring also indicates that mixing and reactive processes evolve extremely slowly at the scale of the decade.

Field Continuous Measurement of Dissolved Gases with a CF-MIMS: Applications to the Physics and Biogeochemistry of Groundwater Flow.

In the perspective of a temporal and spatial exploration of aquatic environments (surface and groundwater), we developed a technique for field continuous measurements of dissolved gases with a precision better than 1% for N2, O2, CO2, He, Ar, 2% for Kr, 8% for Xe, and 3% for CH4, N2O and Ne. With a large resolution (from 1 × 10-9 to 1 × 10-2 ccSTP/g) and a capability of high frequency analysis (1 measure every 2 s), the CF-MIMS (Continuous Flow Membrane Inlet Mass Spectrometer) is an innovative tool allowing the investigation of a large panel of hydrological and biogeochemical processes in aquatic systems. Based on the available MIMS technology, this study introduces the development of the CF-MIMS (conception for field experiments, membrane choices, ionization) and an original calibration procedure allowing the quantification of mass spectral overlaps and temperature effects on membrane permeability. This study also presents two field applications of the CF-MIMS involving the well-logging of dissolved gases and the implementation of groundwater tracer tests with dissolved 4He. The results demonstrate the analytical capabilities of the CF-MIMS in the field. Therefore, the CF-MIMS is a valuable tool for the field characterization of biogeochemical reactivity, aquifer transport properties, groundwater recharge, groundwater residence time and aquifer-river exchanges from few hours to several weeks experiments.

Time-scales of hydrological forcing on the geochemistry and bacterial community structure of temperate peat soils.

Peatlands are an important global carbon reservoir. The continued accumulation of carbon in peatlands depends on the persistence of anoxic conditions, in part induced by water saturation, which prevents oxidation of organic matter, and slows down decomposition. Here we investigate how and over what time scales the hydrological regime impacts the geochemistry and the bacterial community structure of temperate peat soils. Peat cores from two sites having contrasting groundwater budgets were subjected to four controlled drought-rewetting cycles. Pore water geochemistry and metagenomic profiling of bacterial communities showed that frequent water table drawdown induced lower concentrations of dissolved carbon, higher concentrations of sulfate and iron and reduced bacterial richness and diversity in the peat soil and water. Short-term drought cycles (3-9 day frequency) resulted in different communities from continuously saturated environments. Furthermore, the site that has more frequently experienced water table drawdown during the last two decades presented the most striking shifts in bacterial community structure, altering biogeochemical functioning of peat soils. Our results suggest that the increase in frequency and duration of drought conditions under changing climatic conditions or water resource use can induce profound changes in bacterial communities, with potentially severe consequences for carbon storage in temperate peatlands.

Impact of climate changes during the last 5 million years on groundwater in basement aquifers.

Climate change is thought to have major effects on groundwater resources. There is however a limited knowledge of the impacts of past climate changes such as warm or glacial periods on groundwater although marine or glacial fluids may have circulated in basements during these periods. Geochemical investigations of groundwater at shallow depth (80-400 m) in the Armorican basement (western France) revealed three major phases of evolution: (1) Mio-Pliocene transgressions led to marine water introduction in the whole rock porosity through density and then diffusion processes, (2) intensive and rapid recharge after the glacial maximum down to several hundred meters depths, (3) a present-day regime of groundwater circulation limited to shallow depth. This work identifies important constraints regarding the mechanisms responsible for both marine and glacial fluid migrations and their preservation within a basement. It defines the first clear time scales of these processes and thus provides a unique case for understanding the effects of climate changes on hydrogeology in basements. It reveals that glacial water is supplied in significant amounts to deep aquifers even in permafrosted zones. It also emphasizes the vulnerability of modern groundwater hydrosystems to climate change as groundwater active aquifers is restricted to shallow depths.

Origins and processes of groundwater salinization in the urban coastal aquifers of Recife (Pernambuco, Brazil): A multi-isotope approach.

In the coastal multilayer aquifer system of a highly urbanized southern city (Recife, Brazil), where groundwaters are affected by salinization, a multi-isotope approach (Sr, B, O, H) was used to investigate the sources and processes of salinization. The high diversity of the geological bodies, built since the Atlantic opening during the Cretaceous, highly constrains the heterogeneity of the groundwater chemistry, e.g. Sr isotope ratios, and needs to be integrated to explain the salinization processes and groundwater pathways. A paleoseawater intrusion, most probably the 120 kyB.P. Pleistocene marine transgression, and cationic exchange are clearly evidenced in the most salinized parts of the Cabo and Beberibe aquifers. All (87)Sr/(86)Sr values are above the past and present-day seawater signatures, meaning that the Sr isotopic signature is altered due to additional Sr inputs from dilution with different freshwaters, and water-rock interactions. Only the Cabo aquifer presents a well-delimitated area of Na-HCO3 water typical of a freshening process. The two deep aquifers also display a broad range of B concentrations and B isotope ratios with values among the highest known to date (63-68.5‰). This suggests multiple sources and processes affecting B behavior, among which mixing with saline water, B sorption on clays and mixing with wastewater. The highly fractionated B isotopic values were explained by infiltration of relatively salty water with B interacting with clays, pointing out the major role played by (palaeo)-channels for the deep Beberibe aquifer recharge. Based on an increase of salinity at the end of the dry season, a present-day seawater intrusion is identified in the surficial Boa Viagem aquifer. Our conceptual model presents a comprehensive understanding of the major groundwater salinization pathways and processes, and should be of benefit for other southern Atlantic coastal aquifers to better address groundwater management issues.

Groundwater Isolation Governs Chemistry and Microbial Community Structure along Hydrologic Flowpaths.

This study deals with the effects of hydrodynamic functioning of hard-rock aquifers on microbial communities. In hard-rock aquifers, the heterogeneous hydrologic circulation strongly constrains groundwater residence time, hydrochemistry, and nutrient supply. Here, residence time and a wide range of environmental factors were used to test the influence of groundwater circulation on active microbial community composition, assessed by high throughput sequencing of 16S rRNA. Groundwater of different ages was sampled along hydrogeologic paths or loops, in three contrasting hard-rock aquifers in Brittany (France). Microbial community composition was driven by groundwater residence time and hydrogeologic loop position. In recent groundwater, in the upper section of the aquifers or in their recharge zone, surface water inputs caused high nitrate concentration and the predominance of putative denitrifiers. Although denitrification does not seem to fully decrease nitrate concentrations due to low dissolved organic carbon concentrations, nitrate input has a major effect on microbial communities. The occurrence of taxa possibly associated with the application of organic fertilizers was also noticed. In ancient isolated groundwater, an ecosystem based on Fe(II)/Fe(III) and S/SO4 redox cycling was observed down to several 100 of meters below the surface. In this depth section, microbial communities were dominated by iron oxidizing bacteria belonging to Gallionellaceae. The latter were associated to old groundwater with high Fe concentrations mixed to a small but not null percentage of recent groundwater inducing oxygen concentrations below 2.5 mg/L. These two types of microbial community were observed in the three sites, independently of site geology and aquifer geometry, indicating hydrogeologic circulation exercises a major control on microbial communities.

Long-term effects of high nitrogen loads on cation and carbon riverine export in agricultural catchments.

The intensification of agriculture in recent decades has resulted in extremely high nitrogen inputs to ecosystems. One effect has been H(+) release through NH(4)(+) oxidation in soils, which increases rock weathering and leads to acidification processes such as base-cation leaching from the soil exchange complex. This study investigated the evolution of cation concentrations over the past 50 years in rivers from the Armorican crystalline shield (Brittany, western France). On a regional scale, acidification has resulted in increased base-cation riverine exports (Ca(2+), Mg(2+), Na(+), K(+)) correlated with the increased NO(3)(-) concentration. The estimated cation increase is 0.7 mmol(+)/L for Ca(2+) + Mg(2+) and 0.85 mmol(+)/L for total cations. According to mass balance, cation loss represents >30% of the base-cation exchange capacity of soils. Long-term acidification thus contributes to a decline in soil productivity. Estimates of the total organic nitrogen annually produced worldwide indicate that acidification may also constitute an additional carbon source in crystalline catchments if compensated by liming practices.

Influence of depth and time on diversity of free-living microbial community in the variably saturated zone of a granitic aquifer.

We investigated the temporal and vertical changes in the microbial communities related to hydrological variations an aquifer (Brittany, France). Five water samplings were carried out, spanning three hydrological cycles in the variably and the permanently saturated zones. Seasonal variations in the major anion concentrations (NO3 -, SO4 2- and Cl(-) ) indicated that different physical processes occurred during the recharge process in the two zones. The variably saturated zone is mainly dominated by diffusion and advection processes from the soil, whereas the permanently saturated zone is controlled by moderate advective transfer from the variably saturated zone. Bacterial diversity was investigated by flow cytometry, 16S rRNA and narG genes analyses. Part of this diversity was new in that 6 of the 27 16S rRNA gene sequence phylotypes were unknown even at the class or phylum level. The narG gene analysis did not reveal any clear variation in time or depth within the nitrate reducers' community. In contrast, 16S rRNA gene analyses showed modifications of community composition that could be related to the hydrologic and chemical contrast between the two zones. It was concluded that the physical processes of water transfer could influence bacterial diversity at the soil-aquifer interface.