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.     

土壤气相抽提去除土壤中汽油烃污染物柱试验研究

原位处理土壤石油污染对于土壤和地下水的有机污染控制具有极其重要的现实意义.通过砂土柱试验研究了原位物理通风的主要工艺形式及运行参数,并初步分析了汽油烃在砂土柱中的迁移和通风去除机制.结果表明,顶部真空抽提与底部注气两种通风方式相比,顶部真空抽提效果较好,砂土柱汽油烃初始质量浓度为2.937 mg/g时,经过104.5 h通风,砂土柱中汽油烃去除率达80.49%.土壤汽油烃初始浓度影响其在土壤中的迁移和去除,土壤汽油烃初始浓度越大,相同通风条件下,物理通风方法去除土壤中挥发性有机物的效率越低.通风及通风方式对砂土中的汽油烃的去除影响很大,连续通风可在砂土柱中形成稳定的负压环境,在汽油烃初始质量浓度为35.730 mg/g时,连续真空抽提264 h,砂土柱中的汽油烃平均去除率达89.29%;间歇通风在砂土柱中形成的负压环境不稳定,但也可以去除砂土柱中的汽油烃.初步分析认为,汽油烃存在负压作用下的向上挥发和重力作用下的向下迁移两个过程,其综合作用的结果导致汽油烃在砂土柱中的分布状况.     

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.     

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.     

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.     

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     

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.     

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.     

DNAPL source depletion: linking architecture and flux response

The relationship between dense non-aqueous phase liquid (DNAPL) mass reduction and contaminant mass flux was investigated experimentally in four model source zones. The flow cell design for the experiments featured a segmented extraction well that allowed for analysis of spatially resolved flux information. This flux information was coupled with image analysis of the NAPL spatial distribution to investigate the relationship between flux and the up-gradient NAPL architecture. Results indicate that in the systems studied, the relationship between DNAPL mass reduction and contaminant mass flux was primarily controlled by the NAPL architecture. A specific definition of NAPL architecture was employed where the source zone is resolved into a collection of streamtubes with spatial variability in NAPL saturation along each streamtube integrated and transformed into an effective NAPL content for each streamtube. The distribution of NAPL contents among the streamtubes (NAPL architecture) controlled dissolution dynamics. Two simplified models, a streamtube model and an effective Damkohler number model, were investigated for their ability to simulate dissolution dynamics.     

Effect of water content on transient nonequilibrium NAPL-gas mass transfer during soil vapor extraction

The effect of water content on the volatilization of nonaqueous phase liquid (NAPL) in unsaturated soils was characterized by one-dimensional venting experiments conducted to evaluate the lumped mass transfer coefficient. An empirical correlation based upon the modified Sherwood number, Peclet number, and normalized mean grain size was used to estimate initial lumped mass transfer coefficients over a range of water content. The effects of water content on the soil vapor extraction SVE process have been investigated through experimentation and mathematical modeling. The experimental results indicated that a rate-limited NAPL-gas mass transfer occurred in water-wet soils. A severe mass transfer limitation was observed at 61.0% water saturation where the normalized effluent gas concentrations fell below 1.0 almost immediately, declined exponentially from the initiation of venting, and showed long tailing. This result was attributed to the reduction of interfacial area between the NAPL and mobile gas phases due to the increased water content. A transient mathematical model describing the change of the lumped mass transfer coefficient was used. Simulations showed that the nonequilibrium mass transfer process could be characterized by the exponent beta, a parameter which described the reduction of the specific area available for NAPL volatilization. The nonequilibrium mass transfer limitations were controlled by the soil mean grain size and pore gas velocity, were well described by beta values below 1.0 at low water saturation, and were well predicted with beta values greater than 1.0 at high water saturation.     

Assessing performance and closure for soil vapor extraction: integrating vapor discharge and impact to groundwater quality

Soil vapor extraction (SVE) is typically effective for removal of volatile contaminants from higher-permeability portions of the vadose zone. However, contamination in lower-permeability zones can persist due to mass transfer processes that limit the removal effectiveness. After SVE has been operated for a period of time and the remaining contamination is primarily located in lower-permeability zones, the remedy performance needs to be evaluated to determine whether the SVE system should be optimized, terminated, or transitioned to another technology to replace or augment SVE. Numerical modeling of vapor-phase contaminant transport was used to investigate the correlation between measured vapor-phase mass discharge, MF(r), from a persistent, vadose-zone contaminant source and the resulting groundwater contaminant concentrations. This relationship was shown to be linear, and was used to directly assess SVE remediation progress over time and to determine the level of remediation in the vadose zone necessary to protect groundwater. Although site properties and source characteristics must be specified to establish a unique relation between MF(r) and the groundwater contaminant concentration, this correlation provides insight into SVE performance and support for decisions to optimize or terminate the SVE operation or to transition to another type of treatment.     

土壤柴油污染修复的抽气提取去除实验研究

为得到土壤气相抽提(SVE)去除柴油的优化条件,进行了一维土柱抽气提取去除柴油污染物的实验研究,研究不同初始含水率、不同抽气量对污染土壤中柴油去除率的影响及不同深度残留柴油的变化规律。结果表明：在本实验模拟的范围内,抽气量越大,SVE处理效果越好;初始含水率越低,处理效果越好;此外,不同深度去除率变化的规律基本上是随深度的增大而减小。实验结果可为土壤轻油污染实际治理提供实验数据基础。     

The partitioning and modelling of pesticide parathion in a surfactant-assisted soil-washing system

Soil sorption of organic pollutants has long been a problematic in the soil washing process because of its durability and low water solubility. This paper discussed the soil washing phenomena over a wide range of parathion concentrations and several soil samples at various fractions of organic content (foc) levels. When parathion dosage is set below the water solubility, washing performance is stable for surfactant concentrations above critical micelle concentration (cmc) and it is observed that more than 90% of parathion can be washed out when dosage is five times lower than the solubility limit. However, such trends change when non-aqueous phase liquids (NAPL) is present in the system. Parathion extraction depends very much on the surfactant dosage but is not affected by the levels of foc in the system. In between the extreme parathion dosage, a two-stage pattern is observed in these boundary regions. Washing performance is first increased with additional surfactant, but the increase slows down gradually since the sorption sites are believed to be saturated by the huge amount of surfactant in the system. A mathematical model has included foc to demonstrate such behavior and this can be used as a prediction for extraction.     

The effect of entrapped nonaqueous phase liquids on tracer transport in heterogeneous porous media: laboratory experiments at the intermediate scale

This work considers the applicability of conservative tracers for detecting high-saturation nonaqueous-phase liquid (NAPL) entrapment in heterogeneous systems. For this purpose, a series of experiments and simulations was performed using a two-dimensional heterogeneous system (10x1.2 m), which represents an intermediate scale between laboratory and field scales. Tracer tests performed prior to injecting the NAPL provide the baseline response of the heterogeneous porous medium. Two NAPL spill experiments were performed and the entrapped-NAPL saturation distribution measured in detail using a gamma-ray attenuation system. Tracer tests following each of the NAPL spills produced breakthrough curves (BTCs) reflecting the impact of entrapped NAPL on conservative transport. To evaluate significance, the impact of NAPL entrapment on the conservative-tracer breakthrough curves was compared to simulated breakthrough curve variability for different realizations of the heterogeneous distribution. Analysis of the results reveals that the NAPL entrapment has a significant impact on the temporal moments of conservative-tracer breakthrough curves.     

Experimental investigation of pneumatic soil vapor extraction

Soil Vapor Extraction (SVE) is a common remediation technique for removing volatile organic compounds from unsaturated contaminated soils. Soil heterogeneities can however cause serious limitations to the applicability of SVE due to air bypassing low permeable areas of the soil, leading to diffusion limitation of the remediation. To enhance removal from areas subject to diffusion limitation a new remediation technique, pneumatic soil vapor extraction, is proposed. In contrast to traditional SVE, in which soil vapor is extracted continuously by a vacuum pump, pneumatic SVE is based on enforcing a sequence of large pressure drops on the system to enhance the recovery from the low-permeable areas. The pneumatic SVE technique was investigated in the laboratory using TCE as a model contaminant. 2D-laboratory tank experiments were performed on homogeneous and heterogeneous sand packs. The heterogeneous packs consisted of a fine sand lens surrounded by a coarser sand matrix. As expected when using traditional SVE, the removal of TCE from the low permeable lens was extremely slow and subject to diffusion limitation. In contrast when pneumatic venting was used removal rates increased by up to 77%. The enhanced removal was hypothesized to be attributed to mixing of the contaminated air inside the lens and generation of net advective transport out of the lens due to air expansion.     

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.     

Multispectral image analysis method to determine dynamic fluid saturation distribution in two-dimensional three-fluid phase flow laboratory experiments

The need for measuring dynamic fluid saturation distribution in multi-dimensional three-fluid phase flow experiments is hampered by lack of appropriate techniques to monitor full field transient flow phenomena. There is no conventional technique able to measure dynamic three-fluid phase saturation at several array points of the flow field at the same time. A multispectral image analysis technique was developed to determine dynamic NAPL, water and air saturation distribution in two-dimensional three-fluid phase laboratory experiments. Using a digital near-infrared camera, images of sand samples with various degrees of NAPL, water and air saturation were taken, under constant lighting conditions and within three narrow spectral bands of the visible and near-infrared spectrum. It was shown that the optical density defined for the reflected luminous intensity was a linear function of the NAPL and the water saturation for each spectral band and for any two and three-fluid phase systems. This allowed the definition of dimensionless lump reflection coefficients for the NAPL and the water phase within each spectral band. Consequently, at any given time, two images taken within two different spectral bands provided two linear equations which could be solved for the water and the NAPL saturation. The method was applied to two-dimensional three-phase flow experiments, which were conducted to investigate the migration and the distribution of LNAPL in the vadose zone. The method was used to obtain continuous, quantitative and dynamic full field mapping of the NAPL saturation as well as the variation of the water and the air saturation during NAPL flow. The method provides a non-destructive and non-intrusive tool for studying multiphase flow for which rapid changes in fluid saturation in the entire flow domain is difficult to measure using conventional techniques.     

A micromodel analysis of factors influencing NAPL removal by surfactant foam flooding

A methodology to study the trichloroethylene (TCE) and dodecane removal in porous media by surfactant foams (SF) was presented by using etched-glass micromodels. The purpose of this work was to systematically evaluate the impact of various physicochemical factors such as gas fraction (GF), surfactant concentration, pore structure and nonaqueous phase liquid (NAPL) types on NAPL removal during SF flooding. The TCE displacement by SF was dependent on the gas fraction of SF. Low GFs (50% and 66%) were more efficient for TCE removal and sweep efficiencies than a high GF (85%). An increase in TCE removal was observed with increasing surfactant concentration at a fixed GF. TCE removal by SF flooding appeared to be dependent more to the value of Capillary number rather than to the concentration of surfactant solution. The effect of the pore heterogeneity was evaluated by employing two different types of micromodels. The Capillary number is an important parameter in the determination of sweep efficiency or gas saturation of SF in a nonhomogeneous porous medium. However, the TCE removal from a nonhomogeneous porous medium may not be associated with sweep efficiency. The initial configuration of residual TCE blobs in a nonhomogeneous porous medium would also be influential in displacing TCE. Sweep efficiencies and pressure responses of two NAPL systems (TCE and dodecane) were monitored to evaluate foam stability when the foam contacts the NAPLs. Stable foam contacting with TCE is implied, while it appears that dodecane cause the SF to collapse. All results indicate that the Capillary number (a ratio of viscous forces to capillary forces) is the most important parameter for TCE removal by SF flooding. Micromodel visualizations of water, surfactant and SF floods were showed and also discussed.     

有机污染土壤通风去污技术研究进展

土壤气体抽排净化技术是一种安全、高效、经济的有机污染土壤治理技术。本文系统分析了土壤气体抽排净化技术的特点、现状及发展趋势 ;总结了估算污染物净化时间的数学模型 ;概括了影响土壤气体抽排净化技术净化时间的主要因素 ,包括土壤透气率、含水率、有机污染物气体饱和蒸气压、抽排流速、环境温度 ;本文还对目前土壤气体抽排净化技术需要研究的问题及发展前景进行了展望。     

Analysis of a gas-phase partitioning tracer test conducted in an unsaturated fractured-clay formation

The gas-phase partitioning tracer method was used to estimate non-aqueous phase liquid (NAPL), water, and air saturations in the vadose zone at a chlorinated-solvent contaminated field site in Tucson, AZ. The tracer test was conducted in a fractured-clay system that is the confining layer for the underlying regional aquifer. Three suites of three tracers were injected into wells located 14, 24, and 24 m from a single, central extraction well. The tracers comprised noble gases (traditionally thought to be nonsorbing), alkanes (primarily water partitioning), perfluorides (primarily NAPL partitioning), and halons (both NAPL and water partitioning). Observations of vacuum response were consistent with flow in a fractured system. The halon tracers exhibited the greatest amount of retardation, and helium and the perfluoride tracers the least. The alkane tracers were unexpectedly more retarded than the perfluoride tracers, indicating low NAPL saturations and high water saturations. An NAPL saturation of 0.01, water saturation of 0.215, and gas saturation of 0.775 was estimated based on analysis of the suite of tracers comprising helium, perfluoromethylcyclohexane and dibromodifluoromethane, which was considered to be the most robust set. The estimated saturations compare reasonably well to independently determined values.