Anisotropy and depth‐related heterogeneity of hydraulic conductivity in a bog peat. II: modelling the effects on groundwater flow

Anisotropy and heterogeneity of hydraulic conductivity (K) are seldom considered in models of mire hydrology. We investigated the effect of anisotropy and heterogeneity on groundwater flow in bog peat using a steady‐state groundwater model. In five model simulations, four sets of K data were used. The first set comprised measured K values from an anisotropic and heterogeneous bog peat. These data were aggregated to produce the following simplified data sets: an isotropic and heterogeneous distribution of K; an isotropic and homogeneous distribution; and an anisotropic and homogeneous distribution. We demonstrate that, where anisotropy and heterogeneity exist, groundwater flow in bog peat is complex. Fine‐scale variations in K have the potential to influence patterns and rates of groundwater flow. However, for our data at least, it is heterogeneity and not anisotropy that has the greater influence on producing complex patterns of groundwater flow. We also demonstrate that patterns and rates of groundwater flow are simplified and reduced when measured K values are aggregated to create a more uniform distribution of K. For example, when measured K values are aggregated to produce isotropy and homogeneity, the rate of modelled seepage is reduced by 28%. We also show that when measured K values are used, the presence of a drainage ditch can increase seepage through a modelled cross‐section. Our work has implications for the accurate interpretation of hydraulic head data obtained from peat soils, and also the understanding of the effect of drainage ditches on patterns and rates of groundwater flow. Copyright © 2002 John Wiley & Sons, Ltd.     

Does streambed heterogeneity matter for hyporheic residence time distribution in sand-bedded streams?

Stream water residence times within streambed sediments are key values to quantify hyporheic processes including sediment thermal regime, solute transient storage, dilution rates and biogeochemical transformations, such as those controlling degassing nitrous oxide. Heterogeneity of the streambed sediment hydraulic properties has been shown to be potentially an important factor to characterize hyporheic processes. Here, we quantify the importance of streambed heterogeneity on residence times of dune-like bedform induced hyporheic fluxes at the bedform and reach scales. We show that heterogeneity has a net effect of compression of the hyporheic zone (HZ) toward the streambed, changing HZ volume from the homogenous case and thus inducing remarkable differences in the flow field with respect to the homogeneous case. We unravel the physical conditions for which the commonly used homogenous field assumption is applicable for quantifying hyporheic processes thus explaining why predictive measures based on a characteristic residence time, like the Damköhler number, are robust in heterogeneous sand bedded streams.     

Spatial variability in streambed hydraulic conductivity of contrasting stream morphologies: channel bend and straight channel

Streambed hydraulic conductivity is one of the main factors controlling variability in surface water‐groundwater interactions, but only few studies aim at quantifying its spatial and temporal variability in different stream morphologies. Streambed horizontal hydraulic conductivities (K_h) were therefore determined from in‐stream slug tests, vertical hydraulic conductivities (K_v) were calculated with in‐stream permeameter tests and hydraulic heads were measured to obtain vertical head gradients at eight transects, each comprising five test locations, in a groundwater‐dominated stream. Seasonal small‐scale measurements were taken in December 2011 and August 2012, both in a straight stream channel with homogeneous elevation and downstream of a channel meander with heterogeneous elevation. All streambed attributes showed large spatial variability. K_h values were the highest at the depositional inner bend of the stream, whereas high K_v values were observed at the erosional outer bend and near the middle of the channel. Calculated K_v values were related to the thickness of the organic streambed sediment layer and also showed higher temporal variability than K_h because of sedimentation and scouring processes affecting the upper layers of the streambed. Test locations at the channel bend showed a more heterogeneous distribution of streambed properties than test locations in the straight channel, whereas within the channel bend, higher spatial variability in streambed attributes was observed across the stream than along the stream channel. Copyright © 2014 John Wiley & Sons, Ltd.     

Spatial and depth variability of streambed vertical hydraulic conductivity under the regional flow regimes

This study investigated the influence of the regional flow on the streambed vertical hydraulic conductivity (K_v) within the hyporheic zone in three stream reaches of the Weihe River in July 2016. The streambed K_v with two connected depths was investigated at each test reach. Based on the sediment characteristics, the three test reaches could be divided into three categories: a sandy streambed without continuous silt and clay layer, a sandy streambed with continuous silt and clay layer, and a silt–clay streambed. The results demonstrate that the streambed K_v mainly decreases with the depth at the sandy streambed (without continuous silt and clay layer) and increases with the depth at the other two test reaches. At the sandy streambed (with continuous silt and clay layer) where streambed K_v mainly decreases with the depth, the regional upward flux can suspend fine particles and enhance the pore spacing, resulting in the elevated K_v in the upper sediment layers. At another sandy streambed, the continuous silt and clay layer is the main factor that influences the vertical distribution of fine particles and streambed K_v. An increase in streambed K_v with the depth at the silt/clay streambed is attributed to the regional downward movement of water within the sediments that may lead to more fine particles deposited in the pores in the upper sediment layers. The streambed K_v is very close to the bank in the sandy streambed without continuous silt and clay layer and the channel centre in the other two test reaches. Differences in grain size distribution of the sediments at each test reach exercise a strong controlling influence on the streambed K_v. This study promotes the understanding of dynamics influencing the interactions between groundwater and surface water and provides guidelines to scientific water resources management for rivers.     

Heterogeneous characteristics of streambed saturated hydraulic conductivity of the Touchet River,south eastern Washington,USA

Traditionally a streambed is treated as a layer of uniform thickness and low saturated hydraulic conductivity (K) in surface‐ and ground‐water studies. Recent findings have shown a high level of spatial heterogeneity within a streambed and such heterogeneity directly affects surface‐ and ground‐water exchange and can have ecological implications for biogeochemical transformations, nutrient cycling, organic matter decomposition, and reproduction of gravel spawning fish. In this study a detailed field investigation of K was conducted in two selected sites in Touchet River, a typical salmon spawning stream in arid south eastern Washington, USA. In‐stream slug tests were conducted to determine K following the Bouwer and Rice method. For the upper and lower sites, each 50 m long and 9 m wide and roughly 20 m apart, a sampling grid of 5 m longitudinally and 3 m transversely was used. The slug tests were performed for each horizontal coordinate at 0·3–0·45, 0·6–0·75, 0·9–1·05 and 1·2–1·35 m depth intervals unless a shallower impenetrable obstruction was encountered. Additionally, water levels were measured to obtain vertical hydraulic gradient (VHG) between each two adjacent depth intervals. Results indicated that K ranged over three orders of magnitude at both the upper and lower sites and differed between the two sites. At the upper site, K did not differ significantly among different depth intervals based on nonparametric statistical tests for mean, median, and empirical cumulative distribution, but the spatial pattern of K varied among different depth intervals. At the lower site, K for the 0·3–0·45 m depth interval differed statistically from those at other depth intervals, and no similar spatial pattern was found among different depth intervals. Zones of upward and downward water flow based on VHG also varied among different depth intervals, reflecting the complexities of the water flow regime. Detailed characterization of the streambed as attempted in this study should be helpful in providing information on spatial variations of streambed hydraulic properties as well as surface‐ and ground‐water interaction. Copyright © 2009 John Wiley & Sons, Ltd.     

Feasibility of grain-size analysis methods for determination of vertical hydraulic conductivity of streambeds

Accurate estimation of streambed vertical hydraulic conductivity (K_v) is of great importance in the analysis of water quantity exchange and solute transfer between a stream and its sediments. The paper analyzed the inaccuracy of hydraulic conductivity values of sediments derived from grain-size distribution (K_g), which were determined from eight empirical grain-size methods to represent streambed K_v. In this study, the values of K_v for a streambed were derived using falling-head standpipe permeameter tests conducted at eight study sites in the Elkhorn River, Nebraska, and the tested streambed columns were then collected for grain-size analysis by sieving. These empirical methods were used to calculate the K_g values of the streambed from grain-size distribution data of sediments. Unlike many other studies, this study verifies K_g from grain-size distribution with K_v from permeameter tests on the basis of the same samples of streambed sediments. The K_g values derived from the eight empirical methods were larger than the K_v from permeameter tests; there are five methods that give K_g values of about 3–6 times larger than these K_v. The K_g values from the Kozeny formula followed by the Hazen formula give the largest overestimation error if they are used to represent the K_v of the streambed. The USBR and Shepherd formulas generated K_g values close to K_v, but these K_g values are still larger in general than the K_v values. Moreover, the new values of coefficient C for the empirical formulas were revised so that they can be used to calculate the approximate K_v of a streambed. Among the eight methods, the ratios of the original C values to the average new C range from 1.3 to 5.9. It can be hypothesized that smaller C values must be used in the estimation of K_v for general soil samples if these empirical formulas are used to calculate K_v.     

Transient streambed hydraulic conductivity in channel and bar environments,Loup River,Nebraska

Streambed hydraulic conductivity (K) and vertical K (K _v) are key controls on groundwater and surface water exchange and biogeochemical fluxes through the hyporheic zone, but drivers of transient hydraulic properties in different fluvial environments are poorly understood. This study combines hydrogeology, geophysics, and sedimentology to reveal mechanisms of K and K _v transience in the upper 0.5 m of a sandy streambed during low discharge. Hydraulic tests (44 slug tests, 130 falling-head permeameter tests) and 130 grain-size analyses were repeated three times over 8 weeks on a 1,200 m² grid spanning: (a) a channel with continuously flowing water and mobile bed load; (b) an adjacent mid-channel bar that was stationary and infrequently submerged. Aerial photographs and ground-penetrating radar show scour and complete reworking of fluvial sediments in the channel. Bar sediments below the water table remained immobile, but infrequent flows of moderate discharge reworked the uppermost few centimetres of the bar top. Despite differences in sediment mobility and stream flow characteristics across environments, K and K _v exhibited order-of-magnitude differences in spatial heterogeneity and temporal variability in both the channel and bar. Mean K and K _v values in the channel were comparatively stable over time. In the immobile bar, mean K declined 20% and K _v declined 26% after increased discharge temporarily inundated the bar. Grain-size distributions were steady across both environments over time, but repeat geophysical surveys of the bar show a decrease in electrical conductivity, likely from porosity reduction. These findings suggest that sediment dynamics and stream flow characteristics in different streambed environments are important drivers of K transience during low discharge conditions. Specifically, pore clogging can be an important mechanism of transience over short durations (weeks to months) in immobile sediments subject to infrequent flows and minor reworking.     

A heat pulse technique for the determination of small‐scale flow directions and flow velocities in the streambed of sand‐bed streams

The hyporheic zone (HZ) has the capability to eliminate and attenuate nutrients and contaminants in riverine systems. Biogeochemical reactions and the potential elimination of contaminants are strongly controlled by the flow paths and dynamics in the HZ. Nevertheless, an easily applicable method for the field determination of flow patterns in the HZ is still lacking. Therefore, a heat pulse technique, which traces the movement of a short heat pulse in the upper part of the HZ and other sand beds, was developed. Five rods are vertically driven into the sediment of the streambed; one rod with a heater as point source located in about 10‐cm sediment depth and four rods with four temperature sensors in 3 cm distance, arranged concentrically with 7 cm diameter around the heating rod. Subsequently, a heat pulse is applied and the resulting breakthrough curves are indicative of flow velocities and flow directions in the streambed. A rough data analysis procedure is also suggested. In addition, laboratory experiments were performed to test the heat pulse technique. These experiments were validated based on coupled numerical modelling of flow and heat transport. First field tests of the method prove that the method is easily applicable under field conditions. These first field tests showed highly complex flow patterns with flow velocities from 1·8 to 4·9 cm min^?1 and flow directions from parallel to surface flow to opposite to surface flow. This suggests the need for a robust method to quantify hyporheic flow patterns in situ. Copyright © 2011 John Wiley & Sons, Ltd.     

Effect of anisotropy on solute transport in degraded fen peat soils

Peat soils are heterogeneous, anisotropic porous media. Compared to mineral soils, there is still limited understanding of physical and solute transport properties of fen peat soils. In this study, we aimed to explore the effect of soil anisotropy on solute transport in degraded fen peat. Undisturbed soil cores, taken in vertical and horizontal direction, were collected from one drained and one restored fen peatland both in a comparable state of soil degradation. Saturated hydraulic conductivity (K _s) and chemical properties of peat were determined for all soil cores. Miscible displacement experiments were conducted under saturated steady state conditions using potassium bromide as a conservative tracer. The results showed that (1) the K _s in vertical direction (K _sv) was significantly higher than that in horizontal direction (K_sh), indicating that K _s of degraded fen peat behaves anisotropically; (2) pronounced preferential flow occurred in vertical direction with a higher immobile water fraction and a higher pore water velocity; (3) the 5% arrival time (a proxy for the strength of preferential flow) was affected by soil anisotropy as well as study site. A strong correlation was found between 5% arrival time and dispersivity, K _s and mobile water fraction; (4) phosphate release was observed from drained peat only. The impact of soil heterogeneity on phosphate leaching was more pronounced than soil anisotropy. The soil core with the strongest preferential flow released the highest amount of phosphate. We conclude that soil anisotropy is crucial in peatland hydrology but additional research is required to fully understand anisotropy effects on solute transport.     

Evaluation of streambed hydraulic conductivity heterogeneity in an urban watershed

It has long been understood that streambed hydraulic conductivity plays an important role in surface-subsurface solute exchange. Using a portable falling head permeameter in situ, we estimated the horizontal hydraulic conductivity, K, of the near-surface streambed sediments at a total of 85 locations encompassing two depth intervals: 7.5–10 and 10–12.5 cm. The measurements were conducted in an 80 m reach of Indian Creek, a small urban stream in Philadelphia, PA, USA. We found that the ln K data within each sediment layer were Gaussian, but the combined data set was not. The results indicated that while the mean hydraulic conductivity decreased with depth, horizontal heterogeneity (e.g. the variance) increased with depth. This strong contrast between layers suggests that they should be treated as separated entities in modeling studies. Variogram analyses across the stream suggested symmetry with respect to the thalweg in the upper layer and fractality in the lower layer. The variograms along the streams suggested that the K data are random.     

Depth‐dependent hydraulic conductivity distribution patterns of a streambed

Characterization of streambed hydraulic conductivity from the channel surface to a great depth below the channel surface can provide needed information for the determination of stream‐aquifer hydrologic connectedness, and it is also important to river restoration. However, knowledge on the streambed hydraulic conductivity for sediments 1 m below the channel surface is scarce. This study describes a method that was used to determine the distribution patterns of streambed hydraulic conductivity for sediments from channel surface to a depth of 15 m below. The method includes Geoprobe's direct‐push techniques and Permeameter tests. Direct‐push techniques were used to generate the electrical conductivity (EC) logs and to collect sequences of continuous sediment cores from river channels, as well as from the alluvial aquifer connected to the river. Permeameter tests on these sediment cores give the profiles of vertical hydraulic conductivity (K_v) of the channel sediments and the aquifer materials. This method was applied to produce K_v profiles for a streambed and an alluvial aquifer in the Platte River Valley of Nebraska, USA. Comparison and statistical analysis of the K_v profiles from the river channel and from the proximate alluvial aquifer indicates a special pattern of K_v in the channel sediments. This depth‐dependent pattern of K_v distribution for the channel sediments is considered to be produced by hyporheic processes. This K_v‐distribution pattern implied that the effect of hyporheic processes on streambed hydraulic conductivity can reach the sediments about 9 m below the channel surface. Copyright © 2010 John Wiley & Sons, Ltd.     

Two-Dimensional Probabilistic Infiltration Analysis in a Hillslope Using First-Order Moment Approach

A first-order moment analysis method is introduced to evaluate the pore-water pressure variability within a hillslope due to spatial variability in saturated hydraulic conductivity (K_s) during rainfall. The influences of the variance of the natural logarithm of K_s(ln K_s), spatial structure anisotropy of ln K_s, and normalized vertical infiltration flux (q) on the evaluations of the pore-water pressure uncertainty are investigated. Results indicate different responses of pressure head variability in the unsaturated region and the saturated region. In the unsaturated region, a larger variance of ln K_s, a higher spatial structure anisotropy, and a smaller q lead to a larger variability in pressure head, while in the saturated region, the variability in pressure head increases with the increase of variance of ln K_s, the decrease of spatial structure anisotropy, or the increase of q. These variables have great impacts on the range of fluctuation of the phreatic surface within the hillslope. The influences of these three variables on the variance of pressure head within the saturated region are greater than those within the unsaturated region, and the variance of ln K_s has the greatest impact. These results yield useful insight into the effects of heterogeneity on pressure head and uncertainty associated with predicted flow field.     

Effects of variations of river stage and hydraulic conductivity on temporal scaling of groundwater levels: numerical simulations

It was found in previous studies that groundwater levels may fluctuate as a temporal fractal. In this study numerical simulations of groundwater level fluctuations in an unconfined aquifer near a river were conducted to investigate the effects of aquifer heterogeneity and river stage variations on the fractal behavior of the water levels, h(t). Groundwater recharge was taken to be a white-noise process. The aquifer heterogeneity was simulated with a second-order stationary field of hydraulic conductivity (K) with an exponential variogram model. The results showed that groundwater levels fluctuate as a temporal fractal in both homogeneous and heterogeneous aquifers as long as K is less than 10 m/d. Most aquifers may indeed act as a fractal filter which takes a random non-fractal recharge inputs and produces a fractal responses of groundwater level fluctuations. A crossover in temporal scaling of h(t) may appear in more permeable aquifers. Fluctuations of the groundwater level in a homogeneous aquifer are dominated by the recharge process when the river stage is constant or by the river stage variations when the river stage varies in highly permeable aquifers. Heterogeneity plays an important role in the temporal scaling of h(t) in more permeable aquifers: the stronger the heterogeneity, the stronger the temporal scaling of h(t).     

The usefulness of outcrop analogue air permeameter measurements for analyzing aquifer heterogeneity: quantifying outcrop hydraulic conductivity and its spatial variability

Saturated hydraulic conductivity (K) is one of the most important parameters determining groundwater flow and contaminant transport in both unsaturated and saturated porous media. Although several well‐established laboratory methods exist for determining K, in situ measurements of this parameter remain very complex and scale dependent. Often, the limited accessibility of subsurface sediments for sampling means an additional impediment to our ability to quantify subsurface K heterogeneity. One potential solution is the use of outcrops as analogues for subsurface sediments. This paper investigates the use of air permeameter measurements on outcrops of unconsolidated sediments to quantify K and its spatial heterogeneity on a broad range of sediment types. The Neogene aquifer in northern Belgium is used as a case study for this purpose. To characterize the variability in K, 511 small‐scale air permeability measurements were performed on outcrop sediments representative over five of the aquifer's lithostratigraphic units. From these measurements, outcrop‐scale equivalent K tensors were calculated using numerical upscaling techniques. Validation of the air permeameter‐based K values by comparison with laboratory constant head K measurements reveals a correlation of 0.93. Overall, the results indicate that hand‐held air permeameters are very efficient and accurate tools to characterize saturated K, as well as its small‐scale variability and anisotropy on a broad range of unconsolidated sediments. The studied outcrops further provided a qualitative understanding of aquifer hydrostratigraphy and quantitative estimates about K variability at the centimetre‐scale to metre‐scale. Copyright © 2013 John Wiley & Sons, Ltd.     

Bias of Apparent Tracer Ages in Heterogeneous Environments

The interpretation of apparent ages often assumes that a water sample is composed of a single age. In heterogeneous aquifers, apparent ages estimated with environmental tracer methods do not reflect mean water ages because of the mixing of waters from many flow paths with different ages. This is due to nonlinear variations in atmospheric concentrations of the tracer with time resulting in biases of mixed concentrations used to determine apparent ages. The bias of these methods is rarely reported and has not been systematically evaluated in heterogeneous settings. We simulate residence time distributions (RTDs) and environmental tracers CFCs, SF₆, ⁸⁵Kr, and ³⁹Ar in synthetic heterogeneous confined aquifers and compare apparent ages to mean ages. Heterogeneity was simulated as both K‐field variance (σ²) and structure. We demonstrate that an increase in heterogeneity (increase in σ² or structure) results in an increase in the width of the RTD. In low heterogeneity cases, widths were generally on the order of 10 years and biases generally less than 10%. In high heterogeneity cases, widths can reach 100 s of years and biases can reach up to 100%. In cases where the temporal variations of atmospheric concentration of individual tracers vary, different patterns of bias are observed for the same mean age. We show that CFC‐12 and CFC‐113 ages may be used to correct for the mean age if analytical errors are small.     

Field Study of Subsurface Heterogeneity with Steady‐State Hydraulic Tomography

Remediation of subsurface contamination requires an understanding of the contaminant (history, source location, plume extent and concentration, etc.), and, knowledge of the spatial distribution of hydraulic conductivity (K) that governs groundwater flow and solute transport. Many methods exist for characterizing K heterogeneity, but most if not all methods require the collection of a large number of small‐scale data and its interpolation. In this study, we conduct a hydraulic tomography survey at a highly heterogeneous glaciofluvial deposit at the North Campus Research Site (NCRS) located at the University of Waterloo, Waterloo, Ontario, Canada to sequentially interpret four pumping tests using the steady‐state form of the Sequential Successive Linear Estimator (SSLE) ( Yeh and Liu 2000 ). The resulting three‐dimensional (3D) K distribution (or K‐tomogram) is compared against: ( 1 ) K distributions obtained through the inverse modeling of individual pumping tests using SSLE, and ( 2 ) effective hydraulic conductivity (K_eff) estimates obtained by automatically calibrating a groundwater flow model while treating the medium to be homogeneous. Such a K_eff is often used for designing remediation operations, and thus is used as the basis for comparison with the K‐tomogram. Our results clearly show that hydraulic tomography is superior to the inversions of single pumping tests or K_eff estimates. This is particularly significant for contaminated sites where an accurate representation of the flow field is critical for simulating contaminant transport and injection of chemical and biological agents used for active remediation of contaminant source zones and plumes.     

Equivalent diffusion coefficient of clay-rich geological formations: comparison between numerical and analytical estimates

Laboratory experiments in rock samples collected from clay-rich formations indicate that the effective molecular diffusion coefficient (D) is a heterogeneous and anisotropic property. Since laboratory measurements of D are representative of a very small volume, upscaling is necessary in order to incorporate these data in large-scale numerical models of diffusive transport. In this work we address the problem of the estimating the equivalent diffusion coefficient (D _eq ), in terms of total diffusive flux, in a three-dimensional domain characterized by a heterogeneous and anisotropic spatial distribution of D. D _eq was estimated from the results of steady-state diffusive transport simulations through several realizations of the D field. The ensemble averages of D _eq from fields with different degrees of heterogeneity and anisotropy were then compared with estimates from analytical upscaling expressions based on stochastic as well as power-averaging approaches. These expressions are largely based on similar expressions developed for calculating the effective hydraulic conductivity in heterogeneous and anisotropic domains. Comparisons showed that stochastic expressions provide accurate estimates of D _eq only for fields characterized by low heterogeneity. Within the range of heterogeneity and anisotropy considered, our results showed that a power-averaging expression is very accurate in predicting D _eq especially when the parameter p _i is estimated through fitting of the numerical results. Nonetheless, the relationship between this parameter and the anisotropy ratio is linear.     

Numerical modelling of stream–aquifer interaction: Quantifying the impact of transient streambed permeability and aquifer heterogeneity

Stream–aquifer interaction plays a vital role in the water cycle, and a proper study of this interaction is needed for understanding groundwater recharge, contaminants migration, and for managing surface water and groundwater resources. A model‐based investigation of a field experiment in a riparian zone of the Schwarzbach river, a tributary of the Rhine River in Germany, was conducted to understand stream–aquifer interaction under alternative gaining and losing streamflow conditions. An equivalent streambed permeability, estimated by inverting aquifer responses to flood waves, shows that streambed permeability increased during infiltration of stream water to aquifer and decreased during exfiltration. Aquifer permeability realizations generated by multiple‐point geostatistics exhibit a high degree of heterogeneity and anisotropy. A coupled surface water groundwater flow model was developed incorporating the time‐varying streambed permeability and heterogeneous aquifer permeability realizations. The model was able to reproduce varying pressure heads at two observation wells near the stream over a period of 55 days. A Monte Carlo analysis was also carried out to simulate groundwater flow, its age distribution, and the release of a hypothetical wastewater plume into the aquifer from the stream. Results of this uncertainty analysis suggest (a) stream–aquifer exchange flux during the infiltration periods was constrained by aquifer permeability; (b) during exfiltration, this flux was constrained by the reduced streambed permeability; (c) the effect of temporally variable streambed permeability and aquifer heterogeneity were found important to improve the accurate capture of the uncertainty; and (d) probabilistic infiltration paths in the aquifer reveal that such pathways and the associated prediction of the extent of the contaminant plume are highly dependent on aquifer heterogeneity.     

High frequency stream bed mobility of a low‐gradient agricultural stream with implications on the hyporheic zone

Little Kickapoo Creek (LKC), a low‐gradient stream, mobilizes its streambed–fundamentally altering its near‐surface hyporheic zone–more frequently than do higher‐gradient mountain and karst streams. LKC streambed mobility was assessed through streambed surveys, sediment sampling, and theoretical calculations comparing basal shear stress (τ_b) with critical shear stress (τ_c). Baseflow τ_b is capable of entraining a d₅₀ particle; bankfull flow could entrain a 51·2 mm particle. No particle that large occurs in the top 30 cm of the substrate, suggesting that the top 30 cm of the substrate is mobilized and redistributed during bankfull events. Bankfull events occur on average every 7·6 months; flows capable of entraining d₅₀ and d₈₅ particles occur on average every 0·85 and 2·1 months, respectively. Streambed surveys verify streambed mobility at conditions below bankfull. While higher gradient streams have higher potential energy than LKC, they achieve streambed‐mobilization thresholds less frequently. Heterogeneous sediment redistribution creates an environment where substrate hydraulic conductivity (K) varies over four orders of magnitude. The frequency and magnitude of the substrate entrainment has implications on hyporheic zone function in fluid, solute and thermal transport models, interpretations of hyporheic zone stability, and understanding of LKC's aquatic ecosystem. Copyright © 2008 John Wiley & Sons, Ltd.