Mathematical analysis of plankton population dynamics


  • Fatin Nadiah Yussof
  • Normah Maan
  • Nadzri Reba



Stability Analysis, Harmful Algal Blooms, Toxin Producing Phytoplankton, Nutrient Limitation


Harmful algal blooms (HABs) event that causes enormous economic loss and health effect raises concerns among environmentalists. In this paper, a mathematical model of interaction between nutrient, toxin-producing phytoplankton (TPP), non-toxic phytoplankton (NTP), zooplankton, and toxic chemicals is proposed to study on how the process of these HABs occurred. The model of interaction is represented by Ordinary Differential Equations (ODEs) and stability analysis of the model is conducted. Several conditions for the system to be stable around trivial and interior equilibrium point are obtained. From the analysis, it is observed that under nutrient limitation, the amounts of toxic chemicals secreted out by the TPP are increased. As a result, NTP population and zooplankton population are affected by the situation. If this situation is prolonged, this will result in the extinction of both populations. Overall, this study shows that TPP release more toxic chemicals when the nutrient is limited and gives a better understanding on the occurrence of HABs event.



Windust, A. J., Wright, J. L. C., and McLachlan, J. L. 1996. The effects of the diarhetic shellfish poisoning toxins, okadaic acid and dinophysistoxin-1, on the growth of microalgae. Marine biology. 126(1): 19-25.

Hulot, F. D., and Huisman, J. 2004. Allelopathic interactions between phytoplankton species: the roles of heterotrophic bacteria and mixing intensity. Limnology and Oceanography. 49(4part2): 1424-1434.

Lim, H. C., Leaw, C. P., Tan, T. H., Kon, N. F., Yek, L. H., Hii, K. S., Teng, S. T., Razali, R. M., Usup, G., Iwataki, M., and Lim, P. T. 2014. A bloom of Karlodinium australe (Gymnodiniales, Dinophyceae) associated with mass mortality of cage-cultured fishes in West Johor Strait, Malaysia. Harmful Algae. 40: 51-62.

Chakraborty, S., Roy, S., and Chattopadhyay, J. 2008. Nutrient-limited toxin production and the dynamics of two phytoplankton in culture media: A mathematical model. Ecological modelling. 213(2): 191-201.

Roy, S., and Chattopadhyay, J. 2007. Toxin-allelopathy among phytoplankton species prevents competitive exclusion. Journal of Biological Systems. 15(1): 73-93.

Anderson, D. M., Kulis, D. M., Sullivan, J. J., and Hall, S. 1990. Analysis of a dengue fever transmission model. Toxicon. 28(8): 885-893.

Boyer, G. L., Sullivan, J. J., Andersen, R. J., Harrison, P. J., and Taylor, F. J. R. 1987. Effects of nutrient limitation on toxin production and composition in the marine dinoflagellate Protogonyaulax tamarensis. Marine Biology. 96(1): 123-128.

Hall, S. 1982. Toxins and toxicity of Protogonyaulax from the northeast Pacific. PhD thesis. University of Alaska Fairbanks.

McLachlan, J. L., Marr, J. C., ad Conlon-Keily, A., and Adamson, A. 1994. Effects of nitrogen concentration and cold temperature on DSP-toxin concentrations in the dinoflagellate Prorocentrum lima (prorocentrales, dinophyceae). Natural toxins. 2(5): 263-270.

Sohet, K., Pereira, A., and Braekman, J. C., and Houvenaghel, G. 1995. Growth and toxicity of Prorocentrum lima (Ehrenberg) dodge in different culture media: Effect of humic acids and organic phosphorus. in: Lassus, P (Ed.). Harmful marine algal blooms: Proceedings of the Sixth International Conference on toxic marine phytoplankton = Proliférations d'algues marines nuisibles: Sixiéme Conférence Internationale sur le phytoplancton toxique, October 18-22, 1993, Nantes, France. pp. 669-674.

Johansson, N., and Graneli, E. 1999. Influence of different nutrient conditions on cell density, chemical composition and toxicity of Prymnesium parvum (Haptophyta) in semi-continuous cultures. Journal of Experimental Marine Biology and Ecology. 239(2): 243-258.

Rehim, M., Zhang, Z., and Muhammadhaji, A. 2016. Mathematical analysis of a nutrient-plankton system with delay. SpringerPlus. 5(1): 1055.

Chattopadhayay, J., Sarkar, R. R., and Mandal, S. 2002. Toxin-producing plankton may act as a biological control for planktonic blooms field study and mathematical modelling. Journal of Theoretical Biology. 215(3): 333-344.

Das, K., and Ray, S. 2008. Effect of delay on nutrient cycling in phytoplankton-zooplankton interactions in estuarine system. Ecological Modelling. 215(1-3): 69-76.

Dugdale, R. C. J. 1967. Nutrient limitation in the sea: Dynamics, identification, and significance. Limnology and Oceanography. 12(4): 685-695.

Franks, P. J. S. 2002. NPZ models of plankton dynamics: Their construction, coupling to physics and application. Journal of Oceanography. 58(2): 379-387.

Teen, L. P., Pin, L. C., Gires, U. 2012. Harmful algal blooms in Malaysian waters. Sains Malaysiana. 41(12):1509-1515.

Fan, A., Han, P., and Wang, K. 2013. Global dynamics of a nutrient-plankton system in the water ecosystem. Applied Mathematics and Computation. 219(15): 8269-8276.

Chakraborty, K., and Dasb, K. 2015. Modeling and analysis of a two-zooplankton one phytoplankton system in the presence of toxicity. Applied Mathematical Modelling. 39(3-4): 1241-1265.

Khare, S., Dhar, J., and Misra, O. P. 2010. Role of toxin-producing phytoplankton on a plankton ecosystem. Nonlinear Analysis: Hybrid Systems. 4(3): 496-502.

Chatterjee, A., Pal, S., and Chatterjee, S. 2011. Bottom up and top down effect on toxin-producing phytoplankton and its consequence on the formation of plankton bloom. Applied Mathematics and Computation. 218(7): 3387-3398.