Search tips
Search criteria 


Logo of plosonePLoS OneView this ArticleSubmit to PLoSGet E-mail AlertsContact UsPublic Library of Science (PLoS)
PLoS One. 2017; 12(9): e0184241.
Published online 2017 September 5. doi:  10.1371/journal.pone.0184241
PMCID: PMC5584924

Comparing demersal megafaunal species diversity along the depth gradient within the South Aegean and Cretan Seas (Eastern Mediterranean)

Panagiota Peristeraki, Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Supervision, Writing – original draft, Writing – review & editing,1,2,* George Tserpes, Formal analysis, Methodology, Visualization, Writing – review & editing,1 Nikolaos Lampadariou, Methodology, Supervision, Writing – review & editing,3 and Kostantinos I. Stergiou, Supervision, Writing – review & editing4,5
Carlo Nike Bianchi, Editor


Knowledge on biodiversity patterns of demersal megafaunal species in the Mediterranean and particularly in its eastern basin is still very scarce. In the present study, fine-scale diversity patterns in relation to depth were analyzed for three major megafaunal groups (fish, cephalopods and crustaceans) in three subareas of the eastern Mediterranean (Crete, Cyclades and Dodecanese islands). The analysis was based on data from the Mediterranean International Trawl Survey conducted during 2005–2014 and the relationship between depth and two different diversity measures (species richness and Shannon-Weaver) was examined using Generalized Additive Modeling (GAM) techniques. Species richness of fish decreased with depth in two of the three subareas (Cyclades, Dodecanese), while the opposite was true for crustaceans in all subareas. Cephalopods had higher species richness at intermediate depths, near the shelf break. Significant differences among subareas were found, with Crete showing a distinct species richness-depth pattern, which was more obvious for fish and cephalopods. The differences among subareas were also highlighted based on the occurrence of alien species of Indo-Pacific origin, which were more frequent in Crete. Our results suggested that the importance of depth-related factors in structuring communities was higher for cephalopods and less important for fish, and that Crete showed a distinct diversity-depth relationship, a fact that can be attributed to its specific geographical and oceanographic characteristics. These results support the current GFCM/FAO’s characterization of Crete as a unique geographic subarea. The findings of the study contribute to understanding the causes of underlying diversity patterns and would assist various environmental management actions, particularly those related to the establishment of marine-protected areas.


Species diversity has important implications for the functioning of ecosystems. Understanding the causes of underlying diversity patterns, as well as their interaction with environmental factors, is important for ecosystem conservation purposes.

In the past, numerous publications on aquatic sciences, in different geographic areas, have studied faunal zonation, including the distribution of species along environmental gradients. These studies suggest that physical factors, such as depth, temperature, pressure, hydrographic conditions, oxygen content, sediment type, water mass structure and topography, as well as biological interactions, are possible causes of spatial community changes in marine environments [17]. Fishing pressure has been also suggested as a structuring factor for marine communities [812]. Among abiotic parameters, depth is characterised as the most important discriminant factor for faunal communities e.g. [9, 1318, 19].

Various biodiversity aspects have been studied in the Mediterranean Sea e.g. [16, 1928] but such studies are rather scarce for the eastern Mediterranean. Previous studies conducted on a Mediterranean-wide basis have provided some broad information for some eastern Mediterranean areas, such as the Aegean and Cretan Seas, averaging environmental covariates (e.g. depth, temperature) across larger geographic areas [2328]. However, such a coarse-scale approach can mask local species-environment interactions that could provide important information about ecological processes influencing species abundance and distribution [29, 30]. Moreover, fine-scale community studies in the Aegean and Cretan Seas are rather scarce and most of them did not consider any specific community attributes. The few existing studies deal with depth-related distributional patterns of demersal megafaunal species, mostly focusing on the identification of “key-species” by depth zone [14, 3134]. Thus, the lack of fine-scale community studies of demersal megafaunal communities in the eastern Mediterranean is a significant gap in our attempt to better understand species distribution and habitat selection patterns [35].

In the present study, a fine-scale spatial analysis of the changes in demersal megafaunal species diversity in relation to depth is presented. The analysis considers three large faunal groups, fish, cephalopods and crustaceans, and is based on data collected during the Mediterranean International Trawl Survey (MEDITS) carried out between 2005 and 2014 in the South Aegean and Cretan Seas.

Materials and methods

Ethic statement

Data were collected in the framework of the MEDITS project and refer to the South Aegean Sea (eastern Mediterranean). The sampling procedures followed a standardized protocol (for details see [36]), approved by international authorities (EU/DG Mare, FAO/GFCM). The specific surveys (conducted during 2005–2014) comprise a task of the Data Collection Framework for the Greek Fisheries and are accomplished with the permission of the South Aegean and Cretan Prefectures, which are responsible for the protection of the environment in the corresponding areas. The first two authors coordinate and consistently participate in the MEDITS surveys.

Study area

The South Aegean Sea is a geographically and hydrologically complex area, with deep basins interrupted by shallower island complexes. The general large-scale water circulation pattern of the Aegean is anticlockwise. The more saline eastern Mediterranean waters flow northwards in the eastern part of the South Aegean, while the less-saline Black Sea waters flow southwards in the western part [37, 38].

The MEDITS sampling scheme in the South Aegean is rather discontinuous as the complex geomorphology of the area has resulted in large deep-water regions (>800 m) that are not included in the sampling protocol (Fig 1).

Fig 1
Distribution of MEDITS sampling stations (red circles) in the South Aegean and Cretan Seas.

Taking into account the existing oceanographic characteristics and the survey sampling scheme, three subareas were considered:

  1. Crete is an isolated island, one of the biggest in the Mediterranean, with a narrow and relatively steep continental shelf (gradient 1.5°) followed by a steep slope (2°–4°). It is surrounded by the deep Ionian, Libyan and Cretan Seas; the latter is the largest in volume and the deepest (2500 m) basin of the Aegean Sea, separating Crete from the Cycladic and Dodecanese island complexes [38]. The Cretan Sea is one of the most oligotrophic areas of the Mediterranean, with very low primary production [39] and the presence of a strong pelagic microbial loop [40], which reduces particle flux to the sediment [41]. For management purposes, Crete is considered by the General Fisheries Commission for the Mediterranean (GFCM) as a separate Geographical Subarea (GSA).
  2. Cyclades is a complex of islands located in the middle of the Aegean Sea, comprising the Cyclades Plateau. The latter is a complex marginal platform, shallower than 250 m, with numerous islands. It separates the central from the South Aegean, and the deeper parts of the Plateau (150–250 m) are covered by sediments ranging from calcareous muddy sand to sandy mud [42]. The west edge of the Cyclades complex neighbors the Greek mainland and is influenced by the less-saline Black Sea waters [37, 43].
  3. Dodecanese is an island complex, dispersed in the eastern part of Aegean Sea, along the coast of Minor Asia, with a high variability of depths [38] because of the wide continental shelf at the northern part, interrupted by deep waters at the southern part of the area. Most of the Dodecanese islands are situated near the large shallow bays of the Turkish coast. They are mainly influenced from the more saline eastern Mediterranean waters that flow northwards in the eastern part of Aegean. Dodecanese is separated from the Cyclades complex by deep seas, and is considered the gateway for Lessepsian species, many of which establish populations here before their westward expansion [44, 45]. Such species have been considered an important driver of diversity changes in the Mediterranean [46].

Though no specific estimates for the total fishing pressure in each subarea are available, fishing pressure can be considered higher in Cyclades, followed by Dodecanese and lower in Crete, according to the total number of registered commercial fishing vessels in each area. According to the classification of fishing areas by Sylaios et al., [47], Cyclades is characterized as a fishing area with a moderate to high fishery production, while Dodecanese and Crete are characterized as fishing areas with a low to moderate fishery production.


The study utilized data obtained from the MEDITS experimental surveys conducted in the South Aegean and Cretan Seas. The MEDITS survey is carried out annually in various areas of the Mediterranean Sea following a standardized protocol (see [36] for further details) and its primary goal is to monitor changes regarding the abundance and structure of the demersal megafauna community. The survey design is based on a depth-stratified random sampling scheme considering five depth strata: 10–50 m, 50–100 m, 100–200 m, 200–500 m, and 500–800 m. Sampling in the South Aegean and Cretan Seas includes a total of 61 fixed sampling stations, which are distributed from 20 to 800 m according to the surface of the depth strata (Fig 1).

For each station (haul), all fish, cephalopod and crustacean specimens were identified to the species level according to FAO species identification keys. In a few cases, specimens were identified only to the genus or family level due to difficulties in identification on board (fish: some Gobius and Myctophidae specimens; crustaceans: individuals of the genera Munida, Pagurus and Inacus and the family Parthenopidae; cephalopods: specimens of a few sepiolidae species and specimens in bad condition of the genus Alloteuthis, which were assigned to the most abundant Alloteuthis species caught at the station).

Consequently, the total number and biomass by species was estimated for each haul (see MEDITS manual for details [48]). In the present study, we used data from the surveys carried out during 2005–2014. Given that surveys were not conducted in 2007 and 2009–2013 in the study area, the available data refer to 2005, 2006, 2008 and 2014. Both the sampling scheme and the vessel used for the surveys were the same during the study period

Data analysis

For each sampling station, two indices by faunistic category (fish, cephalopods and crustaceans) were estimated: (a) species richness, defined as the total number of species captured at each station and (b) the Shannon-Weaver (H΄) index (S1 Table) [49]:


where, pi = ni/N, ni = number of individuals in the ith species, N = total number of individuals and s = number of species

The index was computed using the “vegan” library under the R language environment [50].

A Generalized Additive Modelling (GAM) approach [51] was used to examine the effect of depth on the above indices by subarea for each faunistic category, separately. The “year” (as a factor) was also tested as a candidate predictor variable in the GAM models to account for accidental sampling inconsistencies among years due, for instance, to bad weather conditions.

Hence, the general form of the nested GAM models was:

Indexαfactor(Year) + s(Depthbysubarea) + e

Indexαs(Depthbysubarea) + e

where, a is the intercept, s indicates the smoother function of the corresponding independent variable and e is a random error term.

The smoother function used was a penalized cubic regression spline and model fitting was accomplished using the “mgcv” library [52, 53] under the R language environment. The procedure automatically selects the degree of smoothing based on the Generalized Cross Validation (GCV) score, which is a proxy for the model predictive performance. However, to avoid dubious relationships, the model was constrained to be at maximum a quartic relationship. Hence, the maximum degrees of freedom for each smoothing term was set to 4 (i.e. k = 5 in the GAM formulation).

Based on the diagnostic residual plots, it was found that a quasi-poisson distribution model and a log link function provided a pertinent fit to the data. The choice of the quasi-poisson distribution allowed for accounting for over-dispersion, which is common in biological data [54]. Best model selection was based on minimizing GCV scores and statistical inference was based on the 95% confidence level.


A total of 245 taxa were recorded in the 244 hauls analyzed. Out of the 245 taxa, 33 were elasmobranchs, 147 osteichthyes, 35 crustaceans and 30 cephalopods. Seven out of the eight alien fish species found were of Indo-Pacific origin, and one, Sphoeroides pachygaster, of Atlantic origin. Five alien species were recorded in Crete, one in Cyclades and four in Dodecanese (Table 1).

Table 1
Alien fish species recorded by subarea.

For all faunistic categories, the total number of recorded taxa was higher in the Dodecanese (S2 Table). However, the median species richness by haul for fish and cephalopods were higher in Cyclades, while the highest median species richness for crustaceans was estimated in the Dodecanese. In all subareas, species richness was higher for fish and lower for crustaceans (Fig 2). Similar patterns were observed for the median Shannon-Weaver index (Fig 3).

Fig 2
Boxplot of species richness by faunistic category and subarea (CR = Crete, CY = Cyclades, DO = Dodecanese).
Fig 3
Boxplot of the Shannon-Weaver index by faunistic category and subarea (CR = Crete, CY = Cyclades, DO = Dodecanese).

In all cases, the inclusion of the “year” factor in the GAMs resulted in models with higher or similar GCV scores and the percentage of the deviance explained was only marginally increased (S3 Table). Thus, models without the “year” effect were considered as the best ones. The analysis of deviance of the applied GAM models on species richness indicated that the effect of depth was always significant (Table 2). Regarding the Shannon-Weaver index, depth was significant in all cases, with the exception of fish in Crete and cephalopods in Cyclades (Table 3). For both indices, the variance explained was higher in the case of cephalopods and lower for fish (Tables (Tables22 and and33).

Table 2
Analysis of deviance for the applied GAM models on species richness.
Table 3
Analysis of deviance for the applied GAM models on the Shannon-Weaver index.

The effect of depth (when significant) on the examined indices is presented in Figs Figs44 and and55.

Fig 4
GAM-derived significant effects of depth on species richness for three megafauna groups (fish, cephalopod, crustacean) in the three subareas (Crete, Cyclades, Dodecanese).
Fig 5
GAM-derived significant effects of depth on the Shannon-Weaver index (H΄) for three megafauna groups (fish, cephalopod, crustacean) in the three subareas (Crete, Cyclades, Dodecanese).

Changes of fish species richness with depth followed similar patterns in Cyclades and Dodecanese. In both subareas, species richness decreased as depth increased. In the case of Crete, the fish species richness pattern was different, remaining rather stable down to 450 m and increasing thereafter (Fig 4).

In all subareas, changes of cephalopod species richness with depth followed a unimodal distribution pattern. However, the depths corresponding to the estimated distribution peaks were not always the same. Although Cyclades and Dodecanese had a similar peak at 250–280 m, the modal peak in Crete was observed at shallower depths (about 150 m) (Fig 4). The model estimates for crustaceans demonstrated that species richness in all subareas increased with depth (Fig 4).

Changes of the Shannon-Weaver index with depth for fish followed different patterns in Cyclades and Dodecanese. In Cyclades, H΄ decreased with an increase in depth. In Dodecanese, H΄ decreased until 300–350 m and increased thereafter. In Crete, the effect of depth on H΄ was not significant (Fig 5).

For cephalopods, the changes in H΄ with depth followed different patterns in Crete and Dodecanese. In Crete, H΄ was rather stable down to 400 m, showing an increasing trend (with relatively high confidence intervals) thereafter. In Dodecanese, H΄ showed a unimodal pattern with a peak at about 250 m. In Cyclades, the effect of depth on H΄ was not significant (Fig 5). For crustaceans, changes in H΄ with depth followed a similar pattern in all subareas, almost like that of species richness (i.e. increasing trend with depth) (Fig 5).


The current study adds important information to our knowledge of depth-related megafauna species diversity, by faunistic category, in three areas of the eastern Mediterranean Sea, where sampling effort is known to be much lower compared to the western Mediterranean [26]. Analysis on a fine spatial scale highlighted differences in the diversity distribution patterns among neighboring areas and faunistic categories. It must be stressed, however, that the results of the present study largely reflect only the summer conditions when the MEDITS surveys take place and thus seasonal effects cannot be taken into account.

Though species richness is assumed to be negatively influenced by low primary production and higher SST values [21, 26, 55], in the case of S. Aegean and Cretan Seas, two highly oligotrophic areas with relatively elevated SST values, a high number of taxa (245), compared to other Mediterranean areas, was recorded. In particular, the recorded numbers of elasmobranch (33), osteichthyes (147) and cephalopod (30) species were higher than the corresponding numbers of species that Colloca et al. [16] record in the western coast of Italy, (14, 128 and 26, respectively) in similar depths. The total number of fish species found (180) in the present study is close to the total number (186) recorded for the whole northern Mediterranean by Granger et al. [27]. These findings, together with previous ones [23, 24, 27], question the belief that species diversity is lower in the eastern Mediterranean [21, 26, 56]. However, the total number of cephalopod species in S. Aegean and Cretan Seas was significantly lower than the one found for the whole northern Mediterranean (58) [28], indicating a high variability in cephalopod species distribution over the Mediterranean Sea. Concerning the total number of crustacean species (35), higher numbers are recorded in different areas of W. Italy in similar depths (37 and 53) [16, 19].

Our results also indicated that alien species of Indo-Pacific origin–at least the ones available to bottom trawling–apart from Dodecanese, which is considered their gateway to the Aegean Sea [44, 45], have also extended their distribution to the Cretan Sea–corroborating findings of previous works [5761]. However, they have not established remarkable populations in Cyclades, probably due to the influence of the Black Sea less-saline cold waters in that area [37].

Our results suggested that the effect of depth on megafauna species richness was always significant, in agreement with previous studies in various Mediterranean areas [9, 10, 16, 28, 33, 34, 6265]. Additionally, species diversity remained stable over the study period, suggesting that environmental or anthropogenic pressure, such as climate change, fishery or touristic development, has not dramatically affected community structure during the last decade. Similar findings are reported from studies based on the MEDITS survey data in different Mediterranean areas [24, 27, 28].

Regarding the pattern of species richness distribution along the depth gradient, clear differences were detected among the three faunistic categories. For fish, several authors commonly report decreasing trends of species richness with depth in various areas of the Mediterranean [33, 34, 65]. A similar pattern was observed in the present study for fish species richness for Dodecanese and Cyclades. However, in Crete, the fish species richness-depth pattern was rather stable until 400 m depth, and increased deeper. The increasing species richness in the middle-slope areas can be related to the narrow shelf and upper slope of Crete, resulting from the strong steepness of the continental margin and the extended surface area of the lower slope [1, 49]. In an area also characterized by a steep and extended lower slope (Latio, W. Italy), Colloca et al. [16] did not find any significant trend in fish species richness with depth. However, other factors derived from anthropogenic disturbance could influence the fish community structure in Crete, such as coastal activities related to the tourist industry and development or the aggregation of fishing pressure on the narrow area of the shelf and upper slope [32].

Concerning the pattern of cephalopod and crustacean species richness along the depth gradient, our results agree with previous studies in other Mediterranean regions [16, 19, 28]. The observed constant increasing trend with depth of crustacean species richness could be attributed to specific adaptations of crustacean species, such as the great variety of feeding strategies, according to the availability of food resources. Indeed, different feeding guilds of plankton, benthos, benthopelagic, scavenging and detritus feeders (the latter mostly in bathyal depths) are described for crustacean species in other Mediterranean regions [6668]. This variety of feeding strategies, which allows them to survive in the harsh environmental conditions of the deep sea, where food resources are limited, is most likely the main reason for the observed high species richness.

Though the interpretation of the Shannon-Weaver (Η΄) index is rather complicated [69], a direct comparison of the diversity-depth patterns derived from both measures used in the present study can provide some useful information about the community structure. For example, in the case of fish, we can assume that, for Cyclades, the species are evenly distributed within each station as both indices show similar distribution patterns. However, at the other two subareas, the observed differences between species richness and Η΄ distribution patterns along the depth gradient may indicate the existence of a few dominant species in the intermediate waters of Dodecanese and in the middle-slope areas of Crete. Comparison of the distribution patterns of both indices for cephalopods may indicate that species are quite evenly distributed, especially in the shallow waters of Dodecanese, while a few dominant species may occur on the shelf and upper slope of Crete and Cyclades. Regarding crustacean species, the similarity of species richness and Η΄ patterns with depth indicates that species are evenly distributed.

The particularly high explained deviance (61.99%) in the GAM models including cephalopod species richness as dependent variable suggests the importance of depth-related factors in structuring cephalopod communities. The observed peak at depths around the shelf break indicates that cephalopod species are favored by higher productivity [6, 70] or other specific environmental and biological factors prevailing in the shelf break zone [16]. As many cephalopods feed in the water column, it is likely that they take advantage of the segregation of planktonic, small pelagic or benthopelagic species, and/or the high diversity of epibenthic communities around the shelf break [16, 70]. In general, the exploitation of several food sources is of great importance in highly oligotrophic environments such as the South Aegean and Cretan Seas [39].

The relatively high deviance explained for crustaceans (49.54%) indicates the importance of depth-related factors in structuring crustacean communities. Effectively, sediment texture, sediment grain size and organic matter, which are directly related to depth [71, 72], play an important role in the distribution of some crustacean species e.g. [7378]. Indeed, previous studies in the same area have reported a clear bathymetric zonation according to sediment type [50, 74].

For fish, the explained deviance was the lowest (28.75%), indicating that, besides depth-related factors, other environmental descriptors, such as habitat type and food availability, may influence community structure [39, 76]. In addition, the particular predatory and competitive abilities of fish, combined with their high mobility, allow them to develop a variety of life strategies and thus occupy larger depth ranges. For example, small-size individuals (recruits) of most species inhabit shallower depths than adults [7983].

Differences in the species richness pattern with depth were also detected among the three subareas, with more similarities existing between Cyclades and Dodecanese (Fig 4). The decreasing trend of fish diversity with increasing depth, observed in Cyclades and Dodecanese, can be attributed to the scarcity of food resources in deeper waters that can mainly support small organisms well adapted to the conditions of deep-sea environments [84]. Previous studies in the central and eastern Mediterranean show that fish assemblages on the slope are characterized by small-sized fish species, which are mainly pelagic feeders (on mesopelagic prey), with relatively high dominance and reduced species diversity [9, 16, 33, 80]. However, our results, as mentioned before, may indicate the existence of a few dominant species in the intermediate waters of Dodecanese and the middle-slope areas of Crete.

The distinct patterns observed in Crete can be partly attributed to its geographic isolation and relatively greater distance from the mainland. In general, the distance of a marine area from the mainland is a critical factor for species diversity [15]. In addition, the specific geomorphological and oceanographic characteristics of Crete [3941, 85, 86] have contributed to the formation of distinct species diversity patterns. Analogous differences in megafauna diversity or species assemblages between adjacent areas are also reported by other authors [12, 34, 87, 88].

The increase in fish species richness in the deeper waters of Crete could be related to the intermediate and deep layers of the Cretan Sea, which are considered a reservoir for heat, salt and dissolved oxygen [82], with the latter positively influencing species diversity [7]. In addition, the narrow, steep shelf and slope of Crete [38] may also contribute to the increase of species richness since it is known that the rate of species change is higher in steeper slopes [87]. Though fishery pressure in Crete is low [47], it is aggregated at the narrow surface of the shelf [32] and this may reduce the megafauna species diversity in the shallower waters.

Cephalopod species richness followed a unimodal distribution pattern along the depth gradient in all areas but differences were observed regarding the observed peaks. In Crete, the peak was observed at depths of around 150 m, while in Cyclades and Dodecanese, the corresponding peaks were observed in deeper waters around 250–280 m. These differences reflect the different depths of the shelf break in the corresponding areas [42, 89], indicating that the specific hydrological conditions prevailing in this transitional area (strong currents, fronts) favor cephalopod species diversity, as mentioned before. The results of a large-scale study on cephalopod diversity in the Mediterranean indicate that the peak of species richness in the whole South Aegean is between 200–300 m [28]. The present fine-scale analysis, however, allowed detecting differences among closely located areas which were masked in the large-scale study.

Regarding crustaceans, no differences were observed between subareas. Higher median values of species richness and H΄ were estimated for Cyclades regarding fish and cephalopod species. For crustaceans, the higher median values were observed in Dodecanese followed by Crete (Figs (Figs22 and and3).3). Such differences can be attributed to the particular characteristics of the examined areas. For instance, the flow of the nutrient-rich Black Sea waters in the western part of Cyclades [37, 90] may influence community structure and favor species diversity. Moreover, it is likely that competition among species for niche occupation is low in the extended continental shelf of the Cyclades Plateau, thus, ranges of species can overlap, a fact that leads to increased species richness [91]. The medium-level fishing pressure in Cyclades, higher in comparison to the other study subareas [47], may also reduce species competition through reduction of the dominant species, favoring a more evenly distribution of species. Additionally, the high heterogeneity of surface sediments on the Cyclades Plateau may produce numerous different microhabitats and thus high a turnover of species, resulting in high species diversity [4, 9294]. The surface sediments of the Cyclades Plateau consist mostly of biogenic sand (50–70%), with large amounts of coralline algae debris and terrigenous sand attributed to relict palaeobeach deposits (30–40%), and the deeper parts (150–250 m) are covered by sediments ranging from calcareous muddy sand to sandy mud [42]. Sediments in the other two areas have been described as terrigenous in origin, with coarser sediments in the coastal zone, calcareous muddy sand with small amounts of terrigenous silt and fine sand at the shelf break, and a more homogeneous fine sediment texture of hemi-pelagic deposition in the middle slope [78]. The sediment characteristics, particularly the high variability in grain size, is probably the main reason for the higher diversity values of crustaceans at Dodecanese and Crete (Fig 2).

Concluding, our results suggest differences in the species diversity–depth patterns depending on the faunal group or subarea studied. Differences were higher between Crete and the other two subareas than between Cyclades and Dodecanese, a fact that could be attributed to the specific geographical and oceanographic characteristics of this particular subarea. This finding supports the characterization by FAO/GFCM that Crete comprises a unique GSA area for fishery management. The analysis of species diversity at a fine spatial scale revealed differences masked in previously contacted large-scale studies. Thus, fine-scale analysis is important for the thorough understanding of factors affecting species diversity and, consequently, the identification of appropriate environmental conservation and management actions. Particularly, the requirements for the establishment of marine-protected areas, as mentioned in article 21 of the Marine Strategy Framework Directive (2008/56/EC), would greatly benefit from such analyses, which would provide essential information for the identification of diversity hot-spot areas.

Supporting information

S1 Table

Species density and Shannon-Weaver index by haul and faunistic category.


S2 Table

List of all species recorded, by subarea and faunistic category.


S3 Table

GCV scores and percentage deviance explained (in parenthesis) of the nested GAM models.



The data used for this study derived from the Data Collection Framework for the Greek Fisheries, funded by the EU and the Greek Ministry of Rural Development and Food. We would like to thank all the colleagues who participated in the MEDITS surveys during 2005–2014 and the crew of the trawler used for the sampling for their work and good collaboration. We also thank the anonymous reviewers for their constructive comments, which helped us improve the quality of our manuscript.

Funding Statement

The data used for this study derived from the Data Collection Framework for the Greek Fisheries, funded by the EU and the Greek Ministry of Rural Development and Food. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Data Availability

Data Availability

All relevant data are within the paper and its Supporting Information files.


1. Rex MA. Biological accommodation in the deep-sea benthos: comparative evidence on the importance of predation and productivity. Deep Sea Res and Oceanog Abstr. Elsevier. 1976; 2310: 975–987.
2. Osman RW, Whitlatch R B. Patterns of species diversity: fact or artefact? Paleobiology, 1978; 401: 41–54.
3. Moreno CA, Jaramillo E. The role of grazers in the zonation of intertidal macroalgae of the Chilean coasts. Oikos. 1983; 41: 73–76.
4. Etter RJ, Grassle JF. Patterns of species diversity in the deep sea as a function of sediment particle size diversity. Nature. 1992; 3606404: 576–578.
5. Gage JD. High benthic species diversity in deep-sea sediments: the importance of hydrodynamics In: Ormond R.E.G., Gage J.D. and Angel M.V. (Editors), Marine Biodiversity: Patterns and Processes. Cambridge University Press, Cambridge, pp. 1997; 148–177
6. Levin L, Etter R, Rex M, Gooday AJ, Smith CR, Pineda J, et al. Environmental influences on regional deep-sea species diversity. Annu Rev Ecol Syst. 2001; 32: 51–93. Available:
7. Levin LA.Oxygen minimum zone benthos: adaptation and community response to hypoxia. Oceanogr. Mar Biol. 2003; 41: 1–45
8. Jennings S, Greenstreet SPR, Reynolds JD. Structural changes in an exploited fish community: a consequence of differential fishing effects on species with contrasting life histories. J of Anim Ecol. 1999; 68: 617–627.
9. Labropoulou A, Papaconstantinou C. Effects of fishing on community structure f demersal fish assemblages. Belg J of Zool. 2005; 135: 191–197
10. Masuti E, Reñones O. Demersal resource assemblages in the trawl fishing grounds of the Balearic Islands (western Mediterranean). Sci. Mar. 2005; 69: 167–366.
11. Kaiser MJ, Clarke K., Hinz H, Austen MCV, Somerfield PJ, Karakassis I. Global analysis of response and recovery of benthic biota to fishing. Mar Ecol Prog Ser. 2006; 311: 1–14.
12. Keskin C, Ordines F, Guijarro B, Massuti E. Comparison of fish assemblages between the Sea of Marmara and the Aegean Sea (north–eastern Mediterranean). J Mar Biol Assoc U K. 2011; 916: 1307–1318.
13. Day DS, Pearcy WG. Species associations on benthic fish on the continental shelf and slope off Oregon. J Fish Res, Board of Canada. 1968; 25: 2665–2675.
14. Tserpes G, Peristeraki P, Potamias G, Tsimenides N. Species distribution in the southern Aegean Sea based on bottom-trawl surveys. Aquat Liv Res. 1999; 12: 167–175.
15. Rex MA, Etter RJ. Deep Sea Biodiversity, Pattern and scale. Harvard University Press, 2010; p. 354.
16. Colloca F, Cardinale M, Belluscio a, Ardizzone G. Pattern of distribution and diversity of demersal assemblages in the central Mediterranean Sea. Est Coast Shelf Sci. 2003; 56: 469–480. doi: doi: 10.1016/S0272-7714(02)00196-8
17. Macpherson E. Species range size distributions for some marine taxa in the Atlantic Ocean. Effect of latitude and depth. Biol J Linn Soc. 2003; 803: 437–455.
18. Reynolds JA. Quantifying habitat associations in marine fisheries: a generalization of the Kolmogorov–Smirnov statistic using commercial logbook records linked to archived environmental data. Can J Fish Aqu Sci. 2003; 604: 370–378.
19. Fanelli E, Colloca F, Ardizzone G. Decapod crustacean assemblages off the west coast of central Italy (western Mediterranean). Sci Mar. 2007; 711: 19–28.
20. Quignard JP, Tomasini JA. Mediterranean fish biodiversity. Biol Mar Med. 2007; 7: 1–66.
21. Ben Rais Lasram F, Guilhaumon F, Mouillot D. Fish diversity patterns in the Mediterranean Sea: deviations from amid-domain model. Mar Ecol Prog Ser. 2009; 376:253–267. doi: 10.3354/meps07786
22. Lasram FBR, Hattab T, Halouani G, Romdhane MS, Le Loc'h F. Modeling of Beta Diversity in Tunisian Waters: Predictions Using Generalized Dissimilarity Modeling and Bioregionalisation Using Fuzzy Clustering. PLoS ONE. 2015; 10(7): e0131728 doi: 10.1371/journal.pone.0131728 [PMC free article] [PubMed]
23. Gaertner J, Bertrand J, Relini G, Papaconstantinou C, Mazouni N, de Sola L, et al. Spatial pattern in species richness of demersal fish assemblages on the continental shelf of the northern Mediterranean Sea: a multiscale analysis. Mar Ecol Prog Ser. 2007; 341: 191–203. doi: 10.3354/meps341191
24. Gaertner JC, Maiorano P, Mérigot B, Colloca F, Politou CY, Gil De Sola L,et al. Large-scale diversity of slope fishes: pattern in consistency between multiple diversity indices. PLoSOne. 2013; 8: doi: 10. 1371/journal.pone.0066753 [PMC free article] [PubMed]
25. Coll M, Piroddi C, Steenbeek J, Kaschner K., Ben RaisLasram F, Aguzzi J, et al. () The biodiversity of the Mediterranean Sea: estimates, patterns, and threats. PLoSOne. 2010; 5:e 11842 doi: 10.1371/ journal.pone.0011842 [PMC free article] [PubMed]
26. Coll M, Piroddi C, Albouy C, Ben RaisLasram F, Cheung WWL, Christensen V, et al. The Mediterranean Sea under siege: spatial overlap between marine biodiversity, cumulative threats and marine reserves. Glob Ecol Biog. 2012; 21:465–480. doi: 10.1111/j.1466-8238.2011.00697.x
27. Granger V, Fromentin JM, Bez N, Relini G, Meynard C, Gaertner JC, et al. Large-scale spatio-temporal monitoring highlights hotspots of demersal fish diversity in the Mediterranean Sea. Prog Oceanogr. 2015; 5–130: 65–74. doi: 10.1016/j.pocean.2014.10.002
28. Keller S, Bartolino V, Hidalgo M, Bittero I, Casciaro L, et al. Spatio -Temporal Patterns of Mediterranean Cephalopod Diversity. PLoSOne. 2016; doi: 10.1371/ journal.pone.0146469 [PMC free article] [PubMed]
29. Bacheler NM, Bailey KM, Ciannelli L, Bartolino V, Chan KS. Density-dependent landscape and climate effects on spawning distribution of walleye Pollock, Theragra chalcogramma. Mar Ecol Prog Ser. 2009; 391: 1–12.
30. Bacheler NM, Ciannelli L, Bailey KM, Bartolino V. Do wall eye pollock exhibit flexibility in where or when they spawn based on variability in water temperature? Deep Sea Res Part II Top Stud Oceanogr. 2012; 65–70:208–216.
31. Tsimenides N, Tserpes G, Machias A, Kallianiotis A. Distribution of fishes on the Cretan shelf. J Fish Biol. 1991; 39: 661–672.
32. Tserpes G, Tzanatos E, Peristeraki P. Spatial management of the Mediterranean bottom-trawl fisheries; the case of the southern Aegean Sea. Hydrob. 2011; 670: 267–274.
33. Kallianiotis A, Sophronidis K, Vidoris P, Tselepides A. Demersal fish and megafaunal assemblages of the Cretan continental shelf and slope (NE Mediterranean): seasonal variation in species density, biomass and diversity. Prog Ocean. 2000; 46: 429–455.
34. Vassilopoulou V, Papaconstantinou C, Bekas P. Demersal fish community patterns in two areas of the Aegean Sea, Greece. It J Zool. 1998; 65-S1, 293–297, doi: 10.1080/11250009809386835
35. Valavanis VD, Pierce GJ, Zuur AF, Palialexis A, Saveliev A, Katara I, et al. Modelling of essential fishhabitat based on remote sensing, spatial analysis and GIS. Hydrobiologia. 2008; 612:5–20.
36. Bertrand J, Gil de Sola L, Papacostantinou C, Relini G, Souplet A. The general specifications of the MEDITS surveys. Sci Mar. 2002; 66: 9–17.
37. Nielsen JN. Hydrography of the Mediterranean and Adjacent Waters, in Report of the Danish Oceanographic Expedition 1908–1910 to the Mediterranean and Adjacent Waters, V. 1, Copenhagen. 1912; pp. 72–191.
38. Lykousis V, Chronis G, Tselepides A, Price NB, Theocharis A, Siokou-Frangou I, et al. Major outputs of the recent multidisciplinary biochemical researches undertaken in the Aegean Sea. J Mar Syst. 2002; 33–34: 313–334
39. Psarra S, Tselepides A, Ignatiades L, Dafnomili E. Primary production estimates in the Cretan Sea. In: Tselepides A, Papadopoulou K-N, Polychronaki T (eds) CINS: Pelagicbenthic coupling in the oligotrophic Cretan Sea, MAST-II Mediterranean Targeted Project. 1996; pp 46–56
40. Danovaro R, Della Croce N, Dell’Anno A, Marrale A, Pusceddu A. Seasonal and spatial distribution of pico-,nano-, and microparticulate matter and bacterioplankton in the Eastern Mediterranean Sea In: Tselepides A, Papadopoulou K-N, Polychronaki T (eds) CINCS: Pelagic-benthic coupling the oligotrophic Cretan Sea, MAST-II Mediterranean Targeted Project; 1996; pp 69–81
41. Chronis G, Stavrakakis S, Danovaro R, Tselepides A, Wassmann P, Duineveld G, et al. Downward fluxes In: Tselepides A, Papadopoulou KN, Polychronaki T (eds) CINCS: Pelagic-benthic coupling the oligotrophic Cretan Sea, MAST-II Mediterranean Targeted Project; 1996; pp 82–111.
42. Volakis S, Anagnostou C. Modern sedimentation at the North part of the Cyclades Plateau. Proc. 4th Nat. Symp. of Oceanography and Fishery, 1993; Rhodes, 101–104.
43. Ovchinnikov I.M. Circulation in the surface and intermediate layers of the Mediterranean. Oceanology. 1966; 6: 48–59.
44. Pancucci-Papadopoulou MA., Raitsos DE, Corsini-Foca M. Biological invasions and climate warming: implications for south-eastern Aegean ecosystem functioning. J Mar Biol Assoc U K. 2012; 92: 777–789.
45. Pancucci-Papadopoulou MA., Raitsos DE, Corsini-Foca M. Biological invasions and climate warming. The Marine Biologist, 2013; 8–11.
46. Azzurro E, Moschella P, Maynou F. Tracking Signals of Change in Mediterranean Fish Diversity Based on Local Ecological Knowledge. PLoS ONE. 2011; 6–9: e24885 doi: 10.1371/journal.pone.0024885 [PMC free article] [PubMed]
47. Sylaios G.K., Koytroumanidis T., Tsicliras A.C. Ranking and classification of fishing areas using fuzzy models and techniques. Fish Manag Ecol, 2010, 17, 240–253
49. Hill M.O. Diversity and evenness: a unifying notation and its consequences. Ecol. 1973; 54: 427–473.
50. R Development Core Team. R: A language and environment for statistical computing R Foundation for Statistical Computing, Vienna, Austria: 2013; ISBN 3-900051-07-0, URL
51. Hastie TJ, Tibshirani RJ. Generalized additive models Chapman and Hall, London, 352 pp. 1990.
52. Wood SN. Modelling and smoothing parameter estimation with multiple quadratic penalties. J R Stat Soc Ser B. 2000; 62: 413–428.
53. Wood SN. Generalized Additive Models: An introduction with R. Florida: Chapman and Hall/CRC, 2006; 391 pp.
54. Harrison XA. Using observation-level random effects to model overdispersion in count data in ecology and evolution. PeerJ. 2014; 2:e616; doi: 10.7717/peerj.616 [PMC free article] [PubMed]
55. Begon M., Townsend C. R., Harper J. L. Ecology: From individuals to ecosystems Malden, MA: Blackwell Pub; 2006
56. Mouillot D, Villéger S, Scherer-Lorenzen M, Mason NWH Functional Structure of Biological Communities Predicts Ecosystem Multifunctionality. PLoS ONE 2011; 6(3): e17476 doi: 10.1371/journal.pone.0017476 [PMC free article] [PubMed]
57. Peristeraki P, Lazarakis G, Skarvelis K, Georgiadis M, Tserpes G. Additional records on the occurrence of alien fish species in the eastern Mediterranean Sea. Med Mar Sci. 2007; 7: 61–67
58. Kasapides P, Peristeraki P, Tserpes G, Magoulas A. First record of the Lessepsian migrant Lagocephalus sceleratus (Gmelin 1789) (Osteichthyes: Tetraodondidae) in the Cretan Sea (Aegean, Greece). Aquat Inv. 2007; 2: 71–73.
59. Kasapides P, Peristeraki P, Tserpes G, Magoulas A. A new record of the Lessepsian invasive fish Etrumeus teres (Osteichthyes: Clupeidae) in the Mediterranean Sea (Aegean, Greece). Aquat Inv. 2007; 2–2: 152–154.
60. Golani D, Appelbaum-Golani B, Peristeraki P. Westward range extension of the Lessepsian migrant the Shrimp scad Alepes djedaba (Forsskål, 1775) in the Mediterranean. Ann Ser Hist Nat. 2013; 232: 115–118.
61. Skarvelis K,Giannakaki A, Kosoglou G, Peristeraki P. Data on frequency of occurrence ofalien species caught on the coastal fisheries of Crete. Proc. 11th Panhell Symp Ocean Fish. 2015.
62. Company JB, Maiorano P, Tselepides A, Politu C-Y, Plaity W, Rotllant G, Sardà F. Deep-sea decapod crustaceans in the western and central Mediterranean Sea: a preliminary aspects of species distribution, biomass and population structure. Sci Mar. 2004; 683:73–86. doi: 10.3989/scimar.2004.68s373
63. Massuti E, Renones O, Carbonell A, Oliver P. Demersal fish communities exploited on the continental shelf and slope off Majorca (Balearic islands, NW Mediterranean). Vie Milieu. 1996; 46: 45–55.
64. Moranta J, Stefanescu C, Massutı E, Morales-Nin B, Lloris D. Fish community structure and depth related trends on the continental slope of the Balearic Islands (Algerian basin, western Mediterranean). Mar Ecol Prog Ser. 1998; 171: 247–259.
65. Tecchio S, Ramirez-Llodra E, Sarda F, Company JB. Biodiversity of deep-sea megafauna in western and central Mediterranean basins. Sci Mar. 2011; 75(2): 341–350. doi: 10.3989/scimar.201175n2341
66. Cartes JE. Feeding strategies and partition of food resources in deep-water decapod crustaceans in relation to depth (between 400 and 2300 m). J Mar Biol Assoc U K 1998; 78:509–524
67. Cartes JE. Diets of deep-water pandalid shrimps on the western Mediterranean slope. Mar. Ecol Prog Ser. 1993; 96: 49–61.
68. Cartes JE,Carrason M. The influence of trophic variables in the depth-range distribution and zonation rates of deep-sea megafauna: the case of the Western Mediterranean assemblages. Deep-Sea Res. 2004; 51: 263–279.
69. Magurran AE. Measuring biological diversity Oxford: Blackwell Science Ltd; 2004; 256 p.
70. Navarro J, Coll M, Somes C J, Olson RJ. Trophic niche of squids: Insights from isotopic data in marin systems worldwide. Deep-Sea Res. 2013; II- 95: 93–102
71. Chronis G, Lykousis V, Georgopoulos D, Zervakis V, Stavrakakis S, Poulos S. Suspended particulate matter and nepheloid layers over the southern margin of the Cretan Sea (N.E. Mediterranean): seasonal distribution and dynamics. Prog Ocean. 2000; 46: 163–185
72. Gray JS. The ecology of marine sediments Cambridge studies in modern biology 2 1981; Cambridge University Press, Cambridge
73. Atkinson RJA, Taylor AC. Physiological ecology of burrowing decapods. In: A.A. Fincham and P.S. Rainbow (eds.), Aspects of decapod crustacean biology, Symp. Zool. Soc. Lond. 1988; 59: 201–226
74. Macpherson E. Biogeography and community structure of the decapod crustacean fauna off Namibia (southeast Atlantic). J. Crust. Biol. 1991; 11: 401–415.
75. Pinheiro MAA, Fransozo A, Negreiros-Fransozo ML. Distribution patterns of Arenaeus cribarius (Lamarck, 1818) (Crustacea, Portunidae) in Fortaleza Bay, Ubatuba (SP), Brazil. Rev Bras Biol. 1996; 56: 705–716.
76. Comeau M, Conan GY, Maynou F, Robichaud G, Therriault JC,Starr M. Growth, spatial distribution, and abundance of benthic stages of the snow crab (Chionoecetes pilio) in Bonne Bay, Newfoundland, Canada. Can J Fish Aquat Sci. 1998; 55: 262–279.
77. Vaz S, Carpentier A, Coppin F. Eastern English Channel fish assemblages: measuring the structuring effect of habitats on distinct sub-communities.–ICES J Mar Sci. 2007; 64: 271–287.
78. Emylianov EM, Schimkus KM. Geochemistry and sedimentology of the Mediterranean Sea. 1986; Dordrecht D. Reidel.
79. Stefanescu C, Rucabado J, Lloris D. Depth -size trends in the western Mediterranean demersal deep-sea fishes. Mar Ecol Prog Ser. 1992; 81: 205–213.
80. Madurell T, Cartes JE, Labropoulou M. Changes in the structure of fish assemblages in a bathyal site oft he Ionian Sea (eastern Mediterranean). Fish Res. 2004; 66: 245–260
81. Bartolino V, Ottavi A, Colloca F, Ardizzone, Domenico Gian, Stefánsson G. Bathymetric preferences of juvenile European hake (Merluccius merluccius). ICES J Mar Sci. 2008; 65: 963–969. Available:
82. Tserpes G, Politou CY, Peristeraki P, Kallianiotis A, Papaconstantinou C. Identification of hake distribution pattern and nursery grounds in the Hellenic seas by means of generalized additive models. Hydrob. 2008; 612:125–133. doi: 10.1007/s10750-008-9486-x
83. Druon JN, Fiorentino F, Murenu M, Knittweis L, Colloca F, Osio C, et al. Modelling of European hake nurseries in the Mediterranean Sea: an ecological niche approach. Prog Ocean. 2015; 188–204.
84. Hutchinson GE. Hommage to Santa Rosalia, or, Why are there so many kinds of animals? Am. Nat. 1959; 93: 154–159.
85. Souvermezoglou E, Krasakopoulou E, Pavlidou A. Modifications of the nutrients and oxygen regime between the Cretan Sea and the Eastern Mediterranean (1986–1995). Proc Internat POEM-BC/MTP Symp. Molitg les Bains-France, 1996; 133–136.
86. Georgopoulos D, Theocharis A, Zodiatis G. Intermediate water formation in the Cretan Sea (S. Aegean Sea). Oceanol. Acta. 1989; 124: 353–359.
87. Kallianiotis A, Vidoris P, Sylaios G. Fish species assemblages and geographical sub-areas in the North Aegean Sea, Greece. Fish Res. 2004; 68: 171–187.
88. Garofalo G, Fiorentino F, Gristina M, Cusumano S, Sinacori G. Stability of spatial pattern of fish species diversity in the Strait of Sicily (central Mediterranean). Hydrob. 2007; 580:117–124 doi: 10.1007/s10750-006-0460-1
89. Theocharis A, Georgopoulos D, Lascaratos A, Nittis K. Water masses and circulation in the central region oft he eastern Mediterranean: eastern Ionian, south Aegean and northern Levantine, 1986–1987. Deep-Sea Res. 1993; II- 40: 1121–1142.
90. Hopkins TS. Physical processes in the Mediterranean basins In Estuarine Transport Processes, Kjerfve B (ed.). University of South Carolina Press: Columbia, SC: 1978; 285–286.
91. Terborgh J. Distribution on Environmental Gradients: Theory and a Preliminary Interpretation of Distributional Patterns in the Avifauna of the Cordillera Vilcabamba, Peru Ecology. 1971; 52–1: 23–40.
92. Levin SA. Dispersion and population interactions. Am. Nat. 1974; 108: 207–228
93. Jorissen F J, De Stigter H C, Widmark JGV. A conceptual model explaining benthic foraminiferal microhabitats. Mar Micropaleont. 1995; 22: 3–15.
94. Hecker B. Variation in megafaunal assemblages on the continental margin south of New England. Deep-Sea Res. 1990; 37: 37–57.

Articles from PLoS ONE are provided here courtesy of Public Library of Science