Plankton occupy a vital place in the marine ecosystem due to their essential role. However small or microscopic, their absence can bring the entire life process to a standstill. In this work, we have proposed a prey–predator ecological model consisting of phytoplankton, zooplankton, and fish, incorporating the cannibalistic nature of zooplankton harvesting the fish population. Due to differences in their feeding habits, zooplankton are divided into two sub-classes: herbivorous and carnivorous. The dynamic behavior of the model is examined for each of the possible steady states. The stability criteria of the model have been analyzed from both local and global perspectives. Hopf bifurcation analysis has been accomplished with the growth rate of carnivorous zooplankton using cannibalism as a bifurcation parameter. To characterize the optimal control, we have used Pontryagin’s maximum principle. Subsequently, the optimal system has been derived and solved numerically using an iterative method with Runge–Kutta fourth-order scheme. Finally, to facilitate the interpretation of our mathematical results, we have proceeded to investigate it using numerical simulations.

It is frequently observed that adult members of prey species sometimes use their predation mechanism on juvenile members of predator species. Ecological literature describes this phenomenon as prey–predator role reversal dynamics.Numerous authors have observed and described the biological development behind this feeding behaviour. However, the dynamics of this role reversal have hardly been illustrated in the literature in a precise way. In this regard, we formulated an ecological model using the standard prey–predator interactions, allowing for a reverse feeding mechanism. The mathematical model consisted of a three-species food-web structure comprising the common prey, intermediate predator, and top predator. Note that a role-reversal mechanism was observed between the intermediate and top predators based on the scarcity of the prey population. However, we observed the most critical parameters had a significant effect on this reverse feeding behaviour. The bifurcation analysis is the primary criterion for this identification. The proposed deterministic model is then extended to its stochastic analogue by allowing for environmental influences on the tri-trophic food web structure. The conditional moment approach is applied to obtain the equilibrium distribution of populations and their conditional moments in the system. The stochastic setup analysis also supports the stability of this food chain structure, with some restricted conditions. Finally, to facilitate the interpretation of our mathematical results, we investigated it using numerical simulations.