Coalescence of interacting cell populations

We analyse the coalescence of invasive cell populations by studying both the temporal and steady behaviour of a system of coupled reaction-diffusion equations. This problem is relevant to recent experimental observations of the dynamics of opposingly directed invasion waves of cells. Two cell types, u and v, are considered with the cell motility governed by linear or nonlinear diffusion. The cells proliferate logistically so that the long-term total cell density, u+v approaches a carrying capacity. The steady-state solutions for u and v are denoted u(s) and v(s). The steady solutions are spatially invariant and satisfy u(s)+v(s)=1. However, this expression is underdetermined so the relative proportion of each cell type u(s) and v(s) cannot be determined a priori. Various properties of this model are studied, such as how the relative proportion of u(s) and v(s) depends on the relative motility and relative proliferation rates. The model is analysed using a combination of numerical simulations and a comparison principle. This investigation unearths some novel outcomes regarding the role of overcrowding and cell death in this type of cell migration assay. These observations have relevance to experimental design and interpretation regarding the identification and parameterisation of mechanisms involved in cell invasion.     

Models for the effect of toxicant in single-species and predator-prey systems

Models of single-species and predator-prey systems in a polluted closed environment are developed and partially analyzed. Three cases are considered: a single influx of toxicant, a constant influx of toxicant, and a periodic pollution of the environment. In the case of single-species growth we are able to determine some local and global dynamics. In the case of predator-prey systems, we investigate the existence of steady states for a small constant influx of toxicant.On leave from Department of Mathematics, Indian Institute of Technology, Kanpur, IndiaResearch partially supported by the Natural Sciences and Engineering Research Council of Canada, Grant No. NSERC A4823     

Scale-free extinction dynamics in spatially structured host-parasitoid systems

Much of the work on extinction events has focused on external perturbations of ecosystems, such as climatic change, or anthropogenic factors. Extinction, however, can also be driven by endogenous factors, such as the ecological interactions between species in an ecosystem. Here we show that endogenously driven extinction events can have a scale-free distribution in simple spatially structured host-parasitoid systems. Due to the properties of this distribution there may be many such simple ecosystems that, although not strictly permanent, persist for arbitrarily long periods of time. We identify a critical phase transition in the parameter space of the host-parasitoid systems, and explain how this is related to the scale-free nature of the extinction process. Based on these results, we conjecture that scale-free extinction processes and critical phase transitions of the type we have found may be a characteristic feature of many spatially structured, multi-species ecosystems in nature. The necessary ingredient appears to be competition between species where the locally inferior type disperses faster in space. If this condition is satisfied then the eventual outcome depends subtly on the strength of local superiority of one species versus the dispersal rate of the other.     

Ratio-Dependent Predator-Prey Models of Interacting Populations

Ratio-dependent predator-prey models are increasingly favored by both the theoretical and experimental ecologists as a more suitable alternative to describe predator-prey interactions when the predators hunt seriously. In this article, the classical Bazykin’s model is modified with ratio-dependent functional response. Stability and bifurcation situations of the system are observed. Since the ratio-dependent model always has difficult dynamics in the vicinity of the origin, the analytical behavior of the system near origin is studied completely. It is found that paradox of enrichment can happen to this system under certain parameter values, although the functional response is ratio-dependent. The parametric space for Turing spatial structure is determined. We also conclude that competition among the predator population might be beneficial for predator species under certain circumstances. Finally, ecological interpretations of our results are presented in the discussion section.     

Approach to equilibrium in age structured populations with an increasing recruitment process

This paper considers a model for the dynamics of an age structured population subject to a density dependent factor which regulates the recruitment. Certain properties of biological interest are obtained and the stability of the equilibrium age distributions is investigated. Finally some applications to known fishery models are considered.Work done under the contract 80.02333.01 of C.N.R.     

Survival in continuous structured populations models

The extended McKendrick-von Foerster structured population model is employed to derive a nonautonomous ordinary differential equation model of a population. The derivation assumes that the individual life history can be delineated into several physiological stages. We study the persistence of the population when the model is autonomous and base the nonautonomous survival analysis on the autonomous case and a comparison principle. A brief excursion into alternate life history strategies is presented.This work was supported in part by the U.S. Environmental Protection Agency under cooperative agreement CR 813353010     

A detailed study of the Beddington-DeAngelis predator-prey model

The present investigation accounts for the influence of intra-specific competition among predators in the original Beddington-DeAngelis predator-prey model. We offer a detailed mathematical analysis of the model to describe some of the significant results that may be expected to arise from the interplay of deterministic and stochastic biological phenomena and processes. In particular, stability (local and global) and bifurcation (Saddle-node, Transcritical, Hopf-Andronov, Bogdanov-Takens) analysis of this model are conducted. Corresponding results from previous well known predator-prey models are compared with the current findings. Nevertheless, we also allow this model in stochastic environment with the influences of both, uncorrelated “white” noise and correlated “coloured” noise. This showing that competition among the predator population is beneficial for a number of predator-prey models by keeping them stable around its positive interior equilibrium (i.e. when both populations co-exist), under environmental stochasticity. Comparisons of these findings with the results of some earlier related investigations allow the general conclusion that predator intra-species competition benefits the predator-prey system under both deterministic and stochastic environments. Finally, an extended discussion of the ecological implications of the analytical and numerical results concludes the paper.     

Stability regions in predator-prey systems with constant-rate prey harvesting

Summary We analyze the global behavior of a predator-prey system, modelled by a pair of non-linear ordinary differential equations, under constant-rate prey harvesting. By methods analogous to those used to study predator harvesting, we characterize the theoretically possible structures and transitions. With the aid of a computer simulation we construct examples to show which of these transitions can be realized in a biologically plausible model.Sponsored by the United States Army under Contract No. DAAG29-75-C-0024 and the National Research Council of Canada, Grant No. 67-3138.     

Persistence of structured populations in random environments

Environmental fluctuations often have different impacts on individuals that differ in size, age, or spatial location. To understand how population structure, environmental fluctuations, and density-dependent interactions influence population dynamics, we provide a general theory for persistence for density-dependent matrix models in random environments. For populations with compensating density dependence, exhibiting “bounded” dynamics, and living in a stationary environment, we show that persistence is determined by the stochastic growth rate (alternatively, dominant Lyapunov exponent) when the population is rare. If this stochastic growth rate is negative, then the total population abundance goes to zero with probability one. If this stochastic growth rate is positive, there is a unique positive stationary distribution. Provided there are initially some individuals in the population, the population converges in distribution to this stationary distribution and the empirical measures almost surely converge to the distribution of the stationary distribution. For models with overcompensating density-dependence, weaker results are proven. Methods to estimate stochastic growth rates are presented. To illustrate the utility of these results, applications to unstructured, spatially structured, and stage-structured population models are given. For instance, we show that diffusively coupled sink populations can persist provided that within patch fitness is sufficiently variable in time but not strongly correlated across space.     

Persistence of predator-prey systems in an uncertain environment

Summary The time derivatives of prey and predator populations are assumed to satisfy a set of inequalities, instead of a precise differential equation, reflecting an uncertain environmental and/or lack of knowledge by the modeler. A system of differential equations is found whose solution gives the boundary of a persistent set, which is positive flow invariant for any system satisfying the inequalities. Conditions are given for the persistent set to be bounded away from both axes, which show that resonance effects cannot drive either predator or prey to extinction if that does not happen for an autonomous system satisfying the inequalities. In general predator-prey systems are more persistent when there is strong asymptotic stability, when there is correlation between prey and predator dynamics, when the effect of perturbations is density dependent, and are more persistent under perturbations of the prey than of the predator.     

Coexistence properties of some predator-prey systems under constant rate harvesting and stocking

The global behaviour of a class of predator-prey systems, modelled by a pair of non-linear ordinary differential equations, under constant rate harvesting and/or stocking of both species, is presented. Theoretically possible structures and transitions are developed and validated by computer simulations. The results are presented as transition loci in the F-G (prey harvest rate-predator harvest rate) plane.Sponsored by the United States Army under Contract No. DAAG29-80-C-0041 and in part by NSERC of Canada, Grant No. 67-3138The authors wish to thank Mr. Al MacKenzie of the Department of Electrical Engineering, University of British Columbia, for preparing the figures in this paper.     

Self-organization of mobile populations in cyclic competition

The formation of out-of-equilibrium patterns is a characteristic feature of spatially extended, biodiverse, ecological systems. Intriguing examples are provided by cyclic competition of species, as metaphorically described by the ‘rock-paper-scissors’ game. Both experimentally and theoretically, such non-transitive interactions have been found to induce self-organization of static individuals into noisy, irregular clusters. However, a profound understanding and characterization of such patterns is still lacking. Here, we theoretically investigate the influence of individuals’ mobility on the spatial structures emerging in rock-paper-scissors games. We devise a quantitative approach to analyze the spatial patterns self-forming in the course of the stochastic time evolution. For a paradigmatic model originally introduced by May and Leonard, within an interacting particle approach, we demonstrate that the system's behavior—in the proper continuum limit—is aptly captured by a set of stochastic partial differential equations. The system's stochastic dynamics is shown to lead to the emergence of entangled rotating spiral waves. While the spirals’ wavelength and spreading velocity is demonstrated to be accurately predicted by a (deterministic) complex Ginzburg-Landau equation, their entanglement results from the inherent stochastic nature of the system. These findings and our methods have important applications for understanding the formation of noisy patterns, e.g. in ecological and evolutionary contexts, and are also of relevance for the kinetics of (bio)-chemical reactions.     

How predator functional responses and Allee effects in prey affect the paradox of enrichment and population collapses

In Rosenzweig-MacArthur models of predator-prey dynamics, Allee effects in prey usually destabilize interior equilibria and can suppress or enhance limit cycles typical of the paradox of enrichment. We re-evaluate these conclusions through a complete classification of a wide range of Allee effects in prey and predator's functional response shapes. We show that abrupt and deterministic system collapses not preceded by fluctuating predator-prey dynamics occur for sufficiently steep type III functional responses and strong Allee effects (with unstable lower equilibrium in prey dynamics). This phenomenon arises as type III functional responses greatly reduce cyclic dynamics and strong Allee effects promote deterministic collapses. These collapses occur with decreasing predator mortality and/or increasing susceptibility of the prey to fall below the threshold Allee density (e.g. due to increased carrying capacity or the Allee threshold itself). On the other hand, weak Allee effects (without unstable equilibrium in prey dynamics) enlarge the range of carrying capacities for which the cycles occur if predators exhibit decelerating functional responses. We discuss the results in the light of conservation strategies, eradication of alien species, and successful introduction of biocontrol agents.     

Ecological principles of species distribution models: the habitat matching rule

Most species distribution models (SDMs) assume that habitats are closed, stable and without competition. In that environmental context, it is ecologically correct to assume that members of a species will be distributed in direct relation to the suitability of the habitat, that is, according to the so‐called habitat matching rule. This paper examines whether it is possible to maintain the assumption of the habitat matching rule in the following circumstances: (1) when habitats are connected and organisms can move between them, (2) when there are disturbances and seasonal cycles that generate instability, and (3) when there is inter‐specific and intra‐specific competition. Here I argue that it is possible as long as the following aspects are taken into account. In open habitats at equilibrium, in which habitat selection and competition operate, the habitat matching rule can be applied in some conditions, while competition tends to homogenize the species distribution in other environmental contexts. In the latter case, two methods can be used to incorporate these effects into SDMs: new parameters can be incorporated into the response functions, or the occurrence of proportions of categories of individuals (adult/young, male/female, or dominant/subordinate species in guilds) can be used instead of the occurrence of organisms. The habitat matching rule is not fulfilled in non‐equilibrium environments. The solution to this problem lies in the design of SDMs with two strategies that depend on scale. Locally, the disequilibrium can be encapsulated using average environmental conditions, with sufficiently large cells (in the case of metapopulations) and/or long enough sampling periods (in the case of seasonal cycles). At coarse scales, the use of presence‐only models can in some cases avoid the destabilizing effect of catastrophic historical processes. The matching law is a strong assumption of SDMs because it is based on population ecology theory and the principle of evolution by natural selection.     

Equilibria in an epistatic viability model under arbitrary strength of selection

A class of viability models that generalize the standard additive model for the case of pairwise additive by additive epistatic interactions is considered. Conditions for existence and stability of steady states in the corresponding two-locus model are analyzed. Using regular perturbation techniques, the case when selection is weaker than recombination and the case when selection is stronger than recombination are investigated. The results derived are used to make conclusions on the dependence of population characteristics on the relation between the strength of selection and the recombination rate.     

On the concept of attractor for community-dynamical processes II: the case of structured populations

In Part I of this paper Jacobs and Metz (2003) extended the concept of the Conley-Ruelle, or chain, attractor in a way relevant to unstructured community ecological models. Their modified theory incorporated the facts that certain parts of the boundary of the state space correspond to the situation of at least one species being extinct and that an extinct species can not be rescued by noise. In this part we extend the theory to communities of physiologically structured populations. One difference between the structured and unstructured cases is that a structured population may be doomed to extinction and not rescuable by any biologically relevant noise before actual extinction has taken place. Another difference is that in the structured case we have to use different topologies to define continuity of orbits and to measure noise. Biologically meaningful noise is furthermore related to the linear structure of the community state space. The construction of extinction preserving chain attractors developed in this paper takes all these points into account.     

Coexistence of two types on a single resource in discrete time

Armstrong and McGehee (1980) have shown that two species modeled in continuous time can coexist on a single resource provided that one species oscillates autonomously. This paper demonstrates the parallel result in discrete time. I consider a deterministic model of two asexual types in a single patch competing for a single resource, and show that such systems generically produce oscillatory coexistence or bistability if one of the types displays periodic or chaotic behavior in isolation. The conditions for coexistence or bistability are derived in terms of the convexity of the functions describing fitness as a function of resource availability. I also analyze whether or not a stable type, a type with a stable equilibrium population size when considered in isolation, can invade a periodic orbit of an unstable type, and show that the same convexity condition distinguishes these two cases. The widely considered exponential or Ricker model for population dynamics lies on the boundary between the two cases and is highly degenerate in this context.     

Competing Effects of Toxin-Producing Phytoplankton on Overall Plankton Populations in the Bay of Bengal

The coexistence of a large number of phytoplankton species on a seemingly limited variety of resources is a classical problem in ecology, known as ‘the paradox of the plankton’. Strong fluctuations in species abundance due to the external factors or competitive interactions leading to oscillations, chaos and short-term equilibria have been cited so far to explain multi-species coexistence and biodiversity of phytoplankton. However, none of the explanations has been universally accepted. The qualitative view and statistical analysis of our field data establish two distinct roles of toxin-producing phytoplankton (TPP): toxin allelopathy weakens the interspecific competition among phytoplankton groups and the inhibition due to ingestion of toxic substances reduces the abundance of the grazer zooplankton. Structuring the overall plankton population as a combination of nontoxic phytoplankton (NTP), toxic phytoplankton, and zooplankton, here we offer a novel solution to the plankton paradox governed by the activity of TPP. We demonstrate our findings through qualitative analysis of our sample data followed by analysis of a mathematical model.     

Weakly coupled hyperbolic systems modeling the circulation of FeLV in structured feline populations

Global existence and regularity results are provided for weakly coupled first order hyperbolic systems modeling the propagation of the Feline Leukemia Virus (FeLV), a retrovirus of domestic cats (Felis catus). In a simple example we find a threshold parameter yielding endemic stationary states.