Understanding the spatial structure of microphytobenthos (MPB) on intertidal flats is necessary to gain insight in the benthic community structure and ecosystem processes. The increasing availability of high resolution satellite sensors provides the opportunity to better understand spatial patterns of MPB on various (meter to km) scales. We tested how MPB patch size (indicated by the range derived from a semi-variogram) and degree of patchiness (indicated by the sill) vary as function of seasons, salinity, tidal flat type (muddy fringing versus sandy mid-channel tidal flats) or ecotopes (defined by hydrodynamics, silt content and elevation), in the Westerschelde estuary, the Netherlands. We used Sentinel-2 imagery (2016–2019) with 10 m spatial resolution to derive (omnidirectional) semi-variogram parameters from the NDVI (used as indicator for MPB biomass) and evaluated (seasonality in) patchiness of MPB in the different categories. We demonstrated that MPB patch size (the range) remains constant from winter to summer, while the sill increased from winter to summer. The location of patches on tidal flats was variable throughout the year and shows a remarkable similarity with seasonality in the spatial heterogeneity of the silt content on tidal flats. The patch size and degree of patchiness is higher on relatively sandy mid-channel tidal flats than on relatively silt rich fringing tidal flats. This implies that spatial patterning of MPB biomass on the meso-scale is likely closely linked to abiotic conditions and that spreading processes or grazing activity play a minor role. We observed visually that some areas with a relatively high MPB biomass (‘patches’) remain visible throughout the year, while other patches were only present during a particular season.@en

The prediction of the erosion of mudflats is hampered by inaccurate estimates of the erodibility distribution of the sediment bed. To investigate how erodibility varies in space and what the vertical distribution over the sediment depth is, comprehensive observations of the sediment properties, hydrodynamics and bed-level changes were conducted on an intertidal flat in the Western Scheldt Estuary, the Netherlands. The erosion potential on a mudflat is determined by the critical shear stress for erosion (τe), erosion rate coefficient (M) and local hydrodynamic conditions. A clear difference in hydrodynamic forcing was observed, leading to significant bed level variations at the low water line, where erosion often occurs during very shallow water condition, and a nearly constant bed level at the upper part. The erosion parameters τe and M could be determined over a sediment bed of 12 cm at the low water line. The erosion coefficient M can be considered constant with depth, although there is a large spreading. A clear vertical variation of τe was found: τe increased significantly downward from 0.10 Pa at the sediment surface to 1.13 Pa at 12 cm below the surface. Additionally, there was a strong indication that the presence of diatoms enhanced τe in the upper 2 mm of sediment by five times of the abiotic τe (from 0.09 Pa to 0.46 Pa). These findings lead to the following improvement for predicting morphological changes of tidal mudflats: (1) very shallow conditions should be better simulated, (2) the vertical distribution of τe should be considered. Otherwise, erosion rates can be overestimated, especially during extreme events, because exposure of the deeper well-consolidated layer likely occurs; and (3) an appropriate description of the effect of diatoms should be considered as part of the bottom boundary condition.@en

The prediction of the erosion of mudflats is hampered by inaccurate estimates of the erodibility distribution of the sediment bed. To investigate how erodibility varies in space and what the vertical distribution over the sediment depth is, comprehensive observations of the sediment properties, hydrodynamics and bed-level changes were conducted on an intertidal flat in the Western Scheldt Estuary, the Netherlands. The erosion potential on a mudflat is determined by the critical shear stress for erosion (τe), erosion rate coefficient (M) and local hydrodynamic conditions. A clear difference in hydrodynamic forcing was observed, leading to significant bed level variations at the low water line, where erosion often occurs during very shallow water condition, and a nearly constant bed level at the upper part. The erosion parameters τe and M could be determined over a sediment bed of 12 cm at the low water line. The erosion coefficient M can be considered constant with depth, although there is a large spreading. A clear vertical variation of τe was found: τe increased significantly downward from 0.10 Pa at the sediment surface to 1.13 Pa at 12 cm below the surface. Additionally, there was a strong indication that the presence of diatoms enhanced τe in the upper 2 mm of sediment by five times of the abiotic τe (from 0.09 Pa to 0.46 Pa). These findings lead to the following improvement for predicting morphological changes of tidal mudflats: (1) very shallow conditions should be better simulated, (2) the vertical distribution of τe should be considered. Otherwise, erosion rates can be overestimated, especially during extreme events, because exposure of the deeper well-consolidated layer likely occurs; and (3) an appropriate description of the effect of diatoms should be considered as part of the bottom boundary condition.@en

Microphytobenthos forms an important part of the diet of macrofauna (macrozoobenthos) in many intertidal ecosystems. It is unclear, however, whether the dependence of macrofauna on microphytobenthos varies spatially within and among tidal systems. We aim (1) to assess the spatial variability in the importance of microphytobenthos in the diet of macrofauna (i.e., between and within two tidal basins and as function of elevation), (2) to quantify grazing pressure of the macrofaunal community on different potential food sources (microphytobenthos, phytoplankton and terrestrial organic material) for several sites in two tidal basins and (3) to compare microphytobenthic production and summer/autumn grazing of the total macrofaunal community and grazing pressure per feeding type, with potential microphytobenthic production estimated from rates in early spring, when grazing was low. Using a natural stable isotope approach, we identified microphytobenthos as a more important food source for macrofauna than phytoplankton and terrestrial organic material. Microphytobenthos dependency differed between tidal basins for the genera Bathyporeia (sand digger shrimp), Macoma (Baltic tellin), and Peringia (mudsnail) and for sampled individuals of all genera combined, and did not vary as function of elevation. We showed that macrofaunal grazing on microphytobenthos is quantitatively important and, in some cases, approached microphytobenthic production rates in early spring. No positive relation between microphytobenthic production in early spring and macrofaunal grazing in summer/autumn was observed. This suggests that the studied consumer-resource interactions are coupled on a larger spatial scale (i.e., mesoscale, ≈10 to 100 km), rather than the fine (mm to m) scale.@en