Interactions of Turing and Hopf bifurcations in an additional food provided diffusive predator-prey model

Complex spatiotemporal dynamics of a diffusive predator-prey system involving additional food supply to predator and intra-specific competition among predator, are investigated. We establish critical conditions of the occurrence of Turing instability, which are necessary and sufficient. Furthermore, we also establish conditions of the occurrence of codimension-2 Turing-Hopf bifurcation and Turing-Turing bifurcation, by exploring interactions of Turing bifurcations and Hopf bifurcation. For Turing-Hopf bifurcation, by analyzing normal form truncated to order 3 which are derived by applying normal form method, it is shown that under proper conditions, diffusive predator-prey system generates interesting spatial, temporal and spatiotemporal patterns, including a pair of spatially inhomogeneous steady states, a spatially homogeneous periodic solution and a pair of spatially inhomogeneous periodic solutions. And numerical simulations are also shown to support theory analysis. Moreover, it is found that proper intra-specific competition among predator helps generate complex spatiotemporal dynamics. And, proper additional food supply to predator helps control the population fluctuations of predator and prey, while large quantity and high quality of additional food supply will lead to the extinction of prey and make predator change the source of food, which finally destroy the ecosystem. These investigations might help understand complex spatiotemporal dynamics of predator-prey system involving additional food supply to predator and intra-specific competition among predator, and help conserve species in an ecosystem via supplying suitable additional food.     

Pattern formation and control of spatiotemporal chaos in a reaction diffusion prey–predator system supplying additional food

Spatiotemporal dynamics of a predator–prey system in presence of spatial diffusion is investigated in presence of additional food exists for predators. Conditions for stability of Hopf as well as Turing patterns in a spatial domain are determined by making use of the linear stability analysis. Impact of additional food is clear from these conditions. Numerical simulation results are presented in order to validate the analytical findings. Finally numerical simulations are carried out around the steady state under zero flux boundary conditions. With the help of numerical simulations, the different types of spatial patterns (including stationary spatial pattern, oscillatory pattern, and spatiotemporal chaos) are identified in this diffusive predator–prey system in presence of additional food, depending on the quantity, quality of the additional food and the spatial domain and other parameters of the model. The key observation is that spatiotemporal chaos can be controlled supplying suitable additional food to predator. These investigations may be useful to understand complex spatiotemporal dynamics of population dynamical models in presence of additional food.     

Dynamical consequences of predator interference in a tri-trophic model food chain

A model food chain involving a specialist and a generalist predator is proposed and studied. One of the salient features of this model food chain is that it combines both the schemes (Volterra and Leslie) of modeling predator–prey interaction in one system in such a way that the demerits of these individual formulations are suppressed and the resulting model system represents a common unit of real world food webs. The stability analysis of the proposed model is carried out. The Hopf bifurcation conditions of the positive equilibrium point are established. Our numerical computations show that chaotic dynamics is sensitive to changes in values of parameters measuring attributes of either interacting populations or their environments. Two dimensional parameter scans suggest that the model food chain displays short-term recurrent chaos. This can be regarded as a plausible explanation for why it has been so difficult to detect deterministic chaos in natural populations.     

Novel dynamics of a predator–prey system with harvesting of the predator guided by its population

A predator–prey model was extended to include nonlinear harvesting of the predator guided by its population, such that harvesting is only implemented if the predator population exceeds an economic threshold. The proposed model is a nonsmooth dynamic system with switches between the original predator-prey model (free subsystem) and a model with nonlinear harvesting (harvesting subsystem). We initially examine the dynamics of both the free and the harvesting subsystems, and then we investigate the dynamics of the switching system using theories of nonsmooth systems. Theoretical results showed that the harvesting subsystem undergoes multiple bifurcations, including saddle-node, supercritical Hopf, Bogdanov–Takens and homoclinic bifurcations. The switching system not only retains all of the complex dynamics of the harvesting system but also exhibits much richer dynamics such as a sliding equilibrium, sliding cycle, boundary node (saddle point) bifurcation, boundary saddle-node bifurcation and buckling bifurcation. Both theoretical and numerical results showed that, by implementing predator population guided harvesting, the predator and prey population could coexist in more scenarios than those in which the predator may go extinct for the continuous harvesting regime. They could either stabilize at an equilibrium or oscillate periodically depending on the value of the economic threshold and the initial value of the system.     

Complex dynamics of an eco-epidemiological model with different competition coefficients and weak Allee in the predator

The paper explores an eco-epidemiological model with weak Allee in predator, and the disease in the prey population. We consider a predator-prey model with type II functional response. The curiosity of this paper is to consider different competition coefficients within the prey population, which leads to the emergent carrying capacity. We perform the local and global stability analysis of the equilibrium points and the Hopf bifurcation analysis around the endemic equilibrium point. Further we pay attention to the chaotic dynamics which is produced by disease. Our numerical simulations reveal that the three species eco-epidemiological system without weak-Allee induced chaos from stable focus for increasing the force of infection, whereas in the presence of the weak-Allee effect, it exhibits stable solution. We conclude that chaotic dynamics can be controlled by the Allee parameter as well as the competition coefficients. We apply basic tools of non-linear dynamics such as Poincare section and maximum Lyapunov exponent to identify chaotic behavior of the system.     

A modified Leslie–Gower predator–prey model with ratio‐dependent functional response and alternative food for the predator

In this work, a modified Leslie–Gower predator–prey model is analyzed, considering an alternative food for the predator and a ratio‐dependent functional response to express the species interaction. The system is well defined in the entire first quadrant except at the origin ( 0 , 0 ) . Given the importance of the origin ( 0 , 0 ) as it represents the extinction of both populations, it is convenient to provide a continuous extension of the system to the origin. By changing variables and a time rescaling, we obtain a polynomial differential equations system, which is topologically equivalent to the original one, obtaining that the non‐hyperbolic equilibrium point ( 0 , 0 ) in the new system is a repellor for all parameter values. Therefore, our novel model presents a remarkable difference with other models using ratio‐dependent functional response. We establish conditions on the parameter values for the existence of up to two positive equilibrium points; when this happen, one of them is always a hyperbolic saddle point, and the other can be either an attractor or a repellor surrounded by at least one limit cycle. We also show the existence of a separatrix curve dividing the behavior of the trajectories in the phase plane. Moreover, we establish parameter sets for which a homoclinic curve exits, and we show the existence of saddle‐node bifurcation, Hopf bifurcation, Bogdanov–Takens bifurcation, and homoclinic bifurcation. An important feature in this model is that the prey population can go to extinction; meanwhile, population of predators can survive because of the consumption of alternative food in the absence of prey. In addition, the prey population can attain their carrying capacity level when predators go to extinction. We demonstrate that the solutions are non‐negatives and bounded (dissipativity and permanence of population in many other works). Furthermore, some simulations to reinforce our mathematical results are shown, and we further discuss their ecological meanings. Copyright © 2017 John Wiley & Sons, Ltd.     

The role of top predator interference on the dynamics of a food chain model

In this paper, the effects of top predator interference on the dynamics of a food chain model involving an intermediate and a top predator are considered. It is assumed that the interaction between the prey and intermediate predator follows the Volterra scheme, while that between the top predator and its favorite food depends on Beddington–DeAngelis type of functional response. The boundedness of the system, existence of an attracting set, local and global stability of non-negative equilibrium points are established. Number of the bifurcation and Lyapunov exponent bifurcation diagrams is established. It is observed that, the model has different types of attracting sets including chaos. Moreover, increasing the top predator interference stabilizes the system, while increasing the normalization of the residual reduction in the top predator population destabilizes the system.     

Dynamics of an intraguild predation food web model with strong Allee effect in the basal prey

Since intraguild predation (IGP) is a ubiquitous and important community module in nature and Allee effect has strong impact on population dynamics, in this paper we propose a three-species IGP food web model consisted of the IG predator, IG prey and basal prey, in which the basal prey follows a logistic growth with strong Allee effect. We investigate the local and global dynamics of the model with emphasis on the impact of strong Allee effect. First, positivity and boundedness of solutions are studied. Then existence and stability of the boundary and interior equilibria are presented and the Hopf bifurcation curve at an interior equilibrium is given. The existence of a Hopf bifurcation curve indicates that if competition between the IG prey and IG predator for the basal resource lies below the curve then the interior equilibrium remains stable, while if it lies above the curve then the interior equilibrium loses its stability. In order to explore the impact of Allee effect, the parameter space is classified into sixteen different regions and, in each region, the number of interior equilibria is determined and the corresponding bifurcation diagrams on the Allee threshold are given. The extinction parameter regions of at least one species and the necessary coexistence parameter regions of all three species are provided. In addition, we explore possible dynamical patterns, i.e., the existence of multiple attractors. By theoretical analysis and numerical simulations, we show that the model can have one (i.e. extinction of all species), two (i.e. bi-stability) or three (i.e. tri-stability) attractors. It is also found by simulations that when there exists a unique stable interior equilibrium, the model may generate multiple attracting periodic orbits and the coexistence of all three species is enhanced as the competition between the IG prey and IG predator for the basal resource is close to the Hopf bifurcation curve from below. Our results indicate that the intraguild predation food web model exhibits rich and complex dynamic behaviors and strong Allee effect in the basal prey increases the extinction risk of not only the basal prey but also the IG prey or/and IG predator.     

A Leslie–Gower predator–prey model with disease in prey incorporating a prey refuge

In this paper, a predator–prey Leslie–Gower model with disease in prey has been developed. The total population has been divided into three classes, namely susceptible prey, infected prey and predator population. We have also incorporated an infected prey refuge in the model. We have studied the positivity and boundedness of the solutions of the system and analyzed the existence of various equilibrium points and stability of the system at those equilibrium points. We have also discussed the influence of the infected prey refuge on each population density. It is observed that a Hopf bifurcation may occur about the interior equilibrium taking refuge parameter as bifurcation parameter. Our analytical findings are illustrated through computer simulation using MATLAB, which show the reliability of our model from the eco-epidemiological point of view.     

Complex dynamics of a discrete predator–prey model with the prey subject to the Allee effect

In this paper, complex dynamics of the discrete predator–prey model with the prey subject to the Allee effect are investigated in detail. Firstly, when the prey intrinsic growth rate is not large, the basins of attraction of the equilibrium points of the single population model are given. Secondly, rigorous results on the existence and stability of the equilibrium points of the model are derived, especially, by analyzing the higher order terms, we obtain that the non-hyperbolic extinction equilibrium point is locally asymptotically stable. The existences and bifurcation directions for the flip bifurcation, the Neimark–Sacker bifurcation and codimension-two bifurcations with 1:2 resonance are derived by using the center manifold theorem and the bifurcation theory. We derive that the model only exhibits a supercritical flip bifurcation and it is possible for the model to exhibit a supercritical or subcritical Neimark–Sacker bifurcation at the larger positive equilibrium point. Chaos in the sense of Marotto is proved by analytical methods. Finally, numerical simulations including bifurcation diagrams, phase portraits, sensitivity dependence on the initial values, Lyapunov exponents display new and rich dynamical behaviour. The analytic results and numerical simulations demonstrate that the Allee effect plays a very important role for dynamical behaviour.     

A stage-structured food chain model with stage dependent predation: Existence of codimension one and codimension two bifurcations

Stage-structured predator–prey models exhibit rich and interesting dynamics compared to homogeneous population models. The objective of this paper is to study the bifurcation behavior of stage-structured prey–predator models that admit stage-restricted predation. It is shown that the model with juvenile-only predation exhibits Hopf bifurcation with the growth rate of the adult prey as the bifurcation parameter; also, depending on parameter values, a stable limit cycle will emerge, that is, the bifurcation will be of supercritical nature. On the other hand, the analysis of the model with adult-stage predation shows that the system admits a fold-Hopf bifurcation with the adult growth rate and the predator mortality rate as the two bifurcation parameters. We also demonstrate the existence of a unique limit cycle arising from this codimension-2 bifurcation. These results reveal far richer dynamics compared to models without stage-structure. Numerical simulations are done to support analytical results.     

Weak Allee effect in a predator–prey model involving memory with a hump

In this paper we consider a predator–prey system which has a factor that allows for a reduction in fitness due to declining population sizes, often termed an Allee effect. We study the influence of the weak Allee effect which is included in the prey equation and we determine conditions for the occurrence of Hopf bifurcation. The prey population is limited by the carrying capacity of the environment, and the predator growth rate depends on past quantities of the prey which is represented by a weight function that specifies a moment in the past when the quantity of food is the most important from the point of view of the present growth of the predator. The stability properties of the system and the biological issues of the memory and Allee effect on the coexistence of the two species are studied. Finally we present some simulations to verify the veracity of the analytical conclusions.     

Conjugate coupling in ecosystems: Cross-predation stabilizes food webs

We study the dynamics of two predator–prey systems that are coupled via cross-predation, in which each predator consumes also the other prey. This setup constitutes a model system in which conjugate coupling emerges naturally and denotes the transition from two separate food chains to a food web. We show that cross-predation of a certain strength leads to amplitude death stabilizing the food web in a new equilibrium.     

Stability and bifurcation analysis for a fractional prey–predator scavenger model

In this study, we consider a fractional prey–predator scavenger model as well as harvesting by a predator and scavenger. We prove the positivity and boundedness of the solutions in this system. The model undergoes a Hopf bifurcation around one of the existing equilibria where the conditions are met for the occurrence of a Hopf bifurcation. The results show that chaos disappears in this biological model. We conclude that the fractional system is more stable compared with the classical case and the stability domain can be extended under fractional order. In addition, a suitable amount of prey harvesting and a fractional order derivative can control the chaotic dynamics and stabilize them. We also present an extended numerical simulation to validate the results.     

Delay-induced Bogdanov–Takens bifurcation in a Leslie–Gower predator–prey model with nonmonotonic functional response

This work is concerned with the dynamics of a Leslie–Gower predator–prey model with nonmonotonic functional response near the Bogdanov–Takens bifurcation point. By analyzing the characteristic equation associated with the nonhyperbolic equilibrium, the critical value of the delay inducing the Bogdanov–Takens bifurcation is obtained. In this case, the dynamics near this nonhyperbolic equilibrium can be reduced to the study of the dynamics of the corresponding normal form restricted to the associated two-dimensional center manifold. The bifurcation diagram near the Bogdanov–Takens bifurcation point is drawn according to the obtained normal form. We show that the change of delay can result in heteroclinic orbit, homoclinic orbit and unstable limit cycle.     

A cannibalistic eco-epidemiological model with disease in predator population

In this present article, we propose and analyze a cannibalistic predator–prey model with disease in the predator population. We consider two important factors for the dynamics of predator population. The first one is governed through cannibalistic interaction, and the second one is governed through the disease in the predator population via cannibalism. The local stability analysis of the model system around the biologically feasible equilibria are investigated. We perform global dynamics of the model using Lyapunov functions. We analyze and compare the community structure of the system in terms of ecological and disease basic reproduction numbers. The existence of Hopf bifurcation around the interior steady state is investigated. We also derive the sufficient conditions for the permanence and impermanence of the system. The study reveals that the cannibalism acts as a self-regulatory mechanism and controls the disease transmission among the predators by stabilizing the predator–prey oscillations.     

Dynamical study of fractional order differential equations of predator‐pest models

To explore the impact of pest‐control strategy through a fractional derivative, we consider three predator‐prey systems by simple modification of Rosenzweig‐MacArthur model. First, we consider fractional‐order Rosenzweig‐MacArthur model. Allee threshold phenomena into pest population is considered for the second case. Finally, we consider additional food to the predator and harvesting in prey population. The main objective of the present investigation is to observe which model is most suitable for the pest control. To achieve this goal, we perform the local stability analysis of the equilibrium points and observe the basic dynamical properties of all the systems. We observe fractional‐order system has the ability to stabilize Rosenzweig‐MacArthur model with low pest density from oscillatory state. In the numerical simulations, we focus on the bistable regions of the second and third model, and we also observe the effect of the fractional order α throughout the stability region of the system. For the third model, we observe a saddle‐node bifurcation due to the additional food and Allee effect to the pest densities. Also, we numerically plot two parameter bifurcation diagram with respect to the harvesting parameter and fractional order of the system. We finally conclude that fractional‐order Rosenzweig‐MacArthur model and the modified Rosenzweig‐MacArthur model with additional food for the predator and harvested pest population are more suitable models for the pest management.     

A discrete-time model with non-monotonic functional response and strong Allee effect in prey

In this paper, we investigate the impact of strong Allee effect on the stability of a discrete-time predator–prey model with a non-monotonic functional response. The dynamics of discrete-time predator–prey models with strong Allee effect is studied earlier. But, the mathematical investigations of predator–prey dynamics in discrete-time set up with Holling type-IV functional response and strong Allee effect in prey are lacking. The proposed model supports the coexistence of two steady states, and the mathematical features of the model are analyzed based on local stability and bifurcation theory. By considering the Allee parameter as the bifurcation parameter, we provide sufficient conditions for the flip and the Neimark–Sacker bifurcations. We observe that Allee parameter plays a significant role in the dynamics of the system.     

Bifurcation and Turing pattern formation in a diffusive ratio-dependent predator–prey model with predator harvesting

The ratio-dependent predator–prey model exhibits rich dynamics due to the singularity of the origin. Harvesting in a ratio-dependent predator–prey model is relatively an important research project from both ecological and mathematical points of view. In this paper, we study the temporal, spatial and spatiotemporal dynamics of a ratio-dependent predator–prey diffusive model where the predator population harvest at catch-per-unit-effort hypothesis. For the spatially homogeneous model, we derive conditions for determining the direction of Hopf bifurcation and the stability of the bifurcating periodic solution by the center manifold and the normal form theory. For the reaction–diffusion model, firstly it is shown that Turing (diffusion-driven) instability occurs, which induces spatial inhomogeneous patterns. Then it is demonstrated that the model exhibit Hopf bifurcation which produces temporal inhomogeneous patterns. Finally, the existence and non-existence of positive non-constant steady-state solutions are established. Moreover, numerical simulations are performed to visualize the complex dynamic behavior.     

Global bifurcation of solutions for a predator–prey model with prey‐taxis

We study pattern formations in a predator–prey model with prey‐taxis. It is proved that a branch of nonconstant solutions can bifurcate from the positive equilibrium only when the chemotactic is repulsive. Furthermore, we find the stable bifurcating solutions near the bifurcation point under suitable conditions. Copyright © 2014 John Wiley & Sons, Ltd.