Prey-taxis destabilizes homogeneous stationary state in spatial Gause–Kolmogorov-type model for predator–prey system

We consider a continuous taxis-diffusion-reaction system of partial-differential equations describing spatiotemporal dynamics of a predator–prey system. The local kinetics of the system is defined by general Gause–Kolmogorov-type model. The predator ability to pursue the prey is modelled by the Patlak–Keller–Segel taxis model, assuming that movement velocities of predators are proportional to the gradients of specific cues emitted by prey (e.g., odour, pheromones, exometabolites). The linear stability analysis of the model showed that the non-trivial homogeneous stationary regime of the model becomes unstable with respect to small heterogeneous perturbations with increase of prey-taxis activity; an Andronov–Hopf bifurcation occurs in the system when the taxis coefficient of predator exceeds its critical bifurcation value that exists for all admissible values of model parameters. These findings generalize earlier results obtained for particular cases of the Gause–Kolmogorov-type model assuming logistic reproduction of the prey population and the Holling types I and II functional responses of the predator population. Numerical simulations with theta-logistic growth of the prey population and the Ivlev functional response of predators illustrate and support results of the analytical study.     

Complexity in a prey-predator model with prey refuge and diffusion

Prey-predator interaction is one of the most commonly observed relationships in ecosystem. In the study of prey-predator models, it is frequently assumed that the changes in population densities are only time-dependent and the dynamics is generally represented by coupled nonlinear ordinary differential equations. In natural system, however, either prey or predator or both move from one place to another for various reasons. In such a case, their dynamic interaction depends both on time and space and requires coupled nonlinear partial differential equations for its dynamic representation. It is also well documented that prey refuges affect the interaction between prey and predator significantly. In this paper, we studied the dynamics of a diffusive prey-predator interaction with prey refuge and type III response function. We have considered both one and two dimensional diffusivity in the model system and presented different stability results under the assumptions that one or both species may be mobile or sedentary. Our results showed that the system may exhibit different spatiotemporal (non-Turing) patterns, like spiral waves, patchy structures, spot pattern, or even spatiotemporal chaos depending on the refuge availability and diffusion rate of species. Another interesting finding was that the dynamic complexity in a prey-predator model increases in case of mobile predator and sedentary prey compare to mobile prey and sedentary predator while refuge availability is varied.     

Dynamical analysis of a prey-predator model incorporating a prey refuge with variable carrying capacity

A prey-predator model incorporating prey refuge with variable carrying capacity and Holling type-II functional response is proposed and analyzed. The model includes a case of increasing carrying capacity as well as a decreasing carrying capacity case. Sufficient conditions are derived to ensure the existence and local stability of the equilibrium points of the proposed model. Moreover, the occurrence of transcritical bifurcation as well as Hopf bifurcation are investigated. The effect of some model parameter related to the prey refuge and the variable carrying capacity on the prey-predator dynamics has been examined. Numerical simulations are presented to demonstrate the theoretical results and to illustrate the effect of these parameters on the model dynamics. Moreover, a comparison with the constant carrying case has been presented.     

一类被开发的HollingⅢ类功能反应模型的定性分析

本文研究了一类捕食种群、食饵种群同时具有收获率的HollingⅢ类功能反应生态系统,其中食饵种群具有非线性密度制约,捕食者无密度制约.应用微分方程定性理论讨论了系统的平衡点,分析了中心焦点的阶数以及稳定性,所给定参数满足一定条件时系统不存在极限环,最后根据细焦点的稳定性判断出极限环的存在性,并验证了极限环的惟一性.     

Multiple Limit Cycles in a Gause Type Predator–Prey Model with Holling Type III Functional Response and Allee Effect on Prey

This work aims to examine the global behavior of a Gause type predator–prey model considering two aspects: (i) the functional response is Holling type III and, (ii) the prey growth is affected by the Allee effect. We prove the origin of the system is an attractor equilibrium point for all parameter values. It has also been shown that it is the ω-limit of a wide set of trajectories of the system, due to the existence of a separatrix curve determined by the stable manifold of the equilibrium point (m,0), which is associated to the Allee effect on prey. When a weak Allee effect on the prey is assumed, an important result is obtained, involving the existence of two limit cycles surrounding a unique positive equilibrium point: the innermost cycle is unstable and the outermost stable. This property, not yet reported in models considering a sigmoid functional response, is an important aspect for ecologists to acknowledge as regards the kind of tristability shown here: (1) the origin; (2) an interior equilibrium; and (3) a limit cycle of large amplitude. These models have undoubtedly been rather sensitive to disturbances and require careful management in applied conservation and renewable resource contexts.     

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.     

Bioeconomic harvesting of a prey-predator fishery

This paper deals with the problem of non-selective harvesting of a prey-predator system by using a reasonable catch-rate function instead of usual catch-per-unit-efforthypothesis. Here both the prey and the predator species obey the law of logistic growth. We have taken the predator functional response to prey density in such a form that each predator's functional response to the prey density approaches a constant as the prey population increases. Boundedness of the exploited system is examined. The existence of its steady states and their stability (local and global) are studied using Eigenvalue analysis. The existence of bionomic equilibria has been illustrated using a numerical example. The problem of determining the optimal harvesting policy is then solved by using Pontryagin's maximum principle.     

On a predator–prey system of Gause type

In this paper a Gause type model of interactions between predator and prey population is considered. We deal with the sufficient condition due to Kuang and Freedman in the generalized form including a kind of weight function. In a previous paper we proved that the existence of such weight function implies the uniqueness of limit cycle. In the present paper we give a new condition equivalent to the existence of a weight function (Theorem 4.4). As a consequence of our result, it is shown that some simple qualitative properties of the trophic function and the prey isocline ensure the uniqueness of limit cycle.     

带有食饵保护的Gause型扩散捕食系统的稳定性分析

考虑了具有扩散项和食饵保护的Gause型捕食系统.该模型带有齐次Neumann边界条件.讨论了系统的全局吸引性以及系统非负常数平衡态的局部稳定性和全局稳定性.其条件依赖于食饵保护参数,表明了食饵保护对系统动力学行为的影响.     

Consequences of refuge for the functional response of Dermestes ater (Coleoptera: Dermestidae) to Musca domestica (Diptera: Muscidae)

It is well known that a predator has the potential to regulate a prey population only if the predator responds to increases in prey density and inflicts greater mortality rates. Predators may cause such density-dependent mortality depending on the nature of the functional and numerical responses. Yet, few studies have examined the relationship between the addition of refuges and the characteristic of functional response fits. We investigated whether addition of a refuge changed the type of functional response exhibited by Dermestes ater on Musca domestica, comparing the inherent ability of D. ater to kill houseflies in the absence and in the presence of refuge. An additional laboratory experiment was also carried out to assess handling and searching times exhibited by D. ater. Logistic regression analyses revealed a type III functional response for predator–prey interaction without refuge, and results were described by the random predator equation. The mean number of prey killed did not differ between experimental habitats, indicating that the addition of refuge did not inhibit predation. However, predators that interacted with prey without refuge spent less time searching for prey at higher densities, increasing predatory interaction. We concluded that this interaction may be weak, because data from experiments with refuge fitted poorly to models. However, the high variability and the nonsignificance of the data from the experiment with refuge show the importance of refuge for promoting spatial heterogeneity, which may prevent prey extinction.     

The role of virus infection in a simple phytoplankton zooplankton system

Many planktonic species show spectacular bursts ("blooms") in population density. Though viral infections are known to cause behavioural and other changes in phytoplankton and other aquatic species, yet their role in regulating the phytoplankton population is still far from being understood. To study the role of viral diseases in the planktonic species, we model the phytoplankton-zooplankton system as a prey-predator system. Here the prey (phytoplankton) species is infected with a viral disease that divides the prey population into susceptible and infected classes, with the infected prey being more vulnerable to predation by the predator (zooplankton). The dynamical behaviour of the system is investigated from the point of view of stability and persistence both analytically and numerically. The model shows that infection can be sustained only above a threshold of force of infection, and, there exists a range in the infection rate where this system shows "bloom"-like stable limit cycle oscillations. The time series of natural "blooms" with different types of irregular oscillations can arise in this model simply from a biologically realistic feature, i.e., by the random variation of the epidemiological parameter (rate of infection) in the infected prey population. The difference in mean strength of infection alone can lead to the different types of patterns observed in natural planktonic blooms.     

具有捕食正效应的捕食-食饵系统

在经典的捕食食饵系统中考虑到由于捕食效应对食饵种群带来的正向调节作用后,提出了具有捕食正效应的捕食-食饵系统.通过对模型的动力学行为的分析,从理论上说明了正向调节作用对系统的影响,并就第一象限内平衡点存在时的相图解释了捕食正效应的作用.结果表明:(1)捕食系统中适当的正向调节作用会增加系统的稳定性;(2)当捕食正效应达到一定的程度后系统拥有一个不稳定的极限环;(3)当捕食正效应过大时会使系统的稳定性发生变化,使捕食者种群与食饵种群同时趋向无穷,出现了调节放纵现象.这些结果在保护生物学中具有重要的意义.     

Population dynamics of thrips prey and their mite predators in a refuge

Prey refuges are expected to affect population dynamics, but direct experimental tests of this hypothesis are scarce. Larvae of western flower thrips Frankliniella occidentalis use the web produced by spider mites as a refuge from predation by the predatory mite Neoseiulus cucumeris. Thrips incur a cost of using the refuge through reduced food quality within the web due to spider mite herbivory, resulting in a reduction of thrips developmental rate. These individual costs and benefits of refuge use were incorporated in a stage-structured predator-prey model developed for this system. The model predicted higher thrips numbers in presence than in absence of the refuge during the initial phase. A greenhouse experiment was carried out to test this prediction: the dynamics of thrips and their predators was followed on plants damaged by spider mites, either with or without web. Thrips densities in presence of predators were higher on plants with web than on unwebbed plants after 3 weeks. Experimental data fitted model predictions, indicating that individual-level measurements of refuge costs and benefits can be extrapolated to the level of interacting populations. Model-derived calculations of thrips population growth rate enable the estimation of the minimum predator density at which thrips benefit from using the web as a refuge. The model also predicted a minor effect of the refuge on the prey density at equilibrium, indicating that the effect of refuges on population dynamics hinges on the temporal scale considered.     

Age structure effects in predator-prey interactions

A mathematical model of predator-prey interactions is proposed which incorporates both age structure in the predators and density dependence in the prey. The properties of the model are investigated by a linearized analysis, which enables the conditions for stability to be formulated. The analysis indicates that for a substantial domain of parameter space, a stable equilibrium is possible with the prey well below its carrying capacity. The effect of violating the stability conditions on the behaviour of the model was investigated by computer simulation. Two further types of behaviour are noted in which coexistence is possible. The first is a two point limit cycle in which young and old predators occur in alternate time periods. The second involves a limit cycle in which the component population trajectories lie on closed curves in phase space.     

Predator size interacts with habitat structure to determine the allometric scaling of the functional response

While both predator body size and prey refuge provided by habitat structure have been established as major factors influencing the functional response (per capita consumption rate as a function of prey density), potential interactions between these factors have rarely been explored. Using a crab predator (Panopeus herbstii) – mussel prey (Brachidontes exustus) system, we examined the allometric scaling of the functional response in oyster (Crassostrea virginica) reef habitat, where crevices within oyster clusters provide mussels refuge from predation. A field survey of mussel distribution showed that mussels attach closer to the cluster periphery at high mussel density, indicating the potential for saturation of the refuge. In functional response experiments, the consumption rate of large crabs was depressed at low prey density relative to small crabs, while at high prey density the reverse was true. Specifically, the attack rate coefficient and handling time both decreased non‐linearly with crab size. An additional manipulation revealed that at low prey densities, the ability of large crabs to maneuver their claws and bodies to extract mussels from crevices was inhibited relative to small crabs by the structured habitat, reducing their attack rate. At high prey densities, crevices were saturated, forcing mussels to the edge of clusters where crabs were only limited by handling time. Our study illuminates a potentially general mechanism where the quality of the prey refuge provided by habitat structure is dependent on the relative size of the predator. Thus anthropogenic influences that alter the natural crab size distribution or degrade reef habitat structure could threaten the long‐term stability of the crab –mussel interaction in reefs.     

Natural selection and population density-feedback II. The evolution of functional response curves

The nature of the functional response may be qualitatively understood as follows. Sigmoid responses to one food type may arise, in the presence of alternate foods, as a result of optimal feeding and foraging behavior. Sigmoid curves resulting from this cause I term class A curves. The same curve may also arise in the absence of alternate foods as a result of learning, individual variations in the level of food density at which predators begin feeding, or training effects. The latter I have termed class B curves. At very high food densities, a drop in food intake per predator might occur because of the tendency for predators to take easily found and captured items first and to become more selective when food is very common. Such “dome-shaped” curves have been found in the laboratory but should be rare in nature. Computer simulation of a three trophic-level system, using the phenotypic selection model of Emlen, indicates that natural selection acting on prey should encourage sigmoidality in the predator's class B functional response, at least in disturbed environments. The opposite force arises from selection acting on predators. However, given the magnitudes of growth efficiencies (see Eq. (8), it appears that at least for terrestrial vertebrates, selection on prey species is more important than selection on predators for determining functional responses. Accordingly, prey-predator systems occupying highly variable environments are expected to show more marked type III (class B) curves than systems in more stable areas. Finally, the role of functional response for prey-predator stability is discussed. Class A (alternate food) responses may result in population control for prey in multiple prey systems. Peterman and Pikitch have modeled systems in which type III functional response by predators, in systems where predation varies independently of prey, may lead to double equilibria. This picture is clouded, however, when predator populations are interactive with their food, though double equilibria are still possible (J. M. Emlen, 1984, “Population Biology: The Coevolution of Population Dynamics and Behavior,” Macmillan, New York, in press).     

The Effect of Refuge on Dermestes ater (Coleoptera: Dermestidae) Predation on Musca domestica (Diptera: Muscidae): Refuge for Prey or the Predator?

It is believed that habitat heterogeneity can change the extent of predator-prey interactions. Therefore, in this study we examined the effect of habitat heterogeneity (characterized here as an addition of refuge) on D. ater predation on M. domestica. Predation of D. ater on M. domestica larvae was carried out in experimental habitats with and without refuge, and examined at different prey densities. The number of prey eaten by beetles over 24 h of predator-prey interaction was recorded, and we investigated the strength of interaction between prey and predator in both experimental habitats by determining predator functional response. The mean number of prey eaten by beetles in the presence of refuge was significantly higher than in the absence of refuge. Females had greater weight gains than males. Logistic regression analyses revealed the type II functional response for both experimental habitats, even though data did not fit well into the random predator model. Results suggest that the addition of refuge in fact enhanced predation, as prey consumption increased in the presence of refuge. Predators kept in the presence of refuge also consumed more prey at high prey densities. Thus, we concluded that the addition of refuge was an important component mediating D. ater-M. domestica population interactions. Refuge actually acted as a refuge for predators from prey, since prey behaviors detrimental to predators were reduced in this case.     

含有食饵自庇护的食饵-捕食模型的模糊滑模控制

针对一类令布食饵自虞护的食饵-捕食系统,当环境变化时,系统参数的扰动处于动态状态,利用其扰动边界的特性,通过设计模糊滑模控制器将失控系统引至平稳轨道,并使系统全局稳定,从理论上实现了失控系统的有效控制,为防治食饵种群灭绝,引起食物链断裂,导致生态平衡破坏提供了一定的理论依据．最后的仿真结果表明了此方法在理论上的有效性．     

Should "handled" prey be considered? Some consequences for functional response, predator-prey dynamics and optimal foraging theory

Predator-prey models consider those prey that are free. They assume that once a prey is captured by a predator it leaves the system. A question arises whether in predator-prey population models the variable describing prey population shall consider only those prey which are free, or both free and handled prey together. In the latter case prey leave the system after they have been handled. The classical Holling type II functional response was derived with respect to free prey. In this article we derive a functional response with respect to prey density which considers also handled prey. This functional response depends on predator density, i.e., it accounts naturally for interference. We study consequences of this functional response for stability of a simple predator-prey model and for optimal foraging theory. We show that, qualitatively, the population dynamics are similar regardless of whether we consider only free or free and handled prey. However, the latter case may change predictions in some other cases. We document this for optimal foraging theory where the functional response which considers both free and handled prey leads to partial preferences which are not observed when only free prey are considered.     

You can’t run but you can hide: refuge use in frog tadpoles elicits density-dependent predation by dragonfly larvae

The potential role of prey refuges in stabilizing predator–prey interactions is of longstanding interest to ecologists, but mechanisms underlying a sigmoidal predator functional response remain to be fully elucidated. Authors have disagreed on whether the stabilizing effect of prey refuges is driven by prey- versus predator-centric mechanisms, but to date few studies have married predator and prey behavioural observations to distinguish between these possibilities. We used a dragonfly nymph–tadpole system to study the effect of a structural refuge (leaf litter) on the predator’s functional response, and paired this with behavioural observations of both predator and prey. Our study confirmed that hyperbolic (type II) functional responses were characteristic of foraging predators when structural cover was low or absent, whereas the functional response was sigmoidal (type III) when prey were provided with sufficient refuge. Prey activity and refuge use were density independent across cover treatments, thereby eliminating a prey-centric mechanism as being the genesis for density-dependent predation. In contrast, the predator’s pursuit length, capture success, and handling time were altered by the amount of structure implying that observed shifts in density-dependent predation likely were related to predator hunting efficiency. Our study advances current theory by revealing that despite fixed-proportion refuge use by prey, presence of a prey refuge can induce density-dependent predation through its effect on predator hunting strategy. Ultimately, responses of predator foraging decisions in response to changes in prey availability and search efficiency may be more important in producing density-dependent predation than the form of prey refuge use.