Transport of Pseudomonas putida in a 3-D Bench Scale Experimental Aquifer

Three-dimensional modeling for colloid-facilitated contaminant transport with the effect of mobile and immobile sorbents

A conceptual equilibrium-based mathematical model for colloid-associated contaminant transport has been developed to study the impact of the subsurface environment on contaminant transport through a three-dimensional, saturated, and homogeneous groundwater flow system with uniform flow. The kinetic model's critical limitation is dealing with the more significant number of parameters utilized upon application to larger scales in three-dimensional regions when a series of transport mechanisms are incorporated. Therefore, the present study is the first attempt to study the equilibrium approach in three-dimensional regions to avoid complexities in the model. The current study, however, shows that the mere existence of colloids does not indicate that contaminants will move more quickly; rather, it also depends on how the aqueous phase interacts with the static solid matrix, captured colloid particles, and mobile colloids as well as how colloids interact with stationary solid matrix phase. We noticed that the affinity of contaminants to immobile sorbents (stationary solid matrix and captured colloids) can reduce the transport even in the presence of colloids. Three-dimensional numerical experiments reveal that contaminants infiltrate more in the downward direction in the absence of colloids and can be distributed more in the longitudinal direction and less in the downward direction when colloids are present. The dual nature of colloids is espied here: first, colloids can remove pollutants from a specific area more quickly, and second, in a similar manner, colloids can pollute a specific region more quickly.

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.

Injection of Zerovalent Iron Gels for Aquifer Nanoremediation: Lab Experiments and Modeling

One of the main technical problems faced during field-scale injections of iron microparticles (mZVI) for groundwater nanoremediation is related to their poor colloidal stability and mobility in porous media. In this study, a shear-thinning gel, composed of a mixture of two environmentally friendly biopolymers, i.e., guar gum and xanthan gum, was employed to overcome these limitations. The slurry rheology and particle mobility were characterized by column transport tests. Then, a radial transport experiment was performed to mimic the particle delivery in more realistic conditions. The gel, even at a low polymeric content (1.75 g/L), proved effective in enhancing the mobility of high concentrated mZVI suspensions (20 g/L) in field-like conditions. The high radius of influence (73 cm) and homogeneous iron distribution were achieved by maintaining a low injection overpressure (<0.4 bar). Based only on the information derived from column tests, the MNMs 2018 software (Micro- and Nanoparticle transport, filtration, and clogging Model-Suite) was able to predict the particle distribution and pressure build-up measured in the radial domain. Experimental and simulated results showed good agreement, thus proving that a simplified experimental-modeling procedure based on 1D column tests could be used to effectively upscale the slurry behavior to more representative scales, e.g., radial domains.

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.

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.

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.

Influence of gravity on transport and retention of representative engineered nanoparticles in quartz sand

Four types of NPs: carbon nanotubes and graphene oxide (carbon-based NPs), titanium dioxide and zinc oxide metal-oxide NPs, were utilized to systematically determine the influence of gravity on the transport of NPs in porous media. Packed column experiments for two types of carbon-based NPs were performed under unfavorable conditions in both up-flow (gravity-negative) and down-flow (gravity-positive) orientations, while for two types of metal-oxide NPs, experiments were performed under both unfavorable and favorable conditions in both up-flow and down-flow orientations. Both breakthrough curves and retained profiles of two types of carbon-based NPs in up-flow orientation were equivalent to those in down-flow orientation, indicating that gravity had negligible effect on the transport and retention of carbon-based NPs under unfavorable conditions. In contrast, under both unfavorable and favorable conditions, the breakthrough curves for two types of metal-oxide NPs in down-flow orientation were lower relative to those in up-flow orientation, indicating that gravity could decrease the transport of metal-oxide NPs in porous media. The distinct effect of gravity on the transport and retention of carbon-based and metal-oxide NPs was mainly attributed to the contribution of gravity to the force balance on the NPs in quartz sand. The contribution of gravity was determined by the interplay of the density and sizes of NP aggregates under examined solution conditions.

Modeling dense-colloid and virus cotransport in three-dimensional porous media

A three-dimensional numerical model was developed to investigate the simultaneous transport (cotransport) of dense colloids and viruses in homogeneous, water saturated, porous media with horizontal uniform flow. The dense colloids are assumed to exist in two different phases: suspended in the aqueous phase, and attached reversibly and/or irreversibly onto the solid matrix. The viruses are assumed to exist in four different phases: suspended in aqueous phase, attached onto the solid matrix, attached onto suspended colloids, and attached onto colloids already attached onto the solid matrix. The viruses in each of the four phases are assumed to undergo inactivation with different rates. Moreover, the suspended dense colloids as well as viruses attached onto suspended dense colloids are assumed to exhibit a "restricted" settling velocity as a consequence of the gravitational force; whereas, viruses due to their small sizes and densities are assumed to have negligible "restricted" settling velocity. The governing differential equations were solved numerically with the finite difference schemes, implicitly or explicitly implemented. Model simulations have shown that the presence of dense colloid particles can either enhance or hinder the horizontal transport of viruses, but also can increase the vertical migration of viruses.

Fringe-controlled biodegradation under dynamic conditions: quasi 2-D flow-through experiments and reactive-transport modeling

Biodegradation in contaminated aquifers has been shown to be most pronounced at the fringe of contaminant plumes, where mixing of contaminated water and ambient groundwater, containing dissolved electron acceptors, stimulates microbial activity. While physical mixing of contaminant and electron acceptor by transverse dispersion has been shown to be the major bottleneck for biodegradation in steady-state plumes, so far little is known on the effect of flow and transport dynamics (caused, e.g., by a seasonally fluctuating groundwater table) on biodegradation in these systems. Towards this end we performed experiments in quasi-two-dimensional flow-through microcosms on aerobic toluene degradation by Pseudomonas putida F1. Plume dynamics were simulated by vertical alteration of the toluene plume position and experimental results were analyzed by reactive-transport modeling. We found that, even after disappearance of the toluene plume for two weeks, the majority of microorganisms stayed attached to the sediment and regained their full biodegradation potential within two days after reappearance of the toluene plume. Our results underline that besides microbial growth, also maintenance and dormancy are important processes that affect biodegradation performance under transient environmental conditions and therefore deserve increased consideration in future reactive-transport modeling.

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.

Bayesian process-identification in bacteria transport in porous media

A Bayesian parameter estimation approach is developed for the estimation of joint probability distribution functions for colloid and bacterial fate and transport model parameters describing breakthrough curves (BTCs) obtained through porous media column studies, and is applied to data involving different ionic strength solutions to fit models of differing complexity. Our approach focuses on the simultaneous fitting of a number of BTCs representing different conditions, and it provides a measure of the goodness of model structure, namely Deviance Information Criteria (DIC). Comparison of DIC per model fit enables the evaluation of the significance of various processes through step-wise increases in complexity due to the addition of process model components. We use the method to investigate the transport of both flagellated and non-flagellated strains of Azotobacter vinelandii in a simulated porous media under three ionic strengths. Three different model structures are considered: one without a detachment process and with Langmuirian blocking function, one with detachment, and one with detachment and a second-order blocking function based on random sequential adsorption. First, the model was applied separately to each single BTC. Next, the model was applied comprehensively to the experiments under various ionic strengths, whereas some transport parameters including dispersivity, detachment coefficient, the fraction of cells undergoing irreversible attachment, and the coefficient of the second-order blocking term were assumed to be the same under different ionic strengths. In most cases, including detachment substantially improved the DIC as expected, whereas using the second-order blocking improved DIC for most of the cases when the method was applied to separate BTCs but not when the method was applied collectively to the three BTCs obtained under various ionic strengths. Also, comparing the outcomes of the separate applications of the parameter estimation algorithm versus the collective application indicates that the uncertainty associated with the estimated parameters is substantially smaller when the collective approach is used and also that the estimated parameters are more consistent with the expectations based on the underlying physical processes.