Dynamic modeling of nitrogen losses in river networks unravels the coupled effects of hydrological and biogeochemical processes

Richard B. Alexander; J.K. Bohlke; Elizabeth W. Boyer; Mark B. David; Judson W. Harvey; Patrick J. Mulholland; Sybil P. Seitzinger; Craig R. Tobias; Christina Tonitto; Wilfred M. Wollheim

doi:10.1007/s10533-008-9274-8

USGS · 70192572

Dynamic modeling of nitrogen losses in river networks unravels the coupled effects of hydrological and biogeochemical processes

Richard B. Alexander·J.K. Bohlke·Elizabeth W. Boyer·Mark B. David·Judson W. Harvey·Patrick J. Mulholland·Sybil P. Seitzinger·Craig R. Tobias·Christina Tonitto·Wilfred M. Wollheim

View original source ↗Export metadata ↓

Abstract

The importance of lotic systems as sinks for nitrogen inputs is well recognized. A fraction of nitrogen in streamflow is removed to the atmosphere via denitrification with the remainder exported in streamflow as nitrogen loads. At the watershed scale, there is a keen interest in understanding the factors that control the fate of nitrogen throughout the stream channel network, with particular attention to the processes that deliver large nitrogen loads to sensitive coastal ecosystems. We use a dynamic stream transport model to assess biogeochemical (nitrate loadings, concentration, temperature) and hydrological (discharge, depth, velocity) effects on reach-scale denitrification and nitrate removal in the river networks of two watersheds having widely differing levels of nitrate enrichment but nearly identical discharges. Stream denitrification is estimated by regression as a nonlinear function of nitrate concentration, streamflow, and temperature, using more than 300 published measurements from a variety of US streams. These relations are used in the stream transport model to characterize nitrate dynamics related to denitrification at a monthly time scale in the stream reaches of the two watersheds. Results indicate that the nitrate removal efficiency of streams, as measured by the percentage of the stream nitrate flux removed via denitrification per unit length of channel, is appreciably reduced during months with high discharge and nitrate flux and increases during months of low-discharge and flux. Biogeochemical factors, including land use, nitrate inputs, and stream concentrations, are a major control on reach-scale denitrification, evidenced by the disproportionately lower nitrate removal efficiency in streams of the highly nitrate-enriched watershed as compared with that in similarly sized streams in the less nitrate-enriched watershed. Sensitivity analyses reveal that these important biogeochemical factors and physical hydrological factors contribute nearly equally to seasonal and stream-size related variations in the percentage of the stream nitrate flux removed in each watershed.

Explore related subjects

Biogeochemistry Article

Keep this discovery

Explore connections, maps & timelines

BibTeX RIS

Richard B. Alexander, J.K. Bohlke, Elizabeth W. Boyer, Mark B. David, Judson W. Harvey, Patrick J. Mulholland, Sybil P. Seitzinger, Craig R. Tobias, Christina Tonitto, Wilfred M. Wollheim. 2009-01-07. Dynamic modeling of nitrogen losses in river networks unravels the coupled effects of hydrological and biogeochemical processes. https://doi.org/10.1007/s10533-008-9274-8

Cite the original work for its findings. Save a collection to share your selection of sources.

KEEP EXPLORING

Related USGS reports

USGS2026

Water table rise sustains carbon release from soils in wetland-dominated landscapes: An intact soil core study

Wetland-dominated landscapes influence carbon cycling through their potential to act as both carbon sinks and sources. Wetlands in low-relief landscapes have dynamic terrestrial-aquatic interfaces that change seasonally with variable surface water and groundwater levels. However, few studies have directly quantified dissolved organic matter (DOM) release and greenhouse gas (CO 2 , CH 4 ) fluxes from wetland soils along terrestrial-aquatic interfaces as they are seasonally re-saturated by groundwater. To estimate groundwater-mediated soil DOM and gas fluxes, we performed laboratory simulations of vertical groundwater rise on intact soil cores collected from four Delmarva bay wetlands located in the Mid-Atlantic United States. At each wetland, one core was collected from within the wetland basin and the other from the transitional zone at the basin edge. Cores were re-saturated with groundwater over 15 days and then kept fully saturated for an additional 25 days. Source groundwater, soil porewater, and exfiltrated surface water samples were collected and analyzed for pH, ORP, DOM concentration, and DOM optical indices. In both the wetland and transition zone cores, porewater DOM concentrations increased over the wet up and were sustained during prolonged saturation. Optical indices shifted from recently produced, microbial-like signatures towards aromatic, terrestrial-like signatures. Fluxes of CO 2 decreased as the duration of soil saturation increased and soil cores switched from CH 4 sinks to sources upon full soil core saturation. Results indicate that groundwater rise sustains carbon mobilization from soils in wetland-dominated landscapes, emphasizing the need to understand how climate-driven changes to groundwater dynamics may affect carbon fluxes along terrestrial-aquatic interfaces.

Delaware, Maryland↗

USGS2025

Seasonal synchronicity and multi-decadal stability of headwater biogeochemistry in the northern temperate zone

Temporal patterns in chemistry of headwater streams reflect responses of water and elemental cycles to perturbations occurring at local to global scales. We evaluated multi-scale temporal patterns in up to 32 y of monthly observations of stream chemistry (ammonium, calcium, dissolved organic carbon, nitrate, total dissolved phosphorus, and sulfate) in 22 reference catchments within the northern temperate zone of North America. Multivariate autoregressive state-space (MARSS) models were applied to quantify patterns at multi-decadal, seasonal, and shorter intervals during a period that encompassed warming climate, seasonal changes in precipitation, and regional declines in atmospheric deposition. Significant long-term trends in solute concentrations within a subset of the catchments were consistent with recovery from atmospheric deposition (e.g., calcium, nitrate, sulfate) and increased precipitation (e.g., dissolved organic carbon). Lack of evidence for multi-decadal trends in most catchments suggests resilience of northern temperate ecosystems or that subtle net effects of simultaneous changes in climate and disturbance regimes do not result in directional trends. Synchronous seasonal oscillations of solute concentrations occurred across many catchments, reflecting shared climate and biotic drivers of seasonality within the northern temperate zone. Despite shared patterns among catchments at a seasonal scale, multi-scale temporal patterns were statistically distinct among even adjacent headwater catchments, implying that local attributes of headwater catchments modify the signals imparted by atmospheric phenomena and regional disturbances. To effectively characterize hydrologic and biogeochemical responses to changing climate and disturbance regimes, catchment monitoring programs could include multiple streams with contributing areas that encompass regional heterogeneity in vegetation, topography, and elevation. Overall, detection of long-term patterns and trends requires monitoring multiple catchments at a frequency that captures periodic variation (e.g., seasonality) and a duration encompassing the perturbations of interest.

Biogeochemistry↗

USGS2025

Rare earth element and phosphorus mobility depend on adsorption to Al-, Fe-, and Mn-oxyhydroxides in a headwater stream in Vermont, USA

We evaluated base cations, Al, Fe, Mn, REEs, DOC, anions, and P mobilization during three discharge events (E-1, E-2, E-3) at the headwater catchment W-9, Sleepers River, Vermont, USA. Peak discharge ranged from 3.696 (E-1) to 0.073 (E-3) mm h −1 . Eight samples from each event were speciated for total (unfiltered-acidified) and dissolved (0.45 µm filtered-acidified). During E-1, total Al, Fe, and Mn increased to maxima of 376 (1790X), 161 (194X), and 38 (45X) µmol L −1 , respectively. Concurrently, total La, Ce, Pr, and Nd increased to maxima of 87 (590X), 114 (671X), 15 (375X), and 53 (408X) nmol L −1 , respectively, greatly exceeding the dissolved fraction. Totals for Er > Yb > Tm > Lu reached comparable enrichments near or at maximum discharge. Discharge ranged from 0.41 to 0.94 mm h −1 during E-2, a snowmelt event. Total Al was comparatively stable; total Fe and Mn increased 15X and 79X, remaining less than total Al. Total La, Ce, Pr, and Nd peaked with total Fe and Mn. E-3, a late summer rain, resembled E-1 but had much lower maxima for all REEs, P, Al, Fe, and Mn. Particulate, total, and dissolved REEs and P correlated with DOC and with total Al, Fe, and Mn over all discharges, with molar Al > Fe > Mn (most samples) for E-1 and E-3. Maximum total and dissolved P declined from 15.4 and 0.89 (E-1) to 0.98 and 0.27 µmol L −1 (E-3), respectively. Particulate REEs correlated strongly (R 2 = 0.95–0.96) with Al, Fe, and Mn particulates eroded from the stream bed and continuously precipitated during and after high discharge of groundwater.

Vermont↗