Modelling of iron cycling and its impact on the electron balance at a petroleum hydrocarbon contaminated site in Hnevice, Czech Republic

Over a period of several decades multiple leaks of large volumes from storage facilities located near Hnevice (Czech Republic) have caused the underlying Quaternary aquifer to be severely contaminated with nonaqueous phase liquid (NAPL) petroleum hydrocarbons. Beginning in the late 1980's the NAPL plume started to shrink as a consequence of NAPL dissolution exceeding replenishment and due to active remediation. The subsurface was classified geochemically into four different zones, (i) a contaminant-free zone never occupied by NAPL or dissolved contaminants, (ii) a re-oxidation zone formerly occupied by NAPL, (iii) a zone currently occupied by NAPL, and (iv) a lower fringe zone between the overlying NAPL and the deeper underlying contaminant-free zone. The study investigated the spatial and temporal variability of the redox zonation at the Hnevice site and quantified the influence of iron-cycling on the overall electron balance. As a first step inverse geochemical modelling was carried out to identify possible reaction models and mass transfer processes. In a subsequent step, two-dimensional (forward) multi-component reactive transport modelling was performed to evaluate and quantify the major processes that control the geochemical evolution at the site. The study explains the observed enrichment of the lower fringe zone with ferrihydrite as a result of the re-oxidation of ferrous iron. It suggests that once the NAPL zone started to shrink the dissolution of previously formed siderite and FeS by oxygen and nitrate consumed a significant part of the oxidation capacity for a considerable time period and therefore limited the penetration of electron acceptors into the NAPL contaminated zone.     

Effect of soil moisture dynamics on dense nonaqueous phase liquid (DNAPL) spill zone architecture in heterogeneous porous media

The amount, location, and form of NAPL in contaminated vadose zones are controlled by the spatial distribution of water saturation and soil permeability, the NAPL spill scenario, water infiltration events, and vapor transport. To evaluate the effects of these processes, we used the three-phase flow simulator STOMP, which includes a new permeability-liquid saturation-capillary pressure (k-S-P) constitutive model. This new constitutive model considers three NAPL forms: free, residual, and trapped. A 2-D vertical cross-section with five stratigraphic layers was assumed, and simulations were performed for seven cases. The conceptual model of the soil heterogeneity was based upon the stratigraphy at the Hanford carbon tetrachloride (CT) spill site. Some cases considered co-disposal of NAPL with large volumes of wastewater, as also occurred at the Hanford CT site. In these cases, the form and location of NAPL were most strongly influenced by high water discharge rates and NAPL evaporation to the atmosphere. In order to investigate the impact of heterogeneity, the hydraulic conductivity within the lower permeability layer was modeled as a realization of a random field having three different classes. For six extreme cases of 100 realizations, the CT mass that reached the water table varied by a factor of two, and was primarily controlled by the degree of lateral connectivity of the low conductivity class within the lowest permeability layer. The grid size at the top boundary had a dramatic impact on NAPL diffusive flux just after the spill event when the NAPL was present near the ground surface. NAPL evaporation with a fine grid spacing at the top boundary decreased CT mass that reached the water table by 74%, compared to the case with a coarse grid spacing, while barometric pumping had a marginal effect for the case of a continuous NAPL spill scenario considered in this work. For low water infiltration rate scenarios, the distribution of water content prior to a NAPL spill event decreased CT mass that reached the water table by 98% and had a significant impact on the formation of trapped NAPL. For all cases simulated, use of the new constitutive model that allows the formation of residual NAPL increased the amount of NAPL retained in the vadose zone. Density-driven advective gas flow from the ground surface controlled vapor migration in strongly anisotropic layers, causing NAPL mass flux to the lower layer to be reduced. These simulations indicate that consideration of the formation of residual and trapped NAPLs and dynamic boundary conditions (e.g., areas, rates, and periods of different NAPL and water discharge and fluctuations of atmospheric pressure) in the context of full three-phase flow are needed, especially for NAPL spill events at the ground surface. In addition, NAPL evaporation, density-driven gas advection, and NAPL vertical movement enhanced by water flow must be considered in order to predict NAPL distribution and migration in the vadose zone.     

A set of constitutive relationships accounting for residual NAPL in the unsaturated zone

Although laboratory experiments show that non-aqueous phase liquid (NAPL) is retained in the unsaturated zone, no existing multiphase flow model has been developed to account for residual NAPL after NAPL drainage in the unsaturated zone. We developed a static constitutive set of saturation-capillary pressure relationships for water, NAPL and air that accounts for both this residual NAPL and entrapped NAPL. The set of constitutive relationships is formulated similarly to the set of scaled relationships that is frequently applied in continuum models. The new set consists of three fluid-phase systems: a three-phase system and a two-phase system, that both comply with the original constitutive model, and a newly introduced residual NAPL system. The new system can be added relatively easily to the original two- and three-phase systems. Entrapment is included in the model. The constitutive relationships of the non-drainable residual NAPL system are based on qualitative fluid behavior derived from a pore scale model. The pore scale model reveals that the amount of residual NAPL depends on the spreading coefficient and the water saturation. Furthermore, residual NAPL is history-dependent. At the continuum scale, a critical NAPL pressure head defines the transition from free, mobile NAPL to residual NAPL. Although the Pc-S relationships for water and total liquid are not independent in case of residual NAPL, two two-phase Pc-S relations can represent a three-phase residual system of Pc-S relations. A newly introduced parameter, referred to as the residual oil pressure head, reflects the mutual dependency of water and oil. Example calculations show consistent behavior of the constitutive model. Entrapment and retention in the unsaturated zone cooperate to retain NAPL. Moreover, the results of our constitutive model are in agreement with experimental observations.     

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

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

Infiltration and redistribution of perchloroethylene in partially saturated, stratified porous media

Contamination of the subsurface by nonaqueous phase liquids (NAPLs) is a widespread problem. To investigate the behavior of a nonspreading, dense NAPL (DNAPL) in the vadose zone, we conducted perchloroethylene (PCE) infiltration experiments in nominally 1- and 2-dimensional (D), stratified porous media. In addition, the usefulness and limitations of a multifluid flow simulator to describe PCE infiltration and redistribution under the experimental conditions were tested. The physical simulations were conducted in a column (1-D) and a flow container (2-D) which were packed with two distinct layers of coarse-grained sand and a fine-grained sand layer in between. Volumetric water and PCE contents were determined with a fully automated dual-energy gamma radiation system. While migrating through the drier parts of the coarse-grained sand layers, PCE appeared to wet the water–air interface rather than displacing any water. In the wetter parts of the porous medium, PCE displaced water and behaved as a true nonwetting fluid. PCE showed a limited response to gradients in capillary pressure and rather high values for the volumetric PCE content were measured in the fine-grained sand layers. This was attributed to the nonspreading nature of PCE. The multifluid flow simulator appeared to predict the initial PCE movement in the vadose zone reasonably well. However, the model was not capable of predicting the final amounts of PCE retained in either the unsaturated or saturated part of the flow domain, mainly because the simulator does not consider the nonspreading flow behavior of NAPLs.     

Diffusive partitioning tracer test for the quantification of nonaqueous phase liquid (NAPL) in the vadose zone: Performance evaluation for heterogeneous NAPL distribution

A partitioning tracer test based on gas-phase diffusion in the vadose zone yields estimates of the residual nonaqueous phase liquid (NAPL) saturation. The present paper investigates this technique further by studying diffusive tracer breakthrough curves in the vadose zone for a heterogeneous NAPL distribution. Tracer experiments were performed in a lysimeter with a horizontal layer of artificial kerosene embedded in unsaturated sand. Tracer disappearance curves at the injection point and tracer breakthrough curves at some distance from the injection point were measured inside and outside of the NAPL layer. A numerical code was used to generate independent model predictions based on the physicochemical sand, NAPL, and tracer properties. The measured and modeled tracer breakthrough curves were in good agreement confirming the validity of important modeling assumptions such as negligible sorption of chlorofluorocarbon (CFC) tracers to the uncontaminated sand and their fast reversible partitioning between the soil air and the NAPL phase. Subsequently, the model was used to investigate different configurations of NAPL contamination. The experimental and model results show that the tracer disappearance curves of a single-well diffusive partitioning tracer test (DPTT) are dominated by the near-field presence of NAPL around the tip of the soil gas probe. In contrast, breakthrough curves of inter-well tracer tests reflect the NAPL saturation in between the probes, although there is no unique interpretation of the tracer signals if the NAPL distribution is heterogeneous. Numerical modeling is useful for the planning of a DPTT application. Simulations suggest that several cubic meters of soil can be investigated with a single inter-well partitioning tracer test of 24-hour duration by placing the injection point in the center of the investigated soil volume and probes at up to 1 m distance for the monitoring of gaseous tracers.     

Kinetic constraints on theIn-situ remediation of soils contaminated with organic chemicals

Cleanup of contaminated soils to comply with soil quality limits currently receives much interest.In-situ remediation of contaminated soils relies on the ability of the techniques employed to enhance the rate of release of contaminants from the soil-sorbed and nonaqueous phase liquid (NAPL) phases into the aqueous or gaseous phases from which they can be more readily removed and treated. Contaminant concentrations in these “environmentally mobile” forms usually decline over time so that the economic efficiency and the overall success of remediation technologies are subject to the “law of diminishing returns”. In this paper we consider the “state of the art” in our understanding of NAPL dissolution and transport, desorption of soilsorbed contaminants and fluid flow in porous media. The extent to which these processes may constrain the success of bioremediation, pump-and-treat remediation and soil venting in relation to established soil quality limits is addressed. Finally, we suggest directions for future research and comment on legislative considerations.     

Impact of nonaqueous phase liquid (NAPL) source zone architecture on mass removal mechanisms in strongly layered heterogeneous porous media during soil vapor extraction

An existing multiphase flow simulator was modified in order to determine the effects of four mechanisms on NAPL mass removal in a strongly layered heterogeneous vadose zone during soil vapor extraction (SVE): a) NAPL flow, b) diffusion and dispersion from low permeability zones, c) slow desorption from sediment grains, and d) rate-limited dissolution of trapped NAPL. The impacts of water and NAPL saturation distribution, NAPL-type (i.e., free, residual, or trapped) distribution, and spatial heterogeneity of the permeability field on these mechanisms were evaluated. Two different initial source zone architectures (one with and one without trapped NAPL) were considered and these architectures were used to evaluate seven different SVE scenarios. For all runs, slow diffusion from low permeability zones that gas flow bypassed was a dominant factor for diminished SVE effectiveness at later times. This effect was more significant at high water saturation due to the decrease of gas-phase relative permeability. Transverse dispersion contributed to fast NAPL mass removal from the low permeability layer in both source zone architectures, but longitudinal dispersion did not affect overall mass removal time. Both slow desorption from sediment grains and rate-limited mass transfer from trapped NAPL only marginally affected removal times. However, mass transfer from trapped NAPL did affect mass removal at later time, as well as the NAPL distribution. NAPL flow from low to high permeability zones contributed to faster mass removal from the low permeability layer, and this effect increased when water infiltration was eliminated. These simulations indicate that if trapped NAPL exists in heterogeneous porous media, mass transfer can be improved by delivering gas directly to zones with trapped NAPL and by lowering the water content, which increases the gas relative permeability and changes trapped NAPL to free NAPL.     

Anionic surfactant remediation of soil columns contaminated by nonaqueous phase liquids

A variety of column experiments have been completed for the purpose of selecting and evaluating suitable surfactants for remediation of nonaqueous phase liquids (NAPLs). The various NAPLs tested in the laboratory experiments were tetrachloroethylene (PCE), trichloroethylene (TCE), jet fuel (JP4) and a dense nonaqueous phase liquid from a site at Hill Air Force Base, UT. Both Ottawa sand and Hill field soil were used in these experiments. Surfactant candidates were first screened using phase behavior experiments and only the best ones were selected for the subsequent column experiments. Surfactants which showed high contaminant solubilization, fast coalescence times, and the absence of liquid crystal phases and gels during the phase behavior experiments were tested in soil column experiments. The primary objective of the soil column experiments was to identify surfactants that recovered at least 99% of the contaminant. The secondary objective was to identify surfactants that show low adsorption and little or no loss of hydraulic conductivity during the column experiments. Results demonstrated that up to 99.9% of the contaminants were removed as a result of surfactant flooding of the soil columns. The addition of xanthan gum polymer to the surfactant solution was shown to increase remediation efficiency as a lower volume of surfactant was required for recovering a given volume of NAPL. Based on these experimental results, guidelines for designing highly efficient and robust surfactant floods have been developed and applied to a field demonstration.     

Migration of contaminants through the unsaturated zone overlying the Hesbaye chalky aquifer in Belgium: a field investigation

This paper presents the results of a detailed field investigation that was performed for studying groundwater recharge processes and solute downward migration mechanisms prevailing in the unsaturated zone overlying a chalk aquifer in Belgium. Various laboratory measurements were performed on core samples collected during the drilling of boreholes in the experimental site. In the field, experiments consisted of well logging, infiltration tests in the unsaturated zone, pumping tests in the saturated zone and tracer tests in both the saturated and unsaturated zones. Results show that gravitational flows govern groundwater recharge and solute migration mechanisms in the unsaturated zone. In the variably saturated chalk, the migration and retardation of solutes is strongly influenced by recharge conditions. Under intense injection conditions, solutes migrate at high speed along the partially saturated fissures, downward to the saturated zone. At the same time, they are temporarily retarded in the almost immobile water located in the chalk matrix. Under normal recharge conditions, fissures are inactive and solutes migrate slowly through the chalk matrix. Results also show that concentration dynamics in the saturated zone are related to fluctuations of groundwater levels in the aquifer. A conceptual model is proposed to explain the hydrodispersive behaviour of the variably saturated chalk. Finally, the vulnerability of the chalk to contamination issues occurring at the land surface is discussed.     

SEAR技术修复土壤和地下水中NAPL污染的研究进展

就SEAR技术修复土壤及地下水中NAPL污染的原理及发展现状进行了综述.SEAR技术可以快速有效地去除土壤和地下水中的NAPL污染源,适于多种污染物.该技术通过增溶和增流2种途径提高NAPL污染物的去除率.表面活性剂的选择和微乳液体系的调配是SEAR技术实施的关键环节.将SEAR技术用于高浓度NAPL污染源的治理,并与生物修复和自然降解相结合,是经济高效的治理方案.     

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.     

Hot water flushing for immiscible displacement of a viscous NAPL

Thermal remediation techniques, such as hot water flooding, are emerging technologies that have been proposed for the removal of nonaqueous phase liquids (NAPLs) from the subsurface. In this study a combined laboratory and modeling investigation was conducted to determine if hot water flooding techniques would improve NAPL mass removal compared to ambient temperature water flushing. Two experiments were conducted in a bench scale two-dimensional sandbox (55 cmx45 cmx1.3 cm) and NAPL saturations were quantified using a light transmission apparatus. In these immiscible displacement experiments the aqueous phase, at 22 degrees C and 50 degrees C, displaced a zone with initial NAPL saturations on the order of 85%. The interfacial tension and viscosity of the selected light NAPL, Voltesso 35, are strongly temperature-dependent. Experimental results suggest that hot water flooding reduced the size of the high NAPL saturation zone, in comparison to the cold water flood, and yielded greater NAPL mass recovery (75% NAPL removal vs. 64%). Hot water flooding did not, however, result in lower residual NAPL saturations. A numerical simulator was modified to include simultaneous flow of water and organic phases, energy transport, temperature and pressure. Model predictions of mass removal and NAPL saturation profiles compared well with observed behavior. A sensitivity analysis indicates that the utility of hot water flooding improves with the increasing temperature dependence of NAPL hydraulic properties.     

Biosurfactant-enhanced solubilization of NAPL mixtures

Remediation of nonaqueous phase liquids (NAPLs) by conventional pump-and-treat methods (i.e., water flushing) is generally considered to be ineffective due to low water solubilities of NAPLs and to mass-transfer constraints. Chemical flushing techniques, such as surfactant flushing, can greatly improve NAPL remediation primarily by increasing the apparent solubility of NAPL contaminants. NAPLs at hazardous waste sites are often complex mixtures. However, the equilibrium and nonequilibrium mass-transfer characteristics between NAPL mixtures and aqueous surfactant solutions are not well understood. This research investigates the equilibrium solubilization behavior of two- and three-component NAPL mixtures (containing akylbenzenes) in biosurfactant solutions. NAPL solubilization is found to be ideal in water (i.e., obeys Raoult's Law), while solubilization in biosurfactant solutions was observed to be nonideal. Specifically, the relatively hydrophobic compounds in the mixture experienced solubility enhancements that were greater than those predicted by ideal enhanced solubilization theory, while the solubility enhancements for the relatively hydrophilic compounds were less than predicted. The degree of nonideality is shown to be a nonlinear function of the NAPL-phase mole fraction. Empirical relationships based on the NAPL-phase mole fraction and/or micelle-aqueous partition coefficients measured in single-component NAPL systems are developed to estimate values for the multicomponent partition coefficients. Empirical relationships that incorporate both the NAPL-phase mole fraction and single-component partition coefficients yield much improved estimates for the multicomponent partition coefficient.     

Plant aided bioremediation in the vadose zone: model development and applications

Phytoremediation has the potential to enhance clean up of land contaminated by various pollutants. A mathematical model that includes a two-fluid phase flow model of water flow as well as a two-region soil model of contaminant reactions was developed and applied to various bioremediation scenarios in the unsaturated zone, especially to plant-aided bioremediation. To investigate model behavior and determine the main parameters and mechanisms that affect bioremediation in unplanted and planted soils, numerical simulations of theoretical scenarios were conducted before applying the model to field data. It is observed from the results that parameters affecting the contaminant concentration in the water phase, such as aqueous solubility, the octanol-water partition coefficient, and organic carbon content of the soil controlled the contaminant fate in the vadose zone. Simulation using the developed model also characterized the fate and transport of the contaminants both in planted and unplanted soils satisfactorily for field applications. Although phytoremediation has the potential for remediation of contaminated soils, results from both modeling and field studies suggested that plants may not always enhance the remediation efficiency when the soil already has a high microbial concentration, when the contaminant bioavailability is low, or when the overall reaction is mass transfer-limited. Therefore, other steps to increase contaminant bioavailability are needed in phytoremediation applications; natural purification mechanisms such as aging, volatilization, and natural bioremediation should be considered to maximize the plant effect and minimize the cost.     

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.     

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.     

Simple screening models of NAPL dissolution in the subsurface

Three simple screening models of nonaqueous phase liquid (NAPL) dissolution in the subsurface are proposed based on the NAPL mass conservation and the assumption of proportionality between the residual NAPL source zone concentration and the remaining residual NAPL mass. The purpose of the proposed models is to predict the solute concentration in the zone of the residual NAPL as a result of dissolution. The predicted source zone concentration decrease is used to simulate and account for the decrease of dissolution rate with time. The proposed simple NAPL dissolution models enable the pseudo-equilibrium formulation to be used and therefore the numerical simulations for field application problems can be simplified compared to the non-equilibrium counterpart. With proper choice of empirical parameters, the proposed simple screening models can work as well as more complex dissolution rate correlation models, such as that of Imhoff et al. [Water Resour. Res. 30 (1994) 307-320]. It is found that the proposed models are very good for quantifying non-equilibrium dissolution, which is characterized by tailing of breakthrough curves. The models are especially useful for situations of small residual NAPL saturation, which are typical for many field applications.     

Non-aqueous-phase-liquid breakthrough during evaporative drying of clay barriers

In this study, an attempt has been made to model a real field scenario, whereby an initially almost saturated clay liner in a waste site is gradually drying, due to evaporation at its lower boundary. A detailed conceptual model that deals with the penetration and breakthrough of non-aqueous-phase-liquid (NAPL) in clay liners is introduced. Water content of clay samples was monitored during ambient evaporation through apertures at the base of sample holders. Clay drying rate served as the primary parameter for the NAPL breakthrough study. The interconnection between drying rates, structural damage formation (cracks and suction) and NAPL penetration is especially addressed. The processes taking place in the clay samples during drying appear to be associated with the capillary effects between the different fluid phases in the vicinity of either the NAPL-clay or the clay-air boundaries. A conceptual model of NAPL penetration and breakthrough of the clay layer has been considered, based on both indirect and direct observations of structural damages produced on either clay boundaries. A mutual interaction between these two boundaries is suggested and discussed. NAPL breakthrough is suggested to take place through cracks initiated on the upper soil surface.     

Three-dimensional visualization and quantification of non-aqueous phase liquid volumes in natural porous media using a medical X-ray Computed Tomography scanner

This study demonstrates the capabilities of a typical medical X-ray Computed Tomography (CT) scanner to non-destructively quantify non-aqueous phase liquid (NAPL) volumes, saturation levels, and three-dimensional spatial distributions in packed soil columns. Columns packed with homogeneous sand, heterogeneous sand, or natural soil, were saturated with water and injected with known quantities of gasoline or tetrachloroethene and scanned. A methodology based on image subtraction was implemented for computing soil porosity and NAPL volumes in each 0.35 mm x 0.35 mm x 1 mm voxel of the columns. Elimination of sample positioning errors and instrument drift artifacts was essential for obtaining reliable estimates of above parameters. The CT data-derived total NAPL volume was in agreement with the measured NAPL volumes injected into the columns. CT data-derived NAPL volume is subject to a 2.6% error for PCE and a 15.5% error for gasoline, at average NAPL saturations as low as 5%, and is mainly due to instrument noise. Non-uniform distributions of NAPL due to preferential flow, and accumulation of NAPL above finer-grained layers could be observed from the data on 3-D distributions of NAPL volume fractions.