Martin Thullner

Influence of Microbial Growth on Hydraulic Properties of Pore Networks

From laboratory experiments it is known that bacterial biomass is able to influence the hydraulic properties of saturated porous media, an effect called bioclogging. To interpret the observations of these experiments and to predict possible bioclogging effects on the field scale it is necessary to use transport models, which are able to include bioclogging. For these models it is necessary to know the relation between the amount of biomass and the hydraulic conductivity of the porous medium. Usually these relations were determined using bundles of parallel pore channels and do not account for interconnections between the pores in more than one dimension. The present study uses two-dimensional pore network models to study the effects of bioclogging on the pore scale. Numerical simulations were done for two different scenarios of the growth of biomass in the pores. Scenario 1 assumes microbial growth in discrete colonies clogging particular pores completely. Scenario 2 assumes microbial growth as a biofilm growing on the wall of each pore. In both scenarios the hydraulic conductivity was reduced by at least two orders of magnitude, but for the colony scenario much less biomass was needed to get a maximal clogging effect and a better agreement with previously published experimental data could be found. For both scenarios it was shown that heterogeneous pore networks could be clogged with less biomass than more homogeneous ones.

Interaction between water flow and spatial distribution of microbial growth in a two-dimensional flow field in saturated porous media.

Bacterial growth and its interaction with water flow was investigated in a two-dimensional flow field in a saturated porous medium. A flow cell (56 x 44 x 1 cm) was filled with glass beads and operated under a continuous flow of a mineral medium containing nitrate as electron acceptor. A glucose solution was injected through an injection port, simulating a point source contamination. Visible light transmission was used to observe the distribution of the growing biomass and water flow during the experiment. At the end of the experiment (on day 31), porous medium samples were destructively collected and analyzed for abundance of total and active bacterial cells, bacterial cell volume and concentration of polysaccharides and proteins. Microbial growth was observed in two stripes along the length of the flow cell, starting at the glucose injection port, where highest biomass concentrations were obtained. The spatial distribution of biomass indicated that microbial activity was limited by transverse mixing between glucose and nitrate media, as only in the mixing zone between the media high biological activities were achieved. The ability of the biomass to change the flow pattern in the flow cell was observed, indicating that the biomass was locally reducing the hydraulic conductivity of the porous medium. This bioclogging effect became evident when the injection of the glucose solution was turned off and water flow still bypassed the area around the glucose injection port, preserving the flow pattern as it was during the injection of the glucose solution. As flow bypass was possible in this system, the average hydraulic properties of the flow cell were not affected by the produced biomass. Even in the vicinity of the injection port, the total volume of the bacterial cells remained below 0.01% of the pore space and was unlikely to be responsible for the bioclogging. However, the bacteria produced large amounts of extracellular polymeric substances (EPS), which likely caused the observed bioclogging effects.

Modeling Microbially Induced Carbon Degradation in Redox-Stratified Subsurface Environments: Concepts and Open Questions

The degradation of organic matter, including organic contaminants, in subsurface environments is controlled by the abundances and functional capabilities of the resident microorganisms. As a consequence, modeling approaches simulating the fate of organics and related changes in redox conditions have to account for the effects of microbial activity on the degradation kinetics, as well as for the spatial and temporal distributions of the chemical species (e.g., terminal electron acceptors, nutrients or toxic substances) that control microbial activity. The present paper reviews the principal modeling approaches that are used to simulate the degradation of organic matter in water-saturated porous media. Special attention is devoted to modeling the bioavailability of chemical substrates of microbial reactions, and the sequential occurrence of terminal electron accepting pathways. While the various model approaches found in the literature are capable of reproducing field data sets from various environmental settings, they are rarely compared in terms of performance and predictive ability. Most approaches incorporate simplifications or empirical rate laws, which limit their range of application. Thus, there remains a need for further development of more general, process-based modeling concepts to represent microbially mediated reactive processes.

Global-scale quantification of mineralization pathways in marine sediments: A reaction-transport modeling approach

[1] The global-scale quantification of organic carbon (Corg) degradation pathways in marine sediments is difficult to achieve experimentally due to the limited availability of field data. In the present study, a numerical modeling approach is used as an alternative to quantify the major metabolic pathways of Corg oxidation (Cox) and associated fluxes of redox-sensitive species fluxes along a global ocean hypsometry, using the seafloor depth (SFD) as the master variable. The SFD dependency of the model parameters and forcing functions is extracted from existing empirical relationships or from the NOAA World Ocean Atlas. Results are in general agreement with estimates from the literature showing that the relative contribution of aerobic respiration to Cox increases from 80% in deep-sea sediments. Sulfate reduction essentially follows an inversed SFD dependency, the other metabolic pathways (denitrification, Mn and Fe reduction) only adding minor contributions to the global-scale mineralization of Corg. The hypsometric analysis allows the establishment of relationships between the individual terminal electron acceptor (TEA) fluxes across the sediment-water interface and their respective contributions to the Corg decomposition process. On a global average, simulation results indicate that sulfate reduction is the dominant metabolic pathway and accounts for approximately 76% of the total Cox, which is higher than reported so far by other authors. The results also demonstrate the importance of bioirrigation for the assessment of global species fluxes. Especially at shallow SFD most of the TEAs enter the sediments via bioirrigation, which complicates the use of concentration profiles for the determination of total TEA fluxes by molecular diffusion. Furthermore, bioirrigation accounts for major losses of reduced species from the sediment to the water column prohibiting their reoxidation inside the sediment. As a result, the total carbon mineralization rate exceeds the total flux of oxygen into the sediment by a factor of 2 globally.

GeoSysBRNS-A flexible multidimensional reactive transport model for simulating biogeochemical subsurface processes

The description of reactive transport processes in subsurface environments requires a sound understanding of both the biogeochemical complexity of the system and the spatially resolved transport of reactive species. However, most existing reactive transport models, for example in the field of contaminant hydrology, are specialized either in the simulation of the reactive or of the flow and transport processes. In this paper, we present and test the coupling of two highly flexible codes for the simulation of reactive transport processes in the subsurface: the Biogeochemical Reaction Network Simulator (BRNS), which contains a solver for kinetically and thermodynamically constrained biogeochemical reactions, and GeoSys/RockFlow, a multidimensional finite element subsurface flow and transport simulator. The new model, named GeoSysBRNS, maintains the full flexibility of the original models. The coupling is handled using an operator splitting scheme, which allows the reactive solver to be compiled into a problem specific library that is accessed by the transport simulator at runtime. The accuracy of the code coupling within GeoSysBRNS is demonstrated using two benchmark problems from the literature: a laboratory experiment on organic carbon degradation in a sand column via multiple microbial degradation pathways, and a dispersive mixing controlled bioreactive transport problem in aquifers, assuming three different reaction kinetics.

Assimilation of toluene carbon along a bacteria–protist food chain determined by 13C-enrichment of biomarker fatty acids

A food chain consisting of toluene, toluene-degrading Pseudomonas sp. PS+ and a bacterivorous flagellated amoebae Vahlkampfia sp. was established in a batch culture. This culture was amended with [U-13C]toluene and served as a model system to elucidate the flux of carbon in the food chain by quantifying bacterial biovolumes and 13C enrichment of phospholipid fatty acid (PLFA) biomarkers of the bacteria and the heterotrophic protists. Major PLFA detected in the batch co-culture included those derived from Pseudomonas sp. PS+ (16:1omega7c and 18:1omega7c) and Vahlkampfia sp. (20:4omega6c and 20:3omega6c). A numerical model including consumption of toluene by the bacteria and predation of the bacteria by the heterotrophic protists was adjusted to the measured toluene carbon, bacterial carbon and delta13C values of bacterial and protist biomass. Using this model, we estimated that 28+/-7% of the consumed toluene carbon was transformed into bacterial biomass, and 12+/-4% of the predated bacterial carbon was incorporated into heterotrophic protist biomass. Our study showed that the 13C enrichment of PLFA biomarkers coupled to biomass determination via biovolume calculations is a suitable method to trace carbon fluxes in protist-inclusive microbial food chains because it does not require the separation of protist cells from bacterial cells and soil particles.

Remediation of groundwater contaminated with MTBE and benzene: the potential of vertical-flow soil filter systems

Field investigations on the treatment of MTBE and benzene from contaminated groundwater in pilot or full-scale constructed wetlands are lacking hugely. The aim of this study was to develop a biological treatment technology that can be operated in an economic, reliable and robust mode over a long period of time. Two pilot-scale vertical-flow soil filter eco-technologies, a roughing filter (RF) and a polishing filter (PF) with plants (willows), were operated independently in a single-stage configuration and coupled together in a multi-stage (RF+PF) configuration to investigate the MTBE and benzene removal performances. Both filters were loaded with groundwater from a refinery site contaminated with MTBE and benzene as the main contaminants, with a mean concentration of 2970±816 and 13,966±1998 μg L(-1), respectively. Four different hydraulic loading rates (HLRs) with a stepwise increment of 60, 120, 240 and 480 L m(-2) d(-1) were applied over a period of 388 days in the single-stage operation. At the highest HLR of 480 L m(-2) d(-1), the mean concentrations of MTBE and benzene were found to be 550±133 and 65±123 μg L(-1) in the effluent of the RF. In the effluent of the PF system, respective mean MTBE and benzene concentrations of 49±77 and 0.5±0.2 μg L(-1) were obtained, which were well below the relevant MTBE and benzene limit values of 200 and 1 μg L(-1) for drinking water quality. But a dynamic fluctuation in the effluent MTBE concentration showed a lack of stability in regards to the increase in the measured values by nearly 10%, which were higher than the limit value. Therefore, both (RF+PF) filters were combined in a multi-stage configuration and the combined system proved to be more stable and effective with a highly efficient reduction of the MTBE and benzene concentrations in the effluent. Nearly 70% of MTBE and 98% of benzene were eliminated from the influent groundwater by the first vertical filter (RF) and the remaining amount was almost completely diminished (∼100% reduction) after passing through the second filter (PF), with a mean MTBE and benzene concentration of 5±10 and 0.6±0.2 μg L(-1) in the final effluent. The emission rate of volatile organic compounds mass into the air from the systems was less than 1% of the inflow mass loading rate. The results obtained in this study not only demonstrate the feasibility of vertical-flow soil filter systems for treating groundwater contaminated with MTBE and benzene, but can also be considered a major step forward towards their application under full-scale conditions for commercial purposes in the oil and gas industries.

A reactive transport modeling approach to simulate biogeochemical processes in pore structures with pore-scale heterogeneities

Redox processes, including degradation of organic contaminants, are often controlled by microorganisms residing in natural porous media like soils or aquifers. These environments are characterized by heterogeneities at various scales which influence the transport of chemical species and the spatial distribution of microorganisms. As a result, the accessibility of the chemical species by the resident microbial populations may be limited, altering the efficiency of the biodegradation process. Hence, the biodegradation rate of contaminants at large scales does not only depend on the degradation capacity of the indigenous microbial population but also on the heterogeneities of the hosting media at smaller scales. It is thus important to establish a link between effective reaction rates and various structural features of porous media which can be directly observed or measured. This link is urgently needed because explicit resolution of heterogeneities within large-scale reactive transport models is still limited by the available computational capacities. The present study introduces a reactive transport modeling approach to determine the influence of pore-scale heterogeneities on biogeochemical processes in porous media. For this purpose, a pore network model, which simulates flow and advective-diffusive transport of chemical species in heterogeneous pore networks is developed and coupled to the Biogeochemical Reaction Network Simulator (BRNS). The resulting coupled model (PNBRNS) is able to simulate the reactive transport of solutes in heterogeneous pore assemblies. The PNBRNS model is applied for the simulation of a test case of bioavailability and effective biodegradation rate of a dissolved contaminant in different pore networks, built using a discrete set of geostatistically derived pore-size or biomass distributions. Results show that the heterogeneity of the pore-size distribution has a significant impact on bioavailability while the heterogeneity of the biomass distribution only leads to minor effects. The model also includes intra-pore bioavailability restrictions using diffusion-limited biodegradation kinetics. The results indicate that intra-pore limitations lead to extra constrains on the biodegradation of contaminants, even in the presence of larger-scale structural heterogeneities.

Dissolved Organic Carbon Enhances the Mass Transfer of Hydrophobic Organic Compounds from Nonaqueous Phase Liquids (NAPLs) into the Aqueous Phase

The hypothesis that dissolved organic carbon (DOC) enhances the mass transfer of hydrophobic organic compounds from nonaqueous phase liquids (NAPLs) into the aqueous phase above that attributable to dissolved molecular diffusion alone was tested. In controlled experiments, mass transfer rates of five NAPL-phase PAHs (log K(OW) 4.15-5.39) into the aqueous phase containing different concentrations of DOC were measured. Mass transfer rates were increased by up to a factor of 4 in the presence of DOC, with the greatest enhancement being observed for more hydrophobic compounds and highest DOC concentrations. These increases could not be explained by dissolved molecular diffusion alone, and point to a parallel DOC-mediated diffusive pathway. The nature of the DOC-mediated diffusion pathway as a function of the DOC concentration and PAH sorption behavior to the DOC was investigated using diffusion-based models. The DOC-enhanced mass transfer of NAPL-phase hydrophobic compounds into the aqueous phase has important implications for their bioremediation as well as bioconcentration and toxicity.

Influence of mass transfer on stable isotope fractionation

Biodegradation of contaminants is a common remediation strategy for subsurface environments. To monitor the success of such remediation means a quantitative assessment of biodegradation at the field scale is required. Nevertheless, the reliable quantification of the in situ biodegradation process it is still a major challenge. Compound-specific stable isotope analysis has become an established method for the qualitative analysis of biodegradation in the field and this method is also proposed for a quantitative analysis. However, to use stable isotope data to obtain quantitative information on in situ biodegradation requires among others knowledge on the influence of mass transfer processes on the observed stable isotope fractionation. This paper reviews recent findings on the influence of mass transfer processes on stable isotope fractionation and on the quantitative interpretation of isotope data. Focus will be given on small-scale mass transfer processes controlling the bioavailability of contaminants. Such bioavailability limitations are known to affect the biodegradation rate and have recently been shown to affect stable isotope fractionation, too. Theoretical as well as experimental studies addressing the link between bioavailability and stable isotope fractionation are reviewed and the implications for assessing biodegradation in the field are discussed.