Kinetic modeling of virus transport at the field scale

Bacteriophage removal by soil passage in two field studies was re-analyzed with the goal to investigate differences between one- and two-dimensional modeling approaches, differences between one- and two-site kinetic sorption models, and the role of heterogeneities in the soil properties. The first study involved removal of bacteriophages MS2 and PRDI by dune recharge, while the second study represented removal of MS2 by deep well injection. In both studies, removal was higher during the first meters of soil passage than thereafter. The software packages HYDRUS-ID and HYDRUS-2D, which simulate water flow and solute transport in one- and two-dimensional variably saturated porous media, respectively, were used. The two codes were modified by incorporating reversible adsorption to two types of kinetic sites. Tracer concentrations were used first to calibrate flow and transport parameters of both models before analyzing transport of bacteriophages. The one-dimensional one-site model did not fully describe the tails of the measured breakthrough curves of MS2 and PRD1 from the dune recharge study. While the one-dimensional one-site model predicted a sudden decrease in virus concentrations immediately after the peaks, measured data displayed much smoother decline and tailing. The one-dimensional two-site model simulated the overall behavior of the breakthrough curves very well. The two-dimensional one-site model predicted a more gradual decrease in virus concentrations after the peaks than the one-dimensional one-site model, but not as good as the one-dimensional two-site model. The dimensionality of the problem hence can partly explain the smooth decrease in concentration after peak breakthrough. The two-dimensional two-site model provided the best results. Values for k(att2) and k(det2) could not be determined at the last two of four monitoring wells, thus suggesting that either a second type of kinetic sites is present in the first few meters of dune passage and not beyond the second monitoring well, or that effects of soil heterogeneity and dimensionality of the problem overshadowed this process. Variations between single collector efficiencies were relatively small, whereas collision efficiencies varied greatly. This implies that the nonlinear removal of MS2 and PRD1 is mainly caused by variations in interactions between grain and virus surfaces rather than by physical heterogeneity of the porous medium. Similarly, a two-site model performed better than the one-site model in describing MS2 concentrations for the deep well injection study. However, the concentration data were too sparse in this study to have much confidence in the fitted parameters.     

Column experiments to study nonlinear removal of bacteriophages by passage through saturated dune sand

In a recent field study on dune recharge, bacteriophages MS2 and PRD1 were found to be removed 3 log10 over the first 2.4 m and only 5 log10 over the next 27 m. To understand the causes of this nonlinear removal, column experiments were carried out under conditions similar to the field: same recharge water, temperature (5 +/- 3 degrees C) and pore water velocity (1.5 m day(-1)). Soil samples were taken along a streamline between the recharge canal and the first monitoring well. Bacteriophage phiX174 was included for comparison. The high initial removal in the field was found not to be due to heterogeneity of phage suspensions but to soil heterogeneity. Phage removal rates correlated strongly positively with soil organic carbon content, and relatively strongly positively with silt content and the presence of ferric oxyhydroxides. Soil organic carbon content, silt content and the presence of ferric oxyhydroxides were found to decrease exponentially with travel distance. Removal rates of phiX174 were found to be 3-10 times higher than those of MS2 and PRD1 due to the lower electrostatic repulsion that the less negatively charged phiX174 experiences. It is suggested that the high initial removal in the field is due to the presence of favorable sites for attachment formed by ferric oxyhydroxides that decrease exponentially with travel distance. Similar removal rates may be found at both laboratory and field scale. However, due to local variations at field scale detailed knowledge on soil heterogeneity may be needed to enable a reliable prediction of removal.     

Interactions between viruses and goethite during saturated flow: effects of solution pH, carbonate, and phosphate

Metal oxides have great potential for controlling the fate and transport of viruses in the subsurface and water-treatment systems. The processes, however, are subject to solution chemistry. In this study, a number of column experiments were conducted to examine the effects of solution pH and anions (carbonate and phosphate) on attachment, transport, and inactivation of two bacteriophages (phiX174 and MS-2) in goethite-coated sand medium. Removal of both viruses on goethite-coated sand increased as solution pH decreased from 9.3 to 7.5, due mostly to virus inactivation. MS-2, a relatively hydrophobic virus with a lower isoelectric point (3.9), was more sensitive to the change of solution pH than phiX174, a relatively hydrophilic virus with a higher isoelectric point (6.6), in terms of their attachment and inactivation on goethite. About 90% of the MS-2 particles removed by goethite (accounting for 81% of the total input) were inactivated at pH 7.5, whereas all of the removed MS-2 particles (accounting for 10% of the total input) still remained infectious at pH 9.3. In comparison, approximately 74% of the goethite-bound phiX174 particles (accounting for 95% of the total input) lost their infectivity at pH 7.5, in contrast to a complete recovery at pH 9.3 (accounting for 65% of the total input) when the columns were eluted using a beef extract solution (pH 9.5). Presence of phosphate (20 mM H(2)PO(4)(-)) in input solution reduced virus attachment and appeared to protect the viruses from being inactivated during transport; this effect was more significant on MS-2 than on phiX174. Specifically, approximately 29% of the phiX174 particles and approximately 49% of MS-2 particles injected into the column were removed during transport. Mass recovery data showed that no phiX174 was inactivated in the presence of phosphate, whereas about 38% of the MS-2 particles attached on goethite lost their infectivity. Conversely, presence of carbonate on goethite increased virus attachment and inactivation due to contribution of additional attachment sites from protonated surface groups of the carbonate ions that were adsorbed on goethite. About 70% of the total input viruses (both phiX174 and MS-2) were removed during transport, of which 35% phiX174 and 85% MS-2 were eventually inactivated.     

A sensibility analysis of model selection in modeling the reactive transport of cesium in crushed granite

We performed a sensibility analysis of model selection in modeling the reactive transport of cesium in crushed granite through model calibration and validation. Based on some solid phase analysis data and kinetic batch experimental results, we hypothesized three two-site sorption models in the LEHGC reactive transport model to fit the breakthrough curves (BTCs) from the corresponding column experiments. The analysis of breakthrough curves shows that both the empirical two-site kinetic linear sorption model and the semi-mechanistic/semi-empirical two-site kinetic surface complexation model, regardless of their complexity, can match our experimental data fairly well under given test conditions. A numerical experiment to further compare the two models shows that they behave differently when the pore velocity is not of the same order of magnitude as our test velocities. This result indicates that further investigations to help determine a better model are needed. We suggest that a multistage column experiment, which tests over the whole range of practical flow velocities, should be conducted to help alleviate inadequate hypothesized models.     

Virus retention and transport through Al-oxide coated sand columns: effects of ionic strength and composition

Knowledge of the factors that influence the fate and transport of viruses in porous media is very important for accurately determining groundwater vulnerability and for developing protective regulations. In this study, six saturated sand column experiments were performed to examine the effects of a positively charged Al-oxide, which was coated on sand particles, on the retention and transport of viruses (phiX174 and MS-2) in background solutions of different ionic strength and composition. We found that the Al-oxide coating on sand significantly removed viruses during their transport in a phosphate buffered saline (PBS) solution. Mass balance calculations showed that 34% of the input MS-2 was inactivated/irreversibly sorbed on the surface of Al-oxide coated sand whereas 100% of phiX174 was recovered. Results from this study also indicated that higher ionic strength facilitated the transport of both phiX174 and MS-2 through the Al-oxide coated sand. This was attributed to the effect of ion shielding, which at higher ionic strength decreased the electrostatic attraction between the viral particles and the sand surface and consequently decreased virus sorption. Strong effect of the ionic strength indicates that an outer-sphere complexation mechanism was responsible for the virus sorption on the Al-oxide coated sand. Ion composition of the background solutions was also found to be a significant factor in influencing virus retention and transport. Virus transport was enhanced in the presence of phosphate (HPO(4)(2-)) as compared to bicarbonate (HCO(3)(-)), and the effect of HPO(4)(2-) was more significant on MS-2 than on phiX174. The presence of bivalent cations (Ca(2+) and Mg(2+)) increased virus transport because the cations partially screened the negative charges on the viruses therefore decreased the electrostatic attraction between the positively charged sand surface and the negatively charged viruses. Mass recovery data indicated that bivalent cations gave rise to a certain degree of inactivation/irreversibly sorption of phiX174 on the surface of Al-oxide coated sand. On the contrary, the bivalent cations appeared to have protected MS-2 from inactivation/irreversibly sorption. This study provides some insights into the mechanisms responsible for virus retention and transport in porous media.     

Effect of goethite coating and humic acid on the transport of bacteriophage PRD1 in columns of saturated sand

The transport of bacteriophage PRD1, a model virus, was studied in columns containing sediment mixtures of quartz sand with goethite-coated sand and using various solutions consisting of monovalent and divalent salts and humic acid (HA). Without HA and in the absence of sand, the inactivation rate of PRD1 was found to be as low as 0.014 day(-1) (at 5+/-3 degrees C), but in the presence of HA it was much lower (0.0009 day(-1)), indicating that HA helps PRD1 to survive. When the fraction of goethite in the sediment was increased, the removal of PRD1 also increased. However, in the presence of HA, C/C0 values of PRD1 increased by as much as 5 log units, thereby almost completely eliminating the effect of addition of goethite. The sticking efficiency was not linearly dependent on the amount of goethite added to the quartz sand; this is apparently due to surface charge heterogeneity of PRD1. Our results imply that, in the presence of dissolved organic matter (DOM), viruses can be transported for long distances thanks to two effects: attachment is poor because DOM has occupied favourable sites for attachment and inactivation of virus may have decreased. This conclusion justifies making conservative assumptions about the attachment of viruses when calculating protection zones for groundwater wells.     

Measurement and analysis of non-Fickian dispersion in heterogeneous porous media

Contaminant breakthrough behavior in a variety of heterogeneous porous media was measured in laboratory experiments, and evaluated in terms of both the classical advection-dispersion equation (ADE) and the continuous time random walk (CTRW) framework. Heterogeneity can give rise to non-Fickian transport patterns, which are distinguished by "anomalous" early arrival and late time tails in breakthrough curves. Experiments were conducted in two mid-scale laboratory flow cells packed with clean, sieved sand of specified grain sizes. Three sets of experiments were performed, using a "homogeneous" packing, a randomly heterogeneous packing using sand of two grain sizes, and an exponentially correlated structure using sand of three grain sizes. Concentrations of sodium chloride tracer were monitored at the inflow reservoir and measured at the outflow reservoir. Breakthrough curves were then analyzed by comparison to fitted solutions from the ADE and CTRW formulations. In all three systems, including the "homogeneous" one, subtle yet measurable differences between Fickian and non-Fickian transport were observed. Quantitative analysis demonstrated that the CTRW theory characterized the full shape of the breakthrough curves far more effectively than the ADE.     

Migration of Cs through quartz sand: experimental results and modeling approaches

We present results from experiments on the migration of ¹³⁷Cs through columns containing quartz sand. Times for ¹³⁷Cs movement through these columns and the quantity of ¹³⁷Cs adsorbed by the sand decreased as the ionic strength of the pore water increased from 0.002 to 0.1 m. The breakthrough curves were characterized by a slow approach towards steady-state concentrations as well as by long tails, indicating that ¹³⁷Cs adsorption to the sand grains was, at least in part, controlled by rate-limited reactions. Various formulations for solute mass transfer were tested for their ability to fit the experimental breakthrough curves. Based on a statistical analysis, a nonlinear, two-site model was identified as the most appropriate for describing the suite of experimental data. Variation in the model parameter that describes the rate of ¹³⁷Cs adsorption to the sand showed no consistent pattern with changes in ionic strength. In contrast, model parameters describing the sorption capacity of the sand grains and the fraction of kinetic sorption sites on the sand decreased with increasing ionic strength. The parameter describing the rate of ¹³⁷Cs desorption varied directly with changes in ionic strength.     

Anomalous transport of colloids and solutes in a shear zone

Transport experiments with colloids and radionuclides in a shear zone were conducted during the Colloid and Radionuclide Retardation experiment (CRR) at Nagra's Grimsel Test Site. Breakthrough curves of bentonite colloids and uranine, a non-sorbing solute, were measured in an asymmetric dipole flow field. The colloid breakthrough is earlier than that of uranine. Both breakthrough curves show anomalously long late time tails and the slope of the late time tails for the colloids is slightly higher. Anomalous late time tails are commonly associated with matrix diffusion processes; the diffusive interaction of solutes transported in open channels with the adjacent porous rock matrix or zones of stagnant water. The breakthrough curves for different colloid size classes are very similar and show no signs of fractionation due to their (size-dependent) diffusivity. It is proposed that tailing of the colloids is mainly caused by the structure of the flow field and that for the colloid transport, matrix diffusion is of minor importance. This has consequences for the interpretation of the uranine breakthrough. Comparisons of experimental results with numerical studies and with the evaluation of the colloid breakthrough with continuous time random theory imply that the tailing in the conservative solute breakthrough in this shear zone is not only caused by matrix diffusion. Part of the tailing can be attributed to advective transport in fracture networks and advection in low velocity regions. Models based on the advection-dispersion equation and matrix diffusion do not properly describe the temporal and spatial evolution of colloid and solute transport in such systems with a consistent set of parameters.     

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.     

A multi-directional tracer test in the fractured Chalk aquifer of E. Yorkshire, UK

A multi-borehole radial tracer test has been conducted in the confined Chalk aquifer of E. Yorkshire, UK. Three different tracer dyes were injected into three injection boreholes and a central borehole, 25 m from the injection boreholes, was pumped at 330 m(3)/d for 8 days. The breakthrough curves show that initial breakthrough and peak times were fairly similar for all dyes but that recoveries varied markedly from 9 to 57%. The breakthrough curves show a steep rise to a peak and long tail, typical of dual porosity aquifers. The breakthrough curves were simulated using a 1D dual porosity model. Model input parameters were constrained to acceptable ranges determined from estimations of matrix porosity and diffusion coefficient, fracture spacing, initial breakthrough times and bulk transmissivity of the aquifer. The model gave equivalent hydraulic apertures for fractures in the range 363-384 microm, dispersivities of 1 to 5 m and matrix block sizes of 6 to 9 cm. Modelling suggests that matrix block size is the primary controlling parameter for solute transport in the aquifer, particularly for recovery. The observed breakthrough curves suggest results from single injection-borehole tracer tests in the Chalk may give initial breakthrough and peak times reasonably representative of the aquifer but that recovery is highly variable and sensitive to injection and abstraction borehole location. Consideration of aquifer heterogeneity suggests that high recoveries may be indicative of a high flow pathway adjacent, but not necessarily connected, to the injection and abstraction boreholes whereas low recoveries may indicate more distributed flow through many fractures of similar aperture.     

Nonequilibrium sorption of phenols onto geosorbents: the impact of pH on intraparticle mass transfer

While the pH effect on sorption equilibrium of weak acids on natural sorbents was investigated in a number of studies, less is known about the pH dependence of sorption kinetics. This paper investigates the impact of pH on sorption kinetics during the transport of some selected phenols through a sandy aquifer material. Breakthrough curves measured in column experiments were analyzed using a mass transfer based nonequilibrium model designated as dispersed flow, film and particle diffusion model (DF-FPDM). In this model, the rate limiting intraparticle diffusion is characterized by the mass transfer coefficient, kSaV, which can be determined from breakthrough curves by curve fitting. The experimental results indicate that the kSaV is pH-dependent and inversely correlated with the pH-dependent distribution coefficient, K(d,app). Regression equations are presented that may be used to estimate approximate values of intraparticle mass transfer coefficients on the basis of experimentally determined or LFER predicted distribution coefficients.     

Filtration and transport of Bacillus subtilis spores and the F-RNA phage MS2 in a coarse alluvial gravel aquifer: implications in the estimation of setback distances

Filtration of Bacillus subtilis spores and the F-RNA phage MS2 (MS2) on a field scale in a coarse alluvial gravel aquifer was evaluated from the authors' previously published data. An advection-dispersion model that is coupled with first-order attachment kinetics was used in this study to interpret microbial concentration vs. time breakthrough curves (BTC) at sampling wells. Based on attachment rates (katt) that were determined by applying the model to the breakthrough data, filter factors (f) were calculated and compared with f values estimated from the slopes of log (cmax/co) vs. distance plots. These two independent approaches resulted in nearly identical filter factors, suggesting that both approaches are useful in determining reductions in microbial concentrations over transport distance. Applying the graphic approach to analyse spatial data, we have also estimated the f values for different aquifers using information provided by some other published field studies. The results show that values of f, in units of log (cmax/co) m(-1), are consistently in the order of 10(-2) for clean coarse gravel aquifers, 10(-3) for contaminated coarse gravel aquifers, and generally 10(-1) for sandy fine gravel aquifers and river and coastal sand aquifers. For each aquifer category, the f values for bacteriophages and bacteria are in the same order-of-magnitude. The f values estimated in this study indicate that for every one-log reduction in microbial concentration in groundwater, it requires a few tens of meters of travel in clean coarse gravel aquifers, but a few hundreds of meters in contaminated coarse gravel aquifers. In contrast, a one-log reduction generally only requires a few meters of travel in sandy fine gravel aquifers and sand aquifers. Considering the highest concentration in human effluent is in the order of 10(4) pfu/l for enteroviruses and 10(6) cfu/100 ml for faecal coliform bacteria, a 7-log reduction in microbial concentration would comply with the drinking water standards for the downgradient wells under natural gradient conditions. Based on the results of this study, a 7-log reduction would require 125-280 m travel in clean coarse gravel aquifers, 1.7-3.9 km travel in contaminated coarse gravel aquifers, 33-61 m travel in clean sandy fine gravel aquifers, 33-129 m travel in contaminated sandy fine gravel aquifers, and 37-44 m travel in contaminated river and coastal sand aquifers. These recommended setback distances are for a worst-case scenario, assuming direct discharge of raw effluent into the saturated zone of an aquifer. Filtration theory was applied to calculate collision efficiency (alpha) from model-derived attachment rates (katt), and the results are compared with those reported in the literature. The calculated alpha values vary by two orders-of-magnitude, depending on whether collision efficiency is estimated from the effective particle size (d10) or the mean particle size (d50). Collision efficiency values for MS-2 are similar to those previously reported in the literature (e.g. ) [DeBorde, D.C., Woessner, W.W., Kiley, QT., Ball, P., 1999. Rapid transport of viruses in a floodplain aquifer. Water Res. 33 (10), 2229-2238]. However, the collision efficiency values calculated for Bacillus subtilis spores were unrealistic, suggesting that filtration theory is not appropriate for theoretically estimating filtration capacity for poorly sorted coarse gravel aquifer media. This is not surprising, as filtration theory was developed for uniform sand filters and does not consider particle size distribution. Thus, we do not recommend the use of filtration theory to estimate the filter factor or setback distances. Either of the methods applied in this work (BTC or concentration vs. distance analyses), which takes into account aquifer heterogeneities and site-specific conditions, appear to be most useful in determining filter factors and setback distances.     

Analysis of the long-term behavior of solute transport with nonlinear equilibrium sorption using breakthrough curves and temporal moments

We analyzed the long-term behavior of breakthrough curves (BTCs) and temporal moments of a solute subjected to Freundlich equilibrium sorption (s = kc(n)). For one-dimensional transport in a homogeneous porous medium, we derived a power-law relation between travel time, tau, and solute displacement, chi, with the exponent being equal to the Freundlich n exponent. The mean solute velocity, derived from the first time moment, was found to change as tau(n-1). For n values larger than 0.66, the second time moment could be related to c chi(2/n), where c is a constant. An approach based on the use of a critical concentration was developed to estimate the presence of the asymptotic regime in the tail of the BTC. This approach was tested successfully using numerical case studies. One-dimensional numerical simulations with varying values of k, n and initial mass were run to verify the closed form analytical expressions for the large time behavior of temporal moments and the tailing part of breakthrough curves. Good agreement between the slope of the tailing part of log-log transformed BTCs and the predicted slope using asymptotic theory was found. Asymptotic theory in general underestimated the magnitude of the concentration in the tail. The quality of the estimated concentrations in the tail improved for small values of the dispersivity. Experimental BTCs of uranin and benazolin were analyzed in combination with sorption/desorption batch experiments using asymptotic theory. A good agreement between the value of n parameter derived from desorption experiment with benazolin and the value of the n parameter derived from the tail of the BTC was found.     

Transport of Escherichia coli through variably saturated sand columns and modeling approaches

A sand column leaching system with well-controlled suction and flow rate was built to investigate the effects on bacterial transport of air-water interface effects (AWI) correlated to water content, particle size, and column length. Adsorption of Escherichia coli strain D to silica sands was measured in batch tests. The average % adsorption for coarse and fine sands was 45.9+/-7.8% and 96.9+/-3.2%, respectively. However, results from static batch adsorption experiments have limited applicability to dynamic bacterial transport in columns. The early breakthrough of E. coli relative to bromide was clear for all columns, namely c. 0.15 to 0.3 pore volume earlier. Column length had no significant effects on the E. coli peak concentration or on total recovery in leachate, indicating retention in the top layer of sands. Tailing of breakthrough curves was more prominent for all fine sand columns than their coarse sand counterparts. Bacterial recovery in leachate from coarse and saturated sand columns was significantly higher than from fine and unsaturated columns. Observed data were fitted by the convection-dispersion model, amended for one-site and two-site adsorption to particles, and for air-water interface (AWI) adsorption. Among all models, the two-site+AWI model achieved consistently high model efficiency for all experiments. Thus it is evident from experimental and modeling results that AWI adsorption plays an important role in E. coli transport in sand columns.     

Breakthrough adsorption study of migratory nickel in fine-grained soil.

The present study was conducted to evaluate the breakthrough curve for nickel adsorption in fine-grained soil from a nearby ash pond site of a thermal power plant. Nickel was found to be the major polluting solute in the ash sluicing wastewater. The adsorption of nickel by vertical soil column batch test and horizontal migration test was carried out in the laboratory. Field investigation was conducted also, by installing test wells around the ash pond site. Experimental results showed a good adsorptive capacity of soil for nickel ions. The breakthrough curves showed a reasonable fitting with a one-dimensional mathematical model. The breakthrough curves yielded from field test results showed good agreement with a two-dimensional mathematical model.     

Cesium migration in saturated silica sand and Hanford sediments as impacted by ionic strength

Large amounts of 137Cs have been accidentally released to the subsurface from the Hanford nuclear site in the state of Washington, USA. The cesium-containing liquids varied in ionic strengths, and often had high electrolyte contents, mainly in the form of NaNO3 and NaOH, reaching concentrations up to several moles per liter. In this study, we investigated the effect of ionic strengths on Cs migration through two types of porous media: silica sand and Hanford sediments. Cesium sorption and transport was studied in 1, 10, 100, and 1000 mM NaCl electrolyte solutions at pH 10. Sorption isotherms were constructed from batch equilibrium experiments and the batch-derived sorption parameters were compared with column breakthrough curves. Column transport experiments were analyzed with a two-site equilibrium-nonequilibrium model. Cesium sorption to the silica sand in batch experiments showed a linear sorption isotherm for all ionic strengths, which matched well with the results from the column experiments at 100 and 1000 mM ionic strength; however, the column experiments at 1 and 10 mM ionic strength indicated a nonlinear sorption behavior of Cs to the silica sand. Transport through silica sand occurred under one-site sorption and equilibrium conditions. Cesium sorption to Hanford sediments in both batch and column experiments was best described with a nonlinear Freundlich isotherm. The column experiments indicated that Cs transport in Hanford sediments occurred under two-site equilibrium and nonequilibrium sorption. The effect of ionic strength on Cs transport was much more pronounced in Hanford sediments than in silica sands. Effective retardation factors of Cs during transport through Hanford sediments were reduced by a factor of 10 when the ionic strength increased from 100 to 1000 mM; for silica sand, the effective retardation was reduced by a factor of 10 when ionic strength increased from 1 to 1000 mM. A two order of magnitude change in ionic strength was needed in the silica sand to observe the same change in Cs retardation as in Hanford sediments.     

Transport of biocolloids in water saturated columns packed with sand: Effect of grain size and pore water velocity

The main objective of this study was to evaluate the combined effects of grain size and pore water velocity on the transport in water saturated porous media of three waterborne fecal indicator organisms (Escherichia coli, MS2, and ΦX174) in laboratory-scale columns packed with clean quartz sand. Three different grain sizes and three pore water velocities were examined and the attachment behavior of Escherichia coli, MS2, and ΦX174 onto quartz sand was evaluated. The mass recoveries of the biocolloids examined were shown to be highest for Escherichia coli and lowest for MS2. However, no obvious relationships between mass recoveries and water velocity or grain size could be established from the experimental results. The observed mean dispersivity values for each sand grain size were smaller for bacteria than coliphages, but higher for MS2 than ΦX174. The single collector removal and collision efficiencies were quantified using the classical colloid filtration theory. Furthermore, theoretical collision efficiencies were estimated only for E. coli by the Interaction-Force-Boundary-Layer, and Maxwell approximations. Better agreement between the experimental and Maxwell theoretical collision efficiencies were observed.     

Continuous time random walk model better describes the tailing of atrazine transport in soil

Contaminant transport in soils is complicated and involves some physical and chemical nonequilibrium processes. In this research, the soil column displacement experiments of Cl⁻ and atrazine under different flow velocities were carried out. The data sets of Cl⁻ transport in sandy loam fitted to the convection dispersion equation (CDE) and the two-region model (TRM) indicated that the effects of physical nonequilibrium process produced by immobile water on the breakthrough curves (BTCs) of Cl⁻ and atrazine transport through the repacking soil columns were negligible. The two-site model (TSM) and the continuous time random walk (CTRW) were also used to fit atrazine transport behavior at the flow rate of 19.86 cm h⁻¹. The CTRW convincingly captured the full evolution of atrazine BTC in the soil column, especially for the part of long tailing. However, the TSM failed to characterize the tailing of atrazine BTC in the soil column. The calculated fraction of equilibrium sorption sites, F, ranging from 0.78 to 0.80 for all flow rates suggested the contribution of nonequilibrium sorption sites to the asymmetry of atrazine BTCs. Furthermore, the data sets for the flow rates of 6.68 cm h⁻¹ and 32.81 cm h⁻¹ were predicted by the TSM and the CTRW. As to the flow rate of 6.68 cm h⁻¹, the CTRW predicted the entire BTC of atrazine transport better than the TSM did. For the flow rate of 32.81 cm h⁻¹, the CTRW characterized the late part of the tail better, while the TSM failed to predict the tailings of atrazine BTC.