Abstract

In this work, our aim is to investigate the impact of a non-Kolmogorov predator-prey-subsidy model incorporating nonlinear prey refuge and the effect of fear with Holling type II functional response. The model arises from the study of a biological system involving arctic foxes, lemmings, and seal carcasses. The positivity and asymptotically uniform boundedness of the solutions of the system have been derived. Analytically, we have studied the criteria for the feasibility and stability of different equilibrium points. In addition, we have derived sufficient conditions for the existence of local bifurcations of codimension 1. It is also observed that there is some time lag between the time of perceiving predator signals through vocal cues and the reduction of prey’s birth rate. So, we have analyzed the dynamical behaviour of the delayed predator-prey-subsidy model. Numerical computations have been performed using MATLAB to validate all the analytical findings. Numerically, it has been observed that the predator, prey, and subsidy can always exist at a nonzero subsidy input rate. But, at a high subsidy input rate, the prey population cannot persist and the predator population has a huge growth due to the availability of food sources.