Next Article in Journal
Hydrogeochemical Characteristics and Human Health Risk Assessment of Fluoride Enrichment in Water in Faulted Basins of Qinghai-Tibet Plateau—A Case Study of Sanhe Plain in Guide Basin
Next Article in Special Issue
An Assessment of Potential Beam Trawling Impact on North-Western Black Sea Benthic Habitats Aiming at a Sustainable Fisheries Management
Previous Article in Journal
How Important for Society Is Recreation Provided by Multi-Purpose Water Reservoirs? Welfare Analysis of the Vltava River Reservoir System
Previous Article in Special Issue
Troubles Never Come Alone: Outcome of Multiple Pressures on a Temperate Rocky Reef
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Benthic Fish Communities Associated with Posidonia oceanica Beds May Reveal the Fishing Impact and Effectiveness of Marine Protected Areas: Two Case Studies in the Southern Tyrrhenian Sea

1
Consorzio Nazionale Interuniversitario per le Scienze del Mare (CoNISMa) Rome, URL-Centro Direzionale—Isola C4, 80143 Napoli, Italy
2
Department of Science and Technology (DiST), Marine Ecology Laboratory, Parthenope University of Naples, Centro Direzionale—Isola C4, 80143 Napoli, Italy
*
Author to whom correspondence should be addressed.
Water 2023, 15(10), 1967; https://doi.org/10.3390/w15101967
Submission received: 5 May 2023 / Revised: 16 May 2023 / Accepted: 20 May 2023 / Published: 22 May 2023
(This article belongs to the Special Issue Anthropogenic Impacts on Benthic Marine Ecosystems)

Abstract

:
The effectiveness of a Marine Protected Area (MPA) is strongly related to its zoning design, management, and surveillance, and fish communities are considered indicators of conservation effects. Posidonia oceanica beds and fish assemblages were studied here to evaluate the fishing impacts and the effectiveness of the Marine Protected Areas (MPAs) of Santa Maria di Castellabate (SMC) and Costa degli Infreschi e della Masseta (CIM) in the South Tyrrhenian Sea, Italy. P. oceanica characteristics were estimated counting and collecting shoots, while fish communities were investigated by visual census. Multivariate analyses were performed to detect differences in fish assemblages among the MPA protection levels and seagrass characteristics. Significant differences in fish communities were observed between the depth and protection levels. The latter are probably due to the zoning design at SMC, which does not allow a spillover effect, and to fishing activities, as showed by the ABC curve results. Trophic analyses suggested the use of longlines at SMC, due to the dominance of planktivorous fish species, and of trammel nets at CIM, due to the great heterogeneity of the fish trophic groups. Shoot density and leaf length were significantly correlated with fish communities. CCA analyses showed omnivorous species, mostly belonging to the Sparidae family, correlated with the highest shoot densities, with their feeding habits mostly directed to crustaceans, while macro- and microinvertivores, mainly belonging to the Serranidae and Labridae families, correlated with the highest leaf lengths, where they sough shelter and feed. In conclusion, the data suggest that a zoning redesign at SMC and an increase in the surveillance of fishing activities at CIM might enhance the effectiveness of these MPAs, highlighting the role of fish community studies in identifying MPA efficiency issues and providing trustful guidelines for their management.

1. Introduction

Marine Protected Areas (MPAs) are management tools aimed at conserving biodiversity, concurrently with an increase in the local economies in a sustainable development scenario. Their success depends on management and enforcement, which are strongly influenced by their social acceptability [1]. Controversies about MPA effectiveness still exist, with MPAs being often exclusively considered as coercive tools established for conservation goals [2]. To smooth out these controversies, the Ecosystem-Based Management approach [3] points out the need to report, by means of investigations, the natural resources enforcement within MPAs [4], recognizable by reserve and spillover effects [5,6,7,8,9,10]. MPA success is primary determined by four issues: (i) engagement level of local populations, (ii) socioeconomic features of the institution areas, (iii) zoning design, and (iv) ecological factors [11]. Leaving aside the first two topics, which are not addressed in the present study, zoning design and ecological factors are the main issues over which researchers, decision makers, and stakeholders argue [1,12]. In particular, zoning design is one of the topics about which stakeholders are more concerned, because it directly influences local population economies. Too often, MPAs have been established only considering population needs or conservation goals, overlooking the meeting point between protection policies and stakeholders that, on the contrary, should be at the base of their establishment process [13,14]. Ecological factors are concretized in the list of priority habitats (Habitats directive 92/43/CEE, Barcelona Convention, SPA/RAC), recognized as those providing ecosystem goods and services for human health and the economy’s well-being [15].
In the Mediterranean Sea, Posidonia oceanica (L.) Delile beds represent the most important and most-studied vegetated habitat [16,17], shaping extensive and extremely productive seagrass beds with high structural complexity [18,19]. P. oceanica leaves can reach over 1 m in height during summer and 40 cm in winter. Usually, shoot density is maximal in shallow waters, where it may attain more than 1000 shoots per m2, and decreases exponentially with depth [20,21]. Nevertheless, edaphic characteristics, such as strong long-shore currents and high water turbidity, can determine highest shoot densities in intermediate areas of the beds, where environmental characteristics are more suitable for their growth [22,23]. Beds can extend from shallow waters to 40 m depth, although in very clear water, they can occur in up to 50 m depth [24]. These habitats play a key role in carbon sequestration and stock, oxygenation of the oceans, prevention of coastal erosion, and provision of food and shelter to many invertebrate and fish species [25]. Indeed, seagrasses are recognized as settlement, spawning, nursery, and feeding areas for many species belonging to all trophic levels, from primary producers, epiphyting P. oceanica leaves and rhizomes, to top predators, finding prey during their different life stages [25].
Fish communities quickly respond to habitat alterations due to anthropic activities [26] and protection actions [7]. Spatial and temporal variations in their structure are influenced by both biological (e.g., predation, competition, larval dynamics, spawning, and recruitment rates) and environmental (e.g., depth, light, nutrient availability, habitat structure, bed complexity, and current) factors [27], as well as by anthropogenic disturbances (e.g., fisheries) [28,29]. Thus, for example, in bordering rock–algal/P. oceanica beds, the Labridae species was observed in higher abundances when compared to other seascape contexts; the zooplanktivorous Oblada melanurus (Linnaeus, 1758) is replaced by Spicara smaris (Linnaeus, 1758), and the species Chromis chromis (Linnaeus, 1758) is less abundant than in mosaic beds [30,31]. Nevertheless, P. oceanica beds are habitats for many fish species, also of commercial interest [32], forming a community structured according to their behaviour. In particular, four behaviours can be recognized: (i) open-water erratic species belonging to the planktivorous Alosidae and piscivorous Carangidae families, (ii) open-water sedentary species Chromis chromis (Pomacentridae) that live on P. oceanica beds, but eat plankton in the water column, far from the canopy, (iii) mesophilic nektobenthic species, primary belonging to the Sparidae family, that live near the canopy and are mainly composed of omnivorous individuals, and (iv) sedentary nektobenthic species, mainly belonging to the Searranidae and Labridae families, that live among the leaves and are mainly composed of macro- and microinvertivore individuals [26,33,34].
Although there are several studies on fish communities and their spatiotemporal variations in P. oceanica beds within Marine Protected Areas [35,36,37,38,39], there is a lack of information on the relationship between P. oceanica bed characteristics and fish community structures. The present study performed at the Santa Maria di Castellabate and Costa degli Infreschi e della Masseta MPAs aims at (i) assessing the relationship between P. oceanica bed characteristics and the associated fish communities and (ii) evaluating fishing impacts and the effectiveness of MPAs as conservation tools through variations in fish assemblages at different MPA protection levels. Finally, possible advice to enhance conservation success and local economies are also provided.

2. Materials and Methods

2.1. Study Area

The MPAs of Santa Maria di Castellabate (SMC) and Costa degli Infreschi e della Masseta (CIM) were established by the 21 October 2009 Ministerial Decree and stretch along the northern and southern parts, respectively, of Cilento Coast (South Tyrrhenian Sea), a vast mountain and coastal area in the Campania region (Italy) [40,41]. They are divided into three protection zones that regulate human activities: the fully protected (A) zone, where any activity that could cause damage or disturbance to the environment is forbidden; the general reserve (B) zone, where some activities are allowed, but regulated and authorized by the managing authority; and the partial reserve (C) zone, where fruition and sustainable use of the sea, even of moderate environmental impact, are allowed.
SMC (Figure 1a) includes 19 km of coastline. The “A” zone extends for 1.7 km2, and it is isolated at the northernmost part of the MPA; the “B” zone is divided in a small area (6.35 km2), surrounding the “A” zone, and in a very large area in front of Punta Licosa (21.35 km2); and the “C” zone includes beaches and deepest areas (37.89 km2). The SMC sea bottoms are characterized by lush macroalgal forests of Fucales in shallow waters [42] and coralligenous bioconstructions in deeper waters. The large acrochorus of Punta Licosa, that gently degrades more than 3 miles offshore, is completely covered by a vast P. oceanica bed, for more than 26 km2, and almost completely included within the “B” zone [43].
CIM (Figure 2b) includes 13 km of coastline. The A zone extends for 0.4 km2, and it is located at the MPA core; the B zone surrounds the A zone and extends for 4.79 km2; the C zone surrounds the B zone and extends for 18.36 km2. Sea bottoms are characterized by photophilic algal and coralligenous habitats at the shallowest and deepest areas, respectively. The steep slope of the sea bottom allows for the settlement of P. oceanica beds (1.41 km2) for a few tens of metres offshore [44].

2.2. Sampling Activities and Data Collection

Sampling was carried out between September and October 2020 on P. oceanica beds at SMC and CIM. In particular, at each MPA, 4 stations were selected for each protection level (“A”, “B”, and “C”) at two different depths: the seagrass bed upper limit (SL) and 15 m depth (15 m) (Figure 1). The present investigation aims at emphasizing differences on fish communities among protection levels and their relations with P. oceanica bed characteristics at each MPA.
P. oceanica beds were sampled following the standardized Marine Strategy Framework Directive protocol [45]. At each station, 6 randomly placed 40 × 40 cm squares were considered to estimate the bed shoot density (shoots/m2). Moreover, 6 shoots were collected to estimate the following phenological parameters: average leaf number per shoot (considering only intermediate and adult leaves), leaf width, green and brown leaf lengths, leaf base length, percentage of intact or broken (by mechanical and biological factors) leaf apices, and total leaf surface.
Total leaf surface T s (m2)/m2 was calculated as in Equation (1):
T s = n ° L ¯ l ¯ w ¯ n ° S
where n ° L ¯ , l ¯ , w ¯ are the average number of leaves, average lengths, and average widths, respectively, and n ° S is the shoot density. Juvenile leaves were excluded because they presented neither photosynthetic nor brown parts.
At each station, fish communities were estimated on 125 m2 by visual census technique replicated 3 times, as described by [5]. Species were classified according to their value as commercial (C), low commercial (LC), and noncommercial (NC), while biomasses were estimated using the population size–structure and length–weight relationships following [46].
Species richness (SR: total species number/125 m2), density (D: total abundance/125 m2), Shannon diversity index (H’), Pielou evenness index (J), and total biomass (B: total kg/125 m2) were also estimated.
Trophic group (TG) of each fish species was assessed following trophic level intervals (TL) defined as [47] pure herbivores (2 ≤ TL ≤ 2.1); omnivores with a preference for vegetable materials (2.1 < TL ≤ 2.9); omnivores with a preference for animals (2.9 < TL ≤ 3.7); and carnivores with a preference for large decapods, cephalopods, and fish (3.7 < TL ≤ 4.5). Finally, the mean sample trophic level T L ¯ was also computed by Equation (2)
T L ¯ = i = 1 n Y i T L i i = 1 n Y i
where Y i is the biomass in 125 m2 of species I; T L i is the trophic level of the species I; and n is the total number of species in the sample [48,49].

2.3. Data Analyses

Distance-based permutation multivariate analyses of variance (PERMANOVAs, [50]) were carried out on P. oceanica data, fish community densities, and biomasses to test for differences within each MPA and among protection levels and different depths. The experimental design involved the “protection level” factor (PL, fixed, three levels: A, B, and C zones), the “depth” factor (De, fixed, two levels: SL and 15 m), and the “station” factor (St, nested in PL and De, 12 levels). Before the analyses, P. oceanica data were standardized, and fish densities and biomasses were square-rooted to reduce the weight of very abundant species. Multivariate analyses were based on Bray–Curtis similarity, and each term was tested by 4999 random permutations [51,52]. Post hoc pairwise comparisons were conducted using the PERMANOVA t-statistic with 4999 permutations. Multivariate patterns were visualized through Canonical Analyses of Principal coordinates (CAP) of PLxDe elements [53]. ABC curves [54] were employed to compare abundance–biomass relationships of commercial species within each MPA protection level. Distance-based linear modelling (DistLM), using stepwise selection procedure and adjusted R2 as selection criterion [55], was performed to find correlations among fish densities, fish biomasses, and P. oceanica bed characteristics, considered as independent variables. Three plot Canonical Correlation Analyses (CCAs, [56]) were also performed to assess correlation among PLxDe elements, trophic levels, and P. oceanica characteristics that significantly affect fish communities resulting from DistLM.
Finally, univariate PERMANOVAs based on Euclidean distance [57] were conducted on species richness, total density, diversity index, evenness, and total biomass of the whole assemblage. The experimental design for these analyses was the same described above for multivariate analyses.

3. Results

3.1. Analysis of Posidonia oceanica Beds

In Table 1, the values of the physiographic parameters detected within MPAs and characterizing P. oceanica beds are reported.
At SMC, the shoot density was 551.5 ± 129.0 shoots/m2 and the leaf length was 39.5 ± 17.18 cm (brown part = 16.5%), with the highest values at the 15 m sites. At CIM, the shoot density was 421.4 ± 93.46 shoots/m2, with the highest values at the SL sites, and the leaf length was 45.5 ± 21.6 cm (brown part = 16.0%), with the highest values at the 15 m sites (Figure 2a).
At SMC, the leaf surface was 10.9 ± 6.8 m2/m2, with the highest values at the 15 m sites, and the mean percentages of the mechanical and biological broken apices were 17.5% and 2.9%, respectively. At CIM, the leaf surface was 10.0 ± 6.6 m2/m2, with the highest values at the SL sites, and the mechanical and biological broken apices were 20.7% and 3.7%, respectively (Figure 2b).

3.2. Descriptive Analysis of Total Assemblages

A total of 36 species of fishes belonging to 11 families were detected at SMC and CIM. In Table 2, commercial value (CV) and trophic group (TG) of each species are reported.
In Table 3, the values of the synecological indices detected within MPAs and characterizing fish communities are reported.
At SMC, the species richness was 8.0 ± 2.1 species/125 m2, with the highest values at the SL sites, and the fish density was 110.1 ± 52.4 individuals/125 m2, with the highest values at the 15 m sites. At CIM, the species richness was 7.5 ± 1.7 species/125 m2 and the fish density was 105.0 ± 51.4 individuals/125 m2, with the highest values at the SL sites (Figure 3a).
At SMC, the Shannon diversity index was 2.0 ± 0.5 and the Pielou evenness index was 0.68 ± 0.1, with the highest values at the SL sites. At CIM, the Shannon diversity index was 2.0 ± 0.4 and the mean Pielou evenness index was 0.69 ± 0.1, with the highest values at the 15 m sites (Figure 3b).
At SMC, the fish biomass was 1.5 ± 0.9, with the highest values at the 15 m sites. At CIM, the fish biomass was 0.9 ± 0.3, with the highest values at the 15 m sites (Figure 3c).
At SMC, the dominant trophic group was the planktivorous one, and the corresponding average trophic level was 3.47 ± 0.17 at SL and 3.57 ± 0.07 at 15 m. At CIM, the trophic groups were very heterogeneous, with a nonevident dominant one; the average trophic level was 3.04 ± 0.17 at SL and 3.49 ± 0.21 at 15 m (Figure 3d).
No significant differences were detected using synecological parameters, except for biomasses. They were different between MPAs (p < 0.001). Within MPAs, at SMC, there was an interaction of PLxDe (p = 0.0084) due to differences between the C zone and the other ones at 15 m (pA-C = 0.0028 and pB-C = 0.0476) and between SL and 15 m (p = 0.0084) within the C zone. At CIM, differences were detected only for the factor PL (p = 0.006) between the B zone and the other ones (pA-B = 0.015 and pB-C = 0.0412).

3.3. Multivariate Analyses on Fish Assemblages

At SMC, significant differences were detected for the PL, due to differences between the “A” and “C” zones, and the De factors. At CIM, significant differences were detected for the interaction of the PLxDe factors. In particular, differences were detected between SL and 15 m in the “A” and “C” zones; differences were also detected between the “A” and “B” zones at SL, while at 15 m, differences were detected between the “C” zone and the other ones for the density data and between the “A” zone and the other ones for the biomass data (Table 4).
At SMC, the CAP analysis on density data showed protection degrees and depth gradients (Figure 4a). Most of the “A” zone elements were distributed on the positive part of the CAP1 axis; the “B” zone elements were located around the 0 value, while the “C” zone elements were polarized on the negative part. Most of the SL elements were distributed on the positive part of the CAP2 axis, while most of the 15 m elements, relating to Serranus scriba (Linnaeus, 1758), Coris julis (Linnaeus, 1758), S. smaris, and Symphodus tinca (Linnaeus, 1758), were distributed on the negative part. A CAP analysis on biomass data (Figure 4b) did not highlight evident gradients, although most of the “A” and “B” zone elements were distributed on the negative part of the CAP1 axis, while most of the “C” zone elements were on the positive part. Most of the SL elements were distributed on the positive part of the CAP2 axis, while the 15 m elements, relating to C. julis, Spondyliosoma cantharus (Linnaeus, 1758), Diplodus sargus (Linnaeus, 1758), S. smaris, and Diplodus vulgaris (Geoffroy Saint-Hilaire, 1817) were distributed on the negative part.
At CIM, the CAP analysis on density data showed protection levels and depth gradients (Figure 4c). “A” zone elements were polarized at the negative part of the CAP 1 axis; “B” zone elements were polarized around the 0 value, while “C” zone elements were distributed on the positive part. Most of the SL elements were distributed on the negative part of the CAP2 axis, relating to D. sargus and D. vulgaris. Most of the 15 m elements were distributed on the positive part of the CAP2 axis, relating to S. cabrilla, C. julis, S. tinca, and S. smaris. A CAP analysis on biomass data showed clear protection levels and depth gradients (Figure 4d). “A” zone elements were polarized on the negative part of the CAP 1 axis, while “B” and “C” zones elements were polarized on the positive part. SL elements were polarized on the negative part of the CAP2 axis, relating to D. sargus and D. vulgaris, while most of the 15 m elements were polarized on the positive part, relating to S. smaris, S. tinca, C. julis, and C. chromis.
At SMC, the ABC curves (Figure 5a) resulted in a negative W statistic in the “A” zone (W = −0.013) and a positive W statistic in the “B” (W = 0.113) and “C” zones (W = 0.038). At CIM, the ABC curves (Figure 5b) resulted in a negative W statistic in the “A” (W = −0.02) and “C” zones (W = −0.106) and a positive W statistic in the “B” zone (W = 0.07).

3.4. Relations between P. oceanica Bed Characteristics and Fish Communities

At SMC, shoot densities and total leaf lengths were the P. oceanica bed characteristics significantly affecting fish communities, while at CIM, fish communities were significantly affected by green and brown leaf lengths and shoot densities (Table 5).
At SMC, the trophic groups strongly related to P. oceanica beds were macroinvertivores, microinvertivores, omnivorous, and planktivorous fishes. They are primary correlated with the 15 m sites that are, in turn, directly related to total leaf lengths and shoot densities (Figure 6a). At CIM, the trophic groups strongly related to P. oceanica beds were microinvertivores, omnivorous, planktivorous, and herbivores fishes. In particular, species were correlated with the 15 m sites, except for the omnivorous ones related to the SL sites of the B zone. The upper sites were correlated with shoot densities, while green and brown leaf lengths were mainly related to the 15 m sites (Figure 6b).

4. Discussion

MPAs are established to protect biodiversity and natural resources for future generations, stemming the decline of marine ecosystems and of the numerous goods and services they provide [59,60]. The urgent need to prove their effectiveness is mandatory in order to take prompt actions to enhance wildlife protection and local economy development. At the same time, it is necessary to know where possible actions, such as the zoning reshaping or the surveillance enhancing, should be implemented to increase the local economic yield (primarily in terms of production by fishing of target species), minimizing potential stakeholder objections. In this study, the effectiveness of the Santa Maria di Castellabate and Costa degli Infreschi e della Masseta MPAs was assessed to evaluate if new conservation actions are necessary. Moreover, the relationship between P. oceanica bed characteristics and fish assemblages was investigated to determine where zoning should be redesigned in order to enhance protection and local economies.

4.1. Posidonia oceanica Beds

P. oceanica beds can be classified as dense prairies [61] in an equilibrium state with environmental constraints [62,63]. In this study, the characteristics of P. oceanica beds were considered as environmental parameters structuring fish communities. Bed density ranged from ~350 to ~600 shoots/m2, with highest values at the 15 m sites of SMC and at the shallow sites of CIM.
Leaf length, both at SMC and CIM, was coherent with previous studies assessing that, and during autumn, average leaf length was ~30–50 cm [19,21], with larger lengths occurring at intermediate depths (e.g., 15 m). Higher lengths at the SMC and CIM intermediate depths corresponded to higher average leaf surfaces that provide a larger available substrate for sessile [64] and vagile [65] organisms. Indeed, leaves sustain many organisms that actively feed on them, eat their remains, or graze on their epiphytes. In turn, these organisms enter the food webs, transferring matter and energy up to the predators of the higher trophic levels [66].
The most frequent cause of broken leaf apices was mechanical damage rather than biological damage (i.e., grazing), also within the no-take (A) zone of both MPAs, where every human activity disturbing or damaging the environment should be forbidden. This result is presumably related to the anchoring activity carried out during summer [67].

4.2. Fish Communities

Most evidence on the MPAs’ effectiveness stems from quantitative assessments of fish populations [68], which rapidly and directly respond to protection instruments in terms of density, body size, and biomass [7]. Compared to other studies, the fish communities investigated here do not have a good status, as evident from the low density and biomass values [9,39,68]. Regarding the synecological parameters, differences in total biomass were detected between the two MPAs and among protection levels and depths. These differences were analysed in more detail considering the density and biomass of the different species, which also showed significant differences among protection levels and depths. Differences between depth levels are coherent with studies identifying different fish assemblages at different depths in P. oceanica beds [69,70]. Differences among protection levels could be related to two possible reasons: MPA zoning design and level of surveillance. In this framework, two well-known reviews written in the early 2000s reported the importance of the positioning of the fully protected (A) zone [2,71]. Indeed, due to its strict level of protection, the “A” zone is the most pristine zone of an MPA, providing refuge, feeding, and nursery areas for fish assemblages [4,12].
SMC is characterized by a fully protected (A) zone confined at the north edge of the MPA, surrounded by a small general reserve (B zone) and a large partial reserve (C zone) at the south. This zoning design does not allow possible spillover effects coming from the potential reserve effect [9]. Moreover, the “B” zone of Punta Licosa, which comprises most of the MPA area, includes the vast majority of the P. oceanica beds, with the most abundant fish populations [31]. Fish assemblage size structure showed a dominance of small- or medium-sized individuals. This result is typical of P. oceanica beds for their role as a nursery area [26], although it is worth pointing out that the lack of large top predators might have also been caused by local fishing [58,72].
At CIM, the fully protected “A” zone is located at the MPA core, surrounded by the general reserve (B zone), which is in turn surrounded by the partial reserve (C zone). Hence, CIM zoning design should allow reserve and spillover effects, i.e., set of MPA positive effects that occur within and outside the fully protected zone [9]. However, the significant interaction of the PLxDe factors, a typical result of these effects [7,69], is not consistent with the very low synecological and biomass values. Weak reserve and spillover effects can be found at the SL sites, as suggested by the biomass CAP analysis, characterized by higher species richness and densities, although with small size classes [38]. Conversely, the reserve and spillover effects completely disappear at the 15 m depth sites, where the very low biomass values suggest a significant negative impact, probably of anthropic origin, on fish communities [73,74,75,76,77]. These findings were also confirmed by the ABC curves. Indeed, the positive W statistic in the B zone may be due to the convergence of incoming individuals from the fully protected “A” zone to the general reserve (B zone) during the adult phase, and the consequent concentration of fishing harvesting here and within the partial reserve (C zone).
In undisturbed communities, the presence of large individuals of the target species results in higher biomasses, whose curves lie completely above the abundance ones, determining a positive W statistic. Conversely, in a highly disturbed community, mainly characterized by small individuals, the abundance curves lie above the biomass ones due to the presence of anthropic impacts, determining a negative W statistic [78,79,80]. SMC showed a negative W statistic only in the fully protected “A” zone, suggesting that despite the strictly protection level, this zone is in an area impacted by anthropic activities, hampering the reserve effect. The latter, on the contrary, is more evident in the other MPA zones (B and C), as showed by the positive W statistics. At CIM, the general reserve B zone is the only one showing a positive W statistic, suggesting that the fully protected “A” and partial reserve “C” zones are subjected to exploitation activities, possibly fishing.
These results are coherent with the trophic analyses that (i) confirm the fishing activities in both MPAs due to the almost complete absence of piscivorous species [58,81] and (ii) suggest the employment of specific and different fishing gear in both MPAs. The selectiveness of fishing gear affects the biomass and trophic structure of fish communities [82,83,84,85,86]. At SMC, the dominance of planktivorous species, such as Oblada melanurus (Linneus 1758), S. smaris, and Sardina pilchardus (Walbaum, 1792), probably depends on the use of selective gear, such as longlines, that act on the trophic structure removing fish predators [80]. At CIM, trophic structure heterogeneity coupled with very low biomasses suggests a larger use of nonselective fishing gear, such as trammel nets [87].

4.3. Relationship between P. oceanica Beds and Fish Communities

The present study aims to enhance the knowledge on the relation between fish communities and P. oceanica bed characteristics, which can influence biodiversity protection and local economies. A clear space distribution of fish species was evident looking at the CAP and CCA analyses. In particular, the water column was dominated by erratic and sedentary fish species, such as S. smaris and C. chromis, both planktivorous species. Nektobenthic fish species were, as expected, more related to P. oceanica rhizomes and leaves, where they seek refuge and feed. In particular, sites characterized by the highest shoot densities, such as those at 15 m in SMC and at SL in CIM, were characterized by omnivores species, such as D. vulgaris and D. sargus, feeding on the most abundant invertebrates in the Posidonia canopy, i.e., decapods and amphipods [64,88]. Moreover, juvenile stages of the Sparidae species were commonly present among P. oceanica shoots, where they can easily find shelter and food, confirming the role of P. oceanica as a nursery habitat [38]. Sites characterized by the highest leaf lengths, such as those at 15 m in both SMC and CIM, were related to seagrass resident species [89], mainly belonging to the Labridae and Serranidae families, such as C. julis, S. tinca, S. scriba, and S. cabrilla, that find shelter and food within the canopy throughout their entire life cycle.

5. Conclusions

The present study has provided useful knowledge about the role of P. oceanica beds in structuring fish communities and in estimating the effectiveness of an MPA. Low fish densities pointed out some issues at both the Santa Maria di Castellabate and Costa degli Infreschi e della Masseta MPAs, probably due to inappropriate zoning design and legal or illegal fishing pressure. It might be essential to reconsider the SMC zoning design, moving the fully protected “A” zone to the wide P. oceanica beds of Punta Licosa, currently included in the “B” zone, in order to adequately protect this key feeding and nursery area, with possible consequent spillover effects in the next years. This new, fully protected (A) zone should be established within the area characterized by the highest Posidonia shoot densities, due to the evident correlations with the Sparidae species, such as D. sargus and D. vulgaris, that are among the most valuable commercial species. As suggested by [80], the results confirm that MPA fish assemblages are threatened also by relevant fishing pressure, primarily performed by longlines at SMC, for the lack of piscivorous species, and trammel nets at CIM, for the high heterogeneity of trophic groups and very low biomasses. Therefore, to enhance conservation and, concurrently, to ensure the survival of local economies, an overall better management and an increased surveillance would be desirable to allow the constant development and maintenance of fish assemblages. These findings highlight the urgent need to carry out up-to-date studies on fish communities at a regional scale to identify possible MPA management issues. Moreover, the idea of a new MPA concept is recently emerging, based on flexible and changeable zoning that considers anthropogenic threats at both global (climatic change) and local (fishery) scales. These studies would also provide useful guidelines to establish new potential MPAs or to enhance the effectiveness of those already existing.

Author Contributions

Conceptualization, L.A.; Methodology, L.A. and F.R.; Software, L.A.; Validation, L.A., F.R. and L.D.; Formal Analysis, L.A.; Investigation, L.A., A.P., F.R. and A.C.d.M.; Resources, F.F. and G.F.R.; Data Curation, F.F. and R.S.; Writing—Original Draft Preparation, L.A.; Writing—Review and Editing, F.R., D.C., A.C.d.M., F.F., F.D.S., L.D. and R.S.; Visualization, R.S.; Supervision, G.F.R.; Project Administration, G.F.R.; Funding Acquisition, G.F.R. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the “Environmental Accounting in the Santa Maria di Castellabate Marine Protected Area” and the “Environmental Accounting in the Costa degli Infreschi e della Masseta Marine Protected Area” projects, financed by the Italian Ministry of the Ecological Transition.

Data Availability Statement

Not applicable.

Acknowledgments

Authors would like to thank the Santa Maria di Castellabate and Costa degli Infreschi e della Masseta Marine Protected Area Directors and Staff for their support and availability during this project.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Gall, S.C.; Rodwell, L.D. Evaluating the social acceptability of Marine Protected Areas. Mar. Policy 2016, 65, 30–38. [Google Scholar] [CrossRef]
  2. Francour, P.; Harmelin, J.G.; Pollard, D.; Sartoretto, S. A review of marine protected areas in the northwestern Mediterranean region: Siting, usage, zonation and management. Aquat. Conserv. Mar. Freshw. Ecosyst. 2001, 11, 155–188. [Google Scholar] [CrossRef]
  3. Katsanevakis, S.; Stelzenmüller, V.; South, A.; Sørensen, T.K.; Jones, P.J.S.; Kerr, S.; Badalamenti, F.; Anagnostou, C.; Breen, P.; Chust, G.; et al. Ecosystem-based marine spatial management: Review of concepts, policies, tools, and critical issues. Ocean Coast. Manag. 2011, 54, 807–820. [Google Scholar] [CrossRef]
  4. Grorud-colvert, K.; Sullivan-stack, J.; Roberts, C.; Constant, V.; Horta e Costa, B.; Pike, E.P.; Kingston, N.; Laffoley, D.; Sala, E.; Claudet, J.; et al. The MPA Guide: A framework to achieve global goals for the ocean. Science 2021, 373, eabf0861. [Google Scholar] [CrossRef] [PubMed]
  5. La Mesa, G.; Vacchi, M. An analysis of the coastal fish assemblage of the Ustica Island marine reserve (Mediterranean Sea). Mar. Ecol. 1999, 20, 147–165. [Google Scholar] [CrossRef]
  6. Valls, A.; Gascuel, D.; Guénette, S.; Francour, P. Modeling trophic interactions to assess the effects of a marine protected area: Case study in the NW Mediterranean Sea. Mar. Ecol. Prog. Ser. 2012, 456, 201–214. [Google Scholar] [CrossRef]
  7. Appolloni, L.; Bevilacqua, S.; Sbrescia, L.; Sandulli, R.; Terlizzi, A.; Russo, G.F. Does full protection count for the maintenance of β-diversity patterns in marine communities? Evidence from Mediterranean fish assemblages. Aquat. Conserv. Mar. Freshw. Ecosyst. 2017, 27, 828–838. [Google Scholar] [CrossRef]
  8. Giakoumi, S.; Scianna, C.; Plass-Johnson, J.; Micheli, F.; Grorud-Colvert, K.; Thiriet, P.; Claudet, J.; Di Carlo, G.; Di Franco, A.; Gaines, S.D.; et al. Ecological effects of full and partial protection in the crowded Mediterranean Sea: A regional meta-analysis. Sci. Rep. 2017, 7, 8940. [Google Scholar] [CrossRef]
  9. Harmelin-Vivien, M.; Le Diréach, L.; Bayle-Sempere, J.; Charbonnel, E.; García-Charton, J.A.; Ody, D.; Pérez-Ruzafa, A.; Reñones, O.; Sánchez-Jerez, P.; Valle, C. Gradients of abundance and biomass across reserve boundaries in six Mediterranean marine protected areas: Evidence of fish spillover? Biol. Conserv. 2008, 141, 1829–1839. [Google Scholar] [CrossRef]
  10. White, J.W.; Scholz, A.J.; Rassweiler, A.; Steinback, C.; Botsford, L.W.; Kruse, S.; Costello, C.; Mitarai, S.; Siegel, D.A.; Drake, P.T.; et al. A comparison of approaches used for economic analysis in marine protected area network planning in California. Ocean Coast. Manag. 2013, 74, 77–89. [Google Scholar] [CrossRef]
  11. Rossiter, J.S.; Levine, A. What makes a “successful” marine protected area? The unique context of Hawaii’s fish replenishment areas. Mar. Policy 2014, 44, 196–203. [Google Scholar] [CrossRef]
  12. Frazão Santos, C.; Ehler, C.N.; Agardy, T.; Andrade, F.; Orbach, M.K.; Crowder, L.B. Marine spatial planning. In World Seas: An Environmental Evaluation Volume III: Ecological Issues and Environmental Impacts; Sheppard, C., Ed.; Candice Janco; Academic Press: Cambridge, MA, USA, 2019; pp. 571–592. ISBN 9780128050521. [Google Scholar] [CrossRef]
  13. Jentoft, S.; Pascual-Fernandez, J.J.; de la Cruz Modino, R.; Gonzalez-Ramallal, M.; Chuenpagdee, R. What Stakeholders Think About Marine Protected Areas: Case Studies from Spain. Hum. Ecol. 2012, 40, 185–197. [Google Scholar] [CrossRef]
  14. Appolloni, L.; Sandulli, R.; Vetrano, G.; Russo, G.F. A new approach to assess marine opportunity costs and monetary values-in-use for spatial planning and conservation; the case study of Gulf of Naples, Mediterranean Sea, Italy. Ocean Coast. Manag. 2018, 152, 135–144. [Google Scholar] [CrossRef]
  15. Buonocore, E.; Grande, U.; Franzese, P.P.; Russo, G.F. Trends and evolution in the concept of marine ecosystem services: An overview. Water 2021, 13, 2060. [Google Scholar] [CrossRef]
  16. Buonocore, E.; Russo, G.F.; Franzese, P.P. Assessing natural capital value in the network of Italian marine protected areas: A comparative approach. Ecol. Quest. 2020, 31, 67–76. [Google Scholar] [CrossRef]
  17. Appolloni, L.; Buonocore, E.; Russo, G.F.; Franzese, P.P. The use of remote sensing for monitoring Posidonia oceanica and Marine Protected Areas: A systemic review. Ecol. Quest. 2020, 31, 7–17. [Google Scholar] [CrossRef]
  18. Duarte, C.M.; Chiscano, C.L. Seagrass biomass and production: A reassessment. Aquat. Bot. 1999, 65, 159–174. [Google Scholar] [CrossRef]
  19. Pergent-Martini, C.; Rico-Raimondino, V.; Pergent, G. Primary production of Posidonia oceanica in the Mediterranean Basin. Mar. Biol. 1994, 120, 9–15. [Google Scholar] [CrossRef]
  20. Hemminga, M.A.; Duarte, C.M. Seagrass Ecology; Cambridge University Press: Cambridge, UK, 2000; ISBN 0521661846. [Google Scholar]
  21. Ott, J.A. Growth and Production in Posidonia oceanica (L.) Delile. Mar. Ecol. 1980, 1, 47–64. [Google Scholar] [CrossRef]
  22. Vacchi, M.; De Falco, G.; Simeone, S.; Montefalcone, M.; Morri, C.; Ferrari, M.; Bianchi, C.N. Biogeomorphology of the Mediterranean Posidonia oceanica seagrass meadows. Earth Surf. Process. Landforms 2017, 42, 42–54. [Google Scholar] [CrossRef]
  23. Montefalcone, M.; Vacchi, M.; Carbone, C.; Cabella, R.; Schiaffino, C.F.; Elter, F.M.; Morri, C.; Bianchi, C.N.; Ferrari, M. Seagrass on the rocks: Posidonia oceanica settled on shallow-water hard substrata withstands wave stress beyond predictions. Estuar. Coast. Shelf Sci. 2016, 180, 114–122. [Google Scholar] [CrossRef]
  24. Borum, J.; Duarte, C.M.; Greve, T.M.; Krause-Jensen, D. European Seagrasses: An Introduction to Monitoring and Management; M & MS Project; 2004; ISBN 8789143213. Available online: http://www.seagrasses.org (accessed on 4 May 2023).
  25. Boudouresque, C.F.; Pergent, G.; Pergent-Martini, C.; Ruitton, S.; Thibaut, T.; Verlaque, M. The necromass of the Posidonia oceanica seagrass meadow: Fate, role, ecosystem services and vulnerability. Hydrobiologia 2015, 781, 25–42. [Google Scholar] [CrossRef]
  26. Guidetti, P. Differences among fish assemblages associated with nearshore Posidonia oceanica seagrass beds, rocky-algal reefs and unvegetated sand habitats in the adriatic sea. Estuar. Coast. Shelf Sci. 2000, 50, 515–529. [Google Scholar] [CrossRef]
  27. Giakoumi, S.; Kokkoris, G.D. Effects of habitat and substrate complexity on shallow sublittoral fish assemblages in the Cyclades Archipelago, North-eastern Mediterranean sea. Mediterr. Mar. Sci. 2013, 14, 58–68. [Google Scholar] [CrossRef]
  28. Di Franco, A.; Bussotti, S.; Navone, A.; Panzalis, P.; Guidetti, P. Evaluating effects of total and partial restrictions to fishing on Mediterranean rocky-reef fish assemblages. Mar. Ecol. Prog. Ser. 2009, 387, 275–285. [Google Scholar] [CrossRef]
  29. La Mesa, G.; Guidetti, P.; Bussotti, S.; Cattaneo-Vietti, R.; Manganaro, A.; Molinari, A.; Russo, G.F.; Spanò, N.; Vetrano, G.; Tunesi, L. Rocky reef fish assemblages at six Mediterranean marine protected areas: Broad-scale patterns in assemblage structure, species richness and composition. Ital. J. Zool. 2013, 80, 90–103. [Google Scholar] [CrossRef]
  30. Karachle, P.K.; Stergiou, K.I. Diet and feeding habits of Spicara maena and S. smaris (Pisces, Osteichthyes, Centracanthidae) in the North Aegean Sea. Acta Adriat. 2014, 55, 75–84. [Google Scholar]
  31. Zubak Čižmek, I.; Stewart, S.T.; Kruschel, C.; Čižmek, H. Seascape Context as a Driver of the Fish Community Structure of Posidonia oceanica Meadows in the Adriatic Sea. Croat. J. Fish. 2021, 79, 89–109. [Google Scholar] [CrossRef]
  32. Diaz-Gil, C.; Grau, A.; Grau, A.M.; Palmer, M.; Cabrera Castro, R.; Jordà, G.; Catalán, I. Changes in the juvenile fish assemblage of a Mediterranean shallow Posidonia oceanica seagrass nursery area after half a century. Mediterr. Mar. Sci. 2019, 20, 603–615. [Google Scholar] [CrossRef]
  33. Harmelin-vivien, M.; Harmelin, J.; Chauvet, C.; Duval, C.; Galzin, R.; Lejeune, P.; Barnabe, G.; Blanc, F.; Chevalier, R.; Duclerc, J.; et al. Evaluation visuelle des peuplements et populations de poissons: Méthodes et problèmes. Rev. d’Écologie 1985, 40, 467–539. [Google Scholar] [CrossRef]
  34. Harmelin, J.G. Structure and Variability of the Ichthyo-fauna in a Mediterranean Protected Rocky Area (National Park of Port-Cros, France). Mar. Ecol. 1987, 8, 263–284. [Google Scholar] [CrossRef]
  35. Bell, J.D.; Harmelin-Vivien, M.L. Fish fauna of French Mediterranean Posidonia oceanica seagrass meadows. II: Feeding habits. Tethys 1983, 11, 1–14. [Google Scholar]
  36. Francour, P. Fish assemblages of Posidonia oceanica beds at port-cros (France, NW Mediterranean): Assessment of composition and long-term fluctuations by visual census. Mar. Ecol. 1997, 18, 157–173. [Google Scholar] [CrossRef]
  37. Moranta, J.; Palmer, M.; Morey, G.; Ruiz, A.; Morales-nin, B. Multi-scale spatial variability in fish assemblages associated with Posidonia oceanica meadows in the Western Mediterranean Sea. Estuar. Coast. Shelf Sci. 2006, 68, 579–592. [Google Scholar] [CrossRef]
  38. Kalogirou, S.; Corsini-Foka, M.; Sioulas, A.; Wennhage, H.; Pihl, L. Diversity, structure and function of fish assemblages associated with Posidonia oceanica beds in an area of the eastern Mediterranean Sea and the role of non-indigenous species. J. Fish Biol. 2010, 77, 2338–2357. [Google Scholar] [CrossRef] [PubMed]
  39. Zubak Čižmek, I.; Kruschel, C.; Schultz, S.T. Predators structure fish communities in Posidonia oceanica meadows: Meta-analysis of available data across the Mediterranean basin. Mar. Ecol. Prog. Ser. 2017, 566, 145–157. [Google Scholar] [CrossRef]
  40. Budillon, F.; Amodio, S.; Alberico, I.; Contestabile, P.; Vacchi, M.; Innangi, S.; Molisso, F. Present-day infralittoral prograding wedges (IPWs) in Central-Eastern Tyrrhenian Sea: Critical issues and challenges to their use as geomorphological indicators of sea level. Mar. Geol. 2022, 450, 106821. [Google Scholar] [CrossRef]
  41. Budillon, F.; Amodio, S.; Contestabile, P.; Alberico, I.; Innangi, S.; Molisso, F. The present-day nearshore submarine depositional terraces off the Campania coast (South-eastern Tyrrhenian Sea): An analysis of their morpho-bathymetric variability. In Proceedings of the International Workshop on Metrology for the Sea, Naples, Italy, 5–7 October 2020; pp. 132–138. [Google Scholar]
  42. Rendina, F.; Falace, A.; Alongi, G.; Buia, M.C.; Neiva, J.; Appolloni, L.; Marletta, G.; Russo, G.F. The Lush Fucales Underwater Forests off the Cilento Coast: An Overlooked Mediterranean Biodiversity Hotspot. Plants 2023, 12, 1497. [Google Scholar] [CrossRef]
  43. Appolloni, L.; Russo, G.F. Nel Regno di Leucosia; Appolloni, L., Russo, G.F., Eds.; Edizioni dell’Ippogrifo: Castellabate, Italy, 2021. [Google Scholar]
  44. Russo, G.F.; Di Donato, R.; Di Stefano, F. Gli habitat sottomarini delle coste della Campania. Ecodinamica 2008, 6, 37–56. [Google Scholar]
  45. Bacci, T.; Penna, M.; Sante Rende, F.; Tomasello, A.; Calvo, S. Scheda Metodologica Posidonia oceanica (L.) Delile Descrittore 1 Biodiversità (Dlgs 190/10) Elemento di Qualità Biologica Angiosperme (Dlgs 152/06) 2020, 14. Available online: https://www.researchgate.net/profile/Marina-Penna/publication/349589723_Nuova_scheda_metodologica_EQB_Angiosperme_Posidonia_dic_2020_-/links/60376b074585158939ca6c97/Nuova-scheda-metodologica-EQB-Angiosperme-Posidonia-dic-2020.pdf (accessed on 20 May 2023).
  46. Froese, R.; Pauly, D. FishBase World Wide Web Electronic Publication, Version (09/2010). 2012. Available online: www.fishbase.org (accessed on 19 May 2023).
  47. Stergiou, K.; Karpouzi, V.S. Feeding Habits and Trophic Levels of Mediterranean Fish. Rev. Fish Biol. Fish. 2002, 11, 217–254. [Google Scholar] [CrossRef]
  48. Pauly, D.; Palomares, M. Fishing down marine food web: It is far more pervasive than we thought. Bull. Mar. Sci. 2005, 76, 197–211. [Google Scholar]
  49. Evagelopoulos, A.; Batjakas, I.E.; Drosos, K. Structure and Diversity of the Demersal Fish Assemblages off Psara Island (Central Aegean Sea) Caught by Experimental Bottom Trawling. Ann. Int. J. Mar. Sci. 2021, 37, 379–391. [Google Scholar] [CrossRef]
  50. Anderson, M.J. A new method for non parametric multivariate analysis of variance. Austral Ecol. 2001, 26, 32–46. [Google Scholar] [CrossRef]
  51. Anderson, M.J. Permutation tests for univariate or multivariate analysis of variance and regression. Can. J. Fish. Aquat. Sci. 2001, 58, 626–639. [Google Scholar] [CrossRef]
  52. Anderson, M.J.; Ter Braak, C.J.F. Permutation Tests for Multi-Factorial Analysis of Variance. J. Stat. Comput. Simul. 2003, 73, 85–113. [Google Scholar] [CrossRef]
  53. Anderson, M.J.; Willis, T.J. Canonical Analysis of Principal Coordinates: A Useful Method of Constrained Ordination for Ecology. Ecology 2003, 84, 511–525. [Google Scholar] [CrossRef]
  54. Clarke, K.R.; Warwick, R.M. Change in Marine Communities: An Approach to Statistical Analysis and Interpretation; Primer-E Ltd: Plymouth, UK, 1994; pp. 1–172. [Google Scholar]
  55. Anderson, M.J. DISTLM v.5: A FORTRAN computer program to calculate a distance-based multivariate analysis for a linear model. Dep. Stat. Univ. Auckl. N. Z. 2004, 10, 2016. [Google Scholar]
  56. Härdle, W.K.; Simar, L. Canonical correlation analysis. In Applied Multivariate Statistical Analysis; Springer: Berlin/Heidelberg, Germany, 2015; pp. 443–454. [Google Scholar]
  57. Terlizzi, A.; Anderson, M.J.; Fraschetti, S.; Benedetti-Cecchi, L. Scales of spatial variation in Mediterranean subtidal sessile assemblages at different depths. Mar. Ecol. Prog. Ser. 2007, 332, 25–39. [Google Scholar] [CrossRef]
  58. Rojo, I.; Anadon, J.D.; Garcia-Charton, J.A. Exceptionally high but still growing predatory reef fish biomass after 23 years of protection in a Marine Protected Area. PLoS ONE 2021, 16, e0246335. [Google Scholar] [CrossRef]
  59. Brondizio, E.S.; Settele, J.; Díaz, S.; Ngo, H.T. Global Assessment Report on Biodiversity and Ecosystem Services of the Intergovernmental Science-Policy Platform on Biodiversity and Ecosystem Services; IPBES Secretariat: Bonn, Germany, 2019. [Google Scholar]
  60. United Nations, E. The First Global Integrated Marine Assessment: World Ocean Assessment; Cambridge University Press: Cambridge, UK, 2017. [Google Scholar]
  61. Giraud, G. Contribution à la Description et à la Phénologie Quantitative des Herbiers de Posidonia oceanica (L.) Delile; Université d’Aix-Marseille II: Marseille, France, 1977. [Google Scholar]
  62. Pergent, G.; Pergent-Martini, C.; Boudouresque, C.-F. Utilisation de l’herbier à Posidonia oceanica comme indicateur biologique de la qualité du milieu littoral en Méditerranée: État des connaissances. Mésogée 1995, 54, 3–27. [Google Scholar]
  63. Pergent-Martini, C.; Pergent, G. Spatio-temporal dynamics of Posidonia oceanica beds near a sewage outfall (Mediterranean, France). In Proceedings of the International Workshop on Seagrass Biology, Rottnest Island, Australia, 25–29 January 1996; pp. 299–306. [Google Scholar]
  64. Donnarumma, L.; Lombardi, C.; Cocito, S.; Gambi, M.C. Settlement pattern of Posidonia oceanica epibionts along a gradient of ocean acidification: An approach with mimics. Mediterr. Mar. Sci. 2014, 15, 498–509. [Google Scholar] [CrossRef]
  65. Bedini, R.; Bedini, M.; Bonechi, L.; Piazzi, L. Patterns of spatial variability of mobile macro-invertebrate assemblages within a Posidonia oceanica meadow. J. Nat. Hist. 2015, 49, 2559–2581. [Google Scholar] [CrossRef]
  66. Peirano, A.; Niccolai, I.; Mauro, R.; Bianchi, C.N. Seasonal grazing and food preference of herbivores in a Posidonia oceanica meadow. Sci. Mar. 2001, 65, 367–374. [Google Scholar] [CrossRef]
  67. Montefalcone, M.; Chiantore, M.; Lanzone, A.; Morri, C.; Albertelli, G.; Bianchi, C.N. BACI design reveals the decline of the seagrass Posidonia oceanica induced by anchoring. Mar. Pollut. Bull. 2008, 56, 1637–1645. [Google Scholar] [CrossRef] [PubMed]
  68. Guidetti, P.; Baiata, P.; Ballesteros, E.; Di Franco, A.; Hereu, B.; Macpherson, E.; Micheli, F.; Pais, A.; Panzalis, P.; Rosenberg, A.A.; et al. Large-scale assessment of mediterranean marine protected areas effects on fish assemblages. PLoS ONE 2014, 9, e91841. [Google Scholar] [CrossRef]
  69. Valle, C.; Bayle-Sempere, J.T. Effects of a marine protected area on fish assemblage associated with Posidonia oceanica seagrass beds: Temporal and depth variations. J. Appl. Ichthyol. 2009, 25, 537–544. [Google Scholar] [CrossRef]
  70. Letourneur, Y.; Ruitton, S.; Sartoretto, S. Environmental and benthic habitat factors structuring the spatial distribution of a summer infralittoral fish assemblage in the north-western Mediterranean Sea. J. Mar. Biol. Assoc. 2003, 83, 193–204. [Google Scholar] [CrossRef]
  71. Halpern, B.S.; Warner, R.R. Matching marine reserve design to reserve objectives. Proc. R. Soc. B Biol. Sci. 2003, 270, 1871–1878. [Google Scholar] [CrossRef]
  72. Soler, G.A.; Edgar, G.J.; Thomson, R.J.; Kininmonth, S.; Campbell, J.; Dawson, T.P.; Barrett, N.S.; Bernard, A.T.F.; Galván, E.; Willis, T.J.; et al. Reef Fishes at All Trophic Levels Respond Positively to Effective Marine Protected Areas. PLoS ONE 2015, 10, e0140270. [Google Scholar] [CrossRef]
  73. Jennings, S.; Kaiser, M.J. The effects of fishing on marine ecosystems. In Advances in Marine Biology; Academic Press: Cambridge, MA, USA, 1998; pp. 201–212. [Google Scholar] [CrossRef]
  74. Mackay, M.; Hardesty, B.D.; Wilcox, C. The Intersection between Illegal Fishing, Crimes at Sea, and Social Well-Being. Front. Mar. Sci. 2020, 7, 589000. [Google Scholar] [CrossRef]
  75. Prato, G.; Barrier, C.; Francour, P.; Cappanera, V.; Markantonatou, V.; Guidetti, P.; Mangialajo, L.; Cattaneo-Vietti, R.; Gascuel, D. Assessing interacting impacts of artisanal and recreational fisheries in a small Marine Protected Area (Portofino, NW Mediterranean Sea). Ecosphere 2016, 7, e01601. [Google Scholar] [CrossRef]
  76. Ravard, D.; Brind’Amour, A.; Trenkel, V.M. Evaluating the potential impact of fishing on demersal species in the Bay of Biscay using simulations and survey data. Fish. Res. 2014, 157, 86–95. [Google Scholar] [CrossRef]
  77. Stagličić, N.; Matić-Skoko, S.; Pallaoro, A.; Grgičević, R.; Kraljević, M.; Tutman, P.; Dragičević, B.; Dulčić, J. Long-term trends in the structure of eastern Adriatic littoral fish assemblages: Consequences for fisheries management. Estuar. Coast. Shelf Sci. 2011, 94, 263–271. [Google Scholar] [CrossRef]
  78. Jouffre, D.; Inejih, C.A. Assessing the impact of fisheries on demersal fish assemblages of the Mauritanian continental shelf, 1987-1999, using dominance curves. ICES J. Mar. Sci. 2005, 62, 380–383. [Google Scholar] [CrossRef]
  79. Dehghan-Madiseh, S.; Nabavi, S.M.B.; Ghofleh-Marammazi, J.; Jahani, N.; Koochaknejad, E. Application of Abundance Biomass Curve in ecological health assessment of Khure-Mussa (northwest of the Persian Gulf). J. Persian Gulf 2012, 3, 1–10. [Google Scholar]
  80. Appolloni, L.; Ciorciaro, D.; Di Stefano, F.; Donnarumma, L.; Ferrigno, F.; Iacono, C.; Miccio, A.; Rendina, F.; Sandulli, R.; Russo, G.F. Different Metiers Affect Fish Catches Accounting in Marine Protected Areas: A Pilot Investigation Method. J. Environ. Account. Manag. 2022, 10, 237–252. [Google Scholar] [CrossRef]
  81. Soler, G.A.; Edgar, G.J.; Stuart-Smith, R.D.; Smith, A.D.M.; Thomson, R.J. Moving beyond trophic groups: Evaluating fishing-induced changes to temperate reef food webs. Mar. Ecol. Prog. Ser. 2018, 587, 175–186. [Google Scholar] [CrossRef]
  82. Zhou, S.; Smith, A.D.M.; Punt, A.E.; Richardson, A.J.; Gibbs, M.; Fulton, E.A.; Pascoe, S.; Bulman, C.; Bayliss, P.; Sainsbury, K. Ecosystem-based fisheries management requires a change to the selective fishing philosophy. Proc. Natl. Acad. Sci. USA 2010, 107, 9485–9489. [Google Scholar] [CrossRef]
  83. Garcia, S.M.; Kolding, J.; Rice, J.; Rochet, M.J.; Zhou, S.; Arimoto, T.; Beyer, J.E.; Borges, L.; Bundy, A.; Dunn, D.; et al. Reconsidering the consequences of selective fisheries. Science 2012, 335, 1045–1047. [Google Scholar] [CrossRef]
  84. Ghosh, B.; Kar, T.K. Sustainable use of prey species in a prey-predator system: Jointly determined ecological thresholds and economic trade-offs. Ecol. Modell. 2014, 272, 49–58. [Google Scholar] [CrossRef]
  85. Walters, C.; Christensen, V.; Fulton, B.; Smith, A.D.M.; Hilborn, R. Predictions from simple predator-prey theory about impacts of harvesting forage fishes. Ecol. Modell. 2016, 337, 272–280. [Google Scholar] [CrossRef]
  86. Moullec, F.; Barrier, N.; Drira, S.; Guilhaumon, F.; Marsaleix, P.; Somot, S.; Ulses, C.; Velez, L.; Shin, Y.J. An end-to-end model reveals losers and winners in a warming Mediterranean Sea. Front. Mar. Sci. 2019, 6, 345. [Google Scholar] [CrossRef]
  87. Zhou, S.; Smith, A.D.M. Effect of fishing intensity and selectivity on trophic structure and fishery production. Mar. Ecol. Prog. Ser. 2017, 585, 185–198. [Google Scholar] [CrossRef]
  88. Gambi, M.C.; Lorenti, M.; Russo, G.F.; Scipione, M.B.; Zupo, V. Depth and Seasonal Distribution of Some Groups of the Vagile Fauna of the Posidonia oceanica Leaf Stratum: Structural and Trophic Analyses. Mar. Ecol. 1992, 13, 17–39. [Google Scholar] [CrossRef]
  89. Bussotti, S.; Guidetti, P. Fish communities associated with different seagrass systems in the Mediterranean Sea. Nat. Sicil. 1999, 23, 245–259. [Google Scholar]
Figure 1. Marine Protected Areas of (a) Santa Maria di Castellabate and (b) Costa degli Infreschi e della Masseta. Different protection levels (A: fully protected zone; B: general reserve zone; and C: partial reserve zone), Posidonia oceanica beds, and sampling stations at the seagrass upper limit and 15 m depth are shown.
Figure 1. Marine Protected Areas of (a) Santa Maria di Castellabate and (b) Costa degli Infreschi e della Masseta. Different protection levels (A: fully protected zone; B: general reserve zone; and C: partial reserve zone), Posidonia oceanica beds, and sampling stations at the seagrass upper limit and 15 m depth are shown.
Water 15 01967 g001
Figure 2. Posidonia oceanica bed characteristics of the Santa Maria di Castellabate (SMC) and Costa degli Infreschi e della Masseta (CIM) MPAs. (a) mean (±SD) shoot densities (n° shoots/m2) and leaf lengths (cm), divided in green and brown parts (%); (b) mean (±SD) leaf surfaces (m2/m2) and ratio of intact or broken (by mechanical or biological damage) apices (%).
Figure 2. Posidonia oceanica bed characteristics of the Santa Maria di Castellabate (SMC) and Costa degli Infreschi e della Masseta (CIM) MPAs. (a) mean (±SD) shoot densities (n° shoots/m2) and leaf lengths (cm), divided in green and brown parts (%); (b) mean (±SD) leaf surfaces (m2/m2) and ratio of intact or broken (by mechanical or biological damage) apices (%).
Water 15 01967 g002
Figure 3. Descriptive fish community analyses for each protection zone and depth of the Santa Maria di Castellabate (SMC) and Costa degli Infreschi e della Masseta (CIM) MPAs. Synecological indices: (a) species richness (SR) and density (D); (b) specific diversity (H’) and evenness (J); (c) total biomass; and (d) percentage of trophic group frequencies and mean trophic level (TL).
Figure 3. Descriptive fish community analyses for each protection zone and depth of the Santa Maria di Castellabate (SMC) and Costa degli Infreschi e della Masseta (CIM) MPAs. Synecological indices: (a) species richness (SR) and density (D); (b) specific diversity (H’) and evenness (J); (c) total biomass; and (d) percentage of trophic group frequencies and mean trophic level (TL).
Water 15 01967 g003
Figure 4. CAP analysis results on fish densities and biomasses at the Santa Maria di Castellabate (a,b) and Costa degli Infreschi e della Masseta MPAs (c,d). All density, biomass datasets, and fish species are correlated with the Pearson coefficient (ρ > 0.55). Straight lines are vectors of species where the orientation and length are proportional to the most-correlated sites. Circles represent the confidence interval (95%).
Figure 4. CAP analysis results on fish densities and biomasses at the Santa Maria di Castellabate (a,b) and Costa degli Infreschi e della Masseta MPAs (c,d). All density, biomass datasets, and fish species are correlated with the Pearson coefficient (ρ > 0.55). Straight lines are vectors of species where the orientation and length are proportional to the most-correlated sites. Circles represent the confidence interval (95%).
Water 15 01967 g004
Figure 5. Fish ABC curves at the (a) Santa Maria di Castellabate (SMC) and (b) Costa degli Infreschi e della Masseta (CIM) MPAs. Filled squares are abundances, and empty squares are biomasses. Red symbols refer to data detected within the “A” zone; yellow symbols refer to the “B” zone; and green symbols refer to the “C” zone.
Figure 5. Fish ABC curves at the (a) Santa Maria di Castellabate (SMC) and (b) Costa degli Infreschi e della Masseta (CIM) MPAs. Filled squares are abundances, and empty squares are biomasses. Red symbols refer to data detected within the “A” zone; yellow symbols refer to the “B” zone; and green symbols refer to the “C” zone.
Water 15 01967 g005
Figure 6. CCA analysis results on trophic group data at (a) Santa Maria di Castellabate and (b) Costa degli Infreschi e della Masseta MPAs. Pisc: piscivorous fishes; Macro: macroinvertivores; Micro: microinvertivores; Omni: omnivorous species; Herb: herbivorous species; Plank: planktivorous fishes; Detr: detritivorous fishes. Straight lines are vectors of P. oceanica bed characteristics, whose orientation and length are proportional to the most-correlated sites.
Figure 6. CCA analysis results on trophic group data at (a) Santa Maria di Castellabate and (b) Costa degli Infreschi e della Masseta MPAs. Pisc: piscivorous fishes; Macro: macroinvertivores; Micro: microinvertivores; Omni: omnivorous species; Herb: herbivorous species; Plank: planktivorous fishes; Detr: detritivorous fishes. Straight lines are vectors of P. oceanica bed characteristics, whose orientation and length are proportional to the most-correlated sites.
Water 15 01967 g006
Table 1. Mean values (±SD) of P. oceanica bed characteristics collected at the Santa Maria di Castellabate (SMC) and Costa degli Infreschi e della Masseta (CIM) MPAs. D: shoot density (shoots/m2), LL: leaf length (cm), Bp: brown leaf part, LS: leaf surface (m2/m2), Mec: apices showing mechanical damage (%), Bio: apices showing biological damage (%).
Table 1. Mean values (±SD) of P. oceanica bed characteristics collected at the Santa Maria di Castellabate (SMC) and Costa degli Infreschi e della Masseta (CIM) MPAs. D: shoot density (shoots/m2), LL: leaf length (cm), Bp: brown leaf part, LS: leaf surface (m2/m2), Mec: apices showing mechanical damage (%), Bio: apices showing biological damage (%).
DLLBpLSMecBio
SMCSLA550.6 ± 75.027.9 ± 11.6 16.60%7.1 ± 3.716.67%4.17%
B599.5 ± 131.0 42.9 ± 18.7 23.00%12.2 ± 7.120.25%2.53%
C427.3 ± 217.2 37.4 ± 17.011.60%7.5 ± 5.3 13.58%6.17%
15 mA599.2 ± 160.739.5 ± 19.916.40%12.3 ± 9.013.51%0.11%
B569.2 ± 93.244.6 ± 17.8 15.80%13.8 ± 8.318.18%2.27%
C563.2 ± 97.245.0 ± 18.115.90%12.6 ± 7.523.26%2.33%
CIMSLA461.9 ± 114.638.4 ± 19.416.00%9.81 ± 6.7920.23%3.37%
B539.4 ± 85.946.9 ± 19.117.80%13.50 ± 7.7815.91%4.55%
C443.2 ± 86.334.0 ± 22.120.70%7.87 ± 5.9525.93%3.71%
15 mA345.9 ± 109.653.9 ± 23.914.70%9.41 ± 6.4420.48%4.82%
B384.2 ± 94.855.9 ± 22.913.00%11.73 ± 7.0323.86%4.55%
C353.6 ± 69.644.2 ± 22.214.20%7.93 ± 5.7617.86%1.19%
Table 2. Species detected at the Santa Maria di Castellabate (SMC) and Costa degli Infreschi e della Masseta (CIM) MPAs. Families, commercial value (CV; NC: noncommercial; LC: low commercial; C: commercial), and trophic group (TG; Pisc: piscivorous fishes; Macro: macroinvertivores; Micro: microinvertivores; Omni: omnivorous species; Herb: herbivorous species; Plan: planktivorous fishes; Detr: detritivorous fishes) are also reported [48,58].
Table 2. Species detected at the Santa Maria di Castellabate (SMC) and Costa degli Infreschi e della Masseta (CIM) MPAs. Families, commercial value (CV; NC: noncommercial; LC: low commercial; C: commercial), and trophic group (TG; Pisc: piscivorous fishes; Macro: macroinvertivores; Micro: microinvertivores; Omni: omnivorous species; Herb: herbivorous species; Plan: planktivorous fishes; Detr: detritivorous fishes) are also reported [48,58].
FamilySpeciesCVTGSMCCIM
ApogonidaeApogon imberbis (Linnaeus, 1758)NCMicroxx
CarangidaeSeriola dumerili (Risso, 1810)CPisc-x
AlosidaeSardina pilchardus (Walbaum, 1792)LCPlanx-
LabridaeCoris julis (Linnaeus, 1758)NCMicroxx
Labrus merula Linnaeus, 1758CMicroxx
Labrus viridis Linnaeus, 1758LCMicroxx
Symphodus doderleini Jordan, 1890NCMicrox-
Centrolabrus melanocercus (Risso, 1810)NCMicro-x
Symphodus ocellatus (Linnaeus, 1758)NCMicroxx
Symphodus roissali (Risso, 1810)NCMicrox-
Symphodus rostratus (Bloch, 1791)NCMicroxx
Symphodus tinca (Linnaeus, 1758)LCMicroxx
Thalassoma pavo (Linnaeus, 1758)NCMicroxx
MoronidaeDicentrarchus labrax (Linnaeus, 1758)CPiscxx
MugilidaeMugil cephalus Linnaeus, 1758CDetr-x
MullidaeMullus surmuletus Linnaeus, 1758CDetrx-
Mullus barbatus Linnaeus, 1758CDetrxx
PomacentridaeChromis chromis (Linnaeus, 1758)NCPlancxx
SerranidaeSerranus cabrilla (Linnaeus, 1758)CMacroxx
Serranus scriba (Linnaeus, 1758)CMacroxx
Epinephelus marginatus (Lowe, 1834)CPiscx-
SparidaeBoops boops (Linnaeus, 1758)CPlanc-x
Dentex dentex (Linnaeus, 1758)CPiscxx
Diplodus annularis (Linnaeus, 1758)COmnixx
Diplodus puntazzo (Walbaum, 1792)COmnix-
Diplodus vulgaris (Geoffroy Saint-Hilaire, 1817)COmnixx
Diplodus sargus (Linnaeus, 1758)COmnixx
Oblada melanurus (Linnaeus, 1758)CPlancxx
Sarpa salpa (Linnaeus, 1758)LCHerbxx
Sparus aurata Linnaeus, 1758COmni-x
Spicara maena (Linnaeus, 1758)LCPlanxx
Spicara smaris (Linnaeus, 1758)LCPlanxx
Spondyliosoma cantharus (Linnaeus, 1758)COmnixx
SphyraenidaeSphyraena viridensis Cuvier, 1829CPisc-x
Table 3. Synecological indices (mean ± SD) at the Santa Maria di Castellabate (SMC) and Costa degli Infreschi e della Masseta (CIM) MPAs. SR: species richness (species/125 m2); D: density (individuals/125 m2); H’: Shannon diversity index; J: Pielou evenness index; B: total biomass (kg/125 m2).
Table 3. Synecological indices (mean ± SD) at the Santa Maria di Castellabate (SMC) and Costa degli Infreschi e della Masseta (CIM) MPAs. SR: species richness (species/125 m2); D: density (individuals/125 m2); H’: Shannon diversity index; J: Pielou evenness index; B: total biomass (kg/125 m2).
MPADePLSRDH’JB
SMCSLA10.5 ± 1.8106.7 ± 48.1 2.3 ± 0.40.71 ± 0.11.4 ± 0.9
B8.3 ± 2.5127.2 ± 41.7 2.2 ± 0.60.74 ± 0.11.3 ± 1.1
C7.2 ± 2.995.3 ± 56.2 1.7 ± 0.60.59 ± 0.21.3 ± 0.1
15 mA7.8 ± 2.0 92.5 ± 56.0 2.0 ± 0.20.68 ± 0.10.9 ± 0.4
B7.5 ± 2.5116.3 ± 57.4 2.0 ± 0.50.69 ± 0.11.7 ± 1.6
C7.0 ± 0.9122.7 ± 54.82.0 ± 0.50.69 ± 0.22.5 ± 1.4
CIMSLA9.0 ± 2.5141.5 ± 76.21.7 ± 0.30.56 ± 0.20.2 ± 0.2
B8.2 ± 1.6111.3 ± 61.22.3 ± 0.30.78 ± 0.10.7 ± 0.6
C7.2 ± 2.2121.3 ± 49.41.9 ± 0.50.68 ± 0.1 0.2 ± 0.1
15 mA7.2 ± 1.966.5 ± 46.02.1 ± 0.4 0.76 ± 0.10.01 ± 0.0
B7.7 ± 1.090.7 ± 34.9 2.0 ± 0.50.69 ± 0.10.9 ± 0.5
C5.8 ± 0.799.0 ± 40.71.8 ± 0.30.71 ± 0.10.1 ± 0.1
Table 4. PERMANOVA and pairwise results on fish densities and biomasses at the Santa Maria di Castellabate and Costa degli Infreschi e della Masseta MPAs; “A”, “B”, and “C” are the levels of the PL and SL factors, and 15 m are the levels of the De factor; pseudo-F and t values are reported; significant p values (<0.05) are in bold.
Table 4. PERMANOVA and pairwise results on fish densities and biomasses at the Santa Maria di Castellabate and Costa degli Infreschi e della Masseta MPAs; “A”, “B”, and “C” are the levels of the PL and SL factors, and 15 m are the levels of the De factor; pseudo-F and t values are reported; significant p values (<0.05) are in bold.
SMCCIM
Density DataBiomass Data Density DataBiomass Data
FactorPseudo-FpPseudo-FpFactorPseudo-FpPseudo-Fp
PL1.890.0211.880.018PL3.310.0022.880.002
De2.750.0052.270.0182De2.770.0042.860.006
PLxDe0.920.55781.190.26PLxDe2.210.0031.950.003
Within PL factorWithin A level
A-B1.150.231.290.11SL-15 m1.690.0051.820.002
A-C1.710.0031.560.006Within B level
B-C1.240.161.280.11SL-15 m1.060.3411.230.108
Within C level
SL-15 m1.710.0091.280.134
Within SL level
A–B1.620.0151.250.103
A–C1.110.3220.930.532
B–C1.190.2281.110.279
Within 15 m level
A–B1.410.0712.010.002
A–C2.210.0042.010.003
B–C1.910.0041.380.051
Table 5. Characteristics of P. oceanica beds selected by distLM analyses affecting denScheme 2. pseudo-F and p values are reported. Significant p values (<0.05) are in bold.
Table 5. Characteristics of P. oceanica beds selected by distLM analyses affecting denScheme 2. pseudo-F and p values are reported. Significant p values (<0.05) are in bold.
SMCCIM
VariableAdj R2Pseudo-FpVariableAdj R2Pseudo-Fp
DensitiesDensities
Total length2.61 × 10−21.93930.0476Brown1.37 × 10−21.48750.0465
Suface5.06 × 10−21.87750.0638Density2.58 × 10−21.42040.1806
BiomassesGreen part2.70 × 10−21.03970.0408
Total length1.83 × 10−21.65130.094Biomasses
Suface5.11 × 10−22.17690.285Density6.10 × 10−31.2150.025
Density5.20 × 10−21.03030.045Green part1.51 × 10−21.31040.0138
Suface2.46 × 10−21.3220.1894
Brown3.21 × 10−21.24780.2356
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Appolloni, L.; Pagliarani, A.; Cocozza di Montanara, A.; Rendina, F.; Donnarumma, L.; Ciorciaro, D.; Ferrigno, F.; Di Stefano, F.; Sandulli, R.; Russo, G.F. Benthic Fish Communities Associated with Posidonia oceanica Beds May Reveal the Fishing Impact and Effectiveness of Marine Protected Areas: Two Case Studies in the Southern Tyrrhenian Sea. Water 2023, 15, 1967. https://doi.org/10.3390/w15101967

AMA Style

Appolloni L, Pagliarani A, Cocozza di Montanara A, Rendina F, Donnarumma L, Ciorciaro D, Ferrigno F, Di Stefano F, Sandulli R, Russo GF. Benthic Fish Communities Associated with Posidonia oceanica Beds May Reveal the Fishing Impact and Effectiveness of Marine Protected Areas: Two Case Studies in the Southern Tyrrhenian Sea. Water. 2023; 15(10):1967. https://doi.org/10.3390/w15101967

Chicago/Turabian Style

Appolloni, Luca, Alberto Pagliarani, Adele Cocozza di Montanara, Francesco Rendina, Luigia Donnarumma, Domenico Ciorciaro, Federica Ferrigno, Floriana Di Stefano, Roberto Sandulli, and Giovanni Fulvio Russo. 2023. "Benthic Fish Communities Associated with Posidonia oceanica Beds May Reveal the Fishing Impact and Effectiveness of Marine Protected Areas: Two Case Studies in the Southern Tyrrhenian Sea" Water 15, no. 10: 1967. https://doi.org/10.3390/w15101967

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop