Energy considerations in groundwater‐ridging mechanism of streamflow generation

Groundwater ridging is the rapid rise of a shallow water table during a rainfall event, in an environment where, in the pre‐event period, the capillary fringe extends to the ground surface. Groundwater ridging is widely cited to account for the observed significant appearance of pre‐event water in a stream stormflow hydrograph. Various hypotheses have been advanced to explain the groundwater‐ridging mechanism; and most recently, from a field study site in South Africa, an energy hypothesis was proposed, which explains that groundwater‐ridging water‐table rise is a result of rapid introduction and transmission of additional pressure head into the capillary fringe from an intense rainfall at the ground surface. However, there is a need for further analysis and evidence from other field study sites to confirm and support this newly proposed energy hypothesis. The objectives of this paper are, therefore, as follows: to review previous observations on groundwater ridging, from other study sites, in order to deduce evidence of the newly proposed energy hypothesis; to present and evaluate a one‐dimensional diffusion mathematical model that can simulate groundwater‐ridging water‐table rise, based on the newly proposed energy hypothesis; and to evaluate the importance of a capillary fringe in streamflow generation. Analysis of previous observations from other study sites generally indicated that the rate of groundwater‐ridging water‐table rise is directly related to the rainfall intensity, hence confirming and agreeing with the newly proposed energy hypothesis. Additionally, theoretical results by the mathematical model agreed fairly well with the field results observed under natural rainfall, confirming that the rapidly rainfall‐induced energy is diffusively transmitted downwards through pore water, elevating the pressure head at every depth. The results in this study also support the concept of a three‐end‐member stream stormflow hydrograph and contribute to the explanation of how catchments can store water for long periods but then release it rapidly during storm events. Copyright © 2015 John Wiley & Sons, Ltd.     

Soil frost effects on soil water and runoff dynamics along a boreal forest transect: 1. Field investigations

To determine how soil frost changes flowpaths of runoff water along a hillslope, a transect consisting of four soil profiles directed towards a small stream in a mature forest stand was investigated at Svartberget, near Vindeln in northern Sweden. Soil temperature, unfrozen water content, groundwater level and snow depth were investigated along the transect, which started at the riparian peat, and extended 30 m upslope into mineral soils. The two, more organic‐rich profiles closest to the stream had higher water retention and wetter autumn conditions than the sandy mineral soils further upslope. The organic content of the soil influenced the variation in frost along the transect. The first winter (1995–96) had abnormally low snow precipitation, which gave a deep frost down to 40–80 cm, whereas the two following winters had frost depths of 5–20 cm. During winter 1995–96, the two organic profiles close to the stream had a shallower frost depth than the mineral soil profile higher upslope, but a considerably larger amount of frozen water. The fraction of water that did not freeze despite several minus degrees in the soil was 5–7 vol.% in the mineral soil and 10–15 vol.% in the organic soil. From the measurements there were no signs of perched water tables during any of the three snowmelt periods, which would have been strong evidence for changed water flowpaths due to soil frost. When shallow soil layers became saturated during snowmelt, especially in 1997 and 1998, it was because of rising groundwater levels. Several rain on frozen ground events during spring 1996 resulted in little runoff, since most of the rain either froze in the soil or filled up the soil water storage. Copyright © 2001 John Wiley & Sons, Ltd.     

Defining hydrochemical evolution of streamflow through flowpath dynamics in Kawakami headwater catchment,Central Japan

The hydrochemical behaviour of catchments is often investigated by inferring stream chemistry through identification of source areas involved in hydrograph separation analysis, yet its dynamic evolution of hydrologic pathways has received little attention. Intensive hydrometric and hydrochemical measurements were performed during two different storms on March 29, 2001 and August 21–22, 2001 to define hydrochemical evolution under the dynamic of flow pathways in a 5·2 ha first‐order drainage of the Kawakami experimental basin (KEB), Central Japan, a forested headwater catchment with various soil depths (1·8 to 5 m) overlying late Neogene of volcanic bedrocks. The hydraulic potential distribution and flow lines data showed that the change in flow direction, which was controlled by rainfall amount and antecedent wetness of the soil profile, agreed well with the hydrochemical change across the slope segment during the storm. Hydrograph separation predicted by end‐member mixing analysis (EMMA) using Ca²⁺ and SiO₂ showed that near surface riparian, hillslope soil water and deep riparian groundwater were important in stream flow generation. The evidence of decrease in solutes concentration at a depth of 1 m in the hillslope and 0·6 m in the near surface riparian during peak storm suggested a flushing of high solutes concentration. Most of the solutes accumulated in the deep riparian groundwater zone, which was due to prominent downward flow and agreed well with the residence time. The distinct flow pathways and chemistry between the near surface riparian and deep riparian groundwater zones and the linkage hillslope aquifer and near surface riparian reservoir, which controls rapid flow and solutes flushing during the storm event, are in conflict with the typical assumption that the whole riparian zone resets flow pathways and chemical signature of hillslope soil water, as has been reported in a previous study. Copyright © 2005 John Wiley & Sons, Ltd.     

A low‐dimensional model of bedrock weathering and lateral flow coevolution in hillslopes: 2. Controls on weathering and permeability profiles,drainage hydraulics,and solute export pathways

The advance of a chemical weathering front into the bedrock of a hillslope is often limited by the rate weathering products that can be carried away, maintaining chemical disequilibrium. If the weathering front is within the saturated zone, groundwater flow downslope may affect the rate of transport and weathering—however, weathering also modifies the rock permeability and the subsurface potential gradient that drives lateral groundwater flow. This feedback may help explain why there tends to be neither “runaway weathering” to great depth nor exposed bedrock covering much of the earth and may provide a mechanism for weathering front advance to keep pace with incision of adjacent streams into bedrock. This is the second of a two‐part paper exploring the coevolution of bedrock weathering and lateral flow in hillslopes using a simple low‐dimensional model based on hydraulic groundwater theory. Here, we show how a simplified kinetic model of 1‐D rock weathering can be extended to consider lateral flow in a 2‐D hillslope. Exact and approximate analytical solutions for the location and thickness of weathering within the hillslope are obtained for a number of cases. A location for the weathering front can be found such that lateral flow is able to export weathering products at the rate required to keep pace with stream incision at steady state. Three pathways of solute export are identified: “diffusing up,” where solutes diffuse up and away from the weathering front into the laterally flowing aquifer; “draining down,” where solutes are advected primarily downward into the unweathered bedrock; and “draining along,” where solutes travel laterally within the weathering zone. For each pathway, a different subsurface topography and overall relief of unweathered bedrock within the hillslope is needed to remove solutes at steady state. The relief each pathway requires depends on the rate of stream incision raised to a different power, such that at a given incision rate, one pathway requires minimal relief and, therefore, likely determines the steady‐state hillslope profile.     

Towards a better understanding of riparian zone water table response to precipitation: surface water infiltration,hillslope contribution or pressure wave processes?

Determining how riparian zone hydrological conditions may change in response to precipitation in various geomorphic settings is critical to determine the occurrence of hot moments of biogeochemical transformations for phosphorus, nitrogen, sulfate, mercury and greenhouse gases in these systems. The author investigate water table response to precipitation at a high temporal resolution (15 min) in a riparian zone located in a deeply incised glacial till valley (20 m) with approximately 2 m of alluvium over a confining layer, in Indiana, USA. During storms, larger water table fluctuations (approximately 100 cm) occurred near the stream than near the toe slope (10–25 cm). A quick rise in water table near the stream occurred for all storms, with partial flow reversals occurring for three of seven storms. The quick rise of the water table near the stream was associated with a decrease in hillslope water contributions to the stream during storms and the development of a water table down valley gradient for most storms. Water table fluctuations, groundwater flow velocities and electrical conductivity data indicated that riparian zone water table response to precipitation was primarily regulated by pressure wave processes. Regardless of the storm, high water tables persisted for at least 2 days after the cessation of precipitation. Although this suggests that high‐resolution precipitation data may be useful to quantify hot moments of biogeochemical transformation associated with high water tables in riparian zones, precipitation data alone are not sufficient to correctly estimate the magnitude of riparian water table level changes during storms. Copyright © 2011 John Wiley & Sons, Ltd.     

Groundwater fluctuations and footslope ferricrete soils in the humid tropical zone of southern Cameroon

This paper discusses the relationship between the differentiation of ferruginous accumulations and the variable water saturation of footslope soil patterns. An analysis of the slope morphology of a typical hill in the forest zone of southern Cameroon and a seasonal survey of the levels of groundwaters, springs and rivers were considered in relation to the petrology of different soil patterns. The study site is a tabular hillock whose slopes present a progressive development from steep to gentle slopes. The variable residence time of water within the soil, creating an alternation of reducing and oxidizing conditions, affects soil chemistry, structure and lateral extension of the soil patterns. The ferruginous soil patterns, being formed on the footslopes, gradually increase in extent with decreasing slope angle and the relative rise of the groundwater level. The steep footslopes, where groundwater has a shorter residence time, show a soft mottled clay pattern, restricted to the bottom part of the slope. The moderate footslopes exhibit a deep permanent and a temporary perched groundwater table. The latter, with its regular capillary fringe, contributes to more reducing conditions within isolated domains in the soil patterns, and thus to the alternation with oxidizing conditions, generating a continuous hard soil pattern (massive carapace). The more gently dipping footslopes exhibit groundwater levels near the surface and also a significant amplitude of groundwater fluctuation. Iron, previously accumulated in moderate footslope patterns, is reduced, remobilized, and leached. The soil patterns formed develop into a variegated carapace, more extended along the slope, containing less iron, but nevertheless more hardened, due to the important fluctuations of the groundwater table. These patterns are limited to the zone of groundwater fluctuation and deteriorate as the water fluctuation zone recedes. Copyright © 2005 John Wiley & Sons, Ltd.     

Throughflow variability during snowmelt in a forested mountain catchment,coastal British Columbia,Canada

This study documented the spatial and temporal variability of outflow from a forested hillslope segment during snowmelt at a small mountain catchment in south coastal British Columbia, Canada. A pit 5 m wide was established just upslope from the stream channel. Outflow from the organic horizon was intercepted and measured by a single trough, and outflow from the mineral horizons was measured separately for three adjacent sections. Throughflow exhibited non‐steady‐state behaviour involving shifting allocations of flow amongst different sections of the outflow pit, as well as threshold effects and hysteresis in the relationship between pit outflow and water table elevation. Most of the pit outflow originated from the mineral horizons, indicating that throughflow was the dominant pathway by which water was delivered to the stream channel. Direct precipitation and snowmelt onto near‐stream saturated areas can account for less than 20% of the total outflow from the hillslope segment. Throughflow from the mineral sections consistently peaked either at the same time as or earlier than stream flow from the 150‐ha catchment during diurnal snowmelt cycles, indicating that throughflow appears to respond rapidly enough to contribute to snowmelt‐induced peak stream flow. Pit outflow cannot be extrapolated reliably to the catchment scale on the basis of simple length‐ or area‐based ratios. Copyright © 2004 John Wiley & Sons, Ltd.     

Tidal effects on groundwater dynamics in unconfined aquifers

The variation of seawater level resulting from tidal fluctuations is usually neglected in regional groundwater flow studies. Although the tidal oscillation is damped near the shoreline, there is a quasi‐steady‐state rise in the mean water‐table position, which may have an influence on regional groundwater flow. In this paper the effects of tidal fluctuations on groundwater hydraulics are investigated using a variably saturated numerical model that includes the effects of a realistic mild beach slope, seepage face and the unsaturated zone. In particular the impact of these factors on the velocity field in the aquifer is assessed. Simulations show that the tidal fluctuation has substantial consequences for the local velocity field in the vicinity of the exit face, which affects the nearshore migration of contaminant in coastal aquifers. An overheight in the water table as a result of the tidal fluctuation is observed and this has a significant effect on groundwater discharge to the sea when the landward boundary condition is a constant water level. The effect of beach slope is very significant and simplifying the problem by considering a vertical beach face causes serious errors in predicting the water‐table position and the groundwater flux. For media with a high effective capillary fringe, the moisture retained above the water table is important in determining the effects of the tidal fluctuations. Copyright © 2001 John Wiley & Sons, Ltd.     

Effects of Airflow Induced by Rainfall on Shallow Groundwater Table Fluctuations

An investigation of groundwater table fluctuations induced by rainfall should consider interactions between the liquid and gas phases in soils. In this study, a water‐air two‐phase flow model was initially verified by simulating an infiltration experiment. It was then employed to model the interactions between liquid and gas phases regarding actions of airflow on the groundwater table and the fluctuations of the phreatic level and water level in the well induced by rainfall. The effects of airflo7w caused by rainfall on phreatic level fluctuations were also studied quantitatively by comparing the results obtained using the proposed model with those obtained from a water single‐phase flow model. The simulation results show that in addition to actual recharge, compressed airflow in unsaturated zones causes the phreatic level to increase, but the rise in the phreatic level is lower than that in the pore‐air pressure head in unsaturated zones due to the mitigation of capillary fringe. The existence of airflow enhances the phreatic level rise during and after rainfall. In addition, the water level in the well, pushed by the phreatic level fluctuations, varies similarly to the phreatic level, but it experiences somewhat delayed and slightly attenuated. The Lisse effect precisely reflects the phreatic level fluctuations before actual recharge. Furthermore, the fluctuations in the phreatic level and water level in the well and the contributions of airflow to phreatic level fluctuations are affected by many factors: rain intensity, initial moisture, overlying aquitard, groundwater table depths, and screen depths of the well.     

A numerical investigation of bedrock groundwater recharge and exfiltration on soil mantled hillslopes

Here we use Richards Equation models of variably saturated soil and bedrock groundwater flow to investigate first-order patterns of the coupling between soil and bedrock flow systems. We utilize a Monte Carlo sensitivity analysis to identify important hillslope parameters controlling bedrock recharge and then model the transient response of bedrock and soil flow to seasonal precipitation. Our results suggest that hillslopes can be divided into three conceptual zones of groundwater interaction, (a) the zone of lateral unsaturated soil moisture accumulation (upper portion of hillslope), (b) the zone of soil saturation and bedrock recharge (middle of hillslope) and (c) the zone of saturated-soil lateral flow and bedrock groundwater exfiltration (bottom of hillslope). Zones of groundwater interaction expand upslope during periods of precipitation and drain downslope during dry periods. The amount of water partitioned to the bedrock groundwater system a can be predicted by the ratio of bedrock to soil saturated hydraulic conductivity across a variety of hillslope configurations. Our modelled processes are qualitatively consistent with observations of shallow subsurface saturation and groundwater fluctuation on hillslopes studied in our two experimental watersheds and support a conceptual model of tightly coupled shallow and deep subsurface circulation where groundwater recharge and discharge continuously stores and releases water from longer residence time storage.     

Groundwater dynamics in a till hillslope: flow directions,gradients and delay

Knowledge of groundwater dynamics is important for the understanding of hydrological controls on chemical processes along the water flow pathways. To increase our knowledge of groundwater dynamics in areas with shallow groundwater, the groundwater dynamics along a hillslope were studied in a boreal catchment in Southern Sweden. The forested hillslope had a 1‐ to 2‐m deep layer of sandy till above bedrock. The groundwater flow direction and slope were calculated under the assumption that the flow followed the slope of the groundwater table, which was computed for different triangles, each defined by three groundwater wells. The flow direction showed considerable variations over time, with a maximum variation of 75°. During periods of high groundwater levels the flow was almost perpendicular to the stream, but as the groundwater level fell, the flow direction became gradually more parallel to the stream, directed in the downstream direction. These findings are of importance for the interpretation of results from hillslope transects, where the flow direction usually is assumed to be invariable and always in the direction of the hillslope. The variations in the groundwater flow direction may also cause an apparent dispersion for groundwater‐based transport. In contrast to findings in several other studies, the groundwater level was most responsive to rainfall and snowmelt in the upper part of the hillslope, while the lower parts of the slope reached their highest groundwater level up to 40 h after the upper parts. This can be explained by the topography with a wetter hollow area in the upper part. Copyright © 2011 John Wiley & Sons, Ltd.     

Impact of Heterogeneity on Oxygen Transfer in a Fluctuating Capillary Fringe

We performed quasi‐two‐dimensional flow through laboratory experiments to study the effect of a coarse‐material inclusion, located in the proximity of the water table, on flow and oxygen transfer in the capillary fringe. The experiments investigate different phases of mass transfer from the unsaturated zone to anoxic groundwater under both steady‐state and transient flow conditions, the latter obtained by fluctuating the water table. Monitoring of flow and transport in the different experimental phases was performed by visual inspection of the complex flow field using a dye tracer solution, measurement of oxygen profiles across the capillary fringe, and determination of oxygen fluxes in the effluent of the flow‐through chamber. Our results show significant effects of the coarse‐material inclusion on oxygen transfer during the different phases of the experiments. At steady state, the oxygen flux across the unsaturated/saturated interface was considerably enhanced due to flow focusing in the fully water‐saturated coarse‐material inclusion. During drainage, a zone of higher water saturation formed in the fine material overlying the coarse lens. The entrapped oxygen‐rich aqueous phase contributed to the total amount of oxygen supplied to the system when the water table was raised back to its initial level. In case of imbibition, pronounced air entrapment occurred in the coarse lens, causing oxygen to partition between the aqueous and gaseous phases. The oxygen mass supplied to the anoxic groundwater following the imbibition event was found to be remarkably higher (approximately seven times) in the heterogeneous system compared with a similar experiment performed in a homogeneous porous medium.     

Flow directions of stream-groundwater exchange in a headwater catchment during the hydrologic year

Understanding near-stream groundwater dynamics and flow directions is important for predicting hillslope-stream connectivity, streamflow generation, and hydrologic controls of streamwater quality. To determine the drivers of groundwater flow in the stream corridor (i.e., the stream channel and the adjacent groundwater in footslopes and riparian areas), we observed the water levels of 36 wells and 7 piezometers along a headwater stream section over a period of 18 months. Groundwater dynamics during events were controlled by the initial position of the groundwater table relative to the subsurface structure. The near-stream groundwater table displayed a fast and pronounced response to precipitation events when lying in fractured bedrock with low storage capacity, and responded less frequently and in a less pronounced way when lying in upper layers with high storage capacity. Precipitation depth, intensity, regolith thickness above the fractured bedrock, and proximity to and elevation above the stream channel also had an effect on the groundwater dynamics, which varied with hydrologic conditions. Our high-frequency and spatially dense measurements highlight the competing influence of groundwater inflow from upslope locations, streamwater level and bedrock properties on the spatiotemporal dynamics of flowpaths in the stream corridor. Near-stream groundwater pointed uniformly towards the stream channel when the stream corridor was hydrologically connected to upslope groundwater. However, local interruptions of the water inflow from upslope locations caused flow reversals towards the footslopes. The direction of near-stream groundwater followed the local fractured bedrock topography during dry hydrologic conditions on a few occasions after events. The outcomes of this research contribute to a better understanding of the drivers controlling spatiotemporal changes in near-stream groundwater dynamics and flow directions in multiple wetness states of the stream corridor.     

Insights on stream temperature processes through development of a coupled hydrologic and stream temperature model for forested coastal headwater catchments

Stream temperature controls a number of biological, chemical, and physical processes occurring in aquatic environments. Transient snow cover and advection associated with lateral throughflow inputs can have a dominant influence on stream thermal regimes for headwater catchments in the rain‐on‐snow zone. Most existing stream temperature models lack the ability to properly simulate these processes. We developed and evaluated a conceptual‐parametric catchment‐scale stream temperature model that includes the role of transient snow cover and lateral advection associated with throughflow. The model consists of routines for simulating canopy interception, snow accumulation and melt, hillslope throughflow runoff and temperature, and stream channel energy exchange processes. The model was used to predict discharge and stream temperature for a small forested headwater catchment near Vancouver, Canada, using long‐term (1963–2013) weather data to compute model forcing variables. The model was evaluated against 4 years of observed stream temperature. The model generally predicted daily mean stream temperature accurately (annual RMSE between 0.57 and 1.24 °C) although it overpredicted daily summer stream temperatures by up to 3 °C during extended low streamflow conditions. Model development and testing provided insights on the roles of advection associated with lateral throughflow, channel interception of snow, and surface–subsurface water interactions on stream thermal regimes. This study shows that a relatively simple but process‐based model can provide reasonable stream temperature predictions for forested headwater catchments located in the rain‐on‐snow zone.     

Hillslope hydrology and wetland response of two small zero‐order boreal catchments on the Precambrian Shield

Two Precambrian Shield zero‐order catchments were monitored from January 2003 to July 2004 to characterize their hydrological and biogeochemical characteristics prior to a forest management experiment. Hydrometric observations were used to examine temporal trends in hillslope‐wetland connectivity and the hillslope runoff processes that control wetland event response. The hillslope groundwater flux from the longer transect (E1) was continuous throughout the study period. Groundwater fluxes from a shorter and steeper hillslope (E0) were intermittent during the study period. Large depression storage elements (termed micro‐basins) located on the upper hillslope of the E1 catchment appeared to be at least partly responsible for the observed rapid wetland runoff responses. These micro‐basins were hydrologically connected to a downslope wetland by a subsurface channel of glacial cobbles that functioned as a macropore channel during episodic runoff events. The runoff response from the hilltop micro‐basins is controlled by antecedent water table position and water is quickly piped to the wetland fringe through the cobble channel during high water table conditions. During periods of low water table position, seepage along the bedrock–soil interface from the hilltop micro‐basin and other hillslopes maintained hillslope–wetland connectivity. The micro‐basins create a dynamic variable source‐area runoff system where the contributing area expands downslope during episodic runoff events. The micro‐basins occupied 30% of the E1 catchment and are a common feature on the Precambrian Shield. Copyright © 2007 John Wiley & Sons, Ltd.     

Modelling the effects of soil type and root distribution on shallow groundwater resources

Vegetated, shallow groundwater environments typically have high environmental and economic value. A sound understanding of the complex interactions and feedbacks between surface vegetation and groundwater resources is crucial to managing and maintaining healthy ecosystems while responding to human needs. A vegetated shallow groundwater environment was modelled using the software HYDRUS 2D to investigate the effects of several combinations of soil type and root distributions on shallow groundwater resources. Three rainfall regimes coupled to both natural and anthropogenically affected groundwater conditions were used to investigate the effect that combinations of four soil types and five root distributions can have on (a) groundwater level drops, (b) groundwater depletion, (c) groundwater recharge and (d) water stress conditions. Vegetation with roots distributed across the whole unsaturated zone and vegetation with dimorphic root systems (i.e. roots having larger concentrations both near the surface and the capillary fringe) behaved differently from vegetation growing roots mainly near the saturated zone. Specifically, vegetation with roots in the unsaturated zone caused water‐table drops and groundwater depletions that were half the amount due to deep‐rooted vegetation. Vegetation with a large portion of roots near the soil surface benefited from rainfall and was less vulnerable to water‐table lowering; as such, the fraction of the total area of roots affected by water stress conditions could be 40% smaller than in the case with deep‐rooted vegetation. However, roots uniformly distributed in the unsaturated zone could halve groundwater recharge rates observed in bare soils. Our analysis provided insights that can enable the formulation of site‐ and purpose‐specific management plans to respond to both human and ecosystem water requirements. Copyright © 2015 John Wiley & Sons, Ltd.     

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.     

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.     

Experimental Sensitivity Analysis of Oxygen Transfer in the Capillary Fringe

Oxygen transfer in the capillary fringe (CF) is of primary importance for a wide variety of biogeochemical processes occurring in shallow groundwater systems. In case of a fluctuating groundwater table two distinct mechanisms of oxygen transfer within the capillary zone can be identified: vertical predominantly diffusive mass flux of oxygen, and mass transfer between entrapped gas and groundwater. In this study, we perform a systematic experimental sensitivity analysis in order to assess the influence of different parameters on oxygen transfer from entrapped air within the CF to underlying anoxic groundwater. We carry out quasi two‐dimensional flow‐through experiments focusing on the transient phase following imbibition to investigate the influence of the horizontal flow velocity, the average grain diameter of the porous medium, as well as the magnitude and the speed of the water table rise. We present a numerical flow and transport model that quantitatively represents the main mechanisms governing oxygen transfer. Assuming local equilibrium between the aqueous and the gaseous phase, the partitioning process from entrapped air can be satisfactorily simulated. The different experiments are monitored by measuring vertical oxygen concentration profiles at high spatial resolution with a noninvasive optode technique as well as by determining oxygen fluxes at the outlet of the flow‐through chamber. The results show that all parameters investigated have a significant effect and determine different amounts of oxygen transferred to the oxygen‐depleted groundwater. Particularly relevant are the magnitude of the water table rise and the grain size of the porous medium.