Hysteresis of colloid retention and release in saturated porous media during transients in solution chemistry

Experimental measurements and numerical simulations of the transport and retention of nanocrystal CdSe/ZnS quantum dots in saturated porous media: Effects of electrolytes, organic ligand, and natural organic matter

This study explores the transport and retention of CdSe/ZnS quantum dot (QD) nanoparticles in water-saturated sand columns as a function of electrolytes (Na+ and Ca2+), ionic strength, organic ligand citrate, and Suwannee River natural organic matter (SRNOM). Numerical simulations were carried out to understand the mechanisms that govern the transport and interactions of QDs in porous media and to assess how environmental parameters impact these mechanisms. An increase in the ionic strength of NaCl and CaCl2 increased QDs retention in porous media. The reduction of the electrostatic interactions screened by dissolved electrolyte ions and the increase of divalent bridging effect are the causes for this enhanced retention behavior. Citrate or SRNOM enhanced QDs transport in NaCl and CaCl2 systems by either increasing the repulsion energy barrier or inducing the steric interactions between QDs and the quartz sand collectors. A non-exponential decay characterized the retention profiles of QDs along the distance to the inlet. The modeling results indicated the four models containing the attachment, detachment, and straining terms - Model 1: M1-attachment, Model 2: M2-attachment and detachment, Model 3: M3-straining, and Model 4: M4-attachment, detachment, and straining - closely simulated the observed breakthrough curves (BTCs) but inadequately described the retention profiles.

Exploring the vertical transport of microplastics in subsurface environments: Lab-scale experiments and field evidence

Microplastics (MPs) defined as smaller 5 mm plastic particles have received increasing attention due to their global occurrence and potential toxicity. This study investigated the effects of environmental factors (rainfall intensity, 13 and 29 mm/h) and MP characteristics (morphology (fiber, flake, and film), polymer type (polypropylene (PP), polyethylene terephthalate (PET), and polystyrene (PS)) and size (100-300, 300-500, and 500-1000 μm)) on the vertical transport of MP in unsaturated soil conditions using lab-scale column experiments. Additionally, the occurrence and characteristics of MP detected in soil/sediment (total 13 samples) and groundwater samples (total 6 samples) were explored in the field study. Laboratory-scale column experiments revealed that heavy rainfall intensity (29 mm/h) increased the degree of MP vertical transport in unsaturated soil conditions and MP fibers showed the greatest vertical mobility among the various morphologies of MPs assessed. For the polymer type and size, the lighter PP polymer or the larger size of MP (500-1000 μm) showed higher mobility. In the field study, a statistical difference in MP abundance was observed depending on the population density and degree of urban development in both soil and groundwater samples. Comparing to the two different types of environmental media samples obtained from the same site, there was a significant difference in the composition of polymer types present while statistically no difference in MP abundance was observed between the two media samples (i.e., soil or sediment and groundwater). In addition, MP fibers and polyethylene (PE) were predominantly detected in our two study areas. These results suggest that various types of MP can pass through the unsaturated zone by water infiltration, even if it takes a long time to reach groundwater. Overall, we found that the degree of vertical transport of the MPs was highly sensitive to environmental conditions and MP characteristics.

Transport of nanoscale zero-valent iron in saturated porous media: Effects of grain size, surface metal oxides, and sulfidation

Knowledge of the fate and transport of nanoscale zero-valent iron (nZVI) in saturated porous media is crucial to the development of in situ remediation technologies. This work systematically compared the retention and transport of carboxymethyl cellulose (CMC) modified nZVI (CMC-nZVI) and sulfidated nZVI (CMC-S-nZVI) particles in saturated columns packed with quartz sand of various grain sizes and different surface metal oxide coatings. Grain size reduction had an inhibitory effect on the transport of CMC-S-nZVI and CMC-nZVI due to increasing immobile zone deposition and straining in the columns. Metal oxide coatings had minor effect on the transport of CMC-S-nZVI and CMC-nZVI because the sand surface was coated by the free CMC in the suspensions, reducing the electrostatic attraction between the nZVI and surface metal oxides. CMC-S-nZVI displayed greater breakthrough (C/C₀ = 0.82-0.90) and higher mass recovery (84.9%-89.3%) than CMC-nZVI (C/C₀ = 0.70-0.80 and mass recovery = 70.9%-79.6%, respectively) under the same experimental conditions. A mathematical model based on the advection-dispersion equation simulated the experimental data of nZVI breakthrough curves very well. Findings of this study suggest sulfidation could enhance the transport of CMC-nZVI in saturated porous media with grain and surface heterogeneities, promoting its application in situ remediation.

Experimental and numerical investigation of the effect of temporal variation in ionic strength on colloid retention and remobilization in saturated porous media

Temporal variations in the chemistry of infiltrating water into the subsurface are known to cause remobilization of colloids from the grain surfaces, thereby increasing the travel distance of the colloidal contaminants. Hence, it is essential to thoroughly understand the transport, deposition, and release mechanisms of colloids in the subsurface, through laboratory experiments and modeling. There are only a few experiments in which the chemistry of inflow water is changed rapidly during colloid transport. Also, although some models have been presented for simulating the effect of transient chemistry on the fate of colloids, there is no consensus in this regard, as the proposed models suffer from shortcomings. In this study, we systematically investigated the effect of temporal variations in ionic strength on the remobilization of deposited colloids in saturated porous media through laboratory column experiments and numerical modeling. Four sets of column experiments were performed, in which we injected carboxylate-modified latex colloids at a given ionic strength for a specified period. After breakthrough of colloids, the ionic strength of inflowing water was decreased in a stepwise manner to 0 mM (DI water). The initial ionic strength values of the four experiments were 100, 50, 25, and 10 mM. We observed partial release of deposited colloids after several steps of ionic strength decrease with significant release observed only when the ionic strength was reduced to below 10 mM. We also found that the fraction of released colloids decreased with increasing value of initial ionic strength of inflow water. We have developed a mathematical model incorporating a novel formulation for ionic strength-dependent deposition and release. The model is found to capture the colloid breakthrough curves reasonably well for all experiments with the same set of parameter values, except the one at the initial ionic strength of 25 mM.

Host-Endosymbiont Relationship Impacts the Retention of Bacteria-Containing Amoeba Spores in Porous Media

Amoebae are protists that are commonly found in water, soil, and other habitats around the world and have complex interactions with other microorganisms. In this work, we investigated how host-endosymbiont interactions between amoebae and bacteria impacted the retention behavior of amoeba spores in porous media. A model amoeba species, Dictyostelium discoideum, and a representative bacterium, Burkholderia agricolaris B1qs70, were used to prepare amoeba spores that carried bacteria. After interacting with B. agricolaris, the retention of D. discoideum spores was enhanced compared to noninfected spores. Diverse proteins, especially proteins contributing to the looser exosporium structure and cell adhesion functionality, are secreted in higher quantities on the exosporium surface of infected spores compared to that of noninfected ones. Comprehensive examinations using a quartz crystal microbalance with dissipation (QCM-D), a parallel plate chamber, and a single-cell force microscope present coherent evidence that changes in the exosporium of D. discoideum spores due to infection by B. agricolaris enhance the connections between spores in the suspension and the spores that were previously deposited on the collector surface, thus resulting in more retention compared to the uninfected ones in porous media. This work provides novel insight into the retention of amoeba spores after bacterial infection in porous media and suggests that the host-endosymbiont relationship regulates the fate of biocolloids in drinking water systems, groundwater, and other porous environments.

Significant Mobility of Novel Heteroaggregates of Montmorillonite Microparticles with Nanoscale Zerovalent Irons in Saturated Porous Media

This study conducted laboratory column experiments to systematically examine the transport of novel heteroaggregates of montmorillonite (Mt) microparticles with nanoscale zerovalent irons (nZVIs) in saturated sand at solution ionic strengths (ISs) ranging from 0.001 to 0.2 M. Spherical nZVIs were synthesized using the liquid phase reduction method and were attached on the plate-shaped Mt surfaces in monolayer. While complete deposition occurred for nZVIs in sand, significant transport was observed for Mt-nZVI heteroaggregates at IS ≤ 0.01 M despite the transport decrease with an increasing loading concentration of nZVIs on Mt. The increased mobility of Mt-nZVI heteroaggregates was because the attractions between nZVIs and sand collectors were reduced by the electrostatic repulsions between the Mt and the collector surfaces, which led to a decreased deposition in the sand columns. Complete deposition occurred for the Mt-nZVI heteroaggregates at IS ≥ 0.1 M due to a favorable deposition at Derjaguin–Landau–Verwey–Overbeek (DLVO) primary energy minima. Interestingly, a large fraction of the deposited heteroaggregates was released by reducing IS because of a monotonic decrease of interaction energy with separation distance for the heteroaggregates at low ISs (resulting in repulsive forces), in contrast to the irreversible deposition of nZVIs. Therefore, the fabricated heteroaggregates could also have high mobility in subsurfaces with saline pore water through continuous capture and release using multiple injections of water with low ISs. Our study was the first to examine the transport of heteroaggregates of a plate-like particle with spherical nanoparticles in porous media; the results have important implications in the use of nanoscale zerovalent iron for in situ soil and groundwater remediation.

Retention of graphene oxide and reduced graphene oxide in porous media: Diffusion-attachment, interception-attachment and straining

The aggregation, deposition and retention of graphene oxide (GO) and reduced graphene oxide (RGO) were investigated systematically to estimate their mobility in the environment. RGO aggregates faster than GO, resulting in weaker diffusive transfer and a lower deposition rate on oxide surfaces. In NaCl, the critical deposition concentration of RGO (CDC_RGO) is smaller than CDC_GO on the SiO₂ surface, indicating that RGO achieves favorable deposition at lower ionic strength. In CaCl₂, Ca²⁺ bridging causes close CDC_GO and CDC_RGO. The retention process was observed in the photolithographic SiO₂ and Al₂O₃ micromodels. GO and RGO particles approach collectors mainly via interception before attachment. The interactive forces have a limited effect on the particle retention. The larger RGO aggregates cause greater extent interception and straining, resulting in lower mobility than GO in porous media. The mobility of GO and RGO show different trends in quartz crystal microbalance with dissipation (QCM-D) and in micromodels because the interception and straining mechanisms exist in pore space. Micromodel observation confirms the processes of interception and straining. The combination of QCM-D and micromodel experiments provides the connection of diffusion-attachment, interception-attachment and straining, which comprehensively explains the higher mobility of GO than RGO in porous media.

Effects of colloids on ammonia nitrogen release under different ion conditions in natural sediments of Lake Taihu, China

Colloids act as vectors for accelerating contaminant movement in natural porous media such as lake sediments. Releasing characteristics of colloids and colloid-adsorbed ammonia nitrogen from lake sediments with the presence of monovalent and divalent cations were studied by an indoor anaerobic flooding incubation experiment. Results show that release of colloids was influenced by valence state and strength of the cation. In the presence of Na⁺ and Ca²⁺, the concentrations of the colloids in the overlying water were 11-163 mg/L and 13-88 mg/L during the entire incubation process and their values increased by 12.7%-122.3% and 1.5%-29.1%, respectively, compared to the control group, which indicated that the promoting effect of monovalent cations on release of colloids was more obvious than that of divalent cations. However, the total mass of colloids release reduced with the increasing ionic strength. Colloid-adsorbed ammonia nitrogen in the overlying water reached 0.15-1.72 mg/L and 0.15-1.12 mg/L with the presence of Na⁺ and Ca²⁺ and was higher 61.7%-161.7% and 21.3%-80.9%, respectively, than in the control group, indicating a consequent effect of ion conditions on the release of ammonium nitrogen from sediments. A significant positive correlation between colloids and ammonia nitrogen concentration further shows that colloidal activity determinately resulted in the increase or decrease in the ammonia nitrogen concentration in the overlying water, which could adsorb ammonia nitrogen and act as vehicles to carry ammonia nitrogen together into the aqueous medium or sink into the sediment. The release of ammonia nitrogen is possibly enhanced by colloidal behavior and varies with spatiotemporal ionic conditions in natural sediments. These findings are essential for improving the understanding of the geological fate of environmental colloids and associated nutritive salts, which provide scientific basis and technical support for the control of endogenous pollution and the comprehensive treatment of water bodies in Lake Taihu.

Review on physical and chemical factors affecting fines migration in porous media

Permeability reduction and formation damage in porous media caused by fines (defined as unconfined solid particles present in the pore spaces) migration is one of the major reasons for productivity decline. It is well accepted that particle detachment occurs under imbalanced torques arising from hydrodynamic and adhesive forces exerted on attached particles. This paper reviewed current understanding on primary factors influencing fines migration as well as mathematical formulations for quantification. We also introduced salinity-related experimental observations that contradict theoretical predictions based on torque balance criteria, such as delayed particle release and attachment-detachment hysteresis. The delay of particle release during low-salinity water injection was successfully explained and formulated by the Nernst-Planck diffusion of ions in a narrow contact area. In addition to the widely recognized explanation by surface heterogeneity and the presence of low-velocity regions, we proposed a hypothesis that accounts for the shifting of equilibrium positions, providing new insight into the interpretation of elusive attachment-detachment hysteresis both physically and mathematically. The review was finalized by discussing the quantification of anomalous salinity effect on adhesion force at low- and high-salinity conditions.

The transport and retention of CQDs-doped TiO2 in packed columns and 3D printed micromodels

CQDs-doped TiO₂ (C-TiO₂) has drawn increased attention in recent because of its excellent catalytic performance. Understanding the transport of C-TiO₂ in porous media is necessary for evaluating the environmental process of this new nanomaterial. Column experiments were used in this study to investigate ionic strength (IS), dissolved organic matter (DOM) and sand grain size on the transport of C-TiO₂. The mobility of C-TiO₂ was inhibited by the increased IS and decreased sand grain size, but was promoted by the increased DOM concentration. The promotion efficiency of DOM ranked as humic acid (HA) > alginate (Alg) > bovine serum albumin (BSA), which was in the same order as their ability to change surface charges. The micromodels of pore network were prepared via 3D printing to further reveal the deposition mechanisms and spatial/temporal distribution of C-TiO₂ in porous space. C-TiO₂ mainly attached to the upstream region of collectors because of interception. The collector ripening was observed after long-time deposition. The existence of DOM caused visible decrease of C-TiO₂ deposition in the pore network. HA caused the most remarkable reduce of deposition in the three types of DOM, which was consistent with the column experiment results. This research is helpful to predict the transport of C-TiO₂ in natural porous media.

Bacteria transport and deposition in an unsaturated aggregated porous medium with dual porosity

Bacterial transport and deposition play an important role in the assessment and prediction of subsurface pollution risks. Bacteria transport experiments were performed under unsaturated flow conditions in an aggregated porous medium at the laboratory column scale, to investigate how the inter- and intra-aggregated pore space of this medium could affect transport and deposition under unsaturated flow conditions, where inter- and intra-pore spaces are not fully activated. The results obtained through experimental observations and numerical simulations showed that some intra- and inter-pore space of this medium was excluded from bacteria transport and retention, as confirmed by the non-uniform transport of bacteria pathways in the aggregated porous media under unsaturated flow conditions. Capillary energy was higher the than other forces acting at bacteria air-water-solid interfaces. If this energy should contribute in increasing bacteria deposition under unsaturated conditions, similar to what has been reported for sandy media, similar overall retention of E. coli and R. rhodochrous was obtained under unsaturated flow conditions, suggesting that capillary energy was not the driving force for bacteria deposition.

Experimental measurements and numerical simulations of the transport and retention of nanocrystal CdSe/ZnS quantum dots in saturated porous media: effects of pH, organic ligand, and natural organic matter

The risks of environmental exposures of quantum dot (QD) nanoparticles are increasing, but these risks are difficult to assess because fundamental questions remain about factors affecting the mobility of QDs. The objective of this study is to help address this shortcoming by evaluating the physico-chemical mechanisms controlling the transport and retention of CdSe/ZnS QDs under various environmental conditions. The approach was to run a series of laboratory-scale column experiments where QDs were transported through saturated porous media with different pH values and concentrations of citrate and Suwannee River natural organic matter (SRNOM). Numerical simulations were then conducted and compared with the laboratory data in order to evaluate parameters controlling transport. QD suspensions were injected into the column in an upward direction and ICP-MS used to analyze Cd²⁺ concentrations (C) in column effluent and sand porous media samples. The increase in the background solution pH values enhanced the QD transport and decreased the QD retention. QD transport recovery percentages obtained from the column effluent samples were 2.6%, 83.2%, 101.7%, 96.5%, and 98.9%, at pH levels of 1.5, 3.5, 5, 7, and 9, respectively. The effects of citrate and SRNOM on the transport and retention of QDs were pH dependent as reflected in the influence of the electrostatic and steric interactions between QDs and sand surfaces. QDs were mobile under unfavorable deposition conditions at environmentally relevant pHs (i.e., 5, 7, and 9). Under favorable pH conditions for deposition (i.e., 1.5), QDs were completely retained within the porous media. The retention profiles of QDs showed a non-exponential decay with distance to the inlet, attributed to multiple deposition rates caused by the QD particles and surface heterogeneities of the quartz silica sand. Results of the diameter ratios of QDs to the median sand grains, in suspensions of DI water at pH 1.5, of citrate at pH 1.5, and of citrate at pH 3.5 indicate straining as the dominating mechanism for QD retention in porous media. The blocking effect and straining were significant under favorable deposition conditions and the detachment effect was non-negligible under unfavorable deposition conditions. Physico-chemical attachment and straining are the governing mechanisms that control the retention of QDs. Overall, experimental results indicate that aggregation, deposition, straining, blocking, and DLVO-type interactions affect the advective transport and retention of QDs in saturated porous media. The simulations were conducted using models that include terms describing attachment, detachment, and straining terms-model 1: M1-attachment, model 2: M2-attachment and detachment, model 3: M3-straining, and model 4: M4-attachment, detachment, and straining. The results from simulations with M2-attachment and detachment and M4-attachment, detachment, and straining matched best the observed breakthrough curves, but all four models inadequately described the retention profiles. Our findings demonstrate that QDs are mobile in porous media under a wide range of physico-chemical conditions representative of the natural environment. The mobility behavior of QDs in porous media indicated the potential risk of soil and groundwater contamination.

Transport and retention of engineered silver nanoparticles in carbonate-rich sediments in the presence and absence of soil organic matter

The transport and retention behavior of polymer- (PVP-AgNP) and surfactant-stabilized (AgPURE) silver nanoparticles in carbonate-dominated saturated and unconsolidated porous media was studied at the laboratory scale. Initial column experiments were conducted to investigate the influence of chemical heterogeneity (CH) and nano-scale surface roughness (NR) arising from mixtures of clean, positively charged calcium carbonate sand (CCS), and negatively charged quartz sands. Additional column experiments were performed to elucidate the impact of CH and NR arising from the presence and absence of soil organic matter (SOM) on a natural carbonate-dominated aquifer material. The role of the nanoparticle capping agent was examined under all conditions tested in the column experiments. Nanoparticle transport was well described using a numerical model that facilitated blocking on one or two retention sites. Results demonstrate that an increase in CCS content in the artificially mixed porous medium leads to delayed breakthrough of the AgNPs, although AgPURE was much less affected by the CCS content than PVP-AgNPs. Interestingly, only a small portion of the solid surface area contributed to AgNP retention, even on positively charged CCS, due to the presence of NR which weakened the adhesive interaction. The presence of SOM enhanced the retention of AgPURE on the natural carbonate-dominated aquifer material, which can be a result of hydrophobic or hydrophilic interactions or due to cation bridging. Surprisingly, SOM had no significant impact on PVP-AgNP retention, which suggests that a reduction in electrostatic repulsion due to the presence of SOM outweighs the relative importance of other binding mechanisms. Our findings are important for future studies related to AgNP transport in shallow unconsolidated calcareous and siliceous sands.

Deposition and release of carboxylated graphene in saturated porous media: Effect of transient solution chemistry

Chemical perturbation of pore-water in porous media may remobilize and release deposited colloids/nanomaterials into bulk flow. This re-entrainment process is important to accurately assessing the fate and transport of colloids/nanomaterials in the subsurface. This study investigated deposition and subsequent release of carboxylated graphene nanomaterials (CG) in water-saturated sand columns by first depositing CG at 100 mM NaCl or 2 mM CaCl₂ (Phase 1), followed by Phase 2 (elution with sequences of 50, 10, and 1 mM NaCl, or sequences of 0.5 and 0.1 mM CaCl₂), and then Phase 3 elution using deionized water. Approximate 89.2%-98.7% of injected CG was retained in sand through Derjaguin-Landau-Verwey-Overbeek (DLVO) interactions, Ca²⁺ bridging, and straining in Phase 1. Sequential reduction of ionic strength in Phases 2 and 3 resulted in increased release of deposited CG mainly due to the expansion of the electrical double layer thickness and thus decreased depth of the attractive secondary minimum. With increasing pulses of flushing solution, unrecoverable CG increased because weakly associated CG via the secondary minimum was likely translated to immobile regions. Significant tailing of CG released in Phase 3 suggests that CG retained in CaCl₂ was more resistant upon detachment than in NaCl. In cation exchange experiment, only 0.7% of applied CG was released, possibly ascribed to the CG remobilized by cation exchange was immediately re-entrained by the secondary minimum in 50 mM NaCl. Our findings indicate that retained nanomaterials (e.g., CG) can be remobilized and transported downward in transient solution chemistries, raising concerns about their potential migration risk to groundwater.

Dynamics and sources of colloids in shallow groundwater in lowland wells and fracture flow in sloping farmland

Field-scale studies of natural colloid mobilization and transport in finely fractured aquifer as well as the source identification of groundwater colloids are of great importance to the safety of shallow groundwater. In this study, the daily monitoring of fracture flow from a sloping farmland plot and the biweekly monitoring of three lowland shallow wells within the same catchment were carried out simultaneously in 2013. The effects of physicochemical perturbations on groundwater colloid dynamics were explored in detail using partial redundancy analysis, structural equation modeling, Pearson correlation and multi-linear regression analyses. The characterization and source identification of groundwater colloids were addressed via multiple parameters. The daily colloid concentration in the fracture flow varied between 0.54 and 31.90 mg/L (1.64 mg/L on average). Unique periods of high colloid concentration (5.59 mg/L on average) occurred during the initially generated flow following the dry season. In comparison, a narrower colloid concentration range of 0.24-11.66 mg/L was observed in the lowland shallow wells, with a smaller temporal variation than that of the fracture flow. A low percentage (2.4-7.0%) of colloids and a high percentage (47.7-92.0%) of coarse particles (2-10 μm) were present in the lowland well water. Hydraulic perturbation by rainwater infiltration in the sloping farmland was the dominant mechanism for colloid mobilization in general; this effect retreated to secondary importance behind chemical perturbations (pH, Mg²⁺ and DOC) at low flow discharges (<1.3 L/min). In contrast, water chemistry (e.g., EC, cations and DOC concentrations) exhibited a major effect on colloid dynamics in the water of the lowland wells, except for the extremely high-salinity water of one well, in which water temperature showed a negative dominant influence on colloid stability. The combined use of multiple parameters (e.g., mineral composition and organic matter, calcium carbonate and δ¹³C contents) traced groundwater colloids to the shallow soil in the upper farmlands. It is strongly advised that in finely fractured aquifers within agricultural catchments, not only the small colloids but also the coarse particles in the size range of 2-10 μm should be monitored in case of colloid-associated contamination from agricultural wastes e.g., N, P, pesticides and/or heavy metals, especially at the early stages of the rainy seasons.

Graphene oxide-facilitated uranium transport and release in saturated medium: Effect of ionic strength and medium structure

Natural subsurface environment is a complex heterogeneous system. To investigate the effect of ionic strength (IS) and heterogeneity on the transport and remobilization of graphene oxide (GO)-facilitated uranium (U(VI)) in saturated porous media, column experiments were performed by the injection of U(VI) alone and U(VI)+GO mixtures into homogeneous and heterogeneous porous media under low and high ionic strength (1 and 50 mM) conditions, and then the columns were successively flushed with background solution and DI water. Results showed that when U(VI) only was introduced into the columns, IS had little effect on the migration of U(VI) alone in both media and the presence of preferential flow in heterogeneous media slightly enhanced the mobility of U(VI). As U(VI)+GO mixtures were injected into the columns, GO showed strong mobility at low IS and high released peak at high IS. The appearance of GO significantly enhanced U(VI) transport in both media. Under low IS condition, the mobility of U(VI) was significantly enhanced at the injection phase, and the medium heterogeneity further promoted the amount of GO-sorbed U(VI) transport. At high IS, less GO-sorbed U(VI) was observed during injection phase, and a large amount of retained GO-sorbed U(VI) were released with GO remobilization during water flushing phase, and the release showed the longer-tailing phenomenon and the release amount was more pronounced in heterogeneous media. The findings in this study showed that the coupled effect of solution chemistry and media heterogeneity played important roles on GO-facilitated U(VI) transport and release in soil and groundwater system.

Impacts of leachate of landfill on the groundwater hydrochemistry and size distributions and heavy metal components of colloids: a case study in NE China

Colloids associated with heavy metals are ubiquitous in contaminated groundwater; waste accumulation at imperfectly sealed landfills can produce large amounts of leachate with colloids and heavy metal contaminants, which can pollute the downstream groundwater. In this study, three sites in a landfill were sampled to reveal heavy metal particle size distributions and their chemical compositions. The > 220 nm particle sizes were the predominant size in the downstream groundwater, while the < 10 nm particle sizes were the predominant size in the upstream groundwater. Total Fe increased from 35.5 μg/L in the upstream groundwater to 107 μg/L in the downstream groundwater. This increase was attributed to the enhanced migration and accumulation of colloids in the aqueous phase. The elements and the colloid size distribution in the landfill indirectly reflected the composition and degradation of the waste. Colloids played a key role in distribution of both solid particles and aqueous contaminants in the landfill. The results of this study will contribute to the knowledge of the effect of different contaminants in the vicinity of landfills without appropriate sealing systems.

Continuum-based models and concepts for the transport of nanoparticles in saturated porous media: A state-of-the-science review

Environmental applications of nanoparticles (NP) increasingly result in widespread NP distribution within porous media where they are subject to various concurrent transport mechanisms including irreversible deposition, attachment/detachment (equilibrium or kinetic), agglomeration, physical straining, site-blocking, ripening, and size exclusion. Fundamental research in NP transport is typically conducted at small scale, and theoretical mechanistic modeling of particle transport in porous media faces challenges when considering the simultaneous effects of transport mechanisms. Continuum modeling approaches, in contrast, are scalable across various scales ranging from column experiments to aquifer. They have also been able to successfully describe the simultaneous occurrence of various transport mechanisms of NP in porous media such as blocking/straining or agglomeration/deposition/detachment. However, the diversity of model equations developed by different authors and the lack of effective approaches for their validation present obstacles to the successful robust application of these models for describing or predicting NP transport phenomena. This review aims to describe consistently all the important NP transport mechanisms along with their representative mathematical continuum models as found in the current scientific literature. Detailed characterizations of each transport phenomenon in regards to their manifestation in the column experiment outcomes, i.e., breakthrough curve (BTC) and residual concentration profile (RCP), are presented to facilitate future interpretations of BTCs and RCPs. The review highlights two NP transport mechanisms, agglomeration and size exclusion, which are potentially of great importance in controlling the fate and transport of NP in the subsurface media yet have been widely neglected in many existing modeling studies. A critical limitation of the continuum modeling approach is the number of parameters used upon application to larger scales and when a series of transport mechanisms are involved. We investigate the use of simplifying assumptions, such as the equilibrium assumption, in modeling the attachment/detachment mechanisms within a continuum modelling framework. While acknowledging criticisms about the use of this assumption for NP deposition on a mechanistic (process) basis, we found that its use as a description of dynamic deposition behavior in a continuum model yields broadly similar results to those arising from a kinetic model. Furthermore, we show that in two dimensional (2-D) continuum models the modeling efficiency based on the Akaike information criterion (AIC) is enhanced for equilibrium vs kinetic with no significant reduction in model performance. This is because fewer parameters are needed for the equilibrium model compared to the kinetic model. Two major transport regimes are identified in the transport of NP within porous media. The first regime is characterized by higher particle-surface attachment affinity than particle-particle attachment affinity, and operative transport mechanisms of physicochemical filtration, blocking, and physical retention. The second regime is characterized by the domination of particle-particle attachment tendency over particle-surface affinity. In this regime although physicochemical filtration as well as straining may still be operative, ripening is predominant together with agglomeration and further subsequent retention. In both regimes careful assessment of NP fate and transport is necessary since certain combinations of concurrent transport phenomena leading to large migration distances are possible in either case.

Role of solution chemistry in the retention and release of graphene oxide nanomaterials in uncoated and iron oxide-coated sand

Understanding the fate and transport including remobilization of graphene oxide nanomaterials (GONMs) in the subsurface would enable us to expedite their benign use and evaluate their environmental impacts and health risks. In this study, the retention and release of GONMs were investigated in water-saturated columns packed with uncoated sand (Un-S) or iron oxide-coated sand (FeS) at environmentally relevant solution chemistries (1-100mM KCl and 0.1-10mM CaCl₂ at pH7 and 11). Our results showed that increasing ionic strength (IS) inhibited GONMs' transport, and the impact of K⁺ was less than Ca²⁺. The positively charged iron oxide coating on sand surfaces immobilized the negatively charged GONMs (pH7) in the primary minimum, yielding hyperexponential retention profiles particularly in Ca²⁺. A stepwise decrease in pore-water IS caused detachment of previously retained GONMs. The mass of GONMs released during each detachment step correlated positively with the difference in secondary minimum depth (ΔΦ_min2) at each IS, indicating that the released GONMs were retained in the secondary minimum. While most retained GONMs were re-entrained upon lowering pore-water IS in Un-S, decreasing IS only released limited GONMs in FeS, which were captured in the primary minimum. Introducing 1mM NaOH (pH11) released most retained GONMs in FeS; and average hydrodynamic diameters of the detached GONMs upon injecting NaOH were significantly smaller than those of GONMs in the influent and retentate, suggesting that NaOH induced GONMs disaggregation. Our findings advance current knowledge to better predict NMs' fate and transport under various solution chemistries such as during rainfall events or in the mixing zones between sea water and fresh water where transient IS changes drastically.

Colloid Mobilization in a Fractured Soil: Effect of Pore-Water Exchange between Preferential Flow Paths and Soil Matrix

Exchange of water and solutes between contaminated soil matrix and bulk solution in preferential flow paths has been shown to contribute to the long-term release of dissolved contaminants in the subsurface, but whether and how this exchange can affect the release of colloids in a soil are unclear. To examine this, we applied rainfall solutions of different ionic strength on an intact soil core and compared the resulting changes in effluent colloid concentration through multiple sampling ports. The exchange of water between soil matrix and the preferential flow paths leading to each port was characterized on the basis of the bromide (conservative tracer) breakthrough time at the port. At individual ports, two rainfalls of a certain ionic strength mobilized different amounts of colloids when the soil was pre-exposed to a solution of lower or higher ionic strength. This result indicates that colloid mobilization depended on rainfall solution history, which is referred as colloid mobilization hysteresis. The extent of hysteresis was increased with increases in exchange of pore water and solutes between preferential flow paths and matrix. The results indicate that the soil matrix exchanged the old water from the previous infiltration with new infiltrating water during successive infiltration and changed the pore water chemistry in the preferential flow paths, which in turn affected the release of soil colloids. Therefore, rainfall solution history and soil heterogeneity must be considered to assess colloid mobilization in the subsurface. These findings have implications for the release of colloids, colloid-associated contaminants, and pathogens from soils.

Critical role of surface roughness on colloid retention and release in porous media

This paper examines the critical role of surface roughness (both nano- and micro-scale) on the processes of colloid retention and release in porous media under steady-state and transient chemical conditions. Nanoscale surface roughness (NSR) in the order of a few nanometers, which is common on natural solid surfaces, was incorporated into extended-DLVO calculations to quantify the magnitudes of interaction energy parameters (e.g. the energy barrier to attachment, ΔΦa , and detachment, ΔΦd , from a primary minimum). This information was subsequently used to explain the behavior of colloid retention and release in column and batch experiments under different ionic strength (IS) and pH conditions. Results demonstrated that the density and height of NSR significantly influenced the interaction energy parameters and consequently the extent and kinetics of colloid retention and release. In particular, values of ΔΦa and ΔΦd significantly decreased in the presence of NSR. Therefore, consistent with findings of column experiments, colloid retention in the primary minimum was predicted to occur at some specific locations on the sand surface, even at low IS conditions. However, NSR yielded a much weaker primary minimum interaction compared with that of smooth surfaces. Colloid release from primary minima upon decreasing IS and increasing pH was attributed to the impact of NSR on the values of ΔΦd . Pronounced differences in the amount of colloid retention in batch and column experiments indicated that primary minimum interactions were weak even at high IS conditions. Negligible colloid retention in batch experiments was attributed to hydrodynamic torques overcoming adhesive torques, whereas significant colloid retention in column experiments was attributed to nano- and micro-scale roughness which would dramatically alter the lever arms associated with hydrodynamic and adhesive torques.

Colloid release and clogging in porous media: Effects of solution ionic strength and flow velocity

The release and retention of in-situ colloids in aquifers play an important role in the sustainable operation of managed aquifer recharge (MAR) schemes. The processes of colloid release, retention, and associated permeability changes in consolidated aquifer sediments were studied by displacing native groundwater with reverse osmosis-treated (RO) water at various flow velocities. Significant amounts of colloid release occurred when: (i) the native groundwater was displaced by RO-water with a low ionic strength (IS), and (ii) the flow velocity was increased in a stepwise manner. The amount of colloid release and associated permeability reduction upon RO-water injection depended on the initial clay content of the core. The concentration of released colloids was relatively low and the permeability reduction was negligible for the core sample with a low clay content of about 1.3%. In contrast, core samples with about 6 and 7.5% clay content exhibited: (i) close to two orders of magnitude increase in effluent colloid concentration and (ii) more than 65% permeability reduction. Incremental improvement in the core permeability was achieved when the flow velocity increased, whereas a short flow interruption provided a considerable increase in the core permeability. This dependence of colloid release and permeability changes on flow velocity and colloid concentration was consistent with colloid retention and release at pore constrictions due to the mechanism of hydrodynamic bridging. A mathematical model was formulated to describe the processes of colloid release, transport, retention at pore constrictions, and subsequent permeability changes. Our experimental and modeling results indicated that only a small fraction of the in-situ colloids was released for any given change in the IS or flow velocity. Comparison of the fitted and experimentally measured effluent colloid concentrations and associated changes in the core permeability showed good agreement, indicating that the essential physics were accurately captured by the model.

Pore-Scale Study of Flow Rate on Colloid Attachment and Remobilization in a Saturated Micromodel

Colloid attachment is an important retention mechanism. It is influenced by colloid size, pore size, and flow rate, among other factors. In this work, we studied colloid attachment experimentally under various flow rates, as well as colloid release in response to a rapid change of flow rate. Colloid transport experiments under saturated conditions and with different flow rates were conducted in a physical micromodel. The micromodel was made of polydimethylsiloxane (PDMS), which is a hydrophobic polymer. Colloids were hydrophilic fluorescent carboxylate-modified polystyrene latex microspheres with a mean diameter of 300 nm. We could directly observe the movement of colloids within the pores using a confocal microscope. We also obtained concentration breakthrough curves by measuring the fluorescence intensity at the outlet of the micromodel. In addition, our experiments were simulated using a pore-network modeling, PoreFlow, based on the pore structure of the micromodel. Local colloid concentrations were calculated by solving local mass balance equations for all network elements and then averaging resulting concentrations over the whole micromodel. The measured breakthrough curves were successfully simulated using PoreFlow. Observed and calculated breakthrough curves showed that colloid attachment rate was smaller for larger flow rate. Temporally enhance colloid release (remobilization of attached colloids) was observed when the flow rate was increased by a factor of 10. But no colloid remobilization was observed when the flow rate decreased by a factor of 10.

In situ measurement and simulation of nano-magnetite mobility in porous media subject to transient salinity

Nanotechnologies have been proposed for a variety of environmental applications, including subsurface characterization, enhanced oil recovery, and in situ contaminant remediation. For such applications, quantitative predictive models will be of great utility for system design and implementation. Electrolyte chemistry, which can vary substantially within subsurface pore waters, has been shown to strongly influence nanoparticle aggregation and deposition in porous media. Thus, it is essential that mathematical models be capable of tracking changes in electrolyte chemistry and predicting its influence on nanoparticle mobility. In this work, a modified version of a multi-dimensional multispecies transport simulator (SEAWAT) was employed to model nanoparticle transport under transient electrolyte conditions. The modeling effort was supported by experimental measurements of paramagnetic magnetite (Fe3O4) nanoparticle, coated with polyacrylamide-methylpropane sulfonic acid - lauryl acrylate (nMag-PAMPS), mobility in columns packed with 40-50 mesh Ottawa sand. Column effluent analyses and magnetic resonance imaging (MRI) were used to quantify nanoparticle breakthrough and in situ aqueous phase concentrations, respectively. Experimental observations revealed that introduction of de-ionized water into the brine saturated column (80 g L(-1) NaCl + 20 g L(-1) CaCl2) promoted release and remobilization of deposited nanoparticles along a diagonal front, coincident with the variable density flow field. This behavior was accurately captured by the simulation results, which indicated that a two-site deposition-release model provided the best fit to experimental observations, suggesting that heterogeneous nanoparticle-surface interactions governed nanoparticle attachment. These findings illustrate the importance of accounting for both physical and chemical processes associated with changes in electrolyte chemistry when predicting nanoparticle transport behavior in subsurface formations.

An explanation for differences in the process of colloid adsorption in batch and column studies

It is essential to understand the mechanisms that control virus and bacteria removal in the subsurface environment to assess the risk of groundwater contamination with fecal microorganisms. This study was conducted to explicitly provide a critical and systematic comparison between batch and column experiments. The aim was to investigate the underlying factors causing the commonly observed discrepancies in colloid adsorption process in column and batch systems. We examined the colloid adsorption behavior of four different sizes of carboxylate-modified latex (CML) microspheres, as surrogates for viruses and bacteria, on quartz sand in batch and column experiments over a wide range of solution ionic strengths (IS). Our results show that adsorption of colloids in batch systems should be considered as an irreversible attachment because the attachment/detachment model was found to be inadequate in describing the batch results. An irreversible attachment-blocking model was found to accurately describe the results of both batch and column experiments. The rate of attachment was found to depend highly on colloid size, solution IS and the fraction of the sand surface area favorable for attachment (Sf). The rate of attachment and Sf values were different in batch and column experiments due to differences in the hydrodynamic of the system, and the role of surface roughness and pore structure on colloid attachment. Results from column and batch experiments were generally not comparable, especially for larger colloids (≥0.5μm). Predictions based on classical DLVO theory were found to inadequately describe interaction energies between colloids and sand surfaces.

Effect of hydrofracking fluid on colloid transport in the unsaturated zone

Hydraulic fracturing is expanding rapidly in the US to meet increasing energy demand and requires high volumes of hydrofracking fluid to displace natural gas from shale. Accidental spills and deliberate land application of hydrofracking fluids, which return to the surface during hydrofracking, are common causes of environmental contamination. Since the chemistry of hydrofracking fluids favors transport of colloids and mineral particles through rock cracks, it may also facilitate transport of in situ colloids and associated pollutants in unsaturated soils. We investigated this by subsequently injecting deionized water and flowback fluid at increasing flow rates into unsaturated sand columns containing colloids. Colloid retention and mobilization was measured in the column effluent and visualized in situ with bright field microscopy. While <5% of initial colloids were released by flushing with deionized water, 32-36% were released by flushing with flowback fluid in two distinct breakthrough peaks. These peaks resulted from 1) surface tension reduction and steric repulsion and 2) slow kinetic disaggregation of colloid flocs. Increasing the flow rate of the flowback fluid mobilized an additional 36% of colloids, due to the expansion of water filled pore space. This study suggests that hydrofracking fluid may also indirectly contaminate groundwater by remobilizing existing colloidal pollutants.

Adhesion of bacterial pathogens to soil colloidal particles: influences of cell type, natural organic matter, and solution chemistry

Bacterial adhesion to granular soil particles is well studied; however, pathogen interactions with naturally occurring colloidal particles (<2 μm) in soil has not been investigated. This study was developed to identify the interaction mechanisms between model bacterial pathogens and soil colloids as a function of cell type, natural organic matter (NOM), and solution chemistry. Specifically, batch adhesion experiments were conducted using NOM-present, NOM-stripped soil colloids, Streptococcus suis SC05 and Escherichia coli WH09 over a wide range of solution pH (4.0-9.0) and ionic strength (IS, 1-100 mM KCl). Cell characterization techniques, Freundlich isotherm, and Derjaguin-Landau-Verwey-Overbeek (DLVO) theory (sphere-sphere model) were utilized to quantitatively determine the interactions between cells and colloids. The adhesion coefficients (Kf) of S. suis SC05 to NOM-present and NOM-stripped soil colloids were significantly higher than E. coli WH09, respectively. Similarly, Kf values of S. suis SC05 and E. coli WH09 adhesion to NOM-stripped soil colloids were greater than those colloids with NOM-present, respectively, suggesting NOM inhibits bacterial adhesion. Cell adhesion to soil colloids declined with increasing pH and enhanced with rising IS (1-50 mM). Interaction energy calculations indicate these adhesion trends can be explained by DLVO-type forces, with S. suis SC05 and E. coli WH09 being weakly adhered in shallow secondary energy minima via polymer bridging and charge heterogeneity. S. suis SC05 adhesion decreased at higher IS 100 mM, which is attributed to the change of hydrophobic effect and steric repulsion resulted from the greater presence of extracellular polymeric substances (EPS) on S. suis SC05 surface as compared to E. coli WH09. Hence, pathogen adhesion to the colloidal material is determined by a combination of DLVO, charge heterogeneity, hydrophobic and polymer interactions as a function of solution chemistry.

Does water content or flow rate control colloid transport in unsaturated porous media?

Mobile colloids can play an important role in contaminant transport in soils: many contaminants exist in colloidal form, and colloids can facilitate transport of otherwise immobile contaminants. In unsaturated soils, colloid transport is, among other factors, affected by water content and flow rate. Our objective was to determine whether water content or flow rate is more important for colloid transport. We passed negatively charged polystyrene colloids (220 nm diameter) through unsaturated sand-filled columns under steady-state flow at different water contents (effective water saturations Se ranging from 0.1 to 1.0, with Se = (θ - θr)/(θs - θr)) and flow rates (pore water velocities v of 5 and 10 cm/min). Water content was the dominant factor in our experiments. Colloid transport decreased with decreasing water content, and below a critical water content (Se < 0.1), colloid transport was inhibited, and colloids were strained in water films. Pendular ring and water film thickness calculations indicated that colloids can move only when pendular rings are interconnected. The flow rate affected retention of colloids in the secondary energy minimum, with less colloids being trapped when the flow rate increased. These results confirm the importance of both water content and flow rate for colloid transport in unsaturated porous media and highlight the dominant role of water content.

Modeling microorganism transport and survival in the subsurface

An understanding of microbial transport and survival in the subsurface is needed for public health, environmental applications, and industrial processes. Much research has therefore been directed to quantify mechanisms influencing microbial fate, and the results demonstrate a complex coupling among many physical, chemical, and biological factors. Mathematical models can be used to help understand and predict the complexities of microbial transport and survival in the subsurface under given assumptions and conditions. This review highlights existing model formulations that can be used for this purpose. In particular, we discuss models based on the advection-dispersion equation, with terms for kinetic retention to solid-water and/or air-water interfaces; blocking and ripening; release that is dependent on the resident time, diffusion, and transients in solution chemistry, water velocity, and water saturation; and microbial decay (first-order and Weibull) and growth (logistic and Monod) that is dependent on temperature, nutrient concentration, and/or microbial concentration. We highlight a two-region model to account for microbe migration in the vicinity of a solid phase and use it to simulate the coupled transport and survival of species under a variety of environmentally relevant scenarios. This review identifies challenges and limitations of models to describe and predict microbial transport and survival. In particular, many model parameters have to be optimized to simulate a diversity of observed transport, retention, and survival behavior at the laboratory scale. Improved theory and models are needed to predict the fate of microorganisms in natural subsurface systems that are highly dynamic and heterogeneous.

Facilitated attachment of nanoparticles at primary minima by nanoscale roughness is susceptible to hydrodynamic drag under unfavorable chemical conditions

This study investigated effects of flow velocity on attachment of nanoparticles at primary minima in the presence of surface roughness under unfavorable chemical conditions. Saturated sand-packed column experiments were conducted at 0.1 and 0.2M NaCl using 30 nm polystyrene latex nanoparticles as model colloids. Particle attachment at primary minima was unambiguously determined by removing particles attached at secondary minima through introducing deionized water and excavating the packed beds. The calculated primary-minimum attachment efficiency was found to decrease with increasing flow velocity, indicating that the fraction of a collector surface that is available for attachment at primary minima decreases with increasing flow velocity. The torque analysis, however, showed that the adhesive torque that the particle experiences at primary minima is much larger than the maximum hydrodynamic drags of a porous medium for the flow velocities used. We attributed the discrepancy to the reason that the sand surface is very rough and the roughness mainly causes the attachment in primary minima under the experimental conditions used in this study. By considering influence of surface roughness in the torque analysis, our calculations show that while particle attachment in primary minima is favored atop of nanoasperities under unfavorable conditions, the adhesive torque that the particle experiences can be greatly reduced and, thus, the attachment is susceptible to flow drag. Whereas the increase of adhesive torque by surface roughness has been widely recognized in the literature, our study indicates that the rough asperities can also decrease adhesive torques for particles attached atop of them.

Release of quantum dot nanoparticles in porous media: role of cation exchange and aging time

Understanding the fate and transport of engineered nanoparticles (ENPs) in subsurface environments is required for developing the best strategy for waste management and disposal of these materials. In this study, the deposition and release of quantum dot (QD) nanoparticles were studied in saturated sand columns. The QDs were first deposited in columns using 100 mM NaCl or 2 mM CaC12 solutions. Deposited QDs were then contacted with deionized (DI) water and/or varying Na(+) concentrations to induce release. QDs deposited in 100 mM Na(+) were easily reversible when the column was rinsed with DI water. Conversely, QDs deposited in the presence of Ca(2+) exhibited resistance to release with DI water. However, significant release occurred when the columns were flushed with NaCl solutions. This release behavior was explained by cation exchange (Ca(2+) in exchange sites were replaced by Na(+)) which resulted in the breakdown of calcium bridging. We also studied the effect of aging time on the QD release. As the aging time increased, smaller amounts of QDs were released following cation exchange. However, deposited QDs were subsequently released when the column was flushed with DI water. The release behavior was modeled using a single first-order kinetic release process and changes in the maximum solid phase concentration of deposited QDs with transition in solution chemistry. The results of this study demonstrate that the presence of carboxyl groups on ENPs and divalent ions in the solution plays a key role in controlling ENP mobility in the subsurface environment.

A theoretical analysis of colloid attachment and straining in chemically heterogeneous porous media

A balance of applied hydrodynamic (T(H)) and resisting adhesive (T(A)) torques was conducted over a chemically heterogeneous porous medium that contained random roughness of height h(r) to determine the fraction of the solid surface area that contributes to colloid immobilization (S(f)*) under unfavorable attachment conditions. This model considers resistance due to deformation and the horizontal component of the adhesive force (F(AT)), spatial variations in the pore scale velocity distribution, and the influence of hr on lever arms for T(H) and T(A). Values of S(f)* were calculated for a wide range of physicochemical properties to gain insight into mechanisms and factors influencing colloid immobilization. Colloid attachment processes were demonstrated to depend on solution ionic strength (IS), the colloid radius (r(c)), the Young's modulus (K), the amount of chemical heterogeneity (P+), and the Darcy velocity (q). Colloid immobilization was also demonstrated to occur on a rough surface in the absence of attachment. In this case, S(f)* depended on IS, r(c), the roughness fraction (f0), h(r), and q. Roughness tended to enhance T(A) and diminish T(H). Consequently, the effect of IS on S(f)* was enhanced by h(r) relative to attachment. In contrast, the effects of r(c) and q on S(f)* were diminished by hr in comparison to attachment. Colloid immobilization adjacent to macroscopic roughness locations shares many similarities to grain-grain contact points and may be viewed as a type of straining process. In general, attachment was more important for higher IS and variance in the secondary minimum, and for smaller r(c), q, and K, but diffusion decreased these values. Conversely, straining was dominant for the opposite conditions. Discrepancies in the literature on mechanisms of colloid retention are likely due to a lack of consideration of all of these factors.

Effects of fluctuations in river water level on virus removal by bank filtration and aquifer passage--a scenario analysis

Riverbank filtration is an effective process for removing pathogenic viruses from river water. Despite indications that changing hydraulic conditions during floods can affect the efficacy of riverbank filtration to remove viruses, the impact on advection and dispersion of viruses in the riverbank is not well understood. We investigated the effects of fluctuations in river water level on virus transport during riverbank filtration, considering 3-D transient groundwater flow and virus transport. Using constant removal rates from published field experiments with bacteriophages, removal of viruses with distance from the riverbank was simulated for coarse gravel, fine gravel and fine sandy gravel. Our simulations showed that, in comparison with steady flow conditions, fluctuations in river water level cause viruses to be transported further at higher concentrations into the riverbank. A 1-5 m increase in river water levels led to a 2- to 4-log (log10 reduction in concentration relative to the initial concentration in the river) increase in virus concentration and to up to 30% shorter travel times. For particular cases during the receding flood, changing groundwater flow conditions caused that pristine groundwater was carried from further inland and that simulated virus concentrations were more diluted in groundwater. Our study suggests that the adverse effect of water level fluctuations on virus transport should be considered in the simulation of safe setback distances for drinking water supplies.

Transport of microbial tracers in clean and organically contaminated silica sand in laboratory columns compared with their transport in the field

Waste disposal on land and the consequent transport of bacterial and viral pathogens in soils and aquifers are of major concern worldwide. Pathogen transport can be enhanced in the presence of organic matter due to occupation of attachment sites in the aquifer materials thus preventing pathogen attachment leading to their faster transport for longer distances. Laboratory column studies were carried out to investigate the effect of organic matter, in the form of dissolved organic carbon (DOC), on the transport of Escherichia coli and MS2 phage in saturated clean silica sand. Transport rates of these microbial tracers were also studied in a contaminated field site. Laboratory column studies showed that low concentrations (0.17 mg L(-1)) of DOC had little effect on E. coli J6-2 removal and slightly reduced the attachment of MS2 phage. After progressive conditioning of the column with DOC (1.7 mg L(-1) and 17 mg L(-1)), neither E. coli J6-2 nor MS2 phage showed any attachment and recovery rates increased dramatically (up to 100%). The results suggest that DOC can affect the transport rates of microbial contaminants. For E. coli J6-2 the predominant effect appeared to be an increase in the secondary energy minimum leading to an increase in E. coli attachment initially. However, after 17 mg L(-1) DOC conditioning of the silica sand no attachment of E. coli was observed as the DOC took up attachment sites in the porous media. MS2 phage appeared to be affected predominantly by out-competition of binding sites in the clean silica sand and a steady reduction in attachment was observed as the DOC conditioning increased. Field study showed a high removal of both E. coli and MS2 phage, although E. coli was removed at a lower rate than MS2 phage. In the field it is likely that a combination of effects are seen as the aquifer material will be heterogeneous in its surface nanoscale properties, demonstrated by the differing removal of E. coli and MS2 phage compared to the laboratory scale experiments. This research demonstrates the importance of combining laboratory scale and field scale studies to fully understand removal of microbes in groundwater aquifers affected by organic matter (DOC).

Application of DLVO energy map to evaluate interactions between spherical colloids and rough surfaces

This study theoretically evaluated interactions between spherical colloids and rough surfaces in three-dimensional space using Derjaguin-Landau-Verwey- Overbeek (DLVO) energy/force map and curve. The rough surfaces were modeled as a flat surface covered by hemispherical protrusions. A modified Derjaguin approach was employed to calculate the interaction energies and forces. Results show that more irreversible attachments in primary minima occur at higher ionic strengths, which theoretically explains the observed hysteresis of colloid attachment and detachment during transients in solution chemistry. Secondary minimum depths can be increased significantly in concave regions (e.g., areas aside of asperities or between asperities) due to sidewall interactions. Through comparing the tangential attractive forces from asperities and the hydrodynamic drag forces in three-dimensional space, we showed that attachment in secondary minima can be located on open collector surfaces of a porous medium. This result challenges the usual belief that the attachment in secondary minima only occurs in stagnation point regions of the porous medium and is absent in shear flow systems such as parallel plate flow chamber and impinging jet apparatus. Despite the argument about the role of secondary minima in colloid attachment remained, our study theoretically justified the existence of attachment in secondary minima in the presence of surface roughness. Further, our study implied that the presence of surface roughness is more favorable for attachment in secondary minima than in primary minima under unfavorable chemical conditions.

Colloid adhesive parameters for chemically heterogeneous porous media

A simple modeling approach was developed to calculate colloid adhesive parameters for chemically heterogeneous porous media. The area of the zone of electrostatic influence between a colloid and solid-water interface (A(z)) was discretized into a number of equally sized grid cells to capture chemical heterogeneity within this region. These cells were divided into fractions having specific zeta potentials (e.g., negative or positive values). Mean colloid adhesive parameters such as the zeta potential, the minimum and maximum in the interaction energy, the colloid sticking efficiency (α), and the fraction of the solid surface area that contributes to colloid immobilization (S(f)) were calculated for possible charge realizations within A(z). The probability of a given charge realization in A(z) was calculated using a binomial mass distribution. Probability density functions (PDFs) for the colloid adhesive parameters on the heterogeneous surface were subsequently calculated at the representative elementary area (REA) scale for a porous medium. This approach was applied separately to the solid-water interface (SWI) and the colloid, or jointly to both the SWI and colloid. To validate the developed model, the mean and standard deviation of the interaction energy distribution on a chemically heterogeneous SWI were calculated and demonstrated to be consistent with published Monte Carlo simulation output using the computationally intensive grid surface integration technique. Our model results show that the PDFs of colloid adhesive parameters at the REA scale were sensitive to the size of the colloid and the heterogeneity, the charge and number of grid cells, and the ionic strength.

Site-specific retention of colloids at rough rock surfaces

The spatial deposition of polystyrene latex colloids (d = 1 μm) at rough mineral and rock surfaces was investigated quantitatively as a function of Eu(III) concentration. Granodiorite samples from Grimsel test site (GTS), Switzerland, were used as collector surfaces for sorption experiments. At a scan area of 300 × 300 μm(2), the surface roughness (rms roughness, Rq) range was 100-2000 nm, including roughness contribution from asperities of several tens of nanometers in height to the sample topography. Although, an increase in both roughness and [Eu(III)] resulted in enhanced colloid deposition on granodiorite surfaces, surface roughness governs colloid deposition mainly at low Eu(III) concentrations (≤5 × 10(-7) M). Highest deposition efficiency on granodiorite has been found at walls of intergranular pores at surface sections with roughness Rq = 500-2000 nm. An about 2 orders of magnitude lower colloid deposition has been observed at granodiorite sections with low surface roughness (Rq < 500 nm), such as large and smooth feldspar or quartz crystal surface sections as well as intragranular pores. The site-specific deposition of colloids at intergranular pores is induced by small scale protrusions (mean height = 0.5 ± 0.3 μm). These protrusions diminish locally the overall DLVO interaction energy at the interface. The protrusions prevent further rolling over the surface by increasing the hydrodynamic drag required for detachment. Moreover, colloid sorption is favored at surface sections with high density of small protrusions (density (D) = 2.6 ± 0.55 μm(-1), asperity diameter (φ) = 0.6 ± 0.2 μm, height (h) = 0.4 ± 0.1 μm) in contrast to surface sections with larger asperities and lower asperity density (D = 1.2 ± 0.6 μm(-1), φ = 1.4 ± 0.4 μm, h = 0.6 ± 0.2 μm). The study elucidates the importance to include surface roughness parameters into predictive colloid-borne contaminant migration calculations.

Causes and implications of colloid and microorganism retention hysteresis

Experiments were designed to better understand the causes and implications of colloid and microorganism retention hysteresis with transients in solution ionic strength (IS). Saturated packed column experiments were conducted using two sizes of carboxyl modified latex (CML) microspheres (0.1 and 1.1 μm) and microorganisms (coliphage φX174 and E. coli D21g) under various transient solution chemistry conditions, and 360 μm Ottawa sand that was subject to different levels of cleaning, namely, a salt cleaning procedure that removed clay particles, and a salt+acid cleaning procedure that removed clay and reduced microscopic heterogeneities due to metal oxides and surface roughness. Comparison of results from the salt and salt+acid treated sand indicated that microscopic heterogeneity was a major contributor to colloid retention hysteresis. The influence of this heterogeneity increased with IS and decreasing colloid/microbe size on salt treated sand. These trends were not consistent with calculated mean interaction energies (the secondary minima), but could be explained by the size of the electrostatic zone of influence (ZOI) near microscopic heterogeneities. In particular, the depth of local minima in the interaction energy has been predicted to increase with a decrease in the ZOI when the colloid size and/or the Debye length decreased (IS increased). The adhesive interaction was therefore largely irreversible for smaller sized 0.1 μm CML colloids, whereas it was reversible for larger 1.1 μm CML colloids. Similarly, the larger E. coli D21g exhibited greater reversibility in retention than φX174. However, direct comparison of CML colloids and microbes was not possible due to differences in size, shape, and surface properties. Retention and release behavior of CML colloids on salt+acid treated sand was much more consistent with mean interaction energies due to reduction in microscopic heterogeneities.