Received 26 February 2016; accepted 8 April 2016; published 11 April 2016
Owing to their position at the interface of land and sea, shallow waters and coastal areas are generally highly productive and ecologically important   . 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  -  . Nevertheless, coastal habitats are threatened by anthropogenic stressors, including coastal development and habitat degradation  -  . 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  -  . As a result, many species and habitats are in jeopardy  , 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   .
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  . 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    . As benthic species are relatively long- lived, 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   .
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  . 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   .
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   , clam harvesting  , bait digging  and pollution due to various anthropogenic activities linked to the phosphate industry  . 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)  , which represent a possible threat due to climate change  .
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-year- cycle. 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.
2. Materials and Μethods
2.1. 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 km2 (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  . At low tide, the Kneiss Islands are surrounded by vast mud and sand flats  , 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  .
Figure 1. Map of study area showing location of sampling stations.
2.2. Sampling and Laboratory Procedures
Samples of the macrobenthic macrofauna were collected at low tide using a corer with a sampling area of 0.0225 m2. Five replicates were carried out at each station, with four samples for biological analysis covering a total surface-area of 0.09 m2, and one sample for sediment analysis. Sampling consisted of collecting the topmost 20 cm of the sediment, being carried out at 34 stations in April 2014, in areas corresponding to fully extended Zostera noltei seagrass beds (Figure 1). Moreover, to identify the seasonal changes in macrofauna diversity, 12 stations around the Bessila islet were sampled four times (in October 2013, and in January, April and July 2014).
Sediment was sieved through a 1-mm mesh; the remaining fraction was fixed in 4% buffered formalin and stained with Rose Bengale. In the laboratory, macrofauna was sorted, identified when possible to the species level under the stereomicroscope, and counted. Species number corresponds to the total number of species recorded in 0.09 m2, and abundance to the number of individuals per 1 m2.
2.2.2. Sediment Analysis
The topmost 3-cm sediment layer was also sampled in each replicate for granulometric analysis (median grain- size). Median grain size was determined after sieving the weighed and dried sediment through a wet column of sieves with decreasing apertures (1000, 500, 250, 125 and 63 μm). Sediment samples were then dried at 80˚C to constant weight and ground to a fine powder. The organic matter (OM) content was determined on the powder samples by “loss on ignition” at 450˚C for 4 hrs. Heavy metal contents (Zn, Pb, Cd, Fe and P) were estimated after digesting the powder sample in aqua regia (HCl−, , H2O) at 95˚C, and analysis by Inductively Coupled Plasma Atomic Emission Spectrometry (ICP-AES) and Mass Spectrometry (ICP-MS)  .
2.3. Data Analysis
2.3.1. Biotic Indices
The original data analyzed here consists of a “stations × species” matrix (34 stations × 125 species), which was obtained after removing rare species (29 species). Species were considered rare when they occurred at less than three stations. Abundances were square-root transformed prior to analysis.
Benthic macrofauna communities are defined quantitatively by selected descriptors which reflect the numerical importance of the species. Specific richness S (number of species in each station), abundance A (ind・m−2), Shannon-Weaver’s diversity index H’  and Pielou’s evenness (J’)  are calculated at each station (Table 1). The values of the index of Shannon-Weaver diversity is interpreted according to the limits set by  used in France as part of the Water Framework Directive (WFD) for coastal areas.
Identified species are classified here into trophic groups according to  and later modified by   -  . Three currently available univariate Biotic Indices (BIs) are tested, namely AMBI. BENTIX and BO2A. AMBI (AZTI Marine Biotic Index) based on the study of  , which considers five ecological groups ranging from sensitive species (EGI) to first-order opportunistic species (EGV)  . BENTIX considers only two groups: sensitive (GS) and tolerant species (GT), which correspond to ecological groups I and II, and to ecological groups III to V of the AMBI, respectively   . The BO2A (Benthic Opportunistic Polychaete Amphipod index) is based on the ratio of opportunistic polychaetes (i.e. polychaetes of ecological groups IV and V of the AMBI) and amphipods (except for the genus Jassa)   .
2.3.2. Univariate Analysis
Regarding the seasonal parameters (Table 2), a Shapiro-Wilk normality test and a Bartlett’s test for homogeneity of variance are performed prior to each ANOVA to check whether the assumptions of ANOVA are met and if data transformation is necessary. Then, ANOVAs are performed to assess the influence of the seasons (winter, spring, summer and autumn) on benthic abundance, taxonomic richness, organic matter, temperature, salinity and pH. A Tukey Honestly Significant Difference test is used to determine differences between the four seasons. The same analyses are performed to assess the difference in organic matter between the three sediment types.
2.3.3. Multivariate Analysis
Multivariate analysis is performed to study spatial variability of the macrozoobenthic communities associated with Zostera (Zosterella) noltei seagrass. A square root transformation is applied to the abundance matrix (data for each station are pooled prior to undertaking further analyses), before calculating the Bray-Curtis similarities using the statistical package PRIMER® 6.0  . 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.
Table 1. Main characteristics of sampling stations and macrofauna collected during spring (April) 2014. S: species richness; A: Abundance (ind・m−2), H’: Shannon-Weaner diversity index; J’: Pielou’s evenness index; AMBI values, BO2A values, BENTIX values; T: Temperature (˚C); Sal: salinity; pH; OM: organic matter (g∙kg−1); Sed: sediment; FS: fine sand; MS: medium sand; mud; Long E: longitude east; Lat N: latitude north.
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 1. Hm: heavy metal contents.
3.1. 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, F2,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.
As regards the heavy metals, Zn varies between 181.49 and 638.75 μg∙g−1, Pb between 18.76 and 44.42 μg∙g−1, Fe between 129.87 and 770.10 μg∙g−1, while Cd levels do not exceed 1 μg∙g−1. Phosphate contents fluctuate between 311.42 (IK3) and 792.01 μg∙g−1 (IK11) (Table 2).
3.2. 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 annelid 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 non- selective 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, F3,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, F3,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).
Figure 2. Spatial and seasonal variability of the principal benthic macrofauna parameters in the Kneiss Islands: specific richness (a), abundance (b), Pielou’s evenness (c) and Shannon-Wiener index (d).
3.3. 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.
SIMPER analysis (Table 3) show that GIII shows the most significant similarity between group members (56.90% contribution to total similarity), this group being characterized by numerous representative species especially Perinereis cultifera (Grube, 1840), Gammarus insensibilis (Stock, 1966), Cerastoderma glaucum (Bruguière, 1789) and Arenicola marina (Linnaeus, 1758). SIMPER analysis also indicates that Melinna palmate (Grube, 1870), Loripes lucinalis (Lamarck, 1818), Euclymene oerstedii (Claparède, 1863), and Cirratulus cirratus (O.F. Müller, 1776) contribute significantly to similarities within groups I and II.
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 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.
This study reports the spatial and temporal variability of benthic macrofauna communities over the whole intertidal zone of the Kneiss Islands in the Gulf of Gabès, i.e. an inventory of species associated with Zostera eelgrass, distributed in the intertidal ecosystem, along with the spatial distribution of macrobenthic fauna assemblages, as well as the seasonal variations according to the main environmental changes and anthropogenic pressures, especially traditional human activities.
4.1. 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  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  . 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  -  .
4.2. 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  , the Bay of Tunis, Ghar El-Melh Lagoon, the Tunis lagoon  , the Venice Lagoon and Valle di Gorino (Adriatic Sea)   . The most abundant species are Gammarus insensibilis (Stock, 1966), Scrobicularia plana (Linnaeus, 1758), Cerithium scabridum (Philippi, 1848), Potamides conicus (Blainville, 1829), Loripes lucinalis (Lamarck, 1818), Cirratulus cirratus (O.F. Müller, 1776), Perinereis cultifera (Grube, 1840) and Euclymene aestuarii (Claparède, 1863). There are scant studies on macrozoobenthic assemblages from seagrass-dominated areas in the intertidal zone, particularly in the Mediterranean region    . However, it is well known that seagrasses increase habitat complexity and provide living space and shelter for a diverse animal community  -  . Seagrass bed communities are usually characterized by a larger number of species and higher abundances than adjacent unvegetated sediments   . 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  . A higher abundance of macrofauna has often been reported in vegetated habitats   and different explanatory mechanisms can be proposed; 1) decreased predation efficiency due to high habitat complexity  ; 2) habitat preference of dense seagrass by prey as an escape mechanism from predation  ; 3) stabilization of sediments leading to accumulation of organic matter, allowing increased settlement and growth of infauna  and 4) high content of organic matter, which attracts a certain type of infauna such as deposit-feeding polychaetes  .
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  and in Atlantic lagoons such as Arcachon Bay   . 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  .
4.3. 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  . 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  . 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   -  . 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   . 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  .
4.4. 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  . 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  . 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  .
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  . 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  . It is not clear whether these currents could have some incidence on the distribution of subtidal macrobenthic assemblages.
The authors thank M. Carpenter for the English revision, as well as the fishers of the Kneiss Islands for their help during the sampling and the Laboratory team of GeoResources Materials, Environment and Global Changes (FSS) for their assistance during sediment analysis.
 Afli, A., Ayari, R. and Brahim, M. (2008) Trophic Organization of the Macro-Zoobenthic Assemblages within Coastal Areas Subjected to Anthropogenic Activities. Journal of the Marine Biological Association of the United Kingdom, 88, 663-674. http://dx.doi.org/10.1017/S0025315408001318
 Morrisey, D.J., Skilleter, G.A., Ellis, J.I., Burns, B.R., Kemp, C.E. and Burt, K. (2003) Differences in Benthic Fauna and Sediment among Mangrove (Avicennia marina var. australasica) Stands of Different Ages in New Zealand. Estuarine, Coastal and Shelf Science, 56, 581-592. http://dx.doi.org/10.1016/S0272-7714(02)00208-1
 Blanchet, H., de Montaudouin, X., Lucas, A. and Chardy, P. (2004) Heterogeneity of Macrozoobenthic Assemblages within a Zostera noltii Seagrass Bed: Diversity, Abundance, Biomass and Structuring Factors. Estuarine, Coastal and Shelf Science, 61, 111-123.
 Bostrom, C., O’Brien, K., Roos, C. and Ekebom, J. (2006) Environmental Variables Explaining Structural and Functional Diversity of Seagrass Macrofauna in an Archipelago Landscape. Journal of Experimental Marine Biology and Ecology, 335, 52-73. http://dx.doi.org/10.1016/j.jembe.2006.02.015
 Kennesh, M.J. (2002) Environmental Threats and Environmental Future of Estuaries. Environmental Conservation, 29, 78-107. http://dx.doi.org/10.1017/s0376892902000061
 Kemp, W.M., Boynton, W.R., Adolf, J.E., Boesch, D.F., Bolcourt, W.C., Brush, G., Cornwell, J.C., et al. (2005) Eutrophication of Chesapeake Bay: Historical Trends and Ecological Interactions. Marine Ecology Progress Series, 303, 1-29. http://dx.doi.org/10.3354/meps303001
 Lotze, H.K., Lenihan, H.S., Bourque, B.J., Bradbury, R.H., Cooke, R.G., Kay, M.C., Kidwell, S.M., et al. (2006) Depletion, Degradation, and Recovery Potential of Estuaries and Coastal Seas. Science, 312, 1806-1809.
 Andersen, J.H., Schuter, L. and ?rtebjerg, G. (2006) Coastal Eutrophication: Recent Developments in Definitions and Implications for Monitoring Strategies. Journal of Plankton Research, 28, 621-628.