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.

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.

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.