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.     

The chaos and control of a food chain model supplying additional food to top-predator

The control and management of chaotic population is one of the main objectives for constructing mathematical model in ecology today. In this paper, we apply a technique of controlling chaotic predator–prey population dynamics by supplying additional food to top-predator. We formulate a three species predator–prey model supplying additional food to top-predator. Existence conditions and local stability criteria of equilibrium points are determined analytically. Persistence conditions for the system are derived. Global stability conditions of interior equilibrium point is calculated. Theoretical results are verified through numerical simulations. Phase diagram is presented for various quality and quantity of additional food. One parameter bifurcation analysis is done with respect to quality and quantity of additional food separately keeping one of them fixed. Using MATCONT package, we derive the bifurcation scenarios when both the parameters quality and quantity of additional food vary together. We predict the existence of Hopf point (H), limit point (LP) and branch point (BP) in the model for suitable supply of additional food. We have computed the regions of different dynamical behaviour in the quantity–quality parametric plane. From our study we conclude that chaotic population dynamics of predator prey system can be controlled to obtain regular population dynamics only by supplying additional food to top predator. This study is aimed to introduce a new non-chemical chaos control mechanism in a predator–prey system with the applications in fishery management and biological conservation of prey predator species.     

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.     

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.     

Stability of a predator–prey model with refuge effect

We consider a predator–prey model, where some prey are completely free from predation within a temporal or spacial refuge. The most common type of spacial refuge, that we investigate here, takes the form where a constant proportion of the prey population is protected. The model is a modification of the classical Nicholson–Bailey host-parasitoid model. In this paper, we study the effect of the presence of refuge on the stability and bifurcation of the system. Moreover, we provide a detailed analysis of the Neimark–Sacker bifurcation of the model.     

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.     

Dynamics analysis of a pest management prey–predator model by means of interval state monitoring and control

In this work, a new pest management strategy by means of interval state monitoring is introduced into a prey–predator model, i.e. when the pest density exceeds the slightly harmful level but is below the damage level, the biological control is adopted in case of the predator density below a maintainable level, once the pest density exceeds the damage level, the chemical control is adopted. In order to determine the frequency of the chemical control and yield of releases of the predator, analysis on the existence of order-1 or order-2 periodic orbit is carried out by the construction of Poincaré map. The results could make the pest control strategy to be a periodic one without real-time monitoring the species. In addition, the stability and attractiveness of the periodic orbit are obtained by geometry approach, which ensures a certain robustness of control, i.e., even though the species densities are detected inaccurately or with a deviation, the system will be eventually stable at the periodic orbit under the control action. Furthermore, to obtain the optimum chemical control strength and yield releases of the predator, an optimization problem is constructed. The analytical results presented in the work are validated by numerical simulations for a specific model.     

Consequences of refuge and diffusion in a spatiotemporal predator–prey model

In this investigation, we offer and examine a predator–prey interacting model with prey refuge in proportion to both the species and Beddington–DeAngelis functional response. We first prove the well-posedness of the temporal and spatiotemporal models which are restricted in a positive invariant region. Then for the temporal model, we analyse its temporal dynamics including uniform boundedness, permanence, stability of all feasible non-negative equilibria and show that refugia can induce periodic oscillation via Hopf bifurcation around the unique positive equilibrium; for the spatiotemporal model, we not only investigate its permanence, stability of non-negative constant steady states and Turing instability but also study the existence and non-existence of non-constant positive steady states by Leray–Schauder degree theory. The key observation is that the coefficient of refuge cooperates a significant part in modifying the dynamics of the current system and mediates the population permanence, stability of coexisting equilibrium and even the Turing instability parameter space. Finally, general numerical simulation consequences are given to illustrate the validity of the theoretical results. Through numerical simulations, one observes that the model dynamics shows prey refugia and self-diffusion control spatiotemporal pattern growth to spots, stripe–spot mixtures and stripes reproduction. The outcomes assign that the dynamics of the model with prey refuge is not simple, but rich and complex. Additionally, numerical simulations show that the other model parameters have an important effect on species’ spatially inhomogeneous distribution, which results in the formation of spots pattern, mixture of spots and stripes pattern, mixture of spots, stripes and rings pattern and anti-spot pattern. This may improve the model dynamics of the prey refuge on the reaction–diffusion predator–prey system.     

Effects of prey refuge on a ratio-dependent predator–prey model with stage-structure of prey population

In this paper, a stage-structured predator–prey model is proposed and analyzed to study how the type of refuges used by prey population influences the dynamic behavior of the model. Two types of refuges: those that protect a fixed number of prey and those that protect a constant proportion of prey are considered. Mathematical analyses with regard to positivity, boundedness, equilibria and their stabilities, and bifurcation are carried out. Persistence condition which brings out the useful relationship between prey refuge parameter and maturation time delay is established. Comparing the conclusions obtained from analyzing properties of two types of refuges using by prey, we observe that value of maturation time at which the prey population and hence predator population go extinct is greater in case of refuges which protect a constant proportion of prey.     

A delayed prey–predator model with Crowley–Martin‐type functional response including prey refuge

In this paper, we have studied a prey–predator model living in a habitat that divided into two regions: an unreserved region and a reserved (refuge) region. The migration between these two regions is allowed. The interaction between unreserved prey and predator is Crowley–Martin‐type functional response. The local and global stability of the system is discussed. Further, the system is extended to incorporate the effect of time delay. Then the dynamical behavior of the system is analyzed, taking delay as a bifurcation parameter. The direction of Hopf bifurcation and the stability of the bifurcated periodic solution are determined with the help of normal form theory and centre manifold theorem. We have also discussed the influence of prey refuge on densities of prey and predator species. The analytical results are supplemented with numerical simulations. Copyright © 2017 John Wiley & Sons, Ltd.     

On Nonautonomous Prey predator Patchy System

§１．ＩｎｔｒｏｄｕｃｔｉｏｎＩｎｔｅｒｅｓｔｈａｓｂｅｅｎｇｒｏｗｉｎｇｉｎｔｈｅｓｔｕｄｙｏｆｍａｔｈｅｍａｔｉｃａｌｍｏｄｅｌｓｏｆｐｏｐｕｌａｔｉｏｎｓｄｉｓｐｅｒｓ－ｉｎｇａｍｏｎｇｐａｔｃｈｅｓｉｎａｈｅｔｅｒｏｇｅｎｅｏｕｓｅｎｖｉｒｏ...     

Bifurcations in a Leslie–Gower model with constant and proportional prey refuge at high and low density

Assuming that the prey refuge is proportional to the prey density if its population size is below a critical threshold, or constant if its size is above the threshold, this paper proposes, and qualitatively analyzes, a Leslie–Gower predator–prey model assuming alternative feeding and harvesting in predators, and a Holling II function as the predator functional response. From the results of the mathematical analysis to the predator–prey models with proportional or constant prey refuge, the proposed model retains the same bifurcation cases obtained for each model analyzed. However, appropriate alterations of the parameters representing the critical threshold of prey population size and harvest in predators allows the formation of at least one limit cycle, stable or unstable, that lives in both vector fields of the proposed model.     

Influence of predator mutual interference and prey refuge on Lotka–Volterra predator–prey dynamics

A Lotka–Volterra predator–prey model incorporating a constant number of prey using refuges and mutual interference for predator species is presented. By applying the divergency criterion and theories on exceptional directions and normal sectors, we show that the interior equilibrium is always globally asymptotically stable and two boundary equilibria are both saddle points. Our results indicate that prey refuge has no influence on the coexistence of predator and prey species of the considered model under the effects of mutual interference for predator species, which differently from the conclusion without predator mutual interference, thus improving some known ones. Numerical simulations are performed to illustrate the validity of our 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.     

Revisited Hastings and Powell model with omnivory and predator switching

The effect of omnivory in predator–prey system is debatable regarding its stabilizing or destabilizing characteristics. Earlier theoretical studies predict that omnivory is stabilizing or destabilizing depending on the condition of the system. The effect of omnivory in the food chain system is not yet properly understood. In the present paper, we study the effect of omnivory in a tri-trophic food chain system on the famous Hastings and Powell model. Omnivory enhances the chance of predator switching between prey and middle predator. The novelty of this paper is to study the effect of predator switching of the top predator which is omnivorous in nature. Our results suggest that in the absence of switching, an increase of omnivory stabilizes the system from chaotic dynamics, however, if we further increase the strength of omnivory, the system becomes unstable and middle predator goes to extinction. It is also observed that the predator switching enhance the stability and persistence of all populations.     

The dynamics of an age structured predator–prey model with disturbing pulse and time delays

Many of the existing predator–prey models on stage structured populations are some ordinary differential equations (ODE) or models without a disturbing effect of human behavior. In reality, death of the juvenile during its immature stage and catching or poisoning for the prey or predator occur continuously. From this basic standpoint, we formulate a general and robust prey-dependent consumption predator–prey model with periodic harvesting (catching or poisoning) for the prey and stage structure for the predator with constant maturation time delay (through-stage time delay) and perform a systematic mathematical and ecological study. We show that the conditions for global attractivity of the ‘predator-extinction’ (‘predator-eradication’) periodic solution and permanence of the population of the model depend on time delay, so, we call it “profitless”. We also show that constant maturation time delay and impulsive catching or poisoning for the prey can bring great effects on the dynamics of system by numerical analysis. In this paper, the main feature is that we introduce time delay and pulse into the predator–prey (natural enemy–pest) model with age structure, exhibit a new modeling method which is applied to investigate impulsive delay differential equations, and give some reasonable suggestions for pest management.     

Extinction and permanence of the predator-prey system with general functional response and impulsive control

Traditional approach for modelling the evolution of populations in the predator-prey ecosystem has commonly been undertaken using specific impulsive response function, and this kind of modelling is applicable only for a specific ecosystem under certain environmental situations only. This paper attempts to fill the gap by modelling the predator-prey ecosystem using a ‘generalized’ impulsive response function for the first time. Different from previous research, the present work develops the modelling for an integrated pest management (IPM) especially when the stocking of predator (natural enemy) and the harvesting of prey (pest) occur impulsively and at different instances of time. The paper firstly establishes the sufficient conditions for the local and the global stabilities of prey eradication periodic solution by applying the Floquet theorem of the Impulsive different equation and small amplitude perturbation under a ‘generalized’ impulsive response function. Subsequently the sufficient condition for the permanence of the system is given through the comparison techniques. The corollaries of the theorems that are established by using the ‘general impulsive response function’ under the locally asymptotically stable condition are found to be in excellent agreement with those reported previously. Theoretical results that are obtained in this work is then validated by using a typical impulsive response function (Holling type-II) as an example, and the outcome is shown to be consistent with the previously reported results. Finally, the implication of the developed theories for practical pest management is illustrated through numerical simulation. It is shown that the elimination of either the preys or the pest can be effectively deployed by making use of the theoretical model established in this work. The developed model is capable to predict the population evolutions of the predator-prey ecosystem to accommodate requirements such as: the combinations of the biological control, chemical control, any functional response function, the moderate impulsive period, the harvest rate for the prey and predator parameter and the incremental stocking of the predator parameter.     

A simple Gause‐type predator–prey model considering social predation

The consumer–resource relationships are among the most fundamental of all ecological relationships and have been the focus of ecology since its beginnings. Usually are described by nonlinear differential equation systems, putting the emphasis in the effect of antipredator behavior (APB) by the prey; nevertheless, a minor quantity of articles has considered the social behavior of predators. In this work, two predator–prey models derived from the Volterra model are analyzed, in which the equation of predators is modified considering cooperation or collaboration among predators. It is well known that competition among predators produces a stabilizing effect on system describing the model, since there exists a wide set in the parameter space where the system has a unique equilibrium point in the phase plane, which is globally asymptotically stable. Meanwhile, the cooperation can originate more complex and unusual dynamics. As we will show, it is possible to prove that for certain subset of parameter values the predator population sizes tend to infinite when the prey population goes to extinct. This apparently contradicts the idea of a realistic model, when it is implicitly assumed that the predators are specialist, ie, the prey is its unique source of food. However, this could be a desirable effect when the prey constitutes a plague. To reinforce the analytical result, numerical simulations are presented.     

The effect of prey refuge and time delay on a diffusive predator‐prey system with hyperbolic mortality

In this article, we investigate the effect of prey refuge and time delay on a diffusive predator‐prey system with Holling II functional response and hyperbolic mortality subject to Neumann boundary condition. More precisely, we study Turing instability of positive equilibrium by using refuge as parameter, instability and Hopf bifurcation induced by time delay. In addition, by the theory of normal form and center manifold, we derive conditions for determining the bifurcation direction and the stability of the bifurcating periodic solution. © 2016 Wiley Periodicals, Inc. Complexity 21: 446–459, 2016     

Global qualitative analysis of a non-smooth Gause predator–prey model with a refuge

The present paper aims to provide a detailed qualitative analysis of a non-smooth Gause predator–prey model. In this model, the saturating functional response function with a discontinuity at a critical prey density was employed to show the effects of a prey refuge on the population dynamic behavior. Analysis of this model revealed rich dynamics including locally (or globally) stable canard cycles, a locally (globally) stable pseudo-equilibrium, unbounded trajectories in which both populations go to infinity or the prey goes to infinity and the predator dies out eventually. The main purpose of the present work is to carry out a completely qualitative analysis for this model. In particular, two sets of sufficient conditions drive both populations to approach infinity and the sufficient and necessary conditions for all of the other main results are presented.