Coupled effects of hydrodynamic and solution chemistry on long-term nanoparticle transport and deposition in saturated porous media

Enhanced sinks of polystyrene nanoplastics (PSNPs) in marine sediment compared to freshwater sediment: Influencing factors and mechanisms

The difference in the transport behaviors of nanoplastics consistently assistant with their toxicities to benthic and other aquatic organisms is still unclear between freshwater and marine sediments. Here, the mobilities of polystyrene nanoplastics (PSNPs) and key environmental factors including salinity and humic acid (HA) were systematically studied. In the sand column experiments, both tested PSNPs in the freshwater system (100 nm NPs (100NPs): 90.15 %; 500 nm NPs (500NPs): 54.22 %) presented much higher penetration ratio than in the marine system (100NPs: 8.09 %; 500NPs: 19.04 %). The addition of marine sediment with a smaller median grain diameter caused a much more apparent decline in NPs mobility (100NPs: from 8.09 % to 1.85 %; 500NPs: from 19.04 % to 3.51 %) than that containing freshwater sediment (100NPs: from 90.15 % to 83.56 %; 500NPs: from 54.22 % to 41.63 %). Interestingly, adding HA obviously led to decreased and slightly increased mobilities for NPs in freshwater systems, but dramatically improved performance for NPs in marine systems. Electrostatic and steric repulsions, corresponding to alteration of zeta potential and hydrodynamic diameter of NPs and sands, as well as minerals owing to adsorption of dissolved organic matter (DOM) and aggregations from varied salinity, are responsible for the mobility difference.

Statuses, shortcomings, and outlooks in studying the fate of nanoplastics and engineered nanoparticles in porous media respectively and borrowable sections from engineered nanoparticles for nanoplastics

This review discussed the research statuses, shortcomings, and outlooks for the fate of nanoplastics (NPs) and engineered nanoparticles (ENPs) in porous media and borrowable sections from ENPs for NPs. Firstly, the most important section was that we reviewed the research statuses on the fate of NPs in porous media and the main influencing factors, and explained the influencing mechanisms. Secondly, in order to give NPs a reference of research ideas and influence mechanisms, we also reviewed the research statuses on the fate of ENPs in porous media and the factors and mechanisms influencing the fate. The main mechanisms affecting the transport of ENPs were summarized (Retention or transport modes: advection, diffusion, dispersion, deposition, adsorption, blocking, ripening, and straining; Main forces and actions: Brownian motion, gravity, electrostatic forces, van der Waals forces, hydration, filtration, bridging; Affecting elements of the forces and actions: the ENP and media grain surface functional groups, size, shape, zeta potential, density, hydrophobicity, and roughness). Instead of using the findings of ENPs, thorough study on NPs was required because NPs and ENPs differed greatly. Based on the limited existing studies on the NP transport in porous media, we found that although the conclusions of ENPs could not be applied to NPs, most of the influencing mechanisms summarized from ENPs were applicable to NPs. Combining the research thoughts of ENPs, the research statuses of NPs, and some of our experiences and reflections, we reviewed the shortcomings of the current studies on the NP fate in porous media as well as the outlooks of future research. This review is very meaningful for clarifying the research statuses and influence mechanisms for the NP fate in porous media, as well as providing a great deal of inspiration for future research directions about the NP fate in porous media.

Particle dispersion through porous media with heterogeneous attractions

Porous media used in many practical applications contain natural spatial variations in composition and surface charge that lead to heterogeneous physicochemical attractions between the media and transported particles. We performed Stokesian dynamics (SD) simulations to examine the effects of heterogeneous attractions on quiescent diffusion and hydrodynamic dispersion of particles within geometrically ordered arrays of nanoposts. We find that transport under quiescent conditions occurs by two mechanisms, diffusion through the void space and intermittent hopping between the attractive wells of different nanoposts. As the attraction heterogeneity increases, the latter mechanism becomes dominant, resulting in an increase in the particle trajectory tortuosity, deviations from Gaussian behavior in the particle displacement distributions, and a decrease in the long-time particle diffusivity. Similarly, under flow conditions corresponding to low Péclet number (Pe), increased attraction heterogeneity leads to transient localization near the nanoposts, resulting in a broadening of the particle distribution and enhanced longitudinal dispersion in the direction of flow. At high Pe where advection strongly dominates, however, the longitudinal dispersion coefficient is insensitive to attraction heterogeneity and exhibits Taylor-Aris dispersion behavior. Our findings provide insight into how heterogeneous interactions may influence particle transport in complex 3-D porous media.

Transformation of organic carbon through medium pressure (polychromatic) UV disinfection of wastewater effluent during wet weather events

An observed decrease in total organic carbon (TOC) and dissolved organic carbon (DOC) concentrations following wastewater disinfection with medium pressure (MP, polychromatic) ultraviolet (UV) irradiation during wet weather flows is investigated. When antecedent rainfall in the previous 7-days was >2 in (5 cm), TOC and DOC concentrations decreased dramatically following MP-UV disinfection. Organic carbon surrogate measurements of biological oxygen demand (BOD), TOC, DOC, turbidity, UVA - 254 nm, SUVA (specific UVA), scanning UV-Visible spectra (200-600 nm), fluorescence excitation-emission matrix (EEM) spectra, and light scattering data are presented for wastewater resource recovery facility (WRRF) influent, secondary effluent (pre-UV-disinfection), and MP-UV-disinfected (final effluent) samples. TOC and DOC in wastewater influent and secondary effluent (i.e., pre-UV disinfection) correlated with antecedent rainfall conditions. The percent TOC and DOC removal through secondary treatment (i.e., from influent to effluent pre-UV) and the percent TOC and DOC removal through MP-UV disinfection (i.e., from effluent pre-UV to effluent post-UV) were compared and the latter approached 90 % through MP-UV disinfection during high antecedent rainfall conditions. Spectroscopy (UV, visible, or fluorescence) was performed on samples after filtration through 0.45 μm filters, i.e., the operationally defined DOC fraction of aquatic carbon. Scanning UV-visible spectra indicated transformation of an unidentified wastewater component into light-scattering entities regardless of antecedent rainfall conditions. The types of organic carbon (diagenetic, biogenic, or anthropogenic) and the significance of wet weather are discussed. An organic carbon contribution via infiltration and inflow was attributed as a source-of-interest in this research.

Prediction of transport behavior of nanoparticles using machine learning algorithm: Physical significance of important features

There is a rising concern related to the possible risk of human exposure to nanoparticles (NPs). Several studies have reported on the transport behavior of NPs in the porous media under varying conditions. Thus, there is a scope to use this information in a predictive model so that the transport behavior of any un-explored NPs could be predicted. The main focus of his study, therefore, is to apply different machine learning (ML) based models to predict the transport efficiency of a wide range of NPs and to identify the important features. To achieve the objective, first, the dataset is prepared by extracting data from published papers for selected NPs [i.e., silver (nAg), titanium dioxide (nTiO2), zinc oxide (nZnO), graphene oxide (nGO), and etc.]. Then, random forest, XGBoost, and CatBoost algorithms combined with synthetic minority oversampling technique (SMOTE) were applied where retention fraction (RF) is considered as the target feature and particle characteristics (i.e., surface charge, size, concentration), solution chemistry [pH, ionic strength (IS]), porous media properties (grain size, porosity) and flow rate are considered as the training features. The outcome of the study indicates that CatBoost combined with SMOTE performed the best in predicting RF for the entire range of NPs (R2 > 0.89 and MSE < 0.007) as well as for individual NPs. Feature importance analysis indicates four features, namely zeta potential, IS, pH, and particle diameter (the entire range of NPs, nGO, nZnO) or grain size (nAg, nTiO2) have significant weightage (>75%). The result suggests that the features overrule the prediction of transport behavior rather than the types of individual NPs. The relative importance of the features depends on the range of the parameter used. The identified important features are in accordance with the underlying physical process, which makes the prediction model more reliable.

Cotransport of microplastics and sulfanilamide antibiotics in groundwater: The impact of MP/SA ratio and aquifer media

The aim of this study was to investigate the effects of the aquifer media, structure type, and initial concentration ratio of contaminants on the cotransport behavior of microplastics (MPs) and sulfanilamide antibiotics (SAs) through a series of one-dimensional column experiments in groundwater. Under a single suspension system, the relative mass recovery rates of fine sand, medium sand, and coarse sand were 25.65%, 37.50%, and 57.91%, respectively. The breakthrough curve of MPs showed a weak and slow upward trend, indicating that the migration of MPs in aqueous media is mainly blocked by the surface. The migration results of different structure type on SAs (ST, SM, SM2, SMX) in a single suspension system indicated that the deposition rate coefficients (k_c) of the four SAs were 1.23 × 10^-1, 9.09 × 10^-2, 1.11 × 10^-1, and 8.87 × 10^-2. Under a binary suspension system (MPs:ST = 1:1), the maximum effluent concentration (MEC) of MPs in fine sand, medium sand, and coarse sand increased to 0.52, 0.64, and 0.88, respectively, and the relative mass recovery rates of ST were 22.79%, 23.59%, 20.25%. This results show that the coexistence of MPs and SAs significantly promotes the migration of MPs and inhibits that of SAs. It is mainly because of their carrier action, adsorption sites and additional deposit sites for MPs through SAs pre-deposition on media. When the initial concentration ratio was 2:1, the particles had the highest Zeta potential (-48.3 mV) and the highest potential barrier (3200 kBT), leading to the formation of complex aggregates (MPs-SAs-MPs) owing to the aggregation of colloidal MPs. The increase in the volume and number of MPs-SAs co-aggregates on the surface of the media as the initial concentration of MPs increases, which was mainly due to the disappearance of surface blocking effect and the occurrence of filtering maturation effect.

Transport of degradable/nondegradable and aged microplastics in porous media: Effects of physicochemical factors

The degradable properties of degradable plastics allow them to form microplastics (MPs) faster. Therefore, degradable MPs may easily be transported in the underground environment. Research on degradable MPs transport in porous media is necessary and urgent. In this study, polylactic acid (PLA) and polyvinyl chloride (PVC) were selected to compare the transport differences between degradable and nondegradable MPs under different factors (flow rates, ionic strengths (ISs), pH, and coexisting cations) through column experiments, and UV irradiation was used to further simulate the effect of aging on different types of MPs. Fourier transform infrared (FT-IR) and X-ray photoelectron spectroscopy (XPS) were used to characterize functional groups and to determine the surface elements of MPs, respectively. The results showed that MPs were more mobile at higher flow rate, lower IS, higher pH, and monovalent cations. The order of transport capacity of MPs was PVC < aged PVC < PLA < aged PLA. This result was mainly attributed to the more negative Zeta potential and higher dispersion stability of aged PLA and PLA, which were caused by abundant O-functional groups. Compared with PVC, the O/C ratio of PLA increased significantly after aging, indicating that PLA was more prone to aging. The advection-dispersion-equation (ADE) fitted the transport data of MPs well. The interaction energy of MPs and quartz sand was accurately predicted by the Derjaguin-Landau-Verwey-Overbeek (DLVO) theory. This work contributes to a comprehensive understanding of the transport of degradable MPs in the environment.

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.

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.

Application of temporal moment analysis to interpret colloid and colloid-facilitated solute transport under varying size exclusion and attachment coefficient

Colloidal particles can attach the contaminants like heavy metals and radionuclides and act as contaminant carriers to provide a faster movement of pollutants through the interconnecting pores of the porous medium. The electrostatic repulsive force between the negatively charged colloids and the solid surface restricts the transport of the mobile colloids to the larger pores of the porous medium and initiates the size exclusion mechanism. The temporal moment analysis is treated as an effective tool to interpret the solute breakthrough curves for analyzing the statistical behavior of the contaminants. In past literature, the temporal moments have not been incorporated with the breakthrough curves of colloids and colloid-facilitated contaminants for statistical interpretation. In this research study, the temporal variations of concentrations of mobile colloids, solute attached to the mobile colloids, and the dissolved solute are obtained numerically in a fully saturated one-dimensional column considering a continuous source for varying size exclusion and colloid attachment coefficient. Utilizing the simulated spatially varying breakthrough curves, the temporal moments are estimated to calculate the mass recovery, average residence time, and the spreading of mobile colloids and dissolved solutes. The temporal moment analysis suggests that the velocity enhancement for higher size exclusion reduces the average residence time of the mobile colloids and the solute adsorbed to the mobile colloids significantly. The mass recovery of mobile colloids and the solute attached to the mobile colloids increases at a specific depth for higher size exclusion. The estimated second central moment attributes that the solute spreading follows the nonlinear trend for low size exclusion. The peaks of the relative concentration of mobile colloids and solute attached to mobile colloids drastically decrease with an increase in attachment coefficient. The peak of the relative concentration of dissolved contaminant enhances with attachment coefficient. The high second temporal moment of the dissolved contaminant at a higher attachment coefficient indicates the slow interaction of dissolved solute and porous medium and that enables a greater spreading of solute through the interconnecting porous medium. The study suggests that the faster movement of mobile colloids and the solute attached to the mobile colloids at higher exclusion imparts a potential risk of groundwater contamination and thorough statistical interpretation is needful to analyze the behavior of colloids and colloid-facilitated contaminants. The research work does not consider the transient flow field and the effect of the presence of air phase in the partially saturated soil column in the groundwater system.

Hydrobiological Mechanism Controlling the Synergistic Effects of Unsaturated Flow and Soil Organic Matter on the Degradation of Emerging Organic Contaminants in Soils

Hydrology is a key factor influencing microbial degradation of emerging organic contaminants (EOCs) in soils, but the underlying mechanisms are not clear. In this study, biotic and abiotic column experiments were performed to investigate the removal and degradation of five EOCs in soils with different soil organic matter (SOM) contents under saturated and unsaturated flow conditions. In biotic experiments, 54-90% of bisphenol A (BPA) and 9-22% of ibuprofen (IBU) were removed from the aqueous phase of saturated columns due to adsorption and biodegradation. The biodegradation removed 26-65% of BPA and 1-22% of IBU. Decreasing soil pore water saturation from 100 to 80% increased BPA removal to 97-100% and IBU removal to 42-43% due to increased biodegradation (67-81% for BPA and 36-39% for IBU). No significant removal of BPA and IBU was observed in SOM-removed soils under saturated and unsaturated flow conditions. The desaturation did not influence sorptive losses of BPA (<27%) and IBU (<7%), suggesting their negligible adsorption at air-water interfaces but increased biodegradation of BPA and IBU sorbed at SOM-water interfaces. The study shows that soil drying and SOM can synergistically degrade BPA and IBU but have no effect on recalcitrant carbamazepine, tetracycline, and ciprofloxacin.

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.

Investigating transport kinetics of polystyrene nanoplastics in saturated porous media

Understanding the fate and transport of polystyrene nanoparticles (PSNPs) in porous media under various conditions is necessary for evaluating and predicting environmental risks caused by microplastics. The transport kinetics of PSNPs are investigated by column experiment and numerical model. The surface of DLVO interaction energy is calculated to analyze and predict the adsorption and aggregation of PSNPs in porous media, which the critical ionic strength of PSNPs can be accurately investigated. The results of the DLVO energy surface suggest that when the concentration of Na⁺ increases from 1 mM to 50 mM, the DLVO energy barrier of PSNPs-silica sand (SS) decreases from 78.37 kT to 5.46 kT. As a result, PSNPs are easily adsorbed on the surface of SS and the mobility of PSNPs is reduced under the condition of a high concentration of Na⁺ (PSNPs recovery rate decreases from 62.16% to 3.65%). When the concentration of Ca²⁺ increases from 0.1 mM to 5 mM, the DLVO energy barrier of PSNPs-SS decreases from 12.10 kT to 1.90 kT, and PSNPs recovery rate decreases from 82.46% to 4.27%. Experimental and model results showed that PSNPs mobility is enhanced by increasing initial concentration, flow velocity and grain size of SS, while the mobility of PSNPs with larger particle diameter is lower. Regression analysis suggests that kinetic parameters related to PSNPs mobility are correlated with DLVO energy barriers. The environmental behavior and mechanism of PSNPs transport in porous media are further investigated in this study, which provides a scientific basis for the systematic and comprehensive evaluation of the environmental risk and ecological safety of nano-plastic particles in the groundwater system.

Fate and removal efficiency of polystyrene nanoplastics in a pilot drinking water treatment plant

Occurrence of microplastics and nanoplastics in aquatic systems, as well as in water compartments used to produce drinking water have become a major concern due to their impact on the environment and public health. Nanoplastics in particular, in regard to their fate and removal efficiency in drinking water treatment plants (DWTP), which ensure water quality and supply drinking water for human consumption have been, by far, rarely investigated. This study investigates the removal efficiency of polystyrene (PS) nanoplastics in a conventional water treatment plant providing drinking water for 500'000 consumers. For that purpose, a pilot-scale DWTP, located within the main treatment plant station, reproducing at a reduced scale the different processes and conditions of the main treatment plant is used. The results show that filtration process through sand and granular activated carbon (GAC) filters in the absence of coagulation achieves an overall nanoplastic removal of 88.1%. The removal efficiency of filtration processes is mainly attributed to physical retention and adsorption mechanisms. On the other hand, it is found that coagulation process greatly improves the removal efficiency of nanoplastics with a global removal efficiency equal to 99.4%. The effective removal efficiency of sand filtration increases considerably from 54.3% to 99.2% in the presence of coagulant, indicating that most of PS nanoplastics are removed during sand filtration process. The higher removal efficiency with the addition of coagulant is related to nanoplastics surface charge reduction and aggregation thus significantly increasing their retention in the filter media.

Transport of nanoparticulate TiO2 UV-filters through a saturated sand column at environmentally relevant concentrations

The fate of sunscreen residues released during bathing activities around recreational areas is an emerging concern regarding the potential ecotoxicity of some of their ingredients, including nanoparticulate TiO₂ UV-filters. To assess the extent of contamination in the natural medium, sand-packed column experiments were carried out with bare TiO₂ engineered nanoparticles (ENPs) and two commercial nano-TiO₂ UV-filters coated with either SiO₂ (hydrophilic) or a combination of Al₂O₃ and simethicone (amphiphilic). The high sensitivity of (single particle)ICPMS online monitoring of the breakthrough curves made it possible to inject the ENPs at trace levels (2-100 μg L^-1) in eluents composed of 10^-3 and 10^-2 M NaCl and pHs of 5.7 and 7.8. The deposition of all ENPs in the sand increased with the ionic strength and decreased with the pH of the carrier. Both bare and SiO₂-coated ENPs showed a clear control by the electrostatic interactions between the particles and the quartz grains surfaces, in partial agreement with classical DLVO theory. Unexpectedly high rates of transfer were observed with the amphiphilic UV-filter, which could be explained by the addition of a contribution to the DLVO model to account for the steric repulsion between the sand collector and the polysiloxane surface layer of this ENP. These results demonstrate the major role played by the coating of UV-filters regarding their fate in porous media like soils, sediments and aquifers.

Effects of physicochemical factors on the transport of aged polystyrene nanoparticles in saturated porous media

Plastic debris, particularly nanoplastics, have attracted substantial attention as an emerging pollutant of global concern. The aging process caused by UV could dramatically alter the physicochemical properties of polystyrene plastics and affect their fate and transport in the subsurface environment. We researched the effects of diverse factors, including flow rate, ionic strength (IS), and cation valence on the transport of aged polystyrene nanoparticles (PSNPs) in saturated porous media and found that their mobility was greatly enhanced by the aging process at all other experimental conditions except coexistence with Al³⁺. Interestingly, we found that the aged PSNPs were polymerized due to the coexistence with Al³⁺, which led the aged PSNPs to exhibit weaker mobility than the pristine. Zeta potential and particle size measurements, FTIR spectroscopy, and XPS were employed to characterize the PSNPs, and the results suggest that UV radiation provides O-containing functional groups for the PSNPs. The experimental results correspond well with the ADR model and the values of Smax and k closely reflect the deposition of the PSNPs in sand columns. Moreover, the Derjaguin-Landau-Verwey-Overbeek (DLVO) theoretical calculation accurately reflects the interaction of the pristine and aged PSNPs and quartz sand. Overall, due to the processes that PSNPs possibly undergo in the environment, their mobility may fluctuate dramatically. These findings help remedy knowledge deficiency regarding nanoplastic mobility being affected by aging processes, further underscore the critical influence of the aging process on the potential risks and environmental fates of nanoplastics.

Band application of flue gas desulfurization gypsum improves sodic soil amelioration

Blending flue gas desulfurization (FGD) gypsum with surface sodic soil is a universally recognized method for the rapid amelioration of sodic soils; however, little information is available on whether other application methods (band application) will reclaim sodic soil. Three FGD gypsum application methods (single-band, dual-band and blend applications) and a control treatment (non-FGD gypsum) were carried out using sodic soil in soil bins to investigate the effects of the application method on the wetting front, major cations in the leachate during the process of water infiltration and soluble and exchangeable cations in the soil profile after infiltration. The results showed that the wetting fronts in the band treatments were denser in the horizontal direction than in the vertical direction, but the blend and control treatments only had vertical migration. The main channel of the stream in the band treatment was concentrated below the application site of FGD gypsum. The orders of desalting capacity were blend treatment, dual-band treatment and single-band treatment for the same volume of outlet water. There was no water outflow in the control treatment even after 115 days of leaching. The dual-band treatment significantly decreased the soil sodicity of the 0-40 cm soil profile, while the single-band treatment only effectively reclaimed (horizontally) half of the soil. In the blend treatment, the exchangeable sodium percentages were 21.3 % and 34.7 % at depths of 30-35 cm and 35-40 cm, respectively, and were close to zero at a depth of 0-30 cm. Compared with blend treatment, band application could be a better way to reclaim sodic soil with FGD gypsum due to its advantages of long-term and efficient amelioration with low consumption.

Mechanisms of bentonite colloid aggregation, retention, and release in saturated porous media: Role of counter ions and humic acid

In the subsurface environment, colloids play an important role in pollutant transport by acting as the carriers. Understanding colloid release, transport, and deposition in porous media is a prerequisite for evaluating the potential role of colloids in subsurface contaminant transport. In this work, the aggregation, retention, and release of bentonite colloid in saturated porous sand media were investigated by kinetic aggregation and column experiments, the correlation and mechanism of these processes were revealed by combining colloid filtration theory, interaction energy calculation and density functional theory. The results showed that the retention and release of colloids were closely related to the dispersion stability and filtration effect. Multivalent cations with higher mineral affinity reduced the colloid stability, and the dispersion stability and mobility of the colloid were greatly improved by humic acid due to the enhancement of electrostatic repulsion and steric hindrance effects. The primary minimum interaction was found to contribute more to irreversible colloid retention in a Ca²⁺ system, while the secondary energy minimum was found to be responsible for colloid release with the occurrence of transient solution chemistry. The deposited colloid aggregates could be redistributed and released when the solution chemistry became favorable towards dispersion. These findings provide essential insight into the environmental colloid fate as well as a vital reference for the risk of colloid-driven transport of contaminants in the subsurface aquifer environment.

Transport of silica encapsulated DNA microparticles in controlled instantaneous injection open channel experiments

Surface water tracing is a widely used technique to investigate in-stream mass transport including contaminant migration. Recently, a microparticle tracer was developed with unique synthetic DNA encapsulated in an environmentally-friendly silica coating (Si-DNA microparticle). Previous tracing applications of such tracers reported detection and quantification, but a massive loss of tracer mass. However, the transport behavior of these DNA-tagged microparticle tracers has not been rigorously quantified and compared with that of solute tracers. Therefore, we compared the transport behavior of Si-DNA microparticles to the behavior of solute NaCl in 6 different, environmentally representative water types using breakthrough curves (BTCs), obtained from laboratory open channel injection experiments, whereby no Si-DNA microparticle tracer mass was lost. Hereafter, we modelled the BTCs using a 1-D advection-dispersion model with one transient storage zone (OTIS) by calibrating the hydrodynamic dispersion coefficient D and a storage zone exchange rate coefficient. We concluded that the transport behavior of Si-DNA microparticles resembled that of NaCl in surface-water relevant conditions, evidenced by BTCs with a similar range of D; however, the Si-DNA microparticle had a more erratic BTC than its solute counterpart, whereby the scatter increased as a function of water quality complexity. The overall larger confidence interval of D_Si-DNA was attributed to the discrete nature of colloidal particles with a certain particle size distribution and possibly minor shear-induced aggregations. This research established a solid methodological foundation for field application of Si-DNA microparticles in surface water tracing, providing insight in transport behavior of equivalent sized and mass particles in rivers.

Porous media transport of iron nanoparticles for site remediation application: A review of lab scale column study, transport modelling and field-scale application

Injection of surface modified zero valent iron nanoparticles for in situ remediation of soil, contaminated with an array of pollutants has attracted great attention due to the high reactivity of zero valent iron towards a broad range of contaminants, its cost effectiveness, minimal physical disruption and low toxicity. The effectiveness of this technology relies on the stability and mobility of injected iron nanoparticles. Hence the development of a modelling tool capable of predicting nZVI transport is indispensable. This review provides state of the art knowledge on the mobility of iron nanoparticles in porous media, mechanisms involved in subsurface retention of nZVI based on continuum models and field scale application. Special attention is given to the identification of the influential parameters controlling the transport potential of iron nanoparticles and the available numerical models for the simulation of laboratory scale transport data.

Cotransport and deposition of colloidal polystyrene microplastic particles and tetracycline in porous media: The impact of ionic strength and cationic types

The cotransport behaviors of colloidal polystyrene microplastic particles (PSMPs) and tetracycline (TC) (20 mg/L) were investigated in saturated porous media in KCl and CaCl₂ solutions of various ionic strengths (1, 10, 50, 100 mM). Furthermore, the effects of TC concentration (0, 1, 5, 10, 20 mg/L) on the cotransport behaviors of PSMPs and TC in 100 mM KCl solution were assessed. The cotransport behaviors were analyzed by comparing the individual transport behaviors of PSMPs or TC. When cotransported, the presence of TC (20 mg/L) slightly inhibited PSMPs mobility in K⁺ solutions (the C/C₀ decreased in the range of 0-5.9%), but facilitated it in Ca²⁺ solutions (the C/C₀ increased in the range of 6.7-42.6%). In KCl solutions, although the presence of TC (PSMPs) did not significantly affect the transport behaviors of PSMPs (TC), the attachment efficiencies of both PSMPs and TC showed a non-linear and non-monotonic change with increase in ionic strength. However, in CaCl₂ solutions, the effects of TC (PSMPs) on the transport behaviors of PSMPs (TC) were remarkable and a non-linear non-monotonic change was observed. The adsorption of TC on PSMPs might play a critical role during the cotransport. Thus, the balance between the transport-inhibiting (e.g., the reduction in electrostatic repulsive force) and transport-facilitating effects (e.g., the effects on hydrophilicity/hydrophobicity of PSMPs due to TC adsorption) may be responsible for the observed changes. Overall, the results demonstrated that the cotransport behaviors of PSMPs and TC were more complicated than their individual transport behaviors in porous media, which might vary considerably with environmental conditions. This work could greatly improve our understanding of complex cotransport behaviors and environmental risk of PSMPs.

Transport of Tl(I) in water-saturated porous media: Role of carbonate, phosphate and macromolecular organic matter

Understanding the transport behaviors of thallium (Tl) in porous media is of considerable interest for both natural soils and artificial filtration removal of Tl. In this context, the transport behaviors of Tl(I) in water-saturated sand columns under different conditions were systematically investigated. It was found that, in addition to the effects of pH and ionic strength (IS), the transport of Tl(I) depended on the carbonate, phosphate and macromolecular organic matter as well. Tl(I) broken the columns more difficultly under higher pH and lower IS conditions. Moreover, the adsorption of carbonate and phosphate on sand surfaces may increase the retention of Tl(I) in columns. As for macromolecular organic matter, humic acid (HA) facilitated Tl(I) transport, especially under neutral and alkaline conditions (7.0 and 9.8), which was possibly associated with Tl-complexes formation and competed adsorption between Tl(I) and HA. However, bovine serum albumin (BSA) impeded Tl(I) transport for the reason that deposited BSA might provide more adsorption sites for Tl(I), though Tl(I) had a slight effect on BSA transport. In order to evaluate the mechanisms of transport, a dual-sites non-equilibrium model was applied to fit the breakthrough curves of Tl(I). Retardation factor (R) values of individual Tl(I) transport from model calculations were found to be higher than that of Tl(I) transport with HA and lower than that of Tl(I) transport with BSA. The fraction of instantaneous sorption sites (β) was found to decrease with increasing pH, implying nonequilibrium sorption is a main sorption mechanism of Tl(I) with pH increasing. The fundamental data obtained herein demonstrated that carbonate, phosphate and macromolecular organic matter significantly influenced the Tl(I) migration and could lead to the leaking or bindings of Tl(I) at Tl-occurring sites.

Transport of polystyrene nanoplastics in natural soils: Effect of soil properties, ionic strength and cation type

Nanoplastics as emerging pollutants have caused growing concerns and posed potential threats to the environment. Nonetheless, only few studies investigated transport behaviors of nanoplastics in natural soils. In this study, column experiments were conducted to investigate the effect of soil properties, ionic strength and cation type on the transport of polystyrene nanoplastics (PSNPs) in a desert soil (DS), a black soil (BS) and a red soil (RS). The effluent recovery of PSNPs in three soils followed the order of DS (0%-96.8%) > BS (0%-87.5%) > RS (0%). The retention of PSNPs was positively correlated with Fe/Al oxides contents (DS: Fe-2.69%, Al-12.6%; BS: Fe-4.04%, Al-15.9%; RS: Fe-6.57%, Al-26.9%), whereas negatively correlated with soil pH (DS: 9.75; BS: 6.57; RS: 4.97). Soil minerals and pH were thus identified as the crucial soil properties determining transport of PSNPs, due to their coupled effects on surface charges to affect electrostatic interactions between soils and PSNPs. In addition, increasing solution ionic strength strongly inhibited the transport of PSNPs in the DS (0%-96.8%) and BS (0%-87.5%). Ca²⁺ (IS: 1-5 mM) was more pronounced in enhancing PSNP retention than Na⁺ (IS: 1-20 mM). Our findings highlight that the transport and fate of PSNPs in natural soils are highly sensitive to soil physicochemical properties, ionic strength and cation type, and reveal that nanoplastics have strong mobility ability in soils with high pH and low Fe/Al oxides contents, which may pose potential risks to the soil and groundwater environment.

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

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

Insight into the Crucial Role of Secondary Mineral Phases in the Transfer of Gold Nanoparticles through a Sand Column Using Online ICP-MS/spICP-MS Monitoring

Investigating the transport of engineered nanoparticles through representative soils is an important issue in assessing their mobility and fate in the environment. In this study, successive injections of gold nanoparticles (AuNPs) were performed in a quartz sand column with an eluent composed of 10^-2 M NaCl at a pH of 7.5. After this series of injections, remobilization of the AuNPs was examined by raising the eluent pH to 10. ¹⁹⁷Au and the conservative ionic tracer ⁷⁹Br were monitored simultaneously by online inductively coupled plasma mass spectrometry (ICP-MS), and the particulate nature of gold eluting from the column was confirmed by setting the ICP-MS in the "single particle" mode. The extent of AuNP attachment was greater than predicted by DLVO theory considering quartz as the sole collector, decreased with the number of injections and with particle size. In contrast with the repulsive interaction energy between the particles and the quartz surface, kaolinite, a secondary mineral of the sand, provided favorable conditions for particle attachment. The superimposed signals of ¹⁹⁷Au and ²⁷Al in the column effluent after pH increase suggest that gold nanoparticles were essentially remobilized as heteroaggregates with the kaolinite colloids they were attached to when favorable conditions for clay detachment from the sand grains were encountered.

Mechanism of Suspended Kaolinite Particle Clogging in Porous Media During Managed Aquifer Recharge

Managed aquifer recharge is an effective strategy for urban stormwater management. Chemical ions are normally retained in stormwater and groundwater and may accelerate clogging during the recharge process. However, the effect of water chemistry on physical clogging has not previously been investigated. In this study, we investigated the hydrogeochemical mechanism of saturated porous media clogging in a series of column experiments. The column was packed with river sand and added suspensions of kaolinite particles. Calcium chloride and sodium chloride are used as representative ions to study chemical effects. We found that an increase in ionic strength resulted in retention of kaolinite solids in the column, with a breakthrough peak of C/C₀ value of 1 to 0.2. The corresponding hydraulic conductivity decreased with increased solids clogging. Divalent cations were also found to have a greater influence on kaolinite particle clogging than monovalent cations. The enhanced hydrochemical-related clogging was caused by kaolinite solids flocculating and increasing the deposition rate coefficient by 1 to 2 times in high ionic strength conditions. Three clogging mechanisms of kaolinite solids are proposed: surface filtration, inner blocking, and attachment. This study further deepens the understanding of the mechanisms of solids clogging during aquifer recharge and demonstrates the significance of ionic strength on recharge clogging risk assessments.

Nanoaggregates of silica with kaolinite and montmorillonite: Sedimentation and transport

Due to a wide range of applications in industrial fields, engineered nanomaterials (ENMs) have a high potential to enter the soil. The soil's major component of clay likely dictates the fate and transport of ENMs in the subsurface. Currently, few studies are available on the fate and transport of nanoparticle silica (nSiO₂) in the presence of clay particles. Therefore, the sedimentation and transport of nSiO₂ with two representative clays (montmorillonite (M) and kaolin (K)) in porous media were investigated in monovalent (Na⁺) and divalent (Ca²⁺) ion solutions with multiple characterizations including SEM/TEM-EDX, zeta potentials, particle sizes and colloid transport modeling. It was shown that nSiO₂-nSiO₂ homoaggregates and nSiO₂-K (or M) heteroaggregates dominated in the nSiO₂-clay nanoaggregate suspension. A distinct decrease in the stability and transport of nSiO₂-M (or K) in NaCl solution and an increase in CaCl₂ occurred when M or K was added to the nSiO₂ suspension at pH 6.0. This was attributed to the faster settlement of the individual M or K in NaCl vs. the better stability in CaCl₂ (compared to nSiO₂ alone). Particularly, more negative individual M platelets occurred in the high NaCl solution until extensive flocculated structures built up, which contributed to the faster deposition of nSiO₂-M compared to nSiO₂-K, even though the nSiO₂-M was more negatively charged. Comparably, the effect of M and K on the fate and transport of nSiO₂ almost disappeared at pH 9.0. The values of the first-order attachment/detachment rate coefficients (k₁/k_1d) and first-order straining coefficient (k₂) obtained from two-site kinetic attachment model fitting are responsible for the deposition of nSiO₂-clay nanoaggregates in sand. This study suggests potential groundwater contamination due to the clay-facilitated transport of ENMs in calcareous soil.

Effects of solution chemistry on the attachment of graphene oxide onto clay minerals

With the increase in production and wide application of graphene oxide (GO), colloidal GO particles are expectantly released into soil and groundwater, where a large number of mineral particles exist. In addition, the porewater chemistry (e.g. organic acid, valence of cation) is a neglected but important aspect to comprehensively investigate the fate of GO. The interactions of GO with three ubiquitous clay minerals (i.e., montmorillonite, kaolinite and diatomite) have been systematically investigated through batch experiments across different solution chemistry conditions. In general, the affinity towards GO is in the order of montmorillonite > kaolinite > diatomite under the same experimental conditions. This observation can be explained by the characteristics of different clay minerals, such as surface charge and surface area. The results indicated that increasing the ionic strength or decreasing the pH enhanced the attachment of GO nanoparticles onto clay minerals as a result of electrostatic interactions. With the increase in concentration of Ca2+, more GO particles were attached onto clay mineral particles. This is caused by complexation between the surface oxygen functional groups of both GO nanoparticles and clay minerals. The presence of 0.1 mM tartaric acid significantly inhibited the attachment of GO onto clay minerals. This is possibly linked to the increased negative charges of the organic acids and the competition between organic acids and GO. The interaction energies were also calculated by applying the classical DLVO theory. The results of this study have helped to understand the behavior and fate of GO in subsurface formations.

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

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

Minimizing Virus Transport in Porous Media by Optimizing Solid Phase Inactivation

The influence of virus type (PRD1 and ΦX174), temperature (flow at 4 and 20°C), a no-flow storage duration (0, 36, 46, and 70 d), and temperature cycling (flow at 20°C and storage at 4°C) on virus transport and fate were investigated in saturated sand-packed columns. The vast majority (84-99.5%) of viruses were irreversibly retained on the sand, even in the presence of deionized water and beef extract at pH = 11. The reversibly retained virus fraction () was small (1.6 × 10 to 0.047) but poses a risk of long-term virus contamination. The value of and associated transport risk was lower at a higher temperature and for increases in the no-flow storage period due to the temperature dependency of the solid phase inactivation. A model that considered advective-dispersive transport, attachment (), detachment (), solid phase inactivation (μ), and liquid phase inactivation (μ) coefficients, and a Langmuirian blocking function provided a good description of the early portion of the breakthrough curve. The removal parameters were found to be in the order of > μ >> μ. Furthermore, μ was an order of magnitude higher than μ for PRD1, whereas μ was two and three orders of magnitude higher than μ for ΦX174 at 4 and 20°C, respectively. Transport modeling with two retention, release, and inactivation sites demonstrated that a small fraction of viruses exhibited a much slower release and solid phase inactivation rate, presumably because variations in the sand and virus surface roughness caused differences in the strength of adhesion. These findings demonstrate the importance of solid phase inactivation, temperature, and storage periods in eliminating virus transport in porous media. This research has potential implications for managed aquifer recharge applications and guidelines to enhance the virus removal by controlling the temperature and aquifer residence time.

Transport of bacteriophage MS2 and PRD1 in saturated dune sand under suboxic conditions

Soil passage of (pretreated) surface water to remove pathogenic microorganisms is a highly efficient process under oxic conditions, reducing microorganism concentrations about 8 log¹⁰ within tens of meters. However, under anoxic conditions, it has been shown that removal of microorganisms can be limited very much. Setback distances for adequate protection of natural groundwater may, therefore, be too short if anoxic conditions apply. Because removal of microorganisms under suboxic conditions is unknown, this research investigated removal of bacteriophage MS2 and PRD1 by soil passage under suboxic conditions at field scale. At the field location (dune area), one injection well and six monitoring wells were installed at different depths along three suboxic flow lines, where oxygen concentrations ranged from 0.4 to 1.7 mg/l and nitrate concentrations ranged from 13 to 16 mg/L. PRD1 and MS2 were injected directly at the corresponding depths and their removal in each flow line was determined. The highest bacteriophage removal was observed in the top layer, with about 9 log removal of MS2, and 7 log removal of PRD1 after 16 meters of aquifer transport. Less removal was observed at 12 m below surface, probably due to a higher groundwater velocity in this coarser grained layer. MS2 was removed more effectively than PRD1 under all conditions. Due to short travel times, inactivation of the phages was limited and the reported log removal was mainly associated with attachment of phages to the aquifer matrix. This study shows that attachment of MS2 and PRD1 is similar for oxic and suboxic sandy aquifers, and, therefore, setback distances used for sandy aquifers under oxic and suboxic conditions provide a similar level of safety. Sticking efficiency and the attachment rate coefficient, as measures for virus attachment, were evaluated as a function of the physico-chemical conditions.

Transport and retention of surfactant- and polymer-stabilized engineered silver nanoparticles in silicate-dominated aquifer material

Packed column experiments were conducted to investigate the transport and blocking behavior of surfactant- and polymer-stabilized engineered silver nanoparticles (Ag-ENPs) in saturated natural aquifer media with varying content of material < 0.063 mm in diameter (silt and clay fraction), background solution chemistry, and flow velocity. Breakthrough curves for Ag-ENPs exhibited blocking behavior that frequently produced a delay in arrival time in comparison to a conservative tracer that was dependent on the physicochemical conditions, and then a rapid increase in the effluent concentration of Ag-ENPs. This breakthrough behavior was accurately described using one or two irreversible retention sites that accounted for Langmuirian blocking on one site. Simulated values for the total retention rate coefficient and the maximum solid phase concentration of Ag-ENPs increased with increasing solution ionic strength, cation valence, clay and silt content, decreasing flow velocity, and for polymer-instead of surfactant-stabilized Ag-ENPs. Increased Ag-ENP retention with ionic strength occurred because of compression of the double layer and lower magnitudes in the zeta potential, whereas lower velocities increased the residence time and decreased the hydrodynamics forces. Enhanced Ag-ENP interactions with cation valence and clay were attributed to the creation of cation bridging in the presence of Ca²⁺. The delay in breakthrough was always more pronounced for polymer-than surfactant-stabilized Ag-ENPs, because of differences in the properties of the stabilizing agents and the magnitude of their zeta-potential was lower. Our results clearly indicate that the long-term transport behavior of Ag-ENPs in natural, silicate dominated aquifer material will be strongly dependent on blocking behavior that changes with the physicochemical conditions and enhanced Ag-ENP transport may occur when retention sites are filled.

Concurrent aggregation and transport of graphene oxide in saturated porous media: Roles of temperature, cation type, and electrolyte concentration

Simultaneous aggregation and retention of nanoparticles can occur during their transport in porous media. In this work, the concurrent aggregation and transport of GO in saturated porous media were investigated under the conditions of different combinations of temperature, cation type (valence), and electrolyte concentration. Increasing temperature (6-24 °C) at a relatively high electrolyte concentration (i.e., 50 mM for Na⁺, 1 mM for Ca²⁺, 1.75 mM for Mg²⁺, and 0.03 and 0.05 mM for Al³⁺) resulted in enhanced GO retention in the porous media. For instance, when the temperature increased from 6 to 24 °C, GO recovery rate decreased from 31.08% to 6.53% for 0.03 mM Al³⁺ and from 27.11% to 0 for 0.05 mM Al³⁺. At the same temperature, increasing cation valence and electrolyte concentration also promoted GO retention. Although GO aggregation occurred in the electrolytes during the transport, the deposition mechanisms of GO retention in the media depended on cation type (valence). For 50 mM Na⁺, surface deposition via secondary minima was the dominant GO retention mechanism. For multivalent cation electrolytes, GO aggregation was rapid and thus other mechanisms such as physical straining and sedimentation also played important roles in controlling GO retention in the media. After passing through the columns, the GO particles in the effluents showed better stability with lower initial aggregation rates. This was probably because less stable GO particles with lower surface charge densities in the porewater were filtered by the porous media, resulting in more stable GO particle with higher surface charge densities in the effluents. An advection-dispersion-reaction model was applied to simulate GO breakthrough curves and the simulations matched all the experimental data well.

Transport of multi-walled carbon nanotubes stabilized by carboxymethyl cellulose and starch in saturated porous media: Influences of electrolyte, clay and humic acid

This study investigated the transport behaviors of carboxymethyl cellulose (CMC) and starch stabilized multi-walled carbon nanotubes (MWNTs) through a saturated quartz sand column in the presence of electrolytes, model clays, and natural organic matter (humic acid) through column breakthrough experiments and model simulations. Both stabilizers, CMC and starch, greatly enhanced the breakthrough of MWNTs, with a full breakthrough plateau (C/C₀) ranging from 0.69 to 0.90 at ionic strength from 0.3 to 10mM. Between the two stabilizers, CMC was more effective in resisting particle deposition, and thus CMC-stabilized MWNTs were more transportable through the medium. While non-stabilized MWNTs were much less transportable and were vulnerable to electrolyte effects (especially Ca²⁺), the stabilized counterparts were much more resistant to the coagulation effects of electrolytes. The presence of colloidal clay particles showed contrasting effects on the transport of bare and stabilized MWNTs. The full breakthrough C/C₀ of bare MWNTs was suppressed by kaolinite and montmorillonite particles from 0.33 to <0.15 with 5mg/L clay, indicating that the presence of both clays enhanced the aggregation and deposition of MWNTs. However, kaolinite particles facilitated the transport of stabilized-MWNTs, while montmorillonite weakened the breakthrough of stabilized MWNTs. Humic acid had less effect on the mobility of stabilized-MWNTs than that of bare MWNTs. The advection-dispersion transport model incorporated with the filtration theory was able to simulate the breakthrough curves and quantitatively interpret the particle deposition. The results can facilitate our understanding of fate and transport of stabilized carbon nanotubes in the environment.

Particle dispersion in porous media: Differentiating effects of geometry and fluid rheology

We investigate the effects of geometric order and fluid rheology on the dispersion of micron-sized particles in two-dimensional microfluidic porous media. Particles suspended in a mixture of glycerol and water or in solutions of partially hydrolyzed polyacrylamide (HPAM) polymers were imaged as they flowed through arrays of microscale posts. From the trajectories of the particles, we calculated the velocity distributions and thereafter obtained the longitudinal and transverse dispersion coefficients. Particles flowed in the shear-thinning HPAM solution through periodic arrays of microposts were more likely to switch between streamlines, due to elastic instabilities. As a result, the distributions of particle velocity were broader in HPAM solutions than in glycerol-water mixtures for ordered geometries. In a disordered array of microposts, however, there was little difference between the velocity distributions obtained in glycerol-water and in HPAM solutions. Correspondingly, particles flowed through ordered post arrays in HPAM solutions exhibited enhanced transverse dispersion. This result suggests that periodic geometric order amplifies the effects of the elasticity-induced velocity fluctuations, whereas geometric disorder of barriers effectively averages out the fluctuations.

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

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

Diffusion of multiwall carbon nanotubes (MWCNTs) through a high density polyethylene (HDPE) geomembrane

The new applications for carbon nanotubes (CNTs) in various fields and consequently their greater production volume have increased their potential release to the environment. Landfills are one of the major locations where carbon nanotubes are expected to be disposed and it is important to ensure that they can limit the release of CNTs. Diffusion of multiwall carbon nanotubes (MWCNTs) dispersed in an aqueous media through a high-density polyethylene (HDPE) geomembrane (as a part of the landfill barrier system) was examined. Based on the laboratory tests, the permeation coefficient was estimated to be less than 5.1×10^-15 m²/s. The potential performance of a HDPE geomembrane and geosynthetic clay liner (GCL) as parts of a composite liner in containing MWCNTs was modelled for six different scenarios. The results suggest that the low value of permeation coefficient of an HDPE geomembrane makes it an effective diffusive barrier for MWCNTs and by keeping the geomembrane defects to minimum during the construction (e.g., number of holes and length of wrinkles) a composite liner commonly used in municipal solid waste landfills will effectively contain MWCNTs.

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

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

In-situ degradation of soil-sorbed 17β-estradiol using carboxymethyl cellulose stabilized manganese oxide nanoparticles: Column studies

This work tested a new remediation technology for in-situ degradation of estrogens by delivering a new class of stabilized manganese oxide (MnO₂) nanoparticles in contaminated soils. The nanoparticles were prepared using a food-grade carboxymethyl cellulose (CMC) as a stabilizer, which was able to facilitate particle delivery into soil. The effectiveness of the technology was tested using 17β-estradiol (E2) as a model estrogen and three sandy loams (SL1, SL2, and SL3) as model soils. Column transport tests showed that the nanoparticles can be delivered in the three soils, though retention of the nanoparticles varied. The nanoparticle retention is strongly dependent on the injection pore velocity. The treatment effectiveness is highly dependent upon the mass transfer rates of both the nanoparticles and contaminants. When the E2-laden soils were treated with 22-130 pore volumes of a 0.174 g/L MnO₂ nanoparticle suspension, up to 88% of water leachable E2 was degraded. The nanoparticles were more effective for soils that offer moderate desorption rates of E2. Decreasing injection velocity or increasing MnO₂ concentration facilitate E2 degradation. The nanoparticles-based technology appears promising for in-situ oxidation of endocrine disruptors in groundwater.

Heteroaggregation of graphene oxide nanoparticles and kaolinite colloids

Graphene oxide (GO) is a material with rapid production growth, and consequently GO nanoparticles are expected to eventually penetrate subsurface formations, where fine mineral particles are in abundance. This study examines the heteroaggregation of GO nanoparticles with kaolinite (KGa-1b) colloids under various conditions. Dynamic batch experiments were conducted in solutions with different pH values (pH=4, 7, and 10), different ionic strengths (I_S=7, 12, and 27mM), and at three controlled temperatures (8, 14, and 25°C). The experimental results showed that a relatively small amount of GO nanoparticles (5-20% of the initial concentration) attached immediately onto KGa-1b colloids, and reached equilibrium in <20min. It was shown that neither temperature nor pH played a significant role in the attachment of GO nanoparticles onto KGa-1b colloids. In contrast, the attachment of GO nanoparticles onto KGa-1b colloids was shown to increase with increasing I_S. Additionally, time-resolved dynamic light scattering (DLS) was used to identify the influence of I_S on heteroaggregation between GO nanoparticles and KGa-1b colloids. The critical coagulation concentration (CCC) for the interaction between GO nanoparticles and KGa-1b colloids was 152mM (NaCl). The interaction energies were calculated, for all experimental conditions, by using measured zeta potentials and applying the classical DLVO theory. The equilibrium experimental data were fitted with a Freundlich isotherm, and the attachment kinetics were described very well with a pseudo-second-order model. Furthermore, thermodynamic analysis revealed that the attachment process was nonspontaneous and exothermic.

Roles of cation valance and exchange on the retention and colloid-facilitated transport of functionalized multi-walled carbon nanotubes in a natural soil

Saturated soil column experiments were conducted to investigate the transport, retention, and release behavior of a low concentration (1 mg L^-1) of functionalized ¹⁴C-labeled multi-walled carbon nanotubes (MWCNTs) in a natural soil under various solution chemistries. Breakthrough curves (BTCs) for MWCNTS exhibited greater amounts of retardation and retention with increasing solution ionic strength (IS) or in the presence of Ca²⁺ in comparison to K⁺, and retention profiles (RPs) for MWCNTs were hyper-exponential in shape. These BTCs and RPs were well described using the advection-dispersion equation with a term for time- and depth-dependent retention. Fitted values of the retention rate coefficient and the maximum retained concentration of MWCNTs were higher with increasing IS and in the presence of Ca²⁺ in comparison to K⁺. Significant amounts of MWCNT and soil colloid release was observed with a reduction of IS due to expansion of the electrical double layer, especially following cation exchange (when K⁺ displaced Ca²⁺) that reduced the zeta potential of MWCNTs and the soil. Analysis of MWCNT concentrations in different soil size fractions revealed that >23.6% of the retained MWCNT mass was associated with water-dispersible colloids (WDCs), even though this fraction was only a minor portion of the total soil mass (2.38%). More MWCNTs were retained on the WDC fraction in the presence of Ca²⁺ than K⁺. These findings indicated that some of the released MWCNTs by IS reduction and cation exchange were associated with the released clay fraction, and suggests the potential for facilitated transport of MWCNT by WDCs.

Temperature dependency of virus and nanoparticle transport and retention in saturated porous media

The influence of temperature on virus (PRD1 and ΦX174) and carboxyl-modified latex nanoparticle (50 and 100nm) attachment was examined in sand-packed columns under various physiochemical conditions. When the solution ionic strength (IS) equaled 10 and 30mM, the attachment rate coefficient (k_att) increased up to 109% (p<0.0002) and the percentage of the sand surface area that contributed to attachment (S_f) increased up to 160% (p<0.002) when the temperature was increased from 4 to 20°C. Temperature effects at IS=10 and 30mM were also dependent on the system hydrodynamics; i.e., enhanced retention at a lower pore water velocity (0.1m/day). Conversely, this same temperature increase had a negligible influence on k_att and S_f values when IS was 1mM or >50mM. An explanation for these observations was obtained from extended interaction energy calculations that considered nanoscale roughness and chemical heterogeneity on the sand surface. Interaction energy calculations demonstrated that the energy barrier to attachment in the primary minimum (∆Φ_a) decreased with increasing IS, chemical heterogeneity, and temperature, especially in the presence of small amounts of nanoscale roughness (e.g., roughness fraction of 0.05 and height of 20nm in the zone of influence). Temperature had a negligible effect on k_att and S_f when the IS=1mM because of the large energy barrier, and at IS=50mM because of the absence of an energy barrier. Conversely, temperature had a large influence on k_att and S_f when the IS was 10 and 30mM because of the presence of a small ∆Φ_a on sand with nanoscale roughness and a chemical (positive zeta potential) heterogeneity. This has large implications for setting parameters for the accurate modeling and transport prediction of virus and nanoparticle contaminants in ground water systems.

Do Goethite Surfaces Really Control the Transport and Retention of Multi-Walled Carbon Nanotubes in Chemically Heterogeneous Porous Media?

Transport and retention behavior of multi-walled carbon nanotubes (MWCNTs) was studied in mixtures of negatively charged quartz sand (QS) and positively charged goethite-coated sand (GQS) to assess the role of chemical heterogeneity. The linear equilibrium sorption model provided a good description of batch results, and the distribution coefficients (K_D) drastically increased with the GQS fraction that was electrostatically favorable for retention. Similarly, retention of MWCNTs increased with the GQS fraction in packed column experiments. However, calculated values of K_D on GQS were around 2 orders of magnitude smaller in batch than packed column experiments due to differences in lever arms associated with hydrodynamic and adhesive torques at microscopic roughness locations. Furthermore, the fraction of the sand surface area that was favorable for retention (S_f) was much smaller than the GQS fraction because nanoscale roughness produced shallow interactions that were susceptible to removal. These observations indicate that only a minor fraction of the GQS was favorable for MWCNT retention. These same observations held for several different sand sizes. Column breakthrough curves were always well described using an advective-dispersive transport model that included retention and blocking. However, depth-dependent retention also needed to be included to accurately describe the retention profile when the GQS fraction was small. Results from this research indicate that roughness primarily controlled the retention of MWCNTs, although goethite surfaces played an important secondary role.