Effects of dispersal on the stability of a gonorrhea endemic model

Using Liapunov's direct method, effects of dispersal on the linear and nonlinear stability of the endemic equilibrium state of the system governing the spread of gonorrhea are investigated. It is noted that the equilibrium state, which is nonlinearly asymptotically stable in the feasible region of the phase plane in the absence of dispersal, remains so with self-dispersal also (cross-dispersal being absent). However, in the presence of both self- and cross-dispersal, the equilibrium state can still remain nonlinearly asymptotically stable in the entire feasible region provided a certain condition involving self- and cross-dispersal coefficients is satisfied. It is also seen in this case that, for the linearly stable equilibrium state, there exists a subregion of the feasible region where it is nonlinearly asymptotically stable.     

Global Stability of a System of Two Interacting Species with Convective and Dispersive Migration

Using Liapunov's direct method, effects of convective and dispersive migration on the global stability of the equilibrium state of a system of two interacting species are investigated. It is shown that the stable equilibrium state without dispersal remains so with dispersal. Further, it is pointed out that stability or instability of the equilibrium state of the system is not affected by convective migration. These results are justified in cases of a system of mutualistic interactions of species and a prey-predator system with functional response.     

Effects of convective and dispersive migration on the linear stability of a two species system with mutualistic interactions and functional response

In this paper, effects of convective and dispersive migration on the linear stability of the equilibrium state of a two species system with mutualistic interactions and functional response have been investigated. In both finite and semi-infinite habitats, it has been shown that the otherwise stable equilibrium state without dispersal remains so with dispersal also, both under flux and reservior conditions. In the case of finite habitat, the degree of stability increases as dispersal coefficients of the two species increase. The effect of convective migration also is to stabilize the equilibrium state in this case.     

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.     

Dispersal delays, predator-prey stability, and the paradox of enrichment

It takes time for individuals to move from place to place. This travel time can be incorporated into metapopulation models via a delay in the interpatch migration term. Such a term has been shown to stabilize the positive equilibrium of the classical Lotka-Volterra predator-prey system with one species (either the predator or the prey) dispersing. We study a more realistic, Rosenzweig-MacArthur, model that includes a carrying capacity for the prey, and saturating functional response for the predator. We show that dispersal delays can stabilize the predator-prey equilibrium point despite the presence of a Type II functional response that is known to be destabilizing. We also show that dispersal delays reduce the amplitude of oscillations when the equilibrium is unstable, and therefore may help resolve the paradox of enrichment.     

Evolution of dispersal in density regulated populations: A haploid model

The evolution of dispersal is explored in a density-dependent framework. Attention is restricted to haploid populations in which the genotypic fitnesses at a single diallelic locus are decreasing functions of the changing number of individuals in the population. It is shown that migration between two populations in which the genotypic response to density is reversed can maintain both alleles when the intermigration rates are constant or nondecreasing functions of the population densities. There is always a unique symmetric interior equilibrium with equal numbers but opposite gene frequencies in the two populations, provided the system is not degenerate. Numerical examples with exponential and hyperbolic fitnesses suggest that this is the only stable equilibrium state under constant positive migration rates (m) less than . Practically speaking, however, there is only convergence after a reasonable number of generations for relatively small migration rates ( ). A migration-modifying mutant at a second, neutral locus, can successfully enter two populations at a stable migration-selection balance if and only if it reduces the intermigration rates of its carriers at the original equilibrium population size. Moreover, migration modification will always result in a higher equilibrium population size, provided the system approaches another symmetric interior equilibrium. The new equilibrium migration rate will be lower than that at the original equilibrium, even when the modified migration rate is a nondecreasing function of the population sizes. Therefore, as in constant viability models, evolution will lead to reduced dispersal.     

Structural sensitivity and resilience in a predator–prey model with density-dependent mortality

Numerous formulations with the same mathematical properties can be relevant to model a biological process. Different formulations can predict different model dynamics like equilibrium vs. oscillations even if they are quantitatively close (structural sensitivity). The question we address in this paper is: does the choice of a formulation affect predictions on the number of stable states? We focus on a predator–prey model with predator competition that exhibits multiple stable states. A bifurcation analysis is realized with respect to prey carrying capacity and species body mass ratio within range of values found in food web models. Bifurcation diagrams built for two type-II functional responses are different in two ways. First, the kind of stable state (equilibrium vs. oscillations) is different for 26.0–49.4% of the parameter values, depending on the parameter space investigated. Using generalized modelling, we highlight the role of functional response slope in this difference. Secondly, the number of stable states is higher with Ivlev's functional response for 0.1–14.3% of the parameter values. These two changes interact to create different model predictions if a parameter value or a state variable is altered. In these two examples of disturbance, Holling's disc equation predicts a higher system resilience. Indeed, Ivlev's functional response predicts that disturbance may trap the system into an alternative stable state that can be escaped from only by a larger alteration (hysteresis phenomena). Two questions arise from this work: (i) how much complex ecological models can be affected by this sensitivity to model formulation? and (ii) how to deal with these uncertainties in model predictions?     

Prey dispersal and predator persistence

To understand how patchiness influences population dynamics of a tri-trophic interaction, a tractable model is formulated in terms of differential equations. Motivated by the structure of systems such as plants, phytophagous mites and predatory mites, the model takes dispersal into account at the middle trophic level. The effect of dispersal for the middle level in a tri-trophic system could be either stabilising or destabilising since the middle level acts both as prey and as predator. First a simple model with logistic growth for the lowest level is formulated. A model with logistic growth for the lowest level and instantaneous dispersal has a globally stable three-species equilibrium, if this equilibrium exists. Addition of a middle level dispersal phase of non-negligible duration influences equilibrium stability. In the absence of the top trophic level a limit cycle can occur, caused by the delay that exists in the reaction of the middle level to the changes in the lowest level. With low predator efficiency, it is also possible to have an unstable three-species equilibrium. So addition of a middle level dispersal phase of non-negligible duration can work destabilising. Next the persistence of the third trophic level is studied. Even when the three-species equilibrium exists, the third trophic level need not be persistent. A two-species limit cycle can keep its stability when a three-species equilibrium exists; the system is then bistable. It is argued that such a bistability can offer an alternative explanation for pesticide-induced outbreaks of spider mites and failure of predator introduction.     

Effect of predator density dependent dispersal of prey on stability of a predator-prey system

This work presents a predator-prey Lotka-Volterra model in a two patch environment. The model is a set of four ordinary differential equations that govern the prey and predator population densities on each patch. Predators disperse with constant migration rates, while prey dispersal is predator density-dependent. When the predator density is large, the dispersal of prey is more likely to occur. We assume that prey and predator dispersal is faster than the local predator-prey interaction on each patch. Thus, we take advantage of two time scales in order to reduce the complete model to a system of two equations governing the total prey and predator densities. The stability analysis of the aggregated model shows that a unique strictly positive equilibrium exists. This equilibrium may be stable or unstable. A Hopf bifurcation may occur, leading the equilibrium to be a centre. If the two patches are similar, the predator density dependent dispersal of prey has a stabilizing effect on the predator-prey system.     

Two-patch population models with adaptive dispersal: the effects of varying dispersal speeds

The population-dispersal dynamics for predator–prey interactions and two competing species in a two patch environment are studied. It is assumed that both species (i.e., either predators and their prey, or the two competing species) are mobile and their dispersal between patches is directed to the higher fitness patch. It is proved that such dispersal, irrespectively of its speed, cannot destabilize a locally stable predator–prey population equilibrium that corresponds to no movement at all. In the case of two competing species, dispersal can destabilize population equilibrium. Conditions are given when this cannot happen, including the case of identical patches.     

Transient dynamics and pattern formation: reactivity is necessary for Turing instabilities.

The theory of spatial pattern formation via Turing bifurcations - wherein an equilibrium of a nonlinear system is asymptotically stable in the absence of dispersal but unstable in the presence of dispersal - plays an important role in biology, chemistry and physics. It is an asymptotic theory, concerned with the long-term behavior of perturbations. In contrast, the concept of reactivity describes the short-term transient behavior of perturbations to an asymptotically stable equilibrium. In this article we show that there is a connection between these two seemingly disparate concepts. In particular, we show that reactivity is necessary for Turing instability in multispecies systems of reaction-diffusion equations, integrodifference equations, coupled map lattices, and systems of ordinary differential equations.     

Global dispersal reduces local diversity

Metapopulation models and stepping-stone models in genetics are based on very different underlying dispersal structures, yet it can be difficult to distinguish the behaviour of the two kinds of models. We demonstrate a striking qualitative difference in the equilibrium behaviour possible with these two kinds of dispersal. If, in a local patch, there are multiple stable equilibria (and consequently an unstable equilibrium), we demonstrate that, for the spatial system with a metapopulation structure, at equilibrium every patch has to be near one of the stable equilibria. This contrasts with the clinal structure possible with a stepping-stone or continuous space model; thus the result can be used to deduce qualitative information about the form of dispersal from observations of allele frequencies.     

Evolutionary games for spatially dispersed populations

An evolutionary game model is developed that incorporates both spatial dispersion and density effects in the evolutionary dynamic. It is shown that a stable equilibrium (e.g. an evolutionarily stable strategy) of the non-dispersed frequency dynamic becomes a stable equilibrium of the larger system if population density stabilizes at these fixed frequencies. It is also shown, by example, that other equilibria, whose frequencies change from one location to another, may appear when dispersal rates are relatively small.Research supported by Natural Sciences and Engineering Research Council of Canada Operating Grant A6187Research supported by Natural Sciences and Engineering Research Council of Canada Operating Grant A7822     

Effects of convective and dispersive interactions on the stability of two species

The evolution and local stability of a system of two interacting species in a finite two-dimensional habitat is investigated by taking into account the effects of self- and cross-dispersion and convection of the species. In absence of cross-dispersion, an equilibrium state which is stable without dispersion is always stable with dispersion provided that the dispersion coefficients of the two species are equal. However, when the dispersion coefficients of the two species are different, the possibility of self-dispersive instability arises. It is also pointed out that the cross-dispersion of species may lead to stability or instability depending upon the nature and the magnitude of the cross-dispersive interactions in comparison to the self-dispersive interactions. The self-convective movement of species increases the stability of the equilibrium state and can stabilize an otherwise unstable equilibrium state. The effect of cross-convection (in absence of self-dispersion and self-convection) is to stabilize the equilibrium state in a prey-predator model with positive cross-dispersion coefficients for the prey species. Finally, it is shown that if the system is stable under homogeneous boundary conditions it remains so under non-homogeneous boundary conditions.     

ASYMMETRIC PATCH SIZE DISTRIBUTION LEADS TO DISRUPTIVE SELECTION ON DISPERSAL

Numerous models have been designed to understand how dispersal ability evolves when organisms live in a fragmented landscape. Most of them predict a single dispersal rate at evolutionary equilibrium, and when diversification of dispersal rates has been predicted, it occurs as a response to perturbation or environmental fluctuation regimes. Yet abundant variation in dispersal ability is observed in natural populations and communities, even in relatively stable environments. We show that this diversification can operate in a simple island model without temporal variability: disruptive selection on dispersal occurs when the environment consists of many small and few large patches, a common feature in natural spatial systems. This heterogeneity in patch size results in a high variability in the number of related patch mates by individual, which, in turn, triggers disruptive selection through a high per capita variance of inclusive fitness. Our study provides a likely, parsimonious and testable explanation for the diversity of dispersal rates encountered in nature. It also suggests that biological conservation policies aiming at preserving ecological communities should strive to keep the distribution of patch size sufficiently asymmetric and variable.     

The evolution of dispersal distance in spatially-structured populations

Most evolutionary models of dispersal have concentrated on dispersal rate, with emigration being either global or restricted to nearest neighbours. Yet most organisms fall into an intermediate region where most dispersal is local but there is a wide range of dispersal distances. We use an individual-based model with 2500 patches each with identical local dynamics and show that the dispersal distance is under selection pressure. The dispersal distance that evolves is critically dependent on the ecological dynamics. When the cost of dispersal increases linearly with distance, selection is for short-distance dispersal under stable and damped local dynamics but longer distance dispersal is favoured as local dynamics become more complex. For the cases of stable, damped and periodic patch dynamics global patch synchrony occurs even with very short-distance dispersal. Increasing the scale of dispersal for chaotic local dynamics increases the scale of synchrony but global synchrony does not neccesarily occur. We discuss these results in the light of other possible causes of dispersal and argue for the importance of incorporating non-equilibrium population dynamics into evolutionary models of dispersal distance.     

Coexistence in competition models with density-dependent mortality

We consider a two-competitor/one-prey model in which both competitors exhibit a general functional response and one of the competitors exhibits a density-dependent mortality rate. It is shown that the two competitors can coexist upon the single prey. As an example, we consider a two-competitor/one-prey model with a Holling II functional response. Our results demonstrate that density-dependent mortality in one of the competitors can prevent competitive exclusion. Moreover, by constructing a Liapunov function, the system has a globally stable positive equilibrium.     

Some circumstances assuring monomorphism in subdivided populations

It is well known that in a subdivided population subject to soft selection with two alleles at one locus, instability of both fixation states (a “protected polymorphism”) entails at least one stable polymorphic equilibrium. Although stable polymorphic and monomorphic equilibria can coexist in general, a stable fixation state (monomorphic equilibrium) precludes the existence of any polymorphic equilibrium under the circumstances of haploid or submultiplicative diploid viabilities. This provides that a stable monomorphism is robust against random fluctuations in allele frequencies. It also increases the known circumstances where there is a unique globally attracting stable equilibrium, i.e., where allele frequencies are determined by the selection-migration structure independent of the history of the system.     

Stabilization through spatial pattern formation in metapopulations with long-range dispersal

Many studies of metapopulation models assume that spatially extended populations occupy a network of identical habitat patches, each coupled to its nearest neighbouring patches by density-independent dispersal. Much previous work has focused on the temporal stability of spatially homogeneous equilibrium states of the metapopulation, and one of the main predictions of such models is that the stability of equilibrium states in the local patches in the absence of migration determines the stability of spatially homogeneous equilibrium states of the whole metapopulation when migration is added. Here, we present classes of examples in which deviations from the usual assumptions lead to different predictions. In particular, heterogeneity in local habitat quality in combination with long-range dispersal can induce a stable equilibrium for the metapopulation dynamics, even when within-patch processes would produce very complex behaviour in each patch in the absence of migration. Thus, when spatially homogeneous equilibria become unstable, the system can often shift to a different, spatially inhomogeneous steady state. This new global equilibrium is characterized by a standing spatial wave of population abundances. Such standing spatial waves can also be observed in metapopulations consisting of identical habitat patches, i.e. without heterogeneity in patch quality, provided that dispersal is density dependent. Spatial pattern formation after destabilization of spatially homogeneous equilibrium states is well known in reaction–diffusion systems and has been observed in various ecological models. However, these models typically require the presence of at least two species, e.g. a predator and a prey. Our results imply that stabilization through spatial pattern formation can also occur in single-species models. However, the opposite effect of destabilization can also occur: if dispersal is short range, and if there is heterogeneity in patch quality, then the metapopulation dynamics can be chaotic despite the patches having stable equilibrium dynamics when isolated. We conclude that more general metapopulation models than those commonly studied are necessary to fully understand how spatial structure can affect spatial and temporal variation in population abundance.     

具Holling-Ⅳ型功能反应函数的捕食者-食饵系统的定性分析(英文)

本文研究一类具Holling-Ⅳ型功能反应函数的捕食者-食饵模型.对模型进行定性分析得知系统正解都是有界的;因此,当平衡点不稳定,系统至少存在一稳定的极限环.本文还运用Poincare形式级数法,得到了正平衡点至多为二阶稳定细焦点的结论.并基于Hopf分支理论得知系统在一定条件下至少存在两个极限环.