Use of tandem circulation wells to measure hydraulic conductivity without groundwater extraction

Conventional methods to measure the hydraulic conductivity of an aquifer on a relatively large scale (10-100 m) require extraction of significant quantities of groundwater. This can be expensive, and otherwise problematic, when investigating a contaminated aquifer. In this study, innovative approaches that make use of tandem circulation wells to measure hydraulic conductivity are proposed. These approaches measure conductivity on a relatively large scale, but do not require extraction of groundwater. Two basic approaches for using circulation wells to measure hydraulic conductivity are presented; one approach is based upon the dipole-flow test method, while the other approach relies on a tracer test to measure the flow of water between two recirculating wells. The approaches are tested in a relatively homogeneous and isotropic artificial aquifer, where the conductivities measured by both approaches are compared to each other and to the previously measured hydraulic conductivity of the aquifer. It was shown that both approaches have the potential to accurately measure horizontal and vertical hydraulic conductivity for a relatively large subsurface volume without the need to pump groundwater to the surface. Future work is recommended to evaluate the ability of these tandem circulation wells to accurately measure hydraulic conductivity when anisotropy and heterogeneity are greater than in the artificial aquifer used for these studies.     

Quantification of groundwater contamination in an urban area using integral pumping tests

In this paper, the integral groundwater investigation method is used for the quantification of PCE and TCE mass flow rates at an industrialized urban area in Linz, Austria. In this approach, pumping wells positioned along control planes perpendicular to the groundwater flow direction are operated for a time period on the order of days and sampled for contaminants. The concentration time series of the contaminants measured during operation of the pumping wells are then used to determine contaminant mass flow rates, mean concentrations and the plume shapes and positions at the control planes. The three control planes used in Linz were positioned downstream of a number of potential source zones, which are distributed over the field site. By use of the integral investigation method, it was possible to identify active contaminant sources, quantify the individual source strength in terms of mass flow rates at the control planes and estimate the contaminant plume position relative to the control planes. The source zones emitting the highest PCE and TCE mass flow rates could be determined, representing the areas where additional investigation and remediation activities will be needed. Additionally, large parts of the area investigated could be excluded from further investigation and remediation activities.     

Natural attenuation of chlorinated ethene compounds: model development and field-scale application at the Dover site

A multi-dimensional and multi-species reactive transport model was developed to aid in the analysis of natural attenuation design at chlorinated solvent sites. The model can simulate several simultaneously occurring attenuation processes including aerobic and anaerobic biological degradation processes. The developed model was applied to analyze field-scale transport and biodegradation processes occurring at the Area-6 site in Dover Air Force Base, Delaware. The model was calibrated to field data collected at this site. The calibrated model reproduced the general groundwater flow patterns, and also, it successfully recreated the observed distribution of tetrachloroethene (PCE), trichloroethene (TCE), dichloroethylene (DCE), vinyl chloride (VC) and chloride plumes. Field-scale decay rates of these contaminant plumes were also estimated. The decay rates are within the range of values that were previously estimated based on lab-scale microcosm and field-scale transect analyses. Model simulation results indicated that the anaerobic degradation rate of TCE, source loading rate, and groundwater transport rate are the important model parameters. Sensitivity analysis of the model indicated that the shape and extent of the predicted TCE plume is most sensitive to transmissivity values. The total mass of the predicted TCE plume is most sensitive to TCE anaerobic degradation rates. The numerical model developed in this study is a useful engineering tool for integrating field-scale natural attenuation data within a rational modeling framework. The model results can be used for quantifying the relative importance of various simultaneously occurring natural attenuation processes.     

Hydraulic constraints on the performance of a groundwater denitrification wall for nitrate removal from shallow groundwater

Denitrification walls are a practical approach for decreasing non-point source pollution of surface waters. They are constructed by digging a trench perpendicular to groundwater flow and mixing the aquifer material with organic matter, such as sawdust, which acts as a carbon source to stimulate denitrification. For efficient functioning, walls need to be permeable to groundwater flow. We examined the functioning of a denitrification wall constructed in an aquifer consisting of coarse sands. Wells were monitored for changes in nitrate concentration as groundwater passed through the wall and soil samples were taken to measure microbial parameters inside the wall. Nitrate concentrations upstream of the wall ranged from 21 to 39 g N m(-3), in the wall from 0 to 2 g N m(-3) and downstream from 19 to 44 g N m(-3). An initial groundwater flow investigation using a salt tracer dilution technique showed that the flow through the wall was less than 4% of the flow occurring in the aquifer. Natural gradient tracer tests using bromide and Rhodamine-WT confirmed groundwater bypass under the wall. Hydraulic conductivity of 0.48 m day(-1) was measured inside the wall, whereas the surrounding aquifer had a hydraulic conductivity of 65.4 m day(-1). This indicated that during construction of the wall, hydraulic conductivity of the aquifer had been greatly reduced, so that most of the groundwater flowed under rather than through the wall. Denitrification rates measured in the center of the wall ranged from 0.020 to 0.13 g N m(-3) day(-1), which did not account for the rates of nitrate removal (0.16-0.29 g N m(-3) day(-1)) calculated from monitoring of groundwater nitrate concentrations. This suggested that the rate of denitrification was greater at the upstream face of the wall than in its center where it was limited by low nitrate concentrations. While denitrification walls can be an inexpensive tool for removing nitrate from groundwater, they may not be suitable in aquifers with coarse textured subsoils where simple inexpensive construction techniques result in major decreases in hydraulic conductivity.     

Importance of considering intraborehole flow in solute transport modeling under highly dynamic flow conditions

Correct interpretation of tracer test data is critical for understanding transport processes in the subsurface. This task can be greatly complicated by the presence of intraborehole flows in a highly dynamic flow environment. At a new tracer test site (Hanford IFRC) a dynamic flow field created by changes in the stage of the adjacent Columbia River, coupled with a heterogeneous hydraulic conductivity distribution, leads to considerable variations in vertical hydraulic gradients. These variations, in turn, create intraborehole flows in fully-screened (6.5m) observation wells with frequently alternating upward and downward movement. This phenomenon, in conjunction with a highly permeable aquifer formation and small horizontal hydraulic gradients, makes modeling analysis and model calibration a formidable challenge. Groundwater head data alone were insufficient to define the flow model boundary conditions, and the movement of the tracer was highly sensitive to the dynamics of the flow field. This study shows that model calibration can be significantly improved by explicitly considering (a) dynamic flow model boundary conditions and (b) intraborehole flow. The findings from this study underscore the difficulties in interpreting tracer tests and understanding solute transport under highly dynamic flow conditions.     

Description of hydrogeologic heterogeneity and evaluation of radionuclide transport at an underground nuclear test

Realistic models of contaminant transport in groundwater demand detailed characterization of the spatial distribution of subsurface hydraulic properties, while at the same time programmatic constraints may limit collection of pertinent hydraulic data. Fortunately, alternate forms of data can be used to improve characterization of spatial variability. We utilize a methodology that augments sparse hydraulic information (hard data) with more widely available hydrogeologic information to generate equiprobable maps of hydrogeologic properties that incorporate patterns of connected permeable zones. Geophysical and lithologic logs are used to identify hydrogeologic categories and to condition stochastic simulations using Sequential Indicator Simulation (SIS). The resulting maps are populated with hydraulic conductivity values using field data and Sequential Gaussian Simulation (SGS). Maps of subsurface hydrogeologic heterogeneity are generated for the purpose of examining groundwater flow and transport processes at the Faultless underground nuclear test, Central Nevada Test Area (CNTA), through large-scale, three-dimensional numerical modeling. The maps provide the basis for simulation of groundwater flow, while transport of radionuclides from the nuclear cavity is modeled using particle tracking methods. Sensitivity analyses focus on model parameters that are most likely to reduce the long travel times observed in the base case. The methods employed in this study have improved our understanding of the spatial distribution of preferential flowpaths at this site and provided the critical foundation on which to build models of groundwater flow and transport. The results emphasize that the impacts of uncertainty in hydraulic and chemical parameters are dependent on the radioactive decay of specific species, with rapid decay magnifying the effects of parameters that change travel time.     

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

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

A controlled field experiment on groundwater contamination by a multicomponent DNAPL: creation of the emplaced-source and overview of dissolved plume development

A unique field experiment has been undertaken at the CFB Borden research site to investigate the development of dissolved chlorinated solvent plumes from a residual dense non-aqueous phase liquid (DNAPL) source. The "emplaced-source" tracer test methodology involved a controlled emplacement of a block-shaped source of sand containing chlorinated solvents below the water table. The gradual dissolution of this residual DNAPL solvent source under natural aquifer conditions caused dissolved solvent plumes of trichloromethane (TCM), trichloroethene (TCE) and perchloroethene (PCE) to continuously develop down gradient. Source dissolution and 3-D plume development were successfully monitored via 173 multilevel samplers over a 475-day tracer test period prior to site remediation research being initiated. Detailed groundwater level and hydraulic conductivity data were collected. Development of plumes with concentrations spanning 1-700,000 micrograms/1 is described and key processes controlling their migration identified. Plumes were observed to be narrow due to the weakness of transverse dispersion processes and long due to advection and significant longitudinal dispersion, very limited sorptive retardation and negligible, if any, attenuation due to biodegradation or abiotic reaction. TCM was shown to be essentially conservative, TCE very nearly conservative and PCE, consistent with its greater hydrophobicity, more retarded yet having a greater mobility than observed in previous Borden field tests. The absence of biodegradation was ascribed to the prevailing aerobic conditions and lack of any additional biodegradable carbon substrates. The transient groundwater flow regime caused significant transverse lateral plume movement, plume asymmetry and was likely responsible for most of the, albeit limited, transverse horizontal plume spreading. In agreement with the widespread incidence of extensive TCE and PCE plumes throughout the industrialized world, the experiment indicates such solvent plumes are likely to be highly mobile and persistent, at least in aquifers that are aerobic and have low sorption potential (low foc content).     

Push-pull partitioning tracer tests using radon-222 to quantify non-aqueous phase liquid contamination

Naturally occurring radon in groundwater can be used as an in situ partitioning tracer for locating and quantifying non-aqueous phase liquid (NAPL) contamination in the subsurface. When combined with the single-well, push-pull test, this methodology has the potential to provide a low-cost alternative to inter-well partitioning tracer tests. During a push-pull test, a known volume of test solution (radon-free water containing a conservative tracer) is first injected ("pushed") into a well; flow is then reversed and the test solution/groundwater mixture is extracted ("pulled") from the same well. In the presence of NAPL radon transport is retarded relative to the conservative tracer. Assuming linear equilibrium partitioning, retardation factors for radon can be used to estimate NAPL saturations. The utility of this methodology was evaluated in laboratory and field settings. Laboratory push-pull tests were conducted in both non-contaminated and trichloroethene NAPL (TCE)-contaminated sediment. The methodology was then applied in wells located in non-contaminated and light non-aqueous phase liquid (LNAPL)-contaminated portions of an aquifer at a former petroleum refinery. The method of temporal moments and an approximate analytical solution to the governing transport equations were used to interpret breakthrough curves and estimate radon retardation factors; estimated retardation factors were then used to calculate TCE saturations. Numerical simulations were used to further investigate the behavior of the breakthrough curves. The laboratory and field push-pull tests demonstrated that radon retardation does occur in the presence of TCE and LNAPL and that radon retardation can be used to calculate TCE saturations. Laboratory injection-phase test results in TCE-contaminated sediment yielded radon retardation factors ranging from 1.1 to 1.5, resulting in calculated TCE saturations ranging from 0.2 to 0.9%. Laboratory extraction-phase test results in the same sediment yielded a radon retardation factor of 5.0, with a calculated TCE saturation of 6.5%. Numerical simulation breakthrough curves provided reasonably good matches to the approximate analytical solution breakthrough curves. However, non-equilibrium radon partitioning and heterogeneous TCE distributions may affect the retardation factors and TCE saturation estimates.     

Design of aquifer remediation systems: (1) describing hydraulic structure and NAPL architecture using tracers

Aquifer heterogeneity (structure) and NAPL distribution (architecture) are described based on tracer data. An inverse modelling approach that estimates the hydraulic structure and NAPL architecture based on a Lagrangian stochastic model where the hydraulic structure is described by one or more populations of lognormally distributed travel times and the NAPL architecture is selected from eight possible assumed distributions. Optimization of the model parameters for each tested realization is based on the minimization of the sum of the square residuals between the log of measured tracer data and model predictions for the same temporal observation. For a given NAPL architecture the error is reduced with each added population. Model selection was based on a fitness which penalized models for increasing complexity. The technique is demonstrated under a range of hydrologic and contaminant settings using data from three small field-scale tracer tests: the first implementation at an LNAPL site using a line-drive flow pattern, the second at a DNAPL site with an inverted five-spot flow pattern, and the third at the same DNAPL site using a vertical circulation flow pattern. The Lagrangian model was capable of accurately duplicating experimentally derived tracer breakthrough curves, with a correlation coefficient of 0.97 or better. Furthermore, the model estimate of the NAPL volume is similar to the estimates based on moment analysis of field data.     

Simulating the fate and transport of TCE from groundwater to indoor air

This work provides an exploratory analysis on the relative importance of various factors controlling the fate and transport of volatile organic contaminants (in this case, TCE) from a DNAPL source zone located below the water table and into the indoor air. The analysis is conducted using the multi-phase compositional model CompFlow Bio, with the base scenario problem geometry reminiscent of a field experiment conducted by Rivett [Rivett, M.O., (1995), Soil–gas signatures from volatile chlorinated solvents: Borden field experiments. Groundwater, 33(1), 84–98.] at the Borden aquifer where groundwater was observed to transport a contaminant plume a substantial distance without vertical mass transport of the contaminant across the capillary fringe and into the vadose zone. Results for the base scenario model indicate that the structure of the permeability field was largely responsible for deflecting the groundwater plume upward towards the capillary fringe, permitting aqueous phase diffusion to transport the TCE into the vadose zone. Alternative permeability realizations, generated as part of a Monte Carlo simulation process, at times deflected the groundwater plume downwards causing the extended thickness of the saturated zone to insulate the vadose zone from exposure to the TCE by upward diffusive transport. Comparison of attenuation coefficients calculated using the CompFlow Bio and Johnson and Ettinger [Johnson, P.C. and Ettinger, R.A., (1991), Heuristic model for predicting the intrusion rate of contaminant vapors into buildings. Environmental Science and Technology, 25, 1445–1452.] heuristic model exhibited fortuitous agreement for the base scenario problem geometry, with this agreement diverging for the alternative permeability realizations as well as when parameters such as the foundation slab fracture aperture, the indoor air pressure drop, the capillary fringe thickness, and the infiltration rate were varied over typical ranges.     

First-passage-time transfer functions for groundwater tracer tests conducted in radially convergent flow

Forced-gradient groundwater tracer tests may be conducted using a variety of hydraulic schemes, so it is useful to have simple semi-analytic models available that can examine various injection/withdrawal scenarios. Models for radially convergent tracer tests are formulated here as transfer functions, which allow complex tracer test designs to be simulated by a series of simple mathematical expressions. These mathematical expressions are given in Laplace space, so that transfer functions may be placed in series by simple multiplication. Predicted breakthrough is found by numerically inverting the composite transfer function to the time-domain, using traditional computer programs or commercial mathematical software. Transport is assumed to be dictated by a radially convergent or uniform flow field, and is based upon an exact first-passage-time solution of the backward Fokker–Planck equation. These methods are demonstrated by simulating a weak-dipole tracer test conducted in a fractured granite formation, where mixing in the injection borehole is non-ideal.     

Design of aquifer remediation systems: (1) Describing hydraulic structure and NAPL architecture using tracers

Aquifer heterogeneity (structure) and NAPL distribution (architecture) are described based on tracer data. An inverse modelling approach that estimates the hydraulic structure and NAPL architecture based on a Lagrangian stochastic model where the hydraulic structure is described by one or more populations of lognormally distributed travel times and the NAPL architecture is selected from eight possible assumed distributions. Optimization of the model parameters for each tested realization is based on the minimization of the sum of the square residuals between the log of measured tracer data and model predictions for the same temporal observation. For a given NAPL architecture the error is reduced with each added population. Model selection was based on a fitness which penalized models for increasing complexity. The technique is demonstrated under a range of hydrologic and contaminant settings using data from three small field-scale tracer tests: the first implementation at an LNAPL site using a line-drive flow pattern, the second at a DNAPL site with an inverted five-spot flow pattern, and the third at the same DNAPL site using a vertical circulation flow pattern. The Lagrangian model was capable of accurately duplicating experimentally derived tracer breakthrough curves, with a correlation coefficient of 0.97 or better. Furthermore, the model estimate of the NAPL volume is similar to the estimates based on moment analysis of field data.     

Measurements of groundwater velocity in discrete rock fractures

Estimating groundwater velocity in fracture networks using a Darcy or cubic law calculation is complicated by the wide distribution of fracture aperture often found in these systems and by the difficulty in measuring hydraulic head in discrete fracture features. Although difficult to conduct in a fractured rock setting, the point dilution method can be utilized to collect direct measurements of groundwater velocity in individual fractures. To compare measured against calculated velocities, more than 100 point dilution experiments were conducted within a 35 x 35 m area of a single fracture and in discrete fracture features within a fracture network at a larger scale. The dilution experiments were conducted by isolating a fracture feature in a borehole, measuring the hydraulic aperture, and measuring the decay of an injected tracer due to the advective groundwater flux across the fracture. Groundwater velocity was estimated using the hydraulic aperture and the rate of decay of the injected tracer. Estimates of the local hydraulic gradient were calculated via the cubic law using the velocity estimate and the hydraulic aperture. The results of the tests conducted in the single fracture show variable (1 to 33 m/day) but on average higher velocities in comparison to that measured during a natural gradient tracer experiment conducted previously (in which the effects of matrix diffusion were accounted for) and to that which would be calculated using the cubic law. Based on these results, it was determined that the best estimate of the average groundwater velocity, at the scale of the measurement area used for the cubic law calculations, could only be obtained using the largest apertures in the aperture distribution. Variability of the velocity measurements was also observed over time. Increases in velocity were attributed to the effect of rainfall although concurrent increases in hydraulic gradient were not detected (likely within the tolerance of the measuring devices). The groundwater velocities measured in the fracture network varied over a wider range than at the scale of the single fracture (from 2 to 388 m/day). No correlation, however, was observed between the size of the fracture aperture and measured velocity.     

Groundwater flow and contaminant transport modeling applications in urban area: scopes and limitations

Background

A three-dimensional groundwater flow model was used to evaluate the groundwater potential and assess the effects of groundwater withdrawal on the regional water level and flow direction in the central Beijing area. A program of groundwater modeling aimed at estimating current contaminant fluxes to the central area and site streams via groundwater was developed.

Results and discussion

The conceptual model developed for the site attempted to incorporate a complex stratigraphic profile in which groundwater flow and contaminant transport is strongly controlled by a shallow aquifer. Here, a conceptual model for groundwater flow and contaminant transport in central Beijing is presented.

Conclusion

Model simulations indicated that a sharp drop in the hydraulic head occurs at the center of the model area, which generates a cone of depression and a continuous decline of head with respect to time as a result of heavy groundwater abstraction.     

Use of principal component analysis to profile temporal and spatial variations of chlorinated solvent concentration in groundwater

This work aims at evaluating spatial distribution patterns of concentration variations for chlorinated solvents in groundwater, based on principal component analysis and geographic information system (GIS) tools. The study investigates long-time series of chlorinated solvent concentrations in groundwater measured for 18 contaminated industrial sites. The characterization of contaminant plumes and delineation of pollutant sources are essential for choosing appropriate monitoring and remediation strategies, as contaminated groundwaters are characterized by complex patterns of spatial and temporal concentration variability, with wide unpredictable fluctuations over time. The present work describes the results of a new exploratory statistical method called the Variability Index Method (VIM) applied to environmental data to assess the performance of using concentration variations as molecular tracers to reveal aquifer dynamics, industrial impacts, and point sources for contamination plumes. The application of this method provides a useful assessment of controls over contaminant concentration variations as well as support for remediation techniques.     

Fluorescent dye imaging of the volume sampled by single well forced-gradient tracer tests evaluated in a laboratory-scale aquifer physical model

This study presents a new method to visualise forced-gradient tracer tests in 2-D using a laboratory-scale aquifer physical model. Experiments were designed to investigate the volume of aquifer sampled in vertical dipole flow tracer tests (DFTT) and push-pull tests (PPT), using a miniature monitoring well and straddle packer arrangement equipped with solute injection and recovery chambers. These tests have previously been used to estimate bulk aquifer hydraulic and transport properties for the evaluation of natural attenuation and other remediation approaches. Experiments were performed in a silica glass bead-filled box, using a fluorescent tracer (fluorescein) to deduce conservative solute transport paths. Digital images of fluorescein transport were captured under ultraviolet light and processed to analyse tracer plume geometry and obtain point-concentration breakthrough histories. Inorganic anion mixtures were also used to obtain conventional tracer breakthrough histories. Concentration data from the conservative tracer breakthrough curves was compared with the digital images and a well characterised numerical model. The results show that the peak tracer breakthrough response in dipole flow tracer tests samples a zone of aquifer close to the well screen, while the sampling volume of push-pull tests is limited by the length of the straddle packers used. The effective sampling volume of these single well forced-gradient tests in isotropic conditions can be estimated with simple equations. The experimental approach offers the opportunity to evaluate under controlled conditions the theoretical basis, design and performance of DFTTs and PPTs in porous media in relation to measured flow and transport properties.     

Dimensionless analysis of two analytical solutions for 3-D solute transport in groundwater

Analytical solutions are widely used as screening tools for estimating the potential for contaminant transport in groundwater, or for interpreting tracer tests or groundwater quality data. A solution for three-dimensional solute migration from a plane-source source that is frequently used in practice is the approximate solution of Domenico [J. Hydrol. 91 (1987) 49-58]. A more rigorous solution to the same problem was provided by Sagar [ASCE, J. Hydraul. Div. 108, no. HY1 (1982) 47-62]. A comprehensive and unambiguous comparison between these two solutions is provided using dimensionless analysis. The solutions are first cast in terms of dimensionless parameters and then used to provide type curves covering a wide range of dimensionless parameter values. Results show that while discrepancies between the two solutions are relatively negligible along the plume centreline (for flow regimes dominated by advection and mechanical dispersion), large concentration differences can be observed as lateral distance from the centreline increases, especially in the presence of solute decay.     

A comparison of numerical solutions of partial differential equations with probabilistic and possibilistic parameters for the quantification of uncertainty in subsurface solute transport

Traditionally, uncertainty in parameters are represented as probabilistic distributions and incorporated into groundwater flow and contaminant transport models. With the advent of newer uncertainty theories, it is now understood that stochastic methods cannot properly represent non random uncertainties. In the groundwater flow and contaminant transport equations, uncertainty in some parameters may be random, whereas those of others may be non random. The objective of this paper is to develop a fuzzy-stochastic partial differential equation (FSPDE) model to simulate conditions where both random and non random uncertainties are involved in groundwater flow and solute transport. Three potential solution techniques namely, (a) transforming a probability distribution to a possibility distribution (Method I) then a FSPDE becomes a fuzzy partial differential equation (FPDE), (b) transforming a possibility distribution to a probability distribution (Method II) and then a FSPDE becomes a stochastic partial differential equation (SPDE), and (c) the combination of Monte Carlo methods and FPDE solution techniques (Method III) are proposed and compared. The effects of these three methods on the predictive results are investigated by using two case studies. The results show that the predictions obtained from Method II is a specific case of that got from Method I. When an exact probabilistic result is needed, Method II is suggested. As the loss or gain of information during a probability–possibility (or vice versa) transformation cannot be quantified, their influences on the predictive results is not known. Thus, Method III should probably be preferred for risk assessments.     

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

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