Ionic strength dependent transport of microparticles in saturated porous media: modeling mobilization and immobilization phenomena under transient chemical conditions

Mobility of crocidolite asbestos in sandy porous media mimicking aquifer systems

Asbestos is widely recognized as being a carcinogen when dispersed in air, but very little is known about its exposure pathways in water and its subsequent effects on human health. Several studies have proved asbestos presence in groundwater but failed to assess its mobility in aquifer systems. This paper aims to fill this gap by studying the transport of crocidolite, an amphibole asbestos, through sandy porous media mimicking different aquifer systems. To this purpose, two sets of column test were performed varying the crocidolite suspension concentration, the quartz sand grain size distribution, and the physicochemical water parameters (i.e., pH). The results proved that crocidolite is mobile in quartz sand due to the repulsive interactions between fibres and porous media. The concentration of fibres at the outlet of the column were found to decrease when decreasing the grain size distribution of the porous medium, with a bigger impact on highly concentrated suspensions. In particular, 5-to-10-µm-long fibres were able to flow through all the tested sands while fibres longer than 10 µm were mobile only through the coarser medium. These results confirm that groundwater migration should be considered a potential exposure pathway while implementing human health risk assessment.

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

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

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.

Nanoplastics removal during drinking water treatment: Laboratory- and pilot-scale experiments and modeling

Microplastics detected in potable water sources and tap water have led to concerns about the efficacy of current drinking water treatment processes to remove these contaminants. It is hypothesized that drinking water resources contain nanoplastics (NPs), but the detection of NPs is challenging. We, therefore, used palladium (Pd)-labeled NPs to investigate the behavior and removal of NPs during conventional drinking water treatment processes including ozonation, sand and activated carbon filtration. Ozone doses typically applied in drinking water treatment plants (DWTPs) hardly affect the NPs transport in the subsequent filtration systems. Amongst the different filtration media, NPs particles were most efficiently retained when aged (i.e. biofilm coated) sand was used with good agreements between laboratory and pilot scale systems. The removal of NPs through multiple filtration steps in a municipal full-scale DWTP was simulated using the MNMs software code. Removal efficiencies exceeding 3-log units were modeled for a combination of three consecutive filtration steps (rapid sand filtration, activated carbon filtration and slow sand filtration with 0.4-, 0.2- and 3.0-log-removal, respectively). According to the results from the model, the removal of NPs during slow sand filtration dominated the overall NPs removal which is also supported by the laboratory-scale and pilot-scale data. The results from this study can be used to estimate the NPs removal efficiency of typical DWTPs with similar water treatment chains.

Measuring and modeling the influence of salinity change on the transport behaviour of Escherichia coli through quartz sand

Pathogenic bacteria can be discharged in the environment through natural as well as anthropogenic activities. Once in the environment, they may contaminate soil and sediments and migrate towards water bodies. Transient chemical conditions may occur in soil/sediments and favor mobilization of bacteria, e.g., upon the reduction of salinity (or ionic strength). However, the magnitude of this phenomenon and its relationship with particle size is not well understood, yet. In this work, we investigated the transport of Escherichia coli under variable salinity conditions (between 1 and 20 part per thousand, ppt) and for different soil grain sizes (between 150 and 710 μm). A model developed in our group was applied in this work. It couples bacteria and salinity transport equations in order to account for transient water composition in the description of bacteria migration. The model was calibrated and validated with laboratory experiments. The tests were monitored continuously with UV-Vis spectroscopy, which allowed to record highly resolved concentration fronts. The results show that salinity increases the retardation of the bacteria. Upon salinity drop, a release of bacteria occurs forming a peak whose magnitude increases with salinity change. This effect becomes more important as the grain size decreases. Simulations suggest that the dominant retention mechanism is attachment for coarse sand and straining for fine sand. The retention can be reversed as the salinity is reduced causing a sudden bacteria mobilization. Such a behaviour may have important implications on microbial contamination of water bodies when soil/sediments undergo transient chemical conditions.

Bentonite colloids immobilization and release in quartz column and its influence on selenite migration

Immobilization and release of colloids are important for colloids-facilitated migrations, and in the safety assessment of geological disposal for high-level radioactive waste, the association between the immobilization and release process of the bentonite colloids with selenite migration has not been well revealed. In this work, the migration of bentonite colloids under different conditions is evaluated, and the effects of colloids immobilization and release on selenite migration are studied. In addition, the cases of in-migration (colloids are immobilized in the quartz sand, and then selenite migrates through the quartz sand with immobilized colloids) and co-migration (colloids bearing selenite are immobilized in the quartz sand) are investigated. The results show that in the systems containing 3.0 mM Mg²⁺, the mobility of the colloids is highly hindered and the colloids are immobilized in the quartz sand mainly by straining effect. The immobilization of bentonite colloids affects selenite migration differently according to the immobilization process (in-migration or co-migration). A more significant retardation effect is observed in the co-migration process than in-migration due to the additional inner-sphere complexed selenite in the co-migration. The immobilized colloids can be more easily released by alkaline DI-water (pH 11.0) than acidic one (pH 6.0) as a result of the more negative surface charges of the immobilized bentonite colloids. The average size of the released colloids is larger than the initial colloids at the same pH. Selenite is found to be released ahead of colloids in either in- or co-migration process, and part of selenite is discovered migrating with released colloids in co-migration process. Since colloids immobilization and release would influence radionuclides migration, further research about colloids immobilization and release with broad range of pH and ionic strength in the host rock and its influence on the migration of other radionuclides are needed.

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.

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 Large-Scale 3D Study on Transport of Humic Acid-Coated Goethite Nanoparticles for Aquifer Remediation

Humic acid-coated goethite nanoparticles (HA-GoeNPs) have been recently proposed as an effective reagent for the in situ nanoremediation of contaminated aquifers. However, the effective dosage of these particles has been studied only at laboratory scale to date. This study investigates the possibility of using HA-GoeNPs in remediation of real field sites by mimicking the injection and transport of HA-GoeNPs under realistic conditions. To this purpose, a three-dimensional (3D) transport experiment was conducted in a large-scale container representing a heterogeneous unconfined aquifer. Monitoring data, including particle size distribution, total iron (Fetot) content and turbidity measurements, revealed a good subsurface mobility of the HA-GoeNP suspension, especially within the higher permeability zones. A radius of influence of 2 m was achieved, proving that HA-GoeNPs delivery is feasible for aquifer restoration. A flow and transport model of the container was built using the numerical code Micro and Nanoparticle transport Model in 3D geometries (MNM3D) to predict the particle behavior during the experiment. The agreement between modeling and experimental results validated the capability of the model to reproduce the HA-GoeNP transport in a 3D heterogeneous aquifer. Such result confirms MNM3D as a valuable tool to support the design of field-scale applications of goethite-based nanoremediation.

Key factors affecting graphene oxide transport in saturated porous media

This study focuses on the transport in porous media of graphene oxide nanoparticles (GONP) under conditions similar to those applied in the generation of in-situ reactive zones for groundwater remediation (i.e. GO concentration of few tens of mg/l, stable suspension in alkaline solution). The experimental tests evaluated the influence on GO transport of three key factors, namely particle size (300-1200 nm), concentration (10-50 mg/L), and sand size (coarse to fine). Three sources of GONP were considered (two commercial and one synthesized in the laboratory). Particles were stably dispersed in water at pH 8.5 and showed a good mobility in the porous medium under all experimental conditions: after injection of 5 pore volumes and flushing, the highest recovery was around 90%, the lowest around 30% (only for largest particles in fine sand). The particle size was by far the most impacting parameter, with increasing mobility with decreasing size, even if sand size and particle concentration were also relevant. The source of GONP showed a minor impact on the mobility. The transport test data were successfully modeled using the advection-dispersion-deposition equations typically applied for spherical colloids. Experimental and modeling results suggested that GONP, under the explored conditions, are retained due to both blocking and straining, the latter being relevant only for large particles and/or fine sand. The findings of this study play a key role in the development of an in-situ groundwater remediation technology based on the injection of GONP for contaminant degradation or sorption. Despite their peculiar shape, GONP behavior in porous media is comparable with spherical colloids, which have been more studied by far. In particular, the possibility of modeling GONP transport using existing models ensures that they can be applied also for the design of field-scale injections of GONP, similarly to other particles already used in nanoremediation.

Studying key processes related to CO2 underground storage at the pore scale using high pressure micromodels

Micromodels experimentation for studying and understanding CO₂ geological storage mechanisms at the pore scale.

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 ᅟ.

Human health risk assessment for nanoparticle-contaminated aquifer systems

Nanosized particles (NPs), such as TiO₂, Silver, graphene NPs, nanoscale zero-valent iron, carbon nanotubes, etc., are increasingly used in industrial processes, and releases at production plants and from landfills are likely scenarios for the next years. As a consequence, appropriate procedures and tools to quantify the risks for human health associated to these releases are needed. The tiered approach of the standard ASTM procedure (ASTM-E2081-00) is today the most applied for human health risk assessment at sites contaminated by chemical substances, but it cannot be directly applied to nanoparticles: NP transport along migration pathways follows mechanisms significantly different from those of chemicals; moreover, also toxicity indicators (namely, reference dose and slope factor) are NP-specific. In this work a risk assessment approach modified for NPs is proposed, with a specific application at Tier 2 to migration in groundwater. The standard ASTM equations are modified to include NP-specific transport mechanisms. NPs in natural environments are typically characterized by a heterogeneous set of NPs having different size, shape, coating, etc. (all properties having a significant impact on both mobility and toxicity). To take into account this heterogeneity, the proposed approach divides the NP population into classes, each having specific transport and toxicity properties, and simulates them as independent species. The approach is finally applied to a test case simulating the release of heterogeneous Silver NPs from a landfill. The results show that taking into account the size-dependent mobility of the particles provides a more accurate result compared to the direct application of the standard ASTM procedure. In particular, the latter tends to underestimate the overall toxic risk associated to the nP release.

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.

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.

A 3-dimensional micro- and nanoparticle transport and filtration model (MNM3D) applied to the migration of carbon-based nanomaterials in porous media

Engineered nanoparticles (NPs) in the environment can act both as contaminants, when they are unintentionally released, and as remediation agents when injected on purpose at contaminated sites. In this work two carbon-based NPs are considered, namely CARBO-IRON®, a new material developed for contaminated site remediation, and single layer graphene oxide (SLGO), a potential contaminant of the next future. Understanding and modeling the transport and deposition of such NPs in aquifer systems is a key aspect in both cases, and numerical models capable to simulate NP transport in groundwater in complex 3D scenarios are necessary. To this aim, this work proposes a modeling approach based on modified advection-dispersion-deposition equations accounting for the coupled influence of flow velocity and ionic strength on particle transport. A new modeling tool (MNM3D - Micro and Nanoparticle transport Model in 3D geometries) is presented for the simulation of NPs injection and transport in 3D scenarios. MNM3D is the result of the integration of the numerical code MNMs (Micro and Nanoparticle transport, filtration and clogging Model - Suite) in the well-known transport model RT3D (Clement et al., 1998). The injection in field-like conditions of CARBO-IRON® (20g/l) amended by CMC (4g/l) in a 2D vertical tank (0.7×1.0×0.12m) was simulated using MNM3D, and compared to experimental results under the same conditions. Column transport tests of SLGO at a concentration (10mg/l) representative of a possible spill of SLGO-containing waste water were performed at different values of ionic strength (0.1 to 35mM), evidencing a strong dependence of SLGO transport on IS, and a reversible blocking deposition. The experimental data were fitted using the numerical code MNMs and the ionic strength-dependent transport was up-scaled for a full scale 3D simulation of SLGO release and long-term transport in a heterogeneous aquifer. MNM3D showed to potentially represent a valid tool for the prediction of the long-term behavior of engineered nanoparticles released in the environment (e.g. from landfills), and the preliminary design of in situ aquifer remediation through injection of suspensions of reactive NPs.

Assessing the capacity of zero valent iron nanofluids to remediate NAPL-polluted porous media

A variety of aqueous suspensions (nanofluids) of zero-valent nano-particles (nZVI) are prepared by wet chemistry techniques, their stability and longevity is evaluated by physic-chemical methods of characterization, and their reactivity toward the dechlorination of per-chloro-ethylene (PCE) is examined with tests in batch reactors. For assessing the mobility, longevity and reactivity of nZVI suspensions (nanofluids), under flow-through conditions, visualization multiphase flow and transport tests are performed on a glass-etched pore network. The nZVI breakthrough curves are constructed by measuring the transient variation of the iron concentration in the effluent with atomic absorption spectroscopy. The capacity of nZVI to remediate the bulk phase of PCE is quantified by detecting the mass loss rate of PCE ganglia trapped in glass-etched pore networks during the continuous injection of nZVI suspension or pure water. The nZVI injection in porous media is simulated as an advection- dispersion process by accounting for the attachment/detachment of nanoparticles on the pore-walls, and describing the kinetics of PCE dissolution and reaction by 1st order equations. Visualization experiments reveal that the gradual elimination of PCE ganglia by the injected nZVI is associated with the preferential "erosion" of the upstream interfacial regions. The step controlling the overall process kinetics might be either (i) the enhanced PCE dissolution or (ii) the direct reaction of bulk PCE with the nZVI deposited upon the ganglia interfaces. Inverse modeling of the experiments under the simplifying assumption of one active mechanism indicates that the estimated kinetic coefficients are increasing functions of the flow rate.

Co-transport of Pesticide Acetamiprid and Silica Nanoparticles in Biochar-Amended Sand Porous Media

The role of biochar as a soil amendment on the transport of acetamiprid, a widely used neonicotinoid pesticide, is little known. We conducted saturated column experiments to examine cotransport of acetamiprid and silica nanoparticles (NPs) in pure and biochar-amended sands. Retention of acetamiprid was minor in the pure sand, whereas application of biochar in the sand significantly increased retention. Retention was greater at lower ionic strengths and near neutral pH values and was attributed to biodegradation and sorption through π-π interaction and pore filling. The convection-diffusion equation with inclusion of first-order sorption, desorption, and degradation well described the transport of acetamiprid in the biochar-amended sand. The simulation results show that the sorption rate did not change with pH. This is because the acetamiprid is nonionic and cannot be bonded with the biochar by protonation or deprotonation. The desorption rate was independent of variation of solution chemistry, indicating that desorption was a physical process (i.e., pore diffusion). Application of biochar in the sand had little influence on the transport of silica NPs in NaCl but caused complete attachment in CaCl. Energy dispersive X-ray spectroscopy suggested that the enhanced attachment was due to cation bridging between silica NPs and functional groups in biochar by the Ca. The co-presence of acetamiprid and silica NPs in the solutions enhanced transport of acetamiprid and NPs in the biochar-amended sand by competing for the binding sites on the biochar surfaces.

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.

Mobilization and transport of metal-rich colloidal particles from mine tailings into soil under transient chemical and physical conditions

Exposed mine tailing wastes with considerable heavy metals can release hazardous colloidal particles into soil under transient chemical and physical conditions. Two-layered packed columns with tailings above and soils below were established to investigate mobilization and transport of colloidal particles from metal-rich mine tailings into soil under transient infiltration ionic strength (IS: 100, 20, 2 mM) and flow rate (FR: 20.7, 41, and 62.3 mm h(-1)), with Cu and Pb as representatives of the heavy metals. Results show that the tailing particles within the colloidal size (below 2 μm) were released from the columns. A step-decrease in infiltration IS and FR enhanced, whereas a step-increase in the IS and FR restrained the release of tailing particles from the column. The effects of step-changing FR were unexpected due to the small size of the released tailing particles (220-342 nm, being not sensitive to hydrodynamic shear force), the diffusion-controlled particle release process and the relatively compact pore structure. The tailing particles present in the solution with tested IS were found negatively charged and more stable than soil particles, which provides favorable conditions for tailing particles to be transported over a long distance in the soil. The mobilization and transport of Cu and Pb from the tailings into soil were mediated by the tailing particles. Therefore, the inherent toxic tailing particles could be considerably introduced into soil under certain conditions (IS reduction or FR decrease), which may result in serious environmental pollution.

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

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

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.

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

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

nZVI injection into variably saturated soils: Field and modeling study

Nano-scale zero valent iron (nZVI) has been used at a number of contaminated sites over the last decade. At most of these sites, significant decreases in contaminant concentrations have resulted from the application of nZVI. However, limited work has been completed investigating nZVI field-scale mobility. In this study, a field test was combined with numerical modeling to examine nZVI reactivity along with transport properties in variably saturated soils. The field test consisted of 142L of carboxymethyle cellulose (CMC) stabilized monometallic nZVI synthesized onsite and injected into a variably saturated zone. Periodic groundwater samples were collected from the injection well, as well as, from two monitoring wells to analyze for chlorinated solvents and other geochemistry indicators. This study showed that CMC stabilized monometallic nZVI was able to decrease tricholorethene (TCE) concentrations in groundwater by more than 99% from the historical TCE concentrations. A three dimensional, three phase, finite difference numerical simulator, (CompSim) was used to further investigate nZVI and polymer transport at the variably saturated site. The model was able to accurately predict the field observed head data without parameter fitting. In addition, the numerical simulator estimated the mass of nZVI delivered to the saturated and unsaturated zones and distinguished the nZVI phase (i.e. aqueous or attached). The simulation results showed that the injected slurry migrated radially outward from the injection well, and therefore nZVI transport was governed by injection velocity and viscosity of the injected solution. A suite of sensitivity analyses was performed to investigate the impact of different injection scenarios (e.g. different volume and injection rate) on nZVI migration. Simulation results showed that injection of a higher nZVI volume delivered more iron particles at a given distance; however, the travel distance was not proportional to the increase in volume. Moreover, simulation results showed that using a 1D transport equation to simulate nZVI migration in the subsurface may overestimate the travel distance. This is because the 1D transport equation assumes a constant velocity while pore water velocity radially decreases from the well during injection. This study suggests that on-site synthesized nZVI particles are mobile in the subsurface and that a numerical simulator can be a valuable tool for optimal design of nZVI field applications.

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.

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

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

Attachment and Detachment Behavior of Human Adenovirus and Surrogates in Fine Granular Limestone Aquifer Material

The transport of human adenovirus, nanoparticles, and PRD1 and MS2 bacteriophages was tested in fine granular limestone aquifer material taken from a borehole at a managed aquifer recharge site in Adelaide, South Australia. Comparison of transport and removal of virus surrogates with the pathogenic virus is necessary to understand the differences between the virus and surrogate. Because experiments using pathogenic viruses cannot be done in the field, laboratory tests using flow-through soil columns were used. Results show that PRD1 is the most appropriate surrogate for adenovirus in an aquifer dominated by calcite material but not under high ionic strength or high pH conditions. It was also found that straining due to size and the charge of the colloid were not dominant removal mechanisms in this system. Implications of this study indicate that a certain surrogate may not represent a specific pathogen solely based on similar size, morphology, and/or surface charge. Moreover, if a particular surrogate is representative of a pathogen in one aquifer system, it may not be the most appropriate surrogate in another porous media system. This was apparent in the inferior performance of MS2 as a surrogate, which is commonly used in virus transport studies.

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.

Functionalization of micromodels with kaolinite for investigation of low salinity oil-recovery processes

Sandstone formations are ubiquitous in both aquifers and petroleum reservoirs, of which clay is a major constituent. The release of clay particles from pore surfaces as a result of reduced injection fluid salinity can greatly modify the recovery of hydrocarbons from subsurface formations by shifting the wettability properties of the rock. In this paper we demonstrate a microfluidic approach whereby kaolinite is deposited into a two-dimensional microfluidic network (micromodel) to enable direct pore-scale, real-time visualization of fluid-solid interactions with representative pore-geometry and realistic surface interactions between the reservoir fluids and the formation rock. Structural characterization of deposited kaolinite particles agrees well with natural modes of occurrence in Berea sandstones; hence, the clay deposition method developed in this work is validated. Specifically, more than 90% of the deposited clay particles formed pore-lining structures and the remainder formed pore bridging structures. Further, regions of highly concentrated clay deposition likely leading to so-called Dalmatian wetting properties were found throughout the micromodel. Two post-deposition treatments are described whereby clay is adhered to the silicon surface reversibly and irreversibly resulting in microfluidic systems that are amenable to studies on (i) the fundamental mechanisms governing the increased oil recovery during low salinity waterfloods and (ii) the effect of a mixed-wet surface on oil recovery, respectively. The reversibly functionalized platform is used to determine the conditions at which stably adhered clay particles detach. Specifically, injection brine salinity below 6000 ppm of NaCl induced kaolinite particle release from the silicon surface. Furthermore, when applied to an aged system with crude oil, the low salinity waterflood recovered an additional 14% of the original oil in place compared to waterflooding with the formation brine.

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.

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

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

Guar gum solutions for improved delivery of iron particles in porous media (part 1): porous medium rheology and guar gum-induced clogging

The present work is the first part of a comprehensive study on the use of guar gum to improve delivery of microscale zero-valent iron particles in contaminated aquifers. Guar gum solutions exhibit peculiar shear thinning properties, with high viscosity in static conditions and lower viscosity in dynamic conditions: this is beneficial both for the storage of MZVI dispersions, and also for the injection in porous media. In the present paper, the processes associated with guar gum injection in porous media are studied performing single-step and multi-step filtration tests in sand-packed columns. The experimental results of single-step tests performed by injecting guar gum solutions prepared at several concentrations and applying different dissolution procedures evidenced that the presence of residual undissolved polymeric particles in the guar gum solution may have a relevant negative impact on the permeability of the porous medium, resulting in evident clogging. The most effective preparation procedure which minimizes the presence of residual particles is dissolution in warm water (60°C) followed by centrifugation (procedure T60C). The multi-step tests (i.e. injection of guar gum at constant concentration with a step increase of flow velocity), performed at three polymer concentrations (1.5, 3 and 4g/l) provided information on the rheological properties of guar gum solutions when flowing through a porous medium at variable discharge rates, which mimic the injection in radial geometry. An experimental protocol was defined for the rheological characterization of the fluids in porous media, and empirical relationships were derived for the quantification of rheological properties and clogging with variable injection rate. These relationships will be implemented in the second companion paper (Part II) in a radial transport model for the simulation of large-scale injection of MZVI-guar gum slurries.

Release of E. coli D21g with transients in water content

Transients in water content are well-known to mobilize microorganisms that are retained in the vadose zone. However, there is no consensus on the relative importance of drainage and imbibition events on microorganism release. To overcome this limitation, we have systematically studied the release of Escherichia coli D21g during cycles of drainage and imbibition under various solution chemistry and initial conditions. Results from these column studies revealed the influence of imbibition and drainage on D21g release. In particular, imbibition efficiently released cells from the air-water interface (AWI) that were initially retained under steady-state unsaturated conditions by expansion of water films and destruction of the AWI. Conversely, significant release and transport of cells during drainage only occurred below a critical water saturation (water film thickness). In this case, a fraction of the cells that were initially retained on the solid-water interface (SWI) partitioned into the mobile aqueous phase and the AWI as the receding water film thickness decreased during drainage. The efficiency of cell release from the SWI during drainage was much less than for the AWI during imbibition. Cycles of drainage and imbibition removed cells from the SWI and the AWI, respectively. However, the peak concentration and amount of cells that were released increased with the number of retained cells and the amount of drainage and imbibition, and decreased with the number of drainage and imbibition cycles. Release of cells during drainage and imbibition was found to be more pronounced in the presence of a weak secondary minimum when the ionic strength (IS) was 5 mM NaCl. Increases in the solution IS decreased the influence of water transients on release, especially during drainage. Complete recovery of the retained cells could be achieved using both IS reduction and cycles of drainage and imbibition, even when the cells were retained under favorable attachment conditions. In general, cell release was more pronounced with transients in water content than transients in IS when the IS ≥ 5 mM.

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

Colloid release from surfaces in response to ionic strength and flow perturbations has been mechanistically simulated. However, these models do not address the mechanism by which colloid attachment occurs, at least in the presence of bulk colloid-collector repulsion (unfavorable conditions), which is a prevalent environmental condition. We test whether a mechanistic model that predicts colloid attachment under unfavorable conditions also predicts colloid release in response to reduced ionic strength (IS) and increased fluid velocity (conditions thought prevalent for mobilization of environmental colloids). The model trades in mean-field colloid-collector interaction for discrete representation of surface heterogeneity, which accounts for a combination of attractive and repulsive interactions simultaneously, and results in an attached colloid population (in primary minimum contact with the surface) having a distribution of strengths of attraction. The model moderates equilibrium separation distance by inclusion of steric interactions. By using the same model parameters to quantitatively predict attachment under unfavorable conditions, simulated release of colloids (for all three sizes) from primary minimum attachment in response to perturbations qualitatively matched experimental results, demonstrating that both attachment and detachment were mechanistically simulated.

Altered transport of lindane caused by the retention of natural particles in saturated porous media

Attachment and straining of colloidal particles in porous media result in their reversible and irreversible retention. The retained particles may either increase the retention of hydrophobic pollutants by sorption onto the particles, or enhance pollutant transport when particles, loaded with the pollutants, are remobilized. The present study examines the effects of retained particles on the transport of the hydrophobic pesticide lindane (gamma-hexachlorocyclohexane) in saturated porous media. The lignite particles used have median diameters of about 3 μm, 1 μm, 0.8 μm, and 0.2 μm, respectively. Laboratory column experiments were analyzed by numerical modeling in order to identify and understand the processes involved in the transport of the particles and of lindane. Four scenarios were considered in which the solution containing lindane is injected either during or after the elution of the particles. The results show that lignite particles retained in a sandy porous medium alter the transport of the invading lindane. Particle retention was high in all scenarios and increased with increasing particle size. Remobilization of particles occurred due to a change in solution chemistry, and continuous particle detachment was observed over time. Numerical modeling of particle transport suggests that both reversible attachment and irreversible straining affected the transport of the particles. Lindane was retarded in all scenarios due to the strong particle retention in conjunction with the sorption of lindane onto the sand and onto retained particles, and the limited number of mobile particles carrying lindane. Moreover, it was found that intra-particle diffusion limited adsorption/desorption of lindane onto/from both limestone fragments of the sand and lignite particles. We assume that retention of lindane is reversible even though lindane recovery was incomplete over the duration of the experiments. The analysis of the effluent concentration suggests that retained particles loaded with lindane may become a secondary source of lindane. Models describing the transport of lindane fitted the experimental data very well and indicated the specific contribution of retained particles to the retardation of lindane. Since the properties of lignite also known as brown coal are similar to those of biochar, the results of the present study could be extended to the potential effects of biochar on lindane and other contaminants in soils, which would include both their retention and their enhanced transport. However, while the transport mechanisms of lindane are similar in water-unsaturated soils and saturated porous media considered here, the behavior of particles is more complex, requiring additional studies.

Particle-facilitated transport of lindane in water-saturated tropical lateritic porous media

The persistent insecticide lindane [(1α,2α,3β,4α,5α,6β)-1,2,3,4,5,6-hexachlorocyclohexane] is still in use in many tropical countries and remains a threat to soil and water quality. We studied the sorption and transport of lindane onto and through lateritic soils in both the absence and presence of lignite particles, onto which lindane may preferably sorb. We determined a linear distribution coefficient of lindane onto the soil matrix of 3.38 ± 0.16 L kg. Soil particles were not released from the porous medium on changing ionic strength, and also transport of lindane was not affected by changes in ionic strength. We fitted coupled transport models for lindane and the particles to the data, revealing that: (i) sorption kinetics of lindane onto the matrix is described best by intraparticle diffusion; (ii) 20% of the total porosity of the lateritic sample is intraparticle porosity; and (iii) only lignite particles with a median diameter <0.45 μm were not retained in the porous medium and thus facilitated the transport of lindane. We conclude that although lindane and similar pollutants may sorb on tropical lateritic porous media, their transport may be facilitated by particles with high organic-C content or dissolved organic C (DOC). This may be of relevance in farmlands and swamp groundwater systems where DOC, produced by leaching or slow biodegradation of surface organic matter, could cause rapid groundwater contamination by sorbing pollutants. Moreover, the results of this study can help to understand nanoparticle behavior in lateritic soils as the size of particles that facilitate lindane transport approaches the nanoparticle size range.

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

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

Modeling microorganism transport and survival in the subsurface

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

Heteroaggregation of microparticles with nanoparticles changes the chemical reversibility of the microparticles' attachment to planar surfaces

•

We theoretically examine detachment of homo- and heteroaggregates from primary minima.

•

Attached homoaggregates in primary minima are irreversible to reduction in ionic strength.

•

The attachment of microparticles coated by nanoparticles is chemically reversible.

•

We explain observed disaggregation of particles from primary minima in the literature.

•

We explain why viruses exhibit most conservative transport behavior in the environment.

This study theoretically investigated detachment of homoaggregates and heteroaggregates attached on the planar surfaces at primary minima during transients in solution chemistry. The homoaggregates were represented as small colloidal clusters with well-defined structures or as clusters generated by randomly packing spheres using Monte Carlo method. The heteroaggregates were modeled as microparticles coated with nanoparticles. Surface element integration technique was adopted to calculate Derjaguin–Landau–Verwey–Overbeek (DLVO) interaction energies for the homoaggregates and heteroaggregates at different ionic strengths. Results show that attached homoaggregates on the planar surface at primary minima are irreversible to reduction in solution ionic strength whether the primary spheres of the homoaggregates are nano- or micro-sized. Heteroaggregation of nanoparticles with a microparticle can cause DLVO interaction energy to decrease monotonically with separation distance at low ionic strengths (e.g., ⩽0.01 M), indicating that the heteroaggregates experience repulsive forces at all separation distances. Therefore, attachment of the heteroaggregates at primary minima can be detached upon reduction in ionic strength. Additionally, we showed that the adhesive forces and torques that the aforementioned heteroaggregates experience can be significantly smaller than those experienced by the microspheres without attaching nanoparticles, thus, the heteroaggregates are readily detached via hydrodynamic drag. Results of study provide plausible explanation for the observations in the literature that attached/aggregated particles can be detached/redispersed from primary minima upon reduction in ionic strength, which challenges the common belief that attachment/aggregation of particles in primary minima is chemically irreversible.

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

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

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

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

Evaluating the effects of variable water chemistry on bacterial transport during infiltration

Bacterial infiltration through the subsurface has been studied experimentally under different conditions of interest and is dependent on a variety of physical, chemical and biological factors. However, most bacterial transport studies fail to adequately represent the complex processes occurring in natural systems. Bacteria are frequently detected in stormwater runoff, and may present risk of microbial contamination during stormwater recharge into groundwater. Mixing of stormwater runoff with groundwater during infiltration results in changes in local solution chemistry, which may lead to changes in both bacterial and collector surface properties and subsequent bacterial attachment rates. This study focuses on quantifying changes in bacterial transport behavior under variable solution chemistry, and on comparing the influences of chemical variability and physical variability on bacterial attachment rates. Bacterial attachment rate at the soil-water interface was predicted analytically using a combined rate equation, which varies temporally and spatially with respect to changes in solution chemistry. Two-phase Monte Carlo analysis was conducted and an overall input-output correlation coefficient was calculated to quantitatively describe the importance of physiochemical variation on the estimates of attachment rate. Among physical variables, soil particle size has the highest correlation coefficient, followed by porosity of the soil media, bacterial size and flow velocity. Among chemical variables, ionic strength has the highest correlation coefficient. A semi-reactive microbial transport model was developed within HP1 (HYDRUS1D-PHREEQC) and applied to column transport experiments with constant and variable solution chemistries. Bacterial attachment rates varied from 9.10×10(-3)min(-1) to 3.71×10(-3)min(-1) due to mixing of synthetic stormwater (SSW) with artificial groundwater (AGW), while bacterial attachment remained constant at 9.10×10(-3)min(-1) in a constant solution chemistry (AGW only). The model matched observed bacterial breakthrough curves well. Although limitations exist in the application of a semi-reactive microbial transport model, this method represents one step towards a more realistic model of bacterial transport in complex microbial-water-soil systems.

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

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

Transport and retention of high concentrated nano-Fe/Cu particles through highly flow-rated packed sand column

The design of an efficient field-scale remediation based on the use of nanoscale zero valent iron (NZVI) requires an accurate assessment of the mobility of such particles in saturated porous media, both during injection in the subsurface (short-term mobility) and later (long-term mobility). In this study, the mobility of highly concentrated dispersions of bimetallic Fe/Cu nanoparticles (d(50) = 70 ± 5 nm) in sand-packed columns (0.5 m length and 0.025 m inner diameter) was studied. In particular, the influence of flow rate (V = 5 × 10(-4), 1 × 10(-3), 2 × 10(-3) m/s) and injected particle concentrations (2, 5, 8, 12 g/l) was addressed. Breakthrough curves and water pressure drop along the column, averaged effective porosity and final distribution of retained particles along the column were measured. Experimental results evidenced a good mobility of the Fe/Cu particles, with significant breakthrough in all explored experimental conditions of flow rate and C(0), without requiring the addition of any stabilizing agent. Clogging phenomenon of the column and also the pore pressure variation during injection period are strongly affected by injected concentration. Clogging due to deposition of particles following a ripening dynamics was observed in particular for C(0) = 8 and 12 g/l. The experimental data were modeled using the E-MNM1D software. The study has implications for field injection of bimetallic nanoparticles, suggesting that particular care is to be devoted when selecting injection concentration, to avoid porous medium clogging and control the radius of influence.