skip to main content


Title: Matrix Diffusion as a Mechanism Contributing to Fractal Stream Chemistry and Long‐Tailed Transit Time Distributions
Abstract

Solute transit or travel time distributions (TTDs) in catchments are relevant to both hydrochemical response and inference of hydrologic mechanisms. Long‐tailed TTDs and fractal scaling behavior of stream concentration power spectra (∼1/frequency, or 1/frequency to a power <2) are widely observed in catchment studies. In several catchments, a significant fraction of streamflow is derived from groundwater in shallow fractured bedrock, where matrix diffusion significantly influences solute transport. I present frequency and time domain theoretical analyses of solute transport to quantify the influence of matrix diffusion on fractal scaling and long‐tailed TTDs. The theoretical concentration power spectra exhibit fractal scaling, and the corresponding TTDs resemble a gamma distribution. The tails of the TTDs are influenced by accessible matrix width, exhibiting a sustained power‐law (rather than exponential) decline for large matrix widths. Application to an experimental catchment shows that theoretical spectra match previously reported power spectral estimates derived from concentration measurements.

 
more » « less
Award ID(s):
2012353 2012264
NSF-PAR ID:
10447319
Author(s) / Creator(s):
 
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Geophysical Research Letters
Volume:
48
Issue:
18
ISSN:
0094-8276
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Spatially integrated transport models have been applied widely to model hydrologic transport. However, we lack simple and process‐based theoretical tools to predict the transport closures—transit time distributions (TTDs) and StorAge Selection (SAS) functions. This limits our ability to infer characteristics of hydrologic systems from tracer observations and to make first‐order estimates of SAS functions in catchments where no tracer data is available. Here we present a theoretical framework linking TTDs and SAS functions to hydraulic groundwater theory at the hillslope scale. For hillslopes where the saturated hydraulic conductivity declines exponentially with depth, analytical solutions for the closures are derived that can be used as hypotheses to test against data. In the simplest form, the hillslope SAS function resembles a uniform or exponential distribution (corresponding to flow pathways in the saturated zone) offset from zero by the storage in the unsaturated zone that does not contribute to discharge. The framework is validated against nine idealized virtual hillslopes constructed using a 2‐D Richards equation‐based model, and against data from tracer experiments in two artificial hillslopes. Modeled internal age, life expectancy, and transit time structures reproduce theoretical predictions. The experimental data also support the theory, though further work is needed to account for the effects of time‐variability. The shape and tailing of TTDs and their power spectra are discussed. The theoretical framework yields several dimensionless numbers that can be used to classify hillslope scale flow and transport dynamics and suggests distinct water age structures for high or low Hillslope number.

     
    more » « less
  2. 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 NO3flux. For the two low‐gradient tundra catchments, we observed consistent DOC enrichment (transport‐limitation) and NO3dilution (source‐limitation) during flow events. Conversely, we found consistent NO3enrichment 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.

     
    more » « less
  3. Abstract

    Processes that drive variability in catchment solute sourcing, transformation, and transport can be investigated using concentration–discharge (C–Q) relationships. These relationships reflect catchment and in‐stream processes operating across nested temporal scales, incorporating both short and long‐term patterns. Scientists can therefore leverage catchment‐scale C–Q datasets to identify and distinguish among the underlying meteorological, biological, and geological processes that drive solute export patterns from catchments and influence the shape of their respective C–Q relationships. We have synthesized current knowledge regarding the influence of biological, geological, and meteorological processes on C–Q patterns for various solute types across diel to decadal time scales. We identify cross‐scale linkages and tools researchers can use to explore these interactions across time scales. Finally, we identify knowledge gaps in our understanding of C–Q temporal dynamics as reflections of catchment and in‐stream processes. We also lay the foundation for developing an integrated approach to investigate cross‐scale linkages in the temporal dynamics of C–Q relationships, reflecting catchment biogeochemical processes and the effects of environmental change on water quality.

    This article is categorized under:

    Science of Water > Hydrological Processes

    Science of Water > Water Quality

    Science of Water > Water and Environmental Change

     
    more » « less
  4. Temporal patterns in stream chemistry provide integrated signals describing the hydrological and ecological state of whole catchments. However, stream chemistry integrates multi-scale signals of processes occurring in both the catchment and stream. Deconvoluting these signals could identify mechanisms of solute transport and transformation and provide a basis for monitoring ecosystem change. We applied trend analysis, wavelet decomposition, multivariate autoregressive state-space modeling, and analysis of concentration–discharge relationships to assess temporal patterns in high-frequency (15 min) stream chemistry from permafrost-influenced boreal catchments in Interior Alaska at diel, storm, and seasonal time scales. We compared catchments that varied in spatial extent of permafrost to identify characteristic biogeochemical signals. Catchments with higher spatial extents of permafrost were characterized by increasing nitrate concentration through the thaw season, an abrupt increase in nitrate and fluorescent dissolved organic matter (fDOM) and declining conductivity in late summer, and flushing of nitrate and fDOM during summer rainstorms. In contrast, these patterns were absent, of lower magnitude, or reversed in catchments with lower permafrost extent. Solute dynamics revealed a positive influence of permafrost on fDOM export and the role of shallow, seasonally dynamic flowpaths in delivering solutes from high-permafrost catchments to streams. Lower spatial extent of permafrost resulted in static delivery of nitrate and limited transport of fDOM to streams. Shifts in concentration–discharge relationships and seasonal trends in stream chemistry toward less temporally dynamic patterns might therefore indicate reorganized catchment hydrology and biogeochemistry due to permafrost thaw. 
    more » « less
  5. Abstract

    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