Biodegradation of hydrocarbons vapors: Comparison of laboratory studies and field investigations in the vadose zone at the emplaced fuel source experiment, Airbase Vaerløse, Denmark

The natural attenuation of volatile organic compounds (VOCs) in the unsaturated zone can only be predicted when information about microbial biodegradation rates and kinetics are known. This study aimed at determining first-order rate coefficients for the aerobic biodegradation of 13 volatile petroleum hydrocarbons which were artificially emplaced as a liquid mixture during a field experiment in an unsaturated sandy soil. Apparent first-order biodegradation rate coefficients were estimated by comparing the spatial evolution of the resulting vapor plumes to an analytical reactive transport model. Two independent reactive numerical model approaches have been used to simulate the diffusive migration of VOC vapors and to estimate degradation rate coefficients. Supplementary laboratory column and microcosm experiments were performed with the sandy soil at room temperature under aerobic conditions. First-order kinetics adequately matched the lab column profiles for most of the compounds. Consistent compound-specific apparent first-order rate coefficients were obtained by the three models and the lab column experiment, except for benzene. Laboratory microcosm experiments lacked of sensitivity for slowly degrading compounds and underestimated degradation rates by up to a factor of 5. Addition of NH3 vapor was shown to increase the degradation rates for some VOCs in the laboratory microcosms. All field models suggested a significantly higher degradation rate for benzene than the rates measured in the lab, suggesting that the field microbial community was superior in developing benzene degrading activity.     

Biodegradation of petroleum hydrocarbon vapors: laboratory studies on rates and kinetics in unsaturated alluvial sand

Predictions of natural attenuation of volatile organic compounds (VOCs) in the unsaturated zone rely critically on information about microbial biodegradation kinetics. This study aims at determining kinetic rate laws for the aerobic biodegradation of a mixture of 12 volatile petroleum hydrocarbons and methyl tert-butyl ether (MTBE) in unsaturated alluvial sand. Laboratory column and batch experiments were performed at room temperature under aerobic conditions, and a reactive transport model for VOC vapors in soil gas coupled to Monod-type degradation kinetics was used for data interpretation. In the column experiment, an acclimatization of 23 days took place before steady-state diffusive vapor transport through the horizontal column was achieved. Monod kinetic parameters Ks and vmax could be derived from the concentration profiles of toluene, m-xylene, n-octane, and n-hexane, because substrate saturation was approached with these compounds under the experimental conditions. The removal of cyclic alkanes, isooctane, and 1,2,4-trimethylbenzene followed first-order kinetics over the whole concentration range applied. MTBE, n-pentane, and chlorofluorocarbons (CFCs) were not visibly degraded. Batch experiments suggested first-order disappearance rate laws for all VOCs except n-octane, which decreased following zero-order kinetics in live batch experiments. For many compounds including MTBE, disappearance rates in abiotic batch experiments were as high as in live batches indicating sorption. It was concluded that the column approach is preferable for determining biodegradation rate parameters to be used in risk assessment models.     

Design and laboratory testing of a chamber device to measure total flux of volatile organic compounds from the unsaturated zone under natural conditions

To determine if an aquifer contaminated with volatile organic compounds (VOCs) has potential for natural remediation, all natural processes affecting the fate and transport of VOCs in the subsurface must be identified and quantified. This research addresses the quantification of air-phase volatile organic compounds (VOCs) leaving the unsaturated zone soil gas and entering the atmosphere-including the additional flux provided by advective soil-gas movement induced by barometric pumping. A simple and easy-to-use device for measuring VOC flux under natural conditions is presented. The vertical flux chamber (VFC) was designed using numerical simulations and evaluated in the laboratory. Mass-balance numerical simulations based on continuously stirred tank reactor equations (CSTR) provided information on flux measurement performance of several sampling configurations with the final chamber configuration measuring greater than 96% of model-simulated fluxes. A laboratory device was constructed to evaluate the flux chamber under both diffusion-only and advection-plus-diffusion transport conditions. The flux chamber measured an average of 82% of 15 diffusion-only fluxes and an average of 95% of 15 additional advection-plus-diffusion flux experiments. The vertical flux chamber has the capability of providing reliable measurement of VOC flux from the unsaturated zone under both diffusion and advection transport conditions.     

Variations in expression of carbon isotope fractionation of chlorinated ethenes during biologically enhanced PCE dissolution close to a source zone

The stable carbon isotope values of tetrachloroethene (PCE) and its degradation products were monitored during studies of biologically enhanced dissolution of PCE dense nonaqueous phase liquid (DNAPL) to determine the effect of PCE dissolution on observed isotope values. The degradation of PCE was monitored in a 2-dimensional model aquifer and in a pilot test cell (PTC) at Dover Air Force Base, both with emplaced PCE DNAPL sources. Within the plume down gradient from the source, the isotopic fractionation of dissolved PCE and its degradation products were consistent with those observed in biodegradation laboratory studies. However, close to the source zone significant shifts in the isotope values of dissolved PCE were not observed in either the model aquifer or PTC due to the constant input of newly dissolved, non fractionated PCE, and the small isotopic fractionation associated with PCE reductive dechlorination by the mixed microbial culture used. Therefore the identification of reductive dechlorination in the presence of PCE DNAPL was based upon the appearance of daughter products and the isotope values of those daughter products. An isotope model was developed to simulate isotope values of PCE during the dissolution and degradation of PCE adjacent to a DNAPL source zone. With the exception of very high degradation rate constants (>1/day) stable carbon isotope values of PCE estimated by the model remained within error of the isotope value of the PCE DNAPL, consistent with measured isotope values in the model aquifer and in the PTC.     

Combined application of conservative transport modelling and compound-specific carbon isotope analyses to assess in situ attenuation of benzene, toluene, and o-xylene

In recent years, compound specific isotope analyses (CSIA) have developed into one of the most powerful tools for the quantification of in situ biodegradation of organic contaminants. In this approach, the calculation of the extent of biodegradation of organic contaminants in aquifers is usually based on the Rayleigh equation, and thus neglects physical transport processes such as dispersion that contribute to contaminant dilution in aquifers. Here we combine compound specific isotope analyses with a conservative transport model to study the attenuation of aromatic hydrocarbons at a former gasworks site. The conservative transport model was first used to simulate concentration reductions caused by dilution at wells downgradient of a BTEX source. In a second step, the diluted concentrations, together with the available stable carbon isotope ratios and carbon fractionation factors for benzene, toluene and o-xylene were applied in the Rayleigh equation to quantify the degree of biodegradation at each of those wells. At the investigated site, where other attenuation processes such as sorption and volatilisation were proven to be negligible, the combined approach is recommended for benzene, which represents a compound for which the effect of biodegradation is comparable to or less than the effect of dilution. As demonstrated for toluene and o-xylene, the application of the Rayleigh equation alone is sufficient if dilution can be proved to be insignificant in comparison to biodegradation. The analysis also suggests that the source width and the position of the observation wells relative to the plume center line are significantly related to the degree of dilution.     

A large-scale experiment on mass transfer of trichloroethylene from the unsaturated zone of a sandy aquifer to its interfaces

A large-scale experiment was conducted to investigate the transport of trichloroethylene (TCE) vapors in the unsaturated zone and to determine the mass transfer to the groundwater and the atmosphere. The experiment involved injection of 5 1 of TCE in the unsaturated zone under controlled conditions, with multidepth sampling of gas and water through the unsaturated zone and across the capillary zone into underlying groundwater. The mass transfer of TCE vapors from the vadose zone to the atmosphere was quantified using a vertical flux chamber. A special soil water sampler was used to monitor transport across the capillary fringe. Experimental data indicated that TCE in the unsaturated zone was mainly transported to the atmosphere and this exchange reduced significantly the potential for groundwater pollution. The maximum measured TCE flux to the atmosphere was about 3 g/m(2)/day. Observed and calculated fluxes based on vertical TCE vapor concentration gradients and Fick's law were in good agreement. This confirms that TCE vapor transport under the experimental conditions was governed essentially by molecular diffusion. TCE vapors also caused a lower, but significant contamination of the underlying groundwater by dispersion across the capillary fringe with a corresponding maximum flux of about 0.1 g/m(2)/day. This mass transfer to groundwater is partly uncertain due to an inadvertent entry of some nonaqueous phase liquid (NAPL) from the source area into the saturated zone. Application of an analytical solution to estimate the TCE flux from the unsaturated zone to the groundwater indicated that this phenomenon is not only influenced by molecular diffusion but also by vertical dispersion. The mass balance indicates that, under the given experimental conditions (e.g. proximity of the source emplacement relative to the soil surface, relatively high permeable porous medium), nearly 95% of the initial TCE mass was transferred to the atmosphere.     

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.     

Microbial in situ degradation of aromatic hydrocarbons in a contaminated aquifer monitored by carbon isotope fractionation

We present an approach for characterizing in situ microbial degradation using the 13C/12C isotope fractionation of contaminants as an indicator of biodegradation. The 13C/12C isotope fractionation of aromatic hydrocarbons was studied in anoxic laboratory soil percolation columns with toluene or o-xylene as the sole carbon and electron source, and sulfate as electron acceptor. After approximately 2 months' of incubation, the soil microbial community degraded 32 mg toluene l(-1) and 44 mg o-xylene l(-1) to less than 0.05 mg l(-1), generating a stable concentration gradient in the column. The 13C/12C isotope ratio in the residual non-degraded fraction of toluene and o-xylene increased significantly, corresponding to isotope fractionation factors (alphaC) of 1.0015 and 1.0011, respectively. When the extent of biodegradation in the soil column was calculated based on the measured isotope ratios (R(t)) and an isotope fractionation factor (alphaC=1.0017) obtained from a sulfate-reducing batch culture the theoretical residual substrate concentrations (C(t)) matched the measured toluene concentrations in the column. This indicated that a calculation of biodegradation based on isotope fractionation could work in systems like soil columns. In a field study, a polluted, anoxic aquifer was analyzed for BTEX and PAH contaminants. These compounds were found to exhibit a significant concentration gradient along an 800-m groundwater flow path downstream of the source of contamination. A distinct increase in the carbon isotope ratio (delta13C) was observed for the residual non-degraded toluene (7.2 per thousand ), o-xylene (8.1 per thousand ) and naphthalene fractions (1.2 per thousand ). Based on the isotope values and the laboratory-derived isotope fractionation factors for toluene and o-xylene, the extent to which the residual substrate fraction in the monitoring wells had been degraded by microorganisms was calculated. The results revealed significant biodegradation along the groundwater flow path. In the wells at the end of the plume, the bioavailable toluene and o-xylene fractions had been almost completely reduced by in situ microbial degradation. Although indane and indene showed decreasing concentrations downstream of the groundwater flow path, suggesting microbial degradation, their carbon isotope ratios remained constant. As the physical properties of these compounds are similar to those of BTEX compounds, the constant isotope values of indane and indene indicated that microbial degradation did not lead to isotope fractionation of all aromatic hydrocarbons. In addition, physical interaction with the aquifer material during the groundwater passage did not significantly alter the carbon isotope composition of aromatic hydrocarbons.     

Transport and degradation of propylene glycol in the vadose zone: model development and sensitivity analysis

Transport and degradation of de-icing chemical (containing propylene glycol, PG) in the vadose zone were studied with a lysimeter experiment and a model, in which transient water flow, kinetic degradation of PG and soil chemistry were combined. The lysimeter experiment indicated that aerobic as well as anaerobic degradation occurs in the vadose zone. Therefore, the model included both types of degradation, which was made possible by assuming advection-controlled (mobile) and diffusion-controlled (immobile) zones. In the mobile zone, oxygen can be transported by diffusion in the gas phase. The immobile zone is always water-saturated, and oxygen only diffuses slowly in the water phase. Therefore, the model is designed in a way that the redox potential can decrease when PG is degraded, and thus, anaerobic degradation can occur. In our model, manganese oxide (MnO₂, which is present in the soil) and NO \(_{3}^{-}\) (applied to enhance biodegradation) can be used as electron acceptors for anaerobic degradation. The application of NO \(_{3}^{-}\) does not result in a lower leaching of PG nor in a slower depletion of MnO₂. The thickness of the snowcover influences the leached fraction of PG, as with a high infiltration rate, transport is fast, there is less time for degradation and thus more PG will leach. The model showed that, in this soil, the effect of the water flow dominates over the effect of the degradation parameters on the leaching at a 1-m depth.     

Reactive transport modeling of processes controlling the distribution and natural attenuation of phenolic compounds in a deep sandstone aquifer

Reactive solute transport modeling was utilized to evaluate the potential for natural attenuation of a contaminant plume containing phenolic compounds at a chemical producer in the West Midlands, UK. The reactive transport simulations consider microbially mediated biodegradation of the phenolic compounds (phenols, cresols, and xylenols) by multiple electron acceptors. Inorganic reactions including hydrolysis, aqueous complexation, dissolution of primary minerals, formation of secondary mineral phases, and ion exchange are considered. One-dimensional (1D) and three-dimensional (3D) simulations were conducted. Mass balance calculations indicate that biodegradation in the saturated zone has degraded approximately 1-5% of the organic contaminant plume over a time period of 47 years. Simulations indicate that denitrification is the most significant degradation process, accounting for approximately 50% of the organic contaminant removal, followed by sulfate reduction and fermentation reactions, each contributing 15-20%. Aerobic respiration accounts for less than 10% of the observed contaminant removal in the saturated zone. Although concentrations of Fe(III) and Mn(IV) mineral phases are high in the aquifer sediment, reductive dissolution is limited, producing only 5% of the observed mass loss. Mass balance calculations suggest that no more than 20-25% of the observed total inorganic carbon (TIC) was generated from biodegradation reactions in the saturated zone. Simulations indicate that aerobic biodegradation in the unsaturated zone, before the contaminant entered the aquifer, may have produced the majority of the TIC observed in the plume. Because long-term degradation is limited to processes within the saturated zone, use of observed TIC concentrations to predict the future natural attenuation may overestimate contaminant degradation by a factor of 4-5.     

In situ biodegradation determined by carbon isotope fractionation of aromatic hydrocarbons in an anaerobic landfill leachate plume (Vejen,Denmark)

Concentrations and isotopic compositions (13C/12C) of aromatic hydrocarbons were determined in eight samples obtained from the strongly anoxic part of the leachate plume downgradient from the Vejen Landfill (Denmark), where methanogenic, sulfate-reducing and iron-reducing conditions were observed. Despite the heterogeneous distribution of the compounds in the plume, the isotope fractionation proved that ethylbenzene and m/p-xylene were subject to significant biodegradation within the strongly anoxic plume. The isotope fractionation factors (alphaC) for the degradation of the m/p-xylene (1.0015) and ethylbenzene (1.0021) obtained from the field observations were similar to factors previously determined for the anaerobic degradation of toluene and o-xylene in laboratory experiments, and suggest that in situ biodegradation is one major process controlling the fate of these contaminants in this aquifer. The isotope fractionation determined for 1,2,4-trimethylbenzene and 2-ethyltoluene suggested in situ biodegradation; however, the isotopic composition did not correlate well with the respective concentration as expressed by the Rayleigh equation. Some other compounds (1,2,3-trimethylbenzene, o-xylene, naphthalene and fenchone) did not show significant enrichments in delta13C values along the flow path. The compound concentrations were too low for accurate isotope analyses of benzene, toluene, 1- and 2-methylnaphthalene, while interferences in the chromatography made it impossible to evaluate the isotopic composition for 4-ethyltoluene, 1,3,5-trimethylbenzene and camphor.In addition to demonstrating the potential of assessing isotopic fractionation as a means for documenting the in situ biodegradation of complex mixtures of aromatic hydrocarbons in leachate plumes, this study also illustrates the difficulties for data interpretation in complex plumes and high analytical uncertainties for isotope analysis of organic compounds in low concentration ranges.     

Stable isotope fractionation analysis as a tool to monitor biodegradation in contaminated acquifers

The assessment of biodegradation in contaminated aquifers has become an issue of increasing importance in the recent years. To some extent, this can be related to the acceptance of intrinsic bioremediation or monitored natural attenuation as a means to manage contaminated sites. Among the few existing methods to detect biodegradation in the subsurface, stable isotope fractionation analysis (SIFA) is one of the most promising approaches which is pronounced by the drastically increasing number of applications. This review covers the recent laboratory and field studies assessing biodegradation of contaminants via stable isotope analysis. Stable isotope enrichment factors have been found that vary from no fractionation for dioxygenase reactions converting aromatic hydrocarbons over moderate fractionation by monooxygenase reactions (epsilon=-3 per thousand) and some anaerobic studies on microbial degradation of aromatic hydrocarbons (epsilon=-1.7 per thousand) to larger fractionations by anaerobic dehalogenation reactions of chlorinated solvents (epsilon=between -5 per thousand and -30 per thousand). The different isotope enrichment factors can be related to the respective biochemical reactions. Based on that knowledge, we discuss under what circumstances SIFA can be used for a qualitative or even a quantitative assessment of biodegradation in the environment. In a steadily increasing number of cases, it was possible to explain biodegradation processes in the field based on isotope enrichment factors obtained from laboratory experiments with pure cultures and measured isotope values from the field. The review will focus on the aerobic and anaerobic degradation of aromatic hydrocarbons and chlorinated solvents as the major contaminants of groundwater. Advances in the instrumental development for stable isotope analysis are only mentioned if it is important for the understanding of the application.     

Hydrochemical and isotopic effects associated with petroleum fuel biodegradation pathways in a chalk aquifer

Hydrochemical data, compound specific carbon isotope analysis and isotopic enrichment trends in dissolved hydrocarbons and residual electron acceptors have been used to deduce BTEX and MTBE degradation pathways in a fractured chalk aquifer. BTEX compounds are mineralised sequentially within specific redox environments, with changes in electron acceptor utilisation being defined by the exhaustion of specific BTEX components. A zone of oxygen and nitrate exhaustion extends approximately 100 m downstream from the plume source, with residual sulphate, toluene, ethylbenzene and xylene. Within this zone complete removal of the TEX components occurs by bacterial sulphate reduction, with sulphur and oxygen isotopic enrichment of residual sulphate (epsilon(s) = -14.4 per thousand to -16.0 per thousand). Towards the plume margins and at greater distance along the plume flow path nitrate concentrations increase with delta15N values of up to +40 per thousand indicating extensive denitrification. Benzene and MTBE persist into the denitrification zone, with carbon isotope enrichment of benzene indicating biodegradation along the flow path. A Rayleigh kinetic isotope enrichment model for 13C-enrichment of residual benzene gives an apparent epsilon value of -0.66 per thousand. MTBE shows no significant isotopic enrichment (delta13C = -29.3 per thousand to -30.7 per thousand) and is isotopically similar to a refinery sample (delta13C = -30.1 per thousand). No significant isotopic variation in dissolved MTBE implies that either the magnitude of any biodegradation-induced isotopic fractionation is small, or that relatively little degradation has taken place in the presence of BTEX hydrocarbons. It is possible, however, that MTBE degradation occurs under aerobic conditions in the absence of BTEX since no groundwater samples were taken with co-existing MTBE and oxygen. Low benzene delta13C values are correlated with high sulphate delta34S, indicating that little benzene degradation has occurred in the sulphate reduction zone. Benzene degradation may be associated with denitrification since increased benzene delta13C is associated with increased delta15N in residual nitrate. Re-supply of electron acceptors by diffusion from the matrix into fractures and dispersive mixing is an important constraint on degradation rates and natural attenuation capacity in this dual-porosity aquifer.     

Methanotrophic activity in a diffusive methane/oxygen counter-gradient in an unsaturated porous medium

Microbial methane (CH4) oxidation is a main control on emissions of this important greenhouse gas from ecosystems such as contaminated aquifers or wetlands under aerobic onditions. Due to a lack of suitable model systems, we designed a laboratory column to study this process in diffusional CH4/O2 counter-gradients in unsaturated porous media. Analysis and simulations of the steady-state CH4, CO2 and O2 gas profiles showed that in a 15-cm-deep active zone, CH4 oxidation followed first-order kinetics with respect to CH4 with a high apparent first-order rate constant of approximately 30 h(-1). Total cell counts obtained using DAPI-staining suggested growth of methanotrophic bacteria, resulting in a high capacity for CH4 oxidation. This together with apparent tolerance to anoxic conditions enabled a rapid response of the methanotrophic community to changing substrate availability, which was induced by changes in O2 concentrations at the top of the column. Microbial oxidation was confirmed by a approximately 7 per thousand enrichment in CH4 stable carbon isotope ratios along profiles. Using a fractionation factor of 1.025+/-0.0005 for microbial oxidation estimated from this shift and the fractionation factor for diffusion, simulations of isotope profiles agreed well with measured data confirming large fractionation associated with microbial oxidation. The designed column should be valuable for investigating response of methanotrophic bacteria to environmental parameters in future studies.     

Tracer test for the measurement of gas diffusion and non-aqueous phase liquid (NAPL) saturation in soil

During soil bioremediation, the diffusion of oxygen into the soil is an important prerequisite for aerobic biodegradation, and the decrease of petroleum products is the ultimate goal. Both processes need to be monitored. The aim of this work was to develop a gas tracer test that yields information on both, gas diffusion and residual saturation with non-aqueous phase liquids (NAPLs) in unsaturated soil heaps. One conservative tracer (methane) and 4 partitioning gas tracers (diethylether, methyl tert-butyl ether, chloroform and n-heptane) were injected as vapors into laboratory columns filled with unsaturated sand with increasing NAPL saturation. Breakthrough curves of gaseous compounds were measured at two points and compared to analytical solutions of an analytical diffusive-reactive transport equation. By fitting of methane data, robust results for effective diffusivity (tortuosity) were obtained. NAPL saturation was most accurately measured by the moderately water soluble tracers (ethers and chloroform). The hydrophobic tracer n-heptane did not partition into water-immersed NAPL. An easy and accurate way to assess air-NAPL partitioning constants from gas chromatography retention times is furthermore reported. It is concluded that gas tracer tests have the potential for measuring two important properties in soil bioremediation systems easily and quickly.     

Modeling of air sparging of VOC-contaminated soil columns

Air sparging is a remediation technology currently being applied for the restoration of sites contaminated with volatile organic compounds (VOCs). Attempts have been made by various researchers to model the fate of VOCs in the gas and liquid phase during air sparging. In this study, a radial diffusion model with an air–water mass transfer boundary condition was developed and applied for the prediction of VOC volatilization from air sparging of contaminated soil columns. The approach taken was to use various parameters such as mass transfer coefficients and tortuosity factors determined previously in separate experiments using a single air channel apparatus and applying these parameters to a complex system with many air channels. Incorporated in the model, is the concept of mass transfer zone (MTZ) where diffusion of VOCs in this zone was impacted by the volatilization of VOCs at the air–water interface but with negligible impact outside the zone. The model predicted fairly well the change in the VOC concentrations in the exhaust air, the final average aqueous VOC concentration, and the total mass removed. The predicted mass removal was within 1% to 20% of the actual experimental mass removed. The results of the model seemed to suggest that air-sparged soil columns may be modeled as a composite of individual air channels surrounded by a MTZ. For a given air flow rate and air saturation, the VOC removal was found to be inversely proportional to the radius of the air channel. The approach taken provided conceptual insights on mass transfer processes during air sparging operations.     

Oxygenated gasoline release in the unsaturated zone, Part 2: Downgradient transport of ethanol and hydrocarbons

In the event of a gasoline spill containing oxygenated compounds such as ethanol and MTBE, it is important to consider the impacts these compounds might have on subsurface contamination. One of the main concerns commonly associated with ethanol is that it might decrease the biodegradation of aromatic hydrocarbon compounds, leading to an increase in the hydrocarbon dissolved plume lengths. The first part of this study (Part 1) showed that when gasoline containing ethanol infiltrates the unsaturated zone, ethanol is likely to partition to and be retained in the unsaturated zone pore water. In this study (Part 2), a controlled field test is combined with a two-dimensional laboratory test and three-dimensional numerical modelling to investigate how ethanol retention in the unsaturated zone affects the downgradient behaviour of ethanol and aromatic hydrocarbon compounds. Ethanol transport downgradient was extremely limited. The appearance of ethanol in downgradient wells was delayed and the concentrations were lower than would be expected based on equilibrium dissolution. Oscillations in the water table resulted in minor flushing of ethanol, but its effect could still be perceived as an increase in the groundwater concentrations downgradient from the source zone. Ethanol partitioning to the unsaturated zone pore water reduced its mass fraction within the NAPL thus reducing its anticipated impact on the fate of the hydrocarbon compounds. A conceptual numerical simulation indicated that the potential ethanol-induced increase in benzene plume length after 20 years could decrease from 136% to 40% when ethanol retention in the unsaturated zone is considered.     

An analytical solution for vertical transport of volatile chemicals in the vadose zone

An analytical solution is presented for one-dimensional vertical transport of volatile chemicals through the vadose zone to groundwater. The solution accounts for the important transport mechanisms of the steady advection of water and gas, diffusion and dispersion in water and gas, as well as adsorption, and first-order degradation. By assuming a linear, equilibrium partitioning between water, gas and the adsorbed chemical phases, the dependent variable in the mathematical model becomes the total resident concentration. The general solution was derived for cases having a constant initial total concentration over a discrete depth interval and a zero initial total concentration elsewhere. A zero concentration gradient is assumed at the groundwater table. Examples are given to demonstrate the application of the new solution for calculating the case of a non-uniform initial source concentration, and estimating the transport of chemicals to the groundwater and the atmosphere. The solution was also used to verify a numerical code called VLEACH. We discovered an error in VLEACH, and found that the new solution agreed very well with the numerical results by corrected VLEACH. A simplified solution to predict the migration of volatile organic chemical due to the gas density effect has shown that a high source concentration may lead to significant downward advective gas-phase transport in a soil with a high air-permeability.     

Investigation of hydrogeologic processes in a dipping layer structure: 2. Transport and biodegradation of organics

Numerical simulation tools have been used to study the dominating processes during transport of aromatic hydrocarbons in the unsaturated soil zone. Simulations were based on field observations at an experimental site located on a glacial delta plain with pronounced layered sedimentary structures. A numerical model for transport in the unsaturated zone, SWMS-3D, has been extended to incorporate coupled multispecies transport, microbial degradation following Monod kinetics and gas diffusive transport of oxygen and hydrocarbons. The flow field parameters were derived from previous work using nonreactive tracers. Breakthrough curves (BTC) from the hydrocarbon field experiment were used to determine sorption parameters and Monod kinetic parameters using a fitting procedure. The numerical simulations revealed that the assumption of homogeneous layers resulted in deviations from the field observations. The deviations were more pronounced with incorporation of reactive transport, compared with earlier work on nonreactive transport. To be able to model reasonable BTC, sorption had to be reduced compared to laboratory experiments. The initial biomass and the maximum utilisation rate could be adjusted to capture both the initial lag phase and the overall degradation rate. Nevertheless, local oxygen limitation is predicted by the model, which was not observed in the field experiment. Incorporation of evaporation and diffusive gas transport of the hydrocarbons did not significantly change the local oxygen demand. The main cause of the observed discrepancies between model and field are attributed to channelling as a result of small-scale heterogeneities such as biopores.