Bifurcation behaviors analysis of a plankton model with multiple delays

A mathematical model describing the dynamics of toxin producing phytoplankton-zooplankton interaction with instantaneous nutrient recycling is proposed. We have explored the dynamics of plankton ecosystem with multiple delays; one due to gestation period in the growth of phytoplankton population and second due to the delay in toxin liberated by TPP. It is established that a sequence of Hopf bifurcations occurs at the interior equilibrium as the delay increases through its critical value. The direction of the Hopf bifurcation and the stability of the bifurcating periodic solutions are determined using the theory of normal form and center manifold. Meanwhile, effect of toxin on the stability of delayed plankton system is also established numerically. Finally, numerical simulations are carried out to support and supplement the analytical findings.