Numerical simulation of a fine-grained denitrification layer for removing septic system nitrate from shallow groundwater

One of the most common methods to dispose of domestic wastewater involves the release of septic effluent from drains located in the unsaturated zone. Nitrogen from such systems is currently of concern because of nitrate contamination of drinking water supplies and eutrophication of coastal waters. It has been proposed that adding labile carbon sources to septic distribution fields could enhance heterotrophic denitrification and thus reduce nitrate concentrations in shallow groundwater. In this study, a numerical model which solves for variably saturated flow and reactive transport of multiple species is employed to investigate the performance of a drain field design that incorporates a fine-grained denitrification layer. The hydrogeological scenario simulated is an unconfined sand aquifer. The model results suggest that the denitrification layer, supplemented with labile organic carbon, may be an effective means to eliminate nitrogen loading to shallow groundwater. It is also shown that in noncalcareous aquifers, the denitrification reaction may provide sufficient buffering capacity to maintain near neutral pH conditions beneath and down gradient of the drain field. Leaching of excess dissolved organic carbon (DOC) from the denitrification layer is problematic, and causes an anaerobic plume to develop in simulations where the water table is less than 5-6 m below ground surface; this anaerobic plume may lead to other down gradient changes in groundwater quality. A drain field and denitrification layer of smaller dimensions is shown to be just as effective for reducing nitrate, but has the benefit of reducing the excess DOC leached from the layer. This configuration will minimize the impact of wastewater disposal in areas where the water table is as shallow as 3.5 m.     

Multicomponent simulation of wastewater-derived nitrogen and carbon in shallow unconfined aquifers. I. Model formulation and performance

One of the most common methods to dispose of domestic wastewater involves the release of septic effluent from drains located in the unsaturated zone. Nitrogen from such systems is currently of concern because of nitrate contamination of drinking water supplies and eutrophication of coastal waters. The objectives of this study are to develop and assess the performance of a mechanistic flow and reactive transport model which couples the most relevant physical, geochemical and biochemical processes involved in wastewater plume evolution in sandy aquifers. The numerical model solves for variably saturated groundwater flow and reactive transport of multiple carbon- and nitrogen-containing species in a three-dimensional porous medium. The reactive transport equations are solved using the Strang splitting method which is shown to be accurate for Monod and first- and second-order kinetic reactions, and two to four times more efficient than sequential iterative splitting. The reaction system is formulated as a fully kinetic chemistry problem, which allows for the use of several special-purpose ordinary differential equation (ODE) solvers. For reaction systems containing both fast and slow kinetic reactions, such as the combined nitrogen-carbon system, it is found that a specialized stiff explicit solver fails to obtain a solution. An implicit solver is more robust and its computational performance is improved by scaling of the fastest reaction rates. The model is used to simulate wastewater migration in a 1-m-long unsaturated column and the results show significant oxidation of dissolved organic carbon (DOC), the generation of nitrate by nitrification, and a slight decrease in pH.     

Investigation of hydrogeologic processes in a dipping layer structure: 2. Transport and biodegradation of organics

Numerical simulation tools have been used to study the dominating processes during transport of aromatic hydrocarbons in the unsaturated soil zone. Simulations were based on field observations at an experimental site located on a glacial delta plain with pronounced layered sedimentary structures. A numerical model for transport in the unsaturated zone, SWMS-3D, has been extended to incorporate coupled multispecies transport, microbial degradation following Monod kinetics and gas diffusive transport of oxygen and hydrocarbons. The flow field parameters were derived from previous work using nonreactive tracers. Breakthrough curves (BTC) from the hydrocarbon field experiment were used to determine sorption parameters and Monod kinetic parameters using a fitting procedure. The numerical simulations revealed that the assumption of homogeneous layers resulted in deviations from the field observations. The deviations were more pronounced with incorporation of reactive transport, compared with earlier work on nonreactive transport. To be able to model reasonable BTC, sorption had to be reduced compared to laboratory experiments. The initial biomass and the maximum utilisation rate could be adjusted to capture both the initial lag phase and the overall degradation rate. Nevertheless, local oxygen limitation is predicted by the model, which was not observed in the field experiment. Incorporation of evaporation and diffusive gas transport of the hydrocarbons did not significantly change the local oxygen demand. The main cause of the observed discrepancies between model and field are attributed to channelling as a result of small-scale heterogeneities such as biopores.     

Electrokinetic ion transport through unsaturated soil: 1. Theory, model development, and testing

An electromigration transport model for non-reactive ion transport in unsaturated soil was developed and tested against laboratory experiments. This model assumed the electric potential field was constant with respect to time, an assumption valid for highly buffered soil, or when the electrode electrolysis reactions are neutralized. The model also assumed constant moisture contents and temperature with respect to time, and that electroosmotic and hydraulic transport of water through the soil was negligible. A functional relationship between ionic mobility and the electrolyte concentration was estimated using the chemical activity coefficient. Tortuosity was calculated from a mathematical relationship fitted to the electrical conductivity of the bulk pore water and soil moisture data. The functional relationship between ionic mobility, pore-water concentration, and tortuosity as a function of moisture content allowed the model to predict ion transport in heterogeneous unsaturated soils. The model was tested against laboratory measurements assessing anionic electromigration as a function of moisture content. In the test cell, a strip of soil was spiked with red dye No 40 and monitored for a 24-h period while a 10-mA current was maintained between the electrodes. Electromigration velocities predicted by the electromigration transport model were in agreement with laboratory experimental results. Both laboratory-measured and model-predicted dye migration results indicated a maximum transport velocity at moisture contents less than saturation due to competing effects between current density and tortuosity as moisture content decreases.     

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.     

Analytical modelling of stable isotope fractionation of volatile organic compounds in the unsaturated zone

Analytical models were developed that simulate stable isotope ratios of volatile organic compounds (VOCs) near a point source contamination in the unsaturated zone. The models describe diffusive transport of VOCs, biodegradation and source ageing. The mass transport is governed by Fick's law for diffusion. The equation for reactive transport of VOCs in the soil gas phase was solved for different source geometries and for different boundary conditions. Model results were compared to experimental data from a one-dimensional laboratory column and a radial-symmetric field experiment. The comparison yielded a satisfying agreement. The model results clearly illustrate the significant isotope fractionation by gas phase diffusion under transient state conditions. This leads to an initial depletion of heavy isotopes with increasing distance from the source. The isotope evolution of the source is governed by the combined effects of isotope fractionation due to vaporisation, diffusion and biodegradation. The net effect can lead to an enrichment or depletion of the heavy isotope in the remaining organic phase, depending on the compound and element considered. Finally, the isotope evolution of molecules migrating away from the source and undergoing degradation is governed by a combined degradation and diffusion isotope effect. This suggests that, in the unsaturated zone, the interpretation of biodegradation of VOC based on isotopic data must always be based on a model combining gas phase diffusion and degradation.     

Characterization of flow and transport processes within the unsaturated zone of Yucca Mountain, Nevada, under current and future climates

This paper presents a large-scale modeling study characterizing fluid flow and tracer transport in the unsaturated zone of Yucca Mountain, Nevada, a potential repository site for storing high-level radioactive waste. The study has been conducted using a three-dimensional numerical model, which incorporates a wide variety of field data and takes into account the coupled processes of flow and transport in the highly heterogeneous, unsaturated fractured porous rock. The modeling approach is based on a dual-continuum formulation of coupled multiphase fluid and tracer transport through fractured porous rock. Various scenarios of current and future climate conditions and their effects on the unsaturated zone are evaluated to aid in the assessment of the proposed repository's system performance using different conceptual models. These models are calibrated against field-measured data. Model-predicted flow and transport processes under current and future climates are discussed.     

Oxygenated gasoline release in the unsaturated zone, Part 2: Downgradient transport of ethanol and hydrocarbons

In the event of a gasoline spill containing oxygenated compounds such as ethanol and MTBE, it is important to consider the impacts these compounds might have on subsurface contamination. One of the main concerns commonly associated with ethanol is that it might decrease the biodegradation of aromatic hydrocarbon compounds, leading to an increase in the hydrocarbon dissolved plume lengths. The first part of this study (Part 1) showed that when gasoline containing ethanol infiltrates the unsaturated zone, ethanol is likely to partition to and be retained in the unsaturated zone pore water. In this study (Part 2), a controlled field test is combined with a two-dimensional laboratory test and three-dimensional numerical modelling to investigate how ethanol retention in the unsaturated zone affects the downgradient behaviour of ethanol and aromatic hydrocarbon compounds. Ethanol transport downgradient was extremely limited. The appearance of ethanol in downgradient wells was delayed and the concentrations were lower than would be expected based on equilibrium dissolution. Oscillations in the water table resulted in minor flushing of ethanol, but its effect could still be perceived as an increase in the groundwater concentrations downgradient from the source zone. Ethanol partitioning to the unsaturated zone pore water reduced its mass fraction within the NAPL thus reducing its anticipated impact on the fate of the hydrocarbon compounds. A conceptual numerical simulation indicated that the potential ethanol-induced increase in benzene plume length after 20 years could decrease from 136% to 40% when ethanol retention in the unsaturated zone is considered.     

Electrokinetic ion transport through unsaturated soil: 2. Application to a heterogeneous field site

Results of a field demonstration of electrokinetic transport of acetate through an unsaturated heterogeneous soil are compared to numerical modeling predictions. The numerical model was based on the groundwater flow and transport codes MODFLOW and MT3D modified to account for electrically induced ion transport. The 6-month field demonstration was conducted in an unsaturated layered soil profile where the soil moisture content ranged from 4% to 28% (m3 m(-3)). Specially designed ceramic-cased electrodes maintained a steady-state moisture content and electric potential field between the electrodes during the field demonstration. Acetate, a byproduct of acetic acid neutralization of the cathode electrolysis reaction, was transported from the cathode to the anode by electromigration. Field demonstration results indicated preferential transport of acetate through soil layers exhibiting higher moisture content/electrical conductivity. These field transport results agree with theoretical predictions that electromigration velocity is proportional to a power function of the effective moisture content. A numerical model using a homogeneous moisture content/electrical conductivity domain did not adequately predict the acetate field results. Numerical model predictions using a three-layer electrical conductivity/moisture content profile agreed qualitatively with the observed acetate distribution. These results suggest that field heterogeneities must be incorporated into electrokinetic models to predict ion transport at the field-scale.     

Simulation of carrier-facilitated transport of phenanthrene in a layered soil profile

The appropriate prediction of the fate of the contaminant is an essential step when evaluating the risk of severe groundwater pollutions-in particular in the context of natural attenuation. We numerically study the reactive transport of phenanthrene at the field scale in a multilayer soil profile based on experimental data. The effect of carrier facilitation by dissolved organic carbon is emphasized and incorporated in the model. Previously published simulations are restricted to the saturated zone and/or to homogeneous soil columns at the laboratory scale. A numerical flow and transport model is extended and applied to understand and quantify the relevant processes in the case of a strongly sorbing hydrophobic organic compound that is subject to carrier facilitation in the unsaturated zone. The contaminant migration is investigated on long- and short-term time scales and compared to predictions without carrier facilitation. The simulations demonstrate the importance of carrier facilitation and suggest strongly to take this aspect into account. By carrier facilitation breakthrough times at the groundwater level decreased from 500 to approximately 8 years and concentration peaks increased by two orders of magnitude in the long-term simulation assuming a temporary spill in an initially unpolluted soil with a non-sorbing carrier.     

Fate of nitrogen for subsurface drip dispersal of effluent from small wastewater systems

Subsurface drip irrigation systems apply effluent from onsite wastewater systems in a more uniform manner at a lower rate than has been possible with other effluent dispersal methods. The effluent is dispersed in a biologically active part of the soil profile for optimal treatment and where the water and nutrients can be utilized by landscape plants. Container tests were performed to determine the fate of water and nitrogen compounds applied to packed loamy sand, sandy loam, and silt loam soils. Nitrogen removal rates measured in the container tests ranged from 63 to 95% despite relatively low levels of available carbon. A Hydrus 2D vadose zone model with nitrification and denitrification rate coefficients calculated as a function of soil moisture content fit the container test results reasonably well. Model results were sensitive to the denitrification rate moisture content function. Two-phase transport parameters were needed to model the preferential flow conditions in the finer soils. Applying the model to generic soil types, the greatest nitrogen losses (30 to 70%) were predicted for medium to fine texture soils and soils with restrictive layers or capillary breaks. The slow transport with subsurface drip irrigation enhanced total nitrogen losses and plant nitrogen uptake opportunity.     

Review of unsaturated-zone transport and attenuation of volatile organic compound (VOC) plumes leached from shallow source zones

Reliable prediction of the unsaturated zone transport and attenuation of dissolved-phase VOC (volatile organic compound) plumes leached from shallow source zones is a complex, multi-process, environmental problem. It is an important problem as sources, which include solid-waste landfills, aqueous-phase liquid discharge lagoons and NAPL releases partially penetrating the unsaturated zone, may persist for decades. Natural attenuation processes operating in the unsaturated zone that, uniquely for VOCs includes volatilisation, may, however, serve to protect underlying groundwater and potentially reduce the need for expensive remedial actions. Review of the literature indicates that only a few studies have focused upon the overall leached VOC source and plume scenario as a whole. These are mostly modelling studies that often involve high strength, non-aqueous phase liquid (NAPL) sources for which density-induced and diffusive vapour transport is significant. Occasional dissolved-phase aromatic hydrocarbon controlled infiltration field studies also exist. Despite this lack of focus on the overall problem, a wide range of process-based unsaturated zone - VOC research has been conducted that may be collated to build good conceptual model understanding of the scenario, particularly for the much studied aromatic hydrocarbons and chlorinated aliphatic hydrocarbons (CAHs). In general, the former group is likely to be attenuated in the unsaturated zone due to their ready aerobic biodegradation, albeit with rate variability across the literature, whereas the fate of the latter is far less likely to be dominated by a single mechanism and dependent upon the relative importance of the various attenuation processes within individual site - VOC scenarios. Analytical and numerical modelling tools permit effective process representation of the whole scenario, albeit with potential for inclusion of additional processes - e.g., multi-mechanistic sorption phase partitioning, and provide good opportunity for further sensitivity analysis and development to practitioner use. There remains a significant need to obtain intermediate laboratory-scale and particularly field-scale (actual site and controlled release) datasets that address the scenario as a whole and permit validation of the available models. Integrated assessment of the range of simultaneous processes that combine to influence leached plume generation, transport and attenuation in the unsaturated zone is required. Component process research needs are required across the problem scenario and include: the simultaneous volatilisation and dissolution of source zones; development of appropriate field-scale dispersion estimates for the unsaturated zone; assessment of transient VOC exchanges between aqueous, vapour and sorbed phases and their influence upon plume attenuation; development of improved field methods to recognise and quantify biodegradation of CAHs; establishment of the influence of co-contaminants; and, finally, translation of research findings into more robust practitioner practice.     

Modelling the geochemical fate and transport of wastewater-derived phosphorus in contrasting groundwater systems

A 1D reactive transport model (RTM) is used to obtain a mechanistic understanding of the fate of phosphorus (P) in the saturated zone of two contrasting aquifer systems. We use the field data from two oxic, electron donor-poor, wastewater-impacted, sandy Canadian aquifers, (Cambridge and Muskoka sites) as an example of a calcareous and non-calcareous groundwater system, respectively, to validate our reaction network. After approximately 10 years of wastewater infiltration, P is effectively attenuated within the first 10 m down-gradient of the source mainly through fast sorption onto calcite and Fe oxides. Slow, kinetic sorption contributes further to P removal, while precipitation of phosphate minerals (strengite, hydroxyapatite) is quantitatively unimportant in the saturated zone. Nitrogen (N) dynamics are also considered, but nitrate behaves essentially as a conservative tracer in both systems. The model-predicted advancement of the P plume upon continued wastewater discharge at the calcareous site is in line with field observations. Model results suggest that, upon removal of the wastewater source, the P plume at both sites will persist for at least 20 years, owing to desorption of P from aquifer solids and the slow rate of P mineral precipitation. Sensitivity analyses for the non-calcareous scenario (Muskoka) illustrate the importance of the sorption capacity of the aquifer solids for P in modulating groundwater N:P ratios in oxic groundwater. The model simulations predict the breakthrough of groundwater with high P concentrations and low N:P ratios after 17 years at 20 m from the source for an aquifer with low sorption capacity (<0.02% w/w Fe(OH)(3)). In this type of system, denitrification plays a minor role in lowering the N:P ratios because it is limited by the availability of labile dissolved organic matter.     

A triple-continuum approach for modeling flow and transport processes in fractured rock

This paper presents a triple-continuum conceptual model for simulating flow and transport processes in fractured rock. Field data collected from the unsaturated zone of Yucca Mountain, a repository site of high-level nuclear waste, show a large number of small-scale fractures. The effect of these small fractures has not been considered in previous modeling investigations within the context of a continuum approach. A new triple-continuum model (consisting of matrix, small-fracture, and large-fracture continua) has been developed to investigate the effect of these small fractures. This paper derives the model formulation and discusses the basic triple-continuum behavior of flow and transport processes under different conditions, using both analytical solutions and numerical approaches. The simulation results from the site-scale model of the unsaturated zone of Yucca Mountain indicate that these small fractures may have an important effect on radionuclide transport within the mountain.     

Reactive transport modeling of processes controlling the distribution and natural attenuation of phenolic compounds in a deep sandstone aquifer

Reactive solute transport modeling was utilized to evaluate the potential for natural attenuation of a contaminant plume containing phenolic compounds at a chemical producer in the West Midlands, UK. The reactive transport simulations consider microbially mediated biodegradation of the phenolic compounds (phenols, cresols, and xylenols) by multiple electron acceptors. Inorganic reactions including hydrolysis, aqueous complexation, dissolution of primary minerals, formation of secondary mineral phases, and ion exchange are considered. One-dimensional (1D) and three-dimensional (3D) simulations were conducted. Mass balance calculations indicate that biodegradation in the saturated zone has degraded approximately 1-5% of the organic contaminant plume over a time period of 47 years. Simulations indicate that denitrification is the most significant degradation process, accounting for approximately 50% of the organic contaminant removal, followed by sulfate reduction and fermentation reactions, each contributing 15-20%. Aerobic respiration accounts for less than 10% of the observed contaminant removal in the saturated zone. Although concentrations of Fe(III) and Mn(IV) mineral phases are high in the aquifer sediment, reductive dissolution is limited, producing only 5% of the observed mass loss. Mass balance calculations suggest that no more than 20-25% of the observed total inorganic carbon (TIC) was generated from biodegradation reactions in the saturated zone. Simulations indicate that aerobic biodegradation in the unsaturated zone, before the contaminant entered the aquifer, may have produced the majority of the TIC observed in the plume. Because long-term degradation is limited to processes within the saturated zone, use of observed TIC concentrations to predict the future natural attenuation may overestimate contaminant degradation by a factor of 4-5.     

A sensitivity analysis of adsorption and degradation parameters in the modeling of pesticide transport in soils

The increasing use of deterministic models in predicting the movement of pesticides in soils, has focused attention on the evaluation of major parameters which represent attenuation factors of organics in the subsurface. These parameters are the degradation rate constant and the adsorption constant for the pesticide. In view of the large in situ variability of these parameters and of the difficulty in obtaining accurate field data, there is a high degree of uncertainty associated with the results obtained from deterministic models. A sensitivity analysis is performed here to quantify the impact of such variation in each of these input parameters on the output results of an unsaturated zone transport model (PRZM). Results show that variations in these parameters about their respective mean values greatly affect the predicted concentration distributions, obtained after three years, of the pesticide aldicarb in all the soil profile. A 15–22% variation in the degradation constant, or a 24% variation in the adsorption constant, lead to a 100% uncertainty in the various simulation results defined as the cumulative quantity of aldicarb or the dissolved aldicarb concentration leached below the root zone (or the unsaturated zone) of the soil. Such a deterministic model presents a high degree of sensitivity to these input parameters. Accurate field data are then needed to obtain reliable model results in predicting pesticide movement inthe unsaturated zone.     

Reactive transport modeling of column experiments on the evolution of saline-alkaline waste solutions

Leakage of saline-alkaline tank waste solutions often creates a serious environmental contamination problem. To better understand the mechanisms controlling the fate of such waste solutions in the Hanford vadose zone, we simulated reactive transport in columns designed to represent local site conditions. The Pitzer ion interaction module was used, with principal geochemical processes considered in the simulation including quartz dissolution, precipitation of brucite, calcite, and portlandite, multi-component cation exchange, and aqueous complexation reactions. Good matches were observed between the simulated and measured column data at ambient temperature ( approximately 21 degrees C). Relatively good agreement was also obtained at high temperature ( approximately 70 degrees C). The decrease of pH at the plume front is examined through formation of secondary mineral phases and/or quartz dissolution. Substantial formation of secondary mineral phases resulting from multi-component cation exchange suggests that these phases are responsible for a decrease in pH within the plume front. In addition, a sensitivity analysis was conducted with respect to cation exchange capacity, selectivity coefficient, mineral assemblage, temperature, and ionic strength. This study could serve as a useful guide to subsequent experimental work, to thermodynamic models developed for the concentrated solutions at high ionic strength and to other types of waste plume studies.     

Dissolved plume attenuation with DNAPL source remediation, aqueous decay and volatilization--analytical solution, model calibration and prediction uncertainty

A vertically-integrated analytical model for dissolved phase transport is described that considers a time-dependent DNAPL source based on the upscaled dissolution kinetics model of Parker and Park with extensions to consider time-dependent source zone biodecay, partial source mass reduction, and remediation-enhanced source dissolution kinetics. The model also considers spatial variability in aqueous plume decay, which is treated as the sum of aqueous biodecay and volatilization due to diffusive transport and barometric pumping through the unsaturated zone. The model is implemented in Excel/VBA coupled with (1) an inverse solution that utilizes prior information on model parameters and their uncertainty to condition the solution, and (2) an error analysis module that computes parameter covariances and total prediction uncertainty due to regression error and parameter uncertainty. A hypothetical case study is presented to evaluate the feasibility of calibrating the model from limited noisy field data. The results indicate that prediction uncertainty increases significantly over time following calibration, primarily due to propagation of parameter uncertainty. However, differences between the predicted performance of source zone partial mass reduction and the known true performance were reasonably small. Furthermore, a clear difference is observed between the predicted performance for the remedial action scenario versus that for a no-action scenario, which is consistent with the true system behavior. The results suggest that the model formulation can be effectively utilized to assess monitored natural attenuation and source remediation options if careful attention is given to model calibration and prediction uncertainty issues.     

Analytical model of solute transport by unsteady unsaturated gravitational infiltration

Penetration of reactive solute into a soil during a cycle of water infiltration and redistribution is investigated by deriving analytical closed form solutions for fluid flux, moisture content and contaminant concentration. The solution is developed for gravitational flow and advective transport and is applied to two scenarios of solute applications encountered in the applications: a finite pulse of solute dissolved in irrigation water and an instantaneous pulse broadcasted onto the soil surface. Through comparison to simulations of Richards' flow, capillary suction is shown to have contrasting effects on the upper and lower boundaries of the fluid pulse, speeding penetration of the wetting front and reducing the rate of drying. This leads to agreement between the analytical and numerical solutions for typical field and experimental conditions. The analytical solution is further incorporated into a stochastic column model of flow and transport to compute mean solute concentration in a heterogeneous field. An unusual phenomenon of plume contraction is observed at long times of solute propagation during the drying stage. The mean concentration profiles match those of the Monte-Carlo simulations for capillary length scales typical of sandy soils.     

Predicting natural attenuation of xylene in groundwater using a numerical model.

The aquifer beneath an abandoned refinery in the Lower Rhine area, Germany, was contaminated with a number of different mineral oil products. Groundwater sampling in the area around the former xylene plant revealed that a xylene plume had developed in the underlying groundwater, and moreover, that there is strong evidence for in situ microbial xylene degradation with oxygen, nitrate, sulfate and ferric iron as electron acceptors. In order to prevent further xylene spreading, three pumping wells extracting contaminated water were installed downgradient of the spill zone. The numerical reactive transport code Transport Biochemisty Chemistry (TBC) was applied to this situation to quantify the relation of microbial degradation to xylene removal by the pumping wells. It could be shown that the unamended in situ degradation was an appreciable xylene removal process that contributed to about one-third to the total xylene removal (degradation plus extraction). A further objective of the model application was to predict xylene spreading under regional flow conditions, i.e. without operation of the three pumping wells, to consider the possible effects of natural xylene attenuation. To accomplish this, the model calibrated for the situation with operating wells was transferred to the hydraulic situation of regional flow while retaining the parameters of the biochemical model. It turned out that the xylene plume that is expected to develop downgradient of the source area will be limited to an extension of not more than 1000 m. An interesting feature of the simulations results was that xylene degradation under iron-reducing conditions, which was of minor importance for the situation with operating pumping wells, becomes the dominant degradation mechanism under regional flow conditions. Moreover, iron reduction will be the key process in controlling plume evolution. The model application illustrates that multi-species reactive transport models are needed to adequately transfer reactive processes from one hydraulic situation to another, while single species models are not suited for this predictive task.