Migration of Petroleum Products Through Sandy Hydrogeologic Systems

Laboratory column experiments were carried out to study the migration of petroleum products through a sandy porous medium. It was found that the oil pressure needed to displace water from the pores of the sand medium increased with depth below the top of the column. While oil under negative pressure displaced water through the variably unsaturated zone, a significant vertical column of oil was needed to displace water from the pores at the water table. These results indicate that oil penetration to and below the water table will occur only if the porous medium is highly conductive and the rate of oil leak is high. For small to moderate leak rates and hydraulic conductivities, oil would preferentially spread laterally through the zones above the water table. This process of spreading could serve as a natural barrier to severe aquifer contamination by petroleum products.
A simplified procedure was developed to give an order-of-magnitude estimate of the preferred subsurface migration pathways of leaked petroleum products. This procedure utilizes the main drainage and wetting curves for oil and water and the interfacial pressure between these two fluids. This approach could be useful in guiding exploratory investigations, reducing both the risk for further spreading of the contaminants and the investigative cost.     

Estimation of groundwater vulnerability to pollution based on DRASTIC in the Niipele sub-basin of the Cuvelai Etosha Basin,Namibia

Surface water is a scarce resource in Namibia with about sixty percent of Namibia's population dependent on groundwater for drinking purposes. With increasing population, the country faces water challenges and thus groundwater resources need to be managed properly. One important aspect of Integrated Water Resources Management is the protection of water resources, including protection of groundwater from contamination and over-exploitation. This study explores vulnerability mapping as a basic tool for protecting groundwater resources from pollution. It estimates groundwater vulnerability to pollution in the upper Niipele sub-basin of the Cuvelai-Etosha in Northern Namibia using the DRASTIC index. The DRASTIC index uses GIS to estimate groundwater vulnerability by overlaying different spatially referenced hydrogeological parameters that affect groundwater contamination. The study assesses the discontinuous perched aquifer (KDP) and the Ohangwena multi-layered aquifer 1 (KOH-1). For perched aquifers, point data was regionalized by a hydrotope approach whereas for KOH-1 aquifer, inverse distance weighting was used. The hydrotope approach categorized different parts of the hydrogeological system with similar properties into five hydrotopes. The result suggests that the discontinuous perched aquifers are more vulnerable than Ohangwena multi-layered aquifer 1. This implies that vulnerability increases with decreasing depth to water table because contaminants have short travel time to reach the aquifer when they are introduced on land surface. The nitrate concentration ranges between 2 and 288 mg/l in perched aquifers while in Ohangwena multi-layered aquifer 1, it ranges between 1 and 133 mg/l. It was observed that perched aquifers have high nitrate concentrations than Ohangwena 1 aquifer, which correlates well with the vulnerability results.     

Bioremediation in Fractured Rock: 1. Modeling to Inform Design,Monitoring, and Expectations

Field characterization of a trichloroethene (TCE) source area in fractured mudstones produced a detailed understanding of the geology, contaminant distribution in fractures and the rock matrix, and hydraulic and transport properties. Groundwater flow and chemical transport modeling that synthesized the field characterization information proved critical for designing bioremediation of the source area. The planned bioremediation involved injecting emulsified vegetable oil and bacteria to enhance the naturally occurring biodegradation of TCE. The flow and transport modeling showed that injection will spread amendments widely over a zone of lower‐permeability fractures, with long residence times expected because of small velocities after injection and sorption of emulsified vegetable oil onto solids. Amendments transported out of this zone will be diluted by groundwater flux from other areas, limiting bioremediation effectiveness downgradient. At nearby pumping wells, further dilution is expected to make bioremediation effects undetectable in the pumped water. The results emphasize that in fracture‐dominated flow regimes, the extent of injected amendments cannot be conceptualized using simple homogeneous models of groundwater flow commonly adopted to design injections in unconsolidated porous media (e.g., radial diverging or dipole flow regimes). Instead, it is important to synthesize site characterization information using a groundwater flow model that includes discrete features representing high‐ and low‐permeability fractures. This type of model accounts for the highly heterogeneous hydraulic conductivity and groundwater fluxes in fractured‐rock aquifers, and facilitates designing injection strategies that target specific volumes of the aquifer and maximize the distribution of amendments over these volumes.     

The pollution of groundwater and unsaturated zone rocks by oil products and PAH in urban areas

An approach to assessing the hazard of unsaturated zone rocks as a source of subsoil-water secondary pollution by oil products and PAH is considered based on environmental-geochemical studies of soils and rocks of the unsaturated zone, groundwater, and snow, which were carried out during the rehabilitation of the water environment of Lefortovskii Park (Moscow). The impact of the primary pollution sources is analyzed, and the formation of groundwater and rock pollution by oil products and PAH is considered. To assess the role of the unsaturated zone as a secondary source of pollution of surface and subsurface waters under the conditions of possible elimination of the primary sources, the transport of oil products (benzol and toluene) and PAH (pyrene) in unsaturated zone was simulated using WHI UnSat Suite Plus software package.     

An approach to catchment‐scale groundwater nitrate risk assessment from diffuse agricultural sources: a case study in the Upper Bann,Northern Ireland

DRASTIC has drawbacks in groundwater risk assessment that are important in guiding activities to prevention agricultural diffuse groundwater pollution. This paper presents an improved and GIS‐based D‐DRASTIC approach for groundwater nitrate risk assessment from diffuse agricultural sources based on DRASTIC. D‐DRASTIC considers the risk concept, nitrate loading, pollutant transport with runoff, depth to water, net recharge, aquifer media, soil media, topography, impact of the vadose zone media, and the hydraulic conductivity of the aquifer. D‐DRASTIC was developed within an ArcGIS environment and applied to the Upper Bann Catchment, Northern Ireland as a case study. D‐DRASTIC shows that ‘very high’ and ‘high’ zones of groundwater nitrate risk occupy 5% and 11% of the case study area, respectively. When considering groundwater pollution sources and pathways, the results using D‐DRASTIC are helpful in guiding the activities of groundwater pollution prevention at the catchment scale in the context of better implementation of the EU Water Framework Directive. Copyright © 2008 John Wiley & Sons, Ltd.     

Impact of the Use of Reclaimed Water on the Quality of Groundwater Resources in the Jordan Valley,Jordan

The use of reclaimed water and its impact on groundwater quality in the middle and southern parts of the Jordan Valley are investigated. The chemical analyses indicate that nitrate and bacteriological pollution is widespread, and thus, seriously affects groundwater use. During the study, 365 water samples were collected from wells and springs to determine the water chemistry and the extent of nitrate pollution. Three hydrochemical facies are identifed, i. e., (Ca–(Mg)–Na–HCO₃), (Ca–Na–SO₄–Cl) and (Ca–Na–Cl). The change of facies is accompanied by a gradual increase in the groundwater total dissolved solids (TDS), which is mainly controlled by evaporates and carbonates dissolution in the aquifer matrix. Water analyses indicate that the shallow aquifer in the study area is affected by non‐point pollution sources, primarily from natural (manure) and chemical nitrogen (N)‐fertilizers and treated wastewater used for agriculture. The concentration of nitrate in the groundwater ranges from 10 to 355 mg/L. Considerable seasonal fluctuations in groundwater quality are observed as a consequence of agricultural practices and other factors such as annual rainfall distribution and the Zarqa River flow. The noticeable levels of total coliform and Escherichia coli in the northern part of the study area may be attributed to contamination from the urban areas, intensive livestock production, and illegal dumping of sewage. Heavy metal concentrations in all samples were found to be significantly lower than the permissible limits for drinking water standards.     

Investigation of Possible Contamination of Shallow Ground Water by Deeply Injected Liquid Industrial Wastes

The objective of this study was to assess the possible impact of deep well disposal operations, conducted between 1958 and 1974, on the ground water quality in a shallow fresh water aquifer beneath Sarnia, Ontario, Canada. Because of the breakout of formation fluids in Sarnia and Port Huron, Michigan, in the early 1970s, it had been hypothesized that liquid waste from the disposal zone in bedrock had leaked through numerous abandoned oil, gas, and salt wells in the area up to the shallow fresh water aquifer and from there to the surface.
A monitoring well network of 29 5cm (2 inch) diameter piezometers was established in the thin sand and shale aquifer system, which exists between 30 and 70m (100 and 230 feet) below ground surface. In addition, a 300m (1000 foot) deep borehole was drilled and instrumented with a Westbay multilevel casing, which permitted sampling of the disposal zone.
Ground water samples from the shallow monitoring wells and the Westbay multilevel casing were analyzed for volatiles by GC/MS. Those volatile aromatics that were conspicuously present in the deep disposal zone, e.g., ethyl toluenes and trimethyl benzene, were not detected in the shallow monitoring wells. Thus, if contaminants from the disposal zone did indeed migrate to the shallow aquifer, contamination was not widespread and probably consisted mostly of displaced chloride-rich formation waters.     

Assessing the Transport of Pharmaceutical Compounds in a Layered Aquifer Discharging to a Stream

A groundwater plume containing high concentrations of pharmaceutical compounds, mainly sulfonamides, barbiturates, and ethyl urethane, in addition to chlorinated ethenes and benzene was investigated. The contamination originating from a former pharmaceutical industry discharges into a multilayered aquifer system and a downgradient stream. In this study, geological and hydrogeological data were integrated into a numerical flow model to examine identified trends using statistical approaches, including principal component analysis and hierarchal cluster analysis. A joint interpretation of the groundwater flow paths and contaminant concentrations in the different compartments (i.e., groundwater and hyporheic zone) provided insight on the transport processes of the different contaminant plumes to the stream. The analysis of historical groundwater concentrations of pharmaceutical compounds at the site suggested these compounds are slowly degrading. The pharmaceutical compounds migrate in both a deep semiconfined aquifer, as well as in the shallow unconfined aquifer, and enter the stream along a 2-km stretch. This contrasted with the chlorinated ethenes, which mainly discharge to the stream as a focused plume from the unconfined aquifer. The integrated approach developed here, combining groundwater flow modeling and statistical analyses of the contaminant concentration data collected in groundwater and the hyporheic zone, lead to an improved understanding of the observed distribution of contaminants in the unconfined and semiconfined aquifers, and thus to their discharge to the stream. This approach is particularly relevant for large and long-lasting contaminant sources and plumes, such as abandoned landfills and industrial production sites, where field investigations may be very expensive.     

Modeling and Remediation of Ground Water Contamination at the Engelse Werk Wellfield

Contaminants have been threatening the Engelse Werk wellfield located between the town of Zwolle and the IJssel River in the Netherlands. Chemical analysis of water samples taken in production wells, both at the IJssel River and near the Zwolle railway station, indicated elevated concentrations of mainly organic contaminants including benzene, bentazon, acenaftene, trichloroethane, and bromacil. Immediate contaminant prevention and remediation measures are needed to safeguard the production wells. Ground water flow and transport models were developed to assist in the design of remediation strategies. Ground water flow models indicated that the IJssel River and a waste disposal ditch at the railway station are within the capture zone of the wellfield. A chloride transport model simulated minimum travel times in the order of four to 13 years for contaminants in the IJssel River to reach the production wells of the wellfield. A transport model for benzene was set up to advise on the remediation measures to be taken at the waste disposal ditch to clean up the contamination in the upper aquifer between this site and the Engelse Werk wellfield. The designed remediation system consists of 12 pumping wells with a combined capacity of 1650 m³/day. The system is capable of reducing the benzene levels at the threatened production wells at the Engelse Werk wellfield to a permissible level below 0.1 μg/L within a period of 5 years.     

Lagged rejuvenation of groundwater indicates internal flow structures and hydrological connectivity

Large proportions of rainwater and snowmelt infiltrate into the subsurface before contributing to stream flow and stream water quality. Subsurface flow dynamics steer the transport and transformation of contaminants, carbon, weathering products and other biogeochemistry. The distribution of groundwater ages with depth is a key feature of these flow dynamics. Predicting these ages are a strong test of hypotheses about subsurface structures and time-varying processes. Chlorofluorocarbon (CFC)-based groundwater ages revealed an unexpected groundwater age stratification in a 0.47 km² forested catchment called Svartberget in northern Sweden. An overall groundwater age stratification, representative for the Svartberget site, was derived by measuring CFCs from nine different wells with depths of 2–18 m close to the stream network. Immediately below the water table, CFC-based groundwater ages of already 30 years that increased with depth were found. Using complementary groundwater flow models, we could reproduce the observed groundwater age stratification and show that the 30 year lag in rejuvenation comes from return flow of groundwater at a subsurface discharge zone that evolves along the interface between two soil types. By comparing the observed groundwater age stratification with a simple analytical approximation, we show that the observed lag in rejuvenation can be a powerful indicator of the extent and structure of the subsurface discharge zone, while the vertical gradient of the age-depth-relationship can still be used as a proxy of the overall aquifer recharge even when sampled in the discharge zone. The single age stratification profile measured in the discharge zone, close to the aquifer outlet, can reveal the main structure of the groundwater flow pattern from recharge to discharge. This groundwater flow pattern provides information on the participation of groundwater in the hydrological cycle and indicates the lower boundary of hydrological connectivity.     

Understanding the dynamics and contamination of an urban aquifer system using groundwater age (14C, 3H,CFCs) and chemistry

The quality of the groundwater supplying drinking water to the Guadalajara metropolitan area has deteriorated due to both endogenic and exogenic processes. Previous studies of this complex neotectonic volcanic environment suggest that the sources of contamination here are underground fluids derived from an active volcanic center and surface wastewater derived from regional land‐use intensification. This study uses isotopic, gaseous, and chemical signatures to more comprehensively characterize this groundwater flow and its contamination paths. Groundwater is mainly recharged at the La Primavera Caldera to the west and is discharged into the Santiago River to the east. The exception to this trend is the Toluquilla area, where groundwater most likely represents rainfall originating from outside the basin limits. Evaporation affects groundwater in these areas, especially waters that have been affected by recycling below urban areas in the Atejamac area and by intensive agricultural activity in the Toluquilla area. Additionally, we present evidence that groundwater flow through alluvial sediments and tuffs in deeper wells mixes with a lower aquifer unit in basaltic‐andesitic rocks, which are in contact with hydrothermal fluids. Groundwater ages range from postbomb in the western and northwestern regions of the study area (i.e., the Atemajac aquifer unit) to Late Pleistocene in the southern and southeastern regions (i.e., the Toluquilla aquifer unit). Recently recharged water records little mixing and is located mostly in or near the La Primavera volcanic system. As groundwater undergoes gravitational flow towards discharge areas, it mixes with older water components. Chloride and sodium concentrations above natural background levels are primarily related to volcanic activity, nitrate is associated with human activities, and sulfate originates from both anthropogenic sources and water–rock interactions. Nitrate originating from land‐use activities (such as sewers, septic tanks, landfills, and agricultural fields) that is introduced into the deeper part of the groundwater system is expected to travel with the groundwater to the discharge areas because oxidizing conditions will prevent microbial reduction. See Supplementary Information.     

Performance Assessments of a Novel Well Design for Reducing Exposure to Bedrock‐Derived Arsenic

Arsenic in groundwater is a serious problem in New England, particularly for domestic well owners drawing water from bedrock aquifers. The overlying glacial aquifer generally has waters with low arsenic concentrations but is less used because of frequent loss of well water during dry periods and the vulnerability to surface‐sourced bacterial contamination. An alternative, novel design for shallow wells in glacial aquifers is intended to draw water primarily from unconsolidated glacial deposits, while being resistant to drought conditions and surface contamination. Its use could greatly reduce exposure to arsenic through drinking water for domestic use. Hypothetical numerical models were used to investigate the potential hydraulic performance of the new well design in reducing arsenic exposure. The aquifer system was divided into two parts, an upper section representing the glacial sediments and a lower section representing the bedrock. The location of the well, recharge conditions, and hydraulic properties were systematically varied in a series of simulations and the potential for arsenic contamination was quantified by analyzing groundwater flow paths to the well. The greatest risk of arsenic contamination occurred when the hydraulic conductivity of the bedrock aquifer was high, or where there was upward flow from the bedrock aquifer because of the position of the well in the flow system.     

Surface water-groundwater exchange dynamics in buried-valley aquifer systems

Groundwater is a primary source of drinking water worldwide, but excess nutrients and emerging contaminants could compromise groundwater quality and limit its usage as a drinking water source. As such contaminants become increasingly prevalent in the biosphere, a fundamental understanding of their fate and transport in groundwater systems is necessary to implement successful remediation strategies. The dynamics of surface water-groundwater (hyporheic) exchange within a glacial, buried-valley aquifer system are examined in the context of their implications for the transport of nutrients and contaminants in riparian sediments. High conductivity facies act as preferential flow pathways which enhance nutrient and contaminant delivery, especially during storm events, but transport throughout the aquifer also depends on subsurface sedimentary architecture (e.g. interbedded high and low conductivity facies). Temperature and specific conductance measurements indicate extensive hyporheic mixing close to the river channel, but surface water influence was also observed far from the stream-aquifer interface. Measurements of river stage and hydraulic head indicate that significant flows during storms (i.e., hot moments) alter groundwater flow patterns, even between consecutive storm events, as riverbed conductivity and, more importantly, the hydraulic connectivity between the river and aquifer change. Given the similar mass transport characteristics among buried-valley aquifers, these findings are likely representative of glacial aquifer systems worldwide. Our results suggest that water resources management decisions based on average (base) flow conditions may inaccurately represent the system being evaluated, and could reduce the effectiveness of remediation strategies for nutrients and emerging contaminants.     

Ground Water Contamination in Kuwait Resulting from the 1991 Gulf War: A Preliminary Assessment

A large number of oil wells in Kuwait were damaged and ignited by the retreating Iraqi troops during the 1991 Gulf War. The resulting spillage of huge volumes of crude oil on the surface gave rise to oil lakes and crude oil–impregnated soil. Moreover, products of crude oil combustion had spread over a large tract of the ground surface, causing widespread contamination of soil. Hydrocarbon contamination of ground water by the infiltrating runoff water carrying the contaminants from the surface soil to the water table and/or through direct contact with the crude oil leaking through the damaged casing in the subsurface was feared. This preliminary study was carried out to investigate the extent and nature of this possible contamination of ground water. The results indicate that the shallow fresh water lenses present under the Umm Al-Aish water field and in the southeastern parts of the Raudhatain water field in North Kuwait were affected by hydrocarbon pollution. Standard methods like the determination of the contents of the total petroleum hydrocarbon by the Fourier transform infrared method and 16 polyaromatic hydrocarbons using the gas chromatography-mass spectrometry method did not work well, possibly due to the environmental degradation of the crude oil over time. The fluorescence methods and the total organic carbon and total organic matter gave better indications of the intensity and the extent of ground water pollution. The brackish water fields of South and Central Kuwait were, however, free of any indications of hydrocarbon contamination.     

Geochemistry of the Springfield Plateau aquifer of the Ozark Plateaus Province in Arkansas,Kansas, Missouri and Oklahoma,USA

Geochemical data indicate that the Springfield Plateau aquifer, a carbonate aquifer of the Ozark Plateaus Province in central USA, has two distinct hydrochemical zones. Within each hydrochemical zone, water from springs is geochemically and isotopically different than water from wells. Geochemical data indicate that spring water generally interacts less with the surrounding rock and has a shorter residence time, probably as a result of flowing along discrete fractures and solution openings, than water from wells. Water type throughout most of the aquifer was calcium bicarbonate, indicating that carbonate‐rock dissolution is the primary geochemical process occurring in the aquifer. Concentrations of calcium, bicarbonate, dissolved oxygen and tritium indicate that most ground water in the aquifer recharged rapidly and is relatively young (less than 40 years). In general, field‐measured properties, concentrations of many chemical constituents, and calcite saturation indices were greater in samples from the northern part of the aquifer (hydrochemical zone A) than in samples from the southern part of the aquifer (hydrochemical zone B). Factors affecting differences in the geochemical composition of ground water between the two zones are difficult to identify, but could be related to differences in chert content and possibly primary porosity, solubility of the limestone, and amount and type of cementation between zone A than in zone B. In addition, specific conductance, pH, alkalinity, concentrations of many chemical constituents and calcite saturation indices were greater in samples from wells than in samples from springs in each hydrochemical zone. In contrast, concentrations of dissolved oxygen, nitrite plus nitrate, and chloride generally were greater in samples from springs than in samples from wells. Water from springs generally flows rapidly through large conduits with minimum water–rock interactions. Water from wells flow through small fractures, which restrict flow and increase water–rock interactions. As a result, springs tend to be more susceptible to surface contamination than wells. The results of this study have important implications for the geochemical and hydrogeological processes of similar carbonate aquifers in other geographical locations. Copyright © 2000 John Wiley & Sons, Ltd.     

Limitations of Monitoring Wells for the Detection and Quantification of Petroleum Products in Soils and Aquifers

Theoretical analysis and laboratory column experiments were carried out to investigate the conditions required for petroleum products (oil) to flow into a well installed through a sandy porous medium contaminated with the oil. The results indicated that oil would flow into a well only after a layer of "free oil" is formed in the adjacent porous medium. Because significant quantities of oil could be stored in the porous medium under the influence of capillary suction prior to the formation of the zone of free oil, the presence of oil in a well would indicate an advanced stage of oil contamination of the subsurface. While monitoring wells could be used to delineate the extent of the free-oil plume and the plume of dissolved petroleum constituents, they are not useful for delineating the extent of capillary held oil.
The experimental results also indicated that the ratio of the oil-layer thickness in the well to that in the porous medium is not a constant as is sometimes assumed in practice. Further, estimates of the oil thickness in the medium based on the oil thickness in wells and on capillary properties measured in the laboratory were sensitive to the values of the parameters used in these estimates. The measured thickness of the oil layer in a monitoring well alone may not yield reliable estimates of the amount of oil in the subsurface, and assuming that the oil-thickness ratio is a constant can lead to inadequate site assessments and inappropriate remedial plans.     

GIS based DRASTIC model for groundwater vulnerability assessment: Case study of the shallow mio-plio-quaternary aquifer (Southeastern Tunisia)

Groundwater is the main source of water in arid regions. Thus, groundwater pollution becomes a major issue due to the increasing contamination, which poses serious and harmful risk to the environment. Groundwater vulnerability maps can be used as a tool to help decision makers to protect groundwater resources from contamination. The vulnerability of the Mio-Plio-Quaternary shallow aquifer (Southeast Tunisia) has been assessed using a DRASTIC model based on Geographic Information System (GIS). The different parameters of the model were collected from several sources and converted into thematic maps using ArcGis©. Each DRASTIC parameter was assigned a weight and rating based on a range of information within the parameter. Groundwater vulnerability map shows a large area (48%) with high risk of pollution. It indicates that the Southern part of the aquifer and the wadi beds are the most susceptible to contamination. The measured nitrate concentration is coherent with the DRASTIC model results.     

222Rn as Natural Tracer for LNAPL Recovery in a Crude Oil‐Contaminated Aquifer

The objective of this study was to investigate whether ²²²Rn in groundwater can be used as a tracer for light non‐aqueous phase liquid (LNAPL) quantification at a field site treated by dual‐phase LNAPL removal. After the break of a pipeline, 5 ha of soil in the nature reserve Coussouls de Crau in southern France was contaminated by 5100 m³ of crude oil. Part of this oil seeped into the underlying gravel aquifer and formed a floating oil body of about 3.9 ha. The remediation consists of plume management by hydraulic groundwater barriers and LNAPL extraction in the source zone. ²²²Rn measurements were performed in 21 wells in and outside the source zone during 15 months. In uncontaminated groundwater, the radon activity was relatively constant and remained always >11 Bq/L. The variability of radon activity measurements in wells affected by the pump‐and‐skim system was consistent with the measurements in wells that were not impacted by the system. The mean activities in wells in the source zone were, in general, significantly lower than in wells upgradient of the source zone, owing to partitioning of ²²²Rn into the oil phase. The lowest activities were found in zones with high non‐aqueous phase liquid (NAPL) recovery. LNAPL saturations around each recovery well were furthermore calculated during a period of high groundwater level, using a laboratory‐determined crude oil–water partitioning coefficient of 38.5 ± 2.9. This yielded an estimated volume of residual crude oil of 309 ± 93 m³ below the capillary fringe. We find that ²²²Rn is a useful and cheap groundwater tracer for finding zones of good LNAPL recovery in an aquifer treated by dual‐phase LNAPL removal, but that quantification of NAPL saturation using Rn is highly uncertain.     

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.