Irreversible sorption of benzene in sandy aquifer materials

Hydrocarbon compounds in aquifers are generally known to show a retardation effect due to sorption onto the surfaces of solid particles. In this study, we investigated the effect of sorption on the transport of benzene in sandy aquifer materials by conducting batch and column tests for both sandy aquifer materials and sandy materials to which had been added 0·5% powdered activated carbon. The batch test was conducted by equilibrating dry materials with benzene solutions of various initial concentrations, and by analysing the concentrations of benzene in the initial and equilibrated solutions using high‐performance liquid chromatography (HPLC). The column test was performed to monitor the concentrations of effluent versus time, known as a breakthrough curve (BTC). We injected KCl and benzene solutions as tracers into the inlet boundary as two different types of square pulse and step, and monitored the effluent concentrations at the exit boundary under a steady‐state condition using an electrical conductivity meter and HPLC. Simulation of benzene transport was performed using the convective–dispersive equation model with the distribution coefficients obtained from the batch test and the transport parameters of the conservative solute KCl from the column test. The observed BTCs of KCl and benzene for pulse injection showed that the arrival times of the peaks of both tracers coincided well, but the relative peak concentration of benzene was much lower than that of KCl. Comparison of the simulated and observed BTCs showed a great discrepancy for all cases of injection mode and material texture, indicating the absence of retardation effect. These results reveal that the predominant process affecting the benzene transport in the sandy aquifer materials is an irreversible sorption rather than retardation. This tentative conclusion was verified by simulation of benzene transport using an irreversible sorption parameter that led to a good agreement between the simulated and observed BTCs. Copyright © 2002 John Wiley & Sons, Ltd.     

Quantification of irreversible benzene sorption in sandy materials

Based on a previous study of the irreversible sorption of benzene in sandy aquifer materials, we further investigated a method to quantify an irreversible sorption coefficient of aqueous benzene. Assuming that the rate of irreversible loss from the solution to the sorption sites followed first‐order kinetics, the irreversible sorption coefficient was derived from a kinetic batch sorption test conducted for an appropriate soil‐to‐solution ratio to reflect the flow conditions imposed on a column test. Simulation results revealed that the irreversible sorption coefficient estimated from the kinetic batch test provided a good agreement with the measured data obtained from the column test, indicating that the method proposed in this study can be used to quantify the irreversible sorption coefficient. Copyright © 2004 John Wiley & Sons, Ltd.     

Effect of NAPL Source Morphology on Mass Transfer in the Vadose Zone

The generation of vapor‐phase contaminant plumes within the vadose zone is of interest for contaminated site management. Therefore, it is important to understand vapor sources such as non‐aqueous‐phase liquids (NAPLs) and processes that govern their volatilization. The distribution of NAPL, gas, and water phases within a source zone is expected to influence the rate of volatilization. However, the effect of this distribution morphology on volatilization has not been thoroughly quantified. Because field quantification of NAPL volatilization is often infeasible, a controlled laboratory experiment was conducted in a two‐dimensional tank (28 cm × 15.5 cm × 2.5 cm) with water‐wet sandy media and an emplaced trichloroethylene (TCE) source. The source was emplaced in two configurations to represent morphologies encountered in field settings: (1) NAPL pools directly exposed to the air phase and (2) NAPLs trapped in water‐saturated zones that were occluded from the air phase. Airflow was passed through the tank and effluent concentrations of TCE were quantified. Models were used to analyze results, which indicated that mass transfer from directly exposed NAPL was fast and controlled by advective‐dispersive‐diffusive transport in the gas phase. However, sources occluded by pore water showed strong rate limitations and slower effective mass transfer. This difference is explained by diffusional resistance within the aqueous phase. Results demonstrate that vapor generation rates from a NAPL source will be influenced by the soil water content distribution within the source. The implications of the NAPL morphology on volatilization in the context of a dynamic water table or climate are discussed.     

Biogeochemical evolution of a landfill leachate plume, Norman, Oklahoma

Leachate from municipal landfills can create groundwater contaminant plumes that may last for decades to centuries. The fate of reactive contaminants in leachate-affected aquifers depends on the sustainability of biogeochemical processes affecting contaminant transport. Temporal variations in the configuration of redox zones downgradient from the Norman Landfill were studied for more than a decade. The leachate plume contained elevated concentrations of nonvolatile dissolved organic carbon (NVDOC) (up to 300 mg/L), methane (16 mg/L), ammonium (650 mg/L as N), iron (23 mg/L), chloride (1030 mg/L), and bicarbonate (4270 mg/L). Chemical and isotopic investigations along a 2D plume transect revealed consumption of solid and aqueous electron acceptors in the aquifer, depleting the natural attenuation capacity. Despite the relative recalcitrance of NVDOC to biodegradation, the center of the plume was depleted in sulfate, which reduces the long-term oxidation capacity of the leachate-affected aquifer. Ammonium and methane were attenuated in the aquifer relative to chloride by different processes: ammonium transport was retarded mainly by physical interaction with aquifer solids, whereas the methane plume was truncated largely by oxidation. Studies near plume boundaries revealed temporal variability in constituent concentrations related in part to hydrologic changes at various time scales. The upper boundary of the plume was a particularly active location where redox reactions responded to recharge events and seasonal water-table fluctuations. Accurately describing the biogeochemical processes that affect the transport of contaminants in this landfill-leachate-affected aquifer required understanding the aquifer's geologic and hydrodynamic framework.     

The Role of Bioattenuation in the Management of Aromatic Hydrocarbon Plumes in Aquifers

Ground water scientists have made significant advances in understanding the soil interactions, hydrogeology, fate and transport, and subsurface microbiology of aromatic hydrocarbons (BTEX) in aquifer systems. It is now generally recognized that a major factor responsible for the attenuation and mass reduction of BTEX in plumes is the widespread occurrence of hydrocarbon biodegradation by indigenous soil microorganisms in aquifer material. Most well-studied BTEX plumes that develop from the accidental release of gasoline fuels contain low levels of soluble hydrocarbons (< 1 to 5000 ppb) and have been shown to be spatially confined because of natural biotransformation mechanisms. These in situ processes are controlled by source and aquifer characteristics, permeability, sorption, and geochemical properties of the aquifer. Many laboratory subsoil-ground water microcosms and field studies (10 to 20 C) have demonstrated the rapid biodecay (1 to SO percent/day for microcosms and 0.5 to 1.5 percent/day for plumes) of these aromatic compounds under primarily aerobic conditions (i.e., those with sufficient dissolved oxygen). The ability to implement ground water bioremediation will depend upon our understanding of source control and aquifer recharge effects on the spatial distribution of plumes. In addition, estimating the biodegradation of sorbed BTEX, determining limits and potential for in situ biostimulation of soluble plumes, and establishing data requirements for predictive modeling of natural attenuation will be useful for this remediation technology. The use of these tools to manage ground water quality appears to represent the most practical alternative, particularly for low-risk ground water supplies.     

Analysis of Cryptosporidium parvum oocyst transport in porous media

Cryptosporidium parvum is a protozoan parasite, transmitted through aqueous environments in the form of an oocyst. In this study, a transport model into which sorption, filtration and inactivation mechanisms are incorporated is applied to simulate laboratory column data, and the suitability of a kinetic model to describe the C. parvum oocyst transport and removal in porous media is compared with an equilibrium model. The kinetic model is applied to simulate previous column experimental data and successfully simulates the concentration peak; the late time tailing effect appeared in the breakthrough curves, indicating that the kinetic model is more suitable than the equilibrium one at simulating the fate and transport of the oocysts in porous media. Simulation illustrates that sorption causes retardation along with a tailing in the breakthrough curve. Additionally, filtration acts as a major mechanism of removing the oocysts from the aqueous phase, whereas the role of inactivation in reducing the viable oocyst concentration is minimal. Copyright © 2004 John Wiley & Sons, Ltd.     

A simple method for calculating growth rates of petroleum hydrocarbon plumes

Consumption of aquifer Fe(III) during biodegradation of ground water contaminants may result in expansion of a contaminant plume, changing the outlook for monitored natural attenuation. Data from two research sites contaminated with petroleum hydrocarbons show that toluene and xylenes degrade under methanogenic conditions, but the benzene and ethylbenzene plumes grow as aquifer Fe(III) supplies are depleted. By considering a one-dimensional reaction front in a constant unidirectional flow field, it is possible to derive a simple expression for the growth rate of a benzene plume. The method balances the mass flux of benzene with the Fe(III) content of the aquifer, assuming that the biodegradation reaction is instantaneous. The resulting expression shows that the benzene front migration is retarded relative to the ground water velocity by a factor that depends on the concentrations of hydrocarbon and bioavailable Fe(III). The method provides good agreement with benzene plumes at a crude oil study site in Minnesota and a gasoline site in South Carolina. Compared to the South Carolina site, the Minnesota site has 25% higher benzene flux but eight times the Fe(III), leading to about one-sixth the expansion rate. Although it was developed for benzene, toluene, ethylbenzene, and xylenes, the growth-rate estimation method may have applications to contaminant plumes from other persistent contaminant sources.     

Sorption Behaviour of Phenols on Natural Sandy Aquifer Material during Flow‐through Column Experiments: The Effect of pH

Flow‐through column experiments were carried out to investigate the influence of pH on the sorption of three phenols (2‐methyl‐4, 6‐dinitrophenol, 2, 4, 6‐trichlorophenol, pentachlorophenol) onto a natural sandy aquifer material collected from a bank filtration site of River Elbe, Germany. For the phenols investigated, an increase in sorption (retardation) with decreasing pH is observed indicating a stronger sorption of the neutral species in comparison to that of the anions formed by dissociation. The anions of 2‐methyl‐4, 6‐dinitrophenol and 2, 4, 6‐trichlorophenol do not show significant sorption. On the contrary, pentachlorophenol showed sorption not only in neutral form but also in ionic form significantly which should be taken into account while assessing the fate and transport of such compound. A linear model based on the degree of protonation (calculated from pH and pK_a) can be used to resolve the apparent (observed) sorption coefficient (K_{d, app}) into its neutral (K_{d, n}) and ionised (K_{d, i}) components. Knowing pK_a, K_{d, n}, and K_{d, i} the apparent sorption coefficient for pH values other than experimentally investigated can be predicted.     

The effects of floods on the temperature of riparian groundwater

The transition area between rivers and their adjacent riparian aquifers, which may comprise the hyporheic zone, hosts important biochemical reactions, which control water quality. The rates of these reactions and metabolic processes are temperature dependent. Yet the thermal dynamics of riparian aquifers, especially during flooding and dynamic groundwater flow conditions, has seldom been studied. Thus, we investigated heat transport in riparian aquifers during 3 flood events of different magnitudes at 2 sites along the same river. River and riparian aquifer temperature and water‐level data along the Lower Colorado River in Central Texas, USA, were monitored across 2‐dimensional vertical sections perpendicular to the bank. At the downstream site, preflood temperature penetration distance into the bank suggested that advective heat transport from lateral hyporheic exchange of river water into the riparian aquifer was occurring during relatively steady low‐flow river conditions. Although a small (20‐cm stage increase) dam‐controlled flood pulse had no observable influence on groundwater temperature, larger floods (40‐cm and >3‐m stage increases) caused lateral movement of distinct heat plumes away from the river during flood stage, which then retreated back towards the river after flood recession. These plumes result from advective heat transport caused by flood waters being forced into the riparian aquifer. These flood‐induced temperature responses were controlled by the size of the flood, river water temperature during the flood, and local factors at the study sites, such as topography and local ambient water table configuration. For the intermediate and large floods, the thermal disturbance in the riparian aquifer lasted days after flood waters receded. Large floods therefore have impacts on the temperature regime of riparian aquifers lasting long beyond the flood's timescale. These persistent thermal disturbances may have a significant impact on biochemical reaction rates, nutrient cycling, and ecological niches in the river corridor.     

Chlorinated Decreasing Solvents: Physical-Chemical Properties Affecting Aquifer Contamination and Remediation

Chlorinated degreasing solvents are multicomponent liquids containing not only the chlorinated hydrocarbons with which their name is associated (e.g., trichloroethylene or |TCE]. perchloroethylene or [PCE], 1,1,1-trichlorocihane [TCA]) but also a number of organic additives included as corrosion inhibitors and antioxidants. The additives, such as 1,4-dioxane, are likely to be of significant public-health importance as ground water contaminants due to their toxicity, solubility, and mobility. Following their use in vapor degreasing systems by industry, chlorinated degreasing solvents will also contain about 25% solubilized oil and grease.
A number of physical-chemical properties become especially important in the light of the multicomponent nature of these solvents. First, the higher aqueous solubility and lower sorption of the additives makes it is reasonable to expect that faster moving plumes of these solvent additives will precede plumes of the chlorinated hydrocarbons. Second, due to high losses of chlorinated hydrocarbons by volatilization from vapor degreasers during years in the middle of the century, it is probable that background concentrations of these hydrocarbons are present in ground water flow systems due to their downwind washout. Finally, the solubilized oil and grease may cause profound changes to the wettability of aquifer materials contacted by the solvents during their subsurface migration. It is argued, therefore, that the wettability of aquifer materials contaminated by chlorinated degreasing solvents needs to be experimentally determined before remediation of DNAPL at each site, rather than being simply assumed as water wet.     

Modeling Arsenic Sorption in the Subsurface with a Dual‐Site Model

Arsenic is a well‐known groundwater contaminant that causes toxicological and carcinogenic effects in humans. Predicting the transport of arsenic in the subsurface is often problematic because of its complex sorption characteristics. Numerous researchers have reported that arsenic sorption on soil material is initially fast and then subsequently slow. A dual‐site numerical sorption model was previously developed to describe arsenic desorption from arsenic‐contaminated soils in batch experiments in terms of two different release mechanisms. Experiments involving synthetic acid rain leaching of four arsenic‐contaminated soil columns were performed to verify the dual‐site numerical sorption model in the context of one‐dimensional vertical transport. The fitted models successfully simulated the signature long tailings and the two‐stage arsenic leaching patterns for all four soil columns. The dual‐site sorption model was incorporated within the general solute transport simulation code Modular Three‐Dimensional Multispecies (MT3DMS), version 5.10. The resulting version was named MT3DDS and is available for public access. This experimental study has shown that MT3DDS is capable of simulating phase redistribution during transport, and thus provides a new numerical tool for simulating arsenic transport in the subsurface.     

Determination of two‐dimensional laboratory‐scale dispersivities

A tracer test was conducted in a laboratory chamber representing a two‐dimensional aquifer to investigate the longitudinal dispersivity (α_L) and the ratio (α_T/α_L) of transverse to longitudinal dispersivity of sandy aquifer materials. Dispersive parameters were obtained by matching the observed chloride plumes at 9 hours and 16 hours after tracer injection with those simulated by a flow and transport model. The best match was found for α_L = 0·2 ? 0·25 cm and α_T/α_L = 0·2. The ratio of α_T/α_L = 0·2 was within the range of laboratory values reported in the literature. Sensitivity analysis revealed that the tracer plume concentration and shape were more sensitive to variations in longitudinal dispersivity than to the ratio of transverse to longitudinal dispersivity. This result contrasted with findings of others, showing that the dispersivity ratio greatly affects contaminant plume shape. However, our experimental boundary conditions restricted expansion of the plume normal to the direction of flow and thus affected the parameter estimation. Copyright © 2004 John Wiley & Sons, Ltd.     

A Coupled Groundwater–Surface Water Modeling Framework for Simulating Transition Zone Processes

There is an identified need for fully representing groundwater–surface water transition zone (i.e., the sediment zone that connects groundwater and surface water) processes in modeling fate and transport of contaminants to assist with management of contaminated sediments. Most existing groundwater and surface water fate and transport models are not dynamically linked and do not consider transition zone processes such as bioturbation and deposition and erosion of sediments. An interface module is developed herein to holistically simulate the fate and transport by coupling two commonly used models, Environmental Fluid Dynamics Code (EFDC) and SEAWAT, to simulate surface water and groundwater hydrodynamics, while providing an enhanced representation of the processes in the transition zone. Transition zone and surface water contaminant processes were represented through an enhanced version of the EFDC model, AQFATE. AQFATE also includes SEDZLJ, a state‐of‐the‐science surface water sediment transport model. The modeling framework was tested on a published test problem and applied to evaluate field‐scale two‐ and three‐dimensional contaminant transport. The model accurately simulated concentrations of salinity from a published test case. For the field‐scale applications, the model showed excellent mass balance closure for the transition zone and provided accurate simulations of all transition zone processes represented in the modeling framework. The model predictions for the two‐dimensional field case were consistent with site‐specific observations of contaminant migration. This modeling framework represents advancement in the simulation of transition zone processes and can help inform risk assessment at sites where contaminant sources from upland areas have the potential to impact sediments and surface water.     

Numerical investigation of road salt impact on an urban wellfield

The impact of road salt on a wellfield in a complex glacial moraine aquifer system is studied by numerical simulation. The moraine underlies an extensive urban and industrial landscape, which draws its water supply from >20 wellfields, several of which are approaching or have exceeded the drinking water limit for chloride. The study investigates the mechanisms of road salt infiltration, storage, and transport in the subsurface and assesses the effectiveness of mitigation measures designed to reduce the impact. The three-dimensional transport model accounts for increases in salt loading, as well as growth of the urbanized area and road network over the past 50 years. The simulations, which focus on one impacted wellfield, show chloride plumes originating mainly at arterial roads and migrating through aquitard windows into the water supply aquifers. The results suggest that the aquifer system contains a large and heterogeneously distributed mass of chloride and that concentrations in the aquifer can be substantially higher than the concentrations in the well water. Future impact scenarios indicate that although the system responds rapidly to reductions in salt loading, the residual chloride mass may take decades to flush out, even if road salting were discontinued. The implications with respect to urban wellfields in typical snow-belt areas are discussed.     

Review of Quantitative Surveys of the Length and Stability of MTBE,TBA, and Benzene Plumes in Groundwater at UST Sites

Quantitative information regarding the length and stability condition of groundwater plumes of benzene, methyl tert‐butyl ether (MTBE), and tert‐butyl alcohol (TBA) has been compiled from thousands of underground storage tank (UST) sites in the United States where gasoline fuel releases have occurred. This paper presents a review and summary of 13 published scientific surveys, of which 10 address benzene and/or MTBE plumes only, and 3 address benzene, MTBE, and TBA plumes. These data show the observed lengths of benzene and MTBE plumes to be relatively consistent among various regions and hydrogeologic settings, with median lengths at a delineation limit of 10 µg/L falling into relatively narrow ranges from 101 to 185 feet for benzene and 110 to 178 feet for MTBE. The observed statistical distributions of MTBE and benzene plumes show the two plume types to be of comparable lengths, with 90th percentile MTBE plume lengths moderately exceeding benzene plume lengths by 16% at a 10‐µg/L delineation limit (400 feet vs. 345 feet) and 25% at a 5‐µg/L delineation limit (530 feet vs. 425 feet). Stability analyses for benzene and MTBE plumes found 94 and 93% of these plumes, respectively, to be in a nonexpanding condition, and over 91% of individual monitoring wells to exhibit nonincreasing concentration trends. Three published studies addressing TBA found TBA plumes to be of comparable length to MTBE and benzene plumes, with 86% of wells in one study showing nonincreasing concentration trends.     

Investigation of Small‐Scale Preferential Flow with a Forced‐Gradient Tracer Test

A new tracer experiment (referred to as MADE‐5) was conducted at the well‐known Macrodispersion Experiment (MADE) site to investigate the influence of small‐scale mass‐transfer and dispersion processes on well‐to‐well transport. The test was performed under dipole forced‐gradient flow conditions and concentrations were monitored in an extraction well and in two multilevel sampler (MLS) wells located at 6, 1.5, and 3.75 m from the source, respectively. The shape of the breakthrough curve (BTC) measured at the extraction well is strongly asymmetric showing a rapidly arriving peak and an extensive late‐time tail. The BTCs measured at seven different depths in the two MLSs are radically different from one another in terms of shape, arrival times, and magnitude of the concentration peaks. All of these characteristics indicate the presence of a complex network of preferential flow pathways controlling solute transport at the test site. Field‐experimental data were also used to evaluate two transport models: a stochastic advection‐dispersion model (ADM) based on conditional multivariate Gaussian realizations of the hydraulic conductivity field and a dual‐domain single‐rate (DDSR) mass‐transfer model based on a deterministic reconstruction of the aquifer heterogeneity. Unlike the stochastic ADM realizations, the DDSR accurately predicted the magnitude of the concentration peak and its arrival time (within a 1.5% error). For the multilevel BTCs between the injection and extraction wells, neither model reproduced the observed values, indicating that a high‐resolution characterization of the aquifer heterogeneity at the subdecimeter scale would be needed to fully capture 3D transport details.     

Effects of rate-limited mass transfer on water sampling with partially penetrating wells

Nonequilibrium concentration type curves are numerically developed and sensitivity analyses are performed to examine the relationships between effluent concentrations in partially penetrating monitoring/extraction wells, the vertical plume shape, and the mass transfer characteristics of the aquifer. The governing two-dimensional, axisymmetric nonequilibrium solute transport equation is solved in three stages using an operator-splitting approach. In the first two stages, the advection and dispersion terms are solved with the Eulerian-Lagrangian method, based on the backward method of characteristics for advection and the standard implicit Galerkin finite element method for dispersion. In the third step, the first-order, immobile-mobile domain mass transfer term is computed analytically for both two-site and lognormally distributed, multirate models. Effluent concentration variations with time and contour plots of the pore water concentration distribution in the aquifer are compared for a wide range of field- and laboratory-measured mass transfer rates, various plume shapes, and relevant physical/chemical parameter values, including pumping rate, vertical anisotropy ratio, retardation factor, and porosity. The simulation results show that rate-limited mass transfer can have a significant impact on sample and aquifer pore water concentrations during three-dimensional transport to a partially penetrating well. An alternative dimensionless form of the nonequilibrium solute transport equation is derived to illustrate the key parameter groupings that quantify rate-limited sorption effects and show the relative importance of individual parameters. A hypothetical field application example demonstrates the fitting of dimensional type curves to discrete-interval sampling data in order to evaluate the mass transfer characteristics of an aquifer and shows how type curve superposition can be used to model complex plume shapes.     

Modeling the effects of nonlinear equilibrium sorption on the transport of solute plumes in saturated heterogeneous porous media

Transport of sorbing solutes in 2D steady and heterogeneous flow fields is modeled using a particle tracking random walk technique. The solute is injected as an instantaneous pulse over a finite area. Cases of linear and Freundlich sorption isotherms are considered. Local pore velocity and mechanical dispersion are used to describe the solute transport mechanisms at the local scale. This paper addresses the impact of the degree of heterogeneity and correlation lengths of the log-hydraulic conductivity field as well as negative correlation between the log-hydraulic conductivity field and the log-sorption affinity field on the behavior of the plume of a sorbing chemical. Behavior of the plume is quantified in terms of longitudinal spatial moments: center-of-mass displacement, variance, 95% range, and skewness. The range appears to be a better measure of the spread in the plumes with Freundlich sorption because of plume asymmetry. It has been found that the range varied linearly with the travelled distance, regardless of the sorption isotherm. This linear relationship is important for extrapolation of results to predict behavior beyond simulated times and distances. It was observed that the flow domain heterogeneity slightly enhanced the spreading of nonlinearly sorbing solutes in comparison to that which occurred for the homogeneous flow domain, whereas the spreading enhancement in the case of linear sorption was much more pronounced. In the case of Freundlich sorption, this enhancement led to further deceleration of the solute plume movement as a result of increased retardation coefficients produced by smaller concentrations. It was also observed that, except for plumes with linear sorption, correlation between the hydraulic conductivity and the sorption affinity fields had minimal effect on the spatial moments of solute plumes with nonlinear sorption.     

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.     

Understanding the Geometry of Connected Fracture Flow with Multiperiod Oscillatory Hydraulic Tests

An understanding of the spatial and hydraulic properties of fast preferential flow pathways in the subsurface is necessary in applications ranging from contaminant fate and transport modeling to design of energy extraction systems. One method for the characterization of fracture properties over interwellbore scales is Multiperiod Oscillatory Hydraulic (MOH) testing, in which the aquifer response to oscillatory pressure stimulations is observed. MOH tests were conducted on isolated intervals of wells in siliciclastic and carbonate aquifers in southern Wisconsin. The goal was to characterize the spatial properties of discrete fractures over interwellbore scales. MOH tests were conducted on two discrete fractured intervals intersecting two boreholes at one field site, and a nest of three piezometers at another field site. Fracture diffusivity estimates were obtained using analytical solutions that relate diffusivity to observed phase lag and amplitude decay. In addition, MOH tests were used to investigate the spatial extent of flow using different conceptual models of fracture geometry. Results indicated that fracture geometry at both field sites can be approximated by permeable two‐dimensional fracture planes, oriented near‐horizontally at one site, and near‐vertically at the other. The technique used on MOH field data to characterize fracture geometry shows promise in revealing fracture network characteristics important to groundwater flow and transport.