A controlled field experiment on groundwater contamination by a multicomponent DNAPL: creation of the emplaced-source and overview of dissolved plume development

A unique field experiment has been undertaken at the CFB Borden research site to investigate the development of dissolved chlorinated solvent plumes from a residual dense non-aqueous phase liquid (DNAPL) source. The "emplaced-source" tracer test methodology involved a controlled emplacement of a block-shaped source of sand containing chlorinated solvents below the water table. The gradual dissolution of this residual DNAPL solvent source under natural aquifer conditions caused dissolved solvent plumes of trichloromethane (TCM), trichloroethene (TCE) and perchloroethene (PCE) to continuously develop down gradient. Source dissolution and 3-D plume development were successfully monitored via 173 multilevel samplers over a 475-day tracer test period prior to site remediation research being initiated. Detailed groundwater level and hydraulic conductivity data were collected. Development of plumes with concentrations spanning 1-700,000 micrograms/1 is described and key processes controlling their migration identified. Plumes were observed to be narrow due to the weakness of transverse dispersion processes and long due to advection and significant longitudinal dispersion, very limited sorptive retardation and negligible, if any, attenuation due to biodegradation or abiotic reaction. TCM was shown to be essentially conservative, TCE very nearly conservative and PCE, consistent with its greater hydrophobicity, more retarded yet having a greater mobility than observed in previous Borden field tests. The absence of biodegradation was ascribed to the prevailing aerobic conditions and lack of any additional biodegradable carbon substrates. The transient groundwater flow regime caused significant transverse lateral plume movement, plume asymmetry and was likely responsible for most of the, albeit limited, transverse horizontal plume spreading. In agreement with the widespread incidence of extensive TCE and PCE plumes throughout the industrialized world, the experiment indicates such solvent plumes are likely to be highly mobile and persistent, at least in aquifers that are aerobic and have low sorption potential (low foc content).     

Effects of sorption on transverse mixing in transient flows

Transverse mixing has been identified as a potentially limiting factor for natural attenuation of plumes originating from continuously emitting sources. Under steady-state flow conditions, dispersion is the only process leading to lateral mixing. This process is very slow and cannot explain the lateral spread of plumes observed in the field. When the flow direction fluctuates with time, transverse dispersion is slightly enhanced, but not very pronounced. Under these flow conditions, however, sorption can contribute to mixing into the mean transverse direction. If the reacting compounds differ in their strength of sorption, chromatographic mixing and separation alternate in time-periodic flows. For instantaneous sorption, the plumes may overlap within a stripe of fixed width. In contrast to sorption in local equilibrium, kinetic sorption contributes to mixing also for compounds with identical sorption strength. I derive an analytical expression for the equivalent transverse dispersion coefficient of a kinetically sorbing compound in a spatially uniform flow field undergoing sinusoidal fluctuations in time. This expression may be used for reactive transport calculations in an equivalent domain with constant flow. The effects are the strongest for compounds with a dimensionless partitioning coefficient of about unity, slow sorption kinetics, and slowly fluctuating velocities. For realistic parameters, kinetic sorption contributes to transverse mixing in the same range as heterogeneity.     

Transverse dispersion of non-reactive tracers in porous media: a new nonlinear relationship to predict dispersion coefficients

Assessing the potential of natural attenuation in groundwater relies on the ability to predict and quantify the processes that occur in contaminant plumes. Transverse dispersion is a significant mass transfer mechanism for mixing of electron acceptors and donors and thus may control the lengths of steady state plumes. Laboratory experiments were carried out using a 2-dimensional acrylic glass tank filled with glass beads, quartz sand and field site material as porous media. Flow velocities and grain sizes were varied in order to cover a large range of Peclet numbers including typical field scenarios. The laboratory study was extended by a comprehensive literature search to compare the new results with earlier work. As a result we propose a new empirical relationship for prediction of transverse dispersion coefficients (Dt) which is based on the Peclet number (Pe). This new relationship indicates a nonlinear dependency on the flow velocity (nu a) and grain size (d), namely a relative decrease of the dispersion coefficient with increasing flow velocity in relatively fast flowing water: Dt/Daq=Dp/Daq+0.28(Pe)0.72 (with Pe=nu a d/Daq; Daq and Dp denote the aqueous and pore diffusion coefficients, resp.).     

Numerical experiments and field results on the size of steady state plumes

Contaminated groundwater poses a serious risk for drinking water supplies. Under certain conditions, however, groundwater contamination remains restricted to a tolerable extent because of natural attenuation processes. We present an innovative approach to evaluate the size of these so-called steady-state plumes by 2-D and 1-D modelling in homogeneous aquifers. If longitudinal mixing is negligible, scenarios can be modelled in a simplified way using a 1-D domain vertical to the direction of flow. We analysed the sensitivity of the plume length with respect to biodegradation kinetics, flow velocity, transverse vertical dispersivity alphat, the source and aquifer geometry and reaction stoichiometry. Our findings indicate that for many readily biodegradable compounds transverse-dispersive mixing rather than reaction kinetics is the limiting factor for natural attenuation. Therefore, if alphat, aquifer and source geometry and concentrations of electron acceptors and donors are known, the length of the steady state contaminant plume can be predicted. The approach is validated under field conditions for an ammonium plume at a former landfill site in SW Germany.     

Choice of dispersion coefficients in reactive transport calculations on smoothed fields

Local dispersion dominates the mixing of compounds that are introduced separately into the subsurface and do not partition into any other than the aqueous phase. Thus, reactions between these compounds are controlled by dispersive mixing if they are not limited by kinetics. I quantify longitudinal dispersive mixing by the longitudinal effective dispersion coefficient of a conservative tracer introduced by a point-like injection [Water Resour. Res. 36 (12) (2000) 3591-3604]. In the upscaling of mixing-controlled reactive transport, I apply the mean velocity and the effective dispersion coefficient to the macroscopic transport calculations, whereas the reactive parameters on the macro-scale are identical to those on the local scale. The applicability of the approach is demonstrated for the transport of compounds undergoing a second-order irreversible bimolecular reaction. Ten realizations of a two-dimensional heterogeneous log-conductivity field are considered. Using the effective dispersion parameters, the overall mass balance is met in the ensemble average, whereas solute spreading is underestimated. I assess the lack of spreading by the difference between the expected macrodispersion and effective dispersion coefficients. I extend the approach to simulations on log-conductivity fields obtained by kriging of regularly spaced conductivity measurements. These fields contain the large-scale features of the true fields but do not resolve the small-scale variability. For the calculations on the kriged fields, the corresponding conditional covariance is substituted into the analytical expressions of effective dispersion, yielding a correction effective dispersion coefficient. The comparison between simulations on the fully resolved fields and on the kriged fields indicates that the approach is valid for wide plumes meeting the ergodicity condition. The high variability of mixing on small scales unresolved by kriging, however, leads to severe uncertainty when mixing-controlled reactions are predicted for narrow plumes.     

Processes controlling the distribution and natural attenuation of dissolved phenolic compounds in a deep sandstone aquifer

Processes controlling the distribution and natural attenuation (NA) of phenol, cresols and xylenols released from a former coal-tar distillation plant in a deep Triassic sandstone aquifer are evaluated from vertical profiles along the plume centerline at 130 and 350 m from the site. Up to four groups of contaminants (phenols, mineral acids, NaOH, NaCl) form discrete and overlapping plumes in the aquifer. Their distribution reflects changing source history with releases of contaminants from different locations. Organic contaminant distribution in the aquifer is determined more by site source history than degradation. Contaminant degradation at total organic carbon (TOC) concentrations up to 6500 mg l(-1) (7500 mg l(-1) total phenolics) is occurring by aerobic respiration NO3-reduction, Mn(IV)-/Fe(III)-reduction, SO4-reduction, methanogenesis and fermentation, with the accumulation of inorganic carbon, organic metabolites (4-hydroxybenzaldehyde, 4-hydroxybenzoic acid), acetate, Mn(II), Fe(II), S(-II), CH4 and H2 in the plume. Aerobic and NO3-reducing processes are restricted to a 2-m-thick plume fringe but Mn(IV)-/Fe(II)-reduction, SO4-reduction, methanogenesis and fermentation occur concomitantly in the plume. Dissolved H2 concentrations in the plume vary from 0.7 to 110 nM and acetate concentrations reach 200 mg l(-1). The occurrence of a mixed redox system and concomitant terminal electron accepting processes (TEAPs) could be explained with a partial equilibrium model based on the potential in situ free energy (deltaGr) yield for oxidation of H2 by specific TEAPs. Respiratory processes rather than fermentation are rate limiting in determining the distribution of H2 and TEAPs and H2 dynamics in this system. Most (min. 90%) contaminant degradation has occurred by aerobic and NO3-reducing processes at the plume fringe. This potential is determined by the supply of aqueous O2 and NO3 from uncontaminated groundwater, as controlled by transverse mixing, which is limited in this aquifer by low dispersion. Consumption to date of mineral oxides and SO4 is, respectively, <0.15% and 0.4% of the available aquifer capacity, and degradation using these oxidants is <10%. Fermentation is a significant process in contaminant turnover, accounting for 21% of degradation products present in the plume, and indicating that microbial respiration rates are slow in comparison with fermentation. Under present conditions, the potential for degradation in the plume is very low due to inhibitory effects of the contaminant matrix. Degradation products correspond to <22% mass loss over the life of the plume, providing a first-order plume scale half-life >140 years. The phenolic compounds are biodegradable under the range of redox conditions in the aquifer and the aquifer is not oxidant limited, but the plume is likely to be long-lived and to expand. Degradation is likely to increase only after contaminant concentrations are reduced and aqueous oxidant inputs are increased by dispersion of the plume. The results imply that transport processes may exert a greater control on the natural attenuation of this plume than aquifer oxidant availability.     

Microbial reduction of sulfate injected to gas condensate plumes in cold groundwater

Despite a rapid expansion over the past decade in the reliance on intrinsic bioremediation to remediate petroleum hydrocarbon plumes in groundwater, significant research gaps remain. Although it has been demonstrated that bacterial sulfate reduction can be a key electron accepting process in many petroleum plumes, little is known about the rate of this reduction process in plumes derived from crude oil and gas condensates at cold-climate sites (mean temperature <10 degrees C), and in complex hydrogeological settings such as silt/clay aquitards. In this field study, sulfate was injected into groundwater contaminated by gas condensate plumes at two petroleum sites in Alberta, Canada to enhance in-situ bioremediation. In both cases the groundwater near the water table had low temperature (6-9 degrees C). Monitoring data had provided strong evidence that bacterial sulfate reduction was a key terminal electron accepting process (TEAP) in the natural attenuation of dissolved hydrocarbons at these sites. At each site, water with approximately 2000 mg/L sulfate and a bromide tracer was injected into a low-sulfate zone within a condensate-contaminant plume. Monitoring data collected over several months yielded conservative estimates for sulfate reduction rates based on zero-order kinetics (4-6 mg/L per day) or first-order kinetics (0.003 and 0.01 day(-1)). These results favor the applicability of in-situ bioremediation techniques in this region, under natural conditions or with enhancement via sulfate injection.     

Concentration field and turbulent fluxes during the mixing of two buoyant plumes

In this work an experimental study of mixing of two identical plumes, carried out in a turbulent neutral boundary layer generated in a wind tunnel, is presented. Measurements have been performed with fast flame ionisation detectors (FFIDs) and a two-component Laser-Doppler Anemometer system. Results allow the study of both the average and the fluctuating concentration field, including the turbulent vertical and longitudinal mass fluxes, in single plumes and during the interaction of two identical plumes. This information gives insight into the details of the mixing phase of the two plumes. Results of trajectories and additional rise (due to plume interactions) have been compared with previous measurements carried out in laminar cross-flows, showing similar behaviour. Concentration distributions in plume cross-sections in turbulent cross-flows differ from those measured in laminar cross-flows. Average vertical and longitudinal velocity measurements into the plume core show the strength of the shielding effect of the upwind plume and some details of interaction between the counter-rotating vortex pairs (CVPs). For large values of the alignment angle φ, between the line joining the stacks and the cross-flow, an average negative vertical velocity is measured in the middle of the plume even if its centre of mass is rising. This downward velocity is induced by the slow interaction of the CVPs and generates a vertical stretching of the plume and negative rise enhancement. Vertical turbulent fluxes change sign on the plume centreline and are of opposite sign with respect to the longitudinal turbulent fluxes. Results indicate a good linearity between vertical turbulent fluxes and concentration gradients, with different proportionality for the top and bottom parts of the plume (especially in the near field) indicating that dispersion could be described by a gradient-transfer model.     

Near-field dispersion modeling for regulatory applications

This paper evaluates the application of dispersion models to estimate near-field pollutant concentrations in two case studies. The Industrial Source Complex Short-Term Model (ISCST3) was evaluated with hexavalent chromium measurements collected within 100 m of two facilities in Barrio Logan, San Diego, CA. ISCST3 provided reasonable estimates for higher pollutant concentrations but underestimated lower concentrations. To understand the observed distribution of concentrations in Barrio Logan, a recently conducted tracer experiment was analyzed. The tracer, sulfur hexafluoride, was released at ambient temperature from an urban facility at the University of California at Riverside, and concentrations were measured within 20 m of the source. Modeling results indicated that Industrial Source Complex-Plume Rise Model Enhancement and American Meteorological Society/U.S. Environmental Protection Agency Regulatory Model-Plume Rise Model Enhancement overestimated high concentrations and underestimated low concentrations. A diagnostic study with a simple Gaussian dispersion model that incorporated site-specific meteorology was used to evaluate model results. This study found that incorporating lateral meandering for nonbuoyant urban plumes in Gaussian dispersion models could improve concentration estimates even when downwash is not considered. Incorporating a meandering component in ISCST3 resulted in improvements in estimating hexavalent chromium concentrations in Barrio Logan. Credible near-source concentration estimates depend on accurate characterization of emissions, onsite micrometeorology, and a method to account for lateral meandering in the near field.     

Near-Field Dispersion Modeling for Regulatory Applications

Abstract

This paper evaluates the application of dispersion models to estimate near-field pollutant concentrations in two case studies. The Industrial Source Complex Short-Term Model (ISCST3) was evaluated with hexavalent chromium measurements collected within 100 m of two facilities in Barrio Logan, San Diego, CA. ISCST3 provided reasonable estimates for higher pollutant concentrations but underestimated lower concentrations. To understand the observed distribution of concentrations in Barrio Logan, a recently conducted tracer experiment was analyzed. The tracer, sulfur hexafluoride, was released at ambient temperature from an urban facility at the University of California at Riverside, and concentrations were measured within 20 m of the source. Modeling results indicated that Industrial Source Complex–Plume Rise Model Enhancement and American Meteorological Society/U.S. Environmental Protection Agency Regulatory Model–Plume Rise Model Enhancement overestimated high concentrations and underestimated low concentrations. A diagnostic study with a simple Gaussian dispersion model that incorporated site-specific meteorology was used to evaluate model results. This study found that incorporating lateral meandering for nonbuoyant urban plumes in Gaussian dispersion models could improve concentration estimates even when downwash is not considered. Incorporating a meandering component in ISCST3 resulted in improvements in estimating hexavalent chromium concentrations in Barrio Logan. Credible near-source concentration estimates depend on accurate characterization of emissions, onsite micrometeorology, and a method to account for lateral meandering in the near field.     

The variability and intrinsic remediation of a BTEX plume in anaerobic sulphate-rich groundwater

Data from long-term groundwater sampling, limited coring, and associated studies are synthesised to assess the variability and intrinsic remediation/natural attenuation of a dissolved hydrocarbon plume in sulphate-rich anaerobic groundwater. Fine vertical scale (0.25- and 0.5-m depth intervals) and horizontal plume-scale (>400 m) characteristics of the plume were mapped over a 5-year period from 1991 to 1996. The plume of dissolved BTEX (benzene, toluene, ethylbenzene, xylene) and other organic compounds originated from leakage of gasoline from a subsurface fuel storage tank. The plume was up to 420 m long, less than 50 m wide and 3 m thick. In the first few years of monitoring, BTEX concentrations near the point of leakage were in approximate equilibrium with non-aqueous phase liquid (NAPL) gasoline. NAPL composition of core material and long-term trends in ratios of BTEX concentrations in groundwater indicated significant depletion (water washing, volatilisation and possibly biodegradation) of benzene from residual NAPL after 1992. Large fluctuations in BTEX concentrations in individual boreholes were shown to be largely attributable to seasonal groundwater flow variations. A combination of temporal and spatial groundwater quality data was required to adequately assess the stationarity of plumes, so as to allow inference of intrinsic remediation. Contoured concentration data for the period 1991 to 1996 indicated that plumes of toluene and o-xylene were, at best, only partially steady state (pseudo-steady state) due to seasonal groundwater flow changes. From this analysis, it was inferred that significant remediation by natural biodegradation was occurring for BTEX component plumes such as toluene and o-xylene, but provided no conclusive evidence of benzene biodegradation. Issues associated with field quantification of intrinsic remediation from groundwater sampling are highlighted. Preferential intrinsic biodegradation of selected organic compounds within the BTEX plume was shown to be occurring, in parallel with sulphate reduction and bicarbonate production. Ratios of average hydrocarbon concentrations to benzene for the period 1991 to 1992 were used to estimate degradation rates (half-lives) at various distances along the plume. The estimates varied with distance, the narrowest range being, for toluene, 110 to 260 days. These estimates were comparable to rates determined previously from an in situ tracer test and from plume-scale modelling.     

Assessing the natural attenuation of organic contaminants in aquifers using plume-scale electron and carbon balances: model development with analysis of uncertainty and parameter sensitivity

A quantitative methodology is described for the field-scale performance assessment of natural attenuation using plume-scale electron and carbon balances. This provides a practical framework for the calculation of global mass balances for contaminant plumes, using mass inputs from the plume source, background groundwater and plume residuals in a simplified box model. Biodegradation processes and reactions included in the analysis are identified from electron acceptors, electron donors and degradation products present in these inputs. Parameter values used in the model are obtained from data acquired during typical site investigation and groundwater monitoring studies for natural attenuation schemes. The approach is evaluated for a UK Permo-Triassic Sandstone aquifer contaminated with a plume of phenolic compounds. Uncertainty in the model predictions and sensitivity to parameter values was assessed by probabilistic modelling using Monte Carlo methods. Sensitivity analyses were compared for different input parameter probability distributions and a base case using fixed parameter values, using an identical conceptual model and data set. Results show that consumption of oxidants by biodegradation is approximately balanced by the production of CH4 and total dissolved inorganic carbon (TDIC) which is conserved in the plume. Under this condition, either the plume electron or carbon balance can be used to determine contaminant mass loss, which is equivalent to only 4% of the estimated source term. This corresponds to a first order, plume-averaged, half-life of > 800 years. The electron balance is particularly sensitive to uncertainty in the source term and dispersive inputs. Reliable historical information on contaminant spillages and detailed site investigation are necessary to accurately characterise the source term. The dispersive influx is sensitive to variability in the plume mixing zone width. Consumption of aqueous oxidants greatly exceeds that of mineral oxidants in the plume, but electron acceptor supply is insufficient to meet the electron donor demand and the plume will grow. The aquifer potential for degradation of these contaminants is limited by high contaminant concentrations and the supply of bioavailable electron acceptors. Natural attenuation will increase only after increased transport and dilution.     

Estimation of mechanical dispersion and dispersivity in a soil-gas system by column experiments and the dusty gas model

In a previous study, column experiments were carried out with Toyoura sand (permeability 2.05×10(-11)m(2)) and Toyoura sand mixed with bentonite (permeability 9.96×10(-13)m(2)) to obtain the molecular diffusion coefficient, the Knudsen diffusion coefficient, the tortuosity for the molecular diffusion coefficient, and the mechanical dispersion coefficient of soil-gas systems. In this study, we conducted column experiments with field soil (permeability 2.0×10(-13)m(2)) and showed that the above parameters can be obtained for both less-permeable and more-permeable soils by using the proposed method for obtaining the parameters and performing column experiments. We then estimated dispersivity from the mechanical dispersion coefficients obtained by the column experiments. We found that the dispersivity depended on the mole fraction of the tracer gas and could be represented by a quadratic equation.     

Natural attenuation of chlorinated ethene compounds: model development and field-scale application at the Dover site

A multi-dimensional and multi-species reactive transport model was developed to aid in the analysis of natural attenuation design at chlorinated solvent sites. The model can simulate several simultaneously occurring attenuation processes including aerobic and anaerobic biological degradation processes. The developed model was applied to analyze field-scale transport and biodegradation processes occurring at the Area-6 site in Dover Air Force Base, Delaware. The model was calibrated to field data collected at this site. The calibrated model reproduced the general groundwater flow patterns, and also, it successfully recreated the observed distribution of tetrachloroethene (PCE), trichloroethene (TCE), dichloroethylene (DCE), vinyl chloride (VC) and chloride plumes. Field-scale decay rates of these contaminant plumes were also estimated. The decay rates are within the range of values that were previously estimated based on lab-scale microcosm and field-scale transect analyses. Model simulation results indicated that the anaerobic degradation rate of TCE, source loading rate, and groundwater transport rate are the important model parameters. Sensitivity analysis of the model indicated that the shape and extent of the predicted TCE plume is most sensitive to transmissivity values. The total mass of the predicted TCE plume is most sensitive to TCE anaerobic degradation rates. The numerical model developed in this study is a useful engineering tool for integrating field-scale natural attenuation data within a rational modeling framework. The model results can be used for quantifying the relative importance of various simultaneously occurring natural attenuation processes.     

Biogeochemical processes at the fringe of a landfill leachate pollution plume: potential for dissolved organic carbon, Fe(II), Mn(II), NH4, and CH4 oxidation

Various redox reactions may occur at the fringe of a landfill leachate plume, involving oxidation of dissolved organic carbon (DOC), CH4, Fe(II), Mn(II), and NH4 from leachate and reduction of O2, NO3 and SO4 from pristine groundwater. Knowledge on the relevance of these processes is essential for the simulation and evaluation of natural attenuation (NA) of pollution plumes. The occurrence of such biogeochemical processes was investigated at the top fringe of a landfill leachate plume (Banisveld, the Netherlands). Hydrochemical depth profiles of the top fringe were captured via installation of a series of multi-level samplers at 18, 39 and 58 m downstream from the landfill. Ten-centimeter vertical resolution was necessary to study NA within a fringe as thin as 0.5 m. Bromide appeared an equally well-conservative tracer as chloride to calculate dilution of landfill leachate, and its ratio to chloride was high compared to other possible sources of salt in groundwater. The plume fringe rose steadily from a depth of around 5 m towards the surface with a few meters in the period 1998-2003. The plume uplift may be caused by enhanced exfiltration to a brook downstream from the landfill, due to increased precipitation over this period and an artificial lowering of the water level of the brook. This rise invoked cation exchange including proton buffering, and triggered degassing of methane. The hydrochemical depth profile was simulated in a 1D vertical reactive transport model using PHREEQC-2. Optimization using the nonlinear optimization program PEST brought forward that solid organic carbon and not clay minerals controlled retardation of cations. Cation exchange resulted in spatial separation of Fe(II), Mn(II) and NH4 fronts from the fringe, and thereby prevented possible oxidation of these secondary redox species. Degradation of DOC may happen in the fringe zone. Re-dissolution of methane escaped from the plume and subsequent oxidation is an explanation for absence of previously present nitrate and anaerobic conditions in pristine groundwater above the plume. Stable carbon isotope (delta13C) values of methane confirm anaerobic methane oxidation immediately below the fringe zone, presumably coupled to reduction of sulfate, desorbed from iron oxide. Methane must be the principle reductant consuming soluble electron-acceptors in pristine groundwater, thereby limiting NA for other solutes including organic micro-pollutants at the fringe of this landfill leachate plume.     

A fluctuating plume model for concentration fluctuations in a plant canopy

A one-particle Lagrangian model for continuous releases in the non-Gaussian inhomogeneous turbulence of a canopy layer is derived based on the fluctuating plume model of Franzese [2003. Lagrangian stochastic modeling of a fluctuating plume in the convective boundary layer. Atmos. Environ. 37, 1691–1701.]. The model equations are filtered by a time-dependent low-pass filter applied to the turbulent kinetic energy in order to obtain a fluctuating plume model able to simulate the vertical meandering of the cloud centroid through non-stationary Lagrangian equations. The model satisfies the well-mixed condition. The relative dispersion of particles and the concentration fluctuation statistics of a passive tracer inside a modeled vegetal canopy are studied. The probability density function of the concentration relative to the plume centroid is parameterized and the mean and variance fields of concentration are simulated and compared with wind-tunnel data and numerical simulations. A skewed, reflected probability density function for the vertical position of the plume centroid is considered.     

Fluorescent dye imaging of the volume sampled by single well forced-gradient tracer tests evaluated in a laboratory-scale aquifer physical model

This study presents a new method to visualise forced-gradient tracer tests in 2-D using a laboratory-scale aquifer physical model. Experiments were designed to investigate the volume of aquifer sampled in vertical dipole flow tracer tests (DFTT) and push-pull tests (PPT), using a miniature monitoring well and straddle packer arrangement equipped with solute injection and recovery chambers. These tests have previously been used to estimate bulk aquifer hydraulic and transport properties for the evaluation of natural attenuation and other remediation approaches. Experiments were performed in a silica glass bead-filled box, using a fluorescent tracer (fluorescein) to deduce conservative solute transport paths. Digital images of fluorescein transport were captured under ultraviolet light and processed to analyse tracer plume geometry and obtain point-concentration breakthrough histories. Inorganic anion mixtures were also used to obtain conventional tracer breakthrough histories. Concentration data from the conservative tracer breakthrough curves was compared with the digital images and a well characterised numerical model. The results show that the peak tracer breakthrough response in dipole flow tracer tests samples a zone of aquifer close to the well screen, while the sampling volume of push-pull tests is limited by the length of the straddle packers used. The effective sampling volume of these single well forced-gradient tests in isotropic conditions can be estimated with simple equations. The experimental approach offers the opportunity to evaluate under controlled conditions the theoretical basis, design and performance of DFTTs and PPTs in porous media in relation to measured flow and transport properties.     

Natural attenuation of xenobiotic organic compounds in a landfill leachate plume (Vejen,Denmark)

Demonstration of natural attenuation of xenobiotic organic compounds (XOCs) in landfill leachate plumes is a difficult task and still an emerging discipline within groundwater remediation. One of the early studies was made at the Vejen Landfill in Denmark in the late 1980s, which suggested that natural attenuation of XOCs took place under strongly anaerobic conditions within the first 150 m of the leachate plume. This paper reports on a revisit to the same plume 10 years later. Within the strongly anaerobic part of the plume, 49 groundwater samples were characterized with respect to redox-sensitive species and XOCs. The analytical procedures have been developed further and more compounds and lower detection limits were observed this time. In addition, the samples were screened for degradation intermediates and for toxicity. The plume showed fairly stationary features over the 10-year period except that the XOC level as well as the level of chloride and nonvolatile organic carbon (NVOC) in the plume had decreased somewhat. Most of the compounds studied were subject to degradation in addition to dilution. Exceptions were benzene, the herbicide Mecoprop (MCPP), and NVOC. In the early study, NVOC seemed to degrade in the first part of the plume, but this was no longer the case. Benzyl succinic acid (BSA) was for the first time identified in a leachate plume as a direct indicator, and as the only intermediate of toluene degradation. Toxicity measurements on solid phase-extracted (SPE) samples revealed that toxic compounds not analytically identified were still present in the plume, suggesting that toxicity measurements could be helpful in assessing natural attenuation in leachate plumes.     

Estimation of local scale dispersion from local breakthrough curves during a tracer test in a heterogeneous aquifer: the Lagrangian approach

The local scale dispersion tensor, Dd, is a controlling parameter for the dilution of concentrations in a solute plume that is displaced by groundwater flow in a heterogeneous aquifer. In this paper, we estimate the local scale dispersion from time series or breakthrough curves, BTCs, of Br concentrations that were measured at several points in a fluvial aquifer during a natural gradient tracer test at Krauthausen. Locally measured BTCs were characterized by equivalent convection dispersion parameters: equivalent velocity, v(eq)(x) and expected equivalent dispersivity, [lambda(eq)(x)]. A Lagrangian framework was used to approximately predict these equivalent parameters in terms of the spatial covariance of log(e) transformed conductivity and the local scale dispersion coefficient. The approximate Lagrangian theory illustrates that [lambda(eq)(x)] increases with increasing travel distance and is much larger than the local scale dispersivity, lambda(d). A sensitivity analysis indicates that [lambda(eq)(x)] is predominantly determined by the transverse component of the local scale dispersion and by the correlation scale of the hydraulic conductivity in the transverse to flow direction whereas it is relatively insensitive to the longitudinal component of the local scale dispersion. By comparing predicted [lambda(eq)(x)] for a range of Dd values with [lambda(eq)(x)] obtained from locally measured BTCs, the transverse component of Dd, DdT, was estimated. The estimated transverse local scale dispersivity, lambda(dT) = DdT/U1 (U1 = mean advection velocity) is in the order of 10(1)-10(2) mm, which is relatively large but realistic for the fluvial gravel sediments at Krauthausen.     

Comment on "Analysis of groundwater contamination using concentration-time series recorded during an integral pumping test: Bias introduced by strong concentration gradients within the plume" by Allelign Zeru and Gerhard Schäfer

We consider the results of a recent paper in this journal [Zeru, A. and Sch?fer, G., 2005. Analysis of groundwater contamination using concentration-time series recorded during an integral pumping test: Bias introduced by strong concentration gradients within the plume. Journal of Contaminant Hydrology 81 (2005) 106-124], which addresses the field-scale characterisation of contaminant plumes in groundwater. There, it is concluded that contaminant concentration gradients can bias Integral Pumping Test (IPT) interpretations considerably, in particular if IPTs are conducted in advective fronts of contaminant plumes. We discuss implications of this setting and also argue that the longitudinal and transverse dispersivities used in the examples of Zeru and Sch?fer (2005) of up to 30 m and 3 m, respectively, are generally very high for the here relevant capture zone scale (<20 m). However, regardless of both longitudinal and transverse concentration gradients, we further show through a counter-example that IPT results are unbiased as long as the concentration attenuation along the flow direction is linear over the capture zone extent.