[1] Advection through hyporheic zones (HZ) consisting of heterogeneous channel bend streambed deposits and their equivalent homogenous medium was investigated using finite difference groundwater flow and transport simulations and forward particle tracking. The top prescribed head boundary was varied in order to mimic various stream channel head distributions resulting from the presence of bed forms and channel curvature. Flux calculations show that heterogeneity causes significant additional HZ flux compared to an equivalent homogenous medium. However, the major cause of HZ flux is a spatially periodic (sinusoidal) head distribution along the boundary, representing the effect of bed forms. The additional influence of heterogeneity on the total channel-bed exchange and the overall HZ geometry are increased when boundary head sinusoidal fluctuation is more subdued. We present dimensionless numbers that summarize these relationships. Heterogeneity's influence is further magnified by considering the effect of channel curvature on boundary heads. The simulations illustrate the dynamic influence of heterogeneity on the hyporheic zone since the various head boundaries employed in our modeling efforts are a proxy for different surface water conditions and bed form states that may occur during a single flood. Furthermore, we show that residence times (total tracking times) of particles originating from the streambed follow a lognormal distribution. In the presence of heterogeneity, residence times can decrease or they can increase compared to residence times for homogeneous conditions depending on the relative positions of the heterogeneities and the bed forms. Hence streambed heterogeneity and stream curvature, factors often neglected in previous modeling efforts, combine with bed form configuration to dynamically determine HZ geometry, fluxes, and residence time distributions.

[1] New semianalytical solutions for evaluation of the drawdown near horizontal and slanted wells with finite length screens in unconfined aquifers are presented. These fully three-dimensional solutions consider instantaneous drainage or delayed yield and aquifer anisotropy. As a basis, solution for the drawdown created by a point source in a uniform anisotropic unconfined aquifer is derived in Laplace domain. Using superposition, the point source solution is extended to the cases of the horizontal and slanted wells. The previous solutions for vertical wells can be described as a special case of the new solutions. Numerical Laplace inversion allows effective evaluation of the drawdown in real time. Examples illustrate the effects of well geometry and the aquifer parameters on drawdown. Results can be used to generate type curves from observations in piezometers and partially or fully penetrating observation wells. The proposed solutions and software are useful for the parameter identification, design of remediation systems, drainage, and mine dewatering.

Commonly used analytical approaches for estimation of pumping-induced drawdown and stream depletion are based on a series of idealistic assumptions about the stream-aquifer system. A new solution has been developed for estimation of drawdown and stream depletion under conditions that are more representative of those in natural systems (finite width stream of shallow penetration adjoining an aquifer of limited lateral extent). This solution shows that the conventional assumption of a fully penetrating stream will lead to a significant overestimation of stream depletion (> 100%) in many practical applications. The degree of overestimation will depend on the value of the stream leakance parameter and the distance from the pumping well to the stream. Although leakance will increase with stream width, a very wide stream will not necessarily be well represented by a model of a fully penetrating stream. The impact of lateral boundaries depends upon the distance from the pumping well to the stream and the stream leakance parameter. In most cases, aquifer width must be on the order of hundreds of stream widths before the assumption of a laterally infinite aquifer is appropriate for stream-depletion calculations. An important assumption underlying this solution is that stream-channel penetration is negligible relative to aquifer thickness. However, an approximate extension to the case of nonnegligible penetration provides reasonable results for the range of relative penetrations found in most natural systems (up to 85 % ). Since this solution allows consideration of a much wider range of conditions than existing analytical approaches, it could prove to be a valuable new tool for water management design and water rights adjudication purposes.

[1] Considerable progress has been made in developing a theoretical framework for modeling slug test responses in formations with high hydraulic conductivity K. However, several questions of practical significance remain unresolved. Given the rapid and often oscillatory nature of test responses, the traditional hydrostatic relationship between the water level and the transducermeasured head in the water column may not be appropriate. A general dynamic interpretation is proposed that describes the relationship between water level response and transducer-measured head. This theory is utilized to develop a procedure for transforming model-generated water level responses to transducer readings. The magnitude of the difference between the actual water level position and the apparent position based on the transducer measurement is a function of the acceleration and velocity of the water column, test geometry, and depth of the transducer. The dynamic approach explains the entire slug test response, including the often-noted discrepancy between the actual initial water level displacement and that measured by a transducer in the water column. Failure to use this approach can lead to a significant underestimation of K when the transducer is a considerable distance below the static water level. Previous investigators have noted a dependence of test responses on the magnitude of the initial water level displacement and have developed various approximate methods for analyzing such data. These methods are re-examined and their limitations clarified. Practical field guidelines are proposed on the basis of findings of this work. The soundness of the dynamic approach is demonstrated through a comparison of K profiles from a series of multilevel slug tests with those from dipole-flow tests performed in the same wells.

[1] We present a three-dimensional model of heterogeneous modern channel bend deposits developed through purely structure-imitating interpolation (kriging) of hydraulic properties. This model, augmented with ground-penetrating radar data and directional variograms, agrees with detailed observations in similar modern environments and leads to a process-based interpretation of the presented hydraulic conductivity structure. Integration of all available information permitted delineation and characterization of the modern streambed as a distinct hydrostratigraphic unit without coring or outcrop studies. Our results imply that the modern streambed is commonly oversimplified in available analytical and numerical models of groundwater-surface water interactions where it is assumed to be homogeneous and isotropic and characterized by a constant width and thickness. This three-dimensional approach that integrates concepts and principles developed in sedimentology, hydrogeology, geophysics, and geostatistics has potential implications on model development of stream-aquifer systems.