Search

Marine phytoplankton generate half of global primary production, making them essential to ecosystem functioning and biogeochemical cycling. Though phytoplankton are phylogenetically diverse, studies rarely designate unique thermal traits to different taxa, resulting in coarse representations of phytoplankton thermal responses. Here we assessed phytoplankton functional responses to temperature using empirically derived thermal growth rates from four principal contributors to marine productivity: diatoms, dinoflagellates, cyanobacteria, and coccolithophores. Using modeled sea surface temperatures for 1950–1970 and 2080–2100, we explored potential alterations to each group’s growth rates and geographical distribution under a future climate change scenario. Contrary to the commonly applied Eppley formulation, our data suggest phytoplankton functional types may be characterized by different temperature coefficients (Q10), growth maxima thermal dependencies, and thermal ranges which would drive dissimilar responses to each degree of temperature change. These differences, when applied in response to global simulations of future temperature, result in taxon-specific projections of growth and geographic distribution, with low-latitude coccolithophores facing considerable decreases and cyanobacteria substantial increases in growth rates. These results suggest that the singular effect of changing temperature may alter phytoplankton global community structure, owing to the significant variability in thermal response between phytoplankton functional types.

Ocean circulation shapes marine phytoplankton communities by setting environmental conditions and dispersing organisms. In addition, processes acting on the water column (e.g., heat fluxes and mixing) affect the community structure by modulating environmental variables that determine in situ growth and loss rates. Understanding the scales over which phytoplankton communities vary in time and space is key to elucidate the relative contributions of local processes and ocean circulation on phytoplankton distributions. Using a global ocean ecosystem model, we quantify temporal and spatial correlation scales for phytoplankton phenotypes with diverse functional traits and cell sizes. Through this analysis, we address these questions: (1) Over what timescales do perturbations in phytoplankton populations persist? and (2) over what distances are variations in phytoplankton populations synchronous? We find that correlation timescales are short in regions of strong currents, such as the Gulf Stream and Antarctic Circumpolar Current. Conversely, in the subtropical gyres, phytoplankton population anomalies persist for relatively long periods. Spatial correlation length scales are elongated near ocean fronts and narrow boundary currents, reflecting flow paths and frontal patterns. In contrast, we find nearly isotropic spatial correlation fields where current speeds are small, or where mixing acts roughly equally in all directions. Phytoplankton timescales and length scales also vary coherently with phytoplankton body size. In addition to aiding understanding of phytoplankton population dynamics, our results provide global insights to guide the design of biological ocean observing networks and to better interpret data collected at long-term monitoring stations.

The mass sinking of phytoplankton cells following blooms is an important source of carbon to the ocean's interior, with some species contributing more to the flux of particulate organic carbon (POC) than others. During the 2008 North Atlantic Bloom Experiment in the Iceland Basin, we examined plankton community composition from surface waters and from sediment traps at depths down to 750m. Samples collected with neutrally buoyant Lagrangian sediment traps captured a major flux event. Diatoms comprised ≥99% of cell flux into the sediment traps, with vegetative cells and resting spores of the genus Chaetoceros contributing 50-95% of cell flux. Resting spores of one species, identified as Chaetoceros aff. diadema, were dominant, comprising 35-92% of cell flux. The flux of resting spores ranged from 2 to 63mgCm-2day-1 and was significantly correlated with POC flux (p=0.003). Over the course of 10 days, the flux of resting spores increased by 26 fold, suggesting that the cells sank en masse, possibly in aggregates. In contrast, vegetative cells of C. aff. diadema sampled from surface waters during the period preceding the flux event generally comprised <1% of the diatom community and never exceeded 5.2%. Resting spores of C. aff. diadema were rarely observed in surface waters but their concentrations increased with depth (to 200m) below the mixed layer. This increase in resting spore abundance, coupled with increased dissolved silicic acid concentrations at depth, suggest that the morphological changes associated with spore formation may have occurred in the mesopelagic zone, while cells were sinking. The values of variable fluorescence (Fv/Fm) measured on sediment trap material dominated by resting spores were among the highest values measured in the study area at any depth. This, in combination with the rapid germination of resting spores in ship-board incubations, suggests that vegetative cells were not physiologically stressed during spore formation. The degradation