• Energy Research

AbstractSubaerial biofilms (SAB) are intricate microbial communities living on terrestrial surfaces, of interest in a variety of contexts including cultural heritage preservation, microbial ecology, biogeochemical cycling, and biotechnology. Here we propose a mathematical model aimed at better understanding the interplay between cyanobacteria and heterotrophic bacteria, common microbial SAB constituents, and their mutual dependence on local environmental conditions. SABs are modeled as thin mixed biofilm-liquid water layers sitting on stone. A system of ordinary differential equations regulates the dynamics of key SAB components: cyanobacteria, heterotrophs, polysaccharides and decayed biomass, as well as cellular levels of organic carbon, nitrogen and energy. These components are interconnected through a network of energetically dominant metabolic pathways, modeled with limitation terms reflecting the impact of biotic and abiotic factors. Daily cylces of temperature, humidity, and light intensity are considered as input model variables that regulate microbial activity by influencing water availability and metabolic kinetics. Relevant physico-chemical processes, including pH regulation, further contribute to a description of the SAB ecology. Numerical simulations explore the dynamics of SABs in a real-world context, revealing distinct daily activity periods shaped by water activity and light availability, as well as longer time scale survivability conditions. Results also suggest that heterotrophs could play a substantial role in decomposing non-volatile carbon compounds and regulating pH, thus influencing the overall composition and stability of the biofilm.

AbstractA multiscale mathematical model is presented to describe de novo granulation, and the evolution of multispecies granular biofilms, in a continuously fed bioreactor. The granule is modelled as a spherical free boundary domain with radial symmetry. The equation governing the free boundary is derived from global mass balance considerations and takes into account the growth of sessile biomass as well as exchange fluxes with the bulk liquid. Starting from a vanishing initial value, the expansion of the free boundary is initiated by the attachment process, which depends on the microbial species concentrations within the bulk liquid and their specific attachment velocity. Nonlinear hyperbolic PDEs model the growth of the sessile microbial species, while quasi-linear parabolic PDEs govern the dynamics of substrates and invading species within the granular biofilm. Nonlinear ODEs govern the evolution of soluble substrates and planktonic biomass within the bulk liquid. The model is applied to an anaerobic, granular-based bioreactor system, and solved numerically to test its qualitative behaviour and explore the main aspects of de novo anaerobic granulation: ecology, biomass distribution, relative abundance, dimensional evolution of the granules and soluble substrates, and planktonic biomass dynamics within the bioreactor. The numerical results confirm that the model accurately describes the ecology and the concentrically layered structure of anaerobic granules observed experimentally, and that it can predict the effects on the process of significant factors, such as influent wastewater composition; granulation properties of planktonic biomass; biomass density; detachment intensity; and number of granules.