A predator–prey model with disease in the predator species only

https://doi.org/10.1016/j.nonrwa.2009.06.012Get rights and content

Abstract

In this paper we propose a predator–prey model with logistic growth in the prey population that includes an SIS parasitic infection in the predator population, with the assumption that the predator has an alternative source of food. For simplicity we initially work with a model involving the fractions of the predator which are susceptible and those infected and then translate the results back to the model with absolute numbers. Important thresholds R01, R02, R03 and R04 are identified and their implications have been explained. Our theoretical study indicates that the absence of prey may be beneficial for predator when a transmissible disease runs among the predator population. One important conclusion is that infection in the predator species may save the prey from extinction even if R02, the basic reproduction number for the prey to be able to invade the predator-only equilibrium, is less than one. Therefore infection in the predator species may be taken as biological control. Finally, analytical results have been supported by numerical simulations with the help of experimental data.

Introduction

Mathematical and statistical models can help us to identify key parameters which determine the rich dynamics of an ecological or epidemiological system. In the development of quantitative theory for interaction of predator and prey, mathematical ecology is also an important factor along with the experimental ecology, see [1] and references therein. Transmissible disease in an ecological situation is fast becoming a major field of study in its own right. The first mathematical description of contagious diseases has been formulated by Kermack and McKendric [2]. Both theoretical and experimental investigations in these two fields namely ecology and epidemiology progressed independently along the years, until the late eighties and early nineties. The ecological literature has increasingly emphasized the importance of parasites in shaping the dynamics of both plant and animal communities, [3], [4], [5].

Some models merging features of the two phenomena, i.e. the demographics of interacting species and an epidemic evolution in such a composite environment, were then considered, see [6], [7], [8], [9], [10], [11], [12], [13], [14], [15], [16], [17] and the references therein.

Most of the above models are based on the assumption that the infected prey is more vulnerable to predation except for those developed by Venturino [11], and Haque and Venturino [15]. In these two papers the authors considered the case when disease is spreading among the predator population. The model considered here also studies the situation where an epidemic runs among the predator population but differs from previous models; it includes the recovery rate and also takes into account also the stability of the positive interior equilibrium point. In addition, we assume that the predator has a logistic growth rate since it has sufficient resources for alternative foods; and it is argued that alternative food sources may have an important role in promoting the persistence of predator–prey systems. van Baalen et al. [18] have observed that the switching of the predator species has a significant contribution to the persistence of species in the predator–prey system. They also have analyzed the conditions for stability as well as long-term behavior of the system under bounded population fluctuations. This leads to specific hypotheses about which types of alternative food (in terms of nutritional quality, availability, and handling time) promote persistence and by which mechanism.

Section snippets

Predator–prey model with disease in the predator species only

Let us consider a predator–prey model, where P(t) represents the number of prey and H(t) denotes the number of predators at the time t. The classical well known model is given by dPdt=aP(1PM)cPH,dHdt=rH(1HK)+ecPH. Here a is the growth rate of the prey and M is the prey carrying capacity of the environment. The total predation rate is cHP and the conversion rate is denoted by e. Again r represents the growth rate of the predator and K is the predator carrying capacity of the environment.

Now

Stability analysis

Here we perform the stability and bifurcation analysis of the system (2.4). System (2.4) can be written in the form Ẋ=F(X) where X=(P,H,I)T, and the Jacobian matrix of the system is J=(j11cP(1I+qI)cPH(q1)ecH(1I+qI)j22ecPH(q1)ecI(1I)(q1)rθIKj33), where j11=a(12P/M)cH(1I+qI), j22=r(12H/K)+ecP(1I+qI) and j33=(β+ecP(q1))(12I)(γ+bθrH/K).

Behavior of the proportional system (2.4) and the absolute system (2.3)

The behavior of the system when one or more of the susceptible predator, infected predator and prey populations are initially absent has already been discussed in Section 3. For starting values for which all three types of individuals are initially present note that for all parameter values there is exactly one LAS equilibrium. The following remarks can be done.

Remark 2

For the proportional system (2.4) with P(0)>0,H(0)>0 and I(0)>0 we have the following results:

(i) For R02<1,R031 and ξ1>0 the system

The role of infection on the predator–prey system

To observe the influence of infection in the predator–prey system we set I=0 in (2.4), and obtain a simple Lotka–Volterra type predator–prey system, with density-dependent logistic growth for the prey population, where the predator has an alternative food source. However it is analyzed by Bazykin et al. [22], but here we report the results in terms of our system parameters for comparison purpose. This system has four biologically relevant equilibria namely : (i) Eˆ0=(0,0), (ii) Eˆ1=(M,0), (iii)

The effect of prey removal on the eco-epidemiological system

To observe the role of the predator in the eco-epidemiological system we put P=0 in our system (2.4) and get a reduced system which is well known and has been analyzed extensively by Diekmann and Heesterbeek [23], Hethcote [24] and Brauer and Castillo-Chavez [25]. Here the local stability results have been reported in terms of our system parameters for comparison purposes. This system has four equilibria, namely (i) Ĕ0=(0,0) which is always admissible, (ii) Ĕ2=(K,0), which is always

Numerical investigations

Many biological examples are given in the Table 2 of [19], here we take the example of rabbit/fox predator–prey system as in the UK, the European rabbit plays a vital role in maintaining bio-diversity and the fox (Vulpus vulpus) has become the main carrier of rabies in Europe. Now let the fox be infected by rabies. The disease is transmitted mostly as a result of an infected animal biting a healthy animal. The birth rate of foxes b, its natural death rate d, transmission coefficient of rabies β

Discussion

In this article we have proposed and analyzed a mathematical model that consists of three non-linear autonomous differential equations for three different populations namely prey (P), susceptible predator (H), infected predator (I). The conditions for existence and stability of the equilibria of the system have been given. The bifurcation situations have also been observed around important equilibrium points.

We have investigated that there are five epidemiological threshold quantities for our

Acknowledgments

I am very grateful to the referees for their comments. I am greatly indebted to John R. King (editor of IMA Journal of Mathematics in Medicine and Biology, associate editor of SIAM Journal of Applied Mathematics), from whom I have learned applications of mathematics to biology and medicine. Finally, I am also thankful to Professor Philip Maini, Mathematical Institute(CMB), University of Oxford, who saw the earlier version of the manuscript. Finally, thanks go to the UK-IERI for financial

References (26)

  • C. Gibson et al.

    The role of the hemiparasitic annual Rhinanthus minor in determining grassland community structure

    Oecologia

    (1992)
  • A. Dobson et al.

    Pathogens and the structure of plant communities

    Trends in Ecology & Evolution

    (1994)
  • K. Hadeler et al.

    Predator–prey populations with parasitic infection

    Journal of Mathematical Biology

    (1989)
  • Cited by (88)

    • Stability, collapse and hyperchaos in a class of tri-trophic predator–prey models

      2023, Physica A: Statistical Mechanics and its Applications
    View all citing articles on Scopus
    View full text