“…This implicit approach has been successfully applied to systems with relatively simple and stable benthic communities (e.g., Soetaert et al, 1996Soetaert et al, , 2000Paraska et al, 2014), but has been insufficient to explain nutrient fluxes in some shallow, near-shore areas (Murray and Parslow, 1997;Holstein and Wirtz, 2009), suggesting that the temporal variability in fauna and flora and associated processes typical for many shallow coastal areas needs to be accounted for by dynamic state variables. However, linking biodiffusion (D b ) to biomass dynamics of functional groups (Ebenhöh et al, 1995;Sohma et al, 2001Sohma et al, , 2004Sohma et al, , 2008Sohma et al, , 2018Vichi, 2002;Fulton et al, 2004a,b;Audzijonyte et al, 2017) may not be enough to account for the great variability in reworking rates measured between coastal areas (Boudreau, 1994;Tromp et al, 1995;Middelburg et al, 1997), and even between replicates from the same site (Morys et al, 2016).…”