# Threats to large brown algal forests in temperate seas: the overlooked role of native herbivorous fish

Fabrizio Gianni, Fabrizio Bartolini, Alexis Pey, Mathieu Laurent, Gustavo M.

Martins, Laura Airoldi, Luisa Mangialajo

► To cite this version:

Fabrizio Gianni, Fabrizio Bartolini, Alexis Pey, Mathieu Laurent, Gustavo M. Martins, et al. Threats to large brown algal forests in temperate seas: the overlooked role of native herbivorous fish. Scientific Reports, 2017, 7, pp.6012. 10.1038/s41598-017-06394-7. hal-01581790

HAL Id: hal-01581790

https://hal.sorbonne-universite.fr/hal-01581790

Submitted on 5 Sep 2017

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
SCIENTIFIC REPORTS

<figure><img src="image_2.png" /></figure>

OPEN

# Threats to large brown algal forests in temperate seas: the overlooked role of native herbivorous fish

Received: 9 March 2017

Accepted: 12 June 2017

Published online: 20 July 2017

Fabrizio Gianni<fnref n="1" />,², Fabrizio Bartolini<fnref n="1" />, Alexis Pey<fnref n="1" />, Mathieu Laurent<fnref n="1" />, Gustavo M. MartinsD<fnref n="1" />, Laura Airoldi<fnref n="4" /> & Luisa Mangialajo<fnref n="1" />,²

Canopy-forming algae are declining globally due to multiple disturbances. This decline has recently been on the increase due to the spread of some tropical herbivorous fishes. This new phenomenon has drawn attention to the effects of fish herbivory in temperate areas, which have been assumed to be negligible compared to that of invertebrates, such as sea urchins. In this study, the impact of a Mediterranean native herbivorous fish (*Sarpa salpa*, salema) was assessed on the canopy-forming seaweed *Cystoseira amentacea* var. *stricta*. *Cystoseira amentacea* forms belts in the infralittoral fringe of wave-exposed shores, which has so far been considered a refuge from fish herbivory. To test the effects of salema feeding on natural *C. amentacea* belts, an innovative herbivore deterrent device was conceived. Salema had a significant effect on *C. amentacea* by decreasing algal size, biomass and fertility, by up to 97%. The results suggest that the contribution of salema feeding to the loss of *Cystoseira* forests in the Mediterranean may have been overlooked. In addition, the analysis of temporal and spatial patterns of salema landings in the Mediterranean Sea suggests that salema abundance may have increased recently. Thus, along with invertebrate herbivory and anthropogenic stressors, fish herbivory may also represent a potential threat to algal forests in temperate areas.

Marine forests of large brown seaweeds are unique habitats that support a great variety of organisms in coastal zones worldwide and are comparable to terrestrial forests for the services they provide<fnref n="1" />,². Several species, mostly belonging to the orders Fucales and Laminariales, can dominate both shallow and deep-water marine ecosystems (up to their light compensation limit)<fnref n="3" />,<fnref n="4" />.

Algal forests are exposed to multiple disturbances that have caused a decline in their abundance in many coastal areas of the world⁵,⁶. A recent survey highlighted that, although global drivers could affect kelp forests at multiple scales, local stressors and regional variations in the effects of these drivers dominate kelp dynamics⁷. Species thriving in shallow-water zones are the most strongly affected because they are located in a boundary environment, exposed to impacts of both terrestrial and marine origin. For instance, contaminants⁸, nutrient enrichment⁹, sediment loads¹⁰, increase in seawater temperature¹¹ and habitat alteration, resulting from coastal urbanisation⁶, are well-known factors responsible for the loss of large brown seaweeds.

Fucoids and kelps can also be highly vulnerable to herbivory¹². Outbreaks of sea urchins, due to natural fluctuations, but more often as an indirect effect of overfishing of their predators, are responsible for the depletion of macroalgal communities and the subsequent formation of extensive barren grounds¹³,¹⁴. This phenomenon has been observed in many coastal areas<fnref n="1" /> and nowadays sea urchin barren grounds are common in many temperate bioregions. As a consequence, sea urchins are considered as the primary herbivorous consumers of algal canopies in the subtidal zone of temperate areas¹⁵. In the intertidal zone, where sea urchins do not usually occur, gastropods and limpets are often considered as the primary herbivores¹⁶,¹⁷. Therefore, in both subtidal and intertidal temperate ecosystems, herbivory is generally associated with invertebrates.

In contrast, the role of herbivorous fishes in regulating macroalgal vegetation is highly variable. In tropical areas, they play an important role in limiting shifts from coral- to macroalgae-dominated communities¹⁸.

<footnote n="1">¹Université Côte d'Azur, CNRS, ECOMERS, Nice, 06108, France. ²Sorbonne Universités, UPMC Univ Paris 06, INSU-CNRS, UMR 7093 Laboratoire d'Océanographie de Villefranche (LOV), Villefranche sur mer, 06230, France.</footnote>

<footnote n="3">³Department of Biology, Faculty of Sciences and Technology & Centre for Ecology, Evolution and Environmental Changes/Azorean Biodiversity Group (cE3c), University of the Azores, 9501-801, Ponta Delgada, Portugal.</footnote>

<footnote n="4">⁴Dipartimento di Scienze Biologiche, Geologiche ed Ambientali BIGEA, University of Bologna, UO CoNISMa, Ravenna, 48123, Italy. Correspondence and requests for materials should be addressed to F.G. (email: fgianni@unice.fr)</footnote>
<figure><img src="image_3.png" /><figcaption>**Figure 1.** Study area. The two sampling sites are located in Villefranche Bay, French Riviera, NW Mediterranean Sea. Maps were made by using Adobe® Illustrator® software CS6.</figcaption></figure>

temperate zones, fish grazing control is assumed to be lower or negligible in comparison to grazing from invertebrates¹⁹⁻²¹. The gradient of control from herbivorous fishes to invertebrates at high latitudes is hypothesized to be driven by the decrease in herbivorous fish species diversity and abundance²².

However, some evidence suggests that the role of herbivorous fishes in temperate regions may have been underestimated. Verlaque observed that the feeding behaviour of *Sarpa salpa* (salema), the only native herbivorous fish in the Mediterranean Sea, is particularly selective, contributing to maintaining high algal biodiversity²³. Salema are also able to determine the distribution of some *Cystoseira* species, restricting grazer-sensitive species to spatial refuges, in either very shallow or deep areas, where fish grazing pressure is lower²⁴. The same observations were reported from central Portugal, where *S. salpa* is responsible for reducing kelp biomass and restricting the survival of recruits to crevices²⁵. Taylor and Schiel did analogous observations on *Odax pullus* in New Zealand, which is able to significantly reduce the cover and biomass of the kelp *Durvillaea antarctica*²⁶. In addition, recent studies have highlighted the impact of tropical herbivorous fishes, by expanding their distribution range into temperate areas, on canopy algae²⁷.

The aim of the present study was to quantify the potential impact of *Sarpa salpa* on large brown algal forests and to provide evidence of a potential increase of salema over the last decades. The herbivory effect of *S. salpa* was investigated on size, biomass and fertility of *Cystoseira amentacea* var. *stricta* (hereafter *C. amentacea*), a species forming narrow belts in the very shallow infralittoral fringe of the Mediterranean Sea, which is considered a refuge from fish herbivory²⁴, ²⁸. In order to experimentally manipulate fish grazing on *C. amentacea*, an innovative herbivorous fish deterrent device was used (DeFish – herbivorous Fish Deterrent for the conservation and restoration of algal forests (Gianni *et al.* in prep.)). The hypotheses of the present study were: i) herbivore exclusion devices are able to effectively reduce *S. salpa* feeding on *C. amentacea*, ii) *C. amentacea* size, biomass and fertility are lower due to fish herbivory if *S. salpa* is not excluded, iii) the magnitude of these effects is greater in low compared to upper level zones of the infralittoral fringe, that are assumed to be less accessible to fish.

To assess potential recent increases of *S. salpa* abundance, temporal patterns in fish landings of salema within the different basins of the Mediterranean Sea were analysed.

## Materials and Methods

**Study area and species.** The experiment was performed between March and June 2014 at two randomly chosen sites in Villefranche Bay, French Riviera (43°41'32"N, 7°18'58"E; Fig. 1). Both sites were characterized by dense intertidal fringing belt populations of *Cystoseira amentacea* (Fig. 2a).

*C. amentacea* is a habitat-forming species, common along North-Western Mediterranean shores²⁹, ³⁰, where it supports complex food webs on rocky-bottoms³¹, ³². Its threatened status is recognised by the Bern and Barcelona Conventions. *C. amentacea* is a long-lived species with a seasonal pattern of growth between March and July. The base is sympodial, formed by a creeping axis from which multiple axes grow. Branches can reach 40 cm length and are lost in late summer, when the dormant season begins. Abundant receptacles, from a few millimetres to 2 cm long, begin to develop at the apical part of all branches at the end of April³³. Since zygotes are heavy and tend to sink close to the parent individual³⁴, *C. amentacea* is characterised by a very limited dispersal ability (dozens of centimetres according to Mangialajo *et al.*³⁵, even if rare events of distant dispersal were suggested by recent studies²⁹, ³⁶).

**Experimental design and set up.** An innovative herbivore deterrent system (DeFish) consisting of a 20 cm long plastic threaded rod with three groups of five cable ties (18 cm long) was designed to deter fish grazing. Each group of cable ties was attached with silicon glue at different heights on the rod and kept straight with a plastic bolt screwed onto the rod (Fig. 2b). To deter fish from passing through the devices, smaller cable ties (10 cm long) were attached perpendicular to the longer ties. Rods were then screwed to drop-in anchors fixed inside holes (2 cm depth) that were drilled into the rock along each side of the protected plots.
<figure><img src="image_4.png" /><figcaption>**Figure 2.** (a) *Cystoseira amentacea* belts in the infralittoral fringe; (b) the deterrent device DeFish used to deter fish grazing on *C. amentacea*; (c) schematic representation of the two vertical zones of the infralittoral fringe characterized by long and short *C. amentacea* branches; a protected plot with DeFish is also shown on the shore; (d) detail of the primary and secondary *C. amentacea* branches with fish bites (white arrows). Photos by Bartolini F. Figure 2c attributions: salema and *Cystoseira* illustrations (Tracey Saxby, Integration and Application Network, University of Maryland Center for Environmental Science (ian.umces.edu/imagelibrary/)).</figcaption></figure>

Compared to the cages that are often used in herbivore exclusion experiments¹², these devices do not require maintenance or cleaning, do not affect light penetration and can be easily removed. However, their installation is time-consuming, logistically challenging for very shallow positions on the shore, and can have a relatively high impact on these fragile environments (involving drilling the rock and fixing by epoxy putty). Therefore a split-plot design³⁷ was applied, at only two sites, to optimise the experimental robustness, while keeping the impact on the forests and logistical constraints at minimum.

Twelve plots (the blocking factor, 40 × 40 cm) were randomly set up in the infralittoral fringe at each of the two sites, spaced several meters apart. Following the split-plot design, two subplots (high-zone and low-zone sub-plots) were attributed to each plot. The high zone subplot corresponds to the upper level of the infralittoral fringe, where belts of *Cystoseira amentacea* are mostly emerged due to wave and tide action and therefore are expected to be less grazed. The low-zone subplot corresponds to the lower level of the infralittoral fringe, where *C. amentacea* belts are almost always immersed and therefore more accessible to fish (Fig. 2c). Then, each plot was randomly assigned to one of three treatments: protected, un-manipulated control and artefact control, for a total of 4 replicated plots per treatment and per site. Protected plots were completely enclosed using at least six deterrent devices fixed on the edges of the plots. In the artefact control plots, three small devices (10 cm high), made with the same materials, were fixed on the edges of the plots, so as to check for a possible effect of the manipulation on *C. amentacea*, and, at the same time, allow fish access. Non-manipulated control plots were marked with epoxy putty.

**Data collection.** The experiment started at the beginning of the *Cystoseira amentacea* growth period, in March 2014. Fish herbivory was estimated by assessing the number of bites on *C. amentacea* branches. *Sarpa salpa* bite marks are evident since the algal branch is clean cut (Fig. 2d). Crabs can also clean cut the branches, but their action is sporadic and limited compared to salema grazing (Gianni et al. under review). Bites cannot be confused with natural broken branches since these are frayed and not precisely broken. The sites selected were virtually inaccessible to people, so that trampling effects could not affect *Cystoseira* branches. A negative effect of the devices on *Cystoseira* branches can be also excluded since they did not swing too much with wave action. To assess the effect of fish feeding on *C. amentacea*, algal size (maximum length in cm), biomass (wet weight in g) and fertility (number of receptacles) were also estimated. Identification of single *C. amentacea* individuals is
<table><caption>**Table 1.** Tables of PERMANOVA calculated on the different parameters: fish herbivory (as number of bites/12.5 cm²), algal size (as mean maximum length in cm), biomass (in g/12.5 cm²) and fertility (as number of receptacles per 12.5 cm²), for both sites and in the different months. Tr: Treatment, Zo: Zone, Pl: plots. Statistically significant values are in bold type. Asterisks indicate that PERMDISP tests were significant (see Supplementary materials).</caption><thead><tr><th rowspan="2"></th><th rowspan="2">df</th><th rowspan="2">Source</th><th colspan="3">March</th><th colspan="3">Cuisse</th><th colspan="3">May</th><th colspan="3">Cuisse</th><th colspan="3">June</th></tr><tr><th>Rubé</th><th>Pseudo-F</th><th>P (perm)</th><th>Rubé</th><th>Pseudo-F</th><th>P (perm)</th><th>Rubé</th><th>Pseudo-F</th><th>P (perm)</th><th>Rubé</th><th>Pseudo-F</th><th>P (perm)</th><th>Rubé</th><th>Pseudo-F</th><th>P (perm)</th></tr></thead><tbody><tr><td rowspan="3">Fish herbivory</td><td>2</td><td>Tr</td><td>1.00</td><td>0.21</td><td>0.82</td><td>6.25</td><td>1.94</td><td>0.17</td><td>204.7</td><td>54.12</td><td><strong>0.005</strong></td><td>97.56</td><td>14.1</td><td><strong>0.006</strong></td><td>184.29</td><td>59.65</td><td><strong>0.004</strong></td><td>66.81</td><td>6.73</td><td><strong>0.013</strong></td></tr><tr><td>1</td><td>Zo</td><td>96.00</td><td>28.51</td><td><strong>0.0003*</strong></td><td>117.04</td><td>27.89</td><td><strong>0.0009*</strong></td><td>31.89</td><td>7.25</td><td><strong>0.02</strong></td><td>109.8</td><td>33.4</td><td><strong>0.0002</strong></td><td>21.40</td><td>9.81</td><td><strong>0.01</strong></td><td>9.37</td><td>1.52</td><td>0.25</td></tr><tr><td>9</td><td>Pl(Tr)</td><td>4.74</td><td>1.40</td><td>0.29*</td><td>3.20</td><td>0.76</td><td>0.64</td><td>3.78</td><td>0.86</td><td>0.59</td><td>6.91</td><td>2.1</td><td>0.14</td><td>3.08</td><td>1.41</td><td>0.30</td><td>9.92</td><td>1.61</td><td>0.23</td></tr><tr><td rowspan="3">Algal size</td><td>2</td><td>TrxZo</td><td>0.68</td><td>0.20</td><td>0.81</td><td>1.62</td><td>0.38</td><td>0.67</td><td>3.25</td><td>0.74</td><td>0.54</td><td>9.03</td><td>2.75</td><td>0.12</td><td>12.08</td><td>5.53</td><td><strong>0.02</strong></td><td>0.84</td><td>0.13</td><td>0.87</td></tr><tr><td>2</td><td>Tr</td><td>1.51</td><td>1.68</td><td>0.23</td><td>0.87</td><td>0.58</td><td>0.57</td><td>154.17</td><td>57.38</td><td><strong>7E-04</strong></td><td>54.66</td><td>4.7</td><td><strong>0.04</strong></td><td>162.21</td><td>65.86</td><td><strong>0.003</strong></td><td>100.2</td><td>21.29</td><td><strong>0.005</strong></td></tr><tr><td>1</td><td>Zo</td><td>32.59</td><td>42.18</td><td><strong>0.0004</strong></td><td>135.3</td><td>67.5</td><td><strong>0.0001*</strong></td><td>57.6</td><td>25.46</td><td><strong>0.001</strong></td><td>276</td><td>28.3</td><td><strong>0.0004</strong></td><td>32.1</td><td>29.37</td><td><strong>0.0007</strong></td><td>207.59</td><td>95.19</td><td><strong>0.0001</strong></td></tr><tr><td rowspan="3">Biomass</td><td>9</td><td>Pl(Tr)</td><td>0.89</td><td>1.16</td><td>0.41</td><td>1.50</td><td>0.75</td><td>0.66</td><td>2.68</td><td>1.18</td><td>0.39</td><td>11.62</td><td>1.19</td><td>0.40</td><td>2.46</td><td>2.25</td><td>0.12</td><td>4.70</td><td>2.15</td><td>0.12</td></tr><tr><td>2</td><td>TrxZo</td><td>0.79</td><td>1.02</td><td>0.38</td><td>0.21</td><td>0.10</td><td>0.90</td><td>4.77</td><td>2.10</td><td>0.17</td><td>26.54</td><td>2.72</td><td>0.11</td><td>0.63</td><td>0.58</td><td>0.57</td><td>0.91</td><td>0.41</td><td>0.65</td></tr><tr><td>2</td><td>Tr</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>1.04E+07</td><td>25.63</td><td><strong>0.005*</strong></td><td>2.11E+07</td><td>13.71</td><td><strong>0.002</strong></td></tr><tr><td rowspan="3">Fertility</td><td>1</td><td>Zo</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>3.77E+06</td><td>32.18</td><td><strong>0.0008</strong></td><td>3.37E+06</td><td>9.85</td><td><strong>0.01</strong></td></tr><tr><td>9</td><td>Pl(Tr)</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>4.05E+05</td><td>3.45</td><td><strong>0.008</strong></td><td>1.54E+06</td><td>0.44</td><td>0.87</td></tr><tr><td>2</td><td>TrxZo</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>1.85E+06</td><td>15.79</td><td><strong>0.001</strong></td><td>4.85E+06</td><td>1.41</td><td>0.29</td></tr><tr><td rowspan="3"></td><td>2</td><td>Tr</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>44779</td><td>8.89</td><td><strong>0.004*</strong></td><td>66686</td><td>40.3</td><td><strong>0.006*</strong></td></tr><tr><td>1</td><td>Zo</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>20817</td><td>8.18</td><td><strong>0.01</strong></td><td>26611</td><td>23.18</td><td><strong>0.001</strong></td></tr><tr><td>9</td><td>Pl(Tr)</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>5031</td><td>1.97</td><td>0.12</td><td>1654</td><td>1.44</td><td>0.26</td></tr><tr><td></td><td>2</td><td>TrxZo</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>18115</td><td>7.11</td><td><strong>0.01</strong></td><td>12901</td><td>11.23</td><td><strong>0.003</strong></td></tr></tbody></table>

difficult, because of its sympodial bearing. Thus, all measurements were completed within a 12.5 cm² reference surface (4 cm diameter circle), assumed to be a proxy of an individual, at the centre of both the high and the low zone of every plot. This approach can mitigate potential biases associated with the split-plot design (i.e. a reduced spatial independence of data). The number of bites was visually assessed in the field together with algal size that was measured with a ruler (accuracy 0.1 cm) in March (before installing the devices), May and June. At the end of June, the experiment was concluded to avoid *C. amentacea* individuals starting to lose their fronds due to high summer temperatures. Therefore, all *C. amentacea* individuals inside the reference surfaces were harvested and transported to the laboratory to estimate biomass and the number of receptacles. The destructive sampling was performed only once, at the end of the experiment, in order to avoid damaging the plots.

Density and size of *Sarpa salpa* were estimated twice, in March and June, during the experiment, using visual counts along ten replicated transects (25 × 5 m) in the experimental area. Transects were performed parallel and perpendicular to the coast on a rocky-sea bottom (maximum depth: 3 m)³⁸. Salema biomass was calculated by using the length-weight relationship reported on FishBase (http://www.fishbase.org).

The datasets analysed during the current study are available from the corresponding author on request.

**Statistical analyses.** Differences in response variables among the treatments and zones were tested by using permutational multivariate analysis of variance (PERMANOVA)³⁹. The use of a multivariate approach on univariate data is justified if a permutational analysis is performed based on Euclidean distance matrices of the measured biological variables, as specified in Anderson *et al.*³⁷. In order to test the homogeneity of dispersions, PERMDISP tests were performed, before permutational analyses, on medians and tables to fit with univariate data³⁷. All PERMDISP tests are reported in the Supplementary Materials. The full split-plot model analysed with PERMANOVA was composed by the fixed factors Treatment with three levels (protect, control, artefact control), Zone with two levels (High and Low Zone) and by the random factor Plot (n = 4), nested in Treatment. The analyses were run using 'Type I Sum of Square', so that the terms were fitted sequentially, according to the split-plot design³⁷. Separate analyses were performed for each time to avoid temporal dependence of data. The sampling sites were also analysed separately because at site 2 a storm in May partially damaged the experimental setup and the deterrent devices did not manage to exclude fish completely. The two sites are therefore referred to hereafter with their local names (i.e. S1: Pointe du Rubé, S2: Pointe de la Cuisse). PERMANOVA p-values were obtained based on 9999 permutations of residuals under the reduced model. Statistical significance was set at the conventional *p* < 0.05 level, except when dispersion tests were significant and, in this case, a more conservative *p*-value (*p* < 0.01) was considered. Pairwise tests were performed on factors or interactions giving statistically significant
<figure><img src="image_5.png" /><figcaption>**Figure 3.** Size class distributions of *Sarpa salpa*. Percentage frequency of salema individuals per size recorded in the two visual census performed during the experiment.</figcaption></figure>

results, calculating Monte Carlo p values in the case of a low number of possible permutations ( < 100). Analyses were completed with the Primer 6 & PERMANOVA + software package.

The net growth potential and reproductive potential of *Cystoseira amentacea* were estimated at both sites. The net growth potential was calculated as the difference in algal branch length between June and March, whereas the reproductive potential was a standardisation in relation to the maximum fertility (number of receptacles multiplied by 100 and divided by the highest number of receptacles recorded at each site). The effect size for protection and zone was calculated on these variables using the log-response ratios⁴⁰ for each zone in each treatment as:

$$R_t = \ln\left(\frac{\bar{x}_{t,z}}{\bar{x}_{c,z}}\right) \qquad (1)$$

where R_t is the log-response ratio for the treatment t, and x̄_{t,z} and x̄_{c,z} are the mean values of net growth potential and reproductive potential, respectively calculated for the treatments protection or artefact control (t) and the treatment control (c), in each zone (z = High or Low).

**Potential patterns of change in salema abundance.** Potential patterns of temporal and spatial variation in salema abundance were investigated in the Mediterranean Sea, using the Mediterranean and the Black Sea capture production quantity (1970–2014) dataset accessed via the FishStatJ (Food and Agriculture Organization, FAO). In accordance with the FAO dataset, the Mediterranean Sea was divided into seven basins (Adriatic Sea, Sardinia Sea, Gulf of Lion, Aegean Sea, Balearic Sea, Ionian Sea, Levant Sea). In order to identify a potential time pattern of the catches (in tonnes), linear trend models were performed for each basin.

## Results

**Sarpa salpa density and size.** The recorded *Sarpa salpa* mean abundance was 0.2 individuals/m² ± 0.06 (SE), corresponding to a mean biomass of 9.88 g/m² ± 2.14. Most individuals (60%) were >14 cm in length (Fig. 3).

**Fish herbivory.** The number of bites recorded on *C. amentacea* ranged between 0/12.5 cm² (in the High Zones of Pointe du Rubé in June) and 27/12.5 cm² (in the Low Zone of an unprotected plot at Pointe de la Cuisse in June). Before setting up the deterrent devices, in March, the number of bites was already significantly greater in the Low Zone compared to the High Zone: on average 0.9 bites ± 0.3 (SE) were recorded in the High Zone and 5.1 ± 1.3 bites in the Low Zone at both sites (Fig. 4 and Table 1). Two months after the setting up of the experiment, the number of fish bites was significantly greater in unprotected than in protected plots (Fig. 4, Table 1 and pairwise tests: Supplementary Tables S1.1–2). The number of bites increased in particular in the High Zone of unprotected plots, even if this number was always lower in comparison to that recorded in the Low Zone (Fig. 4). The experimental protection in the treatment plots effectively limited fish bites at both shore levels and sites (Fig. 4). In June, a significant interaction among the factors Treatment and Zone was recorded at Pointe du Rubé, showing that the two Zones in the protected plots were not statistically different in terms of number of bites due to the effect of protection (Table 1 and pairwise tests: Supplementary Table S1.1). The interaction was not statistically significant at Pointe de la Cuisse, neither the factor Zone (Table 1), potentially due to the fact that the experimental setup was slightly damaged by storms, which may have reduced the ability of the deterrent devices to exclude fish completely.

**Algal size.** The mean maximum length of *Cystoseira amentacea* ranged from 0.5 cm in the Low Zone at both sites in March to 27 cm in the High Zone of a protected plot at Pointe de la Cuisse in May. In March, at the beginning of the growth season, branches of *C. amentacea* were already significantly longer in the High than in the Low
<figure><img src="image_6.png" /><figcaption>**Figure 4.** Fish herbivory. Number of fish bites/12.5 cm² (mean plus SE calculated on all plots, n = 4) for each zone (High and Low) and treatment in the different months and at both sites. P: protected; C: control plots; AC: artefact control. The results of the pairwise tests on the factors Treatment and Zone are reported above the graph (see Table 1 and Supplementary Tables S1.1–2 for more details). Letters above the bars in June at Cap du Rubé indicate significant differences of the pairwise-tests on the interaction TrxZo (Table 1 and Supplementary Table S1.1 for more details).</figcaption></figure>

Zone at both sites (Fig. 5 and Table 1). After two months (and in particular in June), branches of *C. amentacea* were fully developed and significantly longer in the protected plots compared to most unprotected ones and in the High Zone respect to the Low Zone (Fig. 5; Table 1; Supplementary Tables S2.1–2). The length of branches in Control and Artefact Control plots did not differ significantly (Fig. 5 and pairwise tests: Supplementary Tables S2.1–2) and no significant interaction between the main factors was detected (Table 1).

The net growth potential of *C. amentacea*, calculated as the difference in algal length between the beginning and the end of the experiment, responded significantly and positively to herbivore protection in the two zones and sites (Fig. 6). Interestingly, similar effect sizes of the protection against herbivory on *C. amentacea* growth were recorded at the two shore levels at both sites.

**Biomass.** The biomass, expressed as wet weight of *Cystoseira amentacea*, ranged from 0.14 g/12.5 cm² in the Low Zone of an unprotected plot to almost 59 g/12.5 cm² in the High Zone of a protected plot at Pointe du Rubé. A significant Treatment × Zone interaction was highlighted by the analysis of variance, but only at Pointe du Rubé (Table 1). Pairwise tests showed a greater biomass in protected compared to unprotected plots and a consistently
<figure><img src="image_7.png" /><figcaption>**Figure 5.** Algal size. Algal length (mean plus SE calculated on all plots, n = 4) for each zone (High and Low) and treatment in the different months and at both sites. P: protected; C: control plots; AC: artefact control. The results of the pairwise tests on the factors Treatment and Zone are reported above the graph (see Table 1 and Supplementary Tables S2.1–2 for more details).</figcaption></figure>

<figure><img src="image_8.png" /><figcaption>**Figure 6.** Effect sizes of protection and zones based on the net growth potential. Protection (P: protected, AC: artefact control), Zones (H: High, L: Low). Effect sizes are significant if confidence intervals do not overlap zero. R²: log-response ratio for each treatment.</figcaption></figure>
<figure><img src="image_9.png" /><figcaption>**Figure 7.** Biomass. Wet weight (g)/12.5 cm² (mean plus SE calculated on all plots, n = 4) for each Zone (High and Low) and Treatment in June, at the two sites. P: protected; C: control plots; AC: artefact control. Letters above the bars indicate significant differences of the pairwise-test performed on the interaction TrxZo at Pointe du Rubé (see Table 1 and Supplementary Table S3.1 for more details). At Pointe de la Cuisse, the results of the pairwise tests on the factors Treatment and Zone are reported above the graph (see Table 1 and Supplementary Table S3.2 for more details).</figcaption></figure>

<figure><img src="image_10.png" /><figcaption>**Figure 8.** Fertility. Number of receptacles/12.5 cm² (mean plus SE calculated on all plots, n = 4) for each zone and treatment in June. P: protected; C: control plots; AC: artefact control. Letters above the bars indicate significant differences of the pairwise-tests on the interaction TrxZo (see Table 1 and Supplementary Tables S4.1-2).</figcaption></figure>

<figure><img src="image_11.png" /><figcaption>**Figure 9.** Effect sizes of protection and zones based on the reproductive potential. Protection (P: protected, AC: artefact control), Zones (H: High, L: Low). R²: log-response ratio for each treatment.</figcaption></figure>
<figure><img src="image_12.png" /><figcaption>**Figure 10.** Patterns of temporal variation in salema landings. Salema catches (in tonnes) are represented for 7 distinct basins (FAO sub-areas) and the Mediterranean as a whole from 1970 to 2014 (x-axes). Trend lines and the r-squared values are reported on each graph. Note different scales among catch axes.</figcaption></figure>

higher biomass in the High Zone compared to the Low Zone. The High Zone was statistically different from the Low Zone also in the protected plots, but this could be due to the fact that the dispersion test was significant (Fig. 7, Supplementary Table S3.1 and PERMDISP tests). At Pointe de la Cuisse, Treatment and Zone factors were statistically significant (Table 1), showing greater biomass in protected plots and in the High Zone (Fig. 7 and Supplementary Table S3.2). Overall, the fish grazing caused up to 86% of biomass loss, particularly in the Low Zone of the Control plots (Fig. 7).

**Fertility.** The number of receptacles ranged from 0/12.5 cm² in unprotected plots at both sites to 3544/12.5 cm² in the High Zone of a protected plot at Pointe du Rubé. *Cystoseira amentacea* individuals in protected plots and in the High Zone had a greater number of reproductive structures compared to individuals in unprotected plots and in the Low Zone (Fig. 8). In addition, *C. amentacea* individuals in the unprotected plots were often devoid of reproductive structures. The estimated loss of reproductive structures was up to 97%. The analyses of variance showed a significant Treatment × Zone interaction at both sites (Table 1). Pairwise tests showed that the number of receptacles in protected plots were significantly greater than in unprotected ones in both zones and sites (Fig. 8 and Supplementary Tables S4.1–2). However, the differences showed by the pairwise tests on this interaction have
to be taken with caution, because of a significant dispersion of data in the PERMDISP analysis (Supplementary materials).

The reproductive potential, expressed as the ratio of the number of receptacles in each Treatment and Zone compared to the highest number of receptacles recorded at each site, confirmed a clear effect of protection in both zones of the two sites (Fig. 9). Interestingly, at Pointe du Rubé the effect of protection had similar effect sizes in both the High and Low Zones, while at Pointe de la Cuisse effect size of protection was slightly greater in the lower zone, more exposed to grazing.

**Evidence of salema increase.** Salema landings increased since 1970 in four of the seven FAO's Mediterranean sub-areas ( Aegean Sea, Balearic Sea, Ionian Sea and Sardinia Sea) and also in the Gulf of Lion (Fig. 10). In contrast, landings showed a decreasing trend in the Adriatic Sea. In the Levant Sea, salema landings showed a high fluctuation over time (Fig. 10). Despite the spatial variation among FAO sub-areas, overall and across the entire Mediterranean, landings of salema tended to increase with time, showing a 7-fold increase between 1970 and 2014 (Fig. 10). The area investigated in the present study, Villefranche Bay, lies approximately at the intersection between two FAO sub-areas (the Gulf of Lion and the Sardinian Sea, which includes both the Ligurian and the Tyrrhenian Seas). When data from these two sub-areas were combined, there was also an increasing trend in landings of salema over time. The linear trend models showed significant *p*-values for most of the basins, confirming an increase of salema catches with time: Adriatic Sea (n = 44, F = 26.5, *p* = 6.26E-06), Aegean Sea (n = 44, F = 5.23, *p* = 0.02), Balearic Sea (n = 44, F = 24.5, *p* = 1.18E-05), Ionian Sea (n = 44, F = 51.04, *p* = 7.94E-09), Sardinia Sea (n = 44, F = 9.9, *p* = 0.003) and Mediterranean Sea (n = 44, F = 44.3, *p* = 4.07E-08). *P*-values were not significant for the Levant Sea (n = 44, F = 0.0019, *p* = 0.96) and the Gulf of Lion (n = 44, F = 3.31, *p* = 0.07).

**Discussion**

Marine communities are strongly regulated by top-down forces⁴¹. High herbivory rates can lead to regime shifts with a collapse in primary production, biodiversity and ecosystem functioning⁴². The depletion of forests of large brown seaweed and the subsequent creation of extensive barren grounds is typically mediated by sea urchins that are considered the most effective herbivores in temperate regions¹⁵.

The results presented in this study not only confirm that *Sarpa salpa* is able to graze up to the very shallow infralittoral fringe²⁴, but also provide evidence that fish herbivory can be higher than assumed in this zone, considered to be poorly accessible to fish²⁴ and mainly driven by bottom-up processes²⁹, ³⁰. Very high level of herbivory (number of fish bites) on *Cystoseira amentacea* branches was recorded in unprotected plots, resulting in a subsequent reduction in algal features. In fact, during the maximum growth period of *C. amentacea*, from March to June, grazing by salema caused up to 78% reduction in algal size (branch length). Most significantly, fish herbivory affected the biomass of *C. amentacea* up to 86% and fertility up to 97% in the unprotected plots. Interestingly, it was observed that *C. amentacea* branches lacking receptacles seemed to be less intensively grazed than branches bearing reproductive structures, probably because they are rich in fatty acids⁴³. No epiphytes were present on *C. amentacea* branches during the experiment, so they could not have influenced the feeding behaviour of salema⁴⁴.

Differences in grazing along the vertical gradient were observed in the present study, probably because the high level on the shore is accessible to fish only during high tides and with high wave action. Some previous works hypothesised that the physical environment, and in particular the rhythm of emersion/immersion by tides and waves, plays a major role in the development of *Cystoseira* shallow belts in the Mediterranean Sea⁴⁵. Other studies suggest that grazing is a paramount factor affecting the vertical distribution, and in particular the lower limit, of such assemblages²⁴, ³⁵. The present study, based on deterrent devices that were unable to exclude fish completely, did not allow to properly test this hypothesis, as their herbivory-controlling effect may have been dependent on the level of the shore. However, it is worth noting that differences in herbivory and *C. amentacea* features between High and Low Zones were reduced due to the fish deterrent devices, becoming non-significant for the number of bites and the number of receptacles.

The temporal variation recorded in the study was a consequence of the seasonal growth of *C. amentacea* primary branches. Interestingly, at the very beginning of the experiment, when *C. amentacea* biomass was low, herbivory was already very high at the Low Zone, with a particularly clear impact on *C. amentacea* size.

Evidence of the very important role of native herbivorous fishes in temperate areas was also obtained in the South West Pacific. In this region, labrids can remove the entire primary lamina of adult kelp, causing extensive biomass loss and possibly having a significant impact on nutrient cycling²⁶, ⁴⁶. Kelp recruits and juveniles may escape fish herbivory only at wave-exposed sites and under dense canopies that help reduce grazing pressure²⁶. A similar phenomenon has been observed in central Chile, where the blenny *Scartichthys viridis* seems to be able to regulate the presence of a large number of foliose macroalgae in the intertidal zone⁴⁷.

In contrast to the results presented here, Vergés et al. recorded a much lower loss of *C. amentacea* biomass due to *S. salpa* feeding in the Balearic Islands (10% compared to an average of 85% in the present study)²⁴. The recorded difference can be explained by a five times greater density of herbivorous fish at the studied site (0.2 ind./m² ± 0.06, mean ± SE) than at the Balearic Islands site (0.04 ind./m²). In fact, fish density measured in Villefranche Bay is among the highest values reported from the Mediterranean⁴⁸⁻⁵¹. Unfortunately, historical data on salema abundance, useful to assess their potential increase, are not available at this site, to our knowledge.

Natural and human-driven fluctuations of marine organism densities are common, but can often go unnoticed especially for species that are not targeted by fisheries, such as salema. FAO data, even if they need to be interpreted with caution, as they are not standardized regarding fishing effort, provide evidence of a potential increase in the abundance of salema in the Mediterranean since 1970. Increasing salema density over time was also recorded in the Portofino tuna nets (Ligurian Sea, Italy) from the 1950s to the 1970s (data to be taken with
caution due to differences in fishing effort/gear⁵²) and in Portuguese waters in the 1980s⁵³. However, targeted long-time series are scarce, and it is difficult to assess the exact scale and extent of this potential phenomenon. These possible increases may be linked to several factors, most likely in interaction, such as global warming, changes in fishing pressure and decline of predators.

Global warming is driving the rise in seawater temperature⁵⁴, often associated with the spread of invasive species or the proliferation of native species²⁷,⁵⁵. Higher temperatures accelerate metabolic rates⁵⁶ and may favour herbivorous fishes. A striking example is Tosa Bay in southern Japan, where the recent rise in water temperatures enhanced the grazing rate of some tropical fishes, already present in that area, and triggered a shift from kelps to corals⁵⁷. Also in the Mediterranean Sea, the sea surface temperature is rising and this trend is expected to continue to the end of this century⁵⁸; consequently the thermal habitat suitability of salema may increase⁵⁹.

Changes in fishing pressure can favour an increase of *S. salpa* abundance, as it has already been recorded within some Marine Protected Areas (MPAs)⁵⁰,⁶⁰⁻⁶², probably due to fishing regulations. It is worth noting that an increase in salema density, following protection, has been claimed as one of the potential causes of the depletion of large brown seaweed forests in the late nineties in the Portofino MPA⁶³.

Although salema is traditionally fished in several areas, it has never had a high commercial value at fish markets, and data on its abundance are not commonly recorded. Discussions with fishermen in the French Riviera highlighted that the declining demand for Mediterranean fishes, due to the more widespread market distribution of species from more productive environments such as the North West Atlantic, may have caused a corresponding reduction in catches of less valuable species, including salema. At present, at the study area (French Riviera), *S. salpa* is mostly caught accidentally with non-selective gear. The same considerations may apply to other regions, but a coordinated collection of data would be necessary to estimate a possible decrease in herbivorous fishes catches in recent decades.

Another important driver that may have increased salema abundance is the depletion of their predators, such as shark, grouper and leerfish (according to FishBase), due to overfishing⁶⁴,⁶⁵, but relevant data are scarce on this topic, and it is not possible to confirm this cascade effect.

In addition to fluctuations in abundance of native temperate herbivorous fishes due to these factors, several tropical herbivorous fishes have extended their ranges into temperate waters, due to climate change and other anthropogenic factors (e.g. human-facilitated introductions), with a major impact on the structure of benthic assemblages²⁷,⁶⁶. Phase-shifts from large brown macroalgae dominated habitats to less productive barrens were also observed in the Eastern Mediterranean Sea when functionally diverse tropical fishes (*Siganus luridus* and *Siganus rivulatus*) invaded from the Red Sea⁶⁷, replacing the native *S. salpa*, that sharply declined. Even if the North-Western Mediterranean Sea is not currently concerned by this phenomenon, a few individuals of *Siganus* have already been recorded in this basin³⁰,⁶⁸. With the expected increase in sea-water temperature, their potential future spread in the Western Mediterranean Sea over the next decades may further threaten the *Cystoseira* forests.

In conclusion, the results of the present study highlighted the role of *S. salpa* in affecting *Cystoseira* forests and a potential increase in salema populations. Recent studies also suggested a potential role of herbivorous fish on the subtidal *Cystoseira* forests²⁴, ⁴⁸, ⁶⁹, as well as on *Posidonia oceanica* meadows⁷⁰, at depths where sea urchins have been generally considered as the main herbivores¹⁴. The present results support these observations by providing experimental evidence of the role of salema in shaping very shallow *Cystoseira* forests. In the light of these results, we argue that the effects of herbivorous fish on *Cystoseira* forests have been overlooked so far and we cannot exclude that they may have played a major role in the recent loss of algal forests recorded across the Mediterranean Sea, potentially interacting with other well-known causes of loss⁵, ³⁰, ⁷¹⁻⁷³. This is all the more true given that *S. salpa* seems to preferentially graze on reproductive structures, located in the apical part of *Cystoseira* branches. Since most *Cystoseira* species are characterised by limited dispersal ability³⁵, ³⁶, the decrease in fertility may accelerate the loss of forests. The general decrease in complexity of macroalgal communities, from large brown algal forests to turfs of filamentous or coralline algae in other regions⁷⁴, may also be directly or indirectly linked to overlooked fluctuations in abundance of native herbivorous fishes. The present results provide important insights for the conservation of large brown algae forests, not only in the Mediterranean Sea, but also in other geographical areas potentially affected by the same unnoticed phenomenon. The conservation of marine vegetated habitats should take into consideration the role of herbivorous fishes and the assessment of their densities in space and time. The regulation of large herbivorous fish densities should be considered, especially where their abundance is high enough to represent a threat to large brown algae forests conservation or ecological restoration⁷⁵. Devices similar to those proposed in this study, but preferably made of biodegradable materials and using less invasive procedures (without drilling rocks), may represent a solution to protect recently restored or isolated patches of algal forests and promote reproduction and recruitment. A potentially more sustainable measure would be the local regulation of salema populations involving the fishing community and enhancing the species' commercial value.

## References

1. Steneck, R. S. et al. Kelp forest ecosystems: biodiversity, stability, resilience and future. *Environ. Conserv.* **29**, 436–459 (2002).

2. Capdevila, P., Hereu, B., Riera, J. L. & Linares, C. Unravelling the natural dynamics and resilience patterns of underwater Mediterranean forests: insights from the demography of the brown alga *Cystoseira zosteroides*. *J. Ecol.* **104**, 1799–1808 (2016).

3. Nelson, W. et al. Beyond diving depths: deepwater macroalgae in the New Zealand region. *Marine Biodiversity* **45**, 797–818 (2015).

4. Leigh, E. G., Paine, R. T., Quinn, J. F. & Suchanek, T. H. Wave energy and intertidal productivity. *Proc. Natl. Acad. Sci. USA* **84**, 1314–1318 (1987).

5. Mineur, F. et al. European seaweeds under pressure: consequences for communities and ecosystem functioning. *J. Sea Res.* **98**, 91–108 (2015).

6. Coleman, M. A., Kelaher, B. P., Steinberg, P. D. & Millar, A. J. Absence of a large brown macroalga on urbanized rocky reefs around Sydney, Australia, and evidence for historical decline. *J. Phycol.* **44**, 897–901 (2008).

7. Krumhansl, K. A. et al. Global patterns of kelp forest change over the past half-century. *Proc. Natl. Acad. Sci. USA* **113**, 13785–13790 (2016).
8. Sales, M., Cebrian, E., Tomas, F. & Ballesteros, E. Pollution impacts and recovery potential in three species of the genus *Cystoseira* (Fucales, Heterokontophyta). *Estuar. Coast. Shelf Sci.* **92**, 347–357 (2011).

9. Mangialajo, L., Chiantore, M. & Cattaneo-Vietti, R. Loss of fucoid algae along a gradient of urbanisation, and structure of benthic assemblages. *Mar. Ecol. Prog. Ser.* **358**, 63–74 (2008).

10. Eriksson, B. K., Johansson, G. & Snoeijs, P. Long-term changes in the macroalgal vegetation of the inner Gullmar fjord, swedish Skagerrak coast. *J. Phycol.* **38**, 284–296 (2002).

11. Schiel, D. R., Steinbeck, J. R. & Foster, M. S. Ten years of induced ocean warming causes comprehensive changes in marine benthic communities. *Ecology* **85**, 1833–1839 (2004).

12. Poore, A. G. B. et al. Global patterns in the impact of marine herbivores on benthic primary producers. *Ecol. Lett.* **15**, 912–922 (2012).

13. Steneck, R. S. Human influences on coastal ecosystems: does overfishing create trophic cascades? *Trends Ecol. Evol.* **13**, 429–430 (1998).

14. Guidetti, P. & Dulčić, J. Relationships among predatory fish, sea urchins and barrens in Mediterranean rocky reefs across a latitudinal gradient. *Mar. Environ. Res.* **63**, 168–184 (2007).

15. Ling, S. et al. Global regime shift dynamics of catastrophic sea urchin overgrazing. *Philos. T. Roy. Soc. B* **370**, 20130269 (2015).

16. Benedetti-Cecchi, L., Nuti, S. & Cinelli, F. Analysis of spatial and temporal variability in interactions among algae, limpets and mussels in low-shore habitats on the west coast of Italy. *Mar. Ecol. Prog. Ser.* **14**, 87–96 (1996).

17. Hawkins, S. & Hartnoll, R. Factors determining the upper limits of intertidal canopy-forming algae. *Mar. Ecol. Prog. Ser.* **20**, 265–271 (1985).

18. Mumby, P. J. et al. Fishing, trophic cascades, and the process of grazing on coral reefs. *Science* **311**, 98–101 (2006).

19. Gaines, S. D. & Lubchenco, J. A unified approach to marine plant-herbivore interactions. II. Biogeography. *Annu. Rev. Ecol. Syst.* **111**–138 (1982).

20. Choat, J. Fish feeding and the structure of benthic communities in temperate waters. *Annu. Rev. Ecol. Syst.* **13**, 423–449 (1982).

21. Jones, G. P. Ecology of rocky reef fish of north-eastern New Zealand: a review. *N. Z. J. Mar. Freshwat. Res.* **22**, 445–462 (1988).

22. Floeter, S. R., Behrens, M., Ferreira, C., Paddack, M. & Horn, M. Geographical gradients of marine herbivorous fishes: patterns and processes. *Mar. Biol.* **147**, 1435–1447 (2005).

23. Verlaque, M. Relations entre *Sarpa salpa* (Linnaeus, 1758) (Téléostéen, Sparidae), les autres poissons brouteurs et le phytobenthos algal méditerranéen. *Oceanol. Acta* **13**, 373–388 (1990).

24. Vergés, A., Alcoverro, T. & Ballesteros, E. The role of fish herbivory in structuring the vertical distribution of canopy algae (*Cystoseira* spp.) in the Mediterranean Sea. *Mar. Ecol. Prog. Ser.* **375**, 1–11 (2009).

25. Franco, J. et al. Herbivory drives kelp recruits into 'hiding'in a warm ocean climate. *Mar. Ecol. Prog. Ser.* **536**, 1–9 (2015).

26. Taylor, D. I. & Schiel, D. R. Algal populations controlled by fish herbivory across a wave exposure gradient on southern temperate shores. *Ecology* **91**, 201–211 (2010).

27. Vergés, A. et al. The tropicalization of temperate marine ecosystems: climate-mediated changes in herbivory and community phase shifts. *Proc. R. Soc. Lond. B Biol. Sci.* **281**, 20140846 (2014).

28. Thibaut, T. et al. Unexpected temporal stability of *Cystoseira* and *Sargassum* forests in Port-Cros, one of the oldest Mediterranean marine National Parks. *Cryptogamie, Algologie* (2016).

29. Thibaut, T. et al. Unexpected abundance and long-term relative stability of the brown alga *Cystoseira amentacea*, hitherto regarded as a threatened species, in the north-western Mediterranean Sea. *Mar. Pollut. Bull.* **89**, 305–323 (2014).

30. Thibaut, T., Blanfune, A., Boudouresque, C.-F. & Verlaque, M. Decline and local extinction of Fucales in French Riviera: the harbinger of future extinctions? *Mediterr. Mar. Sci.* **16**, 206–224 (2015).

31. Cheminée, A. et al. Nursery value of *Cystoseira* forests for Mediterranean rocky reef fishes. *J. Exp. Mar. Biol. Ecol.* **442**, 70–79 (2013).

32. Thiriet, P. D. et al. Abundance and diversity of crypto- and necto-benthic coastal fish are higher in marine forests than in structurally less complex macroalgal assemblages. *PLoS ONE* **11**, e0164121 (2016).

33. Gómez Garreta, A. et al. Flora Phycologica Iberica Vol. 1. Fucales (2000).

34. Guern, M. Embryologie de quelques espèces du genre *Cystoseira* Agardh 1821 (Fucales). *Vie Milieu* **13**, 649–679 (1962).

35. Mangialajo, L. et al. Zonation patterns and interspecific relationships of fucoids in microtidal environments. *J. Exp. Mar. Biol. Ecol.* **412**, 72–80 (2012).

36. Buonomo, R. et al. Habitat continuity and stepping-stone oceanographic distances explain population genetic connectivity of the brown alga *Cystoseira amentacea*. *Mol. Ecol.* **26**, 766–780 (2017).

37. Anderson, M., Gorley, R. N. & Clarke, R. K. *Permanova+ for Primer: Guide to Software and Statistik Methods* (Plymouth PRIMER-E, UK). 214 (2008).

38. Harmelin-Vivien, M. et al. Evaluation visuelle des peuplements et populations de poissons: méthodes et problèmes. *Revue d'écologie* **40**, 467–539 (1985).

39. Anderson, M. J. Permutation tests for univariate or multivariate analysis of variance and regression. *Can. J. Fish. Aquat. Sci.* **58**, 626–639 (2001).

40. Hedges, L. V., Gurevitch, J. & Curtis, P. S. The meta-analysis of response ratios in experimental ecology. *Ecology* **80**, 1150–1156 (1999).

41. Menge, B. A. Top-down and bottom-up community regulation in marine rocky intertidal habitats. *J. Exp. Mar. Biol. Ecol.* **250**, 257–289 (2000).

42. Pace, M. L., Cole, J. J., Carpenter, S. R. & Kitchell, J. F. Trophic cascades revealed in diverse ecosystems. *Trends Ecol. Evol.* **14**, 483–488 (1999).

43. Vizetto-Duarte, C. et al. Fatty acid profile of different species of algae of the *Cystoseira* genus: a nutraceutical perspective. *Nat. Prod. Res.* 1–7 (2014).

44. Marco-Méndez, C., Ferrero-Vicente, L. M., Prado, P. & Sánchez-Lizaso, J. L. Epiphytes and nutrient contents influence *Sarpa salpa* herbivory on *Caulerpa* spp vs. seagrass species in Mediterranean meadows. *Estuar. Coast. Shelf Sci.* **184**, 54–66 (2017).

45. Epiard-Lahaye, M., Pellegrini, M. & Weiss, H. Influence des rythmes émersion-immersion sur le développement des boutures de *Cystoseira stricta* Sauvageau (Phéophycées, Fucales) en culture. *Bot. Mar.* **30**, 259–266 (1987).

46. Jones, G. & Andrew, N. Herbivory and patch dynamics on rocky reefs in temperate Australasia: the roles of fish and sea urchins. *Aust. J. Ecol.* **15**, 505–520 (1990).

47. Ojeda, F. P. & Munoz, A. A. Feeding selectivity of the herbivorous fish *Scartichthys viridis*: effects on macroalgal community structure in a temperate rocky intertidal coastal zone. *Mar. Ecol. Prog. Ser.* **184**, 219–229 (1999).

48. Hereu, B. Depletion of palatable algae by sea urchins and fishes in a Mediterranean subtidal community. *Mar. Ecol. Prog. Ser.* **313**, 95–103 (2006).

49. Sala, E. et al. The structure of Mediterranean rocky reef ecosystems across environmental and human gradients, and conservation implications. *PLoS ONE* **7**, e32742 (2012).

50. Prado, P., Farina, S., Tomas, F., Romero, J. & Alcoverro, T. Marine protection and meadow size alter fish herbivory in seagrass ecosystems. *Mar. Ecol. Prog. Ser.* **371**, 11–21 (2008).

51. Guidetti, P. et al. Large-scale assessment of Mediterranean marine protected areas effects on fish assemblages. *PLoS ONE* **9**, e91841 (2014).

52. Cattaneo-Vietti, R. a. B. S. *La Tonnarella e la Pesca Tradizionale a Camogli*. 144 (Le Mani-Microart'S 2009).

53. Ribeiro, J. et al. Long-term changes in fish communities of the Ria Formosa coastal lagoon (southern Portugal) based on two studies made 20years apart. *Estuar. Coast. Shelf Sci.* **76**, 57–68 (2008).

54. Doney, S. C. et al. Climate change impacts on marine ecosystems. *Marine Science* **4** (2012).
55. Parravicini, V. et al. Climate change and warm-water species at the north-western boundary of the Mediterranean Sea. *Mar. Ecol.* **36**, 897–909 (2015).

56. Clarke, A. Costs and consequences of evolutionary temperature adaptation. *Trends Ecol. Evol.* **18**, 573–581 (2003).

57. Mezaki, T. & Kubota, S. Changes of hermatypic coral community in coastal sea area of Kochi, high-latitude, Japan. *Aquabiology* **201**, 332–337 (2012).

58. Somot, S., Sevault, F., Déqué, M. & Crépon, M. 21st century climate change scenario for the Mediterranean using a coupled atmosphere-ocean regional climate model. *Global Planet. Change* **63**, 112–126 (2008).

59. Marras, S. et al. Predicting future thermal habitat suitability of competing native and invasive fish species: from metabolic scope to oceanographic modelling. *Conserv. Physiol.* **3**, cou059 (2015).

60. Planes, S., Raventós, N., Ferrari, B. & Alcoverro, T. Fish herbivory leads to shifts in seagrass *Posidonia oceanica* investments in sexual reproduction. *Mar. Ecol. Prog. Ser.* **431**, 205–213 (2011).

61. Raventos, N., Ferrari, B. & Planes, S. Differences in population parameters and behaviour of the herbivorous fish *Sarpa salpa* between protected and unprotected seagrass meadows in the north-western Mediterranean. *J. Mar. Biol. Assoc. U.K.* **89**, 1153–1159 (2009).

62. Guidetti, P. et al. Italian marine reserve effectiveness: does enforcement matter? *Biol. Conserv.* **141**, 699–709 (2008).

63. Parravicini, V. et al. Conserving biodiversity in a human-dominated world: degradation of marine sessile communities within a protected area with conflicting human uses. *PLoS ONE* **8**, e75767 (2013).

64. Ferretti, F., Myers, R. A., Serena, F. & Lotze, H. K. Loss of large predatory sharks from the Mediterranean Sea. *Conserv. Biol.* **22**, 952–964 (2008).

65. Sala, E. The past and present topology and structure of Mediterranean subtidal rocky-shore food webs. *Ecosystems* **7**, 333–340 (2004).

66. Yamano, H., Sugihara, K. & Nomura, K. Rapid poleward range expansion of tropical reef corals in response to rising sea surface temperatures. *Geophys. Res. Lett.* **38** (2011).

67. Sala, E., Kizilkaya, Z., Yildirim, D. & Ballesteros, E. Alien marine fishes deplete algal biomass in the eastern Mediterranean. *PLoS ONE* **6**, e17356 (2011).

68. Daniel, B., Piro, S., Charbonnel, E., Francour, P. & Letourneur, Y. Lessepsian rabbitfish *Siganus luridus* reached the French Mediterranean coasts. *Cybium* **33**, 163–164 (2009).

69. Vergés, A. et al. Tropical rabbitfish and the deforestation of a warming temperate sea. *J. Ecol.* **102**, 1518–1527 (2014).

70. Prado, P., Tomas, F., Alcoverro, T. & Romero, J. Extensive direct measurements of *Posidonia oceanica* consumption confirms the importance of herbivorism in temperate seagrass meadows. *Mar. Ecol. Prog. Ser.* **340**, 63–71 (2007).

71. Blanfuné, A., Boudouresque, C., Verlaque, M. & Thibaut, T. The fate of *Cystoseira crinita*, a forest-forming Fucales (Phaeophyceae, Stramenopiles), in France (North Western Mediterranean Sea). *Estuar. Coast. Shelf Sci.* **181**, 196–208 (2016).

72. Devescovi, M. Effects of bottom topography and anthropogenic pressure on northern Adriatic *Cystoseira* spp.(Phaeophyceae, Fucales). *Aquat. Bot.* **121**, 26–32 (2015).

73. Thibaut, T., Pinedo, S., Torras, X. & Ballesteros, E. Long-term decline of the populations of Fucales (*Cystoseira* spp. and *Sargassum* spp.) in the Albères coast (France, North-western Mediterranean). *Mar. Pollut. Bull.* **50**, 1472–1489 (2005).

74. Wahl, M. et al. The responses of brown macroalgae to environmental change from local to global scales: direct versus ecologically mediated effects. *Perspect. Psychol.* **2**, 11–29 (2015).

75. Gianni, F. et al. Conservation and restoration of marine forests in the Mediterranean Sea and the potential role of Marine Protected Areas. *Adv. Oceanogr. Limnol.* **4**, 83–101 (2013).

## Acknowledgements

This study was carried out under the European Community's Seventh Framework Programme (FP7/2007–2013) [grant number 290056], within the framework of the project MMMPA-Training Network for Monitoring Mediterranean Marine Protected Areas. Funding was also partially provided by Fundação para a Ciência e a Tecnologia (FCT, Portugal) under the project UID/BIA/00329/2013 and in the form of a post-doctoral grant (SFRH/BDP/108114/2015) awarded to GMM. Many thanks to Hugo Dietrich and Dominique Habib for their help in analysing samples, and to Marti J. Anderson and Antonio Di Franco for the precious advice on the experimental design and the statistical analyses. Jeremiah Plass-Johnson undertook a first English revision and offered useful suggestions.

## Author Contributions

F.B. and L.M. conceived the experiment. F.B., P.A., M.L., F.G. and L.M. collected the data and performed the analyses. G.M.M. analysed fish landings data. F.G. and L.M. wrote the manuscript. F.B., P.A., M.L., G.M.M. and L.A. participated in writing and editing the manuscript.

## Additional Information

Supplementary information accompanies this paper at doi:10.1038/s41598-017-06394-7

**Competing Interests:** The authors declare that they have no competing interests.

**Publisher's note:** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_13.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2017