Are we ready for scaling up restoration actions? An insight from Mediterranean macroalgal canopies
Laura Tamburello aff001; Loredana Papa aff003; Giuseppe Guarnieri aff001; Laura Basconi aff004; Serena Zampardi aff001; Maria Beatrice Scipione aff002; Antonio Terlizzi aff001; Valerio Zupo aff002; Simonetta Fraschetti aff001
Authors place of work:
CoNISMa, Roma, Italy
aff001; Stazione Zoologica Anton Dohrn, Napoli, Italy
aff002; Department of Biology, and Environmental Sciences and Technologies, University of Salento, Lecce, Italy
aff003; Department of Environmental Sciences, Informatics and Statistics, Ca’ Foscari University, Venice, Italy
aff004; Department of Life Sciences, University of Trieste, Trieste, Italy
aff005; Department of Biology, University of Naples Federico II, Napoli, Italy
Published in the journal:
PLoS ONE 14(10)
Extensive loss of macroalgal forests advocates for large-scale restoration interventions, to compensate habitat degradation and recover the associated ecological functions and services. Yet, restoration attempts have generally been limited to small spatial extensions, with the principal aim of developing efficient restoration techniques. Here, the success of outplanting Cystoseira amentacea v. stricta germlings cultured in aquaria was experimentally explored at a scale of tens of kms, by means of a multifactorial experimental design. In the intertidal rocky shores of SE Italy, locations with a continuous distribution for hundreds of meters or with few thalli forming patches of few centimeters of C. amentacea canopy were selected. In each location, the effects of adult conspecifics and the exclusion of macrograzers (salema fish and sea urchins) on the survival of germlings were tested. We evaluated the most critical determinants of mortality for germlings, including the overlooked pressure of mesograzers (e.g. amphipods, small mollusks, polychaetes). Despite the high mortality observed during outplanting and early settlement stages, survival of C. amentacea germlings was consistently favored by the exclusion of macrograzers, while the presence of adult conspecifics had no effects. In addition, the cost analysis of the interventions showed the feasibility of the ex-situ method, representing an essential tool for preserving Cystoseira forests. Large scale restoration is possible but requires baseline information with an in-depth knowledge of the species ecology and of the areas to be restored, together with the development of specific cultivation protocols to make consistently efficient restoration interventions.
Algae – Analysis of variance – Death rates – Experimental design – Forests – Grazing – Herbivory – Sea urchins
Coastal ecosystems are globally threatened by multiple and interacting anthropogenic stressors . Shifts from high diversity systems to low diversity ones have been frequently observed, and degraded conditions are often irreversible unless external interventions force the natural recovery of ecosystems . Efforts aimed at mitigating human impacts through management of human activities and conservation of biodiversity, in the attempt of reverting present trajectories of changes, are not always successful , and natural recovery of degraded marine ecosystems rarely occurs in absence of specific conservation strategies . In this scenario, active restoration, alone or in combination with other forms of regulations (e.g. Marine Protected Areas), is considered an effective strategy to assist and speed up the recovery of degraded ecosystems or to foster the backward shift from degraded states [5, 6].
Macroalgal forests (including fucoids and kelps) are severely threatened by several stressors such as direct degradation or destruction of habitat, coastal urbanization, pollution, and herbivores outbreaks, acting in combination with climate change [7–9]. Their decline over vast extensions has been documented in many temperate regions, including the North and South Eastern Pacific [10, 11], the North Eastern and Western Atlantic [12, 13], and the Australian coasts [14, 15]. As many other habitat-forming macroalgal species worldwide, Cystoseira forests have been declining in the whole Mediterranean basin during the last decades [16, 17], and natural recovery has been recorded only occasionally [18, 19]. Although preserving such complex and highly productive habitats represents a priority to maintain the associated biodiversity and ecosystem functioning, restoration of macroalgal forests has been largely neglected compared to other marine habitats [20, 21].
Different methodologies have been attempted for macroalgal reforestation, most frequently testing the efficacy of transplanting adult or juvenile thalli [22–24]. The technique is regarded as successful for several species of Laminariales or Fucales, exhibiting survival rates close to 70% of transplanted adult individuals for Cystoseira species [9, 25, 26]. Yet, the majority of studies revealed spatially variable outcomes, with very low rates or no success of the method, possibly due to dislodgement by intense hydrodynamism [23, 27–29]. Most of these interventions have been carried out at small spatial scales within single locations (never exceeding a few meters) and were limited to a few tens of individuals, with the main purpose to test the efficacy of the method or to apply specific hypotheses within experimental studies (e.g. [26–30]). Yet, the fast growth rate and the high recruitment typical of several kelp species (e.g. Phyllospora comosa ), can favor the successful re-establishment of self-sustaining, expanding adult populations after small-scale interventions, with negligible impact on donor population . On the contrary, the removal of large numbers of adult individuals is hardly sustainable for existing Cystoseira beds and, due to the scarce resilience of compromised canopies, it would result in an irreversible disturbance of donor forests [31, 32].
The outplanting of Cystoseira germlings cultured in laboratory conditions has been suggested as more sustainable than destructive restoration methods [25, 33]. The harvesting of a relatively limited number of fertile apices can provide great amounts of germlings available for large-scale interventions. The natural mortality during the pre- and post-settlement ontogenic shifts (i.e. zygote polarization, adhesion to substrata) can be strongly reduced and the growth to the germling phase can be enhanced through specific cultivation protocols (i.e. optimizing light and nutrient supply), avoiding grazing pressure on eggs and zygotes [30, 33]. Yet, the outplanting process and early settlement stages are likely to be critical for the survival and growth of Cystoseira germlings, representing potential bottlenecks for the realization of large-scale interventions.
Herbivory represents one of the major threats for the survival of Cystoseira juveniles once reintroduced at restoration sites, as it can drastically limit the survival of reintroduced individuals, with estimated losses up to 90% in a few days . In the Mediterranean Sea, grazing by sea urchins and salema fish (Sarpa salpa Linnaeus, 1758) can cause drastic decline of Cystoseira forests, both in shallow subtidal and in intertidal rocky shores [35–37]. In the latter habitat, also limpets and other gastropods may negatively affect the survival of juveniles . Limitation of grazers pressure through caging may represent an efficient strategy to enhance survival of juveniles in the early stages of settlement . However, small-sized recruits can be severely affected also by grazing or bull-dozing disturbance due to mesograzers (e.g. amphipods, small crustaceans and gastropods, polychaetes; sensu Duffy and Hay ), which can be hardly prevented and whose effects are unknown. In addition, the presence of adult Cystoseira conspecifics might ameliorate abiotic conditions and facilitate juvenile survival, by mitigating desiccation stress [40, 41]. Yet, several studies reported lack or negative effects of canopy forests, as juvenile growth might be limited by competitive shading or whiplash [38, 40, 42, 43].
Disentangling the role of the drivers affecting a restoration intervention is a priority to increase the potential of success in reforestation, especially if the challenge is restoring Cystoseira at spatial scales relevant to restore its loss. We focused on Cystoseira amentacea v. stricta Montagne (hereafter C. amentacea) since the genus is considered “of community interest” according to the EU Habitat Directive (92/43/EEC). Cystoseira amentacea has been included among the indicators of environmental quality in Mediterranean coastal waters according to the Water Framework Directive (2000/60/EC), (i.e., EEI and CARLIT). In addition, it is among the species recognized as a priority by the Barcelona Convention and considered vulnerable by several international organizations (e.g. IUCN, RAC/SPA, MedPan), as it is particularly exposed to a suite of human threats. Its reduction/loss has been recorded in several locations in the NW Mediterranean [33, 44], including our study region. However, connectivity among populations is known to be very low, this preventing the recovery of the species once disappeared in one area . Here, the potential of the outplanting of C. amentacea germlings cultured in aquaria was experimentally explored testing the approach for the first time following a hierarchical design, which included the scales of meters, hundreds of meters, and tens of kilometers, allowing to replicate the outplanting technique through space. The manipulation and deployment of tens of thousands of germlings implied facing logistic constraints during the outplanting phase, and contributed to identify the key ecological knowledge and the methodological issues to be addressed for a successful scaling up of Cystoseira restoration.
More specifically we examined 1- the variation in the abundance of germlings during all the outplanting steps; 2- the separate and combined effects of the presence of adult conspecifics and the exclusion of macrograzers (i.e. salema fish and sea urchins) on the abundance of Cystoseira recruits during early settlement phases (i.e. after three months), aiming at identifying which conditions can increase the number of recruits reaching the adult stage; 3- the efficacy of adult transplant as a driver of germlings survival. The final aim was to understand the feasibility of a restoration intervention carried out at a spatial scale of tens of kms also in terms of associated costs, quantifying the mortality associated to each step of the experiment, from the culture in the laboratory to outplanting in the field.
Materials and methods
Field experimental set up
In the intertidal (± 30 cm MSL) rocky shore of Apulia, SE Italy, two alternative conditions were identified: assemblages dominated by a dense canopy of Cystoseira amentacea (> 80% cover), or with sparse and rare individuals of C. amentacea (< 10% cover), representing eligible restoration sites. Two locations for each condition (i.e. donor sites: Marittima and Sant’Isidoro; and restoration sites: Torre Guaceto and Porto Cesareo) were randomly selected among those featured by the requested conditions at a distance of tens of kilometers (Fig 1). On the Adriatic Sea, Marittima (39.997° N, 18.419° E) was characterized by a dense canopy of C. amentacea along the vermetid trottoir with a distribution continuous for hundreds of meters, in an area with limited public access. Both Torre Guaceto and Porto Cesareo are Marine Protected Areas and do not show any evident sign of human threat able to impair this restoration intervention. Within the no-take zone of the MPA of Torre Guaceto (40.716° N, 17.799° E), a rocky platform with sparse individuals of C. amentacea was selected. The historical presence of dense canopy beds is described by Fraschetti et al. , but the species severely declined about ten years ago . On the Ionian Sea, within the MPA of Porto Cesareo, two locations, either characterized by dense canopy beds (Sant’Isidoro; 40.193° N, 17.919° E) or by sparse individuals (Porto Cesareo; 40.247° N, 17.898° E), were identified. Guarnieri et al.  confirmed the historical presence of C. amentacea beds in both locations, although the species has nearly disappeared in Porto Cesareo. Paradoxically, the two locations chosen for restoration, where Cystoseira forests have declined, are within MPAs, both Special Protected Areas for the Mediterranean. The MPA of Torre Guaceto has been recently nominated for a Global Ocean Refuge System Award for its commitment to excellence in marine conservation. The reasons for the decline of the macroalgae are not known, but its conservation status makes the MPA as an excellent area where to implement a restoration intervention. Porto Cesareo is more embedded in a human seascape, but management is effective and water quality is considered very high in the area where the intervention was planned .
At all experimental sites, substrate was nearly horizontal or slightly sloping, and it was either a mixture of vermetid trottoir and rocky platform (Marittima and Sant’Isidoro) or calcarenitic rock (Porto Cesareo and Torre Guaceto). For the duration of the whole experiment, no thermal anomalies were observed in the study region. As far as nutrients, nitrate concentration ranged between 0.700 and 1.545 μmol /l, while fosfate concentration ranged between 0.005 and 0.006 μmol/l in Porto Cesareo and Torre Guaceto, respectively. More information on environmental variables and human pressures for Porto Cesareo, Torre Guaceto and S. Isidoro are available at the website (https://amare.interreg-med.eu/).
Collections and transplantations were carried out under specific permits provided from the MPAs at Torre Guaceto (Prot. 0000995-PM-17, 21/04/2017) and Porto Cesareo (Prot. 0000745-PM-17, 10/04/2017), while no specific permits were required in Marittima, because it is not part of a protected or private area.
Within each location, two sites (tens of meters long) a few hundreds of meters apart were identified. In each site, 18 experimental units (20 x 20 cm quadrats) were permanently marked with epoxy putty and randomly assigned to different treatments.
During the season of maximum vegetative development (June—July) , adult thalli of C. amentacea were transplanted in half of the quadrats at restoration sites, in order to obtain a density comparable to that observed in healthy assemblages (approximately 13 clumps of thalli per quadrat). Specimens were collected from half of the quadrats at donor sites and clumps of adults were dislodged with hammer and chisel, paying attention not to damage their basis. All removed specimens were stored in refreshed coolers with seawater and transported to restoration sites, where they were glued to the substratum with portions of epoxy putty on the bases. Alluminium frames with PVC strings anchored to the substrate facilitated the attachment phase, which was completed within the same day. Specimens’ attachment could require repeated trials before being successful and detached individuals were replaced within the first weeks. Frames were subsequently removed.
As a procedural control of the technique, C. amentacea was cross-transplanted in 3 quadrats between sites of different donor locations (i.e. Marittima and Sant’Isidoro), in order to expose individuals to a similar manipulation and to estimate mortality stress due to handling of C. amentacea adults. Also, in one site for each donor location, specimens were dislodged and relocated in the same position in 3 quadrats, and 3 quadrats were translocated from one site to the other within the same location.
To evaluate the effect of macrograzers on C. amentacea transplants, 30 x 30 x 40 cm exclosure cages were placed in 12 quadrats within each site. Fences made by 0.5 cm double metal mesh were screwed to the substratum and sealed with epoxy putty. Half of the cages had 4 x 4 cm openings on each side, to allow the access of herbivores (e.g. sea urchins, sea snails) and evaluate eventual effects other than herbivory due to the presence of the cage, as control of artifact. For the same reason, 6 quadrats were unfenced in each site, to be accessible to macrograzers (a scheme of the experimental design is presented in Fig 1).
In June 2017, when C. amentacea fronds were fertile at donor sites, apices were collected for fertilization and cultivation of germlings in aquaria. Apices (3–4 cm long) were cut with scissors from about 3000 fronds from the two donor locations, collecting almost 3 fertile apices from each individual, to reduce the risk of sampling highly related specimens due to self‐recruitment, ensuring sufficient genetic variability and avoiding to compromise the reproductive capability of the specimens . Apices packed in aluminum foil were wrapped with seawater-wetted towels and kept in cool, humid and dark conditions during transportation to the nursery facility, to avoid gametes spawning. Transport was completed within 30 hours from collection.
In the nursery, the temperature and photoperiod were selected to reflect typical seasonal conditions at the sampling time during the reproductive season of C. amentacea (June—July). The photoperiod was set to 15:9 h light:dark cycle and the temperature at 20° C; light was provided by LED lamps (AM366 Sicce USA Inc., Knoxville, USA), and irradiance was measured with a LI-COR LI-190/R Photometer (LICOR-Biosciences, Lincoln, NE, USA). Light irradiance was set at 125 μmol photons−2s−1 for the first three days and then at 100 μmol photons−2s−1. Stosch's enriched filtered and autoclaved seawater was used as culture medium following Falace et al. . As the apices were strongly epiphyted, during the first days of culture bacteria and other epiphytes (i.e. diatoms, cianophyta) proliferated. To reduce the growth of these organisms, capable to outcompete germlings, nutrient inputs were reduced to a concentration of 50% of the Stosch's enriched seawater.
A total of 720 clay tiles (ca. 4 cm diameter) were placed in 12 aquaria. Water was continuosly areated by bubbling in the center of each aquarium, and the culture medium was totally renewed every 3 days to minimize possible limiting effect of nutrients depletion.
Once in the laboratory, fertile apices were gently cleaned with a brush and rinsed with sterile seawater, to remove the adhering biofouling and detritus on their surface. Three randomly chosen apices were placed on each clay tile, to guarantee a wide coverage of settled germlings. After 2 hours, gametes were released and were visible on the substrata, so that the apices with receptacles (i.e. swollen ends of the frond containing the reproductive structures) were removed. Cultured germlings grew on clay plates for 24 days, reaching about 0.5 cm of size, before being transported to the storage location and subsequently to the field.
Tiles with germlings were transported from the nursery facility to the experimental locations in July 2017. Tiles were packed to minimize physical damage and maintained in cool conditions during transport. Within 12 hours, they were stored in the laboratory in proximity to restoration sites (storage location, Fig 1), with controlled temperature (22°C). Temperature was slowly increased compared to the nursery value (20°C), in order to acclimatize germlings before outplanting them in the field, where the temperature was 25°C. Outplanting of germling tiles in the field was completed within 3 days. In the storage location, culture medium was replaced with filtered seawater, which was renewed daily. Tiles were maintained in cool conditions during transport from the storage location to the field. In each quadrat, 5 tiles with germlings were fixed with epoxy glue in the inner 20 x 20 cm.
Before transporting germlings from nursery facilities to the storage location (i.e., day 24), 194 randomly chosen tiles from different aquaria were photographycally sampled to quantify the number of recruits. All 720 tiles were sampled in the same way before being fixed to the shore, to estimate the number of germlings that reached the experimental locations during the three days of outplanting (i.e., 180 tiles on day 26, 360 tiles on day 27, 180 tiles on day 28).
The abundance of juveniles was assessed by mean of photographic sampling carried out in October 2017. In the same occasion, the cover of C. amentacea adults in each experimental plot was also visually estimated with a 20 x 20 cm frame subdivided into 25 sub-quadrats. A score between 0 (absence) and 4 (an entire sub-quadrat covered) was assigned to each sub-quadrat, adding up all the 25 estimates to a final value expressed as percentage .
In addition, a destructive sampling was carried out in October 2017 to characterize the benthic assemblage of the study area. Algal assemblages and associated sessile and vagile fauna were collected integrating a pressure vacuum-cleaner connected to a tank and scraping with hammer and chisel. Five 20 x 20 cm randomly placed quadrats were collected at each site. Macroalgal taxa were identified and distinguished into morphological groups, and their abundance was estimated as wet weight (g / 400 cm2). In three randomly-chosen samples per site, invertebrates larger than 1 mm were identified at the species level and their abundance was expressed as number of individuals / 400 cm2. A literature research allowed to identify the functional traits of the identified species with a focus on their trophic group (S1 Table). In May, July and October 2017 the number of sea urchins (Paracentrotus lividus (Lamarck, 1816); and Arbacia lixula (Linnaeus, 1758)) was estimated in 20 randomly placed 50 x 50 cm quadrats at each site.
Linear mixed-effects models (LME) were used to evaluate the abundance of germlings surviving during the three days of outplanting in the field, as the method can handle unbalanced data with a hierarchical temporal structure . The number of germlings per tile was transformed in a logarithmic scale to normalize the data . The model included the fixed factors Day and Halfday (nested within Day, indicating morning and afternoon) and Halfday|Day as random effects. The structure of the random term was selected comparing models with different error structures using the Akaike information criterion (AIC, ) (S4 Table). Significance of the fixed factors was assessed by means of Wald test [52, 53] (S5 Table). To partition and compare variability among factors, variance components of the number of germlings per tile were estimated among days of outplanting, between halfdays within each day, and among tiles within each halfday with a model including 1 | Day/Halfday as random effects . Model assumptions were assessed by visually inspecting plots of residuals vs. fitted values (S2 Fig). The models were fitted using the function lmer in the R package lme4 .
Due to unexpected mortality during the outplanting phase, effects of post-settlement processes (i.e. macrograzing and presence of C. amentacea adults) were evaluated in two locations (Marittima and Torre Guaceto, respectively one donor and one restoration site). To avoid confounding the influence of manipulated factors (i.e., presence of adult conspecifics and herbivores access) with stochastic effects naturally affecting small populations, we a-priori decided to retain only quadrats containing a number of germlings at least equal to 100 at the time of outplanting, ideally guarantying the survival of at least 10 C. amentacea individuals per quadrat (estimates were based on previous experience gained during a small-scale pilot study). At Marittima and Torre Guaceto, nearly all quadrats satisfied this condition (respectively 34 and 32 units out of 36), while at Sant’Isidoro and Porto Cesareo it was rarely fulfilled (respectively 21 and 3 units out of 36, not homogeneously distributed across experimental conditions). Hence, only Marittima and Torre Guaceto were retained for the analysis of post-settlement processes (i.e. macrograzing and presence of C. amentacea adults).
Effects of post-settlement processes on the survival of juveniles after three months were assessed with the analysis of variance (ANOVA) on the total number of germlings per quadrat (i.e. summing the values of the 5 tiles), as sampling units were balanced and all assumptions of the analysis were met. The experimental design included four factors (Fig 1): Location (fixed, with two levels, “Marittima” and “Torre Guaceto”), Site (random, with two levels, nested within Location), Adult transplant (fixed, orthogonal to all other factors, with two levels, “Transplanted adults” and “No adults transplanted”) and Herbivory (fixed, orthogonal to all other factors, with three levels, "Herbivores", "No herbivores", "Artifact control"), with n = 3 replicates for each combination of factors. To account for the variability among quadrats in the number of germlings deployed during outplanting, the analysis was repeated including the initial number of germlings as a covariate. However, as the covariate was not significant, and it did not increase the predictive power of the model, results are not shown.
The efficacy of adult transplant was assessed analyzing C. amentacea cover with ANOVA, hypothesizing that the considered variable would not differ among Conditions (i.e. Donor populations, with unmanipulated C. amentacea, and Restored populations, where the macroalgae had been transplanted). The experimental design included four factors: Condition (fixed, with two levels, "Donor", representing the unmanipulated control, and "Restoration", including transplanted quadrats), Location (random, with two levels, nested within Condition), Site (random, with two levels, nested within Location), and Herbivory (fixed, orthogonal to all other factors, with three levels "Herbivores", "No herbivores", "Control of artifact"), with n = 3 replicates for each combination of factors. Also, the absence of artifact associated to the transplant technique was verified from the lack of differences between control of artifact treatments (i.e. cross transplant between donor locations, translocation between sites within location and dislocation within site) and unmanipulated controls.
ANOVA was also employed to test patterns of abundance of mesograzers. Species were classified according to their trophic group and the total number of grazers and omnivores belonging to each of the most abundant taxonomical groups (mollusks, amphipods, polychaetes) was analyzed by means of two-way ANOVA. As there was no a-priori hypothesis on differences in patterns of distribution of mesograzers among Conditions, the model was simplified to include Location (random, with 4 levels) and Site (random, with 2 levels nested in Location).
In all ANOVAs, homogeneity of variances was checked with Cochran’s C-test and, when necessary, data were log or square-root transformed (respectively abundance of Cystoseira germlings and amphipods). A posteriori Student-Newman-Keuls (SNK) tests were used to compare 1- differences in the abundance of settled Cystoseira germlings among levels of Herbivory within each Site, 2- differences in the abundance of mesograzers among Locations, 3- differences in the cover of Cystoseira adults among Sites. Analyses of variance were performed in R v3.5 (R Core Team 2013) using the package GAD .
To characterize patterns of distribution of macroalgae among study sites, the assemblage structure was analyzed by means of a permutational multivariate analysis of variance (PERMANOVA) [56, 57]. The model included two factors: Location (random, with two levels) and Site (random, with 2 levels nested in Location). Although locations were originally chosen randomly and so treated as a random factor in the analysis, it was nevertheless of interest to assess differences among Locations by using pairwise a posteriori comparison, to better interpret patterns of abundance of Cystoseira germlings after three months. Permutational analysis of multivariate dispersion (PERMDISP) was performed to evaluate the homogeneity of multivariate dispersion among locations. Analyses were done on Bray-Curtis dissimilarity [57, 58] matrix of non-transformed wet weight of macroalgal species in five 20x20 cm quadrats for each site. Similarity percentage analysis (SIMPER) identified the percentage contribution of each macroalgal species to the Bray-Curtis dissimilarity metric between pairs of locations. A 2-dimensional MDS (multidimensional scaling) based on Bray-Curtis measures of dissimilarity was used for a graphical representation of the multivariate structure of macroalgal assemblage. Multivariate analyses have been carried out using PRIMER version 6, including the add-on package PERMANOVA+ .
The cost of restoration intervention was calculated separately for the nursery phase in the laboratory and for the set-up and outplanting in the field. Similar to the approach by Carney et al.  and Verdura et al. , transportation of fertile apices and recruits, travel, equipment and personnel expenses were evaluated. Also, an evaluation of the monitoring cost for one date of sampling was provided. Travel expenses were estimated equal to 0.21 €/km or to 0.27 €/km when highway toll was included, while the personnel cost was estimated equal to 7 € hour/person.
Effect of experimental conditions on germlings
At the end of the nursery phase, tiles were colonized by an average number of 450.65 germlings (SE ± 47.04, n = 194). During transports, the average number of germlings per tile reduced to 34.25 (SE ± 1.45, n = 720). The number of germlings surviving during the outplanting process decreased proportionally to the day of placement in the field and in the afternoon compared to the morning within each day (Fig 2A, S5 Table). Although nearly half of the variability in the number of germlings among tiles (50.71%) was due to factors other than the time of outplanting, the partitioning of variance in Fig 2B indicates that it was partially explained by the time elapsed before outplanting (variability among days 33.41%; within the same day 15.88%).
Effect of Cystoseira adults and herbivory on germlings
Three months after settlement in the field, the mean number of germlings per tile had decreased to 0.69 (SE ± 0.12, n = 307), and 53 out of 360 tiles had been lost. Grazers exclusion differently affected the abundance of C. amentacea germlings among experimental sites (Herbivory X Site (Location) MS = 1.678, F4,48 = 2.718, P < 0.05; Fig 2C, S6 Table). At Marittima, where the overall abundance of germlings was greater (Location MS = 6.357, F1,2 = 39.153, P < 0.05), germlings were favored by grazer exclusion, while at Torre Guaceto a similar trend, although not significant at the SNK test, was clearly observed only at site 2 (Fig 2C, S4 Table). Abundance of germlings was not influenced by the presence of adult conspecifics (Adult Transplant MS = 5.052, F1,2 = 5.657, P = 0.140; S6 Table), although transplant of C. amentacea adults was efficient at maintaining canopy cover comparable between donor and restoration sites (S1 Appendix, S3 Fig and S7 Table).
Multivariate analyses evidenced differences among locations in the structure of macroalgal assemblage (S2 Table, S1 Fig). Beyond C. amentacea, the majority of the macroalgae mostly contributing to differentiate assemblages among locations were species mechanically resistant to grazing (i.e. Ellisolandia elongata, Jania rubens, Corallina officinalis, Titanoderma pustulatum, Halimeda tuna, S3 Table). Also, the analyses on the patterns of distribution of macro and mesograzers evidenced differences among locations. Sea urchins, either P. lividus or A. lixula, were constantly present only in the locations where C. amentacea had naturally disappeared (i.e. restoration sites, Fig 3A). Significant differences among locations were detected also for mesograzers in such taxa as mollusks (Location MS = 668.7, F3,4 = 5.94, P < 0.05, Fig 3B) and polychaetes (Location MS = 598.8, F3,4 = 7.60, P < 0.05, Fig 3D), but not for amphipods (Location MS = 84.9, F3,4 = 5.52, P = 0.06; Fig 3C). SNK tests for mollusks and polychaetes evidenced a lower number of grazers at Marittima compared to other experimental locations, and a similar trend, although without significant differences, was evident also for amphipods. The most abundant herbivores among mollusks were the chitons Callochiton septemvalvis (Montagu, 1803), Lepidochitona monterosatoi Kaas & Van Belle, 1981, and Acanthochitona fascicularis (Linné, 1767), while polychaetes were dominated by Platynereis dumerilii Audouin & Milne-Edwards, 1834 and amphipods by Protohyale schmidtii (Heller), 1866 and Ampithoe ramondi Audouin, 1826.
The total cost of the tested restoration approach and follow-up activity was estimated equal to 19755 € (Table 1). The total number of working hours corresponded to 240 in the laboratory and 1598 in the field (including 48 hours for the collection and processing of fertile apices, 170 hours for cages construction, 876 hours for cages fixing, 292 hours for the transplant of Cystoseira adults and 212 hours for the installation of germlings tiles). Beside the personnel cost, which corresponded to 66%, the equipment consistently represented the highest source of costs, corresponding to the 22% of the total for the intervention. The activities related to herbivory limitation (cages construction and installation) and adult transplant respectively contributed for 45% and 16% of the total cost.
The present study is the first attempt aimed to restore C. amentacea beds over a scale of tens of kms, through extensive outplanting of germlings cultured in aquaria. While this active restoration approach offers several advantages compared to transplantation of adults (e.g., availability of large amounts of individuals with negligible impact for extant macroalgal beds), there is the urgency to identify the procedural steps potentially limiting a tangible scale up of intervention. We separately evaluated how germlings survival was affected by experimental conditions, presence of adult conspecifics and herbivory, including the consideration of the associated costs. The study was featured by a limited number of sites, driven by the challenge of a very intense field work, and the total duration was limited to three months, which is a short time to fully evaluate the success of the restoration actions. However, our intention was to quantitatively describe those steps in the early stages of restoration interventions making restoration particularly vulnerable, as a mandatory knowledge to set future large-scale interventions.
Effect of experimental conditions on germlings mortality
Mortality of germlings during the restoration intervention was very high in each phase: only 7% of individuals grown in the nursery reached restoration locations and only 2% of them survived the initial settlement stage. In the literature, high mortality rates are normally encountered by macroalgal recruits [61, 62], and low survival rates in the first year have been reported also for Cystoseira sp. [63, 64]. Schiel et al.  estimated survival rates from embryos to visible recruits variable from 0.001% to 3.8%. Besides, the first 2–3 months since recruits become visible represent a critical phase, accounting for nearly 40–50% of the mortality occurring within 3 years [40, 41].
Germlings mortality can be density-dependent [40, 41], and intraspecific competition might be exacerbated by limited nutrient availability and elevated temperatures . High temperature itself may reduce survival of juvenile macroalgae [66–68], although undesirable effects have not been detected for the growth of Cystoseira barbata juveniles . Yet, the temperature range experimentally investigated by Irving et al.  simulates naturally occurring spring values (16°), which might result less stressing compared to the extreme temperature values tested by other authors [65–68], better resembling the stressful summer conditions occurring in our restoration experiment. Although germlings were constantly refrigerated during transportation and outplanting, temperature stress might justify the high mortality and in particular the increase in mortality observed for individuals placed in the field in the afternoon compared to those outplanted in the morning, since solar heat can represent a stress during attachment.
Our study evidences increased mortality for germlings stored in the laboratory for prolonged time (24–48 hours before outplanting) and a large percentage of variability in the number of germlings per tile. Experimental conditions maintained in the laboratory are crucial to reduce mortality, representing also a cost as shown by Verdura et al. . Also, abrupt variations of environmental conditions and mechanical damage during transportation or manipulation might have contributed to the mortality of germlings before and during the settlement in the field. Our results show that the availability of nursery facilities in close proximity to restoration sites can significantly contribute to the success of intervention. This is a critical point considering present and future need of restoration actions in areas of the world where facilities are still lacking. In particular, the identification of strategies to minimize mortality during transport and outplanting represents a priority for future interventions. Specific experiments evaluating how vulnerability of recruits to abiotic stress varies with age could indicate the optimal duration of laboratory cultivation phase and the eventual advantage due to postponing the outplanting in the field.
Effect of Cystoseira adults on germlings mortality
Survival of Cystoseira germlings settled in the field was not influenced by the presence of adult conspecifics. Positive and neutral intraspecific interactions have been found common in seaweed communities globally . In particular, recruitment occurring preferentially below canopy fronds or at the edge of canopy beds has been documented for several forest-forming macroalgae, both in the intertidal and in subtidal habitats [23, 70]. The presence of adult conspecifics creates favorable conditions for the settlement of juveniles, by protecting recruits form grazers or competitors [71, 72] or by ameliorating abiotic environmental conditions [73, 74]. In particular, in the intertidal fringe, desiccation represents a major source of stress for macroalgae, as it can impair photosynthetic activity, reduce growth or increase brittleness and mortality due to dislodgement by waves . However, desiccation stress could have exerted a minor role at our study sites, where C. amentacea colonized vermetid-platforms or rocky coastal belts at depths virtually constantly submerged by water. In contrast, the presence of canopy fronds is known to inhibit recruitment of several kelp and fucoids species ( and references therein; ), and neutral or negative effects have been recorded for different Cystoseira species inhabiting subtidal reefs [18, 28, 43, 63]. Algal canopy may negatively affect the survival of germlings through frond abrasion  or by reducing light irradiance up to 95–99%, a resource of primary importance for the development and growth of Cystoseira recruits . Variation in the outcomes of intraspecific interactions is clearly common in nature and can change as a function of the biotic and/or abiotic context.
Effect of herbivory on germlings mortality
Preventing the access of macrograzers by caging successfully guaranteed the survival of germlings in our study, while the absence of protection led to nearly null survival (0.42 ± 0.17 individuals per quadrat). Overgrazing is known to be an important cause of decline of macroalgal forests worldwide and represents the primary cause for unsuccess in macroalgal restoration intervention, with drastic effects described for kelp-forming macroalgae [59, 72] and Cystoseira species [20, 35]. Signs of grazing by salema fish have been observed on transplanted adult individuals of C. amentacea and C. compressa by Susini et al. . Mangialajo et al.  suggested that salema fish may consume also C. amentacea recruits. Salema fish is a very voracious herbivore, capable to graze also on intertidal reefs, to reduce biomass up to 90%, and to decrease fertility by preferentially consuming apices with receptacles [34, 37]. At our study sites, salema schools were frequently observed grazing on the vermetid platform, and signs of bites were noticed on uncaged transplanted adults. Control of artifact cages had windows to allow entrance by grazers, although the complexity of the structure might have prevented grazing by salema fish. Hence, we can suppose that the amount of germlings survived into control of artifacts might represent the number of individuals that survived salema grazing and was impacted only by other herbivores.
Similar to our study, experiments aiming at introducing Cystoseira spp. on artificial reefs successfully adopted caging to prevent grazing of juveniles or displacement and damaging due to non-consumptive biotic interactions (i.e. handling and clipping of thalli; [28, 77]). These experiments evidenced that, in addition to the impact of salema fish, juvenile survival might be threatened by several herbivores and omnivores species, such as mullets, crustaceans (e.g. crabs) or mollusks (e.g. limpets) [34, 77]. The sea urchins P. lividus and A. lixula, although scarcely abundant in intertidal rocky shores, may severely reduce the survival and prevent the recovery of Cystoseira in proximity of sea urchins’ refugees . At our study sites, sea urchins likely exerted a prominent role in determining patterns of distribution of C. amentacea, as they were present only in locations were C. amentacea had historically disappeared. Densities of A. lixula in Porto Cesareo were nearly half compared to those of P. lividus in Torre Guaceto. However, A. lixula is very efficient at maintaining barren habitats, and even low densities can prevent Cystoseira spp. recovery [31, 79].
Other species might be efficient grazers of the small-size, fragile Cystoseira germlings. A high pressure of herbivores at restoration sites (i.e. Porto Cesareo and Torre Guaceto) and in Sant’Isidoro is suggested by the abundance of articulated corallinales (e.g. Jania rubens and Ellisolandia elongata) or of other algal species (e.g. Halimeda tuna) which are mechanically resistant to grazing and scarcely appealing for herbivores . Analyses of the abundance of mesograzers confirmed a reduced presence of herbivore mollusks and polychaetes at Marittima, where we observed the greater survival of Cystoseira germlings. Also, the highest abundance of mollusks at one site in Torre Guaceto corresponded to the lowest survival of Cystoseira germlings, suggesting grazing as a plausible cause of mortality. Effects of mesograzers on macroalgal forests have rarely been studied , possibly due to the difficulties in manipulating them in the field for experimental purposes. However, some studies have shown that small-size species, such as amphipods or gastropods, may be crucial in regulating the success of recruitment of Laminariales or Fucales [82–84]. In some species of macroalgae, juveniles, having lower content of polyphenolic compounds, may be more palatable for mesograzers than adult conspecifics . Ad hoc experimental studies in the laboratory would help elucidating the overlooked impact of mesograzers on early post-settlement stages of Cystoseira, although the small size and high mobility of the herbivorous do not allow to implement strategies for reducing their abundance in the field for restoration purposes. While a characterization of the fish and the benthic assemblages at putative restoration sites may contribute to discard unappropriated locations, the impact of mesograzers should be reckoned as a unavoidable, further cause of mortality for juveniles.
As demonstrated, 135 germlings survived in the optimal conditions (locations with timely outplanting after transport from the nursery and macrograzer exclosure, independently from the presence of adult conspecifics). If resources would have been allocated for all 144 experimental units to this condition, we could have grown nearly 810 germlings, with an estimated cost per individual ~ 23.30 €. This value is nearly double the cost estimated for the culturing and outplanting of Nereocystis luetkeana by Carney et al. . Yet, in their study grazing contributed to the failure of juvenile transplant intervention, and the authors advocated for the use of grazer exclusion devices in sites were grazers are abundant, with consequent increase of the cost of intervention per individual. In our intervention, limitation of herbivory severely conditioned the duration and cost of intervention, requiring 55% of working hours in the field and 45% of total expenses. Recently, Verdura et al.  compared the cost of restoring a population of Cystoseira barbata using the ex situ and in situ methods, considering the travel, transportation, personnel and material expenses. Similar to Carney et al. , they did not use any structure to limit herbivory, as the site did not require them. They concluded that the cost of restoring 25 m2 of this species corresponded to 2665 € using the ex situ seeding technique, with an estimated cost per individual ~ 25.60 €, for a total of ~ 104 plants per site. Although the cost is comparable to our intervention, it is worth noting that their cost-benefit analysis is based on the success of the method in a long-term assessment over 6 years, reaching an adult, self-sustaining population.
Our attempt to operate over large spatial scales required the cultivation of a high number of individuals (approximately 325,000 at the end of the nursery phase), and raised several methodological issues during the outplanting process. Processes deserving attention and still requiring further research when approaching Cystoseira restoration are summarized in Fig 4. Availability of fertile receptacles influences the timing of both ex situ and in situ restoration, as the suitable season might be limited to 2–3 months per year and varies according to the species, to climatic conditions, and to the geographical position of the donor location. The availability of appropriate facilities limits the number of individuals grown, and excessive densities of germlings in the aquaria may result in negative density-dependent effects such as competition or blooms of epiphytes or pathogenic agents (e.g. bacteria). Species-specific culturing conditions, optimal duration of the cultivation phase, and the appropriate substratum need to be identified with proper experimentation, similarly to the approach described by Falace et al. .
One of the most delicate phases is transport from the nursery and storage location to the field, as juveniles are susceptible to variations in abiotic conditions and to mechanical stress. Also, timeliness in the outplanting process is crucial and guidelines for the handling of germlings in the delicate gap between the nursery and the field could be further implemented with new strategies. For example, the temporal desynchronization of germlings cultivation could provide a manageable number of germlings to be outplanted in several subsequent occasions, possibly minimizing germlings mortality. Alternatively, areas for the storage and acclimation of germlings prior to attachment in the field could be individuated in enclosed coastal areas, eventually allowing for a longer time of cultivation in semi-controlled conditions, in order to minimize size-dependent mortality due to stressful conditions (e.g. abiotic and mesograzing pressures). Once in the field, the presence of adult conspecifics does not represent a necessary pre-requisite for the survival of C. amentacea germlings. In contrast, the mitigation of grazing pressure is crucial and could be implemented by setting up protective cages, although these structures cannot prevent the access to mesograzers and their application could be critical over vast extensions (high costs, high probability of dislodgement by waves, use of plastic, esthetic injury). Pilot studies carried out on putative restoration sites should be considered a priority to individuate promising locations according to abiotic (e.g. hydrodynamic conditions) and biotic criteria (e.g. frequentation by macrograzers, structure of mesograzers assemblage). Yet, climate change is expected to severely affect the recruitment and survival of Cystoseira germlings, either directly or indirectly (e.g. increasing grazing rates), eventually setting new priorities for the success of restoration interventions. For example, the negligible presence of Cystoseira adults might become a favourable, necessary condition under warming climate . Also, macrograzers pressure is expected to increase , although no specific predictions have been formulated for the activity of mesograzers.
Taken together, our analyses show that large-scale restoration is possible only if basic information about species, protocols and the areas to be restored is carefully addressed. At the moment, specific details about ecophysiology are available only for a few Cystoseira species (C. amentacea var. stricta, [26, 27, 33]; C. barbata, [9, 18, 25, 28, 29, 42, 60, 77]; C. compressa, [25–27, 77]; C. zosteroides, [63, 64, 68, 87]) out of 45 Cystoseira species existing in the Mediterranean Sea , limiting the potential of selecting suitable species in a restoration framework. Although the Mediterranean Action Plan, adopted within the framework of the Barcelona Convention (1976), identifies, in an amendment of 2009 (Annex IV, SPA/BD Protocol—United Nations Environment Programme [UNEP]), the conservation of all but one Mediterranean Cystoseira species (C. compressa) as a priority, most Cystoseira forests are under severe threat.
Restoring areas depleted by macroalgal forests is an important target given the shifts that we are documenting at global scales. However, at this stage, a priority should be to address efforts and resources devoted to this activity would be better spent in mitigating local stressors and preserving existing forests. A better coordination among complementary strategies of mitigation, conservation and restoration seems to be still the best strategy to halt present trajectories of changes.
1. Halpern BS, Walbridge S, Selkoe KA, Kappel CV, Micheli F, D'Agrosa C, et al. A global map of human impact on marine ecosystems. Science 2008; 319: 948–952. doi: 10.1126/science.1149345 18276889
2. Filbee-Dexter K, Scheibling RE. Sea urchin barrens as alternative stable states of collapsed kelp ecosystems. Mar Ecol Prog Ser. 2014; 495: 1–25. doi: 10.3354/meps10573
3. Lotze HK, Lenihan HS, Bourque BJ, Bradbury RH, Cooke RG, Kay MC, et al. Depletion, degradation, and recovery potential of estuaries and coastal seas. Science 2006; 312: 1806–1809. doi: 10.1126/science.1128035 16794081
4. Lotze HK, Coll M, Magera AM, Ward-Paige C, Airoldi L. Recovery of marine animal populations and ecosystems. Trends Ecol Evol. 2011; 26: 595–605. doi: 10.1016/j.tree.2011.07.008 21852017
5. Abelson A, Halpern BS, Reed DC, Orth RJ, Kendrick GA, Beck MW, et al. Upgrading marine ecosystem restoration using ecological-social concepts. Bioscience 2016; 66: 156–163. doi: 10.1093/biosci/biv171 26977115
6. Bayraktarov E, Saunders MI, Abdullah S, Mills M, Behr J, Possingham HP, et al. The cost and feasibility of marine coastal restoration. Ecol Appl. 2016; 26: 1055–1074. doi: 10.1890/15-1077 27509748
7. Thibaut T, Pinedo S, Torras X, Ballesteros E. Long-term decline of the populations of Fucales (Cystoseira spp. and Sargassum spp.) in the Alberes coast (France, North-western Mediterranean). Mar Pollut Bull. 2005; 50: 1472–1489. doi: 10.1016/j.marpolbul.2005.06.014 16026805
8. Airoldi L, Beck MW. Loss, status and trends for coastal marine habitats of Europe. Oceanogr Mar Biol. 2007; 45: 347–407. doi: 10.1201/9781420050943.ch7
9. 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 S. 2011; 92: 347–357. doi: 10.1016/j.ecss.2011.01.008
10. Steneck RS, Graham MH, Borque BJ, Corbett D, Erlandson JM, Estes JA, et al. Kelp forest ecosystems, biodiversity, stability, resilience and future. Env Cons. 2002; 29: 436–459. doi: 10.1017/S0376892902000322
11. Martínez EA, Cárdenas L. Recovery and genetic diversity of the intertidal kelp Lessonia nigrescens (Phaeophyceae) 20 years after El Niño. J Phycol. 2003; 39: 504–508. doi: 10.1046/j.1529-8817.2003.02191.x
12. Moy FE, Christie H. Large-scale shift from sugar kelp (Saccharina latissima) to ephemeral algae along the south and west coast of Norway. Mar Biol Res. 2012; 8: 309–321 doi: 10.1080/17451000.2011.637561
13. Filbee-Dexter K, Feehan CJ, Scheibling RE. Large-scale degradation of a kelp ecosystem in an ocean warming hotspot. Mar Ecol Prog Ser. 2016; 543: 141–152. doi: 10.3354/meps11554
14. Connell SD, Russell BD, Turner DJ, Shepherd SA, Kildea T, Miller D, Airoldi L, Cheshire A. Recovering a lost baseline: missing kelp forests from a metropolitan coast. Mar Ecol Prog Ser. 2008; 360: 63–72. doi: 10.3354/meps07526
15. Smale DA, Wernberg T. Extreme climatic event drives range contraction of a habitat-forming species. P R Soc B 2013; 280: 20122829. doi: 10.1098/rspb.2012.2829 23325774
16. Bianchi C, Corsini-Foka M, Morri C, Zenetos A. Thirty years after—dramatic change in the coastal marine habitats of Kos Island (Greece), 1981–2013. Mediterr Mar Sci. 2014; 15(3): 482–497. doi: 10.12681/mms.678
17. Thibaut T, Blanfuné A, Boudouresque CF, Verlaque M. Decline and local extinction of Fucales in the French Riviera: the harbinger of future extinctions? Mediterr Mar Sci. 2015; 16/1: 206–224. doi: 10.12681/mms.1032
18. Perkol-Finkel S, Airoldi A. Loss and recovery potential of marine habitats: an experimental study of factors maintaining resilience in subtidal algal forests at the Adriatic Sea. PLoS ONE 2010; 5(5): e10791. doi: 10.1371/journal.pone.0010791 20520726
19. Iveša L, Djakovac T, Devescovi M. Long-term fluctuations in Cystoseira populations along the west Istrian Coast (Croatia) related to eutrophication patterns in the northern Adriatic Sea. Mar Pollut Bull. 2016; 106: 162–173. doi: 10.1016/j.marpolbul.2016.03.010 26975612
20. Gianni F, Bartolini F, Airoldi L, Ballesteros E, Francour P, Guidetti P, et al. Conservation and restoration of marine forests in the Mediterranean Sea and the potential role of Marine Protected Areas. AIOL 2013; 4: 83–101. doi: 10.1080/19475721.2013.845604
21. Jacob C, Buffard A, Pioch S, Thorin S. Marine ecosystem restoration and biodiversity offset. Ecol Eng. 2018; 120: 595–594. doi: 10.1016/j.ecoleng.2017.09.007
22. Correa JA, Lagos NA, Medina MH, Castilla JC, Cerda M, Ramírez M, et al. Experimental transplants of the large kelp Lessonia nigrescens (Phaeophyceae) in high-energy wave exposed rocky intertidal habitats of northern Chile: experimental restoration and management applications. J Exp Mar Biol Ecol. 2006; 335: 13–18. doi: 10.1016/j.jembe.2006.02.010
23. Campbell AH, Marzinelli EM, Vergés A, Coleman MA, Steinberg PD. Towards restoration of missing underwater forests. PLoS ONE 2014; 9(1):e84106. doi: 10.1371/journal.pone.0084106 24416198
24. Westermeir R, Murúa P, Patiño DJ, Muñoz L, Müller DG. Holdfast fragmentation of Macrocystis pyrifera (integrifolia morph) and Lessonia berteroana in Atacama (Chile): a novel approach for kelp bed restoration. J Appl Phycol. 2016; 28: 2969–2977. doi: 10.1007/s10811-016-0827-2
25. Falace A, Zanelli E, Bressan G. Algal transplantation as a potential tool for artificial reef management and environmental mitigation. Bull Mar Sci. 2006; 78: 161–166.
26. Susini ML, Mangialajo L, Thibaut T, Meinesz A. Development of a transplantation technique of Cystoseira amentacea var. stricta and Cystoseira compressa. Hydrobiologia 2007; 580: 241–244. doi: 10.1007/978-1-4020-6156-1_21
27. Mangialajo L, Chiantore M, Susini ML, Meinesz A, Cattaneo-Vietti R, Thibaut T. Zonation patterns and interspecific relationships of fucoids in microtidal environments. J Exp Mar Biol Ecol. 2012; 412: 72–80. doi: 10.1016/j.jembe.2011.10.031
28. Perkol-Finkel S, Ferrario F, Nicotera V, Airoldi L. Conservation challenges in urban seascapes: promoting the growth of threatened species on coastal infrastructures. J Appl Ecol. 2012: 49: 1457–1466. doi: 10.1111/j.1365-2664.2012.02204.x
29. Robvieux P. Conservation des populations de Cystoseira en regions Provence-Alpes-Côte-d’Azur et Corse. 2013 PhD dissertation, Université de Nice-Sophia Antipolis, France
30. De La Fuente G, Chiantore M, Asnaghi V, Kaleb S, Falace A. First ex situ outplanting of the habitat-forming seaweed Cystoseira amentacea var. stricta from a restoration perspective. PeerJ 2019; 7:e7290. doi: 10.7717/peerj.7290 31367482
31. Piazzi L, Ceccherelli G. Concomitance of oligotrophy and low grazing pressure is essential for the resilience of Mediterranean subtidal forests. Mar Pollut Bull. 2017; 123: 197–204. doi: 10.1016/j.marpolbul.2017.08.061 28886921
32. Rindi L, Dal Bello M, Benedetti-Cecchi L. Experimental evidence of spatial signatures of approaching regime shifts in macroalgal canopies. Ecology 2018; 99: 1709–1715. doi: 10.1002/ecy.2391 29797316
33. Falace A, Kaleb S, De La Fuente G, Asnaghi V, Chiantore M. Ex situ cultivation protocol for Cystoseira amentacea var. stricta (Fucales, Phaeophyceae) from a restoration perspective. PLoS ONE 2018; 13(2): e0193011. doi: 10.1371/journal.pone.0193011 29447238
34. Gianni F, Bartolini F, Airoldi L, Mangialajo L. Reduction of herbivorous fish pressure can facilitate focal algal species forestation on artificial structures. Mar Environ Res. 2018; 138: 102–109. doi: 10.1016/j.marenvres.2018.04.007 29706367
35. Hereu B. Depletion of palatable algae by sea urchins and fishes in a Mediterranean subtidal community. Mar Ecol Prog Ser 2006; 313: 95–103. doi: 10.3354/meps313095
36. Agnetta D, Badalamenti F, Ceccherelli G, Di Trapani F, Bonaviri C, Gianguzza P. Role of two co-occurring Mediterranean Sea urchins in the formation of barren from Cystoseira canopy. Estuar Coast Shelf S. 2015; 152: 73–77. doi: 10.1016/j.ecss.2014.11.023
37. Gianni F, Bartolini F, Pey A, Laurent M, Martins GM, Airoldi L, et al. Threats to large brown algal forests in temperate seas: the overlooked role of native herbivorous fish. Sci Rep. 2017; 7: 6012. doi: 10.1038/s41598-017-06394-7 28729633
38. Benedetti-Cecchi L, Cinelli F. Effects of canopy cover, herbivores and substratum type on patterns of Cystoseira spp. settlement and recruitment in littoral rockpools. Mar Ecol Prog Ser. 1992; 90: 183–191. doi: 10.3354/meps090183
39. Duffy JE, Hay ME. Strong impacts of grazing amphipods on the organization of a benthic community. Ecol Monogr 2000; 70:237–263. doi: 10.1890/0012-9615(2000)070[0237:SIOGAO]2.0.CO;2
40. Vadas RL, Johnson S, Norton TA. Recruitment and mortality of early post-settlement stages of benthic algae. Brit Phycol J. 1992; 27: 331–351. doi: 10.1080/00071619200650291
41. Schiel DR, Foster MS. The population biology of large brown seaweeds. Ecological consequences of multiphase life histories in dynamic coastal environments. Annu Rev Ecol Evol S. 2006; 37: 343–72. doi: 10.1146/annurev.ecolsys.37.091305.110251
42. Irving AD, Balata D, Colosio F, Ferrando GA, Airoldi L. Light, sediment, temperature, and the early life-history of the habitat forming alga Cystoseira barbata. Mar Biol. 2009; 156: 1223–1231. doi: 10.1007/s00227-009-1164-7
43. Piazzi L, Balata D, Ceccherelli G. Resilience of Cystoseira beds: lack of canopy effects on recruitment. Bot Mar. 2017; 60: 49–53. doi: 10.1515/bot-2016-0088
44. Buonomo R, Chefaoui RM, Bermejo Lacida R, Engelen AH, Serrão EA, Airoldi L. Predicted extinction of unique genetic diversity in marine forests of Cystoseira spp. Mar Env Res. 2018; 138: 119–128. doi: 10.1016/j.marenvres.2018.04.013 29716751
45. Buonomo R, Assis J, Fernandes F, Engelen AH, Airoldi L, Serrão EA. Habitat continuity and stepping-stone oceanographic distances explain population genetic connectivity of the brown alga Cystoseira amentacea. Mol. Ecol. 2017; 26: 766–780. doi: 10.1111/mec.13960 27997043
46. Fraschetti S, Terlizzi A, Bussotti S, Guarnieri G, D’Ambrosio P, Boero F. Conservation of Mediterranean seascapes: analyses of existing protection schemes. Mar Environ Res. 2005; 59: 309–332. doi: 10.1016/j.marenvres.2004.05.007 15589984
47. Fraschetti S, Guarnieri G, Bevilacqua S, Terlizzi A, Boero F. Protection enhances community and habitat stability: evidence from a Mediterranean Marine Protected Area. PLoS ONE 2013; 8(12): e81838. doi: 10.1371/journal.pone.0081838 24349135
48. Guarnieri G, Bevilacqua S, De Leo F, Farella G, Maffia A, Terlizzi A, et al. The challenge of planning conservation strategies in threatened seascapes: understanding the role of fine scale assessments of community response to cumulative human pressures. PLoS ONE 2016; 11(2): e0149253. doi: 10.1371/journal.pone.0149253 26871942
49. Dethier MN, Graham ES, Cohen S, Tear LM. Visual versus random-point percent cover estimations: “objective” is not always better. Mar Ecol Progr Ser. 1993; 96: 93–100.
50. Bolker BM, Brooks ME, Clark CJ, Geange SW, Poulsen JR, Stevens MHH, White JSS. Generalize linear mixed models: a practical guide for ecology and evolution. Trends Ecol Evol. 2009; 24(3): 127–135. doi: 10.1016/j.tree.2008.10.008 19185386
51. Burnham KP, Anderson DR. Model Selection and Multimodel Inference: A Practical Information-Theoretic Approach. 2002 Springer-Verlag
52. Zuur AF, Ieno EN, Smith GM. Analysing Ecological Data. 2007 Springer-Verlag
53. Crawley MJ. The R Book. 2012 Chichester, West Sussex, United Kingdom, Wiley.
54. Bates D, Maechler M, Bolker B, Walker S. Fitting linear mixed-effects models using lme4. J Stat Softw. 2015; 67(1):1–48. doi: 10.18637/jss.v067.i01
55. Sandrini-Neto L, Camargo MG. GAD: an R package for ANOVA designs from general principles. 2014; Available on CRAN.
56. Anderson MJ, Gorley RN, Clarke KR PERMANOVA+ for PRIMER: guide to software and statistical methods. 2008; PRIMER-E, Plymouth, UK 214 pp.
57. Anderson MJ. Permutational Multivariate Analysis of Variance (PERMANOVA). 2017 Wiley StatsRef doi: 10.1002/9781118445112.stat07841
58. Bray JR, Curtis JT. An ordination of the upland forest communities of Southern Wisconsin. Ecol. Monogr. 1957; 27: 325–349. doi: 10.2307/1942268
59. Carney LT, Waaland JR, Klinger T, Ewing K. Restoration of the bull kelp Nereocystis luetkeana in nearshore rocky habitats. Mar Ecol Prog Ser. 2005; 302: 49–61. doi: 10.3354/meps302049
60. Verdura J, Sales M, Ballesteros E, Cefalì ME, Cebrian E. Restoration of a canopy-forming alga based on recruitment enhancement: methods and long-term success assessment. Front Plant Sci. 2018; 9: 1832. doi: 10.3389/fpls.2018.01832 30619405
61. Gunnill FC. Demography of the intertidal brown alga Pelvetia fastigiata in southern California, USA. Mar Biol. 1980; 59: 169–179. doi: 10.1007/BF00396865
62. Aberg P, Pavia H. Temporal and multiple scale spatial variation in juvenile and adult abundance of the brown alga Ascophyllum nodosum. Mar Ecol Prog Ser. 1997; 158: 111–119. doi: 10.3354/meps158111
63. Capdevila P, Linares C, Aspillaga E, Navarro L, Kersting DK, Hereu B. Recruitment patterns in the Mediterranean deep-water alga Cystoseira zosteroides. Mar Biol. 2015; 162: 1165–1174. doi: 10.1007/s00227-015-2658-0
64. Capdevila P, Linares C, Aspillaga E, Riera JL, Hereu B. Effective dispersal and density-dependence in mesophotic macroalgal forests: insights from the Mediterranean species Cystoseira zosteroides. PLoS ONE 2018; 13(1): e0191346. doi: 10.1371/journal.pone.0191346 29329336
65. Steen H, Scrosati R. Intraspecific competition in Fucus serratus and F. evanescens (Phaeophyceae: Fucales) germlings: effects of settlement density, nutrient concentration, and temperature. Mar Biol. 2004; 144: 61–70. doi: 10.1007/s00227-003-1175-8
66. Lind AC, Konar B. Effects of abiotic stressors on kelp early life-history stages. 2017; Algae 32: 223–233. doi: 10.4490/algae.2017.32.8.7
67. Andrews S, Bennett S, Wernberg T. Reproductive seasonality and early life temperature sensitivity reflect vulnerability of a seaweed undergoing range reduction. 2014; Mar Ecol Prog Ser. 495: 119–129. doi: 10.3354/meps10567
68. Capdevila P, Hereu B, Salguero-Gómez R, la Rovira G, Medrano A, Cebrian E, Garrabou J, Kersting DK, Linares C. Warming impacts on early life stages increase the vulnerability and delay the population recovery of a long-lived habitat-forming macroalga. 2018; J Ecol. 1–12. doi: 10.1111/1365-2745.13090
69. Bennett S, Wernberg T, de Bettignies T, Kendrick GA, Anderson RJ, Bolton JJ, Rodgers KL, Shears NT, Leclerc JC, Leveque L, Davoult D, Christie HC. Canopy interactions and physical stress gradients in subtidal communities. 2015, Ecol Lett. 18: 677–686. doi: 10.1111/ele.12446 25975532
70. Johnson LE, Brawley SH. Dispersal and recruitment of a canopy-forming intertidal alga: the relative roles of propagule availability and post-settlement processes. Oecologia 1998; 117: 517–526. doi: 10.1007/s004420050688 28307677
71. Konar B, Estes JA. The stability of boundary regions between kelp beds and deforested areas. Ecology 2003; 84: 174–185. doi: 10.1890/0012-9658(2003)084[0174:TSOBRB]2.0.CO;2
72. Taylor DI, Schiel DR. Algal populations controlled by fish herbivory across a wave exposure gradient on southern temperate shores. Ecology 2010; 91: 201–211. doi: 10.1890/08-1512.1 20380209
73. Toohey B, Kendrick G, Wernberg T, Phillips J, Malkin S, Prince J. The effects of light and thallus scour from Ecklonia radiata canopy on an associated foliose algal assemblage: the importance of photoacclimation. Mar Biol. 2004; 144: 1019–1027. doi: 10.1007/s00227-003-1267-5
74. Wernberg T, Kendrick GA, Toohey BD. Modification of the physical environment by an Ecklonia radiata (Laminariales) canopy and implications for associated foliose algae. Aquat Ecol. 2005; 39: 419–430. doi: 10.1007/s10452-005-9009-z
75. Haring RN, Dethier MN, Williams SL. Desiccation facilitates wave-induced mortality of the intertidal alga Fucus gardneri. Mar Ecol Prog Ser. 2002; 232: 75–82. doi: 10.3354/meps232075
76. Kendrick GA. Effects of propagule settlement density and adult canopy on survival of recruits of Sargassum spp. (Sargassaceae: Phaeophyta). Mar Ecol Prog Ser. 1994; 103: 129–140. doi: 10.3354/meps103129
77. Ferrario F, Iveša L, Jaklin A, Perkol-Finkel S, Airoldi L. The overlooked role of biotic factors in controlling the ecological performance of artificial marine habitats. J Appl Ecol. 2016; 53: 16–24. doi: 10.1111/1365-2664.12533
78. Benedetti-Cecchi L, Cinelli F. Habitat heterogeneity, sea urchins grazing and the distribution of algae in littoral rock pools on the west coast of Italy (western Mediterranean). Mar Ecol Prog Ser. 1995; 126: 203–212. doi: 10.3354/meps126203
79. Bonaviri C, Vega Fernández T, Fanelli G, Badalamenti F, Gianguzza P. Leading role of the sea urchin Arbacia lixula in maintaining the barren state in southwesters Mediterranean. Mar Biol. 2011; 158: 2505–2513. doi: 10.1007/s00227-011-1751-2
80. Johnson LE, Paine RT. Consistency in a marine algal-grazer interaction over multiple scales. J Phycol. 2016; 52: 942–950. doi: 10.1111/jpy.12475 27711961
81. Zupo V, Alexander T, Edgar G. Relating trophic resources to community structure: a predictive index of food availability. R Soc Open Sci. 2017; 4: 160515. doi: 10.1098/rsos.160515 28386417
82. Sjøtun K, Eggereide SF, Høisaeter T. Grazer-controlled recruitment of the introduced Sargassum muticum (Phaeophyceae, Fucales) in northern Europe. Mar Ecol Prog Ser. 2007; 342: 127–138. doi: 10.3354/meps342127
83. Korpinen S, Jormalainen V. Grazing and nutrients reduce recruitment success of Fucus vesiculosus L. (Fucales: Phaeophyceae). Estuar Coast Shelf S. 2008; 78: 437–444. doi: 10.1016/j.ecss.2008.01.005
84. Henríquez LA, Buschmann AH, Maldonaldo MA, Graham MH, Hernández-González MC, Pereda SV, et al. Grazing on giant kelp microscopic phases and the recruitment success of annual populations of Macrocystis pyrifera (Laminariales, Phaeophyta) in southern Chile. J Phycol. 2011; 47: 252–258. doi: 10.1111/j.1529-8817.2010.00955.x 27021857
85. Chenelot H, Konar B. Lacuna vincta (Mollusca, Neotaenioglossa) herbivory on juvenile and adult Nereocystis luetkeana (Heterokontophyla, Laminariales). Hydrobiologia 2007; 583:107–118. doi: 10.1007/s10750-006-0484-6
86. Vergés A, Doropopoulos C, Malcolm HA, Skye M, Garcia-Pizá M, Marzinelli EM, Campbell AH, Ballesteros E, Hoey AS, Villa-Conceio A, Bozec YM, Steinberg PD. Long-term empirical evidence of ocean warming leading to tropicalization of fish communities, cincreased herbivory and loss of kelp. Proc Natl Acad Sci USA. 2016; 113: 13791–13796. doi: 10.1073/pnas.1610725113 27849585
87. Ballesteros E, Garrabou J, Hereu B, Zabala M, Cebrian E, Sala E. Deep-water stands of Cystoseira zosteroides C. Agardh (Fucales, Ochrophyta) in the Northwestern Mediterranean: insights into assemblage structure and population dynamics. Estuar Coast Shelf S. 2009; 82: 477–484. doi: 10.1016/j.ecss.2009.02.013
88. Guiry MD, Guiry GM. AlgaeBase. Worldwide electronic publication, National University of Ireland, Galway. http://www.algaebase.org; searched on 19 December 2018