Abstract
The patterns of species diversity of plankton functional groups (PFGs) remain poorly understood although they matter greatly for marine ecosystem functioning. Here, we use an ensemble of empirical species distribution models for 845 plankton species to estimate the global species richness of three phytoplankton and 11 zooplankton functional groups as a function of objectively selected environmental predictors. The annual mean species richness of all PFGs decreases from the low to the high latitudes, but the steepness and the shape of this decrease vary significantly across PFGs. Pteropods, small copepods (Oithonids and Poecilostomatoids) and Salps have the steepest latitudinal gradients, whereas Amphipods and the three phytoplankton groups have the weakest ones. Temperature, irradiance and nutrient concentration are the first-order control on the latitudinal richness patterns, whilst the environmental conditions associated to upwelling systems, boundary currents and oxygen minimum zones modulate the position of the peaks and troughs in richness. The species richness of all PFGs increases with net primary production but decreases with particles size and the efficiency of the biological carbon pump. Our study puts forward emergent biodiversity-ecosystem functioning relationships and hypotheses about their underlying drivers for future field-based and modelling research.