Global dynamics of a chemostat competition model with distributed delay

 We study the global dynamics of n-species competition in a chemostat with distributed delay describing the time-lag involved in the conversion of nutrient to viable biomass. The delay phenomenon is modelled by the gamma distribution. The linear chain trick and a fluctuation lemma are applied to obtain the global limiting behavior of the model. When each population can survive if it is cultured alone, we prove that at most one competitor survives. The winner is the population that has the smallest delayed break-even concentration, provided that the orders of the delay kernels are large and the mean delays modified to include the washout rate (which we call the virtual mean delays) are bounded and close to each other, or the delay kernels modified to include the washout factor (which we call the virtual delay kernels) are close in L ¹-norm. Also, when the virtual mean delays are relatively small, it is shown that the predictions of the distributed delay model are identical with the predictions of the corresponding ODEs model without delay. However, since the delayed break-even concentrations are functions of the parameters appearing in the delay kernels, if the delays are sufficiently large, the prediction of which competitor survives, given by the ODEs model, can differ from that given by the delay model. Received: 9 August 1997 / Revised version: 2 July 1998     

Population dynamics and competition in chemostat models with adaptive nutrient uptake

 The standard Monod model for microbial population dynamics in the chemostat is modified to take into consideration that cells can adapt to the change of nutrient concentration in the chemostat by switching between fast and slow nutrient uptake and growing modes with asymmetric thresholds for transition from one mode to another. This is a generalization of a modified Monod model which considers adaptation by transition between active growing and quiescent cells. Global analysis of the model equations is obtained using the theory of asymptotically autonomous systems. Transient oscillatory population density and hysteresis growth pattern observed experimentally, which do not occur for the standard Monod model, can be explained by such adaptive mechanism of the cells. Competition between two species that can switch between fast and slow nutrient uptake and growing modes is also considered. It is shown that generically there is no coexistence steady state, and only one steady state, corresponding to the survival of at most one species in the chemostat, is a local attractor. Numerical simulations reproduce the qualitative feature of some experimental data which show that the population density of the winning species approaches a positive steady state via transient oscillations while that of the losing species approaches the zero steady state monotonically. Received 4 August 1995; received in revised form 15 December 1995     

The asymptotic behavior of a chemostat model with the Beddington-DeAngelis functional response

In this paper, the global asymptotic behavior of a chemostat model with Beddington-DeAngelis functional response is studied. The conditions for the global asymptotical stability of the model with time delays are obtained via monotone dynamical systems. Our results demonstrate that those time delays affect the competitive outcome of the organisms.     

Periodic solutions in a model of competition between plasmid-bearing and plasmid-free organisms in a chemostat with an inhibitor

We obtain necessary and sufficient conditions on the existence of a unique positive equilibrium point and a set of sufficient conditions on the existence of periodic solutions for a 3-dimensional system which arises from a model of competition between plasmid-bearing and plasmid-free organisms in a chemostat with an inhibitor. Our results improve the corresponding results obtained by Hsu, Luo, and Waltman [1]. Received: 20 November 1997 / Revised version: 12 February 1999 / Published online: 20 December 2000     

Periodic response to periodic forcing of the Droop equations for phytoplankton growth

The dynamics of a phytoplankton population growing in a chemostat under a periodic supply of nutrients is investigated with the model proposed by Droop. This model differs from the well-known Monod equations by incorporating nutrient storage by the cells. In spite of its nonlinearity and the time delays introduced by an internal nutrient pool, the model predicts a simple response to a periodic nutrient supply. The population is shown to oscillate with the same frequency as the forcing. To prove the existence of a periodic solution local and global bifurcation results are used. This work establishes a basis on which to evaluate experimental data against the model as a representation of the nutrient-phytoplankton interaction when nutrients fluctuate.     

Modeling aggregation and growth processes in an algal population model: analysis and computations

 Aggregation, the formation of large particles through multiple collision of smaller ones is a highly visible phenomena in oceanic waters which can control material flux to the deep sea. Oceanic aggregates more than 1 cm in diameter have been observed and are frequently described to consist of phytoplankton cells as well as other organic matter such as fecel pellets and mucus nets from pteropods. Division of live phytoplankton cells within an aggregate can also increase the size of aggregate (assuming some daughter cells stay in the aggregate) and hence could be a significant factor in speeding up the formation process of larger aggregate. Due to the difficulty of modeling cell division within aggregates, few efforts have been made in this direction. In this paper, we propose a size structured approach that includes growth of aggregate size due to both cell division and aggregation. We first examine some basic mathematical issues associated with the development of a numerical simulation of the resulting algal aggregation model. The numerical algorithm is then used to examine the basic model behavior and present a comparison between aggregate distribution with and without division in aggregates. Results indicate that the inclusion of a growth term in aggregates, due to cell division, results in higher densities of larger aggregates; hence it has the impact to speed clearance of organic matter from the surface layer of the ocean. Received 1 July 1994; received in revised form 23 February 1996     

Optimal allocation between nutrient uptake and growth in a microbial trichome

 A microbial trichome grows by assimilating nutrients from its environment, and converting these into catalytic macro-molecular machinery. This machinery may be divided into assimilatory machinery and proliferative machinery. The former type is involved in nutrient uptake, whereas the latter type enables the trichome to grow. The cells in the trichome are faced with an allocation problem: given the availability of nutrients in the environment, how many macro-molecular building blocks should be allocated to the synthesis of assimilatory machinery, and how many to the synthesis of proliferative machinery? We answer this question for a particular model, which is a generalization of the Droop quota model. We formulate a two-dimensional non-linear optimal control problem, corresponding to this model. An optimal allocation regime with a singular segment is derived, based on Pontryagin’s maximum principle. We give a direct proof of optimality. We discuss how actual biological cells might implement this optimal regime. Received: 16 December 1996 / Revised version: 14 September 1997     

Multiparametric bifurcations for a model in epidemiology

 In the present paper we make a bifurcation analysis of an SIRS epidemiological model depending on all parameters. In particular we are interested in codimension-2 bifurcations. Received 8 April 1994; received in revised form 29 June 1995     

A model for dengue disease with variable human population

 A model for the transmission of dengue fever with variable human population size is analyzed. We find three threshold parameters which govern the existence of the endemic proportion equilibrium, the increase of the human population size, and the behaviour of the total number of human infectives. We prove the global asymptotic stability of the equilibrium points using the theory of competitive systems, compound matrices, and the center manifold theorem. Received: 3 November 1997 / Revised version: 3 July 1998     

Limit cycles in a chemostat model for a single species with age structure

We model an age-structured population feeding on an abiotic resource by combining the Gurtin-MacCamy [Math. Biosci. 43 (1979) 199] approach with a standard chemostat model. Limit cycles arise by Hopf bifurcations at low values of the chemostat dilution rate, even for simple maternity functions for which the original Gurtin-MacCamy model has no oscillatory solutions. We find the exact location in parameter space of the Hopf bifurcations for special cases of our model. The onset of cycling is largely independent of both the form of the resource uptake function and the shape of the maternity function.     

Nonlinear control for algae growth models in the chemostat

This paper deals with output feedback control of phytoplanktonic algae growth models in the chemostat. The considered class of model is of variable yield type, meaning that the ratio between the environmental nutrient absorption rate and the cells’ growth rate varies, which is different from classical bioprocesses assumptions. On the basis of weak qualitative hypotheses on the analytical expressions of the involved biological phenomena (which guarantee robustness of the procedure toward modeling uncertainties) we propose a nonlinear controller and prove its ability to globally stabilize such processes. Finally, we illustrate our approach with numerical simulations and show its benefits for biological laboratory experiments, especially for ensuring persistence of the culture facing classical experimental problems.     

Analysis of a mathematical model for the growth of tumors

 In this paper we study a recently proposed model for the growth of a nonnecrotic, vascularized tumor. The model is in the form of a free-boundary problem whereby the tumor grows (or shrinks) due to cell proliferation or death according to the level of a diffusing nutrient concentration. The tumor is assumed to be spherically symmetric, and its boundary is an unknown function r=s(t). We concentrate on the case where at the boundary of the tumor the birth rate of cells exceeds their death rate, a necessary condition for the existence of a unique stationary solution with radius r=R ₀ (which depends on the various parameters of the problem). Denoting by c the quotient of the diffusion time scale to the tumor doubling time scale, so that c is small, we rigorously prove that (i) lim inf _t→∞ s(t)>0, i.e. once engendered, tumors persist in time. Indeed, we further show that (ii) If c is sufficiently small then s(t)→R ₀ exponentially fast as t→∞, i.e. the steady state solution is globally asymptotically stable. Further, (iii) If c is not “sufficiently small” but is smaller than some constant γ determined explicitly by the parameters of the problem, then lim sup _t→∞ s(t)<∞; if however c is “somewhat” larger than γ then generally s(t) does not remain bounded and, in fact, s(t)→∞ exponentially fast as t→∞. Received: 25 February 1998 / Revised version: 30 April 1998     

Identifiability of a stochastic model for cell debris in flow cytometry

 One of the most important problems in recovering DNA distribution from flow cytometric DNA measurements is the presence of background noise. In this paper, we analyse a probabilistic model recently proposed for background debris distribution and based on a specific probabilistic mechanism for the DNA fragmentation process of the cell nucleus. In particular, we carry out some sufficient conditions to uniquely identify the original DNA distribution from the flow cytometric data. Received: 15 June 1997 / Revised version: 18 November 1997     

Subharmonic resonance and chaos in forced excitable systems

 Forced excitable systems arise in a number of biological and physiological applications and have been studied analytically and computationally by numerous authors. Existence and stability of harmonic and subharmonic solutions of a forced piecewise-linear Fitzhugh-Nagumo-like system were studied in Othmer ad Watanabe (1994) and in Xie et al. (1996). The results of those papers were for small and moderate amplitude forcing. In this paper we study the existence of subharmonic solutions of this system under large-amplitude forcing. As in the case of intermediate-amplitude forcing, bistability between 1 : 1 and 2 : 1 solutions is possible for some parameters. In the case of large-amplitude forcing, bistability between 2 : 2 and 2 : 1 solutions, which does not occur in the case of intermediate-amplitude forcing, is also possible for some parameters. We identify several new canonical return maps for a singular system, and we show that chaotic dynamics can occur in some regions of parameter space. We also prove that there is a direct transition from 2 : 2 phase-locking to chaos after the first period-doubling bifurcation, rather than via the infinite sequence of period doublings seen in a smooth quadratic interval map. Coexistence of chaotic dynamics and stable phase-locking can also occur. Received: 6 July 1998 / Revised version: 2 October 1998     

Interaction of maturation delay and nonlinear birth in population and epidemic models

 A population with birth rate function B(N) N and linear death rate for the adult stage is assumed to have a maturation delay T>0. Thus the growth equation N′(t)=B(N(t−T)) N(t−T) e⁻ d ₁ T−dN(t) governs the adult population, with the death rate in previous life stages d ₁≧0. Standard assumptions are made on B(N) so that a unique equilibrium N _e exists. When B(N) N is not monotone, the delay T can qualitatively change the dynamics. For some fixed values of the parameters with d ₁>0, as T increases the equilibrium N _e can switch from being stable to unstable (with numerically observed periodic solutions) and then back to stable. When disease that does not cause death is introduced into the population, a threshold parameter R ₀ is identified. When R ₀<1, the disease dies out; when R ₀>1, the disease remains endemic, either tending to an equilibrium value or oscillating about this value. Numerical simulations indicate that oscillations can also be induced by disease related death in a model with maturation delay. Received: 2 November 1998 / Revised version: 26 February 1999     

Global analysis of a model of plasmid-bearing,plasmid-free competition in a chemostat

A model of competition between plasmid-bearing and plasmid-free organisms in a chemostat was proposed in a paper of Stephanopoulis and Lapidus. The model was in the form of a system of nonlinear ordinary differential equations. Such models are relevant to commercial production by genetically altered organisms in continuous culture. The analysis there was local (using index arguments). This paper provides a mathematically rigorous analysis of the global asymptotic behavior of the governing equations in the case of uninhibited specific growth rate.Research supported by the National Council of Science, Republic of ChinaResearch supported by National Science Foundation Grant, DMS-9204490Research supported by the Natural Science and Engineering Council of Canada. This author's contribution was made while on research leave visiting the Department of Ecology and Evolutionary Biology at Princeton University. She would especially like to thank Simon Levin for his guidance as well as for providing an exceptional working environment     

Global stability in consumer-resource cascades

 Models of population growth in consumer-resource cascades (serially arranged containers with a dynamic consumer population, v, receiving a flow of resource, u, from the previous container) with a functional response of the form h(u/v ^b ) are investigated. For b∈[0, 1], it is shown that these models have a globally stable equilibrium. As a result, two conclusions can be drawn: (1) Consumer density dependence in the functional or in the per-capita numerical response can result in persistence of the consumer population in all containers. (2) In the absence of consumer density dependence, the consumer goes extinct in all containers except possibly the first. Several variations of this model are discussed including replacing discrete containers by a spatial continuum and introducing a dynamic resource. Received 25 February 1995 / received in revised form 27 July 1995     

The diffusion model for migration and selection in a plant population

 The diffusion approximation is derived for migration and selection at a multiallelic locus in a partially selfing plant population subdivided into a lattice of colonies. Generations are discrete and nonoverlapping; both pollen and seeds disperse. In the diffusion limit, the genotypic frequencies at each point are those determined at equilibrium by the local rate of selfing and allelic frequencies. If the drift and diffusion coefficients are taken as the appropriate linear combination of the corresponding coefficients for pollen and seeds, then the migration terms in the partial differential equation for the allelic frequencies have the standard form for a monoecious animal population. The selection term describes selection on the local genotypic frequencies. The boundary conditions and the unidimensional transition conditions for a geographical barrier and for coincident discontinuities in the carrying capacity and migration rate have the standard form. In the diallelic case, reparametrization renders the entire theory of clines and of the wave of advance of favorable alleles directly applicable to plant populations. Received 30 August 1995; received in revised form 23 February 1996     

A weakly nonlinear analysis of a model of avascular solid tumour growth

 In this paper we study a mathematical model that describes the growth of an avascular solid tumour. Our analysis concentrates on the stability of steady, radially-symmetric model solutions with respect to perturbations taken from the class of spherical harmonics. Using weakly nonlinear analysis, previous results are extended to show how the amplitudes of the asymmetric modes interact. Attention focuses on a special case for which the model equations simplify. Analysis of the simplified model equations leads to the identification of a two-parameter family of asymmetric steady solutions, the dimensions of whose stable and unstable manifolds depend on the system parameters. The asymmetric steady solutions limit the basin of attraction of the radially-symmetric steady state when it is linearly stable. On the basis of these numerical and analytical results we postulate the existence of fully nonlinear steady solutions which are stable with respect to time-dependent perturbations. Received: 25 October 1998 / Revised version: 20 June 1998