Spatiotemporal Dynamics of Nitrous Oxide Emission Hotspots in Heterogeneous Riparian Sediments. Issue 12 (14th December 2021)

• Record Type: Journal Article
• Authors: Wallace, Corey D.; Tonina, Daniele; McGarr, Jeffrey T.; de Barros, Felipe P. J.; Soltanian, Mohamad Reza
• Is Part Of: Water resources research. Volume 57:Issue 12(2021)
• Page Start: n/a
• Page End: n/a
• Publication Date: 2021-12-14
• Subjects: hyporheic -- nitrous oxide -- denitrification -- sediment heterogeneity -- flow topology -- numerical modeling; Hydrology -- Periodicals; 333.91
• DOI: 10.1029/2021WR030496
• Languages: English
• ISSNs: 0043-1397

Abstract

Abstract: Nitrous oxide (N2 O) is a potent ozone‐depleting greenhouse gas produced by incomplete denitrification. Recent works on riverine N2 O emissions focus mainly on contributions from in‐channel, benthic, and fluvial hyporheic environments under assumptions of steady‐state conditions and homogeneous sediment hydraulic conductivity ( K ). However, riparian floodplains are also a potentially important N2 O source characterized by complex sediment heterogeneity and dynamic surface and groundwater interactions. We use numerical flow and reactive transport models to investigate the influence of complex sedimentary architecture and high‐flow events (e.g., storms) on N2 O production. We interpret the correlation between flow and solute fields with the flow topological Okubo‐Weiss metric (OW) and the scalar dissipation rate weighted by soil organic matter (OM) fraction and soil saturation. We model a heterogeneous riparian floodplain based on field observations from the Theis Environmental Monitoring and Modeling Site, Ohio, USA. N2 O production is greatest within intermediate‐ K sediments (e.g., sands) where denitrification rates are highest, and emissions increase by more than an order of magnitude during storms. Sensitivity analysis reveals that the denitrification rate is most influential for N2 O flux, accounting for nearly 46% of the variance in production rates. Denitrification rates adapt to spatial changes in the flow topology (measured by OW) related to sediment heterogeneity and are strongly influenced by subsurface mixing dynamics. Mixing is greatest in shear strain‐dominated regions, while vorticity promotes OM dissolution and prolongs residence times. Accurate lithologic representation is imperative for characterizing subsurface N2 O production dynamics, especially given growing concern regarding climate change driven hydrologic changes within watersheds worldwide. Key Points: Subsurface N2 O production increases in response to storms, but sediment heterogeneity controls the location and magnitude of emissions Accurate representation of heterogeneity is important for characterizing subsurface nitrogen transformations and N2 O production dynamics Heterogeneity introduces shear into the groundwater flow topology, enhancing mixing and increasing the spatial extent of denitrification