Development of Artificial Neural Network Based Metamodels for Inactivation of Anthrax Spores in Ventilated Spaces Using Computational Fluid Dynamics

ABSTRACT

Linear, quadratic, and artificial neural network (ANN)-based metamodels were developed for predicting the extent of anthrax spore inactivation by chlorine dioxide in a ventilated three-dimensional space over time from computational fluid dynamics model (CFD) simulation data. Dimensionless groups were developed to define the design space of the problem scenario. The Hammersley sequence sampling (HSS) method was used to determine the sampling points for the numerical experiments within the design space. A CFD model, comprised of multiple submodels, was applied to conduct the numerical experiments. Large eddy simulation (LES) with the Smagorinsky subgrid-scale model was applied to compute the airflow. Anthrax spores were modeled as a dispersed solid phase using the Lagrangian treatment. The disinfectant transport was calculated by solving a mass transport equation. Kinetic decay constants were included for spontaneous decay of the disinfectant and for the reaction of the disinfectant with the surfaces of the three-dimensional space. To enhance the mixing of the disinfectant with the room air, a momentum source was included in the simulation. An inactivation rate equation accounted for the reaction between the spores and the disinfectant. The ANN-based metamodels were most successful in predicting the number of viable bioaerosols remaining in an arbitrary enclosed space. Sensitivity analysis showed that the mass fraction of the disinfectant, inactivation rate constant, and contact time had the most influence on the inactivation of the spores.

IMPLICATIONS This investigation presents a framework for the development of user-friendly models; metamodels for the prediction of the number of viable spores remaining in an indoor room during disinfection from accurate but time-consuming CFD studies. During any decontamination event, to know when to stop pumping in the disinfectant and to know what level of log reduction of the spores have been achieved before even starting decontamination would provide valuable guidance. The neural network based metamodels can be applied to obtain quick and relatively accurate answers. This would be necessary when immediate information is required during emergencies.     

Computational fluid dynamics analysis of the effects of reactor configuration on disinfection efficiency.

The efficacy of disinfection processes in water purification systems is governed by several key factors, including reactor hydraulics, disinfectant chemistry, and microbial inactivation kinetics. The objective of this work was to develop a computational fluid dynamics (CFD) model to predict velocity fields, mass transport, chlorine decay, and microbial inactivation in a continuous flow reactor. The CFD model was also used to evaluate disinfection efficiency in alternative reactor designs. The CFD reactor analysis demonstrates that disinfection efficiency is affected by both kinetics and mixing state (i.e., degree of micromixing or segregation). Residence time distributions (RTDs) derived from tracer analysis do not describe intrinsic mixing conditions. The CFD-based disinfection models account for reactor mixing patterns by resolution of the reactor velocity field and thus provide a better prediction of microbial inactivation than models that use an RTD.     

Application of gaseous ozone for inactivation of Bacillus subtilis spores

The effectiveness of gaseous ozone (O3) as a disinfectant was tested on Bacillus subtilis spores, which share the same physiological characteristics as Bacillus anthracis spores that cause the anthrax disease. Spores dried on surfaces of different carrier material were exposed to O3 gas in the range of 500-5000 ppm and at relative humidity (RH) of 70-95%. Gaseous O3 was found to be very effective against the B. subtilis spores, and at O3 concentrations as low as 3 mg/L (1500 ppm), approximately 3-log inactivation was obtained within 4 hr of exposure. The inactivation curves consisted of a short lag phase followed by an exponential decrease in the number of surviving spores. Prehydration of the bacterial spores has eliminated the initial lag phase. The inactivation rate increased with increasing O3 concentration but not >3 mg/L. The inactivation rate also increased with increase in RH. Different survival curves were obtained for various surfaces used to carry spores. Inactivation rates of spores on glass, a vinyl floor tile, and office paper were nearly the same. Whereas cut pile carpet and hardwood flooring surfaces resulted in much lower inactivation rates, another type of carpet (loop pile) showed significant enhancement in the inactivation of the spores.     

One-dimensional simulation of solute transfer in saturated-unsaturated porous media using the discontinuous finite elements method

A one-dimensional transport model for simulating water flow and solute transport in homogeneous-heterogeneous, saturated-unsaturated porous media is presented. The model is composed of a combination of accurate numerical algorithms for solving the nonlinear Richard's and advection-dispersion equations (ADE). The mixed form of Richard's equation is solved using a standard finite element method (FEM) with primary variable switching. The transport equation is solved using operator splitting, with the discontinuous finite element method (DFE) for discretization of the advective term. A slope limiting procedure for DFE avoids numerical instabilities but creates very limited numerical dispersion for high Peclet numbers. An implicit finite differences scheme (FD) is used for the dispersive term. The unsaturated flow and transport model (Wamos-T) is applied to a variety of rigorous problems including transient flow, heterogeneous medium and abrupt variations of velocity in magnitude and direction due to time-varying boundary conditions. It produces accurate and mass-conservative solutions for a very large range of grid Peclet numbers. The Wamos-T model is a good and robust alternative for the simulation of mass transport in unsaturated domain.     

Effect of scale and dimensionality on the surfactant-enhanced solubilization of a residual DNAPL contamination

The mass transfer rate from residual dense non-aqueous phase liquids (DNAPLs) to the mobile aqueous phase is an important parameter for the efficiency of surfactant-enhanced remediation through solubilization of this type of contamination. The mass transfer kinetics are highly dependent on the dimensionality of the system. In this study, irregularly shaped residual TCE saturations in two-dimensional saturated flow fields were flushed with a 2% polyoxyethylene sorbitan (20) monooleate (POESMO) solution until complete removal had been achieved. A numerical model was developed and used for the simulation of the various surfactant-flushing experiments with different initial saturation patterns and flow rates. Through optimization against in situ concentration and saturation data, a phenomenological power-law model for the relationship between the mass transfer rate from the DNAPL to the mobile aqueous phase on the one hand and the residual DNAPL saturation and the flow velocity on the other hand was derived. The obtained mass transfer rate parameters provide a reasonable fit to the experimental data, predicting the cleanup time and the general saturation and concentration pattern quite well but failing to predict the concentration curves at every individual sampling port. The obtained mass transfer rate model gives smaller values for the predicted mass transfer rate but shows a comparable dependence on water flow and saturation as in earlier published one-dimensional column experiments with identical characteristics for porous medium, DNAPL and surfactant. Mass transfer rate predictions were about one order of magnitude lower in the 2-D flow cell experiment than in 1-D column experiments. These results give an indication for the importance of dimensionality during surfactant remediation.     

One-dimensional simulation of solute transfer in saturated–unsaturated porous media using the discontinuous finite elements method

A one-dimensional transport model for simulating water flow and solute transport in homogeneous–heterogeneous, saturated–unsaturated porous media is presented. The model is composed of a combination of accurate numerical algorithms for solving the nonlinear Richard's and advection–dispersion equations (ADE). The mixed form of Richard's equation is solved using a standard finite element method (FEM) with primary variable switching. The transport equation is solved using operator splitting, with the discontinuous finite element method (DFE) for discretization of the advective term. A slope limiting procedure for DFE avoids numerical instabilities but creates very limited numerical dispersion for high Peclet numbers. An implicit finite differences scheme (FD) is used for the dispersive term.The unsaturated flow and transport model (Wamos-T) is applied to a variety of rigorous problems including transient flow, heterogeneous medium and abrupt variations of velocity in magnitude and direction due to time-varying boundary conditions. It produces accurate and mass-conservative solutions for a very large range of grid Peclet numbers. The Wamos-T model is a good and robust alternative for the simulation of mass transport in unsaturated domain.     

Modeling of dispersion near roadways based on the vehicle-induced turbulence concept

A mathematical model is developed for dispersion near roadways by incorporating vehicle-induced turbulence (VIT) into Gaussian dispersion modeling using computational fluid dynamics (CFD). The model is based on the Gaussian plume equation in which roadway is regarded as a series of point sources. The Gaussian dispersion parameters are modified by simulation of the roadway using CFD in order to evaluate turbulent kinetic energy (TKE) as a measure of VIT. The model was evaluated against experimental carbon monoxide concentrations downwind of two major freeways reported in the literature. Good agreements were achieved between model results and the literature data. A significant difference was observed between the model results with and without considering VIT. The difference is rather high for data very close to the freeways. This model, after evaluation with additional data, may be used as a framework for predicting dispersion and deposition from any roadway for different traffic (vehicle type and speed) conditions.     

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.     

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.     

Development of a scalable model for predicting arsenic transport coupled with oxidation and adsorption reactions

Understanding the fundamentals of arsenic adsorption and oxidation reactions is critical for predicting its transport dynamics in groundwater systems. We completed batch experiments to study the interactions of arsenic with a common MnO2(s) mineral, pyrolusite. The reaction kinetics and adsorption isotherm developed from the batch experiments were integrated into a scalable reactive transport model to facilitate column-scale transport predictions. We then completed a set of column experiments to test the predictive capability of the reactive transport model. Our batch results indicated that the commonly used pseudo-first order kinetics for As(III) oxidation reaction neglects the scaling effects with respect to the MnO2(s) concentration. A second order kinetic equation that explicitly includes MnO2(s) concentration dependence is a more appropriate kinetic model to describe arsenic oxidation by MnO2(s) minerals. The arsenic adsorption reaction follows the Langmuir isotherm with the adsorption capacity of 0.053micromol of As(V)/g of MnO2(s) at the tested conditions. The knowledge gained from the batch experiments was used to develop a conceptual model for describing arsenic reactive transport at a column scale. The proposed conceptual model was integrated within a reactive transport code that accurately predicted the breakthrough profiles observed in multiple column experiments. The kinetic and adsorption process details obtained from the batch experiments were valuable data for scaling to predict the column-scale reactive transport of arsenic in MnO2(s)-containing sand columns.     

水质工程用紫外光消毒器的CFD数值模拟

紫外光消毒法是一种重要的污水深度处理方法,消毒器内部辐射强度的分布情况、微生物的停留时间及其运动轨迹对灭菌效果起着至关重要的作用。本文首次在国内系统介绍了采用计算流体动力学(CFD)手段对紫外光消毒器进行数值模拟研究的理论基础和技术路线,指出辐射强度的计算和停留时间的确定是其中的关键。以特定结构的紫外光消毒器为对象,对其内部流场进行了模拟计算,讨论了挡板位置不同所引起的消毒器内部速度场的变化情况;应用离散坐标辐射模型(DO)对紫外光辐射强度分布进行了模拟,同时加入离散粒子模型(DPM)来模拟消毒器内微生物的停留时长和运动轨迹,为最终计算微生物所受紫外光辐射剂量奠定了坚实的基础。     

Transport of coliphage PRD1 in a surface flow constructed wetland.

A tracer study was conducted in a 3-ha surface flow constructed wetland to analyze transport performance of PRD1, an enteric virus model. The convection-dispersion equation (CDE), including a first-order reaction model, adequately simulated transport performance of PRD1 in the wetland under an average hydraulic loading rate of 82 mm/d. Convective velocity (v) and longitudinal dispersion coefficient (D) were estimated by modeling a conservative tracer (bromide) pulse through the wetland. Both PRD1 and bromide were simultaneously added to the entering secondary treated wastewater effluent. The mass of bromide and PRD1 recovered was 76 and 16%, respectively. The PRD1 decay rate was calculated to be 0.3/day. The findings of this study suggest that the CDE model and analytical moment equations represent a suitable option to characterize virus transport performance in surface flow constructed wetlands.     

Peracetic acid disinfection: a feasible alternative to wastewater chlorination.

The paper summarizes the results of a bench-scale study to evaluate the feasibility of using peracetic acid (PAA) as a substitute for sodium hypochlorite both for discharge into surface water and for agricultural reuse. Trials were carried out with increasing doses (1, 2, 3, 5, 10, and 15 mg/L) and contact times (6, 12, 18, 36, 42, and 54 minutes) to study disinfectant decay and bacterial removal and regrowth, using fecal coliform and Escherichia coli (E. coli) as process efficiency indicators. Peracetic acid decay kinetics was evaluated in tap water and wastewater; in both cases, PAA decays according to first-order kinetics with respect to time, and a correlation was found between PAA oxidative initial consumption and wastewater characteristics. The PAA disinfection efficiency was correlated with operating parameters (active concentration and contact time), testing different kinetic models. Two data groups displaying a different behavior on the basis of initial active concentration ranges (1 to 2 mg/L and 5 to 15 mg/L, respectively) can be outlined. Both groups had a "tailing-off" inactivation curve with respect to time, but the second one showed a greater inactivation rate. Moreover, the effect of contact time was greater at the lower doses. Hom's model, used separately for the two data groups, was found to best fit experimental data, and the disinfectant active concentration appears to be the main factor affecting log-survival ratios. Moreover, the S-model better explains the initial resistance of E. coli, especially at low active concentrations (< 2 mg/L) and short contact times (< 12 minutes). Microbial counts, performed by both traditional methods and flow cytometry, immediately and 5 hours after sample collection (both with or without residual PAA inactivation), showed that no appreciable regrowth took place after 5 hours, neither for coliform group bacteria, nor for total heterotrophic bacteria.     

Radial diffusive sampler for the determination of 8-h ambient ozone concentrations

The 8-h ozone radial diffusive sampler was evaluated according to the CEN protocol for the validation of diffusive samplers. All the parameters regarding the sampler characteristics were found to be consistent with the requirements of this protocol apart from the blank value, which must be evaluated and subtracted at each sampling. The nominal uptake rate was determined in laboratory conditions. However, the uptake rate depends on the mass uptake, temperature, humidity and on the combination of temperature and humidity. Based on laboratory experiments, an empirical model has been established which improved the agreement between the radial sampler and the reference method. This improvement was observed under several different meteorological and emission conditions of sampling. By using the model equation of uptake rate, the data quality objective of 30% for the expanded uncertainty included in the O(3) European Directive, is easily attained. Therefore, the sampler represents an appropriate indicative method.     

Response surface modeling for the inactivation of Bacillus subtilis subsp. niger spores by chlorine dioxide gas in an enclosed space

Bacillus subtilis subsp. niger spores are a commonly used biological indicator to evaluate the disinfection of an enclosed space. In the present study, chlorine dioxide (ClO₂) gas was applied to inactivate B. subtilis subsp. niger spores in an enclosed space. The effects of the ClO₂ gas concentration (1-3 mg/l), relative humidity (RH, 30-70%) and exposure time (30-90 min) were investigated using a response surface methodology (RSM). A three-factor Box-Behnken experimental design was used. The obtained data were adequately fitted to a second-order polynomial model with an R_2adj of 0.992. The ClO₂ gas concentration, RH and exposure time all significantly (P<0.05) and positively correlated with the inactivation of B. subtilis subsp. niger spores. The interaction between the ClO₂ gas concentration and RH as well as that between the exposure time and RH indicated significant and synergistic effects (P<0.05). The predictive model was validated by additional eight experiments and proven to be with good accuracy. Overall, this model established by the RSM could show the trend of the inactivation of spores, indicate the interactions between important factors, and provide a reference to determine effective conditions for the disinfection in different enclosed spaces by ClO₂ gas.

Implications: The inactivation of indoor biological contaminants plays an important role in preventing the transmission of pathogens and ensuring human safety. The predictive model using response surface methodology indicates the influence and interaction of the main factors on the inactivation of Bacillus subtilis subsp. niger spores by ClO₂ gas, and can predict a ClO₂ gas treatment condition to achieve an effective sterilization of enclosed spaces. The results in this paper will provide a reference for the application of ClO₂ gas treatments for indoor disinfection.     

Laboratory analog and numerical study of groundwater flow and solute transport in a karst aquifer with conduit and matrix domains

New mathematical and laboratory methods have been developed for simulating groundwater flow and solute transport in karst aquifers having conduits imbedded in a porous medium, such as limestone. The Stokes equations are used to model the flow in the conduits and the Darcy equation is used for the flow in the matrix. The Beavers–Joseph interface boundary conditions are adopted to describe the flow exchange at the interface boundary between the two domains. A laboratory analog is used to simulate the conduit and matrix domains of a karst aquifer. The conduit domain is located at the bottom of the transparent plexiglas laboratory analog and glass beads occupy the remaining space to represent the matrix domain. Water flows into and out of the two domains separately and each has its own supply and outflow reservoirs. Water and solute are exchanged through an interface between the two domains. Pressure transducers located within the matrix and conduit domains of the analog provide data that is processed and stored in digital format. Dye tracing experiments are recorded using time-lapse imaging. The data and images produced are analyzed by a spatial analysis program. The experiments provide not only hydraulic head distribution but also capture solute front images and mass exchange measurements between the conduit and matrix domains. In the experiment, we measure and record pressures, and quantify flow rates and solute transport. The results present a plausible argument that laboratory analogs can characterize groundwater water flow, solute transport, and mass exchange between the conduit and matrix domains in a karst aquifer. The analog validates the predictions of a numerical model and demonstrates the need of laboratory analogs to provide verification of proposed theories and the calibration of mathematical models.     

Kinetics of inorganic carbon utilization by microalgal biofilm in a flat plate photoreactor

A kinetic model was developed to describe inorganic carbon utilization by microalgae biofilm in a flat plate photoreactor. The model incorporates the fundamental mechanisms of diffusive mass transport and biological reaction of inorganic carbon by microalgal biofilm. An advanced numerical technique, the orthogonal collocation method and Gear's method, was employed to solve this kinetic model. The model solutions included the concentration profiles of inorganic carbon in the microalgal biofilm, the growths of suspended microalgae and microalgal biofilm, the effluent concentrations of inorganic carbon, and the flux of inorganic carbon from bulk liquid into biofilm. The batch kinetic test was independently conducted to determine biokinetic parameters used in the microalgal biofilm model simulation while initial thickness of microalgal biofilm were assumed. A laboratory-scale flat plate photoreactor with a high recycle flow rate was setup and conducted to verify the model. The volume of photoreactor is 60 l which yields a hydraulic retention time of 1.67 days. The model-generated inorganic carbon and the suspended microalgae concentration curves agreed well with those obtained in the laboratory-scale test. The fixation efficiencies of HCO(3)(-) and CO(2) are 98.5% and 90% at a steady-state condition, respectively. The concentration of suspended microalgal cell reached up to 12 mg/l at a maximum growth rate while the thickness of microalgal biofilm was estimated to be 104 microm at a steady-state condition. The approaches of experiments and model simulation presented in this study could be employed for the design of a flat plate photoreactor to treat CO(2) by microalgal biofilm in a fossil-fuel power plant.     

Evidence of one-dimensional scale-dependent fractional advection-dispersion

A semi-analytical inverse method and the corresponding program FADEMain for parameter estimation of the fractional advection-dispersion equation (FADE) were developed in this paper. We have analyzed Huang et al.'s [Huang, K., Toride, N., van Genuchten, M.Th., 1995. Experimental investigation of solute transport in large homogeneous and heterogeneous saturated soil columns. Trans. Porous Media 18, 283-302.] laboratory experimental data of conservative solute transport in 12.5-m long homogeneous and heterogeneous soil columns to test the non-Fickian dispersion theory of FADE. The dispersion coefficient was calculated by fitting the analytical solution of FADE to the measured data at different transport scales. We found that the dispersion coefficient increased exponentially with transport scale for the homogeneous column, whereas it increased with transport scale in a power law function for the heterogeneous column. The scale effect of the dispersion coefficient in the heterogeneous soil was much more significant comparing to that in the homogeneous soil. The increasing rate of dispersion coefficient versus transport distance was smaller for FADE than that for the advection-dispersion equation (ADE). Finite difference numerical approximations of the scale-dependent FADE were established to interpret the experimental results. The numerical solutions were found to be adequate for predicting scale-dependent transport in the homogeneous column, while the prediction for the heterogeneous column was less satisfactory.     

Wastewater disinfection by peracetic acid: assessment of models for tracking residual measurements and inactivation.

With its potential for low (if any) disinfection byproduct formation and easy retrofit for chlorine contactors, peracetic acid (PAA) or use of PAA in combination with other disinfectant technologies may be an attractive alternative to chlorine-based disinfection. Examples of systems that might benefit from use of PAA are water reuse schemes or plants discharging to sensitive receiving water bodies. Though PAA is in use in numerous wastewater treatment plants in Europe, its chemical kinetics, microbial inactivation rates, and mode of action against microorganisms are not thoroughly understood. This paper presents results from experimental studies of PAA demand, PAA decay, and microbial inactivation, with a complementary modeling analysis. Model results are used to evaluate techniques for measurement of PAA concentration and to develop hypotheses regarding the mode of action of PAA in bacterial inactivation. Kinetic and microbial inactivation rate data were collected for typical wastewaters and may be useful for engineers in evaluating whether to convert from chlorine to PAA disinfection.     

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.