Linking a one-dimensional pesticide fate model to a three-dimensional groundwater model to simulate pollution risks of shallow and deep groundwater underlying fractured till

The one-dimensional pesticide fate model MACRO was loose-linked to the three-dimensional discrete fracture/matrix diffusion model FRAC3DVS to describe transport of the pesticide mecoprop in a fractured moraine till and local sand aquifer (5-5.5 m depth) overlying a regional limestone aquifer (16 m depth) at Havdrup, Denmark. Alternative approaches to describe the upper boundary in the groundwater model were examined. Field-scale simulations were run to compare a uniform upper boundary condition with a spatially variable upper boundary derived from Monte-Carlo simulations with MACRO. Plot-scale simulations were run to investigate the influence of the temporal resolution of the upper boundary conditions for fluxes in the groundwater model and the effects of different assumptions concerning the macropore/fracture connectivity between the two models. The influence of within-field variability of leaching on simulated mecoprop concentrations in the local aquifer was relatively small. A fully transient simulation with FRAC3DVS gave 20 times larger leaching to the regional aquifer compared to the case with steady-state water flow, assuming full connectivity with respect to macropores/fractures across the boundary between the two models. For fully transient simulations 'disconnecting' the macropores/fractures at the interface between the two models reduced leaching by a factor 24. A fully connected, transient simulation with FRAC3DVS, with spatially uniform upper boundary fluxes derived from a MACRO simulation with 'effective' parameters is therefore recommended for assessing leaching risks to the regional aquifer, at this, and similar sites.     

Multiphase flow and transport in fractured clay/sand sequences.

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.     

Evaluation of TCDD biodegradability under different redox conditions

Polychlorinated dibenzo-p-dioxins have been generated as unwanted by-products in many industrial processes. Although their widespread distribution in different environmental compartments has been recognized, little is known about their fate in the ultimate environment sinks. The highly stable dioxin isomer 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) has been called the most toxic compound known to man. In this laboratory microcosm study, TCDD bioavailability was evaluated under five reduction/oxidation (redox) conditions including aerobic biodegradation, aerobic cometabolism, methanogenesis, iron reduction, and reductive dechlorination. Activated sludge and aquifer sediments from a TCDD and a pentachlorophenol (PCP) contaminated site were used as the inocula. Acetate, sludge cake, and cane molasses were used as the primary substrates (carbon sources) in cometabolism and reductive dechlorination microcosms. After a 90-day incubation period, microcosms constructed under reductive dechlorination conditions were the only treatment showing promising remediation results. The highest TCDD degradation rate [up to 86% of TCDD removal (with an initial concentration of 96 microg/kg of soil)] was observed in the microcosms with anaerobic activated sludge as the microbial inocula and sludge cakes as the primary substrates. Except for reductive dechlorination microcosms, no significant TCDD removal was observed in the microcosms prepared under other conditions. Thus, application of an effective primary substrate to enhance the reductive dechlorination process is a feasible method for TCDD bioremediation. Bioremediation expense can be significantly reduced by the supplement of some less expensive alternative substrates (e.g., sludge cakes, cane molasses). Results would be useful in designing a scale-up in situ or on-site bioremediation system such as bioslurry reactor for field application.     

Two-phase flow and transport in a single fracture-porous medium system

The prognosis for the remediation of contaminated fractured media is much worse than that for more homogeneous units. Fractures act as conduits for the flow of dense non-aqueous phase liquids (DNAPLs), while diffusion is responsible for the storage of dissolved mass in the surrounding matrix. A numerical model incorporating aqueous phase transport in a variable-aperture fracture and its surrounding matrix is developed and coupled with an existing two-phase flow model. The processes of transient two-phase flow, non-equilibrium dissolution, advective–dispersive transport in the fracture, and three-dimensional matrix diffusion are included in the model. Results from various investigations show that the DNAPL distribution is very sensitive to variations in aperture within a single fracture. Diffusion-controlled mass removal from both the matrix and from the hydraulically inaccessible zones within the fracture itself result in extremely large time frames for significant mass removal from these systems. Success in aqueous phase mass removal from the matrix is very sensitive to the effective fracture spacing. The hydraulic gradient in the fracture only affects the amount of water removed from the system, and does not greatly affect the amount of time required to remove the contaminant mass from the matrix. The ability to remove mass is somewhat sensitive to the porosity and effective matrix diffusion coefficient over the range of expected values.     

In situ bioelectrokinetic remediation of phenol-contaminated soil by use of an electrode matrix and a rotational operation mode

In situ bioremediation is a safe and cost-effective technology for the cleanup of contaminated sites, but its remediation rate is usually very slow. This study attempted to accelerate the process of bioremediation by employing non-uniform electrokinetic transport processes to mix organic pollutants and degrading bacteria in soils under in situ conditions (namely, in situ bioelectrokinetic remediation) by use of an electrode matrix and a rotational operation mode. A bench-scale non-uniform electrokinetic system with periodic polarity-reversal was developed for this purpose, and tested by using a sandy loam spiked with phenol as a model organic pollutant. The results demonstrated that non-uniform electrokinetic processes could enhance the in situ biodegradation of phenol in the soil, the efficiency of which depended upon the operational mode of the electric field. Compared with the unidirectional operation and the bidirectional operation, the rotational operation could effectively stimulate the biodegradation of phenol in the soil if adopting appropriate time intervals of polarity-reversal and electrode matrixes. A reversal interval of 3.0 h and a square-shaped electrode matrix with four electrode couples appeared appropriate for the in situ biodegradation of phenol, at which a maximum phenol removal of 58% was achieved in 10d and the bioremediation rate was increased about five times as compared to that with no electric field applied. The results also showed that adopting a small polarity-reversal interval and an appropriate electrode array could produce a high and uniform removal of phenol from the soil. It is believed that in situ bioelectrokinetic remediation holds the potential for field application.     

Numerical analysis of water and solute transport in variably-saturated fractured clayey till

This study numerically investigates the influence of initial water content and rain intensities on the preferential migration of two fluorescent tracers, Acid Yellow 7 (AY7) and Sulforhodamine B (SB), through variably-saturated fractured clayey till. The simulations are based on the numerical model HydroGeoSphere, which solves 3D variably-saturated flow and solute transport in discretely-fractured porous media. Using detailed knowledge of the matrix, fracture, and biopore properties, the numerical model is calibrated and validated against experimental high-resolution tracer images/data collected under dry and wet soil conditions and for three different rain events. The model could reproduce reasonably well the observed preferential migration of AY7 and SB through the fractured till, although it did not capture the exact depth of migration and the negligible impact of the dead-end biopores in a near-saturated matrix. A sensitivity analysis suggests fast flow mechanisms and dynamic surface coating in the biopores, and the presence of a plough pan in the till.     

Migration of colloids in discretely fractured porous media: effect of colloidal matrix diffusion

Numerical simulations of colloid transport in discretely fractured porous media were performed to investigate the importance of matrix diffusion of colloids as well as the filtration and remobilization of colloidal particles in both the fractures and porous matrix. To achieve this objective a finite element numerical code entitled COLDIFF was developed. The processes that COLDIFF takes into account include advective-dispersive transport of colloids, filtration and remobilization of colloidal particles in both fractures and porous matrix, and diffusive interactions of colloids between the fractures and porous matrix. Three sets of simulations were conducted to examine the importance of parameters and processes controlling colloid migration. First, a sensitivity analysis was performed using a porous block containing a single fracture to determine the relative importance of various phenomenological coefficients on colloid transport. The primary result of the analysis showed that the porosity of the matrix and the process of colloid filtration in fractures play important roles in controlling colloid migration. Second, simulations were performed to replicate and examine the results of a laboratory column study using a fractured shale saprolite. Results of this analysis showed that the filtration of colloidal particles in the porous matrix can greatly affect the tailing of colloid concentrations after the colloid source was removed. Finally, field-scale simulations were performed to examine the effect of matrix porosity, fracture filtration and fracture remobilization on long-term colloid concentration and migration distance. The field scale simulations indicated that matrix diffusion and fracture filtration can significantly reduce colloid migration distance. Results of all three analyses indicated that in environments where porosity is relatively high and colloidal particles are small enough to diffuse out of fractures, the characteristics of the porous matrix that affect colloid transport become more important than those of the fracture network. Because the properties of the fracture network tend to have greater uncertainty due to difficulties in their measurement relative to those of the porous matrix, prediction uncertainties associated with colloid transport in discretely fractured porous media may be reduced.     

Experimental investigation of pneumatic soil vapor extraction

Soil Vapor Extraction (SVE) is a common remediation technique for removing volatile organic compounds from unsaturated contaminated soils. Soil heterogeneities can however cause serious limitations to the applicability of SVE due to air bypassing low permeable areas of the soil, leading to diffusion limitation of the remediation. To enhance removal from areas subject to diffusion limitation a new remediation technique, pneumatic soil vapor extraction, is proposed. In contrast to traditional SVE, in which soil vapor is extracted continuously by a vacuum pump, pneumatic SVE is based on enforcing a sequence of large pressure drops on the system to enhance the recovery from the low-permeable areas. The pneumatic SVE technique was investigated in the laboratory using TCE as a model contaminant. 2D-laboratory tank experiments were performed on homogeneous and heterogeneous sand packs. The heterogeneous packs consisted of a fine sand lens surrounded by a coarser sand matrix. As expected when using traditional SVE, the removal of TCE from the low permeable lens was extremely slow and subject to diffusion limitation. In contrast when pneumatic venting was used removal rates increased by up to 77%. The enhanced removal was hypothesized to be attributed to mixing of the contaminated air inside the lens and generation of net advective transport out of the lens due to air expansion.     

Comparative study of microbial dechlorination of chlorinated ethenes in an aquifer and a clayey aquitard

In order to determine whether natural attenuation of chlorinated ethenes by microbial activity occurs in aquitards, sediments at a site contaminated with tetrachloroethene were vertically studied by drilling. The distribution of microbes (Dehalococcoides group and anaerobic hydrogen producers) and the ability of the sediments to sustain microbial dechlorination were determined in an aquitard as well as in an aquifer. Close-spaced sampling revealed the existence of large populations of Dehalococcoides and H(2)-producing bacteria, especially in the organic-rich clayey aquitard rather than in the aquifer. The vinyl chloride reductase gene was also detected in the clay layer. Furthermore, incubation experiments indicated that the clay sediment could sustain transformations of tetrachloroethene at least to vinyl chloride. In contrast, no significant transformation was observed in the aquifer sand. Our results indicate that dechlorination of tetrachloroethene by bacteria can take place in an organic-rich clayey aquitard, and that organic-rich clay may also be important in the natural attenuation in an adjacent aquifer, possibly supplying a carbon source or an electron donor.     

Numerical modeling of coupled variably saturated fluid flow and reactive transport with fast and slow chemical reactions

The couplings among chemical reaction rates, advective and diffusive transport in fractured media or soils, and changes in hydraulic properties due to precipitation and dissolution within fractures and in rock matrix are important for both nuclear waste disposal and remediation of contaminated sites. This paper describes the development and application of LEHGC2.0, a mechanistically based numerical model for simulation of coupled fluid flow and reactive chemical transport, including both fast and slow reactions in variably saturated media. Theoretical bases and numerical implementations are summarized, and two example problems are demonstrated. The first example deals with the effect of precipitation/dissolution on fluid flow and matrix diffusion in a two-dimensional fractured media. Because of the precipitation and decreased diffusion of solute from the fracture into the matrix, retardation in the fractured medium is not as large as the case wherein interactions between chemical reactions and transport are not considered. The second example focuses on a complicated but realistic advective-dispersive-reactive transport problem. This example exemplifies the need for innovative numerical algorithms to solve problems involving stiff geochemical reactions.     

The need for bioaugmentation after thermal treatment of a TCE-contaminated aquifer: Laboratory experiments

A microcosm study was conducted to evaluate the need for bioaugmentation after a thermal treatment to anaerobically dechlorinate trichloroethene (TCE) to ethene. The microcosms were either: heated to 100 degrees C and slowly cooled to simulate thermal remediation while bioaugmenting when the declining temperature reached 10 degrees C; or kept at ambient groundwater temperatures (10 degrees C) and bioaugmented for comparison. Aquifer samples from three sediment locations within a TCE-polluted source zone were investigated in duplicate microcosms. In biostimulated (5 mM lactate) and heated microcosms, no conversion of TCE was observed in 4 out of 6 microcosms, and in the remaining microcosms the dechlorination of TCE was incomplete to cDCE (cis-dichloroethene). By comparison, complete TCE dechlorination to ethene was observed in 4 out of 6 heated microcosms that were bioaugmented with a highly enriched dechlorinating mixed culture, KB-1, but no electron donor, and also in 4 of 6 microcosms that were augmented with KB-1 and an electron donor (5 mM lactate). These data suggest that electron donor released during heating, was capable of promoting complete dechlorination coincident with bioaugmentation. Heated microcosms demonstrated less methanogenesis than unheated microcosms, even with elevated H2 concentrations and addition of KB-1, which contains methanogens. This suggests that the heating process suppressed the native microbial community, which can decrease competition with the bioaugmented culture and increase the effectiveness of dechlorination following a thermal treatment. Specifically, cDCE removal rates were four to six times higher in heated than unheated bioaugmented microcosms. This study confirms the need for bioaugmentation following a laboratory thermal treatment to obtain complete dechlorination of TCE.     

Dechlorination of chloroorganics,decolorization, and simultaneous bioremediation of Cr6+ from real tannery effluent employing indigenous Bacillus cereus isolate

A native Bacillus cereus isolate has been employed, for the first time, for simultaneous decolorization, dechlorination of chloroorganics, and Cr⁶⁺ remediation from the real tannery effluent. Most of the physicochemical variables in 3:1 diluted effluent were well above the standard prescribed limits, which decreased substantially upon microbial treatment. The extent of bioremediation was better in diluted (3:1) as compared to undiluted effluent supplemented with nutrients and augmented with B. cereus isolate. Maximum growth, effluent decolorization (42.5 %), dechlorination (74.1 %), and Cr⁶⁺ remediation (34.2 %) were attained with 4.0 % (v/v) inoculum, 0.8 % glucose, and 0.2 % ammonium chloride in 3:1 diluted effluent at natural pH (8.1) within 72 h of incubation. The efficiency of bioremediation in a bioreactor was higher as compared to a flask trial during 72 h of incubation: decolorization (47.9 %) was enhanced by 5.4 %, dechlorination (77.4 %) by 3.3 %, and Cr⁶⁺ removal (41.7 %) by 7.5 % at an initial color of 286 Pt-Co units and initial concentration of 62 mg chloride ions and 108 mg l^?1 Cr⁶⁺. Immobilized biomass of Pseudomonas putida and B. cereus coculture enhanced the extent of Cr⁶⁺ remediation (51.9 %) by 10.2 % compared to the bioreactor trial. Chromate reductase activity and reduced Cr directly correlated and were mainly associated with soluble fraction of B. cereus plus effluent natural microflora. The GC-MS analyses revealed the formation of metabolites such as acetic acid and 2-butenoic acid in bacterially treated effluent. The supplementation of nutrients along with B. cereus augmentation is required for efficient effluent bioremediation.     

Coupling permanganate oxidation with microbial dechlorination of tetrachloroethene.

For sites contaminated with chloroethene non-aqueous-phase liquids, designing a remediation system that couples in situ chemical oxidation (ISCO) with potassium permanganate (KMnO4) and microbial dechlorination may be complicated because of the potentially adverse effects of ISCO on anaerobic bioremediation processes. Therefore, one-dimensional column studies were conducted to understand the effect of permanganate oxidation on tetrachloroethene (PCE) dechlorination by the anaerobic mixed culture KB-1. Following the confirmation of PCE dechlorination, KMnO4 was applied to all columns at a range of concentrations and application velocities to simulate varied distances from oxidant injection. Immediately following oxidation, reductive dechlorination was inhibited; however, after passing several pore volumes of sterile growth medium through the columns after oxidation, a rebound of PCE dechlorination activity was observed in every inoculated column without the need to reinoculate. The volume of medium required for a rebound of dechlorination activity differed from 1.1 to 8.1 pore volumes (at a groundwater velocity of 4 cm/d), depending on the specific condition of oxidant application.     

Modeling of non-reactive solute transport in fractured clayey till during variable flow rate and time

Fractures and biopores can act as preferential flow paths in clay aquitards and may rapidly transmit contaminants into underlying aquifers. Reliable numerical models for assessment of groundwater contamination from such aquitards are needed for planning, regulatory and remediation purposes. In this investigation, high resolution preferential water-saturated flow and bromide transport data were used to evaluate the suitability of equivalent porous medium (EPM), dual porosity (DP) and discrete fracture/matrix diffusion (DFMD) numerical modeling approaches for assessment of flow and non-reactive solute transport in clayey till. The experimental data were obtained from four large undisturbed soil columns (taken from 1.5 to 3.5 m depth) in which biopores and channels along fractures controlled 96-99% of water-saturated flow. Simulating the transport data with the EPM effective porosity model (FRACTRAN in EPM mode) was not successful because calibrated effective porosity for the same column had to be varied up to 1 order of magnitude in order to simulate solute breakthrough for the applied flow rates between 11 and 49 mm/day. Attempts to simulate the same data with the DP models CXTFIT and MODFLOW/MT3D were also unsuccessful because fitted values for dispersion, mobile zone porosity, and mass transfer coefficient between mobile and immobile zones varied several orders of magnitude for the different flow rates, and because dispersion values were furthermore not physically realistic. Only the DFMD modeling approach (FRACTRAN in DFMD mode) was capable to simulate the observed changes in solute transport behavior during alternating flow rate without changing values of calibrated fracture spacing and fracture aperture to represent the macropores.     

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.     

Enhancing polychlorinated biphenyl dechlorination in fresh water sediment with biostimulation and bioaugmentation

Polychlorinated biphenyls (PCBs) are toxic compounds ubiquitously distributed in ecosystems. Microbial attenuation of these contaminants is a potential means of remediation. Two promising microbial PCB remediation technologies, biostimulation and bioaugmentation, were investigated in different sediments. Biostimulation experiments in which electron donor was supplied (H2 via elemental iron, Fe(0)) resulted in only a marginal improvement in the dechlorination of amended 2,3,4,5-tetrachlorobiphenyl (2,3,4,5-CB), likely because of an inadequate population of indigenous H2-utilizing dechlorinators. Extensive dechlorination was observed, however, after bioaugmenting microcosms with a PCB-dechlorinating enrichment culture. Dechlorination of 2,3,4,5-CB began prior to the 20th day of incubation and proceeded to 2-chlorobiphenyl. This extensive dechlorination activity was maintained in both sediments over 70d at 10 and 25 degrees C. This research demonstrates that although past studies of biostimulation were promising, a great deal must be known about the PCB-dechlorinating organisms present before successful biostimulation is expected. Bioaugmentation, however, appears to be a promising PCB remediation strategy and should be further pursued.     

Plant aided bioremediation in the vadose zone: model development and applications

Phytoremediation has the potential to enhance clean up of land contaminated by various pollutants. A mathematical model that includes a two-fluid phase flow model of water flow as well as a two-region soil model of contaminant reactions was developed and applied to various bioremediation scenarios in the unsaturated zone, especially to plant-aided bioremediation. To investigate model behavior and determine the main parameters and mechanisms that affect bioremediation in unplanted and planted soils, numerical simulations of theoretical scenarios were conducted before applying the model to field data. It is observed from the results that parameters affecting the contaminant concentration in the water phase, such as aqueous solubility, the octanol-water partition coefficient, and organic carbon content of the soil controlled the contaminant fate in the vadose zone. Simulation using the developed model also characterized the fate and transport of the contaminants both in planted and unplanted soils satisfactorily for field applications. Although phytoremediation has the potential for remediation of contaminated soils, results from both modeling and field studies suggested that plants may not always enhance the remediation efficiency when the soil already has a high microbial concentration, when the contaminant bioavailability is low, or when the overall reaction is mass transfer-limited. Therefore, other steps to increase contaminant bioavailability are needed in phytoremediation applications; natural purification mechanisms such as aging, volatilization, and natural bioremediation should be considered to maximize the plant effect and minimize the cost.     

超声波净化石油污染土壤实验研究

为研究超声波技术处理石油污染土壤的影响因素,选取了土壤含油量、土水比、土壤类型、超声处理时间和超声功率5个因素进行了具有交互作用的正交实验,并讨论了各因素对土壤中石油污染物清除率的影响。实验结果表明超声波对土壤中石油污染物的净化效果明显,总石油烃(TPH)的清除率在23.0%到77.3%之间,土壤类型对净化效果的影响最为显著,粘土中的TPH清除率显著低于砂土。由超声空化引起的解吸作用和真空抽滤是实验过程中土壤TPH减少的主要途径。研究结果表明超声波技术用于石油污染土壤修复是可行并且有效果的。     

Miscible fluid displacement stability in unconfined porous media:: Two-dimensional flow experiments and simulations

In situ flushing groundwater remediation technologies, such as cosolvent flushing, rely on the stability of the interface between the resident and displacing fluids for efficient removal of contaminants. Contrasts in density and viscosity between the resident and displacing fluids can adversely affect the stability of the displacement front. Petroleum engineers have developed techniques to describe these types of processes; however, their findings do not necessarily translate directly to aquifer remediation. The purpose of this laboratory study was to investigate how density and viscosity contrasts affected cosolvent displacements in unconfined porous media characterized by the presence of a capillary fringe. Two-dimensional flow laboratory experiments, which were partially scaled to a cosolvent flushing field experiment, were conducted to determine potential implications of flow instabilities in homogeneous sand packs. Numerical simulations were also conducted to investigate the differential impact of fluid property contrasts in unconfined and confined systems. The results from these experiments and simulations indicated that the presence of a capillary fringe was an important factor in the displacement efficiency. Buoyant forces can act to carry a lighter-than-water cosolvent preferentially into the capillary fringe during displacement of the resident groundwater. During subsequent water flooding, buoyancy forces can act to effectively trap the cosolvent in the capillary fringe, contributing to the inefficient removal of cosolvent from the aquifer.     

Axi-symmetric simulation of soil vapor extraction influenced by soil fracturing

Fracturing, either pneumatic or hydraulic, is a method to improve the performance of soil vapor extraction (SVE) in relatively low permeability soils (< 10(-5) cm/s). A two-dimensional model is presented to simulate trichloroethylene (TCE) soil vapor extraction modified by fracturing. Flow and transport is modeled using mobile macropore and micropore networks, which also have been identified in the literature as dual porosity, dual permeability, or heterogeneous flow models. In this model, fluids can flow in both the macropore and micropore networks. This represents a more general model compared to immobile micropore, mobile macropore models presented thus far in the literature for vapor flow and transport in two dimensions. The model considers pressure- and concentration-driven exchange between the macropore and micropore networks, concentration-driven exchange between the gas and sorbed phases within each network, and equilibrium exchange between the gas and water and a sorbed phase within each network. The parameters employed in an example simulation are based on field measurements made at a fractured site. Considered in the simulations were the influence of the volume percentage of fractures, the length of fractures, the relative location of the water table, and the influence of pulsed pumping. For these simulations, internetwork concentration-driven exchange most significantly affected mass removal. The volume percentage of fractures more significantly influence flow and mass removal than the length of fractures. The depth of the water table below the contamination plume only significantly influenced flow and mass removal when the water table was within 60 cm of the bottom of the contaminated soil in the vadose zone for the parameters considered in this study. Pulsed pumping was not found to increase the amount of mass removed in this study.