1. Introduction
The red gorgonian
Paramuricea clavata (Risso, 1827) (Cnidaria: Anthozoa) is among the main habitat-forming species in the Mediterranean circalittoral zone, and its aggregations are known to heavily influence the diversity and structure of the associated assemblage [
1,
2,
3].
P. clavata forests constitute noteworthy seascapes attracting broad fluxes of SCUBA diving tourism, thus resulting in an important economic resource [
4,
5]. Moreover, they host significant populations of commercial fish such as the common dentex
Dentex dentex (Linnaeus, 1758) and the dusky grouper
Epinephelus marginatus (Lowe, 1834), typical targets of artisanal fishing [
6].
P. clavata forests can be impacted by different kinds of pressures, mainly related to global warming and anthropogenic activities; therefore, the species is considered “vulnerable” according to the Red List of the International Union for Conservation of Nature (IUCN), and the reduction of some Mediterranean populations has been assessed in the last three decades [
7]. The first pressure is related to the recent increase of positive thermal anomalies causing cascade effects such as the development of pathogenic bacteria and blooms of filamentous algae and mucilage. In turn, these phenomena affect the gorgonians with extensive diseases that favour the settling of invertebrate epibionts on the damaged colonies [
8,
9,
10,
11,
12,
13,
14,
15,
16]. Marine Heat Waves (MHWs) also result in Mass Mortality Events (MMEs) that have been well-described for many organisms of the coralligenous assemblage [
11,
17,
18,
19,
20,
21,
22,
23], but mainly hit the red gorgonians [
24,
25,
26,
27].
In addition to climate-related impacts, it is well known that also boat anchoring and demersal recreational and artisanal fishing activities threaten
P. clavata due to its large size, branched shape, and skeleton limited flexibility that favour the entanglement of fishing gear [
20,
28,
29]. Bavestrello et al. [
28] described the abrasive action by lost lines in contact with the gorgonian coenenchyme, producing suitable conditions for the settlement of several epibiotic organisms that ultimately cause the total or partial breaking of the colony under the current flow. In this regard, several authors employed the percentage of epibionted colonies as a biological index quantifying the fishing-induced stress of a gorgonian population [
24,
30,
31,
32,
33,
34]. Variable levels of epibiosis, hence the health status of the populations, have been reported in different Mediterranean sites. A study conducted in the Medes Islands (Catalan Sea) indicated that 10–33% of the colonies settled in unprotected areas were partially colonised by epibionts, whereas only 4–10% of the populations inside the borders of the Marine Protected Area (MPA) were involved in the phenomenon [
35]. Along the eastern Adriatic coasts, in an area characterised by low fishing pressure, the colonies with denuded or epibionted branches were less than 10% [
36]. On the contrary, inside the Portofino MPA (Ligurian Sea) and in the Tavolara - Punta Coda Cavallo MPA (Tyrrhenian Sea) colonies with epibiosis exceeded 50% of the total, suggesting a high fishing impact [
6,
37,
38].
Monitoring programs on a broad spatial and temporal scale dedicated to the evaluation of the health status of gorgonian forests are lacking, and the recovery ability of populations after mass mortality episodes was followed only for a few years [
39,
40,
41,
42]. The Marine Strategy Framework Directive (MSFD 2008/56/EC) represents the EU’s Integrated Maritime Policy tool to achieve the Good Environmental Status (GES) of marine waters, based on eleven qualitative descriptors. Within these, the descriptor “biodiversity” states that “
The quality and occurrence of habitats and the distribution and abundance of species are in line with prevailing physiographic, geographic and climatic conditions (MSFD, 2008/56/EC, Annex I). Ecosystem components – groups of species and habitat types, such as the coralligenous – are recognised as indicators of environmental quality, becoming priority habitats target of study and monitoring activities”.
Thanks to the data obtained during the Marine Strategy project (2015–2022), the Ligurian Sea (NW Mediterranean Sea) is nowadays one of the best-known areas along Italian coasts [
43]. During this project, 80 sites were explored (2015–2018) and then subjected to monitoring in the following years. The forests of
P. clavata were widely distributed along the whole Ligurian arc, occurring on sloping outcroppings and coralligenous rocks up to 89 m. The forests dominated by this species represented over 13% of the sampling units considered in the study and were characterised by an average density of 7.4 ± 3.0 colonies m
−2 with maximum values reaching 19 colonies m
−2. This gorgonian also consistently occurred in other communities, such as that dominated by
Corallium rubrum (Linnaeus, 1758) [
43].
Thanks to the large ROV footage collected during the Ligurian Marine Strategy surveys, this paper aims to give a comprehensive view of the level of epibiosis affecting the P. clavata forests along the whole Ligurian coastline and describe a pluriannual trajectory of their health status evolution.
4. Discussion
The diseases hitting the gorgonian populations in the Mediterranean Sea are generally attributed to the acute stress produced by positive thermal anomalies that have progressively increased in frequency and intensity in the last decades [
22,
42]. During the period of our observations, MHW events were recorded in 2015, 2016, 2018, and 2022 [
22,
46].
Our study revealed that the epibiosis affecting
Paramuricea clavata in three macro-areas along the Ligurian coast had an overlapping temporal trajectory that can be generalised at the regional scale. The percentage of involved colonies strongly increased from 2015 until 2019, while in the following period (2019–2022), the level remained unvaried or decreased. This pattern agrees with the trend of MHWs in the western Mediterranean Sea. In fact, the significant sharp increase between 2015 and 2019 could be related to the intense MHWs that took place between 2015 and 2018. Iborra et al. [
42] investigated the gorgonian population of the Gulf of Calvi (NW Corsica Island) over a 15-year period (2004, 2014, and 2019) and the observed changes in the level of necrosis of the colonies were related to the trend of MHWs occurring in that period. According to these authors, the effect of MHWs on the gorgonian integrity reached relatively deep coastal waters (down to 40 m), where temperature increases are generally not recorded [
42]. The shallow water temperature rise induces a strong stratification of the water column with a consequent low availability of trophic resources and a high respiratory demand [
47,
48]. This hypothesis agrees with the bathymetric distribution of the level of epibiosis observed during our study in the Ligurian Sea, where the maximum values were recorded below 50 m. In the same area, previous diseases involving cnidarians and sponges were already observed up to 70 m [
17,
49,
50,
51]. At the same time, the bathymetric distribution of the maximum level of epibiosis and entanglement (50–70 m) can be also related to the bathymetric distribution of
P. clavata forests in the Ligurian Sea, being mainly concentrated between 45 and 65 m [
43].
A remarkable quickness in the evolution of the levels of epibiosis emerges from the present study. In the eastern Ligurian Sea, in four years (2015–2019), the percentage of epibionted colonies increased more than three times, and in the central area increased six times. This evidence supports the damages inflicted by mass diseases following MHWs. In fact, during the 2003 episode, almost 80% of the studied colonies in the Gulf of Genoa developed epibiosis in a few weeks [
21]. A long-term monitoring conducted on the
P. clavata populations in the period 2003–2017 in the Scandola MPA highlighted a substantial decline in density and a slight loss of mean colony biomass after 15 years, and all populations were farther from recovery in 2017–2018 than in 2008. On average, in 2003, values decreased by around 71% and 80%, despite no significant changes in the size structure were observed in any population immediately after the 2003 MHW [
41]. On the other hand, Cupido et al. [
50] observed a strong recovery of the populations in the La Spezia Gulf after the late-summer 1999 and 2003 large mortality events due to thermal stress. The long-term observation (1998 to 2008) confirmed that a positive net recruitment and canopy reestablishment could start some years after a high-mortality event, also under naturally stressed environmental conditions.
Despite this evidence, our data strongly suggest that also human activities are deeply involved in the development of epibiosis in
P. clavata. The correlation, at the regional scale, between the level of epibiosis with the number of fishing vessels and harbours evidences the role of fishing activities. Moreover, the overlapping of the bathymetric trend of epibiosis with that of the entanglement indicates a strong correlation also with this kind of impact. The depth range of epibionted colonies overlaps that of fishing with settled gear targeting spiny lobsters, sparids (common dentex, gilthead seabreams), and groupers [
52]. Finally, the influence of fishing on the level of epibiosis was also supported by the temporal correlation of the two parameters in all the studied areas. In some cases (e.g., A3) it appears evident that a reduction in the entangled colonies anticipates that of epibionted ones.
These data strongly enforce the idea that the observed variations in the Ligurian forests are due to a synergistic effect of natural and anthropogenic causes. Probably, the damages caused by fishing activities pile up on colonies already deeply stressed by thermal diseases. A further indication is that, in the period 2018–2022, despite the occurrence of one of the strongest MHW ever observed [
46], the epibiosis remained stable or decreased concomitantly with a period of strong reduction in the percentage of entanglement supporting a decrease in the fishing effort.
Unfortunately, data about the annual trend of artisanal fishing and, particularly, recreative fishing are virtually impossible to obtain, and therefore, only a hypothesis can be formulated. We suggest that the reduction of the fishing activity speculated for the period 2019–22 could be due to the COVID-19 pandemic lockdown imposed from March to May 2020. Fisheries were limited by the coronavirus pandemic at a global scale. In that period, the lockdown, followed by a strong shrinking in seafood requests, determined a decrease in fishing activities. For example, a recent study demonstrated a reduction of about 50% of the fishing effort for the Adriatic Sea [
53]. In the eastern Mediterranean Sea, from December 2019 to February 2020, the average monthly gross margin for fishermen was 2.5 times less than the usual average [
54]. The fishing effort in an exploited area of the north-western Mediterranean coast (Spain) during the lockdown dropped by 34%, landings were down by 49%, and revenues declined by 39% in comparison with the same period in 2017–2019 [
55]. This scenario well represents also the conditions of the Ligurian fishing economy in that period, an economy based on traditional activities carried out by small artisanal communities, mainly targeting high-value resources [
56].
In addition, during the lockdown period, the media reported the presence of iconic large marine animals, such as marine mammals, elasmobranchs and marine turtles, in unexpected areas, such as very coastal areas or harbours [
55]. Despite the wide interest in the effects of the COVID-19 pandemic on marine habitats and fisheries, no data are available about benthic organisms: our observation could be the first suggestion of an improvement in the health status of
P. clavata forests resulting from a short-term strong reduction of fishing activities. The vulnerability of habitat-forming species is influenced by their low resilience to mechanical impacts, driven by modest to slow growth rates [
57]. This work supports the potential effectiveness of Fisheries Restricted Areas (FRAs) on these complex habitats.
Author Contributions
F Conceptualization, G.B. and M.C.; methodology, M.C., A.D., F.E., M.T.; formal analysis and investigation, M.C., F.E., M.T.; resources, R.B.; data curation, M.C., F.E., M.T.; writing—original draft preparation, G.B., F.B., M.B., M.C., A.D., F.E., M.T.; writing—review and editing, G.B., F.B., M.B., M.C., A.D., F.E., M.T.; visualization, F.B; supervision, G.B., M.B.; project administration, G.B.; funding acquisition, R.B. All authors have read and agreed to the published version of the manuscript.
Figure 2.
(a) Average density (± SE); (b) average height (± SE), and (c) average percentage of epibionted colonies (± SE) in the three macro-areas.
Figure 2.
(a) Average density (± SE); (b) average height (± SE), and (c) average percentage of epibionted colonies (± SE) in the three macro-areas.
Figure 3.
Correlation of the percentage of epibionted colonies in the three macro-areas with the number of fishing harbours ((a), r = 0.98) and fishing boats ((b), r = 0.75) insisting in each macro-area.
Figure 3.
Correlation of the percentage of epibionted colonies in the three macro-areas with the number of fishing harbours ((a), r = 0.98) and fishing boats ((b), r = 0.75) insisting in each macro-area.
Figure 4.
Bathymetric distribution of the percentage (a) of epibionted colonies and (b) of colonies entangled in lost fishing gear.
Figure 4.
Bathymetric distribution of the percentage (a) of epibionted colonies and (b) of colonies entangled in lost fishing gear.
Figure 5.
Eastern macro-area (A1). Differences across time (Before, white bars/After, grey bars) obtained comparing data of (a) epibiosis and (b) entanglement in the same groups of transects replicated after 3–4 years. Inset: Percentage of transects that presented a positive (epibiosis/entanglement increasing, blue bars) or negative (epibiosis/entanglement decreasing, red bars) variation in the value of affected colonies in each considered period.
Figure 5.
Eastern macro-area (A1). Differences across time (Before, white bars/After, grey bars) obtained comparing data of (a) epibiosis and (b) entanglement in the same groups of transects replicated after 3–4 years. Inset: Percentage of transects that presented a positive (epibiosis/entanglement increasing, blue bars) or negative (epibiosis/entanglement decreasing, red bars) variation in the value of affected colonies in each considered period.
Figure 6.
Central macro-area (A2). Differences across time (Before, white bars/After, grey bars) obtained comparing data of (a) epibiosis and (b) entanglement in the same groups of transects replicated after 3–4 years. Inset: Percentage of transects that presented a positive (epibiosis/entanglement increasing, blue bars) or negative (epibiosis/entanglement decreasing, red bars) variation in the value of affected colonies in each considered period.
Figure 6.
Central macro-area (A2). Differences across time (Before, white bars/After, grey bars) obtained comparing data of (a) epibiosis and (b) entanglement in the same groups of transects replicated after 3–4 years. Inset: Percentage of transects that presented a positive (epibiosis/entanglement increasing, blue bars) or negative (epibiosis/entanglement decreasing, red bars) variation in the value of affected colonies in each considered period.
Figure 7.
Western macro-area (A3). Differences across time (Before, white bars/After, grey bars) obtained comparing data of (a) epibiosis and (b) entanglement in the same groups of transects replicated after 3–4 years. Inset: Percentage of transects that presented a positive (epibiosis/entanglement increasing, blue bars) or negative (epibiosis/entanglement decreasing, red bars) variation in the value of affected colonies in each considered period.
Figure 7.
Western macro-area (A3). Differences across time (Before, white bars/After, grey bars) obtained comparing data of (a) epibiosis and (b) entanglement in the same groups of transects replicated after 3–4 years. Inset: Percentage of transects that presented a positive (epibiosis/entanglement increasing, blue bars) or negative (epibiosis/entanglement decreasing, red bars) variation in the value of affected colonies in each considered period.
Table 2.
Results of ANOVA and pair-wise comparisons of Paramuricea clavata, density, colony height, and percentages of epibionted specimens in each macro-area. Bray-Curtis similarity index used for the resemblance matrix construction; permutation n = 9999. Significant effects are in bold.
Table 2.
Results of ANOVA and pair-wise comparisons of Paramuricea clavata, density, colony height, and percentages of epibionted specimens in each macro-area. Bray-Curtis similarity index used for the resemblance matrix construction; permutation n = 9999. Significant effects are in bold.
|
df |
SS |
MS |
Pseudo-F |
P (perm) |
Pair-wises |
T |
P (perm) |
Density |
|
|
|
|
|
|
|
|
Macro-area |
2 |
5201.1 |
2600.5 |
2.716 |
0.3686 |
|
|
|
Res |
107 |
1.0245E+05 |
957.5 |
|
|
|
|
|
Total |
|
|
|
|
|
|
|
|
Height |
|
|
|
|
|
Height |
|
|
Macro-area |
2 |
1360.9 |
680.44 |
35.551 |
0.0267 |
A1 vs A2 |
20.462 |
0.0357 |
Res |
107 |
20480 |
191.4 |
|
|
A1 vs A3 |
0.70001 |
0.5049 |
Total |
109 |
21840 |
|
|
|
A2 vs A3 |
22.188 |
0.0249 |
Epibiosis |
|
|
|
|
|
Epibiosis |
|
|
Macro-area |
2 |
11824 |
5912.2 |
38.355 |
0.006 |
A1 vs A2 |
12.312 |
0.1907 |
Res |
107 |
1.65E+09 |
1541.4 |
|
|
A1 vs A3 |
2.721 |
0.0012 |
Total |
109 |
1.77E+09 |
|
|
|
A2 vs A3 |
16.747 |
0.0545 |
Entanglement |
|
|
|
|
|
Entanglement |
|
|
Macro-area |
2 |
25589 |
12795 |
68.677 |
0.0001 |
A1 vs A2 |
37.941 |
0.0001 |
Res |
107 |
1.99E+09 |
1863 |
|
|
A1 vs A3 |
19.837 |
0.0119 |
Total |
109 |
2.25E+09 |
|
|
|
A2 vs A3 |
16.131 |
0.0553 |
Table 3.
Results of PERMANOVA and pair-wise comparisons for the temporal analysis (Before/After) of percentage of epibiosis and entanglement of Paramuricea clavata in each replicated transect. Bray-Curtis similarity index used for the resemblance matrix construction; permutation n = 9999. Significant effects are in bold.
Table 3.
Results of PERMANOVA and pair-wise comparisons for the temporal analysis (Before/After) of percentage of epibiosis and entanglement of Paramuricea clavata in each replicated transect. Bray-Curtis similarity index used for the resemblance matrix construction; permutation n = 9999. Significant effects are in bold.
|
df |
SS |
MS |
Pseudo-F |
P (perm) |
Pair-wises |
t |
P (perm) |
Unique perms |
P (MC) |
A1 |
|
|
|
|
|
|
|
|
|
|
Epibiosis |
|
|
|
|
|
|
|
|
|
|
Before/After |
5 |
11612 |
2322.3 |
29.303 |
0.0228 |
2015/19 |
40.632 |
0.027 |
35 |
0.0024 |
Res |
20 |
15850 |
792.52 |
|
|
2016/19 |
20.984 |
0.0532 |
35 |
0.0529 |
Total |
25 |
27462 |
|
|
|
2019/22 |
0.79759 |
0.502 |
126 |
0.4923 |
Entanglement |
|
|
|
|
|
|
|
|
|
|
Before/After |
5 |
18518 |
3703.6 |
21.411 |
0.0244 |
2015/19 |
27.952 |
0.0259 |
35 |
0.006 |
Res |
20 |
34595 |
1729.7 |
|
|
2016/19 |
0.4955 |
0.7255 |
35 |
0.7313 |
Total |
25 |
53113 |
|
|
|
2019/22 |
14.721 |
0.0794 |
126 |
0.1193 |
A2 |
|
|
|
|
|
|
|
|
|
|
Epibiosis |
|
|
|
|
|
|
|
|
|
|
Before/After |
5 |
23565 |
4713 |
46.103 |
0.001 |
2015/18 |
27.982 |
0.0135 |
462 |
0.0043 |
Res |
32 |
32713 |
1022.3 |
|
|
2016/20 |
38.243 |
0.0287 |
35 |
0.0029 |
Total |
37 |
56278 |
|
|
|
2018/22 |
0.516 |
0.709 |
8170 |
0.7162 |
Entanglement |
|
|
|
|
|
|
|
|
|
|
Before/After |
5 |
15630 |
3125.9 |
32.804 |
0.0031 |
2015/18 |
2.362 |
0.0203 |
461 |
0.0157 |
Res |
32 |
30493 |
952.92 |
|
|
2016/20 |
0.59803 |
0.8289 |
35 |
0.6893 |
Total |
37 |
46123 |
|
|
|
2018/22 |
30.371 |
0.0046 |
8150 |
0.0043 |
A3 |
|
|
|
|
|
|
|
|
|
|
Epibiosis |
|
|
|
|
|
|
|
|
|
|
Before/After |
5 |
9829.4 |
1965.9 |
12.736 |
0.2349 |
2015/18 |
18.318 |
0.0536 |
35 |
0.0779 |
Res |
22 |
33959 |
1543.6 |
|
|
2016/20 |
11.444 |
0.2483 |
91 |
0.2815 |
Total |
27 |
43788 |
|
|
|
2018/21 |
0.83113 |
0.6803 |
126 |
0.5607 |
Entanglement |
|
|
|
|
|
|
|
|
|
|
Before/After |
5 |
17345 |
3469 |
19.534 |
0.0293 |
2015/18 |
21.685 |
0.0855 |
35 |
0.0435 |
Res |
22 |
39068 |
1775.8 |
|
|
2016/20 |
15.939 |
0.1094 |
91 |
0.1148 |
Total |
27 |
56413 |
|
|
|
2018/21 |
15.329 |
0.0403 |
91 |
0.0872 |