Local‐scale variability of seepage and hydraulic conductivity in a shallow gravel‐bed river

Seepage rate and direction measured with a seepage metre modified for use in flowing water were greatly variable along a 300‐m reach of a shallow, gravel‐bed river and depended primarily on the local‐scale bed topography. The median value of seepage measured at 24 locations was 24 cm/day, but seepage measured at specific sites ranged from ?340 to +237 cm/day. Seepage also varied substantially over periods of hours to days and occasionally reversed direction in response to evolution of the sediment bed. Vertical hydraulic conductivity was related to seepage direction and was larger during upward seepage than during downward seepage; with differences ranging from 4 to 40% in areas of active sediment transport to more than an order of magnitude in areas where current was too slow to mobilize bed sediment. Seepage was poorly related to hydraulic gradient measured over vertical distances of 0·3 m and appeared to be opposite the hydraulic gradient at 18% of the locations where both parameters were measured. Results demonstrate the scale dependence of these measurements in coarse‐grained hyporheic settings and indicate that hydraulic gradients should be determined over a much shorter vertical increment if used to indicate exchange across the sediment–water interface. Published in 2009 by John Wiley & Sons, Ltd.     

In situ quantification of spatial and temporal variability of hyporheic exchange in static and mobile gravel‐bed rivers

Seepage meters modified for use in flowing water were used to directly measure rates of exchange between surface and subsurface water in a gravel‐ and cobble bed river in western Pennsylvania, USA (Allegheny River, Q_mean = 190 m³/s) and a sand‐ and gravel‐bed river in Colorado, USA (South Platte River, Q_mean = 9·7 m³/s). Study reaches at the Allegheny River were located downstream from a dam. The bed was stable with moss, algae, and river grass present in many locations. Median seepage was + 0·28 m/d and seepage was highly variable among measurement locations. Upward and downward seepage greatly exceeded the median seepage rate, ranging from + 2·26 (upward) to ? 3·76 (downward) m/d. At the South Platte River site, substantial local‐scale bed topography as well as mobile bedforms resulted in spatial and temporal variability in seepage greatly in exceedence of the median groundwater discharge rate of 0·24 m/d. Both upward and downward seepage were recorded along every transect across the river with rates ranging from + 2·37 to ? 3·40 m/d. Despite a stable bed, which commonly facilitates clogging by fine‐grained or organic sediments, seepage rates at the Allegheny River were not reduced relative to those at the South Platte River. Seepage rate and direction depended primarily on measurement position relative to local‐ and meso‐scale bed topography at both rivers. Hydraulic gradients were small at nearly all seepage‐measurement locations and commonly were not a good indicator of seepage rate or direction. Therefore, measuring hydraulic gradient and hydraulic conductivity at in‐stream piezometers may be misleading if used to determine seepage flux across the sediment‐water interface. Such a method assumes that flow between the well screen and sediment‐water interface is vertical, which appears to be a poor assumption in coarse‐grained hyporheic settings. Copyright © 2011 John Wiley & Sons, Ltd.     

Investigation of Seepage Meter Measurements in Steady Flow and Wave Conditions

Water exchange between surface water and groundwater can modulate or generate ecologically important fluxes of solutes across the sediment‐water interface. Seepage meters can directly measure fluid flux, but mechanical resistance and surface water dynamics may lead to inaccurate measurements. Tank experiments were conducted to determine effects of mechanical resistance on measurement efficiency and occurrence of directional asymmetry that could lead to erroneous net flux measurements. Seepage meter efficiency was high (average of 93%) and consistent for inflow and outflow under steady flow conditions. Wave effects on seepage meter measurements were investigated in a wave flume. Seepage meter net flux measurements averaged 0.08 cm/h—greater than the expected net‐zero flux, but significantly less than theoretical wave‐driven unidirectional discharge or recharge. Calculations of unidirectional flux from pressure measurements (Darcy flux) and theory matched well for a ratio of wave length to water depth less than 5, but not when this ratio was greater. Both were higher than seepage meter measurements of unidirectional flux made with one‐way valves. Discharge averaged 23% greater than recharge in both seepage meter measurements and Darcy calculations of unidirectional flux. Removal of the collection bag reduced this net discharge. The presence of a seepage meter reduced the amplitude of pressure signals at the bed and resulted in a nearly uniform pressure distribution beneath the seepage meter. These results show that seepage meters may provide accurate measurements of both discharge and recharge under steady flow conditions and illustrate the potential measurement errors associated with dynamic wave environments.     

Influence of Groundwater Hydraulic Gradient on Bank Storage Metrics

The hydraulic gradient between aquifers and rivers is one of the most variable properties in a river/aquifer system. Detailed process understanding of bank storage under hydraulic gradients is obtained from a two‐dimensional numerical model of a variably saturated aquifer slice perpendicular to a river. Exchange between the river and the aquifer occurs first at the interface with the unsaturated zone. The proportion of total water exchanged through the river bank compared to the river bed is a function of aquifer hydraulic conductivity, partial penetration, and hydraulic gradient. Total exchange may be estimated to within 50% using existing analytical solutions provided that unsaturated zone processes do not strongly influence exchange. Model‐calculated bank storage is at a maximum when no hydraulic gradient is present and increases as the hydraulic conductivity increases. However, in the presence of a hydraulic gradient, the largest exchange flux or distance of penetration does not necessarily correspond to the highest hydraulic conductivity, as high hydraulic conductivity increases the components of exchange both into and out of an aquifer. Flood wave characteristics do not influence ambient groundwater discharge, and so in large floods, hydraulic gradients must be high to reduce the volume of bank storage. Practical measurement of bank storage metrics is problematic due to the limitations of available measurement technologies and the nested processes of exchange that occur at the river‐aquifer interface. Proxies, such as time series concentration data in rivers and groundwater, require further development to be representative and quantitative.     

Estimating streambed parameters for a disconnected river

Evaluation of stream–aquifer interaction and water balance for a catchment often requires specific information on streambed parameters, such as streambed hydraulic conductivity, seepage flux across the streambed and so on. This paper describes a simple, inexpensive instrument that is used to measure these streambed parameters under the condition of a stream disconnected from groundwater. Our method includes a seepage cylinder for simulation of river water depth. The proposed method was applied to estimate the vertical hydraulic conductivity of a streambed and the changes in vertical seepage rate from stream to groundwater with varied stream water depth in the Manasi River of Xinjiang Uygur Autonomous Region, China. The vertical hydraulic conductivities of the streambed determined from 12 sites along the Manasi River vary from 1.01 to 29.m/day where the stream disconnects from the groundwater. The experimental results suggest that there are two kinds of relations between the vertical seepage rate and the simulated stream water depth. One is a linear relation between the two variables with low Reynolds numbers (less than 10); the other is a nonlinear relation (exponential relation) between the two variables with larger Reynolds numbers (greater than 10). This second relationship is quite different from the traditional model that usually calculates the vertical seepage rate from stream to groundwater under the condition of disconnection using a linear relation (Darcy's Law). Our results suggest that a linear relation can only be used for a limited range of river water depth. This method gives a convenient tool for rapidly estimating the streambed hydraulic conductivity and the changes in the vertical seepage rate across streambed with varied stream water depths for the case of a stream disconnected from groundwater. Copyright © 2013 John Wiley & Sons, Ltd.     

Effect of rapidly changing river stage on uranium flux through the hyporheic zone

Measurement of ground water/surface water interaction within the hyporheic zone is increasingly recognized as an important aspect of subsurface contaminant fate and transport. Understanding the interaction between ground water and surface water is critical in developing a complete conceptual model of contaminant transport through the hyporheic zone. At the Hanford Site near Richland, Washington, ground water contaminated with uranium discharges to the Columbia River through the hyporheic zone. Ground water flux varies according to changes in hydraulic gradient caused by fluctuating river stage, which changes in response to operation of dams on the Columbia River. Piezometers and continuous water quality monitoring probes were installed in the hyporheic zone to provide long-term, high-frequency measurement of hydraulic gradient and estimated uranium concentrations. Subsequently, the flux of water and uranium was calculated for each half-hour time period over a 15-month study period. In addition, measurement of water levels in the near-shore unconfined aquifer enhanced the understanding of the relationship between river stage, aquifer elevation, and uranium flux. Changing river stage resulted in fluctuating hydraulic gradient within the hyporheic zone. Further, influx of river water caused lower uranium concentrations as a result of dilution. The methods employed in this study provide a better understanding of the interaction between surface and ground water in a situation with a dynamically varying vertical hydraulic gradient and illustrate how the combination of relatively standard methods can be used to derive an accurate estimation of water and contaminant flux through the hyporheic zone.     

Variability and comparison of hyporheic water temperatures and seepage fluxes in a small Atlantic salmon stream

Ground water discharge is often a significant factor in the quality of fish spawning and rearing habitat and for highly biologically productive streams. In the present study, water temperatures (stream and hyporheic) and seepage fluxes were used to characterize shallow ground water discharge and recharge within thestreambed of Catamaran Brook, a small Atlantic salmon (Salmo salar) stream in central New Brunswick, Canada. Three study sites were instrumented using a total of 10 temperature sensors and 18 seepage meters. Highly variable mean seepage fluxes, ranging from 1.7 x 10(-4) to 2.5 cm3 m(-2) sec(-1), and mean hyporheic water temperatures, ranging from 10.5 degrees to 18.0 degrees C, at depths of 20 to 30 cm in the streambed were dependent on streambed location (left versus right stream bank and site location) and time during the summer sampling season. Temperature data were usefulfor determining if an area of the streambed was under discharge (positive flux), recharge (negative flux), or parallel flow (no flux) conditions and seepage meters were used to directly measure the quantity of water flux. Hyporheic water temperature measurements and specific conductance measurements of the seepage meter sample water, mean values ranging from 68.8 to 157.9 microS/cm, provided additional data for determining flux sources. Three stream banks were consistently under discharge conditions, while the other three stream banks showed reversal from discharge to recharge conditions over the sampling season. Results indicate that the majority of the water collected in the seepage meters was composed of surface water. The data obtained suggests that even though a positive seepage flux is often interpreted as ground water discharge, this discharging water may be of stream water origin that has recently entered the hyporheic zone.The measurement of seepage flux in conjunction with hyporheic water temperature or other indicators of water origin should be considered when attempting to quantify the magnitude of exchange and the source of hyporheic water.     

Delineating and quantifying ground water discharge zones using streambed temperatures

Streambed temperature mapping, hydraulic testing using minipiezometers, and geochemical analyses of interstitial water of the streambed were used to delineate the pattern of ground water discharge in a sandy streambed and to develop a flux-based conceptual model for ground water/surface water interactions. A new and simple empirical method was used to relate fluxes obtained from minipiezometer data to streambed temperatures. The relationship allowed flux to be calculated at locations where only streambed temperature measurements were made. Slug testing and potentiomanometer measurements at 34 piezometers indicated ground water discharge ranged from 0.03 to 446 L/m2/day (and possibly as high as 7060 L/m2/day) along a 60 m long by 11 to 14 m wide reach of river. Complex but similar plan-view patterns of flux were calculated for both summer and winter using hundreds of streambed temperatures measured on a 1 by 2 m grid. The reach was dominated by ground water discharge and 5% to 7% of the area accounted for approximately 20% to 24% of the total discharge. < 12% of the total area consisted of recharge zones or no-discharge zones. A conceptual model for ground water/surface water interactions consisting of five different behaviors was developed based on the magnitude and direction of flux across the surface of the streambed. The behaviors include short-circuit discharge (e.g., high-flow springs), high discharge (e.g., preferential flowpaths), low to moderate discharge, no discharge (e.g., horizontal hyporheic or ground water flow), and recharge. Geological variations at depth played a key role in determining which type of flow behavior occurred in the streambed.     

Use of an electromagnetic seepage meter to investigate temporal variability in lake seepage

A commercially available electromagnetic flowmeter is attached to a seepage cylinder to create an electromagnetic seepage meter (ESM) for automating measurement of fluxes across the sediment/water interface between ground water and surface water. The ESM is evaluated through its application at two lakes in New England, one where water seeps into the lake and one where water seeps out of the lake. The electromagnetic flowmeter replaces the seepage-meter bag and provides a continuous series of measurements from which temporal seepage processes can be investigated. It provides flow measurements over a range of three orders of magnitude, and contains no protruding components or moving parts. The ESM was used to evaluate duration of seepage disturbance following meter installation and indicated natural seepage rates resumed approximately one hour following meter insertion in a sandy lakebed. Lakebed seepage also varied considerably in response to lakebed disturbances, near-shore waves, and rainfalls, indicating hydrologic processes are occurring in shallow lakebed settings at time scales that have largely gone unobserved.     

Heat tracing of heterogeneous hyporheic exchange adjacent to in‐stream geomorphic features

There are many field techniques used to quantify rates of hyporheic exchange, which can vary in magnitude and direction spatially over distances of only a few metres, both within and between morphological features. We used in‐stream mini‐piezometers and heat transport modelling of stream and streambed temperatures to quantify the rates and directions of water flux across the streambed interface upstream and downstream of three types of in‐stream geomorphic features: a permanent dam, a beaver dam remnant and a stream meander. We derived hyporheic flux estimates at three different depths at six different sites for a month and then paired those flux rates with measurements of gradient to derive hydraulic conductivity (K) of the streambed sediments. Heat transport modelling provided consistent daily flux estimates that were in agreement directionally with hydraulic gradient measurements and also identified vertical heterogeneities in hydraulic conductivity that led to variable hyporheic exchange. Streambed K varied over an order of magnitude (1·9 × 10^?6 to 5·7 × 10^?5 m/s). Average rates of hyporheic flux ranged from static (q < ±0·02 m/day) to 0·42 m/day. Heat transport modelling results suggest three kinds of flow around the dams and the meander. Copyright © 2010 John Wiley & Sons, Ltd.     

Influence of a thin veneer of low‐hydraulic‐conductivity sediment on modelled exchange between river water and groundwater in response to induced infiltration

A thin layer of fine‐grained sediment commonly is deposited at the sediment–water interface of streams and rivers during low‐flow conditions, and may hinder exchange at the sediment–water interface similar to that observed at many riverbank‐filtration (RBF) sites. Results from a numerical groundwater‐flow model indicate that a low‐permeability veneer reduces the contribution of river water to a pumping well in a riparian aquifer to various degrees, depending on simulated hydraulic gradients, hydrogeological properties, and pumping conditions. Seepage of river water is reduced by 5–10% when a 2‐cm thick, low‐permeability veneer is present on the bed surface. Increasing thickness of the low‐permeability layer to 0·1 m has little effect on distribution of seepage or percentage contribution from the river to the pumping well. A three‐orders‐of‐magnitude reduction in hydraulic conductivity of the veneer is required to reduce seepage from the river to the extent typically associated with clogging at RBF sites. This degree of reduction is much larger than field‐measured values that were on the order of a factor of 20–25. Over 90% of seepage occurs within 12 m of the shoreline closest to the pumping well for most simulations. Virtually no seepage occurs through the thalweg near the shoreline opposite the pumping well, although no low‐permeability sediment was simulated for the thalweg. These results are relevant to natural settings that favour formation of a substantial, low‐permeability sediment veneer, as well as central‐pivot irrigation systems, and municipal water supplies where river seepage is induced via pumping wells. Published in 2011 by John Wiley & Sons, Ltd.     

Analytic solutions to transient groundwater flow under time-dependent sources in a heterogeneous aquifer bounded by fluctuating river stage

Analytical solutions for the water table and lateral discharge in a heterogeneous unconfined aquifer with time-dependent source and fluctuating river stage were derived and compared with those in an equivalent homogeneous aquifer. The heterogeneous aquifer considered consists of a number of sections of different hydraulic conductivity values. The source term and river stage were assumed to be time-dependent but spatially uniform. The solutions derived is useful in studying various groundwater flow problems in a horizontally heterogeneous aquifer since the spatially piecewise-constant hydraulic conductivity and temporally piecewise-constant recharge and lateral discharge can be used to quantify variations in these processes commonly observed in reality. Applying the solutions derived to an aquifer of three sections of different hydraulic conductivity values shown that (1) the aquifer heterogeneity significantly increases the spatial variation of the water table and thus its gradient but it has little effect on lateral discharge in the case of temporally and spatially uniform recharge, (2) the time-dependent but spatially uniform recharge increases the temporal variation of groundwater table over the entire aquifer but its effect on lateral discharge is limited in the zone near the river, and (3) the effect of river stage fluctuation on the water table and lateral discharge is limited in the zone near the river and the effect of the heterogeneity is to increase lateral discharge to or recharge from the river.     

Characterization of complex fluvial systems using remote sensing of spatial and temporal water level variations in the Amazon,Congo, and Brahmaputra Rivers

The Surface Water and Ocean Topography (SWOT) satellite mission will provide global, space‐based estimates of water elevation, its temporal change, and its spatial slope in fluvial environments, as well as across lakes, reservoirs, wetlands, and floodplains. This paper illustrates the utility of existing remote sensing measurements of water temporal changes and spatial slope to characterize two complex fluvial environments. First, repeat‐pass interferometric SAR measurements from the Japanese Earth Resources Satellite are used to compare and contrast floodplain processes in the Amazon and Congo River basins. Measurements of temporal water level changes over the two areas reveal clearly different hydraulic processes at work. The Amazon is highly interconnected by floodplain channels, resulting in complex flow patterns. In contrast, the Congo does not show similar floodplain channels and the flow patterns are not well defined and have diffuse boundaries. During inundation, the Amazon floodplain often shows sharp hydraulic changes across floodplain channels. The Congo, however, does not show similar sharp changes during either infilling or evacuation. Second, Shuttle Radar Topography Mission measurements of water elevation are used to derive water slope over the braided Brahmaputra river system. In combination with in situ bathymetry measurements, water elevation and slope allow one to calculate discharge estimates within 2.3% accuracy. These two studies illustrate the utility of satellite‐based measurements of water elevation for characterizing complex fluvial environments, and highlight the potential of SWOT measurements for fluvial hydrology. Copyright © 2010 John Wiley & Sons, Ltd.     

Groundwater-derived nutrient inputs to the Upper Gulf of Thailand

We report here the first direct measurements of nutrient fluxes via groundwater discharge into the Upper Gulf of Thailand. Nutrient and standard oceanographic surveys were conducted during the wet and dry seasons along the Chao Phraya River, Estuary and out into the Upper Gulf of Thailand. Additional measurements in selected near-shore regions of the Gulf included manual and automatic seepage meter deployments, as well as nutrient evaluations of seepage and coastal waters. The river transects characterized the distribution of biogeochemical parameters in this highly contaminated urban environment. Seepage flux measurements together with nutrient analyses of seepage fluids were used to estimate nutrient fluxes via groundwater pathways for comparison to riverine fluxes.     

Hydraulic visibility: Using satellite altimetry to parameterize a hydraulic model of an ungauged reach of a braided river

What hydraulic information can be gained from remotely sensed observations of a river's surface? In this study, we analyze the relationship between river bed undulations and water surfaces for an ungauged reach of the Xingu River, a first‐order tributary of the Amazon river. This braided reach is crosscut more than 10 times by a ENVISAT (ENVironmental SATellite) track that extends over 100 km. Rating curves based on a modeled discharge series and altimetric measurements are used, including the zero‐flow depth Z ₀ parameter, which describes river's bathymetry. River widths are determined from JERS (Japanese Earth Ressources Satellite) images. Hydrodynamic laws predict that irregularities in the geometry of a river bed produce spatial and temporal variations in the water level, as well as in its slope. Observation of these changes is a goal of the Surface Water and Ocean Topography satellite mission, which has a final objective of determining river discharge. First, the concept of hydraulic visibility is introduced, and the seasonality of water surface slope is highlighted along with different flow regimes and reach behaviors. Then, we propose a new single‐thread effective hydraulic approach for modeling braided rivers flows, based on the observation scales of current satellite altimetry. The effective hydraulic model is able to reproduce water surface elevations derived by satellite altimetry, and it shows that hydrodynamical signatures are more visible in areas where the river bed morphology varies significantly and for reaches with strong downstream control. The results of this study suggest that longitudinal variations of the slope might be an interesting criteria for the analysis of river segmentation into elementary reaches for the Surface Water Ocean Topography mission that will provide continuous measurements of the water surface elevations, the slopes, and the reach widths.     

Vertical Hydraulic Gradient Estimation in Clay Till,Using MiHPT Advanced Direct-Push Technology

The vertical transport of contaminants from source areas is employed in many risk assessment models and screening tools in order to estimate the contaminant mass discharge (CMD) into underlying aquifers. The key parameters for estimating CMD are the contaminant source area and concentration, and the vertical water flux, the latter of which depends on the vertical hydraulic conductivity and the vertical hydraulic gradient in the subsurface. This study focuses on advancing the use of the combined membrane interface probe hydraulic profiling tool (MiHPT) to investigate the vertical hydraulic gradient across a clay till overlying a sandy aquifer at a contaminated site in Denmark. Only the HPT is necessary for the estimate of vertical hydraulic gradient. The hydraulic head, clay till thickness, and resulting vertical hydraulic gradients found using the MiHPT compared well with observations from nearby nested wells. The parameter with the largest discrepancy was the thickness of the clay till. The advantage of the MiHPT is its relatively quick depth discrete access to information regarding subsurface permeability, vertical hydraulic gradients and contaminant distribution across a site. In this case study, performance of additional dissipationtests during the HPT log to acquire determination of the vertical hydraulic gradient increased the cost by 3% compared to standard HPT logs.     

Stochastic evolution of hydraulic geometry relations in the lower Yellow River of China under environmental uncertainties

Hydraulic geometry relations comprise a classic way to understand characteristics of a river. However, environmental changes pose large uncertainties for the reliability of such relations. In the current study, on the basis of the ordinary differential equations (ODEs) formed through linear treatment of the deterministic power-law hydraulic geometry relations, a set of stochastic differential equations (SDEs) driven by Fractional white noise and Poisson noise are developed to simulate the historical dynamic probability distributions of typical hydraulic geometry variables such as slope, width, depth, and velocity with bankfull discharge variation over time in the lower Yellow River of China. One group of possible stochastic average behaviors within the next 50 years are calculated under three different design incoming water-sediment conditions (including 300, 600, and 800 million t of annual average sediment discharge). In each part of the lower reaches, after estimation of the SDE parameters using a nonparametric maximum likelihood estimation (MLE) method, the model is carefully examined using Monte Carlo simulation as compared with the deterministic control models. The results of this comparison reveal the potential responses of hydraulic geometry characteristics to environmental disturbances, and the average trends mainly agree with the measurements. Comparisons among the three different prediction results reveal the stochastic average solution generally is greater than the deterministic solution. The results also confirm the severe negative impacts that result from the condition of 300 million t of incoming sediment, thus, pointing out the need to raise the level of river evolution alert for the lower Yellow River of China in the future. Moreover, with the help of the stochastic computation, the stream power and hydraulic width/depth ratio could be representative of an effective systematic measure for river dynamics. The proposed stochastic approach is not only important to development in the field of fluvial relations, but also beneficial to the practical design and monitoring of a river system according to specified accuracy requirements.     

Scales and magnitude of hyporheic,river–aquifer and bank storage exchange fluxes

Many studies have investigated the exchange processes that occur between rivers and groundwater systems and have successfully quantified the water fluxes involved. Specifically, these exchange processes include hyporheic exchange, river–aquifer exchange (groundwater discharge and river loss) and bank storage exchange. Remarkably, there are relatively few examples of field studies where more than one exchange process is quantified, and as a consequence, the relationships between them are not well understood. To compare the relative magnitudes of these common exchange processes, we have collected data from 54 studies that have quantified one or more of these exchange flux types. Each flux value is plotted against river discharge at the time of measurement to allow the different exchange flux types to be compared. We show that there are positive relationships between the magnitude of each exchange flux type and increasing river discharge across the different studies. For every one order of magnitude increase in river discharge, the hyporheic, river–aquifer and bank storage exchange fluxes increase by factors of 2.7, 2.9 and 2.5, respectively. On average, hyporheic exchange fluxes are almost an order of magnitude greater than river–aquifer exchange fluxes, which are, in turn, approximately four times greater than bank storage exchange fluxes for the same river discharge. Unless measurement approaches that can distinguish between different types of exchange flux are used, there is potential for hyporheic exchange fluxes to be misinterpreted as river–aquifer exchange fluxes, with possible implications for water resource management decisions. Copyright © 2015 John Wiley & Sons, Ltd.     

River-aquifer interactions, geologic heterogeneity, and low-flow management

Low river flows are commonly controlled by river-aquifer exchange, the magnitude of which is governed by hydraulic properties of both aquifer and aquitard materials beneath the river. Low flows are often important ecologically. Numerical simulations were used to assess how textural heterogeneity of an alluvial system influences river seepage and low flows. The Cosumnes River in California was used as a test case. Declining fall flows in the Cosumnes River have threatened Chinook salmon runs. A ground water-surface water model for the lower river basin was developed, which incorporates detailed geostatistical simulations of aquifer heterogeneity. Six different realizations of heterogeneity and a homogenous model were run for a 3-year period. Net annual seepage from the river was found to be similar among the models. However, spatial distribution of seepage along the channel, water table configuration and the level of local connection, and disconnection between the river and aquifer showed strong variations among the different heterogeneous models. Most importantly, the heterogeneous models suggest that river seepage losses can be reduced by local reconnections, even when the regional water table remains well below the riverbed. The percentage of river channel responsible for 50% of total river seepage ranged from 10% to 26% in the heterogeneous models as opposed to 23% in the homogeneous model. Differences in seepage between the models resulted in up to 13 d difference in the number of days the river was open for salmon migration during the critical fall months in one given year.