Transport of Escherichia coli and solutes during waste water infiltration in an urban alluvial aquifer

Recharge of waste water in an unconsolidated poorly sorted alluvial aquifer is a complex process, both physically and hydrochemically. The aim of this paper is to analyse and conceptualise vertical transport mechanisms taking place in an urban area of extensive wastewater infiltration by analysing and combining the water balance, the microbial (Escherichia coli) mass balance, and the mass balance for dissolved solutes. For this, data on sediment characteristics (grain size, organic carbon, reactive iron, and calcite), groundwater levels, and concentrations of E. coli in groundwater and waste water were collected. In the laboratory, data on E. coli decay rate coefficients, and on bacteria retention characteristics of the sediment were collected via column experiments. The results indicated that shallow groundwater, at depths of 50 m below the surface, was contaminated with E. coli concentrations as high as 10(6) CFU/100 mL. In general, E. coli concentrations decreased only 3 log units from the point of infiltration to shallow groundwater. Concentrations were lower at greater depths in the aquifer. In laboratory columns of disturbed sediments, bacteria removal was 2-5 log units/0.5 cm column sediment. Because of the relatively high E. coli concentrations in the shallow aquifer, transport had likely taken place via a connected network of pores with a diameter large enough to allow bacterial transport instead of via the sediment matrix, which was inaccessible for bacteria, as was clear from the column experiments. The decay rate coefficient was determined from laboratory microcosms to be 0.15 d(-1). Assuming that decay in the aquifer was similar to decay in the laboratory, then the pore water flow velocity between the point of infiltration and shallow groundwater, coinciding with a concentration decrease of 3 log units, was 0.38 m/d, and therefore, transport in this connected network of pores was fast. According to the water balance of the alluvial aquifer, determined from transient groundwater modelling, groundwater flow in the aquifer was mainly in vertical downward direction, and therefore, the mass balance for dissolved solutes was simulated using a 1D transport model of a 200 m column of the Quaternary Alluvium aquifer. The model, constructed with PHREEQC, included dual porosity, and was able to adequately simulate removal of E. coli, cation-exchange, and nitrification. The added value of the use of E. coli in this study was the recognition of relatively fast transport velocities occurring in the aquifer, and the necessity to use the dual porosity concept to investigate vertical transport mechanisms. Therefore, in general and if possible, microbial mass balances should be considered more systematically as an integral part of transport studies.     

Concentration statistics for mixing-controlled reactive transport in random heterogeneous media

Uncertainty in the distribution of hydraulic parameters leads to uncertainty in flow and reactive transport. Traditional stochastic analysis of solute transport in heterogeneous media has focused on the ensemble mean of conservative-tracer concentration. Studies in the past years have shown that the mean concentration often is associated with a high variance. Because the range of possible concentration values is bounded, a high variance implies high probability weights on the extreme values. In certain cases of mixing-controlled reactive transport, concentrations of conservative tracers, denoted mixing ratios, can be mapped to those of constituents that react with each other upon mixing. This facilitates mapping entire statistical distributions from mixing ratios to reactive-constituent concentrations. In perturbative approximations, only the mean and variance of the mixing-ratio distribution are used. We demonstrate that the second-order perturbative approximation leads to erroneous or even physically impossible estimates of mean reactive-constituent concentrations when the variance of the mixing ratio is high and the relationship between the mixing ratio and the reactive-constituent concentrations strongly deviates from a quadratic function. The latter might be the case in biokinetic reactions or in equilibrium reactions with small equilibrium constant in comparison to the range of reactive-constituent concentrations. When only the mean and variance of the mixing ratio is known, we recommend assuming a distribution that meets the known bounds of the mixing ratio, such as the beta distribution, and mapping the assumed distribution of the mixing ratio to the distributions of the reactive constituents.     

Upscaling heterogeneity in aquifer reactivity via the exposure-time concept: inverse model

A novel inverse technique is proposed to quantitatively characterize macroscopic variability in aquifer reactivity in a Lagrangian representation. Reactivity heterogeneity is expressed in terms of distributions of flux over cumulative time of exposure of the solution to reactive surface area, termed here 'cumulative reactivity'. In cases involving single aqueous species the combined effects of physical and reactivity heterogeneity on reactive solute transport can often be established and further investigated through joint distributions of flux over travel time and cumulative reactivity. The inverse technique requires the breakthrough curve of a passive tracer to determine the distribution of flux over travel time, and additional breakthrough curves of reactive tracers provide additional moments of the distribution of flux over cumulative reactivity given travel time. Thus breakthroughs of one passive and two reactive tracers can provide the mean and variance of the distribution of flux over cumulative reactivity. This Lagrangian characterization is achieved with knowledge of the types of reactive surfaces present, but not their spatial locations. The distributions can subsequently be applied via forward modeling using the same technique to predict breakthrough curves of other solutes undergoing first-order reactions in similar physically and chemically heterogeneous configurations.     

Simulation of solute transport across low-permeability barrier walls

Low-permeability, non-reactive barrier walls are often used to contain contaminants in an aquifer. Rates of solute transport through such barriers are typically many orders of magnitude slower than rates through the aquifer. Nevertheless, the success of remedial actions may be sensitive to these low rates of transport. Two numerical simulation methods for representing low-permeability barriers in a finite-difference groundwater-flow and transport model were tested. In the first method, the hydraulic properties of the barrier were represented directly on grid cells and in the second method, the intercell hydraulic-conductance values were adjusted to approximate the reduction in horizontal flow, allowing use of a coarser and computationally efficient grid. The alternative methods were tested and evaluated on the basis of hypothetical test problems and a field case involving tetrachloroethylene (PCE) contamination at a Superfund site in New Hampshire. For all cases, advective transport across the barrier was negligible, but preexisting numerical approaches to calculate dispersion yielded dispersive fluxes that were greater than expected. A transport model (MODFLOW-GWT) was modified to (1) allow different dispersive and diffusive properties to be assigned to the barrier than the adjacent aquifer and (2) more accurately calculate dispersion from concentration gradients and solute fluxes near barriers. The new approach yields reasonable and accurate concentrations for the test cases.     

Modelling of diffusion-limited retardation of contaminants in hydraulically and lithologically nonuniform media

A new reactive transport modelling approach and examples of its application are presented, dealing with the impact of sorption/desorption kinetics on the spreading of solutes, e.g. organic contaminants, in groundwater. Slow sorption/desorption is known from the literature to be strongly responsible for the retardation of organic contaminants. The modelling concept applied in this paper quantifies sorption/desorption kinetics by an intra-particle diffusion approach. According to this idea, solute uptake by or release from the aquifer material is modelled at small scale by a "slow" diffusion process where the diffusion coefficient is reduced as compared to the aqueous diffusion coefficient due to (i) the size and shape of intra-particle pores and (ii) retarded transport of solutes within intra-particle pores governed by a nonlinear sorption isotherm. This process-based concept has the advantage of requiring only measurable model parameters, thus avoiding fitting parameters like first-order rate coefficients.In addition, the approach presented here allows for modelling of slow sorption/desorption in lithologically nonuniform media. Therefore, it accounts for well-known experimental findings indicating that sorptive properties depend on (i) the grain size distribution of the aquifer material and (ii) the lithological composition (e.g. percentage of quartz, sandstone, limestone, etc.) of each grain size fraction. The small-scale physico-chemical model describing sorption/desorption is coupled to a large-scale model of groundwater flow and solute transport. Consequently, hydraulic heterogeneities may also be considered by the overall model. This coupling is regarded as an essential prerequisite for simulating field-scale scenarios which will be addressed by a forthcoming publication.This paper focuses on mathematical model formulation, implementation of the numerical code and lab-scale model applications highlighting the sorption and desorption behavior of an organic contaminant (Phenanthrene) with regard to three lithocomponents exhibiting different sorptive properties. In particular, it is shown that breakthrough curves (BTCs) for lithologically nonuniform media cannot be obtained via simple arithmetic averaging of breakthrough curves for lithologically uniform media. In addition, as no analytical solutions are available for model validation purposes, simulation results are compared to measurements from lab-scale column experiments. The model results indicate that the new code can be regarded as a valuable tool for predicting long-term contaminant uptake or release, which may last for several hundreds of years for some lithocomponents. In particular, breakthrough curves simulated by pure forward modelling reproduce experimental data much better than a calibrated standard first-order kinetics reactive transport model, thus indicating that the new approach is of high quality and may be advantageously used for supporting the design of remediation strategies at contaminated sites where some lithocomponents and/or grain size classes may provide a long-term pollutant source.     

Variability of aquifer sorption properties in a field experiment on groundwater transport of organic solutes: Methods and preliminary results

Methods to characterize the organic solute sorption distribution coefficient, organic carbon content, and specific surface area of aquifer solids from the site of a field experiment on solute transport in groundwater were refined for application to small subsamples of 10-cm depth increments taken from 5-cm diameter cores. Initial results indicate that the average sorption characteristics of the Borden aquifer do not vary appreciably along the trajectory of the solute plumes. However, the sorption distribution coefficient of tetrachloroethylene varied over nearly an order of magnitude among 10-cm depth increments in one core sample. Preliminary evidence suggests that the sorption distribution coefficients for four halogenated organic solutes vary proportionally among core strata. However, the distribution coefficients for sorption of tetrachloroethylene on various depth increments are not well correlated with either organic carbon content or specific surface area, suggesting that as yet unidentified mineral phases may play a significant role in sorption of such solutes by the sandy aquifer solids.     

Mineralogical compositions of aquifer matrix as necessary initial conditions in reactive contaminant transport models

Mineralogical compositions and their spatial distributions are important initial conditions for reactive transport modeling. However, popular Kd-based "reactive" transport models only require contaminant concentrations in the pore fluids as initial conditions, and minerals implicitly represent infinite sources and sinks in these models. That situation results in a general neglect of mineralogical characterization in site investigations. This study uses a coupled multi-component reactive mass transport model to predict the natural attenuation of a ground water plume at a uranium mill tailings site in western USA. Numerous ground water geochemistry data are available at this site, but mineralogical data are sketchy. Even given the well-defined pore fluid chemistry, variations of secondary mineral species and mineral abundances in the aquifer resulted in significantly different modeling outcomes. Results show that the amount of calcite in the aquifer determines the distances of plume migration. The possible presence of jurbanite, an aluminum sulfate phase, can store acidity temporarily but cause more severe contamination on a later date. The surfaces of iron oxyhydroxides can store significant amounts of sulfate and protons and serve as a second source for prolonged contamination. These simulations under field conditions illustrate that mineralogical compositions are an essential requirement for accurate prediction of contaminant fate and transport.     

Analytical solution to transport in three-dimensional heterogeneous well capture zones

Solute transport is investigated in a heterogeneous aquifer for combined natural-gradient and well flows. The heterogeneity is associated with the spatially varying hydraulic conductivity K(x, y, z), which is modelled as a log-normal stationary-random function. As such, the conductivity distribution is characterized by four parameters: the arithmetic mean K(A), the variance sigma(Y)(2) (Y=lnK), the horizontal integral scale I of the axisymmetric log-conductivity autocorrelation and the anisotropy ratio e=I(v)/I (I(v) is the vertical integral scale). The well fully penetrates an aquifer of constant thickness B and has given constant discharge QB, while the background aquifer flow is driven by an uniform mean head-gradient, - J. Therefore, for a medium of homogeneous conductivity K(A), the steady-state capture zone has a width 2L=Q/(K(A)|J|) far from the well (herein the term capture zone is used to refer both to the zone from which water is captured by a pumping well and the zone that captures fluid from an injecting well). The main aim is to determine the mean concentration as a function of time in fluid recovered by a pumping well or in a control volume of the aquifer that captures fluid from an injecting well. Relatively simple solutions to these complex problems are achieved by adopting a few assumptions: a thick aquifer B>I(v) of large horizontal extent (so that boundary effects may be neglected), weak heterogeneity sigma(Y)(2)<1, a highly anisotropic formation e<0.2 and neglect of pore-scale dispersion. Transport is analyzed to the first-order in sigma(Y)(2) in terms of the travel time of particles moving from or towards the well along the steady streamlines within the capture zone. Travel-time mean and variance to any point are computed by two quadratures for an exponential log-conductivity two-point covariance. Spreading is reflected by the variance value, which increases with sigma(Y)(2) and I/L. For illustration, the procedure is applied to two particular cases. In the first one, a well continuously injects water at constant concentration. The mean concentration as function of time for different values of the controlling parameters sigma(Y)(2) and I/L is determined within control volumes surrounding the well or in piezometers. In the second case, a solute plume, initially occupying a finite volume Omega(0), is drawn towards a pumping well. The expected solute recovery by the well as a function of time is determined in terms of the previous controlling parameters as well as the location and extent of Omega(0). The methodology is tested against a full three-dimensional simulation of a multi-well forced-gradient flow field test ([Lemke, L., W.B. II, Abriola, L., Goovaerts, P., 2004. Matching solute breakthrough with deterministic and stochastic aquifer models. Ground Water 42], SGS simulations). Although the flow and transport conditions are more complex than the ones pertinent to capture zones in uniform background flow, it was found that after proper adaptation the methodology led to results for the breakthrough curve in good agreement with a full three-dimensional simulation of flow and transport.     

Predictive modelling of dispersion controlled reactive plumes at the laboratory-scale

A model-based interpretation of laboratory-scale experimental data is presented. Hydrolysis experiments carried out using thin glass tanks filled with glass beads to construct a hypothetical and inert, homogeneous porous medium were analysed using a 2D numerical model. A new empirical formula, based upon results for non-reactive (tracer) experiments is used to calculate transversal dispersivity values for a range of grain sizes and any flow velocities. Combined with effective diffusion coefficients calculated from Stokes-Einstein type equations, plume lengths arising from mixing between two solutes can be predicted accurately using numerical modelling techniques. Moreover, pH and ion concentration profiles lateral to the direction of flow of the mixing species can be determined at any given point downstream, without the need for result fitting. In our case, this approach does not lead to overpredictions of lateral mixing, as previously reported when using parameters derived from non-reactive tracer experiments to describe reactive solute transport. The theory is based on the assumption of medium homogeneity.     

A controlled field experiment on groundwater contamination by a multicomponent DNAPL: dissolved-plume retardation

A natural gradient emplaced-source (ES) controlled field experiment was conducted at the Borden aquifer research site, Ontario, to study the transport of dissolved plumes emanating from residual dense nonaqueous-phase liquid (DNAPL) source zones. The specific objective of the work presented here is to determine the effects of solute and co-solute concentrations on sorption and retardation of dissolved chlorinated solvent-contaminant plumes. The ES field experiment comprised a controlled emplacement of a residual multicomponent DNAPL below the groundwater table and intensive monitoring of dissolved-phase plumes of trichloromethane (TCM), trichloroethylene (TCE), and perchloroethylene (PCE) plumes continuously generated in the aquifer down gradient from gradual source dissolution. Estimates of plume retardation (and dispersion) were obtained from 3-D numerical simulations that incorporated transient source input and flow regimes monitored during the test. PCE, the most retarded solute, surprisingly exhibited a retardation factor approximately 3 times lower than observed in a previous Borden tracer test by Mackay et al. [Water Resour. Res. 22 (1986) 2017] conducted approximately 150 m away. Also, an absence of temporal trend in PCE retardation contrasted with the previous Borden test. Supporting laboratory studies on ES site core indicated that sorption was nonlinear and competitive, i.e. reduced sorption of PCE was observed in the presence of TCE. Consideration of the effects of relatively high co-solute (TCE) concentration (competitive sorption) in addition to PCE concentration effects (nonlinear sorption) was necessary to yield laboratory-based PCE retardation estimates consistent with the field plume values. Concentration- and co-solute-based sorption and retardation analysis was also applied to the previous low-concentration pulse injection test of Mackay et al. [Water Resour. Res. 22 (1986) 2017] and was able to successfully predict the temporal field retardation trends observed in that test. While it is acknowledged that other "nonideal transport" effects may contribute, our analysis predicts differences in the PCE retardation magnitude and trend between the two experiments that are consistent with field observations based on the marked solute concentration differences that resulted from contrasting source conditions. Solute and co-solute concentration effects have heretofore received little attention, but may have wide significance in aquifers contaminated by point-source pollutants because many plumes contain mixed solutes over wide concentration ranges in strata that are likely subject to nonlinear sorption.     

Modelling of physical and reactive processes during biodegradation of a hydrocarbon plume under transient groundwater flow conditions

Numerical experiments of non-reactive and reactive transport were carried out to quantify the influence of a seasonally varying, transient flow field on transport and natural attenuation at a hydrocarbon-contaminated field site. Different numerical schemes for solving advective transport were compared to assess their capability to model low transversal dispersivities in transient flow fields. For the field site, it is shown that vertical plume spreading is largely inhibited, particularly if sorption is taken into account. For the reactive simulations, a biodegradation reaction module for the geochemical transport model PHT3D was developed. Results of the reactive transport simulations show that under the site-specific conditions the temporal variations in groundwater flow do, to a modest extent, affect average biodegradation rates and average total (dissolved) contaminant mass in the aquifer. The model simulations demonstrate that the seasonal variability in groundwater flow only results in significantly enhanced biodegradation rates when a differential sorption of electron donor (toluene) and electron acceptor (sulfate) is assumed.     

Interpreting organic solute transport data from a field experiment using physical nonequilibrium models

In a field experiment, two inorganic tracers and five organic solutes were injected into an unconfined sand aquifer. Breakthrough response curves were obtained at several points downgradient of the injection zone. These response curves are analyzed using a model which assumes equilibrium sorption and two models which postulate physical nonequilibrium. The physical nonequilibrium models hypothesize the existence of zones of immobile water, which act as diffusion sources and sinks for the solutes. The physical nonequilibrium models better simulate the sharp breakthrough and extended tailing exhibited by the experimental responses than does the model assuming equilibrium sorption. The reasonableness of parameters obtained from curve-fitting the data is assessed. The two physical nonequilibrium models are compared.     

Pesticide transport in an aerobic aquifer with variable pH--modeling of a field scale injection experiment

Three-dimensional reactive transport simulations were undertaken to study the sorption and degradation dynamics of three herbicides in a shallow aerobic aquifer with spatially variable pH during a 216 days injection experiment. Sorption of two phenoxy acids [(+/-)-2-(4-chloro-2-methylphenoxy) propanoic acid] (MCPP) and [(+/-)-2-(2,4-dichlorophenoxy)propanoic acid] (dichlorprop) was found to be negligible. Degradation of the phenoxy acids was rapid after an initial lag phase. Degradation of the phenoxy acids could only be reproduced satisfactorily by growth-linked microbial degradation. The model fit to the field data was slightly improved if degradation was assumed to be influenced by the local pH that was observed to increase with depth ( approximately 4.5--5.7). In the observed pH-range the nitroaromatic herbicide [2-Methyl-4,6-dinitrophenol] (DNOC) was partly dissociated (pK(a)=4.31) and present in both the neutral and ionized form. The model simulations demonstrated that most of the observed spatial variation in sorption of DNOC could be explained by assuming that only the neutral form of DNOC was subject to sorption. A varying flow field was observed during the injection experiment and the model simulations documented that this most likely resulted in different migration paths for DNOC and the non-sorbing solutes. The model simulations indicated that degradation of DNOC was an important process. The degradation rate of DNOC remained constant over time and was simulated adequately by first-order kinetics. Again, the model fit to field observation was slightly improved if local pH was assumed to influence the degradation rate. Only the maximum utilization rate was estimated from the field data, while the remaining degradation parameters where successfully transferred from the laboratory study.     

An analytical solution for two-dimensional contaminant transport during groundwater extraction

We obtain an analytical solution for two-dimensional steady-state transport of conservative contaminant between injecting and pumping wells. Flow and transport are considered in the vertical cross-section. The Dupuit approximation and conformal mapping onto the complex potential domain are employed to determine the velocity and concentration distributions, respectively. We use this solution to derive a priori conditions under which widely used 1-D analytical solutions with constant velocity and dispersion coefficients provide accurate approximations. These conditions are formulated in terms of aquifer parameters, such as hydraulic conductivity, porosity and dispersivities, and remediation strategy, e.g., well spacing and pumping regimes.     

Rapid transport from the surface to wells in fractured rock: a unique infiltration tracer experiment

A unique infiltration tracer experiment was performed whereby a fluorescent dye was applied to the land surface in an agricultural field, near Perth, Ontario, Canada, to simulate the transport of solutes to two pumped monitoring wells drilled into the granitic gneiss aquifer. This experiment, interpreted using the discrete-fracture capability of the numerical model HydroGeoSphere, showed that solute transport from the surface through thin soil (less than 2m) to wells in fractured bedrock can be extremely rapid (on the order of hours). Also, it was demonstrated that maximum concentrations of contaminants originating from the ground surface will not necessarily be the highest in the shallow aquifer horizon. These are important considerations for both private and government-owned drinking water systems that draw water from shallow fractured bedrock aquifers. This research illustrates the extreme importance of protecting drinking water at the source.     

Reaction capacity characterization of shallow sedimentary deposits in geologically different regions of the Netherlands

Quantitative insight into the reaction capacity of porous media is necessary to assess the buffering capacity of the subsurface against contaminant input via groundwater recharge. Here, reaction capacity is to be considered as a series of geochemical characteristics that control acid/base conditions, redox conditions and sorption intensity. Using existing geochemical analyses, a statistical regional assessment of the reaction capacity was performed for two geologically different areas in the Netherlands. The first area is dominated by Pleistocene aquifer sediments only, in the second area a heterogeneous Holocene confining layer is found on top of the Pleistocene aquifer sediments. Within both areas, two or more regions can be distinguished that have a distinctly different geological build-up of the shallow subsurface. The reactive compounds considered were pyrite, reactive Fe other than pyrite, sedimentary organic matter, carbonate and clay content. This characterization was complemented by the analysis of a dataset of samples newly collected, from two regions within the Pleistocene area, where the sedimentary facies of samples was additionally distinguished. The statistical assessment per area was executed at the levels of region, geological formation and lithology class. For both areas, significant differences in reaction capacities were observed between: 1. different lithology classes within a geological formation in a single region, 2. identical geological formations in different regions and 3. various geological formations within a single region. Here, the reaction capacity is not only controlled by lithostratigraphy, but also by post-depositional diagenesis and paleohydrology. Correlation coefficients among the reactive compounds were generally higher for sand than for clay, but insufficiently high to allow good estimation of reactive compounds from each other. For the sandy Pleistocene aquifer sediments, the content of reactive compounds was frequently observed to be below detection limits. From this, future characterization of sediment reaction capacity is best performed at the sublevel of lithology class, being the geochemically near-uniform unit identifiable for individual geological formations within geographic regions. Additional subdivision on facies provides particular insight in the spatial entity where relatively high reaction capacities may be encountered. To obtain quantitative insight into the reaction capacity of aquifer sediments, non-sandy minor subunits should be well characterised on their reaction capacity as well as their spatial occurrence in the geological formations. A straightforward approach is presented in which the regional statistics on geochemical reactivity become combined with a 3-dimensional geological voxel model. This results into 3-dimensional data fields on reactivity, which are suitable for, for example, groundwater transport modelling. The sedimentological architecture of the deposits becomes well maintained in the geochemical data field, which is an advantage in itself.     

Migration of contaminants through the unsaturated zone overlying the Hesbaye chalky aquifer in Belgium: a field investigation

This paper presents the results of a detailed field investigation that was performed for studying groundwater recharge processes and solute downward migration mechanisms prevailing in the unsaturated zone overlying a chalk aquifer in Belgium. Various laboratory measurements were performed on core samples collected during the drilling of boreholes in the experimental site. In the field, experiments consisted of well logging, infiltration tests in the unsaturated zone, pumping tests in the saturated zone and tracer tests in both the saturated and unsaturated zones. Results show that gravitational flows govern groundwater recharge and solute migration mechanisms in the unsaturated zone. In the variably saturated chalk, the migration and retardation of solutes is strongly influenced by recharge conditions. Under intense injection conditions, solutes migrate at high speed along the partially saturated fissures, downward to the saturated zone. At the same time, they are temporarily retarded in the almost immobile water located in the chalk matrix. Under normal recharge conditions, fissures are inactive and solutes migrate slowly through the chalk matrix. Results also show that concentration dynamics in the saturated zone are related to fluctuations of groundwater levels in the aquifer. A conceptual model is proposed to explain the hydrodispersive behaviour of the variably saturated chalk. Finally, the vulnerability of the chalk to contamination issues occurring at the land surface is discussed.     

Uranium removal from groundwater via in situ biostimulation: Field-scale modeling of transport and biological processes

During 2002 and 2003, bioremediation experiments in the unconfined aquifer of the Old Rifle UMTRA field site in western Colorado provided evidence for the immobilization of hexavalent uranium in groundwater by iron-reducing Geobacter sp. stimulated by acetate amendment. As the bioavailable Fe(III) terminal electron acceptor was depleted in the zone just downgradient of the acetate injection gallery, sulfate-reducing organisms came to dominate the microbial community. In the present study, we use multicomponent reactive transport modeling to analyze data from the 2002 field experiment to identify the dominant transport and biological processes controlling uranium mobility during biostimulation, and determine field-scale parameters for these modeled processes. The coupled process simulation approach was able to establish a quantitative characterization of the principal flow, transport, and reaction processes based on the 2002 field experiment, that could be applied without modification to describe the 2003 field experiment. Insights gained from this analysis include field-scale estimates of the bioavailable Fe(III) mineral threshold for the onset of sulfate reduction, and rates for the Fe(III), U(VI), and sulfate terminal electron accepting processes.     

Prediction of mass-transfer coefficient for solute transport in porous media

Several previously reported laboratory studies related to transport of solutes through packed columns were utilized to develop predictive relationships for mass-transfer rate coefficient. The data were classified into two groups: those obtained under rate-limited mass transfer between mobile and immobile water regions (physical nonequilibrium conditions), and those derived from rate-limited mass transfer between instantaneous and slow sorption sites (sorption nonequilibrium conditions). The mass-transfer coefficient in all these studies was obtained by fitting breakthrough data to a transport model employing a first-order rate limitations with a "constant" mass-transfer coefficient, independent of flow conditions. This study demonstrated that the mass-transfer coefficient in these models is dependent on system parameters including pore-water velocity, length-scale, retardation coefficient, and particle or aggregate size. Predictive relationships were developed, through regression analysis, relating mass-transfer coefficient to residence time. The developed relationships adequately estimated previously reported field mass-transfer values. Successful simulations of field desorption data reported by Bahr [J. Contam. Hydrol. 4 (1989) 205] further demonstrate the potential applicability of the developed relationships.