Use of tracer tests to investigate changes in flow and transport properties due to bioclogging of porous media

Role of Soil Biofilms in Clogging and Fate of Pharmaceuticals: A Laboratory-Scale Column Experiment

Contamination of groundwater with pharmaceutical active compounds (PhACs) increased over the last decades. Potential pathways of PhACs to groundwater include techniques such as irrigation, managed aquifer recharge, or bank filtration as well as natural processes such as losing streams of PhACs-loaded source waters. Usually, these systems are characterized by redox-active zones, where microorganisms grow and become immobilized by the formation of biofilms, structures that colonize the pore space and decrease the infiltration capacities, a phenomenon known as bioclogging. The goal of this work is to gain a deeper understanding of the influence of soil biofilms on hydraulic conductivity reduction and the fate of PhACs in the subsurface. For this purpose, we selected three PhACs with different physicochemical properties (carbamazepine, diclofenac, and metoprolol) and performed batch and column experiments using a natural soil, as it is and with the organic matter removed, under different biological conditions. We observed enhanced sorption and biodegradation for all PhACs in the system with higher biological activity. Bioclogging was more prevalent in the absence of organic matter. Our results differ from works using artificial porous media and thus reveal the importance of utilizing natural soils with organic matter in studies designed to assess the role of soil biofilms in bioclogging and the fate of PhACs in soils.

Cotransport of micro- and nano-plastics with chlortetracycline hydrochloride in saturated porous media: Effects of physicochemical heterogeneities and ionic strength

Global production and use of plastics have resulted in the wide dissemination of micro- and nano-plastics (MNPs) to the natural environment. Potentially acting as a vector, the role of MNPs on the fate and transport of environmental pollutants (e.g., antibiotics such as chlortetracycline hydrochloride; CTC) has garnered global concern recently. Herein, the cotransport of MNPs and CTC in columns packed with uncoated sand or soil colloid-coated sand (SCCS) under different degrees of physicochemical heterogeneity and ionic strength was systematically explored. Our results show that MNPs and CTC inhibit the transport of each other when they coexist. The adsorption of CTC onto sand grains, soil colloids, and MNPs, as well as the aggregation of MNPs in the presence of CTC could be the major contributors to the enhanced retention of CTC and MNPs. In SCCS with different degrees of soil colloid coating, the adsorption of CTC on soil colloids is critical to influence the transport of CTC, and the nonlinear retention of MNPs to soil colloids is mainly attributed to the alteration of collector surface roughness by soil colloids. High ionic strength slightly facilitates CTC transport due to the competition for adsorption sites and the formation of CTC macromolecules, but significantly inhibits MNPs transport by suppressing the electrostatic double layers based on colloid stability theory. Consequently, the cotransport of MNPs and CTC is governed by the coupled interplay of collector surface roughness and chemical heterogeneity, due to the soil colloid coatings and the adsorbed CTC on the surfaces associated with solution chemistries such as ionic strength. Increased cotransport of MNPs and CTC occurred under a higher concentration of MNPs due to a larger number of adsorption sites for CTC. Our findings advance the current understanding of the complex cotransport of MNPs and antibiotics in the environment. This information is valuable for understanding contaminant fate and formulating strategies for environmental remediation due to the contamination of MNPs and co-occurring contaminants.

Contributions of biofilm-induced flow heterogeneities to solute retention and anomalous transport features in porous media

Microbial biofilms are ubiquitous within porous media and the dynamics of their growth influence surface and subsurface flow patterns which impacts the physical properties of porous media and large-scale transport of solutes. A two-dimensional pore-scale numerical model was used to evaluate the impact of biofilm-induced flow heterogeneities on conservative transport. Our study integrates experimental biofilm images of Paenibacillus 300A strain in a microfluidic device packed with cylindrical grains in a hexagonal distribution, with mathematical modeling. Biofilm is represented as a synthetic porous structure with locally varying physical properties that honors the impact of biofilm on the porous medium. We find that biofilm plays a major role in shaping the observed conservative transport dynamics by enhancing anomalous characteristics. More specifically, when biofilm is present, the pore structure in our geometry becomes more spatially correlated. We observe intermittent behavior in the Lagrangian velocities that switches between fast transport periods and long trapping events. Our results suggest that intermittency enhances solute spreading in breakthrough curves which exhibit extreme anomalous slope at intermediate times and very marked late solute arrival due to solute retention. The efficiency of solute retention by the biofilm is controlled by a transport regime which can extend the tailing in the concentration breakthrough curves. These results indicate that solute retention by the biofilm exerts a strong control on conservative solute transport at pore-scale, a role that to date has not received enough attention.

Reactive Transport: A Review of Basic Concepts with Emphasis on Biochemical Processes

Reactive transport (RT) couples bio-geo-chemical reactions and transport. RT is important to understand numerous scientific questions and solve some engineering problems. RT is highly multidisciplinary, which hinders the development of a body of knowledge shared by RT modelers and developers. The goal of this paper is to review the basic conceptual issues shared by all RT problems, so as to facilitate advancement along the current frontier: biochemical reactions. To this end, we review the basic equations to indicate that chemical systems are controlled by the set of equilibrium reactions, which are easy to model, but whose rate is controlled by mixing. Since mixing is not properly represented by the standard advection-dispersion equation (ADE), we conclude that this equation is poor for RT. This leads us to review alternative transport formulations, and the methods to solve RT problems using both the ADE and alternative equations. Since equilibrium is easy, difficulties arise for kinetic reactions, which is especially true for biochemistry, where numerous challenges are open (how to represent microbial communities, impact of genomics, effect of biofilms on flow and transport, etc.). We conclude with the basic eleven conceptual issues that we consider fundamental for any conceptually sound RT effort.

Impact of compost reactive layer on hydraulic transport and C & N cycles: Biogeochemical modeling of infiltration column experiments

Managed Aquifer Recharge (MAR) is a key strategy to increase freshwater resources in many regions facing water scarcity. MAR issues are related to both quantity and quality of the infiltrating water. In most countries, very high quality of the infiltrating water is required, to limit the impact on the aquifer geochemistry. In this paper, the possibility of injecting water of lower quality in the aquifer and letting the biogeochemical reactions take place in order to enhance its quality is explored. Here, we present the fate of nutrients (C, N) in the biogeochemical system of a reactive barrier formed by mixture of different proportions of sand and compost, supplied with treated wastewater to mimic MAR. An integrated conceptual model involving the nutrient cycles and biomass dynamics (auto- and heterotrophic) was developed, and then tested with a number of solute transport experiments in columns with different compost fraction in the column filling. The model incorporated both saturation and inhibition processes (regarding the nutrients and their byproducts) to provide a comprehensive picture of the nutrient dynamics within the column. The model developed (three if considering the 3 column setups) allowed to discriminate the processes that govern the fate of nutrients in relation with the compost enhancing long-term nutrient degradation, yet hindering hydraulic parameters that affect infiltration rates.

High-frequency, in situ sampling of field woodchip bioreactors reveals sources of sampling error and hydraulic inefficiencies

Woodchip bioreactors are a practical, low-cost technology for reducing nitrate (NO₃) loads discharged from agriculture. Traditional methods of quantifying their performance in the field mostly rely on low-frequency, time-based (weekly to monthly sampling interval) or flow-weighted sample collection at the inlet and outlet, creating uncertainty in their performance and design by providing incomplete information on flow and water chemistry. To address this uncertainty, two field bioreactors were monitored in the US and New Zealand using high-frequency, multipoint sampling for in situ monitoring of NO₃-N concentrations. High-frequency monitoring (sub hourly interval) at the inlet and outlet of both bioreactors revealed significant variability in volumetric removal rates and percent reduction, with percent reduction varying by up to 25 percentage points within a single flow event. Time series of inlet and outlet NO₃ showed significant lag in peak concentrations of 1-3 days due to high hydraulic residence time, where calculations from instantaneous measurements produced erroneous estimates of performance and misleading relationships between residence time and removal. Internal porewater sampling wells showed differences in NO₃ concentration between shallow and deep zones, and "hot spot" zones where peak NO₃ removal co-occurred with dissolved oxygen depletion and dissolved organic carbon production. Tracking NO₃ movement through the profile showed preferential flow occurring with slower flow in deeper woodchips, and slower flow further from the most direct flowpath from inlet to outlet. High-frequency, in situ data on inlet and outlet time series and internal porewater solute profiles of this initial work highlight several key areas for future research.

Upscaling of microbially driven first-order reactions in heterogeneous porous media

Reactions mediated by microorganisms determine the fate of many chemicals in natural porous media. At the pore scale, the distribution of chemicals and microorganisms is not homogeneous, leading to heterogeneous distribution of microbial activities at the pore scale. We conducted pore scale reactive transport simulations to investigate the scaling of microbially mediated consumption reaction rates under a range of flow and reaction conditions. The results reveal that the scaling effects largely depended on Péclet and Damköhler numbers. Consumption rate estimates based on volume-averaged concentrations and reaction kinetics overestimated the true volumetric reaction rates, and large-sized biomass aggregates intensified these scaling errors. In contrast, the macroscopic rates estimated via flux-weighted concentrations underestimated the true volumetric reaction rates, with large microbial aggregates reducing scaling errors. This study also demonstrated that macroscopic rate estimates can be improved by combining information on the reaction kinetics with the flux-weighted concentrations.

In-situ microbially induced Ca2+ -alginate polymeric sealant for seepage control in porous materials

This paper presents a novel approach of using in-situ microbially induced Ca2+ -alginate polymeric sealant for seepage control in porous materials. This process comprises two steps: (i) generation of insoluble calcium carbonate inside the pores of porous materials (such as sand) through a microbially induced carbonate precipitation (MICP) process in-situ and (ii) injection of sodium alginate for in-situ gelation via reaction between alginate and Ca2+ ions. The experimental results showed that the hydraulic conductivity/permeability of sand decreased with the increase in alginate concentration. When 5% alginate was used with a CaCO3 concentration of 0.18 g g-1 sand, the permeability of the alginate-treated sand reduced from 5.0 × 10-4 to 2.2 × 10-9  m s-1 . The scanning electron microscopy images revealed that a film-type coating was formed around sand particles with spherical round crystals embedded. Furthermore, the in-situ formed Ca-alginate polymeric sealant can also be used for the removal of Cu2+ ion and suspended particles from contaminated water by more than 90%. Built on the current research, the envisioned practical application of the proposed method may include clogging fractured rock, reducing seepage and prevent piping through dams, excavation dewatering, and forming barriers for remediating specific contaminants.

Removal of micropollutants in biofilters: Hydrodynamic effects on biofilm assembly and functioning

Global water resources contain a variety of micropollutants (MPs), including pharmaceuticals, personal care products, and pesticides. This study investigated the removal of MPs during drinking water production by means of biofiltration. The objective of this work was to investigate the influence of hydrodynamics on biofilm growth and development in a biofiltration process and the consequent effect on MP biotransformation rates. We operated three groups of biofiltration columns continuously for 381 days under three distinct hydrodynamic regimes (superficial velocity: 10, 20, 40 cm h^-1) and fed them a mixture of 29 micropollutants at low concentrations. Total protein concentrations were used as a surrogate measurement for attached biomass and periodic tracer experiments were conducted to estimate dispersivity and assess changes in the depth of the biological zone in each biofilter. These data revealed significant differences in biofilm assembly among the biofilters; higher superficial velocities led to less concentrated surface biomass but a deeper biological zone and more total biomass. Eleven of the 29 MPs were biotransformed and nine of those could be evaluated to estimate biotransformation rates. The second-order rate constants for all nine MPs were not significantly different among the hydrodynamic regimes. However, a depth-based analysis of biotransformation rates revealed significantly greater second-order rate constants for 5 of the MPs at increasing biofilter depths, suggesting that sparse microbial communities found in deeper and more oligotrophic biofilters had a greater activity for the biotransformation of these MPs. The identification of several transformation products at similar relative distributions suggests that the greater activity was not the result of changing metabolic processes under more oligotrophic conditions. These results improve our fundamental understanding of biofilm assembly and functioning in biofiltration processes.

Denitrification in a low-temperature bioreactor system at two different hydraulic residence times: laboratory column studies

Nitrate removal rates in a mixture of pine woodchips and sewage sludge were determined in laboratory column studies at 5°C, 12°C, and 22°C, and at two different hydraulic residence times (HRTs; 58.2-64.0 hours and 18.7-20.6 hours). Baffles installed in the flow path were tested as a measure to reduce preferential flow behavior, and to increase the nitrate removal in the columns. The nitrate removal in the columns was simulated at 5°C and 12°C using a combined Arrhenius-Monod equation controlling the removal rate, and a first-order exchange model for incorporation of stagnant zones. Denitrification in the mixture of pine woodchips and sewage sludge reduced nitrate concentrations of 30 mg N L^-1 at 5°C to below detection limits at a HRT of 58.2-64.0 hours. At a HRT of 18.7-20.6 hours, nitrate removal was incomplete. The Arrhenius frequency factor and activation energy retrieved from the low HRT data supported a biochemically controlled reaction rate; the same parameters, however, could not be used to simulate the nitrate removal at high HRT. The results show an inversely proportional relationship between the advection velocity and the nitrate removal rate, suggesting that bioreactor performance could be enhanced by promoting low advection velocities.

Global sensitivity analysis and Bayesian parameter inference for solute transport in porous media colonized by biofilms

The concept of dual flowing continuum is a promising approach for modeling solute transport in porous media that includes biofilm phases. The highly dispersed transit time distributions often generated by these media are taken into consideration by simply stipulating that advection-dispersion transport occurs through both the porous and the biofilm phases. Both phases are coupled but assigned with contrasting hydrodynamic properties. However, the dual flowing continuum suffers from intrinsic equifinality in the sense that the outlet solute concentration can be the result of several parameter sets of the two flowing phases. To assess the applicability of the dual flowing continuum, we investigate how the model behaves with respect to its parameters. For the purpose of this study, a Global Sensitivity Analysis (GSA) and a Statistical Calibration (SC) of model parameters are performed for two transport scenarios that differ by the strength of interaction between the flowing phases. The GSA is shown to be a valuable tool to understand how the complex system behaves. The results indicate that the rate of mass transfer between the two phases is a key parameter of the model behavior and influences the identifiability of the other parameters. For weak mass exchanges, the output concentration is mainly controlled by the velocity in the porous medium and by the porosity of both flowing phases. In the case of large mass exchanges, the kinetics of this exchange also controls the output concentration. The SC results show that transport with large mass exchange between the flowing phases is more likely affected by equifinality than transport with weak exchange. The SC also indicates that weakly sensitive parameters, such as the dispersion in each phase, can be accurately identified. Removing them from calibration procedures is not recommended because it might result in biased estimations of the highly sensitive parameters.

Dissolved Organic Carbon Mobilisation in a Groundwater System Stressed by Pumping

The concentration and flux of organic carbon in aquifers is influenced by recharge and abstraction, and surface and subsurface processing. In this study groundwater was abstracted from a shallow fractured rock aquifer and dissolved organic carbon (DOC) was measured in observation bores at different distances from the abstraction bore. Groundwater abstraction at rates exceeding the aquifers yield resulted in increased DOC concentration up to 3,500 percent of initial concentrations. Potential sources of this increased DOC were determined using optical fluorescence and absorbance analysis. Groundwater fluorescent dissolved organic material (FDOM) were found to be a combination of terrestrial-derived humic material and microbial or protein sourced material. Relative molecular weight of FDOM within four metres of the abstraction well increased during the experiment, while the relative molecular weight of FDOM between four and ten metres from the abstraction well decreased. When the aquifer is not being pumped, DOC mobilisation in the aquifer is low. We hypothesise that the physical shear stress on aquifer materials caused by intense abstraction significantly increases the temporary release of DOC from sloughing of biofilms and release of otherwise bound colloidal and sedimentary organic carbon (SOC).

Nitrogen removal and spatial distribution of denitrifier and anammox communities in a bioreactor for mine drainage treatment

Mine drainage water may contain high levels of nitrate (NO3(-)) due to undetonated nitrogen-based explosives. The removal of NO3(-) and nitrite (NO2(-)) in cold climates through the microbial process of denitrification was evaluated using a pilot-scale fixed-bed bioreactor (27 m(3)). Surface water was diverted into the above-ground bioreactor filled with sawdust, crushed rock, and sewage sludge. At hydraulic residence times of ca.15 h and with the addition of acetate, NO3(-) and NO2(-) were removed to below detection levels at a NO3(-) removal rate of 5-10 g N m(-3) (bioreactor material) d(-1). The functional groups contributing to nitrogen removal in the bioreactor were studied by quantifying nirS and nirK present in denitrifying bacteria, nosZI and nosZII genes from the nitrous oxide - reducing community, and a taxa-specific part of the16S rRNA gene for the anammox community. The abundances of nirS and nirK were almost 2 orders of magnitude greater than the anammox specific 16S rRNA gene, indicating that denitrification was the main process involved in nitrogen removal. The spatial distribution of the quantified genes was heterogeneous in the bioreactor, with trends observed in gene abundance as a function of depth, distance from the bioreactor inlet, and along specific flowpaths. There was a significant relationship between the abundance of nirS, nirK, and nosZI genes and depth in the bioreactor, such that the abundance of organisms containing these genes may be controlled by oxygen diffusion and substrate supply in the partially or completely water-saturated material. Among the investigated microbial functional groups, nirS and anammox bacterial 16S rRNA genes exhibited a systematic trend of decreasing and increasing abundance, respectively, with distance from the inlet, which suggested that the functional groups respond differently to changing environmental conditions. The greater abundance of nirK along central flowpaths may indicate that the bioreactor design favored preferential flow along these flowpaths, away from the sides of the bioreactor. An improved bioreactor design should consider the role of preferential flowpaths and the heterogeneous distribution of the genetic potential for denitrification, nitrous oxide reduction and anammox on bioreactor function.

Bioclogging in porous media: influence in reduction of hydraulic conductivity and organic contaminants during synthetic leachate permeation

In this study the concept of biofilm accumulation in the sand column was promoted to assess the changes in hydraulic conductivity and concentration of organic contaminants of the synthetic leachate. Four different combinations of column study were carried out using synthetic leachate as a substrate solution. Mixed and stratified mode of experiments with two different sizes (0.3 mm and 0.6 mm) of sand grains were used for column filling. Two columns were acting as a blank, the remaining two columns amended with mixed microbial cultures which were isolated from leachate. The column was operated with continuous synthetic leachate supply for 45 days. The results indicated that the highest hydraulic conductivity reduction occurred in the mixed sand microbial column with 98.8% when compared to stratified sand microbial column. The analysis of organic contaminants of the effluent leachate was also clearly shown that the mixed sand amended with microbes poses a suitable remedial measure when compared to natural and synthetic liners for controlling the leachate migration in the subsurface environment.

Hydraulic conductivity of compacted soils controlled by microbial activity

The hydraulic conductivity defines the displacement of liquids inside porous media and affects the fate and transport of contaminants in the environment. In this research the influence of microbial growth and decay inside soil pores on hydraulic conductivity is analysed. Long-term tests performed in silt-bentonite mixtures permeated with distilled water and a nutrients solution demonstrated that hydraulic conductivity of compacted silt-bentonite samples decreases with time of permeation as a bioclogging mechanism develops. The injection of antibiotics and antifungals in the specimens produces a rebound in the hydraulic conductivity associated with the decay of microbial activity. These results show that biomediated reactions can be used to control the flow rate through compacted soil liners.

Interpreting hydrodynamic behaviour by the model of stirred tanks in series with exchanged zones: preliminary study in lab-scale trickling filters

In trickling filters for wastewater treatment, hydrodynamic behaviour is affected by the growth of biofilm on the porous medium. Therefore, modelling hydrodynamic behaviour is necessary and efficient to predict the biodegradation of pollutants. In this study, laboratory-scale trickling filters were filled with two different porous media (glass beads and plastic rings) and were fed by a synthetic substrate in batch mode. Total organic carbon (TOC) of the effluent was measured and retention time distribution (RTD) was determined by injecting NaCl. Results showed that medium had no significant effect on TOC removal rate (around 80% and 60% respectively for batch time of seven and two days). However, regarding the hydrodynamic behaviour, the effective volume ratio and hydraulic efficiency in the glass beads bed increased remarkably from 28% and 18% to 80% and 70%, respectively, with the reduction of dispersion coefficient (from 4.55 to 1.53). Moreover, the short batch time accelerated this change. Conversely, no variation of hydrodynamic behaviour in plastic rings bed was evident. Along with the feeding of synthetic substrate, biofilm concentration ranged from 1.5 to 10.1 g/L in the glass beads reactor and it achieved around 2.8 g/L in the plastic rings reactor. Hydrodynamic modelling indicated that the model of stirred tanks in series with exchanged zones fitted the experimental results well. These gave values of mobile and immobile volumes of 51 mL and 17 mL, respectively, in the glass beads filter and 25 mL and 15 mL, respectively, in the plastic rings filter.

A dual flowing continuum approach to model denitrification experiments in porous media colonized by biofilms

We present a modeling exercise of solute transport and biodegradation in a coarse porous medium widely colonized by a biofilm phase. Tracer tests in large laboratory columns using both conservative (fluorescein) and biodegradable (nitrate) solutes are simulated by means of a dual flowing continuum approach. The latter clearly distinguishes concentrations in a flowing porous phase from concentrations conveyed in the biofilm. With this conceptual setting, it becomes possible to simulate the sharp front of concentrations at early times and the flat tail of low concentrations at late times observed on the experimental breakthrough curves. Thanks to the separation of flow in two phases at different velocities, dispersion coefficients in both flowing phases keep reasonable values with some physical meaning. This is not the case with simpler models based on a single continuum (eventually concealing dead-ends), for which inferred dispersivity may reach the unphysical value of twice the size of the columns. We also show that the behavior of the dual flowing continuum is mainly controlled by the relative fractions of flow passing in each phase and the rate of mass transfer between phases. These parameters also condition the efficiency of nitrate degradation, the degradation rate in a well-seeded medium being a weakly sensitive parameter. Even though the concept of dual flowing continuum appears promising for simulating transport in complex porous media, its inversion onto experimental data really benefits from attempts with simpler models providing a rough pre-evaluation of parameters such as porosity and mean fluid velocity in the system.

Sand box experiments with bioclogging of porous media: hydraulic conductivity reductions

Tracer experiments during clogging and de-clogging experiments in a 2D sand box were via an image analysis used to establish a data set on the relation between changes in hydraulic conductivity (K) and relative porosity (β). Clogging appears to create a finger-like tracer transport, which could be caused by an initial heterogeneous distribution of biomass in the sand box. De-clogging occurs at a slower rate possibly due to the presence of inert biomass that is not affected by the starvation conditions by sudden removal of the substrate source. The tracer front was observed to get disturbed closer and closer to the substrate source during the experiments suggesting that the zone of clogging moved upstream. Three clogging models, K(β), from the literature were tested for their ability to describe the temporal changes in clogging at the scale of the sand box; the model of Clement et al. (1996) that makes no assumption on biomass distribution, the plug formation model of Thullner et al. (2002a), and the biofilm-plug formation model of Vandevivere (1995). The plug formation and biofilm-plug formation models both match the observed changes between the hydraulic conductivity of the sand box and the relative porosity. Unfortunately our experiments did not reach low relative porosities where the two models predict different behaviors. The model by Clement et al. (1996) underestimates clogging.

Physicochemical heterogeneity controls on uranium bioreduction rates at the field scale

It has been demonstrated in laboratory systems that U(VI) can be reduced to immobile U(IV) by bacteria in natural environments. The ultimate efficacy of bioreduction at the field scale, however, is often challenging to quantify and depends on site characteristics. In this work, uranium bioreduction rates at the field scale are quantified, for the first time, using an integrated approach. The approach combines field data, inverse and forward hydrological and reactive transport modeling, and quantification of reduction rates at different spatial scales. The approach is used to explore the impact of local scale (tens of centimeters) parameters and processes on field scale (tens of meters) system responses to biostimulation treatments and the controls of physicochemical heterogeneity on bioreduction rates. Using the biostimulation experiments at the Department of Energy Old Rifle site, our results show that the spatial distribution of hydraulic conductivity and solid phase mineral (Fe(III)) play a critical role in determining the field-scale bioreduction rates. Due to the dependence on Fe-reducing bacteria, field-scale U(VI) bioreduction rates were found to be largely controlled by the abundance of Fe(III) minerals at the vicinity of the injection wells and by the presence of preferential flow paths connecting injection wells to down gradient Fe(III) abundant areas.

Implications of non-equilibrium transport in heterogeneous reactive barrier systems: evidence from laboratory denitrification experiments

Organic substrates in reactive barrier systems are often heterogeneous material mixtures with relatively large contrasts in hydraulic conductivity and porosity over short distances. These short-range variations in material properties imply that preferential flow paths and diffusion between regions of higher and lower hydraulic conductivity may be important for treatment efficiency. This paper presents the results of a laboratory column experiment where denitrification is investigated using a heterogeneous reactive substrate (sawdust mixed with sewage sludge). Displacement experiments with a non-reactive solute at three different flow rates are used to estimate transport parameters using a dual porosity non-equilibrium model. Parameter estimation from breakthrough curves produced relatively consistent values for the fraction of the porosity consisting of mobile water (β) and the mass transfer coefficient (α), with average values of 0.27 and 0.42 d(-1), respectively. The column system removes >95% of the influent nitrate at low and medium flow, but only 50-75% of the influent nitrate at high flow, suggesting that denitrification kinetics and diffusive mass transfer rates are limiting the degree of treatment at lower hydraulic residence times. Reactive barrier systems containing dual porosity media must therefore consider mass transfer times in their design; this is often most easily accommodated by adjusting flowpath length.

Feedbacks between hydrological heterogeneity and bioremediation induced biogeochemical transformations

For guiding optimal design and interpretation of in situ treatments that strongly perturb subsurface systems, knowledge about the spatial and temporal patterns of mass transport and reaction intensities are important. Here, a procedure was developed and applied to time-lapse concentrations of a conservative tracer (bromide), an injected amendment (acetate) and reactive species (iron(II), uranium(VI) and sulfate) associated with two field scale biostimulation experiments, which were conducted successively at the same field location over two years. The procedure is based on a temporal moment analysis approach that relies on a streamtube approximation. The study shows that biostimulated reactions can be considerably influenced by subsurface hydrological and geochemical heterogeneities: the delivery of bromide and acetate and the intensity of the sulfate reduction is interpreted to be predominantly driven by the hydrological heterogeneity, while the intensity of the iron reduction is interpreted to be primarily controlled by the geochemical heterogeneity. The intensity of the uranium(VI) reduction appears to be impacted by both the hydrological and geochemical heterogeneity. Finally, the study documents the existence of feedbacks between hydrological heterogeneity and remediation-induced biogeochemical transformations at the field scale, particularly the development of precipitates that may cause clogging end flow rerouting.

One-dimensional model for biogeochemical interactions and permeability reduction in soils during leachate permeation

This paper uses the findings from a column study to develop a reactive model for exploring the interactions occurring in leachate-contaminated soils. The changes occurring in the concentrations of acetic acid, sulphate, suspended and attached biomass, Fe(II), Mn(II), calcium, carbonate ions, and pH in the column are assessed. The mathematical model considers geochemical equilibrium, kinetic biodegradation, precipitation-dissolution reactions, bacterial and substrate transport, and permeability reduction arising from bacterial growth and gas production. A two-step sequential operator splitting method is used to solve the coupled transport and biogeochemical reaction equations. The model gives satisfactory fits to experimental data and the simulations show that the transport of metals in soil is controlled by multiple competing biotic and abiotic reactions. These findings suggest that bioaccumulation and gas formation, compared to chemical precipitation, have a larger influence on hydraulic conductivity reduction.