Hydraulic constraints on the performance of a groundwater denitrification wall for nitrate removal from shallow groundwater

Denitrification walls are a practical approach for decreasing non-point source pollution of surface waters. They are constructed by digging a trench perpendicular to groundwater flow and mixing the aquifer material with organic matter, such as sawdust, which acts as a carbon source to stimulate denitrification. For efficient functioning, walls need to be permeable to groundwater flow. We examined the functioning of a denitrification wall constructed in an aquifer consisting of coarse sands. Wells were monitored for changes in nitrate concentration as groundwater passed through the wall and soil samples were taken to measure microbial parameters inside the wall. Nitrate concentrations upstream of the wall ranged from 21 to 39 g N m(-3), in the wall from 0 to 2 g N m(-3) and downstream from 19 to 44 g N m(-3). An initial groundwater flow investigation using a salt tracer dilution technique showed that the flow through the wall was less than 4% of the flow occurring in the aquifer. Natural gradient tracer tests using bromide and Rhodamine-WT confirmed groundwater bypass under the wall. Hydraulic conductivity of 0.48 m day(-1) was measured inside the wall, whereas the surrounding aquifer had a hydraulic conductivity of 65.4 m day(-1). This indicated that during construction of the wall, hydraulic conductivity of the aquifer had been greatly reduced, so that most of the groundwater flowed under rather than through the wall. Denitrification rates measured in the center of the wall ranged from 0.020 to 0.13 g N m(-3) day(-1), which did not account for the rates of nitrate removal (0.16-0.29 g N m(-3) day(-1)) calculated from monitoring of groundwater nitrate concentrations. This suggested that the rate of denitrification was greater at the upstream face of the wall than in its center where it was limited by low nitrate concentrations. While denitrification walls can be an inexpensive tool for removing nitrate from groundwater, they may not be suitable in aquifers with coarse textured subsoils where simple inexpensive construction techniques result in major decreases in hydraulic conductivity.     

Evaluation of groundwater flow patterns around a dual-screened groundwater circulation well

Dual-screened groundwater circulation wells (GCWs) can be used to remove contaminant mass and to mix reagents in situ. GCWs are so named because they force water in a circular pattern between injection and extraction screens. The radial extent, flux and direction of the effective flow of this circulation cell are difficult to measure or predict. The objective of this study is to develop a robust protocol for assessing GCW performance. To accomplish this, groundwater flow patterns surrounding a GCW are assessed using a suite of tools and data, including: hydraulic head, in situ flow velocity, measured hydraulic conductivity data from core samples, chemical tracer tests, contaminant distribution data, and numerical flow and transport models. The hydraulic head data show patterns that are consistent with pumping on a dual-screened well, however, many of the observed changes are smaller than expected. In situ thermal perturbation flow sensors successfully measured horizontal flow, but vertical flow could not be determined with sufficient accuracy to be useful in mapping flow patterns. Two types of chemical tracer tests were utilized at the site and showed that much of the flow occurs within a few meters of the GCW. Flow patterns were also assessed based on changes in contaminant (trichloroethylene, TCE) concentrations over time. The TCE data clearly showed treated water moving away from the GCW at shallow and intermediate depths, but the circulation of that water back to the well, except very close to the well, was less clear. Detailed vertical and horizontal hydraulic conductivities were measured on 0.3 m-long sections from a continuous core from the GCW installation borehole. The measured vertical and horizontal hydraulic conductivity data were used to construct numerical flow and transport models, the results of which were compared to the head, velocity and concentration data. Taken together, the field data and modeling present a fairly consistent picture of flow and transport around the GCW. However, the time and expense associated with conducting all of those tests would be prohibitive for most sites. As a consequence, a sequential protocol for GCW characterization is presented here in which the number of tools used can be adjusted to meet the needs of individual sites. While not perfect, we believe that this approach represents the most efficient means for evaluating GCW performance.     

Analytical solutions for flow fields near continuous wall reactive barriers

Permeable reactive barriers (PRBs) are widely applied for in-situ remediation of contaminant plumes transported by groundwater. Besides the goal of a sufficient contaminant remediation inside the reactive cell (residence time) the width of plume intercepted by a PRB is of critical concern. A 2-dimensional analytical approach is applied to determine the flow fields towards rectangular PRBs of the continuous wall (CW) configuration with and without impermeable side walls (but yet no funnel). The approach is based on the conformal mapping technique and assumes a homogeneous aquifer with a uniform ambient flow field. The hydraulic conductivity of the reactive material is furthermore assumed to exceed the conductivity of the aquifer by at least one order of magnitude as to neglect the hydraulic gradient across the reactor. The flow fields are analyzed regarding the widths and shapes of the respective capture zones as functions of the dimensions (aspect ratio) of the reactive cell and the ambient groundwater flow direction. Presented are an improved characterization of the advantages of impermeable side walls, a convenient approach to improved hydraulic design (including basic cost-optimization) and new concepts for monitoring CW PRBs. Water level data from a CW PRB at the Seneca Army Depot site, NY, are used for field demonstration.     

Plume persistence caused by back diffusion from thin clay layers in a sand aquifer following TCE source-zone hydraulic isolation

This paper concludes that back diffusion from one or a few thin clayey beds in a sand aquifer can cause contaminant persistence above MCLs in a sand aquifer long after the source zone initially causing the plume is isolated or removed. This conclusion is based on an intensive case study of a TCE contaminated site in Florida, with the processes evaluated using numerical modeling. At this site, the TCE DNAPL zone formed decades ago, and was hydraulically isolated by means of an innovative system performing groundwater extraction, treatment and re-injection. Treated water is re-injected in a row of injection wells situated a short distance downgradient of the extraction wells, creating a clean-water displacement front to efficiently flush the downgradient plume. This scheme avoids the creation of stagnation zones typical of most groundwater pump-and-treat systems, thereby minimizing the time for aquifer flushing and therefore downgradient cleanup. The system began operation in August 2002 and although the performance monitoring shows substantial declines in concentrations, detectable levels of TCE and degradation products persist downgradient of the re-injection wells, long after the TCE should have disappeared based on calculations assuming a nearly homogenous sand aquifer. Three hypotheses were assessed for this plume persistence: 1) incomplete source-zone capture, 2) DNAPL occurrence downgradient of the re-injection wells, and 3) back diffusion from one or more thin clay beds in the aquifer. After careful consideration, the first two hypotheses were eliminated, leaving back diffusion as the only plausible hypothesis, supported by detailed measurements of VOC concentrations within and near the clay beds and also by numerical model simulations that closely represent the field site hydrogeologic conditions. The model was also used to simulate a more generalized, hypothetical situation where more thin clayey beds occur in a sand aquifer with an underlying aquitard. While there is no doubt that DNAPL source mass reduction can eventually improve downgradient groundwater quality, the magnitude and time scale over which the improvement occurs is the major uncertainty given current characterization approaches. This study shows that even one thin clay bed, less than 0.2 m thick, can cause plume persistence due to back diffusion for several years or even decades after the flux from the source is completely isolated. Thin clay beds, which have a large storage capacity for dissolved and sorbed contaminant mass, are common in many types of sandy aquifers. However, without careful inspection of continuous cores and sampling, such thin clay beds, and their potential for causing long-term back-diffusion effects, can easily go unnoticed during site characterization.     

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.     

Estimating contaminant mass discharge: a field comparison of the multilevel point measurement and the integral pumping investigation approaches and their uncertainties

In this field study, two approaches to assess contaminant mass discharge were compared: the sampling of multilevel wells (MLS) and the integral groundwater investigation (or integral pumping test, IPT) that makes use of the concentration-time series obtained from pumping wells. The MLS approached used concentrations, hydraulic conductivity and gradient rather than direct chemical flux measurements, while the IPT made use of a simplified analytical inversion. The two approaches were applied at a control plane located approximately 40m downgradient of a gasoline source at Canadian Forces Base Borden, Ontario, Canada. The methods yielded similar estimates of the mass discharging across the control plane. The sources of uncertainties in the mass discharge in each approach were evaluated, including the uncertainties inherent in the underlying assumptions and procedures. The maximum uncertainty of the MLS method was about 67%, and about 28% for the IPT method in this specific field situation. For the MLS method, the largest relative uncertainty (62%) was attributed to the limited sampling density (0.63 points/m(2)), through a novel comparison with a denser sampling grid nearby. A five-fold increase of the sampling grid density would have been required to reduce the overall relative uncertainty for the MLS method to about the same level as that for the IPT method. Uncertainty in the complete coverage of the control plane provided the largest relative uncertainty (37%) in the IPT method. While MLS or IPT methods to assess contaminant mass discharge are attractive assessment tools, the large relative uncertainty in either method found for this reasonable well monitored and simple aquifer suggests that results in more complex plumes in more heterogeneous aquifers should be viewed with caution.     

Influence of temporally variable groundwater flow conditions on point measurements and contaminant mass flux estimations

Monitoring of contaminant concentrations, e.g., for the estimation of mass discharge or contaminant degradation rates, often is based on point measurements at observation wells. In addition to the problem, that point measurements may not be spatially representative, a further complication may arise due to the temporal dynamics of groundwater flow, which may cause a concentration measurement to be not temporally representative. This paper presents results from a numerical modeling study focusing on temporal variations of the groundwater flow direction. “Measurements” are obtained from point information representing observation wells installed along control planes using different well frequencies and configurations. Results of the scenario simulations show that temporally variable flow conditions can lead to significant temporal fluctuations of the concentration and thus are a substantial source of uncertainty for point measurements. Temporal variation of point concentration measurements may be as high as the average concentration determined, especially near the plume fringe, even when assuming a homogeneous distribution of the hydraulic conductivity. If a heterogeneous hydraulic conductivity field is present, the concentration variability due to a fluctuating groundwater flow direction varies significantly within the control plane and between the different realizations. Determination of contaminant mass fluxes is also influenced by the temporal variability of the concentration measurement, especially for large spacings of the observation wells. Passive dosimeter sampling is found to be appropriate for evaluating the stationarity of contaminant plumes as well as for estimating average concentrations over time when the plume has fully developed. Representative sampling has to be performed over several periods of groundwater flow fluctuation. For the determination of mass fluxes at heterogeneous sites, however, local fluxes, which may vary considerably along a control plane, have to be accounted for. Here, dosimeter sampling in combination with time integrated local water flux measurements can improve mass flux estimates under dynamic flow conditions.     

Assessment of Aquifer Vulnerability in an Agricultural Area in Spain Using the DRASTIC Model

The assessment of aquifer vulnerability is a very important task, especially in agricultural areas because the quality and availability of groundwater affects both the sustainability of agriculture and the quality of life. In this study, an integrated approach is considered, with the use of the generic and agricultural DRASTIC models as well as a geographic information system (GIS), to assess groundwater vulnerability in the agricultural area of Barrax, in the province of Albacete, in Spain. Seven parameters—depth to water, net recharge, aquifer media, soil media, topography, impact of vadose zone media, and hydraulic conductivity of the aquifer (DRASTIC)—have been considered as weighted layers to enable an accurate groundwater risk mapping. The results of the generic DRASTIC model indicated very low vulnerability to contamination for Barrax groundwater due to limited urban and industrial development in the wider area. However, agricultural activities impose pressure to groundwater resources and the results of the agricultural DRASTIC model show that 6.86% of the study area is characterized by very high, 2.29% by high, 47.28% by medium, 38.28% by low, and the remaining 5.29% by no vulnerability to groundwater contamination. The distribution of nitrate concentration in groundwater in the area under study is quite well correlated with the agricultural DRASTIC vulnerability index. Sensitivity analysis was also performed to acknowledge statistical uncertainty in the estimation of each parameter used, to assess its impact, and thus to identify the most critical parameters that require further investigation. Depth to water and impact of vadose zone are the parameters that had the most noticeable impact on the generic DRASTIC vulnerability index followed by the soil media and topography. In contrast, the agricultural DRASTIC method is more sensitive to the removal of the depth to water parameter followed by the topography and the soil media parameters.     

Probabilistic study of well capture zones distribution at the Lauswiesen field site

The delineation of well capture zones is of utmost environmental and engineering relevance as pumping wells are commonly used both for drinking water supply needs, where protection zones have to be defined, and for investigation and remediation of contaminated aquifers. We analyze the probabilistic nature of well capture zones within the well field located at the "Lauswiesen" experimental site. The test site is part of an alluvial heterogeneous aquifer located in the Neckar river valley, close to the city of Tübingen in South-West Germany. We explore the effect of different conceptual models of the structure of aquifer heterogeneities on the delineation of three-dimensional probabilistic well catchment and time-related capture zones, in the presence of migration of conservative solutes. The aquifer is modeled as a three-dimensional, doubly stochastic composite medium, where distributions of geo-materials and hydraulic properties are uncertain. We study the relative importance of uncertain facies geometry and uncertain hydraulic conductivity and porosity on predictions of catchment and solute time of travel to the pumping well by focusing on cases in which (1) the facies distribution is random, but the hydraulic properties of each material are fixed, and (2) both facies geometry and material properties vary stochastically. The problem is tackled within a conditional numerical Monte Carlo framework. Results are provided in terms of probabilistic demarcations of the three-dimensional well catchment and time-related capture zones. Our findings suggest that the uncertainty associated with the prediction of the location of the outer boundary of well catchment at the "Lauswiesen" site is significantly affected by the conceptual model adopted to incorporate the heterogeneous nature of the aquifer domain in a predictive framework. Taking into account randomness of both lithofacies distribution and materials hydraulic conductivity allows recognizing the existence of preferential flow paths that influence the extent of the well catchment and the solute travel time distribution at the site.     

Characterisation of microbial activity in the framework of natural attenuation without groundwater monitoring wells?: a new Direct-Push probe

At many contaminated field sites in Europe, monitored natural attenuation is a feasible site remediation option. Natural attenuation includes several processes but only the microbial degradation leads to real contaminant removal and very few methods are accepted by the authorities providing real evidence of microbial contaminant degradation activity. One of those methods is the recently developed in situ microcosm approach (BACTRAP®). These in situ microcosms consist of perforated stainless steel cages or PTFE tubes filled with an activated carbon matrix that is amended with ¹³C-labelled contaminants; the microcosms are then exposed within groundwater monitoring wells. Based on this approach, natural attenuation was accepted by authorities as a site remediation option for the BTEX-polluted site Zeitz in Germany. Currently, the in situ microcosms are restricted to the use inside groundwater monitoring wells at the level of the aquifer. The (classical) system therefore is only applicable on field sites with a network of monitoring wells, and only microbial activity inside the monitoring wells at the level of the aquifer can be assessed. In order to overcome these limitations, a new Direct-Push BACTRAP probe was developed on the basis of the Geoprobe® equipment. With respect to the mechanical boundary conditions of the DP technique, these new probes were constructed in a rugged and segmented manner and are adaptable to various sampling concepts. With this new probe, the approach can be extended to field sites without existing monitoring wells, and microbial activity was demonstrated to be measureable even under very dry conditions inside the vadose zone above the aquifer. In a field test, classical and Direct-Push BACTRAPs were applied in the BTEX-contaminated aquifer at the ModelPROBE reference site Zeitz (Germany). Both types of BACTRAPs were incubated in the centre and at the fringe of the BTEX plume. Analysis of phospholipid fatty acid (PLFA) patterns showed that the bacterial communities on DP-BACTRAPs were more similar to the soil than those found on classical BACTRAPs. During microbial degradation of the ¹³C-labelled substrate on the carrier material of the microcosms, the label was only slightly incorporated into bacterial biomass, as determined by PLFA analysis. This provides clear indication for decreased in situ natural attenuation potential in comparison to earlier sampling campaigns, which is presumably caused by a large-scale source remediation measure in the meantime. In conclusion, Direct-Push-based BACTRAPs offer a promising way to monitor natural attenuation or remediation success at field sites which are currently inaccessible by the technique due to the lack of monitoring wells or due to a main contamination present within the vadose zone.     

Estimating spatially-variable first-order rate constants in groundwater reactive transport systems

Numerical reactive transport models are often used as tools to assess aquifers contaminated with reactive groundwater solutes as well as investigating mitigation scenarios. The ability to accurately simulate the fate and transport of solutes, however, is often impeded by a lack of information regarding the parameters that define chemical reactions. In this study, we employ a steady-state Ensemble Kalman Filter (EnKF), a data assimilation algorithm, to provide improved estimates of a spatially-variable first-order rate constant λ through assimilation of solute concentration measurement data into reactive transport simulation results. The methodology is applied in a steady-state, synthetic aquifer system in which a contaminant is leached to the saturated zone and undergoes first-order decay. Multiple sources of uncertainty are investigated, including hydraulic conductivity of the aquifer and the statistical parameters that define the spatial structure of the parameter field. For the latter scenario, an iterative method is employed to identify the statistical mean of λ of the reference system. Results from all simulations show that the filter scheme is successful in conditioning the λ ensemble to the reference λ field. Sensitivity analyses demonstrate that the estimation of the λ values is dependent on the number of concentration measurements assimilated, the locations from which the measurement data are collected, the error assigned to the measurement values, and the correlation length of the λ fields.     

An analytical solution for two-dimensional contaminant transport during groundwater extraction

We obtain an analytical solution for two-dimensional steady-state transport of conservative contaminant between injecting and pumping wells. Flow and transport are considered in the vertical cross-section. The Dupuit approximation and conformal mapping onto the complex potential domain are employed to determine the velocity and concentration distributions, respectively. We use this solution to derive a priori conditions under which widely used 1-D analytical solutions with constant velocity and dispersion coefficients provide accurate approximations. These conditions are formulated in terms of aquifer parameters, such as hydraulic conductivity, porosity and dispersivities, and remediation strategy, e.g., well spacing and pumping regimes.     

Multi tracer test for the implementation of enhanced in-situ bioremediation at a BTEX-contaminated megasite

At the Centre for Environmental Research Leipzig-Halle (UFZ) research site in Zeitz, Germany, benzene contaminates the lower of two aquifers with concentrations of up to 20 mg/l. Since the benzene plume has a minimum length of approximately 1 km, enhanced natural attenuation measures are being considered as a remediation strategy. This study describes the performance and evaluation of a multi-species reactive tracer test using the tracers fluorescein and bromide as conservative tracers and toluene as reactive tracer. Sampling was performed over a period of six months using a detailed network of multilevel sampling wells. Toluene was only slightly retarded in comparison to bromide, whereas fluorescein was retarded considerably stronger. Therefore, it was not possible to use fluorescein as an in situ tracer for the determination of groundwater velocities. The ionic nature of fluorescein is assumed to be the major reason for its retardation. The results show that the infiltration conditions were suitable to produce a wide spreading of the tracer front along the full thickness of the aquifer. Thus, a large aquifer volume can be treated in future enhanced bioremediation measures. The total quantity of infiltrated toluene (24 l) was degraded under sulfate-reducing conditions over a flow path of 50 m. Benzylsuccinate was identified as a metabolite of toluene degradation under sulfate-reducing conditions at this site. The modelling results show that toluene degradation was described more accurately using Monod kinetics than first-order kinetics. Since toluene was only slightly retarded in comparison to bromide, sorption and desorption processes were considered to be negligible.     

Pulsed pumping process optimization using a potential flow model

A computational model is applied to the optimization of pulsed pumping systems for efficient in situ remediation of groundwater contaminants. In the pulsed pumping mode of operation, periodic rather than continuous pumping is used. During the pump-off or trapping phase, natural gradient flow transports contaminated groundwater into a treatment zone surrounding a line of injection and extraction wells that transect the contaminant plume. Prior to breakthrough of the contaminated water from the treatment zone, the wells are activated and the pump-on or treatment phase ensues, wherein extracted water is augmented to stimulate pollutant degradation and recirculated for a sufficient period of time to achieve mandated levels of contaminant removal. An important design consideration in pulsed pumping groundwater remediation systems is the pumping schedule adopted to best minimize operational costs for the well grid while still satisfying treatment requirements. Using an analytic two-dimensional potential flow model, optimal pumping frequencies and pumping event durations have been investigated for a set of model aquifer-well systems with different well spacings and well-line lengths, and varying aquifer physical properties. The results for homogeneous systems with greater than five wells and moderate to high pumping rates are reduced to a single, dimensionless correlation. Results for heterogeneous systems are presented graphically in terms of dimensionless parameters to serve as an efficient tool for initial design and selection of the pumping regimen best suited for pulsed pumping operation for a particular well configuration and extraction rate. In the absence of significant retardation or degradation during the pump-off phase, average pumping rates for pulsed operation were found to be greater than the continuous pumping rate required to prevent contaminant breakthrough.     

Predicting natural attenuation of xylene in groundwater using a numerical model.

The aquifer beneath an abandoned refinery in the Lower Rhine area, Germany, was contaminated with a number of different mineral oil products. Groundwater sampling in the area around the former xylene plant revealed that a xylene plume had developed in the underlying groundwater, and moreover, that there is strong evidence for in situ microbial xylene degradation with oxygen, nitrate, sulfate and ferric iron as electron acceptors. In order to prevent further xylene spreading, three pumping wells extracting contaminated water were installed downgradient of the spill zone. The numerical reactive transport code Transport Biochemisty Chemistry (TBC) was applied to this situation to quantify the relation of microbial degradation to xylene removal by the pumping wells. It could be shown that the unamended in situ degradation was an appreciable xylene removal process that contributed to about one-third to the total xylene removal (degradation plus extraction). A further objective of the model application was to predict xylene spreading under regional flow conditions, i.e. without operation of the three pumping wells, to consider the possible effects of natural xylene attenuation. To accomplish this, the model calibrated for the situation with operating wells was transferred to the hydraulic situation of regional flow while retaining the parameters of the biochemical model. It turned out that the xylene plume that is expected to develop downgradient of the source area will be limited to an extension of not more than 1000 m. An interesting feature of the simulations results was that xylene degradation under iron-reducing conditions, which was of minor importance for the situation with operating pumping wells, becomes the dominant degradation mechanism under regional flow conditions. Moreover, iron reduction will be the key process in controlling plume evolution. The model application illustrates that multi-species reactive transport models are needed to adequately transfer reactive processes from one hydraulic situation to another, while single species models are not suited for this predictive task.     

Quantification of mass fluxes and natural attenuation rates at an industrial site with a limited monitoring network: a case study

At many "real world" field sites, the number of available monitoring wells is limited due to economic or geological reasons. Under such restricted conditions, it is difficult to perform a reliable field investigation and to quantify primary lines of evidence for natural attenuation (NA), like the documentation of a decrease of contaminant mass flux in flow direction. This study reports the results of a groundwater investigation at a former manufactured gas plant situated in a Quaternary river valley in southwest Germany. The location, infrastructure and aquifer setting are typical of many industrial sites in Germany. Due to difficult drilling conditions (coarse glaciofluvial gravel deposits and an anthropogenic fill above the aquifer), only 12 monitoring wells were available for the investigation and localisation of the contaminant plume. These wells were situated along three control planes (CP) downgradient from the contaminant source, with four wells along each plane. Based on the sparse set of monitoring wells, field scale mass fluxes and first-order natural attenuation rate constants of benzene, toluene, ethylbenzene, and o-xylene and p-xylene (BTEX) and low molecular weight polycyclic aromatic hydrocarbons (PAH) were estimated utilizing different point scale and also a new integral investigation method. The results show that even at a heterogeneous site with a sparse monitoring network point scale investigation methods can provide reliable information on field scale natural attenuation rates, if a dependable flow model or tracer test data is available. If this information is not available, only the new integral investigation method presented can yield adequate results for the quantification of contaminant mass fluxes under sparse monitoring conditions.     

A sequential zero valent iron and aerobic biodegradation treatment system for nitrobenzene

The remediation of nitroaromatic contaminated groundwater is sometimes difficult because nitroaromatic compounds are resistant to biodegradation and, when they do transform, the degradation of the products may also be incomplete. A simple nitroaromatic compound, nitrobenzene, was chosen to assess the feasibility of an in situ multi-zone treatment system at the laboratory scale. The proposed treatment system consists of a zero valent granular iron zone to reduce nitrobenzene to aniline, followed by a passive oxygen release zone for the aerobic biodegradation of the aniline daughter product using pristine aquifer material from Canadian Forces Base (CFB) Borden, Ontario, as an initial microbial source. In laboratory batch experiments, nitrobenzene was found to reduce quickly in the presence of granular iron forming aniline, which was not further degraded but remained partially sorbed onto the granular iron surface. Aniline was found to be readily biodegraded with little metabolic lag under aerobic conditions using the pristine aquifer material. A sequential column experiment, containing a granular iron reducing zone and an aerobic biodegradation zone, successively degraded nitrobenzene and then aniline to below detection limits (0.5 microM) without any noticeable reduction in hydraulic conductivity from biofouling, or through the formation of precipitates.     

The significance of heterogeneity on mass flux from DNAPL source zones: an experimental investigation

Understanding the process of mass transfer from source zones of aquifers contaminated with organic chemicals in the form of dense non-aqueous phase liquids (DNAPL) is of importance in site management and remediation. A series of intermediate-scale tank experiments was conducted to examine the influence of aquifer heterogeneity on DNAPL mass transfer contributing to dissolved mass emission from source zone into groundwater under natural flow before and after remediation. A Tetrachloroethylene (PCE) spill was performed into six source zone models of increasing heterogeneity, and both the spatial distribution of the dissolution behavior and the net effluent mass flux were examined. Experimentally created initial PCE entrapment architecture resulting from the PCE migration was largely influenced by the coarser sand lenses and the PCE occupied between 30 and 60% of the model aquifer depth. The presence of DNAPL had no apparent effect on the bulk hydraulic conductivity of the porous media. Up to 71% of PCE mass in each of the tested source zone was removed during a series of surfactant flushes, with associated induced PCE mobilization responsible for increasing vertical DNAPL distributions. Effluent mass flux due to water dissolution was also found to increase progressively due to the increase in NAPL-water contact area even though the PCE mass was reduced. Doubling of local groundwater flow velocities showed negligible rate-limited effects at the scale of these experiments. Thus, mass transfer behavior was directly controlled by the morphology of DNAPL within each source zone. Effluent mass flux values were normalized by the up-gradient DNAPL distributions. For the suite of aquifer heterogeneities and all remedial stages, normalized flux values fell within a narrow band with mean of 0.39 and showed insensitivity to average source zone saturations.     

Temporal variability of groundwater chemistry in shallow and deep aquifers of Araihazar, Bangladesh

Samples were collected every 2-4 weeks from a set of 37 monitoring wells over a period of 2-3 years in Araihazar, Bangladesh, to evaluate the temporal variability of groundwater composition for As and other constituents. The monitoring wells are grouped in 6 nests and span the 5-91 m depth range. Concentrations of As, Ca, Fe, K, Mg, Mn, Na, P, and S were measured by high-resolution ICPMS with a precision of 5% or better; concentrations of Cl were measured by ion chromatography. In shallow wells <30 m deep, As and P concentrations generally varied by <30%, whereas concentrations of the major ions (Na, K, Mg, Ca and Cl) and the redox-sensitive elements (Fe, Mn, and S) varied over time by up to +/-90%. In wells tapping the deeper aquifers >30 m often below clay layers concentrations of groundwater As were much lower and varied by <10%. The concentrations of major cations also varied by <10% in these deep aquifers. In contrast, the concentration of redox-sensitive constituents Fe, S, and Mn in deep aquifers varied by up to 97% over time. Thus, strong decoupling between variations in As and Fe concentrations is evident in groundwaters from shallow and deep aquifers. Comparison of the time series data with groundwater ages determined by (3)H/(3)He and (14)C dating shows that large seasonal or inter-annual variations in major cation and chloride concentrations are restricted to shallow aquifers and groundwater recharged <5 years ago. There is no corresponding change in As concentrations despite having significant variations of redox sensitive constituents in these very young waters. This is attributed to chemical buffering due to rapid equilibrium between solute and solid As. At two sites where the As content of groundwater in existing shallow wells averages 102 microg/L (range: <5 to 648 microg/L; n=118) and 272 microg/L (range: 10 to 485 microg/L; n=65), respectively, a systematic long-term decline in As concentrations lends support to the notion that flushing may slowly deplete an aquifer of As. Shallow aquifer water with >5 years (3)H/(3)He age show a constant As:P molar ratio of 9.6 over time, suggesting common mechanisms of mobilization.