An integrated sorption-diffusion model for the calculation of consistent distribution and diffusion coefficients in compacted bentonite

A thermodynamic sorption model and a diffusion model based on electric double layer (EDL) theory are integrated to yield a surface chemical model that treats porewater chemistry, surface reactions, and the influence of charged pore walls on diffusing ions in a consistent fashion. The relative contribution of Stern and diffuse layer to the compensation of the permanent surface charge represents a key parameter; it is optimized for the diffusion of Cs in Kunipia-F bentonite, at a dry density of 400 kg/m3. The model is then directly used to predict apparent diffusivities (Da) of Cs, Sr, Cl-, I- and TcO4- and corresponding distribution coefficients (Kd) of Cs and Sr in different bentonites as a function of dry density, without any further adjustment of surface chemical and EDL parameters. Effective diffusivities (De) for Cs, HTO, and TcO4- are also calculated. All calculated values (Da, De, Kd) are fully consistent with each other. A comparison with published, measured data shows that the present model allows a good prediction and consistent explanation of (i) apparent and effective diffusivities for cations, anions, and neutral species in compacted bentonite, and of (ii) Kd values in batch and compacted systems.     

Diffusivity measurement of heavy ions in Wyoming montmorillonite gels by X-ray computed tomography

Medical X-ray computed tomography (CT) was applied to the measurement of the diffusion coefficients of heavy ions in an artificial barrier material for the disposal of nuclear wastes. Cs(+), Sr(2+), I(-), and Br(-) are the heavy ions measured and the barrier used is the water-rich gel of Wyoming montmorillonite (86.5-100 wt.% H(2)O). X-ray CT yields an inevitable artifact (beam-hardening) in the obtained images. Before the diffusion experiments, the polychromatic primary X-ray spectrum of the CT scanner was measured by a CdZnTe detector, and the effects of the artifact were examined for an aqueous CsCl solution sample. The results show that the beam-hardening artifact derived from the polychromatic photon energy distribution can be suppressed by applying a special image reconstruction method assuming the chemical composition of samples. The transient one-dimensional diffusion of heavy ions in a plastic container filled with the gel was imaged nondestructively by the X-ray CT scanner with an in-plane resolution of 0.31 mm and slice thickness of 2 mm. The results show that diffusivities decrease with increasing clay weight fraction. The degree of the diffusivity decrease was high for cations (Cs(+) and Sr(2+)) and low for anions (I(-) and Br(-)). The quantitative decomposition of the contribution of the geometrical tortuosity and of the sorption to the diffusivity was performed by subtracting the diffusivity of nonsorbing I(-) from the measured diffusivities. The results show that the contribution of the sorption is large for Cs(+), Sr(2+) and small for Br(-). Because X-ray CT allows nondestructive and quick measurements of diffusivities, the technique would be useful particularly for measuring the diffusive migration of harmful radioactive elements.     

Role of magnetite and humic acid in radionuclide migration in the environment

Sorption of ¹³⁷Cs, ⁹⁰Sr, ¹⁵⁴Eu and ¹⁴¹Ce by magnetite has been studied at varying pH (4 to 11) in the presence and absence of humic acid. The sorption studies have also been carried out at varying ionic strength (0.01 to 0.2 M NaClO₄) and humic acid concentration (2 to 20 mg/L). Percentage sorption of ¹³⁷Cs and ⁹⁰Sr was found to be pH dependent, with the sorption increasing with increasing pH of the suspension. At any pH, the percentage sorption of ⁹⁰Sr was higher than that of ¹³⁷Cs. The results have been explained in terms of the electrostatic interaction between the positively charged metal ions and the surface charge of the magnetite which becomes increasingly negative with increasing pH. On the other hand, ¹⁵⁴Eu and ¹⁴¹Ce were found to be strongly sorbed by the magnetite at all pH values, with the sorption being independent of pH. The strong sorption of trivalent and tetravalent metal ions suggests the role of complexation reactions during sorption, apart from the electrostatic interactions. However, in the case of ¹⁴¹Ce surface precipitation of Ce(III) formed by reduction of Ce(IV) in the presence of magnetite cannot be ruled out. Presence of humic acid (2 mg/L) was found to have negligible effect on sorption of all metal ions.     

In situ anion diffusion experiments using radiotracers

Diffusion experiments in compacted bentonite have been carried out in situ using the borehole laboratory CHEMLAB. The "ordinary" anion iodide and the redox-sensitive pertechnetate ion have been investigated. In spite of strongly reducing groundwater conditions, technetium was found to diffuse mostly unreduced as TcO4-, although in some spots in the compacted clay, the activity was significantly higher, which may be explained by reduction of some TcO4- by iron-containing minerals in the bentonite. The measured concentration profiles in the clay cannot be accommodated by assuming one single diffusion process. The experimental data are modeled assuming two diffusion paths, intralamellar diffusion and diffusion in external water. The apparent diffusivity for the intralamellar diffusion was found to be 8.6 x 10(-11) m2 s(-1) for iodide with a capacity factor of 0.1, while the apparent diffusivity for the diffusion in external water was found to be 5 x 10(-14) m2 s(-1) with alpha=2.26. The corresponding values for Tc were found to be Da= 6 x 10(-11) m2 s(-1), alpha=0.1 and Da= 1 x 10(-13) m2 s(-1), alpha=0.46, respectively. The diffusion constants and capacity factors obtained in this study are in accordance with data from laboratory experiments.     

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.     

Ground discarded tires remove naphthalene, toluene, and mercury from water

Ground discarded tires adsorb naphthalene, toluene, and mercury ions (Hg2+) from aqueous solutions. Their sorption properties and kinetics were determined by batch equilibration techniques at 20 degrees C. The isotherms were linear for naphthalene and toluene and their sorption coefficients were about 1340 and 255 (ml/g), respectively. Sorption of the organic compounds by the ground rubber particles was relatively fast (within 30 min). However, the mercury isotherms were non-linear, and its sorption was slow as compared to the sorption of the organics. The rubber particles had a strong affinity for Hg2+. These results show that ground discarded tires are effective in removing organic compounds and Hg2+ from wastewater and other contaminated environments. In addition it would be a useful, environmentally friendly use of discarded tires (one tire per year per capita is discarded in the United States).     

Investigation of solution chemistry effects on sorption behavior of Cu(II) on ZSM-5 zeolite

The sorption of Cu(II) from an aqueous solution using ZSM-5 zeolite was investigated by batch technique under ambient conditions. Sorption was investigated as a function of pH, ionic strength, foreign ions, humic substances, and temperature. The results indicate that the sorption of Cu(II) on ZSM-5 zeolite is strongly dependent on pH. Sorption is dependent on ionic strength at low pH, but independent of ionic strength at high pH values. The presence of humic/fluvic acid (HA/FA) enhances the sorption of Cu(II) on ZSM-5 zeolite at low pH values, and reduces Cu(II) sorption at high pH values. Sorption isotherms were well simulated by the Langmuir model. Thermodynamic parameters (i.e., deltaH0, deltaS0 and deltaG0) for the sorption of Cu(II) were determined from temperature-dependent sorption isotherms at 293.15, 313.15, and 333.15 K, respectively. Results indicate that the sorption process of Cu(II) on ZSM-5 zeolite is spontaneous and endothermic.     

Study of sorption kinetics of some ionic liquids on different soil types

In the present contribution sorption kinetics experiments under static conditions were utilized in three selected ionic liquids cations (1-ethyl-3-methylimidazolium, 1-butyl-3-methylimidazolium, 1-hexyl-3-methylimidazolium chlorides) study with five type of soil, differing in total organic carbon (TOC) content. The experimental results indicate the sorption capacity growth with increase in TOC content and hydrophobicity of ionic liquid cation. The obtained kinetic sorption parameters as well as distribution coefficients (K_d) were used to estimate the sorption properties of the soil types towards the ionic liquids in question. The Gibbs free energy values indicate that ionic liquid cations sorption on soils could be generally considered as a physical adsorption with exothermic effect. But the values of −dG for studied cations sorption on soil with very high of TOC content in soil (45%) may testify to nature of chemical adsorption. Sorption of the analyzed compounds occurs probably by means of hydrogen bonds, electrostatic and π  π interaction with the organic matter and the clay minerals of the soils.     

Utilization of Modified Attapulgite for the Removal of Sr(II), Co(II), and Ni(II) Ions from Multicomponent System,Part I: Kinetic Studies

Radionuclide sorption by natural and modified clays is extensively accepted to be an important process from the radioactive waste point of view. This work focused on modification of natural attapulgite with a layered double hydroxide to produce a novel chemisorbent for Sr2+, Ni2+, and Co2+ removal from multicomponent solution. The structural and surface characteristics of both attapulgite (ATP) and modified attapulgite (LDH-ATP) were investigated using XRD, FTIR, SEM, and thermal analysis. Comparison of sorption features of Sr2+, Ni2+, and Co2+ onto ATP and LDH-ATP was achieved; the results indicated that LDH-ATP was the most efficient sorbent for Sr2+, Ni2+, and Co2+. Kinetic studies established that the sorption is fast and reaching >90% within 30 min. The sorption of Sr2+, Ni2+, and Co2+ are well defined by non-linear pseudo-second-order model and controlled by an intra-particle diffusion mechanism. The diffusivity was determined using homogeneous surface diffusion (HSDM) model and found in the order 10−13 m2/min; this confirmed that the sorption of the three ions is chemisorption process. LDH-ATP can be employed as a candidate chemisorbent for the removal of some metal ions from waste solution.     

Sensitivity analysis of radionuclide migration in compacted bentonite: a mechanistic model approach

Mechanistic model calculations for the migration of Cs, Ra, Am and Pb in compacted bentonite have been carried out to evaluate sensitivities with respect to different parameter variations. A surface chemical speciation/electric double layer model is used to calculate: (i) porewater composition and radionuclide speciation in solution and at the bentonite surface, yielding the distribution of mobile and sorbed species and (ii) interaction of diffusing species with negatively charged pore walls to obtain diffusion parameters. The basic scenario considers the interaction of compacted bentonite with a fresh-type groundwater; variations include the presence of bentonite impurities and saline groundwater. It is shown that these scenarios result in significant variations of porewater composition that affect migration via three mechanisms that can partly compensate each other: (1) effects on sorption through radionuclide complexation in solution, and competition of major cations for surface sites; (2) changes in radionuclide solution speciation leading to different diffusing species under different conditions; (3) effects on diffusion through changes in the electric double layer properties of the clay pores as a function of ionic strength.     

Solute transport properties of compacted Ca-bentonite used in FEBEX project

The present Spanish concept of a deep geological high level waste repository includes an engineered clay barrier around the canister. The clay presents a very high sorption capability for radionuclides and a very small hydraulic conductivity, so that the migration process of solutes is limited by sorption and diffusion processes. Therefore, diffusion and distribution coefficients in compacted bentonite (i.e. in "realistic" liquid to solid ratio conditions) are the main parameters that have to be obtained in order to characterise solute transport that could be produced after the canister breakdown. Through-Diffusion (TD) and In-Diffusion (ID) experiments with HTO, Sr, Cs and Se were carried out using compacted FEBEX bentonite, which is the reference material for the Spanish concept of radioactive waste disposal. Experiments were interpreted by means of available analytical solutions that allow the estimation of diffusion coefficients and, in some cases, distribution coefficients. Analytical solutions are simple to use, but rely on hypotheses that do not hold in all the experiments. These experiments were interpreted also using an automatic parameter estimation code that overcomes the limitations of analytical solutions. Numerical interpretation allows the simultaneous estimation of porosity, diffusion and distribution coefficients, accounts for the role of porous sinters and time-varying boundary concentrations, and can use different types of raw concentration data.     

Complexation of the antibiotic tetracycline with humic acid

The effect of solution chemistry and sorbate-to-sorbent ratio on the interaction of the antibiotic tetracycline with Elliott soil humic acid (ESHA) was investigated using equilibrium dialysis and FITEQL modeling. Tetracycline speciation strongly influenced its sorption to ESHA. Sorption was strongly pH-dependent with a maximum around pH 4.3, and competition with H+ and electrolyte cation (Na+) was evident. The pH-dependent trend was consistent with complexation between the cationic/zwitterionic tetracycline species and deprotonated sites in ESHA (mainly carboxylic functional groups). Modification of ESHA by Ca2+ addition increased tetracycline sorption suggesting that ternary complex formation (ESHA-metal-tetracycline) may be important at higher concentrations of multivalent metal cations. The macroscopic data (pH-envelope and sorption isotherms) were successfully modeled using a discrete logK function with the FITEQL 4.0 chemical equilibrium program indicating that ESHA-tetracycline interaction could be reasonably represented as complex formation of a monoacid with discrete sites in humic acid. Sorption-desorption hysteresis was observed; both sorption and desorption isotherms were well described by the Freundlich equation.     

Experimental and modeling studies on sorption and diffusion of radium in bentonite

The sorption and desorption behavior of radium on bentonite and purified smectite was investigated as a function of pH, ionic strength and liquid to solid ratio by batch experiments. The distribution coefficients (Kd) were in the range of 10(2) to > 10(4) ml g-1 and depended on ionic strength and pH. Most of sorbed Ra was desorbed by 1 M KCl. The results for purified smectite indicated that Ra sorption is dominated by ion exchange at layer sites of smectite, and surface complexation at edge sites may increase Ra sorption at higher pH region. Reaction parameters between Ra and smectite were determined based on an interaction model between smectite and groundwater. The reaction parameters were then used to explain the results of bentonite by considering dissolution and precipitation of minerals and soluble impurities. The dependencies of experimental Kd values on pH, ionic strength and liquid to solid ratio were qualitatively explained by the model. The modeling result for bentonite indicated that sorption of Ra on bentonite is dominated by ion exchange with smectite. The observed pH dependency was caused by changes of Ca concentration arising from dissolution and precipitation of calcite. Diffusion behavior of Ra in bentonite was also investigated as a function of dry density and ionic strength. The apparent diffusion coefficients (Da) obtained in compacted bentonite were in the range of 1.1 x 10(-11) to 2.2 x 10(-12) m2 s-1 and decreased with increasing in dry density and ionic strength. The Kd values obtained by measured effective diffusion coefficient (De) and modeled De were consistent with those by the sorption model in a deviation within one order of magnitude.     

Comparative kinetic desorption of 60Co, 85Sr and 134Cs from a contaminated natural silica sand column: influence of varying physicochemical conditions and dissolved organic matter

In order to determine the mechanisms of the retention of 60Co, 85Sr and 134Cs in natural silica sand columns, desorption experiments were performed by changes of pH and ionic strength and by injection of natural organic matter (NOM). Injection of KCl (0.1 M) resulted in a high release of 60Co (60-100%) and 85Sr (72-100%) but a smaller release of 134Cs (31-66%). Only limited release of 60Co (66%) and 85Sr (71%) and no release of 134Cs were observed by injection of NOM. The different percentages of desorption were related to the chemical characteristics of the organic colloids previously retained in columns before the desorption step. The results evidenced different sorption processes on energetically heterogeneous surface sites. According to the initial conditions, the binding of the radionuclides to the solid phase resulted from weak and easily reversible sorption processes to strong association probably by inner sphere complexes. The rather weak release of 134Cs by KCl was attributed to the strong retention of 134Cs by clay coatings on the natural silica sand surfaces.     

A cation exchange model to describe Cs+ sorption at high ionic strength in subsurface sediments at Hanford site, USA

A theoretical and experimental study of cation exchange in high ionic strength electrolytes was performed using pristine subsurface sediments from the U.S. Department of Energy Hanford site. These sediments are representative of the site contaminated sediments impacted by release of high level waste (HLW) solutions containing 137Cs+ in NaNO3 brine. The binary exchange behavior of Cs+-Na+, Cs+-K+, and Na+-K+ was measured over a range in electrolyte concentration. Vanselow selectivity coefficients (Kv) that were calculated from the experimental data using Pitzer model ion activity corrections for aqueous species showed monotonic increases with increasing electrolyte concentrations. The influence of electrolyte concentration was greater on the exchange of Na+-Cs+ than K+-Cs+, an observation consistent with the differences in ion hydration energy of the exchanging cations. A previously developed two-site ion exchange model [Geochimica et Cosmochimica Acta 66 (2002) 193] was modified to include solvent (water) activity changes in the exchanger phase through application of the Gibbs-Duhem equation. This water activity-corrected model well described the ionic strength effect on binary Cs+ exchange, and was extended to the ternary exchange system of Cs+-Na+-K+ on the pristine sediment. The model was also used to predict 137Cs+ distribution between sediment and aqueous phase (Kd) beneath a leaked HLW tank in Hanfordd's S-SX tank using the analytical aqueous data from the field and the binary ion exchange coefficients for the pristine sediment. The Kd predictions closely followed the trend in the field data and were improved by consideration of water activity effects that were considerable in certain regions of the vadose zone plume.     

Simultaneous sorption of lead and chlorobenzene by organobentonite

Clays or organoclays have been used as a barrier to prevent the transport of hazardous contaminants in landfills. However, clays are known to effectively sorb mostly inorganic contaminants, while organoclays are mainly used for organic contaminants. Since the organoclays are basically clay particles modified with cationic surfactants, there might exist an optimal coverage of cationic surfactant on the clay particles to sorb both inorganic and organic contaminants. In order to determine the optimal mass of cationic surfactants on the bentonites, sodium bentonites were treated with various ratios of hexadecyltrimethylammonium (HDTMA) to bentonites. Chlorobenzene and lead were selected as representative contaminants. When either chlorobenzene or lead exists as a single contaminant, chlorobenzene sorption increased with increasing HDTMA to bentonite ratios, and lead sorption decreased with increasing HDTMA to bentonite ratios. Sorption of chlorobenzene was a function of HDTMA coverage on the bentonites, while lead sorption was much more influenced by the initial lead concentration rather than the mass of HDTMA added to the bentonites.     

Sorption of chlorimuron-ethyl on montmorillonite clays: effects of exchangeable cations,pH, and ionic strength

Sorption interaction of chlorimuron-ethyl with montmorillonite clays was investigated under varied types of exchangeable cation, pH, and ionic strength conditions. Chlorimuron-ethyl sorption on bentonites exhibited pronounced cation dependency, and the sorption ability increased as the sequence Ca²⁺- < Na⁺- < Al³⁺- < Fe³⁺-bentonite, due to different sorption mechanisms, whereas the cation dependency was influenced by the clay type and much weaker for montmorillonites. The decrease of pH at the range of 4.0–6.0 prominently increased sorption of chlorimuron-ethyl on all cation-exchanged montmorillonite clays, and nearly a neglected sorption (about 2 %) can be observed at pH over 7.0. In the presence of CaCl₂, sorption of chlorimuron-ethyl on Fe³⁺-bentonite was promoted because of complexion of Ca²⁺ and the surface of Fe³⁺-bentonite. However, as the concentration of CaCl₂ increased, chlorimuron-ethyl sorption on Ca²⁺- and Fe³⁺-exchanged bentonite decreased, suggesting that Ca bridging was not the prevailing mechanism for sorption of chlorimuron-ethyl on these clays. Furthermore, chlorimuron-ethyl sorption was relatively sensitive to pH, and the change of pH may obscure effect of other factors on the sorption, so it was quite necessary to control pH at a constant value when the effect of other factor was being studied.     

A generalised sorption model for the concentration dependent uptake of caesium by argillaceous rocks

A three-site cation exchange model is proposed to describe the concentration dependent uptake of Cs on natural argillaceous rock systems. Major premises in the model are that the sorption of Cs is dominated by the illite mineral component in the rock and that there is a fixed relationship between the site capacities of the three site types denoted as frayed edge, type II and planar sites. The definition of a “reference illite” with a cation exchange capacity of 0.2 equiv. kg^-1 allows the three site capacities to be fixed in the model calculations over the weight fraction of illite in the argillaceous rocks. Up to Cs equilibrium concentrations of 10^-3 M sorption occurs predominantly on the frayed edge and type II sites (higher affinity sites), with the planar site type playing only a minor role. Competition with Cs for sorption on the former two site types arises predominantly from monovalent cations such as K, Rb and NH₄ which have low hydration energies. H and Na (except at high concentrations) are considerably less competitive and bivalent cations such as Mg, Ca and Sr are effectively non-competitive. A consistent set of selectivity coefficients for Cs with respect to K, Rb, NH₄ and Na was derived from analyses and modelling of a wide range of Cs sorption data available in the open literature on pure illites from many different sources. The model was tested against four Cs sorption isotherm data sets determined on argillaceous rocks: Boom clay, Oxford clay, Palfris marl and Opalinus clay. The water chemistries and illite contents given in these experiments allowed the Cs sorption isotherms to be predicted. It is concluded that the Cs sorption model presented here, in which there are no free parameters, can be used to predict the uptake of Cs at equilibrium concentrations below 10^-3 M to within a factor of 2 to 3 in natural argillaceous rock systems.     

Sorption and desorption of Cu and Cd by macroalgae and microalgae

The sorption and desorption of Cu and Cd by two species of brown macroalgae and five species of microalgae were studied. The two brown macroalgae, Laminaria japonica and Sargassum kjellmanianum, were found to have high capacities at pHs between 4.0 and 5.0 while for microalgae, optimum pH lay at 6.7. The presence of other cations in solution was found to reduce the sorption of the target cation, suggesting a competition for sorption sites on organisms. Sorption isotherms obeyed the Freundlich equation, suggesting involvement of a multiplicity of mechanisms and sorption sites. For the microalgae tested, Spirulina platensis had the highest capacity for Cd, followed by Nannochloropsis oculata, Phaeodactylum tricornutum, Platymonas cordifolia and Chaetoceros minutissimus. The reversibility of metal sorption by macroalgae was examined and the results show that both HCl and EDTA solutions were very effective in desorbing sorbed metal ions from macroalgae, with up to 99.5% of metals being recovered. The regenerated biomass showed undiminished sorption performance for the two metals studied, suggesting the potential of such material for use in water and wastewater treatment.     

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.