Release of colloids from primary minimum contact under unfavorable conditions by perturbations in ionic strength and flow rate

Effect of particle density on microplastics transport in artificial and natural porous media

The occurrence and persistence of microplastics (MPs) in natural environments are of increasing concern. Along with this, the transport of MPs in sediments has been investigated mainly focusing on the effect of plastic size and shape, media size effect, and solution chemistry. Yet, the influence of particle density is only partially understood. Therefore, column experiments on the transport of variably buoyant MPs in saturated natural sediments and glass beads were conducted, and transport parameters were quantified using a two-site kinetic transport model with a depth-dependent blocking function (the amount of retained MPs does not decrease at a constant rate with increasing depth, the majority of MPs were retained near the column inlet). Neutral, sinking, and buoyant MPs within the same size range were selected, with stable water isotope applied as conservative tracer to explore water and MP movement in the tested sediments. The results showed that 95.5 ± 1.4% of sinking MPs remained in columns packed with gravel, followed by buoyant and neutral MPs, thus indicating that particle density does affect MP mobility. Similar recovered amounts of MPs were found in columns packed with glass beads, indicating that tested sediment types do not affect the deposition behavior of MPs. The breakthrough curves of MPs were accurately described by the selected model. However, the simulated retention profiles overestimated the observed MP amount in layers closest to the column inlet. The coupled experimental and modeled results suggest an enhanced retention of sinking MPs, while neutrally and buoyant MPs exhibit a higher mobility in comparison. Thus, neutral or buoyant MPs can potentially pose a higher contamination risk to subsurface porous media environments compared to sinking MPs. Discrepancies between observed and simulated retention profiles indicate that future model development is needed for advancing the MP deposition as affected by particle density.

Polymer Coatings Affect Transport and Remobilization of Colloidal Activated Carbon in Saturated Sand Columns: Implications for In Situ Groundwater Remediation

Colloidal activated carbon (CAC) is an emerging technology for the in situ remediation of groundwater impacted by per- and polyfluoroalkyl substances (PFAS). In assessing the long-term effectiveness of a CAC barrier, it is crucial to evaluate the potential of emplaced CAC particles to be remobilized and migrate away from the sorptive barrier. We examine the effect of two polymer stabilizers, carboxymethyl cellulose (CMC) and polydiallyldimethylammonium chloride (PolyDM), on CAC deposition and remobilization in saturated sand columns. CMC-modified CAC showed high mobility in a wide ionic strength (IS) range from 0.1 to 100 mM, which is favorable for CAC delivery at a sufficient scale. Interestingly, the mobility of PolyDM-modified CAC was high at low IS (0.1 mM) but greatly reduced at high IS (100 mM). Notably, significant remobilization (release) of deposited CMC-CAC particles occurred upon the introduction of solution with low IS following deposition at high IS. In contrast, PolyDM-CAC did not undergo any remobilization following deposition due to its favorable interactions with the quartz sand. We further elucidated the CAC deposition and remobilization behaviors by analyzing colloid-collector interactions through the application of Derjaguin-Landau-Verwey-Overbeek theory, and the inclusion of a discrete representation of charge heterogeneity on the quartz sand surface. The classical colloid filtration theory was also employed to estimate the travel distance of CAC in saturated columns. Our results underscore the roles of polymer coatings and solution chemistry in CAC transport, providing valuable guidelines for the design of in situ CAC remediation with maximized delivery efficiency and barrier longevity.

3D pore-scale characterization of colloid aggregation and retention by confocal microscopy: Effects of fluid structure and ionic strength

Understanding the mechanisms of colloid transport and retention as well as the spatial distribution of colloids in porous media is an important topic for contamination transport and remediation in subsurface environments. Utilizing advanced three-dimensional visualization experiments, we effectively capture the intricate distribution characteristics of colloids in the 3D pore space and quantify the size of colloid clusters that aggregate at fluid-fluid interfaces and solid surfaces during two-phase flow. Our experimental results reveal the influence of pore-scale events, such as Haines jumps and pinch-off, on colloid retention. Our results also indicate that large drainage rates can facilitate colloid retention on solid surfaces, especially under the condition of high ionic strength. This can be attributed to the migration of colloids from the fluid-fluid interface to the solid surface, propelled by transients in the local fluid structure. The findings reveal a synergistic effect of the ionic strength and hydrodynamic conditions on colloid transport and retention during two-phase flow and provide important insights for predicting the fate and transport of contaminants in soil and groundwater environments involving multiple fluid phases.

An extended colloid filtration theory for modeling Escherichia coli transport in 3-D fracture networks

Microbial transport in fractured carbonate rock using enhanced solutions is a significant and neglected research topic in the literature. We propose an extended colloid filtration theory (CFT) combined with a particle-tracking following streamlines (PTFS) model for the rapid prediction of breakthrough curves (BTCs) and plumes of pathogens in three-dimensional (3-D) discrete fracture networks (DFNs). We adapted CFT in porous media to pathogen transport in fractures containing Terra Rossa (soil) deposits. As an example of the model capability, a simulation was used to predict the 3-D motion field and Escherichia coli count in groundwater originating from the Forcatella managed aquifer recharge (MAR) Facility (Brindisi, Italy) using a DFN composed of 3,900 fractures. In arid regions, MAR facilities are significant for sustaining basic human needs, such as freshwater supply for drinking and crop production. The Markov chain Monte Carlo (MCMC) technique was applied to E. coli counts in the collected water samples to increase data representativeness. The pathogen transport coefficients were further supported by batch filtration tests carried out in the CNR/IRSA Laboratory (Bari, Italy). The mean E. coli attachment rate coefficient of 0.15 × 10-8 m2 d-1 (sticking efficiency = 1.1 × 10-8 m) resulted in a 2.1 log10 removal in 600 m of reclaimed water filtration. The simulation output visualized the E. coli 3-D pathways in groundwater and the positions of contaminated groundwater spring outflows on Forcatella Beach. The simulation results agreed with the mean MCMC output of E. coli concentrations in bathing water under unperturbed geochemical and environmental flow and transport conditions. However, results indicate that concentrations of pathogenic strains, parasites, and enteric viruses may enter the marine environment of MAR sites during flood periods.

QMRA of Ralstonia solanacearum in potato cultivation: Risks associated with irrigation water recycled through managed aquifer recharge

Agricultural aquifer storage recovery and transfer (ASTR) stores excess fresh water for later reuse in irrigation. Moreover, water quality improves because chemical pollutants and pathogens will be removed by degradation and attachment to the aquifer material. The source water may contain the bacterial plant pathogen Ralstonia solanacearum which causes plant infections and high yield losses. We used quantitative microbial risk assessment (QMRA) to investigate the removal of R. solanacearum during ASTR to predict infection risks of potato plants after irrigation with the recovered water. Laboratory experiments analyzed the ASTR treatment by investigating the bacterial die-off in the water phase and the removal by attachment to the aquifer sediment. Die-off in the water phase depends on the residence time and ranged between 1.3 and 2.7 log10 after 10 or 60 days water storage, respectively. A subpopulation of the bacteria persisted for a prolonged time at low concentrations which may pose a risk if the water is recovered too early. However, the natural aquifer sand filtration proofed to be highly effective in removing R. solanacearum by attachment which depends on the distance between injection and abstraction well. The high removal by attachment alone (18 log10 after 1 m) would reduce bacterial concentrations to negligible numbers. Upscaling to longer soil passages is discussed in the paper. Infection risks of potato plants were calculated using a dose-response model and ASTR treatment resulted in negligible infection risks of a single plant, but also when simulating the irrigation of a 5 ha potato field. This is the first QMRA that analyzed an agricultural ASTR and the fate of a plant pathogen focusing on plant health. QMRA is a useful (water) management tool to evaluate the treatment steps of water reclamation technologies with the aim to provide safe irrigation water and reduce risks disseminating plant diseases.

Retention and release of black phosphorus nanoparticles in porous media under various physicochemical conditions

Black phosphorus nanosheets/nanoparticles (BPNs) are widely applied in many fields. However, the transport of BPNs in the subsurface still has not yet been reported and there is increasing concern about potential adverse impacts on ecosystems. Roles of median grain size and surface roughness, BPN concentration, and solution chemistries (pH, ionic strength, and cation types) on the retention and release of BPNs in column experiments were therefore investigated. The mobility of BPNs significantly increased with increasing grain size and decreasing surface roughness due to their influence on the mass transfer rate, number of deposition sites and retention capacity, and straining processes. Transport of BPNs was enhanced with an increase in pH and a decrease in ionic strength because of surface deprotonation and stronger repulsion that tends to reduce aggregation. The BPN transport was significantly sensitive to ionic strength, compared with other engineered nanoparticles. Additionally, charge heterogeneity and cation-bridging played a critical role in the retention of BPNs in the presence of divalent cations. Higher input concentrations increased the retention of BPNs, probably because collisions, aggregation at pore throat locations, and hydrodynamic bridging were more pronounced. Small fractions of BPNs can be released under decreasing IS and increasing pH due to the expansion of the electrical double layer and increased repulsion at convex roughness locations. A mathematical model that includes provisions for advective dispersive transport and time-dependent retention with blocking or ripening terms well described the retention and release of BPNs. These findings provide fundamental information that helps to understand the transport of BPNs in the subsurface environments.

Relating mechanistic fate with spatial positioning for colloid transport in surface heterogeneous porous media

HYPOTHESES

The transport behavior of colloids in subsurface porous media is altered by surface chemical and physical heterogeneities. Understanding the mechanisms involved and distribution outcomes is crucial to assess and control groundwater contamination. The multi-scale processes that broaden residence time distribution for particles in the medium are here succinctly described with an upscaling model. Experiments/model: The spatial distribution of silver particles along glass bead-packed columns obtained from X-ray micro-computed tomography and a mechanistic upscaling model were used to study colloid retention across interface-, collector-, pore-, and Darcy-scales. Simulated energy profiles considering variable colloid-grain interactions were used to determine collector efficiencies from particle trajectories via full force-torque balance. Rate coefficients were determined from collector efficiencies to parameterize the advective-dispersive-reactive model that reports breakthrough curves and depth profiles.

FINDINGS

Our results indicate that: (i) with surface heterogeneity, individual colloid-grain interactions are non-unique and span from repulsive to attractive extremes; (ii) experimentally observed spatial positioning of retention at grain-water interfaces and grain-to-grain contacts is governed respectively by mechanistic attachment to the grain surface and retention without contact at rear-flow stagnation zones, and (iii) experimentally observed non-monotonic retention profiles and heavy-tailed breakthrough curves can be modeled with explicit implementation of heterogeneity at smaller scales.

Influence of biochar incorporation on the collector surface properties and the transport of silver nanoparticles in porous media

Biochar is widely used as a soil amendment due to its environmental friendliness and convenient availability. It is believed that the presence of biochar in porous media can influence the transport of colloidal and solute contaminants. In this study, different mass ratios of biochar were added to packed sand with a rough or smooth surface to determine the significance of biochar on the retention and release of silver nanoparticles (AgNPs). The results showed that biochar reduced the transport of AgNPs in rough and smooth sands under different solution conditions. A small amount of biochar (0.1-1% in mass percentage) can significantly enhance the retention of AgNPs due to the alteration in collector surface roughness and chemical heterogeneity that potentially reduce the energy barrier for retention. Furthermore, the retention of AgNPs in rough sand was always higher than that in smooth sand under the same experimental conditions. The presence of biochar also produced nonmonotonic retention of AgNPs mainly due to the changes in collector surface roughness. Additionally, the AgNPs retention associated with biochar tended to be irreversible due to the charge heterogeneity, while the reversible retention could mainly occur on a rough sand surface via shallow primary minima. This work highlights the significance of collector surface roughness that needs to be considered in the process of biochar amendment for practical applications to effectively immobilize colloidal contaminants in soil or groundwater.

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.

Effects of ionic strength and particle size on transport of microplastic and humic acid in porous media

As an emerging pollutant, the transport behavior of colloidal microplastic particles (CMPs) in saturated porous media may be affected by the simultaneous presence of other substances in the natural environment. In this study, colloidal polystyrene microplastic particles (PSMPs) were selected as the representative of CMPs to investigate the cotransport behaviors of CMPs in the presence of humic acid (HA) under varied environmental conditions (ionic strength: 1, 100 mM KCl; HA concentration: 0, 5, 10, 20 mg⋅L^-1) in porous media. The presence of HA with different concentrations was found to increase the mobility of 1.0-μm and 0.2-μm CMPs in porous media in a non-linear and non-monotonic manner. Furthermore, the HA-facilitated transport of CMPs occurred under both electrostatically unfavorable and favorable attachment conditions (limited to the conditions examined in this study, corresponding to 1 and 100 mM KCl, respectively). The transport behavior of the smaller-sized CMPs (0.2-μm CMPs) was more sensitive to the change of ionic strength and the presence of HA than that of the larger-sized CMPs (1.0-μm CMPs). The cotransport process of CMPs and HA was affected by many factors. Modeling results showed that a small amount of competitive blocking occurred during the cotransport process. Moreover, both the presence of HA and change in ionic strength could affect the surface properties of CMPs. Thus, the cotransport behavior of CMPs with HA was different from the transport of individual CMPs in porous media. Experimental results revealed that HA induced complexity in the transport behavior of CMPs in the aqueous environment. Therefore, undeniably, a lot more systematic explorations are further demanded to better comprehend the CMPs cotransport mechanism in the presence of other substances.

Review of Particle Detachment and Attachment in Porous Media

Particle detachment, migration and attachment are common processes in porous media, especially in unconsolidated formations. In this review, the processes are discussed and equations describing the processes are presented. Two particle detachment processes analyzed are the hydrodynamic forces and electric double-layer forces. The particle detachments equations were critically examined to determine if they reflect crucial factors that trigger particle detachment in porous media. Essential factors that are missing in the equation are the effect of pressure and the level of rock consolidation. Incorporating the level of rock cementation and the effect of pressure in the equations will make the models more empirical and less theoretical. For particle attachment, Van der Waals forces, adhesion, particle attachment efficiency, and straining processes and their equations are considered. The colloidal forces are all embraced in terms of capturing important elements that mobilize particles in porous media, however, the practical application of the models can pose a challenge. For particle adsorption on grain surfaces, it is recommended that the effect of pressure and temperature be studied.

Micro- and nanoplastics retention in porous media exhibits different dependence on grain surface roughness and clay coating with particle size

The presence and/or coating of natural colloids (e.g., clays and metal oxides or hydroxides) on collector surfaces has frequently been demonstrated to enhance the retention of engineered colloids that are negatively charged due to favorable electrostatic interactions. However, this work demonstrates that the presence of natural clay coating can lead to reduced or nonmonotonic retention of micro- and nanoplastics (MNPs). Column experiments were carried out to systematically investigate the transport of MNPs with different sizes in relatively smooth and rough sands that had various clay coating fractions. These coating fractions on the collector were found to significantly influence MNP retention in a complex manner that changed with the colloid size and the roughness properties of the sand. This observation was attributed to the impact of clay coatings on the roughness and morphology properties of collector surfaces that were dominant over surface charge. Scanning electron microscopy and interaction energy calculations on surfaces with pillars or valleys indicate that mechanisms that contributed to MNP retention changed with the colloid size. In particular, retention of nanosized plastics was mainly controlled by interactions on convex/concave locations that changed with the solution chemistry, whereas microsized plastics were also strongly influenced by the applied hydrodynamic torque and straining processes. Additionally, the significant sensitivity of MNP retention under a low-level ionic strength also reflects the importance of roughness and charge heterogeneities. These observations are important for investigating the mechanisms of colloid transport in natural systems that ubiquitously exhibit clay coating on their surfaces.

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.

Important Role of Concave Surfaces in Deposition of Colloids under Favorable Conditions as Revealed by Microscale Visualization

This study conducted saturated column experiments to systematically investigate deposition of 1 μm positively charged polystyrene latex micro-colloids (representing microplastic particles) on negatively charged rough sand, glass beads, and soil with pore water velocities (PWV) from 4.9 × 10^-5 to 8.8 × 10^-4 m/s. A critical value of PWV was found below which colloidal attachment efficiency (AE) increased with increasing PWV. The increase in AE with PWV was attributed to enhanced delivery of the colloids and subsequent attachment at concave locations of rough collector surfaces. The AE decreased with further increasing PWV beyond the threshold because the convex sites became unavailable for colloid attachment. By simulating the rough surfaces using the Weierstrass-Mandelbrot equation, the extended Derjaguin-Landau-Verwey-Overbeek (XDLVO) interaction energy calculations and torque analysis revealed that the adhesive torques could be reduced to be comparable or smaller than hydrodynamic torques even under the favorable conditions. Interestingly, scanning electron microscopic experiments showed that blocking occurred at convex sites at all ionic strengths (ISs) (e.g., even when the colloid-colloid interaction was attractive), whereas at concave sites, blocking and ripening (i.e., attached colloids favor subsequent attachment) occurred at low and high ISs, respectively. To our knowledge, our work was the first to show coexistence of blocking and ripening at high ISs due to variation of the collector surface morphology.

Favorable and unfavorable attachment of colloids in a discrete sandstone fracture

The transport of cationic amine-modified latex (AML) and anionic carboxylate-modified latex (CML) microspheres through a discrete sandstone fracture with mineralogical heterogeneity and roughness was studied. Two microsphere sizes (200 nm and 1000 nm), two ionic strengths (5 mM and 10 mM), and two specific discharges (0.35 mm.s-1 and 0.70 mm.s-1) were tested to observe the impact on transport under favorable and unfavorable conditions. The difference in retention between AML (net favorable) and CML (net unfavorable) microsphere attachment was 25% for the 200 nm microspheres and 13% for the 1000 nm microspheres. Less than 50% of the AML microspheres were retained in the fracture, postulated to be due to the effects of mineralogical heterogeneity and fracture surface roughness. The effect of an increase in ionic strength in increasing retention was significant for unfavorable attachment, but was not significant for favorable attachment conditions. The effect of specific discharge was minor for all but the 200 nm CML microspheres at 10 mM ionic strength. When flushing the fracture first with cationic microspheres, then with anionic microspheres, the recovery of anionic microspheres resembled favorable attachment presumably due to interaction with cationic microspheres that remained attached to the sandstone surface. Colloid breakthrough curves could be fit well with a two site attachment model, with reversible and irreversible sites.

Managed aquifer recharge implementation criteria to achieve water sustainability

Depletion of groundwater is accelerated due to an increase in water demand for applications in urbanized areas, agriculture sectors, and energy extraction, and dwindling surface water during changing climate. Managed aquifer recharge (MAR) is one of the several methods that can help achieve long-term water sustainability by increasing the natural recharge of groundwater reservoirs with water from non-traditional supplies such as excess surface water, stormwater, and treated wastewater. Despite the multiple benefits of MAR, the wide-scale implementation of MAR is lacking, partly because of challenges to select the location for MAR implementation and identify the MAR type based on site conditions and needs. In this review, we provide an overview of MAR types with a basic framework to select and implement specific MAR at a site based on water availability and quality, land use, source type, soil, and aquifer properties. Our analysis of 1127 MAR projects shows that MAR has been predominantly implemented in sites with sandy clay loam soil (soil group C) and with access to river water for recharge. Spatial analysis reveals that many regions with depleting water storage have opportunities to implement MAR projects. Analyzing data from 34 studies where stormwater was used for recharge, we show that MAR can remove dissolved organic carbon, most metals, E. coli but not efficient at removing most trace organics, and enterococci. Removal efficiency depends on the type of MAR. In the end, we highlight potential challenges for implementing MAR at a site and additional benefits such as minimizing land subsidence, flood risk, augmenting low dry-season flow, and minimizing salt-water intrusion. These results could help identify locations in the water-stressed regions to implement specific MAR for water sustainability.

A modeling approach for quantitative assessment of interfacial interaction between two rough particles in colloidal systems

HYPOTHESIS AND BACKGROUND

The simulation of rough particle surface is important to understand and control the interface behavior of particles in colloidal systems. Literature analysis suggested a lack of information for an accurate model simulating the interfacial interaction between two rough particles. It is hypothesized that the total interfacial energy developed between two rough particles would depend on the surface morphologies of particles, and it could be predicted if a mathematical model to represent the interaction of two rough particles were created accurately.

EXPERIMENTS

In this study, mathematical models were developed to determine the interfacial energy created between two particles according to the XDLVO theory by considering the rippled particle theory and surface element integral (SEI) method. Three different scenarios of particle interactions were assumed in the simulation. The present study provides deep insights into particle interactions via considering aspect ratio, size, and surface roughness of two particles in colloidal systems.

FINDINGS

The assessment of the interfacial interaction revealed that an increase in the aspect ratio, surface roughness, and relative surface roughness of particles would weaken the total interaction energy generated between particles and promote particle aggregation. Increased interaction energy was predicted for the interaction of particles by increasing the particle size. The asperity ratio was more effective than the asperity number in controlling the interfacial energy between two particles. The results of this study could be used for foreseeing the interaction of rough particles, which has a significant application in particle coagulation or dispersion in colloidal systems.

Influence of agricultural organic inputs and their aging on the transport of ferrihydrite nanoparticles: From enhancement to inhibition

Organic matter effectively regulates nanoparticles transport. However, little is known about the effect of agricultural organic inputs on the transport of ferrihydrite nanoparticles (FHNPs) during aging. In this study, columns were filled with sand mixed with varying proportions of pristine, water-processing, or alkali-processing biochar or swine manure and used to simulate the release of organic matter and changes in surface roughness of sand grains during field aging. The influence of these factors on FHNPs transport was investigated using column experiments. The dissolved organic matter (DOM) (0.008-24.8 mg L^-1) released from agricultural organic inputs decreased the zeta potential of the FHNPs from 30.8 mV to 14.6--48.9 mV and further caused electrostatic repulsion, osmotic repulsion, and elastic-steric repulsion between FHNPs and mixed sand, thus enhancing FHNPs transport. Ferrihydrite nanoparticles transport increased with increasing content of biochar and swine manure due to the increased amount of DOM. However, with the presence of organic inputs, surface roughness up to a certain degree (the increase in specific surface area up to 4.6 m²) became the dominant inhibition factor affecting FHNPs transport. After DOM release, agricultural organic inputs decreased the enhancement of FHNPs transport; with the increase input, their rougher surface gradually increased inhibition of FHNPs transport. The strongest FHNPs retention in the alkali-processing biochar (0.2-10%) or swine manure (1-2%) mixed sand columns indicated that fully aged agricultural organic inputs strongly inhibited FHNPs transport. Our findings provided novel insights into the critical influence of agricultural organic inputs and their aging on FHNPs transport, which changed gradually from enhancement to inhibition gradually.

Humic acid induced weak attachment of fullerene nC60 nanoparticles and subsequent detachment upon reduction of solution ionic strength in saturated porous media

Sand column experiments were performed under saturated conditions to investigate impact of humic acid (HA) on attachment of nC₆₀ nanoparticles (NPs) in NaCl and CaCl₂ at ionic strengths (ISs) from 1 mM to 100 mM and subsequent detachment via reducing solution IS. The attachment increased with increasing IS due to reduced repulsive Derjaguin-Landau-Verwey-Overbeek (DLVO) interaction energy and accordingly increased retention in primary energy wells. More attachments occurred in CaCl₂ compared to NaCl because Ca²⁺ exhibited greater charge screen ability and served as a bridging agent between the NPs and sand surfaces. The presence of HA significantly reduced nC₆₀ NPs attachment on sand surfaces (especially on nanoscale physical heterogeneities) in 10 mM NaCl and 1 mM CaCl₂ because of enhanced electrostatic and steric repulsions. Interestingly, although the HA did not cause reduction of attachment in 100 mM NaCl and 10 mM CaCl₂ compared to the case in absence of HA, the HA caused weak attachment of nC₆₀ on sand surfaces and then much more significant detachment by decreasing IS. The HA did not alter both attachment and detachment in 100 mM CaCl₂, because the Ca²⁺ at the high concentration caused formation of very stable complex of HA and NPs, and strong interaction of the complex with the sand surfaces via cation bridge. Our study highlighted that the HA can not only enhance the transport of NPs by inhibiting attachment as revealed in the literature, but also by the continuous capture and release of the NPs from surfaces in subsurface environments.

Release of colloidal biochar during transient chemical conditions: The humic acid effect

Our understanding of colloidal biochar (CB) transport and release is largely unknown in environments with transient chemical conditions, e.g., ionic strength (IS), pH, and especially humic acid (HA). In this study, column experiments were conducted to investigate CB transport and retention in the presence and absence of HA, and CB release under transient IS and pH conditions in saturated sand. Step reductions in solution IS from 25 to 0.01 mM produced significant release peaks of CB due to a reduction in the depth of the primary minima on rough surfaces with small energy barriers. In contrast, step increases of solution pH from 4 to 10 only slightly increased CB release presumably due to the strong buffering capacity of CB. The CB retention was diminished by HA during the deposition phase. However, the release of CB with transients in IS and pH was not influenced much when deposition occurred in the presence of HA. These observations indicate that HA increased the energy barrier during deposition but did not have a large influence on the depth of the interacting minimum during transient release. Potential explanations for these effects of HA on CB retention and transient release include enhanced repulsive electrostatic interactions and/or altering of surface roughness properties. Our findings indicated that the release of retained CB is sensitive to transient IS conditions, but less dependent on pH increases and CB deposition in the presence of HA. This information is needed to quantify potential benefits and/or adverse risks of mobile CB in natural environments.

A Unified Force and Torque Balance for Colloid Transport: Predicting Attachment and Mobilization under Favorable and Unfavorable Conditions

Colloid attachment and detachment behaviors concern a wide range of environmental contexts but have typically been mechanistically predicted exclusive of one another despite their obvious coupling. Furthermore, previous mechanistic prediction often addressed packed column contexts, wherein specific forces and torques on the colloid could not be well-constrained, preventing robust predictions. These weaknesses were addressed through direct observation experiments under conditions where perfect sink assumptions fail and allow calibration of the contact between the colloid and collector. Attachment and flow perturbation experiments in the presence of colloid-collector attraction (favorable conditions) permitted calibration of contact parameters without the complexity that comes with colloid-collector repulsion (unfavorable conditions). Combining calibrated contact parameters with discrete representative nanoscale heterogeneity, developed to predict unfavorable attachment, provided an independent means to predict unfavorable detachment. The result was mechanistic prediction of colloid attachment and detachment that quantitatively agreed with experimental observation for both ionic strength and flow perturbation results, improving significantly upon previous qualitative prediction.

An experimental study of cotransport of heavy metals with kaolinite colloids

Cotransport of heavy metals, Pb, Cu and Zn (multi-metal system), and transport of those metals (single-metal system) were investigated by performing laboratory soil column experiment under the presence of kaolinite colloids. Preequilibrated kaolinite colloids with heavy metal solution was injected to the column until 10 pore volumes under two different flow rates and three different concentration of kaolinite colloids. Heavy metal concentration in effluent showed that the mobility of Pb was facilitated as kaolinite colloids concentration (C_c0) increases under high flow rate while the mobility of Pb and Cu were retarded as C_c0 increases under low flow rate. In addition, optimized first order rate coefficient related to sand-heavy metal interaction and estimated bed efficiency of experimental breakthrough curves demonstrated that the presence of mobile kaolinite colloids delayed the adsorption of heavy metals to the sand and facilitated the transport. Colloid associated contaminant transport model used in this study was found to be well fitted to the experimental breakthrough curves with the parameters associated with observed heavy metal transport without kaolinite colloids and adsorption/desorption between the heavy metals and the mobile kaolinite colloids.

Interactions between nanoparticles and fractal surfaces

This study evaluated attachment of a 30-nm nanoparticle to and detachment from fractal surfaces by calculating Derjaguin-Landau-Verwey-Overbeek (DLVO) interaction energies in three-dimensional space using the surface element integration technique. The fractal surfaces were generated using the Weierstass-Mandelbrot function with varying values of fractal dimension D (2.3 ≤ D ≤ 2.7) and fractal roughness G (0.000136 ≤ G ≤ 0.136). Results show that maximum energy barrier is reduced at peak areas of a fractal surface, and hence attachment in primary minima is favored. Some nanoparticles attached in primary minima at the peak areas can be detached by decreasing ionic strength (IS) due to monotonic decrease of interaction energy with increasing separation distance at low ISs. While the attachment in primary minima at valley areas is irreversible to IS reduction, the attachment is inhibited due to enhanced maximum energy barrier at these areas. A nonmonotonic variation of attachment efficiency in primary minimum (AEPM) with IS is present at high fractal dimension (D ≥ 2.4) or low fractal roughness (G < 0.00136), whereas the AEPM decreases monotonically with decreasing IS at low fractal dimension (D < 2.4) or high fractal roughness (G ≥ 0.00136). The AEPM decreases monotonically with increasing D or decreasing G at ISs from 1 mM to 200 mM. The decrease of AEPM with D or G is much slower at 10 mM compared to other ISs. These theoretical findings can explain various experimental observations in the literature, and can have important utility to development of water filtration techniques in engineered systems and to assessment of environmental risks of nanoparticles.

Hydrodynamic versus Surface Interaction Impacts of Roughness in Closing the Gap between Favorable and Unfavorable Colloid Transport Conditions

Recent experiments revealed that roughness decreases the gap in colloid attachment between favorable (repulsion absent) and unfavorable (repulsion present) conditions through a combination of hydrodynamic slip and surface interactions with asperities. Hydrodynamic slip was calibrated to experimentally observed tangential colloid velocities, demonstrating that slip length was equal to maximum asperity relief, thereby providing a functional relationship between slip and roughness metrics. Incorporation of the slip length in mechanistic particle trajectory simulations yielded the observed modest decrease in attachment over rough surfaces under favorable conditions, with the observed decreased attachment being due to reduced colloid delivery rather than decreased attraction. Cumulative interactions with multiple asperities acting within the zone of colloid-surface interaction were unable to produce the observed dramatic increased attachment and decreased reversibility with increased roughness under unfavorable conditions, necessitating inclusion of nanoscale attractive heterogeneity that was inferred to have codeveloped with roughness. Simulated attachment matched experimental observations when the spatial frequency of larger heterodomains (nanoscale zones of attraction) increased disproportionately relative to smaller heterodomains as roughness increased, whereas attachment was insensitive to asperity properties, including the number of interactions per asperity and asperity height; colloid detachment simulations were highly sensitive to these parameters. These cumulative findings reveal that hydrodynamic slip moderately decreases colloid bulk delivery, nanoscale heterogeneity dramatically enhances colloid attachment, and multiple interactions among asperities decrease detachment from rough surfaces.

Deposition of engineered nanoparticles (ENPs) on surfaces in aquatic systems: a review of interaction forces, experimental approaches, and influencing factors

The growing development of nanotechnology has promoted the wide application of engineered nanomaterials, raising immense concern over the toxicological impacts of nanoparticles on the ecological environment during their transport processes. Nanoparticles in aquatic systems may undergo deposition onto environmental surfaces, which affects the corresponding interactions of engineered nanoparticles (ENPs) with other contaminants and their environmental fate to a certain extent. In this review, the most common ENPs, i.e., carbonaceous, metallic, and nonmetallic nanoparticles, and their potential ecotoxicological impacts on the environment are summarized. Colloidal interactions, including Derjaguin-Landau-Verwey-Overbeek (DLVO) and non-DLVO forces, involved in governing the depositional behavior of these nanoparticles in aquatic systems are outlined in this work. Moreover, laboratory approaches for examining the deposition of ENPs on collector surfaces, such as the packed-bed column and quartz crystal microbalance (QCM) method, and the limitations of their applications are outlined. In addition, the deposition kinetics of nanoparticles on different types of surfaces are critically discussed as well, with emphasis on other influencing factors, including particle-specific properties, particle aggregation, ionic strength, pH, and natural organic matter. Finally, the future outlook and challenges of estimating the environmental transport of ENPs are presented. This review will be helpful for better understanding the effects and transport fate of ENPs in aquatic systems. Graphical abstract ᅟ.

Why Variant Colloid Transport Behaviors Emerge among Identical Individuals in Porous Media When Colloid-Surface Repulsion Exists

We herein demonstrate the cause of well-observed variant transport behaviors for apparently identical colloids in porous media under conditions of colloid-collector repulsion (unfavorable attachment conditions). We demonstrate that variant colloid transport behavior under unfavorable conditions can be explained by inherently variable colloid residence times prior to arrest on grains (collectors). We demonstrate that the residence time distributions derived from particle trajectory simulations incorporating representative nanoscale heterogeneity provide quantitative prediction of colloid transport under unfavorable conditions. We quantitatively predict hyper-exponential retention profiles in glass beads from representative nanoscale heterogeneity determined for glass, and we qualitatively predict nonmonotonic retention profiles in quartz sand from an estimated representative nanoscale heterogeneity for quartz. We also demonstrate that the transition from hyper-exponential to nonmonotonic profiles among glass beads versus quartz sand under otherwise equivalent conditions is primarily driven by greater grain angularity and consequent greater length and number of grain to grain contacts in quartz sand relative to glass beads. That continuum-scale transport behaviors emerge from upscaling of simulated pore-scale colloid residence times corroborates the utility of representative nanoscale heterogeneity for quantitative prediction of colloid transport under unfavorable conditions.

Enhanced transport of ferrihydrite colloid by chain-shaped humic acid colloid in saturated porous media

Both humic acid and colloid particle size effectively regulate colloid transport. However, little is known about effect of particle size and configuration of humic acid colloid (HA_colloid) on enhanced-transport of ferrihydrite colloid (FH_colloid) in porous media. Co-transport of HA_colloid and FH_colloid at different pH was systematically investigated by monitoring breakthrough curves (BTCs) in saturated sand columns. The colloid transport model and the (X)DLVO theory were used to reveal the mechanism of HA_colloid-enhanced FH_colloid transport in the columns. Results showed that HA_colloid enhanced FH_colloid transport in neutral and alkaline conditions. In neutral conditions, small HA_colloid (F-HA_colloid) with chain-shaped structure enhanced FH_colloid transport more prominently than pristine granular HA_colloid. The chain-shaped F-HA_colloid caused osmotic repulsion and elastic-steric repulsion between colloids and sand, leading to enhanced transport. However, the granular HA_colloid readily occurred as deposition due to attachment and straining, which decreased the enhanced transport of FH_colloid. In alkaline conditions, both HA_colloid and F-HA_colloid were chain-shaped, with longer chains of HA_colloid than F-HA_colloid. Ferrihydrite colloid transport was enhanced by HA_colloid more significantly than F-HA_colloid due to stronger repulsion between mixed HA_colloid-FH_colloid and sand. It suggested that regulation of particle size and morphology of HA_colloid would enhance FH_colloid transport and further help in understanding FH_colloid-facilitated contaminants transport in porous media.

A binomial modeling approach for upscaling colloid transport under unfavorable attachment conditions: Emergent prediction of non-monotonic retention profiles

We used a recently developed simple mathematical network model to upscale pore-scale colloid transport information determined under unfavorable attachment conditions. Classical log-linear and non-monotonic retention profiles, both well-reported under favorable and unfavorable attachment conditions, respectively, emerged from our upscaling. The primary attribute of the network is colloid transfer between bulk pore fluid, the near surface fluid domain (NSFD), and attachment (treated as irreversible). The network model accounts for colloid transfer to the NSFD of down-gradient grains and for reentrainment to bulk pore fluid via diffusion or via expulsion at rear flow stagnation zones (RFSZs). The model describes colloid transport by a sequence of random trials in a 1D network of Happel cells, which contain a grain and a pore. Using combinatorial analysis that capitalizes on the binomial coefficient, we derived from the pore-scale information the theoretical residence time distribution of colloids in the network. The transition from log-linear to non-monotonic retention profiles occurs when the conditions underlying classical filtration theory are not fulfilled, i.e., when a NSFD colloid population is maintained. Then, nonmonotonic retention profiles result, potentially both for attached and NSFD colloids. The concentration maxima shift downgradient depending on specific parameter choice. The concentration maxima were also shown to shift downgradient temporally (with continued elution) under conditions where attachment is negligible, explaining experimentally-observed down-gradient transport of retained concentration maxima of adhesion-deficient bacteria. For the case of zero reentrainment, we develop closed form, analytical expressions for the shape and the maximum of the colloid retention profile.

Strong release of viruses in fracture flow in response to a perturbation in ionic strength: Filtration/retention tests and modeling

Effluents derived from a municipal wastewater treatment plant were used for virus filtration/retention experiments by using a horizontal laboratory filter. Filtration tests were performed to examine how soil geochemical heterogeneity and fracture patterns affected the transport of viruses in groundwater in order to model the influence of reductive perturbations in ionic strength (IS) during wastewater filtration. Although perturbations of IS and velocity are known to result in resuspension of colloids, we found that the effect of soil geochemical heterogeneity can produce strong and instantaneous virus releases in fractured aquifers, likely an internal additional source of viruses. Sixteen limestone slabs were packed in a PVC box filter at the Bari Laboratory (South Italy) to replicate wastewater filtration throughout a fractured medium similar to the Bari carbonate aquifer. Terra rossa, which is an aggregate of sand, silt and clay, was unevenly spread on the surface of each limestone slab within the filter. Since the mineralogical composition of terra rossa includes iron (hematite, magnetite, and goethite) oxides, the soil exhibited localized unfavorable colloid/collector interactions for attachment. In contrast, soil-free parts of the fracture surfaces maintained favorable colloid/collector interactions. We found in our experiments that the lowering of IS due to the reduction of water salt content, which could occur during runoff injections after rainfall, might be sufficient to cause strong detachment of viruses from fracture surfaces, allowing further migration into the groundwater. The model in this work can predict the count and pathways of released viruses in groundwater fractures under soil geochemical heterogeneity and originated by reductions of IS, by using analytical solutions.

Influence of Biochar on Deposition and Release of Clay Colloids in Saturated Porous Media

Although the potential application of biochar in soil remediation has been recognized, the effect of biochar on the transport of clay colloids, and accordingly the fate of colloid-associated contaminants, is unclear to date. This study conducted saturated column experiments to systematically examine transport of clay colloids in biochar-amended sand porous media in different electrolytes at different ionic strengths. The obtained breakthrough curves were simulated by the convection-diffusion equation, which included a first-order deposition and release terms. The deposition mechanisms were interpreted by calculating Derjaguin-Landau-Verwey-Overbeek interaction energies. A linear relationship between the simulated deposition rate or the attachment efficiency and the fraction of biochar was observed ( ≥ 0.91), indicating more favorable deposition in biochar than in sand. The interaction energy calculations show that the greater deposition in biochar occurs because the half-tube-like cavities on the biochar surfaces favor deposition in secondary minima and the nanoscale physical and chemical heterogeneities on the biochar surfaces increase deposition in primary minima. The deposited clay colloids in NaCl can be released by reduction of ionic strength, whereas the presence of a bivalent cation (Ca) results in irreversible deposition due to the formation of cation bridging between the colloids and biochar surfaces. The deposition and release of clay colloids on or from biochar surfaces not only change their mobilizations in the soil but also influence the efficiency of the biochar for removal of pollutants. Therefore, the influence of biochar on clay colloid transport must be considered before application of the biochar in soil remediation.

Modeling of retention and re-entrainment of mono- and poly-disperse particles: Effects of hydrodynamics, particle size and interplay of different-sized particles retention

In this paper, numerical simulations of experimental data were performed with kinetic rate coefficients to characterize the retention and re-entrainment dynamics under different hydrodynamic conditions for monodisperse and polydisperse latex particles (3, 10, 16μm and the mixture). The results show that drastic increase in fluid velocity provokes hardly any remarkable decrease in retention in the presence of large energy barriers (>2000kT). Systematical increases in deposition and re-entrainment dynamic rates were observed with fluid velocity and/or particle size. Increased irreversible deposition rate indicates straining and wedging dominate deposition in this study. Excess retention of 3μm particle in the polydisperse particle suspension was observed. The origins are reckoned that deposited larger particles may hinder the re-entrainment of smaller particles near the grain-to-grain contact and can provide additional sites of attachment.

Colloid characterization and in situ release in shallow groundwater under different hydrogeology conditions

Changes to groundwater hydrodynamics and chemistry can lead to colloid release that can have a major impact on the groundwater environment. To analyze the effects of colloid release caused by artificial groundwater recharge, field and laboratory tests on colloid characterization and colloid release were conducted. The field tests were carried out at an artificial recharge test site in Shandong Province. In the field investigation, one recharge water sample and five groundwater samples were collected and filtered through three levels of ultrafiltration membranes, with pore sizes of 0.45 μm, 100 kDa, and 50 kDa. The field results indicated that the colloid mass concentrations in groundwater retained between membranes with pore sizes of 100 kDa-0.45 μm and 50 kDa-100 kDa were 19 and 62 mg/L, respectively. In recharge water, the colloid mass concentrations retained by 100-kDa-0.45-μm and 50-kDa-100-kDa membranes were 3 and 99 mg/L, respectively. Colloids detected on the ultrafiltration membranes were mainly inorganic between 100 kDa and 0.45 μm, and mainly organic between 50 and 100 kDa. Based on the field colloid investigation results, the organic colloid was chosen in the laboratory experiments to reveal its release behavior under different conditions. Porous media diameter, flux, ionic strength (IS), and ion valence were changed to determine their influences on organic colloid concentration outflow from undisturbed porous media. The experiment's results indicate that decreasing the diameter, and increasing the flux, ionic strength, and the number of divalent cations, can promote organic colloid release. The organic colloid release rate in the early stage was high and is thus likely to affect the quality of groundwater. The results provide a useful scientific basis for minimizing changes to hydrodynamic and hydrochemical conditions during artificial recharge, thus safeguarding groundwater quality.

Synergies of media surface roughness and ionic strength on particle deposition during filtration

Although it is widely believed that media/collector roughness can enhance particle deposition on surfaces, this effect has not been consistently observed nor systematically described. Here, column tests were conducted to: 1) evaluate media roughness impacts on particle deposition in the presence of an energy barrier (i.e., at low ionic strength conditions), and 2) describe the concurrent impacts of collector surface roughness and suspension fluid ionic strength on particle deposition in packed beds. This work presents a first, systematic demonstration that media/collector surface roughness consistently influences particle deposition in a non-linear, non-monotonic manner, irrespective of the presence of an energy barrier. Notably, ionic strength-associated changes in DLVO interaction energy could not solely explain observed differences in particle deposition associated with collector surface roughness. Particle-to-roughness element and particle-to-smooth/bottom surface interactions contributed to a critical roughness size associated with a minimum DLVO interaction energy; however, that critical size is not necessarily the same as the critical size associated with minimal particle deposition rates. Surface roughness and ionic strength concurrently affected particle deposition in a manner that is not simply additive; rather, particle deposition rates were highly correlated with inverse Debye-Hückel length (i.e., ln [κ^-1]) using second-order polynomial functions. Notably, the secondary energy minimum alone appears inadequate for explaining the observed particle deposition behavior. These relationships may provide insight for further development of physico-chemical filtration models for describing particle deposition on surfaces.

Observed Dependence of Colloid Detachment on the Concentration of Initially Attached Colloids and Collector Surface Heterogeneity in Porous Media

Sand column experiments were conducted to examine the effects of the concentration of attached colloids (CAC) on their subsequent detachment upon decreasing solution ionic strength (IS). Different pore volumes of latex microparticle suspensions were injected into the columns to allow different amounts of colloids to attach at ISs of 0.001, 0.01, and 0.2 M. Then, deionized water was introduced to release the attached colloids. Results show that the fraction of attachments that were reversible to reduction of IS (FRA) increased with increasing CAC at a given IS if the sand was extensively treated using acids to reduce surface charge heterogeneity. This indicates that colloids were preferentially immobilized in sites favoring irreversible attachment and then gradually occupied reversible sites. In contrast, the FRA decreased with increasing CAC at 0.001 M in sand without the acid treatment, illustrating the opposite attachment sequence. Scanning electron microscope examinations reveal that the concave regions favored irreversible colloid attachment. Reversible attachment is likely due to immobilization on flat surfaces with charge heterogeneities, retention in stagnation point regions via secondary minimum association, ripening in the acid-treated sand, and capture of colloids by protruding asperities with charge heterogeneity in the untreated sand. At ISs of 0.01 and 0.2 M, the FRA was essentially independent of CAC in the untreated sand because the colloids were randomly attached on the sand surfaces over time.

Contribution of Nano- to Microscale Roughness to Heterogeneity: Closing the Gap between Unfavorable and Favorable Colloid Attachment Conditions

Surface roughness has been reported to both increase as well as decrease colloid retention. In order to better understand the boundaries within which roughness operates, attachment of a range of colloid sizes to glass with three levels of roughness was examined under both favorable (energy barrier absent) and unfavorable (energy barrier present) conditions in an impinging jet system. Smooth glass was found to provide the upper and lower bounds for attachment under favorable and unfavorable conditions, respectively. Surface roughness decreased, or even eliminated, the gap between favorable and unfavorable attachment and did so by two mechanisms: (1) under favorable conditions attachment decreased via increased hydrodynamic slip length and reduced attraction and (2) under unfavorable conditions attachment increased via reduced colloid-collector repulsion (reduced radius of curvature) and increased attraction (multiple points of contact, and possibly increased surface charge heterogeneity). Absence of a gap where these forces most strongly operate for smaller (<200 nm) and larger (>2 μm) colloids was observed and discussed. These observations elucidate the role of roughness in colloid attachment under both favorable and unfavorable conditions.

Detachment of fullerene nC60 nanoparticles in saturated porous media under flow/stop-flow conditions: Column experiments and mechanistic explanations

This study was aimed at investigating the detachment of fullerene nC60 nanoparticles (NPs) in saturated sand porous media under transient and static conditions. The nC60 NPs were first attached at primary minima of Derjaguin-Landau-Verwey-Overbeek (DLVO) interaction energy profiles in electrolyte solutions with different ionic strengths (ISs). The columns were then eluted with deionized water to initiate nC60 NP detachment by decreasing solution IS. Finally, the flow of the columns was periodically interrupted to investigate nC60 NP detachment under static condition. Our results show that the detachment of nC60 NPs occurred under both transient and static conditions. The detachment under transient conditions was attributed to the fact that the attractions acting on the nC60 NPs at primary minima were weakened by nanoscale physical heterogeneities and overcome by hydrodynamic drags at lower ISs. However, a fraction of nC60 NPs remained at shallow primary minima in low flow regions, and detached via Brownian diffusion during flow interruptions. Greater detachment of nC60 NPs occurred under both transient and static conditions if the NPs were initially retained in electrolyte solutions with lower valent cations due to lower attractions between the NPs and collectors. Decrease in collector surface chemical heterogeneities and addition of dissolved organic matter also increased the extent of detachment by increasing electrostatic and steric repulsions, respectively. While particle attachment in and subsequent detachment from secondary minima occur in the same electrolyte solution, our results indicate that perturbation in solution chemistry is necessary to lower the primary minimum depths to initiate spontaneous detachment from the primary minima. These findings have important implications for predicting the fate and transport of nC60 NPs in subsurface environments during multiple rainfall events and accordingly for accurately assessing their environmental risks.

Concurrent Modeling of Hydrodynamics and Interaction Forces Improves Particle Deposition Predictions

It is widely believed that media surface roughness enhances particle deposition-numerous, but inconsistent, examples of this effect have been reported. Here, a new mathematical framework describing the effects of hydrodynamics and interaction forces on particle deposition on rough spherical collectors in absence of an energy barrier was developed and validated. In addition to quantifying DLVO force, the model includes improved descriptions of flow field profiles and hydrodynamic retardation functions. This work demonstrates that hydrodynamic effects can significantly alter particle deposition relative to expectations when only the DLVO force is considered. Moreover, the combined effects of hydrodynamics and interaction forces on particle deposition on rough, spherical media are not additive, but synergistic. Notably, the developed model's particle deposition predictions are in closer agreement with experimental observations than those from current models, demonstrating the importance of inclusion of roughness impacts in particle deposition description/simulation. Consideration of hydrodynamic contributions to particle deposition may help to explain discrepancies between model-based expectations and experimental outcomes and improve descriptions of particle deposition during physicochemical filtration in systems with nonsmooth collector surfaces.

Spontaneous Detachment of Colloids from Primary Energy Minima by Brownian Diffusion

The Derjaguin-Landau-Verwey-Overbeek (DLVO) interaction energy profile has been frequently used to interpret the mechanisms controlling colloid attachment/detachment and aggregation/disaggregation behavior. This study highlighted a type of energy profile that is characterized by a shallow primary energy well (i.e., comparable to the average kinetic energy of a colloid) at a small separation distance and a monotonic decrease of interaction energy with separation distance beyond the primary energy well. This energy profile is present due to variations of height, curvature, and density of discrete physical heterogeneities on collector surfaces. The energy profile indicates that colloids can be spontaneously detached from the shallow primary energy well by Brownian diffusion. The spontaneous detachment from primary minima was unambiguously confirmed by conducting laboratory column transport experiments involving flow interruptions for two model colloids (polystyrene latex microspheres) and engineered nanoparticles (fullerene C60 aggregates). Whereas the spontaneous detachment has been frequently attributed to attachment in secondary minima in the literature, our study indicates that the detached colloids could be initially attached at primary minima. Our study further suggests that the spontaneous disaggregation from primary minima is more significant than spontaneous detachment because the primary minimum depth between colloid themselves is lower than that between a colloid and a collector surface.

Determining Parameters and Mechanisms of Colloid Retention and Release in Porous Media

A modeling framework is presented to determine fundamental parameters and controlling mechanisms of colloid (microbes, clays, and nanoparticles) retention and release on surfaces of porous media that exhibit wide distributions of nanoscale chemical heterogeneity, nano- to microscale roughness, and pore water velocity. Primary and/or secondary minimum interactions in the zone of electrostatic influence were determined over the heterogeneous solid surface. The Maxwellian kinetic energy model was subsequently employed to determine the probability of immobilization and diffusive release of colloids from each of these minima. In addition, a balance of applied hydrodynamic and resisting adhesive torques was conducted to determine locations of immobilization and hydrodynamic release in the presence of spatially variable water flow and microscopic roughness. Locations for retention had to satisfy both energy and torque balance conditions for immobilization, whereas release could occur either due to diffusion or hydrodynamics. Summation of energy and torque balance results over the elementary surface area of the porous medium provided estimates for colloid retention and release parameters that are critical to predicting environmental fate, including the sticking and release efficiencies and the maximum concentration of retained colloids on the solid phase. Nanoscale roughness and chemical heterogeneity produced localized primary minimum interactions that controlled long-term retention, even when mean chemical conditions were unfavorable. Microscopic roughness played a dominant role in colloid retention under low ionic strength and high hydrodynamic conditions, especially for larger colloids.

Equilibrium and kinetic models for colloid release under transient solution chemistry conditions

We present continuum models to describe colloid release in the subsurface during transient physicochemical conditions. Our modeling approach relates the amount of colloid release to changes in the fraction of the solid surface area that contributes to retention. Equilibrium, kinetic, equilibrium and kinetic, and two-site kinetic models were developed to describe various rates of colloid release. These models were subsequently applied to experimental colloid release datasets to investigate the influence of variations in ionic strength (IS), pH, cation exchange, colloid size, and water velocity on release. Various combinations of equilibrium and/or kinetic release models were needed to describe the experimental data depending on the transient conditions and colloid type. Release of Escherichia coli D21g was promoted by a decrease in solution IS and an increase in pH, similar to expected trends for a reduction in the secondary minimum and nanoscale chemical heterogeneity. The retention and release of 20nm carboxyl modified latex nanoparticles (NPs) were demonstrated to be more sensitive to the presence of Ca(2+) than D21g. Specifically, retention of NPs was greater than D21g in the presence of 2mM CaCl2 solution, and release of NPs only occurred after exchange of Ca(2+) by Na(+) and then a reduction in the solution IS. These findings highlight the limitations of conventional interaction energy calculations to describe colloid retention and release, and point to the need to consider other interactions (e.g., Born, steric, and/or hydration forces) and/or nanoscale heterogeneity. Temporal changes in the water velocity did not have a large influence on the release of D21g for the examined conditions. This insensitivity was likely due to factors that reduce the applied hydrodynamic torque and/or increase the resisting adhesive torque; e.g., macroscopic roughness and grain-grain contacts. Our analysis and models improve our understanding and ability to describe the amounts and rates of colloid release and indicate that episodic colloid transport is expected under transient physicochemical conditions.

Prediction of Nanoparticle and Colloid Attachment on Unfavorable Mineral Surfaces Using Representative Discrete Heterogeneity

Despite several decades of research there currently exists no mechanistic theory to predict colloid attachment in porous media under environmental conditions where colloid-collector repulsion exists (unfavorable conditions for attachment). It has long been inferred that nano- to microscale surface heterogeneity (herein called discrete heterogeneity) drives colloid attachment under unfavorable conditions. Incorporating discrete heterogeneity into colloid-collector interaction calculations in particle trajectory simulations predicts colloid attachment under unfavorable conditions. As yet, discrete heterogeneity cannot be independently measured by spectroscopic or other approaches in ways directly relevant to colloid-surface interaction. This, combined with the fact that a given discrete heterogeneity representation will interact differently with differently sized colloids as well as different ionic strengths for a given sized colloid, suggests a strategy to back out representative discrete heterogeneity by a comparison of simulations to experiments performed across a range of colloid size, solution IS, and fluid velocity. This has recently been performed for interaction of carboxylate-modified polystyrene latex (CML) microsphere attachment to soda lime glass at pH 6.7 with NaCl electrolyte. However, extension to other surfaces, pH values, and electrolytes is needed. For this reason, the attachment of CML (0.25, 1.1, and 2.0 μm diameters) from aqueous suspension onto a variety of unfavorable mineral surfaces (soda lime glass, muscovite, and albite) was examined for pH values of 6.7 and 8.0), fluid velocities (1.71 × 10(-3) and 5.94 × 10(-3) m s(-1)), IS (6.0 and 20 mM), and electrolytes (NaCl, CaSO4, and multivalent mixtures). The resulting representative heterogeneities (heterodomain size and surface coverage, where heterodomain refers to nano- to microscale attractive domains) yielded colloid attachment predictions that were compared to predictions from existing applicable semiempirical expressions in order to examine the strengths and weaknesses of the discrete heterogeneity approach and opportunities for improvement.

Influence of surface heterogeneities on reversibility of fullerene (nC60) nanoparticle attachment in saturated porous media

This study systematically investigated influence of surface roughness and surface chemical heterogeneity on attachment and detachment of nC60 nanoparticles in saturated porous media by conducting laboratory column experiments. Sand and glass beads were employed as a model collectors to represent a different surface roughness. The two collectors were treated by washing with only deionized water or by using acids to extensively remove chemical heterogeneities. Results show that both attachment and detachment were more in the acid-treated sand than those in the acid-treated glass beads. The greater attachment and detachment were attributed to the reason that sand surfaces have much more nanoscale asperities, which facilitates particle attachment atop of them at primary minima and subsequent detachment upon reduction of ionic strength. No detachment was observed if the water-washed collectors were employed, demonstrating that the couple of chemical heterogeneity with nanoscale roughness causes irreversible attachment in primary minima. Whereas existing studies frequently represented surface rough asperities as regular geometries (e.g., hemisphere, cone, pillar) for estimating influence of surface roughness on Derjaguin-Landau-Verwey-Overbeek (DLVO) interaction energies, our theoretical calculations indicate that the assumptions could underestimate both attachment and detachment because these geometries cannot account for surface curvature effects.