Emergent ecology in a microscale model of the surface ocean

Microbial processes operate at the microscale, which is not resolved by existing ecosystem models. Here, we present a novel model that simulates a 1 mL three-dimensional cube using a hybrid Lagrangian–Eulerian approach, at ecologically relevant timescales. The model simulates individual microbes, including three phytoplankton size classes with healthy, senescent, and dead lifecycle stages; copiotrophic and oligotrophic heterotrophic bacteria; and dissolved organic matter at 50 µm resolution. Diffusion, shear, sedimentation, chemotaxis, and attachment processes are explicitly resolved. The emerging quantitative representation of the ecosystem shows that (1) copiotrophs grow mostly attached to eukaryotic phytoplankters and get almost all of their carbon from them vs. oligotrophs that grow on exudates and lysates of cyanobacteria; (2) contrasting diel patterns in substrate appearance in the phycosphere vs. ambient water and growth of particle-associated copiotrophs vs. free-living oligotrophs; (3) attached bacteria reduce carbon flux from the phycosphere, lowering chemotactic efficiency toward eukaryotes below that toward cyanobacteria; (4) shear reduces chemotactic efficiency and fitness of the copiotroph; and (5) the main benefit of chemotaxis is to locate attachment partners. These patterns are consistent with available observations. Our study provides insights into the microscale ecology of marine bacteria, and the open-source code is a tool for further research in this area.