Modeling Ion Transport across Thin-Film Composite Membranes During Saltwater Electrolysis.

Affordable thin-film composite (TFC) membranes are a potential alternative to more expensive ion exchange membranes in saltwater electrolyzers used for hydrogen gas production. We used a solution-friction transport model to study how the induced potential gradient controls ion transport across the polyamide (PA) active layer and support layers of TFC membranes during electrolysis. The set of parameters was simplified by assigning the same size-related partition and friction coefficients for all salt ions through the membrane active layer. The model was fit to experimental ion transport data from saltwater electrolysis with 600 mM electrolytes at a current density of 10 mA cm-2. When the electrolyte concentration and current density were increased, the transport of major charge carriers was successfully predicted by the model. Ion transport calculated using the model only minimally changed when the negative active layer charge density was varied from 0 to 600 mM, indicating active layer charge was not largely responsible for controlling ion crossover during electrolysis. Based on model simulations, a sharp pH gradient was predicted to occur within the supporting layer of the membrane. These results can help guide membrane design and operation conditions in water electrolyzers using TFC membranes.

Relative Insignificance of Polyamide Layer Selectivity for Seawater Electrolysis Applications.

Low-cost polyamide thin-film composite (TFC) membranes are being explored as alternatives to cation exchange membranes for seawater electrolysis. An optimal membrane should have a low electrical resistance to minimize applied potentials needed for water electrolysis and be able to block chloride ions present in a seawater catholyte from reaching the anode. The largest energy loss associated with a TFC membrane was the Nernstian overpotential of 0.74 V (equivalent to 37 Ω cm2 at 20 mA cm-2), derived from the pH difference between the anolyte and catholyte and not the membrane ohmic overpotential. Based on analysis using electrochemical impedance spectroscopy, the pristine TFC membrane contributed only 5.00 Ω cm2 to the ohmic resistance. Removing the polyester support layer reduced the resistance by 79% to only 1.04 Ω cm2, without altering the salt ion transport between the electrolytes. Enlarging the pore size (∼5 times) in the polyamide active layer minimally impacted counterion transport across the membrane during electrolysis, but it increased the total concentration of chloride transported by 60%. Overall, this study suggests that TFC membranes with thinner but mechanically strong supporting layers and size-selective active layers should reduce energy consumption and the potential for chlorine generation for seawater electrolyzers.

Downscaling Industrial-Scale Syngas Fermentation to Simulate Frequent and Irregular Dissolved Gas Concentration Shocks.

In large-scale syngas fermentation, strong gradients in dissolved gas (CO, H2) concentrations are very likely to occur due to locally varying mass transfer and convection rates. Using Euler-Lagrangian CFD simulations, we analyzed these gradients in an industrial-scale external-loop gas-lift reactor (EL-GLR) for a wide range of biomass concentrations, considering CO inhibition for both CO and H2 uptake. Lifeline analyses showed that micro-organisms are likely to experience frequent (5 to 30 s) oscillations in dissolved gas concentrations with one order of magnitude. From the lifeline analyses, we developed a conceptual scale-down simulator (stirred-tank reactor with varying stirrer speed) to replicate industrial-scale environmental fluctuations at bench scale. The configuration of the scale-down simulator can be adjusted to match a broad range of environmental fluctuations. Our results suggest a preference for industrial operation at high biomass concentrations, as this would strongly reduce inhibitory effects, provide operational flexibility and enhance the product yield. The peaks in dissolved gas concentration were hypothesized to increase the syngas-to-ethanol yield due to the fast uptake mechanisms in C. autoethanogenum. The proposed scale-down simulator can be used to validate such results and to obtain data for parametrizing lumped kinetic metabolic models that describe such short-term responses.

Unique stratification of biofilm density in heterotrophic membrane-aerated biofilms: An experimental and modeling study.

We consistently find a band of high cell density develop within heterotrophic membrane-aerated biofilms. This study reports and attempts to explain this unique behavior. Biofilm density affects volumetric reaction rates, biofilm growth rates, substrate diffusion, and mechanical behavior. Yet the mechanisms and dynamics of biofilm density development are poorly understood. In this study, a membrane-aerated biofilm, where O2 was supplied from the base of the biofilm and acetate from the bulk liquid, was used to explore spatial and temporal patterns of density development. Biofilm density was assessed by optical coherence tomography. After inoculation, the biofilm quickly increased in thickness, with a low density throughout. However, as the biofilm reached a stable thickness of around 1000 µm, a high-density layer developed in the biofilm interior. The layer slowly expanded over time. Oxygen microprofiles in the biofilm showed this layer coincided with the most metabolically active zone, resulting from counter-diffusing O2 and acetate. The formation of this dense layer appeared to be related to changes in growth rates. Initially, high growth rates throughout the biofilm presumably led to fast-growing, low-density biofilms. As the biofilm became thicker, and as substrates became limiting in the biofilm interior, growth rates decreased, resulting in new growth at a higher density. A 1-D mathematical model with variable biofilm density was developed by linking the rates of extracellular polymeric substances (EPS) production to the growth rate. The model captured the initial fast growth at a low density, followed by a slower, substrate-limited growth in the biofilm interior, producing a dense band within the biofilm. Together, these results suggest that low growth rates can lead to high-density zones within the interior of counter-diffusional biofilms. These findings should also be relevant to conventional, co-diffusional biofilms, although differences in density may be less obvious.

No model is perfect, but some are useful.

Agent-based model should inform the action plan to curb algal blooms in Lake Erie.

Controlling the hydraulic resistance of membrane biofilms by engineering biofilm physical structure.

The application of membrane technology for water treatment and reuse is hampered by the development of a microbial biofilm. Biofilm growth in micro-and ultrafiltration (MF/UF) membrane modules, on both the membrane surface and feed spacer, can form a secondary membrane and exert resistance to permeation and crossflow, increasing energy demand and decreasing permeate quantity and quality. In recent years, exhaustive efforts were made to understand the chemical, structural and hydraulic characteristics of membrane biofilms. In this review, we critically assess which specific structural features of membrane biofilms exert resistance to forced water passage in MF/UF membranes systems applied to water and wastewater treatment, and how biofilm physical structure can be engineered by process operation to impose less hydraulic resistance ("below-the-pain threshold"). Counter-intuitively, biofilms with greater thickness do not always cause a higher hydraulic resistance than thinner biofilms. Dense biofilms, however, had consistently higher hydraulic resistances compared to less dense biofilms. The mechanism by which density exerts hydraulic resistance is reported in the literature to be dependant on the biofilms' internal packing structure and EPS chemical composition (e.g., porosity, polymer concentration). Current reports of internal porosity in membrane biofilms are not supported by adequate experimental evidence or by a reliable methodology, limiting a unified understanding of biofilm internal structure. Identifying the dependency of hydraulic resistance on biofilm density invites efforts to control the hydraulic resistance of membrane biofilms by engineering internal biofilm structure. Regulation of biofilm internal structure is possible by alteration of key determinants such as feed water nutrient composition/concentration, hydraulic shear stress and resistance and can engineer biofilm structural development to decrease density and therein hydraulic resistance. Future efforts should seek to determine the extent to which the concept of "biofilm engineering" can be extended to other biofilm parameters such as mechanical stability and the implication for biofilm control/removal in engineered water systems (e.g., pipelines and/or, cooling towers) susceptible to biofouling.

Multiphysics modelling of photon, mass and heat transfer in coral microenvironments.

Coral reefs are constructed by calcifying coral animals that engage in a symbiosis with dinoflagellate microalgae harboured in their tissue. The symbiosis takes place in the presence of steep and dynamic gradients of light, temperature and chemical species that are affected by the structural and optical properties of the coral and their interaction with incident irradiance and water flow. Microenvironmental analyses have enabled quantification of such gradients and bulk coral tissue and skeleton optical properties, but the multi-layered nature of corals and its implications for the optical, thermal and chemical microenvironment remains to be studied in more detail. Here, we present a multiphysics modelling approach, where three-dimensional Monte Carlo simulations of the light field in a simple coral slab morphology with multiple tissue layers were used as input for modelling the heat dissipation and photosynthetic oxygen production driven by photon absorption. By coupling photon, heat and mass transfer, the model predicts light, temperature and O2 gradients in the coral tissue and skeleton, under environmental conditions simulating, for example, tissue contraction/expansion, symbiont loss via coral bleaching or different distributions of coral host pigments. The model reveals basic structure-function mechanisms that shape the microenvironment and ecophysiology of the coral symbiosis in response to environmental change.

Microbial competition reduces metabolic interaction distances to the low µm-range.

Metabolic interactions between cells affect microbial community compositions and hence their function in ecosystems. It is well-known that under competition for the exchanged metabolite, concentration gradients constrain the distances over which interactions can occur. However, interaction distances are typically quantified in two-dimensional systems or without accounting for competition or other metabolite-removal, conditions which may not very often match natural ecosystems. We here analyze the impact of cell-to-cell distance on unidirectional cross-feeding in a three-dimensional aqueous system with competition for the exchanged metabolite. Effective interaction distances were computed with a reaction-diffusion model and experimentally verified by growing a synthetic consortium of 1 µm-sized metabolite producer, receiver, and competitor cells in different spatial structures. We show that receivers cannot interact with producers located on average 15 µm away from them, as product concentration gradients flatten close to producer cells. We developed an aggregation protocol and varied the receiver cells' product affinity, to show that within producer-receiver aggregates even low-affinity receiver cells could interact with producers. These results show that competition or other metabolite-removal of a public good in a three-dimensional system reduces metabolic interaction distances to the low µm-range, highlighting the importance of concentration gradients as physical constraint for cellular interactions.

New insight into CO2-mediated denitrification process in H2-based membrane biofilm reactor: An experimental and modeling study.

The H2-based membrane biofilm reactor (H2-MBfR) is an emerging technology for removal of nitrate (NO3-) in water supplies. In this research, a lab-scale H2-MBfR equipped with a separated CO2 providing system and a microsensor measuring unit was developed for NO3- removal from synthetic groundwater. Experimental results show that efficient NO3- reduction with a flux of 1.46 g/(m2⋅d) was achieved at the optimal operating conditions of hydraulic retention time (HRT) 80 min, influent NO3- concentration 20 mg N/L, H2 pressure 5 psig and CO2 addition 50 mg/L. Given the complex counter-diffusion of substrates in the H2-MBfR, mathematical modeling is a key tool to both understand its behavior and optimize its performance. A sophisticated model was successfully established, calibrated and validated via comparing the measured and simulated system performance and/or substrate gradients within biofilm. Model results indicate that i) even under the optimal operating conditions, denitrifying bacteria (DNB) in the interior and exterior of biofilm suffered low growth rate, attributed to CO2 and H2 limitation, respectively; ii) appropriate operating parameters are essential to maintaining high activity of DNB in the biofilm; iii) CO2 concentration was the decisive factor which matters its dominant role in mediating hydrogenotrophic denitrification process; iv) the predicted optimum biofilm thickness was 650 µm that can maximize the denitrification flux and prevent loss of H2.

Noise-driven cell differentiation and the emergence of spatiotemporal patterns.

The emergence of phenotypic diversity in a population of cells and their arrangement in space and time is one of the most fascinating features of living systems. In fact, understanding multicellularity is unthinkable without explaining the proximate and the ultimate causes of cell differentiation in time and space. Simpler forms of cell differentiation can be found in unicellular organisms, such as bacterial biofilm, where reversible cell differentiation results in phenotypically diverse populations. In this manuscript, we attempt to start with the simple case of reversible nongenetic phenotypic to construct a model of differentiation and pattern formation. Our model, which we refer to as noise-driven differentiation (NDD) model, is an attempt to consider the prevalence of noise in biological systems, alongside what is known about genetic switches and signaling, to create a simple model which generates spatiotemporal patterns from bottom-up. Our simulations indicate that the presence of noise in cells can lead to reversible differentiation and the addition of signaling can create spatiotemporal pattern.

Computational and Experimental Investigation of Biofilm Disruption Dynamics Induced by High-Velocity Gas Jet Impingement.

Experimental data showed that high-speed microsprays can effectively disrupt biofilms on their support substratum, producing a variety of dynamic reactions such as elongation, displacement, ripple formation, and fluidization. However, the mechanics underlying the impact of high-speed turbulent flows on biofilm structure is complex under such extreme conditions, since direct measurements of viscosity at these high shear rates are not possible using dynamic testing instruments. Here, we used computational fluid dynamics simulations to assess the complex fluid interactions of ripple patterning produced by high-speed turbulent air jets impacting perpendicular to the surface of Streptococcus mutans biofilms, a dental pathogen causing caries, captured by high-speed imaging. The numerical model involved a two-phase flow of air over a non-Newtonian biofilm, whose viscosity as a function of shear rate was estimated using the Herschel-Bulkley model. The simulation suggested that inertial, shear, and interfacial tension forces governed biofilm disruption by the air jet. Additionally, the high shear rates generated by the jet impacts coupled with shear-thinning biofilm property resulted in rapid liquefaction (within milliseconds) of the biofilm, followed by surface instability and traveling waves from the impact site. Our findings suggest that rapid shear thinning under very high shear flows causes the biofilm to behave like a fluid and elasticity can be neglected. A parametric sensitivity study confirmed that both applied force intensity (i.e., high jet nozzle air velocity) and biofilm properties (i.e., low viscosity and low air-biofilm surface tension and thickness) intensify biofilm disruption by generating large interfacial instabilities.IMPORTANCE Knowledge of mechanisms promoting disruption though mechanical forces is essential in optimizing biofilm control strategies which rely on fluid shear. Our results provide insight into how biofilm disruption dynamics is governed by applied forces and fluid properties, revealing a mechanism for ripple formation and fluid-biofilm mixing. These findings have important implications for the rational design of new biofilm cleaning strategies with fluid jets, such as determining optimal parameters (e.g., jet velocity and position) to remove the biofilm from a certain zone (e.g., in dental hygiene or debridement of surgical site infections) or using antimicrobial agents which could increase the interfacial area available for exchange, as well as causing internal mixing within the biofilm matrix, thus disrupting the localized microenvironment which is associated with antimicrobial tolerance. The developed model also has potential application in predicting drag and pressure drop caused by biofilms on bioreactor, pipeline, and ship hull surfaces.

Biofilm compressibility in ultrafiltration: A relation between biofilm morphology, mechanics and hydraulic resistance.

Poroelastic fluid-structure interaction models were coupled to experimental data to determine the effects of biofilm spatial distribution of mechanical and hydraulic properties on the biofilm hydraulic resistance and compressibility in membrane filtration processes. Biofilms were cultivated on ultrafiltration membranes for 20 and 30 days under high (0.28 bar) and low (0.06 bar) transmembrane pressure (TMP), in dead-end filtration mode. Subsequently, biofilms were subjected to a compression/relaxation cycles by step-wise TMP changes. Structural deformation of biofilms during compression was observed in-situ using optical coherence tomography. Experimental results show that the observed increase in the biofilm hydraulic resistance during compression is not necessarily accompanied by a detectable biofilm thickness reduction. A dual-layer biofilm model with a dense base and porous top layer could explain these observed results. Because porosity controls indirectly the mechanical response of biofilms under compression, results could be described without assuming a gradient in mechanical properties within the biofilm. The biofilm surface roughness did not significantly influence the water flux in this study. However, the fraction of biofilm base layer directly exposed to bulk liquid could be a good indicator in the determination of water flux. The main implications of this study for the design and operation of low-pressure membrane systems (e.g., MF and UF with fouling layer being the main filtration resistance) lays in the selection of favorable operational TMP and biofilm morphology.

Determination of mechanical properties of biofilms by modelling the deformation measured using optical coherence tomography.

The advantage of using non-invasive imaging such as optical coherence tomography (OCT) to asses material properties from deformed biofilm geometries can be compromised by the assumptions made on fluid forces acting on the biofilm. This study developed a method for the determination of elastic properties of biofilms by modelling the biofilm deformation recorded by OCT imaging with poroelastic fluid-structure interaction computations. Two-dimensional biofilm geometries were extracted from OCT scans of non-deformed and deformed structures as a result of hydrodynamic loading. The biofilm geometries were implemented in a model coupling fluid dynamics with elastic solid mechanics and Darcy flow in the biofilm. The simulation results were compared with real deformed geometries and a fitting procedure allowed estimation of the Young's modulus in given flow conditions. The present method considerably improves the estimation of elastic moduli of biofilms grown in mini-fluidic rectangular channels. This superior prediction is based on the relaxation of several simplifying assumptions made in past studies: shear stress is not anymore taken constant over the biofilm surface, total stress including also pressure is accounted for, any biofilm shape can be used in the determinations, and non-linear behavior of mechanical properties can be estimated. Biofilm elastic moduli between 70 and 700 Pa were obtained and biofilm hardening at large applied stress due to increasing flow velocity was quantified. The work performed here opens the way for in-situ determination of other mechanical properties (e.g., viscoelastic properties, relaxation times, plastic yields) and provides data for modelling biofilm deformation and detachment with eventual applications in biofilm control and removal strategies.

Effect of biofilm structural deformation on hydraulic resistance during ultrafiltration: A numerical and experimental study.

Biofilm formation in membrane systems negatively impacts the filtration system performances. This study evaluated how biofilm compression driven by permeate flow increases the hydraulic resistance and leads to reduction in permeate flux. We analysed the effect of biofilm compression on hydraulic resistance and permeate flux through computational models supported by experimental data. Biofilms with homogeneous surface structure were subjected to step-wise changes in flux and transmembrane pressure during compression and relaxation tests. Biofilm thickness under applied forces was measured non-invasively in-situ using optical coherence tomography (OCT). A numerical model of poroelasticity, which couples water flow through the biofilm with biofilm mechanics, was developed to correlate the structural deformation with biofilm hydraulics (permeability and resistance). The computational model enabled extracting mechanical and hydrological parameters corresponding to the experimental data. Homogeneous biofilms under elevated compression forces experienced a significant reduction in thickness while only a slight increase in resistance was observed. This shows that hydraulic resistance of homogeneous biofilms was affected more by permeability decrease due to pore closure than by a decrease in thickness. Both viscoelastic and elastoplastic models could describe well the permanent biofilm deformation. However, for biofilms under study, a simpler elastic model could also be used due to the small irreversible deformations. The elastic moduli fitting the measured data were in agreement with other reported values for biofilm under compression. Biofilm stiffening under larger flow-driven compression forces was observed and described numerically by correlating inversely the elastic modulus with biofilm porosity. The importance of this newly developed method lies in estimation of accurate biofilm mechanical parameters to be used in numerical models for both membrane filtration system and biofouling cleaning strategies. Such model can ultimately be used to identify optimal operating conditions for membrane system subjected to biofouling.

Quantitative chemical biosensing by bacterial chemotaxis in microfluidic chips.

Whole-cell bacterial bioreporters are proposed as alternatives to chemical analysis of, for example, pollutants in environmental compartments. Commonly based on reporter gene induction, bioreporters produce a detectable signal within 30 min to a few hours after exposure to the chemical target, which is impractical for applications aiming at a fast response. In an attempt to attain faster readout but maintain flexibility of chemical targeting, we explored the concept for quantitative chemical sensing by bacterial chemotaxis. Chemotaxis was quantified from enrichment of cells across a 600 µm-wide chemical gradient stabilized by parallel flow in a microfluidic chip, further supported by transport and chemotaxis steady state and kinetic modelling. As proof-of-concept, we quantified Escherichia coli chemotaxis towards serine, aspartate and methylaspartate as a function of attractant concentration and exposure time. E. coli chemotaxis enrichment increased sharply between 0 and 10 µM serine, before saturating at 100 µM. The chemotaxis accumulation rate was maximal at 10 µM serine, leading to observable cell enrichment within 5 min. The potential application for biosensing of environmental toxicants was investigated by quantifying chemotaxis of Cupriavidus pinatubonensis JMP134 towards the herbicide 2,4-dichlorophenoxyacetate. Our results show that bacterial chemotaxis can be quantified on a scale of minutes and may be used for developing faster bioreporter assays.

Mechanisms of nitrous oxide (N2 O) formation and reduction in denitrifying biofilms.

Nitrous oxide (N2 O) is a potent greenhouse gas that can be formed in wastewater treatment processes by ammonium oxidizing and denitrifying microorganisms. While N2 O emissions from suspended growth systems have been extensively studied, and some recent studies have addressed emissions from nitrifying biofilms, much less is known about N2 O emissions from denitrifying biofilm processes. This research used modeling to evaluate the mechanisms of N2 O formation and reduction in denitrifying biofilms. The kinetic model included formation and consumption of key denitrification species, including nitrate (NO3-), nitrite (NO2-), nitric oxide (NO), and N2 O. The model showed that, in presence of excess of electron donor, denitrifying biofilms have two distinct layers of activity: an outer layer where there is net production of N2 O and an inner layer where there is net consumption. The presence of oxygen (O2 ) had an important effect on N2 O emission from suspended growth systems, but a smaller effect on biofilm systems. The effects of NO3- and O2 differed significantly based on the biofilm thickness. Overall, the effects of biofilm thickness and bulk substrate concentrations on N2 O emissions are complex and not always intuitive. A key mechanism for denitrifying biofilms is the diffusion of N2 O and other intermediates from one zone of the biofilm to another. This leads to zones of N2 O formation or consumption transformations that would not exist in suspended growth systems.

Predicting N2O emissions from nitrifying and denitrifying biofilms: a modeling study.

Wastewater treatment plants can be significant sources of nitrous oxide (N2O), a potent greenhouse gas. While our understanding of N2O emissions from suspended-growth processes has advanced significantly, less is known about emissions from biofilm processes. Biofilms may behave differently due to their substrate gradients and microbial stratification. In this study, we used mathematical modeling to explore the mechanisms of N2O emissions from nitrifying and denitrifying biofilms. Our ammonia-oxidizing bacteria biofilm model suggests that N2O emissions from biofilm can be significantly greater than from suspended-growth systems. The driving factor is the diffusion of hydroxylamine, a nitrification intermediate, from the aerobic to the anoxic regions of the biofilm. The presence of nitrite-oxidizing bacteria further increased emissions. For denitrifying biofilms, our results suggest that emissions are generally greater than for suspended-growth systems. However, the magnitude of the difference depends on the bulk dissolved oxygen, chemical oxygen demand, and nitrate concentrations, as well as the biofilm thickness. Overall, the accumulation and diffusion of key intermediates, i.e. hydroxylamine and nitrite, distinguish biofilms from suspended-growth systems. Our research suggests that the mechanisms of N2O emissions from biofilms are much more complex than suspended-growth systems, and that emissions may be higher in many cases.

Impact of cell cluster size on apparent half-saturation coefficients for oxygen in nitrifying sludge and biofilms.

A three-dimensional (3-D) diffusion-reaction model was used to assess the effects of nitrifiers growing in cell clusters on the apparent oxygen half-saturation coefficients in activated sludge flocs. The model allows conciliation of seemingly contradictory reports by several research groups. Although intrinsic half-saturation coefficients (i.e., not affected by diffusion) show a better affinity for oxygen for ammonia oxidizing (AOB) than for nitrite oxidizing bacteria (NOB) (KO,AOB < KO,NOB), measurements in flocs often produced reversed apparent values (KO,AOB,app > KO,NOB,app), which can now be explained by the 3-D model with AOB and NOB microcolonies. This effect cannot be described with a conventional 1-D homogeneous model because the reversion of the AOB/NOB apparent KO is caused by the high biomass density and resulting concentration gradients inside the microcolonies. Two main factors explain the reversion of the half-saturation coefficients: the difference in oxygen yields (for NOB lower than for AOB) and the difference in colony size (NOB colonies are smaller than those of AOB). The strongest increase in the apparent half-saturation coefficients is linked to the colony size, rather than to the floc size. For high-density microbial aggregates (i.e., granular sludge), the need for a stratified population (AOB outer shell, NOB inner layers) was revealed in order to outcompete NOB. This study stresses the need for a more detailed description of the biomass distribution in activated sludge, granular sludge and biofilm reactors when elucidating the mechanisms for NOB repression.

Biofouling in forward osmosis systems: An experimental and numerical study.

This study evaluates with numerical simulations supported by experimental data the impact of biofouling on membrane performance in a cross-flow forward osmosis (FO) system. The two-dimensional numerical model couples liquid flow with solute transport in the FO feed and draw channels, in the FO membrane support layer and in the biofilm developed on one or both sides of the membrane. The developed model was tested against experimental measurements at various osmotic pressure differences and in batch operation without and with the presence of biofilm on the membrane active layer. Numerical studies explored the effect of biofilm properties (thickness, hydraulic permeability and porosity), biofilm membrane surface coverage, and biofilm location on salt external concentration polarization and on the permeation flux. The numerical simulations revealed that (i) when biofouling occurs, external concentration polarization became important, (ii) the biofilm hydraulic permeability and membrane surface coverage have the highest impact on water flux, and (iii) the biofilm formed in the draw channel impacts the process performance more than when formed in the feed channel. The proposed mathematical model helps to understand the impact of biofouling in FO membrane systems and to develop possible strategies to reduce and control biofouling.

Challenges in microbial ecology: building predictive understanding of community function and dynamics.

The importance of microbial communities (MCs) cannot be overstated. MCs underpin the biogeochemical cycles of the earth's soil, oceans and the atmosphere, and perform ecosystem functions that impact plants, animals and humans. Yet our ability to predict and manage the function of these highly complex, dynamically changing communities is limited. Building predictive models that link MC composition to function is a key emerging challenge in microbial ecology. Here, we argue that addressing this challenge requires close coordination of experimental data collection and method development with mathematical model building. We discuss specific examples where model-experiment integration has already resulted in important insights into MC function and structure. We also highlight key research questions that still demand better integration of experiments and models. We argue that such integration is needed to achieve significant progress in our understanding of MC dynamics and function, and we make specific practical suggestions as to how this could be achieved.