Temporal evolution of DNAPL source and contaminant flux distribution: impacts of source mass depletion

We investigated, using model simulations, the changes occurring in the distribution of dense non-aqueous phase liquid (DNAPL) mass (Sn) within the source zone during depletion through dissolution, and the resulting changes in the contaminant flux distribution (J) at the source control plane (CP). Two numerical codes (ISCO3D and T2VOC) were used to simulate selected scenarios of DNAPL dissolution and transport in three-dimensional, heterogeneous, spatially correlated, random permeability fields with emplaced sources. Data from the model simulations were interpreted based on population statistics (mean, standard deviation, coefficient of variation) and spatial statistics (centroid, second moments, variograms). The mean and standard deviation of the Sn and J distributions decreased with source mass depletion by dissolution. The decrease in mean and standard deviation was proportional for the J distribution resulting in a constant coefficient of variation (CV), while for the Sn distribution, the mean decreased faster than the standard deviation. The spatial distributions exhibited similar behavior as the population distribution, i.e., the CP flux distribution was more stable (defined by temporally constant second moments and range of variograms) than the Sn distribution. These observations appeared to be independent of the heterogeneity of the permeability (k) field (variance of the log permeability field=1 and 2.45), correlation structure (positive vs. negative correlation between the k and Sn domains) and the DNAPL dissolution model (equilibrium vs. rate-limited), for the cases studied. Analysis of data from a flux monitoring field study (Hill Air Force Base, Utah) at a DNAPL source CP before and after source remediation also revealed temporal invariance of the contaminant flux distribution. These modeling and field observations suggest that the temporal evolution of the contaminant flux distribution can be estimated if the initial distribution is known. However, the findings are preliminary and broader implications to sampling strategies for remediation performance assessment need to be evaluated in additional modeling and experimental studies.     

Changes in contaminant mass discharge from DNAPL source mass depletion: evaluation at two field sites

Changes in contaminant fluxes resulting from aggressive remediation of dense nonaqueous phase liquid (DNAPL) source zone were investigated at two sites, one at Hill Air Force Base (AFB), Utah, and the other at Ft. Lewis Military Reservation, Washington. Passive Flux Meters (PFM) and a variation of the Integral Pumping Test (IPT) were used to measure fluxes in ten wells installed along a transect down-gradient of the trichloroethylene (TCE) source zone, and perpendicular to the mean groundwater flow direction. At both sites, groundwater and contaminant fluxes were measured before and after the source-zone treatment. The measured contaminant fluxes (J; ML(-2)T(-1)) were integrated across the well transect to estimate contaminant mass discharge (M(D); MT(-1)) from the source zone. Estimated M(D) before source treatment, based on both PFM and IPT methods, were approximately 76 g/day for TCE at the Hill AFB site; and approximately 640 g/day for TCE, and approximately 206 g/day for cis-dichloroethylene (DCE) at the Ft. Lewis site. TCE flux measurements made 1 year after source treatment at the Hill AFB site decreased to approximately 5 g/day. On the other hand, increased fluxes of DCE, a degradation byproduct of TCE, in tests subsequent to remediation at the Hill AFB site suggest enhanced microbial degradation after surfactant flooding. At the Ft. Lewis site, TCE mass discharge rates subsequent to remediation decreased to approximately 3 g/day for TCE and approximately 3 g/day for DCE approximately 1.8 years after remediation. At both field sites, PFM and IPT approaches provided comparable results for contaminant mass discharge rates, and show significant reductions (>90%) in TCE mass discharge as a result of DNAPL mass depletion from the source 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.     

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.     

Assessing impacts of partial mass depletion in DNAPL source zones: II. Coupling source strength functions to plume evolution

Analytical solutions, describing the time-dependent DNAPL source-zone mass and contaminant discharge rate, derived previously in Part I [Falta, R.W., Rao, P.S., Basu, N., this issue. Assessing the impacts of partial mass depletion in DNAPL source zones: I. Analytical modeling of source strength functions and plume response. J. Contam. Hydrol.] are used as a flux-boundary condition in a semi-analytical contaminant transport model. These analytical solutions assume a power relationship between the flow-averaged source concentration, and the source DNAPL mass; the empirical exponent (gamma) is a function of the flow field heterogeneity, DNAPL architecture, and the correlation between them. The DNAPL source strength terms can account for partial source remediation, either at time zero, or at some later time after the DNAPL release. The transport model considers advection, retardation, three-dimensional dispersion, and sequential first-order decay/production of several species. A separate solution is used to compute the time-dependent mass of each contaminant in the plume. A series of examples using different values of gamma shows how the benefits of partial DNAPL source remediation can vary with site conditions. In general, when gamma>1, relatively large short-term reductions in the plume concentrations and mass occur, but the source longevity is not strongly affected. Conversely, when gamma<1, the short-term reductions in the plume concentrations and mass are smaller, but the source longevity can be greatly reduced. In either case, the source remediation effort is much more effective if it is undertaken at an early time, before much contaminant mass has entered the plume. If the remediation effort is significantly delayed, the leading parts of the plume are not affected by the source remediation, and additional control or remediation of the plume itself is required.     

Assessing the impacts of partial mass depletion in DNAPL source zones I. Analytical modeling of source strength functions and plume response

Analytical solutions are developed for approximating the time-dependent contaminant discharge from DNAPL source zones undergoing dissolution and other decay processes. The source functions assume a power relationship between source mass and chemical discharge and can consider partial DNAPL source remediation (depletion) at any time after the initial DNAPL release. The source functions are used as a time-dependent boundary condition in an idealized chemical transport model to develop leading order approximations of the plume response to DNAPL source removal. The results suggest that partial DNAPL remediation does not tend to have a dramatic impact on the maximum extent of the plume if very low concentration values are used to define the plume boundaries. However, the solutions show that partial DNAPL removal from the source zone is likely to lead to large reductions in plume concentrations and mass, and it reduces the longevity of the plume. When the mass discharge from the source zone is linearly related to the DNAPL mass, it is shown that partial DNAPL depletion leads to linearly proportional reductions in the plume mass and concentrations.     

A model for contaminant mass flux in capped sediment under consolidation

The paper presents a model for contaminant transport and flux through a consolidating subaqueous sediment and overlying cap. The formulation is based on the effect of consolidation and excess pore pressure dissipation on transient, nonlinear advective component of transport through sediment and the cap. The consolidation is induced by the buoyant weight of the cap when it is placed on the contaminated sediments. One equation is presented for advective-diffusive transport through the sediment that is dependent upon soil/contaminant properties and transient advective velocity, which is calculated from a second equation based on the Terzaghi consolidation theory. A third equation is provided to describe the transport of contaminants in the cap. The parameters, including advective velocity, and boundary conditions used for contaminant transport through the cap are derived from the solution of the first two equations. The finite difference method is used to solve the system of equations for consolidation and contaminant transport. A hypothetical case is analyzed to demonstrate the formulation, and the results show that advection due to consolidation can accelerate breakthrough of contaminant through the cap by orders of magnitude. The derivation and results show that consolidation should be included for cap design, and that reactive caps are essential for delaying and reducing dissolved contaminant flux.     

Influence of temporally variable groundwater flow conditions on point measurements and contaminant mass flux estimations

Monitoring of contaminant concentrations, e.g., for the estimation of mass discharge or contaminant degradation rates, often is based on point measurements at observation wells. In addition to the problem, that point measurements may not be spatially representative, a further complication may arise due to the temporal dynamics of groundwater flow, which may cause a concentration measurement to be not temporally representative. This paper presents results from a numerical modeling study focusing on temporal variations of the groundwater flow direction. “Measurements” are obtained from point information representing observation wells installed along control planes using different well frequencies and configurations. Results of the scenario simulations show that temporally variable flow conditions can lead to significant temporal fluctuations of the concentration and thus are a substantial source of uncertainty for point measurements. Temporal variation of point concentration measurements may be as high as the average concentration determined, especially near the plume fringe, even when assuming a homogeneous distribution of the hydraulic conductivity. If a heterogeneous hydraulic conductivity field is present, the concentration variability due to a fluctuating groundwater flow direction varies significantly within the control plane and between the different realizations. Determination of contaminant mass fluxes is also influenced by the temporal variability of the concentration measurement, especially for large spacings of the observation wells. Passive dosimeter sampling is found to be appropriate for evaluating the stationarity of contaminant plumes as well as for estimating average concentrations over time when the plume has fully developed. Representative sampling has to be performed over several periods of groundwater flow fluctuation. For the determination of mass fluxes at heterogeneous sites, however, local fluxes, which may vary considerably along a control plane, have to be accounted for. Here, dosimeter sampling in combination with time integrated local water flux measurements can improve mass flux estimates under dynamic flow conditions.     

Simultaneous optimization of dense non-aqueous phase liquid (DNAPL) source and contaminant plume remediation

A framework is developed for simultaneous, optimal design of groundwater contaminant source removal and plume remediation strategies. The framework allows for varying degrees of effort and cost to be dedicated to source removal versus plume remediation. We have accounted for the presence of physical heterogeneity in the DNAPL source, since source heterogeneity controls mass release into the plume and the efficiency of source removal efforts. We considered high and low estimates of capital and operating costs for chemical flushing removal of the source, since these are expected to vary form site to site. Using the lower chemical flushing cost estimates, it is found that the optimal allocation of funds to source removal or plume remediation is sensitive to the degree of heterogeneity in the source. When the time elapsed between the source release and the implementation of remediation was varied, it was found that, except for the longest elapsed time (50,000 days), a combination of partial source removal and plume remediation was most efficient. When first-order, dissolved contaminant degradation was allowed, source removal was found to be unnecessary for the cases where the degradation rate exceeded intermediate values of the first-order rate constant. Finally, it was found that source removal became more necessary as the degree of aquifer heterogeneity increased.     

Backward location and travel time probabilities for a decaying contaminant in an aquifer

Backward location and travel time probabilities can be used to determine the prior position of contamination in an aquifer. These probabilities, which are related to adjoint states of concentration, can be used to improve characterization of known sources of groundwater contamination, to identify previously unknown contamination sources, and to delineate capture zones. The first contribution of this paper is to extend the adjoint model to the case of a decaying solute (first-order decay), and to describe two different interpretations of backward probabilities. The conventional interpretation accounts for the probability that a contaminant particle could decay before reaching the detection location. The other interpretation is conditioned on the fact that the detected contaminant particle actually reached the detection location, despite this possibility of decay. In either case, travel time probabilities are skewed toward earlier travel times, relative to a conservative solute. The second contribution of this paper is to verify the load term for a monitoring well observation. We provide examples using one-dimensional models and hypothetical aquifers. We employ an infinite domain in order to verify the monitoring well load. This new but simple one-dimensional adjoint solution can also be used to verify higher-dimensional numerical models of backward location and travel time probabilities. We employ a semi-infinite domain to illustrate the effect of decay on backward models of pumping well probabilistic capture zones. Decay causes the capture zones to fall closer to the well.     

Interphase mass transfer during chemical oxidation of TCE DNAPL in an aqueous system

The chemical oxidation of trichloroethene dense non-aqueous phase liquid by permanganate was studied in an aqueous system using micro-reaction/extraction vessels in a novel approach. Experiments were conducted at ambient temperature ( approximately 20 degrees C) under static and mixed conditions to evaluate the rate of TCE(DNAPL) dissolution as a function of permanganate concentration. Chemical oxidation by permanganate was shown to increase the rate of TCE(DNAPL) dissolution under static conditions and decrease the rate of dissolution under mixed conditions. The apparent inconsistency in results appears to result from the local deposition of a film at the DNAPL interface composed of manganese oxide solids as discovered through visual observation with the aid of a Goniometer. Data from interfacial deposition experiments suggest that the film formed rapidly and reached maturation within approximately 2 h with little or no growth occurring thereafter. A conceptual model of the reaction and mass transfer processes occurring at the DNAPL interface was proposed based on the experimental results.     

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.     

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.     

Analytical solutions of contaminant transport in homogeneous and isotropic aquifer in three-dimensional groundwater flow

Environmental Science and Pollution Research - Modeling three-dimensional contaminant transport released from arbitrary shape source geometries is useful in hydrological and environmental sciences....     

A large-scale experiment on mass transfer of trichloroethylene from the unsaturated zone of a sandy aquifer to its interfaces

A large-scale experiment was conducted to investigate the transport of trichloroethylene (TCE) vapors in the unsaturated zone and to determine the mass transfer to the groundwater and the atmosphere. The experiment involved injection of 5 1 of TCE in the unsaturated zone under controlled conditions, with multidepth sampling of gas and water through the unsaturated zone and across the capillary zone into underlying groundwater. The mass transfer of TCE vapors from the vadose zone to the atmosphere was quantified using a vertical flux chamber. A special soil water sampler was used to monitor transport across the capillary fringe. Experimental data indicated that TCE in the unsaturated zone was mainly transported to the atmosphere and this exchange reduced significantly the potential for groundwater pollution. The maximum measured TCE flux to the atmosphere was about 3 g/m(2)/day. Observed and calculated fluxes based on vertical TCE vapor concentration gradients and Fick's law were in good agreement. This confirms that TCE vapor transport under the experimental conditions was governed essentially by molecular diffusion. TCE vapors also caused a lower, but significant contamination of the underlying groundwater by dispersion across the capillary fringe with a corresponding maximum flux of about 0.1 g/m(2)/day. This mass transfer to groundwater is partly uncertain due to an inadvertent entry of some nonaqueous phase liquid (NAPL) from the source area into the saturated zone. Application of an analytical solution to estimate the TCE flux from the unsaturated zone to the groundwater indicated that this phenomenon is not only influenced by molecular diffusion but also by vertical dispersion. The mass balance indicates that, under the given experimental conditions (e.g. proximity of the source emplacement relative to the soil surface, relatively high permeable porous medium), nearly 95% of the initial TCE mass was transferred to the atmosphere.     

Sublethal contaminant exposure alters behavior in a common insect: Important implications for trophic transfer

This study examined the effects of sub-lethal exposure of the ubiquitous pesticide malathion on the behavior of the model orthopteran species, the house cricket (Acheta domesticus). Increasing concentrations of malathion caused male crickets to increase periods of non-directional movement, such as twitching and grooming, directional movement, and to seek out shelter less. These are all behavioral alterations that may increase the cricket's chances of being preyed upon, and thus have the potential for serious ecological consequences through trophic transfer. This study also revealed that female crickets appeared to be less affected by malathion than their male conspecifics, indicating a potential sex-bias in both susceptibility and possible predator attack.     

Field scale characterization and modeling of contaminant release from a coal tar source zone

A coal tar contaminated site was characterized using traditional and innovative investigation methods. A careful interpretation of hydrogeological and hydrogeochemical data allowed for the conceptualization of the heterogeneous coal tar distribution in the subsurface. Past and future contaminant release from the source zone was calculated using a modeling framework consisting of a three-dimensional steady-state groundwater flow model (MODFLOW) and two hydrogeochemical models (MIN3P). Computational time of long-term simulations was reduced by simplifying the coal tar composition using 3 composite and 2 individual constituents and sequential application of a 2D centerline model (for calibration and predictions) and a 3D model (only for predictions). Predictions were carried out for a period of 1000 years. The results reveal that contaminant mass flux is governed by the geometry of zones containing residual coal tar, amount of coal tar, its composition and the physicochemical properties of the constituents. The long-term predictions made using the 2D model show that even after 1000 years, source depletion will be small with respect to phenanthrene, 89% of initial mass will be still available, and for the moderately and sparingly soluble composite constituents, 60% and 98%, respectively.     

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.     

DNAPL distribution in the source zone: effect of soil structure and uncertainty reduction with increased sampling density

This paper focuses on parameters describing the distribution of dense nonaqueous phase liquid (DNAPL) contaminants and investigates the variability of these parameters that results from soil heterogeneity. In addition, it quantifies the uncertainty reduction that can be achieved with increased density of soil sampling. Numerical simulations of DNAPL releases were performed using stochastic realizations of hydraulic conductivity fields generated with the same geostatistical parameters and conditioning data at two sampling densities, thus generating two simulation ensembles of low and high density (three-fold increase) of soil sampling. The results showed that DNAPL plumes in aquifers identical in a statistical sense exhibit qualitatively different patterns, ranging from compact to finger-like. The corresponding quantitative differences were expressed by defining several alternative measures that describe the DNAPL plume and computing these measures for each simulation of the two ensembles. The uncertainty in the plume features under study was affected to different degrees by the variability of the soil, with coefficients of variation ranging from about 20% to 90%, for the low-density sampling. Meanwhile, the increased soil sampling frequency resulted in reductions of uncertainty varying from 7% to 69%, for low- and high-uncertainty variables, respectively. In view of the varying uncertainty in the characteristics of a DNAPL plume, remedial designs that require estimates of the less uncertain features of the plume may be preferred over others that need a more detailed characterization of the source zone architecture.