Relationship between mass-flux reduction and source-zone mass removal: analysis of field data

The magnitude of contaminant mass-flux reduction associated with a specific amount of contaminant mass removed is a key consideration for evaluating the effectiveness of a source-zone remediation effort. Thus, there is great interest in characterizing, estimating, and predicting relationships between mass-flux reduction and mass removal. Published data collected for several field studies were examined to evaluate relationships between mass-flux reduction and source-zone mass removal. The studies analyzed herein represent a variety of source-zone architectures, immiscible-liquid compositions, and implemented remediation technologies. There are two general approaches to characterizing the mass-flux-reduction/mass-removal relationship, end-point analysis and time-continuous analysis. End-point analysis, based on comparing masses and mass fluxes measured before and after a source-zone remediation effort, was conducted for 21 remediation projects. Mass removals were greater than 60% for all but three of the studies. Mass-flux reductions ranging from slightly less than to slightly greater than one-to-one were observed for the majority of the sites. However, these single-snapshot characterizations are limited in that the antecedent behavior is indeterminate. Time-continuous analysis, based on continuous monitoring of mass removal and mass flux, was performed for two sites, both for which data were obtained under water-flushing conditions. The reductions in mass flux were significantly different for the two sites (90% vs. approximately 8%) for similar mass removals ( approximately 40%). These results illustrate the dependence of the mass-flux-reduction/mass-removal relationship on source-zone architecture and associated mass-transfer processes. Minimal mass-flux reduction was observed for a system wherein mass removal was relatively efficient (ideal mass-transfer and displacement). Conversely, a significant degree of mass-flux reduction was observed for a site wherein mass removal was inefficient (non-ideal mass-transfer and displacement). The mass-flux-reduction/mass-removal relationship for the latter site exhibited a multi-step behavior, which cannot be predicted using some of the available simple estimation functions.     

Assessment of a simple function to evaluate the relationship between mass flux reduction and mass removal for organic-liquid contaminated source zones

The efficacy of a simple mass-removal function for characterizing mass-flux-reduction/mass-removal behavior for organic-liquid contaminated source zones was evaluated using the data obtained from a series of flow-cell experiments. The standard function, which employs a constant exponent, could not adequately reproduce the non-singular (multi-step) behavior exhibited by the measured data. Allowing the exponent to change as a function of mass removal (as the organic-liquid distribution and relative permeability change) produced non-singular relationships similar to those exhibited by the measured data. Four methods were developed to dynamically inform the exponent through use of measurable system-indicator parameters. Key factors that mediate the magnitude of mass flux (dilution and source accessibility) were accounted for using measures of source zone cross-sectional area, ganglia-to-pool (GTP) ratio, and relative permeability. The two methods that incorporated only the ganglia-to-pool ratio produced adequate simulations of the observed behavior for early stages of mass removal, but not for later stages. The method that incorporated parameters accounting for the source zone cross-sectional area (i.e., measure of system dilution) and source accessibility (GTP ratio and relative permeability) provided the most representative simulations of the observed data.     

Source-zone characterization of a chlorinated-solvent contaminated Superfund site in Tucson, AZ

An extensive site-characterization project was conducted at a large chlorinated-solvent contaminated Superfund site in Tucson, AZ. The project consisted of several components, including traditional site-characterization activities, tracer tests, laboratory experiments conducted with core material collected from the site, and mathematical modeling. The primary focus of the work presented herein is the analysis of induced-gradient contaminant elution tests conducted in a source zone at the site, investigation of the potential occurrence of immiscible liquid in the saturated zone, characterization of the relationship between mass flux reduction and mass removal, and evaluation of the impact of source-zone management on site remediation. The results of the present study, along with those of prior work, indicate that immiscible liquid is likely present in the saturated zone at the site source zones. Extensive tailing and rebound was observed for the contaminant-elution tests, indicating nonideal transport and mass-transfer behavior. The elution data were analyzed with a source-zone-scale mathematical model, and the results indicated that nonideal immiscible-liquid dissolution was the primary cause of the observed behavior. The time-continuous relationship between mass flux reduction and mass removal associated with the plume-scale pump-and-treat operation exhibited an initial large drop in mass flux with minimal mass removed, followed by a period of minimal mass flux reduction and a second period of large reduction. This behavior reflects the impact of both source-zone and aqueous-plume mass removal dynamics. Ultimately, a greater than 90% reduction in mass flux was achieved for a mass removal of approximately 50%. The influence of source-zone management on site remediation was evaluated by conducting two predictive simulations, one for which the source zones were controlled and one for which they were not. A plume-scale model was used to simulate the composite contaminant concentrations associated with groundwater extracted with the pump-and-treat system, which were compared to measured data. The information generated from this study was used to enhance the site conceptual model, help optimize operation of the pump-and-treat system, and evaluate the utility of source-zone remediation.     

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.     

Evaporative mass transfer behavior of a complex immiscible liquid

A series of laboratory experiments was conducted with a multiple-component immiscible liquid, collected from the Picillo Farm Superfund Site in Rhode Island, to examine liquid-vapor mass-transfer behavior. The immiscible liquid, which comprises solvents, oils, pesticides, PCBs, paint sludges, explosives, and other compounds, was characterized using gas chromatography and gas chromatography/mass spectrometry to determine mole fractions of selected constituents. Batch experiments were conducted to evaluate equilibrium phase-partitioning behavior. Two sets of air-stripping column studies were conducted to examine the mass-transfer dynamics of five selected target compounds present in the immiscible-liquid mixture. One set of column experiments was designed to represent a system with free-phase immiscible liquid present; the other was designed to represent a system with a residual phase of immiscible liquid. Initial elution behavior of all target components generally appeared to be ideal for both systems, as the initial vapor-phase concentrations were similar to vapor-phase concentrations measured for the batch experiment and those estimated using Raoult's law (incorporating the immiscible-liquid composition data). Later-stage removal of 1,2-dichlorobenzene appeared to be rate limited for the columns containing free-phase immiscible liquid and no porous medium. Conversely, evaporative mass transfer appeared to be ideal throughout the experiment conducted with immiscible liquid distributed relatively uniformly as a residual phase within a sandy porous medium.     

Design of aquifer remediation systems: (2) estimating site-specific performance and benefits of partial source removal

A Lagrangian stochastic model is proposed as a tool that can be utilized in forecasting remedial performance and estimating the benefits (in terms of flux and mass reduction) derived from a source zone remedial effort. The stochastic functional relationships that describe the hydraulic "structure" and non-aqueous phase liquid (NAPL) "architecture" have been described in a companion paper (Enfield, C.G., Wood, A.L., Espinoza, F.P., Brooks, M.C., Annable, M., Rao, P.S.C., this issue. Design of aquifer remediation systems: (1) describing hydraulic structure and NAPL architecture using tracers. J. Contam. Hydrol.). The previously defined functions were used along with the properties of the remedial fluids to describe remedial performance. There are two objectives for this paper. First, is to show that a simple analytic element model can be used to give a reasonable estimate of system performance. This is accomplished by comparing forecast performance to observed performance. The second objective is to display the model output in terms of change in mass flux and mass removal as a function of pore volumes of remedial fluid injected. The modelling results suggest that short term benefits are obtained and related to mass reduction at the sites where the model was tested.     

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.     

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.     

Assessing the impact of source-zone remediation efforts at the contaminant-plume scale through analysis of contaminant mass discharge

The long-term impact of source-zone remediation efforts was assessed for a large site contaminated by trichloroethene. The impact of the remediation efforts (soil vapor extraction and in-situ chemical oxidation) was assessed through analysis of plume-scale contaminant mass discharge, which was measured using a high-resolution data set obtained from 23 years of operation of a large pump-and-treat system. The initial contaminant mass discharge peaked at approximately 7kg/d, and then declined to approximately 2kg/d. This latter value was sustained for several years prior to the initiation of source-zone remediation efforts. The contaminant mass discharge in 2010, measured several years after completion of the two source-zone remediation actions, was approximately 0.2kg/d, which is ten times lower than the value prior to source-zone remediation. The time-continuous contaminant mass discharge data can be used to evaluate the impact of the source-zone remediation efforts on reducing the time required to operate the pump-and-treat system, and to estimate the cost savings associated with the decreased operational period. While significant reductions have been achieved, it is evident that the remediation efforts have not completely eliminated contaminant mass discharge and associated risk. Remaining contaminant mass contributing to the current mass discharge is hypothesized to comprise poorly accessible mass in the source zones, as well as aqueous (and sorbed) mass present in the extensive lower-permeability units located within and adjacent to the contaminant plume. The fate of these sources is an issue of critical import to the remediation of chlorinated-solvent contaminated sites, and development of methods to address these sources will be required to achieve successful long-term management of such sites and to ultimately transition them to closure.     

Operation of a two-phase partitioning bioreactor for the oxidation of anthracene by the enzyme manganese peroxidase

A study was conducted to determine the potential of a two-phase partitioning bioreactor (TPPB) for the treatment of a poorly soluble compound, anthracene, by the enzyme manganese peroxidase (MnP) from the fungus Bjerkandera sp. BOS55. Silicone oil was used as the immiscible solvent, which contained anthracene at high concentrations. The optimization of the oxidation process was conducted taking into account the factors which may directly affect the MnP catalytic cycle (the concentration of H(2)O(2) and malonic acid) and those that affect the mass transfer of anthracene between the organic and the aqueous phase (solvent and agitation speed). The main objective was carried out in terms of improved efficiency, i.e., maximizing the anthracene oxidized per unit of enzyme used. The TPPB reached nearly complete oxidation of anthracene at a conversion rate of 1.8mgl(-1)h(-1) in 56h, which suggests the application of enzymatic TPPBs for the removal of poorly soluble compounds.     

A stochastic-advective transport model for NAPL dissolution and degradation in non-uniform flows in porous media

Remediation schemes for contaminated sites are often evaluated to assess their potential for source zone reduction of mass, or treatment of the contaminant between the source and a control plane (CP) to achieve regulatory limits. In this study, we utilize a stochastic stream tube model to explain the behavior of breakthrough curves (BTCs) across a CP. At the local scale, mass dissolution at the source is combined with an advection model with first-order decay for the dissolved plume. Field-scale averaging is then employed to account for spatial variation in mass within the source zone, and variation in the velocity field. Under the assumption of instantaneous mass transfer from the source to the moving liquid, semi-analytical expressions for the BTC and temporal moments are developed, followed by derivation of expressions for effective velocity, dispersion, and degradation coefficients using the method of moments. It is found that degradation strongly influences the behavior of moments and the effective parameters. While increased heterogeneity in the velocity field results in increased dispersion, degradation causes the center of mass of the plume to shift to earlier times, and reduces the dispersion of the BTC by lowering the concentrations in the tail. Modified definitions of effective parameters are presented for degrading solutes to account for the normalization constant (zeroth moment) that keeps changing with time or distance to the CP. It is shown that anomalous dispersion can result for high degradation rates combined with wide variation in velocity fluctuations. Implications of model results on estimating cleanup times and fulfillment of regulatory limits are discussed. Relating mass removal at the source to flux reductions past a control plane is confounded by many factors. Increased heterogeneity in velocity fields causes mass fluxes past a control plane to persist, however, aggressive remediation between the source and CP can reduce these fluxes.     

A pragmatic approach for estimation of source-zone emissions at LNAPL contaminated sites

When considering natural attenuation as a remediation strategy at a site contaminated by a light non-aqueous phase liquid (LNAPL), it is important to consider the emission of contaminants from the source zone. A quantification of source-zone emissions is essential both for comparison with down-gradient mass fluxes to provide an estimate of fractional mass flux reduction, as well as for estimating the source lifetime. Because the spatial distribution of LNAPL at a field site is strongly dependent on both the spill circumstances and the heterogeneity of the geologic materials, which can be problematic for in-situ determination, alternative methods for estimating source-zone emissions are needed. In this work, a three-dimensional multiphase flow and transport modelling approach is used to investigate the relationship between the lateral extent of an LNAPL body and the emission of contaminants to groundwater at a contaminated site. For simulations involving an LNAPL release in an aquifer comprised of heterogeneous porosity and permeability distributions that were generated geostatistically, it is shown that a simple linear relationship exists between the lateral extent of the LNAPL body in the capillary fringe and the emission to the aqueous phase. The parameters describing the relationship are found to be linear functions of the groundwater flow velocity and the vertical infiltration rate. This site-specific relationship provides a simple method to estimate contaminant emissions to groundwater at LNAPL contaminated sites.     

U.S. Environmental Protection Agency’s Super-fund Innovative Technology Evaluation of Pneumatic Fracturing ExtractionSM

The U.S. Environmental Protection Agency (EPA), in cooperation with Accutech Remedial Systems (ARS) and the New Jersey Institute of Technology (NJIT), performed a field demonstration of Pneumatic Fracturing Extraction (PFE)^SM for the removal of chlorinated-volatile organic compounds (VOCs) from vadose zones of low permeability. The demonstration was conducted in the fall of 1992 at an industrial park in Somerville, New Jersey, where removal of VOC contamination in shale bedrock was required to comply with New Jersey’s Environmental Cleanup Responsibility Act (ECRA). During the demonstration, airflow and contaminant concentrations were monitored to establish a database against which the developer’s claims about the technology were evaluated. The developer contended that PFE would increase extracted airfkiw rates from the subsurface formation by at least 100 percent and would increase the mass removal rate for the key contaminant, trichloroethene (TCE), by at least 50 percent. Also, during the demonstration hot-gas injection was evaluated. Based on comparisons of four-hour test results before and after fracturing, airflow rates increased more than 600 percent, and TCE mass-removal rates increased about 675 percent. The increase in TCE mass-removal rates appeared to be a result, primarily, of the increased airflow. In addition, the extracted air contained significantly higher concentrations of other VOCs after fracturing. Using data developed in the four-hour postfracture test, the estimated cost for a hypothetical one-year clean-up is $140 per pound of TCE removed, or $140 per ton of soil contaminated with one pound of TCE. Experiments to evaluate the effects of injecting heated air, at 200 to 250°F, into the vadose zone gave inconclusive results.     

Nonideal transport of reactive contaminants in heterogeneous porous media: 7. distributed-domain model incorporating immiscible-liquid dissolution and rate-limited sorption/desorption

The purpose of this work is to present a distributed-domain mathematical model incorporating the primary mass-transfer processes that mediate the transport of immiscible organic liquid constituents in water-saturated, locally heterogeneous porous media. Specifically, the impact of grain/pore-scale heterogeneity on immiscible-liquid dissolution and sorption/desorption is represented in the model by describing the system as comprising a continuous distribution of mass-transfer domains. With this conceptualization, the distributions of the initial dissolution rate coefficient and the sorption/desorption rate coefficient are represented as probability density functions. Several sets of numerical experiments are conducted to examine the effects of heterogeneous dissolution and sorption/desorption on contaminant transport and elution. Four scenarios with different combinations of uniform/heterogeneous rate-limited dissolution and uniform/heterogeneous rate-limited sorption/desorption are evaluated. The results show that both heterogeneous rate-limited sorption/desorption and heterogeneous rate-limited dissolution can significantly increase the time or pore volumes required to elute immiscible-liquid constituents from a contaminated porous medium. However, sorption/desorption has minimal influence on elution behavior until essentially all of the immiscible liquid has been removed. For typical immiscible-liquid constituents that have relatively low sorption, the asymptotic elution tailing produced by heterogeneous rate-limited sorption/desorption begins at effluent concentrations that are several orders of magnitude below the initial steady-state concentrations associated with dissolution of the immiscible liquid. Conversely, the enhanced elution tailing associated with heterogeneous rate-limited dissolution begins at concentrations that are approximately one-tenth of the initial steady-state concentrations. Hence, dissolution may generally control elution behavior of immiscible-liquid constituents in cases wherein grain/pore-scale heterogeneity significantly influences both dissolution and sorption/desorption.     

Two-dimensional laboratory simulation of LNAPL infiltration and redistribution in the vadose zone

A quantitative two-dimensional laboratory experiment was conducted to investigate the immiscible flow of a light non-aqueous phase liquid (LNAPL) in the vadose zone. An image analysis technique was used to determine the two-dimensional saturation distribution of LNAPL, water and air during LNAPL infiltration and redistribution. Vertical water saturation variations were also continuously monitored with miniature resistivity probes. LNAPL and water pressures were measured using hydrophobic and hydrophilic tensiometers. This study is limited to homogeneous geological conditions, but the unique experimental methods developed will be used to examine more complex systems. The pressure measurements and the quantification of the saturation distribution of all the fluids in the entire flow domain under transient conditions provide quantitative data essential for testing the predictive capability of numerical models. The data are used to examine the adequacy of the constitutive pressure-saturation relations that are used in multiphase flow models. The results indicate that refinement of these commonly used hydraulic relations is needed for accurate model prediction. It is noted in particular that, in three-fluid phase systems, models should account for the existence of a residual NAPL saturation occurring after NAPL drainage. This is of notable importance because residual NAPL can act as a non negligible persistent source of contamination.     

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).     

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.     

液桥凝聚用于办公废纸脱墨的研究

研究了混合办公废纸脱墨的选择性凝聚技术,在废纸浆中添加分散在水中并不与水相溶的长链烷作凝聚剂,这些凝聚剂选择性地将印刷墨粉通过液桥连结在一起,且形成大颗粒的球形聚团,再分选脱除.研究表明只需要加入占纸重量4%的廉价长链烷,即可获得95%以上的脱墨率.因此在混合办公废纸脱墨时加入凝聚剂进行液桥凝聚,可降低脱墨浆中的尘埃度并增加产品的白度.     

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.