arxiv: v1 [q-bio.pe] 10 Oct 2015

Transcription

1 PREY CANNIBALISM ALTERS THE DYNAMICS OF HOLLING-TANNER TYPE PREDATOR-PREY MODELS arxiv:5.2989v [q-bio.pe] Oct 25 Abstract. Cannibalism, which is the act of killing and at least partial consmption of conspecifics, is biqitos in natre. Mathematical models have considered cannibalism in the predator primarily, and show that predator cannibalism in two species ODE models provides a strong stabilizing effect. There is strong ecological evidence that cannibalism exists among prey as well, yet this phenomenon has been mch less investigated. In the crrent manscript, we investigate both the ODE and spatially explicit forms of a Holling-Tanner model, with ratio dependent fnctional response. We show that cannibalism in the predator provides a stabilizing inflence as expected. However, when cannibalism in the prey is considered, we show that it cannot stabilise the nstable interior eqilibrim in the ODE case, bt can destabilise the stable interior eqilibrim. In the spatially explicit case, we show that in certain parameter regime, prey cannibalism can lead to pattern forming Tring dynamics, which is an impossibility withot it. Lastly we consider a stochastic prey cannibalism rate, and find that it can alter both spatial patterns, as well as limit cycle dynamics. Aladeen Basheer, Emmanel Qansah, Schman Bhowmick and Rana D. Parshad Department of Mathematics, Clarkson University, Potsdam, New York 3699, USA.. Introdction Cannibalism is the act of killing and at least partial consmption of conspecifics. This has been docmented in hman poplations all arond the world. It is/has been practiced in varios indigenos Soth American and New Ginean tribes as a social norm [6], and also reported in Africa [23]. It is also practiced among Aghoris or witch doctors in Northern India, in the hope of gaining immortality [25]. However, if one steps away from these seemingly macabre and gholish hman practices, and observes real predator-prey systems, then cannibalism is actally commonplace [8]. It is seen to actively occr in more than 3 species in natre [24]. Apart from being a fascinating socio-anthropological/ecological sbject in its own right, cannibalism can lead to varios conter intitive effects as concerns poplation dynamics. These inclde the life boat mechanism, where the act of cannibalism cases persistence, in a poplation doomed to go extinct [2]. It can also indces stability, in an otherwise cycling poplation [7]. A detailed srvey of the crrent and past mathematical literatre on cannibalism, shows that varios types of predator-prey models, with the inclsion of cannibalism, 99 Mathematics Sbject Classification. Primary: 35B36, 37C75, 6H35; Secondary: 92D25, 92D4. Key words and phrases. Holling-Tanner model, prey cannibalism, stability, Tring instability, white noise.

2 2 BASHEER, QUANSAH, BHOWMICK, PARSHAD have been investigated. These inclde two species ODE models, three species ODE models incorporating stage strctre, two species PDE models, discrete models and integro-differential eqation models, as well as recent models of cannibalism that inclde diseased predators [7, 4, 3, 5, 2]. The two species nstrctred ODE models that have been investigated, are of the general type d dt = fv, K 2 dv = dv + ɛfv cc, v, dt where is the prey, and v the predator. Here cc, v models cannibalism in the predator, with c being the rate of cannibalism. In sch models, cannibalism is seen to have a strong stabilising inflence. Kohlmeir [7] was the first to show this, by considering cannibalism in the predator, in the classic Rosenzweig-McArthr model. That is when f is of Holling type II. They showed that predator cannibalism can srprisingly stabilise the nstable interior eqilibrim of the model. We note that cannibalism, has been primarily modeled as an act in the predator species only [, 34, 7]. Note, there is strong evidence in ecology that cannibalism often occrs in the prey species [27, 26]. Recent experiments with prey cannibalism, show that they can have a significant impact on the poplation dynamics of predator-prey systems [28]. However, there is not a single model in the mathematical literatre, that considers the effect of prey cannibalism in two species nstrctred ODE or nstrctred PDE models. To the best of or knowledge there are very few works overall, that have considered prey canibalism, in some form or the other. Solis et al. considered the effect of prey cannibalism in a mlti component strctred model, via integer differential eqations [3]. Chow and Jang investigated a strctred discrete model with prey cannibalism finding that it may either stabilize or destabilize, depending on cannibalism rates and varios other parameters [7]. Also, Zhang et al. considered a three species strctred ODE model, with strctre and cannibalism in the prey. They showed that cannibalism strengthens the conditions of both global asymptotical stability, and persistence of the interior eqilibrim point [37]. Many models of cannibalism with size and stage strctre effects have been investigated. The approach in the literatre has been to strctre the predator class into adlts and jveniles and allow predation of the jvenile predator by the adlt predator. This essentially yields a three species ODE model. The reader is referred to [5, 4, 37], and again therein cannibalism is seen to provide a stabilising inflence, depending on parameter regime. However, depending on model strctre and parameter regimes, cannibalism is also seen to destabilize [8] in some cases. Varios other classes of models sch as integro-differential eqation models, as well as discrete models have also been investigated [7, 3, 2]. Since species disperse in space in search of food, shelter, mates and to avoid predators, spatially dispersing poplations are often modeled via partial differential eqations PDE/spatially explicit models of interacting species [36, 2, 22, 29]. Of particlar interest to many is the phenomenon of Tring instability, first introdced by Alan Tring [35], which shows that diffsion in a system can lead to pattern forming instabilities. In the context of cannibalism, there are very few works in the

3 PREY CANNIBALISM 3 literatre that deal with spatially explicit models and/or Tring instability. Sn [34] was the first to show that cannibalism in the classical diffsive Rosenzweig- McArthr model, can bring abot Tring patterns. This is an impossibility in the system withot cannibalism. Fasani et al. [] considered the same model, bt with a general form for the cannibalism term of which the model of Sn et al. is a special case, and also fond that cannibalism can bring abot Tring patterns. The other direction does not seemed to have been stdied. That is when and nder what conditions, can cannibalism remove Tring instability? In the crrent manscript we arge that these above approaches, do not captre the fll breadth of cannibalism as it exists in natre. We cite the following reasons: Many mathematical models of predator-prey systems, do not have a symmetric fnctional responses between predator and prey. For example, the models of Leslie- Gower type [8], and Holling-Tanner type [4]. These models are biqitos in the mathematical literatre see [2, 4] and the references therein, however, the effect of cannibalism, in particlar the effect of prey cannibalism in sch models, has not been considered. Also, In all of the investigations in the literatre, only a deterministic rate of cannibalism, c, has been considered. However, predator and prey species are often sbject to ncertainties in the environment, modeled best via environmental noise. For these reasons it is possible that the varios rates, birth rate, death rate, predation rate etc., one considers for modeling prposes, may not be known precisely, and may always involve ncertainty. Yet there are no stdies where a model with a noisy cannibalism rate has been considered, that is when c c + noise. In the crrent manscript we consider a Holling-Tanner predator-prey model, with a ratio dependent fnctional response [], that has been well stdied in the literatre [2]. Next we introdce into this model a cannibalism occrring only in the predator species, b cannibalism occring only in the prey species. We investigate in detail: ODE versions modeling a and b. 2 The PDE, stochastic ODE and stochastic PDE versions modeling b. Or contribtions therein are the following: We show that predator cannibalism, in certain parameter regime, can indeed stabilize the nstable interior eqilibrim, in the ODE case, as expected from the known reslts on other models, with cannibalism in the predator. 2 We show that prey cannibalism cannot stabilize the nstable interior eqilibrim, in the ODE case. This is contrary to what is known abot models with cannibalism in the predator. 3 We show that prey cannibalism can destabilize the spatially homogenos eqilibrim in the PDE case. In particlar, we show that in certain parameter regime sch that there is no Tring instability in the Holling-Tanner model, introdcing prey cannibalism can actally bring abot Tring instability. Alternatively, we show that in certain parameter regime, introdcing predator cannibalism cannot bring abot Tring instability. 4 We show that in certain parameter regime, if Tring instability does not exist in

4 4 BASHEER, QUANSAH, BHOWMICK, PARSHAD the model with prey cannibalism, then removing prey cannibalism cannot indce Tring instability. 5 The ODE version of the Holling-Tanner model we consider is known to possess a niqe limit cycle. We investigate the effects of a stochastic prey cannibalism rate c, on the limit cycle dynamics, finding that classical white noise cannot indce stability. 6 We also consider the effects of a stochastic prey cannibalism rate c, on the spatial dynamics, finding that space-time white noise, can significantly alter the spatial dynamics of the model. 2. The model formlation 2.. The Holling-Tanner model. We first present the following Holling-Tanner model, with ratio dependent fnctional response in prey, [2]. d dt = v 3 + αv dv dt = δv v nc model We will refer to the above model as the nc-model no cannibalism henceforth, and recap certain dynamical aspects of the model. Note, 3 is a nondimensionalised version of the model presented in [2]. The reader is referred to [2], for a complete derivation on the parameters, their biological significance and the nondimensionalization process. Here and v are the poplations of the prey and predator species respectively. The dynamics of the prey species follows a ratio dependent fnctional response []. The predator is modeled according to the Leslie-Gower scheme [8]. This assmes that the carrying capacity of the predator, is not constant bt rather proportional to its food, the prey. Hence a reasonable choice for the carrying capacity K = K is K = + d 3. If d 3 >, the predator is thoght of as a generalist, that is it can change its food sorce in the absence of its favorite prey. If d 3 =, the predator is thoght of as a specialist. In this case the reslting model is of Holling-Tanner type [4]. Here is the intrinsic growth rate prey handling time of the predator v, and α is essentially given by α = prey captre rate. The non-trivial steady state nc, v nc is given by: since we reqire nc = + α, v nc = + α nc >, = > + α Ths we obtain the following restriction on de to a feasibility reqirement, < α For stability we reqire that the T race J nc < which implies α + + α + 2 < δ

5 PREY CANNIBALISM 5 As the determinant Det J nc >, is always positive, as we have < α, which is always tre as this is a feasibility condition The Holling-Tanner model with predator cannibalism. We now consider the inclsion of predator cannibalism in 3. This yields the following model 4 d = v dt + αv dv dt = δv v γ + cv predator cannibalism model In deriving the above model we make the following assmptions The predator v is depredating on the prey species, as well as on its own species v. Hence the food sorce of v is γ + cv, where cv is the cannibalism term, with cannibalism rate c. Here we assme c+γ, with γ <, c <. Ths the food intake of the predator stays roghly the same, even if it trns into a cannibal. It starts eating less prey, so γ <, and now also eats some conspecifics, so < c <. There is a clear gain of energy to the cannibalistic predator, from the act of cannibalism. This gain reslts in an increase in reprodction in the predator. This in trn leads to a gain in predator poplation, modeled via adding a v term to the predator eqation. We therefore assme >. The non-trivial steady state c, v c is given by 5 c = γ γ, v c = c α γ c + c since c > we obtain γ > and de to the assmptions and feasibility conditions, γ <, α γ c + c > we have αγ γ + > c The Jacobian Matrix at c, v c : c + c + αγ 2 J c = γ c c 2 α γ c + 2 γδ 2 δ c Notice that J 2 and J 22 are negative since γ < and c >. Det J c = c + c + αγ 2 δ c + γδ 2 c 2 γ c α γ c + 2 > T c = c + c + αγ 2 + δ c γ c We can ths state the following lemma

6 6 BASHEER, QUANSAH, BHOWMICK, PARSHAD Lemma 2.. Consider a parameter set s.t. T nc >, that is the interior eqilibrim point nc, v nc of the no cannibalism model 3 is nstable. If γ < then there exists a cannibalism rate c s.t T c <. This implies that the nstable interior eqilibrim point nc, v nc of the no cannibalism model 3 can be stabilised via predator cannibalism. However If γ > then the nstable interior eqilibrim point nc, v nc of the no cannibalism model 3 cannot be stabilised via predator cannibalism. Proof. We take a simple geometric approach and visalise the nllclines below in fig.. In order to shift the nllcline to the right de to predator cannibalism, we will γ reqire that the slope of the predator nllcline is less than. That is c <. If γ = ɛ <, we can always choose c small enogh s.t ɛ < c. In this case the slope of the predator nllcline with predator cannibalism will be less than. If on the other hand γ γ > > c. <, this is clearly an impossibility as J nc > J c < Prey Nllcline Predator Nllcline Withot Cannibalsm Nllcline With Predator Cannibalsm v > > Figre. Here we see the nllclines of the model 3 vs the nllclines of model 4. We see that if γ < predator cannibalism, decreases the slope of the predator nllcline, shifting it to the right, and making it possible for J c < and hence for T c <, while T nc >. In this figre the parameters chosen are α =.89, = 3.5, c =.8, δ =., γ =.2. Remark. The cannibals feeding rate of prey γ and the ratio, are critical in terms of stability. The stability condition tells s that if the cannibal eats less prey

7 PREY CANNIBALISM 7 so γ is small in comparison to, then only a little bit of feeding of conspecifics is enogh to stabilize the dynamics The Holling-Tanner model with prey cannibalism. We now consider the inclsion of prey cannibalism in 3. This yields the following model d 6 dt = + c v + αv c 2 + d dv dt = δv v prey cannibalism model In deriving the above model we make the following assmptions The prey is depredating on its own species. Here the generic cannibalism term C, is added in the prey eqation, and is given by C = c + d. That is the fnctional response of the cannibalistic prey is of Holling type II, with cannibalism rate c. There is a clear gain of energy to the cannibalistic prey, from the act of cannibalism. This gain reslts in an increase in reprodction in the prey. This in trn leads to a gain in prey poplation, modeled via adding a c term to the prey eqation. We restrict c < c, as it takes depredation of a nmber of prey by the cannibal to prodce one new offspring. The non-trivial steady state c, v c is given by c = m + d + c c + m + d + c c 2 + 4d m + 4dc 2 where m = α + <. The feasibility conditions are divided into the following cases, if m < + c and m + d + c + c 2 if m < + c and m + d + c < + c The Jacobian matrix of system 6 at the point c, v c given by: J c = + α 2 c cd c c + d 2 2 δ The conditions for stability are: + α 2 δ The Trace J c < if Trace J c <, det J c > + α 2 c cd c c + d 2 < δ

8 8 BASHEER, QUANSAH, BHOWMICK, PARSHAD The above implies Note + α 2 c + d 2 c δ < c d c Det J c = δ c[d 2 + 2d c + cd + c 2 ] d 2 + 2d c + c 2 > This is always tre Effect of prey cannibalism on stability. We now ask the following qestion. Consider the nc-model in a parameter regime, sch that the interior eqilibrim nc, v nc is stable. Can prey cannibalism destabilize this eqilibrim? If so, nder what conditions? This is settled via the following lemma, Lemma 2.2. For any selection of parameters, there always exists a cannibalism rate c s.t the interior eqilibrim point nc, v nc of the nc-model is stable, while the interior eqilibrim point c, v c of the prey cannibalism model is nstable. Proof. Note that the condition nder which the interior eqilibrim point nc, v nc of the nc-model is stable, while the interior eqilibrim point c, v c of the prey cannibalism model is nstable, is T nc < < T c The above is eqivalent to showing that c s.t Note c + cd c c + d 2 < nc. c+d 2 2 c + cd c c + d 2 < c + c c + d 2 = c + c 2 < nc. Ths we want to show c + c 2 < nc. Now c = m + d + c c + m + d + c c 2 + 4d m + 4dc 2 It is easily comptable that d c dc = 2 + m + d + c c <, 2 m + d + c c 2 + 4d m + 4dc that is c is monotonically decreasing in c, and so decreasing c to a small enogh level, will eventally case c + c 2 < nc. This proves the lemma. We next ask the following qestion. Consider the nc-model in a parameter regime, sch that the interior eqilibrim nc, v nc is nstable. Can prey cannibalism stabilize this eqilibrim? If so, nder what conditions? In order to answer this qestion, we mst frther analyse the dynamical behavior of the prey cannibalism model. We state or reslt via the following lemma:

9 PREY CANNIBALISM 9 Lemma 2.3. If m + d + c < + c then then there does not exist a cannibalism rate c s.t c + cd c c + d 2 > nc. This implies that no amont of prey cannibalism can stabilize the nstable interior eqilibrim point nc, v nc of the prey cannibalism model. Proof. The explicit condition nder which the eqilibrim point nc, v nc of the nc-model is nstable and the eqilibrim point c, v c of the prey cannibalism model is stable is given by T c < < T nc This is eqivalent to c + cd c c + d 2 < nc. Now, c + cd c c + d 2 = c + c 2d c 2 c + d 2 < c + c c + d 2 2 c + d 2 = c + c 2. Ths if we can show c + c 2 < nc then this will imply that c + We let fc = c + c 2 nc, note f =. cd c c + d 2 < c + c 2 < nc. We next compte the derivative of f w.r.t. c. fc = m + d + c c + m + d + c c 2 + 4d m + 4dc + c 2 2 nc f c = 2 + m + d + c c + 2 m + d + c c d m + 4dc m + d + c c = 2 m + d + c c 2 + 4d m + 4dc < 7 since m + d + c < + c. This implies that f c <. Ths fc is decreasing, and since f =, it mst be that fc = c + c 2 nc <, for all c. c + Ths cd c c + d 2 < c + c 2 < nc. This proves the lemma. Remark 2. The above proof is nder the first feasibility condition m + d + c < + c, that is small cannibalism parameter c. Nmerically, we cold not find a cannibalism rate c sch that the eqilibrim point nc, v nc of the nc-model is nstable while the eqilibrim point c, v c of the prey cannibalism model is stable, nder the second feasibility condition m + d + c > + c, that is when the cannibalism parameter c can be very large. We tried a parameter sweep of all parameters in the interval [, ]. We ths beleive that prey cannibalism cannot

10 BASHEER, QUANSAH, BHOWMICK, PARSHAD stabilise, whatever the range of the cannibalism parameter c be, large or small. However, the large c case is nproven at this point. Nonetheless, we make the following conjectre Conjectre. There does not exist a cannibalism rate c, sch that the eqilibrim point nc, v nc of the nc-model is nstable while the eqilibrim point c, v c of the prey cannibalism model is stable. 3. The effect of Noisy cannibalism rate on limit cycle dynamics In this section we aim to investigate the effect of a noisy cannibalistic parameter on the limit cycle dynamics in the prey cannibalism model. We proceed by stating the following lemma, Lemma 3.. The soltions of system 6 are always positive and bonded, frthermore there exists T sch that < t < + c, < vt < + c, for all t > T. Proof. The positivity easily follows by the qasi monotonicity of the reaction terms. For bondedness note that for system 6, we have d dt + c A standard comparison argment shows that lim t t + c.it follows that there exists T > sch that t < + c, for t > T. From the second eqation of system 6, we see that, for t > T, A standard comparison argment shows that lim t vt + c. Now, it can be shown that nder certain parametric restrictions, there exists a limit cycle in the prey cannibalism model 6. We state this via the following lemma, Lemma 3.2. Consider the prey cannibalism model 6. Assme that the following conditions hold 8 α + >, + α 2 c + d 2 c δ > c, d c then the system possesses at least one limit cycle. Here c is the prey eqilibrim in model 6. Proof. From the local stability analysis, we have that if the above conditions hold, then the positive eqilibrim E c = c, v c of system 6 is an nstable node or focs. The existence of a limit cycle now follows from lemma 3. and Poincairé- Bendixson theorem [32]. Remark 3. Althogh the existence of a niqe limit cycle for the nc-model is well known, and the existence of a limit cycle with the addition of prey cannibalism is seen via lemma 3.2, we are nable to prove niqeness, in this case. We trn or attention to investigating the effect that a stochastic cannibalism rate has on the limit cycle dynamics.

11 V V PREY CANNIBALISM Biological species are always sbject to random environmental flctations, be it via climatic factors, or ecological. Also, many times the parameters we se for modeling sch processes many not be known with certainty. In this present paper, we only assme that the cannibalism rate c is affected by environmental flctations, which is define as mortality rate de to cannibalism. We introdce a stochastic factor in time ηt, into the cannibalism rate c, and so the parameter c is treated as a Gassian random variable with mean ĉ and intensity α. Here ηt = d dt W t is a classical white noise. 9 c ĉ + α ηt. Under these assmptions we consider the following SDE model, where forsimplicity, we redefine mean ĉ to be the parameter c. d = + c v + αv c 2 2 dt + α + d + d dw t dv = δv v dt We observe that as α is increased, the limit cycle dynamics becomes aperiodic, bt some cyclicity remains. The noise is nable to stabilise, that is drive the poplation cycle to a steady state. This is shown via the following plots, U U Figre 2. Here we see a limit cycle in nc-model in the left panel with α =.. We observe that the effect of a noisy prey cannibalism rate is to distort the limit cycle, bt some cyclicity remains with α =.. For the list of parameters chosen please refer to Table. Remark 4. Ths it seems that even with some stochasticity in the cannibalism parameter, prey cannibalism is nable to stabilise dynamics. It is well known that large enogh noise intensity will drive everything to extinction, bt this is an ninteresting case. This leads s to make the following conjectre Conjectre 2. There does not exist a stochastic cannibalism rate c, sch that the non trivial eqilibrim point nc, v nc of the nc-model is nstable, while the nontrivial eqilibrim point c, v c of the stochastic prey cannibalism model is stable.

12 2 BASHEER, QUANSAH, BHOWMICK, PARSHAD Time Time Figre 3. The left panel in the above is the time series corresponding to the phase plot in the left panel in fig. 2 with α =.. The right panel corresponds to a slightly greater noise intensity with α =.. For the list of parameters chosen please refer to Table Time Time Figre 4. These figres show the effect of frther increasing the noise intensity, on the limit cycle dynamics where the left panel shows α =.5 and right panel shows α = 2.. For the list of parameters chosen please refer to Table. 4. The effect of prey cannibalism on Tring instability Tring instability occrs when the positive interior eqilibrim point for an ODE system is stable in the absence of diffsion, bt becomes nstable de to the addition of diffsion, that is in the PDE model. This phenomena is also referred to as diffsion driven instability, [35]. To formlate the spatially explicit form of the earlier described ODE models, we assme that the predator and prey poplations move actively in space. Random movement of animals occrs becase of varios reqirements and necessities like, search for better food, better opportnity for social interactions sch as finding mates [22]. Food availability and living conditions demand that these animals migrate to other spatial locations. In the proposed model, we have inclded diffsion assming that the animal movements are niformly distribted in all directions. The models with diffsion are described as follows t = D xx + v +αv, } v t = D v v xx + δv v no cannibalism model

13 PREY CANNIBALISM 3 and the diffsive model with prey cannibalism 2 t = D xx + + c v +αv c2 +d, } v t = D v v xx + δv v prey cannibalism model The spatial domain for the above is Ω R. Ω is assmed bonded, and we prescribe Nemann bondary conditions x = v x =, for x Ω, and sitable initial conditions. The main reason for choosing sch bondary conditions is that we are interested in the self-organization of patterns, and zero-flx bondary conditions imply that no external inpt is imposed [22, 2]. In the nmerical simlations we take Ω = [, π]. Here x, t is the density of the prey species at any given time t, that is logistically controlled, and predated on by a predator species with density vx, t. Remark 5. For varios models, that fall into the above general category, it has been shown that Tring instability exists. However, this is only for certain parameters that lie in the so called Tring space. For parameters otside of that space Tring patterns are an impossibility. 4.. General criterion for no Tring instability in nc-model. Or goal in the crrent section is to investigate the conditions nder which there is no Tring instability in, bt there are Tring patterns in 2. Ths we aim to qantify the conditions nder which cannibalism in the prey species, can indce Tring patterns. We begin by linearizing model system abot its homogenos steady state nc, v nc, by introdcing both space and time-dependent flctations. In order to derive the precise conditions for Tring instability there is a well known procedre. We refer the reader to [36]. The calclations therein can be stated via the following theorem, Theorem 4. Tring Instability Condition. Consider nc, v nc which is a spatially homogenos eqilibrim point of the model [ system ]. For a given set of J J parameters, if at nc, v nc, the Jacobian J = 2, of the reaction terms, J 2 J 22 and the diffsion coefficients D, D v satisfy J + J 22 < J J 22 J 2 J 2 >, D v J + D J 22 >. D v J + D J D D v J J 22 J 2 J 2 >, then nc, v nc is said to be linearly stable in the absence of diffsion and linearly nstable in the presence of diffsion. And easy corollary derived from the above imposes an additional necessary condition for Tring instability Corollary Necessary condition for Tring Instability. If nc, v nc is a spatially homogenos eqilibrim point of the model system which is linearly stable in the absence of diffsion, and linearly nstable in the presence of diffsion, then

14 4 BASHEER, QUANSAH, BHOWMICK, PARSHAD [ ] J J either J <, J 22 > or J >, J 22 <, where J = 2, is the Jacobian J 2 J 22 matrix of the reaction terms. We begin by focsing on, investigated in detail by Banerjee [2]. The reslts therein reveal that tring instability in model occrs, if α <, for a range of parameters in the so called Tring space. Note that the jacobin matrix of the reaction terms of at the point nc, v nc is given by J = α + + α + 2 α δ δ Since the parameters >, δ >, we always have that J 22 <. Ths an application of corollary yields the following lemma: Lemma 4.2. The eqilibrim point nc, v nc cannot be driven nstable via diffsion if J <, or eqivalently if α + 2 < α + 2. It is a matter of simple algebra to convert the above. That is: 3 α + 2 < α + 2 α + α Remark 6. Let s investigate this ratio frther. Note, is the intrinsic growth prey handling time rate of the predator, and α = prey captre rate. These are nondimensionalised parameters, and de to feasibility constraints α <. Roghly speaking, 3 is saying that if predator growth rate prey handling time prey captre rate then Tring instability is an impossibility in the nc model. Now, given a predator growth rate, If prey handling time is close enogh to the prey captre rate, then 3 will hold, and Tring instability cannot occr in. We now pase and ask the following qestion: Can cannibalism in the prey species alter this scenario? Specifically, can prey cannibalism, bring abot Tring instability in 2, even if 3 holds? To answer this qestion, we investigate the Tring instability conditions for 2. The calclations therein yield that for the spatially homogenos eqilibrim soltion c, v c to 2 the Jacobian matrix at c, v c is given by: J = + α 2 c cd c c + d 2 + α 2 2 δ δ For cannibalism in the prey to alter Tring dynamics via lemma 4.2, at the very least we reqire that J >. Clearly there is no Tring, if J <. Now, if J > it does not garantee there is Tring, as this condition is only necessary, bt there is a chance for Tring instability to occr. Ths In order for there to be a

15 PREY CANNIBALISM 5 possibility of Tring instability in, while there is no Tring instability in 2, we wold reqire J c >, that is 4 + α 2 c + d 2 c > cd, whilst < c α + α. We can also take a geometric approach similar to Fasani [] to proceed. We investigate the nllclines of, shown in fig Nllcline Withot Prey Cannibalsm Predator Nllcline Nllcline With Prey Cannibalsm J nc < v --> J c > > Figre 5. Here we see the nllclines of the model vs the nllclines of model 2. We see prey cannibalism drives down the prey nllcline. Remark 7. Since the sitation in fig. 5 is clearly achievable geometrically, there is a chance that prey cannibalism cold bring abot Tring instability, even if < α+. We then perform a parameter sweep based on the Tring α conditions that we derive, in order to see if the appropriate parameters can be fond Tring instability de to prey cannibalism. For details on the Tring calclations we refer the reader to [36]. The calclations therein yield varios necessary and sfficient conditions for Tring instability to occr, if prey cannibalism is present. We state or reslt via the following theorem.

16 t t 6 BASHEER, QUANSAH, BHOWMICK, PARSHAD Theorem 4.3 Prey cannibalism indced Tring Instability. Consider c, v c which is a spatially homogenos eqilibrim point of the model system 2. If there[ exist parameters ] α,, c, d, c all strictly positive, s.t at c, v c, the Jacobian J J J = 2, and the diffsion coefficients D J 2 J, D v satisfy 22 + α 2 c + d 2 c > c > d c + α 2 c + d 2 c δ d c D v + α 2 c cd c c + d 2 D δ > D v + α 2 c cd c δd c + d 2 D δ > 2 D D 2 c+2δd c 2 +cδd c+δ c 3 v d 2 +2d c+ c, 2 α + α, then c, v c is linearly stable in the absence of diffsion and linearly nstable in the presence of diffsion, whilst nc, v nc which is a spatially homogenos eqilibrim point of the model system is linearly stable in the presence of diffsion, for the same parameters α,, d, and c = c =. That is, in the absence of prey cannibalism. Remark 8. Note for system 2 and for we have that J J 22 J 2 J 2 >. A parameter sweep sing the above conditions yields that sch a parametric space is definitely possible. For chosen parameters in this space, we see Tring patterns as in fig v x x.22 Figre 6. Here the Tring patterns in predator and prey, are shown for model 2, that is when there is prey cannibalism. Spatial patterns are formed. In these simlations Ω = [, π], and the spatially homogenos steady state is given a small pertrbation by a fnction of the form ɛ cos 2 x, ɛ <<. For the parameters sed please see Table. Now we move in the opposite direction. Spposing that the conditions of theorem 4.3 do not hold, and we do not have Tring instability in the nc model 2. Can removing the prey cannibalism bring abot Tring instability? We next provide a

17 PREY CANNIBALISM x x a b Figre 7. Here the model 2 that is when there is prey cannibalism, is contrasted with model. We see in Fig a model 2 the density of the prey at time t=5. There is a spatial pattern. In Fig b model there is no pattern. For the parameters sed please see Table. sfficient condition nder which this is not possible. This is done via the following lemma Lemma 4.4. Consider c, v c which is a spatially homogenos eqilibrim point of the model system 2 which is linearly stable both in the absence of diffsion and in the presence of diffsion. If J c < or + α 2 c c + d 2 c < cd, and if m+d+c +c, and +α < 2+α. Then nc, v nc which is a spatially homogenos eqilibrim point of the model system, for the same parameter set, bt when c = c =, is also linearly stable in the presence of diffsion Proof. The proof follows via observing that if 5 + α 2 c c + d 2 c > cd => + α 2 < c + cdc c + d 2 < nc = + α < 2 + α This implies that α + 2 < + α 2, which in trn implies that <, and so Tring instability is an impossibility in, from lemma α+ α 4.2. Note c + cdc < c +d 2 nc follows from the assmption m + d + c + c, and lemma 2.3. Remark 9. Essentially removing prey cannibalism, shifts the prey nll cline p, still leaving J nc <, at the new eqilibrim nc, v nc. This is easily seen in fig. 5.

18 8 BASHEER, QUANSAH, BHOWMICK, PARSHAD We lastly state a reslt concerning Tring pattern formation and predator cannibalism Lemma 4.5. The eqilibrim point nc, v nc cannot be driven nstable via diffsion, in the presence of predator cannibalism, if and γ <. However, if γ >, then nc, v nc can be driven nstable via diffsion, in the presence of predator cannibalism. Proof. The proof is an immediate conseqence of the nll cline analysis. If J nc = α+ <, Tring instability is an impossibility. However, predator α cannibalism will only psh the nll cline frther to the right if γ <, and so J c <, ths Tring patterns are an impossibility, in this case. However, if γ >, then predator cannibalism will psh the predator nll cline to the left, and it is possible that J c >, so that Tring patterns cold occr. 5. Dynamics of the PDE model with noise In mathematical modeling, deterministic systems of PDEs have been relied pon to captre the relevant dynamics of some biological systems. However, these systems are freqently sbject to noisy environments, inpts and signalling. It was reported that songbirds experience decreased predation rate in the environments with noise [], and hermit crabs were distracted by boat motor noise and hence less vigilant against approaching predators [6, 3]. Ths stochastic pertrbations are important when considering the ability of sch models to reprodce reslts consistently. In particlar, owing to the Tring mechanism being very sensitive to pertrbations, stochasticity is able to affect evolving patterns in ways not seen in deterministic simlations [2]. Previos reslt in stochastic predator-prey models focsed on internal noise[9, 5, 33], where it is assmed that noise affect the model system itself. This allows to accont for sch flctations even in the framework of a continm model [5]. However, predator-prey models with Holling type II external noise had been generally overlooked despite its potential ecological reality and theoretical interest. Recent theoretical work indicates that cannibalism can strongly promote coexistence between a prey and a predator, especially if the mortality rate de to cannibalism is sfficiently larger than the heterospecific predation rate [27], as a reslt, we investigate the effect of external spatial temporal noise on Tring patterns in a predator-prey model with diffsion and Holling type II fnctional response, with prey cannibalism. We model these pertrbations as said, by means of independent Gassian white noises. In this present paper, we only assme some control parameter of the system is affected by environmental flctations sch as epidemics, weather and natre disasters. We consider flctations presented in by the coefficient of the cannibalistic term of the prey, which is define as mortality rate de to cannibalism. Hence by introdcing a stochastic factor both in space and time ηx, t into the coefficient of the cannibalistic term c, the parameter c is treated as a Gassian random variable with mean ĉ and intensity ɛ. 6 c ĉ + ɛηx, t,

19 PREY CANNIBALISM 9 We therefore introdce eqation 6 in order to extend model system 2 sch that we obtain a system of Itô s stochastic differential eqations in the form, where for simplicity, we redefine mean ĉ to be the parameter c. : 7a 7b t = D + + c v t = D v v + δv v. v + αv c2 + d ɛ2 ηx, t, + d 5.. Nmerical Experiment. In this section, we discretize nmerically the stochastic predator-prey model 7 with white noise. We redefine the noise term ηx, t as 2 W x, t/ x t, where W x, t is a one-dimensional Brownian motion. The spatial approximation is constrcted from a Chebychev collocation scheme[3, 9]. Ths the spatial approximation is constrcted as a linear combination of the interpolating splines on the Gass-Lobatto qadratre x j, with a spatial domain size of N=256. This redces model system 2 to a system of Itô s stochastic ordinary differential eqation. The reslting system is then integrated in time sing an implicit Milsten scheme for approximating Itô s stochastic ordinary differential eqation, with a time stepsize of t =.. This therefore redce the Itô s stochastic ordinary differential eqation to the form 8 n+ j n j = 9 v n+ j v n j = D N i= D v N i= D 2 ji n+ i + n+ j + c n+ D 2 jiv n+ i + δv n+ j ɛn j 2 n j + d W j ɛ2 2 j v n+ j n+ j j n+ vn+ j n+ t, j j + αv n+ n j 3 n j + 2d n j + d3 cn+ j 2 n+ j + d W j 2 t where D 2 ji is known as the 2 nd order spectral differentiation matrix and n j = x j, t n, whereas vj n = vx j, t n, for j =,..., N. Also W j = z j t, where z j N, is normally distribted random nmbers with mean and standard deviation. The implicit Milsten scheme sed to advance the initial vale problem in time necessitates solving a nonlinear system of eqations at each time step. This is achieved throgh a Newton-Raphson method The effect of Noisy cannibalism rate on Tring Instability. In this section we investigate nmerically by simlation the effects of noise on Tring patterns at time t = 5. Extensive testing was performed throgh nmerical integration to describe model 7, and the reslts are shown in this section. This is done by sing an initial condition defined in a small pertrbation arond the positive homogenos steady state given, which is as = + α cos 2 nxx > x < π, v = v + α 2 cos 2 nxx > x < π,, t

20 2 BASHEER, QUANSAH, BHOWMICK, PARSHAD where α =.2 and α2 =.2. Dring the simlations, different types of dynamics are obtained for both the prey and predator. Conseqently in this section, we can restrict or analysis of pattern formation to prey poplation since mch of the noise in model system 7 affect the prey poplation. In or nmerical simlation with a specific parameter set in the Tring domain as show in Fig.8, we show the evoltion of fixed spatial pattern of prey poplation from deterministic model model system 2 to stochastic model model system 7 where with noise we consider =. and.4 with small random pertrbation of the stationary soltion, v as given above. It is nmerically evidenced from Fig.8 first plot =, middle plot =. and last plot =.4 that large variety of distinct patterns can be obtained by making small changes to the level of intensity in the stochastic model, which can also case the disappears of Tring patterns as well as the pattern formed can be enhance by small intensity levels. Therefore in or experiments it was realized that for a given parameter set in the Tring domain, there exist a specific intensity level of the noise that can case the average prey poplation species to decrease significantly as well as the disappearance of pattern formation. We next investigate whether noise can indce pattern formation. This achieved t t x x Figre 8 by sing a non-tring parameter set. Fig.9 shows the dynamic behavior of model 7 with non-tring parameter set. The first plot of Fig.9 shows the plot of the deterministic model where with diffsion the prey poplation still goes into a steady statehence no Tring Patterns formed, the proceeding plot of Fig.9 shows the evoltion from non-tring patterns to some form of weak patterns as changes from =. to = x Figre 9 Clearly from above, it is comprehensible that noise can indce instability, and the effect of noise on model system 7 is drastic. Ths noise can enhance, as well as case the disappearance of patterns formed in the deterministic model. These behaviors confirm how sensitive model system 7 is to different noise intensities.

21 PREY CANNIBALISM 2 Table. Parameters sed in the simlations. In all figres. Fig. α c c d δ D D v 2,3, ,7a b Discssion and Conclsion Cannibalism, althogh biqitos in both predator and prey commnities, has primarily been stdied as a phenomenon in the predator species. In the crrent manscript we have considered the effect of both predator and prey cannibalism on a Holling Tanner model, with ratio dependent fnctional response. Or reslts show that in the ODE version of the model, in certain parameter regime, predator cannibalism is able to stabilize the nstable interior eqilibrim, that is the eqilibrim is nstable if there is no cannibalism. This stability reslt depends critically on the cannibals feeding rate of prey γ, and the ratio. That is if the cannibalistic predator eats less prey so γ is small in comparison to, then only a little bit of feeding of conspecifics is enogh to stabilize the dynamics. This is seen via lemma 2., and is in accordance with known facts abot predator cannibalism in other models. However, when we consider prey cannibalism, we see that it is not able to stabilize the nstable interior eqilibrim. This is rigorosly proven for small prey cannibalism rate c, via lemma 2.3. We conjectre this is tre even for large cannibalism rate. This reslt tells s that for parameters sch that the no cannibalism model has an nstable eqilibrim, or essentially cyclical dynamics, no amont of prey cannibalism can stabilize this. From an ecological perspective, this tells s that in a predator-prey commnity with cyclical flctation in the predator and prey poplations so essentially limit cycle dynamics, cannibalism in the prey wold jst maintain the cyclically, and the oscillations in poplation cannot be driven to a stable steady state. Interestingly, this is seen to hold even if we assme some stochasticity in the prey cannibalism rate, see figs We also see that prey cannibalism can lead to limit cycle dynamics, bt the niqeness or non niqeness of limit cycles is an open qestion. It wold be extremely interesting to investigate this frther. In particlar one may ask, if changing the form of the prey cannibalism term C may lead to mltiple limit cycles, or render the limit cycle niqe. Another interesting open problem, wold be to consider the effect of both predator and prey cannibalism together, on models of Holling Tanner type. In particlar, one might try to classify regions, where they have conflicting effects, that is say predator cannibalism stabilises while prey cannibalism destabilises, and to then consider the net effect, in terms of the parameters concerned. To the best of or knowledge, this has not been considered in any work in the literatre. Also, of mch interest lately is the problem of a diseased predator

22 22 BASHEER, QUANSAH, BHOWMICK, PARSHAD in predator-prey models with cannibalistic predators. It is seen that cannibalism in the predator, can prevent disease in predator from spreading, as the cannibal cold wipe ot the diseased predator [3]. It wold be interesting to consider sch a system, with cannibalism now in the prey, and also a diseased prey. Here one can ask if prey cannibalism cold stop disease in the prey species, as is known from similar models in predator cannibalism. In case of the spatially explicit model, we ncover certain interesting reslts as well. What we observe is that for small predator growth rate, or if prey handling time and prey captre rate are close, Tring patterns are an impossibility in the no cannibalism model. That is spatial segregation cannot occr in. However, even if one maintains this small predator growth rate, or if prey handling time and prey captre rate are close, if the prey cannibalism rate is chosen in the right regime, see theorem 4.3, then one can still have spatial patterns in the prey cannibalism model 2. From a spatial ecology perspective this is qite interesting, as this is telling s that if the predator growth rate is crtailed for some reason, in the Holling Tanner model, spatial segregation cannot occr, nless prey cannibalism takes place, in which case spatial segregation can occr. Interestingly, if one compares the destabilizing inflence of prey cannibalism in the ODE and PDE models, we find that if the cannibalism rate c < 2 nc c, then it will destabilize in the ODE case. However, to destabilize in the PDE case, it cannot be to small, and mst be sch that + α 2 c + d 2 c > c > d c + α 2 c δ c + d 2 as seen in theorem 4.3. Alternatively, we show that predator cannibalism cannot case Tring patterns if the predator feeding rate is crtailed, according to the parametric restriction γ <, as seen in lemma 4.5. Lastly, ecosystems are always sbject to random forcing. However there have been no stdies with a random cannibalism rate, even as far as predator cannibalism goes, let alone prey cannibalsim. Or reslts with a random cannibalism rate show that it can affect both limit cycle dynamics, as well as completely alter the spatial Tring patterns, that form in the deterministic case, see figs In case of limit cycle dynamics, noisy cannibalism makes things aperiodic, bt cannot remove the poplation cycles completely. In the case of the spatially explicit model, the noise is seen to both be able to create patterns, as well as to destroy them. All in all, we hope that the crrent work leads to many frther investigations, into the fascinating effects of prey cannibalism on predator prey dynamics. Sch endeavors will bring s closer to nderstanding the fll extent of cannibalism as it occrs in natre, and the role it plays in shaping predator-prey commnities. 7. Acknowledgements We wold like to acknowledge very helpfl conversations with Professor Volker Rdolf, in the department of Biosciences at Rice University, on varios ecological concepts of prey cannibalism, and its mathematical modeling. References [] Peter A Abrams and Lev R Ginzbrg. The natre of predation: prey dependent, ratio dependent or neither? Trends in Ecology & Evoltion, 58:337 34, 2. d c

23 PREY CANNIBALISM 23 [2] Malay Banerjee and Santo Banerjee. Tring instabilities and spatio-temporal chaos in ratiodependent holling tanner model. Mathematical biosciences, 236:64 76, 22. [3] Santosh Biswas, Samrat Chatterjee, and Joydev Chattopadhyay. Cannibalism may control disease in predator poplation: reslt drawn from a model based stdy. Mathematical Methods in the Applied Sciences, 24. [4] B Bonomo, D Lacitignola, and S Rionero. Effect of prey growth and predator cannibalism rate on the stability of a strctred poplation model. Nonlinear Analysis: Real World Applications, 2:7 8, 2. [5] Brno Bonomo and Deborah Lacitignola. On the stabilizing effect of cannibalism in stagestrctred poplation models. Mathematical biosciences and engineering: MBE, 34:77 73, 26. [6] Alvin Aaden Yim-Hol Chan, Palina Giraldo-Perez, Sonja Smith, and Daniel T Blmstein. Anthropogenic noise affects risk assessment and attention: the distracted prey hypothesis. Biology Letters, 64:458 46, 2. [7] Ynshyong Chow and Sophia R-J Jang. Cannibalism in discrete-time predator prey systems. Jornal of biological dynamics, 6:38 62, 22. [8] David Claessen, Andre M De Roos, and Lennart Persson. Poplation dynamic theory of size dependent cannibalism. Proceedings of the Royal Society of London B: Biological Sciences, 27537:333 34, 24. [9] Wai Sn Don and Alex Solomonoff. Accracy and speed in compting the chebyshev collocation derivative. SIAM Jornal on Scientific Compting, 66: , 995. [] Stefano Fasani and Sergio Rinaldi. Remarks on cannibalism and pattern formation in spatially extended prey predator systems. Nonlinear Dynamics, 674: , 22. [] Clinton D Francis, Catherine P Ortega, and Alexander Crz. Noise polltion changes avian commnities and species interactions. Crrent biology, 96:45 49, 29. [2] Ph Getto, O Diekmann, and AM De Roos. On the dis advantages of cannibalism. Jornal of mathematical biology, 56:695 72, 25. [3] David Gottlieb and Livi Lstman. The spectrm of the chebyshev collocation operator for the heat eqation. SIAM jornal on nmerical analysis, 25:99 92, 983. [4] Sze-Bi Hs and Tzy-Wei Hang. Global stability for a class of predator-prey systems. SIAM Jornal on Applied Mathematics, 553: , 995. [5] Jan Kelkel and Christina Srlesc. On a stochastic reaction diffsion system modeling pattern formation on seashells. Jornal of mathematical biology, 66: , 2. [6] Debra Kelly. The cltre that still practices cannibalism, October 23. [7] C Kohlmeier and W Ebenhöh. The stabilizing role of cannibalism in a predator-prey system. Blletin of Mathematical Biology, 573:4 4, 995. [8] PH Leslie and JC Gower. The properties of a stochastic model for the predator-prey type of interaction between two species. Biometrika, pages , 96. [9] Li Li and Zhi-Jn Wang. Dynamics in a predator prey model with space and noise. Applied Mathematics and Comptation, 292: , 23. [2] Philip K Maini, Thomas E Woolley, Rth E Baker, Eamonn A Gaffney, and S Seirin Lee. Tring s model for biological pattern formation and the robstness problem. Interface focs, page rsfs23, 22. [2] James D Mrray. Mathematical biology i: An introdction, vol. 7 of interdisciplinary applied mathematics, 22. [22] Akira Okbo and Smon A Levin. Diffsion and ecological problems: modern perspectives, volme 4. Springer Science & Bsiness Media, 23. [23] CR Pennell. Cannibalism in early modern north africa. British Jornal of Middle Eastern Stdies, 82:69 85, 99. [24] Gary A Polis. The evoltion and dynamics of intraspecific predation. Annal Review of Ecology and Systematics, pages , 98. [25] Associated Press. Indian doc focses on hind cannibal sect, October 25. [26] Volker HW Rdolf. Conseqences of stage-strctred predators: cannibalism, behavioral effects, and trophic cascades. Ecology, 882:299 33, 27. [27] Volker HW Rdolf. The interaction of cannibalism and omnivory: conseqences for commnity dynamics. Ecology, 88: , 27. [28] Volker HW Rdolf. The impact of cannibalism in the prey on predator-prey systems. Ecology, 89:36 327, 28.