Time scales of organic contaminant dissolution from complex source zones: coal tar pools vs. blobs

Groundwater contamination due to complex organic mixtures such as coal tar, creosote and fuels is a widespread problem in industrialized regions. Although most compounds in these mixtures are biodegradable, the contaminant sources are very persistent for many decades after the contamination occurred (e.g., more than 100 years ago at gasworks sites). This limited bioavailability is due to slow dissolution processes. This study presents results from a large scale tank experiment (8 m long) on the long-term (354 days) dissolution kinetics of BTEX and PAHs from a 2.5 m long coal tar pool and 0.5 m long (smear) zone containing coal tar blobs distributed in a coarse sand. The results inidicate (1) that Raoult's law holds for estimation of the saturation aqueous concentrations of the coal tar constituents, (2) that for the dissolution of smear zones longer than approximately 0.1 m and with more than 3-5% residual saturation, the local equilibrium assumption is valid and (3) that although very small (< 0.1 mm), the transverse vertical dispersivity dominates the pool dissolution processes. Typical time scales for removal of the pollutants from the blob zone and the pool are in the order of a few weeks to more than 10,000 years, respectively.     

Rebound of a coal tar creosote plume following partial source zone treatment with permanganate

The long-term management of dissolved plumes originating from a coal tar creosote source is a technical challenge. For some sites stabilization of the source may be the best practical solution to decrease the contaminant mass loading to the plume and associated off-site migration. At the bench-scale, the deposition of manganese oxides, a permanganate reaction byproduct, has been shown to cause pore plugging and the formation of a manganese oxide layer adjacent to the non-aqueous phase liquid creosote which reduces post-treatment mass transfer and hence mass loading from the source. The objective of this study was to investigate the potential of partial permanganate treatment to reduce the ability of a coal tar creosote source zone to generate a multi-component plume at the pilot-scale over both the short-term (weeks to months) and the long-term (years) at a site where there is >10 years of comprehensive synoptic plume baseline data available. A series of preliminary bench-scale experiments were conducted to support this pilot-scale investigation. The results from the bench-scale experiments indicated that if sufficient mass removal of the reactive compounds is achieved then the effective solubility, aqueous concentration and rate of mass removal of the more abundant non-reactive coal tar creosote compounds such as biphenyl and dibenzofuran can be increased. Manganese oxide formation and deposition caused an order-of-magnitude decrease in hydraulic conductivity. Approximately 125 kg of permanganate were delivered into the pilot-scale source zone over 35 days, and based on mass balance estimates <10% of the initial reactive coal tar creosote mass in the source zone was oxidized. Mass discharge estimated at a down-gradient fence line indicated >35% reduction for all monitored compounds except for biphenyl, dibenzofuran and fluoranthene 150 days after treatment, which is consistent with the bench-scale experimental results. Pre- and post-treatment soil core data indicated a highly variable and random spatial distribution of mass within the source zone and provided no insight into the mass removed of any of the monitored species. The down-gradient plume was monitored approximately 1, 2 and 4 years following treatment. The data collected at 1 and 2 years post-treatment showed a decrease in mass discharge (10 to 60%) and/or total plume mass (0 to 55%); however, by 4 years post-treatment there was a rebound in both mass discharge and total plume mass for all monitored compounds to pre-treatment values or higher. The variability of the data collected was too large to resolve subtle changes in plume morphology, particularly near the source zone, that would provide insight into the impact of the formation and deposition of manganese oxides that occurred during treatment on mass transfer and/or flow by-passing. Overall, the results from this pilot-scale investigation indicate that there was a significant but short-term (months) reduction of mass emanating from the source zone as a result of permanganate treatment but there was no long-term (years) impact on the ability of this coal tar creosote source zone to generate a multi-component plume.     

Natural attenuation of a plume from an emplaced coal tar creosote source over 14 years

An emplaced source of coal tar creosote within the sandy Borden research aquifer has documented the long-term (5140 days) natural attenuation for this complex mixture. Plumes of dissolved chemicals were produced by the essentially horizontal groundwater flowing at about 9 cm/day. Eleven chemicals have been extensively sampled seven times using a monitoring network of approximately 280, 14-point multilevel samplers. A model of source dissolution using Raoult's Law adequately predicted the dissolution of 9 of 11 compounds. Mass transformation has limited the extent of the plumes as groundwater has flowed more than 500 m, yet the plumes are no longer than 50 m. Phenol and xylenes have been removed and naphthalene has attenuated from its maximum extent on day 1357. Some compound plumes have reached an apparent steady state and the plumes of other compounds (dibenzofuran and phenanthrene) are expected to continue to expand due to an increasing mass flux and limited degradation potential. Biotransformation is the major process controlling natural attenuation at the site. The greatest organic mass lost is associated with the high solubility compounds. However, the majority of the mass loss for most compounds has occurred in the source zone. Oxygen is the main electron acceptor, yet the amount of organics lost cannot be accounted for by aerobic mineralization or partial mineralization alone. The complex evolution of these plumes has been well documented but understanding the controlling biotransformation processes is still elusive. This study has shown that anticipating bioattenuation patterns should only be considered at the broadest scale. Generally, the greatest mass loss is associated with those compounds that have a high solubility and low partitioning coefficients.     

A robust approach for iterative contaminant source location and release history recovery

Contamination source identification is a crucial step in environmental remediation. The exact contaminant source locations and release histories are often unknown due to lack of records and therefore must be identified through inversion. Coupled source location and release history identification is a complex nonlinear optimization problem. Existing strategies for contaminant source identification have important practical limitations. In many studies, analytical solutions for point sources are used; the problem is often formulated and solved via nonlinear optimization; and model uncertainty is seldom considered. In practice, model uncertainty can be significant because of the uncertainty in model structure and parameters, and the error in numerical solutions. An inaccurate model can lead to erroneous inversion of contaminant sources. In this work, a constrained robust least squares (CRLS) estimator is combined with a branch-and-bound global optimization solver for iteratively identifying source release histories and source locations. CRLS is used for source release history recovery and the global optimization solver is used for location search. CRLS is a robust estimator that was developed to incorporate directly a modeler's prior knowledge of model uncertainty and measurement error. The robustness of CRLS is essential for systems that are ill-conditioned. Because of this decoupling, the total solution time can be reduced significantly. Our numerical experiments show that the combination of CRLS with the global optimization solver achieved better performance than the combination of a non-robust estimator, i.e., the nonnegative least squares (NNLS) method, with the same solver.     

Biogeochemical characterisation of a coal tar distillate plume

The distillation of acidified coal tars for up to 50 years has given rise to a phenol plume approximately 500 m long, 50 m deep and containing up to 15 g l(-1) dissolved organic carbon (DOC) in the Triassic Sandstones aquifer. A conceptual biogeochemical model based on chemical and microbiological analysis of groundwater samples has been developed as a preliminary to more detailed studies of the controls on natural attenuation. While the development of redox zones and the production of methane and carbon dioxide provide evidence of natural attenuation, it appears that degradation is slow. The existence of sulphate in the plume indicates that this electron acceptor has not been depleted and that consequently methanogenesis is probably limited. Based on a simple estimate of sulphate input concentration, a half-life of about 15 years has been estimated for sulphate reduction. Geochemical modelling predicts that increased alkalinity within the plume has not led to carbonate precipitation, and thus within the limits of accuracy of the measurement, alkalinity may reflect the degree of biodegradation. This implies a loss of around 18% of the DOC over a 30-year period. Despite limited degradation, microbial studies show that there are diverse microbial communities in the aquifer with the potential for both anaerobic and aerobic biodegradation. Microbial activity was found to be greatest at the leading edge of the plume where DOC concentrations are 60 mg l(-1) or less, but activity could still be observed in more contaminated samples even though cells could not be cultured. The study suggests that degradation may be limited by the high phenol concentrations within the core of the plume, but that once diluted by dispersion, natural attenuation may proceed. More detailed studies to confirm these initial findings are identified and form the basis of associated papers.     

Microcosm studies of microbial degradation in a coal tar distillate plume

Investigation of a groundwater plume containing up to 24 g l(-1) phenolic compounds suggested that over a period of nearly 50 years, little degradation had occurred despite the presence of a microbial community and electron acceptors within the core of the plume. In order to study the effect of contaminant concentration on degradation behaviour, laboratory microcosm experiments were performed under aerobic and anaerobic conditions at four different concentrations obtained by diluting contaminated with uncontaminated groundwater. The microcosms contained groundwater with total phenols at ca. 200, 250, 660 and 5000 mg l(-1), and aquifer sediment that had been acclimatised within the plume for several months. The microcosms were operated for a period of 390-400 days along with sterile controls to ascertain whether degradation was microbially mediated or abiotic. Under aerobic conditions, degradation only occurred at concentrations up to 660 mg l(-1) total phenols. At phenol concentrations below 250 mg l(-1) a benzoquinone intermediate, thought to originate from the degradation of 2,5-dimethylphenol, was isolated and identified. This suggested an unusual degradative pathway for this compound; its aerobic degradation more commonly proceeding via catecholic intermediates. Under anaerobic conditions, degradation only occurred in the most dilute microcosm (total phenols 195 mg l(-1)) with a loss of p-cresol accompanied by a nonstoichiometric decrease in nitrate and sulphate. By inference, iron(III) from the sediment may also have been used as a terminal electron acceptor, in which case the amount of biologically available iron released was calculated as 1.07 mg Fe(III)/g of sediment. The study shows that natural attenuation is likely to be stimulated by dilution of the plume.     

Temporal evolution of DNAPL source and contaminant flux distribution: impacts of source mass depletion

We investigated, using model simulations, the changes occurring in the distribution of dense non-aqueous phase liquid (DNAPL) mass (Sn) within the source zone during depletion through dissolution, and the resulting changes in the contaminant flux distribution (J) at the source control plane (CP). Two numerical codes (ISCO3D and T2VOC) were used to simulate selected scenarios of DNAPL dissolution and transport in three-dimensional, heterogeneous, spatially correlated, random permeability fields with emplaced sources. Data from the model simulations were interpreted based on population statistics (mean, standard deviation, coefficient of variation) and spatial statistics (centroid, second moments, variograms). The mean and standard deviation of the Sn and J distributions decreased with source mass depletion by dissolution. The decrease in mean and standard deviation was proportional for the J distribution resulting in a constant coefficient of variation (CV), while for the Sn distribution, the mean decreased faster than the standard deviation. The spatial distributions exhibited similar behavior as the population distribution, i.e., the CP flux distribution was more stable (defined by temporally constant second moments and range of variograms) than the Sn distribution. These observations appeared to be independent of the heterogeneity of the permeability (k) field (variance of the log permeability field=1 and 2.45), correlation structure (positive vs. negative correlation between the k and Sn domains) and the DNAPL dissolution model (equilibrium vs. rate-limited), for the cases studied. Analysis of data from a flux monitoring field study (Hill Air Force Base, Utah) at a DNAPL source CP before and after source remediation also revealed temporal invariance of the contaminant flux distribution. These modeling and field observations suggest that the temporal evolution of the contaminant flux distribution can be estimated if the initial distribution is known. However, the findings are preliminary and broader implications to sampling strategies for remediation performance assessment need to be evaluated in additional modeling and experimental studies.     

Evolution of bacterial community during bioremediation of PAHs in a coal tar contaminated soil

The monitoring of a windrow treatment applied to soil contaminated by mostly 2-, 3- and 4-ring PAHs produced by coal tar distillation was performed by following the evolution of both PAH concentration and the bacterial community. Total and PAH-degrading bacterial community structures were followed by 16S rRNA PCR-DGGE in parallel with quantification by bacterial counts and 16 PAH measurements. Six months of biological treatment led to a strong decrease in 2-, 3- and 4-ring PAH concentrations (98, 97 and 82% respectively). This result was associated with the activity of bacterial PAH-degraders belonging mainly to the Gamma-proteobacteria, in particular, the Enterobacteria and Pseudomonas genera, which were detected over the course of the treatment. This group was considered to be a good bioindicator to determine the potential PAH biodegradation of contaminated soil. Conversely, other species, like the Beta-proteobacteria, were detected after 3months, when 2-, 3- and 4-ring PAHs were almost completely degraded. Thus, presence of the Beta-proteobacteria group could be considered a good candidate indicator to estimate the endpoint of biotreatment of this type of PAH-contaminated soil.     

A laboratory study of sediment and contaminant release during gas ebullition

Significant quantities of gas are generated from labile organic matter in contaminated sediments. The implications for the gas generation and subsequent release of contaminants from sediments are unknown but may include enhanced direct transport such as pore water advection and diffusion. The behavior of gas in sediments and the resulting migration of a polyaromatic hydrocarbon, viz phenanthrene, were investigated in an experimental system with methane injection at the base of a sediment column. Hexane above the overlying water layer was used to trap any phenanthrene migrating out of the sediment layer. The rate of suspension of solid particulate matter from the sediment bed into the overlying water layer was also monitored. The experiments indicated that significant amounts of both solid particulate matter and contaminant can be released from a sediment bed by gas movement with the amount of release related to the volume of gas released. The effective mass transfer coefficient of gas bubble-facilitated contaminant release was estimated under field conditions, being around three orders of magnitude smaller than that of bioturbation. A thin sand-capping layer (2 cm) was found to dramatically reduce the amount of contaminant or particles released with the gas because it could prevent or at least reduce sediment suspension. Based on the experimental observations, gas bubble-facilitated contaminant transport pathways for both uncapped and capped systems were proposed. Sediment cores were sliced to obtain phenanthrene concentration. X-ray computed tomography (CT) was used to investigate the void space distribution in the sediment penetrated by gas bubbles. The results showed that gas bubble migration could redistribute the sediment void spaces and may facilitate pore water circulation in the sediment.     

Sampling and analysis of particulate and gaseous polycyclic aromatic hydrocarbons from coal tar sources in the working environment

A sampling method for polycyclic aromatic hydrocarbons (PAH) in workplace atmospheres using a standard glass fiber filter with a back-up section of Amberlite XAD-2 is described. Filters and XAD are solvent-desorbed and, without further clean-up, the solutions are submitted to high performance liquid chromatographic analysis using fluorimetric detection. Recovery of PAH from XAD-2 was studied, and was found to be in the range 80–100%. The sampling method was evaluated in a coke plant, an aluminium plant and in a creosote impregnating plant. The method described is rapid and sensitive for the determination of the most important workplace PAH, and accounts for both particulate and gaseous PAH.     

Migration and natural fate of a coal tar creosote plume: 1. Overview and plume development

A volume of sand containing coal tar creosote was emplaced below the water table at CFB Borden to investigate natural attenuation processes for complex biodegradable mixtures. Coal tar creosote is a mixture of more than 200 polycyclic aromatic hydrocarbons, heterocyclic compounds and phenolic compounds. A representative group of seven compounds was selected for detailed study: phenol, m-xylene, naphthalene, phenanthrene, 1-methylnaphthalene, dibenzofuran and carbazole. Movement of groundwater through the source led to the development of a dissolved organic plume, which was studied over a 4-year period. Qualitative plume observations and mass balance calculations indicated two key conclusions: (1) compounds from the same source can display distinctly different patterns of plume development and (2) mass transformation was a major influence on plume behaviour for all observed compounds.     

Using polymer solutions to enhance recovery of mobile coal tar and creosote DNAPLs

Direct pumping and enhanced recovery of coal tar and creosote dense, non-aqueous phase liquids (DNAPLs) from the subsurface have had mixed results because these DNAPLs are viscous fluids that can potentially alter aquifer wettability. To improve the inefficiencies associated with waterflooding, the research presented here considered the use of a polymer solution that can be added to the injected flood solution to increase the viscosity and decrease the velocity of the flooding solution. Results from one-dimensional, vertically oriented laboratory column experiments that evaluate the recovery of coal-derived DNAPL with both water and polymer flooding solutions are presented. The final DNAPL saturation remaining in the column was assessed in water and oil-wet systems for three viscous DNAPLs. Adding polymer to increase the aqueous solution viscosity did not have a significant impact in water-wet systems. A final DNAPL saturation of approximately 19% was achieved for both water and polymer floods. In contrast, the addition of polymer significantly improved recovery in oil-wet systems. The final saturation was over 40% in oil-wet systems after waterflooding, but approximately 19% with a polymer flushing solution. Although the final saturation produced with polymer flooding was similar between the oil- and water-wet systems, differences in the relative permeability and distribution of DNAPL in the porous matrix caused the DNAPL recovery to be much slower in the oil-wet system.     

扬州古运河沉积物污染物释放强度与特征研究

研究了扬州古运河沉积物中氮、磷、铁和锰等污染物释放强度与特征,着重考察了上覆水中DO和污染物负荷对污染物释放强度的影响,分析了氮、磷、铁和锰等污染物释放条件.结果表明:扬州古运河有一定的自净能力,且呈现初期释放NH4+ -N、后期转化为NO3-- N的特征;随着水环境条件改变,沉积物NH4+ -N、PO34-、铁和锰平均释放强度变化范围是123.14～445.68、13.50～32.45、25.20～198.47、19.97～25.02 mg/(m2·d);上覆水DO浓度越低,污染物平均释放强度越大;上覆水污染物浓度越小,污染物(除铁外)平均释放强度越大;在DO＜5 mg/L条件下,沉积物向上覆水释放NH4+ -N;在DO＜1 mg/L条件下,沉积物向上覆水释放PO34-、铁和锰;在DO为2～5 mg/L条件下,上覆水PO34-和铁向沉积物中汇入,但锰则呈现释放和汇入的交替变化.     

Spatial and temporal features of density-dependent contaminant transport: Experimental investigation and numerical modeling

We investigate the spatial and temporal features of variable-density contaminant plumes migration in porous materials. Our analysis is supported by novel experimental results concerning concentration profiles inside a vertical column setup that has been conceived at CEA to this aim. The experimental method relies on X-ray spectrometry, which allows determining solute profiles as a function of time at several positions along the column. The salient outcomes of the measurements are elucidated, with focus on miscible fluids in homogeneous saturated media. The role of the injected solution molarity is evidenced. As molarity increases, the solutes plume transport progressively deviates from the usual Fickian behavior, and pollutants distribution becomes skewed in the direction dictated by gravity. By resorting to a finite elements approach, we numerically solve the nonlinear equations that rule the pollutants migration: a good agreement is found between the simulated profiles and the experimental data. At high molarity, a strong dependence on initial conditions is found. Finally, we qualitatively explore the (unstable) interfacial dynamics between the dense contaminant plume and the lighter resident fluid that saturates the column, and detail its evolution for finite-duration contaminant injections.     

Influence of mass transfer characteristics for DNAPL source depletion and contaminant flux in a highly characterized glaciofluvial aquifer

The transfer of contaminant mass between the nonaqueous- and aqueous-phases is a process of central importance for the remediation of sites contaminated by dense nonaqueous-phase liquids (DNAPLs). This paper describes a comparison of the results obtained with various alternative DNAPL-aqueous-phase mass transfer models contained in the literature for predicting DNAPL source-zone depletion times in groundwater systems. These dissolution models were largely developed through laboratory column experiments. To gain insight into the implications of various representations of the local-scale kinetic as well as equilibrium DNAPL dissolution processes, aquifer heterogeneity and the complex architecture of a DNAPL source-zone, the aqueous-phase contaminant concentrations and mass fluxes arriving at a down-gradient compliance boundary are analyzed in a conditional stochastic framework. The hydrogeologic setting is a heterogeneous fluvial aquifer in Southwest Germany, referred to as the aquifer analog dataset, that was intensively characterized in three dimensions for hydrogeological parameters that include permeability, effective porosity, grain size, mineralogy and sorption coefficients. By embedding the various dissolution models into the compositional, multiphase flow model, CompFlow, the relative times predicted for complete depletion of a released DNAPL source due to natural dissolution are explored. Issues related to achieving environmental benefits through, for example, partial DNAPL-zone source removal via enhanced remedial technologies are also discussed. In this context, performance metrics in the form of peak aqueous-phase contaminant concentrations and mass fluxes arriving at a down-gradient compliance boundary are compared to each other. This is done for each of the alternative mass transfer models. A significant reduction in the fractional flux at a downstream location from the DNAPL source can be achieved by partial source-zone mass reduction; however, peak concentration levels at the same location remain much higher than the United States Environment Protection Agency (US-EPA) drinking water limits. Although groundwater quality was found to improve more rapidly for the equilibrium dissolution model, it is also shown that dissolution models that promote rapid DNAPL disappearance produce greater prediction uncertainty in the aqueous-phase flux reduction.     

Geochemical modeling of arsenic release from a deep natural solid matrix under alternated redox conditions

Dissolved arsenic (As) concentrations detected in groundwater bodies of the Emilia-Romagna Region (Italy) exhibit values which are above the regulation limit and could be related to the natural composition of the host porous matrix. To support this hypothesis, we present the results of a geochemical modeling study reproducing the main trends of the dynamics of As, Fe, and Mn concentrations as well as redox potential and pH observed during batch tests performed under alternating redox conditions. The tests were performed on a natural matrix extracted from a deep aquifer located in the Emilia-Romagna Region (Italy). The solid phases implemented in the model were selected from the results of selective sequential extractions performed on the tested matrix. The calibrated model showed that large As concentrations have to be expected in the solution for low crystallinity phases subject to dissolution. The role of Mn oxides on As concentration dynamics appears significant in strongly reducing environments, particularly for large water–solid matrix interaction times. Modeled data evidenced that As is released firstly from the outer surface of Fe oxihydroxides minerals exhibiting large concentrations in water when persistent reducing conditions trigger the dissolution of the crystalline structure of the binding minerals. The presence of organic matter was found to strongly affect pH and redox conditions, thus influencing As mobility.

Assessing chlorinated ethene degradation in a large scale contaminant plume by dual carbon-chlorine isotope analysis and quantitative PCR

The fate of chlorinated ethenes in a large contaminant plume originating from a tetrachloroethene (PCE) source in a sandy aquifer in Denmark was investigated using novel methods including compound-specific carbon and chlorine isotope analysis and quantitative real-time polymerase chain reaction (qPCR) methods targeting Dehaloccocoides sp. and vcrA genes. Redox conditions were characterized as well based on concentrations of dissolved redox sensitive compounds and sulfur isotopes in SO(4)(2-). In the first 400 m downgradient of the source, the plume was confined to the upper 20 m of the aquifer. Further downgradient it widened in vertical direction due to diverging groundwater flow reaching a depth of up to 50 m. As the plume dipped downward and moved away from the source, O(2) and NO(3)(-) decreased to below detection levels, while dissolved Fe(2+) and SO(4)(2-) increased above detectable concentrations, likely due to pyrite oxidation as confirmed by the depleted sulfur isotope signature of SO(4)(2-). In the same zone, PCE and trichloroethene (TCE) disappeared and cis-1,2-dichloroethene (cDCE) became the dominant chlorinated ethene. PCE and TCE were likely transformed by reductive dechlorination rather than abiotic reduction by pyrite as indicated by the formation of cDCE and stable carbon isotope data. TCE and cDCE showed carbon isotope trends typical for reductive dechlorination with an initial depletion of (13)C in the daughter products followed by an enrichment of (13)C as degradation proceeded. At 1000 m downgradient of the source, cDCE was the dominant chlorinated ethene and had reached the source δ(13)C value confirming that cDCE was not affected by abiotic or biotic degradation. Further downgradient (up to 1900 m), cDCE became enriched in (13)C by up to 8 ‰ demonstrating its further transformation while vinylchloride (VC) concentrations remained low (<1 μg/L) and ethene was not observed. The correlated shift of carbon and chlorine isotope ratios of cDCE by 8 and 3.9 ‰, respectively, the detection of Dehaloccocides sp genes, and strongly reducing conditions in this zone provide strong evidence for reductive dechlorination of cDCE. The significant enrichment of (13)C in VC indicates that VC was transformed further, although the mechanism could not be determined. The transformation of cDCE was the rate limiting step as no accumulation of VC occurred. In summary, the study demonstrates that carbon-chlorine isotope analysis and qPCR combined with traditional approaches can be used to gain detailed insight into the processes that control the fate of chlorinated ethenes in large scale plumes.     

A PCE groundwater plume discharging to a river: influence of the streambed and near-river zone on contaminant distributions

An investigation of a tetrachloroethene (PCE) groundwater plume originating at a dry cleaning facility on a sand aquifer and discharging to a river showed that the near-river zone strongly modified the distribution, concentration, and composition of the plume prior to discharging into the surface water. The plume, streambed concentration, and hydrogeology were extensively characterized using the Waterloo profiler, mini-profiler, conventional and driveable multilevel samplers (MLS), Ground Penetrating Radar (GPR) surveys, streambed temperature mapping (to identify discharge zones), drivepoint piezometers, and soil coring and testing. The plume observed in the shallow streambed deposits was significantly different from what would have been predicted based on the characteristics of the upgradient plume. Spatial and temporal variations in the plume entering the near-river zone contributed to the complex contaminant distribution observed in the streambed where concentrations varied by factors of 100 to 5000 over lateral distances of less than 1 to 3.5 m. Low hydraulic conductivity semi-confining deposits and geological heterogeneities at depth below the streambed controlled the pattern of groundwater discharge through the streambed and influenced where the plume discharged into the river (even causing the plume to spread out over the full width of the streambed at some locations). The most important effect of the near-river zone on the plume was the extensive anaerobic biodegradation that occurred in the top 2.5 m of the streambed, even though essentially no biodegradation of the PCE plume was observed in the upgradient aquifer. Approximately 54% of the area of the plume in the streambed consisted solely of PCE transformation products, primarily cis-1,2-dichloroethene (cDCE) and vinyl chloride (VC). High concentrations in the interstitial water of the streambed did not correspond to high groundwater-discharge zones, but instead occurred in low discharge zones and are likely sorbed or retarded remnants of past high-concentration plume discharges. The high-concentration areas (up to 5529 microg/l of total volatile organics) in the streambed are of ecological concern and represent potential adverse exposure locations for benthic and hyporheic zone aquatic life, but the effect of these exposures on the overall health of the river has yet to be determined. Even if the upgradient source of PCE is remediated and additional PCE is prevented from reaching the streambed, the high-concentration deposits in the streambed will likely take decades to hundreds of years to flush completely clean under natural conditions because these areas have low vertical groundwater flow velocities and high retardation factors. Despite high concentrations of contaminants in the streambed, PCE was detected in the surface water only rarely due to rapid dilution in the river and no cDCE or VC was detected. Neither the sampling of surface water nor the sampling of the groundwater from the aquifer immediately adjacent to the river gave an accurate indication of the high concentrations of PCE biodegradation products present in the streambed. Sampling of the interstitial water of the shallow streambed deposits is necessary to accurately characterize the nature of plumes discharging to rivers.     

Effectiveness of nanoscale zero-valent iron for treatment of a PCE-DNAPL source zone

Nanoscale zero-valent iron (nZVI) has received considerable attention as a potential in situ remediation technology for treating chlorinated solvent source zones. Experimental and mathematical modeling studies were conducted to investigate the performance of nZVI in the transformation of tetrachloroethene (PCE) entrapped as a dense nonaqueous phase liquid (DNAPL). Injection of a 60 g/L suspension of nZVI into a column containing 20-30 mesh Ottawa sand and PCE-DNAPL at a residual saturation of 5.5% resulted in a uniform distribution of nZVI and minimal displacement of PCE. Subsequent flushing with 267 pore volumes of water containing 3mM CaCl(2) at a Darcy velocity of 0.75 m/day resulted in steady-state effluent concentrations of PCE near the solubility limit (ca. 200mg/L) and production of dissolved-phase ethene (10-30 mg/L). Over the duration of the experiment, approximately 30% of the initial PCE-DNAPL mass reacted to form ethene, 50% was eluted as dissolved-phase PCE, and 20% remained in the column as PCE-DNAPL. To further explore the implications of the nZVI column results, a multiphase transport model was developed that incorporated rate-limited PCE-DNAPL dissolution and reactions with nZVI. Using a fitted pseudo first-order transformation rate coefficient of 1.421/h, the model accurately captured observed trends in effluent concentrations of PCE and ethene and overall mass balance. A model sensitivity study reveals a strong dependence of treatment effectiveness on system characteristics. The sensitivity analysis suggests that an increase in the extent of PCE transformation is facilitated by decreasing flow rate, emplacement of nZVI down-gradient of the DNAPL source zone, and decreasing length of the DNAPL source zone. These findings indicate that, although emplacement of high concentrations of nZVI within a PCE-DNAPL source zone can result in substantial transformation of the parent compound, careful attention to design parameters (e.g. flow rate, location and amount nZVI delivered) will be required to achieve complete conversion to benign reaction products.     

A comparison of contaminant dynamics in arctic and temperate fish: A modeling approach

In order to compare the abilities of arctic and temperate fish to accumulate PCBs we conduct a metabolic analysis to determine how process rates in a mathematical fish contaminant model change with temperature. We evaluate the model by applying the original and adapted models to estimate PCB concentrations in lake trout (Salvelinus namaycush) in Trout Lake, Ontario, Canada, and in arctic char (Salvelinus alphinus) in Lake ?yangen, in the Norwegian high arctic. Modeled concentrations are, for the most part, within 50% of mean measured values and are comparable to the error associated with the fish data. In order to evaluate differences in fish bioaccumulation processes, the model is applied to hypothetical arctic and temperate systems, assuming the same contaminant input values in water and diet. The model predicts that temperate salmonids are able to biomagnify PCBs 6-60% more than arctic salmonids. For all congeners, the lower BMF(MAX) of arctic fish contribute to their lower concentrations. For congeners with log K(ow) < 6.0, the lower concentrations in arctic fish are also attributed to faster loss due to gill ventilation. Faster growth rates for temperate fish reduce the difference in bioaccumulation for congeners with log K(ow) > 7.0. These processes are controlled by the influence of lipid in the fish and their diet as well as the dependence of growth on temperature. We suggest that fish models originally calibrated for temperate systems may be directly applied to arctic lakes after accounting for the lipid content of the fish and their diet as well as water temperature.