This study presents an assessment of the diversity and spatial distribution of benthic macrofauna communities along the Moulay Bousselham lagoon and discusses the environmental factors contributing to observed patterns. In the autumn of 2018, 68 stations were sampled with three replicates per station in subtidal and intertidal areas. Environmental conditions showed that the range of water temperature was from 25.0 °C to 12.3 °C, the salinity varied between 38.7 and 3.7, while the average of pH values fluctuated between 7.3 and 8.0. In vegetated habitats, biomass values of the seagrass Zostera noltei Hornemann ranged between 31.7 gDW/m² and 170.2 gDW/m² while the biomass of the seagrass Ruppia cirrhosa (Petagna) Grande between 54.2 gDW/m² and 84.7 gDW/m². Sediment analyses showed that the lagoon is mainly composed of sandy and silty sediments. We recorded 37,165 individuals of macrofauna distributed in 63 taxa belonging to 50 families, with a mean abundance value of 4582.8 ind/m² and biomass average of 22.2 g/m². Distance-based linear modeling analysis (DISTLM) identified sediment characteristics, water parameters and habitat type (biomass of Z. noltei) as the major environmental drivers influencing macrozoobenthos patterns. Our results clearly revealed that the hydrographic regime (marine and terrestrial freshwater), sediment distribution and characteristics and the type of habitat (vegetated vs. unvegetated substrate) are the key factors determining the species composition and patterns of macrozoobenthos assemblages.

Large-scale research on seagrass-associated benthic fauna is very important for future regional marine conservation. In our study, we investigated spatial and latitudinal variation of benthic macroinvertebrate assemblages associated to Zostera noltei Hornemann, 1832 beds from five semi-enclosed coastal systems (SECSs) ranging from 23°N to 34°N along the Atlantic coast of Morocco. Overall, 17,320 individuals were reported as belonging to 96 taxa. The ecological community descriptors differ significantly at the level of the site. Specific richness showed an inconsistent significant pattern with latitude. The multivariate analyses of the assemblages composition showed 57 of total variation observed in benthic assemblages, while the PERMANOVA analysis confirmed that this variation is significant at the level of the site. According to DistLM results, variations in belowground biomass, and percentage of mud, were the important predictor variables explaining this variation along the large scale of the studied SECS. However, such patterns could be related to other factors such as habitat heterogeneity and regional, biogeographic, and anthropogenic factors. The present study marked the first attempt on broad-scale ecological research of seagrass beds in Morocco and offers baseline data for planning the broad-scale conservation of biodiversity in seagrass beds that remain suffering from multiple human-induced threats such as coastal developments and climate change.