Search

Laboratory batch experiments were performed to assess the reduction of trichloroethene (TCE) and oxygen via natural ferrous minerals. TCE reduction under anoxic conditions was measured via the generation of reduced gases, while oxygen reduction via the generation of hydroxyl radicals was measured as a surrogate for potential TCE oxidation. Results showed that TCE reduction under anoxic conditions was observed for ankerite, siderite, and illite, but not for biotite; acetylene was the primary identified dechlorination product. With the exception of biotite, first-order dechlorination rate constants increased with increasing ferrous content of the mineral, with rate constants ranging from 3.1 × 10−8 to 4.8 10−7 L g−1 d−1. Measured reduction potentials (mV vs SHE) ranged from −104 for illite to +84 for biotite. When normalizing measured first-order dechlorination rate constants to the estimated ferrous iron mineral specific surface area (where surface area was based on nitrogen adsorption analysis of the minerals), TCE dechlorination rate constants increased with increasing reduction potentials. Under oxic conditions, hydroxyl radicals were generated with each of the four minerals. However, mineral activity showed no readily apparent correlation to ferrous content or mineral surface area. In terms of TCE and oxygen reduced per mole of ferrous iron initially present in each mineral, illite was the most reactive of the four minerals. Together, these results suggest that several ferrous minerals may contribute to abiotic dechlorination in the natural environment, and (at least for TCE reduction under anoxic conditions) measurement of ferrous mineral content and reduction potential may serve as useful tools for estimating TCE first-order abiotic dechlorination rate constants.

Low permeability zones (LPZs) are major sources of groundwater contamination after active remediation to remove pollutants in adjacent high permeability zones (HPZs). Slow back diffusion from LPZs to HPZs can extend management of polluted sites by decades. Numerical models are often used to simulate back diffusion, estimate cleanup times, and develop site management strategies. Sharp concentration gradients of pollutants are present at the interface between HPZs and LPZs, and hence accurate simulation requires fine grid sizes resulting in high computational burden. Since the MODFLOW family of codes is widely used in practice, we develop a new approach for modeling pollutant back diffusion using MODFLOW/RT3D that eliminates the need for fine discretization of the LPZ. Instead, the LPZ is treated as an impermeable region in MODFLOW, while in RT3D the LPZ is conceptualized as a series of immobile zones coupled with a mobile zone at the HPZ/LPZ interface. Finite volume discretization of diffusion and reaction within the LPZ is then modeled as mass transfer and reaction among several immobile species. This results in a simulation domain with significantly fewer grid cells compared to that required if all LPZs are discretized, providing potential for improved computational efficiency. Cases, including a layer of HPZ over an LPZ, a thin/thick lens of LPZ embedded in HPZ, and multiple lens of LPZs embedded in HPZ are tested by the new approach for tracer and reactive scenarios.

Laboratory batch experiments were performed to assess the impacts of temperature and mineralogy on the abiotic dechlorination of tetrachloroethene (PCE) or trichloroethene (TCE) due to the presence of ferrous minerals in natural aquifer clayey soils under anaerobic conditions. A combination of x-ray diffraction (XRD), magnetic susceptibility, and ferrous mineral content were used to characterize each of the 3 natural soils tested in this study, and dechlorination at temperatures ranging from 20 to 55°C were examined. Results showed that abiotic dechlorination occurred in all 3 soils examined, yielding reduced gas abiotic dechlorination products acetylene, butane, ethene, and/or propane. Bulk first-order dechlorination rate constants (kbulk), scaled to the soil:water ratio expected for in situ conditions, ranged from 2.0×10-5day-1 at 20°C, to 32×10-5day-1 at 55°C in the soil with the greatest ferrous mineral content. For the generation of acetylene and ethene from PCE, the reaction was well described by Arrhenius kinetics, with an activation energy of 91kJ/mol. For the generation of coupling products butane and propane, the Arrhenius equation did not provide a satisfactory description of the data, likely owing to the complex reaction mechanisms associated with these products and/or diffusional mass transfer processes associated with the ferrous minerals likely responsible for these coupling reactions. Although the data set was too limited to determine a definitive correlation, the two soils with elevated ferrous mineral contents had elevated abiotic dechlorination rate constants, while the one soil with a low ferrous mineral content had a relatively low abiotic dechlorination rate constant. Overall, results suggest intrinsic abiotic dechlorination rates may be an important long-term natural attenuation component in site conceptual models for clays that have the appropriate iron mineralogy.

Physical, chemical, and biological interactions between groundwater and sedimentary rock directly control the fundamental subsurface properties such as porosity, permeability, and flow. This is true for a variety of subsurface scenarios, ranging from shallow groundwater aquifers to deeply buried hydrocarbon reservoirs. Microfluidic flow cells are now commonly being used to study these processes at the pore scale in simplified pore structures meant to mimic subsurface reservoirs. However, these micromodels are typically fabricated from glass, silicon, or polydimethylsiloxane (PDMS), and are therefore incapable of replicating the geochemical reactivity and complex three-dimensional pore networks present in subsurface lithologies. To address these limitations, we developed a new microfluidic experimental test bed, herein called the Real Rock-Microfluidic Flow Cell (RR-MFC). A porous 500μm-thick real rock sample of the Clair Group sandstone from a subsurface hydrocarbon reservoir of the North Sea was prepared and mounted inside a PDMS microfluidic channel, creating a dynamic flow-through experimental platform for real-time tracking of subsurface reactive transport. Transmitted and reflected microscopy, cathodoluminescence microscopy, Raman spectroscopy, and confocal laser microscopy techniques were used to (1) determine the mineralogy, geochemistry, and pore networks within the sandstone inserted in the RR-MFC, (2) analyze non-reactive tracer breakthrough in two- and (depth-limited) three-dimensions, and (3) characterize multiphase flow. The RR-MFC is the first microfluidic experimental platform that allows direct visualization of flow and transport in the pore space of a real subsurface reservoir rock sample, and holds potential to advance our understandings of reactive transport and other subsurface processes relevant to pollutant transport and cleanup in groundwater, as well as energy recovery.

Low permeability source zones sustain long-term trichloroethene (TCE) groundwater contamination. In anaerobic environments, TCE is transformed by both biological reductive dechlorination and abiotic reactions with reactive minerals. Little is known about the relative contribution of these two pathways as TCE diffuses from low permeability zones (LPZs) into high permeability zones (HPZs). This study combines a flow cell experiment, batch experiments, and a diffusion-reaction model to evaluate the contributions of biotic and abiotic TCE transformation in LPZs. Natural clay (LPZ) and sand (HPZ) from a former Air Force base were used in all experiments. In batch, the LPZ material transformed TCE and cis-1,2-dichloroethene (cis-DCE) to acetylene with pseudo first-order rate constants of 8.57 × 10−6 day−1 and 1.02 × 10−6 day−1, respectively. Biotic and abiotic pathways were then evaluated together in a bench-scale flow cell (16.5 cm × 2 cm × 16.5 cm) that contained a LPZ layer, with a source of TCE at the base, overlain by a HPZ continuously purged with lactate-amended groundwater. Diffusion controlled mass transfer in the LPZ, while advection controlled migration in the HPZ. The mass discharge rate of TCE and its biotic (cis-DCE and vinyl chloride) and abiotic (acetylene) transformation products were measured over 180 days in the flow cell effluent. Depth profiles of these compounds through the LPZ were determined after terminating the experiment. A one-dimensional diffusion-reaction model was used to interpret the effluent and depth profile data and constrain reaction parameters. Abiotic transformation rate constants for TCE to acetylene, normalized to in situ solids loading, were approximately 13 times greater in batch than in the flow cell. Slower transformation rates in the flow cell indicate elevated TCE concentration and/or further degradation of acetylene to other reduced gas compounds in the flow cell LPZ (thereby partially masking TCE abiotic transformation). Biotic and abiotic parameters used to interpret the flow cell data were then used to simulate a field site with a 300 cm thick LPZ. Abiotic processes contributed to a 2% reduction in TCE flux after 730 days. When abiotic rate constants were changed to that observed in batch, or to rate constants previously reported for a pyrite rich mudstone, the TCE flux reduction was 21% and 53%, respectively, after 730 days. Though biotic processes dominated TCE transformation in the flow cell experiment, the simulations indicate that abiotic processes have potential to significantly contribute to TCE attenuation in electron donor limited environments provided suitable reactive minerals are present.

Predicting the longevity of non-aqueous phase liquid (NAPL) source zones has proven to be a difficult modeling problem that has yet to be resolved. Research efforts towards understanding NAPL depletion have focused on developing empirical models that relate lumped mass transfer rates to velocities and organic saturations. These empirical models are often unable to predict NAPL dissolution for systems different from those used to calibrate them, indicating that system-specific factors important for dissolution are not considered. This introduces the need for a calibration step before these models can be reliably used to predict NAPL dissolution for systems of arbitrary characteristics. In this paper, five published Sherwood–Gilland models are evaluated using experimental observations from the dissolution of two laboratory-scale complex three-dimensional NAPL source zones. It is shown that the relative behavior of the five models depends on the system and source zone characteristics. Through a theoretical analysis, comparing Sherwood–Gilland type models to a process-based, thermodynamic dissolution model, it is shown that the coefficients of the Sherwood–Gilland models can be related to measurable soil properties. The derived dissolution model with soil-dependent coefficients predicts concentrations identical to those predicted by the thermodynamic dissolution model for cases with negligible hysteresis. This correspondence breaks down when hysteresis has a significant impact on interfacial areas. In such cases, the derived dissolution model will slightly underestimate dissolved concentrations at later times, but is more likely to capture system-specific dissolution rates than Sherwood–Gilland models.

Contaminant hydrogeological processes occurring in porous media are typically not amenable to direct observation. As a result, indirect measurements (e.g., contaminant breakthrough at a fixed location) are often used to infer processes occurring at different scales, locations, or times. To overcome this limitation, non-invasive imaging methods are increasingly being used in contaminant hydrogeology research. Four of the most common methods, and the subjects of this review, are optical imaging using UV or visible light, dual-energy gamma radiation, X-ray microtomography, and magnetic resonance imaging (MRI). Non-invasive imaging techniques have provided valuable insights into a variety of complex systems and processes, including porous media characterization, multiphase fluid distribution, fluid flow, solute transport and mixing, colloidal transport and deposition, and reactions. In this paper we review the theory underlying these methods, applications of these methods to contaminant hydrogeology research, and methods' advantages and disadvantages. As expected, there is no perfect method or tool for non-invasive imaging. However, optical methods generally present the least expensive and easiest options for imaging fluid distribution, solute and fluid flow, colloid transport, and reactions in artificial two-dimensional (2D) porous media. Gamma radiation methods present the best opportunity for characterization of fluid distributions in 2D at the Darcy scale. X-ray methods present the highest resolution and flexibility for three-dimensional (3D) natural porous media characterization, and 3D characterization of fluid distributions in natural porous media. And MRI presents the best option for 3D characterization of fluid distribution, fluid flow, colloid transport, and reaction in artificial porous media. Obvious deficiencies ripe for method development are the ability to image transient processes such as fluid flow and colloid transport in natural porous media in three dimensions, the ability to image many reactions of environmental interest in artificial and natural porous media, and the ability to image selected processes over a range of scales in artificial and natural porous media.

Nonaqueous phase liquid (NAPL) dissolution was studied in three-dimensional (3D) heterogeneous experimental aquifers (25.5 cm x 9 cm x 8.5 cm) with two different longitudinal correlation lengths (2.1 cm and 1.1 cm) and initial spill volumes (22.5 ml and 10.5 ml). Spatial and temporal distributions of NAPL during dissolution were measured using magnetic resonance imaging (MRI). At high NAPL spill volume, average effluent concentrations initially increased during dissolution, as NAPL pools transitioned to NAPL ganglia, and then decreased as the total NAPL-water interfacial area decreased over time. Experimental results were used to test six dissolution models: (i and ii) a one-dimensional (1D) model using either specific NAPL-water interfacial area values estimated from MR images at each time step (i.e., 1D quasi-steady state model), or an empirical mass transfer (Sh') correlation (i.e., 1D transient model), (iii and iv) a multiple analytical source superposition technique (MASST) using either the NAPL distribution determined from MR images at each time step (i.e., MASST steady state model), or the NAPL distribution determined from mass balance calculations (i.e., MASST transient model), (v) an equilibrium streamtube model, and (vi) a 3D grid-scale pool dissolution model (PDM) with a dispersive mass flux term. The 1D quasi-steady state model and 3D PDM captured effluent concentration values most closely, including some concentration fluctuations due to changes in the extent of flow reduction. The 1D transient, MASST steady state and transient, and streamtube models all showed a monotonic decrease in effluent concentration values over time, and the streamtube model was the most computationally efficient. Changes during dissolution of the effective NAPL-water interfacial area estimated from imaging data are similar to changes in effluent concentration values. The 1D steady state model incorporates estimates of the effective NAPL-water interfacial area directly at each time point; the 3D PDM does so indirectly through mass balance and a relative permeability function, which causes reduced water flow through high saturation NAPL regions. Hence, when model accuracy is required, the results indicate that a surrogate of this effective interfacial area is required. Approaches to include this surrogate in the MASST and streamtube models are recommended.

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.

Mass transfer between aquifer material and groundwater is often modeled as first-order rate-limited sorption or diffusive exchange between mobile zones and immobile zones with idealized geometries. Recent improvements in experimental techniques and advances in our understanding of pore-scale heterogeneity demonstrate that two (or even a few) rate coefficients are insufficient in many cases. Here, we investigate a piece-wise linear model for a continuous distribution of rate coefficients, that has several advantages over previously used `statistical' distribution models (with functional form from gamma or lognormal PDF's): (1) distributions of arbitrary, even bimodal, shapes can be represented; (2) linear estimation methods can be applied to determine the distribution from experimental data; (3) the uncertainty in the distribution can be determined for each of its sections; and (4) the relationship between the time scales of available data and those of estimatable mass transfer processes can be investigated. A statistical model refinement algorithm is presented that reduces the number of parameters (sections of the piece-wise linear model) to the admissible minimum. We show that purging experiments allow estimation of a wider zone of the rate distribution than do batch experiments, and hence will provide predictions that are accurate over a wider range of time scales. Finally, in an application to TCE gas-purging desorption data, the piece-wise linear rate-distribution model has a higher probability of being adequate than those using a gamma or lognormal distribution or a single rate coefficient.

The purpose of this work is to identify the mechanisms that govern the removal of carbon tetrachloride (CT) during soil vapor extraction (SVE) by comparing numerical and analytical model simulations with a detailed data set from a well-defined intermediate-scale flow cell experiment. The flow cell was packed with a fine-grained sand layer embedded in a coarse-grained sand matrix. A total of 499 mL CT was injected at the top of the flow cell and allowed to redistribute in the variably saturated system. A dual-energy gamma radiation system was used to determine the initial NAPL saturation profile in the fine-grained sand layer. Gas concentrations at the outlet of the flow cell and 15 sampling ports inside the flow cell were measured during subsequent CT removal using SVE. Results show that CT mass was removed quickly in coarse-grained sand, followed by a slow removal from the fine-grained sand layer. Consequently, effluent gas concentrations decreased quickly at first, and then started to decrease gradually, resulting in long-term tailing. The long-term tailing was mainly due to diffusion from the fine-grained sand layer to the coarse-grained sand zone. An analytical solution for a one-dimensional advection and a first-order mass transfer model matched the tailing well with two fitting parameters. Given detailed knowledge of the permeability field and initial CT distribution, we were also able to predict the effluent concentration tailing and gas concentration profiles at sampling ports using a numerical simulator assuming equilibrium CT evaporation. The numerical model predictions were accurate within the uncertainty of independently measured or literature derived parameters. This study demonstrates that proper numerical modeling of CT removal through SVE can be achieved using equilibrium evaporation of NAPL if detailed fine-scale knowledge of the CT distribution and physical heterogeneity is incorporated into the model. However, CT removal could also be fit by a first-order mass transfer analytical model, potentially leading to an erroneous conclusion that the long-term tailing in the experiment was kinetically controlled due to rate-limited NAPL evaporation.

The amount, location, and form of NAPL in contaminated vadose zones are controlled by the spatial distribution of water saturation and soil permeability, the NAPL spill scenario, water infiltration events, and vapor transport. To evaluate the effects of these processes, we used the three-phase flow simulator STOMP, which includes a new permeability-liquid saturation-capillary pressure (k-S-P) constitutive model. This new constitutive model considers three NAPL forms: free, residual, and trapped. A 2-D vertical cross-section with five stratigraphic layers was assumed, and simulations were performed for seven cases. The conceptual model of the soil heterogeneity was based upon the stratigraphy at the Hanford carbon tetrachloride (CT) spill site. Some cases considered co-disposal of NAPL with large volumes of wastewater, as also occurred at the Hanford CT site. In these cases, the form and location of NAPL were most strongly influenced by high water discharge rates and NAPL evaporation to the atmosphere. In order to investigate the impact of heterogeneity, the hydraulic conductivity within the lower permeability layer was modeled as a realization of a random field having three different classes. For six extreme cases of 100 realizations, the CT mass that reached the water table varied by a factor of two, and was primarily controlled by the degree of lateral connectivity of the low conductivity class within the lowest permeability layer. The grid size at the top boundary had a dramatic impact on NAPL diffusive flux just after the spill event when the NAPL was present near the ground surface. NAPL evaporation with a fine grid spacing at the top boundary decreased CT mass that reached the water table by 74%, compared to the case with a coarse grid spacing, while barometric pumping had a marginal effect for the case of a continuous NAPL spill scenario considered in this work. For low water infiltration rate scenarios, the distribution of water content prior to a NAPL spill event decreased CT mass that reached the water table by 98% and had a significant impact on the formation of trapped NAPL. For all cases simulated, use of the new constitutive model that allows the formation of residual NAPL increased the amount of NAPL retained in the vadose zone. Density-driven advective gas flow from the ground surface controlled vapor migration in strongly anisotropic layers, causing NAPL mass flux to the lower layer to be reduced. These simulations indicate that consideration of the formation of residual and trapped NAPLs and dynamic boundary conditions (e.g., areas, rates, and periods of different NAPL and water discharge and fluctuations of atmospheric pressure) in the context of full three-phase flow are needed, especially for NAPL spill events at the ground surface. In addition, NAPL evaporation, density-driven gas advection, and NAPL vertical movement enhanced by water flow must be considered in order to predict NAPL distribution and migration in the vadose zone.

Soil vapor extraction (SVE) is commonly used to remediate nonaqueous phase liquids (NAPLs) from the vadose zone. This paper aims to determine the effect of grain size heterogeneity on the removal of NAPL in porous media during SVE. Magnetic resonance imaging (MRI) was used to observe and quantify the amount and location of NAPL in flow-through columns filled with silica gel grains. MRI is unique because it is nondestructive, allowing three-dimensional images to be taken of the phases as a function of space and time. Columns were packed with silica gel in three ways: coarse grains (250-550 microm) only, fine grains (32-63 microm) only, and a core of fine grains surrounded by a shell of coarse grains. Columns saturated with water were drained under a constant suction head, contaminated with decane, and then drained to different decane saturations. Each column was then continuously purged with water-saturated nitrogen gas and images were taken intermittently. Results showed that at residual saturation, a sharp volatilization front moved through the columns filled with either coarse-grain or fine-grain silica gel. In the heterogeneous columns, the volatilization front in the core lagged just behind the shell because gas flow was greater through the shell and decane in the core diffused outward to the shell. When decane saturation in the core was above residual saturation, decane volatilization occurred near the inlet, the relative decane saturation throughout the core dropped uniformly, and decane in the core flowed in the liquid phase to the shell to replenish volatilized decane. These results indicate that NAPL trapped in low-permeability zones can flow to replenish areas where NAPL is lost due to SVE. However, when residual NAPL saturation is reached, NAPL flow no longer occurs and diffusion limits removal from low-permeability zones.

To determine the effects of concentration history on slow sorption processes, desorption kinetic profiles for trichloroethene (TCE) were measured for a soil at 100% relative humidity subject to different exposure concentrations and exposure times. Exposure concentrations ranged from 1% to 80% of the saturation vapor pressure ( P s ) for TCE, and exposure times (i.e., time allowed for sorption before desorption begins) ranged from 1 to 96 days. A spherical diffusion model based on a γ distribution of sorption rates and a γ distribution of desorption rates was developed and applied to the data. At 80% P / P s , the entire γ distributions of sorption and desorption rates were available for TCE. In accordance with a micropore filling mechanism, the fraction of these distributions available for TCE sorption decreased with decreasing P / P s . Experimental results are consistent with a micropore-filling mechanism, where the amount of slow desorbing mass decreased with decreasing exposure time, and the fraction of slow desorbing sites filled decreased with decreasing exposure concentration. Simulation results suggest that diffusion limits the rates that micropores fill, and that rates of sorption and desorption for soil contaminated at smaller values of P / P s are, on average, less than those at larger values of P / P s (i.e., slow desorption rates are a function of exposure concentration). Simulation results also suggest that the model adequately describes the effects of exposure concentration and exposure time on the rates of sorption and desorption, but not on the capacity of the slow sites for TCE. This work is important because contaminant concentrations in the subsurface vary in space and time, and the proposed model represents a new and mechanistically based approach to capture the effects of this heterogeneity on slow desorption.