Bistability in a predator–prey model characterized by the Crowley–Martin functional response: Effects of fear, hunting cooperation, additional foods and nonlinear harvesting

This study focuses on unraveling key factors influencing predator–prey interactions with Crowley–Martin functional response. Specifically, it explores the roles of additional food sources, harvesting practices, hunting cooperation, fear and its carry-over effects. We analyze equilibrium points and their stability properties through rigorous mathematical methods. Numerical illustrations showcase a diverse range of bifurcations including Hopf, saddle–node, and transcritical, providing a comprehensive understanding of the system’s dynamics. We find that the collaboration among predators during hunting induces instability in the system, leading to the emergence of population cycles from a stable state. Further, we place emphasis on investigating the impact of seasonal forcing by introducing time-varying parameters into our model. We reveal the emergence of periodic solutions, higher periodic solutions and chaotic dynamics due to the seasonal variations of the prey’s birth rate and the degree of hunting cooperation. We also emphasize the significance of incorporating different periodicity of seasonally forced parameters, leading to a more precise understanding of predator–prey dynamics.