The model estimates sulfide availability to the aggregation as a whole by summing the fluxes separately determined for each 2-cm disc composing the rhizosphere. Depth-dependant boundary conditions were set for each disc separately based on the sediment profiles (Figure 5). Diffusional fluxes into each disc were calculated from the shape of the concentration profiles according to the following function [48]: