A New Mathematical Model to Explore Microbial Processes and their Constraints in Phytoplankton Colonies and Sinking Marine Aggregates

N₂-fixing colonies of cyanobacteria and aggregates of phytoplankton and detritus sinking hundreds of meters per day are instrumental for the ocean’s sequestration of CO₂ from the atmosphere. Understanding of small-scale microbial processes associated with phytoplankton colonies and aggregates is therefore crucial for understanding large-scale biogeochemical processes in the ocean. Phytoplankton colonies and sinking aggregates are characterized by steep concentration gradients of gases and nutrients in their interior. Here, we present a mechanistic mathematical model designed to perform modeling of small-scale fluxes and evaluate the physical, chemical, and biological constraints of processes that co-occur in phytoplankton colonies and sinking porous aggregates. The model accurately reproduced empirical measurements of O₂ concentrations and fluxes measured in sinking aggregates. Common theoretical assumptions of either constant concentration or constant flux over the entire surface did not apply to sinking aggregates. Consequently, previous theoretical models over estimate O₂ flux in these aggregates by as high as 15-fold.