Global stability and persistence in LG–Holling type II diseased predator ecosystems

A Leslie–Gower–Holling type II model is modified to introduce a contagious disease in the predator population, assuming that disease cannot propagate to the prey. All the system’s equilibria are determined and the behaviour of the system near them is investigated. The main mathematical issues are global stability and bifurcations for some of the equilibria, together with sufficient conditions for persistence of the ecosystem. Counterintuitive results on the role played by intraspecific competition are highlighted.     

Effect of delay in a Lotka-Volterra type predator-prey model with a transmissible disease in the predator species

We consider a system of delay differential equations modeling the predator-prey ecoepidemic dynamics with a transmissible disease in the predator population. The time lag in the delay terms represents the predator gestation period. We analyze essential mathematical features of the proposed model such as local and global stability and in addition study the bifurcations arising in some selected situations. Threshold values for a few parameters determining the feasibility and stability conditions of some equilibria are discovered and similarly a threshold is identified for the disease to die out. The parameter thresholds under which the system admits a Hopf bifurcation are investigated both in the presence of zero and non-zero time lag. Numerical simulations support our theoretical analysis.     

The role of transmissible diseases in the Holling-Tanner predator-prey model

Transmissible diseases are known to induce remarkable major behavioral changes in predator-prey systems. However, little attention has been paid to model such situations. The latter would allow to predict useful applications in both dynamics and control. Here the Holling-Tanner model is revisited to account for the influence of a transmissible disease, under the assumption that it spreads among the prey species only. We have found the equilibria and analyzed the behavior of the system around each one of them. A threshold result determining when the disease dies out has been identified. We also investigated the parametric space under which the system enters into Hopf and transcritical bifurcations, around the disease free equilibrium. The system is shown to experience neither saddle-node nor pitch-fork bifurcation. Global stability results are obtained by constructing suitable Lyapunov functions.     

一类具分段常数变量的捕食-食饵系统的Neimark-Sacker分支

讨论了具时滞与分段常数变量的捕食-食饵生态模型的稳定性及Neimark-Sacker分支;通过计算得到连续模型对应的差分模型,基于特征值理论和Schur-Cohn判据得到正平衡态局部渐进稳定的充分条件;以食饵的内禀增长率为分支参数,运用分支理论和中心流形定理分析了Neimark-Sacker分支的存在性与稳定性条件;通过举例和数值模拟验证了理论的正确性。     

Incorporating prey refuge in a prey–predator model with a Holling type I functional response: random dynamics and population outbreaks

A prey–predator discrete-time model with a Holling type I functional response is investigated by incorporating a prey refuge. It is shown that a refuge does not always stabilize prey–predator interactions. A prey refuge in some cases produces even more chaotic, random-like dynamics than without a refuge and prey population outbreaks appear. Stability analysis was performed in order to investigate the local stability of fixed points as well as the several local bifurcations they undergo. Numerical simulations such as parametric basins of attraction, bifurcation diagrams, phase plots and largest Lyapunov exponent diagrams are executed in order to illustrate the complex dynamical behavior of the system.     

A predator-prey reaction-diffusion system with nonlocal effects

We consider a predator-prey system in the form of a coupled system of reaction-diffusion equations with an integral term representing a weighted average of the values of the prey density function, both in past time and space. In a limiting case the system reduces to the Lotka Volterra diffusion system with logistic growth of the prey. We investigate the linear stability of the coexistence steady state and bifurcations occurring from it, and expressions for some of the bifurcating solutions are constructed. None of these bifurcations can occur in the degenerate case when the nonlocal term is in fact local.     

The dynamics of a Lotka–Volterra predator–prey model with state dependent impulsive harvest for predator

According to the economic and biological aspects of renewable resources management, we propose a Lotka–Volterra predator–prey model with state dependent impulsive harvest. By using the Poincaré map, some conditions for the existence and stability of positive periodic solution are obtained. Moreover, we show that there is no periodic solution with order larger than or equal to three under some conditions. Numerical results are carried out to illustrate the feasibility of our main results. The bifurcation diagrams of periodic solutions are obtained by using the numerical simulations, and it is shown that a chaotic solution is generated via a cascade of period-doubling bifurcations, which implies that the presence of pulses makes the dynamic behavior more complex.     

The dynamics of two-species allelopathic competition with optimal harvesting

This paper analyses a bionomic model of two competitive species in the presence of toxicity with different harvesting efforts. An interesting dynamics in the first quadrant is analysed and two saddle-node bifurcations are detected for different bifurcation parameters. It is noted that under certain parametric restrictions, the model has a unique positive equilibrium point that is globally asymptotically stable whenever it is locally stable. It is also noted that the model can have zero, one or two feasible equilibria appearing through saddle-node bifurcations. The non-existence of a limit cycle in the interior of the first quadrant is also discussed using the Poincare–Dulac criteria. The saddle-node bifurcations are studied using Sotomayor's theorem. Numerical simulations are carried out to validate the analytical findings. The conditions for the existence of bionomic equilibria are discussed and an optimal harvesting policy is derived using Pontryagin's maximum principle.     

Global stability and optimisation of a general impulsive biological control model

An impulsive model of augmentative biological control consisting of a general continuous predator-prey model in ordinary differential equations, i.e. a meta-model, augmented by a discrete part describing periodic introductions of predators is considered. The existence of an invariant periodic solution that corresponds to prey eradication is shown and a condition ensuring its global asymptotic stability is given. An optimisation problem related to the preemptive use of augmentative biological control is then considered. It is assumed that the per time unit budget of biological control (i.e. the number of predators to be released) is fixed and the best deployment of this budget is sought in terms of release frequency. The cost function to be minimised is the time needed to reduce an unforeseen prey (pest) invasion occurring at a worst time instant under some harmless level. The analysis shows that the optimisation problem admits a countable infinite number of solutions. An argumentation considering the required robustness of the optimisation result with respect to the invasive prey population level and to the model parameters is then conducted. It is shown that the cost function is decreasing in the predator release frequency so that the best deployment of the biocontrol agents is to carry out as frequent introductions as possible.     

Predator-prey relationships: Indiscriminate predation

The Gurtin and Levine model⁵ is studied in this paper under the assumption that the fecundity of prey depends on age as well as on the total population sizes of prey and predators. The purpose of this study is to see the effect of this density dependence on the stability criteria for the equilibria of the model equations. It is shown that there are cases when, due to density dependence, the model which is originally neutrally stable becomes stable.     

Long-term variation of link strength in a simple benthic food web

1. The predatory isopod Saduria entomon (L.) and its amphipod prey Monoporeia affinis (Lindstr?m) are key components of the food web in the northern Baltic Sea, together representing 80-90% of the macrobenthic biomass. We use 20 years of stomach content data for Saduria to investigate how diet dynamics affect the stability of the interaction between Saduria and Monoporeia. 2. Consumption of the main prey, Monoporeia, fitted a type III functional response. Consumption rates of the most important alternative prey, mysids, were found to be unrelated to mysid densities but negatively related to the density of Monoporeia. The fit of consumption data to a model that assumes passive prey selection was poor. Thus we conclude that some form of active choice is involved. 3. The effect of consumption of mysids, the alternative prey, on the stability of this system was investigated using a 'one predator-two prey' model with stochastic environmental variation. Analysis of the model suggests that feeding on mysids leads to a decreased extinction risk for the predator, Saduria, and reduced density oscillations for both Saduria and its main prey, Monoporeia.     

The complex dynamics of a diffusive prey–predator model with an Allee effect in prey

This paper investigates complex dynamics of a predator–prey interaction model that incorporates: (a) an Allee effect in prey; (b) the Michaelis–Menten type functional response between prey and predator; and (c) diffusion in both prey and predator. We provide rigorous mathematical results of the proposed model including: (1) the stability of non-negative constant steady states; (2) sufficient conditions that lead to Hopf/Turing bifurcations; (3) a prior estimates of positive steady states; (4) the non-existence and existence of non-constant positive steady states when the model is under zero-flux boundary condition. We also perform completed analysis of the corresponding ODE model to obtain a better understanding on effects of diffusion on the stability. Our analytical results show that the small values of the ratio of the prey's diffusion rate to the predator's diffusion rate are more likely to destabilize the system, thus generate Hopf-bifurcation and Turing instability that can lead to different spatial patterns. Through numerical simulations, we observe that our model, with or without Allee effect, can exhibit extremely rich pattern formations that include but not limit to strips, spotted patterns, symmetric patterns. In addition, the strength of Allee effects also plays an important role in generating distinct spatial patterns.     

具有年龄结构的食蚜蝇-蚜虫捕食模型

现有的具有年龄结构的捕食-食饵模型总是假设只有成年捕食者捕食猎物,这与实际情况不符。建立了一个幼年捕食者捕食食饵的具有年龄结构的食蚜蝇-蚜虫模型,应用微分方程定性理论,讨论了系统平衡点及其稳定性:其中平衡点E1(0,0,0)为不稳定的;满足一定条件时,边界平衡点E2(K,0,0)及正平衡点E3(x*,y1*,y2*)为局部渐近稳定的;且应用一致持续生存理论得到了系统永久持续生存的条件,为有害生物综合治理提供了理论依据。     

Effects of density-dependent migrations on stability of a two-patch predator-prey model

We consider a predator-prey model in a two-patch environment and assume that migration between patches is faster than prey growth, predator mortality and predator-prey interactions. Prey (resp. predator) migration rates are considered to be predator (resp. prey) density-dependent. Prey leave a patch at a migration rate proportional to the local predator density. Predators leave a patch at a migration rate inversely proportional to local prey population density. Taking advantage of the two different time scales, we use aggregation methods to obtain a reduced (aggregated) model governing the total prey and predator densities. First, we show that for a large class of density-dependent migration rules for predators and prey there exists a unique and stable equilibrium for migration. Second, a numerical bifurcation analysis is presented. We show that bifurcation diagrams obtained from the complete and aggregated models are consistent with each other for reasonable values of the ratio between the two time scales, fast for migration and slow for local demography. Our results show that, under some particular conditions, the density dependence of migrations can generate a limit cycle. Also a co-dim two Bautin bifurcation point is observed in some range of migration parameters and this implies that bistability of an equilibrium and limit cycle is possible.     

A comparison of two predator-prey models with Holling's type I functional response

In this paper, we analyze a laissez-faire predator-prey model and a Leslie-type predator-prey model with type I functional responses. We study the stability of the equilibrium where the predator and prey coexist by both performing a linearized stability analysis and by constructing a Lyapunov function. For the Leslie-type model, we use a generalized Jacobian to determine how eigenvalues jump at the corner of the functional response. We show, numerically, that our two models can both possess two limit cycles that surround a stable equilibrium and that these cycles arise through global cyclic-fold bifurcations. The Leslie-type model may also exhibit super-critical and discontinuous Hopf bifurcations. We then present and analyze a new functional response, built around the arctangent, that smoothes the sharp corner in a type I functional response. For this new functional response, both models undergo Hopf, cyclic-fold, and Bautin bifurcations. We use our analyses to characterize predator-prey systems that may exhibit bistability.     

Hopf bifurcation in three-species food chain models with group defense.

Three-species food-chain models, in which the prey population exhibits group defense, are considered. Using the carrying capacity of the environment as the bifurcation parameter, it is shown that the model without delay undergoes a sequence of Hopf bifurcations. In the model with delay it is shown that using a delay as a bifurcation parameter, a Hopf bifurcation can also occur in this case. These occurrences may be interpreted as showing that a region of local stability (survival) may exist even though the positive steady states are unstable. A computer code BIFDD is used to determine the stability of the bifurcation solutions of a delay model.     

Biomass and biodiversity in species-rich tritrophic communities

We consider a tritrophic system with one basal and one top species and a large number of primary consumers, and derive upper and lower bounds for the total biomass of the middle trophic level. These estimates do not depend on dynamical regime, holding for fixed point, periodic, or chaotic dynamics. We have two kinds of estimates, depending on whether the predator abundance is zero. All these results are uniform in a self-limitation parameter, which regulates prey diversity in the system. For strong self-limitation, diversity is large; for weak self-limitation, it is small. Diversity depends on the variance of species’ parameter values. The larger this variance, the lower the diversity, and vice versa. Moreover, variation in the parameters of the Holling type II functional response changes the bifurcation character, with the equilibrium state with nonzero predator abundance losing stability. If that variation is small then the bifurcation can lead to oscillations (the Hopf bifurcation). Under certain conditions, there exists a supercritical Hopf bifurcation. We then find a connection between diversity and Hopf bifurcations. We also show that the system exhibits top-down regulation and a hump-shaped diversity-productivity curve.We then extend the model by allowing species to experience self-regulation. For this extended model, explicit estimates of prey diversity are obtained. We study the dynamics of this system and find the following. First, diversity and system dynamics crucially depend on variation in species parameters. We show that under certain conditions, the system undergoes a supercritical Hopf bifurcation. We also establish a connection between diversity and Hopf bifurcations. For strong self-limitation, diversity is large and complex dynamics are absent. For weak self-limitation, diversity is small and the equilibrium with non-zero predator abundance is unstable.     

食饵庇护所对斑块环境下Leslie-Gower捕食系统的影响

建立了一个显式含有空间庇护所的两斑块Leslie-Gower捕食者-食饵系统。假设只有食饵种群在斑块间以常数迁移率迁移,且在每个斑块上食饵间的迁移比局部捕食者-食饵相互作用发生的时间尺度要快。利用两个时间尺度,可以构建用来描述所有斑块总的食饵和捕食者密度的综合系统。数学分析表明,在一定条件下,存在唯一的正平衡点,并且此平衡点全局稳定。进一步,捕食者的数量随着食饵庇护所数量增加而降低;在一定条件下,食饵的数量随着食饵庇护所数量增加先增加后降低,在足够强的庇护所强度下,两物种出现灭绝。对比以往研究,利用显式含有和隐含空间庇护所的数学模型所得结论不一致,这意味着在研究庇护所对捕食系统种群动态影响时,空间结构可能起着重要作用。     

具有时滞和密度制约基于比率的捕食-被捕食系统的稳定性

考虑了一类具有时滞和密度制约基于比率的捕食-被捕食系统的稳定性.证明了系统在一定条件下的一致持久性,并且得到了系统正平衡位置局部稳定和全局稳定的充分条件     

Some results on global stability of a predator-prey system

In this paper we derive some results to ensure the global stability of a predator-prey system. The results cover most of the models which have been proposed in the ecological literature for predator-prey systems. The first result is very geometric and it is very easy to check from the graph of prey and predator isoclines. The second one is purely algebraic, however, it covers the defects of the first one especially in dealing with Holling's type-3 functional response in some sense. We also discuss the global stability of Kolmogorov's model. Some examples are presented in the discussion section.Works partially supported by the National Science Council of the Republic of China