The large-scale distribution pattern of megafauna communities along the Mediterranean middle slope was explored. The study was conducted between 500 and 800 m depth where deep-water fishery occurs. Although community studies carried out deeper than 500 m are partly available for some geographic areas, few large-scale comparative studies have been carried out. Within the framework of the MEDITS survey programme, we compared the megafauna community structure in ten geographical sub-areas (GSAs) along the Mediterranean coasts. Additionally, the spatial distribution of fishing was analysed using vessel monitoring by satellite information. Overall, the community showed a significant difference between sub-areas, with a decreasing eastward pattern in abundance and biomass. Longitude was the main factor explaining variation among sub-areas (by generalized additive models). However, we found a region which did not follow the general pattern. GSA 6 (northern Spain) showed significantly lower abundance and a different composition structure to the adjacent areas. The decrease in community descriptors (i.e. abundance and biomass) in this area is probably a symptom of population changes induced by intense fishery exploitation. Overall, a combination of environmental variables and human-induced impacts appears to influence the bentho-pelagic communities along the slope areas of the Mediterranean.
RESUMEN
En este estudio se describe la estructura y patrones de distribución de la comunidad de megafauna que habita en el margen continental medio a lo largo del Mediterráneo. El estudio se realizó entre los 500 y 800 m, coincidiendo espacialmente con las pesquerías de profundidad. A pesar de que se conoce parcialmente la estructura de las comunidades que habitan por debajo de 500 m, existe la necesidad de estudiar estas comunidades a una escala espacial más amplia. Dentro del marco del proyecto internacional MEDITS, se comparó la estructura de las comunidades en diez sub-áreas geográficas (GSAs) a lo largo de las costas mediterráneas. Además se analizó la distribución espacial del esfuerzo pesquero utilizando la información de los datos de seguimiento de buques. En general los resultados mostraron diferencias significativas entre subáreas mostrando un patrón decreciente en los valores de biomasa hacia el este, siendo la longitud el principal factor explicativo del modelo de distribución (GAMs). Sin embargo, encontramos una subárea que no seguía el patrón general, la GSA6 (norte de España). La GSA6 mostró una biomasa y estructura de la comunidad diferente a las áreas adyacentes. La disminución de la biomasa en esta área parece ser un síntoma de los cambios poblaciones causados por la elevada intensidad de pesca en la zona. Los resultados sugieren que la distribución y estructura de las comunidades bento-pelágicas parecen estar moduladas por la combinación de las variables ambientales y los impactos producidos por la actividad humana.
The continental slopes are the connecting habitat between the shelf (200 m) and the deep basins (4000 m). They play an important role in the global carbon cycle, cycling nutrient and providing valuable food and energy resources (Levin and Dayton 2009). These biological systems cannot support high levels of exploitation due to their low productivity and the presence of productive but not yet well-identified habitats (de Juan and Lleonart 2010, Norse et al. 2012, Watson and Morato 2013). Several deep-water species have biological characteristics that make them more vulnerable to fishing exploitation: K-type life-history traits, low fecundity and aggregation behaviour in restricted topographic areas (Merrett and Haedrich 1997, Clark et al., 2016, Fernandez-Arcaya et al. 2016). However, several major sources of impacts are threatening continental slopes around the world (i.e. fishing, marine litter, chemical contamination, species invasion, ocean acidification and climate change–related stressors) (Galil 2009, Pham et al. 2014, Fernandez-Arcaya et al. 2017), including those of the Mediterranean Sea, whose physical characteristics make it particularly sensitive to these impacts (Coll et al. 2012, Ramirez-Llodra et al. 2010, 2013). The Mediterranean is a semi-enclosed area characterized by oligotrophic conditions (Margalef 1985). The biological production decreases from north to south and west to east and is inversely related to the increase in temperature and salinity, being the eastern basin one of the poorest marine areas in the world (Azov 1991, Danovaro et al. 2010).
Among the human activities that are impacting the Mediterranean continental slope and its associated biological communities, bottom trawling is the most important and the cause of most concern (Guijarro et al. 2017, Piroddi et al. 2017). Mediterranean deep-water fisheries began in the first decades of the last century, coinciding with the development of new technologies that made fisheries in deeper waters possible (Relini and Orsi Relini 1987, Spedicato et al. 1995, Sardà and Demestre 1987). Larger vessels with more powerful winches and stronger cables were increasingly used from Mediterranean harbours (Farrugio et al. 1993, Guijarro et al. 2011, Quetglas et al. 2017). Two main deep-sea bottom fisheries have extended along the Mediterranean continental slope: (i) an upper slope (~200-600 m) fishery targeting Nephrops norvegicus and Parapenaeus longirostris, and ii) a middle slope (~400-800 m) fishery targeting Aristaeomorpha foliacea and Aristeus antennatus (Sardà et al. 2004a).
The scientific assessment conducted on deep-sea Mediterranean shrimp fisheries is giving rise to alarm about overexploitation and major damage to the associated fauna (Gorelli et al. 2016, Colloca et al. 2017, Vielmini et al. 2017). However, the community associated with red shrimp is not well known and large spatial scale analysis is still lacking. Mediterranean Sea community studies carried out in continental margin areas deeper than 500 m are partly available for some geographic areas. Most of these studies have been conducted in the NW Mediterranean Sea, which is currently one of the most studied deep-sea regions of the word (Sardà et al. 2004b, Massutí and Reñones 2005, Fernandez-Arcaya et al. 2016, among others). Additionally, the published data focus on specific taxonomic groups such as crustaceans (Company et al. 2004, Follesa et al. 2009), cephalopods (Keller et al. 2016) and fishes (D’Onghia et al. 2004, Moranta et al. 2008, Granger et al. 2015), while studies of non-crustacean invertebrate communities are still very limited (Cartes et al. 2009, Gori et al. 2013, Mecho et al. 2014). Thus, there is an increasing need for better understanding of the community distribution patterns as a whole, in particular in the context of an ecosystem management approach that assesses the whole community rather than evaluating single resources (Quetglas et al. 2017). Additionally, it is essential to identify the main drivers affecting the spatial distribution of deep-sea communities in order to evaluate their distribution patterns (Navarro et al. 2016, Rowden et al. 2017, Allen et al. 2018). This information is also needed in order to identify the major impacts and implement effective management plans within deep-sea habitats. Previous studies suggest that, in addition to the environmental variables, human-induced impacts also affect species composition, abundance and distribution (Worm et al. 2006, Colloca et al. 2017, Navarro et al. 2016). However, these patterns are far from being well understood.
The objective of the present work is to describe and compare abundance, biomass and assemblage structure of the whole megafauna community along the Mediterranean continental slope, using independent MEDITS fishing cruise data. A second objective is to analyse the relationship between environmental variables, fishing intensity and the megabenthic community distribution in deep-water shrimp fishing grounds in order to better understand the dominant forces driving the distribution patterns of deep-sea assemblages.
MATERIALS AND METHODS
The study area comprised a wide area of the Mediterranean Sea, from the Alboran Sea to the western Ionian Sea (Fig. 1). The study covers ten geographical sub-areas (GSAs) with specific environmental and biological characteristics defined by the General Fisheries Commission for the Mediterranean (GFCM): the northern Alboran Sea (GSA 1), Alboran Island (GSA2), the Balearic Islands (GSA 5), northern Spain (GSA 6), the Gulf of Lions (GSA 7), eastern Corsica (GSA 8), the central-southern Tyrrhenian Sea (GSA10), the Strait of Sicily (GSA 16), the southern Adriatic Sea (GSA 18) and the western Ionian Sea (GSA 19).
Map of the Mediterranean Sea showing the sampling stations in each GSA. The number of hauls in each GSA is shown in parenthesis. GSA 1 (48), GSA 2 (16), GSA 5 (40), GSA 6 (35), GSA 7 (17), GSA 8 (11), GSA 10 (92), GSA 16 (122), GSA 18 (60), GSA 19 (119).
The sampling of benthic megafauna was carried out with an experimental bottom trawl net (IFREMER reference GOC73), the standard sampler for study of deep-sea megafauna (Cartes et al. 2017, Foveau et al. 2017), following a standardized survey design and sampling methodology. The average haul rate for all the GSAs was proportional to the total surface area (one station per 60 km²). The details of sampling methodology can be found on Bertrand et al. 2002a, b and in the MEDITS Handbook (http://www.sibm.it/MEDITS%202011/principaledownload.htm).
Only hauls between 500 and 800 m (559) were selected, in order to analyse the preferred habitat of deep shrimp fishing grounds (Sardà et al. 2004a) from 2012 to 2015. The total catch from each haul was sorted, identified to the lowest taxonomical level possible, counted and weighed following the MEDITS survey protocol. Only taxonomical groups that were appropriately sampled and identified in all GSAs were included. The species were grouped by phylum, except for Actinopterygii and Mollusca, which were grouped by subgroups because of their different life history and moving capacity. As a result, nine groups were used for the analysis: Osteichthyes, Chondrichthyes, Crustacea, Cephalopoda, other Mollusca, Echinodermata, Cnidaria, Porifera and Tunicata.
Total abundance (n km–2) and biomass (g km–2) indexes of each species were standardized by swept area (in km2). Net horizontal opening and recorded distances were used to calculate the swept area. After standardization, specific data were placed in the above groups and their percentage and biomass index descriptors were box-plotted by GSAs.
One-way ANOVA tests on log10-transformed data of abundance and biomass indexes of the entire community, followed by pairwise Tukey honestly significant difference tests, were applied to test for significant (Ndifferences among GSAs. Previously, all data were tested for normality using the Shapiro-Wilk test. In this and forward analyses, the selected index was biomass because for several groups it can be complicated to enumerate the individuals caught, as is the case of colonial organisms (like corals or tunicates) or sponges (see Results section for more details).
In addition, data matrices (species biomass vs haul) were square-root transformed and similarity between all samples was calculated using the Bray-Curtis measure. A group average cluster analysis on transformed biomass indices was carried out, and a similarity profile test (SIMPROF) permutation routine (1000 restarts) was applied to test for the significance of genuine clustering. Similarities between sampled hauls were also visualized using a multidimensional scaling plot (MDS) and similarity percentage (SIMPER) analyses were run to determine species contributions to the observed communities. Analysis of similarity (ANOSIM) was used to test for differences in community composition by zones and years. Ecological analyses were performed with the PRIMER-E 6 and PERMANOVA software (Clarke and Gorley 2006).
Generalized additive models (GAMs) were used to determine significant relationships between biomass index (response variables) and latitude and longitude of the hauls (explanatory variables) for each group and for the entire community. The models were defined as follows: (y) ~β0 + s1(x1) + s2(x2), where β0 is the intercept, s represents a polynomial smooth function, substituting the slope parameter in a linear regression, and x represents explanatory variables (Wood 2006). The response variable (biomass index) was log-transformed to approximate normal distribution. The most parsimonious models were identified for the forward selection using the GCV scores (Burnham and Anderson 2003). The mgcv statistical package in R statistical software (R3.3.2) (http://www.r-project.org/) was used to fit the models.
Fishing effort
The fishing effort in the area was assessed by means of the Vessel Monitoring System (VMS) . To reduce the spatial derived effect on community distribution (e.g. the longitude effect), we used the available VMS data from the Spanish Mediterranean Area (GSAs 1, 2, 5 and 6) for the two-years study period (2013 and 2014). VMS data were provided by the Spanish Ministry of Agriculture, Food and Environment and consisted of records containing data on the geographic position, date, time and instantaneous velocity for each boat, approximately every two hours. Bottom trawlers usually fish at velocities of between 2 and 3.5 knots, so only signals showing these velocities were included in the analysis so as to remove VMS signals from boats transiting to fishing grounds or ports. To estimate the fishing effort, a grid of 0.03×0.03 degrees was created around each haul to calculate the number of VMS signals associated with it.
GAM models were applied to analyse the effect of fishing effort and spatial distribution (latitude and longitude) on the community biomass (see above for more details).
RESULTS
A total of 432 species belonging to 9 different taxonomical groups were selected. The overall taxonomic composition showed that the most diverse groups were Osteichthyes (127 sp.) and Crustacea (91 sp.), followed by the other 6 groups (54 Echinodermata, 38 Cephalopoda, 36 Cnidaria, 26 Chondrichthyes, 36 Mollusca (non-cephalopod), 10 Tunicata and 8 Porifera). The percentages of abundance and biomass for each taxonomic group for all GSAs are shown in Figure 2. In abundance, Osteichthyes and Crustacea dominated the fauna overall, except in Alboran Sea (GSAs 1 and 2) where Chondrichthyes were the dominant taxa (33.7% and 43%, respectively, considering the biomass indexes). In terms of biomass Osteichthyes and Chondrichthyes (mostly Galeus melastomus species) dominated the community in these GSAs. The more sessile taxonomical groups (benthic species) showed a lower contribution to the overall community. They were therefore plotted separately in order to visualize their patterns of dominance (Fig. 2C, D). Echinodermata and Cnidaria dominated the benthic fauna in terms of abundance and biomass. However, abundance and biomass showed different patterns, possibly because of the high weight of some hauls conducted over mud facies of corals, such as Isidella. In fact,biomass was composed of Cnidaria in the central basin (GSA 8, 99.5%, GSA 10, 79%, GSA 16, 73.7%, GSA 18, 95%), mostly the coral Isidella elongata, and of the actinia Actinauge richardi in GSA 16 (Fig. 2D). In contrast, the biomass in the western GSAs was dominated by several species of holothurians (Molpadia musculus and Parastichopus regalis), ophiouras (Leptometra sp. and Hymenodiscus coronata) and sea-urchins (Cidaris cidaris) (all the mentioned species together: GSA 1, 69.7%; GSA 2, 66.6%; GSA 5, 56.8%; and GSA 7, 64.5%) (Fig. 2D). In terms of abundance, in GSA 6 (northern Spain) other Mollusca represent more than 50% of the community. This percentage was mostly caused by the high densities of the species Aporrhais serresianus (Fig. 2C).
Percentage of abundance (N km–2) (A) and biomass (g km–2) indexes (B) for each taxonomic group by GSAs. Below, detailed figure of abundance (C) and biomass indexes (D) of the most sessile taxonomic groups. Ostei, Osteichthyes; Chon, Chondrichthyes; Crus, Crustacea; Cepha, Cephalopoda; Cnid, Cnidaria; Echi, Echinodermata; OtherMol, other Mollusca (non-cephalopod); Porif, Porifera; Tunic, Tunicata.
The overall community box-plot (Fig. 3) showed a significant difference between areas (one-way ANOVA, F=227.9, p<0.001), with a decreasing eastward pattern in abundance and biomass indices, except in GSA 6 (northern Spain). GSA 6 showed lower levels than the rest of the sub-areas of the western Mediterranean (GSAs 1, 2, 5, 7 and 8, with no significant differences), and equivalent ones to those in the more eastern areas. The sampled sub-areas in the eastern area (GSAs 10, 16, 18 and 19) showed no significant difference among them (Table 1, Fig. 3). The spatial pattern was consistent within all taxa, showing significantly lower values for demersal (e.g. Osteichthyes, Chondrichthyes Crustacea and Cephalopoda) and benthic taxa (Cnidaria, Echinodermata, other Mollusca, Porifera and Tunicata) in GSAs 6, 10, 16, 18 and 19 (Fig. 4). A haul conducted over the crinoid beds (Leptometra sp.) in 2013 was responsible for the abundance peak (and high standard error) observed in GSA 7 (Fig 4B). A close-up of Figure 4D showed a peak of biomass values in GSAs 8 and 18 caused by several hauls conducted over an octocoralline community of I. elongata and a peak of Echinodermata Holothuria tubulosa (GSA 2).
Total catch in the ten sampled GSAs, considering standardized abundance (N km–2) and biomass (g km–2) indexes. Boxes are the lower to the upper quartile and the centre line is the median with whiskers from minimum to maximum values.
P values of Tukey post hoc tests after one-way ANOVA for biomass with GSAs as fixed factors, F=227.9. * significant differences p<0.01.
GSA
1
2
5
6
7
8
10
16
18
1
2
0.29
5
0.88
0.02
6
*
*
*
7
*
*
0.20
*
8
0.99
0.39
0.99
*
*
10
*
*
*
1.0
*
0.25
16
*
*
*
0.59
*
*
0.24
18
*
*
*
0.70
*
*
0.34
0.96
19
*
*
*
0.16
*
*
0.16
0.98
*
Abundances (N km–2) and biomass (g km–2) indexes by GSAs of demersal taxonomic groups (A and B) and more sessile taxonomic groups (C and D) (values are mean ±SE). Ostei, Osteichthyes; Chon, Chondrichthyes; Crus, Crustacea; Cepha, Cephalopoda; Cnid, Cnidaria; Echi, Echinodermata; OtherMol, other Mollusca (non-cephalopod); Porif, Porifera; Tunic, Tunicata.
A group average sorting dendrogram showing the percentage similarity of the megabenthic community and the corresponding MDS ordination at each station sampled is shown in Figure 5. Grouping appeared to be strongly influenced by zone and suggests a longitudinal pattern of species distribution. The analysis delineated five clusters in two main community-level groups at the 25% level of similarity (Fig. 5A). These two main groups delineated species from the western and central basin as follows: Group a) a more heterogeneous group formed by samples from the western Mediterranean Basin (GSAs 1, 2, 5, 7 and 8); and Group b) a more homogeneous group formed by samples from the central Mediterranean Basin (GSA 10, 16, 18 and 19) and northern Spain (GSA 6). Three subgroups were formed by western basin samples, based again on longitude. These subgroups overlapped with GSAs, except GSA 7, which is widely distributed in space and showed a transitional zone between the Balearic basin and Corsica (Fig. 5B, C). From Group b, all hauls from the central basin were clustered together and separated by samples from northern Spain (Fig. 5A).
Cluster sorting dendrogram showing the percentage similarity of the hauls in relation to the fauna collected (A) and the corresponding two-dimensional MDS ordination (B) by the sorting cluster based upon an arbitrary level of similarity of 0.35; a and b are the two main groups. C, a map of the Mediterranean Sea showing the haul distribution by sorting cluster. ALBO, Alboran Sea and Alboran Island; ESPN, northern Spain; GL-COR, Gulf of Lions and Corsica; CENTM, central Mediterranean Basin; BAL-GL, Balearic basin and Gulf of Lions.
Ordination of data showed good correspondence with clusters, with the samples being adequately represented (stress=0.16) in two-dimensional space (Fig. 5B). The ANOSIM test performed on the groupings showed that the clusters were significantly different from one another (global R=0.86, p<0.01; Table 2). Analysis of similarities of the complete data set showed an overall significant difference between GSAs (R=0.67, p<0.01; Table 3). In contrast, there were no significant temporal differences among samples (R=0.03, p=0.1).
Pairwise test conducted by ANOSIM using the factor Cluster. p<0.01 for all pairwise comparisons. Global R=0.85.
Clusters
CENTM
ALBO
BAL-GL
COR-GL
CENTM
ALBO
1.00
BAL-GL
0.92
0.76
COR-GL
0.94
0.95
0.67
ESPN
0.57
1.00
0.62
0.99
Pairwise test conducted by ANOSIM using the factor GSA. p<0.01 for all pairwise comparisons. ´ no significant difference (p=36.1). Global R=0.67.
GSA
10
1
16
18
19
2
5
6
7
10
1
1.00
16
0.40
0.99
18
0.73
0.98
0.31
19
0.43
0.99
0.33
0.24
2
1.00
0.02´
1.00
1.00
1.00
5
0.97
0.82
0.96
0.96
0.94
0.80
6
0.83
0.97
0.76
0.84
0.69
1.00
0.74
7
1.00
0.66
0.96
0.96
0.96
0.69
0.69
0.92
8
1.00
0.77
0.92
0.93
0.97
0.99
0.83
1.00
0.57
SIMPER analysis was used to identify the species that contributed to the highest similarity within clusters (Supplementary material, Table S1). The clusters formed exclusively by areas, such as ALBO (Alboran Sea and Alboran Island) and ESPN (northern Spain) showed the highest values of average similarity, 55.6% and 48. 3% respectively, and three other clusters showed a similar percentage of similarity (close to 40%). According to SIMPER, although 43 species were listed, there were four dominant species: Phycis blennoides Galeus melastomus, Plesionika martia and Todarodes sagittatus. These four species showed a wide range of distribution, appearing within all clusters. However, they showed a decreasing eastward pattern in their average biomass, except Plesionika martia, which showed an opposite trend (Supplementary material Table S1). The analysis also showed the presence of species exclusive to certain zones: for example, Nezumia sclerorhyncus appeared only in Corsica and the central basin (clusters GL-COR and CENTM), while Nezumia aequalis was limited to the western Mediterranean (ALBO, ESPN and BAL-GL). Aristaeomorpha foliacea contributed greatly to the community of the central Mediterranean (CENTM), while Aristeus antennatus did so to the community of northern Spain and the Balearic Islands (BAL-GL and ESPN) and Nephrops norvegicus to the cluster of the Gulf of Lions and Corsica (GL-COR). Six species of Chondrichthyes (G. melastomus, Galeus atlanticus, Etmopterus spinax, Chimaera monstrosa, Dipturus oxyrinchus and Scyliorhinus canicula) contributed to the overall similarity within groups. G. melastomus and S. canicula showed the highest percentage of contribution in the northern Spain cluster (ESPN). In contrast, E. spinax contributed to all clustersexcept ESPN. G. atlanticus contributed exclusively to ALBO, and D. oxyrinchus exclusively toGL-COR. Additionally, one species of Cnidaria, Isidella elongata,and the Gastropoda Galeodea rugosa contributed to the overall similarity within GL-COR and ALBO, respectively.
Based on generalized cross-validation and Aikake information criterion values, the best models explained 85.2% of the deviance variability of total biomass, including the additive effect of latitude and longitude (Table 4). The GAM models identified a significant statistical decrease in biomass eastward, but with minimum values in the longitude corresponding to GSA6 (i.e. –1° to 3.16°, Fig. 6B). Longitude was the single variable which explained the highest variability for biomass. These results were observed for the overall community as well as for all species groups analysed (Table 4). For latitude, the highest values seem to be related to higher latitude, although there was also a maximum around 36-37°, around which both low and high values were found (Fig. 6B), as they correspond to different trends found in the Alboran area (GSAs 1 and 2) and the Strait of Sicily (GSA 16).
Generalized additive models (GAM) fitted to assess the effect of latitude and longitude on biomass index for the overall sampled community (Total) and by species group. GCV, Generalized cross-validation; AIC, Akaike Information Criterion; Num, Number of samples; %Dev, percentage of explained deviance. Ostei, Osteichthyes; Chon, Chondrichthyes; Crus, Crustacea; Cepha, Cephalopoda; Cnid, Cnidaria; Echi, Echinodermata; OtherMol, other Mollusca (non-cephalopod).
Species group
Response variable
Factor
GCV
AIC
Num
%Dev
p
Total
LogBiom (g km–2)
~ lat
0.29
854
534
24.4
<0.001
LogBiom (g km–2)
~ long
0.08
214
534
77.2
<0.001
LogBiom (g km–2)
~ (lat, long )
0.06
17
534
85.2
<0.001
Cnid
LogBiom (g km–2)
~ lat
0.92
816
296
23.9
<0.001
LogBiom (g km–2)
~ long
0.51
643
296
57.8
<0.001
LogBiom (g km–2)
~ (lat, long )
0.45
603
296
66.5
<0.001
Cepha
LogBiom (g km–2)
~ lat
0.46
984
476
14.2
<0.001
LogBiom (g km–2)
~ long
0.34
848
476
36.2
<0.001
LogBiom (g km–2)
~ (lat, long )
0.31
801
476
45.5
<0.001
Crus
LogBiom (g km–2)
~ lat
0.17
568
534
13.8
<0.001
LogBiom (g km–2)
~ long
0.11
346
534
43.6
<0.001
LogBiom (g km–2)
~ (lat, long )
0.08
213
534
58.7
<0.001
Echi
LogBiom (g km–2)
~ lat
0.71
592
237
16.3
<0.001
LogBiom (g km–2)
~ long
0.36
429
237
58.7
<0.001
LogBiom (g km–2)
~ (lat, long )
0.34
417
237
64.5
<0.001
Chon
LogBiom (g km–2)
~ lat
0.44
1072
532
21.2
<0.001
LogBiom (g km–2)
~ long
0.18
623
532
66.3
<0.001
LogBiom (g km–2)
~ (lat, long )
0.15
512
532
74.4
<0.001
OtherMoll
LogBiom (g km–2)
~ lat
0.64
510
212
5.08
<0.05
LogBiom (g km–2)
~ long
0.21
277
212
69.9
<0.001
LogBiom (g km–2)
~ (lat, long )
0.2
268
212
74.2
<0.001
Oste
LogBiom (g km–2)
~ lat
0.64
993
534
21.3
<0.001
LogBiom (g km–2)
~ long
0.21
455
534
71.5
<0.001
LogBiom (g km–2)
~ (lat, long )
0.2
310
534
79.8
<0.001
Results of significant regression models fitted for overall community to assess the connection of log-transformed biomass (g km–2) index with latitude (A), longitude (B; negative values: West, positive values: East), and with the interaction of latitude and longitude (C). Shadowed areas represent 95% confidence intervals. The colour gradient of the geographical effect codes from low values (darker colours) to high values (lighter colours). The statistical summary of all models for the entire community and for different groups is summarized in Table 4.
Regarding the assessment of the possible effect of fishing effort on biomass values, the best models explained 90.7% of the deviance variability of total biomass, including latitude, longitude and fishing effort (Table 5). Spatial distribution explained higher variability than fishing effort. However, a significant decreasing pattern between biomass index and fishing effort was found, explaining 16.4% of the variability (Table 5, Fig 7A). The spatial effect in biomass and fishing effort showed opposite patterns, with the highest biomass in the areas with the lowest fishing effort (GSA 1), and the lowest biomass overlapped with the highest fishing effort (north of GSA 6) (Fig. 7B, C). An exception was found in GSA2, which was relatively highly impacted but showed high biomass values.
Summary of regression models results assessing the effects of fishing effort (FE), latitude and longitude on mean biomass values for the Spanish area.
Response variable
Factor
GCV
AIC
Num
%Dev
p
LogBiom (g km–2)
~ lat
0.11
46.3
68
84.5
<0.001
LogBiom (g km–2)
~ long
0.09
34.5
68
86.7
<0.001
LogBiom (g km–2)
~ FE
0.53
150.6
68
16.4
<0.001
LogBiom (g km–2)
~ s(FE) + s(lat, long )
0.06
10.4
68
90.7
<0.001
Partial effects of fishing effort (FE), and spatial distribution on the biomass index (g km–2) in the western Mediterranean (GSAs 1, 2, 5 and 6) using logarithmic transformation of the biomass index (A). The colour gradient of the geographical effect codes from low values (darker colours) to high values (lighter colours) (B). The statistical summary of plots is summarized in Table 5.
DISCUSSION
The present study explores long-term distribution patterns of megafauna in the Mediterranean Basin and how it may be affected by fishing. The results, overall, showed an eastward decreasing pattern in the abundance of megafauna assemblages collected with a trawl net of 500 to 800 m on the Mediterranean continental slope. This finding was observed for all the analysed taxonomic groups (Osteichthyes, Chondrichthyes, Crustacea, Echinodermata, Cephalopoda, other Mollusca, Tunicata and Porifera).
The decreasing eastward pattern observed is consistent with the available information limited to a single studied taxon (Cartes et al. 2002, Company et al. 2004, D’Onghia et al. 2004) and with results reported for deeper areas (below 800 m depth) (Danovaro et al. 2010, Tecchio et al. 2011a), and is primarily related to a decreasing west-east gradient of biological production (Coll et al. 2010). The Mediterranean Sea is characterized by a well-known west to east decrease in nutrient concentration, phytoplankton biomass and surface primary production (Moutin and Raimbault 2002) that seems to be mirrored by a similar decrease in benthic megafaunal biomass. The contrasting results in community descriptors obtained between the western and central Mediterranean basins in this study suggest that differences in biological production (or more generally, in food availability) between the two basins might be the strongest driver for community abundance and biomass. However, other drivers such as fishing effort may also modulate the pattern.
We observed that communities from northern Spain (GSA 6) do not follow the general eastward decreasing trend, showing different compositions and lower abundance values in all the studied taxonomic groups in comparison with adjacent areas. It is well known that the deep NW Mediterranean Basin is an area of high productivity in comparison with the other parts of the Mediterranean Sea, because of higher surface productivity, coastal inputs and exposure to massive inputs of organic matter caused by cyclic cascading events (Margalef 1985, Canals et al. 2006, 2009). However, deep-water bottoms in this area have been intensively trawled during the last century (Gorelli et al. 2016, Pitcher et al. 2017), and thus the presence of the lowest biomass values in these communities inhabiting more productive waters suggest that fishing exploitation might have an important impact on communities in this area. In fact, our results showed that deep trawl-fishing exploitation in GSA 6 is two to four times higher than that in adjacent areas of the Spanish continental slope margin, with the exception of Alboran Island (Fig. 7B). These results suggest that in this area human activity is an important factor in describing species community and partially modulates community descriptors. This is particularly true on the continental margins, where the life-history characteristics of several deep-sea species make them potentially more vulnerable to human impacts. These results are in accordance with those of previous studies showing that trawl-fishing exploitation in GSA 6 is very intense and much higher, for instance, than in the Balearic Islands (GSA 5) and the Alboran Sea (GSA 1), two sites that were also found to have healthier communities (Ramírez-Amaro et al. 2016, Quetglas et al. 2017). Following the present study trends but in a whole Mediterranean context, Colloca et al. (2017) identified GSA 6 (with GSA 9 and 17) as the fisheries with the lowest ecosystem sustainability. Moreover, Company et al. (2004) in a comparative study of Mediterranean decapods, found the same general eastward decreasing pattern in species dwelling below 1000 m depth, but not in species from the middle slope (between 600 and 800 m depth), also potentially due to the intense fishing activity targeting A. antennatus (western Mediterranean) and A. foliacea (in Greek waters). Note that in the present study the results of the fishing pressure in deep-water communities were skewed by the high effort values found at Alboran Island: in fact, if we deleted the eight trawls conducted in this area, the fishing effort explained 34% of the biomass variability. Alboran Island is a seamount of volcanic origin characterized by high productivity (Maldonado and Comas 1992) that has a special fishing regulation (BOE 1998). While it is becoming increasingly clear that human activities are impacting deep-sea Mediterranean communities (Cartes et al. 2004, Navarro et al. 2015, Piroddi et al. 2017), in order to set up effective mitigation strategies, the effects of fishing impacts need to be better quantified. Thus, these results are still preliminary and an analysis of a wider spatial scale will be conducted.
The longitudinal gradient found in this study for biomass and abundance was also observed in the assemblage composition. The MDS results showed samples clearly grouped among geographic sub-areas, except again for samples from northern Spain (GSA 6). The cluster analysis showed different species in the community composition between the western and central Mediterranean (the most significant examples were the two species of Nezumia spp. and Aristeids: N. aequalis was exclusively found in the western basin and N. sclerorhynchus in the central basin of Mediterranean; A. antennatus dominated the Balearic Sea, while A. foliacea was more abundant in Adriatic and Ionian Sea communities). These species could be preferentially distributed due to their biology and fishing could also potentially alter the community compositions in highly exploited areas (e.g. GSA 6). Additionally, we found higher homogeneity between samples from the central basin than those from the western basin. The more complex oceanographic dynamism of the western Mediterranean Basin (Béranger et al. 2005) and the human-induced changes may have driven the observed results and may indicate that each sub-area is characterized by the presence of a specific assemblage with high biogeographic complexity.
An increase in small and fast-growing species has been observed to be a general community response to trawling (Hiddink et al. 2006). The dominance of crustaceans (which are show the highest percentage of the community) throughout the Mediterranean middle slope suggests that they may be more competitive in a overfishing context than other megabenthos groups (e.g. elasmobranchs and non-crustacean invertebrates). High values of crustacean density have been related to medium and high levels of fishing effort (Guijarro et al. 2011) and, in fact, high catches of A. antennatus have been related to the ability of fishermen to remove competitors at a higher rate (Sardà and Maynou 1998). Interestingly, the only exception was the Alboran Basin (GSA 1 and 2) where the community was dominated by Chondrichthyes. The dominance of species characterized by slow population growth, such as sharks and rays, might be due to the fact that that the open slope has remained almost unexploited at depths lower than 500 m (Torres et al. 2001, Rey et al. 2005, Ramírez-Amaro et al. 2016), in addition to the particular oceanographic conditions, river discharges and the presence of submarine canyon which result in very productive areas (Millot 1999, Moranta et al. 2007). Moreover, sessile fauna has been suggested to be the first to undergo the negative effects of fishing (Jennings and Kaiser 1998, Roberts 2002). In accordance, the abundance and biomass of invertebrate non-crustaceans was much lower than that of more mobile fauna, such as fish and crustacean decapods, in the whole Mediterranean continental slope, except for punctual peaks in the values of Echinodermata (GSA 7) and Cnidaria (GSAs 8 and 18). The high densities and biomasses of these two taxonomic groups were caused by isolated hauls conducted over Leptometra sp. grounds and deep coral facies of Isidella elongata. The compact mud facies with Isidella elongata have been listed by the GFCM as sensitive habitats because of their importance as an essential habitat for certain crustacean species (Aristeus antennatus and Aristaeomorpha foliacea) (GFCM 2009, Fabri et al. 2014). Recent studies have shown that I. elongata is only present in areas where the fishing effort is low or absent (Maynou and Cartes 2012, Lauria et al. 2017), a finding that may explain the isolated presence of this species in the present study. This apparent impact of fishing on the Mediterranean continental slope is consistent with earlier descriptions of deep-sea communities in some areas of the Mediterranean (e.g. the Balearic Islands, Farriols et al. 2017) and recent analyses of sediment cores from the NW Mediterranean continental slope (Cartes et al. 2017), which found significant quantities of deep-water coral aggregates that have almost disappeared, indicating their previous/historic presence in these areas. However, further analysis, including the combined use of different gears (e.g. OTMS and the Agassiz dredge or benthic sledge) would be required in order to provide an integrated image of the megafaunal compartment, particularly regarding non-crustacean invertebrates (i.e. echinoderms), for which little information is currently available on their distribution (Tecchio et al. 2011b, Mecho et al. 2015).
In conclusion, the present study is the first attempt to describe the middle slope community along a wide spatial scale. It is in agreement with previous works which report a general eastward decreasing pattern in community descriptors (e.g. Coll et al. 2010, Keller et al. 2016) potentially related to environmental drivers (primary production and temperature). In addition, fishing activity also potentially modulates the deep Mediterranean community. This statement is particularly reflected in highly impacted areas, such as the northwestern Spanish middle slope, where the community does not and may be no longer follow the general eastward decreasing pattern, suggesting that the human activity has largely modified the community structure in this area.
Further investigations would require high-spatial resolution data (i.e. data from VMS or the Automatic Information System coupled to information on gear impact on the seabed) to determine the potential impact of fishing on a large spatial scales (Eigaard et al. 2015, 2017). Other alternative approaches need to be developed in order to assess the impact on ecosystem functioning and the response to fishing pressure and other human stressors. These could be based on life-history and biological traits of impacted species (de Juan and Demestre 2012, Rijnsdorp et al. 2016) or on community functional diversity indicators (Mouillot et al. 2013, Beauchard et al. 2017). They could be based on modelling approaches that estimate the reduction in community biomass or the shift in community composition corresponding to the estimated fishing intensity relative to system carrying capacity and natural disturbances (Pitcher et al. 2017). Additionally, the use of habitat suitability modelling, which predicts the suitability of a location for a species or group of species, based on their observed relationship with environmental conditions (Davies et al. 2015, Rowden et al. 2017), could help in spatially managing the impact of fishing on Mediterranean continental margins.
ACKNOWLEDGEMENTS
This paper was prepared in the framework of the MEDITS survey programme within the Data Collection Framework (DCF). The European Union and the Member States involved in the DCF are thankfully acknowledged. Fernandez-Arcaya was funded by a post-doctoral grant co-funded by the Regional Government of the Balearic Islands and the European Social Fund 2014–2020. The present study could not have been done without the projects supporting this research (ECLIPSAME; CTM2012-37701; CLIFISH: CTM2015-66400-C3-1-R MINECO/FEDER) and the FPI Fellowship granted to Maria Teresa Farriols (BES-2013-065112), both funded by the Spanish Ministry of Economy and Competitiveness. Sandrine Vaz and Angélique Jadaud received additional financial support from EU DG Environment funded the IDEM project (Implementation of the MSFD to the Deep Mediterranean Sea) under grant agreement 11.0661/2017/750680/SUB/ENV.C2.
REFERENCESAllenR.M.MetaxasA.SnelgroveP.V.Applying movement ecology to marine animals with complex life cycles2018101942https://doi.org/10.1146/annurev-marine-121916-063134AzovY.Eastern Mediterranean a marine desert?199123225232https://doi.org/10.1016/0025-326X(91)90679-MBeauchardO.VeríssimoH.QueirósA.M.The use of multiple biological traits in marine community ecology and its potential in ecological indicator development2017768196https://doi.org/10.1016/j.ecolind.2017.01.011BérangerK.MortierL.CréponM.Seasonal variability of water transport through the Straits of Gibraltar, Sicily and Corsica, derived from a high-resolution model of the Mediterranean circulation200566341364https://doi.org/10.1016/j.pocean.2004.07.013BertrandJ.LeonoriI.DremièreP.Y.Depth trajectory and performance of a trawl used for an international bottom trawl survey in the Mediterranean200266Suppl. 2169182BertrandJ.A.Gil de SolaL.PapaconstantinouC.The general specifications of the MEDITS surveys200266917Boletín Oficial del Estado (BOE)22319983254132541BurnhamK.P.AndersonD.R.2003Springer ScienceBusiness MediaCanalsM.PuigP.de MadronX.D.Flushing submarine canyons2006444354357https://doi.org/10.1038/nature05271CanalsM.DanovaroR.HeussnerS.Cascades in Mediterranean submarine grand canyons2009222643https://doi.org/10.5670/oceanog.2009.03CartesJ.AbellóP.LlorisD.Analysis of feeding guilds of fish and decapod crustaceans during the MEDITS-99 cruise along the Iberian Peninsula Mediterranean coasts200266209220https://doi.org/10.3989/scimar.2002.66s2209CartesJ.MaynouF.SardàF.The Mediterranean deep-sea ecosystems: an overview of their diversity, structure, functioning and anthropogenic impactsPart 12004WWFIUCN938CartesJ.MaynouF.FanelliE.Long-term changes in the composition and diversity of deep-slope megabenthos and trophic webs off Catalonia (western Mediterranean): Are trends related to climatic oscillations?2009823246https://doi.org/10.1016/j.pocean.2009.03.003CartesJ.E.SchironeA.BarsantiM.Recent reconstruction of deep-water macrofaunal communities recorded in Continental Margin sediments in the Balearic Basin20171255264https://doi.org/10.1016/j.dsr.2017.04.016ClarkM.R.AlthausF.SchlacherT.A.The impacts of deep-sea fisheries on benthic communities: a review201673i51i69https://doi.org/10.1093/icesjms/fsv123ClarkeK.GorleyR.2006PlymouthCollM.PiroddiC.SteenbeekJ.The Biodiversity of the Mediterranean Sea: Estimates, Patterns, and Threats20105e11842https://doi.org/10.1371/journal.pone.0011842CollM.PiroddiC.AlbouyC.The Mediterranean Sea under siege: spatial overlap between marine biodiversity, cumulative threats and marine reserves201221465480https://doi.org/10.1111/j.1466-8238.2011.00697.xCollocaF.ScarcellaG.LibralatoS.Recent Trends and Impacts of Fisheries Exploitation on Mediterranean Stocks and Ecosystems20174244https://doi.org/10.3389/fmars.2017.00244CompanyJ.B.MaioranoP.PlaityW.Deep-sea decapod crustaceans in the western and central Mediterranean Sea: preliminary aspects of species distribution, biomass and population structure200468Suppl. 37386https://doi.org/10.3989/scimar.2004.68s373D'OnghiaG.PolitouC.Y.BozzanoA.Deep-water fish assemblages in the Mediterranean Sea.200468Suppl. 38799https://doi.org/10.3989/scimar.2004.68s387DanovaroR.CompanyJ.B.CorinaldesiC.Deep-Sea Biodiversity in the Mediterranean Sea: The Known, the Unknown, and the Unknowable20105e11832https://doi.org/10.1371/journal.pone.0011832DaviesJ.S.StewartH.A.NarayanaswamyB.E.Benthic Assemblages of the Anton Dohrn Seamount (NE Atlantic): Defining Deep-Sea Biotopes to Support Habitat Mapping and Management Efforts with a Focus on Vulnerable Marine Ecosystems201510e0124815https://doi.org/10.1371/journal.pone.0124815de JuanS.DemestreM.A Trawl Disturbance Indicator to quantify large scale fishing impact on benthic ecosystems201218183190https://doi.org/10.1016/j.ecolind.2011.11.020de JuanS.LleonartJ.A conceptual framework for the protection of vulnerable habitats impacted by fishing activities in the Mediterranean high seas201053717723https://doi.org/10.1016/j.ocecoaman.2010.10.005EigaardO.R.BastardieF.BreenM.Estimating seabed pressure from demersal trawls, seines, and dredges based on gear design and dimensions201573i27i43https://doi.org/10.1093/icesjms/fsv099EigaardO.R.BastardieF.HintzenN.T.The footprint of bottom trawling in European waters: distribution, intensity, and seabed integrity201774847865FabriM.C.PedelL.BeuckL.Megafauna of vulnerable marine ecosystems in French mediterranean submarine canyons: Spatial distribution and anthropogenic impacts2014104184207https://doi.org/10.1016/j.dsr2.2013.06.016FarriolsM.T.OrdinesF.SomerfieldP.J.Bottom trawl impacts on Mediterranean demersal fish diversity: Not so obvious or are we too late?201713784102https://doi.org/10.1016/j.csr.2016.11.011FarrugioH.OliverP.BiagiF.An overview of the history, knowledge, recent and future research trends in Mediterranean fisheries199357105119Fernandez-ArcayaU.DrazenJ.C.MuruaH.Bathymetric gradients of fecundity and egg size in fishes: A Mediterranean case study2016116106117https://doi.org/10.1016/j.dsr.2016.08.005Fernandez-ArcayaU.Ramirez-LlodraE.AguzziJ.Ecological role of submarine canyons and need for canyon conservation: a review201745https://doi.org/10.3389/fmars.2017.00005FollesaM.C.PorcuC.GastoniA.Community structure of bathyal decapod crustaceans off South‐Eastern Sardinian deep‐waters (Central‐Western Mediterranean)200930s1188199https://doi.org/10.1111/j.1439-0485.2009.00323.xFoveauA.VazS.DesroyN.Process-driven and biological characterisation and mapping of seabed habitats sensitive to trawling201712e0184486https://doi.org/10.1371/journal.pone.0184486GalilB.Taking stock: inventory of alien species in the Mediterranean Sea200911359372https://doi.org/10.1007/s10530-008-9253-yGeneral Fisheries Commission for the Mediterranean (GFCM), SAC32009Málaga3GorelliG.SardàF.CompanyJ.B.Fishing Effort Increase and Resource Status of the Deep-Sea Red Shrimp Aristeus antennatus (Risso 1816) in the Northwest Mediterranean Sea Since the 1950s.201624192202https://doi.org/10.1080/23308249.2015.1119799GoriA.OrejasC.MadurellT.Bathymetrical distribution and size structure of cold-water coral populations in the Cap de Creus and Lacaze-Duthiers canyons (northwestern Mediterranean)2013102049https://doi.org/10.5194/bg-10-2049-2013GrangerV.FromentinJ.-M.BezN.Large-scale spatio-temporal monitoring highlights hotspots of demersal fish diversity in the Mediterranean Sea20151306574https://doi.org/10.1016/j.pocean.2014.10.002GuijarroB.TserpesG.MorantaJ.Assessment of the deep water trawl fishery off the Balearic Islands (western Mediterranean): from single to multi-species approach20116706785https://doi.org/10.1007/s10750-011-0670-zGuijarroB.OrdinesF.MassutíE.Improving the ecological efficiency of the bottom trawl fishery in the Western Mediterranean: It's about time!201783204214https://doi.org/10.1016/j.marpol.2017.06.007HiddinkJ.HuttonT.JenningsS.Predicting the effects of area closures and fishing effort restrictions on the production, biomass, and species richness of benthic invertebrate communities200663822830https://doi.org/10.1016/j.icesjms.2006.02.006JenningsS.KaiserM.J.The effects of fishing on marine ecosystems199834201352https://doi.org/10.1016/S0065-2881(08)60212-6KellerS.BartolinoV.HidalgoM.Large-scale spatio-temporal patterns of Mediterranean cephalopod diversity201611e0146469https://doi.org/10.1371/journal.pone.0146469LauriaV.GarofaloG.FiorentinoF.Species distribution models of two critically endangered deep-sea octocorals reveal fishing impacts on vulnerable marine ecosystems in central Mediterranean Sea201778049https://doi.org/10.1038/s41598-017-08386-zLevinL.A.DaytonP.K.Ecological theory and continental margins: where shallow meets deep200924606617https://doi.org/10.1016/j.tree.2009.04.012MaldonadoA.ComasM.Geology and geophysics of the Alboran Sea: An introduction1992126165https://doi.org/10.1007/BF02084913MargalefR.1985OxfordPergamon Press363MassutíE.ReñonesO.Demersal resource assemblages in the trawl fishing grounds off the Balearic Islands (western Mediterranean)200569167181https://doi.org/10.3989/scimar.2005.69n1167MaynouF.CartesJ.E.Effects of trawling on fish and invertebrates from deep-sea coral facies of Isidella elongata in the western Mediterranean20129215011507https://doi.org/10.1017/S0025315411001603MechoA.BillettD.S.Ramírez-LlodraE.First records, rediscovery and compilation of deep-sea echinoderms in the middle and lower continental slope of the Mediterranean Sea201478281302https://doi.org/10.3989/scimar.03983.30CMechoA.Fernandez-ArcayaU.AguzziJ.Reproductive biology of the seastar Ceramaster grenadensis from the deep north-western Mediterranean Sea201595805815https://doi.org/10.1017/S0025315415000065MerrettN.R.HaedrichR.L.1997LondonChapman and HallMillotC.Circulation in the Western Mediterranean Sea199920423442https://doi.org/10.1016/S0924-7963(98)00078-5MorantaJ.MassutíE.PalmerM.Geographic and bathymetric trends in abundance, biomass and body size of four grenadier fishes along the Iberian coast in the western Mediterranean2007726383https://doi.org/10.1016/j.pocean.2006.09.003MorantaJ.QuetglasA.MassutíE.Spatio-temporal variations in deep-sea demersal communities off the Balearic Islands (western Mediterranean)200871346366https://doi.org/10.1016/j.jmarsys.2007.02.029MouillotD.GrahamN.A.J.VillégerS.A functional approach reveals community responses to disturbances201328167177https://doi.org/10.1016/j.tree.2012.10.004MoutinT.RaimbaultP.Primary production, carbon export and nutrients availability in western and eastern Mediterranean Sea in early summer 1996 (MINOS cruise)200233273288https://doi.org/10.1016/S0924-7963(02)00062-3NavarroJ.CollM.CardadorL.The relative roles of the environment, human activities and spatial factors in the spatial distribution of marine biodiversity in the Western Mediterranean Sea2015131126137https://doi.org/10.1016/j.pocean.2014.12.004NavarroJ.CardadorL.FernándezÁ.M.Differences in the relative roles of environment, prey availability and human activity in the spatial distribution of two marine mesopredators living in highly exploited ecosystems201643440450https://doi.org/10.1111/jbi.12648NorseE.A.BrookeS.CheungW.W.Sustainability of deep-sea fisheries201236307320https://doi.org/10.1016/j.marpol.2011.06.008PhamC.K.Ramirez-LlodraE.AltC.H.Marine litter distribution and density in European seas, from the shelves to deep basins20149e95839https://doi.org/10.1371/journal.pone.0095839PiroddiC.CollM.LiqueteC.Historical changes of the Mediterranean Sea ecosystem: modelling the role and impact of primary productivity and fisheries changes over time2017744491https://doi.org/10.1038/srep44491PitcherC.R.EllisN.JenningsS.Estimating the sustainability of towed fishing-gear impacts on seabed habitats: a simple quantitative risk assessment method applicable to data-limited fisheries20178472480https://doi.org/10.1111/2041-210X.12705QuetglasA.MerinoG.GonzálezJ.Harvest strategies for an ecosystem approach to fisheries management in western Mediterranean demersal fisheries20174106Ramírez-AmaroS.OrdinesF.TerrasaB.Demersal chondrichthyans in the western Mediterranean: assemblages and biological parameters of their main species201667636652https://doi.org/10.1071/MF15093Ramirez-LlodraE.BrandtA.DanovaroR.Deep, diverse and definitely different: unique attributes of the world's largest ecosystem2010728512899https://doi.org/10.5194/bg-7-2851-2010Ramirez-LlodraE.De MolB.CompanyJ.B.Effects of natural and anthropogenic processes in the distribution of marine litter in the deep Mediterranean Sea2013118273287https://doi.org/10.1016/j.pocean.2013.07.027ReliniG.Orsi ReliniL.The decline of blue red shrimps stocks in the gulf of Genoa198751245260ReyJ.Gil de SolaL.MassutíE.Distribution and biology of the blackmouth catshark Galeus melastomus in the Alboran Sea (Southwestern Mediterranean)200535215223RijnsdorpA.D.BastardieF.BolamS.G.Towards a framework for the quantitative assessment of trawling impact on the seabed and benthic ecosystem201673i127i138https://doi.org/10.1093/icesjms/fsv207RobertsC.M.Deep impact: the rising toll of fishing in the deep sea200217242245https://doi.org/10.1016/S0169-5347(02)02492-8RowdenA.A.AndersonO.F.GeorgianS.E.High-Resolution Habitat Suitability Models for the Conservation and Management of Vulnerable Marine Ecosystems on the Louisville Seamount Chain, South Pacific Ocean20174335https://doi.org/10.3389/fmars.2017.00335SardàF.DemestreM.Estudio biológico de la gamba Aristeus antennatus (Risso, 1816) en el Mar Catalán (NE de España)198751213232SardàF.MaynouF.Assessing perceptions: do Catalan fishermen catch more shrimp on Fridays?199836149157https://doi.org/10.1016/S0165-7836(98)00102-7SardàF.D'OnghiaG.PolitouC.Y.Deep-sea distribution, biological and ecological aspects of Aristeus antennatus (Risso, 1816) in the western and central Mediterranean Sea200468117127https://doi.org/10.3989/scimar.2004.68s3117SardàF.CalafatA.FlexasM.M.An introduction to Mediterranean deep-sea biology200468738https://doi.org/10.3989/scimar.2004.68s37SpedicatoM.T.GrecoS.LemboG.Prime valutazioni sulla struttura dello stock di Aristeus antennatus (Risso, 1816) nel Tirreno Centro Meridionale19952239244TecchioS.Ramírez-LlodraE.SardàF.Biodiversity of deep-sea demersal megafauna in western and central Mediterranean basins201175341350https://doi.org/10.3989/scimar.201175n2341TecchioS.Ramírez-LlodraE.SardàF.Drivers of deep Mediterranean megabenthos communities along longitudinal and bathymetric gradients2011439181192https://doi.org/10.3354/meps09333TorresP.GonzálezM.ReyJ.Rose shrimp fishery’s associated fauna in not exploited grounds on the Alboran Sea slope (Western Mediterranean Sea)200136330330VielminiI.PerryA.L.CornaxM.J.Untying the Mediterranean Gordian Knot: A Twenty First Century Challenge for Fisheries Management20174195WormB.BarbierE.B.BeaumontN.Impacts of Biodiversity Loss on Ocean Ecosystem Services2006314787790https://doi.org/10.1126/science.1132294WatsonR.A.MoratoT.Fishing down the deep: Accounting for within-species changes in depth of fishing20131406365https://doi.org/10.1016/j.fishres.2012.12.004WoodS.N.On confidence intervals for generalized additive models based on penalized regression splines200648445464https://doi.org/10.1111/j.1467-842X.2006.00450.xSupplementary MateriAL
The following supplementary material is available through the online version of this article and at the following link: http://scimar.icm.csic.es/scimar/supplm/sm04852esm.pdf
Table S1. – Results of similarity percentage analysis (SIMPER) of the overall community on each GSA. Av.Abund (average abundance), Av.Sim (average similarity), Sim/SD (ratio of the average similarity and the standard deviation), Contrib% (percentage contribution), Cum% (cumulative percentages).