Internal Biofilm Heterogeneities Enhance Solute Mixing and Chemical Reactions in Porous Media.

Bacterial biofilms can form in porous media that are of interest in industrial applications ranging from medical implants to biofilters as well as in environmental applications such as in situ groundwater remediation, where they can be critical locations for biogeochemical reactions. The presence of biofilms modifies porous media topology and hydrodynamics by clogging pores and consequently solutes transport and reactions kinetics. The interplay between highly heterogeneous flow fields found in porous media and microbial behavior, including biofilm growth, results in a spatially heterogeneous biofilm distribution in the porous media as well as internal heterogeneity across the thickness of the biofilm. Our study leverages highly resolved three-dimensional X-ray computed microtomography images of bacterial biofilms in a tubular reactor to numerically compute pore-scale fluid flow and solute transport by considering multiple equivalent stochastically generated internal permeability fields for the biofilm. We show that the internal heterogeneous permeability mainly impacts intermediate velocities when compared with homogeneous biofilm permeability. While the equivalent internal permeability fields of the biofilm do not impact fluid-fluid mixing, they significantly control a fast reaction. For biologically driven reactions such as nutrient or contaminant uptake by the biofilm, its internal permeability field controls the efficiency of the process. This study highlights the importance of considering the internal heterogeneity of biofilms to better predict reactivity in industrial and environmental bioclogged porous systems.

Direct Measurements of the Forces between Silver and Mica in Humic Substance-Rich Solutions.

Deposition of engineered nanoparticles onto porous media from flowing suspensions is important for soil and groundwater quality. The deposition mechanism is controlled by interaction forces between particles and collectors. We investigated the origin and magnitude of opposing forces between silver and mica surfaces (representing nanosilver and sand grains) in solutions relevant to agricultural soils with direct measurements using a surface force apparatus. Solutions of variable NaNO3, Ca(NO3)2, and humic acid (HA) concentrations were used to differentiate individual contributing forces and quantify surface properties. The measured Hamaker constant for silver-water-mica was consistent with Lifshitz theory. Our results indicate that HA forms an adsorbed surface layer, but its charge, thicknesses, compressibility, and mass are significantly larger on mica than silver. Ca2+ primarily reduced the differences between the initially adsorbed HA layer properties on each surface, making them more similar. Force-distance profiles indicate that, when silver-mica systems were exposed to HA, osmotic-steric, electrostatic, and van der Waals forces dominate. Soft particle theory was deemed inappropriate for this system. Derjaguin's approximation was utilized to translate force measurements into interaction energy between nanosilver particles and mica collectors. We propose attachment efficiency estimates from measured surface properties, which suggest high particle mobility when nanosilver is applied to HA-rich agricultural soils with modest ionic strength.

Morphology of Shear-Induced Colloidal Aggregates in Porous Media: Consequences for Transport, Deposition, and Re-entrainment.

Colloid deposition in granular media is relevant to numerous environmental problems. Classic filtration models assume a homogeneous pore space and largely ignore colloid aggregation. However, substantial evidence exists on the ubiquity of aggregation within porous media, suggesting that deposition is enhanced by it. This work studies the deposition process in relation to aggregate size and structure. We demonstrate that aggregation is induced at typical groundwater velocities by comparing the repulsive DLVO force between particle pairs to the hydrodynamic shear force opposing it. Column experiments imaged with high-resolution X-ray computed tomography are used to measure aggregate structure and describe their morphology probability distribution and spatial distribution. Aggregate volume and surface area are found to be power-law distributed, while Feret diameter is exponentially distributed with some flow rate dependencies caused by erosion and restructuring by the fluid shear. Furthermore, size and shape of aggregates are heterogeneous in depth, where a small number of large aggregates control the concentration versus depth profile shape. The range of aggregate fractal dimensions found (2.22-2.42) implies a high potential for restructuring or breaking during transport. Shear-induced aggregation is not currently considered in macroscopic models for particle filtration, yet is critical to consider in the processes that control deposition.

Effect of hydrofracking fluid on colloid transport in the unsaturated zone.

Hydraulic fracturing is expanding rapidly in the US to meet increasing energy demand and requires high volumes of hydrofracking fluid to displace natural gas from shale. Accidental spills and deliberate land application of hydrofracking fluids, which return to the surface during hydrofracking, are common causes of environmental contamination. Since the chemistry of hydrofracking fluids favors transport of colloids and mineral particles through rock cracks, it may also facilitate transport of in situ colloids and associated pollutants in unsaturated soils. We investigated this by subsequently injecting deionized water and flowback fluid at increasing flow rates into unsaturated sand columns containing colloids. Colloid retention and mobilization was measured in the column effluent and visualized in situ with bright field microscopy. While <5% of initial colloids were released by flushing with deionized water, 32-36% were released by flushing with flowback fluid in two distinct breakthrough peaks. These peaks resulted from 1) surface tension reduction and steric repulsion and 2) slow kinetic disaggregation of colloid flocs. Increasing the flow rate of the flowback fluid mobilized an additional 36% of colloids, due to the expansion of water filled pore space. This study suggests that hydrofracking fluid may also indirectly contaminate groundwater by remobilizing existing colloidal pollutants.

Surfactant-mediated control of colloid pattern assembly and attachment strength in evaporating droplets.

This study demonstrates that the pattern assembly and attachment strength of colloids in an evaporating sessile droplet resting on a smooth substrate can be controlled by adding nonionic solutes (surfactant) to the solution. As expected, increasing the surfactant concentration leads to a decrease in initial surface tension of the drop, σ(0). For the range of initial surface tensions investigated (39-72 mN m(-1)), three distinct deposition patterns were produced: amorphous stains (σ(0) = 63-72 mN m(-1)), coffee-ring stains (σ(0) = 48-53 mN m(-1)), and concentric rings (σ(0) = 39-45 mN m(-1)). A flow-displacement system was used to measure the attachment strength of the dried colloids. Characteristic drying regimes associated with the three unique pattern formations are attributed to abrupt transitions of contact line dynamics during evaporation. The first transition from slipping- to pinned-contact line was found to be a direct result of the competition between mechanical instability of the droplet and the friction generated by pinned colloids at the contact line. The second transition from pinned- to recurrent-stick-rip-slip-contact line was caused by repeated liquid film rupturing from evaporation-intensified surfactant concentration. Data from flow-displacement tests indicate that attachment strength of dried particles is strongest for amorphous stains (lowest surfactant concentration) and weakest for concentric rings (highest surfactant concentration). The mechanism behind these observations was ascribed to the formation and adsorption of micelles onto colloid and substrate surfaces as the droplet solution evaporates. The range of attachment forces observed between the colloids and the solid substrate were well captured by extended-DLVO interactions accounting for van der Waals attraction, electric double layer repulsion, and micelle-protrusion repulsion. Both empirical and theoretical results suggest that an increasingly dense layer of adsorbed micellar-protrusions on colloid and substrate surfaces acts as a physical barrier that hinders strong van der Waals attractive interactions at close proximity. Thereby, colloid stains dried at higher surfactant concentrations are more easily detached from the substrate when dislodging forces are applied than stains dried at lower surfactant concentrations.

Quantification of colloid retention and release by straining and energy minima in variably saturated porous media.

The prediction of colloid transport in unsaturated porous media in the presence of large energy barrier is hampered by scant information of the proportional retention by straining and attractive interactions at surface energy minima. This study aims to fill this gap by performing saturated and unsaturated column experiments in which colloid pulses were added at various ionic strengths (ISs) from 0.1 to 50 mM. Subsequent flushing with deionized water released colloids held at the secondary minimum. Next, destruction of the column freed colloids held by straining. Colloids not recovered at the end of the experiment were quantified as retained at the primary minimum. Results showed that net colloid retention increased with IS and was independent of saturation degree under identical IS and Darcian velocity. Attachment rates were greater in unsaturated columns, despite an over 3-fold increase in pore water velocity relative to saturated columns, because additional retention at the readily available air-associated interfaces (e.g., the air-water-solid [AWS] interfaces) is highly efficient. Complementary visual data showed heavy retention at the AWS interfaces. Retention by secondary minima ranged between 8% and 46% as IS increased, and was greater for saturated conditions. Straining accounted for an average of 57% of the retained colloids with insignificant differences among the treatments. Finally, retention by primary minima ranged between 14% and 35% with increasing IS, and was greater for unsaturated conditions due to capillary pinning.

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

HYPOTHESES: The transport behavior of colloids in subsurface porous media is altered by surface chemical and physical heterogeneities. Understanding the mechanisms involved and distribution outcomes is crucial to assess and control groundwater contamination. The multi-scale processes that broaden residence time distribution for particles in the medium are here succinctly described with an upscaling model. Experiments/model: The spatial distribution of silver particles along glass bead-packed columns obtained from X-ray micro-computed tomography and a mechanistic upscaling model were used to study colloid retention across interface-, collector-, pore-, and Darcy-scales. Simulated energy profiles considering variable colloid-grain interactions were used to determine collector efficiencies from particle trajectories via full force-torque balance. Rate coefficients were determined from collector efficiencies to parameterize the advective-dispersive-reactive model that reports breakthrough curves and depth profiles. FINDINGS: Our results indicate that: (i) with surface heterogeneity, individual colloid-grain interactions are non-unique and span from repulsive to attractive extremes; (ii) experimentally observed spatial positioning of retention at grain-water interfaces and grain-to-grain contacts is governed respectively by mechanistic attachment to the grain surface and retention without contact at rear-flow stagnation zones, and (iii) experimentally observed non-monotonic retention profiles and heavy-tailed breakthrough curves can be modeled with explicit implementation of heterogeneity at smaller scales.

Correlation equation for predicting attachment efficiency (α) of organic matter-colloid complexes in unsaturated porous media.

Naturally occurring polymers such as organic matter have been known to inhibit aggregation and promote mobility of suspensions in soil environments by imparting steric stability. This increase in mobility can significantly reduce the water filtering capacity of soils, thus jeopardizing a primary function of the vadose zone. Improvements to classic filtration theory have been made to account for the known decrease in attachment efficiency of electrostatically stabilized particles, and more recently, of sterically stabilized particles traveling through simple and saturated porous media. In the absence of an established unsaturated transport expression, and in the absence of applicable theoretical approaches for suspensions with asymmetric and nonindifferent electrolytes, this study presents an empirical correlation to predict attachment efficiency (α) for electrosterically stabilized suspensions in unsaturated systems in the presence of nonideal electrolytes. We show that existing models fall short in estimating polymer-coated colloid deposition in unsaturated media. This deficiency is expected given that the models were developed for saturated conditions where the mechanisms controlling colloid deposition are significantly different. A new correlation is derived from unsaturated transport data and direct characterization of microspheres coated with natural organic matter over a range of pH and CaCl(2) concentrations. The improvements to existing transport models include the following: adjustment for a restricted liquid-phase in the medium, development of a quantitative term to account for unsaturated transport phenomena, and adjustments in the relative contribution of steric stability parameters based on direct measurements of the adsorbed polymer layer characteristics. Differences in model formulation for correlations designed for saturated systems and the newly proposed correlation for unsaturated systems are discussed, and the performance of the new model against a comprehensive set of experimental observations is evaluated.

Biofilms in 3D porous media: Delineating the influence of the pore network geometry, flow and mass transfer on biofilm development.

This study investigates the functional correspondence between porescale hydrodynamics, mass transfer, pore structure and biofilm morphology during progressive biofilm colonization of a porous medium. Hydrodynamics and the structure of both the porous medium and the biofilm are experimentally measured with 3D particle tracking velocimetry and micro X-ray Computed Tomography, respectively. The analysis focuses on data obtained in a clean porous medium after 36 h of biofilm growth. Registration of the particle tracking and X-ray data sets allows to delineate the interplay between porous medium geometry, hydrodynamic and mass transfer processes on the morphology of the developing biofilm. A local analysis revealed wide distributions of wall shear stresses and concentration boundary layer thicknesses. The spatial distribution of the biofilm patches uncovered that the wall shear stresses controlled the biofilm development. Neither external nor internal mass transfer limitations were noticeable in the considered system, consistent with the excess supply of nutrient and electron acceptors. The wall shear stress remained constant in the vicinity of the biofilm but increased substantially elsewhere.

Correction: Biofilm imaging in porous media by laboratory X-Ray tomography: Combining a non-destructive contrast agent with propagation-based phase-contrast imaging tools.

[This corrects the article DOI: 10.1371/journal.pone.0180374.].

Biofilm imaging in porous media by laboratory X-Ray tomography: Combining a non-destructive contrast agent with propagation-based phase-contrast imaging tools.

X-ray tomography is a powerful tool giving access to the morphology of biofilms, in 3D porous media, at the mesoscale. Due to the high water content of biofilms, the attenuation coefficient of biofilms and water are very close, hindering the distinction between biofilms and water without the use of contrast agents. Until now, the use of contrast agents such as barium sulfate, silver-coated micro-particles or 1-chloronaphtalene added to the liquid phase allowed imaging the biofilm 3D morphology. However, these contrast agents are not passive and potentially interact with the biofilm when injected into the sample. Here, we use a natural inorganic compound, namely iron sulfate, as a contrast agent progressively bounded in dilute or colloidal form into the EPS matrix during biofilm growth. By combining a very long source-to-detector distance on a X-ray laboratory source with a Lorentzian filter implemented prior to tomographic reconstruction, we substantially increase the contrast between the biofilm and the surrounding liquid, which allows revealing the 3D biofilm morphology. A comparison of this new method with the method proposed by Davit et al (Davit et al., 2011), which uses barium sulfate as a contrast agent to mark the liquid phase was performed. Quantitative evaluations between the methods revealed substantial differences for the volumetric fractions obtained from both methods. Namely, contrast agent-biofilm interactions (e.g. biofilm detachment) occurring during barium sulfate injection caused a reduction of the biofilm volumetric fraction of more than 50% and displacement of biofilm patches elsewhere in the column. Two key advantages of the newly proposed method are that passive addition of iron sulfate maintains the integrity of the biofilm prior to imaging, and that the biofilm itself is marked by the contrast agent, rather than the liquid phase as in other available methods. The iron sulfate method presented can be applied to understand biofilm development and bioclogging mechanisms in porous materials and the obtained biofilm morphology could be an ideal basis for 3D numerical calculations of hydrodynamic conditions to investigate biofilm-flow coupling.

Reverse engineering of biochar.

This study underpins quantitative relationships that account for the combined effects that starting biomass and peak pyrolysis temperature have on physico-chemical properties of biochar. Meta-data was assembled from published data of diverse biochar samples (n=102) to (i) obtain networks of intercorrelated properties and (ii) derive models that predict biochar properties. Assembled correlation networks provide a qualitative overview of the combinations of biochar properties likely to occur in a sample. Generalized Linear Models are constructed to account for situations of varying complexity, including: dependence of biochar properties on single or multiple predictor variables, where dependence on multiple variables can have additive and/or interactive effects; non-linear relation between the response and predictors; and non-Gaussian data distributions. The web-tool Biochar Engineering implements the derived models to maximize their utility and distribution. Provided examples illustrate the practical use of the networks, models and web-tool to engineer biochars with prescribed properties desirable for hypothetical scenarios.

Removal of sulfamethoxazole and sulfapyridine by carbon nanotubes in fixed-bed columns.

Sulfamethoxazole (SMX) and sulfapyridine (SPY), two representative sulfonamide antibiotics, have gained increasing attention because of the ecological risks these substances pose to plants, animals, and humans. This work systematically investigated the removal of SMX and SPY by carbon nanotubes (CNTs) in fixed-bed columns under a broad range of conditions including: CNT incorporation method, solution pH, bed depth, adsorbent dosage, adsorbate initial concentration, and flow rate. Fixed-bed experiments showed that pH is a key factor that affects the adsorption capacity of antibiotics to CNTs. The Bed Depth Service Time model describes well the relationship between service time and bed depth and can be used to design appropriate column parameters. During fixed-bed regeneration, small amounts of SMX (3%) and SPY (9%) were irreversibly bonded to the CNT/sand porous media, thus reducing the column capacity for subsequent reuse from 67.9 to 50.4 mg g(-1) for SMX and from 91.9 to 72.9 mg g(-1) for SPY. The reduced column capacity resulted from the decrease in available adsorption sites and resulting repulsion (i.e., blocking) of incoming antibiotics from those previously adsorbed. Findings from this study demonstrate that fixed-bed columns packed with CNTs can be efficiently used and regenerated to remove antibiotics from water.

Effect of surface modification on single-walled carbon nanotube retention and transport in saturated and unsaturated porous media.

This work investigated the effect of different surface modification methods, including oxidization, surfactant coating, and humic acid coating, on single-walled carbon nanotube (SWNT) stability and their mobility in granular porous media under various conditions. Characterization and stability studies demonstrated that the three surface modification methods were all effective in solubilizing and stabilizing the SWNTs in aqueous solutions. Packed sand column experiments showed that although the three surface medication methods showed different effect on the retention and transport of SWNTs in the columns, all the modified SWNTs were highly mobile. Compared with the other two surface modification methods, the humic acid coating method introduced the highest mobility to the SWNTs. While reductions in moisture content in the porous media could promote the retention of the surface modified SWNTs in some sand columns, results from bubble column experiment suggested that only oxidized SWNTs were retention in unsaturated porous media through attachment on air-water interfaces. Other mechanisms such as grain surface attachment and thin-water film straining could also be responsible for the retention of the SWNTs in unsaturated porous media. An advection-dispersion model was successfully applied to simulate the experimental data of surface modified SWNT retention and transport in porous media.

Deposition and transport of functionalized carbon nanotubes in water-saturated sand columns.

Knowledge of the fate and transport of functionalized carbon nanotubes (CNTs) in porous media is crucial to understand their environmental impacts. In this study, laboratory column and modeling experiments were conducted to mechanistically compare the retention and transport of two types of functionalized CNTs (i.e., single-walled nanotubes and multi-walled nanotubes) in acid-cleaned, baked, and natural sand under unfavorable conditions. The CNTs were highly mobile in the acid-cleaned sand columns but showed little transport in the both natural and baked sand columns. In addition, the retention of the CNTs in the both baked and natural sand was strong and almost irreversible even after reverse, high-velocity, or surfactant flow flushing. Both experimental and modeling results showed that pH is one of the factors dominating CNT retention and transport in natural and baked sand. Retention of the functionalized CNTs in the natural and baked sand columns reduced dramatically when the system pH increased. Our results suggest that the retention and transport of the functionalized CNTs in natural sand porous media were mainly controlled by strong surface deposition through the electrostatic and/or hydrogen-bonding attractions between surface function groups of the CNTs and metal oxyhydroxide impurities on the sand surfaces.

Colloid retention at the meniscus-wall contact line in an open microchannel.

Colloid retention mechanisms in partially saturated porous media are currently being researched with an array of visualization techniques. These visualization techniques have refined our understanding of colloid movement and retention at the pore scale beyond what can be obtained from breakthrough experiments. One of the remaining questions is what mechanisms are responsible for colloid immobilization at the triple point where air, water, and soil grain meet. The objective of this study was to investigate how colloids are transported to the air-water-solid (AWS) contact line in an open triangular microchannel, and then retained as a function of meniscus contact angle with the wall and solution ionic strength. Colloid flow path, meniscus shape and meniscus-wall contact angle, and colloid retention at the AWS contact line were visualized and quantified with a confocal microscope. Experimental results demonstrated that colloid retention at the AWS contact line was significant when the meniscus-wall contact angle was less than 16°, but was minimal for the meniscus-wall contact angles exceeding 20°. Tracking of individual colloids and computational hydrodynamic simulation both revealed that for small contact angles (e.g., 12.5°), counter flow and flow vortices formed near the AWS contact line, but not for large contact angles (e.g., 28°). This counter flow helped deliver the colloids to the wall surface just below the contact line. In accordance with DLVO and hydrodynamic torque calculations, colloid movement may be stopped when the colloid reached the secondary minimum at the wall near the contact line. However, contradictory to the prediction of the torque analysis, colloid retention at the AWS contact line decreased with increasing ionic strength for contact angles of 10-20°, indicating that the air-water interface was involved through both counter flow and capillary force. We hypothesized that capillary force pushed the colloid through the primary energy barrier to the primary minimum to become immobilized, when small fluctuations in water level stretched the meniscus over the colloid. For large meniscus-wall contact angles counter flow was not observed, resulting in less colloid retention, because a smaller number of colloids were transported to the contact line.

Impact of dissolved organic matter on colloid transport in the vadose zone: deterministic approximation of transport deposition coefficients from polymeric coating characteristics.

Although numerous studies have been conducted to discern colloid transport and stability processes, the mechanistic understanding of how dissolved organic matter (DOM) affects colloid fate in unsaturated soils (i.e., the vadose zone) remains unclear. This study aims to bridge the gap between the physicochemical responses of colloid complexes and porous media interfaces to solution chemistry, and the effect these changes have on colloid transport and fate. Measurements of adsorbed layer thickness, density, and charge of DOM-colloid complexes and transport experiments with tandem internal process visualization were conducted for key constituents of DOM, humic (HA) and fulvic acids (FA), at acidic, neutral and basic pH and two CaCl(2) concentrations. Polymeric characteristics reveal that, of the two tested DOM constituents, only HA electrosterically stabilizes colloids. This stabilization is highly dependent on solution pH which controls DOM polymer adsorption affinity, and on the presence of Ca(+2) which promotes charge neutralization and inter-particle bridging. Transport experiments indicate that HA improved colloid transport significantly, while FA only marginally affected transport despite having a large effect on particle charge. A transport model with deposition and pore-exclusion parameters fit experimental breakthrough curves well. Trends in deposition coefficients are correlated to the changes in colloid surface potential for bare colloids, but must include adsorbed layer thickness and density for sterically stabilized colloids. Additionally, internal process observations with bright field microscopy reveal that, under optimal conditions for retention, experiments with FA or no DOM promoted colloid retention at solid-water interfaces, while experiments with HA enhanced colloid retention at air-water interfaces, presumably due to partitioning of HA at the air-water interface and/or increased hydrophobic characteristics of HA-colloid complexes.

Colloid transport and retention in unsaturated porous media: effect of colloid input concentration.

Colloids play an important role in facilitating transport of adsorbed contaminants in soils. Recent studies showed that under saturated conditions colloid retention was a function of its concentration. It is unknown if this is the case under unsaturated conditions. In this study, the effect of colloid concentration on colloid retention was investigated in unsaturated columns by increasing concentrations of colloid influents with varying ionic strength. Colloid retention was observed in situ by bright field microscopy and quantified by measuring colloid breakthrough curves. In our unsaturated experiments, greater input concentrations resulted in increased colloid retention at ionic strength above 0.1 mM, but not in deionized water (i.e., 0 mM ionic strength). Bright field microscope images showed that colloid retention mainly occurred at the solid-water interface and wedge-shaped air-water-solid interfaces, whereas the retention at the grain-grain contacts was minor. Some colloids at the air-water-solid interfaces were rotating and oscillating and thus trapped. Computational hydrodynamic simulation confirmed that the wedge-shaped air-water-solid interface could form a "hydrodynamic trap" by retaining colloids in its low velocity vortices. Direct visualization also revealed that colloids once retained acted as new retention sites for other suspended colloids at ionic strength greater than 0.1 mM and thereby could explain the greater retention with increased input concentrations. Derjaguin-Landau-Verwey-Overbeek (DLVO) energy calculations support this concept. Finally, the results of unsaturated experiments were in agreement with limited saturated experiments under otherwise the same conditions.