Skip Navigation Links
Journal of Environmental Accounting and Management
António Mendes Lopes (editor), Jiazhong Zhang(editor)
António Mendes Lopes (editor)

University of Porto, Portugal


Jiazhong Zhang (editor)

School of Energy and Power Engineering, Xi'an Jiaotong University, Xi'an, Shaanxi Province 710049, China

Fax: +86 29 82668723 Email:

Global Bifurcation Analysis of Higher Codimension in a Ratio Dependent Predator-Prey Model with Prey Refuge and Constant Yield-Harvesting of Both Species

Journal of Environmental Accounting and Management 11(4) (2023) 375--401 | DOI:10.5890/JEAM.2023.12.002

Nawaj Sarif, Sahabuddin Sarwardi

Download Full Text PDF



In this paper, a modified ratio-dependent Bazykin prey-predator model with prey refuge and constant harvesting of both species has been considered. The significant mathematical attributes of the proposed model are studied with the assistance of equilibria, local stability analysis, and bifurcation theory. The model displays a complex dynamics within the prey-predator plane. The existence and stability of possible equilibria are examined. The existence of Hopf bifurcation is obtained by analyzing the associated characteristic equation. The direction of the Hopf-bifurcating periodic solution and its stability are determined by calculating the first Lyapunov coefficient. Also, it has been proven that the proposed model exhibits Bogdanov-Takens bifurcation of codimension two for an appropriate choice of bifurcation parameters and calculating the conventional form. Furthermore, some numerical examples are provided to illustrate the theoretical finding. Lastly, discuss their ecological interpretations.


  1. [1]  Anderson, R.M. and May, R.M. (1979), Population biology of infectious diseases: Part I, Nature, 280, 361-367.
  2. [2]  Hastings, A. and Powell, T. (1991), Chaos in three-species food chain, Ecology, 72, 896-903.
  3. [3]  Turchin, P. (2003), Complex Population Dynamics: A Theoretical/Empirical Synthesis, Princeton University Press, Princeton.
  4. [4]  Lotka, A. (1925), Elements of Physical Biology, Baltimore, Williams and Williams.
  5. [5]  Volterra, V. (1926), Fluctuation in the abundance of species considered mathematically, Nature CXVII, 558-560.
  6. [6]  Lu, M. and Huang, J. (2021), Global analysis in Bazykin's model with Holling II functional response and predator competition, Journal of Differential Equations 280, 99-138.
  7. [7]  Xiang, C., Huang, J., Ruan, S., and Xiao, D. (2020), Bifurcation analysis in a host-generalist parasitoid model with Holling II functional response, Journal of Differential Equations, 268, 4618-4662.
  8. [8]  Molla, H., Rahman, Md.S., and Sarwardi, S. (2019) Dynamics of a predator-prey model with Holling type II functional response incorporating a prey refuge depending on both the species, International Journal of Nonlinear Sciences and Numerical Simulation, 20(1), 89-104.
  9. [9]  Benjamin, D.D., Thomann, E., Medlock, J., and De Leenheer, P. (2020), Global analysis of a predator-prey model with variable predator search rate, Journal of Mathematical Biology, 81, 159-183.
  10. [10]  Arditi,R. and Ginzburg, L. R. (1989), Coupling in predator-prey Dynamics: Ratio-depedence, Journal of theoretical Biology, 139, 311-326.
  11. [11]  Arditi, R. and Akcakaya, H.R. (1990) Underestimation of the mutual interference of predators, Oecologia (Berlin), 83, 358-361.
  12. [12]  Abrams, P.A. (1994), The fallacies of ``ratio-dependent'' predation, Ecology, 75, 1842-1850.
  13. [13]  Sarnelle, O., and Wilson, A.E. (2008), Type III functional response in Daphnia, Ecology, 89, 1723-1732.
  14. [14]  Bishop, M.J., Kelaher, B.P., Smith, M., York, P.H., and Booth, D.J. (2006), Ratio-dependent response of a temperate Australian estuarine system to sustained nitrogen loading, Oecologia, 149, 701-708.
  15. [15]  Arditi, R. and Berryman, A.A. (1991), The biological control paradox, Trends in Ecology and Evolution, 6.
  16. [16]  Kuang, Y. and Beretta, E. (1998), Global qualitative analysis of a ratio-dependent predator-prey system, Journal of Mathematical Biology, 36, 389-406.
  17. [17]  Sen, M., Banerjee, M. and Morozov, A., (2012), Bifurcation analysis of a ratio-dependent prey-predator model with Allee effect, Ecological Complexity, 11, 12-27.
  18. [18]  Berezovskaya, F., Karev, G., and Arditi, R. (2001), Parametric analysis of the ratio-dependent predator-prey model, Journal of Mathematical Biology, 43, 221-246.
  19. [19]  Makinde, O.D. (2007), Solving ratio-dependent predator-prey system with constant effort harvesting using Adomian decomposition method, Applied Mathematics and Computation, 186, 17-22.
  20. [20]  Clark, C.W. (1990), The optimal management of renewable resources, Mathematical Bioeconomics, John Wiley and Sons, New york, Toronto.
  21. [21]  Biswas, Md.H.A., Hossain, Md.R., and Mondal, M.K. (2017), Mathematical modelling applied to sustainable management of marine resources, Procedia Engineering, 194, 337-344.
  22. [22]  Huang, J.C., Gong, Y.J., and Ruan, S.G. (2013), Bifurcation analysis in a predator-prey model with constant-yield predator harvesting, Discrete and Continuous Dynamical Systems - Series B, 18, 2101-2121.
  23. [23]  Zhu, Q. and Peng, H.Q. (2019) Bifurcation analysis of a stage-structured predator-prey model with prey refuge, Discrete and Continuous Dynamical Systems - Series B, 12, 2195-2209.
  24. [24]  Adak, D., Bairagi, N., and Hakl, R. (2020), Chaos in delay-induced Leslie-Gower prey-predator-parasite model and its control through prey harvesting, Nonlinear Analysis: Real World Applications, 51, 102998.
  25. [25]  Choi, J. and Kim, Y.J. (2018), Predator-prey equations with constant harvesting and planting, Journal of Theoretical Biology, 458, 47-57.
  26. [26]  Luo, J. and Zhao, Y. (2017), Stability and bifurcation analysis in a predator-prey system with constant harvesting and prey group defense, International Journal of Bifurcation and Chaos, 27, 750179.
  27. [27]  Lee, J. and Baek, H. (2017), Dynamics of a Beddington-DeAngelis-type predator-prey system with constant rate harvesting, Electronic Journal of Qualitative Theory of Differential Equations, 1, 1-20.
  28. [28]  Xiao, D. and Ruan, S. (1999), Bogdanov-Takens bifurcations in predator-prey systems with constant rate harvesting, Fields Institute Communications, 21, 493-506.
  29. [29]  Haque, M. and Sarwardi, S. (2016), Effect of toxicity on a harvested fishery model, Modeling Earth Systems and Environment, 2, 1-10.
  30. [30]  Haque, M.M. and Sarwardi, S. (2018), Dynamics of a harvested prey-predator model with prey refuge dependent on both species, International Journal of Bifurcation and Chaos, 28(12), 1-16.
  31. [31]  Liu, X., and Huang, O. (2020), Analysis of optimal harvesting of a predator-prey model with Holling type IV functional response, Ecological complexity, 42, 100816.
  32. [32]  Sun, G. and Dai, B. (2020), Stability and bifurcation of a delayed diffusive predator-prey system with food-limited and nonlinear harvesting, Mathematical Biosciences and Engineering, 17, 3520-3552.
  33. [33]  Hu, D.P. and Cao, H.J. (2017), Stability and bifurcation analysis in a predator-prey system with Michaelis-Menten type predator harvesting, Nonlinear Analysis: Real World Applications, 33, 58-82.
  34. [34]  Bazykin, A.D. (1998), Non-linear dynamics of interacting populations, World Scientific Nonlinear Science Series A.
  35. [35]  May, R.M., Beddington, J.R., Clark, C.W., Holt, S.J. and Laws, R.M. (1979), Management of multispecies fisheries, Science, 205, 267-277.
  36. [36]  Ross, S.L. (1984), Differential Equations, John Wiley \& Sons, U.K.
  37. [37]  Murray, J.D. (1993), Mathematical Biology, Springer, New York.
  38. [38]  Perko, L. (1996), Differential equations and dynamical systems, Springer, New York.
  39. [39]  Bogdanov, R.I. (1976), Bifurcations of a limit cycle for a family of vector fields on the plane, Trudy Sem. Petrovsk. Vyp., 2, 23-35.
  40. [40]  Bogdanov, R.I. (1976), The versal deformations of a singular point on the plane in the case of zero eigenvalues, Trudy Sem. Petrovsk. Vyp., 2, 37-65.
  41. [41]  Takens, F. (1974), Forced oscillations and bifurcation, in \textquotedblleft{Applications of Global Analysis, I}\textquotedblright (Sympos., Utrecht State Univ., Utrecht), Math. Inst. Rijksuniv. Utrecht., 1-59.
  42. [42]  Kuznetsov, Y.A. (2013), Elements of Applied Bifurcation Theory, Springer Science \& Business Media.
  43. [43]  Sarwardi, S., Haque, Md. M. and Hossain, S. (2020), Analysis of Bogdanov-Takens bifurcations in a spatiotemporal harvested-predator and prey system with Beddington-DeAngelis type response function, Nonlinear Dynamics, 100, 1755--1778.
  44. [44]  Chow, S.N. and Hale, J.K. (1983), Methods of Bifurcation Theory, Springer, New York.
  45. [45]  Chow, S.N., Li, C., and Wang, D. (1994), Normal Forms and Bifurcation of Planar Vector Fields, Cambridge University Press, New York.
  46. [46]  Guckenheimer, J. and Holmes, P. (1983), Nonlinear Oscillations, Dynamical Systems, and Bifurcation of Vector Fields, Springer-Verlag, New York.
  47. [47]  Roy, M., Mandal, S., and Ray, S. (2008), Detrital ontogenic model including decomposer diversity, Ecological Modelling 215(1-3), 200-206.
  48. [48]  Sarwardi, S., Haque, M., and Mandal, P.K. (2012), Ratio-dependent predator-prey model of interacting population with delay effect, Nonlinear Dynamics, 69, 817-836.
  49. [49]  Sarif, N. and Sarwardi, S. (2021), Analysis of Bogdanov-Takens bifurcation of codimension two in a Gause-type model with constant harvesting of both species and delay effect, Journal of Biological Systems, 29, 741-771.
  50. [50]  Sarwardi, S., Haque, M. and Mandal, P.K. (2014), Persistence and global stability of Bazykin predator-prey model with Beddington-DeAngelis response function, Communications in Nonlinear Science and Numerical Simulation, 19(1), 189-209.