Geothermal gradients and ground water circulation in fissured and karstic rocks: The role played by the structure of the permeable network

In fissured and karstic rocks the general movement of underground waters (forced convection) can modify geothermic gradients. This depends both on the discontinuous structure (channels and fissures) and on hydrodynamic conditions which can vary with the weather, e.g. during the recharging of reserviors in rainy periods.An experimental analysis has been carried out in the broken and karstified Mesozoic limestone in the South of France, on shallow boreholes (60 m) grouped in a closely-spaced network. Nearly a hundred thermal loggings have been measured in the homothermic zone below 25 m. The gradients in dry periods, varying from one drilling to another, are between 0.01 and 0.03°C m⁻¹ for an average thermal conductivity of rock of 2.56 Wm⁻¹ °C⁻¹. During recharging of the aquifer by rain, the gradients do not change in some drillings. This always occurs in those which cut through networks of slightly karstified fissures with low hydraulic conductivity. The slow circulation allows the water to be in thermal quasi-equilibrium with the rock. In other drillings, however, recharging causes local and sometimes very significant modifications of the gradients. Disturbances are temporary and appear directly over well-developed karstic channels which rapidly draw down the infiltrated cold water to the bottom. Thermal profiles, either stable or disturbed, can be surveyed simultaneously in drillings situated at least 10 m from each other. The position and nature of the karstic channels in which the forced convection is most active can be identified through observations by videologging and flow velocity tests.     

Climate,vegetation, and soil controls on hydraulic redistribution in shallow tree roots

Hydraulic redistribution defined as the translocation of soil moisture by plant root systems in response to water potential gradients is a phenomenon widely documented in different climate, vegetation, and soil conditions. Past research has largely focused on hydraulic redistribution in deep tree roots with access to groundwater and/or winter rainfall, while the case of relatively shallow (i.e., ≈1–2 m deep) tree roots has remained poorly investigated. In fact, it is not clear how hydraulic redistribution in shallow root zones is affected by climate, vegetation, and soil properties. In this study, we developed a model to investigate the climate, vegetation, and soil controls on the net direction and magnitude of hydraulic redistribution in shallow tree root systems at the growing season to yearly timescale. We used the model to evaluate the effect of hydraulic redistribution on the water stress of trees and grasses. We found that hydraulic lift increases with decreasing rainfall frequency, depth of the rooting zone, root density in the deep soil and tree leaf area index; at the same time for a given rainfall frequency, hydraulic lift increases with increasing average rainstorm depth and soil hydraulic conductivity. We propose that water drainage into deeper soil layers can lead to the emergence of vertical water potential gradients sufficient to explain the occurrence of hydraulic lift in shallow tree roots without invoking the presence of a shallow water table or winter precipitation. We also found that hydraulic descent reduces the water stress of trees and hydraulic lift reduces the water stress of grass with important implications on tree–grass interactions.     

Effects of eastern redcedar encroachment on soil hydraulic properties along Oklahoma's grassland‐forest ecotone

In north‐central Oklahoma eastern redcedar (Juniperus virginiana), encroachment into grassland is widespread and is suspected of reducing streamflow, but the effects of this encroachment on soil hydraulic properties are unknown. This knowledge gap creates uncertainty in understanding the hydrologic effects of eastern redcedar encroachment and obstructs fact‐based management of encroached systems. The objective of this study was to quantify the effects of eastern redcedar encroachment into tallgrass prairie on soil hydraulic properties. Leaf litter depth, soil organic matter, soil water repellency, soil water content, sorptivity, and unsaturated hydraulic conductivity were measured near Stillwater, OK, along 12 radial transects from eastern redcedar trunks to the center of the grassy intercanopy space. Eastern redcedar encroachment in the second half of the 20th century caused the accumulation of 3 cm of hydrophobic leaf litter near the trunks of eastern redcedar trees. This leaf litter was associated with increased soil organic matter in the upper 6 cm of soil under eastern redcedar trees (5.96% by mass) relative to the grass‐dominated intercanopy area (3.99% by mass). Water repellency was more prevalent under eastern redcedar than under grass, and sorptivity under eastern redcedar was 0.10 mm s^?1/2, one seventh the sorptivity under adjacent prairie grasses (0.68 mm s^?1/2). Median unsaturated hydraulic conductivity under grass was 2.52 cm h^?1, four times greater than under eastern redcedar canopies (0.57 cm h^?1). Lower sorptivity and unsaturated hydraulic conductivity would tend to decrease infiltration and increase runoff, but other factors such as rainfall interception by the eastern redcedar canopy and litter layer, and preferential flow induced by hydrophobicity must be examined before the effects of encroachment on streamflow can be predicted. Copyright © 2011 John Wiley & Sons, Ltd.     

Effects of experimental uncertainty on the calculation of hillslope flow paths

Measurement uncertainty is a key hindrance to the quantification of water fluxes at all scales of investigation. Predictions of soil‐water flux rely on accurate or representative measurements of hydraulic gradients and field‐state hydraulic conductivity. We quantified the potential magnitude of errors associated with the parameters and variables used directly and indirectly within the Darcy – Buckingham soil‐water‐flux equation. These potential errors were applied to a field hydrometric data set collected from a forested hillslope in central Singapore, and their effect on flow pathway predictions was assessed. Potential errors in the hydraulic gradient calculations were small, approximately one order of magnitude less than the absolute magnitude of the hydraulic gradients. However, errors associated with field‐state hydraulic conductivity derivation were very large. Borehole (Guelph permeameter) and core‐based (Talsma ring permeameter) techniques were used to measure field‐saturated hydraulic conductivity. Measurements using these two approaches differed by up to 3\9 orders of magnitude, with the difference becoming increasingly marked within the B horizon. The sensitivity of the shape of the predicted unsaturated hydraulic conductivity curve to ±5% moisture content error on the moisture release curve was also assessed. Applied moisture release curve error resulted in hydraulic conductivity predictions of less than ±0\2 orders of magnitude deviation from the apparent conductivity. The flow pathways derived from the borehole saturated hydraulic conductivity approach suggested a dominant near‐surface flow pathway, whereas pathways calculated from the core‐based measurements indicated vertical percolation to depth. Direct tracer evidence supported the latter flow pathway, although tracer velocities were approximately two orders of magnitude smaller than the Darcy predictions. We conclude that saturated hydraulic conductivity is the critical hillslope hydrological parameter, and there is an urgent need to address the issues regarding its measurement further. Copyright © 2000 John Wiley & Sons, Ltd.     

Dynamic river–groundwater exchange in the presence of a saline,semi‐confined aquifer

Understanding groundwater–surface water exchange in river banks is crucial for effective water management and a range of scientific disciplines. While there has been much research on bank storage, many studies assume idealized aquifer systems. This paper presents a field‐based study of the Tambo Catchment (southeast Australia) where the Tambo River interacts with both an unconfined aquifer containing relatively young and fresh groundwater (<500 μS/cm and <100 years old) and a semi‐confined artesian aquifer containing old and saline groundwater (electrical conductivity > 2500 μS/cm and >10 000 years old). Continuous groundwater elevation and electrical conductivity monitoring within the different aquifers and the river suggest that the degree of mixing between the two aquifers and the river varies significantly in response to changing hydrological conditions. Numerical modelling using MODFLOW and the solute transport package MT3DMS indicates that saline water in the river bank moves away from the river during flooding as hydraulic gradients reverse. This water then returns during flood recession as baseflow hydraulic gradients are re‐established. Modelling also indicates that the concentration of a simulated conservative groundwater solute can increase for up to ~34 days at distances of 20 and 40 m from the river in response to flood events approximately 10 m in height. For the same flood event, simulated solute concentrations within 10 m of the river increase for only ~15 days as the infiltrating low‐salinity river water drives groundwater dilution. Average groundwater fluxes to the river stretch estimated using Darcy's law were 7 m³/m/day compared with 26 and 3 m³/m/day for the same periods via mass balance using Radon (²²²Rn) and chloride (Cl), respectively. The study shows that by coupling numerical modelling with continuous groundwater–surface water monitoring, the transient nature of bank storage can be evaluated, leading to a better understanding of the hydrological system and better interpretation of hydrochemical data. Copyright © 2015 John Wiley & Sons, Ltd.     

Processes of water movement through a chalk coombe deposit in Southeast England

The processes of water movement through the Coombe Deposit in a chalk dry valley near Eastbourne in Southeast England were investigated using simple methods based on regular weekly measurements of rainfall, soil water content, and soil water potential. The drainage flux (recharge) through the soil was determined using the water balance method during the winter and the zero flux plane (ZFP) method after the appearance of the ZFP in the spring. The unsaturated hydraulic conductivity was derived applying Darcy's Law in a novel way using the measured potential gradients and weekly drainage fluxes. The derived conductivity characteristics were adequate to identify the flow mechanisms, to illustrate the difference in behaviour between the horizons of the soil profile, and to give some indication of pore water velocities. The mean daily drainage flux at 2.85 m depth during the recharge period from 10 October 1980 to 29 May 1981 was 1.6 mm d^?1. Weekly mean rates of up to 3.7 mm d^?1 were observed, but peak short term rates must have considerably exceeded this figure. It was shown that, in the lower part of the Coombe Deposit, when drainage fluxes are large, much of the flux passes through a very small proportion of the wetted cross-sectional area of the soil. This gives rise to pore water velocities of at least 3 m d^?1 at a depth of 2.85 m and 0.5 m d^?1 between 0.5 m and 2.5 m depth. These results show that pollutants may be moved very rapidly through the profile in this, and similar, material. The core sampling techniques normally used to monitor pollutant movement in the chalk are unlikely to succeed in detecting this movement, not only because it is transient but also because it occupies only a very small proportion of the water filled pores.     

A Simple Field Method for Assessing Near‐Surface Saturated Hydraulic Conductivity

This technical note describes an effective and inexpensive field technique for measuring the saturated hydraulic conductivity of both undisturbed cores and repacked soil samples. The method requires no specialized equipment; everything that is required can be obtained in a hardware store. The method is a straightforward field implementation of the widely used falling‐head laboratory analysis directly derived from Darcy's law. As such, it sidesteps the need for empirical assumptions about soil texture and the relationship between saturated and unsaturated flow components which many permeameter‐based methods rely upon. The method is shown to produce results that are consistent with K values obtained elsewhere in the same homogeneous sand formation. Furthermore, the proposed method is useful for measuring hydraulic conductivity in drill cuttings obtained from direct push or auguring drill techniques, which cannot be done with any other field method. The range of hydraulic conductivity values that this test is appropriate for is on the order of 1E ? 7 m/s to 1E ? 3 m/s.     

Influence of Degree of Saturation in the Electric Resistivity–Hydraulic Conductivity Relationship

The relationship between aquifer hydraulic conductivity and aquifer resistivity, either measured on the ground surface by vertical electrical sounding (VES) or from resistivity logs, or measured in core samples have been published for different types of aquifers in different locations. Generally, these relationships are empirical and semi-empirical, and confined in few locations. This relation has a positive correlation in some studies and negative in others. So far, there is no potentially physical law controlling this relation, which is not completely understood. Electric current follows the path of least resistance, as does water. Within and around pores, the model of conduction of electricity is ionic and thus the resistivity of the medium is controlled more by porosity and water conductivity than by the resistivity of the rock matrix. Thus, at the pore level, the electrical path is similar to the hydraulic path and the resistivity should reflect hydraulic conductivity. We tried in this paper to study the effect of degree of groundwater saturation in the relation between hydraulic conductivity and bulk resistivity via a simple numerical analysis of Archie’s second law and a simplified Kozeny-Carmen equation. The study reached three characteristic non-linear relations between hydraulic conductivity and resistivity depending on the degree of saturation. These relations are: (1) An inverse power relation in fully saturated aquifers and when porosity equals water saturation, (2) An inverse polynomial relation in unsaturated aquifers, when water saturation is higher than 50%, higher than porosity, and (3) A direct polynomial relation in poorly saturated aquifers, when water saturation is lower than 50%, lower than porosity. These results are supported by some field scale relationships.     

SEAWAT‐Based Simulation of Axisymmetric Heat Transport

Simulation of heat transport has its applications in geothermal exploitation of aquifers and the analysis of temperature dependent chemical reactions. Under homogeneous conditions and in the absence of a regional hydraulic gradient, groundwater flow and heat transport from or to a well exhibit radial symmetry, and governing equations are reduced by one dimension (1D) which increases computational efficiency importantly. Solute transport codes can simulate heat transport and input parameters may be modified such that the Cartesian geometry can handle radial flow. In this article, SEAWAT is evaluated as simulator for heat transport under radial flow conditions. The 1971, 1D analytical solution of Gelhar and Collins is used to compare axisymmetric transport with retardation (i.e., as a result of thermal equilibrium between fluid and solid) and a large diffusion (conduction). It is shown that an axisymmetric simulation compares well with a fully three dimensional (3D) simulation of an aquifer thermal energy storage systems. The influence of grid discretization, solver parameters, and advection solution is illustrated. Because of the high diffusion to simulate conduction, convergence criterion for heat transport must be set much smaller (10^?10) than for solute transport (10^?6). Grid discretization should be considered carefully, in particular the subdivision of the screen interval. On the other hand, different methods to calculate the pumping or injection rate distribution over different nodes of a multilayer well lead to small differences only.     

How Does Subsurface Characterization Affect Simulations of Hyporheic Exchange?

We investigated the role of increasingly well‐constrained geologic structures in the subsurface (i.e., subsurface architecture) in predicting streambed flux and hyporheic residence time distribution (RTD) for a headwater stream. Five subsurface realizations with increasingly resolved lithological boundaries were simulated in which model geometries were based on increasing information about flow and transport using soil and geologic maps, surface observations, probing to depth to refusal, seismic refraction, electrical resistivity (ER) imaging of subsurface architecture, and time‐lapse ER imaging during a solute tracer study. Particle tracking was used to generate RTDs for each model run. We demonstrate how improved characterization of complex lithological boundaries and calibration of porosity and hydraulic conductivity affect model prediction of hyporheic flow and transport. Models using hydraulic conductivity calibrated using transient ER data yield estimates of streambed flux that are three orders of magnitude larger than uncalibrated models using estimated values for hydraulic conductivity based on values published for nearby hillslopes (10^?4 vs. 10^?7 m²/s, respectively). Median residence times for uncalibrated and calibrated models are 10³ and 10⁰ h, respectively. Increasingly well‐resolved subsurface architectures yield wider hyporheic RTDs, indicative of more complex hyporheic flowpath networks and potentially important to biogeochemical cycling. The use of ER imaging to monitor solute tracers informs subsurface structure not apparent from other techniques, and helps to define transport properties of the subsurface (i.e., hydraulic conductivity). Results of this study demonstrate the value of geophysical measurements to more realistically simulate flow and transport along hyporheic flowpaths.     

Identification of optimal soil hydraulic functions and parameters for predicting soil moisture

Abstract

The accuracy of six combined methods formed by three commonly-used soil hydraulic functions and two methods to determine soil hydraulic parameters based on a soil hydraulic parameter look-up table and soil pedotransfer functions was examined for simulating soil moisture. A novel data analysis and modelling approach was used that eliminated the effects of evapotranspiration so that specific sources of error among the six combined methods could be identified and quantified. By comparing simulated and observed soil moisture at six sites of the USDA Soil Climate Analysis Network, we identified the optimal soil hydraulic functions and parameters for predicting soil moisture. Through sensitivity tests, we also showed that adjusting only the soil saturated hydraulic conductivity, K_s , is insufficient for representing important effects of macropores on soil hydraulic conductivity. Our analysis illustrates that, in general, soil hydraulic conductivity is less sensitive to K_s than to the soil pore-size distribution parameter.

Editor D. Koutsoyiannis; Associate editor D. Hughes

Citation Pan, F., McKane, R.B. and Stieglitz, M., 2012. Identification of optimal soil hydraulic functions and parameters for predicting soil moisture. Hydrological Sciences Journal, 57 (4), 723–737.     

Determining hydraulic properties of a loam soil by alternative infiltrometer techniques

Testing infiltrometer techniques to determine soil hydraulic properties is necessary for specific soils. For a loam soil, the water retention and hydraulic conductivity predicted by the BEST (Beerkan Estimation of Soil Transfer parameters) procedure of soil hydraulic characterization was compared with data collected by more standard laboratory and field techniques. Six infiltrometer techniques were also compared in terms of saturated soil hydraulic conductivity, K_s. BEST yielded water retention values statistically similar to those obtained in the laboratory and K_s values practically coinciding with those determined in the field with the pressure infiltrometer (PI). The unsaturated soil hydraulic conductivity measured with the tension infiltrometer (TI) was reproduced satisfactorily by BEST only close to saturation. BEST, the PI, one‐potential experiments with both the TI and the mini disk infiltrometer (MDI), the simplified falling head (SFH) technique and the bottomless bucket (BB) method yielded statistically similar estimates of K_s, differing at the most by a factor of three. Smaller values were obtained with longer and more soil‐disturbing infiltration runs. Any of the tested infiltration techniques appears usable to obtain the order of magnitude of K_s at the field site, but the BEST, BB and PI data appear more appropriate to characterize the soil at some stage during a rainfall event. Additional investigations on both similar and different soils would allow development of more general procedures to apply infiltrometer techniques for soil hydraulic characterization. Copyright © 2015 John Wiley & Sons, Ltd.     

Matrix Diffusion Effects on Contaminant Migration from an Injection Well in Fractured Sandstone

Deep-well injection into fractured sandstone is an option for the disposal of contaminated mine dewatering discharge from an open pit uranium mine. As part of the assessment of potential contaminant migration from deep-well injection, the effect of matrix diffusion was evaluated. An analytical mathematical model was developed for the simulation of the radial movement of a contaminant front away from an injection point under steady flow conditions in a planar fracture with uniform properties. The model includes the effects of advection in the fracture, diffusion of contaminants from the fracture into the rock matrix, and equilibrium adsorption on the fracture surface as well as in the rock matrix. Effective diffusion coefficients obtained from laboratory experiments on 11 intact core samples varied from 3.4 × 10⁻⁸ to 3.2 × 10⁻⁷ cm²/s. Model simulations were made with diffusion coefficient values in this range and with single-fracture injection rates estimated from fracture frequencies in boreholes, and from bulk hydraulic conductivity values obtained from field tests. Because of matrix diffusion, the rate of outward movement of the front of the nonreactive contaminants from the injection well is much slower than the rate of water flow in the fractures. Simulations of the movement of contaminants that undergo adsorption indicate that even a small distribution coefficient for the rock matrix causes the contaminants to remain very close to the injection well during the one-year period. The results of the simplified model demonstrate that matrix diffusion is an important process that cannot be neglected in the assessment of a waste disposal scheme located in fractured porous rock. However, in order to make a definitive assessment of the capability of matrix diffusion and associated matrix adsorption to significantly limit the extent of contaminant migration around injection wells, it would be necessary to conduct field tests such as a preliminary or experimental injection.     

Hydraulic conductivity in frozen soils

An apparatus has been developed for investigation of hydraulic conductivity of frozen soils. The test procedure is isothermal and involves the passage of water from one reservoir into the frozen sample and out of the frozen sample into a second reservoir. The water in the reservoirs remains unfrozen because it contains dissolved lactose. The concentration of lactose is such that, initially, the water in the reservoirs is in thermodynamic equilibrium with the water in the soil. On application of pressure to one reservoir a known hydraulic gradient is established and flow takes place. Flow is shown to vary linearly with hydraulic gradient. The hydraulic conductivity coefficient depends on soil type and temperature and is related to the unfrozen water content. At temperatures within a few tenths of 0°C the coefficient apparently ranges from 10^?5 to 10^?9 cm sec^?1, and decreases only slowly below about ?0·5°C. Soils known to be susceptible to frost heave are shown to have significant hydraulic conductivities well below 0°C.     

The exponential decline in saturated hydraulic conductivity with depth: a novel method for exploring its effect on water flow paths and transit time distribution

The strong vertical gradient in soil and subsoil saturated hydraulic conductivity is characteristic feature of the hydrology of catchments. Despite the potential importance of these strong gradients, they have proven difficult to model using robust physically based schemes. This has hampered the testing of hypotheses about the implications of such vertical gradients for subsurface flow paths, residence times and transit time distribution. Here we present a general semi‐analytical solution for the simulation of 2D steady‐state saturated‐unsaturated flow in hillslopes with saturated hydraulic conductivity that declines exponentially with depth. The grid‐free solution satisfies mass balance exactly over the entire saturated and unsaturated zones. The new method provides continuous solutions for head, flow and velocity in both saturated and unsaturated zones without any interpolation process as is common in discrete numerical schemes. This solution efficiently generates flow pathlines and transit time distributions in hillslopes with the assumption of depth‐varying saturated hydraulic conductivity. The model outputs reveal the pronounced effect that changing the strength of the exponential decline in saturated hydraulic conductivity has on the flow pathlines, residence time and transit time distribution. This new steady‐state model may be useful to others for posing hypotheses about how different depth functions for hydraulic conductivity influence catchment hydrological response. Copyright © 2016 John Wiley & Sons, Ltd.     

Impact of land use on the hydraulic properties of the topsoil in a small French catchment

The hydraulic properties of the topsoil control the partition of rainfall into infiltration and runoff at the soil surface. They must be characterized for distributed hydrological modelling. This study presents the results of a field campaign documenting topsoil hydraulic properties in a small French suburban catchment (7 km²) located near Lyon, France. Two types of infiltration tests were performed: single ring infiltration tests under positive head and tension‐disk infiltration using a mini‐disk. Both categories were processed using the BEST—Beerkan Estimation of Soil Transfer parameters—method to derive parameters describing the retention and hydraulic conductivity curves. Dry bulk density and particle size data were also sampled. Almost all the topsoils were found to belong to the sandy loam soil class. No significant differences in hydraulic properties were found in terms of pedologic units, but the results showed a high impact of land use on these properties. The lowest dry bulk density values were obtained in forested soils with the highest organic matter content. Permanent pasture soils showed intermediate values, whereas the highest values were encountered in cultivated lands. For saturated hydraulic conductivity, the highest values were found in broad‐leaved forests and small woods. The complementary use of tension‐disk and positive head infiltration tests highlighted a sharp increase of hydraulic conductivity between near saturation and saturated conditions, attributed to macroporosity effect. The ratio of median saturated hydraulic conductivity to median hydraulic conductivity at a pressure of − 20 mm of water was about 50. The study suggests that soil texture, such as used in most pedo‐transfer functions, might not be sufficient to properly map the variability of soil hydraulic properties. Land use information should be considered in the parameterizations of topsoil within hydrological models to better represent in situ conditions, as illustrated in the paper. Copyright © 2010 John Wiley & Sons, Ltd.     

Modeling changes in the permeability in a gypsified fractured-porous rock massif

A method is proposed for evaluating rock hydraulic conductivity (k _x;y , LT^?1) in grid model blocks to be used in finite-difference simulation of mass exchange in gypcified fractured-porous rocks containing groundwater. The values of k _x and k _y were determined taking into account the dominating form of transport (advection or transverse dispersion), contributing most to gypsum dissolution. The method was used to simulate groundwater flow in a vertical section of gypsified rocks in the base of the dam at the Lower-Kafirniganskii Hydropower System. The changes in rock permeability in the base of the dam were used to evaluate whether there is a need to take seepage-control measures (cement-grout curtain, blanket, and a gypsum fill over the blanket). A leaching zone was found to form near the curtain and blanket during the simulation period (100 years) under background seepage conditions. In this period, the total water flow in the lower pool can increase by no more than 10–15%. Under such conditions, the “seepage-control” measures can enhance the reliability of the dam at the Lower-Kafirniganskii Hydropower System.     

ALLUVIAL CHARACTERISTICS,GROUNDWATER–SURFACE WATER EXCHANGE AND HYDROLOGICAL RETENTION IN HEADWATER STREAMS

Conservative solute injections were conducted in three first-order montane streams of different geological composition to assess the influence of parent lithology and alluvial characteristics on the hydrological retention of nutrients. Three study sites were established: (1) Aspen Creek, in a sandstone–siltstone catchment with a fine-grained alluvium of low hydraulic conductivity (1·3×10⁻⁴ cm/s), (2) Rio Calaveras, which flows through volcanic tuff with alluvium of intermediate grain size and hydraulic conductivity (1·2×10⁻³ cm/s), and (3) Gallina Creek, located in a granite/gneiss catchment of coarse, poorly sorted alluvium with high hydraulic conductivity (4·1×10⁻³ cm/s). All sites were instrumented with networks of shallow groundwater wells to monitor interstitial solute transport. The rate and extent of groundwater–surface water exchange, determined by the solute response in wells, increased with increasing hydraulic conductivity. The direction of surface water–groundwater interaction within a stream was related to local variation in vertical and horizontal hydraulic gradients. Experimental tracer responses in the surface stream were simulated with a one-dimensional solute transport model with inflow and storage components (OTIS). Model-derived measures of hydrological retention showed a corresponding increase with increasing hydraulic conductivity. To assess the temporal variability of hydrological retention, solute injection experiments were conducted in Gallina Creek under four seasonal flow regimes during which surface discharge ranged from baseflow (0·75 l/s in October) to high (75 l/s during spring snowmelt). Model-derived hydrological retention decreased with increasing discharge. The results of our intersite comparison suggest that hydrological retention is strongly influenced by the geologic setting and alluvial characteristics of the stream catchment. Temporal variation in hydrological retention at Gallina Creek is related to seasonal changes in discharge, highlighting the need for temporal resolution in studies of the dynamics of surface water–groundwater interactions in stream ecosystems. © 1997 by John Wiley & Sons, Ltd.