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

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

Locally-calibrated light transmission visualization methods to quantify nonaqueous phase liquid mass in porous media

Five locally-calibrated light transmission visualization (LTV) methods were tested to quantify nonaqueous phase liquid (NAPL) mass and mass reduction in porous media. Tetrachloroethylene (PCE) was released into a two-dimensional laboratory flow chamber packed with water-saturated sand which was then flushed with a surfactant solution (2% Tween 80) until all of the PCE had been dissolved. In all the LTV methods employed here, the water phase was dyed, rather than the more common approach of dyeing the NAPL phase, such that the light adsorption characteristics of NAPL did not change as dissolution progressed. Also, none of the methods used here required the use of external calibration chambers. The five visualization approaches evaluated included three methods developed from previously published models, a binary method, and a novel multiple wavelength method that has the advantage of not requiring any assumptions about the intra-pore interface structure between the various phases (sand/water/NAPL). The new multiple wavelength method is also expected to be applicable to any translucent porous media containing two immiscible fluids (e.g., water-air, water-NAPL). Results from the sand-water-PCE system evaluated here showed that the model that assumes wetting media of uniform pore size (Model C of Niemet and Selker, 2001) and the multiple wavelength model with no interface structure assumptions were able to accurately quantify PCE mass reduction during surfactant flushing. The average mass recoveries from these two imaging methods were greater than 95% for domain-average NAPL saturations of approximately 2.6x10(-2), and were approximately 90% during seven cycles of surfactant flushing that sequentially reduced the average NAPL saturation to 7.5x10(-4).     

Influence of cation type, ionic strength, and pH on solubilization and mobilization of residual hydrocarbon by a biosurfactant

This study investigated the effect of cation type, ionic strength, and pH on the performance of an anionic monorhamnolipid biosurfactant for solubilization and removal of residual hexadecane from sand. Three common soil cations, Na⁺, Mg²⁺, and Ca²⁺, were used in these experiments and hexadecane was chosen to represent a nonaqueous phase liquid (NAPL) less dense than water. Results showed that hexadecane solubility in rhamnolipid solution was significantly increased by the addition of Na⁺ and Mg²⁺. Addition of up to 0.2 mM Ca²⁺ also increased hexadecane solubility. For Ca²⁺ concentrations greater than 0.2 mM there was little effect on hexadecane solubility due to competing effects of calcium-induced rhamnolipid precipitation and enhanced hexadecane solubilization. Efficiency of NAPL solubilization can be expressed in terms of molar solubilization ratios (MSR). The results showed that MSR values for hexadecane in rhamnolipid solutions increased 7.5-fold in the presence of 500 mM Na⁺, and 25-fold in the presence of 1 mM Mg²⁺. The presence of cations also reduced the interfacial tension between rhamnolipid solutions and hexadecane. For example, an increase in Na⁺ from 0 to 800 mM caused a decrease in interfacial tension from 2.2 to 0.89 dyn cm⁻¹. Similarly, decreasing pH caused a reduction in interfacial tension. The lowest interfacial tension value observed in this study was 0.02 dyn cm⁻¹ at pH 6 in the presence of 320 mM Na⁺. These conditions were also found to be optimal for removal of hexadecane residual from sand columns, with 58% of residual removed within three pore volumes. The removal of residual NAPL from the packed columns was primarily by mobilization, even though solubilization was significantly increased in the presence of Na⁺.     

Investigation of surfactant-enhanced mass removal and flux reduction in 3D correlated permeability fields using magnetic resonance imaging

Magnetic resonance imaging (MRI) was used to visualize the NAPL source zone architecture before and after surfactant-enhanced NAPL dissolution in three-dimensional (3D) heterogeneously packed flowcells characterized by different longitudinal correlation lengths: 2.1 cm (aquifer 1) and 1.1 cm (aquifer 2). Surfactant flowpaths were determined by imaging the breakthrough of a paramagnetic tracer (MnCl(2)) analyzed by the method of moments. In both experimental aquifers, preferential flow occurred in high permeability materials with low NAPL saturations, and NAPL was preferentially removed from the top of the aquifers with low saturation. Alternate flushing with water and two surfactant pulses (5-6 pore volumes each) resulted in approximately 63% of NAPL mass removal from both aquifers. However, overall reduction in mass flux (Mass Flux 1) exiting the flowcell was lower in aquifer 2 (68%) than in aquifer 1 (81%), and local effluent concentrations were found to increase by as high as 120 times at local sampling ports from aquifer 2 after surfactant flushing. 3D MRI images of NAPL revealed that NAPL migrated downward and created additional NAPL source zones in previously uncontaminated areas at the bottom of the aquifers. The additional NAPL source zones were created in the direction transverse to flow in aquifer 2, which explains the higher mass flux relative to aquifer 1. Analysis using a total trapping number indicates that mobilization of NAPL trapped in the two coarsest sand fractions is possible when saturation is below 0.5 and 0.4, respectively. Results from this study highlight the potential impacts of porous media heterogeneity and NAPL source zone architecture on advanced in-situ flushing technologies.     

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.     

Effect of soil amendments on sorption and mobility of metribuzin in soils

Metribuzin (4-amino-6-tert-butyl-4,5-dihydro-3-methylthio-1,2,4-triazin-5-one), is weakly sorbed to soil therefore, leaches easily to lower soil profiles. Soil amendments play a significant role in the management of leaching losses of pesticides. Therefore, present study reports the effect of organic manure and fly ash amendments on metribuzin downward mobility in sandy loam soil columns. Application of animal manure [T-1(OM) and T-2(OM)] and fly ash [T-1(FA) and T-2(FA)] at 2.5% and 5.0% levels increased the metribuzin retention in the soil. Freundlich constant [K(f)(1/n)] values of metribuzin for treatments T-1(OM) and T-2(OM) were 0.70 and 1.11, respectively, which were significantly higher than the value (0.27) in natural soil (T-0). The respective values for treatments T-1(FA) and T-2(FA) were 1.80 and 4.61. Downward mobility of metribuzin was studied in packed soil columns [300 mm (l)x59 mm (i.d.)]. Both the amendments significantly reduced the downward mobility of metribuzin and affected breakthrough time and maximum concentration of metribuzin in the leachate. Leaching losses of metribuzin were decreased from 97% in natural soil (T-0) column to 64% [T-1(OM)] and 42% [T-2(OM)] for animal manure-amended columns and 26% [T-1(FA)] to 100% [T-2(FA)] for fly ash-amended columns, as metribuzin did not leach out of 5% fly ash-amended column. Study indicates that both animal manure and fly ash were quite effective in reducing the downward mobility of metribuzin in packed soil columns of a sandy loam soil.     

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.     

Impact of edible oil injection on the permeability of aquifer sands

Recent laboratory and field studies have shown that food-grade edible oils can be injected into the subsurface for installation of in-situ permeable reactive barriers. However to be effective, the oil must be distributed out away from the oil injection points without excessive permeability loss. In this work, we examine the distribution of soybean oil in representative aquifer sediments as non-aqueous phase liquid oil (NAPL oil) or as an oil-in-water emulsion. Laboratory columns packed with sands or clayey sands were flushed with either NAPL oil or a soybean emulsion followed by plain water, while monitoring permeability loss and the final oil residual saturation. NAPL oil can be injected into coarse-grained sands. However NAPL injection into finer grained sediments requires high injection pressures which may not be feasible at some sites. In addition, NAPL injection results in high oil residual saturations and moderate permeability losses. In contrast, properly prepared emulsions can be distributed through sands with varying clay content without excessive pressure buildup, low oil retention and very low to moderate permeability loss. For effective transport, the emulsion must be stable, the oil droplets must be significantly smaller than the mean pore size of the sediment and the oil droplets should have a low to moderate tendency to stick to each other and the aquifer sediments. In our work, oil retention and associated permeability loss increased with sediment clay content and with the ratio of droplet size to pore size. For sandy sediments, the permeability loss is modest (0-40% loss) and is proportional to the oil residual saturation.     

Surfactant enhanced recovery of tetrachloroethylene from a porous medium containing low permeability lenses. 2. Numerical simulation

A numerical model of surfactant enhanced solubilization was developed and applied to the simulation of nonaqueous phase liquid recovery in two-dimensional heterogeneous laboratory sand tank systems. Model parameters were derived from independent, small-scale, batch and column experiments. These parameters included viscosity, density, solubilization capacity, surfactant sorption, interfacial tension, permeability, capillary retention functions, and interphase mass transfer correlations. Model predictive capability was assessed for the evaluation of the micellar solubilization of tetrachloroethylene (PCE) in the two-dimensional systems. Predicted effluent concentrations and mass recovery agreed reasonably well with measured values. Accurate prediction of enhanced solubilization behavior in the sand tanks was found to require the incorporation of pore-scale, system-dependent, interphase mass transfer limitations, including an explicit representation of specific interfacial contact area. Predicted effluent concentrations and mass recovery were also found to depend strongly upon the initial NAPL entrapment configuration. Numerical results collectively indicate that enhanced solubilization processes in heterogeneous, laboratory sand tank systems can be successfully simulated using independently measured soil parameters and column-measured mass transfer coefficients, provided that permeability and NAPL distributions are accurately known. This implies that the accuracy of model predictions at the field scale will be constrained by our ability to quantify soil heterogeneity and NAPL distribution.     

Magnetic resonance imaging of nonaqueous phase liquid during soil vapor extraction in heterogeneous porous media

Soil vapor extraction (SVE) is commonly used to remediate nonaqueous phase liquids (NAPLs) from the vadose zone. This paper aims to determine the effect of grain size heterogeneity on the removal of NAPL in porous media during SVE. Magnetic resonance imaging (MRI) was used to observe and quantify the amount and location of NAPL in flow-through columns filled with silica gel grains. MRI is unique because it is nondestructive, allowing three-dimensional images to be taken of the phases as a function of space and time. Columns were packed with silica gel in three ways: coarse grains (250-550 microm) only, fine grains (32-63 microm) only, and a core of fine grains surrounded by a shell of coarse grains. Columns saturated with water were drained under a constant suction head, contaminated with decane, and then drained to different decane saturations. Each column was then continuously purged with water-saturated nitrogen gas and images were taken intermittently. Results showed that at residual saturation, a sharp volatilization front moved through the columns filled with either coarse-grain or fine-grain silica gel. In the heterogeneous columns, the volatilization front in the core lagged just behind the shell because gas flow was greater through the shell and decane in the core diffused outward to the shell. When decane saturation in the core was above residual saturation, decane volatilization occurred near the inlet, the relative decane saturation throughout the core dropped uniformly, and decane in the core flowed in the liquid phase to the shell to replenish volatilized decane. These results indicate that NAPL trapped in low-permeability zones can flow to replenish areas where NAPL is lost due to SVE. However, when residual NAPL saturation is reached, NAPL flow no longer occurs and diffusion limits removal from low-permeability zones.     

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.     

Performance assessment of NAPL remediation in heterogeneous alluvium

Over the last few years, more than 40 partitioning interwell tracer tests (PITTs) have been conducted at many different sites to measure nonaqueous phase liquid (NAPL) saturations in the subsurface. While the main goal of these PITTs was to estimate the NAPL volume in the subsurface, some were specifically conducted to assess the performance of remedial actions involving NAPL removal. In this paper, we present a quantitative approach to assess the performance of remedial actions to recover NAPL that can be used to assess any NAPL removal technology. It combines the use of PITTs (to estimate the NAPL volume in the swept pore volume between injection and extraction wells of a test area) with the use of several cores to determine the vertical NAPL distribution in the subsurface. We illustrate the effectiveness of such an approach by assessing the performance of a surfactant/foam flood conducted at Hill Air Force Base, UT, to remove a TCE-rich NAPL from alluvium with permeability contrasts as high as one order of magnitude. In addition, we compare the NAPL volumes determined by the PITTs with volumes estimated through geostatistical interpolation of aquifer sediment core data collected with a vertical frequency of 5-10 cm and a lateral borehole spacing of 0.15 m. We demonstrate the use of several innovations including the explicit estimation of not only the errors associated with NAPL volumes and saturations derived from PITTs but also the heterogeneity of the aquifer sediments based upon permeability estimates. Most importantly, we demonstrate the reliability of the     

Remediation of NAPL below the water table by steam-induced heat conduction

Previous experimental studies have shown that NAPL will be removed when it is contacted by steam. However, in full-scale operations, steam may not contact the NAPL directly and this is the situation addressed in this study. A two-dimensional intermediate scale sand box experiment was performed where an organic contaminant was emplaced below the water table at the interface between a coarse and a fine sand layer. Steam was injected above the water table and after an initial heating period the contaminant was recovered at the outlet. The experiment was successfully modeled using the numerical code T2VOC and the dominant removal mechanism was identified to be heat conduction induced boiling of the separate phase contaminant. Subsequent numerical modeling showed that this mechanism was insensitive to the porous medium properties and that it could be evaluated by considering only one-dimensional heat conduction.     

Fingerprinting of Gasoline and Coal Tar NAPL Volatile Hydrocarbons Dissolved in Groundwater

Accidental spills and chronic leaks of fuel oil or other hydrocarbon material (e.g., coal tar) often result in subsurface accumulation of nonaqueous phase liquid (NAPL), which can be a subsequent source of contamination in groundwater. Linking hydrocarbons in groundwater to a source NAPL has been difficult when using standard target analytes (e.g., BTEX) because of differences in partitioning properties of the analytes between the source NAPL and groundwater. Because aqueous solubility is predicted to be the controlling influence in the partitioning of hydrocarbons from NAPL to groundwater, a solubility-based approach to matching dissolved hydrocarbons in groundwater to their source NAPL has been developed and validated for two sites with commonly encountered types of NAPL contamination. Specifically, a gasoline LNAPL and a coal tar DNAPL from two separate sites (West Virginia and California) and groundwater interfaced with these NAPLs were analyzed for approximately 50 gasoline-range hydrocarbons consisting of paraffin, isoparaffin, (mono-) aromatic, naphthene, and olefin compounds (PIANO). Solubility characteristics of selected alkyl aromatic hydrocarbons from the PIANO analysis were used to identify a set of diagnostic hydrocarbons, expressed as hydrocarbon ratios, which were found to be useful in distinguishing the source(s) of hydrocarbons in groundwater. At the West Virginia site, the diagnostic ratios in a downgradient groundwater sample were similar to those of a gasoline NAPL at that site, indicating the source of hydrocarbons to the groundwater was the upgradient gasoline NAPL. The diagnostic ratios of the groundwater in contact with the gasoline NAPL and the remote groundwater were also similar, providing evidence that the diagnostic ratios were retained during transport in the aquifer. At the California site, diagnostic ratios in a cross-gradient groundwater sample differed from those of the coal tar NAPL at that site, indicating that the remote groundwater hydrocarbons did not originate from the coal tar contamination. Environmental factors such as selective degradation of specific isomers and various geological conditions (e.g., soil mineralogy, and organic content) may confound the application of this solubility-based fingerprinting approach. Thus, it is recommended that multiple diagnostic pairs be simultaneously evaluated when considering this fingerprinting approach for specific sites and product types.     

Movement of metolachlor and terbuthylazine in core and packed soil columns

Movement of metolachlor and terbuthylazine including a bromide tracer was studied in core and packed soil columns in PVC pipes (80 mm diameter, 15 mm depth) with two German soil types viz: silt loam and loamy silt. The breakthrough curves (BTCs) for bromide indicated some preferential flow of water both under conventional tillage (CN) and no-till (NT) simulation with silt loam soil. The herbicides leached to a greater extent in NT columns than in CN columns. Leaching was higher in loamy silt soil than in silt loam soil under CN conditions. This result is in agreement with the higher sorption capacity of silt loam having higher organic carbon compared to loamy silt having low organic carbon. Adsorption strength of the herbicides did not affect their breakthrough time, but was reflected in the slope and maximum height of the BTCs. The BTCs showed the expected inverse relationship between leaching and adsorption with greater mobility of the weakly-sorbed metolachlor than the more strongly sorbed terbuthylazine. Maximum amounts of the applied herbicides were recovered from the top soil layer in intact columns. Metolachlor was more mobile in packed columns than in core columns.     

Flow behavior and residual saturation formation of liquid carbon tetrachloride in unsaturated heterogeneous porous media

The formation of residual, discontinuous nonaqueous phase liquids (NAPLs) in the vadose zone is a process that is not well understood. To obtain data that can be used to study the development of a residual NAPL saturation in the vadose zone and to test current corresponding models, detailed transient experiments were conducted in intermediate-scale columns and flow cell. The column experiments were conducted to determine residual carbon tetrachloride (CCl(4)) saturations of two sands and to evaluate the effect of CCl(4) vapors on the water distribution. In the intermediate-scale flow cell experiment, a rectangular zone of the fine-grained sand was packed in an otherwise medium-grained matrix. A limited amount of CCl(4) was injected from a small source and allowed to redistribute until a pseudo steady state situation had developed. A dual-energy gamma radiation system was used to determine fluid saturations at numerous locations. The experiments clearly demonstrated the formation of residual CCl(4) saturations in both sands. Simulations with an established multifluid flow simulator show the shortcomings of current relative permeability-saturation-capillary pressure (k-S-P) models. The results indicate that nonspreading behavior of NAPLs should be implemented in simulators to account for the formation of residual saturations.     

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.     

NAPL migration and ecotoxicity of conventional and renewable fuels in accidental spill scenarios

Fuels derived from non-petroleum renewable resources have raised interest due to their potential in replacing petroleum-based fuels, but information on their fate and effects in the terrestrial and aquatic environments in accidental spill scenario is limited. In this study, migration of four fuels (conventional diesel, conventional gasoline, renewable diesel NExBTL, and ethanol-blended gasoline RE85 containing maximum 85 % ethanol) as non-aqueous phase liquids (NAPL) in soil was demonstrated in a laboratory-scale experiment. Ecotoxicity data was produced for the same fuels. There was no significant difference in migration of conventional and renewable diesel, but gasoline migrated 1.5 times deeper and 7–9 times faster in sand than diesel. RE85 spread horizontally wider but not as deep (p?Eisenia fetida followed by ethanol-blended gasoline (LC50 1,643 mg/kg THC) and conventional diesel (LC50 2,432 mg/kg THC), although gasoline evaporated fast from soil. For comparison, the toxicity of the water-accommodated fractions (WAF) of the fuels was tested with water flea Daphnia magna and Vibrio fischeri, also demonstrating groundwater toxicity. The WAF of conventional gasoline and RE85 showed almost similar toxicity to both the aquatic test species. EC50 values of 1:10 (by volume) WAF were 9.9 %WAF (gasoline) and 9.3 %WAF (RE85) to D. magna and 9.3 %WAF (gasoline) and 12.3 %WAF (RE85) to V. fischeri. Low solubility decreased toxicity potential of conventional diesel in aquatic environment, but direct physical effects of oil phase pose a threat to organisms in nature. Renewable diesel NExBTL did not show clear toxicity to any test species.     

Evaluation of simplified mass transfer models to simulate the impacts of source zone architecture on nonaqueous phase liquid dissolution in heterogeneous porous media

Nonaqueous phase liquid (NAPL) dissolution was studied in three-dimensional (3D) heterogeneous experimental aquifers (25.5 cm x 9 cm x 8.5 cm) with two different longitudinal correlation lengths (2.1 cm and 1.1 cm) and initial spill volumes (22.5 ml and 10.5 ml). Spatial and temporal distributions of NAPL during dissolution were measured using magnetic resonance imaging (MRI). At high NAPL spill volume, average effluent concentrations initially increased during dissolution, as NAPL pools transitioned to NAPL ganglia, and then decreased as the total NAPL-water interfacial area decreased over time. Experimental results were used to test six dissolution models: (i and ii) a one-dimensional (1D) model using either specific NAPL-water interfacial area values estimated from MR images at each time step (i.e., 1D quasi-steady state model), or an empirical mass transfer (Sh') correlation (i.e., 1D transient model), (iii and iv) a multiple analytical source superposition technique (MASST) using either the NAPL distribution determined from MR images at each time step (i.e., MASST steady state model), or the NAPL distribution determined from mass balance calculations (i.e., MASST transient model), (v) an equilibrium streamtube model, and (vi) a 3D grid-scale pool dissolution model (PDM) with a dispersive mass flux term. The 1D quasi-steady state model and 3D PDM captured effluent concentration values most closely, including some concentration fluctuations due to changes in the extent of flow reduction. The 1D transient, MASST steady state and transient, and streamtube models all showed a monotonic decrease in effluent concentration values over time, and the streamtube model was the most computationally efficient. Changes during dissolution of the effective NAPL-water interfacial area estimated from imaging data are similar to changes in effluent concentration values. The 1D steady state model incorporates estimates of the effective NAPL-water interfacial area directly at each time point; the 3D PDM does so indirectly through mass balance and a relative permeability function, which causes reduced water flow through high saturation NAPL regions. Hence, when model accuracy is required, the results indicate that a surrogate of this effective interfacial area is required. Approaches to include this surrogate in the MASST and streamtube models are recommended.     

Design of aquifer remediation systems: (1) describing hydraulic structure and NAPL architecture using tracers

Aquifer heterogeneity (structure) and NAPL distribution (architecture) are described based on tracer data. An inverse modelling approach that estimates the hydraulic structure and NAPL architecture based on a Lagrangian stochastic model where the hydraulic structure is described by one or more populations of lognormally distributed travel times and the NAPL architecture is selected from eight possible assumed distributions. Optimization of the model parameters for each tested realization is based on the minimization of the sum of the square residuals between the log of measured tracer data and model predictions for the same temporal observation. For a given NAPL architecture the error is reduced with each added population. Model selection was based on a fitness which penalized models for increasing complexity. The technique is demonstrated under a range of hydrologic and contaminant settings using data from three small field-scale tracer tests: the first implementation at an LNAPL site using a line-drive flow pattern, the second at a DNAPL site with an inverted five-spot flow pattern, and the third at the same DNAPL site using a vertical circulation flow pattern. The Lagrangian model was capable of accurately duplicating experimentally derived tracer breakthrough curves, with a correlation coefficient of 0.97 or better. Furthermore, the model estimate of the NAPL volume is similar to the estimates based on moment analysis of field data.