Remediation Endpoints Research Articles

A screening model for predicting the potential of soil colloids-enhanced leaching of hydrophobic organic contaminants to groundwater at contaminated sites

Modeling the fate and transport of organic pollutants at contaminated sites is critical for risk assessment and management practices, such as establishing realistic cleanup standards or remediation endpoints. Against the conventional wisdom that highly hydrophobic persistent organic pollutants (POPs) (e.g., polybrominated diphenyl ethers and polycyclic aromatic hydrocarbons) in surface soils are essentially immobile, mounting evidence has demonstrated the potential of these contaminants leaching into the groundwater, due to enhanced transport by soil colloids. Here, we develop a Colloids-Enhanced Transport (CET) model, which can be used as a simple screening tool to predict the leaching potential of POPs into groundwater, as mediated by soil colloids. The CET model incorporates several processes, including the release of POPs-bearing colloids into the porewater, the vertical transport of colloids and associated POPs in the vadose zone, the mixing of POPs-containing soil leachate with groundwater, and the migration of POPs-bearing colloids in saturated zone. Thus, using parameters that can be easily obtained (e.g., annual rainfall, soil type, and common hydrogeological properties of the subsurface porous media), the CET model can estimate the concentrations of POPs in the saturated zone from the observed POPs concentrations in surface or shallow subsurface zones. The CET model can also be used to derive soil quality standards or cleanup endpoints by back-calculating soil concentrations based on groundwater protection limits.

Bioengineering remediation of former industrial sites contaminated with chemical mixtures

Former gasworks sites are known to be contaminated with complex chemical mixtures that require remediation before redevelopment. Bioamendments such as biochar and spent mushroom compost (SMC) offer a green and sustainable remediation approach to help tackle this issue. However, the effectiveness of different biochar types and their interactions with the soil microbial community is still not well understood. To address this, a full factorial microcosm experiment was carried out using biochar derived from rice husk (RHB) and wheat straw (WSB) mixed with soil from a former gasworks site at varying concentrations (0%, 2.5%, and 5%), with and without SMC. The experiment aimed to evaluate the fate of contaminants including alkanes, PAHs, and metals, and their effect on the soil microbial community, as well as the implications for remediation endpoints. The results showed that the bioamendments had an average TPH reduction of 92%, with SMC and WSB-SMC having the highest degradation rates at 93%. While the bioamendments did not significantly affect the extent of TPH removal compared to the control, they did improve the degradation of high molecular weight (HMW) PAHs, particularly in RHB-SMC for EC17-20 (60%) and EC21-35 (62%) of total PAH concentration, and in WSB-SMC for HMW bioavailable PAH concentration (89%). The bioamendments also affected the partitioning and distribution of metals after 120 days of treatment, leading to decreased available phase fractions. The treatments increased microbial abundance in the soil, with Gram positives, Gram negatives, and fungi increasing by 4%, 8%, and 38%, respectively, after 120 days, particularly in SMC and mixed treatments. This was mirrored in increased microbial soil respiration. After 120 days, low metal (178 ± 5 mg/kg) and TPH (21 ± 7 mg/kg) bioavailability translated into higher EC50 (10624 ± 710 mg/L), indicating lower toxicity. There was a strong correlation between bioavailability and toxicity of TPH and metals with microbial relative abundance and activity. In summary, while green and sustainable remediation may accelerate the remediation process, monitored natural attenuation may be sufficient for site reclamation. However, this strategy, as demonstrated here, can reduce metal bioavailability, and promote the biodegradation of HMW PAHs.

In situ (bio)remediation treatment options for U and Sr contaminated land: a comparison of radionuclide retention and remobilisation

This is the first comprehensive study to compare the long term stability of U and Sr (bio)remediation end-points following oxidative remobilisation (via oxygen/nitrate) yielding insights into retention capabilities and removal mechanisms.

Effects of thermal desorption on ecotoxicological characteristics of heavy petroleum-contaminated soil.

This study comprehensively evaluates the ecotoxicity of high-concentration heavy petroleum (HCHP)-contaminated soil before and after thermal desorption (TD) remediation at different temperatures and times. The results showed that the detoxification of contaminated soil was effectively achieved by extending the remediation duration at 400-600 °C. After treatment at 400 °C for 60 min, the toxicological indicators including bioluminescence EC50 (acute toxicity), seed germination ratio (Gr) and plant biomass of Brassica juncea (subacute toxicity), and diversity of the microbial community (chronic toxicity) reached a maximum. The value of the SOS-Induction Factor (SOSIF), characterizing genotoxicity was below 1.5, indicating that it was non-toxic. Pearson's correlation analysis illustrated that the water-soluble fraction (WSF), ALK1-3 and ARO1-3 of petroleum hydrocarbons were the primary sources of ecotoxicity. Notably, although the total ratio of petroleum removed from the soil reached 87.26 ± 4.38 %-98.69 ± 1.61 % under high-temperature thermal desorption (HTTD, 500-600 °C), the ecotoxicity was not lower than that at 400 °C. The pyrolysis products of petroleum macromolecules and extreme changes in soil properties were the leading causes of soil ecotoxicity following HTTD. The inconsistency between the removal of petroleum pollutants and ecological health risks reveals the significance of soil ecotoxicological assessments for identifying TD remediation endpoints and process optimization.

Phytotoxicity complements chemical assessment for re-use and re-purposing of refinery wastes for soil amendment purposes after bioremediation

The objective of this study was to assess the suitability of onsite re-use of mature compost for landscaping and tree mulching, produced from the bioremediation of oily sludge from the refinery. Compost samples from the co-composting process were analysed for a range of contaminants, including a human health risk assessment fractionation (HRAF) of the remaining petroleum hydrocarbons, as well as a phytotoxicity test. The chemical characterisation demonstrated that the process removed more than 94% of the original petroleum hydrocarbons from the sludge, and the removal rates were high at 1155 mg/kg/day. The HRAF demonstrated no residual risks, posed by the petroleum hydrocarbons present in the compost to human health if used on-site when compared to the relevant Australian environmental investigation levels (ILs). However, the phytotoxicity assessment demonstrated that the compost was toxic to germinating lettuce. The gap in the literature this study addressed was to provide an estimate of the LD50 and no effect concentration (NEC) for the compost using a standard plant bioassay containing a range of residual (aged) and bioremediated refinery process wastes, including petroleum hydrocarbons. The values estimated for LD50 and NEC were approximately 125 and 43 mg/kg, respectively for compost containing residual petroleum hydrocarbon fractions, filling a gap in the current literature which has limited data on standard toxicity values that can be used in determining and informing commercial remediation strategies and their outcomes with aged sludges. Phytotoxicity was shown to be an important complement to conventional analyses and HRAF data when characterising the sludge, and understanding its potential for re-use. The novelty of the study is that it highlighted a gap in the complementary use of chemical and bioassay analyses for evaluating refinery waste remediation endpoints, which has potential for broader application to other projects.

Toward Optimizing LNAPL Remediation

AbstractOptimizing the remediation of light nonaqueous phase liquids (LNAPLs) to achieve an acceptable endpoint status for a site is not trivial. Recently, Sookhak Lari, Johnston, et al. (2018, https://doi.org/10.1016/j.jhazmat.2017.11.006), Sookhak Lari, Rayner, and Davis (2018, https://doi.org/10.1016/j.jenvman.2018.07.041) conducted three‐dimensional multiphase, multicomponent simulations to address LNAPL remediation endpoints for a single recovery well. However, optimized LNAPL remediation for multiple wells is not addressed in the literature. In the first part of this paper, we establish a matrix of 10 simulation scenarios to show the sensitivity of the remedial endpoint (i.e., what is feasibly achieved) to several parameters including viscosity and partitioning attributes of the LNAPL, heterogeneity of the formation, and the location and number of the recovery wells. While this addresses the variability of LNAPL removal from the subsurface and is valuable in its own right, it does not address the optimal removal of LNAPL. We address this in the second part of the paper by linking a genetic algorithm to TMVOC‐MP to allow, for example, the assessment of the optimal number and location of LNAPL recovery wells in a field‐scale problem. Using supercomputing facilities and within 49 genetic algorithm generations, each including 150 members, highly optimized answers to different objective functions were obtained. For the first time, such a multiphase, multicomponent optimization tool promises the possibility of optimizing LNAPL remediation at field scale to achieve practicable endpoint conditions, within computational affordability.

Kinetics of phytoremediation of petroleum hydrocarbon contaminated soil

Earthmaster and the University of Waterloo have successfully developed plant growth promoting rhizobacteria (PGPR) – Enhanced Phytoremediation Systems (PEPSystems™) which have been deployed across Canada for the treatment of soil contaminated with petroleum hydrocarbons (PHCs), including CCME fractions F2 and F3. A challenge with phytoremediation is to predict the length of time to remediate a site so that site owners will be inclined to use the technology. In previous field trials of PEPSystems, it was determined that PHC was mostly degraded by microbes in the rhizosphere, following first-order exponential decay kinetics. Using new PEPSystems data collected from multiple commercial remediation sites across Western Canada, the kinetic equations of PHC decay were tested to determine if remediation time was accurately predicted. In general, when compared to the predicted time to remediation endpoint, data from recent commercial field applications showed that 35% and 20% less time was needed to reach remediation endpoints for fractions F2 and F3, respectively. As a result, the predictive kinetic equation for fraction F2 degradation was updated to reflect current remediation outcomes. Insufficient data were available to update the F3 equation. Being able to more accurately predict remediation timelines will enhance the value and utilization of PEPSystems.

Towards characterizing LNAPL remediation endpoints

Remediating sites contaminated with light non-aqueous phase liquids (LNAPLs) is a demanding and often prolonged task. It is vital to determine when it is appropriate to cease engineered remedial efforts based on the long-term effectiveness of remediation technology options. For the first time, the long term effectiveness of a range of LNAPL remediation approaches including skimming and vacuum-enhanced skimming each with and without water table drawdown was simulated through a multi-phase and multi-component approach. LNAPL components of gasoline were simulated to show how component changes affect the LNAPL's multi-phase behaviour and to inform the risk profile of the LNAPL. The four remediation approaches along with five types of soils, two states of the LNAPL specific mass and finite and infinite LNAPL plumes resulted in 80 simulation scenarios. Effective conservative mass removal endpoints for all the simulations were determined. As a key driver of risk, the persistence and mass removal of benzene was investigated across the scenarios. The time to effectively achieve a technology endpoint varied from 2 to 6 years. The recovered LNAPL in the liquid phase varied from 5% to 53% of the initial mass. The recovered LNAPL mass as extracted vapour was also quantified. Additional mass loss through induced biodegradation was not determined. Across numerous field conditions and release incidents, graphical outcomes provide conservative (i.e. more prolonged or greater mass recovery potential) LNAPL remediation endpoints for use in discussing the halting or continuance of engineered remedial efforts.

Field-scale multi-phase LNAPL remediation: Validating a new computational framework against sequential field pilot trials

Remediation of subsurface systems, including groundwater, soil and soil gas, contaminated with light non-aqueous phase liquids (LNAPLs) is challenging. Field-scale pilot trials of multi-phase remediation were undertaken at a site to determine the effectiveness of recovery options. Sequential LNAPL skimming and vacuum-enhanced skimming, with and without water table drawdown were trialled over 78days; in total extracting over 5m3 of LNAPL. For the first time, a multi-component simulation framework (including the multi-phase multi-component code TMVOC-MP and processing codes) was developed and applied to simulate the broad range of multi-phase remediation and recovery methods used in the field trials. This framework was validated against the sequential pilot trials by comparing predicted and measured LNAPL mass removal rates and compositional changes. The framework was tested on both a Cray supercomputer and a cluster. Simulations mimicked trends in LNAPL recovery rates (from 0.14 to 3mL/s) across all remediation techniques each operating over periods of 4–14days over the 78day trial. The code also approximated order of magnitude compositional changes of hazardous chemical concentrations in extracted gas during vacuum-enhanced recovery. The verified framework enables longer term prediction of the effectiveness of remediation approaches allowing better determination of remediation endpoints and long-term risks.

Characteristics of PAH tar oil contaminated soils—Black particles, resins and implications for treatment strategies

Tar oil contamination is a major environmental concern due to health impacts of polycyclic aromatic hydrocarbons (PAH) and the difficulty of reaching acceptable remediation end-points. Six tar oil-contaminated soils with different industrial histories were compared to investigate contamination characteristics by black particles. Here we provide a simple method tested on 6 soils to visualize and identify large amounts of black particles (BP) as either solid aggregates of resinified and weathered tar oil or various wood/coke/coal-like materials derived from the contamination history. These materials contain 2–10 times higher PAH concentrations than the average soil and were dominantly found in the sand fraction containing 42–86% of the total PAH. The PAH contamination in the different granulometric fractions was directly proportional to the respective total organic carbon content, since the PAH were associated to the carbonaceous particulate materials. Significantly lower (bio)availability of PAH associated to these carbonaceous phases is widely recognized, thus limiting the efficiency of remediation techniques. We provide a conceptual model of the limited mass transfer of PAH from resinated tar oil phases to the water phase and emphasize the options to physically separate BP based on their lower bulk density and slower settling velocity.

Mixed contaminant interactions in soil: Implications for bioavailability, risk assessment and remediation

Majority of contaminated sites in the world contain complex mixtures of heavy metals and organic contaminants from diverse natural processes and anthropogenic activities. Mixed interactions of heavy metals and organic contaminants may affect their bioavailability and accumulation in soil and biota through synergistic or antagonistic processes. Evaluation of contaminant bioavailability is a necessary component of the overall site assessment process for establishing either bioavailability-based or risk-based, site-specific remedial options. However, contemporary approaches aimed at the effective characterisation of contaminated soils for risk assessment, remedial and regulatory purposes are frequently challenged by knowledge gaps in contaminant bioavailability, mixed contaminant effects and emerging contaminants. Understanding mixed contaminant interactions at the elemental and molecular levels is, therefore, imperative not only to explain the underlying mechanisms controlling the fate and transport of these contaminants in soils, but also predict their bioavailability, ecotoxicological effects on natural communities under realistic exposure conditions and remediation endpoints. In this paper, scattered literature is harnessed to review specific soil-contaminant interactions, inter-contaminant (metal-metal, organic-organic, metal-organic) interactions and their implications for bioavailability, risk assessment and soil remediation. Key words: Heavy metals, organic contaminants, mixed contaminant interactions, co-contaminated soil, bioavailability, risk assessment, soil remediation.

Ozonation of diesel–fuel contaminated sand and the implications for remediation end-points

In this study, we investigate specifically the influence of soil grain size and water content on the degradation of n-alkane fractions and the formation of aldehydes and carboxylic acid during ozonation. 15g of quartz sand spiked with diesel (25gkg−1) were exposed to ozone for 20h at concentrations of 10, 30 and 50mgL−1, respectively. Results indicated that ozonation of the n-alkanes in fine grain size sand (0.15–0.25mm) was 1.2 times faster than coarse sand due to higher surface contact area between O3 and sand particles. Soil moisture below 18% w/w did not influence the ozonation efficiency. In contrast the ozonation led to an increase of acidity of the sand samples (pH=3.0) after 20h treatment. This was due to the formation of carboxylic acid. Formaldehyde, one of the key by-products of ozonation, was always <13mgkg−1 after the treatment which is below the industrial soil clean-up target level. While the aldehydes and carboxylic acid further reacted with O3 and their ozonation rate were slower than those of the alkanes suggesting that the hydroxylated by-products accumulated in the sand during the process. Overall the findings demonstrated that not only the alkanes but also aldehydes and carboxylic acid should be considered when defining remediation end-points.

Machine learning models for predicting PAHs bioavailability in compost amended soils

Compost addition to polluted soils is a strategy for waste reuse and soil remediation, while bioavailability is a key parameter for environmental assessment. Empirical data from an 8-month microcosm experiment were used to assess the ability and performance of six machine learning (ML) models to predict temporal bioavailability changes of 16 polycyclic aromatic hydrocarbons (PAHs) in contaminated soils amended with compost. The models included multilayer perceptrons (MLPs), radial basis function (RBF), support vector regression (SVR), M5 model tree (M5P), M5 rule (M5R) and linear regression (LR). Overall, the performance of the six models, determined by 10-fold cross validation method, was ranked as follows: RBF>M5P>SVR>MLP>M5R>LR. Results further demonstrated that the ML models successfully identified the relative importance of each variable (i.e. incubation time, organic carbon content, soil moisture content, nutrient levels) on the temporal bioavailability change of individual PAH. Such models can potentially be useful for predicting the concentration of a wide range of pollutants in soils, which could contribute to reduce chemical monitoring at site and help decision making for remediation end points and risk assessment.

Mineralogical study and leaching behavior of a stabilized harbor sediment with hydraulic binder

The environmental assessment of potential effects of contaminated harbor sediments stabilized with hydraulic binders and the determination of remediation endpoints require the determination of pollutants leaching potentials. Moreover, little information about the speciation and mobility of inorganic contaminants in these specific solid matrices is available in the literature. The objective of this paper is to investigate the relationship between mineralogy and leachability of contaminants (copper, lead, and zinc) present in a French harbor sediment stabilized with quicklime and Portland cement. Batch equilibrium leaching tests at various pH, chemical analysis of leachates, and mineralogical studies (X-ray diffraction, scanning electron microscopy-energy dispersive spectroscopy, and diffuse reflectance infrared Fourier transform) have been combined in the present investigation. The acid neutralization capacity of the stabilized matrix studied is first controlled by the dissolution of portlandite (pH ~12), followed by the dissolution of C-S-H (pH ~11) and the dissolution of ettringite (pH ~10). Finally, a very high buffering capacity of this stabilized sediment is observed for pH values around 6. This equilibrium is mainly controlled by the dissolution of iron sulfides and carbonate minerals. Consequently, the mobilization of inorganic contaminants as a function of pH remains very low (<0.1 wt%) for pH values above 6 and significantly increases for pH below these values. This research confirms the importance of a combined methodology for the intrinsic characterization of potential mobilization of contaminants in a stabilized sediment and for a better understanding of geochemical processes that affect contaminant fate, transformation, and transport in the subsurface environment.

A Bayesian approach to integrate temporal data into probabilistic risk analysis of monitored NAPL remediation

Upon their release into the subsurface, non-aqueous phase liquids (NAPLs) dissolve slowly in groundwater and/or volatilize in the vadose zone threatening the environment and public health over extended periods of time. The failure of a treatment technology at any given site is often due to the unnoticed presence of dissolved NAPL trapped in low permeability areas and/or the remaining presence of substantial amounts of pure phase NAPL after remediation efforts. The design of remediation strategies and the determination of remediation endpoints are traditionally carried out without quantifying risks associated with the failure of such efforts. We conduct a probabilistic risk analysis (PRA) to estimate the likelihood of failure of an on-site NAPL treatment technology. The PRA integrates all aspects of the problem (causes, pathways, and receptors) without resorting to extensive modeling. It accounts for a combination of multiple mechanisms of failure of a monitoring system, such as bypassing, insufficient sampling frequency and malfunctioning of the observation wells. We use a Bayesian framework to update the likelihood of failure of the treatment technology with observed measurements of concentrations at nearby monitoring wells.

Hybridized Design Concepts and Their Application to ERD Systems

Optimization of large‐scale injection‐based remedial systems requires engineering to intentionally capitalize on the biological, chemical, and physical mechanisms that occur within and between the zones of reagent application. These types of systems can be called hybrid designs as they employ multiple processes to achieve remediation endpoints (Figure 1), resulting in optimized system performance and a reduction in the overall life‐cycle cost. While all remedial applications incorporate these mechanisms to some extent, the importance of each of these processes is magnified in large‐scale applications. This column discusses the dominant mechanisms responsible for mass reduction within both source and distal plume footprints, with a focus on the application of “Hybridized Design” for enhanced reductive dechlorination (ERD) systems. Diagram showing the hybrid design approach which encompasses physical (sorption, advection, diffusion), chemical (mass flux, abiotic degradation) and biological (metabolic and cometabolic degradation) processes. image

Verification of Site Cleanup Using a Dual-Equilibrium Desorption Model

Abstract The dual-equilibrium desorption (DED) model was introduced in 2002 illustrating the apparent biphasic nature of desorption of contaminants from soil. This model was applied to a gasoline fueling station investigation in northwestern California in which it was tested against empirical data for dissolved-phase gasoline constituents in groundwater in the attempt to predict remediation end points. In this work, we are further testing the applicability of DED to field-based investigations by comparing its model predictions with flow modeling and advective-dispersive contaminant transport empirical results. The model was tested in two parts. First, groundwater flow was modeled using the Laplace equation for groundwater flow in two dimensions. This was compared to observations and used to describe the distribution of hydraulic heads across the flow field. This distribution was then used as input for the contaminant transport model. The equations for advective-dispersive transport were modified by substituting the distribution coefficient with that predicted in the DED. Predicted concentrations were then compared to concentrations observed in a downgradient monitoring well. Results indicate that this procedure is useful in predicting the timing of spikes in contaminant transport to within six months to one year and that the magnitude of those peaks were under predicted by approximately 2 %. The implication is that contaminant spikes are not transported in coherent packets but in a more chaotic manner, possible due to fractionation of gasoline. While better precision may be desired, this is one of the few methods available for use with chemical mixtures, such as total petroleum hydrocarbons.

Biodegradation of Organic Chemicals in Soil/Water Microcosms System: Model Development

The chemical interactions of hydrophobic organic contaminants with soils and sediments may result in strong binding and slow subsequent release rates that significantly affect remediation rates and endpoints. In order to illustrate the recalcitrance of chemical to degradation on sites, a sorption mechanism of intraparticle sequestration was postulated to operate on chemical remediation sites. Pseudo-first order sequestration kinetics is used in the study with the hypothesis that sequestration is an irreversibly surface-mediated process. A mathematical model based on mass balance equations was developed to describe the fate of chemical degradation in soil/water microcosm systems. In the model, diffusion was represented by Fick’s second law, local sorption-desorption by a linear isotherm, irreversible sequestration by a pseudo-first order kinetics and biodegradation by Monod kinetics. Solutions were obtained to provide estimates of chemical concentrations. The mathematical model was applied to a benzene biodegradation batch test and simulated model responses correlated well compared to measurements of biodegradation of benzene in the batch soil/water microcosm system. A sensitivity analysis was performed to assess the effects of several parameters on model behavior. Overall chemical removal rate decreased and sequestration increased quickly with an increase in the sorption partition coefficient. When soil particle radius, a, was greater than 1 mm, an increase in radius produced a significant decrease in overall chemical removal rate as well as an increase in sequestration. However, when soil particle radius was less than 0.1 mm, an increase in radius resulted in small changes in the removal rate and sequestration. As pseudo-first order sequestration rate increased, both chemical removal rate and sequestration increased slightly. Model simulation results showed that desorption resistance played an important role in the bioavailability of organic chemicals in porous media. Complete biostabilization of chemicals on remediation sites can be achieved when the concentration of the reversibly sorbed chemical reduces to zero (i.e., undetectable), with a certain amount of irreversibly sequestrated chemical left inside the soil particle solid phase.

Effects of source zone heterogeneity on surfactant-enhanced NAPL dissolution and resulting remediation end-points

The effectiveness of removal of nonaqueous phase liquids (NAPLs) from the entrapment source zone of the subsurface has been limited by soil heterogeneity and the inability to locate all entrapped sources. The goal of this study was to demonstrate the uncertainty of degree of source removal associated with aquifer heterogeneity. In this demonstration, source zone NAPL removal using surfactant-enhanced dissolution was considered. Model components that simulate the processes of natural dissolution in aqueous phase and surfactant-enhanced dissolution were incorporated into an existing code of contaminant transport. The dissolution modules of the simulator used previously developed Gilland–Sherwood type phenomenological models of NAPL dissolution to estimate mass transfer coefficients that are upscaleable to multidimensional flow conditions found at field sites. The model was used to simulate the mass removal from 10 NAPL entrapment zone configurations based on previously conducted two-dimensional tank experiments. These entrapment zones represent the NAPL distribution in spatially correlated random fields of aquifer hydraulic conductivity. The numerical simulations representing two-dimensional conditions show that effectiveness of mass removal depends on the aquifer heterogeneity that controls the NAPL entrapment and delivery of the surfactant to the locations of entrapped NAPLs. Flow bypassing resulting from heterogeneity and the reduction of relative permeability due to NAPL entrapment reduces the delivery efficiency of the surfactant, thus prolonging the remediation time to achieve desired end-point NAPL saturations and downstream dissolved concentrations. In some extreme cases, the injected surfactant completely bypassed the NAPL source zones. It was also found that mass depletion rates for different NAPL source configurations vary significantly. The study shows that heterogeneity result in uncertainties in the mass removal and achievable end-points that are directly related to dissolved contaminant plume development downstream of the NAPL entrapment zone.

Long-Term Bioremediation of a Subsurface Plume in Silty Soil

In northern Virginia, a loss from a tank farm has produced two plumes, containing about 200,000 gal of diesel fuel, jet-A fuel, and gasoline. Evidence suggests that the longest part of the contamination plume moved to its present length of 2,500 ft in less than 5 years. Since natural biodegradation would require about 2,500 years to reduce the hydrocarbon contamination to the remediation endpoints, other methods have been considered. Excavation of the plumes would take an estimated 5 years. However, the tank farm is surrounded by commercial buildings and expensive homes, and many of these buildings would have to be removed to reach the plumes. Enhanced natural bioremediation would require about 200 years at a start-up cost of about $1 million dollars and recurring costs of approximately $500,000/year. Infiltration galleries and enhanced subsurface permeability could reduce the remediation time to as little as 20 years.