The effect of spontaneous gas expansion and mobilization on the aqueous-phase concentrations above a dense non-aqueous phase liquid pool

The spontaneous expansion and mobilization of discontinuous gas above dense non-aqueous-phase liquid (DNAPL) pools can affect the aqueous-phase concentrations of the DNAPL constituents above the pool. The results of an intermediate-scale, two-dimensional flow cell experiment showed that the discontinuous gas flow produced by spontaneous expansion, driven by the partitioning of 1,1,1-TCA from the surface of a DNAPL pool, resulted in detectable aqueous-phase concentrations of 1,1,1-TCA well above the pool surface. In comparison to a conventional model for DNAPL pool dissolution in the absence of a discontinuous gas phase, these concentrations were greater than expected, and were present at greater than expected elevations. Additionally, this study showed that the discontinuous gas flow produced transient behavior in the aqueous-phase concentrations, where the elevated concentrations occurred as short-term, pulse-like events. These results suggest that the spontaneous expansion and mobilization of discontinuous gas in DNAPL source zones could lead to the misdiagnosis of source zone architecture using aqueous concentration data, and that the transient nature of the elevated concentrations could further complicate the difficult task of source zone characterization.     

Comparison of alternative upscaled model formulations for simulating DNAPL source dissolution and biodecay

This paper compares the performance of analytical and numerical approaches for modeling DNAPL dissolution with biodecay. A solution derived from a 1-D advective transport formulation (“Parker” model) is shown to agree very closely with high resolution numerical solutions. A simple lumped source mass balance solution in which with decay is assumed proportional to DNAPL mass (“Falta1” model) over- or underpredicts aqueous phase biodecay depending on the magnitude of the exponential factor governing the relationship between dissolution rate and DNAPL mass. A modification of the Falta model that assumes decay proportional to the source exit concentration is capable of accurately simulating source behavior with strong aqueous phase biodecay if model parameters are appropriately selected or calibrated (“Falta2” model). However, parameters in the lumped models exhibit complex interdependencies that cannot be quantified without consideration of transport processes within the source zone. Combining the Falta2 solution with relationships derived from the Parker model was found to resolve these limitations and track the numerical model results. A method is presented to generalize the analytical solutions to enable simulation of partial mass removal with changes in source parameters over time due to various remedial actions. The algorithm is verified by comparison with numerical simulation results. An example application is presented that demonstrates the interactions of partial mass removal, enhanced biodecay, enhanced mass transfer and source zone flow reduction applied at various time periods on contaminant flux reduction. Increasing errors that arise in numerical solutions with coarse discretization and high decay rates are shown to be controlled by using an adjusted decay coefficient derived from the Parker analytical solution.     

Experimental Study of the Effects of DNAPL Distribution on Mass Rebound

The release of stored dissolved contaminants from low permeability zones contributes to plume persistence beyond the time when dense nonaqueous phase liquid (DNAPL) has completely dissolved. This is fundamental to successfully meeting acceptable low concentrations in groundwater that are driven by site‐specific cleanup goals. The study goals were to assess the role of DNAPL entrapment morphology on mass storage and plume longevity. As controlled field studies are not feasible, two‐dimensional (2D) test tanks were used to quantify the significance of mass loading processes from source dissolution and stored mass rebound. A simple two‐layer soil domain representing a high permeable formation sand overlying a zone of lower permeability sand was used in the tests. DNAPL mass depletion through dissolution was monitored via X‐ray photon attenuation, and effluent samples were used to monitor the plume. These data enabled analysis of the DNAPL distribution, the dissolved plume, and the dissolved phase distribution within the low permeability layer. Tests in an intermediate tank showed that mass storage contributes substantially to plume longevity. Detectable effluent concentrations persisted long after DNAPL depletion. The small tank results indicated that the DNAPL morphology influenced the flow field and caused distinctive transport mechanisms contributing to mass storage. Zones of high DNAPL saturation at the interface between the low and high permeability layers exhibited flow bypassing and diffusion dominated transport into the low permeability layer. In the absence of a highly saturated DNAPL zone near the soil interface the contaminant penetrated deeper into the low permeability layer caused by a combination of advection and diffusion.     

Mass Flux Analysis of Abiotic Tetrachloroethene Dense Nonaqueous Phase Liquid Pool Dissolution in a Heterogeneous Flow Environment

Porous aquifer materials are often characterized by layered heterogeneities that influence groundwater flow and present complexities in contaminant transport modeling. Such flow variations also have the potential to impact the dissolution flux from dense nonaqueous phase liquid (DNAPL) pools. This study examined how these heterogeneous flow conditions affected the dissolution of a tetrachloroethene (PCE) pool in a two-dimensional intermediate-scale flow cell containing coarse sand. A steady-state mass-balance approach was used to calculate the PCE dissolution rate at three different flow rates. As expected, aqueous PCE concentrations increased along the length of the PCE pool and higher flow rates decreased the aqueous PCE concentration in the effluent. Nonreactive tracer studies at two flow rates confirmed the presence of a vertical flow gradient, with the most rapid velocity located at the bottom of the tank. These results suggest that flow focusing occurred near the DNAPL pool. Effluent PCE concentrations and pool dissolution flux rates were compared to model predictions assuming local equilibrium (LE) conditions at the DNAPL pool/aqueous phase interface and a uniform distribution of flow. The LE model did not describe the data well, even over a wide range of PCE solubility and macroscopic transverse dispersivity values. Model predictions assuming nonequilibrium mass-transfer-limited conditions and accounting for vertical flow gradients, however, resulted in a better fit to the data. These results have important implications for evaluating DNAPL pool dissolution in the field where subsurface heterogeneities are likely to be present.     

Transport and fate of nonaqueous phase liquid (NAPL) in variably saturated porous media with evolving scales of heterogeneity

 A stochastic simulation is performed to study multiphase flow and contaminant transport in fractal porous media with evolving scales of heterogeneity. Numerical simulations of residual NAPL mass transfer and subsequent transport of dissolved and/or volatilized NAPL mass in variably saturated media are carried out in conjunction with Monte Carlo techniques. The impact of fractal dimension, plume scale and anisotropy (stratification) of fractal media on relative dispersivities is investigated and discussed. The results indicate the significance of evolving scale of porous media heterogeneity to the NAPL transport in the subsurface. In general, the fractal porous media enhance the dispersivities of NAPL mass plume transport in both the water phase and the gas phase while the influence on the water phase is more significant. The porous media with larger fractal dimension have larger relative dispersivities. The aqueous horizontal dispersivity exhibits a most significant increase against the plume scale.     

Hydraulic displacement of dense nonaqueous phase liquids for source zone stabilization

Hydraulic displacement is a mass removal technology suitable for stabilization of a dense, nonaqueous phase liquid (DNAPL) source zone, where stabilization is defined as reducing DNAPL saturations and reducing the risk of future pool mobilization. High resolution three-dimensional multiphase flow simulations incorporating a spatially correlated, heterogeneous porous medium illustrate that hydraulic displacement results in an increase in the amount of residual DNAPL present, which in turn results in increased solute concentrations in groundwater, an increase in the rate of DNAPL dissolution, and an increase in the solute mass flux. A higher percentage of DNAPL recovery is associated with higher initial DNAPL release volumes, lower density DNAPLs, more heterogeneous porous media, and increased drawdown of groundwater at extraction wells. The fact that higher rates of recovery are associated with more heterogeneous porous media stems from the fact that larger contrasts in permeability provide for a higher proportion of capillary barriers upon which DNAPL pooling and lateral migration can occur. Across all scenarios evaluated in this study, the ganglia-to-pool (GTP) ratio generally increased from approximately 0.1 to between approximately 0.3 and 0.7 depending on the type of DNAPL, the degree of heterogeneity, and the imposed hydraulic gradient. The volume of DNAPL recovered as a result of implementing hydraulic displacement ranged from between 9.4% and 45.2% of the initial release volume, with the largest percentage recovery associated with 1,1,1 trichloroethane, the least dense of the three DNAPLs considered.     

Bioaugmentation in a Well‐Characterized Fractured Rock DNAPL Source Area

A field demonstration was performed at Edwards Air Force Base to assess bioaugmentation for treatment of a well‐characterized tetrachloroethene (PCE) dense nonaqueous phase liquid (DNAPL) source area in fractured rock. Groundwater recirculation was employed to deliver remedial amendments, including bacteria, to facilitate reductive dechlorination and enhance DNAPL dissolution. An active treatment period of 9 months was followed by a 10‐month posttreatment rebound evaluation. Dechlorination daughter products were observed in both the shallow and deep fracture zones following treatment. In the shallow fracture zone, the calculated DNAPL mass removed was approximately equal to the DNAPL mass estimated using partitioning tracer testing, and no rebound in chlorinated ethenes or ethene was observed during the posttreatment period. A maximum DNAPL dissolution enhancement factor of 5 was observed in the shallow fracture zone. In the deep fracture zone, only approximately 45% of the DNAPL mass—as estimated via partitioning tracer testing—was removed and rebound in the total molar chlorinated ethenes + ethene was observed. The difference in behavior between the shallow and deep fracture zones was attributed to DNAPL architecture and the fracture flow field.     

222Rn as Natural Tracer for LNAPL Recovery in a Crude Oil‐Contaminated Aquifer

The objective of this study was to investigate whether ²²²Rn in groundwater can be used as a tracer for light non‐aqueous phase liquid (LNAPL) quantification at a field site treated by dual‐phase LNAPL removal. After the break of a pipeline, 5 ha of soil in the nature reserve Coussouls de Crau in southern France was contaminated by 5100 m³ of crude oil. Part of this oil seeped into the underlying gravel aquifer and formed a floating oil body of about 3.9 ha. The remediation consists of plume management by hydraulic groundwater barriers and LNAPL extraction in the source zone. ²²²Rn measurements were performed in 21 wells in and outside the source zone during 15 months. In uncontaminated groundwater, the radon activity was relatively constant and remained always >11 Bq/L. The variability of radon activity measurements in wells affected by the pump‐and‐skim system was consistent with the measurements in wells that were not impacted by the system. The mean activities in wells in the source zone were, in general, significantly lower than in wells upgradient of the source zone, owing to partitioning of ²²²Rn into the oil phase. The lowest activities were found in zones with high non‐aqueous phase liquid (NAPL) recovery. LNAPL saturations around each recovery well were furthermore calculated during a period of high groundwater level, using a laboratory‐determined crude oil–water partitioning coefficient of 38.5 ± 2.9. This yielded an estimated volume of residual crude oil of 309 ± 93 m³ below the capillary fringe. We find that ²²²Rn is a useful and cheap groundwater tracer for finding zones of good LNAPL recovery in an aquifer treated by dual‐phase LNAPL removal, but that quantification of NAPL saturation using Rn is highly uncertain.     

Transition probability/Markov chain analyses of DNAPL source zones and plumes

At sites where a dense nonaqueous phase liquid (DNAPL) was spilled or released into the subsurface, estimates of the mass of DNAPL contained in the subsurface from core or monitoring well data, either in the nonaqueous or aqueous phase, can be highly uncertain because of the erratic distribution of the DNAPL due to geologic heterogeneity. In this paper, a multiphase compositional model is applied to simulate, in detail, the DNAPL saturations and aqueous-phase plume migration in a highly characterized, heterogeneous glaciofluvial aquifer, the permeability and porosity data of which were collected by researchers at the University of Tübingen, Germany. The DNAPL saturation distribution and the aqueous-phase contaminant mole fractions are then reconstructed by sampling the data from the forward simulation results using two alternate approaches, each with different degrees of sampling conditioning. To reconstruct the DNAPL source zone architecture, the aqueous-phase plume configuration, and the contaminant mass in each phase, one method employs the novel transition probability/Markov chain approach (TP/MC), while the other involves a traditional variogram analysis of the sampled data followed by ordinary kriging. The TP/MC method is typically used for facies and/or hydraulic conductivity reconstruction, but here we explore the applicability of the TP/MC method for the reconstruction of DNAPL source zones and aqueous-phase plumes. The reconstructed geometry of the DNAPL source zone, the dissolved contaminant plume, and the estimated mass in each phase are compared using the two different geostatistical modeling approaches and for various degrees of data sampling from the results of the forward simulation. It is demonstrated that the TP/MC modeling technique is robust and accurate and is a preferable alternative compared to ordinary kriging for the reconstruction of DNAPL saturation patterns and dissolved-phase contaminant plumes.     

Dependence of lumped mass transfer coefficient on scale and reactions kinetics for biologically enhanced NAPL dissolution

Residual dense non-aqueous liquids (NAPLs) in aquifers constitute a great challenge for groundwater cleanup. Active engineered treatment of regions that contain residual NAPLs is often required to shorten the long-term impact of NAPLs on groundwater quality. Enhanced residual NAPL cleanup can be achieved by promoting biodegradation of NAPL components in the aqueous phase, thereby increasing contaminant fluxes from the NAPL phase. Reaction-enhanced NAPL dissolution is often mathematically simulated under the assumption that lumped mass transfer coefficients, used to describe the dissolution behavior of the NAPL phase, are independent of the reactions. However, this assumption is not warranted because reactions occurring near the water–NAPL interface can reduce characteristic mass transfer lengths, which tend to enhance mass transfer over the no-reaction case.In this study, we mathematically investigated the connections between lumped mass transfer coefficients and reaction kinetics over an idealized residual NAPL domain. Since mass transfer is frequently a scale-dependent process, we also examined the influence of system extent on mass transfer coefficients. For our idealized domain with an assumed first-order decay reaction, the results show that lumped mass transfer coefficients depend on reaction kinetics and system scale. The mass transfer coefficient derived from the non-reactive case cannot properly represent the mass transfer process under the reactive conditions. When the advection time scale is long in comparison to the transverse dispersion time scale in the system, a fast reaction can increase significantly the lumped mass transfer coefficient. The mass transfer coefficient used for simulation was also found to be affected by the nature of the numerical scheme used.     

Integrity of Clay Till Aquitards to DNAPL Migration: Assessment Using Current and Emerging Characterization Tools

Field investigations were carried out to determine the occurrence of tetrachloroethene (PCE) dense nonaqueous phase liquid (DNAPL), the source zone architecture and the aquitard integrity at a 30‐ to 50‐year old DNAPL release site. The DNAPL source zone is located in the clay till unit overlying a limestone aquifer. The DNAPL source zone architecture was investigated through a multiple‐lines‐of‐evidence approach using various characterization tools; the most favorable combination of tools for the DNAPL characterization was geophysical investigations, membrane interface probe, core subsampling with quantification of chlorinated solvents, hydrophobic dye test with Sudan IV, and Flexible Liner Underground Technologies (FLUTe) NAPL liners with activated carbon felt (FACT). While the occurrence of DNAPL was best determined by quantification of chlorinated solvents in soil samples supported by the hydrophobic dye tests (Sudan IV and NAPL FLUTe), the conceptual understanding of source zone architecture was greatly assisted by the indirect continuous characterization tools. Although mobile or high residual DNAPL (S _t > 1%) only occurred in 11% of the source zone samples (intact cores), they comprised 86% of the total PCE mass. The dataset, and associated data analysis, supported vertical migration of DNAPL through fractures in the upper part of the clay till, horizontal migration along high permeability features around the redox boundary in the clay till, and to some extent vertical migration through the fractures in the reduced part of the clay till aquitard to the underlying limestone aquifer. The aquitard integrity to DNAPL migration was found to be compromised at a thickness of reduced clay till of less than 2 m.     

Estimating the Impact of Vadose Zone Sources on Groundwater to Support Performance Assessment of Soil Vapor Extraction

Soil vapor extraction (SVE) is a prevalent remediation remedy for volatile organic compound (VOC) contaminants in the vadose zone. To support selection of an appropriate condition at which SVE may be terminated for site closure or for transition to another remedy, an evaluation is needed to determine whether vadose zone VOC contamination has been diminished sufficiently to keep groundwater concentrations below threshold values. A conceptual model for this evaluation was developed for VOC fate and transport from a vadose zone source to groundwater when vapor‐phase diffusive transport is the dominant transport process. A numerical analysis showed that, for these conditions, the groundwater concentration is controlled by a limited set of parameters, including site‐specific dimensions, vadose zone properties, and source characteristics. On the basis of these findings, a procedure was then developed for estimating groundwater concentrations using results from the three‐dimensional multiphase transport simulations for a matrix of parameter value combinations and covering a range of potential site conditions. Interpolation and scaling processes are applied to estimate groundwater concentrations at compliance (monitoring) wells for specific site conditions of interest using the data from the simulation results. The interpolation and scaling methodology using these simulation results provides a far less computationally intensive alternative to site‐specific three‐dimensional multiphase site modeling, while still allowing for parameter sensitivity and uncertainty analyses. With iterative application, the approach can be used to consider the effect of a diminishing vadose zone source over time on future groundwater concentrations. This novel approach and related simulation results have been incorporated into a user‐friendly Microsoft® Excel®‐based spreadsheet tool entitled SVEET (Soil Vapor Extraction Endstate Tool), which has been made available to the public.     

Dissolution of nonaqueous phase liquid pools in anisotropic aquifers

A two-dimensional numerical transport model is developed to determine the effect of aquifer anisotropy and heterogeneity on mass transfer from a dense nonaqueous phase liquid (DNAPL) pool. The appropriate steady state groundwater flow equation is solved implicitly whereas the equation describing the transport of a sorbing contaminant in a confined aquifer is solved by the alternating direction implicit method. Statistical anisotropy in the aquifer is introduced by two-dimensional, random log-normal hydraulic conductivity field realizations with different directional correlation lengths. Model simulations indicate that DNAPL pool dissolution is enhanced by increasing the mean log-transformed hydraulic conductivity, groundwater flow velocity, and/or anisotropy ratio. The variance of the log-transformed hydraulic conductivity distribution is shown to be inversely proportional to the average mass transfer coefficient.     

DNAPL to LNAPL Transitions During Horizontal Cosolvent Flooding

Cosolvent flooding is a technology with the potential to remove nonaqueous phase liquid (NAPL) sources from the subsurface. It can be used to initiate separate phase mobilization, which allows removal of NAPL within very few pore volumes. Mobilization may result in a sinking DNAPL bank during horizontal flooding of NAPLs denser than water. Reversal of phase density difference between aqueous and DNAPL phases could potentially avoid this downward migration of mobilized DNAPLs. We achieved phase density difference reversal and made DNAPLs float using two components in the cosolvent flooding solution. A low-density cosolvent partitions preferentially into the DNAPL and swells it, which causes a reduction in density of the DNAPL and reversal of the density difference between the NAPL and aqueous phases. A highdensity additive that remains in the aqueous phase allows the cosolvent flooding solution overall to have a density greater than that of water and permits control of the flooding instability. This study focused on tert-butanol as the swelling cosolvent and tetrachloroethylene as the contaminant. In batch tests with sucrose and glycerol as dense additives, phase density difference reversal occurred. To investigate the applicability of phase density difference reversal as a remediation technology, horizontal column and sandbox experiments were performed. These experiments demonstrated the occurrence of phase density difference reversal and effective remediation in horizontal cosolvent floods.     

Modeling unstable alcohol flooding of DNAPL-contaminated columns

Alcohol flooding, consisting of injection of a mixture of alcohol and water, is one source removal technology for dense non-aqueous phase liquids (DNAPLs) currently under investigation. An existing compositional multiphase flow simulator (UTCHEM) was adapted to accurately represent the equilibrium phase behavior of ternary and quaternary alcohol/DNAPL systems. Simulator predictions were compared to laboratory column experiments and the results are presented here. It was found that several experiments involved unstable displacements of the NAPL bank by the alcohol flood or of the alcohol flood by the following water flood. Unstable displacement led to additional mixing compared to ideal displacement. This mixing was approximated by a large dispersion in one-dimensional simulations and or by including permeability heterogeneities on a very small scale in three-dimensional simulations. Three-dimensional simulations provided the best match. Simulations of unstable displacements require either high-resolution grids, or need to consider the mixing of fluids in a different manner to capture the resulting effects on NAPL recovery.     

Modeling metabolic reductive dechlorination in dense non-aqueous phase liquid source-zones

Recent laboratory experimental evidence has suggested that bioremediation may be an attractive management strategy for dense non-aqueous phase liquid (DNAPL) source-zones. In particular, metabolic reductive dechlorination has been shown to reduce aqueous phase chlorinated ethene contaminant concentrations and enhance DNAPL dissolution, reducing source longevity. Transitioning this technology from the laboratory to the field will be facilitated by tools capable of simulating bioenhanced dissolution. This work presents a mathematical model for metabolic reductive dechlorination in a macroscale two-phase (aqueous-organic) environment. The model is implemented through adaptation of an existing multi-phase compositional simulator, which has been modified to incorporate eight chemical components and four microbial populations: a fermentative population, two dechlorinating populations, and a competitor population (e.g., methanogens). Monod kinetics, modified to incorporate electron donor thresholds, electron acceptor competition, and competitor inhibition, are used to simulate microbial growth and component degradation. The developed model is numerically verified and demonstrated through comparisons with published column-scale dechlorination data. Dechlorination kinetics, electron donor concentrations, and DNAPL saturation and distribution are all found to affect the extent of dissolution enhancement, with enhancements ranging from 1.0 to ∼1.9. Comparison of simulation results with those from a simplified analytic modeling approach suggest that the analytical model may tend to over-predict dissolution enhancement and fail to account for the transient nature of dissolution enhancement, leading to significant (70%) under-prediction of source longevity.     

Experimental Investigation of Surfactant-Enhanced Dissolution of Residual NAPL in Saturated Soil

Nonaqueous phase liquid (NAPL) is a long-term source of ground water contamination as the pollutant slowly partitions into the air and water phases. The objective of this work was to study the efficacy of aqueous surfactant solution to enhance the dissolution of a residual NAPL below the capillary fringe, hence reducing the time needed for aquifer restoration. An analytical technique was developed to measure the concentration of NAPL in a nonionic surfactant. Soil column experiments simulated conditions in the saturated soil where a NAPL may become trapped as a discontinuous immobile phase. Experimental results indicate that dissolution was a rate-limited process, approaching equilibrium concentrations after 24 hours. The relative permeability of the aqueous phase initially decreased as surfactant was injected, but increased over time as the saturation of residual NAPL was reduced through mass transfer into the surfactant-enhanced aqueous phase. These findings suggest that enhancing the aqueous phase with a nonionic surfactant may significantly enhance the in situ recovery or residual NAPL.     

A practical model for mobile, residual, and entrapped NAPL in water-wet porous media

Flow of nonvolatile nonaqueous phase liquid (NAPL) and aqueous phases that account for mobile, entrapped, and residual NAPL in variably saturated water-wet porous media is modeled and compared against results from detailed laboratory experiments. Residual saturation formation in the vadose zone is a process that is often ignored in multifluid flow simulators, which might cause an overestimation of the volume of NAPL that reaches the ground water. Mobile NAPL is defined as being continuous in the pore space and flows under a pressure gradient or gravitational body force. Entrapped NAPL is defined as being occluded by the aqueous phase, occurring as immobile ganglia surrounded by aqueous phase in the pore space and formed when NAPL is replaced by the aqueous phase. Residual NAPL is defined as immobile, nonwater entrapped NAPL that does not drain from the pore spaces and is conceptualized as being either continuous or discontinuous. Free NAPL comprises mobile and residual NAPL. The numerical model is formulated on mass conservation equations for oil and water, transported via NAPL and aqueous phases through variably saturated porous media. To account for phase transitions, a primary variable switching scheme is implemented for the oil-mass conservation equation over three phase conditions: (1) aqueous or aqueous-gas with dissolved oil, (2) aqueous or aqueous-gas with entrapped NAPL, and (3) aqueous or aqueous gas with free NAPL. Two laboratory-scale column experiments are modeled to verify the numerical model. Comparisons between the numerical simulations and experiments demonstrate the necessity to include the residual NAPL formation process in multifluid flow simulators.     

Effect of NAPL Source Morphology on Mass Transfer in the Vadose Zone

The generation of vapor‐phase contaminant plumes within the vadose zone is of interest for contaminated site management. Therefore, it is important to understand vapor sources such as non‐aqueous‐phase liquids (NAPLs) and processes that govern their volatilization. The distribution of NAPL, gas, and water phases within a source zone is expected to influence the rate of volatilization. However, the effect of this distribution morphology on volatilization has not been thoroughly quantified. Because field quantification of NAPL volatilization is often infeasible, a controlled laboratory experiment was conducted in a two‐dimensional tank (28 cm × 15.5 cm × 2.5 cm) with water‐wet sandy media and an emplaced trichloroethylene (TCE) source. The source was emplaced in two configurations to represent morphologies encountered in field settings: (1) NAPL pools directly exposed to the air phase and (2) NAPLs trapped in water‐saturated zones that were occluded from the air phase. Airflow was passed through the tank and effluent concentrations of TCE were quantified. Models were used to analyze results, which indicated that mass transfer from directly exposed NAPL was fast and controlled by advective‐dispersive‐diffusive transport in the gas phase. However, sources occluded by pore water showed strong rate limitations and slower effective mass transfer. This difference is explained by diffusional resistance within the aqueous phase. Results demonstrate that vapor generation rates from a NAPL source will be influenced by the soil water content distribution within the source. The implications of the NAPL morphology on volatilization in the context of a dynamic water table or climate are discussed.     

Impact of seasonal variations in hydrological stresses and spatial variations in geologic conditions on a TCE plume at an industrial complex in Wonju,Korea

Spatial and temporal variations in a trichloroethylene (TCE) plume at an industrial complex in Wonju, Korea, were examined based on hydrogeological data and seven rounds of groundwater quality data collected over a year. The site has considerable vertical heterogeneities; the top layer of soil is covered by impermeable paving material at several locations, followed by a series of reclaimed or residual soil layers, and with weathered rocks to the crystalline biotite granite at the bottom. Areal heterogeneity in the surface conditions plays an important role in controlling groundwater recharge. The heterogeneity structure is influenced by complex surface conditions paved with asphalt and concrete. Owing to the presence of limited recharge area and concentrated summer precipitation events, the effects of seasonal variations on groundwater hydraulics tend to diminish with distance from the recharge area. This result was established by analysing the influence of the contrasting surface recharge conditions between the near‐source zone and the far zone, and the seasonally concentrated precipitation on the transport patterns of a TCE plume. In addition, variations in the plume's downstream contaminant flux levels were also analysed along a transect line near the source zone. The results show that the general tendency of the TCE plume contaminant concentration and mass discharges were reproducible if we account for seasonal recharge variations and the associated changes in the groundwater level. During recharge events, the TCE concentration variations appear to be influenced by leaching of the residual dense non‐aqueous‐phase liquid (DNAPL) TCE trapped in the unsaturated zone. This result shows that seasonal variations in contaminant plume near the source zone is inevitable at this site, and that these variations indicate the presence of residual DNAPL at or above the water table, at least in some isolated locations. Copyright © 2011 John Wiley & Sons, Ltd.