USGS Science⌕ Search

SEARCH · USGS Science

Results for “Journal of Geophysical Research - Atmospheres”

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

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

At least 793 records · Page 44Linked to original sources

Hydrodynamics of a tidally‐forced coral reef atoll

The hydrodynamics of a tidally forced semi‐enclosed coral reef atoll (North Scott) at the edge of the continental shelf of northwestern Australia were investigated by combining field observations and numerical modeling. The observations revealed that the spring tidal range outside the atoll reaches 4 m, and as the water level drops below mean sea level, the reef rim surrounding the shallow (~10–15 m) lagoon becomes exposed. During this time, the lagoon can only exchange with the open ocean through two narrow channels, resulting in highly asymmetric water levels and velocities that were most pronounced during spring tide. On average, the ebb tide duration was ~2 hr longer than the flood, with rapid flood velocities in the channel reaching 2 m/s. We applied an unstructured grid model Delft3D‐Flexible Mesh to simulate the atoll hydrodynamics and were able to replicate the asymmetric water levels and complex velocities in the lagoon. The results revealed that at higher tidal stages, a dominant momentum balance exists between the pressure gradient (established by the propagation of the tide on the shelf) and the local flow acceleration of water throughout the interior of the atoll. At lower tidal stages, which coincided with a reversal of the offshore tidal pressure gradient, the lagoon became isolated from offshore dynamics and all momentum terms were negligible. This resulted in a tidally averaged residual westward flow within the lagoon that drove an asymmetric flushing pattern within the atoll, which we propose would be a common flushing mechanism within other tide‐dominated atolls worldwide.

Journal of Geophysical Research C: Oceans↗

Field observations of alongshore runup variability under dissipative conditions in presence of a shoreline sandwave

Video measurements of runup were collected at low tide along several profiles covering an alongshore distance of 500 m. The morphology displayed a complex shape with a shoreline sandwave in the lower beach face of about 250 m long mirrored in the inner sandbar. Wave conditions were stationary and moderate (offshore height of 2 m and peak period of nearly 13 s) but yet dissipative. Runup energy was dominated by infragravity frequencies. Alongshore variations in runup (by a factor up to 3) observed both in the incident and infragravity bands were much higher than reported previously (e.g., Guedes et al., 2012, https://doi.org/10.1016/j.csr.2012.08.022 ; Ruggiero et al., 2004, https://doi.org/10.1029/2003JC002160 ) while the alongshore variations in other environmental parameters (e.g., foreshore beach slope) appear to be much lower. Our data suggest that the beach morphology in the inner surf zone plays a crucial role by inducing rapid and significant modification in the incident wave pattern and the alongshore coherence length scales were consistent with the typical alongshore length scale of the morphology.

Journal of Geophysical Research C: Oceans↗

Submarine permafrost map in the arctic modelled using 1D transient heat flux (SuPerMAP)

Offshore permafrost plays a role in the global climate system, but observations of permafrost thickness, state, and composition are limited to specific regions. The current global permafrost map shows potential offshore permafrost distribution based on bathymetry and global sea level rise. As a first‐order estimate, we employ a heat transfer model to calculate the subsurface temperature field. Our model uses dynamic upper boundary conditions that synthesize Earth System Model air temperature, ice mass distribution and thickness, and global sea level reconstruction and applies globally distributed geothermal heat flux as a lower boundary condition. Sea level reconstruction accounts for differences between marine and terrestrial sedimentation history. Sediment composition and pore water salinity are integrated in the model. Model runs for 450 ka for cross‐shelf transects were used to initialize the model for circumarctic modeling for the past 50 ka. Preindustrial submarine permafrost (i.e., cryotic sediment), modeled at 12.5‐km spatial resolution, lies beneath almost 2.5 ×10 6 km 2 of the Arctic shelf. Our simple modeling approach results in estimates of distribution of cryotic sediment that are similar to the current global map and recent seismically delineated permafrost distributions for the Beaufort and Kara seas, suggesting that sea level is a first‐order determinant for submarine permafrost distribution. Ice content and sediment thermal conductivity are also important for determining rates of permafrost thickness change. The model provides a consistent circumarctic approach to map submarine permafrost and to estimate the dynamics of permafrost in the past.

Journal of Geophysical Research C: Oceans↗

Seasonal, spring-neap, and tidal variation in cohesive sediment transport parameters in estuarine shallows

Numerical models for predicting sediment concentrations and transport rely on parameters such as settling velocity and bed erodibility that describe sediment characteristics, yet these parameters are rarely probed directly. We investigated temporal and spatial variation in sediment parameters in the shallows of San Pablo Bay, CA. Flow, turbulence, and suspended sediment data were measured at sites located at 1 and 2 m below mean lower low water (MLLW) from November 2013 through April 2015, supplemented by monthlong periods in 2011, 2012, and 2016. Maximum current velocities were 0.40–0.47 m s -1 at these depths; the strongest currents decreased to 0.27–0.34 m s -1 during neap periods. Winters 2013–2014 and 2014–2015 experienced strong drought conditions, limiting the potential for seasonal impact on sediment conditions during this experiment. Despite this, the more storm‐influenced site showed clear changes during the winter: the roughness parameter decreased from 10 −4 to 10 −5 m, from hydrodynamically rough to smooth conditions, and bed erodibility increased by an order of magnitude. Median settling velocity was 2.05·10 −4 m s -1 ; it varied twofold within a tidal cycle, decreasing as current velocity grew during flood and ebb. This tidal control on floc size affected settling velocity on the spring‐neap timescale, possibly driving a spring‐neap oscillation in erodibility. Our findings highlight variation in sediment dynamics that is commonly ignored in numerical models and the need for field observations to ground truth ongoing modeling efforts.

California↗

Degradation of 100‐m‐scale rocky ejecta craters at the InSight Landing Site on Mars and implications for surface processes and erosion rates in the hesperian and amazonian

Rocky ejecta craters (RECs) at the Interior Exploration using Seismic Investigations, Geodesy and Heat Transport (InSight) landing site on Elysium Planitia, Mars, provide constraints on crater modification and rates for the Hesperian and Amazonian. The RECs are between 10 m and 1.2 km in diameter and exhibit five classes of preservation. Class 1 represents pristine craters with sharp rims and abundant ejected rocks. From Classes 2 to 5, rims become more subdued, craters are infilled, and the ejecta become discontinuously distributed. High‐Resolution Imaging Science Experiment digital elevation models indicate a maximum depth to diameter ratio of ~0.15, which is lower than pristine models for craters of similar size. The low ratio is related to the presence of a loosely consolidated regolith and early‐stage eolian infill. Rim heights have an average height to diameter ratio of ~0.03 for the most pristine class. The size‐frequency distribution of RECs, plotted using cumulative and differential methods, indicates that crater classes within the diameter range of 200 m to 1.2 km are separated by ~100 to 200 Myr. Smaller craters degrade faster, with classes separated by <100 Myr. Rim erosion can be entirely modeled by nonlinear diffusional processes using the calculated timescales and a constant diffusivity of 8 × 10 −7 m 2 /year for craters 200 to 500 m in diameter. Diffusion models only partly capture depth‐related degradation, which requires eolian infill. Depth degradation and rim erosion rates are 10 −2 to 10 −3 m/Myr, respectively. The rates are consistent with relatively slow modification that is typical of the last two epochs of Martian history.

Journal of Geophysical Research - Planets↗

Floor-fractured craters on Ceres and implications for interior processes

Several of the impact craters on Ceres have sets of fractures on their floors. These fractures appear similar to those found within a class of lunar craters referred to as floor-fractured craters (FFCs). We have cataloged the Ceres FFCs according to the classification scheme designed for the Moon. An analysis of the depth to diameter ratio for Ceres craters shows that, like lunar FFCs, the Ceres FFCs are anomalously shallow. Large (>50 km) Ceres FFCs are most consistent with Class 1 lunar FFCs, while smaller craters on Ceres are more consistent with Class 4 lunar FFCs. This suggests that Ceres FFCs may similarly be undergoing fracturing due to the intrusion of a low-density material below the craters. While on the Moon (and Mars) the intrusive material is hypothesized to be silicate magma, cryomagmatic intrusions are more likely responsible for the formation of the Ceres FFCs. However, new models suggest that at least some of the FFC fractures may have formed due to the solid-state flow of a low-viscosity and low-density material into the crater wall.

Journal of Geophysical Research↗

Overview of spirit microscopic imager results

This paper provides an overview of Mars Exploration Rover Spirit Microscopic Imager (MI) operations and the calibration, processing, and analysis of MI data. The focus of this overview is on the last five Earth years (2005–2010) of Spirit 's mission in Gusev crater, supplementing the previous overview of the first 450 sols of the Spirit MI investigation. Updates to radiometric calibration using in‐flight data and improvements in high‐level processing are summarized. Released data products are described, and a table of MI observations, including target/feature names and associated data sets, is appended. The MI observed natural and disturbed exposures of rocks and soils as well as magnets and other rover hardware. These hand‐lens‐scale observations have provided key constraints on interpretations of the formation and geologic history of features, rocks, and soils examined by Spirit . MI images complement observations by other Spirit instruments, and together show that impact and volcanic processes have dominated the origin and evolution of the rocks in Gusev crater, with aqueous activity indicated by the presence of silica‐rich rocks and sulfate‐rich soils. The textures of some of the silica‐rich rocks are similar to terrestrial hot spring deposits, and observations of subsurface cemented layers indicate recent aqueous mobilization of sulfates in places. Wind action has recently modified soils and abraded many of the rocks imaged by the MI, as observed at other Mars landing sites.

Journal of Geophysical Research E: Planets↗

Submarine deposition of a subaerial landslide in Taan Fiord, Alaska

A large subaerial landslide entered Taan Fiord, Alaska, on 17 October 2015 producing a tsunami with runup to 193 m. We use LiDAR data to show the slide volume to be 76 + 3/−4 million cubic meters and that 51,000,000 m 3 entered Taan Fiord. In 2016, we mapped the fjord with multibeam bathymetry and high‐resolution seismic data. Landslide and postlandslide deposits extend 6 km downfjord, are up to 70 ± 11 m thick, and have a total volume of ~147,000,000 m 3 . Seismic data image a blocky landslide unit and two units deposited immediately after the landslide. The blocky landslide unit is ~65,000,000 m 3 . We infer it consists dominantly of subaerially derived material and secondarily of fjord floor sediment. The overlying units are likely megaturbidites presumably deposited within minutes to days after the landslide. We infer that these deposits dominantly consist of fjord floor material mobilized and suspended as the slide entered and traveled downfjord. The lower postlandslide unit is up to 35 ± 6 m thick, and the upper unit is up to 12 ± 3 m thick. These deposits are distinctive and will leave a lasting record of the event. This subaerial‐to‐submarine landslide deposit is distinct from other submarine landslide deposits studied in Alaskan fjords because it has a much greater thickness, larger and more angular blocks, distinctive postlandslide megaturbidites, and a higher‐amplitude acoustic signature of the blocky deposit. The tight constraints on the landslide source and deposit volumes, topography, bathymetry, and tsunami runup heights and flow directions should make this a benchmark site for landslide‐tsunami models.

Alaska↗

Geomorphic evolution of a gravel‐bed river under sediment‐starved vs. sediment‐rich conditions: River response to the world's largest dam removal

Understanding river response to sediment pulses is a fundamental problem in geomorphic process studies, with myriad implications for river management. However, because large sediment pulses are rare and usually unanticipated, they are seldom studied at field scale. We examine fluvial response to a massive (~20 Mt) sediment pulse released by the largest dam removal globally, on the Elwha River, Washington, United States, in an 11‐year before‐after/control‐impact study of channel morphology and grain size. We test the hypothesis that for a given flow magnitude, greater geomorphic change occurs under sediment‐rich conditions than under sediment‐starved conditions. Channel response to flow forcing was significantly different during the sediment‐pulse peak, 1–2 years after dam removal began, than earlier or later. During peak sediment supply our hypothesis was supported; major geomorphic change occurred under low flows and unit stream power ≤60 W/m 2 . However, by 4–6 years after dam removal began, rates of geomorphic change and sensitivity to stream power had decreased substantially such that our hypothesis was no longer unequivocally supported. These findings are consistent with a two‐phase conceptual model of dam‐removal response, involving a transport‐limited state followed by a more supply‐limited state. From comparisons with other dam removals and natural sediment pulses, we infer that the longevity of sediment‐pulse signals in gravel‐bed rivers depends upon gradient, river discharge, valley morphology, and sediment grain size. Stream power associated with substantial geomorphic change varies with sediment supply, such that assigning a general threshold stream power to gravel‐bed rivers may be untenable.

Washington↗

Downstream‐propagating channel responses to decadal‐scale climate variability in a glaciated river basin

Regional climate is an important control on the rate of coarse sediment mobilization and transport in alpine river systems. Changes in climate are then expected to cause a cascade of geomorphic responses, including adjustments in downstream channel morphology. However, the mechanics and sensitivity of channel response to short‐term climate variability remain poorly documented. In the Nooksack River, which drains a glaciated stratovolcano in Washington State, bed elevation changes were inferred from shifting stage–discharge relations at seven USGS stream gages. Decadal‐scale elevation trends at most sites can be explained as a downstream‐propagating channel response to regional climate variability, where periods of persistent warm, dry [cool, wet] conditions corresponded to periods of aggradation [incision]. The channel elevation response propagated downstream at a rate of one to four kilometers per year; propagation rate scaled closely with channel slope. Historical trends in glacier extent and flood intensity both show some potential to explain climate–sediment linkages, though assessing causation is complicated by the shared climate signal in both records. Results show the influence of the Pacific Decadal Oscillation, with relatively high coarse sediment yields prior to 1950 and since 1980, and notably lower sediment yields from 1950 to 1980. Measured sediment yields from nearby glaciated basins corroborate this history, suggesting a regional coherence to these climate–sediment linkages. These results document consistent relations between climate, sediment supply and downstream channel response at the basin‐scale, with channel responses propagating downstream over periods of decades with little apparent attenuation.

Washington↗

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↗

Multidecadal geomorphic evolution of a profoundly disturbed gravel-bed river system—a complex, nonlinear response and its impact on sediment delivery

A 2.5-km3 debris avalanche during the 1980 eruption of Mount St. Helens reset the fluvial landscape of upper North Fork Toutle River valley. Since then, a new drainage network has formed and evolved. Cross-section surveys repeated over nearly 40 years at 16 locations along a 20-km reach of river valley document channel evolution, geomorphic processes, and their impacts on sediment delivery. We analyzed spatial and temporal changes in channel morphology using two new metrics: 1) a shape index that defines the degree of U-shaped or V-shaped valley geometry; and 2) an alluvial phase-space diagram that relates bed degradation or aggradation between consecutive surveys to increases or decreases in cross-section area. Metric relations reveal more diverse channel evolution than originally described by a simple, linear-response model of sequential channel initiation and incision; aggradation and widening; and subsequent episodic scour and fill with little change in bed elevation. Instead, vertical and lateral adjustments have been crucial processes intertwined throughout channel evolution. Channel evolution has followed a distinctly nonlinear and non-sequential trajectory, migrating through several phase spaces and involving varied combinations of (1) degradation and aggradation with widening and narrowing, (2) bed-level fluctuations with little change in cross-section area, and (3) changes in cross-section area with little change of bed elevation. Persistent channel widening and reworking of the channel bed presently drive elevated sediment delivery from this basin. Elevated sediment delivery is likely to persist until valley-floor widths greatly exceed that of the channel-migration corridor, and/or channel banks and valley walls stabilize.

Washington↗

Ice wedge degradation and stabilization impacts water budgets and nutrient cycling in Arctic trough ponds

Trough ponds are ubiquitous features of Arctic landscapes and an important component of freshwater aquatic ecosystems. Permafrost thaw causes ground subsidence, creating depressions that gather water, creating ponds. Permafrost thaw also releases solutes and nutrients, which may fertilize these newly formed ponds. We measured water budget elements and chloride, ammonium, and dissolved organic nitrogen (DON) across a chronosequence of trough ponds representing different stages of ice wedge degradation and stabilization. We developed a coupled hydrologic and biogeochemical model to explore how ice wedge degradation affects hydrology and nutrient availability in trough ponds in the advanced degradation stages (DAs), which are characterized by deep troughs with warmer temperatures relative to the other stages. DAs experienced greater evaporation than the other stages, and subsurface inflows entered the DAs from a wide area. Chloride accumulated in the ponds with time since thaw, implying that subsurface fluxes are delivering solutes from the thawing permafrost. Ammonium accumulated at high rates in the initial degradation stage and was seasonally depleted over the summer in all degradation stages. Ammonium trends in the DAs were consistent with high concentration inflows and in‐pond assimilation at rates between 0.37 and 2.0 mg N m −2 day −1 . Seasonal DON trends indicated that the accumulation of recalcitrant organic matter may eventually limit aquatic ecosystem production and foster pond infilling. These results provide direct evidence of nutrient release from thawing permafrost and the utilization of these nutrients by Arctic trough pond ecosystems and highlight infilling as a mechanism by which Arctic surface waters may be lost

Journal of Geophysical Research: Biogeosciences↗

Sediment trapping and carbon sequestration in floodplains of the lower Atchafalaya Basin, LA: Allochthonous vs. autochthonous carbon sources

Recent studies suggest that about 2 Pg of organic C is stored on floodplains worldwide. The present study indicates the Atchafalaya River, fifth largest river in the United States in terms of discharge, traps 30 mm/y of sediment on average within its floodplain, which is the highest average non‐episodic rate of fluvial deposition on the U.S. Coastal Plain. We installed sediment sampling stations at 23 sites, normally in transect, in the Atchafalaya Basin; these sites represent the range of hydrogeomorphic conditions on the floodplain based on hydrologic connectivity with the river main stem. The rate of sedimentation translates into about 12.5 Tg/y and includes 694 Mg/y of organic C. Highest sedimentation rates are associated with areas of high connectivity to channels and prograding deltaic processes. The δ 13 C content suggests that 35% of deposited C is derived from river‐suspended sediment compared to litterfall in the Basin. Thus, much of the organic C sequestered is allochthonous material. However, isolated interior sites with limited connectivity to the channel may generate and sequester large amounts of autochthonous C. The substantial trapping of both auto‐ and allochthonous C (392 Mg/y) make this freshwater‐forested floodplain critical in storage of material before reaching the coastal delta and estuary. This C deposition rate (340 g C m ‐2 y ‐1 ) exceeds all other rates reported in recent Blue Carbon and Tidal Freshwater Forested Wetland studies. Atchafalaya C sequestration occurs in/near areas with tidal influence and like other coastal systems is an important site for trapping mineral and organic sediment and in global C cycling.

Louisiana↗

Estimating soil respiration in a subalpine landscape using point, terrain, climate and greenness data

Landscape carbon (C) flux estimates are necessary for assessing the ability of terrestrial ecosystems to buffer further increases in anthropogenic carbon dioxide (CO2) emissions. Advances in remote sensing have allowed for coarse-scale estimates of gross primary productivity (GPP) (e.g., MODIS 17), yet efforts to assess spatial patterns in respiration lag behind those of GPP. Here, we demonstrate a method to predict growing season soil respiration at a regional scale in a forested ecosystem. We related field measurements (n=144) of growing season soil respiration across subalpine forests in the Southern Rocky Mountains ecoregion to a suite of biophysical predictors with a Random Forest model (30 m pixel size). We found that Landsat Enhanced Vegetation Index (EVI), growing season AI, temperature, precipitation, elevation, and slope aspect explained spatiotemporal variability in soil respiration. Our model had a psuedo-r2 of 0.45 and root mean squared error (RMSE) of roughly one-quarter of the mean value of respiration. Predicted growing season soil respiration across the region was remarkably consistent across 2004, 2005 and 2006 (150-d averages of 542.8, 544.3, and 536.5 g C m-2, respectively). Yet, we observed substantial variability in spatial patterns of soil respiration predictions that varied between years, suggesting that our method is sensitive to changes in respiration drivers. We compared our estimates to MODIS GPP and nocturnal net ecosystem exchange (NEE) derived from eddy covariance towers as a proxy for ecosystem respiration. Averaged across the predictive region, mean predicted growing season soil respiration was 73% of MODIS GPP, while predicted soil respiration was generally within 20% of nocturnal NEE from eddy covariance towers. This study demonstrated that geospatial and remotely-sensed datasets can be used in a statistical modeling framework to estimate soil respiration at landscape scales.

Colorado, Wyoming↗

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↗

El Niño increases high‐tide flooding in tidal wetlands along the U.S. Pacific coast.

Periodic oscillations between El Niño and La Niña conditions in the Pacific Basin affect oceanographic and meteorological phenomena globally, with impacts on the abundance and distribution of marine species. However, El Niño effects on estuarine hydrology and tidal wetland processes have seldom been examined rigorously. We used detailed wetland elevation and local inundation data from 10 tidal wetlands located along the Pacific coast of the United States to assess changes in flooding during the 2015–2016 El Niño and to determine decadal‐scale relationships between estuarine sea‐level anomalies and Pacific Basin climate indices for this region. During the 2015–2016 El Niño all sites experienced significant increases in high‐tide water levels exceeding those predicted by astronomical tides, and increased flooding frequency during at least one of the El Niño subperiods relative to pre‐El Niño conditions. The magnitude of positive sea‐level anomalies varied by site (4–15 cm), with local hot spots of high water in southern Oregon, northern California, and Pt. Mugu lagoon in the Southern California Bight. Furthermore, over the last three decades of historic tide records, there were positive relationships between high‐tide sea‐level anomalies and equatorial Pacific Basin sea surface temperature anomalies across the region, and negative relationships with the Northern Oscillation Index. Increases of 1 °C in equatorial sea surface temperature were associated with 3–5 cm of increased high‐tide flooding at the sites. Elevated estuarine flooding associated with future El Niños could impact important tidal wetland processes and could be an additive stressor for wetlands facing accelerating sea‐level rise.

California, Oregon, Washington↗

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↗