Temporal changes in kerosene content and composition in field soil as a result of leaching

A field experiment was designed to determine the combined effect of leaching and natural attenuation on the redistribution dynamics of kerosene--a volatile petroleum hydrocarbon mixture (VPHM)--and of its selected individual components in the soil subsurface. Variables included the composition of contaminant spilled, the soil water content before contamination and the leaching pattern. Temporal changes in the residual kerosene concentration and composition in the soil subsurface of the experimental field during 39 days and leaching by 500 mm of irrigation water were determined to a depth of 100 cm. The main processes controlling contaminant attenuation were volatilization and redistribution with depth. Soil hydration status was found to affect the attenuation, redistribution and composition of VPHM in the porous media. An initial relative increase of n-alkanes in the subsurface compared with the total VPHM in the first leaching period was a result of the volatilization of low vapor pressure compounds. The redistribution of individual components in the soil profile during leaching was in accordance with their physico-chemical properties.     

非饱和土壤轻油污染多相流研究进展

轻油泄漏到地下后 ,在非饱和土壤中水、气、油三相相互作用、相互影响 ,共同构成多相流体系统。在简要分析轻油污染问题的产生、危害以及轻油泄漏在地表以下运移方式的基础上 ,介绍并总结了近些年来国内外土壤轻油污染有关的多相流研究方面的进展 ,包括国外的模型研究成果、国内外试验研究方法和研究成果 ,并指出目前研究中存在的问题和需要进一步探讨的方向     

Oxygenated gasoline release in the unsaturated zone - Part 1: Source zone behavior

Oxygenates present in gasoline, such as ethanol and MTBE, are a concern in subsurface contamination related to accidental spills. While gasoline hydrocarbon compounds have low solubility, MTBE and ethanol are more soluble, ethanol being completely miscible with water. Consequently, their fate in the subsurface is likely to differ from that of gasoline. To evaluate the fate of gasoline containing oxygenates following a release in the unsaturated zone shielded from rainfall/recharge, a controlled field test was performed at Canadian Forces Base Borden, in Ontario. 200L of a mixture composed of gasoline with 10% ethanol and 4.5% MTBE was released in the unsaturated zone, into a trench 20cm deep, about 32cm above the water table. Based on soil cores, most of the ethanol was retained in the source, above the capillary fringe, and remained there for more than 100 days. Ethanol partitioned from the gasoline to the unsaturated pore-water and was retained, despite the thin unsaturated zone at the site (~35cm from the top of the capillary fringe to ground surface). Due to its lower solubility, most of the MTBE remained within the NAPL as it infiltrated deeper into the unsaturated zone and accumulated with the gasoline on top of the depressed capillary fringe. Only minor changes in the distribution of ethanol were noted following oscillations in the water table. Two methods to estimate the capacity of the unsaturated zone to retain ethanol are explored. It is clear that conceptual models for sites impacted by ethanol-fuels must consider the unsaturated zone.     

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.     

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.     

地下水中轻质有机污染物(LNAPL)透镜体研究

在二维砂槽模型中模拟了轻质油在均匀多孔介质地下水非饱和区中的运移过程。模拟结果表明,地下水毛细区是轻质油污染的重点区,除了ＬＮＡＰＬ的残留以外,进入地下水饱和的ＬＮＡＰＬ终将被地下水顶托回到毛细区中,毛细区以上的约大多民将在重力作用下进入毛细区中,试验中观察到达稳定状态时ＬＮＡＰＬ透镜体的上边缘略微高出毛细区。利用多孔介质毛细管模型,建立了利用界面张力、接触角、介质特征孔隙直径等物理量估算不同位置     

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.     

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.     

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.     

A sensitivity analysis of adsorption and degradation parameters in the modeling of pesticide transport in soils

The increasing use of deterministic models in predicting the movement of pesticides in soils, has focused attention on the evaluation of major parameters which represent attenuation factors of organics in the subsurface. These parameters are the degradation rate constant and the adsorption constant for the pesticide. In view of the large in situ variability of these parameters and of the difficulty in obtaining accurate field data, there is a high degree of uncertainty associated with the results obtained from deterministic models. A sensitivity analysis is performed here to quantify the impact of such variation in each of these input parameters on the output results of an unsaturated zone transport model (PRZM). Results show that variations in these parameters about their respective mean values greatly affect the predicted concentration distributions, obtained after three years, of the pesticide aldicarb in all the soil profile. A 15–22% variation in the degradation constant, or a 24% variation in the adsorption constant, lead to a 100% uncertainty in the various simulation results defined as the cumulative quantity of aldicarb or the dissolved aldicarb concentration leached below the root zone (or the unsaturated zone) of the soil. Such a deterministic model presents a high degree of sensitivity to these input parameters. Accurate field data are then needed to obtain reliable model results in predicting pesticide movement inthe unsaturated zone.     

Numerical methods for describing chemical transport in the unsaturated zone of the subsurface

Finite-difference and finite-element methods of approximation have been extended to solve the one-dimensional nonlinear partial differential equations that describe the simultaneous transport of heat, moisture and chemical in the unsaturated zone. Especially for chemical transport, nodal spacing criteria are required to minimize numerical dispersion and oscillatory behavior in the solution vector for chemical concentration. Conservative criteria for nodal spacing for saturated flow can be used to set nodal spacing for unsaturated zone transport. When nodal spacing criteria are satisfied, for the same set of transport and boundary conditions, chemical concentration profiles calculated by the two numerical methods will be almost the same. A situation that is simulated very well with one-dimensional models, is the application of chemicals to land surfaces. To compare and contrast the characteristics of solutions given by the two numerical methods, moisture content, temperature and chemical concentration profiles for a 75-day period after application in the unsaturated zone are calculated for two representative types of organic chemicals. In the first, the chemical is very slowly degraded in the subsurface environment but strongly sorbed to soil surfaces. In the second, the chemical is rapidly degraded but weakly sorbed to soil surfaces. Because of differences in sorption coefficients and mechanisms of degradation, for the same set of hydrodynamic properties of the subsurface, the weakly sorbed chemical is more widely distributed throughout the unsaturated zone, whereas the strongly sorbed chemical stays very close to where it is put initially with little penetration into the subsurface. Satisfying nodal spacing criteria minimizes the impact of the method of approximation on the calculated solutions of the transport equations. For better model predictive performance, however, there are needs for more fundamental information on processes governing transport in the subsurface.     

Modelling assessment of regional groundwater contamination due to historic smelter emissions of heavy metals

Historic emissions from ore smelters typically cause regional soil contamination. We developed a modelling approach to assess the impact of such contamination on groundwater and surface water load, coupling unsaturated zone leaching modelling with 3D groundwater transport modelling. Both historic and predictive modelling were performed, using a mass balance approach for three different catchments in the vicinity of three smelters. The catchments differ in their hydrology and geochemistry. The historic modelling results indicate that leaching to groundwater is spatially very heterogeneous due to variation in soil characteristics, in particular soil pH. In the saturated zone, cadmium is becoming strongly retarded due to strong sorption at neutral pH, even though the reactivity of the sandy sediments is low. A comparison between two datasets (from 1990 to 2002) on shallow groundwater and modelled concentrations provided a useful verification on the level of statistics of "homogeneous areas" (areas with comparable land use, soil type and geohydrological situation) instead of comparison at individual locations. While at individual locations observations and the model varies up to two orders of magnitude, for homogeneous areas, medians and ranges of measured concentrations and the model results are similar. A sensitivity analysis on metal input loads, groundwater composition and sediment geochemistry reveals that the best available information scenario based on the median value of input parameters for the model predicts the range in observed concentrations very well. However, the model results are sensitive to the sediment contents of the reactive components (organic matter, clay minerals and iron oxides). Uncertainty in metal input loads and groundwater chemistry are of lesser importance. Predictive modelling reveals a remarkable difference in geochemical and hydrological controls on subsurface metal transport at catchment-scale. Whether the surface water load will peak within a few decades or continue to increase until after 2050 depends on the dominant land use functions in the areas, their hydrology and geochemical build-up.     

Non-aqueous phase liquid spreading during soil vapor extraction

Many non-aqueous phase liquids (NAPLs) are expected to spread at the air-water interface, particularly under non-equilibrium conditions. In the vadose zone, this spreading should increase the surface area for mass transfer and the efficiency of volatile NAPL recovery by soil vapor extraction (SVE). Observations of spreading on water wet surfaces led to a conceptual model of oil spreading vertically above a NAPL pool in the vadose zone. Analysis of this model predicts that spreading can enhance the SVE contaminant recovery compared to conditions where the liquid does not spread. Experiments were conducted with spreading volatile oils hexane and heptane in wet porous media and capillary tubes, where spreading was observed at the scale of centimeters. Within porous medium columns up to a meter in height containing stagnant gas, spreading was less than ten centimeters and did not contribute significantly to hexane volatilization. Water film thinning and oil film pinning may have prevented significant oil film spreading, and thus did not enhance SVE at the scale of a meter. The experiments performed indicate that volatile oil spreading at the field scale is unlikely to contribute significantly to the efficiency of SVE.     

Effects of degree of water saturation on dispersivity and immobile water in sandy soil columns

Three natural nonaggregated soil samples, with similar grain-size distributions, have been used to determine the dispersive behavior of porous media under steady, saturated and unsaturated flow conditions. Tritium was used as a tracer and was found to have no sorption on the solid matrix. Generated breakthrough curves (BTCs) for the unsaturated experiments were symmetrical with no evidence of tailing. The unsaturated experiments for two of the soils were adequately described by considering all the water in the pore volume as mobile. However, about 10% of the pore water, independent of the degree of saturation, was found to be immobile in the case of the third soil during unsaturated flow. For this soil, there was no mass transfer between the two water regions, indicating that the immobile water is essentially isolated from the flowing water fraction. For all three soils, dispersivity under unsaturated conditions was found to be higher, independent of the degree of water saturation, than the value determined for the saturated experiments. This is inconsistent with what would be expected from the simple bundle-of-capillary-tubes model and does not agree well with a more sophisticated conceptualization of the porous medium. The data, however, clearly indicate a wider range in pore-water velocities when these soils are desaturated.     

Metal ion and other key water quality constituent attenuation in soil columns

A set of soil columns was constructed to simulate discharge of disinfected tertiary treated wastewater to a river via nearby land application or indirect discharge. The system was primarily designed to observe the fate of metal ions and nutrients. The following three experiments were conducted: (1) flow through saturated soils only, which simulates indirect discharge where water is directly applied to groundwater; (2) flow through unsaturated soil followed by saturated flow, which simulates vadose then saturated zone transport; and (3) saturated flow only using ethylene diamine tetraacetic acid-metal chelates, which determined effects of metal organic complexes on metal mobility through the soil. Metal ion attenuation was substantial but not complete in experiments 1 and 2 (removal: 68% Cu2+, 43% Ni2+, 98% Pb2+, and 96% Hg2+), which was somewhat contrary to modeling results. Cyanide attenuation was also monitored (92% removal). In experiment 3, lead attenuation was somewhat reduced (92% removal) and delayed (requiring additional residence time); copper attenuation was significantly reduced (38% removal) and delayed; and nickel concentrations were higher in the 28-day sample (> 80 microg/L) than in the column feed water (58 microg/L). Near-complete denitrification and total phosphorus attenuation were observed. For the water quality constituents studied, unsaturated (vadose zone) transport did not appear to add additional benefit.     

Factors affecting the distribution of hydrocarbon contaminants and hydrogeochemical parameters in a shallow sand aquifer

The distributions of hydrocarbon contaminants and hydrogeochemical parameters were investigated in a shallow sand aquifer highly contaminated with petroleum hydrocarbons leaked from solvent storage tanks. For these purposes, a variety of field investigations and studies were performed, which included installation of over 100 groundwater monitoring wells and piezometers at various depths, soil logging and analyses during well and piezometer installation, chemical analysis of groundwater, pump tests, and slug tests. Continuous water level monitoring at three selected wells using automatic data-logger and manual measuring at other wells were also conducted. Based on analyses of the various investigations and tests, a number of factors were identified to explain the distribution of the hydrocarbon contaminants and hydrogeochemical parameters. These factors include indigenous biodegradation, hydrostratigraphy, preliminary pump-and-treat remedy, recharge by rainfall, and subsequent water level fluctuation. The permeable sandy layer, in which the mean water table elevation is maintained, provided a dominant pathway for contaminant transport. The preliminary pump-and-treat action accelerated the movement of the hydrocarbon contaminants and affected the redox evolution pattern. Seasonal recharge by rain, together with indigenous biodegradation, played an important role in the natural attenuation of the petroleum hydrocarbons via mixing/dilution and biodegradation. The water level fluctuations redistributed the hydrocarbon contaminants by partitioning them into the soil and groundwater. The identified factors are not independent but closely inter-correlated.     

Characterization of flow and transport processes within the unsaturated zone of Yucca Mountain, Nevada, under current and future climates

This paper presents a large-scale modeling study characterizing fluid flow and tracer transport in the unsaturated zone of Yucca Mountain, Nevada, a potential repository site for storing high-level radioactive waste. The study has been conducted using a three-dimensional numerical model, which incorporates a wide variety of field data and takes into account the coupled processes of flow and transport in the highly heterogeneous, unsaturated fractured porous rock. The modeling approach is based on a dual-continuum formulation of coupled multiphase fluid and tracer transport through fractured porous rock. Various scenarios of current and future climate conditions and their effects on the unsaturated zone are evaluated to aid in the assessment of the proposed repository's system performance using different conceptual models. These models are calibrated against field-measured data. Model-predicted flow and transport processes under current and future climates are discussed.     

Estimation of ecotoxicity of petroleum hydrocarbon mixtures in soil based on HPLC–GCXGC analysis

Detailed HPLC–GCXGC/FID (high performance liquid chromatography followed by comprehensive two-dimensional gas chromatography with flame-ionization detection) analysis of oil-contaminated soils was performed to interpret results of selected acute ecotoxicity assays. For the five ecotoxicity assays tested, plant seed germination and Microtox^® were selected as most sensitive for evaluating ecotoxicity of the oil in the soil phase and in the leaching water, respectively. The measured toxicity for cress when testing the soil samples did not correspond to TPH concentration in the soil. A detailed chemical composition analysis of the oil contamination using HPLC–GCXGC/FID allows to better predict the ecotoxicological risk and leaching potential of petroleum hydrocarbons in soil. Cress biomass production per plant was well correlated to the total aromatic hydrocarbon concentration (R² = 0.79, n = 6), while cress seed germination was correlated (R² = 0.82, n = 6) with total concentration of “highly water-soluble aromatic hydrocarbons” (HSaromatics). The observed ecotoxicity of the leaching water for Microtox-bacteria related well to calculated (based on the HPLC–GCXGC/FID results) petroleum hydrocarbon equilibrium concentrations in water.     

Effect of biochar on the fate of volatile petroleum hydrocarbons in an aerobic sandy soil

Biochar addition to soil is currently being investigated as a novel technology to remediate polluted sites. A critical consideration is the impact of biochar on the intrinsic microbial pollutant degradation, in particular at sites polluted with a mixture of readily biodegradable and more persistent organic pollutants. We therefore studied the impact of biochar (2% on dry weight basis) on the fate of volatile petroleum hydrocarbons in an aerobic sandy soil with batch and column studies. The soil-water partitioning coefficient, K(d), was enhanced in the biochar-amended soil up to a factor 36, and petroleum hydrocarbon vapor migration was retarded accordingly. Despite increased sorption, in particular of monoaromatic hydrocarbons, the overall microbial respiration was comparable in the biochar-amended and unamended soil. This was due to more rapid biodegradation of linear, cyclic and branched alkanes in the biochar amended soil. We concluded that the total petroleum hydrocarbon degradation rate was controlled by a factor other than substrate availability and the reduced availability of monoaromatic hydrocarbons in the biochar amended soil led to greater biodegradation of the other petroleum compounds.