Transport of hydrocarbons from an emplaced fuel source experiment in the vadose zone at Airbase Vaerløse, Denmark

Evaluation of seasonal factors and rainfall on temporal variation of molecular composition of two distinct hydrocarbon soil vapor plumes

Seasonal changes and significant rainfall event can significantly alter temperatures and rainfall infiltration which greatly affects the vadose zone gas composition. This study reports a high-resolution analysis of temporal variations in vadose zone gas compositions of two distinct vapor plumes emanating from two different light non-aqueous phase liquid (LNAPL) sources (aliphatic-rich naphtha for Zone #1vs aromatic-rich pyrolysis gasoline for Zone #2) at a large petrochemical site. Increase in soil moisture caused by significant rainfall reduced the influx of atmospheric oxygen into the soil layer and inhibited aerobic biodegradation of hydrocarbon vapor, thus resulting in the accumulation of VOCs and CH4 in the shallow soil vapor. The impacts of temperature and rainfall infiltration on the molecular composition of shallow soil vapor were significantly greater than that on deep soil vapor. During the hot and rainy summer, the O2 concentration in the shallow soil gas was at its lowest level, while the detection rate of VOCs in shallow soil gas was the highest. The largest volume (9.20 × 104-9.36 × 104 m3 for Zone #1 and 1.30 × 105-1.37 × 105 m3 for Zone #2) and mass (1181-1222 kg for Zone #1 and 221-237 kg for Zone #2) of vapor plume were observed during the summer. The concentration and composition of shallow soil vapor fluctuate due to seasonal factors, but the diagnostic ratios of conservative markers "alkylcyclopentanes" consistently maintained the distinctive chemical characteristics of the source LNAPL in all sampling campaigns. This field data offers important implications for practitioners regarding the timing of site investigation and vapor intrusion risk assessment.

Influence of dissolved and non-aqueous phase toluene on spectral induced polarization signatures of soils

Fifty-two laboratory experiments are undertaken to analyze the sensitivity of spectral induced polarization (SIP) to the presence of toluene in soils. Among these experiments, four experiments are conducted to collect SIP responses of soils containing dissolved phase toluene within the pore water using columns. The results demonstrate that SIP is not sensitive to the presence of dissolved phase toluene in soils. The remaining forty-eight experiments are undertaken with four types of soils mixed with non-aqueous phase toluene. The experimental results prove that SIP is sensitive to toluene saturation under varying salinity conditions. These observations are well-explained by a published petrophysical model accounting for the effects of water saturation on complex conductivity. The water saturation exponent n and quadrature conductivity exponent p in this model are obtained by fitting complex conductivity data versus saturation at different saturation levels. The petrophysical model is tested where in-phase and quadrature conductivity responses are predicted from water saturation, soil cation exchange capacity (CEC), and pore water conductivity. The petrophysical model provides satisfactory predictions for non-aqueous phase toluene saturation. Overall, this study contributes to our understanding of SIP as a non-intrusive tool for characterizing toluene contamination in soils with applications to the field.

Understanding petroleum vapor fate and transport through high resolution analysis of two distinct vapor plumes

No field study has provided a detailed characterization of the molecular composition and spatial distribution of a vadose zone plume of petroleum volatile organic compounds (VOCs), which is critical to improve the current understanding of petroleum VOC transport and fate. This is study reports a high-resolution analysis of two distinct vapor plumes emanating from two different light non-aqueous phase liquid (LNAPL) sources (an aliphatic-rich LNAPL for Zone #1vs an aromatic-rich LNAPL for Zone #2) at a large petrochemical site. Although deep soil vapor signatures were similar to the source zone LNAPL signatures, the composition of the shallow soil vapors reflected preferential attenuation of certain hydrocarbons over others during upward transport in the vadose zone. Between deeper and shallower soil gas samples, attenuation of aromatics was observed under all conditions, but important differences were observed in attenuation to aliphatic compound classes. Attenuation of all aliphatic compounds was observed under aerobic conditions but little attenuation of any aliphatics was observed under anoxic conditions without methane. In contrast, under methanogenic conditions, paraffins attenuated more than isoparaffins and naphthenes. These results suggest that isoparafins and naphthenes may present more of a vapor intrusion risk than benzene or other aromatic hydrocarbons commonly considered to be petroleum vapor intrusion risk drivers. While the overall vapor composition changed significantly within the vadose zone, diagnostic ratios of relatively recalcitrant alkylcyclopentanes were preserved in shallow soil vapor samples. These alkylcyclopentanes may be useful for distinguishing between petroleum vapor intrusion and other sources of petroleum VOCs detected in indoor air.

Transport and natural attenuation of benzene vapor from a point source in the vadose zone

The vadose zone is a very dynamic and active environment that directly affects natural attenuation and vapor intrusion of volatile organic compounds (VOCs). Therefore, it is important to understand the fate and transport of VOCs in the vadose zone. A column experiment combined with model study was conducted to investigate the influence of soil type, vadose zone thickness, and soil moisture content on benzene vapor transport and natural attenuation in the vadose zone. Vapor-phase biodegradation and volatilization to atmosphere for benzene are two main natural attenuation mechanism in the vadose zone. Our data showed that biodegradation in black soil is the main natural attenuation mechanism (82.8%) while volatilization is the main natural attenuation mechanism in quartz sand, floodplain soil, lateritic red earth and yellow earth (>71.9%). The R-UNSAT model-predicted soil gas concentration profile and flux were close with four soil column data except for yellow earth. Increasing the vadose zone thickness and soil moisture content significantly reduced the volatilization contribution while increased biodegradation contribution. The volatilization loss decreased from 89.3% to 45.8% when the vadose zone thickness increased from 30 cm to 150 cm. The volatilization loss decreased from 71.9% to 10.1% when the soil moisture content increased from 6.4% to 25.4%. Overall, this study provided valuable insights into clarifying the roles of soil type, moisture, and other environmental conditions in vadose zone natural attenuation mechanism and vapor concentration.

Assessing toluene biodegradation under temporally varying redox conditions in a fractured bedrock aquifer using stable isotope methods

In complex hydrogeological settings little is known about the extent of temporally varying redox conditions and their effect on aromatic hydrocarbon biodegradation. This study aims to assess the impact of changing redox conditions over time on aromatic hydrocarbon biodegradation in a fractured bedrock aquifer using stable isotope methods. To that end, four snapshots of highly spatio-temporally resolved contaminant and redox sensitive species concentrations, as well as stable isotope ratio profiles, were determined over a two-years time period in summer 2016, spring 2017, fall 2017 and summer 2018 in a toluene contaminated fractured bedrock aquifer. The concentration profiles of redox sensitive species and stable isotope ratio profiles for dissolved inorganic carbon (DIC) and sulfate (δ¹³C_DIC, δ³⁴S_SO4, δ¹⁸O_SO4) revealed that the aquifer alternates between oxidising (spring 2017/summer 2018) and reducing conditions (summer 2016/fall 2017). This alternation was attributed to a stronger aquifer recharge with oxygen-rich meltwater in spring 2017/summer 2018 compared to summer 2016/fall 2017. The temporally varying redox conditions coincided with various extents of toluene biodegradation revealed by the different magnitude of heavy carbon (¹³C) and hydrogen (²H) isotope enrichment in toluene. This indicated that the extent of toluene biodegradation and its contribution to plume attenuation was controlled by the temporally changing redox conditions. The highest toluene biodegradation was observed in summer 2016, followed by spring 2017 and fall 2017, whereby these temporal changes in biodegradation occurred throughout the whole plume. Thus, under temporally varying recharge conditions both the core and the fringe of a contaminant plume can be replenished with terminal electron acceptors causing biodegradation in the whole plume and not only at its distal end as previously suggested by the plume fringe concept. Overall, this study highlights the importance of highly temporally resolved groundwater monitoring to capture temporally varying biodegradation rates and to accurately predict biodegradation-induced contaminant attenuation in fractured bedrock aquifers.

A simple contaminant fate and transport modelling tool for management and risk assessment of groundwater pollution from contaminated sites

Contaminated sites pose a significant threat to groundwater resources. The resources that can be allocated by water regulators for site investigation and cleanup are limited compared to the large number of contaminated sites. Numerical transport models of individual sites require large amounts of data and are labor intensive to set up, and thus they are likely to be too expensive to be useful in the management of thousands of contaminated sites. Therefore, simple tools based on analytical solutions of contaminant transport models are widely used to assess (at an early stage) whether a site might pose a threat to groundwater. We present a tool consisting of five different models, representing common geological settings, contaminant pathways, and transport processes. The tool employs a simplified approach for preliminary, conservative, fast and inexpensive estimation of the contamination levels of aquifers. This is useful for risk assessment applications or to select and prioritize the sites, which should be targeted for further investigation. The tool is based on steady-state semi-analytical models simulating different contaminant transport scenarios from the source to downstream groundwater, and includes both unsaturated and saturated transport processes. The models combine existing analytical solutions from the literature for vertical (from the source to the top of the aquifer) and horizontal (within the aquifer) transport. The effect of net recharge causing a downward migration and an increase of vertical dispersion and dilution of the plume is also considered. Finally, we illustrate the application of the tool for a preliminary assessment of two contaminated sites in Denmark and compare the model results with field data. The comparison shows that a first preliminary assessment with conservative, and often non-site specific parameter selection, is qualitatively consistent with broad trends in observations and provides a conservative estimate of contamination.

Toluene biodegradation in the vadose zone of a poplar phytoremediation system identified using metagenomics and toluene-specific stable carbon isotope analysis

Biodegradation is an important mechanism of action of phytoremediation systems, but performance evaluation is challenging. We applied metagenomic molecular approaches and compound-specific stable carbon isotope analysis to assess biodegradation of toluene in the vadose zone at an urban pilot field system where hybrid poplars were planted to remediate legacy impacts to an underlying shallow fractured bedrock aquifer. Carbon isotope ratios were compared spatio-temporally between toluene dissolved in groundwater and in the vapor phase. Enrichment of ¹³C from toluene in the vapor phase compared to groundwater provided evidence for biodegradation in the vadose zone. Total bacterial abundance (16S rRNA) and abundance and expression of degradation genes were determined in rhizosphere soil (DNA and RNA) and roots (DNA) using quantitative PCR. Relative abundances of degraders in the rhizosphere were on average higher at greater depths, except for enrichment of PHE-encoding communities that more strongly followed patterns of toluene concentrations detected. Quantification of RMO and PHE gene transcripts supported observations of active aerobic toluene degradation. Finally, spatially-variable numbers of toluene degraders were detected in poplar roots. We present multiple lines of evidence for biodegradation in the vadose zone at this site, contributing to our understanding of mechanisms of action of the phytoremediation system.

Role of the source to building lateral separation distance in petroleum vapor intrusion

The adoption of source to building separation distances to screen sites that need further field investigation is becoming a common practice for the evaluation of the vapor intrusion pathway at sites contaminated by petroleum hydrocarbons. Namely, for the source to building vertical distance, the screening criteria for petroleum vapor intrusion have been deeply investigated in the recent literature and fully addressed in the recent guidelines issued by ITRC and U.S.EPA. Conversely, due to the lack of field and modeling studies, the source to building lateral distance received relatively low attention. To address this issue, in this work we present a steady-state vapor intrusion analytical model incorporating a piecewise first-order aerobic biodegradation limited by oxygen availability that accounts for lateral source to building separation. The developed model can be used to evaluate the role and relevance of lateral vapor attenuation as well as to provide a site-specific assessment of the lateral screening distances needed to attenuate vapor concentrations to risk-based values. The simulation outcomes showed to be consistent with field data and 3-D numerical modeling results reported in previous studies and, for shallow sources, with the screening criteria recommended by U.S.EPA for the vertical separation distance. Indeed, although petroleum vapors can cover maximum lateral distances up to 25-30m, as highlighted by the comparison of model outputs with field evidences of vapor migration in the subsurface, simulation results by this new model indicated that, regardless of the source concentration and depth, 6m and 7m lateral distances are sufficient to attenuate petroleum vapors below risk-based values for groundwater and soil sources, respectively. However, for deep sources (>5m) and for low to moderate source concentrations (benzene concentrations lower than 5mg/L in groundwater and 0.5mg/kg in soil) the above criteria were found extremely conservative as the model results indicated that for such scenarios the lateral screening distance may be set equal to zero.

Characterization and fingerprinting of soil and groundwater contamination sources around a fuel distribution station in Galicia (NW Spain)

Soil and groundwater contamination around a fuel distribution station in Tomiño (NW Spain) was evaluated. For this purpose, top and subsoil (up to 6.4 m) and groundwater were sampled around the station, approximately in a 60-m radius. Samples were analysed by HS-SPME-GC-MS to identify and quantify volatile fuel organic compounds (VFOC) (MTBE, ETBE and BTEX) and diesel range organics (DRO). Analysis and fingerprinting data suggested that the contamination of soil and groundwater was provoked by a fuel leak from underground storage tanks. This was reflected by hydrocarbon indices and principal component analysis, which discriminated a direct source of contamination of the subsoil samples around the station. The contaminants probably migrated from tank nearby soils to surrounding soils and leached to groundwater, following a SW direction. Irrigation with contaminated groundwater provoked a severe contamination of topsoils, which were enriched with the lightest components of gasoline and diesel. Fingerprinting also revealed the continuity of the leak, reflected by the presence of volatiles in some samples, which principally appeared in fresh leaks. MTBE was detected in a very high concentration in groundwater samples (up to 690 μg L(-1)), but it was not detected in fresh gasoline. This also evidenced an old source of contamination, probably starting in the mid-1990s, when the use of MTBE in gasoline was regulated.

Advanced multivariate analysis to assess remediation of hydrocarbons in soils

Accurate monitoring of degradation levels in soils is essential in order to understand and achieve complete degradation of petroleum hydrocarbons in contaminated soils. We aimed to develop the use of multivariate methods for the monitoring of biodegradation of diesel in soils and to determine if diesel contaminated soils could be remediated to a chemical composition similar to that of an uncontaminated soil. An incubation experiment was set up with three contrasting soil types. Each soil was exposed to diesel at varying stages of degradation and then analysed for key hydrocarbons throughout 161 days of incubation. Hydrocarbon distributions were analysed by Principal Coordinate Analysis and similar samples grouped by cluster analysis. Variation and differences between samples were determined using permutational multivariate analysis of variance. It was found that all soils followed trajectories approaching the chemical composition of the unpolluted soil. Some contaminated soils were no longer significantly different to that of uncontaminated soil after 161 days of incubation. The use of cluster analysis allows the assignment of a percentage chemical similarity of a diesel contaminated soil to an uncontaminated soil sample. This will aid in the monitoring of hydrocarbon contaminated sites and the establishment of potential endpoints for successful remediation.

Can soil gas VOCs be related to groundwater plumes based on their isotope signature?

The isotope evolution of tetrachloroethene (PCE) during its transport from groundwater toward the soil surface was investigated using laboratory studies and numerical modeling. During air-water partitioning, carbon and chlorine isotope ratios evolved in opposite directions, with a normal isotope effect for chlorine (ε = -0.20‰) and an inverse effect for carbon (ε = +0.46‰). During the migration of PCE from groundwater to the unsaturated zone in a 2D laboratory system, small shifts of carbon and chlorine isotope ratios (+0.8‰) were observed across the capillary fringe. Numerical modeling showed that these shifts are due to isotope fractionation associated with air-water partitioning and gas-phase diffusion. Carbon and chlorine isotope profiles were constant throughout the unsaturated zone once a steady state was reached. However, depending on the thickness of the unsaturated zone and its lithology, depletion in heavy isotopes may occur with distance during the transient migration of contaminants. Additionally, variations of up to +1.5‰ were observed in the unsaturated zone for chlorine isotopes during water table fluctuations. However, at steady state, it is possible to link a groundwater plume to gas-phase contamination and/or to differentiate sources of contamination based on isotope ratios.

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.

Modeling volatilization of residual VOCs in unsaturated zones: a moving boundary problem

It is of practical interest to investigate the natural evaporation of volatile organic compounds (VOCs) after the removal of a leaking tank situated on the top of the soil. This study aims to develop a mathematical model to predict mole fraction distributions and migration of evaporation front for two VOCs emanating from residual non-aqueous phase liquid (NAPL) due to the leak from the tank in a homogeneous soil. Considering the location of the front and the regions above and below the front, a numerical model for the diffusive transport of VOCs in unsaturated soils was developed using the finite difference method with a moving grid approach. The model was further simplified to the case of single VOC and solved analytically by Boltzmann's transformation with a moving boundary. Analytical expressions for the depth and moving speed of the front for a single VOC were then obtained for practical use. Finally, the developed model was used to predict the concentration distributions of VOCs below the land surface and examine the factors affecting the location and moving speed of the evaporation front.

Photoacoustic infrared spectroscopy for conducting gas tracer tests and measuring water saturations in landfills

Gas tracer tests can be used to determine gas flow patterns within landfills, quantify volatile contaminant residence time, and measure water within refuse. While gas chromatography (GC) has been traditionally used to analyze gas tracers in refuse, photoacoustic spectroscopy (PAS) might allow real-time measurements with reduced personnel costs and greater mobility and ease of use. Laboratory and field experiments were conducted to evaluate the efficacy of PAS for conducting gas tracer tests in landfills. Two tracer gases, difluoromethane (DFM) and sulfur hexafluoride (SF(6)), were measured with a commercial PAS instrument. Relative measurement errors were invariant with tracer concentration but influenced by background gas: errors were 1-3% in landfill gas but 4-5% in air. Two partitioning gas tracer tests were conducted in an aerobic landfill, and limits of detection (LODs) were 3-4 times larger for DFM with PAS versus GC due to temporal changes in background signals. While higher LODs can be compensated by injecting larger tracer mass, changes in background signals increased the uncertainty in measured water saturations by up to 25% over comparable GC methods. PAS has distinct advantages over GC with respect to personnel costs and ease of use, although for field applications GC analyses of select samples are recommended to quantify instrument interferences.

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 growth with vapor-phase substrate

Limited information exists on influences of the diffusive transport of volatile organic contaminants (VOC) on bacterial activity in the unsaturated zone of the terrestrial subsurface. Diffusion of VOC in the vapor-phase is much more efficient than in water and results in effective VOC transport and high bioavailability despite restricted mobility of bacteria in the vadose zone. Since many bacteria tend to accumulate at solid-water, solid-air and air-water interfaces, such phase boundaries are of a special interest for VOC-biodegradation. In an attempt to evaluate microbial activity toward air-borne substrates, this study investigated the spatio-temporal interplay between growth of Pseudomonas putida (NAH7) on vapor-phase naphthalene (NAPH) and its repercussion on vapor-phase NAPH concentrations. Our data demonstrate that growth rates of strain PpG7 were inversely correlated to the distance from the source of vapor-phase NAPH. Despite the high gas phase diffusivity of NAPH, microbial growth was absent at distances above 5 cm from the source when sufficient biomass was located in between. This indicates a high efficiency of suspended bacteria to acquire vapor-phase compounds and influence headspace concentration gradients at the centimeter-scale. It further suggests a crucial role of microorganisms as biofilters for gas-phase VOC emanating from contaminated groundwater or soil.

Analytical modelling of stable isotope fractionation of volatile organic compounds in the unsaturated zone

Analytical models were developed that simulate stable isotope ratios of volatile organic compounds (VOCs) near a point source contamination in the unsaturated zone. The models describe diffusive transport of VOCs, biodegradation and source ageing. The mass transport is governed by Fick's law for diffusion. The equation for reactive transport of VOCs in the soil gas phase was solved for different source geometries and for different boundary conditions. Model results were compared to experimental data from a one-dimensional laboratory column and a radial-symmetric field experiment. The comparison yielded a satisfying agreement. The model results clearly illustrate the significant isotope fractionation by gas phase diffusion under transient state conditions. This leads to an initial depletion of heavy isotopes with increasing distance from the source. The isotope evolution of the source is governed by the combined effects of isotope fractionation due to vaporisation, diffusion and biodegradation. The net effect can lead to an enrichment or depletion of the heavy isotope in the remaining organic phase, depending on the compound and element considered. Finally, the isotope evolution of molecules migrating away from the source and undergoing degradation is governed by a combined degradation and diffusion isotope effect. This suggests that, in the unsaturated zone, the interpretation of biodegradation of VOC based on isotopic data must always be based on a model combining gas phase diffusion and degradation.

Variability of soil potential for biodegradation of petroleum hydrocarbons in a heterogeneous subsurface

Quantifying the spatial variability of factors affecting natural attenuation of hydrocarbons in the unsaturated zone is important to (i) performing a reliable risk assessment and (ii) evaluating the possibility for bioremediation of petroleum-polluted sites. Most studies to date have focused on the shallow unsaturated zone. Based on a data set comprising analysis of about 100 soil samples taken in a 16 m-deep unsaturated zone polluted with volatile petroleum compounds, we statistically and geostatistically analysed values of essential soil properties. The subsurface of the site was highly layered, resulting in an accumulation of pollution within coarse sandy lenses. Air-filled porosity, readily available phosphorous, and the first-order rate constant (k(1)) of benzene obtained from slurry biodegradation experiments were found to depend on geologic sample characterization (P<0.05), while inorganic nitrogen was homogenously distributed across the soil stratigraphy. Semivariogram analysis showed a spatial continuity of 4-8.6 m in the vertical direction, while it was 2-5 times greater in the horizontal direction. Values of k(1) displayed strong spatial autocorrelation. Even so, the soil potential for biodegradation was highly variable, which from autoregressive state-space modeling was partly explained by changes in soil air-filled porosity and gravimetric water content. The results suggest considering biological heterogeneity when evaluating the fate of contaminants in the subsurface.

Gas-phase diffusivity and tortuosity of structured soils

Modeling gas-phase diffusion of volatile contaminants in the unsaturated zone relies on soil-gas diffusivity models often developed for repacked and structureless soil columns. These suffer from the flaw of not reflecting preferential diffusion through voids and fractures in the soil, thus possibly causing an underestimation of vapor migration towards building foundations and vapor intrusion to indoor environments. We measured the ratio of the gas diffusion coefficient in soil and in free air (D(p)/D(0)) for 42 variously structured, intact, and unsaturated soil cores taken from 6 Danish sites. Whilst the results from structureless fine sand were adequately described using previously proposed models, results that were obtained from glacial clay till and limestone exhibited a dual-porosity behavior. Instead, these data were successfully described using a dual-porosity model for gas-phase diffusivity, considering a presence of drained fractures surrounded by a lower diffusivity matrix. Based on individual model fits, the tortuosity of fractures in till and limestone was found to be highest in samples with a total porosity <40%, suggesting soil compaction to affect the geometry of the fractures. In summary, this study highlights a potential order of magnitude underestimation associated in the use of classical models for prediction of subsurface gas-phase diffusion coefficients in heterogeneous and fractured soils.

Unsaturated zone leaching models for assessing risk to groundwater of contaminated sites

Risk assessments of sites contaminated with organic contaminants are typically conducted using models that ignore gas phase transport in the unsaturated zone. Here a general approach to developing analytical solutions to multiphase transport is presented. The approach is based on a combined gas and aqueous phase contaminant transport equation. The equation has the same general form as the standard advection-diffusion equation for which many analytical solutions have been derived. Four new analytical solutions are developed using this approach: a three-dimensional solution accounting for infiltration, lateral gas diffusion, sorption and degradation; a simple one-dimensional screening model, and two one-dimensional radial gas diffusion models for use in simulating volatile organic contaminant diffusion in unsaturated soils with an impermeable cover. The models show that both degradation and diffusion are important mechanisms for attenuation of contaminant concentrations at the water table. Finally, model results are compared with field data to illustrate the applicability of the solutions in risk assessment.

Carbon isotope fractionation during diffusion and biodegradation of petroleum hydrocarbons in the unsaturated zone: field experiment at Vaerløse Airbase, Denmark, and modeling

A field experiment was conducted in Denmark in order to evaluate the fate of 13 volatile organic compounds (VOCs) that were buried as an artificial fuel source in the unsaturated zone. Compound-specific isotope analysis showed distinct phases in the 13C/12C ratio evolution in VOC vapors within 3 m from the source over 114 days. At day 3 and to a lesser extent at day 6, the compounds were depleted in 13C by up to -5.7% per hundred with increasing distance from the source compared to the initial source values. This trend can be explained by faster outward diffusion of the molecules with 12C only compared to molecules with a 13C. Then, the isotope profile leveled out, and several compounds started to become enriched in 13C by up to 9.5% per hundred with increasing distance from the source, due to preferential removal of the molecules with 12C only, through biodegradation. Finally, as the amount of a compound diminished in the source, a 13C enrichment was also observed close to the source. The magnitude of isotope fractionation tended to be larger the smaller the mass of the molecule was. This study demonstrates that, in the unsaturated zone, carbon isotope ratios of hydrocarbons are affected by gas-phase diffusion in addition to biodegradation, which was confirmed using a numerical model. Gas-phase diffusion led to shifts in delta(13)C >1% per hundred during the initial days after the spill, and again during the final stages of source volatilization after >75% of a compound had been removed. In between, diffusion has less of an effect, and thus isotope data can be used as an indicator for hydrocarbon biodegradation.

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.