Lagrange stability and ecological systems

This paper considers the ecological and mathematical relationships among various stability concepts. It is shown that diagonally dominant nonlinear systems have positive solutions for positive initial conditions. The boundedness of such positive solutions leads to a total systems concept of stability: Lagrange stability. This concept is discussed as a compatible method for determining the stability of ecosystems and their mathematical models. Stability in the sense of Lagrange is a total systems concept which does not directly address the problem of the stability of equilibria. A topological presentation is used to describe the Lagrange stability of nonlinear ecosystem compartment models. Finally, the concept of practical stability is developed and discussed with respect to maintaining the desired state of ecosystems in the presence of perturbations.     

Global stability in Lotka-Volterra systems with diffusion

Summary Sufficient conditions for global stability inn-species Lotka-Volterra systems with diffusion are derived. Both continuous environments with 0-flux boundary conditions and environments consisting of discrete patches are considered.     

Chemotaxis-induced spatio-temporal heterogeneity in multi-species host-parasitoid systems

When searching for hosts, parasitoids are observed to aggregate in response to chemical signalling cues emitted by plants during host feeding. In this paper we model aggregative parasitoid behaviour in a multi-species host-parasitoid community using a system of reaction-diffusion-chemotaxis equations. The stability properties of the steady-states of the model system are studied using linear stability analysis which highlights the possibility of interesting dynamical behaviour when the chemotactic response is above a certain threshold. We observe quasi-chaotic dynamic heterogeneous spatio-temporal patterns, quasi-stationary heterogeneous patterns and a destabilisation of the steady-states of the system. The generation of heterogeneous spatio-temporal patterns and destabilisation of the steady state are due to parasitoid chemotactic response to hosts. The dynamical behaviour of our system has both mathematical and ecological implications and the concepts of chemotaxis-driven instability and coexistence and ecological change are discussed. I. G. Pearce gratefully acknowledges the financial support of the NERC.     

Preferential mating in symmetric multilocus systems: stability conditions of the central equilibrium

Several multilocus models that incorporate both preferential mating and viability selection are studied. Specifically, a class of symmetric heterozygosity models are considered that assign individuals to phenotypic classes according to which loci are in heterozygous state regardless of the actual allelic content. Otherwise, an arbitrary number of loci, number of alleles per locus, and arbitrary recombination scheme, viability parameters and preferential mating pattern based on phenotypes are allowed. The conditions for the stability of a central polymorphism are indicated and interpreted. The effects of viability and preference selection may be summarized in a single quantity for each phenotypic class, a generalized fitness. Preferential assortative mating alone can produce stability for a central polymorphism as in the case of viability selection when sexual attractiveness or general fitness increases with higher levels of heterozygosity. The situation is more complex with sexual selection.     

Effects of dispersion on stability of multispecies prey-predator systems

The general multispecies prey-predator system with Gompertz's antisymmetric interactions is nonlinearly stable in the absence of dispersion and continues to remain stable with dispersion under both homogeneous reservoir and zero flux boundary conditions in a region containing the equilibrium state. It is proved that a general multispecies food-web model without antisymmetric interactions is stable in the absence of dispersion and remains stable with dispersion in the above-mentioned region.     

A predator–prey refuge system: Evolutionary stability in ecological systems

A refuge model is developed for a single predator species and either one or two prey species where no predators are present in the prey refuge. An individual’s fitness depends on its strategy choice or ecotype (predators decide which prey species to pursue and prey decide what proportion of their time to spend in the refuge) as well as on the population sizes of all three species. It is shown that, when there is a single prey species with a refuge or two prey species with no refuge compete only indirectly (i.e. there is only apparent competition between prey species), that stable resident systems where all individuals in each species have the same ecotype cannot be destabilized by the introduction of mutant ecotypes that are initially selectively neutral. In game-theoretic terms, this means that stable monomorphic resident systems, with ecotypes given by a Nash equilibrium, are both ecologically and evolutionarily stable. However, we show that this is no longer the case when the two indirectly-competing prey species have a refuge. This illustrates theoretically that two ecological factors, that are separately stabilizing (apparent competition and refuge use), may have a combined destabilizing effect from the evolutionary perspective. These results generalize the concept of an evolutionarily stable strategy (ESS) to models in evolutionary ecology. Several biological examples of predator–prey systems are discussed from this perspective.     

Dynamics and stability in coevolutionary ecological systems. I. Community stability and coevolutionarily stable states

An extension of J. Roughgarden's [1979, Theor. Pop. Biol. 9, 388; 1979, "An Introduction to Evolutionary Ecology and Population Genetic Theory," Macmillan, New York] formalism for investigating the effects of coevolution on community structure is presented. The extension assumes the result that a coevolved community is asymptotically stable when coevolution takes place at a genetically noninvasible boundary. This is proved for the general case of n interacting species. From this a community persistence function, phi (P), is defined that allows measuring the domain of attraction for the community as well as the resilience time, that is, the time taken for a perturbation to decay to 1-1/e (63%) of its initial value.     

Effects of diffusion on total biomass in heterogeneous continuous and discrete-patch systems

Theoretical models of populations on a system of two connected patches previously have shown that when the two patches differ in maximum growth rate and carrying capacity, and in the limit of high diffusion, conditions exist for which the total population size at equilibrium exceeds that of the ideal free distribution, which predicts that the total population would equal the total carrying capacity of the two patches. However, this result has only been shown for the Pearl-Verhulst growth function on two patches and for a single-parameter growth function in continuous space. Here, we provide a general criterion for total population size to exceed total carrying capacity for three commonly used population growth rates for both heterogeneous continuous and multi-patch heterogeneous landscapes with high population diffusion. We show that a sufficient condition for this situation is that there is a convex positive relationship between the maximum growth rate and the parameter that, by itself or together with the maximum growth rate, determines the carrying capacity, as both vary across a spatial region. This relationship occurs in some biological populations, though not in others, so the result has ecological implications.     

Modelling the spatio-temporal dynamics of multi-species host-parasitoid interactions: heterogeneous patterns and ecological implications

A mathematical model of the spatio-temporal dynamics of a two host, two parasitoid system is presented. There is a coupling of the four species through parasitism of both hosts by one of the parasitoids. The model comprises a system of four reaction-diffusion equations. The underlying system of ordinary differential equations, modelling the host-parasitoid population dynamics, has a unique positive steady state and is shown to be capable of undergoing Hopf bifurcations, leading to limit cycle kinetics which give rise to oscillatory temporal dynamics. The stability of the positive steady state has a fundamental impact on the spatio-temporal dynamics: stable travelling waves of parasitoid invasion exhibit increasingly irregular periodic travelling wave behaviour when key parameter values are increased beyond their Hopf bifurcation point. These irregular periodic travelling waves give rise to heterogeneous spatio-temporal patterns of host and parasitoid abundance. The generation of heterogeneous patterns has ecological implications and the concepts of temporary host refuge and niche formation are considered.     

Effects of diffusion on the electrophoretic behavior of associating systems: the Gilbert-Jenkins theory revisited

The Gilbert-Jenkins theory predicts the asymptotic shape of moving-boundary sedimentation and electrophoretic patterns and broad zone molecular sieve chromatographic elution profiles for the class of interacting systems, A + B in equilibrium C, in which two dissimilar macromolecules react reversibly to form a complex. A particularly provocative case is the one in which the complex has a greater migration velocity than that of either reactant, each of which has a different velocity. Depending upon conditions, this case predicts, for example, that in the asymptotic limit an ascending electrophoretic pattern or a frontal gel chromatographic elution profile can show two hypersharp reaction boundaries separated by a plateau. This prediction is now confirmed by numerical solution of transport equations which retain the second-order diffusional term and extrapolation of the computed patterns to zero diffusion coefficient. For finite diffusion coefficient, however, the two hypersharp reaction boundaries are separated by a weak negative gradient. These calculations are extended to an examination of the transitions between the three types of patterns admitted by the case under consideration in order to gain physical understanding and to define criteria for recognizing the transitions. Studies of this kind not only establish confidence in the Gilbert-Jenkins theory, but, in addition, they provide new insights which make for more effective application of the theory to real systems.     

Convergence to the equilibrium state in the Volterra-Lotka diffusion equations

Summary We introduce a diffusion term in the Volterra-Lotka model for two interacting species. For certain simple boundary conditions there exists a Lyapunov functional which allows to investigate the asymptotic behavior. Either the solution converges to equilibrium or sustained oscillations occur.     

Effects of diffusion on the electrophoretic behavior of associating systems: The Gilbert-Jenkins theory revisited

The Gilbert-Jenkins theory predicts the asymptotic shape of moving-boundary sedimentation and electrophoretic patterns and broad zone molecular sieve chromatographic elution profiles for the class of interacting systems, A + B ⇄ C, in which two dissimilar macromolecules react reversibly to form a complex. A particularly provocative case is the one in which the complex has a greater migration velocity than that of either reactant, each of which has a different velocity. Depending upon conditions, this case predicts, for example, that in the asymptotic limit an ascending electrophoretic pattern or a frontal gel chromatographic elution profile can show two hypersharp reaction boundaries separated by a plateau. This prediction is now confirmed by numerical solution of transport equations which retain the second-order diffusional term and extrapolation of the computed patterns to zero diffusion coefficient. For finite diffusion coefficient, however, the two hypersharp reaction boundaries are separated by a weak negative gradient. These calculations are extended to an examination of the transitions between the three types of patterns admitted by the case under consideration in order to gain physical understanding and to define criteria for recognizing the transitions. Studies of this kind not only establish confidence in the Gilbert-Jenkins theory, but, in addition, they provide new insights which make for more effective application of the theory to real systems.     

The effects of interference competition on stability,structure and invasion of a multi-species system

Summary An interference competition model for a many species system is presented, based on Lotka-Volterra equations in which some restrictions are imposed on the parameters. The competition coefficients of the Lotka-Volterra equations are assumed to be expressed as products of two factors: the intrinsic interference to other individuals and the defensive ability against such interference. All the equilibrium points of the model are obtained explicitly in terms of its parameters, and these equilibria are classified according to the concept of sector stability. Thus survival or extinction of species at a stable equilibrium point can be determined analytically.The result of the analysis is extended to the successional processes of a community. A criterion for invasion of a new species is obtained and it is also shown that there are some characteristic quantities which show directional changes as succession proceeds.     

Effects of asymmetric dispersal and environmental gradients on the stability of host–parasitoid systems

We consider host–parasitoid systems spatially distributed on a row of patches connected by dispersal. We analyze the effects of dispersal frequency, dispersal asymmetry, number of patches and environmental gradients on the stability of the host–parasitoid interactions. To take into account dispersal frequency, the hosts and parasitoids are allowed to move from one patch to a neighboring patch a certain number of times within a generation. When this number is high, aggregation methods can be used to simplify the proposed initial model into an aggregated model describing the dynamics of both the total host and parasitoid populations. We show that as the number of patches increases less asymmetric parasitoid dispersal rates are required for stability. We found that the 'CV²>1 rule' is a valid approximation for stability if host growth rate is low, otherwise the general condition of stability we establish should be preferred. Environmental variability along the row of patches is introduced as gradients on host growth rate and parasitoid searching efficiency. We show that stability is more likely when parasitoids move preferentially towards patches where they have high searching efficiency or when hosts go mainly to patches where they have a low growth rate.