Sudden shifts in marine plankton communities in response to environmental changes are of special concern because of their low predictability and high potential impacts on ocean ecosystems. We explored how anthropogenic climate change influences the spatial extent and frequency of changepoints in plankton populations by comparing the behavior of a plankton community in a coupled Earth system model under pre-industrial, historical 20th century, and projected 21st century forcing. The ocean areas where surface ocean temperature, nutrient concentrations, and different plankton types exhibited changepoints expanded over time. In contrast, regional hotspots where changepoints occur frequently largely disappeared. Heterotrophy and larger organism sizes were associated with more changepoints. In the pre-industrial and 20th century, plankton changepoints were associated with shifts in physical fronts, and more often with changepoints for iron and silicate than for nitrate and phosphate. In the 21st century, climate change disrupts these interannual-variability-driven changepoint patterns. Together, our results suggest anthropogenic climate change may drive less frequent but more widespread changepoints simultaneously affecting several components of pelagic food webs.
Negrete-García, Gabriela, Jessica Y Luo, Matthew C Long, Keith Lindsay, Michael Levy, and Andrew D Barton, December 2022: Plankton energy flows using a global size-structured and trait-based model. Progress in Oceanography, 209, 102898, DOI:10.1016/j.pocean.2022.102898. Abstract
Plankton community models are critical tools for understanding the processes that shape marine plankton communities, how plankton communities impact biogeochemical cycles, and the feedbacks between community structure and function. Here, using the flexible Marine Biogeochemistry Library (MARBL), we present the Size-based Plankton ECological TRAits (MARBL-SPECTRA) model, which is designed to represent a diverse plankton community while remaining computationally tractable. MARBL-SPECTRA is composed of nine phytoplankton and six zooplankton size classes represented using allometric scaling relationships for physiological traits and interactions within multiple functional types. MARBL-SPECTRA is embedded within the global ocean component of the Community Earth System Model (CESM) and simulates large-scale, emergent patterns in phytoplankton growth limitation, plankton phenology, plankton generation time, and trophic transfer efficiency. The model qualitatively reproduces observed global patterns of surface nutrients, chlorophyll biomass, net primary production, and the biogeographies of a range of plankton size classes. In addition, the model simulates how predator:prey dynamics and trophic efficiency vary across gradients in total ecosystem productivity. Shorter food chains that export proportionally more carbon from the surface to the ocean interior occur in productive, eutrophic regions, whereas in oligotrophic regions, the food chains are relatively long and export less organic matter from the surface. The union of functional type modeling with size-resolved, trait-based modeling approaches allows MARBL-SPECTRA to capture both large-scale elemental cycles and the structure of planktonic food webs affecting trophic transfer efficiency.
McGinty, Niall, Andrew D Barton, Nicholas R Record, Zoe V Finkel, David G Johns, Charles A Stock, and Andrew J Irwin, April 2021: Anthropogenic climate change impacts on copepod trait biogeography. Global Change Biology, 27(7), DOI:10.1111/gcb.154991431-1442. Abstract
Copepods are among the most abundant marine metazoans and form a key link between marine primary producers, higher trophic levels, and carbon sequestration pathways. Climate change is projected to change surface ocean temperature by up to 4°C in the North Atlantic with many associated changes including slowing of the overturning circulation, areas of regional freshening, and increased salinity and reductions in nutrients available in the euphotic zone over the next century. These changes will lead to a restructuring of phytoplankton and zooplankton communities with cascading effects throughout the food web. Here we employ observations of copepods, projected changes in ocean climate, and species distribution models to show how climate change may affect the distribution of copepod species in the North Atlantic. On average species move northeast at a rate of 14.1 km decade−1. Species turnover in copepod communities will range from 5% to 75% with the highest turnover rates concentrated in regions of pronounced temperature increase and decrease. The changes in species range vary according to copepod traits with the largest effects found to occur in the cooling, freshening area in the subpolar North Atlantic south of Greenland and in an area of significant warming along the Scotian shelf. Large diapausing copepods (>2.5 mm) which are higher in lipids and a crucial food source for whales, may have an advantage in the cooling waters due to their life‐history strategy that facilitates their survival in the arctic environment. Carnivorous copepods show a basin wide increase in species richness and show significant habitat area increases when their distribution moves poleward while herbivores see significant habitat area losses. The trait‐specific effects highlight the complex consequences of climate change for the marine food web.
Barton, Andrew D., Fernando Gonzalez Taboada, Angus Atkinson, Claire E Widdicombe, and Charles A Stock, August 2020: Integration of temporal environmental variation by the marine plankton community. Marine Ecology Progress Series, 647, DOI:10.3354/meps13432. Abstract
Theory and observations suggest that low frequency variation in marine plankton populations, or red noise, may arise through cumulative integration of white noise atmospheric forcing by the ocean and its amplification within food webs. Here, we revisit evidence for the integration of stochastic atmospheric variations by comparing the power spectra of time series of atmospheric and oceanographic conditions to the population dynamics of 150 plankton taxa at Station L4 in the Western English Channel. The power spectra of oceanographic conditions (sea surface temperature, surface nitrate) are redder than those of atmospheric forcing (surface wind stress, net heat fluxes) at Station L4. However, plankton populations have power spectral slopes across trophic levels and body sizes that are redder than atmospheric forcing but whiter than oceanographic conditions. While zooplankton have redder spectral slopes than phytoplankton, there is no significant relationship between power spectral slope and body size or generation length. Using a predator-prey model, we show that the whitening of plankton time series relative to oceanographic conditions arises from noisy plankton bloom dynamics in this strongly seasonal system. The model indicates that, for typical predator-prey interactions, where the predator is on average 10 times longer than the prey, grazing leads to a modest reddening of phytoplankton variability by their larger and longer lived zooplankton consumers. Our findings suggest that, beyond extrinsic forcing by the environment, predator-prey interactions play a role in influencing the power spectra of time series of plankton populations.
Seasonal to interannual predictions of ecosystem dynamics have the potential to improve the management of living marine resources. Prediction of oceanic net primary production (NPP), the foundation of marine food webs and the biological carbon pump, is particularly promising, with recent analysis suggesting that ecosystem feedback processes may lead to higher predictability of NPP at interannual scales than for physical variables like sea surface temperature (SST). Here, we assessed the potential predictability of oceanic NPP and SST across seasonal to interannual lead times using reduced dimension, linear dynamical spatio-temporal models (rDSTM). This approach combines empirical orthogonal function (EOF) analysis with vector autoregressive (VAR) modeling to simplify the analysis of spatio-temporal data. The rDSTMs were fitted to monthly NPP and SST anomalies derived from 20 years of remote sensing data (1997-2017), considering two alternative algorithms commonly used to estimate NPP (VGPM and Eppley-VGPM) and optimally analyzed SST fields (AVHRR OISST). The local decay of anomalies provided high predictability up to three months, and subsequent interactions with remote forcing significantly extended predictability beyond the initial anomaly decay. Indeed, interactions among spatial modes associated with the propagation of major climate modes, particularly the El Niño-Southern Oscillation (ENSO), extended the predictability horizon above one year in some regions. Patterns of enhanced NPP predictability matched the location of oligotrophic gyres and transition regions between ocean biomes, where fluctuations in biome boundaries generate large biogeochemical perturbations that leave lasting imprints on NPP. In these areas, the predictability horizon for NPP was longer than for SST, although SST was more predictable over large areas of the equatorial and northeast Pacific. Our results support the potential for extending seasonal to interannual physical climate predictions to predict ocean productivity.
The measured concentration of chlorophyll a in the surface ocean spans four orders of magnitude, from ∼0.01 mg m-3 in the oligotrophic gyres to >10 mg m-3 in coastal zones. Productive regions encompass only a small fraction of the global ocean area yet they contribute disproportionately to marine resources and biogeochemical processes, such as fish catch and coastal hypoxia. These regions and/or the full observed range of chlorophyll concentration, however, are often poorly represented in global earth system models (ESMs) used to project climate change impacts on marine ecosystems. Furthermore, recent high resolution (∼10 km) global earth system simulations suggest that this shortfall is not solely due to coarse resolution (∼100 km) of most global ESMs. By integrating a global biogeochemical model that includes two phytoplankton size classes (typical of many ESMs) into a regional simulation of the California Current System (CCS) we test the hypothesis that a combination of higher spatial resolution and enhanced resolution of phytoplankton size classes and grazer linkages may enable global ESMs to better capture the full range of observed chlorophyll. The CCS is notable for encompassing both oligotrophic (<0.1 mg m-3) and productive (>10 mg m-3) endpoints of the global chlorophyll distribution. As was the case for global high-resolution simulations, the regional high-resolution implementation with two size classes fails to capture the productive endpoint. The addition of a third phytoplankton size class representing a chain-forming coastal diatom enables such models to capture the full range of chlorophyll concentration along a nutrient supply gradient, from highly productive coastal upwelling systems to oligotrophic gyres. Weaker ‘top-down’ control on coastal diatoms results in stronger trophic decoupling and increased phytoplankton biomass, following the introduction of new nutrients to the photic zone. The enhanced representation of near-shore chlorophyll maxima allows the model to better capture coastal hypoxia along the continental shelf of the North American west coast and may improve the representation of living marine resources.
Barton, Andrew D., Andrew J Irwin, Zoe V Finkel, and Charles A Stock, March 2016: Anthropogenic climate change drives shift and shuffle in North Atlantic phytoplankton communities. Proceedings of the National Academy of Sciences, 113(11), DOI:10.1073/pnas.1519080113. Abstract
Anthropogenic climate change has shifted the biogeography and phenology of many terrestrial and marine species. Marine phytoplankton communities appear sensitive to climate change, yet understanding of how individual species may respond to anthropogenic climate change remains limited. Here, using historical environmental and phytoplankton observations, we characterize the realized ecological niches for 87 North Atlantic diatom and dinoflagellate taxa and project changes in species biogeography between mean historical (1951–2000) and future (2051–2100) ocean conditions. We find that the central positions of the core range of 74% of taxa shift poleward at a median rate of 12.9 km per decade (km⋅dec−1), and 90% of taxa shift eastward at a median rate of 42.7 km⋅dec−1. The poleward shift is faster than previously reported for marine taxa, and the predominance of longitudinal shifts is driven by dynamic changes in multiple environmental drivers, rather than a strictly poleward, temperature-driven redistribution of ocean habitats. A century of climate change significantly shuffles community composition by a basin-wide median value of 16%, compared with seasonal variations of 46%. The North Atlantic phytoplankton community appears poised for marked shift and shuffle, which may have broad effects on food webs and biogeochemical cycles.