SEARCH · Geology Reports

Results for “Interfaces”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 865 records · Page 48Linked to original sources

USGS2008

Not all water becomes wine: Sulfur inputs as an opportune tracer of hydrochemical losses from vineyards

California's widespread and economically important vineyards offer substantial opportunities to understand the interface between hydrology and biogeochemistry in agricultural soils. The common use of native sulfur (S) as a fumigant or soil additive provides a novel way to isotopically differentiate among sulfate (SO 4 2− ) pools, allowing the estimation of water and SO 4 2− budgets. The objectives of this study were (1) to characterize the near‐surface hydrological flow paths in a vineyard during irrigation and storm events and (2) to determine how those flow paths affect the fate and transport of SO 4 2− across seasons. Integrating hydrological theory with measurements of SO 4 2− concentration and sulfate‐S isotopic ratios (expressed as [SO 4 2− ] and δ 34 S, respectively) in inputs, soil water, and leachate provided a means of determining flow paths. Low [SO 4 2− ] and δ 34 S in leachate during 4‐h irrigation events reflect minimal engagement of the soil matrix, indicating that preferential flow was the dominant path for water in the near surface. In contrast, high [SO 4 2− ] and δ 34 S values during 8‐h irrigation and storm events reflect near‐complete engagement of the soil matrix, indicating that lateral flow was the dominant pathway. Because hydrologic response and SO 4 2− mobility are tightly coupled in these soils, the magnitude of water fluxes through the near surface controls S cycling both on and off site. These results indicate that preferential flow is an important loss pathway to consider in managing both water resources and water quality (reactive elements) in vineyard land use systems.

Water Resources Research↗

USGS2009

Three‐dimensional model of Hellenic Arc deformation and origin of the Cretan uplift

[1] The Hellenic Arc of Greece is the most seismically active part of Europe, but little is know about its mechanics. We modeled deformation along the arc using a finite element model. The model was intended to capture large‐scale 3‐D structure of Nubian plate subduction beneath the Aegean block and its deformational consequences. The shape of the interface was developed using mapped traces at the surface and earthquake hypocenters at depth. Model block motions were constrained by recent compilations of GPS velocity vectors. We simulated a 10 ka period of convergence between Nubia and the Aegean and calculated the strain field in the overriding plate as well as the spatial distribution and orientation of differential stress (∣ σ 1 − σ 3 ∣). From these calculations we derived testable quantities such as the expected seismic moment rate on the interplate contact, uplift pattern, and distribution of strain modes. Our relatively simple model broadly reproduced observed uplift patterns, earthquake activity, and loci of extension and contraction. The model showed a localization of uplift near the island of Crete, where the fastest Aegean uplift rates are well documented. Comparison of calculated expected seismic moment and observed earthquake catalogs implies a nearly fully coupled interplate contact. On the basis of our modeling results, we suggest that south Aegean deformation is driven primarily by the fast moving (∼33 mm a −1 ) Aegean upper plate overriding a nearly stalled (∼5 mm a −1 ) Nubian lower plate. This tectonic setting thus more closely resembles a continental thrust than it does a typical oceanic subduction zone.

Aegean Sea, Hellenic Arc↗

USGS2008

Scaling hyporheic exchange and its influence on biogeochemical reactions in aquatic ecosystems

Hyporheic exchange and biogeochemical reactions are difficult to quantify because of the range in fluid‐flow and sediment conditions inherent to streams, wetlands, and nearshore marine ecosystems. Field measurements of biogeochemical reactions in aquatic systems are impeded by the difficulty of measuring hyporheic flow simultaneously with chemical gradients in sediments. Simplified models of hyporheic exchange have been developed using Darcy's law generated by flow and bed topography at the sediment‐water interface. However, many modes of transport are potentially involved (molecular diffusion, bioturbation, advection, shear, bed mobility, and turbulence) with even simple models being difficult to apply in complex natural systems characterized by variable sediment sizes and irregular bed geometries. In this study, we synthesize information from published hyporheic exchange investigations to develop a scaling relationship for estimating mass transfer in near‐surface sediments across a range in fluid‐flow and sediment conditions. Net hyporheic exchange was quantified using an effective diffusion coefficient ( D e ) that integrates all of the various transport processes that occur simultaneously in sediments, and dimensional analysis was used to scale D e to shear stress velocity, roughness height, and permeability that describe fluid‐flow and sediment characteristics. We demonstrated the value of the derived scaling relationship by using it to quantify dissolved oxygen (DO) uptake rates on the basis of DO profiles in sediments and compared them to independent flux measurements. The results support a broad application of the D e scaling relationship for quantifying coupled hyporheic exchange and biogeochemical reaction rates in streams and other aquatic ecosystems characterized by complex fluid‐flow and sediment conditions.

Water Resources Research↗

USGS2010

Cyclic ground tilt associated with the 2004–2008 eruption of Mount St. Helens

The long‐term behavior of the 2004–2008 effusive eruption of Mount St. Helens was characterized by a gradual decline in the rates of seismicity, dome growth, and broad‐scale ground deformation, but shallow near‐periodic “drumbeat” earthquakes over timescales of minutes indicated episodic short‐term behavior. In part to better characterize this behavior and any associated ground deformation, a network of tiltmeters was installed and recorded thousands of cyclic tilt events within the crater. The duration of these events, from minutes to hours, was too long for them to be directly caused by the drumbeat seismicity. Tilt events were characterized by a recoverable, asymmetric pattern of rapid tilt away from the vent followed by a more gradual reversal, were highly correlated between different tiltmeters in the crater, sometimes occurred in association with volcanic seismicity, and ceased at the end of the eruption. Tilt vectors converged on a point just south of the center of the preexisting 1980s lava dome, and the absence of detectable tilt outside the crater suggests a shallow source (<1 km). We examine several models, including cycles of conduit pressurization and plug slip or gas loss and stick‐slip behavior on the interface between the ascending plug and the 1980s lava dome. The small number of stations within the crater prevents a unique determination of source type or geometry, but results are consistent with a mechanism involving extrusion of the semisolid dacite plug and/or cycles of conduit pressurization.

Washington↗

USGS2011

Flow structures and sandbar dynamics in a canyon river during a controlled flood, Colorado River, Arizona

In canyon rivers, debris fan constrictions create rapids and downstream pools characterized by secondary flow structures that are closely linked to channel morphology. In this paper we describe detailed measurements of the three‐dimensional flow structure and sandbar dynamics of two pools along the Colorado River in the Grand Canyon during a controlled flood release from Glen Canyon Dam. Results indicate that the pools are characterized by large lateral recirculation zones (eddies) resulting from flow separation downstream from the channel constrictions, as well as helical flow structures in the main channel and eddy. The lateral recirculation zones are low‐velocity areas conducive to fine sediment deposition, particularly in the vicinity of the separation and reattachment points and are thus the dominant flow structures controlling sandbar dynamics. The helical flow structures also affect morphology but appear secondary in importance to the lateral eddies. During the controlled flood, sandbars in the separation and reattachment zones at both sites tended to build gradually during the rising limb and peak flow. Deposition in shallow water on the sandbars was accompanied by erosion in deeper water along the sandbar slope at the interface with the main channel. Erosion occurred via rapid mass failures as well as by gradual boundary shear stress driven processes. The flow structures and morphologic links at our study sites are similar to those identified in other river environments, in particular sharply curved meanders and channel confluences where the coexistence of lateral recirculation and helical flows has been documented.

Journal of Geophysical Research F: Earth Surface↗

USGS2010

Use of electrical imaging and distributed temperature sensing methods to characterize surface water–groundwater exchange regulating uranium transport at the Hanford 300 Area, Washington

We explored the use of continuous waterborne electrical imaging (CWEI), in conjunction with fiber‐optic distributed temperature sensor (FO‐DTS) monitoring, to improve the conceptual model for uranium transport within the Columbia River corridor at the Hanford 300 Area, Washington. We first inverted resistivity and induced polarization CWEI data sets for distributions of electrical resistivity and polarizability, from which the spatial complexity of the primary hydrogeologic units was reconstructed. Variations in the depth to the interface between the overlying coarse‐grained, high‐permeability Hanford Formation and the underlying finer‐grained, less permeable Ringold Formation, an important contact that limits vertical migration of contaminants, were resolved along ∼3 km of the river corridor centered on the 300 Area. Polarizability images were translated into lithologic images using established relationships between polarizability and surface area normalized to pore volume ( S por ). The FO‐DTS data recorded along 1.5 km of cable with a 1 m spatial resolution and 5 min sampling interval revealed subreaches showing (1) temperature anomalies (relatively warm in winter and cool in summer) and (2) a strong correlation between temperature and river stage (negative in winter and positive in summer), both indicative of reaches of enhanced surface water–groundwater exchange. The FO‐DTS data sets confirm the hydrologic significance of the variability identified in the CWEI and reveal a pattern of highly focused exchange, concentrated at springs where the Hanford Formation is thickest. Our findings illustrate how the combination of CWEI and FO‐DTS technologies can characterize surface water–groundwater exchange in a complex, coupled river‐aquifer system.

Washington↗

USGS2012

Slab1.0: A three-dimensional model of global subduction zone geometries

We describe and present a new model of global subduction zone geometries, called Slab1.0. An extension of previous efforts to constrain the two-dimensional non-planar geometry of subduction zones around the focus of large earthquakes, Slab1.0 describes the detailed, non-planar, three-dimensional geometry of approximately 85% of subduction zones worldwide. While the model focuses on the detailed form of each slab from their trenches through the seismogenic zone, where it combines data sets from active source and passive seismology, it also continues to the limits of their seismic extent in the upper-mid mantle, providing a uniform approach to the definition of the entire seismically active slab geometry. Examples are shown for two well-constrained global locations; models for many other regions are available and can be freely downloaded in several formats from our new Slab1.0 website, http://on.doi.gov/d9ARbS . We describe improvements in our two-dimensional geometry constraint inversion, including the use of ‘average’ active source seismic data profiles in the shallow trench regions where data are otherwise lacking, derived from the interpolation between other active source seismic data along-strike in the same subduction zone. We include several analyses of the uncertainty and robustness of our three-dimensional interpolation methods. In addition, we use the filtered, subduction-related earthquake data sets compiled to build Slab1.0 in a reassessment of previous analyses of the deep limit of the thrust interface seismogenic zone for all subduction zones included in our global model thus far, concluding that the width of these seismogenic zones is on average 30% larger than previous studies have suggested.

Journal of Geophysical Research B: Solid Earth↗

USGS2012

Sediment entrainment by debris flows: In situ measurements from the headwaters of a steep catchment

Debris flows can dramatically increase their volume, and hence their destructive potential, by entraining sediment. Yet quantitative constraints on rates and mechanics of sediment entrainment by debris flows are limited. Using an in situ sensor network in the headwaters of a natural catchment we measured flow and bed properties during six erosive debris-flow events. Despite similar flow properties and thicknesses of bed sediment entrained across all events, time-averaged entrainment rates were significantly faster for bed sediment that was saturated prior to flow arrival compared with rates for sediment that was dry. Bed sediment was entrained from the sediment-surface downward in a progressive fashion and occurred during passage of dense granular fronts as well as water-rich, inter-surge flow. En masse failure of bed sediment along the sediment-bedrock interface was never observed. Large-magnitude, high-frequency fluctuations in total normal basal stress were dissipated within the upper 5 cm of bed sediment. Within this near surface layer, concomitant fluctuations in Coulomb frictional resistance are expected, irrespective of the influence of pore fluid pressure or fluctuations in shear stress. If the near-surface sediment was wet as it was overridden by a flow, additional large-magnitude, high-frequency pore pressure fluctuations were measured in the near-surface bed sediment. These pore pressure fluctuations propagated to depth at subsonic rates and in a diffusive manner. The depth to which large excess pore pressures propagated was typically less than 10 cm, but scaled as ( D / f i ) 0.5 , in which D is the hydraulic diffusivity and f i is the frequency of a particular pore pressure fluctuation. Shallow penetration depths of granular-normal-stress fluctuations and excess pore pressures demonstrate that only near-surface bed sediment experiences the full dynamic range of effective-stress fluctuations, and as a result, can be more easily entrained than deeper sediment. These data provide robust tests for mechanical models of entrainment and demonstrate that a debris flow over wet bed sediment will be larger than the same flow over dry bed sediment.

Colorado↗

USGS2012

Hydrologic conditions controlling runoff generation immediately after wildfire

We investigated the control of postwildfire runoff by physical and hydraulic properties of soil, hydrologic states, and an ash layer immediately following wildfire. The field site is within the area burned by the 2010 Fourmile Canyon Fire in Colorado, USA. Physical and hydraulic property characterization included ash thickness, particle size distribution, hydraulic conductivity, and soil water retention curves. Soil water content and matric potential were measured indirectly at several depths below the soil surface to document hydrologic states underneath the ash layer in the unsaturated zone, whereas precipitation and surface runoff were measured directly. Measurements of soil water content showed that almost no water infiltrated below the ash layer into the near-surface soil in the burned site at the storm time scale (i.e., minutes to hours). Runoff generation processes were controlled by and highly sensitive to ash thickness and ash hydraulic properties. The ash layer stored from 97% to 99% of rainfall, which was critical for reducing runoff amounts. The hydrologic response to two rain storms with different rainfall amounts, rainfall intensity, and durations, only ten days apart, indicated that runoff generation was predominantly by the saturation-excess mechanism perched at the ash-soil interface during the first storm and predominantly by the infiltration-excess mechanism at the ash surface during the second storm. Contributing area was not static for the two storms and was 4% (saturation excess) to 68% (infiltration excess) of the catchment area. Our results showed the importance of including hydrologic conditions and hydraulic properties of the ash layer in postwildfire runoff generation models.

Colorado↗

USGS2012

A meeting of the waters: interdisciplinary challenges and opportunities in tidal rivers

At the interface of estuarine tides and freshwater rivers lie wetland and aquatic ecosystems, which experience dramatic effects of sea level rise. There, nontidal channels and riparian floodplains are transforming into tidal ecosystems, and tidal freshwater ecosystems are receiving increasing salinity. These river-floodplain systems have both fluvial characteristics, including meandering channels and expansive floodplain forests, and estuarine characteristics, including tides and intertidal wetlands [see Barendregt et al., 2009; Conner et al., 2007, and references therein]. Because tidal rivers lie at the disciplinary divide between fluvial and estuarine science, a knowledge gap has developed in scientists' understanding of the geomorphic and biogeochemical response of these environments to sea level rise, climate change, and anthropogenically driven variations in watershed exports.

Eos, Transactions, American Geophysical Union↗

USGS2012

Carbon dioxide and methane emissions from the Yukon River system

Carbon dioxide (CO 2 ) and methane (CH 4 ) emissions are important, but poorly quantified, components of riverine carbon (C) budgets. This is largely because the data needed for gas flux calculations are sparse and are spatially and temporally variable. Additionally, the importance of C gas emissions relative to lateral C exports is not well known because gaseous and aqueous fluxes are not commonly measured on the same rivers. We couple measurements of aqueous CO 2 and CH 4 partial pressures ( p CO 2 , p CH 4 ) and flux across the water-air interface with gas transfer models to calculate subbasin distributions of gas flux density. We then combine those flux densities with remote and direct observations of stream and river water surface area and ice duration, to calculate C gas emissions from flowing waters throughout the Yukon River basin. CO 2 emissions were 7.68 Tg C yr −1 (95% CI: 5.84 −10.46), averaging 750 g C m −2 yr −1 normalized to water surface area, and 9.0 g C m −2 yr −1 normalized to river basin area. River CH 4 emissions totaled 55 Gg C yr −1 or 0.7% of the total mass of C emitted as CO 2 plus CH 4 and ∼6.4% of their combined radiative forcing. When combined with lateral inorganic plus organic C exports to below head of tide, C gas emissions comprised 50% of total C exported by the Yukon River and its tributaries. River CO 2 and CH 4 derive from multiple sources, including groundwater, surface water runoff, carbonate equilibrium reactions, and benthic and water column microbial processing of organic C. The exact role of each of these processes is not yet quantified in the overall river C budget.

Yukon River Basin↗

USGS2012

Plate interaction in the NE Caribbean subduction zone from continuous GPS observations

Kinematic similarities between the Sumatra and Puerto Rico Trenches highlight the potential for a mega-earthquake along the Puerto Rico Trench and the generation of local and trans-Atlantic tsunamis. We used the horizontal components of continuous GPS (cGPS) measurements from 10 sites on NE Caribbean islands to evaluate strain accumulation along the North American (NA) - Caribbean (CA) plate boundary. These sites move westward and slightly northward relative to CA interior at rates ≤2.5 mm/y. Provided this motion originates in the subduction interface, the northward motion suggests little or no trench-perpendicular thrust accumulation and may in fact indicate divergence north of Puerto Rico, where abnormal subsidence, bathymetry, and gravity are observed. The Puerto Rico Trench, thus, appears unable to generate mega-earthquakes, but damaging smaller earthquakes cannot be discounted. The westward motion, characterized by decreasing rate with distance from the trench, is probably due to eastward motion of CA plate impeded at the plate boundary by the Bahamas platform. Two additional cGPS sites in Mona Passage and SW Puerto Rico move to the SW similar to Hispaniola and unlike the other 10 sites. That motion relative to the rest of Puerto Rico may have given rise to seismicity and normal faults in Mona Rift, Mona Passage, and SW Puerto Rico.

Geophysical Research Letters↗

USGS2012

Juan de Fuca slab geometry and its relation to Wadati-Benioff zone seismicity

A new model of the subducted Juan de Fuca plate beneath western North America allows first-order correlations between the occurrence of Wadati-Benioff zone earthquakes and slab geometry, temperature, and hydration state. The geo-referenced 3D model, constructed from weighted control points, integrates depth information from earthquake locations and regional seismic velocity studies. We use the model to separate earthquakes that occur in the Cascadia forearc from those that occur within the underlying Juan de Fuca plate and thereby reveal previously obscured details regarding the spatial distribution of earthquakes. Seismicity within the slab is most prevalent where the slab is warped beneath northwestern California and western Washington suggesting that slab flexure, in addition to expected metamorphic dehydration processes, promotes earthquake occurrence within the subducted oceanic plate. Earthquake patterns beneath western Vancouver Island are consistent with slab dehydration processes. Conversely, the lack of slab earthquakes beneath western Oregon is consistent with an anhydrous slab. Double-differenced relocated seismicity resolves a double seismic zone within the slab beneath northwestern California that strongly constrains the location of the plate interface and delineates a cluster of seismicity 10 km above the surface that includes the 1992 M7.1 Mendocino earthquake. We infer that this earthquake ruptured a surface within the Cascadia accretionary margin above the Juan de Fuca plate. We further speculate that this earthquake is associated with a detached fragment of former Farallon plate. Other subsurface tectonic elements within the forearc may have the potential to generate similar damaging earthquakes.

British Columbia;California;Oregon;Washington↗

USGS2012

Drivers of change in estuarine-coastal ecosystems: Discoveries from four decades of study in San Francisco Bay

Poised at the interface of rivers, ocean, atmosphere and dense human settlement, estuaries are driven by a large array of natural and anthropogenic forces. San Francisco Bay exemplifies the fast-paced change occurring in many of the world's estuaries, bays and inland seas in response to these diverse forces. We use observations from this particularly well-studied estuary to illustrate responses to six drivers that are common agents of change where land and sea meet: water consumption and diversion; human modification of sediment supply; introduction of non-native species; sewage input; environmental policy; and climate shifts. In San Francisco Bay, responses to these drivers include, respectively, shifts in the timing and extent of freshwater inflow and salinity intrusion; decreasing turbidity; restructuring of plankton communities; nutrient enrichment; elimination of hypoxia and reduced metal contamination of biota; and food web changes that decrease resistance of the estuary to nutrient pollution. Detection of these changes and discovery of their causes through environmental monitoring have been essential for establishing and measuring outcomes of environmental policies that aim to maintain high water quality and sustain services provided by estuarine-coastal ecosystems. The wide range of variability time scales and the multiplicity of interacting drivers place heavy demands on estuarine monitoring programs. But the San Francisco Bay case study illustrates why the imperative for monitoring has never been greater.

California↗

USGS2019

The 12 November 2017 Mw 7.3 Ezgeleh–Sarpolzahab (Iran) earthquake and active tectonics of the Lurestan arc

The 12 November 2017 M w 7.3 Ezgeleh‐Sarpolzahab earthquake is the largest instrumentally recorded earthquake in the Zagros Simply Folded Belt by a factor of ∼10 in seismic moment. Exploiting local, regional, and teleseismic data and synthetic aperture radar interferometry imagery, we characterize the rupture, its aftershock sequence, background seismicity, and regional tectonics. The mainshock ruptured slowly (∼2 km/s), unilaterally southward, for ∼40 km along an oblique (dextral‐thrust) fault that dips ∼14°E beneath the northwestern Lurestan arc. Slip is confined to basement depths of ∼12–18 km, resolvably beneath the sedimentary cover which is ∼8 km thick in this area. The gentle dip angle and basement location allow for a broad slip area, explaining the large magnitude relative to earthquakes in the main Fars arc of the Zagros, where shallower, steeper faults are limited in rupture extent by weak sedimentary layers. Early aftershocks concentrate around the southern and western edges of the mainshock slip area and therefore cluster in the direction of rupture propagation, implying a contribution from dynamic triggering. A cluster of events ∼100 km to the south near Mandali (Iraq) reactivated the ∼50° dipping Zagros Foredeep Fault. The basement fault responsible for the Ezgeleh‐Sarpolzahab earthquake probably accounts for the ∼1 km elevation contrast between the Lurestan arc and the Kirkuk embayment but is distinct from sections of the Mountain Front Fault that define frontal escarpments elsewhere in the Zagros. It may be related to a seismic interface underlying the central and southern Lurestan arc, and a key concern is whether or not the more extensive regional structure is also seismogenic.

Journal of Geophysical Research: Solid Earth↗

USGS2019

Prediction of ice‐free conditions for a perennially ice‐covered Antarctic lake

Although perennially ice‐covered Antarctic lakes have experienced variable ice thicknesses over the past several decades, future ice thickness trends and associated aquatic biological responses under projected global warming remain unknown. Heat stored in the water column in chemically stratified Antarctic lakes that have middepth temperature maxima can significantly influence the ice thickness trends via upward heat flux to the ice/water interface. We modeled the ice thickness of the west lobe of Lake Bonney, Antarctica, based on possible future climate scenarios utilizing a 1D thermodynamic model that accounts for surface radiative fluxes as well as the heat flux associated with the temperature evolution of the water column. Model results predict that the ice cover of Lake Bonney will shift from perennial to seasonal within one to four decades, a change that will drastically influence ecosystem processes within the lake.

Journal of Geophysical Research F: Earth Surface↗

USGS2019

Seasonal and spatial variation in the location and reactivity of a nitrate-contaminated groundwater discharge zone in a lakebed

Groundwater discharge delivering anthropogenic N from surrounding watersheds can impact lake nutrient budgets. However, upgradient groundwater processes and changing dynamics in N biogeochemistry at the groundwater-lake interface are complex and difficult to resolve. In this study, hydrograph variations in a groundwater flow-through lake altered discharge patterns of a wastewater-derived, groundwater contaminant plume, thereby affecting biogeochemical processes controlling N transport. Groundwater geochemistry 15 cm under the lakebed along transects perpendicular to shore varied from oxic to anoxic with increasing nitrate concentrations (10-75 M) and corresponding gradients in nitrite and nitrous oxide. Porewater depth profiles of nitrate concentrations and stable isotope compositions largely reflected upgradient groundwater N sources and N-cycle processes, with minor additional nitrate reduction in the shallowest lakebed sediments. Potential denitrification rates determined in laboratory microcosms were 10-100 fold higher in near-surface sediments (0-5 cm) than in deeper sediments (5-30 cm) and were correlated with sediment carbon content and abundance of denitrification genes (nirS, nosZI, and nosZII). Potential anammox-driven N2 production was highest in deeper anoxic sediments. Injection of bromide and nitrite in the lake sediments indicated a vertical porewater velocity of 4-5 cm hr-1, with highest nitrite consumption rates above 10 cm. However, short residence times in the shallow sediments allowed only a small fraction of the contaminant nitrate to be removed before discharging into the lake. Results demonstrate the importance of resolving local versus upgradient biogeochemical processes affecting contaminant distribution in discharge areas, and transient migration of local gradients and processes in response to changing lake levels and groundwater flow paths.

Journal of Geophysical Research: Biogeosciences↗

USGS2019

Residence time controls on the fate of nitrogen in flow‐through lakebed sediments

For many glacial lakes with highly permeable sediments, water exchange rates control hydrologic residence times within the sediment‐water interface (SWI) and the removal of reactive compounds such as nitrate, a common pollutant in lakes and groundwater. Here we conducted a series of focused tracer injection experiments in the upper 20 cm of the naturally downwelling SWI in a flow‐through lake on Cape Cod, MA. We systematically varied residence time and reactant controls on nitrate processing, using isotopically labeled 15 N nitrate to monitor the effect of these changes on nitrate removal via denitrification. The addition of acetate, a labile carbon compound, triggered the lake SWI to switch from net production to net removal of nitrate. When acetate was combined with increased residence time created by controlled reductions in water flux, we observed a fivefold increase in nitrate removal, a 26‐fold increase in N 2 production, and a 42‐fold increase in N 2 O production. We demonstrate that water residence time is an important control on the fate of nitrate in these lake SWIs and illustrate that seasonal conditions that alter lake exchange rates and variability in lake carbon may predict dynamic nitrate removal across the SWI. Additionally, observed N 2 O production during the oxic pore water experiments paired with geophysical characterization of the sediment porosity revealed that the lake SWI has less mobile pores occupying upward of 50% of the total porosity volume, which function as reactive microzones for nitrate processing.

Journal of Geophysical Research: Biogeosciences↗