Bifurcation Analysis in a Harvested Modified Leslie–Gower Model Incorporated with the Fear Factor and Prey Refuge

Fear and prey refuges are two significant topics in the ecological community because they are closely associated with the connectivity of natural resources. The effect of fear on prey populations and prey refuges (proportional to both the prey and predator) is investigated in the nonlinear-type predator-harvested Leslie–Gower model. This type of prey refuge is much more sensible and realistic than the constant prey refuge model. Because there is less research on the dynamics of this type of prey refuge, the current study has been considered to strengthen the existing literature. The number and stability properties of all positive equilibria are examined. Since the calculations for the determinant and trace of the Jacobian matrix are quite complicated at these equilibria, the stability of certain positive equilibria is evaluated using a numerical simulation process. Sotomayor’s theorem is used to derive a precise mathematical confirmation of the appearance of saddle-node bifurcation and transcritical bifurcation. Furthermore, numerical simulations are provided to visually demonstrate the dynamics of the system and the stability of the limit cycle is discussed with the help of the first Lyapunov number. We perform some sensitivity investigations on our model solutions in relation to three key model parameters: the fear impact, prey refuges, and harvesting. Our findings could facilitate some biological understanding of the interactions between predators and prey.