Abstract. Solute concentrations in stream water vary with discharge in patterns that record complex feedbacks between hydrologic and biogeochemical processes. In a comparison of three shale-underlain headwater catchments located in Pennsylvania, USA (the forested Shale Hills Critical Zone Observatory), and Wales, UK (the peatland-dominated Upper Hafren and forest-dominated Upper Hore catchments in the Plynlimon forest), dissimilar concentration–discharge (C–Q) behaviors are best explained by contrasting landscape distributions of soil solution chemistry – especially dissolved organic carbon (DOC) – that have been established by patterns of vegetation and soil organic matter (SOM). Specifically, elements that are concentrated in organic-rich soils due to biotic cycling (Mn, Ca, K) or that form strong complexes with DOC (Fe, Al) are spatially heterogeneous in pore waters because organic matter is heterogeneously distributed across the catchments. These solutes exhibit non-chemostatic behavior in the streams, and solute concentrations either decrease (Shale Hills) or increase (Plynlimon) with increasing discharge. In contrast, solutes that are concentrated in soil minerals and form only weak complexes with DOC (Na, Mg, Si) are spatially homogeneous in pore waters across each catchment. These solutes are chemostatic in that their stream concentrations vary little with stream discharge, likely because these solutes are released quickly from exchange sites in the soils during rainfall events. Furthermore, concentration–discharge relationships of non-chemostatic solutes changed following tree harvest in the Upper Hore catchment in Plynlimon, while no changes were observed for chemostatic solutes, underscoring the role of vegetation in regulating the concentrations of certain elements in the stream. These results indicate that differences in the hydrologic connectivity of organic-rich soils to the stream drive differences in concentration behavior between catchments. As such, in catchments where SOM is dominantly in lowlands (e.g., Shale Hills), we infer that non-chemostatic elements associated with organic matter are released to the stream early during rainfall events, whereas in catchments where SOM is dominantly in uplands (e.g., Plynlimon), these non-chemostatic elements are released later during rainfall events. The distribution of SOM across the landscape is thus a key component for predictive models of solute transport in headwater catchments.
Concentration–discharge (C‐Q) relationships reflect material sources, storage, reaction, proximity, and transport in catchments. Differences in hydrologic pathways and connectivity influence observed C‐Q patterns at the catchment outlet. We examined solute and sediment C‐Q relationships at event and interannual timescales in a small mid‐Atlantic (USA) catchment. We found systematic differences in the C‐Q behaviour of geogenic/exogenous solutes (e.g., calcium and nitrate), biologically associated solutes (e.g., dissolved organic carbon), and particulate materials (e.g., total suspended solids). Negative log(C)–log(Q) regression slopes, indicating dilution, were common for geogenic solutes whereas positive slopes, indicating concentration increase, were common for biologically associated solutes. Biologically associated solutes often exhibited counterclockwise hysteresis during events whereas geogenic solutes exhibited clockwise hysteresis. Across event and interannual timescales, solute C‐Q patterns are linked to the spatial distribution of hydrologic sources and the timing and sequence of hydro‐biogeochemical source contributions to the stream. Groundwater is the primary source of stormflow during the earliest and latest stages of events, whereas precipitation and soil water become increasingly connected to the stream near peakflow. This sequence and timing of flowpath connectivity results in dilution and clockwise hysteresis for geogenic/exogenous solutes and concentration increase and counterclockwise hysteresis for biologically associated solutes. Particulate materials demonstrated positive C‐Q slopes over the long‐term and clockwise hysteresis during individual events. Drivers of particulate and solute C‐Q relationships differ, based on longitudinal and lateral expansion of active channels and changing shear stresses with increasing flows. Although important distinctions exist between the drivers of solute and sediment C‐Q relationships, overall solute and sediment C‐Q patterns at event and interannual timescales reflect consistent catchment hydro‐biogeochemical processes.
more » « less- NSF-PAR ID:
- 10066269
- Publisher / Repository:
- Wiley Blackwell (John Wiley & Sons)
- Date Published:
- Journal Name:
- Hydrological Processes
- Volume:
- 32
- Issue:
- 18
- ISSN:
- 0885-6087
- Page Range / eLocation ID:
- p. 2829-2844
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
-
more » « less
Abstract High‐frequency in situ sensors have enabled researchers to measure solute concentrations at a time scale that captures the variability in stream discharge. We analyzed discrete samples and high‐frequency time series of solutes to characterize how nitrate (NO3−) and fluorescent dissolved organic matter (fDOM; a proxy for dissolved organic carbon) respond to changes in discharge at annual and intra‐annual timescales across a stream network in New Hampshire, USA. NO3−and fDOM exhibited highly variable concentration‐discharge (c‐Q) behavior at intra‐annual scales. Transport limitation, source limitation, and chemostatic behavior were observed to occur within and among years in all our study watersheds. Annual assessment of c‐Q misclassified streams 31% of the time, as the annual time step missed seasonal and event‐induced shifts in c‐Q dynamics. In some instances, anomalous events lasting less than 5% of the year determine the annual c‐Q behavior for a site. Catchment land use appeared to drive some of the variability among watersheds in c‐Q relationships and their temporal variability. Forested streams had highly variable NO3−c‐Q behavior and streams draining watersheds with more development had greater variability in fDOM c‐Q behavior. Sample frequency impacts how hydrologic systems are characterized and extrapolating c‐Q behavior from discrete samples alone can bias interpretations of c‐Q dynamics and our understanding of solute transport.
-
High elevation mountain watersheds are undergoing rapid warming and declining snow fractions worldwide, causing earlier and quicker snowmelt. Understanding how this hydrologic shift affects subsurface flow paths, biogeochemical reactions, and solute export has been challenging due to the entanglement of hydrological and biogeochemical processes. Coal Creek, a high-elevation catchment (2,700 3,700 m, 53 km2) in Colorado, is experiencing a higher rate of warming than surrounding low-lying areas. This warming corresponds with dynamic and increased responses from biogenic solutes and dissolved organic carbon (DOC), whereas the behavior of geogenic solutes and dissolved inorganic carbon (DIC) has remained relatively unchanged. DOC has experienced the largest concentration increase (>3x), with annual average flow weighted concentrations positively correlated to average annual temperature. This suggests temperature is the main driver of increasing DOC levels. Although DOC and DIC response to warming is influenced by many drivers, the relative contribution of each remains unknown. DOC and DIC were analyzed to incorporate both carbon component products of soil respiration (DOC and CO2) and to represent high solute concentrations transported by shallow (DOC) versus deep (DIC) subsurface flow. The contrasting behavior of these carbon solutes indicates climate change and warming are driving changes in organic matter decomposition and soil respiration. Modeling results from the process-based model HBV-BioRT show increased temperatures cause earlier snowmelt and streamflow generation and lower peak discharge. As stream flow generation occurs earlier, so do DOC flushing and DIC dilution events. Additionally, post-snowmelt periods show greater DOC production and concentrations under warming scenarios. Results indicated increased production of DOC in post-snowmelt periods. DOC is then flushed out by earlier snowmelt partitioned through the shallow soil zone. Most process-based studies lack a watershed-scale understanding of carbon transformation and flow path alterations. This work demonstrates complex hydrologic and biogeochemical coupling at the watershed scale to illustrate how water flow paths and chemistry are responding to a changing climate in highelevation mountain watersheds.more » « less
-
more » « less
Abstract Climate change is intensifying the Arctic hydrologic cycle, potentially accelerating the release of carbon and nutrients from permafrost landscapes to rivers. However, there are limited riverine flow and solute data of adequate frequency and duration to test how seasonality and catchment landscape characteristics influence production and transport of carbon and nutrients in Arctic river networks. We measured high frequency hydrochemical dynamics at the outlets of three headwater catchments in Arctic Alaska over 3 years. The catchments represent common Arctic landscapes: low‐gradient tundra, low‐gradient and lake‐influenced tundra, and high‐gradient alpine tundra. Using in‐situ spectrophotometers, we measured dissolved organic carbon (DOC) and nitrate (NO3−) concentrations at 15‐min intervals through the flow seasons of 2017, 2018, and 2019. These high‐frequency data allowed us to quantify concentration–discharge (C‐Q) responses during individual storm events across the flow season. Differences in C‐Q responses among catchments indicated strong landscape and seasonal controls on lateral DOC and NO3−flux. For the two low‐gradient tundra catchments, we observed consistent DOC enrichment (transport‐limitation) and NO3−dilution (source‐limitation) during flow events. Conversely, we found consistent NO3−enrichment and DOC dilution in the high‐gradient alpine catchment. Our analysis revealed how high flow events may contribute disproportionately to downstream export in these Arctic streams. Because the duration of the flow season is expected to lengthen and the intensity of Arctic storms are expected to increase, understanding how discharge and solute concentration are coupled is crucial to understanding carbon and nutrient dynamics in rapidly changing permafrost ecosystems.
-
This study examines dissolved rhenium (Re) concentrations as a function of water runoff using river samples from two contrasting mountainous watersheds, the Eel and Umpqua Rivers in the Pacific Northwest, USA. These watersheds share many key characteristics in terms of size, discharge, climate, and vegetation, but they have a 15-fold difference in sediment yield due to differences in their tectonic setting and uplift and erosion rates. We evaluate concentration-runoff (C-R) relationships and ratios of coefficients of variation (CVC/CVR) for major cations, anions, dissolved inorganic carbon, selected trace elements including Re, and 87Sr/86Sr ratios. Recent research outlines the potential of Re to serve as a tracer for the oxidation of ancient/fossil organic matter because of its close association with petrogenic carbon (OCpetro) in rocks. In both the Eel and Umpqua Rivers, our measurements show that Re behaves similarly to major weathering derived-solutes corrected for atmospheric input, such as Ca2+*, Mg2+*, and Na+* with modest dilution across all tributaries with increasing runoff. Rhenium behaves dissimilarly from other trace elements, such as Mo and U, and is also dissimilar to biologically-cycled nutrients, such as NO3 – , PO4 3 , and K+*, suggesting differences in sources, solute generation mechanisms, and flowpaths. Rhenium behavior is also distinct from that of colloids, which have increasing concentrations with increasing runoff. We find that Re and sulfate corrected for atmospheric input (SO4 2 *) have distinct CR relationships, in which SO4 2 * undergoes greater dilution with increasing runoff. This implies that Re is not dominantly sourced from sulfide weathering, which leaves primary bedrock minerals and OCpetro hosted in bedrock of these watersheds as the likely dominant sources of dissolved Re release. At mean discharge, Re concentration in the Eel river (3.5 pmol L-1) is more than two times greater than Re concentrations in the Umpqua River (1.5 pmol L-1). Furthermore, comparison of two tributary watersheds with similar bedrock but marked differences in erosion rates show higher Re concentrations in Bull Creek (erosion rate of 0.5 mm yr 1) relative to Elder Creek (erosion rate of 0.2 mm yr 1). The results of this study suggest that dissolved Re in the Eel and Umpqua River basins is likely derived from primary mineral dissolution or OCpetro oxidation, and Re fluxes are higher in areas with higher erosion rates, suggesting that tectonic setting is one factor that controls Re release and therefore OCpetro oxidation.more » « less
