Prochlorococcus, the most abundant photosynthetic organism on Earth, dominates the phytoplankton of the ocean's subtropical gyres, yet statistical and dynamical models long gave opposing predictions of its future abundance. Bian, Cai and Follett (2023) traced the conflict to the fact that the statistical model's predictive power resides almost entirely in sharp, temperature-thresholded range transitions rather than in concentration changes within already-warm waters, and they left open how the hypothesized gyre decline should be tested, what sets the realized as opposed to thermal-maximal range (the observed abundance transition spans more than 15 degrees C across regions), how species distribution models should be validated across spatial-temporal scales, and how statistical and dynamical approaches should be integrated with feedbacks. Ribalet et al. (2025) have since supplied the requested direct test of the decline hypothesis: in situ division rates collapse above roughly 28 degrees C, and an ecosystem model embedding this response projects 17-51% declines of tropical Prochlorococcus production driven by thermal stress rather than nutrient limitation. The surviving open problem is to quantitatively disentangle the bottom-up and top-down controls that set the realized range and its more-than-15-degrees-C transition-temperature spread, to build integrated statistical-dynamical projections of abundance under transient warming, and to formalize multi-scale, pattern-versus-magnitude validation protocols for marine ecological forecasting.

Background

Prochlorococcus is a unicellular marine cyanobacterium roughly 0.5-1 micrometers in diameter and the most abundant photosynthetic organism on Earth. It dominates the phytoplankton of the subtropical gyres — the vast, permanently stratified, nutrient-poor (oligotrophic) central regions of the ocean basins that together cover about half the surface ocean — and contributes a substantial share of global ocean primary production, so how its abundance and geographic range shift as the ocean warms is a first-order question for future marine food webs and carbon cycling. Two modeling traditions gave opposing answers. Species distribution models (SDMs), which are statistical regressions of observed abundance on environmental predictors, were applied by Flombaum et al. (2013), who fitted Prochlorococcus abundance to sea-surface temperature and light (photosynthetically active radiation, PAR) and, combined with warming projections from Earth system models, forecast large global increases by 2100. Dynamical ecosystem models (such as the MIT-Darwin family), which explicitly resolve plankton growth, nutrient cycling, grazing, and viral mortality within an ocean circulation, instead projected decreasing concentrations inside the warming, increasingly stratified gyre interiors, where the supply of deep nutrients declines. Bian, Cai and Follett (2023) diagnosed this conflict. Probing the SDM across spatial-temporal scales — monthly fluctuations at time-series stations, sharp spatial transitions along ocean transects, and basin-scale patterns — they found that temperature and light carry almost no predictive power for fluctuations where Prochlorococcus is already abundant; the SDM's skill comes almost entirely from the large abundance jump between cold, nutrient-rich waters with few cells and warm, oligotrophic waters with many. A two-state model with a single thermal threshold reproduces most of that skill, meaning the SDM captures pattern (where the population is present versus absent) but not magnitude (how many cells). The Flombaum range projection is therefore best read as an estimate of the maximal thermally viable range, while its predicted concentration increases within the warm gyres are unsupported — consistent with the dynamical-model expectation of gyre declines. Bian et al. left several issues open: (i) the hypothesized decline of concentration in the subtropical gyres needed direct experimental and modeling tests, with both its magnitude and its dominant mechanism (nutrient limitation versus grazing or viral lysis) unquantified; (ii) the realized range, as opposed to the thermal maximal range, is set by additional controls — top-down mortality agents (grazing by protists and zooplankton, viral lysis, interactions with heterotrophic bacteria) and bottom-up resources (nutrients, light, seasonality) — and the observed sharp abundance transition occurs across a spread of more than 15 degrees C in sea-surface temperature between regions and seasons, which remains to be quantitatively explained; (iii) the multi-scale validation idea had not been formalized into a systematic protocol for SDMs used in marine ecological forecasting; (iv) integrated statistical-dynamical models including feedbacks between temperature, nutrient cycles, and community structure under transient climate change had not been developed. Subsequently, Ribalet et al. (2025) supplied a direct test of item (i): a decade of continuous-flow-cytometry measurements with the SeaFlow instrument (which counts, sizes, and measures chlorophyll fluorescence of individual phytoplankton cells continuously from a ship's seawater intake) across the tropical and subtropical Pacific showed that in situ division rates of Prochlorococcus rise exponentially with temperature up to about 28 degrees C and then decline sharply. Embedding this thermal response in a global ocean ecosystem model projected 17-51% declines of Prochlorococcus production in the tropical ocean by 2090-2100 under moderate-to-high emissions scenarios and near-collapse in the Western Pacific Warm Pool (the warmest open-ocean region), with direct thermal stress — not reduced nutrient supply — identified as the dominant driver, and hypothetical warm-adapted strains delaying but not preventing the decline. What remains open is the subject of this problem.

Problem Statement

The problem is to quantitatively explain and predict the realized abundance and biogeography of Prochlorococcus in a warming ocean, and to place such prediction on a validated footing. It comprises three linked components that jointly define this single objective: (1) Realized-range controls: identify and quantify, among the candidate bottom-up processes (thermal physiology, nutrient supply, light) and top-down mortality agents (grazing, viral lysis, interactions with heterotrophic bacteria, and programmed cell death — regulated self-destruction of individual cells), the controls that set the realized range of Prochlorococcus as distinct from its thermal maximal range, and give a quantitative, testable account of why the observed sharp transition between high- and low-abundance waters occurs across a spread of more than 15 °C in sea-surface temperature between ocean regions and seasons. (2) Integrated prediction: develop integrated statistical-dynamical models that embed these controls — including feedbacks between temperature, nutrient cycles, and community structure under transient, non-equilibrium warming — and project Prochlorococcus abundance (not only potential range) under specified climate scenarios, consistent with the now-established in-situ thermal response (division rates declining sharply above roughly 28 °C), and determining how realized concentration in the gyre interiors responds given mortality compensation and possible in-situ warm adaptation. (3) Validated forecasting methodology: formalize the multi-scale validation principle of Bian, Cai and Follett (2023) — testing a model's skill separately across spatial-temporal scales, and decomposing its errors into pattern (where populations occur) versus magnitude (how many) components — into a systematic, reproducible validation protocol for species distribution models used in marine ecological forecasting, and demonstrate, against a stated skill or error-separation measure evaluated on a stated held-out benchmark, that it discriminates reliable from unreliable predictions for Prochlorococcus.

References

  1. Bian, Vincent; Cai, Merrick; Follett, Christopher L. (2023). 'Understanding opposing predictions of Prochlorococcus in a changing climate.' Nature Communications 14:1445. doi:10.1038/s41467-023-36928-9. https://doi.org/10.1038/s41467-023-36928-9
  2. Ribalet, François; Dutkiewicz, Stephanie; Monier, Erwan; Armbrust, E. Virginia (2025). 'Future ocean warming may cause large reductions in Prochlorococcus biomass and productivity.' Nature Microbiology 10(10):2441-2453. doi:10.1038/s41564-025-02106-4. https://doi.org/10.1038/s41564-025-02106-4
  3. Flombaum, Pedro; Gallegos, José L.; Gordillo, Rodolfo A.; Rincón, José; Zabala, Lina L.; Jiao, Nianzhi; Karl, David M.; Li, William K. W.; Lomas, Michael W.; Veneziano, Daniele; Vera, Carolina S.; Vrugt, Jasper A.; Martiny, Adam C. (2013). 'Present and future global distributions of the marine cyanobacteria Prochlorococcus and Synechococcus.' Proceedings of the National Academy of Sciences USA 110(24):9824-9829. doi:10.1073/pnas.1307701110. https://doi.org/10.1073/pnas.1307701110
  4. Follett, Christopher L.; Dutkiewicz, Stephanie; Ribalet, François; Zakem, Emily; Caron, David; Armbrust, E. Virginia; Follows, Michael J. (2022). 'Trophic interactions with heterotrophic bacteria limit the range of Prochlorococcus.' Proceedings of the National Academy of Sciences USA 119(2):e2110993118. doi:10.1073/pnas.2110993118. https://doi.org/10.1073/pnas.2110993118
  5. Carlson, Michael C. G.; Ribalet, François; Maidanik, Ilia; Durham, Bryndan P.; Hulata, Yotam; Ferrón, Sara; Weissenbach, Julia; Shamir, Nitzan; Goldin, Svetlana; Baran, Nava; Cael, B. B.; Karl, David M.; White, Angelicque E.; Armbrust, E. Victoria; Lindell, Debbie (2022). 'Viruses affect picocyanobacterial abundance and biogeography in the North Pacific Ocean.' Nature Microbiology 7(4):570-580. doi:10.1038/s41564-022-01088-x. https://doi.org/10.1038/s41564-022-01088-x
  6. Beckett, Stephen J.; Demory, David; Coenen, Ashley R.; Casey, John R.; Dugenne, Mathilde; Follett, Christopher L.; Connell, Paige; Carlson, Michael C. G.; Hu, Sarah K.; Wilson, Samuel T.; Muratore, Daniel; Rodriguez-Gonzalez, Rogelio A.; Peng, Shengyun; Becker, Kevin W.; Mende, Daniel R.; Armbrust, E. Virginia; Caron, David A.; Lindell, Debbie; White, Angelicque E.; Ribalet, François; Weitz, Joshua S. (2024). 'Disentangling top-down drivers of mortality underlying diel population dynamics of Prochlorococcus in the North Pacific Subtropical Gyre.' Nature Communications 15:2105. doi:10.1038/s41467-024-46165-3. https://doi.org/10.1038/s41467-024-46165-3
  7. Frémont, Paul; Beckett, Stephen J.; Demory, David; Carr, Eric; Follett, Christopher L.; Lindell, Debbie; Talmy, David; Dutkiewicz, Stephanie; Weitz, Joshua S. (2026). 'Coexistence of photosynthetic marine microorganisms, viruses and grazers: towards integration in ocean ecosystem models.' Environmental Microbiology 28(4):e70295. doi:10.1111/1462-2920.70295. https://doi.org/10.1111/1462-2920.70295
  8. Pantiukhin, Dmitrii; Verhaegen, Gerlien; Havermans, Charlotte (2024). 'Pan-Arctic distribution modeling reveals climate-change-driven poleward shifts of major gelatinous zooplankton species.' Limnology and Oceanography 69(6):1316-1334. doi:10.1002/lno.12568. https://doi.org/10.1002/lno.12568
  9. Gaudin, Marinna; Eveillard, Damien; Chaffron, Samuel (2024). 'Ecological associations distribution modelling of marine plankton at a global scale.' Philosophical Transactions of the Royal Society B 379(1909):20230169. doi:10.1098/rstb.2023.0169. https://doi.org/10.1098/rstb.2023.0169
  10. Du Toit, Andrea (2025). 'Warming waters, declining Prochlorococcus.' Nature Reviews Microbiology 23(11):682. doi:10.1038/s41579-025-01244-0. https://doi.org/10.1038/s41579-025-01244-0
  11. Wang, P.; Laws, E. A.; Liang, C.; Liu, X.; Chen, J.; Zhong, Y.; Li, C.; Huang, B. (2025). 'Cell death helps to stabilize populations of Prochlorococcus.' Journal of Geophysical Research: Oceans 130(8):e2024JC020962. doi:10.1029/2024JC020962. https://doi.org/10.1029/2024JC020962
  12. Sweeney, Carey P.; Harman, Avery; Kerkar, Anvita U.; Aasjord, Anne E.; Chourrout, Daniel M.; Hiebert, Terra; Sutherland, Kelly; Thompson, Anne W. (2026). 'Prochlorococcus predation by a globally abundant filter feeder.' Marine Ecology Progress Series 788:1-13. doi:10.3354/meps15163. https://doi.org/10.3354/meps15163
  13. Hetharua, Buce; Xu, Min; Sun, Shan; Zhang, Kaidian; Yang, Huidan; Liu, Hongbing; Kao, Shuh-Ji (2025). 'Temperature-driven nitrogen mixotrophy shapes marine cyanobacteria Prochlorococcus and Synechococcus latitudinal distribution pattern.' Communications Earth & Environment 6:149. doi:10.1038/s43247-025-02102-w. https://doi.org/10.1038/s43247-025-02102-w