  The coexistence of multiple phytoplankton species despite their reliance on
similar resources is often explained with mean-field models assuming mixed
populations. In reality, observations of phytoplankton indicate spatial
aggregation at all scales, including at the scale of a few individuals. Local
spatial aggregation can hinder competitive exclusion since individuals then
interact mostly with other individuals of their own species, rather than
competitors from different species. To evaluate how microscale spatial
aggregation might explain phytoplankton diversity maintenance, an
individual-based, multispecies representation of cells in a hydrodynamic
environment is required. We formulate a three-dimensional and multispecies
individual-based model of phytoplankton population dynamics at the Kolmogorov
scale. The model is studied through both simulations and the derivation of
spatial moment equations, in connection with point process theory. The spatial
moment equations show a good match between theory and simulations. We
parameterized the model based on phytoplankters' ecological and physical
characteristics, for both large and small phytoplankton. Defining a zone of
potential interactions as the overlap between nutrient depletion volumes, we
show that local species composition-within the range of possible
interactions-depends on the size class of phytoplankton. In large
phytoplankton, individuals are surrounded by cells from other species, while in
small phytoplankton, individuals remain in mostly monospecific clusters.
Spatial structure therefore favours intra- over inter-specific interactions for
small phytoplankton, which likely contributes to coexistence mechanisms. Other
factors behind diversity maintenance must be examined for large phytoplankton.
