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.     

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

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

Dynamics of a stochastic predator-prey model with habitat complexity and prey aggregation

Interplay between predator and prey is a complex process in ecosystems due to its nature. The population dynamics can be affected by many extrinsic and intrinsic factors. In this paper, we make an attempt to uncover the effects from environmental disturbances when populations are subject to habitat complexity and aggregation effect. We firstly propose a stochastic predator-prey model with habitat complexity and aggregation efficiency for prey. We then mathematically analyze the model, to demonstrate the existence, uniqueness and the stochastically ultimately boundedness of the global positive solution, and to establish sufficient conditions for the existence of ergodic stationary distribution of the solution. We also establish sufficient conditions under which either only predator population dies out or the entire predator-prey model becomes extinct. Our theoretical and numerical results indicate that: (1) the environmental noises are disadvantage for the survival of biological populations; (2) when the density of prey is greater than one, prey aggregation can heighten the capability of predator species to capture prey and reduce the effect of environmental fluctuations, while when the density of prey is less than one, the results are opposite; (3) habitat complexity is propitious to the survival of prey population and may seriously threaten the persistence of the predator population.     

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.     

Sigmoid functional responses and population stability

Sigmoid functional responses are known to stabilize the differential Lotka-Volterra predator-prey model. However, we have found that they have no such effect in a comparable discrete generation model. The difficulty in stabilizing this model results from the one-generation time delay between changes in predator population density and the level of prey mortality. By contrast, sigmoid functional responses can stabilize the system if the predator population remains relatively constant, as is more likely of generalist predators.     

Predation across spatial scales in heterogeneous environments

A hierarchy of scales is introduced to the spatially heterogeneous Lotka-Volterra predator-prey diffusion model, and its effects on the model's spatial and temporal behavior are studied. When predators move on a large scale relative to prey, local coupling of the predator-prey interaction is replaced by global coupling. Prey with low dispersal ability become narrowly confined to the most productive habitats, strongly amplifying the underlying spatial pattern of the environment. As prey diffusion rate increases, the prey distribution spreads out and predator abundance declines. The model retains neutrally stable Lotka-Volterra temporal dynamics: different scales of predator and prey dispersal do not stabilize the interaction. The model predicts that, for prey populations that are limited by widely ranging predators, species with low dispersal ability should be restricted to discrete high density patches, and those with greater mobility should be more uniformly distributed at lower density.     

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.     

The Lotka-Volterra predator-prey model with foraging-predation risk trade-offs

This article studies the effects of adaptive changes in predator and/or prey activities on the Lotka-Volterra predator-prey population dynamics. The model assumes the classical foraging-predation risk trade-offs: increased activity increases population growth rate, but it also increases mortality rate. The model considers three scenarios: prey only are adaptive, predators only are adaptive, and both species are adaptive. Under all these scenarios, the neutral stability of the classical Lotka-Volterra model is partially lost because the amplitude of maximum oscillation in species numbers is bounded, and the bound is independent of the initial population numbers. Moreover, if both prey and predators behave adaptively, the neutral stability can be completely lost, and a globally stable equilibrium would appear. This is because prey and/or predator switching leads to a piecewise constant prey (predator) isocline with a vertical (horizontal) part that limits the amplitude of oscillations in prey and predator numbers, exactly as suggested by Rosenzweig and MacArthur in their seminal work on graphical stability analysis of predator-prey systems. Prey and predator activities in a long-term run are calculated explicitly. This article shows that predictions based on short-term behavioral experiments may not correspond to long-term predictions when population dynamics are considered.     

Local spatial structure and predator-prey dynamics: counterintuitive effects of prey enrichment

The Lotka-Volterra predator-prey model with prey density dependence shows the final prey density to be independent of its vital rates. This result assumes the community to be well mixed so that encounters between predators and prey occur as a product of the landscape densities, yet empirical evidence suggests that over small spatial scales this may not be the normal pattern. Starting from an individual-based model with neighborhood interactions and movements, a deterministic approximation is derived, and the effect of local spatial structure on equilibrium densities is investigated. Incorporating local movements and local interactions has important consequences for the community dynamics. Now the final prey density is very much dependent on its birth, death, and movement rates and in ways that seem counterintuitive. Increasing prey fecundity or mobility and decreasing the coefficient of competition can all lead to decreases in the final density of prey if the predator is also relatively immobile. However, analysis of the deterministic approximation makes the mechanism for these results clear; each of these changes subtly alters the emergent spatial structure, leading to an increase in the predator-prey spatial covariance at short distances and hence to a higher predation pressure on the prey.     

Effects of Density Dependent Migrations on the Dynamics of a Predator Prey Model

We study the effects of density dependent migrations on the stability of a predator-prey model in a patchy environment which is composed with two sites connected by migration. The two patches are different. On the first patch, preys can find resource but can be captured by predators. The second patch is a refuge for the prey and thus predators do not have access to this patch. We assume a repulsive effect of predator on prey on the resource patch. Therefore, when the predator density is large on that patch, preys are more likely to leave it to return to the refuge. We consider two models. In the first model, preys leave the refuge to go to the resource patch at constant migration rates. In the second model, preys are assumed to be in competition for the resource and leave the refuge to the resource patch according to the prey density. We assume two different time scales, a fast time scale for migration and a slow time scale for population growth, mortality and predation. We take advantage of the two time scales to apply aggregation of variables methods and to obtain a reduced model governing the total prey and predator densities. In the case of the first model, we show that the repulsive effect of predator on prey has a stabilizing effect on the predator-prey community. In the case of the second model, we show that there exists a window for the prey proportion on the resource patch to ensure stability.     

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

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

Bacterial Predator-Prey Interaction at Low Prey Density

A bacterial predator-prey interaction was studied using Bdellovibrio and bioluminescent prey bacteria. The attacking bdellovibrio causes decay of bioluminescence, which is correlated with bdellovibrio penetration into the prey. The behavior of the prey and predator populations over time was found to be well described by a Lotka-Volterra model. By using this model, the probability of bdellovibrio penetration after encountering a prey cell was found to be approximately 3.0%. The prey density required to give the bdellovibrios a 50% chance of survival was calculated to be at least 3.0 × 10⁶ cells per ml, and the density required for population equilibria was calculated to be about 7 × 10⁵ prey bacteria per ml. These values, not generally characteristic of natural habitats, suggest that the existence of Bdellovibrio in nature is limited to special ecological niches.     

Examination of the effects of toxicity and nutrition on a two-prey one-predator system with a metabolomics-inspired model

Mercury (Hg) sequestration by phytoplankton results in intracellular concentrations that are multiple times greater than ambient water levels, and therefore the consumption of contaminated phytoplankton by herbivorous zooplankton, such as Daphnia, and their inefficient excretion of methylmercury (MeHg) can mediate the transfer to higher trophic levels. Employing a modified version of a metabolomics-inspired Daphnia ecophysiological model, the present study introduces two prey species to a simple Lotka-Volterra predator-prey system in order to shed light on the implications for the integrity of zooplankton assemblages, when experiencing multiple prey items of different toxicity and nutritional quality. We also examine the capacity of adaptive strategies of the predator (homeostatic rigidity, energetic investments to cope with toxicity) to shape predator-prey interactions. Our analysis suggests that the degree of nutritional quality of the prey items is a predominant driver of the predator-prey relationships, shifting from prey- to predator-dominated food webs with increasing nutritional quality. Increasing prey nutritional content leads to the emergence of oscillatory behaviour, which can be further modulated by the growth rates and degree of toxicity of different prey species. Severe exposure to contamination could lead to a decline of the predator biomass with faster growth rates of low nutritional quality prey, even though the increase of its MeHg somatic quota is only modest. In stark contrast, when a prey assemblage of superior nutritional quality prevails in an environment of elevated toxicity, faster prey growth rates are conducive to higher predator biomass levels, albeit its distinctly higher internal contaminant content. Owing to the heightened somatic growth dilution, the ingestion of carbon and nutritional metabolites is significantly higher relative to the MeHg intake rates, which leads to faster net growth of the predator and thus reinforces the benefits brought about by the nutritional value of their diet. Our results suggest that the homeostatic rigidity of the predator can assist in coping with toxic exposure. With a tighter range between the minimum and optimum somatic quotas, the predator population appears to be more resilient and its decline begins at higher levels of MeHg exposure. The predator-prey system displays a greater propensity for oscillatory behaviour, with their amplitude being driven by the interplay between the degree of saturation for nutritionally beneficial metabolites, and the energetic investments allotted to cope with toxicity and/or the excretion of excess metabolic by-products. We conclude by highlighting the prospect of our modelling work to guide new directions of research, to test a multitude of hypotheses pertaining to various ecophysiological facets of predator-prey systems, and extend its use to other contexts, such as the implications of toxin-producing algae for the predator physiology.     

Cannibalism in an age-structured predator-prey system

Recently, Kohlmeier and Ebenhöh showed that cannibalism can stabilize population cycles in a Lotka-Volterra type predator-prey model. Population cycles in their model are due to the interaction between logistic population growth of the prey and a hyperbolic functional response. In this paper, we consider a predator-prey system where cyclic population fluctuations are due to the age structure in the predator species. It is shown that cannibalism is also a stabilizing mechanism when population oscillations are due to this age structure. We conclude that in predator-prey systems, cannibalism by predators can stabilize both externally generated (consumer-resource) as well as internally generated (agestructure) fluctuations.     

A powerful truncated tail strength method for testing multiple null hypotheses in one dataset

This article re-analyses a prey-predator model with a refuge introduced by one of the founders of population ecology Gause and his co-workers to explain discrepancies between their observations and predictions of the Lotka-Volterra prey-predator model. They replaced the linear functional response used by Lotka and Volterra by a saturating functional response with a discontinuity at a critical prey density. At concentrations below this critical density prey were effectively in a refuge while at a higher densities they were available to predators. Thus, their functional response was of the Holling type III. They analyzed this model and predicted existence of a limit cycle in predator-prey dynamics. In this article I show that their model is ill posed, because trajectories are not well defined. Using the Filippov method, I define and analyze solutions of the Gause model. I show that depending on parameter values, there are three possibilities: (1) trajectories converge to a limit cycle, as predicted by Gause, (2) trajectories converge to an equilibrium, or (3) the prey population escapes predator control and grows to infinity.     

The stability of the Boubaker polynomials expansion scheme (BPES)-based solution to Lotka-Volterra problem

This comment addresses critics on the claimed stability of solution to the accelerated-predator-satiety Lotka-Volterra predator-prey problem, proposed by Dubey al. (2010. A solution to the accelerated-predator-satiety Lotka-Volterra predator-prey problem using Boubaker polynomial expansion scheme. Journal of Theoretical Biology 264, 154-160). Critics are based on incompatibilities between the claimed asymptotic behavior and the presumed Malthusian growth of prey population in absence of predator.     

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.     

Explicit model for searching behavior of predator

The authors present an approach for explicit modeling of spatio-temporal dynamics of predator-prey community. This approach is based on a reaction-diffusion-adjection PD (prey dependent) system. Local kinetics of population is determined by logistic reproduction function of prey, constant natural mortality of predator and Holling type 2 trophic function. Searching behavior of predator is described by the advective term in predator balance equation assuming the predator acceleration to be proportional to the prey density gradient. The model was studied with zero-flux boundary conditions. The influence of predator searching activity on the community dynamics, in particular, on the emergence of spatial heterogeneity, has been investigated by linear analysis and numerical simulations. It has been shown how searching activity may effect the persistence of species, stabilizing predator-prey interactions at very low level of pest density. It has been demonstrated that obtaining of such dynamic regimes does not require the use of complex trophic functions.     

Cannibalism in discrete-time predator-prey systems

In this study, we propose and investigate a two-stage population model with cannibalism. It is shown that cannibalism can destabilize and lower the magnitude of the interior steady state. However, it is proved that cannibalism has no effect on the persistence of the population. Based on this model, we study two systems of predator-prey interactions where the prey population is cannibalistic. A sufficient condition based on the nontrivial boundary steady state for which both populations can coexist is derived. It is found via numerical simulations that introduction of the predator population may either stabilize or destabilize the prey dynamics, depending on cannibalism coefficients and other vital parameters.     

A functional response model of a predator population foraging in a patchy habitat

1. Functional response models (e.g. Holling's disc equation) that do not take the spatial distributions of prey and predators into account are likely to produce biased estimates of predation rates. 2. To investigate the consequences of ignoring prey distribution and predator aggregation, a general analytical model of a predator population occupying a patchy environment with a single species of prey is developed. 3. The model includes the density and the spatial distribution of the prey population, the aggregative response of the predators and their mutual interference. 4. The model provides explicit solutions to a number of scenarios that can be independently combined: the prey has an even, random or clumped distribution, and the predators show a convex, sigmoid, linear or no aggregative response. 5. The model is parameterized with data from an acarine predator-prey system consisting of Phytoseiulus persimis and Tetranychus urticae inhabiting greenhouse cucumbers. 6. The model fits empirical data quite well and much better than if prey and predators were assumed to be evenly distributed among patches, or if the predators were distributed independently of the prey. 7. The analyses show that if the predators do not show an aggregative response it will always be an advantage to the prey to adopt a patchy distribution. On the other hand, if the predators are capable of responding to the distribution of prey, then it will be an advantage to the prey to be evenly distributed when its density is low and switch to a more patchy distribution when its density increases. The effect of mutual interference is negligible unless predator density is very high. 8. The model shows that prey patchiness and predator aggregation in combination can change the functional response at the population level from type II to type III, indicating that these factors may contribute to stabilization of predator-prey dynamics.