N2-fixing colonies of cyanobacteria and aggregates of phytoplankton and detritus sinking hundreds of meters per day are instrumental for the ocean's sequestration of CO2 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 O2 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 overestimate O2 flux in these aggregates by as high as 15-fold.