Abstract
In this paper, we study the dynamics of a time fractional order toxic-phytoplankton-phytoplankton-zooplankton system (TPPZS). Routh-Hurwitz criteria is applied to discuss the stability analysis of biologically feasible equilibrium points for the given system in terms of reproduction numbers (associated with ecological as well as disease phenomena). Local stability properties of the toxic-phytoplankton-free equilibrium are also investigated. Numerical simulations are performed for a hypothetical set of parameter values in support of our analytic results. (C) 2015 Elsevier B.V. All rights reserved.