Fig. 5

PCO ordination plot of the distance between centroids calculated by taxa and aggregating the stations within the levels of the “Season × Area” interaction (WiMar, Winter Margherita di Savoia; WiBar, Winter Barletta; WiLes, Winter Lesina; WiVar, Winter Varano; SuMar, Summer Margherita; SuBar, Summer Barletta; SuLes, Summer Lesina; SuVar, Summer Varano). Vectors indicate the correlation of taxa to the multidimensional space according to Pearson’s correlation coefficient (ρ > 0.75). Taxa are coded as: Acanthocardia tuberculata (Act), Acteon tornatilis (Acto), Amphiura spp. (Ams), Astropecten spp. (Ass), Bolinus brandaris (Bob), Diogenes pugilator (Dip), Dosinia lupinus (Dol), Donax semistriatus (Dos), Echinocardium cordatum (Ecc), Glycera rouxii (Glr), Glycimeris glycimeris (Glg), Hyalinoecia tubicola (Hyt), Holothuria poli (Hop), Liocarcinus depurator (Lid), Mactra glauca (Mag), Maldanidae (Mald), Moerella donacina (Mod), Oestergrenia digitata (Oed), Owenia fusiformis (Owf), Sipunculus nududs (Sin), Spisula subtruncata (Sps), Tritia mutabilis (Trm)