Vadose zone attenuation of organic compounds at a crude oil spill site - interactions between biogeochemical reactions and multicomponent gas transport

Sorption and attenuation of petroleum VOCs in five unsaturated soils: Microcosms and column experiments

The fate of volatile organic compounds (VOC) vapors in the unsaturated zone is the basis for evaluating the natural attenuation potential and vapor intrusion risk. Microcosm and column experiments were conducted to study the effects chemical speciation and soil types/properties on the fate of petroleum VOCs in unsaturated zone. The biodegradation and total attenuation rates of the seven VOCs obtained by microcosm experiments in black soil and yellow earth were also generally higher than those in floodplain soil, lateritic red earth, and quartz sand. The VOC vapors in floodplain soil, lateritic red earth, and quartz sand showed slow total attenuation rates (<0.3 d-1). N-pentane, methylcyclopentane, and methylcyclohexane showed lower biodegradation rates than octane and three monoaromatic hydrocarbons. Volatilization into the atmosphere and biodegradation are two important natural attenuation paths for VOCs in unsaturated soil columns. The volatilization loss fractions of different volatile hydrocarbons in all five unsaturated soils were generally in the order: n-pentane (93.5%-97.8%) > methylcyclopentane (77.2%-85.5%) > methylcyclohexane (53.5%-69.2%) > benzene (17.1%-73.3%) > toluene (0-45.7%) > octane (1.9%-34.2%) > m-xylene (0-5.7%). The fractions by volatilization into the atmosphere of all seven hydrocarbons in quartz sand, lateritic red earth, and floodplain soil were close and higher compared to the yellow earth and black soil. Overall, this study illustrated the important roles chemical speciation and soil properties in determining the vapor-phase transport and natural attenuation of VOCs in the unsaturated zone.

Advective and diffusive gas phase transport in vadose zones: Importance for defining vapour risks and natural source zone depletion of petroleum hydrocarbons

Quantifying the interlinked behaviour of the soil microbiome, fluid flow, multi-component transport and partitioning, and biodegradation is key to characterising vapour risks and natural source zone depletion (NSZD) of light non-aqueous phase liquid (LNAPL) petroleum hydrocarbons. Critical to vapour transport and NSZD is transport of gases through the vadose zone (oxygen from the atmosphere, volatile organic compounds (VOCs), methane and carbon dioxide from the zone of LNAPL biodegradation). Volatilisation of VOCs from LNAPL, aerobic biodegradation, methanogenesis and heat production all generate gas pressure changes that may lead to enhanced gas fluxes apart from diffusion. Despite the importance of the gaseous phase dynamics in the vadose zone processes, the relative pressure changes and consequent scales of advective (buoyancy and pressure driven) / diffusive transport is less studied. We use a validated multi-phase multi-component non-isothermal modelling framework to differentiate gas transport mechanisms. We simulate a multicomponent unweathered gasoline LNAPL with high VOC content to maximise the potential for pressure changes due to volatilisation and to enable the joint effects of methanogenesis and shallower aerobic biodegradation of vapours to be assessed, along with heat production. Considering a uniform fine sand profile with LNAPL resident in the water table capillary zone, results suggest that biodegradation plays the key role in gas phase formation and consequent pressure build-up. Results suggest that advection is the main transport mechanism over a thin zone inside the LNAPL/capillary region, where the effective gaseous diffusion is very low. In the bulk of the vadose zone above the LNAPL region, the pressure change is minimal, and gaseous diffusion is dominant. Even for high biodegradation rate cases, pressure build-up due to heat generation (inducing buoyancy effects) is smaller than the contribution of gas formation due to biodegradation. The findings are critical to support broader assumptions of diffusive transport being dominant in vapour transport and NSZD assessments.

Influence of advection on the soil gas radon deficit technique for the quantification of LNAPL

The Radon (Rn) deficit technique is a rapid, low-cost, and non-invasive method to identify and quantify light non-aqueous phase liquids (LNAPL) in the soil. LNAPL saturation is typically estimated from Rn deficit using Rn partition coefficients, assuming equilibrium conditions. This work examines the applicability of this method in the presence of local advective fluxes that can be generated by groundwater fluctuations or biodegradation processes in the source zone. To this end, a one-dimensional analytical model was developed to simulate the steady-state diffusive-advective transport of soil gas Rn in the presence of LNAPL. The analytical solution was first validated against an existing numerical model adapted to include advection. Then a series of simulations to study the effect of advection on Rn profiles were carried out. It was found that in high-permeability soils (such as sandy soils), advective phenomena can significantly affect Rn deficit curves in the subsurface compared with those expected, assuming either equilibrium conditions or a diffusion-dominated transport. Namely, in the presence of pressure gradients generated by groundwater fluctuations, applying the traditional Rn deficit technique (assuming equilibrium conditions) can lead to an underestimation of LNAPL saturation. Furthermore, in the presence of methanogenesis processes (e.g., in the case of a fresh LNAPL of petroleum hydrocarbons), local advective fluxes can be expected above the source zone. In such cases, Rn concentrations above the source zone can be higher than those above background areas without advective phenomena, resulting in Rn deficits higher than 1 (i.e., Rn excess), and thus leading to a wrong interpretation regarding the presence of LNAPL in the subsurface if advection is not considered. Overall, the results obtained suggest that advection should be considered in the presence of pressure gradients in the subsurface to ensure an effective application of the soil gas Rn-deficit technique for quantitative estimation of LNAPL saturation.

Effects of freeze-thaw cycles on methanogenic hydrocarbon degradation: Experiment and modeling

Cold regions are warming much faster than the global average, resulting in more frequent and intense freeze-thaw cycles (FTCs) in soils. In hydrocarbon-contaminated soils, FTCs modify the biogeochemical and physical processes controlling petroleum hydrocarbon (PHC) biodegradation and the associated generation of methane (CH₄) and carbon dioxide (CO₂). Thus, understanding the effects of FTCs on the biodegradation of PHCs is critical for environmental risk assessment and the design of remediation strategies for contaminated soils in cold regions. In this study, we developed a diffusion-reaction model that accounts for the effects of FTCs on toluene biodegradation, including methanogenic biodegradation. The model is verified against data generated in a 215 day-long batch experiment with soil collected from a PHC contaminated site in Ontario, Canada. The fully saturated soil incubations with six different treatments were exposed to successive 4-week FTCs, with temperatures oscillating between -10 °C and +15 °C, under anoxic conditions to stimulate methanogenic biodegradation. We measured the headspace concentrations and ¹³C isotope compositions of CH₄ and CO₂ and analyzed the porewater for pH, acetate, dissolved organic and inorganic carbon, and toluene. The numerical model represents solute diffusion, volatilization, sorption, as well as a reaction network of 13 biogeochemical processes. The model successfully simulates the soil porewater and headspace concentration time series data by representing the temperature dependencies of microbial reaction and gas diffusion rates during FTCs. According to the model results, the observed increases in the headspace concentrations of CH₄ and CO₂ by 87% and 136%, respectively, following toluene addition are explained by toluene fermentation and subsequent methanogenesis reactions. The experiment and the numerical simulation show that methanogenic degradation is the primary toluene attenuation mechanism under the electron acceptor-limited conditions experienced by the soil samples, representing 74% of the attenuation, with sorption contributing to 11%, and evaporation contributing to 15%. Also, the model-predicted contribution of acetate-based methanogenesis to total produced CH₄ agrees with that derived from the ¹³C isotope data. The freezing-induced soil matrix organic carbon release is considered as an important process causing DOC increase following each freezing period according to the calculations of carbon balance and SUVA index. The simulation results of a no FTC scenario indicate that, in the absence of FTCs, CO₂ and CH₄ generation would decrease by 29% and 26%, respectively, and that toluene would be biodegraded 23% faster than in the FTC scenario. Because our modeling approach represents the dominant processes controlling PHC biodegradation and the associated CH₄ and CO₂ fluxes, it can be used to analyze the sensitivity of these processes to FTC frequency and duration driven by temperature fluctuations.

Natural Source Zone Depletion (NSZD) Quantification Techniques: Innovations and Future Directions

Natural source zone depletion (NSZD) is an emerging technique for sustainable and cost-effective bioremediation of light non-aqueous phase liquid (LNAPL) in oil spill sites. Depending on regulatory objectives, NSZD has the potential to be used as either the primary or sole LNAPL management technique. To achieve this goal, NSZD rate (i.e., rate of bulk LNAPL mass depletion) should be quantified accurately and precisely. NSZD has certain characteristic features that have been used as surrogates to quantify the NSZD rates. This review highlights the most recent trends in technology development for NSZD data collection and rate estimation, with a focus on the operational and technical advantages and limitations of the associated techniques. So far, four principal techniques are developed, including concentration gradient (CG), dynamic closed chamber (DCC), CO2 trap and thermal monitoring. Discussions revolving around two techniques, “CO2 trap” and “thermal monitoring”, are expanded due to the particular attention to them in the current industry. The gaps of knowledge relevant to the NSZD monitoring techniques are identified and the issues which merit further research are outlined. It is hoped that this review can provide researchers and practitioners with sufficient information to opt the best practice for the research and application of NSZD for the management of LNAPL impacted sites.

Tracking NSZD mass removal rates over decades: Site-wide and local scale assessment of mass removal at a legacy petroleum site

Long-term estimates of natural source zone depletion (NSZD) rates for petroleum LNAPL (light non-aqueous phase liquid) sites are not available. One-off measurements are often thought valid over the lifetime of LNAPL sites. In the context of site-wide LNAPL mass estimates, we report site-specific gasoline and diesel NSZD rates spanning 21-26 years. Using depth profiles of soil gases (oxygen, carbon dioxide, methane, volatiles) above LNAPL, NSZD rates were estimated in 1994, 2006 and 2020 for diesel and 1999, 2009 and 2020 for gasoline. Each date also had soil-core mass estimates, which together with NSZD rates allow estimation of the longevity for LNAPL presence. Site-wide coring (in 1992, 2002, 2007) estimated LNAPL mass reductions of 12,000 t. For diesel NSZD, the ratio of NSZD rates for 2006 (16,000-49,000 L/ha/y) to those in 2020 (2600-14,000 L/ha/y) was ~3-6. By 2020, the 1994 diesel NSZD rates would have predicted the entire removal of measured mass (16-42 kg/m²). For gasoline, NSZD rates in 1999 were extremely high (50,000-270,000 L/ha/y) but 9-27 times lower (5800-10,000 L/ha/y) a decade later. The gasoline NSZD rates in 1999 predicted near complete mass removal in 2-12 years, but 10-11 kg/m² was measured 10 and 21 years later which is 26% of the initial mass in 1999. The outcomes substantiate the need to understand NSZD rate changes over the lifetime of LNAPL-impacted sites.

Towards a digital twin for characterising natural source zone depletion: A feasibility study based on the Bemidji site

Natural source zone depletion (NSZD) of light non-aqueous phase liquids (LNAPLs) may be a valid long-term management option at petroleum impacted sites. However, its future long-term reliability needs to be established. NSZD includes partitioning, biotic and abiotic degradation of LNAPL components plus multiphase fluid dynamics in the subsurface. Over time, LNAPL components are depleted and those partitioning to various phases change, as do those available for biodegradation. To accommodate these processes and predict trends and NSZD over decades to centuries, for the first time, we incorporated a multi-phase multi-component multi-microbe non-isothermal approach to representatively simulate NSZD at field scale. To validate the approach we successfully mimic data from the LNAPL release at the Bemidji site. We simulate the entire depth of saturated and unsaturated zones over the 27 years of post-release measurements. The study progresses the idea of creating a generic digital twin of NSZD processes and future trends. Outcomes show the feasibility and affordability of such detailed computational approaches to improve decision-making for site management and restoration strategies. The study provided a basis to progress a computational digital twin for complex subsurface systems.

Refinement of the gradient method for the estimation of natural source zone depletion at petroleum contaminated sites

Rates of natural source zone depletion (NSZD) are increasingly being used to aid remedial decision making and light non-aqueous phase liquid (LNAPL) longevity estimates at petroleum release sites. Current NSZD estimate methods, based on analyses of carbon dioxide (CO₂) and oxygen (O₂) soil-gas concentration gradients ("gradient method") assume linear concentration profiles with depth. This assumption can underestimate the concentration gradients especially above LNAPL sources that are typically characterized by curvilinear or semi-curvilinear O₂ and CO₂ concentration profiles. In this work, we proposed a new method that relies on calculating the O₂ and CO₂ concentration gradient using a first-order reaction model. The method requires an estimate of the diffusive reaction length that can be easily derived from soil-gas concentration data. A simple step-by-step guide for applying the new method is provided. Nomographs were also developed to facilitate method application. Application of the nomographs using field data from published literature showed that NSZD rates could be underestimated by nearly an order of magnitude if reactivity in the vadose zone is not accounted for. The new method helps refine NSZD rates estimation and improve risk-based decision making at certain petroleum contaminated sites.

Quantifying natural source zone depletion at petroleum hydrocarbon contaminated sites: A comparison of 14C methods

Surficial CO₂ efflux surveys have been used to delineate hydrocarbon source zones in contaminated aquifers and provide estimates of hydrocarbon biodegradation rates. This approach requires distinguishing between CO₂ derived from petroleum degradation and CO₂ produced from natural soil respiration. To this end, radiocarbon has been used to differentiate between ¹⁴C-depleted CO₂ from hydrocarbon degradation and ¹⁴C-enriched CO₂ from natural soil respiration to effectively quantify the contribution of each source to total CO₂ efflux, and by deduction natural source zone depletion (NSZD) rates. In this study, a systematic method comparison has been conducted to evaluate available approaches for collecting CO₂ gas samples for radiocarbon analysis used to correct total CO₂ efflux measurements for quantifying natural source zone depletion rates. Gas samples for radiocarbon analysis were sampled from (i) dynamic closed chambers (located at ground surface), (ii) static chambers (also at ground surface), (iii) shallow soil gas probes (0.3 m bgs), and (iv) soil gas monitoring wells (~0.6 m below ground surface) during a CO₂ efflux survey conducted at the site of a historical pipeline rupture near Bemidji, MN. The mean fraction of radiocarbon (F¹⁴C) obtained from samples overlying the source zone were (i) 0.93 ± 0.01, (ii) 0.73 ± 0.03, (iii) 0.71 ± 0.04, and (iv) 0.41 ± 0.06, for the four methods respectively. These F¹⁴C values were used to apportion total CO₂ efflux measurements into contributions of contaminant-derived CO₂ efflux and natural soil respiration to evaluate natural source zone depletion processes. Results suggest that the method of radiocarbon sampling has a significant effect on the calculated fraction of the CO₂ efflux originating from contaminant-related soil respiration, with contributions ranging between 27% and 59% of total soil respiration. Results indicate that radiocarbon sampled from static chambers and shallow soil gas probes methods offer the best compromise between CO₂ sample yield and sample representativeness, providing the most reliable estimates of CO₂ effluxes originating from contaminant degradation. However, the results also show that for this study, all methods agree within a factor of <2.3 regarding the inferred NSZD rates.

Impacts of water table fluctuations on actual and perceived natural source zone depletion rates

A viable means of quantifying the rate of natural source zone depletion (NSZD) at hydrocarbon contaminated sites is by the measurement of carbon dioxide (CO₂) and methane (CH₄) effluxes at the surface. This methodology assumes that gas effluxes are reflective of actual contaminant degradation rates in the subsurface, which is only accurate for quasi-steady state conditions. However, in reality, subsurface systems are highly dynamic, often resulting in fluctuations of the water table. To quantify the effects of water table fluctuations on NSZD rates, a simulated biodiesel spill in a 400 cm long, 100 cm wide and 150 cm tall sandtank was subjected to lowering and raising the water table, while soil-gas chemistry and surface CO₂ and CH₄ effluxes were measured. Results show that water table fluctuations have both short-term (perceived) and long-term (actual) effects on NSZD rates, interpreted using surface efflux measurements. When the water table was lowered, surface effluxes immediately increased up to 3 and 344 times higher than baseline for CO₂ and CH₄ effluxes, respectively, due to the liberation of anaerobically produced gas accumulated below the water table. After this immediate release, the system then reached quasi-steady state conditions 1.4 to 1.6 times higher for CO₂ than baseline conditions, attributed to increased aerobic degradation in the broadened and exposed smear zone. When the water table was raised, quasi-steady state CO₂ and CH₄ effluxes declined to values of 0.9 and 0.4 times baseline effluxes, respectively, implying that contaminant degradation rates were reduced due to submergence of the smear zone. The findings of this study show that the dynamic effects of water table fluctuations and redistribution of the contaminants affect surface effluxes as well as short-term (perceived) and long-term (actual) contaminant degradation rates. Therefore, water table fluctuations need to be considered when quantifying NSZD at contaminated sites using sparse temporal rate measurements to estimate NSZD rates for extended periods of time (e.g., annual rates).

Benzene homologues contaminants in a former herbicide factory site: distribution, attenuation, risk, and remediation implication

Benzene homologues often used as organic raw materials or as detergents in chemical industry are prone to accidental release into the environment which can cause serious long-term soil pollutions. In a large former herbicide factory site, we investigated 43 locations for benzene homologues contaminations in soil, soil gas, and groundwater and studied the hydrogeological conditions. An inverse distance weighted interpolation method was employed to determine the pollutants three-dimensional spatial distribution in the soils. Results showed that benzene homologues residues were mainly originated from the herbicide production workshop and that the pollution had horizontally expanded at the deeper soil layer. Contaminants had already migrated 15 m downward from ground surface. Contaminant phase distribution study showed that NAPL was the primary phase (> 99%) for the pollutants accumulated in the unsaturated zone, while it had not migrated to groundwater. The primary mechanism for contaminant transport and attenuation included dissolution of "occluded" NAPL into pore water and pollutant volatilization into soil pore space. Risk assessment revealed that the pollutants brought unacceptable high carcinogenic and non-carcinogenic risks to public health. In order to convert this former chemical processing factory site into a residential area, a remediation to the polluted production workshop sites is urgently required.

Laboratory-scale experimental and modelling investigations of 222Rn profiles in chemically heterogeneous LNAPL contaminated vadose zones

The potential of LNAPL delineation by ²²²Rn soil-gas monitoring in a chemically heterogeneous vadose zone was investigated in this study based on laboratory (batch and columns) experiments and numerical modelling. An enhanced version of the MIN3P reactive transport code was used to simulate Rn transport in both uncontaminated and NAPL-contaminated vadose zones and results were validated against analytical solutions and laboratory experiments. Results show that ²²²Rn activity profiles are mainly controlled by porous media ²²²Rn production, vadose zone fluid saturations and especially the type and distribution of NAPL in contaminated areas. The results also show that decreases in ²²²Rn activity and variations in activity gradients provide evidence for the presence and saturation of NAPL. This study demonstrates that LNAPL delineation via ²²²Rn gas surveys at contaminated sites works best, if gas measurements extend as deep as possible and include regions where ²²²Rn activity decreases due to elevated NAPL content. In addition, collection and analysis of depth-discrete gas samples allows the characterization of vertical NAPL distribution based on the ²²²Rn activity gradient. The determination of ²²²Rn production in the unsaturated zone, as well as water capillary pressure curves are of key importance in enabling the discrimination of an uncontaminated from a NAPL-contaminated area.

Natural source zone depletion of LNAPL: A critical review supporting modelling approaches

Natural source zone depletion (NSZD) of light non-aqueous phase liquids (LNAPLs) includes partitioning, transport and degradation of LNAPL components. NSZD is being considered as a site closure option during later stages of active remediation of LNAPL contaminated sites, and where LNAPL mass removal is limiting. To ensure NSZD meets compliance criteria and to design enhanced NSZD actions if required, residual risks posed by LNAPL and its long term behaviour require estimation. Prediction of long-term NSZD trends requires linking physicochemical partitioning and transport processes with bioprocesses at multiple scales within a modelling framework. Here we expand and build on the knowledge base of a recent review of NSZD, to establish the key processes and understanding required to model NSZD long term. We describe key challenges to our understanding, inclusive of the dominance of methanogenic or aerobic biodegradation processes, the potentially changeability of rates due to the weathering profile of LNAPL product types and ages, and linkages to underlying bioprocesses. We critically discuss different scales in subsurface simulation and modelling of NSZD. Focusing on processes at Darcy scale, 36 models addressing processes of importance to NSZD are investigated. We investigate the capabilities of models to accommodate more than 20 subsurface transport and transformation phenomena and present comparisons in several tables. We discuss the applicability of each group of models for specific site conditions.

Multi-year CO2 efflux measurements for assessing natural source zone depletion at a large hydrocarbon-impacted site

The changing landscape of fuel consumption related, in part, to increased engine efficiency and the inexpensive supply of natural gas, has led to the closure of multiple refineries. As the operational lifetime of many refineries exceeds 100 years, historical releases of oil and refined products is common. To evaluate remediation and rehabilitation options, there is a need to understand the rate and distribution of natural hydrocarbon mass losses across these large properties. Here, surficial CO₂ flux measurements were used to evaluate naturally occurring hydrocarbon mass losses at a large-scale former refinery that has been closed since 1982. Natural source zone depletion (NSZD) rates over a five-year period (2012-2016) were derived from surficial CO₂ efflux measurements on a high-resolution grid (N > 80). Results demonstrate substantial variations of mass loss rates across the site. Average site-wide mass loss rates ranged from 1.1-5.4 g TPH m^-2 d^-1 as C₁₀H₂₂ with a multi-year average of 4.0 g TPH m^-2 d^-1 as decane (C₁₀H₂₂), consistent with observations at other sites. Statistical analysis demonstrated that the same average mass loss rates would have been obtained with fewer measurement locations (N = 20-30). Comparing NSZD rates to site metadata show CO₂ fluxes to be a reasonably good proxy for zones of subsurface hydrocarbon contamination - particularly with respect to vadose zone impacts. It is hypothesized that the observed decline of NSZD rates over the study period is related to rise of groundwater levels, leading to increased submergence of the smear zone. Overall, mass loss rates calculated from CO₂ fluxes show NSZD can result in substantial contaminant removal, which may rival that obtained from engineered remediation, under some conditions.

Weathering of Oil in a Surficial Aquifer

The composition of crude oil in a surficial aquifer was determined in two locations at the Bemidji, MN, spill site. The abundances of 71 individual hydrocarbons varied within 16 locations sampled. Little depletion of these hydrocarbons (relative to the pipeline oil) occurred in the first 10 years after the spill, whereas losses of 25% to 85% of the total measured hydrocarbons occurred after 30 years. The C_6-30 n-alkanes, toluene, and o-xylene were the most depleted hydrocarbons. Some hydrocarbons, such as the n-C_10-24 cyclohexanes, tri- and tetra- methylbenzenes, acyclic isoprenoids, and naphthalenes were the least depleted. Benzene was detected at every sampling location 30 years after the spill. Degradation of the oil led to increases in the percent organic carbon and in the δ ¹³ C of the oil. Another method of determining hydrocarbon loss was by normalizing the total measured hydrocarbon concentrations to that of the most conservative analytes. This method indicated that the total measured hydrocarbons were depleted by 47% to 77% and loss of the oil mass over 30 years was 18% to 31%. Differences in hydrocarbon depletion were related to the depth of the oil in the aquifer, local topography, amount of recharge reaching the oil, availability of electron acceptors, and the presence of less permeable soils above the oil. The results from this study indicate that once crude oil has been in the subsurface for a number of years there is no longer a "starting oil concentration" that can be used to understand processes that affect its fate and the transport of hydrocarbons in groundwater.

Applying the Rayleigh Approach for Stable Isotope-Based Analysis of VOC Biodegradation in Diffusion-Dominated Systems

Compound-specific stable isotope analysis (CSIA) has become an established tool for assessing biodegradation in the subsurface. Diffusion-dominated vapor phase transport thereby is often excluded from quantitative assessments due to the problem of diffusive mixing of concentrations with different isotopic signatures for CSIA interpretation. In soils and other unsaturated porous media volatile organic compounds (VOCs) however, are mainly transported via gas-phase diffusion and may thus prohibit a CSIA-based quantitative assessment of the fate of VOCs. The present study presents and verifies a concept for the assessment of biodegradation-induced stable isotope fractionation along a diffusive transport path of VOCs in unsaturated porous media. For this purpose data from batch and column toluene biodegradation experiments in unsaturated porous media were combined with numerical reactive transport simulations; both addressing changes of concentration and stable isotope fractionation of toluene. The numerical simulations are in good agreement with the experiment data, and our results show that the presented analytically derived assessment concept allows using the slope of the Rayleigh plot to obtain reasonable estimates of effective in situ fractionation factors in spite of diffusion-dominated transport. This enlarges the application range of CSIA and provides a mean for a better understanding of VOC fate in the unsaturated subsurface.

Relative contributions of microbial and infrastructure heat at a crude oil-contaminated site

Biodegradation of contaminants can increase the temperature in the subsurface due to heat generated from exothermic reactions, making temperature observations a potentially low-cost approach for determining microbial activity. For this technique to gain more widespread acceptance, it is necessary to better understand all the factors affecting the measured temperatures. Biodegradation has been occurring at a crude oil-contaminated site near Bemidji, Minnesota for 39 years, creating a quasi-steady-state plume of contaminants and degradation products. A model of subsurface heat generation and transport helps elucidate the contribution of microbial and infrastructure heating to observed temperature increases at this site. We created a steady-state, two-dimensional, heat transport model using previous-published parameter values for physical, chemical and biodegradation properties. Simulated temperature distributions closely match the observed average annual temperatures measured in the contaminated area at the site within less than 0.2 °C in the unsaturated zone and 0.4 °C in the saturated zone. The model results confirm that the observed subsurface heat from microbial activity is due primarily to methane oxidation in the unsaturated zone resulting in a 3.6 °C increase in average annual temperature. Another important source of subsurface heat is from the active, crude-oil pipelines crossing the site. The pipelines impact temperatures for a distance of 200 m and contribute half the heat. Model results show that not accounting for the heat from the pipelines leads to overestimating the degradation rates by a factor of 1.7, demonstrating the importance of identifying and quantifying all heat sources. The model results also highlighted a zone where previously unknown microbial activity is occurring at the site.

Understanding complex LNAPL sites: Illustrated handbook of LNAPL transport and fate in the subsurface

The goal of the paper is to highlight the management of the complexities and risks for light non-aqueous phase liquid (LNAPL) sites, and how the Illustrated Handbook of LNAPL Transport and Fate in the Subsurface (CL:AIRE, London. ISBN 978-1-905046-24-9. http://www.claire.co.uk/LNAPL; "LNAPL illustrated handbook") is useful guidance and a tool for professionals to understand these complexities and risks. The LNAPL illustrated handbook provides a clear and concise best-practice guidance document, which is a valuable decision support tool for use in discussions and negotiations regarding LNAPL impacted sites with respect to the risks of LNAPL. The LNAPL illustrated handbook is a user-friendly overview of the nature of LNAPL contamination in various geological settings including unconsolidated, consolidated, and fractured rock environments to best understand its fate and behavior leading to the appropriate management and/or remedial approach of the two major risks associated with a LNAPL source. As a source term, LNAPL has chemicals that form dissolved- and vapor-phase plumes, which are referred to as composition-based risks; and being a liquid there is the risk that the source may expand impacting a greater volume of the aquifer, which are referred to as saturation-based risks. There have been significant developments in recent years on the understanding of the complex behavior of LNAPL and associated groundwater and vapor plumes; however, the state of practice has often lagged these improvements in knowledge. The LNAPL illustrated handbook aids the site investigator, site owners, and regulators to understand these risks, and understand how these risks behave through better conceptual understanding of LNAPL transport and fate in the subsurface.

Bubble-Facilitated VOC Transport from LNAPL Smear Zones and Its Potential Effect on Vapor Intrusion

Most conceptual and mathematical models of soil vapor intrusion assume that the transport of volatile organic compounds (VOCs) from a source toward a building is limited by diffusion through the soil gas. Under conditions where advection occurs, transport rates are higher and can lead to higher indoor air concentrations. Advection-dominated conditions can be created by gas bubble flow in the saturated zone. A series of laboratory column experiments were conducted to measure mass flux due to bubble-facilitated VOC transport from light nonaqueous phase liquid (LNAPL) smear zones. Smear zones that contained both LNAPL residual and trapped gas, as well as those that contained only LNAPL residual, were investigated. Results showed that the VOC mass flux due to bubble-facilitated transport was orders-of-magnitude higher than under diffusion-limited conditions. Results also showed that the mass flux due to bubble-facilitated transport was intermittent, and increased with an increased supply of dissolved gases.

Measuring Carbon-based Contaminant Mineralization Using Combined CO2 Flux and Radiocarbon Analyses

A method is described which uses the absence of radiocarbon in industrial chemicals and fuels made from petroleum feedstocks which frequently contaminate the environment. This radiocarbon signal - or rather the absence of signal - is evenly distributed throughout a contaminant source pool (unlike an added tracer) and is not impacted by biological, chemical or physical processes (e.g., the ¹⁴C radioactive decay rate is immutable). If the fossil-derived contaminant is fully degraded to CO2, a harmless end-product, that CO2 will contain no radiocarbon. CO2 derived from natural organic matter (NOM) degradation will reflect the NOM radiocarbon content (usually <30,000 years old). Given a known radiocarbon content for NOM (a site background), a two end-member mixing model can be used to determine the CO2 derived from a fossil source in a given soil gas or groundwater sample. Coupling the percent CO2 derived from the contaminant with the CO2 respiration rate provides an estimate for the total amount of contaminant degraded per unit time. Finally, determining a zone of influence (ZOI) representing the volume from which site CO2 is collected allows determining the contaminant degradation per unit time and volume. Along with estimates for total contaminant mass, this can ultimately be used to calculate time-to-remediate or otherwise used by site managers for decision-making.

Biodegradation of vapor-phase toluene in unsaturated porous media: Column experiments

Biodegradation of organic chemicals in the vapor phase of soils and vertical flow filters has gained attention as promising approach to clean up volatile organic compounds (VOC). The drivers of VOC biodegradation in unsaturated systems however still remain poorly understood. Here, we analyzed the processes controlling aerobic VOC biodegradation in a laboratory setup mimicking the unsaturated zone above a shallow aquifer. The setup allowed for diffusive vapor-phase transport and biodegradation of three VOC: non-deuterated and deuterated toluene as two compounds of highly differing biodegradability but (nearly) identical physical and chemical properties, and MTBE as (at the applied experimental conditions) non-biodegradable tracer and internal control. Our results showed for toluene an effective microbial degradation within centimeter VOC transport distances despite high gas-phase diffusivity. Degradation rates were controlled by the reactivity of the compounds while oxic conditions were found everywhere in the system. This confirms hypotheses that vadose zone biodegradation rates can be extremely high and are able to prevent the outgassing of VOC to the atmosphere within a centimeter range if compound properties and site conditions allow for sufficiently high degradation rates.

PHT3D-UZF: A Reactive Transport Model for Variably-Saturated Porous Media

A modified version of the MODFLOW/MT3DMS-based reactive transport model PHT3D was developed to extend current reactive transport capabilities to the variably-saturated component of the subsurface system and incorporate diffusive reactive transport of gaseous species. Referred to as PHT3D-UZF, this code incorporates flux terms calculated by MODFLOW's unsaturated-zone flow (UZF1) package. A volume-averaged approach similar to the method used in UZF-MT3DMS was adopted. The PHREEQC-based computation of chemical processes within PHT3D-UZF in combination with the analytical solution method of UZF1 allows for comprehensive reactive transport investigations (i.e., biogeochemical transformations) that jointly involve saturated and unsaturated zone processes. Intended for regional-scale applications, UZF1 simulates downward-only flux within the unsaturated zone. The model was tested by comparing simulation results with those of existing numerical models. The comparison was performed for several benchmark problems that cover a range of important hydrological and reactive transport processes. A 2D simulation scenario was defined to illustrate the geochemical evolution following dewatering in a sandy acid sulfate soil environment. Other potential applications include the simulation of biogeochemical processes in variably-saturated systems that track the transport and fate of agricultural pollutants, nutrients, natural and xenobiotic organic compounds and micropollutants such as pharmaceuticals, as well as the evolution of isotope patterns.

Relating subsurface temperature changes to microbial activity at a crude oil-contaminated site

Crude oil at a spill site near Bemidji, Minnesota has been undergoing aerobic and anaerobic biodegradation for over 30 years, creating a 150-200 m plume of primary and secondary contaminants. Microbial degradation generates heat that should be measurable under the right conditions. To measure this heat, thermistors were installed in wells in the saturated zone and in water-filled monitoring tubes in the unsaturated zone. In the saturated zone, a thermal groundwater plume originates near the residual oil body with temperatures ranging from 2.9°C above background near the oil to 1.2°C down gradient. Temperatures in the unsaturated zone above the oil body were up to 2.7°C more than background temperatures. Previous work at this site has shown that methane produced from biodegradation of the oil migrates upward and is oxidized in a methanotrophic zone midway between the water table and the surface. Enthalpy calculations and observations demonstrate that the temperature increases primarily result from aerobic methane oxidation in the unsaturated zone above the oil. Methane oxidation rates at the site independently estimated from surface CO2 efflux data are comparable to rates estimated from the observed temperature increases. The results indicate that temperature may be useful as a low-cost measure of activity but care is required to account for the correct heat-generating reactions, other heat sources and the effects of focused recharge.

A Petroleum Vapor Intrusion Model Involving Upward Advective Soil Gas Flow Due to Methane Generation

At petroleum vapor intrusion (PVI) sites at which there is significant methane generation, upward advective soil gas transport may be observed. To evaluate the health and explosion risks that may exist under such scenarios, a one-dimensional analytical model describing these processes is introduced in this study. This new model accounts for both advective and diffusive transport in soil gas and couples this with a piecewise first-order aerobic biodegradation model, limited by oxygen availability. The predicted results from the new model are shown to be in good agreement with the simulation results obtained from a three-dimensional numerical model. These results suggest that this analytical model is suitable for describing cases involving open ground surface beyond the foundation edge, serving as the primary oxygen source. This new analytical model indicates that the major contribution of upward advection to indoor air concentration could be limited to the increase of soil gas entry rate, since the oxygen in soil might already be depleted owing to the associated high methane source vapor concentration.

Measurement of Natural Losses of LNAPL Using CO₂ Traps

Efflux of CO2 above releases of petroleum light nonaqueous phase liquids (LNAPLs) has emerged as a critical parameter for resolving natural losses of LNAPLs and managing LNAPL sites. Current approaches for resolving CO2 efflux include gradient, flux chamber, and mass balance methods. Herein a new method for measuring CO2 efflux above LNAPL bodies, referred to as CO2 traps, is introduced. CO2 traps involve an upper and a lower solid phase sorbent elements that convert CO2 gas into solid phase carbonates. The sorbent is placed in an open vertical section of 10 cm ID polyvinyl chloride (PVC) pipe located at grade. The lower sorbent element captures CO2 released from the subsurface via diffusion and advection. The upper sorbent element prevents atmospheric CO2 from reaching the lower sorbent element. CO2 traps provide integral measurement of CO2 efflux based over the period of deployment, typically 2 to 4 weeks. Favorable attributes of CO2 traps include simplicity, generation of integral (time averaged) measurement, and a simple means of capturing CO2 for carbon isotope analysis. Results from open and closed laboratory experiments indicate that CO2 traps quantitatively capture CO2 . Results from the deployment of 23 CO2 traps at a former refinery indicate natural loss rates of LNAPL (measured in the fall, likely concurrent with high soil temperatures and consequently high degradation rates) ranging from 13,400 to 130,000 liters per hectare per year (L/Ha/year). A set of field triplicates indicates a coefficient of variation of 18% (resulting from local spatial variations and issues with measurement accuracy).

Documentation of time-scales for onset of natural attenuation in an aquifer treated by a crude-oil recovery system

A pipeline transporting crude-oil broke in a nature reserve in 2009 and spilled 5100 m(3) of oil that partly reached the aquifer and formed progressively a floating oil lens. Groundwater monitoring started immediately after the spill and crude-oil recovery by dual pump-and-skim technology was operated after oil lens formation. This study aimed at documenting the implementation of redox-specific natural attenuation processes in the saturated zone and at assessing whether dissolved compounds were degraded. Seven targeted water sampling campaigns were done during four years in addition to a routine monitoring of hydrocarbon concentrations. Liquid oil reached the aquifer within 2.5 months, and anaerobic processes, from denitrification to reduction of sulfate, were observable after 8 months. Methanogenesis appeared on site after 28 months. Stable carbon isotope analyses after 16 months showed maximum shifts in δ(13)C of +4.9±0.22‰ for toluene, +2.4±0.19‰ for benzene and +0.9±0.51‰ for ethylbenzene, suggesting anaerobic degradation of these compounds in the source zone. Estimations of fluxes of inorganic carbon produced by biodegradation revealed that, in average, 60% of inorganic carbon production was attributable to sulfate reduction. This percentage tended to decrease with time while the production of carbon attributable to methanogenesis was increasing. Within the investigation time frame, mass balance estimations showed that biodegradation is a more efficient process for control of dissolved concentrations compared to pumping and filtration on an activated charcoal filter.

Combined radiocarbon and CO2 flux measurements used to determine in situ chlorinated solvent mineralization rate

A series of combined measurements was made at the Naval Air Station North Island (NASNI) Installation Restoration Site 5, Unit 2 during July and August 2013. Combined measurements included CO2 respiration rate, CO2 radiocarbon content to estimate chlorinated hydrocarbon (CH) mineralization and a zone of influence (ZOI) model. CO2 was collected continuously over 2 two-week periods by recirculating monitoring well headspace gas through NaOH traps. A series of 12 wells in the main CH plume zone and a background well with no known historical contamination were sampled. The background well CO2 was used to determine radiocarbon content derived from respired natural organic matter. A two end-member mixing model was then used to determine the amount of CH-derived carbon present in the CO2 collected from plume region wells. The ZOI model provided an estimate for the soil volume sampled at each well. CH mineralization rates were highest upgradient and at the plume fringe for areas of high historical contamination and ranged from 0.02 to 5.6 mg CH carbon per day. Using the ZOI model volume estimates, CH-carbon removal ranged from 0.2 to 32 mg CH-carbon m(-3) per day. Because the rate estimates were based on a limited sampling (temporally), they were not further extrapolated to long-term contaminant degradation estimates. However, if the site manager or regulators required them, estimates - subject to long-term variability uncertainties - could be made using volume and rate data determined over short timescales. A more comprehensive seasonal sampling is needed to constrain long-term remediation models for the entire impacted area and identify environmental conditions related to more rapid turnover times amongst the wells.

A mass balance approach to investigating geochemical controls on secondary water quality impacts at a crude oil spill site near Bemidji, MN

Secondary water quality impacts can result from a broad range of coupled reactions triggered by primary groundwater contaminants. Data from a crude-oil spill research site near Bemidji, MN provide an ideal test case for investigating the complex interactions controlling secondary impacts, including depleted dissolved oxygen and elevated organic carbon, inorganic carbon, CH4, Mn, Fe, and other dissolved ions. To better understand these secondary impacts, this study began with an extensive data compilation of various data types, comprising aqueous, sediment, gas, and oil phases, covering a 260m cross-sectional domain over 30years. Mass balance calculations are used to quantify pathways that control secondary components, by using the data to constrain the sources and sinks for the important redox processes. The results show that oil constituents other than BTEX (benzene, toluene, ethylbenzene, o-, m- and p-xylenes), including n-alkanes and other aromatic compounds, play significant roles in plume evolution and secondary water quality impacts. The analysis underscores previous results on the importance of non-aqueous phases. Over 99.9% of the Fe(2+) plume is attenuated by immobilization on sediments as Fe(II) and 85-95% of the carbon biodegradation products are outgassed. Gaps identified in carbon and Fe mass balances and in pH buffering mechanisms are used to formulate a new conceptual model. This new model includes direct out-gassing of CH4 and CO2 from organic carbon biodegradation, dissolution of directly produced CO2, and sorption with H(+) exchange to improve pH buffering. The identification of these mechanisms extends understanding of natural attenuation of potential secondary impacts at enhanced reductive dechlorination sites, particularly for reduced Fe plumes, produced CH4, and pH perturbations.

Comparison of surficial CO2 efflux to other measures of subsurface crude oil degradation

At a spill site near Bemidji, Minnesota, crude oil at the water table has been undergoing anaerobic biodegradation for over 30years. Previous work at this site has shown that methane produced from biodegradation of the oil migrates upward and is oxidized in a methanotrophic zone midway between the water table and the surface. To compare microbial activity measurement methods from multiple locations in the oil body, surficial carbon dioxide efflux, methanogen and methanotroph concentrations, and oil degradation state were collected. Carbon dioxide effluxes over the oil body averaged more than four times those at the background site. Methanotrophic bacteria concentrations measured using pmoA were over 10(5) times higher above the oil-contaminated sediments compared with the background site. Methanogenic archaea measured using mcrA ranged from 10(5) to over 10(7) in the oil and were below detection in the background. Methanogens correlated very well with methanotroph concentrations (r=0.99), n-alkylcyclohexane losses as a proxy for degradation state (r=-0.96), and somewhat less well with carbon dioxide efflux (r=0.92). Carbon dioxide efflux similarly correlated to methanotroph concentrations (r=0.90) and n-alkylcyclohexane losses (r=-0.91).

Hexadecane and pristane degradation potential at the level of the aquifer--evidence from sediment incubations compared to in situ microcosms

Monitored natural attenuation is widely accepted as a sustainable remediation method. However, methods providing proof of proceeding natural attenuation within the water-unsaturated (vadose) zone are still relying on proxies such as measurements of reactive and non-reactive gases, or sediment sampling and subsequent mineralisation assays, under artificial conditions in the laboratory. In particular, at field sites contaminated with hydrophobic compounds, e.g. crude oil spills, an in situ evaluation of natural attenuation is needed, because in situ methods are assumed to provide less bias than investigations applying either proxies for biodegradation or off-site microcosm experiments. In order to compare the current toolbox of methods with the recently developed in situ microcosms, incubations with direct push-sampled sediments from the vadose and the aquifer zones of a site contaminated with crude oil were carried out in conventional microcosms and in situ microcosms. The results demonstrate the applicability of the in situ microcosm approach also outside water-saturated aquifer conditions in the vadose zone. The sediment incubation experiments demonstrated turnover rates in a similar range (vadose, 4.7 mg/kg*day; aquifer, 6.4 mghexadecane/kgsoil/day) of hexadecane degradation in the vadose zone and the aquifer, although mediated by slightly different microbial communities according to the analysis of fatty acid patterns and amounts. Additional experiments had the task of evaluating the degradation potential for the branched-chain alkane pristane (2,6,10,14-tetramethylpentadecane). Although this compound is regarded to be hardly degradable in comparison to n-alkanes and is thus frequently used as a reference parameter for indexing the extent of biodegradation of crude oils, it could be shown to be degraded by means of the incubation experiments. Thus, the site had a high inherent potential for natural attenuation of crude oils both in the vadose zone and the aquifer.

Numerical model investigation for potential methane explosion and benzene vapor intrusion associated with high-ethanol blend releases

Ethanol-blended fuel releases usually stimulate methanogenesis in the subsurface, which could pose an explosion risk if methane accumulates in a confined space above the ground where ignitable conditions exist. Ethanol-derived methane may also increase the vapor intrusion potential of toxic fuel hydrocarbons by stimulating the depletion of oxygen by methanotrophs, and thus inhibiting aerobic biodegradation of hydrocarbon vapors. To assess these processes, a three-dimensional numerical vapor intrusion model was used to simulate the degradation, migration, and intrusion pathway of methane and benzene under different site conditions. Simulations show that methane is unlikely to build up to pose an explosion hazard (5% v/v) if diffusion is the only mass transport mechanism through the deeper vadose zone. However, if methanogenic activity near the source zone is sufficiently high to cause advective gas transport, then the methane indoor concentration may exceed the flammable threshold under simulated conditions. During subsurface migration, methane biodegradation could consume soil oxygen that would otherwise be available to support hydrocarbon degradation, and increase the vapor intrusion potential for benzene. Vapor intrusion would also be exacerbated if methanogenic activity results in sufficiently high pressure to cause advective gas transport in the unsaturated zone. Overall, our simulations show that current approaches to manage the vapor intrusion risk for conventional fuel released might need to be modified when dealing with some high ethanol blend fuel (i.e., E20 up to E95) releases.

Numerical modeling analysis of VOC removal processes in different aerobic vertical flow systems for groundwater remediation

Vertical flow systems filled with porous medium have been shown to efficiently remove volatile organic contaminants (VOCs) from contaminated groundwater. To apply this semi-natural remediation strategy it is however necessary to distinguish between removal due to biodegradation and due to volatile losses to the atmosphere. Especially for (potentially) toxic VOCs, the latter needs to be minimized to limit atmospheric emissions. In this study, numerical simulation was used to investigate quantitatively the removal of volatile organic compounds in two pilot-scale water treatment systems: an unplanted vertical flow filter and a planted one, which could also be called a vertical flow constructed wetland, both used for the treatment of contaminated groundwater. These systems were intermittently loaded with contaminated water containing benzene and MTBE as main VOCs. The highly dynamic but permanently unsaturated conditions in the porous medium facilitated aerobic biodegradation but could lead to volatile emissions of the contaminants. Experimental data from porous material analyses, flow rate measurements, solute tracer and gas tracer test, as well as contaminant concentration measurements at the boundaries of the systems were used to constrain a numerical reactive transport modeling approach. Numerical simulations considered unsaturated water flow, transport of species in the aqueous and the gas phase as well as aerobic degradation processes, which made it possible to quantify the rates of biodegradation and volatile emissions and calculating their contribution to total contaminant removal. A range of degradation rates was determined using experimental results of both systems under two operation modes and validated by field data obtained at different operation modes applied to the filters. For both filters, simulations and experimental data point to high biodegradation rates, if the flow filters have had time to build up their removal capacity. For this case volatile emissions are negligible and total removal can be attributed to biodegradation, only. The simulation study thus supports the use of both of these vertical flow systems for the treatment of groundwater contaminated with VOCs and the use of reactive transport modeling for the assessment of VOCs removal and operation modes in these high performance treatment systems.

In situ vadose zone bioremediation

Contamination of the vadose zone with various pollutants is a world-wide problem, and often technical or economic constraints impose remediation without excavation. In situ bioremediation in the vadose zone by bioventing has become a standard remediation technology for light spilled petroleum products. In this review, focus is given on new in situ bioremediation strategies in the vadose zone targeting a variety of other pollutants such as perchlorate, nitrate, uranium, chromium, halogenated solvents, explosives and pesticides. The techniques for biostimulation of either oxidative or reductive degradation pathways are presented, and biotransformations to immobile pollutants are discussed in cases of non-degradable pollutants. Furthermore, research on natural attenuation in the vadose zone is presented.

Methane emissions and contaminant degradation rates at sites affected by accidental releases of denatured fuel-grade ethanol

The recent increase in the use of denatured fuel-grade ethanol (DFE) has enhanced the probability of its environmental release. Due to the highly labile nature of ethanol (EtOH), it is expected to rapidly biodegrade, increasing the potential for inducing methanogenic conditions in the subsurface. As environmental releases of DFE can be expected to occur at the ground surface or in the vadose zone (e.g., due to surficial spills from rail lines or tanker trucks and leaking underground storage tanks), the potential for methane (CH4) generation at DFE spill sites requires evaluation. An assessment is needed because high CH4 generation rates may lead to CH4 fluxes towards the ground surface, which is of particular concern if spills are located close to human habitation-related to concerns of soil vapor intrusion (SVI). This work demonstrates, for the first time, the measurement of surficial gas release rates at large volume DFE spill sites. Two study sites, near Cambria and Balaton, in MN are investigated. Total carbon emissions at the ground surface (summing carbon dioxide (CO2) and CH4 emissions) are used to quantify depth-integrated DFE degradation rates. Results from both sites demonstrate that substantial CO2 and CH4 emissions do occur-even years after a spill. However, large total carbon fluxes, and CH4 emissions in particular, were restricted to a localized area within the DFE source zone. At the Balaton site, estimates of total DFE carbon losses in the source zone ranged between 5 and 174 μmol m(-2) s(-1), and CH4 effluxes ranged between non-detect and 9 μmol m(-2) s(-1). At the Cambria site estimates of total DFE carbon losses in the source zone ranged between 8 and 500 μmol m(-2) s(-1), and CH4 effluxes ranged between non-detect and 393 μmol m(-2) s(-1). Substantial CH4 accumulation, coupled with oxygen (O2) depletion, measured in samples collected from custom-designed gas collection chambers at the Cambria site suggests that the development of explosion or asphyxiation hazards is possible in confined spaces above a rapidly degrading DFE release. However, the results also indicate that the development of such hazards is locally constrained, will require a high degree of soil moisture, close proximity to the source zone, a good connection between the soil and the confined space, and poorly aerated conditions.

Biodegradation of chlorobenzene, 1,2-dichlorobenzene, and 1,4-dichlorobenzene in the vadose zone

Much of the microbial activity in nature takes place at interfaces, which are often associated with redox discontinuities. One example is the oxic/anoxic interface where polluted groundwater interacts with the overlying vadose zone. We tested whether microbes in the vadose zone can use synthetic chemicals as electron donors and thus protect the overlying air and buildings from groundwater pollutants. Samples from the vadose zone of a site contaminated with chlorobenzene (CB), 1,2-dichlorobenzene (12DCB), and 1,4-dichlorobenzene (14DCB) were packed in a multiport column to simulate the interface of the vadose zone with an underlying groundwater plume. A mixture of CB, 12DCB, and 14DCB in anoxic water was pumped continuously through the bottom of column to an outlet below the first sampling port to create an oxic/anoxic interface and a capillary fringe. Removal to below the detection limits by rapid biodegradation with rates of 21 ± 1 mg of CB • m(-2) • d(-1), 3.7 ± 0.5 mg of 12DCB • m(-2) • d(-1), and 7.4 ± 0.7 mg of 1.4 DCB • m(-2) • d(-1) indicated that natural attenuation in the capillary fringe can prevent the migration of CB, 12DCB, and 14DCB vapors. Enumeration of bacteria capable of degrading chlorobenzenes suggested that most of the biodegradation takes place within the first 10 cm above the saturated zone. Biodegradation also increased the upward flux of contaminants and thus enhanced their elimination from the underlying water. The results revealed a substantial biodegradation capacity for chlorinated aromatic compounds at the oxic/anoxic interface and illustrate the role of microbes in creating steep redox gradients.

Role of natural attenuation in modeling the leaching of contaminants in the risk analysis framework

Natural attenuation (NA) processes occurring in the subsurface can significantly affect the impact on groundwater from contamination sources located in the vadose zone, especially when mobile and readily biodegradable compounds, such as BTEX, are present. Besides, in the last decades several studies have shown natural attenuation to take place also for more persistent compounds, such as Polycyclic Aromatic Hydrocarbons (PAHs). Nevertheless, common risk analysis frameworks, based on the ASTM RBCA (Risk Based Corrective Action) approach, do not include NA pathways in the fate and transport models, thus possibly leading to an overestimation of the calculated risk. The aim of this study was to provide an insight on the relevance of the different key natural attenuation processes usually taking place in the subsurface and to highlight for which contamination scenarios their inclusion in the risk-analysis framework could provide a more realistic risk assessment. To this end, an analytical model accounting for source depletion and biodegradation, dispersion and diffusion during leaching was developed and applied to several contamination scenarios. These scenarios included contamination by BTEX, characterized by relatively high mobility and biodegradation rate, and PAHs, i.e. a more persistent class of compounds. The obtained results showed that BTEX are likely to be attenuated in the source zone due to their mobility and ready biodegradation (assuming biodegradation constant rates in the order of 0.01-1 d(-1)). Instead, attenuation along transport through the vadose zone was found to be less important, as the residence time of the contaminant in the unsaturated zone is often too low with respect to the time required to get a relevant biodegradation of BTEX. On the other hand, heavier compounds such as PAHs, were found to be attenuated during leaching since the residence time in the vadose zone can reach values up to thousands of years. In these cases, even with the relatively slow biodegradation rate of PAHs, in the order of 0.0001-0.001 d(-1), attenuation can result significant. These conclusions were also confirmed by comparing the model results with experimental data collected at an hydrocarbon-contaminated site. The proposed model, that neglects the transport of NAPLs, could be easily included in the risk-analysis framework, allowing to get a more realistic assessment of risks, while keeping the intrinsic simplicity of the ASTM-RBCA approach.

Methane bioattenuation and implications for explosion risk reduction along the groundwater to soil surface pathway above a plume of dissolved ethanol

Fuel ethanol releases can stimulate methanogenesis in impacted aquifers, which could pose an explosion risk if methane migrates into enclosed spaces where ignitable conditions exist. To assess this potential risk, a flux chamber was emplaced on a pilot-scale aquifer exposed to continuous release (21 months) of an ethanol solution (10% v:v) that was introduced 22.5 cm below the water table. Despite methane concentrations within the ethanol plume reaching saturated levels (20-23 mg/L), the maximum methane concentration reaching the chamber (21 ppm(v)) was far below the lower explosion limit in air (50,000 ppm(v)). The low concentrations of methane observed in the chamber are attributed to methanotrophic activity, which was highest in the capillary fringe. This was indicated by methane degradation assays in microcosms prepared with soil samples from different depths, as well as by PCR measurements of pmoA, which is a widely used functional gene biomarker for methanotrophs. Simulations with the analytical vapor intrusion model "Biovapor" corroborated the low explosion risk associated with ethanol fuel releases under more generic conditions. Model simulations also indicated that depending on site-specific conditions, methane oxidation in the unsaturated zone could deplete the available oxygen and hinder aerobic benzene biodegradation, thus increasing benzene vapor intrusion potential. Overall, this study shows the importance of methanotrophic activity near the water table to attenuate methane generated from dissolved ethanol plumes and reduce its potential to migrate and accumulate at the surface.

Loss of volatile hydrocarbons from an LNAPL oil source

The light nonaqueous phase liquid (LNAPL) oil pool in an aquifer that resulted from a pipeline spill near Bemidji, Minnesota, was analyzed for volatile hydrocarbons (VHCs) to determine if the composition of the oil remains constant over time. Oil samples were obtained from wells at five locations in the oil pool in an anaerobic part of the glacial outwash aquifer. Samples covering a 21-year period were analyzed for 25 VHCs. Compared to the composition of oil from the pipeline source, VHCs identified in oil from wells sampled in 2008 were 13 to 64% depleted. The magnitude of loss for the VHCs analyzed was toluene≫o-xylene, benzene, C(6) and C(10-12)n-alkanes>C(7)-C(9)n-alkanes>m-xylene, cyclohexane, and 1- and 2-methylnaphthalene>1,2,4-trimethylbenzene and ethylbenzene. Other VHCs including p-xylene, 1,3,5- and 1,2,3-trimethylbenzenes, the tetramethylbenzenes, methyl- and ethyl-cyclohexane, and naphthalene were not depleted during the time of the study. Water-oil and air-water batch equilibration simulations indicate that volatilization and biodegradation is most important for the C(6)-C(9)n-alkanes and cyclohexanes; dissolution and biodegradation is important for most of the other hydrocarbons. Depletion of the hydrocarbons in the oil pool is controlled by: the lack of oxygen and nutrients, differing rates of recharge, and the spatial distribution of oil in the aquifer. The mass loss of these VHCs in the 5 wells is between 1.6 and 7.4% in 29years or an average annual loss of 0.06-0.26%/year. The present study shows that the composition of LNAPL changes over time and that these changes are spatially variable. This highlights the importance of characterizing the temporal and spatial variabilities of the source term in solute-transport models.

Modeling of vapor intrusion from hydrocarbon-contaminated sources accounting for aerobic and anaerobic biodegradation

A one-dimensional steady state vapor intrusion model including both anaerobic and oxygen-limited aerobic biodegradation was developed. The aerobic and anaerobic layer thickness are calculated by stoichiometrically coupling the reactive transport of vapors with oxygen transport and consumption. The model accounts for the different oxygen demand in the subsurface required to sustain the aerobic biodegradation of the compound(s) of concern and for the baseline soil oxygen respiration. In the case of anaerobic reaction under methanogenic conditions, the model accounts for the generation of methane which leads to a further oxygen demand, due to methane oxidation, in the aerobic zone. The model was solved analytically and applied, using representative parameter ranges and values, to identify under which site conditions the attenuation of hydrocarbons migrating into indoor environments is likely to be significant. Simulations were performed assuming a soil contaminated by toluene only, by a BTEX mixture, by Fresh Gasoline and by Weathered Gasoline. The obtained results have shown that for several site conditions oxygen concentration below the building is sufficient to sustain aerobic biodegradation. For these scenarios the aerobic biodegradation is the primary mechanism of attenuation, i.e. anaerobic contribution is negligible and a model accounting just for aerobic biodegradation can be used. On the contrary, in all cases where oxygen is not sufficient to sustain aerobic biodegradation alone (e.g. highly contaminated sources), anaerobic biodegradation can significantly contribute to the overall attenuation depending on the site specific conditions.

Nitrogen as an indicator of mass transfer during in-situ gas sparging

Aiming at the stimulation of intrinsic microbial activity, pulses of pure oxygen or pressurized air were recurrently injected into groundwater polluted with chlorobenzene. To achieve well-controlled conditions and intensive sampling, a large, vertical underground tank was filled with the local unconfined sandy aquifer material. In the course of two individual gas injections, one using pure oxygen and one using pressurized air, the mass transfer of individual gas species between trapped gas phase and groundwater was studied. Field data on the dissolved gas composition in the groundwater were combined with a kinetic model on gas dissolution and transport in porous media. Phase mass transfer of individual gas components caused a temporary enrichment of nitrogen, and to a lower degree of methane, in trapped gas leading to the formation of excess dissolved nitrogen levels downgradient from the dissolving gas phase. By applying a novel gas sampling method for dissolved gases in groundwater it was shown that dissolved nitrogen can be used as a partitioning tracer to indicate complete gas dissolution in porous media.

Review of unsaturated-zone transport and attenuation of volatile organic compound (VOC) plumes leached from shallow source zones

Reliable prediction of the unsaturated zone transport and attenuation of dissolved-phase VOC (volatile organic compound) plumes leached from shallow source zones is a complex, multi-process, environmental problem. It is an important problem as sources, which include solid-waste landfills, aqueous-phase liquid discharge lagoons and NAPL releases partially penetrating the unsaturated zone, may persist for decades. Natural attenuation processes operating in the unsaturated zone that, uniquely for VOCs includes volatilisation, may, however, serve to protect underlying groundwater and potentially reduce the need for expensive remedial actions. Review of the literature indicates that only a few studies have focused upon the overall leached VOC source and plume scenario as a whole. These are mostly modelling studies that often involve high strength, non-aqueous phase liquid (NAPL) sources for which density-induced and diffusive vapour transport is significant. Occasional dissolved-phase aromatic hydrocarbon controlled infiltration field studies also exist. Despite this lack of focus on the overall problem, a wide range of process-based unsaturated zone - VOC research has been conducted that may be collated to build good conceptual model understanding of the scenario, particularly for the much studied aromatic hydrocarbons and chlorinated aliphatic hydrocarbons (CAHs). In general, the former group is likely to be attenuated in the unsaturated zone due to their ready aerobic biodegradation, albeit with rate variability across the literature, whereas the fate of the latter is far less likely to be dominated by a single mechanism and dependent upon the relative importance of the various attenuation processes within individual site - VOC scenarios. Analytical and numerical modelling tools permit effective process representation of the whole scenario, albeit with potential for inclusion of additional processes - e.g., multi-mechanistic sorption phase partitioning, and provide good opportunity for further sensitivity analysis and development to practitioner use. There remains a significant need to obtain intermediate laboratory-scale and particularly field-scale (actual site and controlled release) datasets that address the scenario as a whole and permit validation of the available models. Integrated assessment of the range of simultaneous processes that combine to influence leached plume generation, transport and attenuation in the unsaturated zone is required. Component process research needs are required across the problem scenario and include: the simultaneous volatilisation and dissolution of source zones; development of appropriate field-scale dispersion estimates for the unsaturated zone; assessment of transient VOC exchanges between aqueous, vapour and sorbed phases and their influence upon plume attenuation; development of improved field methods to recognise and quantify biodegradation of CAHs; establishment of the influence of co-contaminants; and, finally, translation of research findings into more robust practitioner practice.

CO2-efflux measurements for evaluating source zone natural attenuation rates in a petroleum hydrocarbon contaminated aquifer

In order to gain regulatory approval for source zone natural attenuation (SZNA) at hydrocarbon-contaminated sites, knowledge regarding the extent of the contamination, its tendency to spread, and its longevity is required. However, reliable quantification of biodegradation rates, an important component of SZNA, remains a challenge. If the rate of CO(2) gas generation associated with contaminant degradation can be determined, it may be used as a proxy for the overall rate of subsurface biodegradation. Here, the CO(2)-efflux at the ground surface is measured using a dynamic closed chamber (DCC) method to evaluate whether this technique can be used to assess the areal extent of the contaminant source zone and the depth-integrated rate of contaminant mineralization. To this end, a field test was conducted at the Bemidji, MN, crude oil spill site. Results indicate that at the Bemidji site the CO(2)-efflux method is able to both delineate the source zone and distinguish between the rates of natural soil respiration and contaminant mineralization. The average CO(2)-efflux associated with contaminant degradation in the source zone is estimated at 2.6 μmol m(-2) s(-1), corresponding to a total petroleum hydrocarbon mineralization rate (expressed as C(10)H(22)) of 3.3 g m(-2) day(-1).