Cotransport of clay colloids and viruses in water saturated porous media

Upscaling of nanoparticle transport in porous media under unfavorable conditions: Pore scale to Darcy scale

Transport and deposition of nanoparticles in porous media is a multi-scale problem governed by several pore-scale processes, and hence, it is critical to link the processes at pore scale to the Darcy-scale behavior. In this study, using pore network modeling, we develop correlation equations for deposition rate coefficients for nanoparticle transport under unfavorable conditions at the Darcy scale based on pore-scale mechanisms. The upscaling tool is a multi-directional pore-network model consisting of an interconnected network of pores with variable connectivities. Correlation equations describing the pore-averaged deposition rate coefficients under unfavorable conditions in a cylindrical pore, developed in our earlier studies, are employed for each pore element. Pore-network simulations are performed for a wide range of parameter values to obtain the breakthrough curves of nanoparticle concentration. The latter is fitted with macroscopic 1-D advection-dispersion equation with a two-site linear reversible deposition accounting for both equilibrium and kinetic sorption. This leads to the estimation of three Darcy-scale deposition coefficients: distribution coefficient, kinetic rate constant, and the fraction of equilibrium sites. The correlation equations for the Darcy-scale deposition coefficients, under unfavorable conditions, are provided as a function of measurable Darcy-scale parameters, including: porosity, mean pore throat radius, mean pore water velocity, nanoparticle radius, ionic strength, dielectric constant, viscosity, temperature, and surface potentials of the particle and grain surfaces. The correlation equations are found to be consistent with the available experimental results, and in qualitative agreement with Colloid Filtration Theory for all parameters, except for the mean pore water velocity and nanoparticle radius.

Porous media grain size distribution and hydrodynamic forces effects on transport and deposition of suspended particles

The effects of porous media grain size distribution on the transport and deposition of polydisperse suspended particles under different flow velocities were investigated. Selected Kaolinite particles (2-30μm) and Fluorescein (dissolved tracer) were injected in the porous media by step input injection technique. Three sands filled columns were used: Fine sand, Coarse sand, and a third sand (Mixture) obtained by mixing the two last sands in equal weight proportion. The porous media performance on the particle removal was evaluated by analysing particles breakthrough curves, hydro-dispersive parameters determined using the analytical solution of convection-dispersion equation with a first order deposition kinetics, particles deposition profiles, and particle-size distribution of the recovered and the deposited particles. The deposition kinetics and the longitudinal hydrodynamic dispersion coefficients are controlled by the porous media grain size distribution. Mixture sand is more dispersive than Fine and Coarse sands. More the uniformity coefficient of the porous medium is large, higher is the filtration efficiency. At low velocities, porous media capture all sizes of suspended particles injected with larger ones mainly captured at the entrance. A high flow velocity carries the particles deeper into the porous media, producing more gradual changes in the deposition profile. The median diameter of the deposited particles at different depth increases with flow velocity. The large grain size distribution leads to build narrow pores enhancing the deposition of the particles by straining.

Experimental and Numerical Investigations of Silver Nanoparticle Transport under Variable Flow and Ionic Strength in Soil

Unsaturated column experiments were conducted with an undisturbed loamy sand soil to investigate the influence of flow interruption (FI) and ionic strength (IS) on the transport and retention of surfactant-stabilized silver nanoparticles (AgNP) and the results were compared to those obtained under continuous flow conditions. AgNP concentrations for breakthrough curves (BTCs) and retention profiles (RPs) were analyzed by ICP-MS. Experimental results were simulated by the numerical code HP1 (Hydrus-PhreeqC) with the DLVO theory, extended colloid filtration theory and colloid release model. BTCs of AgNP showed a dramatic drop after FI compared to continuous flow conditions. Evaporation increased due to FI, resulting in increased electrical conductivity of the soil solution, which led to a totally reduced mobility of AgNP. A reduction of IS after FI enhanced AgNP mobility slightly. Here the strongly increased Al and Fe concentration in the effluent suggested that soil colloids facilitated the release of AgNP (cotransport). The numerical model reproduced the measured AgNP BTCs and indicated that attachment to the air-water interface (AWI) occurring during FI was the key process for AgNP retention.

Interactions of infectious F-specific RNA bacteriophages with suspended matter and sediment: Towards an understanding of FRNAPH distribution in a river water system

The association of viruses with settling particles is certainly a major process controlling the spread of viral pollution in surface water and sediment. To better understand the viral distribution in a river system, the behavior of F-specific RNA bacteriophages (FRNAPHs) was investigated in relationship with the suspended solids and sediment. The partitioning of phage particles (free or associated with solids) in surface water and the attachment capabilities of eight distinct strains of phages to sediment were studied in lab experiments. In situ observations were also performed with the genotyping of 166 individual plaques of FRNAPHs isolated from surface water and sediment. The results reported here demonstrate a variation of the status of infectious phages as a function of the hydro-climatological conditions. Phage-solid association seems to mainly occur during the peak of rainfall-runoff events but also to a certain extent during the recession phase compared to low flow conditions. The transfer of phages from the water column to sediment may occur at this time. Furthermore, the ability of FRNAPHs to interact with sediment was established for six strains out of eight, belonging to genogroups II, III and IV. A similar dynamic was observed for strains within a same genogroup despite different intensity of attachment and inactivation rates for strains of genogroups III and IV. The latter results match the in situ observations in the water and sediment compartments of the studied area. Infectious FRNAPH genogroup II was more abundant in sediment than in surface water. Its capability to sorb to sediment and its higher persistence in the environment compared to genogroups III and IV were the two main explanations. Together, lab and in situ experiments produce an overall vision of the mechanisms governing FRNAPH distribution among the water column and riverbed sediment.

Impact of manure-related DOM on sulfonamide transport in arable soils

Field application of livestock manure introduces colloids and veterinary antibiotics, e.g. sulfonamides (SAs), into farmland. The presence of manure colloids may potentially intensify the SAs-pollution to soils and groundwater by colloid-facilitated transport. Transport of three SAs, sulfadiazine (SDZ), sulfamethoxypyridazine (SMPD), and sulfamoxole (SMOX), was investigated in saturated soil columns with and without manure colloids from sows and farrows, weaners, and fattening pigs. Experimental results showed that colloid-facilitated transport of SMOX was significant in the presence of manure colloids from fattening pigs with low C/N ratio, high SUVA280nm and protein C, while manure colloids from sows and farrows and weaners had little effect on SMOX transport. In contrast, only retardation was observed for SDZ and SMPD when manure colloids were present. Breakthrough curves (BTCs) of colloids and SAs were replicated well by a newly developed numerical model that considers colloid-filtration theory, competitive kinetic sorption, and co-transport processes. Model results demonstrate that mobile colloids act as carriers for SMOX, while immobile colloids block SMOX from sorbing onto the soil. The low affinity of SMOX to sorb on immobile colloids prevents aggregation and also promotes SMOX's colloid-facilitated transport. Conversely, the high affinity of SDZ and SMPD to sorb on all types of immobile colloids retarded their transport. Thus, manure properties play a fundamental role in increasing the leaching risk of hydrophobic sulfonamides.

Cotransport of clay colloids and viruses through water-saturated vertically oriented columns packed with glass beads: Gravity effects

The cotransport of clay colloids and viruses in vertically oriented laboratory columns packed with glass beads was investigated. Bacteriophages MS2 and ΦX174 were used as model viruses, and kaolinite (ΚGa-1b) and montmorillonite (STx-1b) as model clay colloids. A steady flow rate of Q=1.5 mL/min was applied in both vertical up (VU) and vertical down (VD) flow directions. In the presence of KGa-1b, estimated mass recovery values for both viruses were higher for VD than VU flow direction, while in the presence of STx-1b the opposite was observed. However, for all cases examined, the produced mass of viruses attached onto suspended clay particles were higher for VD than VU flow direction, suggesting that the flow direction significantly influences virus attachment onto clays, as well as packed column retention of viruses attached onto suspended clays. KGa-1b hindered the transport of ΦX174 under VD flow, while STx-1b facilitated the transport of ΦX174 under both VU and VD flow directions. Moreover, KGa-1b and STx-1b facilitated the transport of MS2 in most of the cases examined except of the case where KGa-1b was present under VD flow. Also, the experimental data were used for the estimation of virus surface-coverages and virus surface concentrations generated by virus diffusion-limited attachment, as well as virus attachment due to sedimentation. Both sedimentation and diffusion limited virus attachment were higher for VD than VU flow, except the case of MS2 and STx-1b cotransport. The diffusion-limited attachment was higher for MS2 than ΦΧ174 for all cases examined.

Interaction Between Graphene Oxide Nanoparticles and Quartz Sand

In this study, the influence of pH, ionic strength (IS), and temperature on graphene oxide (GO) nanoparticles attachment onto quartz sand were investigated. Batch experiments were conducted at three controlled temperatures (4, 12, and 25 °C) in solutions with different pH values (pH 4, 7, and 10), and ionic strengths (IS = 1.4, 6.4, and 21.4 mM), under static and dynamic conditions. The surface properties of GO nanoparticles and quartz sand were evaluated by electrophoretic mobility measurements. Derjaguin-Landau-Verwey-Overbeek (DLVO) potential energy profiles were constructed for the experimental conditions, using measured zeta potentials. The experimental results showed that GO nanoparticles were very stable under the experimental conditions. Both temperature and pH did not play a significant role in the attachment of GO nanoparticles onto quartz sand. In contrast, IS was shown to influence attachment. The attachment of GO particles onto quartz sand increased significantly with increasing IS. The experimental data were fitted nicely with a Freundlich isotherm, and the attachment kinetics were satisfactorily described with a pseudo-second-order model, which implies that the quartz sand exhibited substantial surface heterogeneity and that GO retention was governed by chemisorption. Furthermore, thermodynamic analysis revealed that the attachment process was nonspontaneous and endothermic, which may be associated with structural changes of the sand surfaces due to chemisorption. Therefore, secondary minimum interaction may not be the dominant mechanism for GO attachment onto the quartz sand under the experimental conditions.

Modeling the co-transport of viruses and colloids in unsaturated porous media

A mathematical model is developed to simulate the co-transport of viruses and colloids in unsaturated porous media under steady-state flow conditions. The virus attachment to the mobile and immobile colloids is described using a linear reversible kinetic model. Colloid transport is assumed to be decoupled from virus transport; that is, we assume that colloids are not affected by the presence of attached viruses on their surface. The governing equations are solved numerically using an alternating three-step operator splitting approach. The model is verified by fitting three sets of experimental data published in the literature: (1) Syngouna and Chrysikopoulos (2013) and (2) Walshe et al. (2010), both on the co-transport of viruses and clay colloids under saturated conditions, and (3) Syngouna and Chrysikopoulos (2015) for the co-transport of viruses and clay colloids under unsaturated conditions. We found a good agreement between observed and fitted breakthrough curves (BTCs) under both saturated and unsaturated conditions. Then, the developed model was used to simulate the co-transport of viruses and colloids in porous media under unsaturated conditions, with the aim of understanding the relative importance of various processes on the co-transport of viruses and colloids in unsaturated porous media. The virus retention in porous media in the presence of colloids is greater during unsaturated conditions as compared to the saturated conditions due to: (1) virus attachment to the air-water interface (AWI), and (2) co-deposition of colloids with attached viruses on its surface to the AWI. A sensitivity analysis of the model to various parameters showed that the virus attachment to AWI is the most sensitive parameter affecting the BTCs of both free viruses and total mobile viruses and has a significant effect on all parts of the BTC. The free and the total mobile viruses BTCs are mainly influenced by parameters describing virus attachment to the AWI, virus interaction with mobile and immobile colloids, virus attachment to solid-water interface (SWI), and colloid interaction with SWI and AWI. The virus BTC is relatively insensitive to parameters describing the maximum adsorption capacity of the AWI for colloids, inlet colloid concentration, virus detachment rate coefficient from the SWI, maximum adsorption capacity of the AWI for viruses and inlet virus concentration.

DLVO and XDLVO calculations for bacteriophage MS2 adhesion to iron oxide particles

In this study, batch experiments were performed to examine the adhesion of bacteriophage MS2 to three iron oxide particles (IOP1, IOP2 and IOP3) with different particle properties. The characteristics of MS2 and iron oxides were analyzed using various techniques to construct the classical DLVO and XDLVO potential energy profiles between MS2 and iron oxides. X-ray diffractometry peaks indicated that IOP1 was mainly composed of maghemite (γ-Fe2O3), but also contained some goethite (α-FeOOH). IOP2 was composed of hematite (α-Fe2O3) and IOP3 was composed of iron (Fe), magnetite (Fe3O4) and iron oxide (FeO). Transmission electron microscope images showed that the primary particle size of IOP1 (γ-Fe2O3) was 12.3±4.1nm. IOP2 and IOP3 had primary particle sizes of 167±35nm and 484±192nm, respectively. A surface angle analyzer demonstrated that water contact angles of IOP1, IOP2, IOP3 and MS2 were 44.83, 64.00, 34.33 and 33.00°, respectively. A vibrating sample magnetometer showed that the magnetic saturations of IOP1, IOP2 and IOP3 were 176.87, 17.02 and 946.85kA/m, respectively. Surface potentials measured in artificial ground water (AGW; 0.075mM CaCl2, 0.082mM MgCl2, 0.051mM KCl, and 1.5mM NaHCO3; pH7.6) indicated that iron oxides and MS2 were negatively charged in AGW (IOP1=-0.0185V; IOP2=-0.0194V; IOP3=-0.0301V; MS2=-0.0245V). Batch experiments demonstrated that MS2 adhesion to iron oxides was favorable in the order of IOP1>IOP2>IOP3. This tendency was well predicted by the classical DLVO model. In the DLVO calculations, both the sphere-plate and sphere-sphere geometries predicted the same trend of MS2 adhesion to iron oxides. Additionally, noticeable differences were not found between the DLVO and XDLVO interaction energy profiles, indicating that hydrophobic interactions did not play a major role; electrostatic interactions, however, did influence MS2 adhesion to iron oxides. Furthermore, the aggregation of iron oxides was investigated with a modified XDLVO model. This model included magnetic interactions between the particles in order to predict the aggregation of iron oxides. Even though iron oxide particle aggregation could occur under experimental conditions, the DLVO model results using primary particle size were more suitable for the interactions between MS2 and the iron oxides because of fast sorption of MS2 onto the surfaces of iron oxides.

Modeling dense-colloid and virus cotransport in three-dimensional porous media

A three-dimensional numerical model was developed to investigate the simultaneous transport (cotransport) of dense colloids and viruses in homogeneous, water saturated, porous media with horizontal uniform flow. The dense colloids are assumed to exist in two different phases: suspended in the aqueous phase, and attached reversibly and/or irreversibly onto the solid matrix. The viruses are assumed to exist in four different phases: suspended in aqueous phase, attached onto the solid matrix, attached onto suspended colloids, and attached onto colloids already attached onto the solid matrix. The viruses in each of the four phases are assumed to undergo inactivation with different rates. Moreover, the suspended dense colloids as well as viruses attached onto suspended dense colloids are assumed to exhibit a "restricted" settling velocity as a consequence of the gravitational force; whereas, viruses due to their small sizes and densities are assumed to have negligible "restricted" settling velocity. The governing differential equations were solved numerically with the finite difference schemes, implicitly or explicitly implemented. Model simulations have shown that the presence of dense colloid particles can either enhance or hinder the horizontal transport of viruses, but also can increase the vertical migration of viruses.

Individual and Co Transport Study of Titanium Dioxide NPs and Zinc Oxide NPs in Porous Media

The impact of pH and ionic strength on the mobility (individual and co-transport) and deposition kinetics of TiO2 and ZnO NPs in porous media was systematically investigated in this study. Packed column experiments were performed over a series of environmentally relevant ionic strengths with both NaCl (0.1-10 mM) and CaCl2 (0.01-0.1mM) solutions and at pH 5, 7, and 9. The transport of TiO2 NPs at pH 5 was not significantly affected by ZnO NPs in solution. At pH 7, a decrease in TiO2 NP transport was noted with co-existence of ZnO NPs, while at pH 9 an increase in the transport was observed. At pH 5 and 7, the transport of ZnO NPs was decreased when TiO2 NPs was present in the solution, and at pH 9, an increase was noted. The breakthrough curves (BTC) were noted to be sensitive to the solution chemistries; the decrease in the breakthrough plateau with increasing ionic strength was observed under all examined pH (5, 7, and 9). The retention profiles were the inverse of the plateaus of BTCs, as expected from mass balance considerations. Overall, the results from this study suggest that solution chemistries (ionic strength and pH) are likely the key factors that govern the individual and co-transport behavior of TiO2 and ZnO NPs in sand.

Interaction of human adenoviruses and coliphages with kaolinite and bentonite

Human adenoviruses (hAdVs) are pathogenic viruses responsible for public health problems worldwide. They have also been used as viral indicators in environmental systems. Coliphages (e.g., MS2, ΦX174) have also been studied as indicators of viral pollution in fecally contaminated water. Our objective was to evaluate the distribution of three viral fecal indicators (hAdVs, MS2, and ΦΧ174), between two different phyllosilicate clays (kaolinite and bentonite) and the aqueous phase. A series of static and dynamic experiments were conducted under two different temperatures (4, 25°C) for a time period of seven days. HAdV adsorption was examined in DNase I reaction buffer (pH=7.6, and ionic strength (IS)=1.4mM), whereas coliphage adsorption in phosphate buffered saline solution (pH=7, IS=2mM). Moreover, the effect of IS on hAdV adsorption under static conditions was evaluated. The adsorption of hAdV was assessed by real-time PCR and its infectivity was tested by cultivation methods. The coliphages MS2 and ΦΧ174 were assayed by the double-layer overlay method. The experimental results have shown that coliphage adsorption onto both kaolinite and bentonite was higher for the dynamic than the static experiments; whereas hAdV adsorption was lower under dynamic conditions. The adsorption of hAdV increased with decreasing temperature, contrary to the results obtained for the coliphages. This study examines the combined effect of temperature, agitation, clay type, and IS on hAdV adsorption onto clays. The results provide useful new information on the effective removal of viral fecal indicators (MS2, ΦX174 and hAdV) from dilute aqueous solutions by adsorption onto kaolinite and bentonite. Factors enabling enteric viruses to penetrate soils, groundwater and travel long distances within aquifers are important public health issues. Because the observed adsorption behavior of surrogate coliphages MS2 and ΦΧ174 is substantially different to that of hAdV, neither MS2 nor ΦΧ174 is recommended as a suitable model for adenovirus.

Transport of Human Adenoviruses in Water Saturated Laboratory Columns

Groundwater may be contaminated with infective human enteric viruses from various wastewater discharges, sanitary landfills, septic tanks, agricultural practices, and artificial groundwater recharge. Coliphages have been widely used as surrogates of enteric viruses, because they share many fundamental properties and features. Although a large number of studies focusing on various factors (i.e. pore water solution chemistry, fluid velocity, moisture content, temperature, and grain size) that affect biocolloid (bacteria, viruses) transport have been published over the past two decades, little attention has been given toward human adenoviruses (hAdVs). The main objective of this study was to evaluate the effect of pore water velocity on hAdV transport in water saturated laboratory-scale columns packed with glass beads. The effects of pore water velocity on virus transport and retention in porous media was examined at three pore water velocities (0.39, 0.75, and 1.22 cm/min). The results indicated that all estimated average mass recovery values for hAdV were lower than those of coliphages, which were previously reported in the literature by others for experiments conducted under similar experimental conditions. However, no obvious relationship between hAdV mass recovery and water velocity could be established from the experimental results. The collision efficiencies were quantified using the classical colloid filtration theory. Average collision efficiency, α, values decreased with decreasing flow rate, Q, and pore water velocity, U, but no significant effect of U on α was observed. Furthermore, the surface properties of viruses and glass beads were used to construct classical DLVO potential energy profiles. The results revealed that the experimental conditions of this study were unfavorable to deposition and that no aggregation between virus particles is expected to occur. A thorough understanding of the key processes governing virus transport is pivotal for public health protection.

Influence of clay particles on the transport and retention of titanium dioxide nanoparticles in quartz sand

This study investigated the influence of two representative suspended clay particles, bentonite and kaolinite, on the transport of titanium dioxide nanoparticles (nTiO2) in saturated quartz sand in both NaCl (1 and 10 mM ionic strength) and CaCl2 solutions (0.1 and 1 mM ionic strength) at pH 7. The breakthrough curves of nTiO2 with bentonite or kaolinite were higher than those without the presence of clay particles in NaCl solutions, indicating that both types of clay particles increased nTiO2 transport in NaCl solutions. Moreover, the enhancement of nTiO2 transport was more significant when bentonite was present in nTiO2 suspensions relative to kaolinite. Similar to NaCl solutions, in CaCl2 solutions, the breakthrough curves of nTiO2 with bentonite were also higher than those without clay particles, while the breakthrough curves of nTiO2 with kaolinite were lower than those without clay particles. Clearly, in CaCl2 solutions, the presence of bentonite in suspensions increased nTiO2 transport, whereas, kaolinite decreased nTiO2 transport in quartz sand. The attachment of nTiO2 onto clay particles (both bentonite and kaolinite) were observed under all experimental conditions. The increased transport of nTiO2 in most experimental conditions (except for kaolinite in CaCl2 solutions) was attributed mainly to the clay-facilitated nTiO2 transport. The straining of larger nTiO2-kaolinite clusters yet contributed to the decreased transport (enhanced retention) of nTiO2 in divalent CaCl2 solutions when kaolinite particles were copresent in suspensions.

Effect of gravity on colloid transport through water-saturated columns packed with glass beads: modeling and experiments

The role of gravitational force on colloid transport in water-saturated columns packed with glass beads was investigated. Transport experiments were performed with colloids (clays: kaolinite KGa-1b, montmorillonite STx-1b). The packed columns were placed in various orientations (horizontal, vertical, and diagonal) and a steady flow rate of Q = 1.5 mL/min was applied in both up-flow and down-flow modes. All experiments were conducted under electrostatically unfavorable conditions. The experimental data were fitted with a newly developed, analytical, one-dimensional, colloid transport model. The effect of gravity is incorporated in the mathematical model by combining the interstitial velocity (advection) with the settling velocity (gravity effect). The results revealed that flow direction influences colloid transport in porous media. The rate of particle deposition was shown to be greater for up-flow than for down-flow direction, suggesting that gravity was a significant driving force for colloid deposition.

Cotransport of titanium dioxide and fullerene nanoparticles in saturated porous media

This study investigated the cotransport of titanium dioxide nanoparticles (nTiO2) and fullerene nanoparticles (nC60), two of the most widely utilized nanoparticles, in saturated quartz sand under a series of ionic strengths in NaCl solutions (0.1-10 mM) at both pH 5 and 7. Under all examined ionic strengths at pH 5, both breakthrough curves and retained profiles of nTiO2 in the copresence of nC60 were similar to those without nC60, indicating that nC60 nanoparticles copresent in suspensions did not significantly affect the transport and retention of nTiO2 in quartz sand at pH 5. In contrast, under all examined ionic strengths at pH 7, the breakthrough curves of nTiO2 in the copresence of nC60 in suspensions were higher and the retained profiles were lower than those without nC60, which demonstrated that the presence of nC60 in suspensions increased the rate of transport (decreased retention) of nTiO2 in quartz sand at pH 7. Competition of deposition sites on quartz sand surfaces by the copresence of nC60 was found to contribute to the increased nTiO2 transport at pH 7. Under all examined ionic strength conditions at both pH 5 and 7, the breakthrough curves of nC60 were reduced in the copresence of nTiO2, and the corresponding retained profiles were higher than those without nTiO2, indicating that the presence of nTiO2 decreased the transport of nC60 in quartz sand. Co-deposition of nC60 with nTiO2 in the form of nTiO2-nC60 clusters as well as the deposition of nC60 onto previously deposited nTiO2 were responsible for the increased nC60 deposition in the presence of nTiO2 at pH 5, whereas deposition of nC60 onto surfaces of predeposited nTiO2 was found to be responsible for the increased nC60 deposition at pH 7.