Abstract
In this paper, we investigate the global Hopf bifurcation of a diffusive modified Leslie–Gower predator–prey model with delay and Michaelis–Menten type prey harvesting. First, we obtain the stability of positive steady state and the existence of local Hopf bifurcation under certain conditions. Second, we get the permanence of the system by using the comparison theorem. Moreover, by constructing a suitable Lyapunov function, we derive sufficient conditions for the global attractivity of the unique positive steady state for the system without delay. Then, the global existence of positive periodic solutions is established by using the global Hopf bifurcation result of Wu. Finally, the results are verified by numerical simulation.
Similar content being viewed by others
Data Availability
We do not analyse or generate any datasets, because our work proceeds within a theoretical and mathematical approach.
References
Leslie, P., Gower, J.: The properties of a stochastic model for the predator–prey type of interaction between two species. Biometrica 47, 219–234 (1960)
Yadav, R., Mukherjee, N., Sen, M.: Spatiotemporal dynamics of a prey-predator model with Allee effect in prey and hunting cooperation in a Holling type III functional response. Nonlinear Dyn. 107(1), 1397–1410 (2021)
Lamontagne, Y., Coutu, C., Rousseau, C.: Bifurcation analysis of a predator–prey system with generalised holling type III functional response. J. Dyn. Differ. Equ. 20(3), 535–571 (2008)
Yang, R., Wei, J.: Stability and bifurcation analysis of a diffusive prey predator system in Holling type III with a prey refuge. Nonlinear Dyn. 79(1), 631–646 (2015)
Liu, Z., Zhong, S., Yin, C., Chen, W.: Dynamics of impulsive reaction–diffusion predator–prey system with Holling type III functional response. Appl. Math. Model. 35, 5564–5578 (2011)
Yang, W., Li, Y.: Dynamics of a diffusive predator–prey model with modified Leslie–Gower and Holling-type III schemes. Comput. Math. Appl. 65, 1727–1737 (2013)
Huang, J., Ruan, S., Song, J.: Bifurcations in a predator–prey system of Leslie type with generalized Holling type III functional response. J. Differ. Equ. 257, 1721–1752 (2014)
Zhang, Y., Meng, X.: Dynamic study of a stochastic Holling III predator-prey system with a prey refuge. J. Differ. Equ. 55–3, 73–78 (2022)
Liu, C., Chang, L., Huang, Y., Wang, Z.: Turing patterns in a predator–prey model on complex networks. Nonlinear Dyn. 99(4), 3313–3322 (2020)
Blyuss, K., Kyrychko, S., Kyrychko, Y.: Time-delayed and stochastic effects in a predator–prey model with ratio dependence and Holling type III functional response. Chaos 31(7), 073141 (2021)
Tiwari, V., Tripathi, J.P., Abbas, S., Wang, J.S., Sun, G.Q., Jin, Z.: Qualitative analysis of a diffusive Crowley–Martin predator-prey model: the role of nonlinear predator harvesting. Nonlinear Dyn. 98, 1169–1189 (2019)
Tripathi, J.P., Abbas, S., Sun, G.Q., Jana, D., Wang, C.H.: Interaction between prey and mutually interfering predator in prey reserve habitat: pattern formation and the Turing–Hopf bifurcation. J. Franklin. I(355), 7466–7489 (2018)
May, R.M., Beddington, J.R., Clark, C.W., Holt, S.J., Laws, R.M.: Management of multispecies fisheries. Science 205(4403), 267–277 (1979)
Clark, C.W.: Mathematical Bioeconomics: The Optimal Management of Renewable Resources. Wiley, New York (1976)
Clark, C.W., Mangel, M.: Aggregation and fishery dynamics: a theoretical study of schooling and the purse seine tuna fisheries. Fish. Bull. 77, 317–337 (1979)
Wang, W., Liu, Q., Jin, Z.: Spatiotemporal complexity of a ratio-dependent predator–prey system. Phys. Rev. E 75, 051913 (2007)
Zhang, L., Wang, W., Xue, Y.: Spatiotemporal complexity of a predator–prey system with constant harvest rate. Chaos Soliton. Fract. 41, 38–46 (2009)
Yuan, R., Jiang, W., Wang, Y.: Saddle-node-Hopf bifurcation in a modified Leslie–Gower predator–prey model with time-delay and prey harvesting. J. Math. Anal. Appl. 422(2), 1072–1090 (2014)
Yan, X., Zhang, C.: Global stability of a delayed diffusive predator–prey model with prey harvesting of Michaelis–Menten type. Appl. Math. Lett. 114, 106904 (2021)
Liu, B., Wu, R., Chen, L.: Patterns induced by super cross-diffusion in a predator–prey system with Michaelis–Menten type harvesting. Math. Biosci. 298, 71–79 (2018)
Hu, D., Cao, H.: Stability and bifurcation analysis in a predator–prey system with Michaelis–Menten type predator harvesting. Nonlinear Anal. Real. 33, 58–82 (2017)
Liu, W., Jiang, Y.: Bifurcation of a delayed Gause predator–prey model with Michaelis–Menten type harvesting. J. Theor. Biol. 438, 116–132 (2018)
Gupta, R., Peeyush, C.: Bifurcation analysis of modified Leslie–Gower predator–prey model with Michaelis–Menten type prey harvesting. J. Math. Anal. Appl. 398, 278–295 (2013)
Zhang, X., Zhao, H.: Stability and bifurcation of a reaction–diffusion predator–prey model with non-local delay and Michaelis-=Menten-type prey-harvesting. Int. J. Comput. Math. 93(9), 1447–1469 (2016)
Jiang, J., Li, X., Wu, X.: The dynamics of a bioeconomic model with Michaelis–Menten type prey harvesting. B. Malays Math. Sci. So. 46(2), 57 (2023)
Zuo, W., Ma, Z., Cheng, Z.: Spatiotemporal dynamics induced by Michaelis–Menten type prey harvesting in a diffusive Leslie–Gower predator–prey model. Int. J. Bifurc. Chaos 30(14), 2050204 (2020)
Liu, M., Hu, D., Meng, F.: Stability and bifuration analysis in a delay-induced predator–prey model with Michaelis–Menten type predator harvesting. Discrete Contin. Dyn. Syst. 14(9), 3197–3222 (2021)
Liu, W., Jiang, Y.: Bifurcation of a delayed Gause predator–prey model with Michaelis–Menten type harvesting. J. Theor. Biol. 438, 116–132 (2018)
Tiwari, V., Tripathi, J.P., Mishra, S., Upadhyay, R.K.: Modeling the fear effect and stability of non-equilibrium patterns in mutually interfering predator–prey systems. Appl. Math. Comput. 124948, 371 (2020)
Tiwari, V., Tripathi, J.P., Jana, D., Tiwari, S.K., Upadhyay, R.K:. Exploring complex dynamics of spatial predator-prey system: role of predator interference and additional food. Int. J. Bifurc. Chaos. 30(7), 2050102 (2020)
Kohlmeier, C., Ebenhoh, W.: The stabilizing role of cannibalism in a predator–prey system. B. Math. Biol. 57, 401–411 (1995)
Vaidyaa, N.K., Wu, J.H.: Modeling Spruce Budworm population revisited: impact of physiological structure on outbreak control. B. Math. Biol. 70, 769–784 (2008)
Yang, R., Nie, C., Jin, D.: Spatiotemporal dynamics induced by nonlocal competition in a diffusive predator–prey system with habitat complexity. Nonlinear Dyn. 110, 879–900 (2022)
Yang, R., Wang, F., Jin, D.: Spatially inhomogeneous bifurcating periodic solutions induced by nonlocal competition in a predator–prey system with additional food. Math. Methods Appl. Sci. 45(16), 9967–997 (2022)
Xu, C., Zhang, W., Aouiti, C., Liu, Z., Yao, Y.: Bifurcation insight for a fractional-order stage-structured predator–prey system incorporating mixed time delays. Math. Methods Appl. Sci. 46(8), 9103–9118 (2023)
Yang, R., Zhang, C.: Dynamics in a diffusive modified Leslie–Gower predator–prey model with time delay and prey harvesting. Nonlinear Dyn. 87(2), 863–878 (2016)
Yang, R., Wei, J.: The effect of delay on a diffusive predator–prey system with modified Leslie–Gower functional response. Bull. Malays. Math. Sci. Soc. 40(1), 51–73 (2015)
Tian, Y., Weng, P.: Asymptotic stability of diffusive predator-prey model with modified Leslie–Gower and Holling-type II schemes and nonlocal time delays. Ima. J. Math. Control. I(31), 1–14 (2014)
Chen, S., Shi, J., Wei, J.: Global stability and Hopf bifurcation in a delayed diffusive Leslie–Gower predator–prey system. Int. J. Bifurc. Chaos 22(3), 1250061 (2012)
Song, Y., Peng, Y., Wei, J.: Bifurcations for a predator–prey system with two delays. J. Math. Anal. Appl. 337(1), 466–479 (2008)
Li, Y., Wang, M.: Hopf bifurcation and global stability of a delayed predator–prey model with prey harvesting. Comput. Math. Appl. 69(5), 398–410 (2015)
Yan, X., Zhang, C.: Global stability of a delayed diffusive predator–prey model with prey harvesting of Michaelis–Menten type. Appl. Math. Lett. 114, 106904 (2021)
Xu, X., Liu, M.: Global Hopf bifurcation of a general predator-prey system with diffusion and stage structures. J. Differ. Equ. 269, 8370–8386 (2020)
Liu, M., Cao, J.: Global Hopf bifurcation in a phytoplankton–zooplankton system with delay and diffusion. Int. J. Bifurc. Chaos. 31(8), 2150114 (2021)
Xu, X., Wei, J.: Bifurcation analysis of a spruce budworm model with diffusion and physiological structures. J. Differ. Equ. 262(10), 5206–5230 (2017)
Su, Y., Wei, J., Shi, J.: Hopf bifurcation in a diffusive logistic equation with mixed delayed and instantaneous density dependence. J. Dyn. Differ. Equ. 24, 897–925 (2012)
Wu, J.: Theory and Applications of Partial Functional Differential Equations. Springer-Verlag, New York (1996)
Hsu, S., Hwang, T.: Global stability for a class of predator–prey systems. SIAM J. Appl. Math. 55, 763–783 (1995)
Du, Y., Hsu, S.: A diffusive predator–prey model: in heterogeneous environment. J. Differ. Equ. 203, 331–364 (2004)
Acknowledgements
This work is supported by the Natural Science Foundation for Heilongjiang Province (No. LH2021A018) and the National Natural Science Foundation of China (No. 12071115).
Author information
Authors and Affiliations
Contributions
All authors participated in the writing and coordination of the manuscript and all authors read and approved the final manuscript.
Corresponding author
Ethics declarations
Conflict of interest
The authors declare that they have no competing interests.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
This work is supported by the Natural Science Foundation of Heilongjiang Province (No. LH2021A018) and the National Natural Science Foundation of China (No. 12071115).
Appendix
Appendix
In this appendix, we provide the details of the proof of Theorem 2.1. To prove the existence and uniqueness of positive constant steady state, we only need to prove that \(\phi (u)=0\) has a unique positive root. If \(\mathbb {E}_1>0\), we have \(\phi (0)>0\). Since \(\phi (\infty )=-\infty \), \(\phi (u)=0\) has at least one positive root. From \(\phi '(0)=\mathbb {D}_{1}\) and \(\phi ''(0)=\mathbb {C}_{1}\), we know that if one of the following cases holds for \(u\in [0,+\infty )\), \(\phi (u)=0\) has a unique positive root.
- \((s_{1})\):
-
: \(\phi (u)\) monotonically decreasing;
- \((s_{2})\):
-
: \(\phi (u)\) increases and then decreases;
- \((s_{3})\):
-
: \(\phi (u)\) decreases and then increases and then decreases and either \(\phi (u_{1}^{2})>0\) or \(\phi (u_{1}^{1})<0\);
- \((s_{4})\):
-
: \(\phi (u)\) increases and then decreases and then increases and then decreases and either \(\phi (u_{1}^{2})>0\) or \(\phi (u_{1}^{3})>0\) and \(\phi (u_{1}^{1})<0\).
Next, we analyze \((H_{1a})\)-\((H_{1l})\) in turn. The following discussion satisfies \(u\in [0,+\infty )\).
\((H_{1a})\). If \(\triangle \le 0\), than we have \(\phi ''(u)<0\), which implies that \(\phi '(u)\) monotonically decreasing, and \(\phi (u)\) monotonically decreasing or increases and then decreases.
\((H_{1b})\). If \(\triangle >0\), \(\mathbb {B}_{1}\le 0\), \(\mathbb {C}_1>0\) and \(\mathbb {D}_{1}\ge 0\), this means \(\phi ''(u)=0\) has a positive root, \(\phi '(u)\) increases and then decreases and \(\phi (u)\) increases and then decreases.
\((H_{1c})\). If \(\triangle >0\), \(\mathbb {B}_{1}\le 0\) and \(\mathbb {C}_1>0\), than \(\phi ''(u)=0\) has a positive root, which implies \(\phi '(u)\) increases and then decreases. Under the conditions that \(\mathbb {D}_{1}<0\) and \(\phi '(u_{2}^{+})\le 0\), we have \(\phi '(u)<0\) and \(\phi (u)\) monotonically decreasing.
\((H_{1d})\). Under the conditions of \(\triangle >0\), \(\mathbb {B}_{1}\le 0\) and \(\mathbb {C}_1>0\), we obtain \(\phi ''(u)=0\) has a positive root, \(\phi '(u)\) increases and then decreases. If \(\mathbb {D}_{1}<0\) and \(\phi '(u_{2}^{+})>0\) hold, the maximum point of \(\phi '(u)\) to be greater than 0, \(\phi (u)\) decreases and then increases and then decreases.
\((H_{1e})\). In the case of \(\triangle >0\), \(\mathbb {B}_{1}\le 0\) and \(\mathbb {C}_1\le 0\), we have \(\phi ''(u)=0\) has no positive root, \(\phi '(u)\) monotonically decreasing and \(\phi (u)\) monotonically decreasing or increases and then decreases.
\((H_{1f})\). If \(\triangle >0\), \(\mathbb {B}_{1}>0\) and \(\mathbb {C}_1\ge 0\), than \(\phi ''(u)=0\) has a positive root and \(\phi '(u)\) increases and then decreases. Under the conditions of \(\mathbb {D}_1\ge 0\), \(\phi (u)\) increases and then decreases.
\((H_{1g})\). Under the conditions that \(\triangle >0\), \(\mathbb {B}_{1}>0\), \(\mathbb {C}_1\ge 0\), \(\mathbb {D}_{1}<0\) and \(\phi '(u_{2}^{+})\le 0\), we obtain \(\phi '(u)\le 0\), and \(\phi (u)\) monotonically decreasing.
\((H_{1h})\). If \(\triangle >0\), \(\mathbb {B}_{1}>0\), \(\mathbb {C}_1\ge 0\), \(\mathbb {D}_{1}<0\) and \(\phi '(u_{2}^{+})>0\), than the maximum point of \(\phi '(u)\) to be greater than 0. So \(\phi (u)\) decreases and then increases and then decreases.
\((H_{1i})\). In the case of \(\triangle >0\), \(\mathbb {B}_{1}>0\) and \(\mathbb {C}_1<0\), we have \(\phi '(u)\) decreases and then increases and then decreases. If \(\mathbb {D}_{1}>0\) and \(\phi '(u_{2}^{-})\ge 0\), the minimum point of \(\phi '(u)\) to be greater than 0, which implies that \(\phi (u)\) increases and then decreases.
\((H_{1j})\). Under the conditions that \(\triangle >0\), \(\mathbb {B}_{1}>0\), \(\mathbb {C}_1<0\), \(\mathbb {D}_{1}>0\), \(\phi '(u_{2}^{+})>0\) and \(\phi '(u_{2}^{-})<0\), we obtain the maximum point of \(\phi '(u)\) to be greater than 0, and the minimum point of \(\phi '(u)\) to be less than 0, this means \(\phi (u)\) increases and then decreases and then increases and then decreases.
\((H_{1k})\). If \(\triangle >0\), \(\mathbb {B}_{1}>0\) and \(\mathbb {C}_1<0\), this means \(\phi '(u)\) decreases and then increases and then decreases. In the case of \(\phi '(u_{2}^{+})\le 0\), we have the maximum point of \(\phi '(u)\) to be less than 0, and \(\phi (u)\) monotonically decreasing or increases and then decreases.
\((H_{1l})\). Under the conditions of \(\triangle >0\), \(\mathbb {B}_{1}>0\), \(\mathbb {C}_1<0\), \(\mathbb {D}_{1}\le 0\) and \(\phi '(u_{2}^{+})>0\), which inplies the maximum point of \(\phi '(u)\) to be greater than 0, and \(\phi (u)\) decreases and then increases and then decreases.
The proof is complete.
Rights and permissions
Springer Nature or its licensor (e.g. a society or other partner) holds exclusive rights to this article under a publishing agreement with the author(s) or other rightsholder(s); author self-archiving of the accepted manuscript version of this article is solely governed by the terms of such publishing agreement and applicable law.
About this article
Cite this article
Wang, K., Xu, X. & Liu, M. Global Hopf Bifurcation of a Diffusive Modified Leslie–Gower Predator–Prey Model with Delay and Michaelis–Menten Type Prey Harvesting. Qual. Theory Dyn. Syst. 23, 81 (2024). https://doi.org/10.1007/s12346-023-00939-0
Received:
Accepted:
Published:
DOI: https://doi.org/10.1007/s12346-023-00939-0