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

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.

In Situ Measurements of Dynamic Bacteria Transport and Attachment in Heterogeneous Sand-Packed Columns

Prevention, mitigation, and regulation of bacterial contaminants in groundwater require a fundamental understanding of the mechanisms of transport and attachment in complex geological materials. Discrepancies in bacterial transport behaviors observed between field studies and laboratory experiments indicate an incomplete understanding of dynamic bacterial transport and immobilization processes in realistic heterogeneous geologic systems. Here, we develop a new experimental approach for in situ quantification of dynamic bacterial transport and attachment distribution in geologic media that relies on radiolabelingEscherichia coliwith positron-emitting radioisotopes and quantifying transport with three-dimensional (3D) positron emission tomography (PET) imaging. Our results indicate that the highest bacterial attachment occurred at the interfaces between sand layers oriented orthogonal to the direction of flow. The predicted bacterial attachment from a 3D numerical model matched the experimental PET results, highlighting that the experimentally observed bacterial transport behavior can be accurately captured with a distribution of a first-order irreversible attachment model. This is the first demonstration of the direct measurement of attachment coefficient distributions from bacterial transport experiments in geologic media and provides a transformational approach to better understand bacterial transport mechanisms, improve model parametrization, and accurately predict how local geologic conditions can influence bacterial fate and transport in groundwater.

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.

Using colloidal AFM probe technique and XDLVO theory to predict the transport of nanoplastics in porous media

The plastic concentration in terrestrial systems is orders of magnitude higher than that found in marine ecosystems, which has raised global concerns about their potential risk to agricultural sustainability. Previous research on the transport of nanoplastics in soil relied heavily on the qualitative prediction of the mean-field extended Derjaguin-Landau-Verwey-Overbeek theory (XDLVO), but direct and quantitative measurements of the interfacial forces between single nanoplastics and porous media are lacking. In this study, we conducted multiscale investigations ranging from column transport experiments to single particle measurements. The maximum effluent concentration (C/C0) of amino-modified nanoplastics (PS-NH₂) was 0.94, whereas that of the carboxyl-modified nanoplastics (PS-COOH) was only 0.33, indicating PS-NH₂ were more mobile than PS-COOH at different ionic strengths (1-50 mM) and pH values (5-9). This phenomenon was mainly attributed to the homogeneous aggregation of PS-COOH. In addition, the transport of PS-NH₂ in the quartz sand column was inhibited with the increase of ionic strength and pH, and pH was the major factor governing their mobility. The transport of PS-COOH was inhibited with increasing ionic strength and decreasing pH. Hydrophilicity/hydrophobicity-mediated interactions and particle heterogeneity strongly interfered with interfacial forces, leading to the qualitative prediction of XDLVO, contrary to experimental observations. Through the combination of XDLVO and colloidal atomic force microscopy, accurate and quantitative interfacial forces can provide compelling insight into the fate of nanoparticles in the soil environment.

Mechanism comparisons of transport-deposition-reentrainment between microplastics and natural mineral particles in porous media: A theoretical and experimental study

The migration and distribution of microplastic particles (MPs) in the natural environment has attracted global attention in recent years. However, little is known about the transport-deposition-reentrainment differences between MPs and natural mineral particles in porous media. In this study, polystyrene (PS) and silica (SiO₂) particles, representing model MPs and natural mineral particles, respectively, were selected to study the responses of different particle types to changes in specific particle size and flow velocity. Three typical particle sizes and various flow velocities were chosen to compare and delineate the transport-deposition-reentrainment characteristics of PS and SiO₂ in a packed-bed laboratory column. Collector efficiency was calculated using Tufenkji and Elimelech (TE) equation. The particle fractions released from the collector surfaces were predicted using DLVO theory and force analysis. Two types of particles were attached in the secondary minimum, which were either retained on the collector surface or reentrained to the fluid. The staged elution experiment wherein the flow velocity was increased experienced a period of flow shock, thus breaking the force balance of the particle. An increase in the flow velocity resulted in various degrees of particle elution. The breakthrough experiment at a specific flow velocity showed that the corresponding velocity alteration in staged elution experiment contributed to reentrainment to varying extents. When the effect of gravity on particle deposition was negligible, the particle size was larger, and the lower the velocity for releasing the particles. However, the opposite tendency was observed when considering the effect of gravity on particle deposition. Moreover, the deposition, mainly due to gravity, easily causes particle reentrainment as the flow velocity increases. This study further predicts and reveals the nature of transport and deposition differences between MPs and natural mineral particles, which helps to further assess the risk and potential of groundwater contamination with MPs of different sizes.

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.

Permanganate oxidation of polycyclic aromatic compounds (PAHs and polar PACs): column experiments with DNAPL at residual saturation

Permanganate is an oxidant usually applied for in situ soil remediation due to its persistence underground. It has already shown great efficiency for dense nonaqueous phase liquid (DNAPL) degradation under batch experiment conditions. In the present study, experimental permanganate oxidation of a DNAPL - coal tar - sampled in the groundwater of a former coking plant was carried out in a glass bead column. Several glass bead columns were spiked with coal tar using the drainage-imbibition method to mimic on-site pollution spread at residual saturation as best as possible. The leaching of organic pollutants was monitored as the columns were flushed by successive sequences: successive injections of hot water, permanganate solution for oxidation, and ambient temperature water, completed by two injections of a tracer before and after oxidation. Sixteen conventional US-EPA PAHs and selected polar PACs were analyzed in the DNAPL remaining in the columns at the end of the experiment and in the particles collected at several steps of the flushing sequences. Permanganate oxidation of the pollutants was rapidly limited by interfacial aging of the DNAPL drops. Moreover, at the applied flow rate chosen to be representative of in situ injections and groundwater velocities, the reaction time was not sufficient to reach high degradation yields but induced the formation and the leaching of oxygenated PACs.

Delineating the Relationship between Nanoparticle Attachment Efficiency and Fluid Flow Velocity

The ability to fundamentally describe nanoparticle (NP) transport in the subsurface underpins environmental risk assessment and successful material applications, including advanced remediation and sensing technologies. Despite considerable progress, our understanding of NP deposition behavior remains incomplete as there are conflicting reports regarding the effect of fluid flow velocity on attachment efficiency. To directly address this and more accurately describe NP attachment behavior, we have developed a novel protocol using a quartz crystal microbalance with dissipation monitoring (QCM-D) to separate and individually observe deposition mechanisms (diffusion and sedimentation), providing in situ, real-time information about particle diffusion (from the bulk liquid to solid surface). Through this technique, we have verified that the approaching velocity of NPs via diffusion increases (0.8-6.7 μm/s) with increasing flow velocity (6.1-106.0 μm/s), leading to an increased NP kinetic energy, thus affecting deposition processes. Further, in the presence of a secondary energy minimum associated with organic surface coatings, secondary minimum deposition decreases and primary minimum deposition increases with the flow velocity. NPs deposited at the primary minimum are relatively more resistant to hydrodynamic energies (including detachment associated energies), resulting in an increase of observed attachment efficiencies. Taken together, this work not only describes a novel method to delineate and quantify physical processes underpinning particle behavior but also provides direct measurements regarding key factors defining the relationship(s) of flow velocity and particle attachment. Such insight is valuable for next-generation fate and transport model accuracy, especially under unfavorable attachment regimes, which is a current and critical need for subsurface material applications and implication paradigms.

Multiscale dynamics of colloidal deposition and erosion in porous media

Diverse processes-e.g., environmental pollution, groundwater remediation, oil recovery, filtration, and drug delivery-involve the transport of colloidal particles in porous media. Using confocal microscopy, we directly visualize this process in situ and thereby identify the fundamental mechanisms by which particles are distributed throughout a medium. At high injection pressures, hydrodynamic stresses cause particles to be continually deposited on and eroded from the solid matrix-notably, forcing them to be distributed throughout the entire medium. By contrast, at low injection pressures, the relative influence of erosion is suppressed, causing particles to localize near the inlet of the medium. Unexpectedly, these macroscopic distribution behaviors depend on imposed pressure in similar ways for particles of different charges, although the pore-scale distribution of deposition is sensitive to particle charge. These results reveal how the multiscale interactions between fluid, particles, and the solid matrix control how colloids are distributed in a porous medium.

Synergies of surface roughness and hydration on colloid detachment in saturated porous media: Column and atomic force microscopy studies

Saturated column experiments were conducted to systematically examine the influence of hydration on the detachment of nano- and micro-sized latex colloids (35 nm and 1 μm, respectively) from sand. The colloids were attached on the sand in primary minima (PM) using high ionic strength (IS) NaCl solutions. The PM were predicted to be shallower and located farther from sand surfaces with increasing IS due to the hydration force. Consequently, a greater amount of colloid detachment occurred in deionized water when the colloids were initially deposited at a higher IS. Atomic force microscopy (AFM) examinations showed that both nanoscale protruding asperities and large wedge-like valleys existed on the sand surface. The influence of these surface features on the interaction energies/forces was modeled by approximating the roughness as cosinoidal waves and two intersecting half planes, respectively. The PM were deep and attachment was irreversible at concave regions for all ISs, even if the hydration force was included. Conversely, colloids were weakly attached at protruding asperities due to a reduced PM depth, and thus were responsible for the detachment upon IS reduction. The AFM examinations confirmed that the adhesive forces were enhanced and reduced (or even completely eliminated) at concave and convex locations of sand surfaces, respectively. These results have important implications for surface cleaning and prediction of the transport and fate of hazardous colloids and colloid-associated contaminants in subsurface environments.

Quantitative Linking of Nanoscale Interactions to Continuum-Scale Nanoparticle and Microplastic Transport in Environmental Granular Media

Quantitative linkage of fundamental physicochemical characteristics to rate coefficients used in simulations of experimentally observed transport behaviors of nanoparticles and microplastics (colloids) in environmental granular media is an active area of research. Quantitative linkage is herein demonstrated for (i) colloids ranging from nano- to microscale; in two field-based granular media of contrasting grain size, (ii) natural fine sand at the column scale; and (ii) streambed-equilibrated commercial pea gravel at the field scale. Continuum-scale rate coefficients were linked to nanoscale interactions via mechanistic pore-scale colloid trajectory simulations that predicted and defined fast- and slow-attaching subpopulations, as well as nonattaching subpopulations that either remained in the near-surface pore water or re-entrained to bulk pore water. These subfractions of the classic collector efficiency were upscaled to continuum-scale rate coefficients that produced experimentally observed colloid breakthrough-elution concentration histories and nonexponential colloid distributions from the source. The simulations explained transition from hyperexponential to nonmonotonic colloid distributions from the source as driven accumulation of mobile near-surface colloids due to relatively strong secondary minimum interaction and weak diffusion for microscale colloids. The assumption of depletion of the fast-attaching colloid subpopulation by attachment to grain surfaces produced the experimentally observed contrasting distances across which nonexponential colloid distribution from the source occurred in the fine sand versus pea gravel. Rate coefficients were quantitatively calculated from physicochemical parameters and the following three fit parameters: (i) fractional coverage by nanoscale heterogeneity; (ii) efficiency of return to the near-surface domain; and (iii) in explicit simulations, characteristic velocity for scaling transfer to near-surface pore water.

Investigation for Synergies of Ionic Strength and Flow Velocity on Colloidal-Sized Microplastic Transport and Deposition in Porous Media Using the Colloidal-AFM Probe

Studies that explore the transport and retention behavior of colloidal-sized microplastic (MP) with focusing on the governing mechanisms for their attachment and detachment process using colloidal-atomic force microscopy (C-AFM) were still limited. In the present study, multiscale investigations ranging from pore-scale column test to microscale visualization and eventually to nanoscale interfacial and adhesive force measurement were conducted. Pore- and microscale tests were conducted at various flow velocity and over a broad range of IS values and found that IS and flow velocity could synergically impact the deposition of MPs during filtration, in particular under unfavorable condition at small flow velocity. The net difference between the highest and lowest deposition conditions became smaller while flow velocity was decreasing in porous media. However, the net difference between the high and low IS conditions in parallel plate chamber were not sensitive to the change of flow velocity. The measurement from C-AFM suggested that not only the interfacial force but also the adhesive forces changed while MP was approaching/retracting to the collector surface. Information related to the magnitude, location, and occurrence of interfacial/adhesive forces were analyzed. Comparisons of the interaction energy determined from the measured force and ones derived from surface energy components using DLVO theory were conducted to explain the synergies of IS and flow velocity on pathogenic size MPs transport and deposition during filtration.