Global Stability and Sensitivity Analysis of a Delayed Nutrient-Plankton Model Under Climate Change and Oceans Deoxygenation

Ocean deoxygenation driven by climate change and nutrient enrichment threatens marine ecosystems. This study analyzes a six-compartment nutrient-plankton-oxygen model incorporating three time-delays. Global asymptotic stability of the coexistence equilibrium is proved using the Lyapunov–Krasovskii method, indicating intrinsic ecosystem resilience. A novel threshold metric, the concentration number is derived via next-generation matrix techniques, signaling significant bloom-induced deoxygenation risk. Sensitivity analysis identifies nutrient-to-phytoplankton conversion ( ), nutrient uptake by phytoplankton ( ), and external nutrient input ( ) as strong positive drivers (+1.0 indices), while washout ( = −0.8) and phytoplankton mortality ( = −0.2) exert mitigating effects. Results show that plankton-mediated oxygen production and greenhouse gas removal provide stabilizing feedback sufficient to offset warming-induced oxygen losses under the parameter values considered. These findings highlight nutrient management and plankton conservation as critical strategies to mitigate climate-driven ocean deoxygenation and preserve marine ecosystem sustainability.

Climate change, Deoxygenation, Hypoxic, Ocean, Plankton.

1. Achimugwu, P. U., Kinyanjui, M. N., & Malonza, D. M. (2023). Mitigation of climate change due to excessive carbon dioxide emission and accumulation: a mathematical model approach. Commun. Math. Biol. Neurosci., 2023, 70. [Google Scholar] [Crossref]

2. Akagha, S. C., Nwankwo, D. I., & Yin, K. (2020). Dynamics of nutrient and phytoplankton in Epe Lagoon, Nigeria: possible causes and consequences of reoccurring cyanobacterial blooms. Applied Water Science, 10(5), 109. [Google Scholar] [Crossref]

3. Carpenter, S. R., Caraco, N. F., Correll, D. L., Howarth, R. W., Sharpley, A. N., &Smith, V. H. (1998). Nonpoint pollution of surface waters with phosphorus and nitrogen. Ecological Applications, 8(3), 559–568. [Google Scholar] [Crossref]

4. Cheng, L., Abraham, J., Trenberth, K. E., Boyer, T., Mann, M. E., & Zhu, J. (2024). New record ocean temperatures and related climate indicators in 2023. Advances in Atmospheric Sciences, 41(6), 1068–1082. [Google Scholar] [Crossref]

5. Chowdhury, P. R., Banerjee, M., & Petrovskii, S. (2024). A two-timescale model of plankton–oxygen dynamics predicts formation of oxygen minimum zones and global anoxia. Journal of Mathematical Biology, 89(1), 8. [Google Scholar] [Crossref]

6. Devlin, M., & Brodie, J. (2023). Nutrients and eutrophication. In Marine pollution– monitoring, management and mitigation (pp. 75–100). Springer. [Google Scholar] [Crossref]

7. Gruber, N., Clement, D., Carter, B. R., Feely, R. A., Van Heuven, S., & Hoppema, M. (2019). The oceanic sink for anthropogenic CO2 from 1994 to 2007. Science, 363(6432), 1193–1199. [Google Scholar] [Crossref]

8. Laffoley, D., & Baxter, J. M. (2019). Ocean deoxygenation: Everyone’s problem-causes, impacts, consequences and solutions. IUCN Gland, Switzerland. [Google Scholar] [Crossref]

9. Macdonald, J. C., & Gulbudak, H. (2023). Forward hysteresis and hopf bifurcation in an npzd model with application to harmful algal blooms. arXiv preprint arXiv:2307.06234. [Google Scholar] [Crossref]

10. Misra, A. K., Singh, R. K., Tiwari, P. K., Khajanchi, S., & Kang, Y. (2020). Dynamics of algae blooming: effects of budget allocation and time delay. Nonlinear Dynamics, 100(2), 1779–1807. [Google Scholar] [Crossref]

11. Mukherjee, N., Ghorai, S., & Banerjee, M. (2019). Detection of turing patterns in a three species food chain model via amplitude equation. Communications in Nonlinear Science and Numerical Simulation, 69, 219–236. [Google Scholar] [Crossref]

12. NOAA. (2024, June 16). Ocean oxygen. https://oceanservice.noaa.gov/facts/ ocean-oxygen.html. (National Oceanic and Atmospheric Administration) [Google Scholar] [Crossref]

13. Ochieng, F. O. (2024). A novel data-driven dynamical model for predicting future climate trends. Modeling Earth Systems and Environment, 10(4), 4599–4610. [Google Scholar] [Crossref]

14. Octavina, C., Maghfirah, M., Lacroix, G., Suciati, N., Fajar, M., Syufi, N. H., & Haditiar, Y. (2025). Spatial analysis of plankton distribution in northern waters of aceh: An indicator of marine environmental quality. In Bio web of conferences (Vol. 156, p. 02005). [Google Scholar] [Crossref]

15. Pachiappan, P., Santhanam, P., Begum, A., & Balaji Prasath, B. (2019). An introduction to plankton. Basic and applied phytoplankton biology, 1–24. [Google Scholar] [Crossref]

16. Prentice, I. C., Farquhar, G., Fasham, M., Goulden, M. L., Heimann, M., & Jaramillo, V. (2001). The carbon cycle and atmospheric carbon dioxide. Climate change 2001: the scientific basis, Intergovernmental panel on climate change. [Google Scholar] [Crossref]

17. Rehim, M., Zhang, Z., & Muhammadhaji, A. (2016). Mathematical analysis of a nutrient– plankton system with delay. SpringerPlus, 5, 1–22. [Google Scholar] [Crossref]

18. Rhein, M., Rintoul, S. R., Aoki, S., Campos, E., Chambers, D., Feely, R. A. & (2013). Observations: ocean. [Google Scholar] [Crossref]

19. Roy, S., Broomhead, D. S., Platt, T., Sathyendranath, S., & Ciavatta, S. (2012). Sequential variations of phytoplankton growth and mortality in an npz model: A remote-sensingbased assessment. Journal of Marine Systems, 92(1), 16–29. [Google Scholar] [Crossref]

20. Schmidtko, S., Stramma, L., & Visbeck, M. (2017). Decline in global oceanic oxygen content during the past five decades. Nature, 542(7641), 335–339. [Google Scholar] [Crossref]

21. Sekerci, Y., & Ozarslan, R. (2020). Respiration effect on plankton–oxygen dynamics in view of non-singular time fractional derivatives. Physica A: Statistical Mechanics and its Applications, 553, 123942. [Google Scholar] [Crossref]

22. Sterpu, M., Rocs，oreanu, C., Efrem, R., & Campbell, S. A. (2023). Stability and bifurcations in a nutrient–phytoplankton–zooplankton model with delayed nutrient recycling with gamma distribution. Mathematics, 11(13), 2911. [Google Scholar] [Crossref]

23. Solomon U. A., Adamu I. I., Benshak S. M. (2026). Mathematical Model for The Dynamics of Nutrient-Phytoplankton System with Delay Under Climate Change and Deoxygination of Oceans. Iconic Research And Engineering Journals, 10(1). doi: https://doi.org/10.64388/IREV10I1-1719581. [Google Scholar] [Crossref]

24. Yang, H. (2022). Global dynamics of a diffusive phytoplankton-zooplankton model with toxic substances effect and delay. Mathematical Biosciences and Engineering, 19(7), 6712–6730. [Google Scholar] [Crossref]

25. Zhuang, K., Jia, G., & Liu, D. (2019). Stability and hopf bifurcation in a three-component planktonic model with spatial diffusion and time delay. Complexity, 2019(1), 4590915. [Google Scholar] [Crossref]

• Assessment of the Role of Artificial Intelligence in Repositioning TVET for Economic Development in Nigeria
• Teachers’ Use of Assure Model Instructional Design on Learners’ Problem Solving Efficacy in Secondary Schools in Bungoma County, Kenya
• “E-Booksan Ang Kaalaman”: Development, Validation, and Utilization of Electronic Book in Academic Performance of Grade 9 Students in Social Studies
• Analyzing EFL University Students’ Academic Speaking Skills Through Self-Recorded Video Presentation
• Major Findings of The Study on Total Quality Management in Teachers’ Education Institutions (TEIs) In Assam – An Evaluative Study