How flocculation can explain coexistence in the chemostat

Local intraspecific aggregation in phytoplankton model communities: spatial scales of occurrence and implications for coexistence

The coexistence of multiple phytoplankton species despite their reliance on similar resources is often explained with mean-field models assuming mixed populations. In reality, observations of phytoplankton indicate spatial aggregation at all scales, including at the scale of a few individuals. Local spatial aggregation can hinder competitive exclusion since individuals then interact mostly with other individuals of their own species, rather than competitors from different species. To evaluate how microscale spatial aggregation might explain phytoplankton diversity maintenance, an individual-based, multispecies representation of cells in a hydrodynamic environment is required. We formulate a three-dimensional and multispecies individual-based model of phytoplankton population dynamics at the Kolmogorov scale. The model is studied through both simulations and the derivation of spatial moment equations, in connection with point process theory. The spatial moment equations show a good match between theory and simulations. We parameterized the model based on phytoplankters' ecological and physical characteristics, for both large and small phytoplankton. Defining a zone of potential interactions as the overlap between nutrient depletion volumes, we show that local species composition-within the range of possible interactions-depends on the size class of phytoplankton. In small phytoplankton, individuals remain in mostly monospecific clusters. Spatial structure therefore favours intra- over inter-specific interactions for small phytoplankton, contributing to coexistence. Large phytoplankton cell neighbourhoods appear more mixed. Although some small-scale self-organizing spatial structure remains and could influence coexistence mechanisms, other factors may need to be explored to explain diversity maintenance in large phytoplankton.

Global stability of a continuous bioreactor model under persistent variation of the dilution rate

In this work, the global stability of a continuous bioreactor model is studied, with the concentrations of biomass and substrate as state variables, a general non-monotonic function of substrate concentration for the specific growth rate, and constant inlet substrate concentration. Also, the dilution rate is time varying but bounded, thus leading to state convergence to a compact set instead of an equilibrium point. Based on the Lyapunov function theory with dead-zone modification, the convergence of the substrate and biomass concentrations is studied. The main contributions with respect to closely related studies are: i) The convergence regions of the substrate and biomass concentrations are determined as function of the variation region of the dilution rate (D) and the global convergence to these compact sets is proved, considering monotonic and non-monotonic growth functions separately; ii) several improvements are proposed in the stability analysis, including the definition of a new dead zone Lyapunov function and the properties of its gradient. These improvements allow proving convergence of substrate and biomass concentrations to their compact sets, while tackling the interwoven and nonlinear nature of the dynamics of biomass and substrate concentrations, the non-monotonic nature of the specific growth rate, and the time-varying nature of the dilution rate. The proposed modifications are a basis for further global stability analysis of bioreactor models exhibiting convergence to a compact set instead of an equilibrium point. Finally, the theoretical results are illustrated through numerical simulation, showing the convergence of the states under varying dilution rate.

Interspecific density-dependent model of predator–prey relationship in the chemostat

The objective of this study is to analyze a model of competition for one resource in the chemostat with general interspecific density-dependent growth rates, taking into account the predator–prey relationship. This relationship is characterized by the fact that the prey species promotes the growth of the predator species which in turn inhibits the growth of the first species. The model is a three-dimensional system of ordinary differential equations. With the same dilution rates, the model can be reduced to a planar system where the two models have the same local and even global behavior. The existence and stability conditions of all steady states of the reduced model in the plane are determined according to the operating parameters. Using the nullcline method, we present a geometric characterization of the existence and stability of all equilibria showing the multiplicity of coexistence steady states. The bifurcation diagrams illustrate that the steady states can appear or disappear only through saddle-node or transcritical bifurcations. Moreover, the operating diagrams describe the asymptotic behavior of this system by varying the control parameters and show the effect of the inhibition of predation on the emergence of the bistability region and the reduction until the disappearance of the coexistence region by increasing this inhibition parameter.

A density-dependent model of competition for one resource in the chemostat

This paper deals with a two-microbial species model in competition for a single-resource in the chemostat including general intra- and interspecific density-dependent growth rates with distinct removal rates for each species. In order to understand the effects of intra- and interspecific interference, this general model is first studied by determining the conditions of existence and local stability of steady states. With the same removal rate, the model can be reduced to a planar system and then the global stability results for each steady state are derived. The bifurcations of steady states according to interspecific interference parameters are analyzed in a particular case of density-dependent growth rates which are usually used in the literature. The operating diagrams show how the model behaves by varying the operating parameters and illustrate the effect of the intra- and interspecific interference on the disappearance of coexistence region and the occurrence of bi-stability region. Concerning the small enough interspecific interference terms, we would shed light on the global convergence towards the coexistence steady state for any positive initial condition. When the interspecific interference pressure is large enough this system exhibits bi-stability where the issue of the competition depends on the initial condition.

Competition for a single resource and coexistence of several species in the chemostat

We study a model of the chemostat with several species in competition for a single resource. We take into account the intra-specific interactions between individuals of the same population of micro-organisms and we assume that the growth rates are increasing and the dilution rates are distinct. Using the concept of steady-state characteristics, we present a geometric characterization of the existence and stability of all equilibria. Moreover, we provide necessary and sufficient conditions on the control parameters of the system to have a positive equilibrium. Using a Lyapunov function, we give a global asymptotic stability result for the competition model of several species. The operating diagram describes the asymptotic behavior of this model with respect to control parameters and illustrates the effect of the intra-specific competition on the coexistence region of the species.

Practical coexistence of two species in the chemostat - a slow-fast characterization

We show that the chemostat model with two species having different but close break-even concentrations exhibits a slow-fast dynamics. Considering small perturbations about the dilution rate for which break-even concentrations are identical, we use the Fenichel theory to show the coexistence of species for large times. Then we determine the reduced dynamics, which is non-trivial and characterized by the slopes of the growth functions about their break-even concentrations.