Explosive tritrophic food chain models with interference: A comparative study
Introduction
Ecological models are increasingly used to emphasize the linkage and relationship between social and physical environments. They demonstrate our responses to such factors over entire life cycles. They are easy to understand and interpret, they are often based on causality and they can include time delay, age/stage structure, functional response, etc. The effect of these factors has illustrated complicated and rich dynamics such as stability, bifurcation, periodic solutions, persistence, etc [1]. Research results herein are helpful to predict the developing/decaying tendency of specific populations to determine the key factors behind extinction and to seek optimum strategies of preventing and controlling extinction of various species [2].
In population dynamics, a functional response refers to the change in the density of prey consumed per unit time by per predator, as the prey density changes [3]. Based on different backgrounds and experimental data, various forms of functional responses have been developed mimicing various processes of energy transfer in predator-prey systems [4]. Categories which specify the functional response: prey dependent (function of only prey density, viz. Holling type I - IV [5]) and prey-predator dependent (function of both prey-predator density, viz. Beddington-DeAngelis [6], Crowley-Martin [7], Hassel-Varley [8] etc) those are describing the predator’s interference on their hunting procedure. The inherent feature of these functions means that the more prey in the environment the better off the predator, which is true in many predator-prey interactions [9]. Also, non-monotonic response occurs at the microbial level when the nutrient concentration reaches a high level an inhibitory effect on the specific growth rate may occur, Holling type IV functional response are used to model such inhibitory effect. Some researchers suggest that a more suitable general predator-prey model should have a prey-predator dependent functional response [10], [11], [12], [13] and all the prey-predator dependent functional responses can provide better description of predator feeding over a range of predator-prey abundances present [14]. Feeding procedure of higher trophic level species is more complicated than its below trophic level species [15], [16]. In general it is shown in the nature that if lower trophic level species is specialist type, then there is high chance to be generalist type of the higher trophic level species [15], [16]. Sexual reproduction is very common process for the growth of that kind of higher trophic level generalist species [17], [18], [19], [20], [21]. Also gestation period is very important factor in population dynamics and assuming that reproduction of predator after consuming prey is not instantaneous but mediated by some time lag required for gestation [22], [23]. It has been shown that many models in the class of where the top predator is modeled via the modified Leslie-Gower scheme for sexually reproductive generalist predator [22], [24].
We first provide some ecological background of invasive species. We are motivated primarily to control non-native species, which is a cental problem in spatial ecology. Data on species such as the invasive Burmese python (Python bivittatus) in the Florida everglades, show an exponential increase in python population, which have resulted in local prey population reducing severely [25]. An invasive species is formally defined as any species capable of propagating itself in a nonnative environment and thus establish a self-sustained population. In the United States alone damages caused by invasive species to agriculture, forests, fisheries and businesses, have been estimated to be $120 billion a year [26]. Another potential issue is that the environment may turn favorable for a certain species while becoming unfavorable for its competitors or natural enemies. This can result in the favored species to outbreak [27]. For example, in the European Alps certain seasonal environmental conditions enable the population of the larch budmoth to become large enough to defoliate entire forests [28]. then for sufficiently large initial data of top predator density, a potential to blow-up/explode in finite time can occur [24]. Biological control is an adopted strategy to limit harmful populations [29]. The objective of a biological control is to establish a management strategy that best controls and decreases the harmful population to healthy levels as opposed to high and risky levels. Naturally, how does one define high level, and further, how well does the biological control actually work, at various high levels? We have recently started investigating this question via the mathematical property of finite time blow-up [30], [31], [32].
This ideas have also been applied to model invasive populations that seems to be “exploding”, under a variety of ecological scenarios. In the current manuscript we ask:
- •
What is the foraging pressure of intermediate specialist predator on prey if it is modeled by Lotka–Volterra scheme and generalist sexually reproductive top predator on intermediate predator if the top predator is modeled according to the modified Leslie–Gower scheme?
- •
What are the dynamical outcomes when the predation process of intermediate and top predator are considered by different combinations of prey dependent and prey-predator dependent functional responses? With respect to the different combinations of functional responses different models will come for analysis.
- •
We ask what would the effects be, if there is predator interference as well.
- •
What will be periodic effect of gestation delay of intermediate predator on the steady state of the system?
- •
In particular, we ask what would the finite time blow up be, if the top generalist predator grows by sexual reproduction as well.
- •
How would the blow-up/blow up-prevention be effected in the case of various functional responses/feeding rates of the preys.
In this study, the effect of top predator interference on the dynamics of a general tri-trophic food chain model with five different combinations of functional responses is investigated. The manuscript is organized as follows: Models are formulated in Section 2 and stability properties of these models are studied in Section 3. Also, existence of Hopf-bifurcation for coexisting equilibrium point of all models are investigated in this section. In Section 4, finite time blow up behavior is investigated and numerical simulation has been carried out in Section 5. A comparative study between different models and a brief conclusion are in Sections 6 and 7 respectively.
Section snippets
Model formulation
First, we represent prey, intermediate specialist predator and sexually reproductive generalist top predator as X, Y and Z respectively in a general model, the system is:where all parameters are positive constants. a1 is the intrinsic growth rate and K is the environmental carrying capacity of X. a2 and θ are the death rate and food conversion efficiency of Y. c denotes the growth rate of Z. D represents the residual loss in Z
Stability properties of the systems
Feasibility or biologically positivity studies aim to objectively and rationally uncover the strengths and weaknesses of an existing proposed model in a given environment. Ecological stability can refer to types of stability in a continuum ranging from regeneration via resilience (returning quickly to a previous state), to constancy to persistence. The precise definition depends on the ecosystem in question, the variable or variables of interest and the overall context. In the context of
Motivation, control mechanisms and “Ecological” damping
In this section, we present results on finite time blow-up of the five delayed systems (2)–(6) and their non-delayed counterparts. From an ecological viewpoint, it is critical to derive the sharp conditions of population explosion/blow-up for any ecological model system. Parshad et al. [34] established that the solution to the three species model studied by Aziz-Alaoui [35], can blow-up in finite time even under the restriction derived in [35]. It is also useful then to study ecological
Numerical simulation
To demonstrate the analytic results of models (2)–(6) and their non-delayed counter part, we choose five different parameter sets for each model which are given in Table 9. For each model, we show the stability of the coexistence equilibrium point (E*) by phase-space diagram (Fig. 1.n(i), n=a,b,c,d,e) and bifurcation diagram (Fig. 1.n(ii), n=a,b,c,d,e). In Fig. 1.n(i) (n=a,b,c,d,e), if we choose (black trajectory) and τ < τ0 (green trajectory), then the coexistence equilibrium point (E*) is
Comparative study between different models
From the model formulation it is clear that the common issues of all the models are prey (X) grows logistically and top predator (Z) grows by sexual reproduction (for so mating of male and female happen and represented by term in third equation), also Z is generalist type and it is represented by modified Leslie-Grower scheme (Third equation of each model). Only differences between the systems (2)-(6) come out with respect to the different combinations of functional responses of
Conclusions and discussions
In this paper, a general framework has been structured to study the dynamical behavior of hybrid tritrophic food chain model consisting specialist middle predator and sexually reproductive generalist top predator. The different combinations of prey dependent and prey-predator dependent functional responses are used to formulate five differnt model systems. We discussed the stability properties of each equilibrium points for non-delayed couter-part of each model system. Existence of
Acknowledgments
The authors are grateful to the anonymous referees and Editor in Chief of this journal for their careful reading, valuable comments and helpful suggestions, which have helped them to improve the presentation of this work significantly.
References (59)
- et al.
Stability and Hopf bifurcation analysis of a ratio-dependent predatorprey model with two time delays and Holling type III functional response
Appl. Math. Comput.
(2015) - et al.
Diseased prey predator model with general Holling type interactions
Appl. Math. Comput.
(2014) - et al.
A food chain model with impulsive perturbations and Holling IV functional response
Chaos, Solitons Fractals
(2005) - et al.
Analysis of a diffusive Leslie–Gower predator-prey model with nonmonotonic functional response
Chaos Solitons Fractals
(2014) - et al.
Complex dynamics of a discrete-time predator-prey system with Holling IV functional response
Chaos Solitons Fractals
(2016) - et al.
Persistence and global stability of Bazykin predator-prey model with Beddington–DeAngelis response function
Commun. Nonlinear Sci. Numer. Simul.
(2014) - et al.
Species extinction problem: genetic vsecological factors
Appl. Math. Modell.
(2001) - et al.
Do prey handling predators really matter: Subtle effects of a Crowley–Martin functional response
Chaos Solitons Fractals
(2017) - et al.
Dynamics of generalist predator in a stochastic environment: effect of delayed growth and prey refuge
Appl. Math. Comput.
(2015) The theory and classification of outbreaks
Insect Outbreaks
(1987)
Finite time blowup in a realistic food chain model
ISRN Biomath.
A remark on “study of a Leslie-gower-type tritrophic population model”
Chaos Solitons Fractals
Uniform blow-up rates and asymptotic estimates of solutions for diffusion systems with weighted localized sources
J. Appl. Math. Comput.
Global boundedness and blow-up for a parabolic system with positive Dirichlet boundary value
J. Appl. Math. Comput.
Blow-up in a three species cooperating model
Appl. Math. Lett.
On diffusion induced blowups in a mutualistic model
Nonlinear Anal.
Dynamics of a three species competition model
Discr. Contin. Dyn. Syst. A
A users guide to PDE models for chemotaxis
J. Math. Biol.
Predator community structure and trophic linkage strength to a focal prey
Mol. Ecol.
Intraguild predation among biological-control agents: theory and practice
Biol. Control
Nonfragile asynchronous control for uncertain chaotic lurie network systems with bernoulli stochastic process
Int. J. Robust Nonlinear Control
The functional response of predators to prey density and its role in mimicry and population dynamics
Mem. Entomol. Soc. Can.
Mutual interference between parasites or predators and its effect on searching efficiency
J. Anim. Ecol.
Functional responses and interference within and between year classes of a dragonfly population
J. North Am. Benthol. Soc.
New inductive population model for insect parasites and its bearing on biological control
Nature
Variation in plankton densities among lakes: a case for ratio-dependent predation models
Am. Nat.
Physiological basis of ratio-dependent predator–prey theory: the metabolic pool model as a paradigm
Ecology
Dynamical analysis of a prey–predator model with Beddington-DeAngelis type function response incorporating a prey refuge
Nonlinear Dyn.
Functional responses with predator interference: viable alternatives to the Holling type II model
Ecology
Cited by (3)
Controlling chaos and pattern formation study in a tritrophic food chain model with cannibalistic intermediate predator
2022, European Physical Journal PlusModeling the plankton–fish dynamics with top predator interference and multiple gestation delays
2020, Nonlinear Dynamics