ABSTRACT: Diatom aggregates contribute significantly to the vertical sinking flux of particulate matter in the ocean. These fragile structures form a specific microhabitat for the aggregated cells, but their internal chemical and physical characteristics remain largely unknown. Studies on the impact of aggregation on the Si cycle led to apparent inconsistency. Despite a lower biogenic silica (bSiO₂) dissolution rate and diffusion of the silicic acid (dSi) being similar in aggregates and in seawater, dSi surprisingly accumulates in aggregates. A reaction-diffusion model helps to clarify this incoherence by reconstructing dSi accumulation measured during batch experiments with aggregated and non-aggregated Skeletonema marinoi and Chaetoceros decipiens. The model calculates the effective bSiO₂ dissolution rate as opposed to the experimental apparent bSiO₂ dissolution rate, which is the results of the effective dissolution of bSiO₂ and transport of dSi out of the aggregate. In the model, dSi transport out of the aggregate is modulated by alternatively considering retention (decrease of the dSi diffusion constant) and adsorption (reversible chemical bonds between dSi and the aggregate matrix) processes. Modelled bSiO₂ dissolution is modulated by the impact of dSi concentration inside aggregates and diatom viability, as enhanced persistence of metabolically active diatoms has been observed in aggregates. Adsorption better explains dSi accumulation within and outside aggregates, raising the possible importance of dSi travelling within aggregates to the deep sea (potentially representing 20% of the total silica flux). The model indicates that bSiO₂ dissolution is effectively decreased in aggregates mainly due to higher diatom viability but also to other parameters discussed herein.