Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Contrasts in the marine ecosystem of two Macaronesian islands: A comparison between the remote Selvagens Reserve and Madeira Island

  • Alan M. Friedlander ,

    Roles Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Project administration, Resources, Software, Supervision, Validation, Visualization, Writing – original draft, Writing – review & editing

    alan.friedlander@hawaii.edu

    Affiliations Pristine Seas, National Geographic Society, Washington, DC, United States of America, Fisheries Ecology Research Lab, University of Hawai‘i at Mānoa, Honolulu, Hawai‘i, United States of America

  • Enric Ballesteros,

    Roles Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Resources, Software, Validation, Visualization, Writing – original draft, Writing – review & editing

    Affiliation Departamento Ecología Marina, Centre d'Estudis Avançats de Blanes-CSIC, Blanes, Spain

  • Sabrina Clemente,

    Roles Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Resources, Software, Validation, Visualization, Writing – original draft, Writing – review & editing

    Affiliations Departamento Ecología Marina, Centre d'Estudis Avançats de Blanes-CSIC, Blanes, Spain, Departamento de Biología Animal, Universidad de La Laguna, Tenerife, Canary Islands, Spain

  • Emanuel J. Gonçalves,

    Roles Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Resources, Software, Validation, Visualization, Writing – original draft, Writing – review & editing

    Affiliation MARE – Marine and Environmental Sciences Centre, ISPA – Instituto Universitário, Lisbon, Portugal

  • Andrew Estep,

    Roles Conceptualization, Investigation, Methodology, Resources, Software, Writing – review & editing

    Affiliation Waitt Institute, La Jolla, California, United States of America

  • Paul Rose,

    Roles Conceptualization, Investigation, Project administration, Supervision

    Affiliations Pristine Seas, National Geographic Society, Washington, DC, United States of America, Royal Geographical Society, London, United Kingdom

  • Enric Sala

    Roles Conceptualization, Funding acquisition, Investigation, Project administration, Resources, Supervision, Validation, Writing – original draft, Writing – review & editing

    Affiliation Pristine Seas, National Geographic Society, Washington, DC, United States of America

Abstract

The islands of Madeira and Selvagens are less than 300 km apart but offer a clear contrast between a densely populated and highly developed island (Madeira), and a largely uninhabited and remote archipelago (Selvagens) within Macaronesia in the eastern Atlantic. The Madeira Archipelago has ~260,000 inhabitants and receives over six million visitor days annually. The Selvagens Islands Reserve is one of the oldest nature reserves in Portugal and comprises two islands and several islets, including the surrounding shelf to a depth of 200 m. Only reserve rangers and a small unit of the maritime police inhabit these islands. The benthic community around Selvagens was dominated by erect and turf algae, while the community at Madeira was comprised of crustose coralline and turf algae, sessile invertebrates, and sea urchin barrens. The sea urchin Diadema africanum was 65% more abundant at Madeira than at Selvagens. Total fish biomass was 3.2 times larger at Selvagens than at Madeira, and biomass of top predators was more than 10 times larger at Selvagens. Several commercially important species (e.g., groupers, jacks), which have been overfished throughout the region, were more common and of larger size at Selvagens than at Madeira. Important sea urchin predators (e.g., hogfishes, triggerfishes) were also in higher abundance at Selvagens compared to Madeira. The effects of fishing and other anthropogenic influences are evident around Madeira. This is in stark contrast to Selvagens, which harbors healthy benthic communities with diverse algal assemblages and high fish biomass, including an abundance of large commercially important species. The clear differences between these two island groups highlights the importance of expanding and strengthening the protection around Selvagens, which harbors one of the last intact marine ecosystems in the North Atlantic, and the need to increase management and protection around Madeira.

Introduction

Macaronesia is a collection of four archipelagos (Madeira, Selvagens, Azores, and Canaries) located in the North Atlantic Ocean off the coasts of Europe and Africa [13], although Cape Verde is also included when considering terrestrial ecosystems [4]. The Macaronesian biogeographic region has long been noted for its terrestrial biodiversity and high endemism [5]. It also has a unique marine fauna, which has been influenced by west Africa, the Mediterranean Sea, and continental western Europe [13, 67].

Madeira is a Portuguese archipelago within Macaronesia and is situated ~ 650 km west of Morocco at its closest point (Fig 1). The archipelago includes the islands of Madeira (pop. 262,456), Porto Santo (pop. 5,483), and the Desertas, which consists of three uninhabited islands: Deserta Grande, Bugio and Ilhéu de Chão. The island of Madeira is a large shield volcano that rises > 6,000 m and is the largest island of the group (741 km2, [8]). It is home to Funchal, which is the capital and principal city of the Autonomous Region of Madeira. Tourism is the most important sector of the economy of the archipelago, with > 6.5 million overnight stays reported in 2014 [9].

thumbnail
Fig 1. Sampling locations around Madeira and Selvagens in the eastern Atlantic Ocean.

https://doi.org/10.1371/journal.pone.0187935.g001

The Selvagens Islands are located ~ 290 km southeast of Madeira and 170 km north of Tenerife, Canary Islands. They consist of two major islands (Selvagem Grande and Selvagem Pequena), and several small islets of varying sizes. The islands are administered by the Portuguese municipality of Funchal, and are part of the Madeira civil parish of Sé. In 1971, the Portuguese government acquired these islands, converting them into a strict nature reserve [10]. The Selvagens Islands Nature Reserve includes the surrounding shelf to a depth of 200 m (total marine area—34 km2) and is one of the oldest nature reserves in Portugal [10]. On Selvagem Grande, there is a permanent station with two year-around rangers and two members of the maritime police, while Selvagem Pequena is typically manned by two rangers between May and October. These are the only human inhabitants on these islands. The Selvagens Islands Nature Reserve is listed as an IUCN Protected Area Category Ia. Human activity is strictly controlled, with visits on land restricted to paths and a small area at the administrative facility. The average number of permits per year to visit the reserve is 122, involving ~ 550 people. The Selvagens are protected under European Commission legislation, included within the Natura 2000 Network as a Special Conservation Area, and designated a Special Protection Area (SPA) under the Birds Directive of the European Union [11]. It is the world’s largest breeding area for Cory’s Shearwaters (Calonectris diomedea, [12]), and is important to a number of other endangered and threatened sea bird species that are known to breed there [13].

Since its discovery and colonization in the 15th century, Madeira has had a long history of fishing [14]. The narrow and steep insular platforms, along with low productivity waters, limit the available habitat for demersal fisheries species, and restricts the fishing methods that are suitable for these habitats [15]. Intensive overfishing of many prized nearshore species has been a concern throughout the region for decades [1618]. This has not only reduced nearshore fish stocks, but has also had a cascading effect on the entire ecosystem by removing the predators of sea urchins, which have proliferated in the absence of these predators. This hyper-abundance of sea urchins has subsequently grazed down the erect vegetative framework of the reefs, creating barrens that have further reduced the productivity of nearshore ecosystems [1921].

In contrast to most locations in the region, the Selvagens Islands are remote and nearly uninhabited, allowing for comparisons of marine ecosystems between two distinct levels of anthropogenic influence, while still within the same biogeographic region. Because of its remoteness and strict protection, the Selvagens Islands may represent one of the last remaining intact marine ecosystems in the eastern Atlantic. With this in mind, we set out to examine the value of Selvagens as a baseline for ecosystem health that can be used to assess ecosystem state not only within the region, but also globally. Our main objective was to compare the benthic communities and fish assemblages around a densely populated and highly developed island (Madeira) with the Selvagens Islands Nature Reserve, a strict marine protected area with almost no direct human impacts. An additional goal of our research was to examine the efficacy of marine protected areas around Madeira, comparing them to the Selvagens Islands Nature Reserve, and fished areas around Madeira.

Materials and methods

Ethics statement

Data were collected by all authors in a collaborative effort. Non-invasive research was conducted, which included photographs and visual estimates described in the methods. The Serviço do Parque Natural da Madeira granted all necessary permission to conduct this research. No vertebrate sampling was conducted and therefore no approval was required by the Institutional Animal Care and Use Committee. Our data are available at Data Dryad: doi:10.5061/dryad.322q2.

Sampling design

We conducted in situ surveys of fishes, algae, and benthic macro-invertebrates within two depth strata (10 and 20 m) at 29 locations at Selvagens in September 2015 and 36 locations around Madeira in July 2016 (Fig 1). Oceanographic conditions (e.g., sea surface temperature, Chlorophyll-a) are similar around both locations during the time periods when sampling was conducted [22], thus reducing the effects of seasonality. Previous studies in the region have not shown seasonal differences in reef fish abundance [16, 23], and we therefore feel that our sampling effort in not confounded by potential temporal variability in these reef fish assemblages. Sampling locations were restricted to rocky boulder habitat at both locations to reduce the influence of habitat on fish assemblage structure, which is known to exist around Madeira [24]. Samples were allocated haphazardly between the 10 and 20 m bathymetric contour, with sites spaced ~ 1 km apart around Selvagens and ~ 2.5 km apart around Madeira. Weather conditions precluded sampling along parts of the north coast of Madeira.

At Madeira Island, the Garajau Marine Protected Area (MPA) prohibits the take of marine life except bait fish for the tuna fishery, which consists of small coastal pelagic species [e.g., mainly Atlantic chub mackerel (Scomber colias), European pilchard (Sardina pilchardus) and bogue (Boops boops)]. The Rocha do Navio MPA permits extraction of nearshore resources using spearfishing and line fishing, as well as collecting of invertebrates and tuna bait fishing. Samples from the Garajau MPA were excluded in overall comparisons between Madeira and Selvagens as it was the only MPA with a high level of protection at Madeira. The Garajau MPA accounted for only 4.4% of the coastline around Madeira and its exclusion allowed for a better comparison of the two island groups without the confounding effect of this small (3.76 km2) MPA.

All surveys were conducted on open-circuit scuba between the hours of 9:00 and 16:00 to reduce the influence of crepuscular variation in the fish assemblages. Dive duration ranged from 60–90 min depending on habitat and environmental conditions.

Benthos

Characterization of the benthos was conducted along one 50 m-long transect run parallel to the shoreline at each of the two depth strata. For algae and sessile invertebrates, we used a line-point intercept methodology along each transect, recording the species or taxa found every 20 cm on a fiberglass measuring tape. Sessile benthic cover was grouped into turf algae (< 3 cm in height, e.g., Caulerpa webbiana, Cottoniella filamentosa), erect algae (e.g., Halopteris scoparia, Padina pavonica), canopy algae (e.g., Cystoseira abies-marina, Sargassum spp.), erect non-canopy algae (e.g., Lobophora variegata), crustose coralline algae (CCA, e.g., Neogoniolithon spp.), sessile invertebrates (e.g., Balanus sp.), and barrens, which consisted of bare rock. Sponges accounted for <5% of the sessile cover at both island groups and were included as sessile invertebrates. For mobile invertebrates, we counted individuals in twenty-five 50 x 50 cm quadrats randomly placed along each of the 50-m transects.

Fishes

At each depth stratum within a site, divers counted and estimated lengths of all fishes encountered within fixed-length (25-m) belt transects whose widths differed depending on the direction of swim. All fish ≥ 20 cm total length (TL) were tallied within a 4-m wide strip surveyed on an initial “swim-out” as the transect line was laid out (transect area = 100 m2). All fishes < 20 cm TL were tallied within a 2-m wide strip surveyed on the return swim back along the laid transect line (transect area = 50 m2). Swimming duration per transect varied from 10–15 min, depending on habitat complexity and fish abundance. Three replicate transects were performed per site at each depth stratum.

Fishes were identified to species level in all cases. Fish total length (TL) was estimated to the nearest cm. Fishes were tallied by length and individual-specific lengths were converted to body weights. Numerical density (abundance) was expressed as number of individuals per m2 and biomass density was expressed as g per m2. The biomass of individual fishes was estimated using the allometric length-weight conversion: W = aTLb, where parameters a and b are species-specific constants, TL is total length in cm, and W is weight in grams. Length-weight fitting parameters were obtained from FishBase [25]. The sum of all individual weights and numerical densities was used to estimate biomass and numerical density by species. Fish species diversity were calculated from the Shannon-Weaver diversity index: , where pi is the proportion of all individuals counted that were of taxa i. Fishes were categorized into four trophic groups (top predators, herbivores, secondary consumers, and planktivores) based on information from FishBase [25]. Resource (commercially targeted) species were designated based on expert opinion of dive operators, fishers, and resource managers.

Statistical analysis

Principal Coordinate Analysis (PCO) was used to compare sessile benthic functional groups between islands. Data were arcsine square root transformed prior to analysis. Drivers of sessile benthic community structure were investigated using permutation-based multivariate analysis of variance (PERMANOVA). Similarity percentages analysis (SIMPER) was used to examine differences in sessile benthic functional group cover between islands and depths. PCO, SIMPER, and PERMANOVA were also used to compare mobile invertebrates and fishes between island groups. Densities (individuals m-2) of the top 10 mobile invertebrate species, which together accounted for > 90% of the individuals in each island group, were square root transformed prior to all analyses. Fish assemblage structure by biomass (g m-2) was square root transformed, while numerical abundance (individuals m-2) was ln(x+1) transformed prior to multivariate analyses. All PERMANOVA, PCOs, and SIMPER analyses were conducted using Primer v6 [26].

Fish assemblage characteristics and trophic biomass between islands were compared using linear mixed models (LMMs) with island group, depth strata, and their interaction treated as fixed factors. Stations were treated as a random effect to account for spatial autocorrelation in these data. Unplanned comparisons between pairs were examined using the Tukey-Kramer HSD (honestly significant difference) test (α = 0.05). Numerical abundance, total biomass, resource fish biomass, and trophic group biomass were all ln(x+1) transformed prior to conducting the LMMs.

Fish trophic group biomass was tested for differences between islands and depths using multivariate analysis of variance (MANOVA). The multivariate test statistic Pillai’s Trace was used because it is robust to heterogeneity of variance and is less likely to involve type I errors than are comparable tests [27]. Canonical discriminate analysis was used to identify and display the nature of the significant differences among islands and depths found by the MANOVA. Trends in the trophic groups were represented as vectors given by correlations of these variables with the canonical variates. These vectors were plotted on the first two canonical axes, together with the treatment centroids and 95% confidence clouds. The strength of each variable in discriminating among groups was displayed graphically as the length of these vectors. Trophic group biomasses were ln(x+1) transformed to conform to the assumption of multivariate normality prior to analysis. We performed univariate GLMs, as described above, if the MANOVA was significant.

Although our main focus was on comparisons between Selvagens and Madeira, the few small MPAs around Madeira provided insights into the potential benefits of protection within this highly developed island. Owing to small sample sizes, comparisons of total fish biomass among MPAs on Madeira with Selvagens and fished areas around Madeira were conducted using a Kruskal-Wallis rank-sum test, with Dunn’s test for unplanned multiple comparisons [28]. All LMM, MANOVA and non-parametric analyses were performed using JMP Pro 12.2 [29].

Overall community structure between islands and depths were compared using Principal Components Analysis on dominant sessile benthic cover, mobile invertebrates, and fish biomass by trophic group with supplemental variables (island, depth strata) projected onto the unconstrained ordination using the ordination program CANOCO version 5.0 [30]. Sessile benthic cover data were arcsine square root transformed, mobile invertebrates and fish biomass were ln(x+1) transformed prior to analysis. Observations in the input matrix were unique sampling stations (n = 122). Functional groups were centered (subtraction of the average of the column values) and standardized (division by the standard deviation of the column values), resulting in each column having zero mean and unit variance.

Results

Benthic community

There was a significant difference in the assemblages of sessile benthic functional groups between islands (Pseudo-F1,121 = 35.6, p < 0.001) and between depths (Pseudo-F1,121 = 3.6, p = 0.02), but not in their interaction (Pseudo-F1,121 = 0.1, p = 0.91). Turf algae cover (35%) and erect algae (30%) accounted for the majority of the sessile benthic cover around Selvagens, while sessile cover at Madeira was dominated by crustose coralline algae (CCA, 30%), turf algae (21%), sessile invertebrates (16%, primarily Balanus sp.), erect algae (12%), and barrens (11%) (Table 1). Average dissimilarity in sessile benthic functional group cover between islands was 58.6%, with turf algae contributing the most to this dissimilarity (22%), followed by CCA (22%), and erect algae (20%) (Table 1). Average dissimilarity of sessile benthic cover between depths was 47.8% (Table 1). Turf algal cover was higher at 20 m and contributed 24% to this dissimilarity, followed by CCA (18%), which was also higher at 20 m. Erect algae was higher at 10 m and contributed an additional 15% to the dissimilarity between islands.

thumbnail
Table 1.

A. Similarity of Percentages (SIMPER) for sessile benthic cover by functional group most responsible for the percent dissimilarities between islands using Bray-Curtis similarity analysis of hierarchical agglomerative group average clustering. Average dissimilarity = 58.6%, with one standard deviation of the mean in parentheses. B. SIMPER for sessile benthic cover by functional group most responsible for the percent dissimilarities between depths. Average dissimilarity = 47.8%. Abundance values are means and one standard deviation of the mean in parentheses.

https://doi.org/10.1371/journal.pone.0187935.t001

Analysis of sessile benthic functional assemblage structure showed clear separation in ordination space, with sites at Selvagens showing higher concordance relative to Madeira (Fig 2). PCO1 explained 71% of the variation in sessile benthic functional cover between islands, with CCA, barrens, and invertebrates trending in the direction of Madeira and erect algae and turf in the direction of Selvagens. PCO2 only explained 12% of the variation in sessile benthic functional cover, with encrusting algae trending towards shallower sites.

thumbnail
Fig 2. Principal Coordinate Analysis of sessile benthic functional groups.

Percent cover data were arcsine square root transformed prior to analysis.

https://doi.org/10.1371/journal.pone.0187935.g002

Mobile invertebrates

There was a significant difference in the assemblages of mobile benthic functional groups between islands (Pseudo-F1,121 = 22.5, p < 0.001) and between depths (Pseudo-F1,121 = 14.2, p < 0.001) but not in their interaction (Pseudo-F1,121 = 0.5, p = 0.79). The sea urchin Diadema africanum was the most common mobile invertebrate observed at both islands (Table 2). The average density at Madeira (4.4 individuals m-2, ± 4.4 sd) was however 65% higher when compared with Selvagens (2.7 individuals m-2, ± 2.7 sd), and accounted for 45.5% of the dissimilarity between island groups. The hermit crab, Calcinus tubularis, was the second most abundant mobile invertebrate at Madeira, while at Selvagens the purple sea urchin, Paracentrotus lividus was most abundant. Densities of D. africanum were 35% higher at the 20 m sites compared with the 10 m sites and accounted for 45.0% of the dissimilarity between depths (Table 2). In contrast, densities of the black (Arbacia lixula) and purple (P. lividus) sea urchins were > 4 times higher at 10 m compared with 20 m, although these species were much more abundant at Selvagens compared with Madeira.

thumbnail
Table 2.

A. SIMPER for mobile invertebrates most responsible for the percent dissimilarities between islands using Bray-Curtis similarity analysis of hierarchical agglomerative group average clustering. Average dissimilarity = 71.5%, with one standard deviation of the mean in parentheses. B. SIMPER for mobile invertebrates most responsible for the percent dissimilarities between depths. Average dissimilarity = 66.8%. Abundance values are means and one standard deviation of the mean in parentheses.

https://doi.org/10.1371/journal.pone.0187935.t002

Analysis of mobile invertebrate densities showed clear separation in ordination space between islands, with Selvagens showing higher concordance among sites (Fig 3). PCO1 explained 32.7% of the variation in mobile invertebrate assemblage structure and was correlated with the combination of island and depth. PCO2 explained an additional 24.5% of the variation. D. africanum, along with Stramonita haemastoma and Percnon gibbesi drove much of the separation along PCO1, while P. lividus and C. tubularis were orthogonal to these species, with the former explaining variation at Selvagens and the latter at Madeira.

thumbnail
Fig 3. Principal Coordinate Analysis of mobile invertebrate communities between islands and depths.

Top 10 species together accounted for > 90% of the individuals in each island group. Data square root transformed prior to analysis.

https://doi.org/10.1371/journal.pone.0187935.g003

Fishes

Fish assemblage characteristics.

Overall species richness per transect was similar between islands, but there was a significant interaction between depth and island, with the 20 m sites at Madeira having significantly higher species richness compared with the other three depth x island combinations (Table 3, Fig 4). Density (number of individuals m-2) was significantly greater at Selvagens compared with Madeira, and significantly higher in the 10 vs. 20 m depth stratum. Total fish biomass (g m-2) was 3.2 times higher at Selvagens compared to Madeira, but there was a significant interaction between island and depth, with deep sites at Madeira not significantly different from deep sites at Selvagens despite an 80% higher biomass at the latter. Biomass of resource species followed a similar pattern. Diversity was higher at Madeira compared with Selvagens and higher in 20 m compared with 10 m.

thumbnail
Fig 4. Comparisons of fish assemblage characteristics between islands and depths.

Box plots showing median (black line), mean (red dashed line), upper and lower quartiles, and 5th and 95th percentiles. A. Species richness transect, B. Number of individuals m-2, C. Grams m-2, and D. Shannon-Weiner Diversity.

https://doi.org/10.1371/journal.pone.0187935.g004

thumbnail
Table 3. Comparisons of fish assemblage characteristics between islands and depth strata.

Results of Linear Mixed Models with station as a random effect, and Tukey-Kramer HSD (honestly significant difference) pairwise comparisons among fixed factors. Underlined island x depth combinations are not significantly different (α = 0.05). Mad.–Madeira, Sel.–Selvagens.

https://doi.org/10.1371/journal.pone.0187935.t003

Fish assemblage structure.

Fish biomass assemblage structure was significantly different between Madeira and Selvagens (pseudo-F1,121 = 24.5, p < 0.001), between depths (pseudo-F1,121 = 4.4, p < 0.001), as well as their interaction (pseudo-F1,121 = 4.1, p < 0.001). Fish assemblage structure based on numerical abundance showed similar patterns with significant differences between Madeira and Selvagens (pseudo-F1,121 = 21.7, p < 0.001), between depths (pseudo-F1,121 = 5.9, p < 0.001), as well as their interaction (pseudo-F1,121 = 4.0, p = 0.004). Fish assemblage structure based on species biomass showed clear separation in ordination space between Selvagens and Madeira based on Principal Coordinates Analysis (Fig 5). The 1st PCO axis explained ~23% of the variation in assemblage structure. This separation was driven by Diplodus vulgaris towards Madeira and Kyphosus sectatrix, Serranus atricauda, Thalassoma pavo, Sparisoma cretense, and Bodianus scrofa towards Selvagens. Boops boops was orthogonal to this primary axis. Average dissimilarity between islands was 79.2% and was driven by K. sectatrix (22.6%), S. cretense (8.6%), B. boops (8.3%), and T. pavo (7.8%) (Table 4). Average dissimilarity between depth strata was 72.3% and was driven by K. sectatrix (17.7%), B. boops (10.8%), Chromis limbata (9.3%), and S. cretense (7.4%) (Table 4).

thumbnail
Fig 5. Principal Coordinates Analysis of fish species biomass.

Square root transformation. Bray Curtis Similarity matrix.

https://doi.org/10.1371/journal.pone.0187935.g005

thumbnail
Table 4.

A. SIMPER for fish species by biomass (g m-2) most responsible for the percent dissimilarities between islands using Bray-Curtis similarity analysis of hierarchical agglomerative group average clustering. Average dissimilarity = 79.2%, with one standard deviation of the mean in parentheses. B. SIMPER for fish species by biomass most responsible for the percent dissimilarities between depths. Average dissimilarity = 72.3%. Abundance values are means and one standard deviation of the mean in parentheses.

https://doi.org/10.1371/journal.pone.0187935.t004

There were large differences in the biomass of important resource fish species between Selvagens and Madeira (Table 4). Chubs (K. sectatrix) accounted for 34% of the biomass at Selvagens, with absolute biomass 182% higher than at Madeira. The amberjack Seriola dumerili was one of the most important species by weight at Selvagens (16% of total biomass), but was not observed on transects around Madeira. Another amberjack, S. rivoliana, was 7 times more abundant at Selvagens compared to Madeira. Barred hogfish (B. scrofa), a sea urchin predator, was 3.4 times more abundant by weight at Selvagens compared to Madeira. Another predator of small sea urchins, S. cretense [20], was 6.8 times more abundant by weight at Selvagens compared with Madeira. A total of fifteen Dusky Grouper (Epinephelus marginatus), which is listed as endangered by IUCN, were observed around Selvagens, ranging in size from 40 to 120 cm. Only four dusky groupers were observed around Madeira and these were all exclusively found within the Garajau MPA.

Fish trophic structure.

Trophic structure based on biomass was significantly different between island groups (F1,118 = 26.85, p < 0.001) but not depth (F1,118 = 0.10, p = 0.756); however, the interaction between island and depth was significant (F1,118 = 5.56, p = 0.020). The 10-m stratum at Selvagens was distinct from the other three island x depth combinations, with herbivores most responsible for this separation (Fig 6). The Selvagens 20 m stratum was distinct from the 20 m Madeira stratum but overlapped with the Madeira 10 m stratum. Top predators accounted for this separation and was orthogonal to herbivores.

thumbnail
Fig 6. Canonical discriminate analysis of the difference in trophic group biomass between island groups and depths.

Treatment (island group x depth) centroids and 95% confidence clouds are plotted together with the direction and importance of trends in trophic group.

https://doi.org/10.1371/journal.pone.0187935.g006

Biomass of top predators was an order of magnitude higher at Selvagens compared to Madeira, and accounted for 22% of total biomass at Selvagens but < 4% at Madeira (Table 5, Fig 7). Herbivore biomass was seven times higher at Selvagens compared with Madeira, but there was a significant interaction with depth. Herbivores comprised 45% of the total biomass at Selvagens and < 20% at Madeira. Secondary consumers had significantly higher biomass at Selvagens, as well as in the deep depth stratum. Planktivores were not significantly different between the two island groups or between depth strata.

thumbnail
Table 5. Comparisons of fish trophic biomass between islands and depth strata.

Results of Linear Mixed Models with station as a random effect, and Tukey-Kramer HSD (honestly significant difference) pairwise comparisons among fixed factors. Underlined island x depth combinations are not significantly different (α = 0.05). Mad.–Madeira, Sel.–Selvagens.

https://doi.org/10.1371/journal.pone.0187935.t005

thumbnail
Fig 7. Comparisons of trophic biomass between islands and depths.

Box plots showing median (black line), mean (red line), upper and lower quartiles, and 5th and 95th percentiles. A. Top predators, B. herbivores, C. secondary consumers, and D. planktivores. Values are in g m-2.

https://doi.org/10.1371/journal.pone.0187935.g007

Marine protected areas (MPAs).

Total fish biomass was significantly different among management regimes (X2 = 46.3, p < 0.001). Biomass in the Garajau MPA was 2.4 times higher than areas open to fishing around Madeira (p < 0.001) and 2.2 times higher than the partially protected Rocha do Navio MPA (Fig 8). The latter comparison, while not significant, is suggestive of differences between the two MPAs (p = 0.07). Total biomass at Garajau was only 1/3rd lower than Selvagens, while biomass at partially protected Rocha do Navio was 2.8 times lower than Selvagens. Top predators comprised 15% of the total biomass at Garajau but only < 4% at both Rocha do Navio and areas open to fishing around Madeira. Top predator biomass at Garajau was 90% lower than Selvagens, compared to Rocha do Navio, where top predator biomass was more than 15 times lower than Selvagens.

thumbnail
Fig 8. Comparisons of fish trophic biomass (g m-2) by management regime.

Values are means and standard error of the mean. Kruskal-Wallis Rank Sum comparisons of total fish biomass among management regimes was statistically different (X2 = 46.3, p < 0.001). Management regimes with the same letter are not significantly different from one another (Dunn’s unplanned multiple comparisons procedures, α = 0.05).

https://doi.org/10.1371/journal.pone.0187935.g008

Whole community comparisons between islands.

There was clear separation between Madeira and Selvagens based on dominant sessile benthic cover, mobile invertebrates, and fish biomass by trophic group (Fig 9). Depth strata within islands were close to one another in ordination space, with the two depth strata at Selvagens having higher concordance compared to Madeira. The first two axes of the PCA biplot explained 46.5% of the variance in community structure and 63% of the community-island x depth relationships (Table 6). The major drivers of this separation for Selvagens were fish biomass of top predators and canopy algae. The separation for Madeira was driven by the hermit crab C. tubularis, secondary consumers and CCA. Erect algae, turf, and herbivore biomass were orthogonal to the primary vectors of island separation towards Selvagens, while barrens, invertebrates, and D. africanum were orthogonal in the direction of Madeira.

thumbnail
Table 6. A. Results of Principal Components Analysis with supplemental variables on dominant sessile benthic cover, mobile invertebrates, and fish biomass by trophic group.

https://doi.org/10.1371/journal.pone.0187935.t006

thumbnail
Fig 9. Biplot of results of Principal Components Analysis with supplemental variables on dominant sessile benthic cover, mobile invertebrates, and fish biomass by trophic group.

Sessile benthic cover data were arcsine square root transformed, mobile invertebrates and fish biomass were ln(x+1) transformed prior to analysis. All data were centered and standardized. Statistical results are shown in Table 6.

https://doi.org/10.1371/journal.pone.0187935.g009

Discussion

The Selvagens and Madeira islands present a sharp contrast between a highly urbanized and developed island (Madeira) and a remote and virtually uninhabited marine reserve. The effects of fishing and land-based activities around Madeira have been well documented [3134], and while we did not directly evaluate the anthropogenic stressors between Madeira and Selvagens, the extensive development and human activity at the former and the remoteness and near absence of human habitation at the latter clearly highlights the disparities in human pressure experienced between these two locations. As a result of these differences, the Selvagens Islands may represent one of the last remaining intact marine ecosystems in the eastern Atlantic.

Remote locations with limited fishing pressure and few land-based stressors are some of the few remaining examples of marine ecosystems without major anthropogenic influences [3536]. Although the North Atlantic has been overexploited for centuries [3739], the establishment of the Selvagens Islands Reserve in 1971 has resulted in the maintenance of a healthy ecosystem in a region where the surrounding seas are intensely fished and largely degraded. While illegal fishing has been reported around Selvagens, the presence of park rangers and the recent addition of a Portuguese maritime police force seems to have limited the impact of these activities as shown by the high abundance of large fisheries species, including the dusky groupers—one of the main targets of recreational and commercial fisheries in the nearby islands and in the Mediterranean.

The intertidal community around Selvagens was noteworthy for its abundance of large sun limpets (Patella candei), particularly on Selvagem Pequena (Fig 10). This species has become rare throughout much of its range due to overfishing and is listed as in danger of extinction in the Canary Islands [40]. In addition, top-shell snails (Phorcus atratus) and other limpets, mainly Patella aspera and Siphonaria pectinata, were also common. These extremely high densities of intertidal grazers appear to limit macroalgal growth, which was restricted to small patches of turf algae, mainly Jania cf. rubens in the lower intertidal. The intertidal around Selvagens likely represents one of the few remaining intact ecosystems of its kind in the region.

thumbnail
Fig 10. The intertidal habitats of Ilhas Selvagens are some of the least disturbed within the region, with the highly prized sun limpet (Patella candei) common and of large size, particularly around Selvagem Pequena.

https://doi.org/10.1371/journal.pone.0187935.g010

Turf algae and erect algae accounted for the majority of the sessile benthic cover around Selvagens, while Madeira was dominated by CCA, turf algae, sessile invertebrates, and sea urchin barrens. Our results from Madeira are consistent with previous studies of the benthos around this island [3133]. Erect macroalgae provides the main biological substrate for many organisms in the eastern Atlantic [17, 4142]. In the nearby Canary Islands, upright seaweeds are the principal engineering organisms on shallow rocky bottoms, providing complex habitats that support highly diverse communities [17, 43].

The sea urchin D. africanum was the most common mobile invertebrate observed at both islands, but densities were 65% higher at Madeira compared with Selvagens. The densities of D. africanum that we observed around Madeira are similar to those documented by researchers back in the 1990s and early 2000s [3133]. In the eastern Atlantic, removal of top predators because of overfishing has been linked to hyperabundances of sea urchins (predominantly Diadema), with the subsequent elimination of erect vegetative frameworks and the creation of barrens as an alternate stable state [17, 20, 4347]. These upright seaweed beds exist due to the balance between seaweeds, herbivores, and predators [19, 48], but human-derived stressors (e.g., overfishing, pollution) can reduce the resilience of these desirable macroalgal beds, while strengthening the influence of urchin barrens, and thus exacerbating the threat, spatial extent, and irreversibility of an unwanted regime shift [1718, 21]. Previous studies around Madeira have shown similar negative relationships between sea urchin densities and macroalgal cover [3233].

The identity of predatory fishes on the sea urchin D. africanum is well known from previous studies in the Canary Islands by directly assessing the frequency of predation events and indirectly through examination of fish stomach contents [20]. These authors also showed that the depletion of sea urchin predators in fished areas compared with MPAs resulted in increased sea urchin populations and cascading effects that reduces benthic diversity in areas open to fishing. Species known to feed on and control sea urchins (e.g., Balistes capriscus, B. scrofa, Canthidermis sufflamen, and S. cretense, [20]) were in high abundance at Selvagens compared to Madeira, where they are also important fisheries species.

Total fish biomass was more than three times higher at Selvagens compared to Madeira, and biomass of top predators was an order of magnitude greater. Several commercial species (e.g., groupers, jacks, triggerfishes) that have been overfished elsewhere in Macaronesia [1618, 49], but were common and of large size at Selvagens compared with Madeira. Previous studies of the fish assemblages around Madeira dating back to the mid-1990s documented few large predatory fishes or other species with high resource value [24, 5051], suggesting that these species may have been subjected to overfishing for many years.

The IUCN endangered Dusky Grouper (E. marginatus) is the best-known grouper of the Mediterranean Sea and adjacent European and North African coasts, but it has been overexploited throughout much of its range [5253]. The prevalence and large sizes of this species at Selvagens is striking compared to Madeira, where it was only found within the highly restricted Garajau MPA. A similar pattern has been observed in other Macaronesian regions such as in the Canary Islands, where groupers are much more abundant around islands where fishing activities are restricted by MPAs and human population is low [16]. The Garajau MPA was established in 1986 and was the first exclusively marine protected area in Portugal. Overall biomass and that of top predators and secondary consumers was higher within this MPA when compared with areas open to fishing around Madeira and was comparable to Selvagens; however, its small size (3.8 km2) likely results in limited benefits beyond its boundaries.

The effects of fishing and other anthropogenic influences associated with urbanization are evident around Madeira [15, 34, 5455]. This is in stark contrast to Selvagens, which harbors rich benthic communities with diverse algal assemblages and high fish biomass, along with an abundance of large resource species. The clear differences between these two islands highlights the importance of protecting the Selvagens, which harbors one of the last intact marine ecosystems in the North Atlantic, and the need to increase management and protection around Madeira if the ecosystem is to recover and provide the ecosystem services essential to the island community. No-take areas have been identified as an effective tool to restore erect macroalgal beds in other Macaronesian islands [56], and are known to provide overall ecosystem benefits both within and beyond their borders [57].

Benthic and fish communities at Selvagens resemble those of some MPAs around other areas of Macaronesia. Densities of the most common mobile invertebrate D. africanum at Selvagens are within the range recorded in the no-take reserve of La Graciosa in the Canary Islands [45]. However, sea urchin abundances are still an order of magnitude higher than densities recorded at the most effective MPA in the Canary Islands, Mar de Las Calmas, on El Hierro Island [45]. Similarly, fish biomass at Selvagens were similar to those in the La Graciosa no-take reserve, but lower than those in the Mar de Las Calmas [58]. These differences may be driven by inherent characteristics of both Selvagens and the Canary Islands, with the former showing a higher proportion of warm-temperate species and the latest having more complex and variable assemblages, likely due to the larger sizes of the islands and more heterogeneous habitats [59].

Threats to the Selvagens include illegal fishing within the reserve and unregulated or weakly monitored fishing for tuna and other target species surrounding the reserve. The current 200 m depth limit for the reserve means that fishing can occur very close to the islands, with potential impacts to nearshore species. Expansion of the reserve would provide protection for coastal and pelagic species, as well as reducing by-catch of sea birds, marine mammals, and large pelagic fishes (e.g., tuna and billfishes) that frequent the area. In addition, the Selvagens may be an important biogeographical link between Madeira and the Canary Islands [60].

The results of this study reinforce patterns observed at smaller spatial scales within the Canary Islands [1617, 4546], and are consistent with observations contrasting remote vs. inhabited marine ecosystems around the world [6163]. The Selvagens harbor one of the last intact marine ecosystems in the North Atlantic, and maintains high coastal fish species diversity within a relatively small area [60]. The Selvagens Nature Reserve serves as a global model for what can be achieved elsewhere if species are protected and allowed to recover within their borders. Increased protection for this unique area is a precautionary bulwark against the degradation and decline of marine ecosystems throughout the region. Our results also identify the need for better fisheries and coastal zone management, as well as the need for larger and more effective marine protected areas around Madeira if the ecosystem is to recover and provide the ecosystem services essential to the island community.

Acknowledgments

We would like to thank the captain and crew of the MY Plan B, as well as Pedro Vasconcelos and Hugo Silva of Focus Natura Dive Center for field support during the expeditions. Fundação para a Ciência e a Tecnologia (FCT) supported this work under the strategic project UID/MAR/04292/2013. We are grateful to the Serviço do Parque Natural da Madeira, especially Paulo Oliveira and Carolina Jardim Santos, for their assistance with permits, logistics, and local knowledge. The Hydrographic Institute of Portugal graciously provided bathymetric information for Selvagens.

References

  1. 1. Floeter SR, Rocha LA, Robertson DR, Joyeux JC, Smith-Vaniz WF, Wirtz P, et al. Atlantic reef fish biogeography and evolution. J Biogeogr. 2008;35: 22–47.
  2. 2. Briggs JC, Bowen BW. A realignment of marine biogeographic provinces with particular reference to fish distributions. J Biogeogr. 2012;39: 12–30.
  3. 3. Almada VC, Toledo JF, Brito A, Levy A, Floeter SR, Robalo JI, et al. Complex origins of the Lusitania biogeographic province and northeastern Atlantic fishes. Front Biogeogr. 2013;5: 20–28.
  4. 4. Fernández-Palacios JM, de Nascimento L, Otto R, Delgado JD, García-del-Rey E, Arévalo JR, et al. A reconstruction of Palaeo-Macaronesia, with particular reference to the long-term biogeography of the Atlantic island laurel forests. J Biogeogr. 2011;38: 226–246.
  5. 5. Kondraskov P, Schütz N, Schüßler C, de Sequeira MM, Guerra AS, Caujapé-Castells J, et al. Biogeography of Mediterranean Hotspot Biodiversity: Re-Evaluating the 'Tertiary Relict' Hypothesis of Macaronesian Laurel Forests. PloS ONE. 2015;10: e0132091. pmid:26173113
  6. 6. Wirtz P. Twelve invertebrate and eight fish species new to the marine fauna of Madeira, and a discussion of the zoogeography of the area. Helgolander Meeresuntersuchungen. 1998;52: 197–208.
  7. 7. Wirtz P, Fricke R, Biscoito MJ. The coastal fishes of Madeira Island–new records and an annotated check-list. Zootaxa. 2008;1715: 1–26.
  8. 8. Nguyen HT, Fernandez-Steeger TM, Wiatr T, Rodrigues D, Azzam R. Use of terrestrial laser scanning for engineering geological applications on volcanic rock slopes-an example from Madeira Island (Portugal). Nat Hazard Earth Sys 2011;11: 807–817.
  9. 9. Brandmüller T, Önnerfors A. Eurostat Regional Yearbook 2016 edition. Luxembourg: Publications Office of the European Union. http://ec.europa.eu/eurostat/documents/3217494/7604195/KS-HA-16-001-EN-N.pdf/76c007e9-6c1d-435a-97f8-e5ea700aa149
  10. 10. Menezes DM, Cunha R, Freitas CM, Mendes JP, Oliveira M, Parrot C. As Ilhas Selvagens. Governo Regional da Madeira, Secretaria Regional do Ambiente e dos Recursos Naturais, Servico do Parque Natural da Madeira. 2004.
  11. 11. Grimmett RF, Jones T. Important Bird Areas in Europe. Cambridge, UK: International Council for Bird Preservation Technical Publication No. 9; 1989.
  12. 12. Granadeiro JP, Dias MP, Rebelo R, Santos CD, Catry P. Numbers and populations trends of Cory’s shearwater Calonectris diomedea at Salvagem Grande, Northeast Atlantic. Waterbirds. 2006;29: 56–60.
  13. 13. Olivera P, Menezes D, Trout R, Buckle A, Geraldes P, Jesus J. Successful eradication of the European rabbit (Oryctolagus cuniculus) and house mouse (Mus musculus) from the island of Selvagem Grande (Macaronesian archipelago), in the Eastern Atlantic. Integr Zool. 2010;5: 70–83. pmid:21392324
  14. 14. Santos F. A pesca na História da Madeira: Estado dos conhecimentos e problemáticas futuras de análise. Anuário 2010, Centro de Estudos de História do Atlântico, 354–397.
  15. 15. Hermida M, Delgado J. High trophic level and low diversity: Would Madeira benefit from fishing down?. Mar Policy. 2016;73: 130–137.
  16. 16. Tuya F, Sanchez-Jerez P, Haroun RJ. Populations of inshore serranids across the Canarian Archipelago: relationships with human pressure and implications for conservation. Biol Conserv. 2006;128: 13–24
  17. 17. Hernández JC, Clemente S, Sangil C, Brito A. The key role of the sea urchin Diadema aff. antillarum in controlling macroalgae assemblages throughout the Canary Islands (eastern subtropical Atlantic): an spatio-temporal approach. Mar Environ Res. 2008;66: 259–70. pmid:18479745
  18. 18. Sangil C, Martín-García L, Hernández JC, Concepción L, Fernández R, Clemente S. Impacts of fishing and environmental factors driving changes on littoral fish assemblages in a subtropical oceanic island. Estuar Coast Shelf Sci. 2013;128: 22–32.
  19. 19. Sala E, Boudouresque CF, Harmelin-Vivien M. Fishing, trophic cascades, and the structure of algal assemblages: evaluation of an old but untested paradigm. Oikos. 1998;82: 425–439.
  20. 20. Clemente S, Hernández JC, Rodríguez A, Brito A. Identifying keystone predators and the importance of preserving functional diversity in sublittoral rocky-bottom areas. Mar Ecol Prog Ser. 2010;413: 55–67.
  21. 21. Ling SD, Scheibling RE, Rassweiler A, Johnson CR, Shears N, Connell SD, et al. Global regime shift dynamics of catastrophic sea urchin overgrazing. Phil Trans R Soc B. 2015;370: 20130269.18.
  22. 22. Martins AM, Amorim AS, Figueiredo MP, Sousa RJ, Mendonça A, Bashmachnikov IL, et al. Sea Surface Temperature (AVHRR, MODIS) and Ocean Colour (MODIS) seasonal and interannual variability in the Macaronesian islands of Azores, Madeira, and Canaries. Proc SPIE. 2007;6743: 67430A1–67430A15.
  23. 23. Tuya F, Ortega-Borges L, Sanchez-Jerez P, Haroun RJ. Effect of fishing pressure on the spatio-temporal variability of the parrotfish, Sparisoma cretense (Pisces: Scaridae), across the Canarian Archipelago (eastern Atlantic). Fish Res. 2006;77: 24–33.
  24. 24. Ribeiro C, Almeida AJ, Araújo R, Biscoito M, Freitas M. Fish assemblages of Cais do Carvao Bay (Madeira Island) determined by the visual census technique. J Fish Biol. 2005;67: 1568–84.
  25. 25. Froese R, Pauly D, 2012. FishBase. [cited 2017, Mar 1]. www.Fishbase.org.
  26. 26. Anderson MJ, Gorley RN, Clarke KR. PERMANOVA+ for PRIMER: Guide to Software and Statistical Methods. Plymouth, UK. PRIMER-E, 2008.
  27. 27. Green RH. Sampling Design and Statistical Methods for Environmental Biologists. New York: Wiley Interscience; 1979.
  28. 28. Zar JH (1999) Biostatistical analysis, 4th edition. India: Pearson Education. 662 p.
  29. 29. JMP Pro 12.2, SAS Institute Inc., Cary, NC, 1989–2007.
  30. 30. ter Braak CJF, Šmilauer P. Canoco Reference Manual and User’s Guide: Software for Ordination, Version 5.0. Ithaca, NY, USA. Microcomputer Power: 2012.
  31. 31. Bianchi CN, Morri C, Sartoni G, Wirtz P. Sublittoral epibenthic communities around Funchal (Ilha da Madeira, NE Atlantic). Bol Mus Mun Funchal. 1998;5: 59–80.
  32. 32. Alves F, Chícharo L, Serrao E, Abreu AD. Algal cover and sea urchin spatial distribution at Madeira Island (NE Atlantic). Sci Mar. 2001;65: 383–92.
  33. 33. Alves F, Chícharo L, Serrao E, Abreu AD. Grazing by Diadema antillarum (Philippi) upon algal communities on rocky substrates. Sci Mar. 2007;67: 307–311.
  34. 34. Baioni D. Human activity and damaging landslides and floods on Madeira Island. Nat Hazard Earth Sys. 2011;11: 3035–3046.
  35. 35. D'agata S, Mouillot D, Wantiez L, Friedlander AM, Kulbicki M, Vigliola L. Marine reserves lag behind wilderness in the conservation of key functional roles. Nat Commun. 2016;7: 12000. pmid:27354026
  36. 36. Knowlton N, Jackson JB. Shifting baselines, local impacts, and global change on coral reefs. PLoS Biol. 2008;6: e54. pmid:18303956
  37. 37. Wolff WJ. The south-eastern North Sea: losses of vertebrate fauna during the past 2000 years. Biol Conserv. 2000;95: 209–217.
  38. 38. Jackson JB, Kirby MX, Berger WH, Bjorndal KA, Botsford LW, Bourque BJ, et al. Historical overfishing and the recent collapse of coastal ecosystems. Science. 2001;293: 629–638. pmid:11474098
  39. 39. Barrett JH, Locker AM, Roberts CM. The origins of intensive marine fishing in medieval Europe: the English evidence. Proc R Soc Lond B Biol Sci. 2004;271: 2417–21.
  40. 40. González-Lorenzo G, Hernández EM, Pérez-Dionis G, Hernández AB, Santos BG, Diez JB. Ineffective conservation threatens Patella candei, an endangered limpet endemic to the Macaronesian islands. Biol Conserv. 2015;192: 428–435.
  41. 41. Clemente S, Hernández JC. Influence of wave exposure and habitat complexity in determining spatial variation of the sea urchin Diadema aff. antillarum (Echinoidea: Diadematidae) populations and macroalgal cover (Canary Islands-Eastern Atlantic Ocean). Rev Biol Trop. 2008;56: 229–254.
  42. 42. Sangil C, Sansón M, Clemente S, Afonso-Carrillo J, Hernández JC. Contrasting the species abundance, species density and diversity of seaweed assemblages in alternative states: Urchin density as a driver of biotic homogenization. J Sea Res. 2014;85: 92–103.
  43. 43. Sangil C, Sansón M, Afonso-Carrillo J. Spatial variation patterns of subtidal seaweed assemblages along a subtropical oceanic archipelago: thermal gradient vs herbivore pressure. Estuar Coast Shelf Sci. 2011; 94: 322–33.
  44. 44. Tuya F, Haroun RJ, Boyra A, Sanchez-Jerez P. Sea urchin Diadema antillarum: different functions in the structure and dynamics of reefs on both sides of the Atlantic. Mar Ecol Prog Ser. 2005;302: 307–10.
  45. 45. Hernández JC, Clemente S, Sangil C, Brito A. Actual status of the sea urchin Diadema aff. antillarum populations and macroalgal cover in marine protected areas compared to a highly fished area (Canary Islands—eastern Atlantic Ocean). Aquat Conserv. 2008;18: 1091–1108.
  46. 46. Clemente S, Hernández JC, Brito A. Context-dependent effects of marine protected areas on predatory interactions. Mar Ecol Prog Ser. 2011;437: 119–133.
  47. 47. Tuya F, Boyra A, Sanchez-Jerez P, Barberá C, Haroun RJ. Relationships between rocky-reef fish assemblages, the sea urchin Diadema antillarum and macroalgae throughout the Canarian Archipelago. Mar Ecol Prog Ser. 2004;278: 157–169.
  48. 48. Pinnegar JK, Polunin NV, Francour P, Badalamenti F, Chemello R, Harmelin-Vivien ML, et al. Trophic cascades in benthic marine ecosystems: lessons for fisheries and protected-area management. Environ Conserv. 2000;27: 179–200.
  49. 49. Falcón JM, Bortone SA, Brito A, Bundrick CM. Structure and relationships within and between the littoral rock-substrate fish communities off four islands in the Canarian Archipelago. Mar Biol. 1996; 125: 215–231.
  50. 50. Andrade CA, Albuquerque FM. Fish assemblages associated with bottom habitats on the south coast of Madeira. Boletim do Museu Municipal do Funchal Supplement 4. 1995; 9–20.
  51. 51. Delgado C. Caracterização faunística, batimétrica e geomorfológica da reserva natural do Garajau Parque Natural da Madeira. Relatório de estágio de licenciatura, Universidade do Porto. 1998.
  52. 52. Cornish A, Harmelin-Vivien M. (Grouper & Wrasse Specialist Group). 2004. Epinephelus marginatus. The IUCN Red List of Threatened Species 2004: e.T7859A12857009. http://dx.doi.org/10.2305/IUCN.UK.2004.RLTS.T7859A12857009.en. Downloaded on 02 June 2017.
  53. 53. Schunter C, Carreras-Carbonell J, Planes S, Sala E, Ballesteros E, Zabala M, et al. Genetic connectivity patterns in an endangered species: The dusky grouper (Epinephelus marginatus). J Exp Mar Biol Ecol. 2011;401: 126–133.
  54. 54. Baeta F, Costa MJ, Cabral H. Changes in the trophic level of Portuguese landings and fish market price variation in the last decades. Fish Res. 2009;97: 216–222.
  55. 55. Pereira E, Figueira C, Aguiar N, Vasconcelos R, Vasconcelos S, Calado G, et al. Microbiological and mycological beach sand quality in a volcanic environment: Madeira archipelago, Portugal. Sci. Total Environ. 2013;461: 469–79. pmid:23747562
  56. 56. Sangil C, Clemente S, Martín-García L, Hernández JC. No-take areas as an effective tool to restore urchin barrens on subtropical rocky reefs. Estuar Coast Shelf Sci. 2012;112: 207–215.
  57. 57. Gaines SD, White C, Carr MH, Palumbi SR. Designing marine reserve networks for both conservation and fisheries management. Proc Natl Acad Sci USA. 2010;107: 18286–18293. pmid:20200311
  58. 58. Brito A, Falcón JM, González G, Dorta C, Rodríguez A, Martín-Sosa Rodríguez, et al. Estudio científico para el seguimiento de las Reservas Marinas del entorno de la Isla de La Graciosa e Islotes del Norte de Lanzarote y del entorno de Punta Restinga-Mar de las Calmas. Viceconsejería de Pesca del Gobierno de Canarias. 2008. 205 pp.
  59. 59. Falcón JM. Ictiofauna de las Islas Canarias: Análisis biogeográfico [dissertation]. San Cristóbal de La Laguna, Tenerife (Canary Islands): Universidad de La Laguna; 2015.
  60. 60. Almada F, Abecasis D, Villegas-Rios D, Henriques S, Pais MP, Batista M, et al. Ichthyofauna of the Selvagens Islands. Do small coastal areas show high species richness in the northeastern Atlantic?. Mar Biol Res. 2015;11: 49–61.
  61. 61. Friedlander AM, DeMartini EE. Contrasts in density, size, and biomass of reef fishes between the northwestern and the main Hawaiian islands: the effects of fishing down apex predators. Mar Ecol Prog Ser. 2002;230: 253–264.
  62. 62. Friedlander AM, Ballesteros E, Beets J, Berkenpas E, Gaymer CF, Gorny M, et al. Effects of isolation and fishing on the marine ecosystems of Easter Island and Salas y Gómez, Chile. Aquat Conserv. 2013;23: 515–531.
  63. 63. Friedlander AM, Obura D, Aumeeruddy R, Ballesteros E, Church J, Cebrian E, et al. Coexistence of low coral cover and high fish biomass at Farquhar Atoll, Seychelles. PloS ONE. 2014;9: e87359. pmid:24489903