Tay CJ, Mohd MH, Teh SY, Koh HL. Internal phosphorus recycling promotes rich and complex dynamics in an algae-phosphorus model: Implications for eutrophication management.
J Theor Biol 2022;
532:110913. [PMID:
34562459 DOI:
10.1016/j.jtbi.2021.110913]
[Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/16/2021] [Revised: 08/30/2021] [Accepted: 09/15/2021] [Indexed: 10/20/2022]
Abstract
Internal phosphorus recycling in lakes is an important nutrient source that promotes algal growth. Its persistence impedes the effort to improve water quality and thus poses a challenge to the management of eutrophication in lakes, especially in shallow lakes where the occurrence of internal phosphorus recycling is reportedly more common. This paper aims to provide crucial insights on the effects of internal phosphorus recycling on eutrophication dynamics for effective management of lake eutrophication. For this purpose, a mathematical model for lake eutrophication, comprising two compartments of algae and phosphorus, is first formulated for application to a eutrophic tropical lake named Tasik Harapan in Universiti Sains Malaysia. Numerical bifurcation analysis of the model is then performed to assess the combined influences of internal phosphorus recycling, algal mortality and external phosphorus loading on Tasik Harapan eutrophication dynamics. Specifically, co-dimension one bifurcation analysis of algal mortality rate is carried out by means of XPPAUT for various external phosphorus loading rates. The emergence of limit cycle for a certain range of algal mortality rate could be related to the hydra effect (i.e., algal concentration increases in response to greater algal mortality) and the paradox of enrichment (i.e., destabilization of algae in nutrient rich environment). To trace the locus of co-dimension one bifurcation, co-dimension two bifurcation analysis is performed by means of MatCont. The analysis demonstrated that the inclusion of the internal phosphorus recycling term induces rich and complex dynamics of the model. These dynamics include saddle-node bifurcation, cusp, Bogdanov-Takens bifurcation, Generalized Hopf bifurcation and limit point bifurcation of cycles. The results suggest that high internal phosphorus recycling rate promotes bistability and catastrophic shift in a shallow and tropical lake like Tasik Harapan. Hence, the key to effective management of eutrophication in shallow and tropical lakes is the control of internal phosphorus recycling.
Collapse