Multi-nutrient, multi-group model of present and future oceanic phytoplankton communities

• Main Authors: E. Litchman, C. A. Klausmeier, J. R. Miller, O. M. Schofield, P. G. Falkowski
• Format: Article in Journal/Newspaper
• Language: English
• Published: Copernicus Publications 2006
• Subjects: Ecology; QH540-549.5; Life; QH501-531; Geology; QE1-996.5; Pacific; North Atlantic; Subarctic
• Online Access: https://doaj.org/article/92c2b7f5b7db4deabaa23116dda4defb

Phytoplankton community composition profoundly affects patterns of nutrient cycling and the dynamics of marine food webs; therefore predicting present and future phytoplankton community structure is crucial to understand how ocean ecosystems respond to physical forcing and nutrient limitations. We develop a mechanistic model of phytoplankton communities that includes multiple taxonomic groups (diatoms, coccolithophores and prasinophytes), nutrients (nitrate, ammonium, phosphate, silicate and iron), light, and a generalist zooplankton grazer. Each taxonomic group was parameterized based on an extensive literature survey. We test the model at two contrasting sites in the modern ocean, the North Atlantic (North Atlantic Bloom Experiment, NABE) and subarctic North Pacific (ocean station Papa, OSP). The model successfully predicts general patterns of community composition and succession at both sites: In the North Atlantic, the model predicts a spring diatom bloom, followed by coccolithophore and prasinophyte blooms later in the season. In the North Pacific, the model reproduces the low chlorophyll community dominated by prasinophytes and coccolithophores, with low total biomass variability and high nutrient concentrations throughout the year. Sensitivity analysis revealed that the identity of the most sensitive parameters and the range of acceptable parameters differed between the two sites. We then use the model to predict community reorganization under different global change scenarios: a later onset and extended duration of stratification, with shallower mixed layer depths due to increased greenhouse gas concentrations; increase in deep water nitrogen; decrease in deep water phosphorus and increase or decrease in iron concentration. To estimate uncertainty in our predictions, we used a Monte Carlo sampling of the parameter space where future scenarios were run using parameter combinations that produced acceptable modern day outcomes and the robustness of the predictions was determined. Change in the onset and ...