Search

Ball milling has emerged as a promising destructive technique for treating per- and polyfluoroalkyl substances (PFAS)-impacted soils. Environmental media properties such as reactive species generated upon ball milling and particle size are postulated to influence the effectiveness of the technology. In this study, four media types amended with perfluorooctanoic acid (PFOA) and perfluorooctane sulfonate (PFOS) were planetary ball milled to investigate destruction, fluoride recovery without additional co-milling reagents and the relationship between PFOA and PFOS destruction, particle size during milling, and electron generation. Silica sand, nepheline syenite sand, calcite and marble were sieved to achieve similar initial particle sizes (6/35 distribution), amended with PFOA and PFOS, and milled for 4 h. Particle size analysis was conducted throughout milling and 2,2-diphenyl-1-picrylhydrazyl (DPPH•) was used as a radical scavenger to assess electron generation from the four media types. Particle size reduction was observed to be positively correlated to PFOA and PFOS destruction and DPPH• neutralization (demonstrating electron generation by milling) in silica sand and nepheline syenite sand. Milling of a fine fraction (< 500 μm) of silica sand revealed less destruction compared to the 6/35 distribution suggesting the ability to fracture grains in silicate media is integral to PFOA and PFOS destruction. DPPH• neutralization was demonstrated in all four amended media types, confirming silicate sands and calcium carbonates generate electrons as a reactive species during ball milling. Fluoride loss as a function of milling time was observed in all amended media types. A sodium fluoride (NaF) spiked was used to quantify fluoride loss in the media independent of PFAS. A method was developed using the NaF-amended media fluoride concentrations to estimate the total fluorine liberated from PFOA and PFOS by ball milling. Estimates produced suggest complete recovery of theoretical fluorine yield is obtained. Data from this study was used to propose a reductive destruction mechanism for PFOA and PFOS., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2023 Elsevier B.V. All rights reserved.)

Thermal conductive heating (TCH) is an in-situ thermal treatment (ISTT) technology for treating non-aqueous phase liquid (NAPL) source zones. Numerical models can be useful tools for improving remedial performance, but traditional multiphase flow models are rarely used to simulate mass recovery during ISTT applications at the field scale due to their computational expense. This study developed a 3D model based on macroscopic invasion percolation to simulate the vaporization of NAPL, and the subsequent vapor migration and potential condensation at the field scale. The model was used to simulate the mass recovery of trichloroethene (TCE) from a NAPL source zone under seven scenarios of different heater placements, including three scenarios with an undersized target treatment zone (TTZ). Simulation results showed that TCH was effective in removing NAPL within the TTZ, but the treatment zone did not extend far from the perimeter heaters. In addition, during heating, NAPL condensation outside the TTZ due to the escaping vapor was observed in all scenarios. Overall, the resulting mass recovery was lower in the three scenarios with an undersized TTZ (91-95%) than in the other four scenarios (≈ 99%). Moreover, the locations of unrecovered/condensed NAPL could be inferred by monitoring mass recovery tailing at individual extraction wells., Competing Interests: Declaration of Competing Interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2023 Elsevier B.V. All rights reserved.)

There is a need for destructive technologies for per- and polyfluoroalkyl substances (PFAS) in soil. While planetary ball mill have been shown successful degradation of PFAS, there are issues surrounding scale up (maximum size is typically 0.5 L cylinders). While having lower energy outputs, horizontal ball mills, for which scale up is not a limiting factor, already exist at commercial/industrial sizes from the mining, metallurgic and agricultural industries, which could be re-purposed. This study evaluated the effectiveness of horizontal ball mills in degrading perfluorooctanesulfonate (PFOS), 6:2 fluorotelomer sulfonate (6:2 FTSA), and aqueous film forming foam (AFFF) spiked on nepheline syenite sand. Horizontal ball milling was also applied to two different soil types (sand dominant and clay dominant) collected from a firefighting training area (FFTA). Liquid chromatography tandem mass spectrometry was used to track 21 target PFAS throughout the milling process. High-resolution accurate mass spectrometry was also used to identify the presence and degradation of 19 non-target fluorotelomer substances, including 6:2 fluorotelomer sulfonamido betaine (FtSaB), 7:3 fluorotelomer betaine (FtB), and 6:2 fluorotelomer thioether amido sulfonate (FtTAoS). In the presence of potassium hydroxide (KOH), used as a co-milling reagent, PFOS, 6:2 FTSA, and the non-target fluorotelomer substances in the AFFF were found to undergo upwards of 81%, 97%, and 100% degradation, respectively. Despite the inherent added complexity associated with field soils, better PFAS degradation was observed on the FFTA soils over the spiked NSS, and more specifically, on the FFTA clay over the FFTA sand. These results held through scale-up, going from the 1 L to the 25 L cylinders. The results of this study support further scale-up in preparation for on-site pilot tests., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2022 Elsevier B.V. All rights reserved.)

The removal of dissolved volatile organic compounds (VOCs) from low-permeability lenses is important to limit back diffusion at sites impacted by dense non-aqueous phase liquids (DNAPLs). In situ thermal treatment (ISTT) technologies have the potential to treat DNAPL-impacted sites by enhancing diffusion from low-permeability lenses during heating. A series of two-dimensional laboratory tank experiments was conducted to investigate heating, gas formation, and trichloroethene (TCE) removal from a clay lens surrounded by sand. Results showed preferential heating of the clay and substantial TCE removal, with post-heating relative concentrations less than 0.06. The extent of TCE removal was not explained by only an increase in the aqueous TCE diffusion coefficient with increased temperature. Modelling estimates based on 1D diffusion from the lens showed that diffusion through both gas and water phases was required to match observations. Gas formation in the interior of the lens was also indicated by measured changes in bulk electrical conductivity of the clay during cool down, with gas saturations estimated to be greater than 0.21 at the end of heating. These estimates were larger than those needed to match the observed removal by diffusion, and suggest that connected gas pathways were created in the lens during heating, but that not all of the gas produced was part of those pathways. These results suggest that ISTT technologies may be effective in removing dissolved VOCs from thin clay lenses, and that gas formation within the clay should be considered when predicting the extent and rate of removal., (Copyright © 2021 Elsevier B.V. All rights reserved.)

Per- and polyfluoroalkyl substances (PFAS) are manmade, fluorinated organic chemicals which have been identified as persistent organic pollutants. PFAS have surface active properties that have made them suitable for applications in oil- and water-resistant products, as well as many firefighting foams. No on-site remediation strategies exist to treat PFAS impacted soils. Mechanochemical remediation of PFOS- and PFOA-amended sand via a planetary ball mill was studied. The effect of sand mass, KOH as a co-milling reagent, and water saturation on the degradation of PFOA and PFOS was evaluated. By 4 h of milling concentrations were reduced by up to 98% for PFOS-amended dry sand and 99% for PFOA-amended dry sand without the addition of a co-milling reagent. Water saturation was determined to be a significant hindrance on the mechanochemical destruction of PFOS and PFOA. A maximum of 89% of fluoride was recovered from PFOS-amended sand when KOH was used as a co-milling reagent. It is hypothesized that reactive particles generated from the fracture of sand grains react with PFAS molecules to initiate destruction, which can result in full defluorination. Milling experiments were also conducted on soils from a Canadian firefighting training area (FFTA), demonstrating that PFOS concentrations can be reduced by up to 96% in site soils. For the first time, ball milling for the remediation of PFAS in environmental media has been demonstrated using amended sand and legacy soils from a FFTA., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2020 Elsevier B.V. All rights reserved.)

In situ thermal treatment (ISTT) technologies can be used to remove mass from non-aqueous phase liquid (NAPL) source zones. Ensuring the vaporization of NAPL and the capture of vapors are crucial, and numerical models are useful for understanding the processes that affect performance to help improve design and operation. In this paper, a two-dimensional model that combines a continuum approach based on finite difference for heat transfer with a macroscopic invasion percolation (macro-IP) approach for gas migration was developed to simulate thermal conductive heating (TCH) applications at the field-scale. This approach simulates heat transport and gas migration, but is different than a traditional continuum multiphase approach. Mass recovery for 60 randomly generated realizations under three degrees of heterogeneity of the permeability field were simulated. The mass recovery curves had an overall similar shape for the various permeability fields. However, a wider range of completion times was observed for domains with a higher permeability variance. Results also showed that NAPL pools that were highly saturated, deep, and away from the heaters needed more heating time to be depleted, and that total NAPL mass was not a good indicator of completion time. The completion time was positively correlated with the maximum value of the mixed spatial moment of NAPL saturation about the heaters in the lateral and vertical direction, and the NAPL pool with the highest moment could increase the heating time by as much as 35%. This effect was most notable in simulations with a high permeability variance and suggests the potential to reduce heating time by locating the largest NAPL pools and placing TCH heaters accordingly., Competing Interests: Declaration of Competing Interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2020 Elsevier B.V. All rights reserved.)

To address the challenge of trichloroethene (TCE) remediation in low permeability zone, an inexpensive Cu-Ni bimetallic cathode was proposed in electrokinetic (EK) remediation system to couple electrokinetic migration with in situ electrochemical hydrodechlorination. Aqueous phase TCE was originally added into the anolyte so that breakthrough curves through the low permeability porous soil compartment could be obtained to better understand TCE migration driven by electroosmosis flow using different cathodes. The Cu-Ni cathode resulted in more TCE migration of 7.64 mg compared to that of 5.99 mg with Ni and 4.22 mg with mixed metal oxide (MMO) cathode, suggesting that the Cu-Ni cathode was capable of driving more TCE flux out of the contaminated soil. With the Cu-Ni cathode, 98.4% of TCE flux that reached the cathode was electrochemically reduced on the cathode, which was much higher than that with MMO cathode (77.9%) or Ni cathode (59.6%). TCE mass that was transported by electroosmosis flow increased from 2.04 to 6.68 mg when the voltage gradient increased from 1 to 4 V cm -1 , with the normalized energy consumption increasing from 0.06 to 0.16 kWh kg -1 per unit water movement, and from 0.54 to 2.55 kWh g -1 per unit TCE transport. For TCE that did reach the cathode compartment, > 98% degradation maintained at the Cu-Ni cathode with various voltage gradients. The coupled electrokinetic and electrochemical hydrodechlorination technology appears to be a promising strategy for the remediation of low permeability porous media., (Copyright © 2020 Elsevier Ltd. All rights reserved.)

Treatability tests can be carried out to assess the potential effectiveness of thermal treatment technologies under different site conditions and are important for specific technology selection and design. In order to reduce the costs for laboratory tests and expand the insights from previous treatability studies, a one-dimensional (1D) radial finite difference model was developed to simulate the removal of volatile organic compounds (VOCs) in laboratory thermal treatability tests. The processes considered in the model include heat conduction, co-boiling of single-component or multi-component NAPLs with water, and water boiling. An explicit approach is used to simulate the evolution of NAPL composition for multi-component NAPLs during heating. The developed model adopts only two fitting parameters and was calibrated and validated using previous laboratory experiments. In this paper, the developed model was first calibrated to three laboratory experiments using temperature measurements, which resulted in matches to the NAPL and gas saturations. After calibration, the model was able to predict the temperature, NAPL and gas saturations for the remaining seven experiments, including those with single and multi-component NAPLs, using the average value of each fitting parameter., (Copyright © 2019 Elsevier B.V. All rights reserved.)

In situ thermal treatment (ISTT) technologies have been applied at sites impacted by non-aqueous phase liquids (NAPLs). There is a need to establish expectations for the treatment of semi-volatile NAPLs, including those consisting primarily of polycyclic aromatic hydrocarbons (PAHs), and the potential benefits and limitations of partial NAPL removal. A series of laboratory experiments was conducted to investigate NAPL removal and soil concentrations during the heating of creosote-impacted sand, as well as aqueous concentrations during post-heating dissolution. The results showed co-boiling near the water boiling temperature due to the low volatility of most creosote components, with limited decreases in NAPL saturation (from 30% to 21% of the pore space). Decreases in soil concentration were more substantial than decreases in NAPL saturation (by a factor of 2-180), with greater removal for higher-volatility components at higher treatment temperatures. Results of the dissolution experiments showed mixed results, with decreases in the aqueous concentrations for 12 of 15 components, but increases in aqueous concentrations for phenanthrene, fluoranthene and pyrene after heating to 205 °C or 320 °C. Overall, the results illustrate the utility of bench-scale treatability tests in helping to establish ISTT goals and expectations., (Copyright © 2018 Elsevier B.V. All rights reserved.)

A series of intermediate-scale laboratory experiments were completed in a two-dimensional flow cell to investigate gas production and migration during the application of electrical resistance heating (ERH) for the removal of dense non-aqueous phase liquids (DNAPLs). Experiments consisted of heating water in homogeneous silica sand and heating 270 mL of trichloroethene (TCE) and chloroform (CF) DNAPL pools in heterogeneous silica sands, both under flowing groundwater conditions. Spatial and temporal distributions of temperature were measured using thermocouples and observations of gas production and migration were collected using front-face image capture throughout the experiments. Post-treatment soil samples were collected and analyzed to assess DNAPL removal. Results of experiments performed in homogeneous sand subject to different groundwater flow rates showed that high groundwater velocities can limit subsurface heating rates. In the DNAPL pool experiments, temperatures increased to achieve DNAPL-water co-boiling, creating estimated gas volumes of 131 and 114 L that originated from the TCE and CF pools, respectively. Produced gas migrated vertically, entered a coarse sand lens and subsequently migrated laterally beneath an overlying capillary barrier to outside the heated treatment zone where 31-56% of the original DNAPL condensed back into a DNAPL phase. These findings demonstrate that layered heterogeneity can potentially facilitate the transport of contaminants outside the treatment zone by mobilization and condensation of gas phases during ERH applications. This underscores the need for vapor phase recovery and/or control mechanisms below the water table during application of ERH in heterogeneous porous media during the co-boiling stage, which occurs prior to reaching the boiling point of water., (Copyright © 2015 Elsevier B.V. All rights reserved.)

In situ thermal treatment technologies, such as electrical resistance heating and thermal conductive heating, use subsurface temperature measurements in addition to the analysis of soil and groundwater samples to monitor remediation performance. One potential indication of non-aqueous phase liquid (NAPL) removal is an increase in temperature following observations of a co-boiling plateau, during which subsurface temperatures remain constant as NAPL and water co-boil. However, observed co-boiling temperatures can be affected by the composition of the NAPL and the proximity of the NAPL to the temperature measurement location. Results of laboratory heating experiments using single-component and multi-component NAPLs showed that local-scale temperature measurements can be mistakenly interpreted as an indication of the end of NAPL-water co-boiling, and that significant NAPL saturations (1% to 9%) remain despite observed increases in temperature. Furthermore, co-boiling of multi-component NAPL results in gradually increasing temperature, rather than a co-boiling plateau. Measurements of gas production can serve as a complementary metric for assessing NAPL removal by providing a larger-scale measurement integrated over multiple smaller-scale NAPL locations. Measurements of the composition of the NAPL condensate can provide ISTT operators with information regarding the progress of NAPL removal for multi-component sources., (Copyright © 2014 Elsevier B.V. All rights reserved.)

Semianalytical transient solutions have been developed to evaluate what level of fractured porous media (e.g., bedrock or clay) matrix cleanup must be achieved in order to achieve compliance of fracture pore water concentrations within a specified time at specified locations of interest. The developed mathematical solutions account for forward and backward diffusion in a fractured porous medium where the initial condition comprises a spatially uniform, nonzero matrix concentration throughout the domain. Illustrative simulations incorporating the properties of mudstone fractured bedrock demonstrate that the time required to reach a desired fracture pore water concentration is a function of the distance between the point of compliance and the upgradient face of the domain where clean groundwater is inflowing. Shorter distances correspond to reduced times required to reach compliance, implying that shorter treatment zones will respond more favorably to remediation than longer treatment zones in which back-diffusion dominates the fracture pore water response. For a specified matrix cleanup goal, compliance of fracture pore water concentrations will be reached sooner for decreased fracture spacing, increased fracture aperture, higher matrix fraction organic carbon, lower matrix porosity, shorter aqueous phase decay half-life, and a higher hydraulic gradient. The parameters dominating the response of the system can be measured using standard field and laboratory techniques., (© 2012, The Author(s). Ground Water © 2012, National Ground Water Association.)

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., (© 2012, The Author(s). Ground Water © 2012, National Ground Water Association.)

A two-dimensional experiment employing a heterogeneous sand pack incorporating two pools of trichloroethylene (TCE) was performed to assess the efficacy of electrical resistance heating (ERH) under passive venting conditions. Temperature monitoring displayed the existence of a TCE-water co-boiling plateau at 73.4°C, followed by continued heating to 100°C. A 5cm thick gas accumulation formed beneath a fine-grained capillary barrier during and after co-boiling. The capillary barrier did not desaturate during the course of the experiment; the only pathway for gas escape being through perforated wells traversing the barrier. The thickness of the accumulation was dictated by the entry pressure of the perforated well. The theoretical maximum TCE soil concentration within the region of gas accumulation, following gas collapse, was estimated to be 888mg/kg. Post-heating soil sampling revealed TCE concentrations in this region ranging from 27mg/kg to 96.7mg/kg, indicating removal of aqueous and gas phase TCE following co-boiling as a result of subsequent boiling of water. The equilibrium concentrations of TCE in water corresponding to the range of post-treatment concentrations in soil (6.11mg/kg to 136mg/kg) are calculated to range from 19.8mg/l to 440mg/l. The results of this experiment illustrate the importance of providing gas phase venting during the application of ERH in heterogeneous porous media., (Copyright © 2011 Elsevier B.V. All rights reserved.)

Numerical modeling was employed to study the performance of thermal conductive heating (TCH) in fractured shale under a variety of hydrogeological conditions. Model results show that groundwater flow in fractures does not significantly affect the minimum treatment zone temperature, except near the beginning of heating or when groundwater influx is high. However, fracture and rock matrix properties can significantly influence the time necessary to remove all liquid water (i.e., reach superheated steam conditions) in the treatment area. Low matrix permeability, high matrix porosity, and wide fracture spacing can contribute to boiling point elevation in the rock matrix. Consequently, knowledge of these properties is important for the estimation of treatment times. Because of the variability in boiling point throughout a fractured rock treatment zone and the absence of a well-defined constant temperature boiling plateau in the rock matrix, it may be difficult to monitor the progress of thermal treatment using temperature measurements alone., (Copyright © 2010 The Author(s). Journal compilation © 2010 National Ground Water Association.)

The influence of source zone concentration reduction on solute plume detachment and recession times in fractured rock was investigated using new semianalytical solutions to transient solute transport in the presence of advection, dispersion, sorption, matrix diffusion, and first-order decay. Novel aspects of these solutions are: (1) the source zone concentration behavior is simulated using a constant concentration with the option for either an instantaneous reduction to zero concentration or an exponentially decaying source zone concentration initiated at some time (t*) after the source is introduced, and (2) different biodegradation rates in the fracture and rock matrix. These solutions were applied for sandstone bedrock and revealed that biodegradation in the matrix, not the fracture, may be the most significant attenuation mechanism and therefore may dictate remediation time scales. Also, instantaneous and complete source concentration reduction in aged plumes may not be beneficial with respect to plume response because back-diffusion can sustain plume migration for long periods of time. Moderate source zone concentration reduction has a similar impact on the rate of advance of the leading edge of the plume as aggressive concentration reduction. If the source zone concentration reduction half-life is less than the plume decay half-life, then volatile organic compound (VOC) mass sequestered in the rock matrix will ultimately dictate plume persistence and not the presence of the source zone.

A two-dimensional, transient-flow, and transport numerical model was developed to simulate in situ chemical oxidation (ISCO) of trichloroethylene and tetrachloroethylene by potassium permanganate in fractured clay. This computer model incorporates dense, nonaqueous phase liquid dissolution, reactive aquifer material, multispecies matrix diffusion, and kinetic formulations for the oxidation reactions. A sensitivity analysis for two types of parameters, hydrogeological and engineering, including matrix porosity, matrix organic carbon, fracture aperture, potassium permanganate dosage, and hydraulic gradient, was conducted. Remediation metrics investigated were the relative rebound concentrations arising from back diffusion and percent mass destroyed. No well-defined correlation was found between the magnitude of rebound concentrations during postremedy monitoring and the amount of contaminant mass destroyed during the application. Results indicate that all investigated parameters affect ISCO remediation in some form. Results indicate that when advective transport through the fracture is dominant relative to diffusive transport into the clay matrix (large System Peclet Number), permanganate is more likely to be flushed out of the system and treatment is not optimal. If the System Peclet Number is too small, indicating that diffusion into the matrix is dominant relative to advection through the fracture, permanganate does not traverse the entire fracture, leading to postremediation concentration rebound. Optimal application of ISCO requires balancing advective transport through the fracture with diffusive transport into the clay matrix.

A fixed-volume release of 1,2-DCE, tracked in space and time with a light transmission/image analysis system, provided a data set for the infiltration, redistribution, and immobilisation of a dense non-aqueous phase liquid (DNAPL) in a heterogeneous porous medium. The two-dimensional bench scale flow cell was packed with a spatially correlated, random heterogeneous distribution of six sand types. In order to provide the necessary modelling parameters, detailed constitutive relationships were measured at the local scale for the six sands. These experiments revealed that nonwetting phase (NWP) relative permeability-saturation (k(rN)-S(W)) relationships are strongly correlated to sand type. Trends in the best-fit k(rN)-S(W) parameters reflected a positive correlation between mean grain diameter and the maximum NWP relative permeability, k(rN)(max). Multiphase flow simulations of the bench scale experiment best reproduced the experimental observations, producing excellent matches in both time and space, when the measured, correlated local scale k(rN)-S(W) relationships were employed.

A mathematical solution for solute transport in a three-dimensional porous medium with a patch source under steady-state, uniform ground water flow conditions was developed by Domenico (1987). The solution derivation strategy used an approximate approach to solve the boundary value problem, resulting in a nonexact solution. Variations of the Domenico (1987) solution are incorporated into the software programs BIOSCREEN and BIOCHLOR, which are frequently used to evaluate subsurface contaminant transport problems. This article mathematically elucidates the error in the approximation and presents simulations that compare different versions of the Domenico (1987) solution to an exact analytical solution to demonstrate the potential error inherent in the approximate expressions. Results suggest that the accuracy of the approximate solutions is highly variable and dependent on the selection of input parameters. For solute transport in a medium-grained sand aquifer, the Domenico (1987) solution underpredicts solute concentrations along the centerline of the plume by as much as 80% depending on the case of interest. Increasing the dispersivity, time, or dimensionality of the system leads to increased error. Because more accurate exact analytical solutions exist, we suggest that the Domenico (1987) solution, and its predecessor and successor approximate solutions, need not be employed as the basis for screening tools at contaminated sites.

The time required for dense nonaqueous phase liquid (DNAPL) to cease migrating following release to the subsurface is a valuable component of a site conceptual model. This study uses numerical simulation to investigate the migration of six different DNAPLs in sandy aquifers. The most influential parameters governing migration cessation time are the density and viscosity of the DNAPL and the mean hydraulic conductivity of the aquifer. Releases of between 1 and 40 drums of chlorinated solvent DNAPLs, characterized by relatively high density and low viscosity, require on the order of months to a few years to cease migrating in a heterogeneous medium sand aquifer having an average hydraulic conductivity of 7.4 x 10(-3) cm/s. In contrast to this, the release of 20 drums of coal tar (rho(D)= 1061 kg/m(3), micro(D)= 0.161 Pa.s) requires more than 100 years to cease migrating in the same aquifer. Altering the mean hydraulic conductivity of the aquifer results in a proportional change in cessation times. Parameters that exhibit relatively little influence on migration time scales are the DNAPL-water interfacial tension, release volume, source capillary pressure, mean aquifer porosity, and ambient ground water hydraulic gradient. This study also demonstrates that low-density DNAPLs (e.g., coal tar) give rise to greater amounts of lateral spreading and greater amounts of pooling on capillary barriers than high-density DNAPLs such as trichloroethylene or tetrachloroethylene.

Numerical simulation is used to examine the relative velocities of DNAPL and aqueous phase plumes in sandy aquifers where lateral spreading of DNAPL has occurred at the base of the aquifer. The scenario being modeled is one where a permeable aquifer is underlain by a sloping aquitard, which results in lateral migration of the DNAPL down the slope, in addition to lateral migration of an aqueous phase plume subject to a specified hydraulic gradient. A sensitivity analysis is presented to the impacts of both DNAPL properties and geologic properties. The most important chemical properties governing the relative velocities of the DNAPL and the shallow aqueous phase plume are the DNAPL viscosity and the aqueous component soil-water partition coefficient (Kd). The dip of the underlying aquitard was found to be relatively unimportant, at least for the range of values studied. The scenario under consideration can be important in conceptual model development and remedial design, as in certain cases DNAPL could be migrating in areas without the evidence of a well-developed aqueous phase plume. The implication of this work is that the absence of a shallow aqueous phase plume directly downgradient of a DNAPL source zone does not rule out the possibility of deep occurrences of DNAPL beyond the shallow monitoring well network. A further finding of this study is that the occurrence of a highly sorbing compound in groundwater at virtually any concentration may indicate the immediate upgradient presence of residual or pooled DNAPL.

Phase behaviour experiments employing PCB (Aroclor 1242)/alcohol/water systems were conducted with ethanol (EtOH) and n-propanol (nPA). Both exhibited an affinity for the aqueous phase within the entire two-phase region. As much as 88% by volume (88% vol.) EtOH and 80% vol. nPA were necessary to achieve full miscibility of the PCB in the aqueous phase. DNAPL-water interfacial tension (IFT) was reduced from 38.9 dyn/cm to 4.7 dyn/cm and 2.4 dyn/cm with 80% vol. EtOH and 76% vol. nPA. The addition of alcohol brought about 41% and 54% reductions in DNAPL viscosity at maximal concentrations of EtOH and nPA. Density of the PCB-DNAPL was relatively unaffected by the presence of alcohol. A series of seven experiments were conducted where successive slugs of nPA and xanthan gum polymer solutions were injected into a fractured shale sample. A 30% vol. nPA solution injected under a hydraulic gradient of 0.36 allowed enhanced PCB removal primarily through reduction of IFT and resulted in 72% DNAPL recovery. Several pore volumes of alcohol solution were necessary to displace all the potentially mobile non-wetting phase since the high-viscosity DNAPL was mobilized at a lower flow rate than the overall fluid velocity, illustrating non-piston displacement. The injection of a 95% vol. nPA alcohol solution, theoretically at a sufficient concentration to produce fully miscible displacement of the residual DNAPL at equilibrium, resulted in non-equilibrium partitioning of the PCB into the flushing solution, likely due to the high fluid velocities in the fracture. The injection of 200 pore volumes of 95% vol. nPA solution resulted in 94% DNAPL recovery. Alcohol floods operated below the miscibility envelope appear to be a valuable source zone remedial alternative where the objective is to reduce DNAPL mobility to zero, but it should be noted that DNAPL mobility is increased during the application of the technology and steps may need to be taken to prevent unwanted vertical mobilization.

An emulsion stabilized by colloidal suspensions of finely divided solids is known as a Pickering emulsion. The potential for polychlorinated biphenyls (PCBs) to form Pickering emulsions ex situ when in contact with powdered solids, such as clays and metal oxides, is investigated here. Bentonite, iron oxide and magnesium oxide dispersions proved to be robust Pickering emulsion stabilizers, whereas manganese oxide dispersions were not. Batch experiments revealed that emulsions can be formed using a moderately low energy input and can be stabilized with solid concentrations as low as 0.5 wt.%. For the base conditions (volumetric oil fraction (phi(oil))=30 vol.%; solid concentration (chi)=2 wt.%), the formed emulsions were indefinitely stable and the initial average droplet diameters varied from 80 to 258 mum, depending on the solid used in the colloidal dispersion. The average droplet size varied at early time, but for most conditions stabilized to a steady-state value 1 week after preparation. The effect of Ostwald ripening was limited. At greater than 0.5 wt.% concentration, the efficiency of the solid dispersion as a stabilizer was dependant on the volumetric oil fraction but not on the solid concentration. Generally, systems with volumetric oil fractions outside of the 20-70 vol.% range were unstable. The emulsions' droplet stability, average droplet size and size distribution were observed to vary as a function of the amount of energy provided to the system, the volumetric oil fraction, and the concentration of the solid in the aqueous dispersion. It is hypothesized that drilling through fractured rock in the immediate vicinity of dense, non-aqueous phase liquid (DNAPL) PCBs may provide both the energy and solid material necessary to form Pickering emulsions.

Matrix diffusion can attenuate the rate of plume migration in fractured bedrock relative to the rate of ground water flow for both conservative and nonconservative solutes of interest. In a system of parallel, equally spaced constant aperture fractures subject to steady-state ground water flow and an infinite source width, the degree of plume attenuation increases with time and travel distance, eventually reaching an asymptotic level. The asymptotic degree of plume attenuation in the absence of degradation can be predicted by a plume attenuation factor, beta, which is readily estimated as R' (phi(m)/phi(f)), where R' is the retardation factor in the matrix, phi(m) is the matrix porosity, and phi(f) is the fracture porosity. This dual-porosity relationship can also be thought of as the ratio of primary to secondary porosity. Beta represents the rate of ground water flow in fractures relative to the rate of plume advance. For the conditions examined in this study, beta increases with greater matrix porosity, greater matrix fraction organic carbon, larger fracture spacing, and smaller fracture aperture. These concepts are illustrated using a case study where dense nonaqueous phase liquid in fractured sandstone produced a dissolved-phase trichloroethylene (TCE) plume approximately 300 m in length. Transport parameters such as matrix porosity, fracture porosity, hydraulic gradient, and the matrix retardation factor were characterized at the site through field investigations. In the fractured sandstone bedrock examined in this study, the asymptotic plume attenuation factors (beta values) for conservative and nonconservative solutes (i.e., chloride and TCE) were predicted to be approximately 800 and 12,210, respectively. Quantitative analyses demonstrate that a porous media (single-porosity) solute transport model is not appropriate for simulating contaminant transport in fractured sandstone where matrix diffusion occurs. Rather, simulations need to be conducted with either a discrete fracture model that explicitly incorporates matrix diffusion, or a dual-continuum model that accounts for mass transfer between mobile and immobile zones. Simulations also demonstrate that back diffusion from the matrix to fractures will likely be the time-limiting factor in reaching ground water cleanup goals in some fractured bedrock environments.

Flow-through column tests were conducted to investigate the performance of iron wall remediation systems for the degradation of aqueous-phase trichloroethylene (TCE). Concentration profiles under steady-state transport conditions were generated by measuring TCE concentrations at sample ports located at various locations along the length of the column. The results indicated that a pseudo-first-order model is adequate at describing degradation kinetics for low initial TCE concentrations, but not for higher initial concentrations. The deviation from pseudo-first-order kinetics can be explained by interspecies competition for reactive sites between TCE and a dominant reaction product. A modification of the pseudo-first-order model that accounts for product interference predicts laboratory data for high initial concentration profiles, but deviates slightly as initial concentrations approach the solubility of TCE. The data clearly demonstrate the importance of accurately describing reaction kinetics for the purpose of designing iron wall treatment systems.

The migration of Dense, Non-Aqueous Phase Liquid (DNAPL) and dissolved phase contamination through a fractured heterogeneous porous medium has been investigated through the use of a multiphase compositional model. The sensitivity of the timescales of migration and the distribution of contaminant in the subsurface to the mean permeability, the variance of the permeability, and the degree of fracturing of the domain were examined. It was found that increasing the mean permeability of the domain allowed the DNAPL to penetrate deeper into the subsurface, while decreasing the mean permeability caused the DNAPL to pool at shallower depths. The presence of fractures within the system was found to control the infiltration only in the most fractured domain. Moment analysis of the nonwetting phase showed that large-scale movement had ceased after approximately 9 years (maximum duration of the source-on condition was approximately 4.5 years). This tended to be due to a redistribution of the DNAPL towards a residual configuration, as was evidenced by the gradual trending of average nonwetting phase saturations within the domain to a static value. The dissolved phase plume was found to migrate at essentially the same rate as the nonwetting phase, due to the reduced relative permeability of lenses containing DNAPL, and due to diffusive losses of mass to the matrix of fractured clay and silty-clay lenses. Some exceptions to this were found when the DNAPL could not overcome the displacement pressure of a lens, and could not by-pass the lens due to the lack of available driving force after the source had been shut off., (Copyright 2003 Elsevier B.V.)

The migration of five dense nonaqueous phase liquids (DNAPLs) through a single fracture in a clay aquitard was numerically simulated with the use of a compositional simulator. The effects of fracture aperture, fracture dip, matrix porosity, and matrix organic carbon content on the migration of chlorobenzene, 1,2-dichloroethylene, trichloroethylene, tetra-chloroethylene, and 1,2-dibromoethane were examined. Boundary conditions were chosen such that DNAPL entry into the system was allowed to vary according to the stresses applied. The aperture is the most important factor of those studied controlling the migration rate of DNAPL through a single fracture embedded in a clay matrix. Loss of mass to the matrix through diffusion does not significantly retard the migration rate of the DNAPL, particularly in larger aperture fractures (e.g., 50 microm). With time, the ratio of diffusive loss to the matrix to DNAPL flux into the fracture approaches an asymptotic value lower than unity. The implication is that matrix diffusion cannot arrest the migration of DNAPL in a single fracture. The complex relationships between density, viscosity, and solubility that, to some extent, govern the migration of DNAPL through these systems prevent accurate predictions without the use of numerical models. The contamination potential of the migrating DNAPL is significantly increased through the transfer of mass to the matrix. The occurrence of opposite concentration gradients within the matrix can cause dissolved phase contamination to exist in the system for more than 1000 years after the DNAPL has been completely removed from the fracture.

A numerical model (Queen's University Multi-Phase Flow Simulator, QUMPFS) was used to assess the rate of trichloroethylene (TCE) dense, non-aqueous phase liquid (DNAPL) migration through fractured clay, with special attention focused on the influence of interbedded sand lenses. The presence of these sand lenses was found to increase the time required for the non-wetting phase to migrate through the full 30 m vertical extent of the clay sequence from a few days to several years. Applied vertical hydraulic gradients were found to be moderately influential in systems consisting solely of fractured clays, yet one of the dominant factors controlling speed of vertical migration when sand lenses were present. Larger displacement pressure of the sands relative to that of the fractures leads to slower DNAPL migration rates, due to the delays that occur during build-up of capillary pressures. Dissolution of DNAPL and subsequent matrix diffusion of the aqueous phase has little effect on the rate of DNAPL migration through systems consisting of fractured clay only, yet slows the rate of migration in systems containing sand lenses. In all cases examined, the rate of DNAPL loading to the lower aquifer far exceeded the rate of aqueous phase mass loading. It was also found that DNAPL reaches the lower aquifer at approximately the same time as the aqueous phase plumes even for systems experiencing downward groundwater flow due to the attenuation of the aqueous phase through matrix diffusion.