Abstract:
The interplay of phytoplankton, zooplankton, and fish is one of the most important aspects of the aquatic environment. In this paper, we propose to explore the dynamics of a phytoplankton–zooplankton–fish system, with fear-induced birth rate reduction in the middle predator by the top predator and an additional food source for the top predator fish. Phytoplankton–zooplankton and zooplankton–fish interactions are handled using Holling type IV and II responses, respectively. First, we prove the well-posedness of the system, followed by results related to the existence of possible equilibrium points. Conditions under which a different number of interior equilibria exist are also derived here. We also show this existence numerically by varying the intrinsic growth rate of phytoplankton species, which demonstrates the model’s vibrant nature from a mathematical point of view. Furthermore, we performed the local and global stability analysis around the above equilibrium points, and the transversality conditions for the occurrence of Hopf bifurcations and transcritical bifurcations are established. We observe numerically that for low levels of fear, the system behaves chaotically, and as we increase the fear parameter, the solution approaches a stable equilibrium by the route of period-halving. The chaotic behavior of the system at low levels of fear can also be controlled by increasing the quality of additional food. To corroborate our findings, we constructed several phase portraits, time-series graphs, and one- and two-parametric bifurcation diagrams. The computation of the largest Lyapunov exponent and a sketch of Poincaré maps verify the chaotic character of the proposed system. On varying the parametric values, the system exhibits phenomena like multistability and the enrichment paradox, which are the basic qualities of non-linear models. Thus, the current study can also help ecologists to estimate the parameters to study and obtain such important findings related to non-linear PZF systems. Therefore, from a biological and mathematical perspective, the analysis and the corresponding results of this article appear to be rich and interesting.