No abstract

Fifty years of hyporheic zone research have shown the important role played by the hyporheic zone as an interface between groundwater and surface waters. However, it is only in the last two decades that what began as an empirical science has become a mechanistic science devoted to modeling studies of the complex fluid dynamical and biogeochemical mechanisms occurring in the hyporheic zone. These efforts have led to the picture of surface-subsurface water interactions as regulators of the form and function of fluvial ecosystems. Rather than being isolated systems, surface water bodies continuously interact with the subsurface. Exploration of hyporheic zone processes has led to a new appreciation of their wide reaching consequences for water quality and stream ecology. Modern research aims toward a unified approach, in which processes occurring in the hyporheic zone are key elements for the appreciation, management, and restoration of the whole river environment. In this unifying context, this review summarizes results from modeling studies and field observations about flow and transport processes in the hyporheic zone and describes the theories proposed in hydrology and fluid dynamics developed to quantitatively model and predict the hyporheic transport of water, heat, and dissolved and suspended compounds from sediment grain scale up to the watershed scale. The implications of these processes for stream biogeochemistry and ecology are also discussed.

SUMMARY Biofilm cells are less susceptible to antimicrobials than their planktonic counterparts. While this phenomenon is multifactorial, the ability of the biofilm matrix to reduce antibiotic penetration into the biofilm is thought to be of limited importance, as previous studies suggest that antibiotics move fairly rapidly through biofilms. In this study, we monitored the transport of two clinically relevant antibiotics, tobramycin and ciprofloxacin, into non-mucoid P. aeruginosa biofilms. To our surprise, we showed that the positively charged antibiotic tobramycin is sequestered to the biofilm periphery, while the neutral antibiotic ciprofloxacin readily penetrated. We provide evidence that tobramycin in the biofilm periphery both stimulated a localized stress response and killed bacteria in these regions, but not in the underlying biofilm. Although it is unclear which matrix component binds tobramycin, its penetration was increased by the addition of cations in a dose-dependent manner, which led to increased biofilm death. These data suggest that ionic interactions of tobramycin with the biofilm matrix limit its penetration. We propose that tobramycin sequestration at the biofilm periphery is an important mechanism in protecting metabolically active cells that lie just below the zone of sequestration.

No abstract

1] Temporary storage of solutes in streams is often controlled by flow-induced uptake in hyporheic zones. This phenomenon accounts for the tails that are generally observed following the passage of a solute pulse, and such exchange is particularly important for the transport of reactive substances that can be subject to various biogeochemical processes in the subsurface. Advective pumping, induced by streamflow over an irregular permeable bed, leads to a distribution of pore water flow paths in the streambed and a corresponding distribution of subsurface solute residence times. This paper describes a modeling framework that couples longitudinal solute transport in streams with solute advection along a continuous distribution of hyporheic flow paths. Moment methods are used to calculate the shape of solute breakthrough curves in the stream based on various representations of hyporheic exchange, including both advective pumping and several idealized formulations. Basic hydrodynamic principles are used to derive the distribution of solute residence times due to pumping. The model provides an accurate representation of the breakthrough curves of tritium along a 30 km reach of Säva Brook in Uppland County in Sweden. Both hydrodynamic theory for pumping exchange and pore water samples obtained from the bed during the tracer experiment suggest that the residence time for solutes in the hyporheic zone is characterized by a log normal probability density function. Closed-form solutions of the central temporal moments of solute breakthrough curves in the stream reveal a significant similarity between this new model and existing models of hyporheic exchange, including the Transient Storage Model. The new model is advantageous because its fundamentally derived exchange parameters can be expressed as functions of basic hydrodynamic quantities, which allows the model results to be generalized to conditions beyond those directly observed during tracer experiments. The utility of this approach is demonstrated by using the pumping theory to relate the spatial variation of hyporheic exchange rate along Säva Brook with the local Froude number, hydraulic conductivity and water depth.

[1] It is necessary to improve our understanding of the exchange of dissolved constituents between surface and subsurface waters in river systems in order to better evaluate the fate of water-borne contaminants and nutrients and their effects on water quality and aquatic ecosystems. Here we present a model that can predict hyporheic exchange at the bed-form-to-reach scale using readily measurable system characteristics. The objective of this effort was to compare subsurface flow induced at scales ranging from very small scale bed forms up to much larger planform geomorphic features such as meanders. In order to compare exchange consistently over this range of scales, we employed a spectral scaling approach as the basis for a generalized analysis of topography-induced stream-subsurface exchange. The spectral model involves a first-order approximation for local flow-boundary interactions but is fully three-dimensional and includes the lateral hyporheic zone in addition to the flow directly beneath the streambed. The primary model input parameters are stream velocity and slope, sediment permeability and porosity, and detailed measurements of the stream channel topography. The primary outputs are the distribution of water flux across the stream channel boundary, the resulting pore water flow paths, and the subsurface residence time distribution. We tested the bed-formexchange component of the model using a highly detailed two-dimensional data set for exchange with ripples and dunes and then applied the model to a three-dimensional meandering stream in a laboratory flume. Having spatially explicit information allowed us to evaluate the contributions of both gravitational and current-driven hyporheic flow through various classes of stream channel features including ripples, dunes, bars, and meanders. The model simulations indicate that all scales of topography between ripples and meanders have a significant effect on pore water flow fields and residence time distributions. Furthermore, complex interactions across the spectrum of topographic features play an important role in controlling the net interfacial flux and spatial distribution of hyporheic exchange. For example, shallow exchange induced by current-driven interactions with small bed forms dominates the interfacial flux, but local pore water flows are modified significantly by larger-scale surface-groundwater interactions. As a result, simplified representations of the stream topography do not adequately characterize patterns and rates of hyporheic exchange.

[1] Hyporheic exchange is generally analyzed with the assumption of a homogeneous hyporheic zone. In reality, streambed sediments have a heterogeneous structure, and this natural heterogeneity produces spatially variable interfacial fluxes and complex hyporheic exchange patterns. To assess the basic effects of sediment structure on hyporheic exchange, we performed salt and dye injection experiments in a recirculating laboratory flume with two heterogeneous sediment beds characterized by negativeexponential correlated random hydraulic conductivity fields. Dye injections showed that the hyporheic flow structure was controlled by the spatial relationship of bed forms to high-and low-permeability regions of the streambed. As no existing model could represent these effects, we developed a new finite element model to calculate the pore water flow field resulting from the interaction of the bed form-induced boundary head distribution and the heterogeneous sediment structure. A numerical particle-tracking approach was then used to assess the resulting hyporheic exchange. The combined flowtransport model did an excellent job of predicting the complex hyporheic flow pathways in the heterogeneous bed and the net hyporheic exchange up to t % 30 hours. The heterogeneous hydraulic conductivity field caused both greater spatial variability in the water flux through the bed surface and a greater average interfacial flux than would have occurred with a homogeneous bed. The layered correlation structure of the streambed produced an effective anisotropy that favored longitudinal pore water flow and caused a relatively rapid decrease of the mean pore water velocity with depth. As a result, solute penetration into the bed was confined to a more shallow region than would have occurred with a homogeneous bed. The combination of faster near-surface transport and shallower solute penetration produced a shorter mean hyporheic residence time. On the basis of the combination of experimental results and model simulations we conclude that the structural heterogeneity of streambed sediments produces more spatially limited hyporheic exchange that occurs with greater spatial variability and at a higher overall rate.

Abstract. Fine sediment exchange between a stream and the surrounding subsurface influences downstream contaminant transport and stream ecology. Fundamental models for this exchange were developed on the basis of (1) the hydraulics of bed form-driven advective pore water flow and (2) subsurface colloid transport processes. First, a model was developed to predict the advective flow induced in a sand bed by stream flow over bedforms. The resulting "pumping" exchange rate was calculated based on the streamflow conditions, bed form geometry, and bed depth. The pumping exchange of suspended sediment was then calculated by superimposing advective transport and particle settling in the bed and including the effect of physicochemical filtration by bed sediment. The filtration coefficient approach was used to predict the reduction in the concentration of transported particles. Both settling and filtration cause colloids to be trapped in stream beds, producing a higher net exchange rate relative to conservative solutes. When transported particles are completely trapped in a single pass through the bed, the exchange calculation is simplified because only the particle flux to the bed must be considered. In this case, the net exchange rate may be adequately represented by an effective piston velocity (flux/concentration) or loss rate to the bed in the advectiondispersion equation for the stream. Solute and colloid exchanges are predicted by the models without the use of fitting coefficients; only measurable hydraulic and particle parameters were used as model inputs. Simulations are presented which show the effect of stream parameters, settling, and filtration on net particle exchange. This fundamental approach to modeling stream-subsurface exchange potentially has great utility for understanding and predicting the transport and fate of reactive substances in streams.