Spatial and Temporal Structures of the Macrozoobenthos from the Intertidal Zone of the Kneiss Islands ( Central Mediterranean Sea )

The present study analyses the spatio-temporal variability of the macrozoobenthos from the intertidal zone of the Kneiss Islands (Gulf of Gabès, Tunisia). Samples were collected once from 34 stations, while seasonal variations were studied by sampling four times at 12 stations over the period 2013-2014. A total of 159 macrobenthos taxa associated with intertidal Zostera noltei beds are identified from the 34 stations, with a taxonomic dominance of crustaceans (32%), molluscs (29%) and annelids (27%). Abundance varies from 9244 to 36,844 ind·m−2 with a mean value of 14,346 ind·m−2. Analysis of the trophic structure shows that the majority of stations are strongly represented by carnivores (41%), followed by the non-selective deposit feeders (16%). Cluster analysis and multidimensional scaling allow identification of three main benthic assemblages based on species abundance, corresponding to different sediment types and organic matter contents. The seasonal variability in abundance, diversity and community structure is mainly due to spring and summer recruitment. The biotic indices (i.e. AMBI, BO2A and BENTIX) show that the intertidal area of Kneiss Islands exhibits a good ecological status.


Introduction
Owing to their position at the interface of land and sea, shallow waters and coastal areas are generally highly productive and ecologically important [1] [2].However, coastal ecosystems often form a mosaic of interlinked habitats that should not be considered in isolation.Each habitat has its own characteristics and hosts particular benthic assemblages [3]- [5].Nevertheless, coastal habitats are threatened by anthropogenic stressors, including coastal development and habitat degradation [6]- [9].Indeed, many human activities (pollution, tourism, clam harvesting, bait digging, commercial fisheries, eutrophication, sediment discharge, sand extraction and marine transportation) have directly and indirectly affected the biodiversity of these ecosystems on a worldwide scale [10]- [13].As a result, many species and habitats are in jeopardy [14], considering that degradation has modified coastal habitats to such a degree that they no longer fulfill their functional roles as nurseries, in feeding or in reproduction [15] [16].
In Tunisia, very few studies have been carried out on intertidal benthic macrofauna.Nevertheless, the macrobenthos of the Gulf of Gabès has been relatively more studied [17].It is known that the benthic macrofauna is an important component in marine ecosystems, playing a vital role in nutrient cycling, detrital decomposition and as a food source for higher trophic levels.Due to their relatively sessile habit, which implies the inability to avoid unfavourable conditions, macrobenthic species are sensitive indicators of changes in marine environments caused by natural stress or anthropogenic disturbances [13] [18] [19].As benthic species are relatively longlived, they integrate water and sediment quality conditions over time, thus indicating temporal as well as chronic disturbances; consequently, macrofauna can be used as a powerful tool to detect even slight environmental changes [20] [21].
Moreover, the Kneiss Islands, located in the Gulf of Gabès (central Mediterranean Sea), represent a site of international interest in terms of their ornithological diversity (Important Bird area, SPAMI, and Ramsar site), with over 330,000 water birds counted on this wetland in spring [22].This avifaunal biodiversity is closely related to the benthic communities which play a crucial role in the food web either by feeding on detritus, or as food for aquatic birds [23] [24].
Likewise, the marine ecosystems of the Kneiss Islands in the Gulf of Gabès are subject to the loss and degradation of biodiversity caused by bottom trawling [25] [26], clam harvesting [27], bait digging [28] and pollution due to various anthropogenic activities linked to the phosphate industry [29].In addition, discharges to the sea are rendering the Skhira harbour environment hostile for native species and favouring the proliferation of invasive species (both native and non-native) [30], which represent a possible threat due to climate change [31].Mosbahi et al. (2015) presented a first description of the intertidal macrozoobenthic communities in a small part of the Kneiss Islands.By extending this previous research, the present study investigates the spatial and seasonal structure of the macrobenthic communities of a large area around the Kneiss Islands over a 2-yearcycle.The main aims are: 1) to investigate the structural diversity of the benthic macrofauna based on taxonomic, trophic groups and the biodiversity parameters of the benthic macrofauna; 2) to identify the role of the main environmental factors that determine the benthic community structure and functioning; and 3) to provide a reliable assessment of the general ecological status of the macrobenthos.

Study Area
Located in the north-western part of the Gulf of Gabès, between latitudes 34˚10'N -34˚30'N and longitudes 10˚E -10˚30'E, the intertidal zone around the Kneiss Islands extends over an area of 220 km 2 (Figure 1).This area is characterized by the highest tidal range in the western Mediterranean Sea; the tide is semi-diurnal with an amplitude varying from 0.8 to 2.3 m [32].At low tide, the Kneiss Islands are surrounded by vast mud and sand flats [33], making this site the most important area for wintering of migratory waders in the Mediterranean region.
Due to the high biodiversity of benthic communities, the foreshore of the Kneiss Islands represents the largest area in terms of water bird conservation in the Mediterranean zone (Important Bird Area in 2003) and has been designated a Ramsar site since 2007.Moreover, the tidal flats of Kneiss Islands are colonized by seagrass meadows composed of Zostera noltei Hornemann [17].for each station are pooled prior to undertaking further analyses), before calculating the Bray-Curtis similarities using the statistical package PRIMER ® 6.0 [48].A dendrogram is thus created, with group averages expressed in the cluster mode.Then, a non-parametric multi-dimensional scaling (MDS) ordination, using the Bray-Curtis similarity measure, is applied to the abundance matrix (after square-root transformation), with the objective of examining the spatial variability of macrozoobenthic communities.The SIMilarity PERcentages (SIMPER) routine is used to establish which species contribute most to the observed differences in the data.

Abiotic Parameters
The granulometric analyses reveal that there are three main types of sediment: medium sand, fine sand, and mud (Table 1).Moreover, all sampled stations are colonized by Zostera noltei eelgrass beds.
For the physical and chemical parameters, only one measurement was carried out for the whole set of 34 sampled stations in April 2014 (Table 1).The temperatures range from 22.6 to 24.4˚C (mean = 23.3 ± SD 0.46) for temperature, the salinity from 39.1 to 40.5 (mean = 39.6 ± SD 0.32), the pH from 8.05 to 8.95 (mean = 8.7 ± SD 0.18) and organic matter content from 1.32% to 19.12% (mean = 6.16 ± SD 5.78).Temperature, salinity and pH appear homogeneous between the 34 stations.The organic matter shows significant differences between the three sediment types (ANOVA, F 2,31 = 134.9;p < 0.01); the percentage of organic matter is higher in the mud than in the fine or medium sands.No difference is observed between the medium and fine sand (Tukey test).
ANOVA shows that temperature, salinity and pH differ significantly between seasons, with temperatures varying between 12.3˚C in winter (IK8) and 31˚C in summer (station 6), salinity between 38.2 in winter (IK2) and 41.8 in autumn (IK5), pH between 7.84 (IK2) in winter and 8.99 in autumn (IK9) and organic matter between 1.32 (IK4) and 5.28% (IK12).The Tukey test reveals that the temperature and salinity differ significantly between the four seasons, whereas the pH differs only in winter.

Macrofaunal Characteristics
The taxonomic identification of the collected invertebrates produced a list of 159 taxa belonging to six zoological groups.Crustaceans are dominant (32% of the total number of taxa), followed by molluscs (29%) and anne-lid polychaetes (27%).The three other phyla identified (echinoderms, cnidarians and tunicates) account for only 12% of the total number of taxa.The faunistic parameters show a wide range of variability, from 9244 to 36,844 ind•m −2 for abundance (with a mean abundance of 14,346 ± S.D. 3540 ind•m −2 ), 22 to 63 taxa for specific richness, 0.68 to 0.90 for evenness and 3.53 to 5.18 bits./ind.for the Shannon index (Table 1).Trophic structure analysis shows that the majority of stations are strongly represented by carnivores (41%), followed by nonselective deposit feeders NSDF (16%), herbivores (16%), selective deposit feeders SDF (15%) and suspension feeders SF (12%).
For the 12 stations monitored seasonally over one year (Figure 2), the specific richness (a) and the abundance (b) display temporal fluctuations.The number of taxa varies significantly between 16 in winter (IK11 and IK12) and 40 in summer and autumn (IK5 and IK8) (ANOVA, F 3,44 = 14.98; p < 0.001).The number of individuals varies significantly between 3111 ind•m −2 in winter (IK12) and 21,022 ind•m −2 in spring (IK4) (ANOVA, F 3,44 = 48.6;p < 0.001).The Tukey test reveals the absence of any difference in abundance between summer and the autumn; whereas for the taxonomic richness, the winter period differs from the three other seasons with a depletion of species.The evenness (J') and Shannon-Wiener diversity (H') also show fluctuations, with the highest values being observed in summer and autumn.The Shannon-Weaver diversity shows a diversified population across the seasons and is highly diversified during the summer (except for IK1) (Figure 2).Pielou's evenness (J') reflects the absence of species dominance in the population between the four seasons.
AMBI, BO2A and BENTIX indices are consistent with the results classifying all stations as having good ecological status (Table 1).In fact, the majority of stations were strongly dominated by sensitive (EGI) and indifferent species (EGII).ANOVA did not indicate any significant differences for these biotic indices between the all stations (p < 0.05).

Multivariate Analysis
The dendrogram and MDS ordination (Figure 3) allow us to separate the 34 stations into three main groups: the first group (GI) corresponds to the twelve stations sampled around Bessila Islet (IK1 to IK12), mainly made up of fine sand, except for IK2 and IK11 which are composed of medium sand; the second group (GII) corresponding to the twelve stations sampled on Laboua Islet (IK13 to IK24), composed solely of medium sand; and the third and last group (GIII) comprising the other 10 stations (IK25 to IK34) located near the port of embarkation and along a tidal channel characterized by mud sediments.

Characteristics of the Different Faunal Assemblages
In the present study, three distinct macrofaunal assemblages can be recognized in the apparently homogeneous Zostera noltei meadow.The assemblages are characterized by spatial and temporal changes in the population, and their distribution pattern seems to be entirely governed by the physico-chemical and edaphic characteristics of the environment.The faunal assemblages of groups II (12 stations) and III (10 stations) make up the two most representative aspects of the Z. noltei macrobenthic communities.They comprise 65% of stations, being densely vegetated compared to the poorly vegetated stations of group I. Granulometric analyses show that the sediments on the foreshore of the Kneiss Islands are generally fine grained.Seagrass cover and different sediment characteristics indicate that these two factors play an important role in the establishment of different macrozoobenthic assemblages within a Zostera noltei habitat, thus confirming the results found by [4] for the benthic assemblage patterns in Arcachon Bay, France.Furthermore, many ecological studies on benthic fauna note the importance of sediment grain-size in controlling the distribution of the communities.Indeed, certain species colonize very fine sediments in which they build tubes or other protecting biogenic constructions of mud, such as the tubicolous polychaete Melinna palmata Grube, 1870.Other organisms prefer sand or coarse sediment with high porosities [49].Thus, sedimentary texture is a major factor affecting the distribution of benthic communities, and represents a paramount parameter that should not be ignored in the study of benthic fauna [50]- [53].

Comparison with Other Zostera noltei Beds
In this study, we identify 159 taxa, associated with intertidal Zostera noltei seagrass beds, which are unequally distributed among the sampling stations.The benthic community in the Kneiss Islands seems to be characterized by a structure and organization closely similar to ecosystems in other parts of the Mediterranean, dominated mainly by crustaceans, molluscs, polychaetes, cnidarians and echinoderms.This pattern has also been observed in other Mediterranean coastal sites, such as in the Kerkennah Islands [54], the Bay of Tunis, Ghar El-Melh Lagoon, the Tunis lagoon [55], the Venice Lagoon and Valle di Gorino (Adriatic Sea) [56] , 1863).There are scant studies on macrozoobenthic assemblages from seagrass-dominated areas in the intertidal zone, particularly in the Mediterranean region [4] [58] [59].However, it is well known that seagrasses increase habitat complexity and provide living space and shelter for a diverse animal community [60]- [62].Seagrass bed communities are usually characterized by a larger number of species and higher abundances than adjacent unvegetated sediments [63] [64].The effect on these meadows is both structural, since it increases the complexity of the habitat, allowing different species to occupy various ecological niches within a given area, and trophic because it supports epiphytes, a resource for many grazers [65].A higher abundance of macrofauna has often been reported in vegetated habitats [65] [66] and different explanatory mechanisms can be proposed; 1) decreased predation efficiency due to high habitat complexity [67]; 2) habitat preference of dense seagrass by prey as an escape mechanism from predation [5]; 3) stabilization of sediments leading to accumulation of organic matter, allowing increased settlement and growth of infauna [68] and 4) high content of organic matter, which attracts a certain type of infauna such as deposit-feeding polychaetes [62].
The dominant macrozoobenthic assemblages within the intertidal area colonized by the Zostera noltei seagrass beds of the Kneiss Islands (Table 3) are made up of species that are dominant in the Venice lagoon [57] and in Atlantic lagoons such as Arcachon Bay [4] [69].Indeed, a large number of species and high abundances are observed in these study areas, associated with the probable positive effect of seagrass which is commonly related to the structural complexity of the canopy.This favors the successful recruitment and colonization of animals, while the interlacing rhizome layer and roots increase sediment stability [70].

Seasonal Variability of Macrozoobenthic Communities
During the study period, the macrozoobenthic composition was subject to predictable seasonal changes: the main community structure parameters fluctuated according to the typical seasonal cycle of coastal temperate waters [71].The benthic population density and species richness shows seasonal variation in such a way that, at all stations, the maximum is recorded in spring and summer and the minimum during winter.The seasonal variations in mean abundance, diversity and community structure of the macrofauna benthic are mainly caused by recruitment in spring and summer [20].Recruitment is reflected by a general increase in abundance and number of species during spring and summer and a decline during winter.Recruitment is known to be highly variable in space and time and is mainly influenced by numerous factors during the seasonal cycle such as food availability, water temperature, predation and the hydrodynamic regime [20] [72]- [74].Thus, larval availability, larval settlement and the consequent post-recruitment processes such as growth and mortality will have a strong effect on the seasonal variability of the benthic communities [73] [75].The seasonal variability of intertidal macrofauna benthic communities in the Kneiss Islands could be due to the impact of clam harvesting, which is intensively carried out during the winter, leading to a significant decrease in the taxonomic richness and abundance of the surrounding benthic macrofauna [27].

Ecological Status of the Kneiss Islands
Taken together, the calculated values of the various biotic indices show a certain degree of intercorrelation among themselves.Indeed, the majority of these indices classify all stations colonized by intertidal eelgrass meadows as having a good ecological status, with some having a high quality status (Table 1).An analysis of the ecological groups shows that the area is mostly represented by sensitive and indifferent species (EGI and EGII), while there is a lack of first-order opportunistic species (EGV) indicative of polluted environments.An ecological disturbance indicator is useful if it can be easily calculated, while remaining sensitive, anticipatory and integrative across key environmental gradients [76].However, although a single index may provide a good overview of the gradient status of a given benthic environment, it is difficult to choose a particular index that will definitively establish the true status of a specific site [77].For example, the AMBI index sometimes appears to be more appropriate for Atlantic ecosystems and estuarine areas with low biodiversity, few species and high densities, whereas the BENTIX index seems better suited to determining EcoQ in Mediterranean coastal ecosystems with high biodiversity [78].

Conclusions
Long-term studies are essential for understanding the natural variability of species composition, dominance structure and the functional diversity of benthic communities in relation to natural environmental drivers, all of which are affected by climate change [79].This study focuses on the spatial and temporal variability of the taxonomic diversity and structure of the benthic macrofauna in the intertidal zone of the Kneiss Islands, which represents an important area in the Mediterranean as regards tidal range and avian biodiversity.The present study contributes an inventory of 159 macrobenthic species dominated by crustaceans and molluscs, which appears similar in terms of richness to the taxonomic list of other Mediterranean ecosystems.Moreover, we find that sediment type and organic matter content are the fundamental factors accounting for the distribution of assemblages of macrofauna on the tidal flats of the Kneiss Islands.
Finally, it would be of great interest to carry out the same type of study on the subtidal zones of the Kneiss Islands, characterized by the presence of submerged tidal channels.This particular Mediterranean ecosystem is subject to a range of anthropogenic pressures, which lead to a decrease in sediment grain-size from downstream to upstream under the action of tidal currents [80].It is not clear whether these currents could have some incidence on the distribution of subtidal macrobenthic assemblages.

Figure 1 .
Figure 1.Map of study area showing location of sampling stations.

Figure 3 .
Figure 3. Dendrogram and MDS ordination of Bray-Curtis similarities from abundance data (square root transformation) for 34 sampling stations in the intertidal zone of the Kneiss Islands.

Table 2 .
Spatial-temporal variations of main physico-chemical parameters across three seasons (A: autumn; W: winter; S: summer), the parameters for spring are reported in

Table 3 .
MDS formed groups, with indication of similarities within each group (%) and the most representative species (%) contributing to the similarity within the group, determined by SIMPER analysis.