Bifurcation and global stability in an eco-epidemic model with refuge

Abstract

In this work, we formulate a predator–prey–pathogen model in which the predator is specialist in nature and infected prey can undergo refugia of constant size to avoid predator attack. To investigate the predation effect on the epidemics, we take a situation where the predator eats infected prey only. This is in accordance with the fact that the infected individuals are less active and can be caught more easily. Though it is a well-known fact that consumption of infected prey may harm predator population, the opposite holds in few cases. This leads to a controlled measure of disease prevalence. As predator consumes a particular type of prey species, Holling type II functional response is appropriate. This corroborates to the specialist type of predator. For biological validity of the model, boundedness of the system is studied. The dynamical behavior of the model has been analyzed throughly. Model analysis shows that all the population remains in coexistence when predator consumes the infected prey rather than the susceptible one. The results establish the fact that the effects of refuge used by prey decrease the equilibrium density of susceptible prey population, whereas the opposite holds for infected prey population. However, equilibrium density of predator may decrease or increase by increasing the amount of prey refuge. Global stability of the coexistence equilibrium point is developed by using Li and Muldowney’s high-dimensional Bendixson’s criterion. Numerical simulations are performed to validate our theoretical results.

Keywords

Eco-epidemic model / Prey refuge / Global stability / Hopf bifurcation / Persistence

References

[1] Bera SP, Maiti A, Samanta GP. A prey–predator model with infection in both prey and predator. Filomat, 2015, 29(8): 1753-1767

[2] Berryman AA, Lima M, Hawkins BA. Population regulation, emergent properties, and a requiem for density dependence. Oikos, 2002, 99: 600-606

[3] Birkhoff G, Rota CG. Ordinary differential equation, 1982 Boston Ginn and Co.

[4] Collings JB. Bifurcation and stability analysis of a temperature-dependent mite predator-prey interaction model incorporating a prey refuge. Bull Math Biol, 1995, 57: 63-76

[5] Dobson AP. The population biology of parasite-induced changes in host behavior. Q Rev Biol, 1988, 63: 139-165

[6] Duffy MA, Sivars-Becker L. Rapid evolution and ecological host–parasite dynamics. Ecol Lett, 2007, 10: 44-53

[7] Dutta H. Insights into the impacts of four current environmental problems on flying birds. Energy Ecol Environ, 2017, 2(5): 329-349

[8] Fan M, Michael Y, Wang K. Global stability of an SEIS epidemic model with recruitment and a varying total population size. Math Biosci, 2001, 170: 199-208

[9] Flick AJ, Acevedo MA, Elderd BD. The negative effects of pathogen-infected prey on predators: a meta-analysis. Oikos, 2016, 125(11): 1554-1560

[10] Gause GF, Smaragdova NP, Witt AA. Further studies of interaction between predators and prey. J Anim Ecol, 1936, 5: 1-18

[11] Gonzalez-Olivares E, Ramos-Jiliberto R. Dynamic consequences of prey refuges in a simple model system: more prey, fewer predators and enhanced stability. Ecol Model, 2003, 166: 135-146

[12] Hadeler KP, Freedman HI. Predator–prey populations with parasitic infection. J Math Biol, 1989, 27: 609-631

[13] Haque M, Venturino E. Modelling disease spreading in symbiotic communities, wildlife destruction, conservation and biodiversity, 2009 New York Nova Science Publishers

[14] Hassell MP. The dynamics of arthropod predator–prey systems, 1978 Princeton Princeton University Press

[15] Hassell MP, May RM. Stability in insect host–parasite models. J Anim Ecol, 1973, 42: 693-726

[16] Hethcote H, Wang W, Han L, Ma Z. A predator–prey model with infected prey. Theor Popul Biol, 2004, 66: 259-268

[17] Hochberg ME, Holt RD. Refuge evolution and the population dynamics of coupled host–parasitoid associations. Evol Ecol, 1995, 9: 633-661

[18] Holt RD, Roy M. Predation can increase the prevalence of infectious disease. Am Nat, 2007, 169: 690-699

[19] Holmes JC, Bethal WM (1972) Modifications of intermediate host behavior by parasites. In: Canning EV, Wright CA (eds) Behavioural aspects of parasite transmission. Suppl I to Zool. f Linnean Soc vol 51, pp 123–149

[20] Hoy MA (1985) Almonds (California). In: Helle W, Sabelis MW (eds) Spider mites: their biology, natural enemies and control, World Crop Pests, vol 1B. Elesvier, Amsterdam, pp 229–310

[21] Hsieh YH, Hsiao CK. Predator–prey model with disease infection in both populations. Math Med Biol, 2008, 25: 247-266

[22] Huang Y, Chen F, Zhong L. Stability analysis of a prey–predator model with Holling type III response function incorporating a prey refuge. Appl Math Comput, 2006, 182: 672-683

[23] Hutson V. A theorem on average Lyapunov function. Monatsh Math, 1984, 98: 267-275

[24] Hudson PJ, Dobson AP, Newborn D. Do parasite make prey vulnerable to predation Red grouse and parasite. J Anim Ecol, 1992, 61: 681-692

[25] Johnson PTJ, Stanton DE, Preu ER, Forshay KJ, Carpentar SR. Dining on disease: how interactions between infection and environment affect predation risk. Ecology, 2006, 87: 1973-1980

[26] Kabata Z. Parasites and disease of fish cultured in the tropics, 1985 London Taylor and Francis

[27] Kar TK. Modelling and analysis of a harvested prey–predator system incorporating a prey refuge. J Comput Appl Math, 2006, 185: 19-33

[28] Kermack W, McKendrick A. A contribution to the mathematical theory of epidemics. Proc R Soc A, 1927, 115: 700-721

[29] Krivan V. Effects of optimal anti predator behavior of prey on predator–prey dynamics: the role of refuges. Theor Popul Biol, 1998, 53: 131-142

[30] Krivan V. On the Gause predator–prey model with a refuge: a fresh look at the history. J Theor Biol, 2011, 274: 67-73

[31] Lafferty KD. Foraging on prey that are modified by parasites. Am Nat, 1992, 140: 854-867

[32] Liao CM, Yeh CH, Chen SC. Predation affects the susceptibility of hard clam M e r e t r i x l u s o r i a \documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$Meretrix \, lusoria$$\end{document} to Hg-stressed birnavirus. Ecol Model, 2008, 210: 253-262

[33] Li M, Muldowney J. On Bendixson’s criterion. J Differ Equ, 1993, 106: 27-39

[34] Maynard Smith J. Models in ecology, 1974 Cambridge Cambridge University

[35] McNair JN. The effects of refuges on predator–prey interactions: a reconsideration. Theor Popul Biol, 1986, 29: 38-63

[36] Mondal R, Kesh D, Mukherjee D. Influence of induced plant volatile and refuge in tritrophic model. Ener Ecol Environ, 2018, 3(3): 171-184

[37] Moore J. Parasites and the behaviour of animals, 2002 Oxford Oxford University Press

[38] Mukherjee D. The effect of prey refuges on a three species food chain model. Differ Equ Dyn Syst, 2014, 22: 413-426

[39] Mukherjee D. The effect of refuge and immigration in a predator–prey systems in the presence of a competitor for the prey. Nonlinear Anal Real World Appl, 2016, 31: 277-287

[40] Mukherjee D. Study of refuge use on a predator-prey system with a competitor for the prey. Int J Biol Math, 2017, 10(2): 1750023

[41] Packer C, Holt RD, Hudson PJ, Lafferty KD, Dobson AP. Keeping the herds healthy and alert: implications of predator control for infectious disease. Ecol Lett, 2003, 6: 797-802

[42] Pal AK, Samanta GP. A ratio-dependent eco-epidemiological model incorporating a prey refuge. Univ J Appl Math, 2013, 1(2): 86-100

[43] Ruxton GD. Short term refuge use and stability of predator–prey models. Theor Popul Biol, 1995, 47: 1-17

[44] Sharma S, Samanta GP. A Leslie–Gower predator–prey model with disease in prey incorporating a prey refuge. Chaos Slitons Fractals, 2015, 70: 69-84

[45] Sih A. Prey refuges and predator–prey stability. Theor Popul Biol, 1987, 31: 1-12

[46] Sotomayor J. Peixoto MM. Generic bifurcations of dynamical systems. Dynamical systems, 1973 New York Academic Press 549-560

[47] Sun C, Loreau M. Dynamics of a three-species food chain model with adaptive traits. Chaos Solitons Fractals, 2009, 41: 2812-2819

[48] Tewa JJ, Djeumen VY, Bowong S. Predator–prey model with Holling response function of type II and SIS infectious disease. Appl Math Model, 2013, 37: 4825-41

[49] Van Dobben WH. Sublethal effects of three ectoparasites on fish. J Fish Biol, 1952, 7: 283-294

[50] Venturino E. The influence of diseases on Lotka-Volterra systems. Rocky Mt J Math, 1994, 24: 381-402

[51] Xiao YN, Chen LS. Analysis of a three species eco-epidemiological model. J Math Anal Appl, 2001, 258(2): 733-754

[52] Zhou XY, Shi XY, Song XY. Analysis of a delay prey–predator model with disease in the prey species only. J Korean Math Soc, 2009, 46(4): 713-731