An experimental study of water table recharge by seepage losses from a ditch with intermittent flow

Farmed catchments in the Mediterranean area often exhibit dense networks of ditches which are also preferential zones of water table recharge, and thereby of groundwater contamination. This study presents an experimental analysis of seepage losses and related groundwater recharge patterns during a typical Mediterranean runoff event at the scale of a ditch located above a shallow water table. The objectives were (i) to evaluate the patterns of water table recharge by seepage in a ditch, (ii) to study the main flow processes occurring during recharge, and (iii) to estimate solute propagation in case of contaminated flow in the ditch. The field observation indicated three major points. Firstly, they showed that seepage losses during a runoff event in a ditch can rapidly lead to a significant recharge of a shallow water table. Secondly, the recharge induces a groundwater mound much larger than the event plume. The infiltrated water and the accompanying solutes remained in the vicinity of the ditch. The patterns of groundwater recharge and contamination appeared very different. Lastly, both unsaturated and saturated‐piston flow processes were observed which suggests that a variably‐saturated flow modelling approach ought to be used to simulate the ditch‐water shallow table interaction. Finally, the study indicates that the patterns of water table recharge and contamination in Mediterranean catchments with dense ditches network vary largely in space and time, and will require dense monitoring networks to estimate the evolution of the average contamination levels. Copyright © 2008 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.     

Modelling subsurface flow conditions in a salinized catchment in south‐western Australia,with a view to improving management practices

Finite element modelling of the saturated–unsaturated surface–subsurface flow mechanisms operative in a small salinized catchment in south‐western Australia was used to help define the flow system and explain the causes of waterlogging and salinization there. Data available at the site from a previous study were used to obtain a first approximation to the flow system. Altering the properties of some of the strata gave a closer calibration. It was found that the modelled saturated hydraulic conductivity of the B horizon in the duplex soil zone needed to be at least an order of magnitude lower than that measured in order to reproduce the perching conditions observed in the field. Also, the model indicated the influence of a doleritic dyke, whose presence was confirmed by field measurement. Our analysis showed that there were two main flow systems operating in the hillslope. The first, and most dominant, was the recharge occurring through the upslope gradational soil zone and percolating down to both the deeply weathered regolith and the basal aquifer. The second flow system is an unsaturated flow system operating in the high permeability A horizon in the downslope duplex soil zone. The first system is primarily responsible for the saline seepage zone in the valley bottom. The second contributes to the waterlogging and perching occurring upslope of the seepage zone. Vertical flow through the higher permeability B horizon in the gradational soil zone in the upper slopes is a major contributor of recharge. Recharge by flow through macropores occurs where, but only where, perched aquifers develop and allow the macropores to be activated. Areas with perched aquifers occurred in downslope locations and near a doleritic dyke located upslope. Thus, the area where macropore recharge occurred was not large. The recharge rate required to maintain the piezometric levels at present values is only about 30 mm/yr (about 5% of the annual rainfall). The piezometric levels under the upper part of the catchment varied greatly with only small changes in recharge rate. A 50% reduction in recharge rate had the effect of reducing the length of the seepage zone at the end of winter by 40%. Changes in recharge rate had little effect on the extent of the perched aquifer at the end of winter. Deep‐rooted perennial forages, shrubs or trees on the gradational soil zone in the upper part of the catchment and on the zones upslope of geological barriers to flow would be required to reduce the recharge and to allow for rehabilitation of the saline valley floor. Waterlogging associated with the perched water table in the bottom part of the catchment would be best addressed by tree plantations located just upslope of the salinized zone in the valley floor. Copyright © 1999 John Wiley & Sons, Ltd.     

Aquifer response to stream-stage and recharge variations. II. Convolution method and applications

In this second of two papers, analytical step-response functions, developed in the companion paper for several cases of transient hydraulic interaction between a fully penetrating stream and a confined, leaky, or water-table aquifer, are used in the convolution integral to calculate aquifer heads, streambank seepage rates, and bank storage that occur in response to stream-stage fluctuations and basinwide recharge or evapotranspiration. Two computer programs developed on the basis of these step-response functions and the convolution integral are applied to the analysis of hydraulic interaction of two alluvial stream–aquifer systems in the northeastern and central United States. These applications demonstrate the utility of the analytical functions and computer programs for estimating aquifer and streambank hydraulic properties, recharge rates, streambank seepage rates, and bank storage. Analysis of the water-table aquifer adjacent to the Blackstone River in Massachusetts suggests that the very shallow depth of water table and associated thin unsaturated zone at the site cause the aquifer to behave like a confined aquifer (negligible specific yield). This finding is consistent with previous studies that have shown that the effective specific yield of an unconfined aquifer approaches zero when the capillary fringe, where sediment pores are saturated by tension, extends to land surface. Under this condition, the aquifer's response is determined by elastic storage only. Estimates of horizontal and vertical hydraulic conductivity, specific yield, specific storage, and recharge for a water-table aquifer adjacent to the Cedar River in eastern Iowa, determined by the use of analytical methods, are in close agreement with those estimated by use of a more complex, multilayer numerical model of the aquifer. Streambank leakance of the semipervious streambank materials also was estimated for the site. The streambank-leakance parameter may be considered to be a general (or lumped) parameter that accounts not only for the resistance of flow at the river–aquifer boundary, but also for the effects of partial penetration of the river and other near-stream flow phenomena not included in the theoretical development of the step-response functions.     

Quantifying River‐Groundwater Interactions of New Zealand's Gravel‐Bed Rivers: The Wairau Plain

New Zealand's gravel‐bed rivers have deposited coarse, highly conductive gravel aquifers that are predominantly fed by river water. Managing their groundwater resources is challenging because the recharge mechanisms in these rivers are poorly understood and recharge rates are difficult to predict, particularly under a more variable future climate. To understand the river‐groundwater exchange processes in gravel‐bed rivers, we investigate the Wairau Plain Aquifer using a three‐dimensional groundwater flow model which was calibrated using targeted field observations, “soft” information from experts of the local water authority, parameter regularization techniques, and the model‐independent parameter estimation software PEST. The uncertainty of simulated river‐aquifer exchange flows, groundwater heads, spring flows, and mean transit times were evaluated using Null‐space Monte‐Carlo methods. Our analysis suggests that the river is hydraulically perched (losing) above the regional water table in its upper reaches and is gaining downstream where marine sediments overlay unconfined gravels. River recharge rates are on average 7.3 m³/s, but are highly dynamic in time and variable in space. Although the river discharge regularly hits 1000 m³/s, the net exchange flow rarely exceeds 12 m³/s and seems to be limited by the physical constraints of unit‐gradient flux under disconnected rivers. An important finding for the management of the aquifer is that changes in aquifer storage are mainly affected by the frequency and duration of low‐flow periods in the river. We hypothesize that the new insights into the river‐groundwater exchange mechanisms of the presented case study are transferable to other rivers with similar characteristics.     

Interaction of Aquifer and River‐Canal Network near Well Field

The article presents semi‐analytical mathematical models to asses (1) enhancements of seepage from a canal and (2) induced flow from a partially penetrating river in an unconfined aquifer consequent to groundwater withdrawal in a well field in the vicinity of the river and canal. The nonlinear exponential relation between seepage from a canal reach and hydraulic head in the aquifer beneath the canal reach is used for quantifying seepage from the canal reach. Hantush's (1967) basic solution for water table rise due to recharge from a rectangular spreading basin in absence of pumping well is used for generating unit pulse response function coefficients for water table rise in the aquifer. Duhamel's convolution theory and method of superposition are applied to obtain water table position due to pumping and recharge from different canal reaches. Hunt's (1999) basic solution for river depletion due to constant pumping from a well in the vicinity of a partially penetrating river is used to generate unit pulse response function coefficients. Applying convolution technique and superposition, treating the recharge from canal reaches as recharge through conceptual injection wells, river depletion consequent to variable pumping and recharge is quantified. The integrated model is applied to a case study in Haridwar (India). The well field consists of 22 pumping wells located in the vicinity of a perennial river and a canal network. The river bank filtrate portion consequent to pumping is quantified.     

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.     

Three-dimensional flow patterns at the river–aquifer interface — a case study at the Danube

The three-dimensional groundwater flow patterns in a gravel bar at the Danube east of Vienna were investigated and are discussed in this paper. The observed groundwater level gradients are highly dynamic and respond very quickly to changes in the river water levels. A variably saturated groundwater model was calibrated to the data to describe the complex dynamics of flow in the gravel bar. The model results suggest that short-term (6–48 h) fluctuations of river water levels cause variations in the exchange flow rates from − 35 l/s to 82 l/s. The highest rates occur during brief infiltration after rapidly rising river water levels. Simulations of different scenarios indicate that riverbank clogging will decrease the exchange fluxes by up to 80%, while clogging of both riverbank and riverbed essentially stops the flow exchange. The groundwater model is also used to simulate the transport of a conservative tracer. The variation of river water levels over time is shown to increase the extent of the active river–aquifer mixing zone in the gravel bar. These dynamic factors significantly enhance the dilution of conservative tracer concentrations in this zone.     

Groundwater–surface water interactions in a large semi‐arid floodplain: implications for salinity management

Flow regulation and water diversion for irrigation have considerably impacted the exchange of surface water between the Murray River and its floodplains. However, the way in which river regulation has impacted groundwater–surface water interactions is not completely understood, especially in regards to the salinization and accompanying vegetation dieback currently occurring in many of the floodplains. Groundwater–surface water interactions were studied over a 2 year period in the riparian area of a large floodplain (Hattah–Kulkyne, Victoria) using a combination of piezometric surface monitoring and environmental tracers (Cl⁻, δ²H, and δ¹⁸O). Despite being located in a local and regional groundwater discharge zone, the Murray River is a losing stream under low flow conditions at Hattah–Kulkyne. The discharge zone for local groundwater, regional groundwater and bank recharge is in the floodplain within ∼1 km of the river and is probably driven by high rates of transpiration by the riparian Eucalyptus camaldulensis woodland. Environmental tracers data suggest that the origin of groundwater is principally bank recharge in the riparian zone and a combination of diffuse rainfall recharge and localized floodwater recharge elsewhere in the floodplain. Although the Murray River was losing under low flows, bank discharge occurred during some flood recession periods. The way in which the water table responded to changes in river level was a function of the type of stream bank present, with point bars providing a better connection to the alluvial aquifer than the more common clay‐lined banks. Understanding the spatial variability in the hydraulic connection with the river channel and in vertical recharge following inundations will be critical to design effective salinity remediation strategies for large semi‐arid floodplains. Copyright © 2005 John Wiley & Sons, Ltd.     

Artificial recharge through a thick, heterogeneous unsaturated zone

Thick, heterogeneous unsaturated zones away from large streams in desert areas have not previously been considered suitable for artificial recharge from ponds. To test the potential for recharge in these settings, 1.3 x 10(6) m(3) of water was infiltrated through a 0.36-ha pond along Oro Grande Wash near Victorville, California, between October 2002 and January 2006. The pond overlies a regional pumping depression 117 m below land surface and is located where thickness and permeability of unsaturated deposits allowed infiltration and saturated alluvial deposits were sufficiently permeable to allow recovery of water. Because large changes in water levels caused by nearby pumping would obscure arrival of water at the water table, downward movement of water was measured using sensors in the unsaturated zone. The downward rate of water movement was initially as high as 6 m/d and decreased with depth to 0.07 m/d; the initial time to reach the water table was 3 years. After the unsaturated zone was wetted, water reached the water table in 1 year. Soluble salts and nitrate moved readily with the infiltrated water, whereas arsenic and chromium were less mobile. Numerical simulations done using the computer program TOUGH2 duplicated the downward rate of water movement, accumulation of water on perched zones, and its arrival at the water table. Assuming 10 x 10(6) m(3) of recharge annually for 20 years, a regional ground water flow model predicted water level rises of 30 m beneath the ponds, and rises exceeding 3 m in most wells serving the nearby urban area.     

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.     

Heterogeneity influences on stream water–groundwater interactions in a gravel-dominated floodplain

ABSTRACT

Floodplains are composed of complex depositional patterns of ancient and recent stream sediments, and research is needed to address the manner in which coarse floodplain materials affect stream–groundwater exchange patterns. Efforts to understand the heterogeneity of aquifers have utilized numerous techniques typically focused on point-scale measurements; however, in highly heterogeneous settings, the ability to model heterogeneity is dependent on the data density and spatial distribution. The objective of this research was to investigate the correlation between broad-scale methodologies for detecting heterogeneity and the observed spatial variability in stream/groundwater interactions of gravel-dominated alluvial floodplains. More specifically, this study examined the correlation between electrical resistivity (ER) and alluvial groundwater patterns during a flood event at a site on Barren Fork Creek, in the Ozark ecoregion of Oklahoma, USA, where chert gravels were common both as streambed and as floodplain material. Water table elevations from groundwater monitoring wells for a flood event on 1–5 May 2009 were compared to ER maps at various elevations. Areas with high ER matched areas with lower water table slope at the same elevation. This research demonstrated that ER approaches were capable of indicating heterogeneity in surface water–groundwater interactions, and that these heterogeneities were present even in an aquifer matrix characterized as highly conductive. Portions of gravel-dominated floodplain vadose zones characterized by high hydraulic conductivity features can result in heterogeneous flow patterns when the vadose zone of alluvial floodplains activates during storm events.

EDITOR D. Koutsoyiannis; ASSOCIATE EDITOR X. Chen     

Stream water infiltration, bank storage, and storage zone changes due to stream-stage fluctuations

During a flood period, stream-stage increases induce infiltration of stream water into an aquifer; subsequent declines in stream stage cause a reverse motion of the infiltrated water. This paper presents the results of the water exchange rate between a stream and aquifer, the storage volume of the infiltrated stream water in the surrounding aquifer (bank storage), and the storage zone. The storage zone is the part of aquifer where groundwater is replaced by stream water during the flood. MODFLOW was used to simulate stream–aquifer interactions and to quantify rates of stream infiltration and return flow. MODPATH was used to trace the pathlines of the infiltrated stream water and to determine the size of the storage zone. Simulations were focused on the analyses of the effects of the stream-stage fluctuation, aquifer properties, the hydraulic conductivity of streambed sediments, regional hydraulic gradients, and recharge and evapotranspiration (ET) rates on stream–aquifer interactions. Generally, for a given stream–aquifer system, larger flow rates result from larger stream-stage fluctuations; larger storage volumes and storage zones are produced by larger and longer-lasting fluctuations. For a given stream-stage hydrograph, a lower-permeable streambed, an aquitard, or an anisotropic aquifer of low vertical hydraulic conductivity can significantly reduce the rate of infiltration and limit the size of the storage zone. The bank storage solely caused by the stage fluctuation differs slightly between gaining and losing streams. Short-term rainfall recharge and ET loss in the shallow groundwater slightly influence on the flow rate, but their effects on bank storage in a larger area for a longer period can be considerable.     

Assessment of the quantitative and qualitative buffer function of an alluvial wetland: hydrological modelling of a large floodplain (Garonne River,France)

For 2 years, water flow‐patterns in the Garonne floodplain of south‐western France were studied in the field and through hydrodynamic modelling (MARTHE Hydrodynamic Software developed by BRGM). Water flow‐paths and the transport of dissolved elements between river and aquifer have been investigated and modelled. In order to quantify the buffer function of the alluvial floodplain, we focused our work on the effect of a major flood on the water flow‐direction, and on nitrate transport. Thus, we showed that the effect of a large flood in the river was rapidly lost with increasing distance from the river. During the observation period, a hydrologically active strip only 300 m wide on either side of the riverbed played a buffering role in absorbing the flood crest. It was also found that meanders favour the exchange between river and alluvial aquifer, shown by the creation of bypasses between the upstream and downstream parts of meanders. This, in turn, contributes to a dilution of nitrates in the phreatic aquifer, which here has higher nitrate content than the surface water; such dilution may result in an overestimation of the denitrification process in the wooded riverbanks. The coupling of chemical measurements—especially of chlorides and nitrate—with modelling of the dissolved‐element transport allows us to establish the water balance for the riparian wetland, and to separate the effect of dilution and denitrification on nitrate concentration. This indicated the existence of areas in the riparian wetlands where denitrification is particularly strong, leading to reductions in nitrate concentrations of 10 to 30 mg/l NO₃^? during the flood. Copyright © 2003 John Wiley & Sons, Ltd.     

Studies of a regulated dryland river: surface–groundwater interactions

Unlike rivers in humid regions, dryland rivers typically exhibit reduced flow in the downstream direction as a result of transmission losses, which include seepage of streamflow into the aquifer, evaporation, and transpiration. However, much remains to be learned about the nature of the exchange between surface water and groundwater in these landscapes, especially in terms of spatial and temporal variability. Our study focused on streambank seepage and groundwater flow in the alluvial aquifer, specifically on answering questions such as: Is there seasonal variability in seepage losses? Is seepage permanently lost? Can losses be reduced by killing riparian vegetation? To better understand the magnitude, variability, and fate of streambank seepage, we assessed river stages, groundwater hydraulic gradients, and groundwater flow paths at two sites along a reach of the Pecos River, a dryland perennial river in West Texas. We found that along this reach the river was losing water to the aquifer even under low‐flow conditions; but seepage was controlled by a number of different mechanisms. Seepage increased not only during high‐flow events but also when the groundwater level was declining owing to long periods of no irrigation release. Tamarix (saltcedar) control did not affect hydraulic gradients nor reduce streambank seepage and given that this reach of the Pecos River is a losing one, streamflow will not be enhanced by controlling saltcedar. These findings can be used to improve basic conceptual models of dryland river systems and to predict hydrologic responses to changes in the timing and magnitude of streamflows and to riparian vegetation management. Copyright © 2012 John Wiley & Sons, Ltd.     

Hydrogeological interactions between fault zones and alluvial aquifers in regional flow systems

The hydrological influence of fault zones in tectonic areas is usually difficult to depict from field data. Numerical simulation allows representation of such flow systems and an estimation of flow lines and rates. This paper reports on simulations of the groundwater flow in a range‐and‐basin area affected by a regional fault zone, which may drain or recharge an overlaying alluvial aquifer. Different hydraulic conductivity values for the range rocks, the fault‐zone, and the sedimentary infill of the basin are considered, as well as different fault‐zone widths and boundary conditions. Results show that upward and downward fluxes develop in the upper part of the fault zone, controlled by the action of the alluvial aquifer, influencing the recharge of the sedimentary basin. This paper shows the hydrological efficiency of fault zones as preferential flow; it also analyses the constraints that determine groundwater recharge to the surrounding basins. These results contribute to the understanding of hydrogeological dynamics in tectonic areas. Copyright © 2008 John Wiley & Sons, Ltd.     

Reducing monitoring gaps at the aquifer–river interface by modelling groundwater–surface water exchange flow patterns

This study investigates spatial patterns and temporal dynamics of aquifer–river exchange flow at a reach of the River Leith, UK. Observations of sub‐channel vertical hydraulic gradients at the field site indicate the dominance of groundwater up‐welling into the river and the absence of groundwater recharge from surface water. However, observed hydraulic heads do not provide information on potential surface water infiltration into the top 0–15 cm of the streambed as these depths are not covered by the existing experimental infrastructure. In order to evaluate whether surface water infiltration is likely to occur outside the ‘window of detection’, i.e. the shallow streambed, a numerical groundwater model is used to simulate hydrological exchanges between the aquifer and the river. Transient simulations of the successfully validated model (Nash and Sutcliff efficiency of 0·91) suggest that surface water infiltration is marginal and that the possibility of significant volumes of surface water infiltrating into non‐monitored shallow streambed sediments can be excluded for the simulation period. Furthermore, the simulation results show that with increasing head differences between river and aquifer towards the end of the simulation period, the impact of streambed topography and hydraulic conductivity on spatial patterns of exchange flow rates decreases. A set of peak flow scenarios with altered groundwater‐surface water head gradients is simulated in order to quantify the potential for surface water infiltration during characteristic winter flow conditions following the observation period. The results indicate that, particularly at the beginning of peak flow conditions, head gradients are likely to cause substantial increase in surface water infiltration into the streambed. The study highlights the potential for the improvement of process understanding of hyporheic exchange flow patterns at the stream reach scale by simulating aquifer‐river exchange fluxes with a standard numerical groundwater model and a simple but robust model structure and parameterization. Copyright © 2011 John Wiley & Sons, Ltd.     

Development and application of an analytical model of stream/aquifer interaction

A mathematical model to simulate stream/aquifer interactions in an unconfined aquifer subjected to time varying river stage was developed from the linearized Boussinesq equation using the principle of superposition and the concept of semigroups. The mathematical model requires an estimate of three parameters to simulate ground-water elevations; transmissivity, specific yield, and recharge. The solution has physical significance and includes terms for the steady-state water level, the steady-state water level as influenced by a change in river stage, a transient redistribution of water levels in the aquifer from the previous day, and a transient change in water level caused by a change in river stage. The mathematical model was tested using observed water table elevations at three locations across a 2-km-wide alluvial valley aquifer. The average absolute deviation between observed and simulated daily water levels was 0.09 m. The difference in river stage over the test year was 4.9 m.     

Geomorphology,hydrology, and aquatic vegetation drive seasonal hyporheic flow patterns across a gravel‐dominated floodplain

Across 1·7 km² of the Umatilla River floodplain (Oregon, USA), we investigated the influences of an ephemeral tributary and perennial ‘spring channel’ (fed only by upwelling groundwater) on hyporheic hydrology. We derived maps of winter and summer water‐table elevations from data collected at 46 monitoring wells and 19 stage gauges and used resulting maps to infer groundwater flow direction. Groundwater flow direction varied seasonally across the floodplain and was influenced by main channel stage, flooding, the tributary creek, and the location and direction of hyporheic exchange in the spring channel. Hyporheic exchange in the spring channel was evaluated with a geochemical mixing model, which confirmed patterns of floodplain groundwater movement inferred from water‐table maps and showed that the spring channel was fed predominantly by hyporheic water from the floodplain aquifer (87% during winter, 80% during summer), with its remaining flow supplied by upslope groundwater from the adjacent catchment aquifer. Summertime growth of aquatic macrophytes in the spring channel also influenced patterns of hyporheic exchange and groundwater flow direction in the alluvial aquifer by increasing flow resistance in the spring channel, locally raising surface water stage and adjacent water‐table elevation, and thereby altering the slope of the water‐table in the hyporheic zone. The Umatilla River floodplain is larger than most sites where hyporheic hydrology has been investigated in detail. Yet, our results corroborate other research that has identified off‐channel geomorphic features as important drivers of hyporheic hydrology, including previously published modeling efforts from a similar river and field observations from smaller streams. Copyright © 2007 John Wiley & Sons, Ltd.     

Aquifer modelling of the Ganga–Mahawa sub‐basin,a part of the Central Ganga Plain,Uttar Pradesh,India

The Ganga–Mahawa sub‐basin, which has an area of 1280 km² forms the western part of the Central Ganga Plain in the Moradabad and Badaun districts of western Uttar Pradesh, India. The Bundelkhand granite forms the basement complex, overlain unconformably by the upper Vindhyan sequence, which is further overlain by the Neogene (Middle and Upper) Siwaliks and finally by Quaternary alluvium. Four geomorphological units, the Varanasi older alluvial plain, Aligarh older alluvial plain, terrace zones and the Ganga recent floodplain, abandoned channels, channel scars and meander scars represent various landforms. The hydrogeological cross‐sections indicate the occurrence of a single aquifer down to 120 m. Some influent seepage from the River Ganga could be seen around Gangeswari, but the rest of the River Ganga is effluent. Groundwater‐flow modelling was carried out to assess the degree of Ganga river and aquifer interaction. The River Ganga marks the western boundary; boundaries to the northeast and southeast are set as fixed heads to simulate lateral inflow into and outflow from the sub‐basin respectively. The eastern boundary is simulated as a no‐flow condition. The Mahawa and Badmar rivers are considered to be effluent. The area modelled is covered by a grid of 34 rows×46 columns with three layers, viz., an unconfined aquifer, an aquitard which is underlain by a semi‐confined to confined aquifer. The permeability distribution was inferred from morphometric analysis and pumping tests. Natural recharge due to monsoon rainfall forms the main input. The River Ganga stage data at Ahar, Naora and Ramghat has been used for assigning surface water levels and river bed elevations in the model. Abstraction from all existing deep and shallow tube wells has been assigned as output at various cells. A steady state flow simulation was carried out and calibrated against the June 1986 water level; subsequent transient conditions were calibrated up to May 1995. The computed groundwater balance was comparable to that estimated from field investigations. The aquifer modelling study has attempted to integrate all available information and provided a tool that could be used for predictive simulation. Copyright © 2000 John Wiley & Sons, Ltd.