PCE dissolution and simultaneous dechlorination by nanoscale zero-valent iron particles in a DNAPL source zone

While the capability of nanoscale zero-valent iron (NZVI) to dechlorinate organic compounds in aqueous solutions has been demonstrated, the ability of NZVI to remove dense non-aqueous phase liquid (DNAPL) from source zones under flow-through conditions similar to a field scale application has not yet been thoroughly investigated. To gain insight on simultaneous DNAPL dissolution and NZVI-mediated dechlorination reactions after direct placement of NZVI into a DNAPL source zone, a combined experimental and modeling study was performed. First, a DNAPL tetrachloroethene (PCE) source zone with emplaced NZVI was built inside a small custom-made flow cell and the effluent PCE and dechlorination byproducts were monitored over time. Second, a model for rate-limited DNAPL dissolution and NZVI-mediated dechlorination of PCE to its three main reaction byproducts with a possibility for partitioning of these byproducts back into the DNAPL was formulated. The coupled processes occurring in the flow cell were simulated and analyzed using a detailed three-dimensional numerical model. It was found that subsurface emplacement of NZVI did not markedly accelerate DNAPL dissolution or the DNAPL mass-depletion rate, when NZVI at a particle concentration of 10g/L was directly emplaced in the DNAPL source zone. To react with NZVI the DNAPL PCE must first dissolve into the groundwater and the rate of dissolution controls the longevity of the DNAPL source. The modeling study further indicated that faster reacting particles would decrease aqueous contaminant concentrations but there is a limit to how much the mass removal rate can be increased by increasing the dechlorination reaction rate. To ensure reduction of aqueous contaminant concentrations, remediation of DNAPL contaminants with NZVI should include emplacement in a capture zone down-gradient of the DNAPL source.     

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.     

Convergent-flow tracer tests in heterogeneous media: combined experimental-numerical analysis for determination of equivalent transport parameters

In modeling transport within naturally heterogeneous aquifers, it is usually assumed that the transport equations valid at local scales can also be applied at larger scales. At larger scales, the heterogeneous domain is represented by an equivalent homogeneous medium. Convergent-flow tracer tests constitute one of the most frequently used field tests to estimate effective input parameters of equivalent homogeneous aquifers. Traditionally, statistical approaches applied to groundwater flow and solute transport have provided tools to estimate these equivalent parameters. These approaches are based on a number of simplifications including the assumption that the point transmissivity values follow a multilog-normal random function. Several investigators have found that this assumption may not be valid in many field cases. In order to study the applicability of the equivalent homogeneous formulation in a nontraditional stochastic field, a number of experimental and numerical studies were conducted. The results are used to determine the apparent values of porosity and dispersivity that would be obtained if convergent-flow tracer tests were conducted in a deterministically generated heterogeneous transmissivity field displaying anisotropy in the correlation structure. It is shown that in this particular heterogeneous media, apparent porosity strongly depends on connectivity rather than on transmissivity. This dependence on connectivity questions the theoretical results obtained in continuum equivalent fields to estimate effective porosity.     

Investigation of surfactant-enhanced dissolution of entrapped nonaqueous phase liquid chemicals in a two-dimensional groundwater flow field.

Because of their low solubility, waste chemicals in the form of nonaqueous phase liquids (NAPLs) that are entrapped in subsurface formations act as long-term sources of groundwater contamination. In the design of remediation schemes that use surfactants, it is necessary to estimate the mass transfer rate coefficients under multi-dimensional flow fields that exit at field sites. In this study, we investigate mass transfer under a two-dimensional flow field to obtain an understanding of the basic mechanisms of surfactant-enhanced dissolution and to quantify the mass transfer rates. Enhanced dissolution experiments in a two-dimensional test cell were conducted to measure rates of mass depletion from entrapped NAPLs to a flowing aqueous phase containing a surfactant. In situ measurement of transient saturation changes using a gamma attenuation system revealed dissolution patterns that are affected by the dimensionality of the groundwater flow field. Numerical modeling of local flow fields that changed with time, due to depletion of NAPL sources, enabled the examination of the basic mechanisms of NAPL dissolution in complex groundwater systems. Through nonlinear regression analysis, mass transfer rates were correlated to porous media properties, NAPL saturation and aqueous phase velocity. Results from the experiments and numerical analyses were used to identify deficiencies in existing methods of analysis that uses assumptions of one-dimensional flow, homogeneity of aquifer properties, local equilibrium and idealized transient mass transfer.     

Empirical correlations to estimate agglomerate size and deposition during injection of a polyelectrolyte-modified Fe0 nanoparticle at high particle concentration in saturated sand

Controlled emplacement of polyelectrolyte-modified nanoscale zerovalent iron (NZVI) particles at high particle concentration (1-10 g/L) is needed for effective in situ subsurface remediation using NZVI. Deep bed filtration theory cannot be used to estimate the transport and deposition of concentrated polyelectrolyte-modified NZVI dispersions (>0.03 g/L) because particles agglomerate during transport which violates a fundamental assumption of the theory. Here we develop two empirical correlations for estimating the deposition and transport of concentrated polyelectrolyte-modified NZVI dispersions in saturated porous media when NZVI agglomeration in porous media is assumed to reach steady state quickly. The first correlation determines the apparent stable agglomerate size formed during NZVI transport in porous media for a fixed hydrogeochemical condition. The second correlation estimates the attachment efficiency (sticking coefficient) of the stable agglomerates. Both correlations are described using dimensionless numbers derived from parameters affecting deposition and agglomeration in porous media. The exponents for the dimensionless numbers are determined from statistical analysis of breakthrough data for polyelectrolyte-modified NZVI dispersions collected in laboratory scale column experiments for a range of ionic strength (1, 10, and 50mM Na(+) and 0.25, 1, and 1.25 mM Ca(2+)), approach velocity (0.8 to 55 × 10(-4)m/s), average collector sizes (d(50)=99 μm, 300 μm, and 880 μm), and polyelectrolyte surface modifier properties. Attachment efficiency depended on approach velocity and was inversely related to collector size, which is contrary to that predicted from classic filtration models. High ionic strength, the presence of divalent cations, lower extended adsorbed polyelectrolyte layer thickness, decreased approach velocity, and a larger collector size promoted NZVI agglomeration and deposition and thus limited its mobility in porous media. These effects are captured quantitatively in the two correlations developed. The application and limitations of using the correlations for preliminary design of in situ NZVI emplacement strategies is discussed.     

Intermediate-scale 2D experimental investigation of in situ chemical oxidation using potassium permanganate for remediation of complex DNAPL source zones

In situ chemical oxidation is a technology that has been applied to speed up remediation of a contaminant source zone by inducing increased mass transfer from DNAPL sources into the aqueous phase for subsequent destruction. The DNAPL source zone can consist of one or more individual sources that may be present as an interconnected pool of high saturation, as a region of disconnected ganglia at residual saturation, or as combinations of these two morphologies. Potassium permanganate (KMnO(4)) is a commonly employed oxidant that has been shown to rapidly destroy DNAPL compounds like PCE and TCE following second-order kinetics in an aqueous system. During the oxidation of a target DNAPL compound, or naturally occurring reduced species in the subsurface, manganese oxide (MnO(2)) solids are produced. Research has shown that these manganese oxide solids may result in permeability reductions in the porous media thus reducing the ability for oxidant to be transported to individual DNAPL sources. It can also occur at the DNAPL-water interface, decreasing contact of the oxidant with the DNAPL. Additionally, MnO(2) formation at the DNAPL-water interface, and/or flow-bypassing as a result of permeability reductions around the source, may alter the mass transfer from the DNAPL into the aqueous phase, potentially diminishing the magnitude of any DNAPL mass depletion rate increase induced by oxidation. An experiment was performed in a two-dimensional (2D) sand-filled tank that included several discrete DNAPL source zones. Spatial and temporal monitoring of aqueous PCE, chloride, and permanganate concentrations was used to relate changes in mass depletion of, and mass flux, from DNAPL residual and pool source zones to chemical oxidation performance and MnO(2) formation. During the experiment, permeability changes were monitored throughout the 2D tank and these were related to MnO(2) deposition as measured through post-oxidation soil coring. Under the conditions of this experiment, MnO(2) formation was found to reduce permeability in and around DNAPL source zones resulting in changes to the overall flow pattern, with the effects depending on source zone configuration. A pool with little or no residual around it, in a relatively homogeneous flow field, appeared to benefit from resulting MnO(2) pore-blocking that substantially reduced mass transfer from the pool even though there was relatively little PCE mass removed from the pool. In contrast, a pool with residual around it (in a more typical heterogeneous flow field) appeared to undergo increased mass transfer as MnO(2) reduced permeability, altering the water flow and increasing the mixing at the DNAPL-water interface. Further, the magnitude of increased PCE mass depletion during oxidation appeared to depend on the PCE source configuration (pool versus ganglia) and decreased as MnO(2) was formed and deposited at the DNAPL-water interface. Overall, the oxidation of PCE mass appeared to be rate-limited by the mass transfer from the DNAPL to aqueous phase.     

The significance of heterogeneity on mass flux from DNAPL source zones: an experimental investigation

Understanding the process of mass transfer from source zones of aquifers contaminated with organic chemicals in the form of dense non-aqueous phase liquids (DNAPL) is of importance in site management and remediation. A series of intermediate-scale tank experiments was conducted to examine the influence of aquifer heterogeneity on DNAPL mass transfer contributing to dissolved mass emission from source zone into groundwater under natural flow before and after remediation. A Tetrachloroethylene (PCE) spill was performed into six source zone models of increasing heterogeneity, and both the spatial distribution of the dissolution behavior and the net effluent mass flux were examined. Experimentally created initial PCE entrapment architecture resulting from the PCE migration was largely influenced by the coarser sand lenses and the PCE occupied between 30 and 60% of the model aquifer depth. The presence of DNAPL had no apparent effect on the bulk hydraulic conductivity of the porous media. Up to 71% of PCE mass in each of the tested source zone was removed during a series of surfactant flushes, with associated induced PCE mobilization responsible for increasing vertical DNAPL distributions. Effluent mass flux due to water dissolution was also found to increase progressively due to the increase in NAPL-water contact area even though the PCE mass was reduced. Doubling of local groundwater flow velocities showed negligible rate-limited effects at the scale of these experiments. Thus, mass transfer behavior was directly controlled by the morphology of DNAPL within each source zone. Effluent mass flux values were normalized by the up-gradient DNAPL distributions. For the suite of aquifer heterogeneities and all remedial stages, normalized flux values fell within a narrow band with mean of 0.39 and showed insensitivity to average source zone saturations.     

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.     

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.