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

Nearshore marine biodiversity of Osa Peninsula, Costa Rica: Where the ocean meets the rainforest

  • 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, Hawaiʿi Institute of Marine Biology, University of Hawaiʿi, Kāneʻohe, Hawaiʿi, United States of America

  • Enric Ballesteros,

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

    Affiliation Centre d’Estudis Avancats de Blanes-CSIC, Blanes, Girona, Spain

  • Odalisca Breedy,

    Roles Data curation, Formal analysis, Investigation, Methodology, Resources, Validation, Writing – review & editing

    Affiliation Centro de Investigación en Ciencias del Mar y Limnología, Universidad de Costa Rica, San José, Costa Rica

  • Beatriz Naranjo-Elizondo,

    Roles Data curation, Formal analysis, Investigation, Methodology, Resources, Validation, Writing – review & editing

    Affiliation Centro de Investigación en Ciencias del Mar y Limnología, Universidad de Costa Rica, San José, Costa Rica

  • Noelia Hernández,

    Roles Investigation, Project administration, Resources, Validation, Writing – review & editing

    Affiliation Asociación Conservación Osa, Puerto Jiménez, Golfito, Costa Rica

  • Pelayo Salinas-de-León,

    Roles Data curation, Investigation, Methodology, Resources, Writing – review & editing

    Affiliation Charles Darwin Research Station, Charles Darwin Foundation, Puerto Ayora, Galápagos Islands, Ecuador

  • Enric Sala,

    Roles Conceptualization, Funding acquisition, Resources, Supervision, Validation, Writing – review & editing

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

  • Jorge Cortés

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

    Affiliation Centro de Investigación en Ciencias del Mar y Limnología, Universidad de Costa Rica, San José, Costa Rica

Abstract

Osa Peninsula in remote southwest Costa Rica harbors 2.5% of global terrestrial biodiversity in only 1,200 km2 and has the largest remaining tract of Pacific lowland wet forest in Mesoamerica. However, little is known about the marine ecosystems of this diverse region. Much of the coastline consists of soft sediment exposed to strong wave action. Three major hard bottom habitat types define this region, including: 1) coral reefs around Isla del Caño Biological Reserve, a no-take marine protected area (MPA) of 52 km2, 2) coastal rocky reefs and islets along the peninsula, including Corcovado National Park, and 3) submerged pinnacles just outside the Isla del Caño MPA. Average coral cover at Isla del Caño was 21%, composed primarily of Porites lobata and Pocillopora elegans. In contrast, coastal rocky reefs were dominated by turf algae (39.8%) and macroalgae (20.7%) with low coral cover (1.1%). Submerged pinnacles were dominated by crustose coralline algae (33.3%) and erect coralline algae (25.7%). Fish assemblage characteristics (species richness, abundance, biomass) were significantly higher at the pinnacles compared to the other habitats and was dominated by schooling species such as Haemulon steindachneri, and the herbivores Kyphosus ocyurus, and Acanthurus xanthopterus. Top predators, primarily Triaenodon obesus, Caranx sexfasciatus, and Lutjanus argentimaculatus, were also most abundant at these pinnacles and accounted for the largest differences in fish trophic structure among habitats. Despite Isla del Caño being fully protected from fishing, biomass was similar to fished areas along the coast and lower than the adjacent submerged pinnacles outside the reserve. Similarly, Corcovado National Park includes 20.3 km2 of no-take MPAs; however, there is limited enforcement, and we noted several instances of fishing within the park. The unique configuration of healthy offshore coral reefs and pinnacles connected to coastal habitats provides corridors for many species including large predators such as sharks and other marine megafauna, which warrants additional protection.

Introduction

Osa Peninsula in the remote southwest corner of Costa Rica is recognized as one of the most biodiverse places on earth relative to its size (~ 1,200 km2), with 2.5% of the world’s known terrestrial species of plants and animals [1]. Once an island in the Pacific, the peninsula evolved in isolation until it merged with mainland Costa Rica two million years ago [2]. One of the last places in Costa Rica to be settled by large human groups, the Osa Peninsula has the largest remaining tract of Pacific lowland wet forest in Mesoamerica [1], which extends all the way to the Pacific Ocean. High species richness has also been reported from Osa’s marine-coastal ecosystems [3]. Within this area are the largest mangrove areas in the country, protected by the Térraba-Sierpe National Wetland [4], and the only tropical fjord with anoxic waters, Golfo Dulce, which separates Osa Peninsula from the Costa Rican mainland [5]. This embayment of tectonic origin, despite being anoxic at depths > 60 m, has a high diversity of marine fauna [68], and is important for marine mammals, sharks, and seabirds [3, 9, 10].

In addition, there are a diversity of marine ecosystems from shallow to deep waters, including rocky outcrops scattered along the coast, and submerged pinnacles and various types of hard and soft bottom habitats in a wide range of depths, making this region one of the most diverse in the country [11]. The Osa Peninsula is an important breeding habitat for highly migratory species such as scalloped hammerhead sharks (Sphyrna lewini) [12] and humpback whales (Megaptera novaeangliae) [3]. However, the biodiversity and marine-coastal environments have been degraded by various human activities, mainly due to the rapid expansion of agricultural and tourism development, land contamination, and overexploitation of natural resources [13].

Covering much of the western portion of the Osa Peninsula is the Corcovado National Park (422.3 km2 including the marine portion). Corcovado is an icon of terrestrial conservation since its creation in 1975, consisting of lowland rainforest, highland cloud forest, jolillo palm forest, mangrove swamps, marine coasts, and beach habitats [14]. The marine boundaries of this park extend ~ 40 km along the coast and out to 500 m (covering 20.4 km2), which were added in 1980 mainly to protect sea turtles nesting sites although fishing is also prohibited [15]. Several sections along the shallow coast have well developed coral reefs; however, much of nearshore ecosystem has been poorly explored [16, 17].

The Isla del Caño Biological Reserve, located ~ 15 km off the coast, is a no-take marine protected area (MPA, 52 km2), which is home to the most extensive coral reefs on the Pacific coast of Costa Rica [1821]. The corals of Isla del Caño were severely impacted by the 1982–1983 El Niño event, resulting in losses of up to 50% of live coral cover [18, 19, 22]; however, subsequent El Niño events in the 1990s impacted these reefs to a lesser extent [23]. Isla del Caño has a high richness of octocorals with at least 12 species known from 10 to 30 m. Between 30 and 50 m, octocoral species richness declines; however, there is a unique faunal assemblage in these deeper depths and a new species to science was recently described that suggests a higher number of octocoral species are still to be discovered [24].

According to the marine conservation gap analysis conducted by the Costa Rican National System of Conservation Areas (SINAC), the Osa region is one of the main areas of the country where significant gaps were found [11]. To date, these conservation gaps have not been effectively addressed, and marine conservation efforts occur in only two small MPAs. Because of the high terrestrial biodiversity of Osa Peninsula, it is important to also understand the adjacent marine environment and its importance in land-sea linkages. The marine corridor adjacent to Osa Peninsula is largely unknown and the major goals of this study were to describe the marine habitats and ecosystems of the Osa region, identify key ecosystem elements, and help inform proposals for increased marine protection for the region.

Materials and methods

Ethics statement

Data were collected by all authors in a collaborative effort. Non-invasive research was conducted, which included photographs, video, and visual estimates described in the methods below. The Costa Rica National System of Conservation Areas (SINAC), Osa Conservation Area office (ACOSA), and Costa Rican Fisheries and Aquaculture Institute (INCOPESCA) granted all necessary permissions to conduct this research. No vertebrate sampling was conducted and therefore no approval was required by any Animal Care and Use Committee. Our data are publicly available at Data Dryad: https://doi.org/10.5061/dryad.np5hqbzws

Habitats and sampling locations

Precipitation in the region averages 4,000 to 6,000 mm y−1 with distinct rainy (May–December) and dry seasons (January–April) [25]. Sampling was conducted from 11 to 22 March 2019, during the dry season. Habitat strata were defined a priori based on previous knowledge of the area [16, 17, 21] and included: 1) coral reefs around Isla del Caño Biological Reserve (IC), 2) coastal rocky reefs and islets along the Osa Peninsula, including Corcovado National Park (COAST), and 3) submerged pinnacles just outside the Isla del Caño MPA (PIN) (Fig 1). Much of the coastline along the Osa Peninsula consists of soft sediment exposed to strong wave action and large river inputs [26], which limited sampling in this region. The coastal rocky reefs and islets that do exist are mostly shallow (< 20 m deep) and dominated by barnacles, crustose coralline algae with few scleractinian corals [16]. In contrast, Isla del Caño, which is located ~ 15 km of the coast, has minimal coastal influences and harbors some of the best developed coral reefs in Costa Rica [19, 21]. The pinnacles just outside the Isla del Caño MPA come to within 10–15 m of the surface and represent a unique habitat within the region.

thumbnail
Fig 1. Major habitat types observed around the Osa Peninsula.

A. Coastal rocky reef, B. Roca Corcovado, a coastal islet, C. Coral reefs around Isla del Caño Biological Reserve, and D. Bajo Diablo, a submerged pinnacle outside the Isla del Caño Biological Reserve. All photos courtesy of Manu San Félix.

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

Sampling was stratified by depth strata (shallow = ~ 10 and deep = ~ 20 m). Efforts were made to sample both strata at each site when possible. However, deep habitat was not always present at some of the COAST sites and shallow sites were not always present at PIN sites, resulting in a slightly unbalanced sampling design.

Underwater visual censuses

For algae, corals, and other sessile invertebrates we used a line-point intercept methodology along each transect, recording the species or taxa found every 20 cm on the measuring tape. Five 10-m long transects were conducted at each depth strata (~ 10 and 20 m). Mobile invertebrate densities were quantified in 50 x 50 cm quadrats, with 25 quadrats per depth stratum. All individuals were identified to the lowest possible taxa.

Within each depth strata fish transects were conducted along three 25-m long transects laid along isobaths within a homogeneous habitat, with 5 m separating each transect. At each survey site, a scuba diver counted and sized all fishes of total length (TL) ≥20 cm within a 4 m wide strip on an initial ‘‘swim out” while the transect line was laid (transect area = 100 m2), and all fishes of TL <20 cm were counted in a 2 m wide strip along the transect line on the way back (transect area = 50 m2). Total fish lengths were estimated to the nearest cm. All fishes were identified to the species level based on Robertson and Allen [27] and categorized into four trophic groups: top predators (trophic level ≥ 4.0), herbivores, secondary consumers, and planktivores using FishBase (http://www.FishBase.org). Valid scientific names were verified according to Fricke et al. [28].

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 (www.fishbase.org). The sum of all individual weights and numerical densities was used to estimate biomass density by species. Species diversity was calculated using the Shannon-Wiener diversity index [29]: H’ = Σ (pi × ln pi), where pi is the proportion of all individuals counted that were of species i.

Statistical analyses

Sessile benthic taxa were grouped into functional groups for higher level analyses. These included: crustose coralline algae (CCA), scleractinian coral, soft corals (Alcyonacea), erect coralline algae, cyanobacteria, macroalgae, sponges, turf, other invertebrates (e.g., mollusks, ascidians, hydroids). Drivers of sessile benthic communities, mobile invertebrate communities, and fish assemblage structures were investigated using permutation-based multivariate analysis of variance (PERMANOVA). Bray–Curtis similarity matrices were created from percent cover of sessile benthic taxa, mobile invertebrate abundance, and biomass of fish taxa. Prior to analyses, sessile benthic taxa cover was arcsine-sqrt transformed, mobile invertebrate community structure was ln(x+1) transformed, and fish biomass was 4th root transformed. Fixed factors for the PERMANOVA were depth (shallow = 10 m, deep = 20 m) and habitat (IC, PIN, COAST).

Interpretation of PERMANOVA results was aided using individual analysis of similarities (ANOSIM), and similarity of percentages analysis (SIMPER) of species responsible for such patterns [30]. ANOSIM is a permutation-based hypothesis testing analysis of similarities that generates an R statistic that is on a scale from 0 or negative value (identical assemblages) to 1 (completely dissimilar assemblages). The resulting P value indicates the probability that the two assemblages come from a similar distribution [31]. Pairwise ANOSIM R statistics represent comparisons that are well separated (R> 0.75), overlapping but clearly different (R > 0.5), or barely separable at all (R< 0.25). SIMPER was used to determine the sessile benthic taxa, mobile invertebrate taxa, and fish species most responsible for the percentage dissimilarities between habitats using Bray-Curtis similarity analysis of hierarchical agglomerative group average clustering. Principal Coordinate Analysis (PCO) was used to compare sessile benthic, mobiles invertebrate, and fish assemblage structures among habitat types.

The number of sessile benthic taxa and fish assemblage characteristics (e.g., species richness, number of individuals, biomass, diversity) were compared using generalized linear mixed models (GLMMs) with a Poisson distribution and log-link function with habitat, depth strata, and their interaction treated as fixed factors. Site was treated as a random effect in the models. Contrasts between fixed factors were conducted using likelihood ratio tests (α = 0.05). The overdispersion parameter was estimated by Pearson’s χ2/df.

Results

Nearshore ecosystems

We conducted a total of 25 quantitative sampling surveys among three habitat types (IC, PIN, COAST) and at shallow (~ 10 m) and deep (~ 20 m) depth strata in most cases (N = 48), with 4 additional qualitative surveys along Osa Peninsula and Isla del Caño (Fig 2, S1 Table).

thumbnail
Fig 2. Locations visited on the expedition to the Osa Peninsula.

Basemap derived from GEBCO Compilation Group (2020) GEBCO 2020 Grid (doi:10.5285/a29c5465-b138-234de053-6c86abc040b9). Processing and assembly of the Global Self-consistent, Hierarchical, High-resolution Geography Database for shoreline data from [32].

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

Benthic communities

We identified a total of 95 different benthic taxa from 11 phyla, 17 classes, 35 orders, and 53 families around the Osa Peninsula (S2 Table). Turf algae accounted for 33.2% of total sessile benthic cover, followed by crustose coralline algae (CCA, 15.1%), the crustose calcareous alga Peyssonnelia sp. (7.1%), and the brown algae Dictyota humifusa (5.4%). Scleractinian hard corals accounted for 13% of total bottom cover, with Porites lobata (4.5%), Pocillopora elegans (3.9%), Porites panamensis (1.4%), and Psammocora stellata (1.2%) accounting for most of this assemblage. Soft corals (Alcyonacea) comprised 2.6% of total benthic cover with Leptogorgia alba and Carijoa riisei accounting for 78.7% of the total cover within this group. There was no significant difference in the number of benthic taxa among habitats (GLMM, χ2 = 3.483, p = 0.175), between depth strata (χ2 = 3.435, p = 0.064), or their interaction (χ2 = 3.273, p = 0.195).

Sessile benthic community structure based on cover by taxa was significantly different among major habitat types (Pseudo-F2,47−12.718, p = 0.001), with the differences between IC and PIN the most distinct (ANOSIM R = 0.856, p = 0.001), followed by PIN and COAST (ANOSIM R = 0.637, p = 0.001), with IC and COAST being most similar (ANOSIM R = 0.516, p = 0.001). Sessile benthic community structure was not significantly different between depth strata (Pseudo-F1,47−1.177, p = 0.300) or the interaction of depth and habitats (Pseudo-F2,47−1.021, p = 0.431). Sessile benthic cover by functional group also showed significant differences in community structure among habitats (Pseudo-F2,47−21.848, p = 0.001), but not between depths (Pseudo-F2,47−1.054, p = 0.358), or their interaction (Pseudo-F2,47−0.979, p = 0.423). Similar to the taxa level analysis, functional group assemblage structure was most distinct between IC and PIN (ANOSIM R = 0.860, p = 0.001), followed by COAST and PIN (ANOSIM R = 0.573, p = 0.001), and COAST and IC (ANOSIM R = 0.500, p = 0.001).

The average SIMPER dissimilarity of sessile benthic community structure between PIN and IC was 79.3%, with turf algae accounting for 21.6% of this dissimilarity and driven largely by the high cover of turf at IC ( = 38.7% ± 16.7) compared to the nearby PIN ( = 8.0% ± 17.5) (Table 1). The cover of CCA at PIN ( = 33.3 ± 16.7) was 3.5 times higher than at IC ( = 9.5 ± 8.1), contributing 16.0% to the total dissimilarity between these two habitats. Peyssonnelia sp. accounted for an additional 14.3% of the dissimilarity between these two habitats, with cover of this red calcareous alga at PIN ( = 24.9 ± 15.1) more than 10 times higher than at IC ( = 2.4 ± 3.1).

thumbnail
Table 1. Similarity of percentages (SIMPER) for benthic taxa most responsible for the percent dissimilarities among habitats using Bray-Curtis similarity analysis of hierarchical agglomerative group average clustering.

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

Benthic community structure between COAST and PIN also had high dissimilarity (68.0%), with turf algae nearly 5 times higher along the coast ( = 39.8% ± 19.1%) compared to PIN and accounted for 26.2% of the dissimilarity between these two habitats. CCA accounted for an additional 15.7% of the dissimilarity between these two habitats, with COAST cover ( = 14.5 ± 7.6) 2.3 times lower than at PIN. Dissimilarity between IC and COAST was 54.7%, with turf accounting for 18.2% of this difference although cover of this group was only 2.8% higher within COAST compared to IC. Two scleractinian hard corals (Porites lobata and Pocillopora elegans), accounted for 6.7% and 5.6%, respectively, of the dissimilarity between COAST and IC. Cover of Porites lobata was 21.0 times higher at IC ( = 7.6 ± 7.6) compared to COAST and Pocillopora elegans was 24.9 times higher at IC ( = 6.2 ± 9.0) compared to COAST.

Sampling stations within habitat types based on benthic taxa were concordant with one another and separated in ordination space by habitat (Fig 3A). PCO1 accounted for 36.1% of total variation with locations around IC clustered towards the higher end of this axis and PIN clustering towards the lower end of the axis. Peyssonnelia sp. (r = -0.836), CCA (r = -0.809) and the orange cup coral Tubastraea coccinea (r = -0.604) were all highly correlated with PIN sites, while the scleractinian coral Porites lobata (r = 0.728), turf algae (r = 0.713), and P. panamensis (r = 0.575) were highly correlated with IC sites. PCO2 explained an addition 15.7% of the variation in the benthic community. The calcareous green alga Halimeda discoidea (r = 0.595), the purple colonial tunicate Eudistoma angolanum (r = 0.583), and the sponge Aplysina cf. revillagigedoi (r = 0.575) were most closely correlated with COAST sites. Pocillopora elegans (r = -0.533) was correlated with sites around IC along PCO2 and Balanus sp. (r = -0.488) was most closely correlated with PIN stations along this same axis.

thumbnail
Fig 3. Principal coordinates analysis of benthic community composition based on percent cover.

Data were arcsine-sqrt transformed prior to analyses. Vectors are the relative contribution and direction of influence of taxa to the observed variation among sites (Pearson product-moment correlations ≥ 0.5). A. Benthic community composition by taxa, B. Benthic community composition by functional group.

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

Sampling stations within habitat types based on benthic functional groups showed greater separated in ordination space by habitat compared to examination at the taxa level (Fig 3B). PCO1 accounted for 57.5% of the total variation and was positively correlated with turf algae (r = 0.739), and scleractinian coral (r = 0.674) towards IC. PCO1 was negatively correlated with crustose coralline algae (r = -0.901), CCA (r = -0.807), and other inverts (r = -0.785). An additional 19.7% of total variation was explained by PCO2 and negatively correlated with macroalgae (r = -0.709) towards shallow sites and positively correlated with turf algae (r = 0.588) towards deeper sites.

Mobile invertebrates

There were 22 unique mobile invertebrate taxa recorded on transects from 3 phyla, 6 classes, 9 orders, and 15 families (S3 Table). Mollusca accounted for 59.1% of all taxa, followed by Echinodermata (36.4%), and Arthropoda (4.5%). The sea urchin Hesperocidaris asteriscus had the highest density among all mobile invertebrates ( = 73.3 ± 158.3 individuals 100 m-2) and occurred on 42% of transects. The long-spined sea urchin Diadema mexicanum was the most frequently encountered mobile invertebrate, occurring on 79% of transects with a density of 59.6 ± 58.6 indiv. 100 m-2. All other mobile invertebrates were in relatively low densities compared with these two sea urchin species.

Mobile invertebrate community structure was significantly different among major habitat types (Pseudo-F2,43−9.21, p = 0.001), between depths (Pseudo-F1,43−2.92, p = 0.022), and their interaction (Pseudo-F2,43−3.37, p = 0.002). The differences between IC and COAST were the most distinct (ANOSIM R = 0.530, p = 0.001), followed by IC and PIN (ANOSIM R = 0.315, p = 0.004), and PIN and COAST, which were indistinguishable from one another (ANOSIM R = 0.151, p = 0.051).

The average SIMPER dissimilarity of sessile benthic community structure between IC and COAST was 85.83%, with D. mexicanum accounting for 40.65% of this dissimilarity and H. asteriscus comprising an additional 31.58% (Table 2). Densities of D. mexicanum were 55.8% higher at IC compared with COAST while densities of H. asteriscus were 3.6 times higher at COAST sites. PIN and COAST showed 80.3% dissimilarity with H. asteriscus accounting for 44.91% of this dissimilarity, with densities at PIN sites 53.7% higher than COAST sites. IC and PIN were the least dissimilar (72.1%), but was also driven by H. asteriscus, which was 5.6 times more abundant at PIN sites compared to IC. D. mexicanum accounted for 37.7% of the dissimilarity between depths with densities in the shallow stratum 35.3% higher than deep stratum. H. asteriscus was slightly more abundant in the deep stratum ( 72.9 ± 204.4) compared to the shallow stratum ( 71.5 ± 123.1) and contributed an additional 33.2% to the dissimilarity between depth strata.

thumbnail
Table 2. Similarity of percentages (SIMPER) for mobile invertebrates most responsible for the percent dissimilarities among habitats using Bray-Curtis similarity analysis of hierarchical agglomerative group average clustering.

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

Sampling stations within PIN where concordant with one another, while IC stations separated out by depth for the most part (Fig 4). COAST stations showed the most variability in ordination space. PCO1 accounted for 39.4% of total variation with PCO2 comprising an additional 25.1. Pharia pyramidata was negatively correlated with PCO1(r = -0.413) and positively corrected with PCO2 (r = 0.675). H. asteriscus was also negatively correlated with PCO1 (r = -0.340) and PCO2 (r = -0.300). D. mexicanum was strongly negatively correlated with PCO2 (r = -0.423). H. asteriscus was orthogonal to P. pyramidata and D. mexicanum, which were opposite one another in ordination space.

thumbnail
Fig 4. Principal coordinates analysis of mobile invertebrate assemblage structure composition based on numerical density (number of individuals 100 m-2).

Data were ln(x+1)-transformed prior to analyses. Vectors are the relative contribution and direction of influence of taxa to the observed variation among sites (Pearson product-moment correlations ≥ 0.5).

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

Fishes

We observed a total of 129 different fish taxa during the expedition from 49 families (S4 Table). Of these, 112 taxa were recorded on quantitative transects. The most taxa rich families were Labridae (N = 10), Carangidae (N = 10), Lutjanidae (N = 8), Epinephelidae (N = 7), and Haemulidae (N = 7). Bluestriped chubs (Kyphosus ocyurus) were the most important species by weight, accounting for 16% of the total biomass but only occurring at 15% of the stations (S4 Table). Yellowfin surgeonfish, Acanthurus xanthopterus, was the next most important species by weight and was present at 58% of the stations. The razor surgeonfish (Prionurus laticlavius) and bigeye trevally (Caranx sexfasciatus) each comprised ~ 7% of total biomass. However, the razor surgeonfish was observed much more frequently (69% of the stations) compared with the bigeye trevally, which only was present at 23% of the stations.

Overall mean fish taxa richness per transect was 20.4 (± 5.3) and was significantly higher at PIN compared to COAST and IC, which were not significantly different from one another (Table 3, Fig 5). Numerical abundance was significantly higher at PIN ( = 7.1 ± 2.5), compared to IC ( = 4.2 ± 2.2) and COAST (3.0 ± 1.6), which were not significantly different from one another. Fish biomass was significantly higher at PIN by 4.6 times ( = 7.9 ± 11.4) compared to IC ( = 1.7 ± 2.6) and 3.6 times higher than COAST ( = 2.2 ± 5.9).

thumbnail
Fig 5. Comparisons of fish assemblage characteristics among habitats.

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

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

thumbnail
Table 3. Comparisons of fish assemblage characteristics among habitats and depth strata.

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

Biomass between COAST and IC did not differ significantly from each other. Owing to its small size, only two stations were sampled at Islote Matapalo along the coast. However, this isolated rock had the highest biomass (8.9 ± 13.1 g m-2) observed during our surveys and was comprised of large schools of bluestriped chub (K. ocyurus), which accounted for 27% of the total biomass at this site. This was followed by the razor surgeonfish (P. laticlavius) and bigeye trevally (C. sexfasciatus), which accounted for an additional 22% and 18%, respectively, of the biomass at Islote Matapalo. Although fish biomass was more than 2 times higher in the deep depth stratum (4.5 ± 8.3) compared to the shallow depth stratum (2.2 ± 2.0) these differences were not significant. Diversity did not differ significantly among habitats, between depths, or their interaction. However, this metric showed the highest variability within PIN habitats and the least within COAST habitats.

The highest dissimilarity in fish assemblages based on taxa biomass was between COAST and PIN (SIMPER average dissimilarity = 92.0), with the greatest contribution to this dissimilarity coming from C. sexfasciatus (8.4%), followed by Haemulon steindachneri (8.0%), and K. ocyurus (7.9%) (Table 4). All these schooling species had much higher biomass at PIN stations compared to COAST. IC and COAST had an average dissimilarity of 88.1% and was driven by Scarus rubroviolaceus (9.6%), which had biomass 14.3 times higher at IC compared to COAST. Two species of schooling surgeonfishes, P. laticlavius and A. xanthopterus, accounted for an additional 7.3% and 7.0% of the dissimilarity between groups, respectively, with A. xanthopterus 15.2 times more abundant by weight at IC and P. laticlavius 3.5 times more abundant by weight at IC compared to COAST. Dissimilarity between PIN and IC was 85.7% and was driven by H. steindachneri (8.4%), which had a biomass of 70.5 g m-2 (± 176.1) at PIN but was not present on transects at IC. Triaenodon obesus was 3.7 times more abundant by weight at PIN vs. IC and contributed 7.9% to the dissimilarity between these groups. Caranx sexfasciatus accounted for an additional 6.6% of the dissimilarity between these two habitats and had a biomass of 63.4 g m-2 (± 155.5) at PIN stations but was not observed at IC. Similarly, biomass of K. ocyurus was three orders of magnitude greater at PIN 159.2 (± 551.4) vs. IC (0.8 ± 0.8) stations.

thumbnail
Table 4. Similarity of percentages (SIMPER) for fish biomass most responsible for the percent dissimilarities among habitats using Bray-Curtis similarity analysis of hierarchical agglomerative group average clustering.

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

Fish taxa biomass assemblage structure was significantly different among habitats (pseudo-F2,54 = 5.64, p < 0.001) and between depths (pseudo-F1,42 = 2.14, p = 0.007), but not their interaction (pseudo-F2,54 = 1.26, p = 0.138). Fish assemblage structure between COAST and IC was most distinct (ANOSIM R = 0.613, p = 0.001), while IC and PIN were most similar (ANOSIM R = 0.381, p = 0.001). Fish assemblage structure based on taxa biomass showed separation in ordination space among habitats based on Principal Coordinates Analysis (Fig 6A). The 1st PCO axis explained 17.9% of the variation in assemblage structure. This separation was driven by Arothron meleagris (r = -0.594) and Bodianus diplotaenia (r = -0.492), which were negatively correlated with PCO1 and trended towards PIN sites. PCO2 accounted for an additional 13.5% of the variation in fish assemblage structure with Chaetodon humeralis (r = 0.591), Anisotremus taeniatus (r = 0.577), and Caranx sexfasciatus (r = 0.470) the most highly correlated with this axis.

thumbnail
Fig 6. Principal coordinates analysis of fish assemblage structure based on biomass (g m-2).

Data were 4th root-transformed prior to analyses. Vectors are the relative contribution and direction of influence of taxa to the observed variation among sites (Pearson product-moment correlations ≥ 0.5). A. Fish assemblage structure based on taxa; B. Fish assemblage structure based on trophic group.

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

Fish trophic biomass assemblage structure was significantly different among habitats (pseudo-F2,54 = 7.42, p = 0.001), but not between depths (pseudo-F1,54 = 1.82, p = 0.141) or their interaction (pseudo-F2,54 = 0.82, p = 0.553). COAST vs. PIN trophic assemblages were overlapping but clearly different (R = 0.352), while COAST vs. IC and IC vs. PIN assemblages were barely separable (both R< 0.25). Trophic biomass assemblage structure showed separation in ordination space by habitats and depth based on Principal Coordinates Analysis (Fig 6B). The 1st PCO axis explained 58.1% of the variation in assemblage structure and was positively correlated with all four trophic groups. Top predators had the strongest correlation with PCO1 (r = 0.561), followed by secondary consumers (r = 0.556), herbivores (r = 0.477), and planktivores (r = 0.400). PIN stations had high concordance and trends towards the higher end of PCO1. COAST stations trends towards the lower end of PCO1 but were more widely dispersed in ordination space.

Discussion

The Osa Peninsula possesses the last large stand of tropical moist forest on the Pacific Coast of Mesoamerica and this nearly virgin rainforest extends all the way to the Pacific Ocean. Due to the high rainfall and steep terrain of Osa there are strong land-sea linkages that can have both positive (e.g., nutrient enrichment) and negative (e.g., sedimentation) effects on the marine environment. Much of the coastline around Osa is soft sediment and the extensive beach habitat provides critical nesting habitat for three species of threatened and endangered sea turtles–the Pacific green (Chelonia mydas), leatherback (Dermochelys coriacea), and olive ridley (Lepidochelys olivacea). The coastal rocky reefs of Osa have evolved in a high sediment and nutrient rich environment, which are dominated by turf, macroalgae, and crustose corallines, with low coral cover. The coastal islets and emergent rocks are exposed to high wave with benthic communities dominated by encrusting organisms and a relatively high cover of gorgonians, which are adapted to these harsh environmental conditions.

The lower reaches of the streams and rivers of Corcovado NP harbor a diverse assemblage of fishes, including a number of marine species [33], many of which were found on the coastal and offshore reefs during our study. These species feed on the abundant crustaceans, mollusks, and small fishes found in these habitats, which also serve as important nurseries for the juveniles of a number of these species. This strong connection between freshwater and marine environments of this region highlights the importance of healthy land/sea connections for ecosystem productivity in the area.

Approximately 15 km offshore from the peninsula is Isla del Caño, which has some of the most extensive coral reefs in the region [34] and the highest species richness of corals in Costa Rica [17, 21]. These reefs, like many in the region, suffered high mortality during the 1982–1983 El Niño event, but appear to have recovered to a large extent [21, 23, 35]. This unique ecosystem is a reservoir of biodiversity that may be resilient to near-term climate change.

The submerged pinnacles just outside the Isla del Caño Biological Reserve harbored unique invertebrate communities that includes the recent discovery of a new gorgonian [24], with high fish biomass, taxonomic richness, and abundance compared to nearby Isla del Caño and the coast. Fish biomass at these pinnacles was higher than most unprotected locations in the ETP [36] and comparable to iconic highly protected offshore MPAs in the region like Cocos Island National Park [37] and Malpelo [36, 38], as well as remote unfished reefs in the Pacific and Indian oceans [39]. Overall fish species richness from our quantitative surveys (N = 129) is higher than those reported from other parts of Costa Rica [40, 41] and similar to some well-studied areas in the region [42, 43].

Reefs in Golfo Dulce, a better studied nearby area, have been impacted by natural and anthropogenic disturbances. The decline of coral reefs in Golfo Dulce has resulted from increased sediments loads due to deforestation, slash-burn agriculture, road construction, and mining [3, 13]. There have also been increases in freshwater input, extreme ENSO events, and intense dinoflagellate blooms in the area [44].

The main source of sedimentation to the western portion of Osa and Isla del Caño is the Térraba River, which runs through the Térraba- Sierpe National Wetland—the largest mangrove estuary in Costa Rica [4, 16]. This ecosystem is an important nursery habitat for fishes and invertebrates that migrate to adjacent habitats along the coast as they mature [45, 46]. The considerable contribution of terrestrial sediments to the nearshore marine ecosystem has accelerated in recent years due to erosion of exposed lands, as well as coastal alteration [47]. These fine-grained sediments readily sticks to coral tissue and are more easily resuspended due to their lower density, contributing to reductions in water clarity for prolonged periods, locally reducing pH and oxygen conditions [48]. Coral cover along the coast of the Osa Peninsula is low and some locations appear to have suffered mortality from sedimentation [16]. Also, levels of bioerosion (destruction of the coral skeleton by organisms) are high in this region [13] and these coral reefs are being destroyed faster than they are accreting [49, 50].

Despite Isla del Caño Biological Reserve being fully protected from fishing by law, biomass was similar to fished areas along the coast and lower than the adjacent pinnacles, which are outside of the protected area. Illegal fishing and derelict fishing gear have been observed inside the reserve, particularly on the east side of the island, away from the ranger station. Owing to its high biodiversity and uniqueness, Isla del Caño experiences intensive ecotourism and these activities have the potential to negatively impact the reefs and associated fauna if not managed correctly.

Currently, the Corcovado National Park only include 500 m of ocean along the coast (20 km2 vs 422 km2 of protected land). Some of the main conservation interests of the national park are the high density of bull sharks (Carcharhinus leucas) in the estuarine environments and the cetaceans that transit the area [15]. Several rocks, islets, and islands (e.g., Islote Matapalo, Roca Corcovado) are currently not included within the park boundaries and these unique features represent important biodiversity hotspots, which should be included within the park and protected from fishing. There is limited enforcement of the marine portion of Corcovado NP, and we noted several instances of fishing within the park during our short expedition. Rough sea conditions limit fishing by boat to some extent so increased enforcement of the coastal zone could benefit local fish and invertebrate populations.

Conservation

The unique configuration of healthy offshore coral reefs and pinnacles connected to coastal habitats around the Osa Peninsula provides corridors for many species including large predators such as sharks and other marine megafauna, which warrants additional protection. This is part of a larger corridor of MPAs that includes Isla del Coco National Park, the Galápagos Marine Reserve, Coiba National Park, and Malpelo Fauna and Flora Sanctuary. Oviedo and Solís [51] emphasize the need to ensure connectivity between key sites for marine mammals like humpback whales (Megaptera novaengliae), protecting a corridor between Isla del Caño and Drake Bay, as well as expanding the protected waters of Corcovado National Park. Also, the importance of both Golfo Dulce and Térraba-Sierpe National Wetland as reproductive sites for the scalloped hammerhead (S. lewini) have been demonstrated [12], and the southern Pacific coast of Costa Rica in general has been recognized for its important to endangered and threatened sea turtles [52].

The submerged pinnacles just outside the Isla del Caño Biological Reserve harbor extremely high fish biomass and unique biodiversity yet they are currently not protected within the Isla del Caño Biological Reserve, and fishing was evident as we found fishing lines on some of our dives. A logical conservation measure would be to expand the reserve to include these important ecological features. In addition, there is an urgent need for integrated ridge to reef management to reduce the impacts of coastal and road construction, discharge of effluent and pollutants into the coastal environment, illegal logging, and mining adjacent to Corcovado National Park.

Osa Peninsula is at a critical juncture. It is the last remaining section of Costa Rica’s Pacific coast where eco-tourism is the dominant type of tourism and a significant sector of the local economy [53]. Therefore, appropriate small-scale development must be balanced with sound environmental practices to provide sustainable social and economic opportunities. The area is relatively lightly populated, but previous land use and resource use practices have shifted these ecosystems into a compromised state, which from a water quality and marine resources perspective are far less productive than they might be under different future planning scenarios. A better understanding of ridge-to-reef connectivity is critically needed to ensure better management of nearshore ecosystem health [5456]

Fisheries in the region are mainly artisanal, using hook-and-line, long-line, and trammel nets to target snappers, mackerel, and grunts [57]. However, there is a high bycatch of sharks, bony fishes, and turtles associated with the Costa Rican longline fishery and high catch rates of juvenile blacktip sharks (Carcharhinus limbatus) near the Osa Peninsula indicates that Osa is an important nursery area for this species [58]. Protection is urgently needed from this fishery in this area, which would also protect green and olive ridley turtles that nest on the Osa Peninsula. Shrimp trawling in Costa Rica has been shown to have significant environmental impacts on marine habitats and populations, including the bycatch of sub-reproductive sized fishes, as well as threatened and endangered species such as sea turtles [59, 60]. Shrimp trawling in Costa Rica is extremely unpopular with artisanal fishers and the creation of an MPA in the region that excluded this practice could benefit local fisheries [61]. Climate change and ENSO events have already impacted this region and while MPAs and integrated ridge to reef management cannot directly affect these processes, they can reduce local stressors and increase ecosystem resilience and resistance to climate change.

Considered one of the twenty-five most biodiverse countries of the planet, Costa Rica has earned global recognition for its extraordinary national parks and conservation of nature. Costa Rica is a founding member and co-chair of the High Ambition Coalition for Nature and People, whose objective is to set a global target to protect 30% of the planet’s sea and land by year 2030. Currently ~11% of the EEZ of Costa Rica is fully protected from fishing and other extractive activities, mostly within the Cocos Island National Park, which was recently expanded on December 17, 2021. The creation of a large MPA in the Osa region that includes the existing MPAs with greater protection and enforcement would benefit the rich biodiversity of this part of the country as well as replenishing nearby overexploited important fisheries resources.

Supporting information

S1 Table. Metadata for surveys conducted around Osa Peninsula, south Pacific Costa Rica.

Geomor.–Geomorphic habitat types: sites around Isla del Caño (IC), 2) pinnacles around Isla del Caño (PIN), and 3) coastal rocky reefs and nearshore islets (COAST). Mangt.–Management type (no-take marine protected area, open to fishing). Stat.–stations (99# are qualitative surveys). CR–coral reef, P–pinnacle, IS–coastal islet, RR–rocky reef. No. benthic–number of benthic transects. No. fish–number of fish transects.

https://doi.org/10.1371/journal.pone.0271731.s001

(DOCX)

S2 Table. Benthic taxa identified during the expedition to Osa Peninsula.

https://doi.org/10.1371/journal.pone.0271731.s002

(DOCX)

S3 Table. Macro invertebrate taxa observed on quantitative transects around Osa Peninsula.

https://doi.org/10.1371/journal.pone.0271731.s003

(DOCX)

S4 Table. Fish species list from Osa Peninsula.

Biomass based on quantitative surveys. See methods for a description. %Freq.–percent frequency of occurrence (n = 26 stations). Data type–Quan–quantitative, Qual–qualitative.

https://doi.org/10.1371/journal.pone.0271731.s004

(DOCX)

Acknowledgments

The expedition to the Osa Peninsula was made possible thanks to the support of the Ministry of Environment and Energy, particularly Minister Carlos Manuel Rodríguez. We would also like to thank Vice Minister of Water and Seas Haydée Rodríguez Romero for additional support. Numerous government and non-government organizations provided assistance during our expedition, including: the Centro de Investigación en Ciencias del Mar y Limnología, Universidad de Costa Rica, the National System of Conservation Areas (SINAC), Osa Conservation Area (ACOSA), Caño Island Biological Reserve Park Rangers (RBIC), Costa Rican Institute of Fisheries and Aquaculture (INCOPESCA), the Municipality of Osa, Lodge Punta Marenco, the Center for Research in Marine Sciences and Limnology (CIMAR) of the University of Costa Rica, the captain and crew of the Argo, Corcovado National Park, the Wyss Foundation, and the artisanal fishing communities of the Osa Peninsula. Arturo Angulo, Juan José Alvarado Barrientos, and Cindy Fernández reviewed the taxonomic lists for accuracy. This is contribution no. 1896from the Hawai‘i Institute of Marine Biology and no. 11537from the School of Ocean and Earth Science and Technology at the University of Hawai‘i.

References

  1. 1. Gutiérrez BL, Zambrano AMA, Zambrano SLA, Gil CAQ, Bohlman S, Arias EA, et al. An island of wildlife in a human-dominated landscape: The last fragment of primary forest on the Osa Peninsula’s Golfo Dulce coastline, Costa Rica. PLOS ONE. 2019;14: e0214390. pmid:30913255
  2. 2. Sak PB, Fisher DM, Gardner TW. Effects of subducting seafloor roughness on upper plate vertical tectonism: Osa Peninsula, Costa Rica. Tectonics. 2004;23.
  3. 3. Quesada-Alpízar MA, Cortés J. Los ecosistemas marinos del Pacífico sur de Costa Rica: estado del conocimiento y perspectivas de manejo. Revista de Biología Tropical. 2006;54(Suppl 1): 101–145.
  4. 4. Silva-Benavides AM. Mangroves. In: Wehrtmann IS, Cortés J, editors. Marine Biodiversity of Costa Rica, Central America. Dordrecht: Springer Netherlands; 2009. pp. 123–130. https://doi.org/10.1007/978-1-4020-8278-8_7
  5. 5. Hebbeln D, Cortés J. Sedimentation in a tropical fjord: Golfo Dulce, Costa Rica. Geo-Marine Letters. 2001;20: 142–148.
  6. 6. Córdoba-Muñoz R, Vargas-Zamora JA. Temperature, salinity, oxygen, and nutrient profiles at a 200 m deep station in Golfo Dulce, Pacific coast of Costa Rica. Revista de Biología Tropical. 1996;44(Suppl 3): 233–236.
  7. 7. Morales-Ramírez Á. La diversidad marina del Golfo Dulce, Pacífico sur de Costa Rica: amenazas de conservación. Biocenosis. 2011;24. Available: https://revistas.uned.ac.cr/index.php/biocenosis/article/view/1198
  8. 8. Stief P, Lundgaard ASB, Morales-Ramírez Á, Thamdrup B, Glud RN. Fixed-Nitrogen loss associated with sinking zooplankton carcasses in a coastal oxygen minimum zone (Golfo Dulce, Costa Rica). Frontiers in Marine Science. 2017;4. Available: https://www.frontiersin.org/article/10.3389/fmars.2017.00152
  9. 9. Wolff M, Hartmann HJ, Koch V. A pilot trophic model for Golfo Dulce, a fjord-like tropical embayment, Costa Rica. Revista de Biología Tropical. 1996;44(Suppl 3): 215–231.
  10. 10. Bessesen BL, González-Suárez M. The value and limitations of local ecological knowledge: Longitudinal and retrospective assessment of flagship species in Golfo Dulce, Costa Rica. People and Nature. 2021;3: 627–638.
  11. 11. Sistema Nacional de Áreas de Conservación (SINAC), Ministerio de Ambiente, Energía y Telecomunicaciones (MINAET). GRUAS II: Propuesta de Ordenamiento Territorial para la conservación de la biodiversidad de Costa Rica. Volumen 3: Análisis de vacíos en la representatividad e integridad de la biodiversidad marina y costera. 2008 p. 60.
  12. 12. Zanella I, López-Garro A, Cure K. Golfo Dulce: critical habitat and nursery area for juvenile scalloped hammerhead sharks Sphyrna lewini in the Eastern Tropical Pacific Seascape. Environmental Biology of Fishes. 2019;102: 1291–1300.
  13. 13. Fonseca AC, Alvarado JJ, Cortés J. Planificación marino-costera del Área de Conservación Osa (ACOSA), Costa Rica. San José, Costa Rica; 2009. Report No.: Serie Técnica: Apoyando los esfuerzos en el manejo y protección de la biodiversidad tropical, No.: 12.
  14. 14. Rivera CJ. Facing the 2013 gold rush: a population viability analysis for the endangered white-lipped peccary (Tayassu pecari) in Corcovado National Park, Costa Rica. Natural Resources. 2014;05: 1007.
  15. 15. Salas E, Ross Salazar E, Arias A. Diagnóstico de áreas marinas protegidas y áreas marinas para la pesca responsable en el Pacífico costarricense. San José, Costa Rica: Fundación MarViva; 2012 p. 174.
  16. 16. Cortés J, Jiménez CE. Coastal-marine environments of Parque Nacional Corcovado, Puntarenas, Costa Rica. Revista de Biología Tropical. 1996;44(Suppl 3): 35–40.
  17. 17. Cortés J. The Pacific Coastal and Marine Ecosystems. In: Kappelle M, editor. Costa Rican Ecosystems. University of Chicago Press; 2016. pp. 97–138. https://doi.org/10.7208/9780226121642-010
  18. 18. Guzmán HM. Estructura de la comunidad arrecifal de la Isla del Caño, Costa Rica, y el efecto de perturbaciones naturales severas. M.Sc. Thesis, Universidad de Costa Rica. 1986.
  19. 19. Guzmán HM, Cortés J. Coral reef community structure at Caño Island, Pacific Costa Rica. Marine Ecology. 1989;10: 23–41.
  20. 20. Cortés J, Jiménez C. Corals and coral reefs of the Pacific of Costa Rica: history, research and status. In: Cortés J, editor. Latin American Coral Reefs. Amsterdam: Elsevier Science; 2003. pp. 361–385. https://doi.org/10.1016/B978-044451388-5/50017-5
  21. 21. Cortés J, Jiménez CE, Fonseca AC, Alvarado JJ. Status and conservation of coral reefs in Costa Rica. Revista de Biología Tropical. 2010;58(Suppl 1): 33–50. pmid:20873039
  22. 22. Glynn P, Cortés J, Guzmán H, Richmond RH. El Niño (1982–83) associated coral mortality and relationship to sea surface temperature deviations in the tropical eastern Pacific. Proceedings of the 6th International Coral Reef Symposium. Townsville, Australia; 1988. pp. 237–243.
  23. 23. Guzmán HM, Cortés J. Changes in reef community structure after fifteen years of natural disturbances in the eastern Pacific (Costa Rica). Bulletin of Marine Science. 2001;69: 133–149.
  24. 24. Breedy O, Guzman HM, Murillo-Cruz C, Vargas S. Discovery of a new species of Leptogorgia Milne Edwards & Haime, 1857 (Anthozoa: Octocorallia: Gorgoniidae) from eastern tropical Pacific mesophotic reefs. Marine Biodiversity. 2021;51: 95.
  25. 25. Matlaga TJH. Mechanisms underlying the occurrence of species in complex modified tropical landscapes: a case study of amphibians in the Osa Peninsula, Costa Rica. Journal of Tropical Ecology. 2018;34: 32–40.
  26. 26. Marshall JS. The Geomorphology and Physiographic Provinces of Central America. In: Bundschuh J, Alvarado G, editors. Central America: Geology, Resources and Hazards. 2007. pp. 75–122.
  27. 27. Robertson DR, Allen GR. Shorefishes of the Tropical Eastern Pacific: online information system. Version 2.0. In: Shorefishes of the Tropical Eastern Pacific: online information system. Version 2.0 [Internet]. 2015 [cited 27 Jan 2022]. Available: https://biogeodb.stri.si.edu/sftep/en/pages
  28. 28. Fricke R, Eschmeyer WN, Van Der Laan R. Eschmeyer’s Catalog of Fishes: Genera, Species, References. In: Eschmeyer’s Catalog of Fishes: Genera, Species, References [Internet]. 11 Nov 2014 [cited 27 Jan 2022]. Available: http://researcharchive.calacademy.org/research/ichthyology/catalog/fishcatmain.asp
  29. 29. Ludwig JA, Reynolds JF. Statistical ecology: a primer on methods and computing. New York, New York: John Wiley and Sons; 1988.
  30. 30. Clarke KR. Non-parametric multivariate analyses of changes in community structure. Australian Journal of Ecology. 1993;18: 117–143.
  31. 31. Clarke KR, Gorley RN, Somerfield PJ, Warwick RM. Change in marine communities: an approach to statistical analysis and interpretation, 3rd edn. Plymouth: Primer-E Ltd; 2014. Available: http://plymsea.ac.uk/id/eprint/7656/
  32. 32. Wessel P, Smith WHF. A global, self-consistent, hierarchical, high-resolution shoreline database. Journal of Geophysical Research B: Solid Earth. 1996;101: 8741–8743.
  33. 33. Winemiller KO. An introduction to the freshwater fish communities of Corcovado National Park, Costa Rica. Brenesia. 1983;21: 47–66.
  34. 34. Glynn PW, Alvarado JJ, Banks S, Cortés J, Feingold JS, Jiménez C, et al. Eastern Pacific Coral Reef Provinces, Coral Community Structure and Composition: An Overview. In: Glynn PW, Manzello DP, Enochs IC, editors. Coral Reefs of the Eastern Tropical Pacific: Persistence and Loss in a Dynamic Environment. Dordrecht: Springer Netherlands; 2017. pp. 107–176. https://doi.org/10.1007/978-94-017-7499-4_5
  35. 35. Cortés J. A history of marine biodiversity scientific research in Costa Rica. In: Wehrtmann IS, Cortés J, editors. Marine Biodiversity of Costa Rica, Central America. Dordrecht: Springer Netherlands; 2009. pp. 47–80. https://doi.org/10.1007/978-1-4020-8278-8_2
  36. 36. Edgar GJ, Banks SA, Bessudo S, Cortés J, Guzmán HM, Henderson S, et al. Variation in reef fish and invertebrate communities with level of protection from fishing across the Eastern Tropical Pacific seascape. Global Ecology and Biogeography. 2011;20: 730–743.
  37. 37. Friedlander AM, Zgliczynski BJ, Ballesteros E, Aburto-Oropeza O, Bolaños A, Sala E. The shallow-water fish assemblage of Isla del Coco National Park, Costa Rica: Structure and patterns in an isolated, predator-dominated ecosystem. Revista de Biologia Tropical. 2012;60(Suppl 3): 321–338.
  38. 38. Quimbayo JP, Mendes TC, Kulbicki M, Floeter SR, Zapata FA. Unusual reef fish biomass and functional richness at Malpelo, a remote island in the Tropical Eastern Pacific. EnvironmentL Biology of Fishes. 2017;100: 149–162.
  39. 39. McClanahan TR, Schroeder RE, Friedlander AM, Vigliola L, Wantiez L, Caselle JE, et al. Global baselines and benchmarks for fish biomass: Comparing remote reefs and fisheries closures. Marine Ecology Progress Series. 2019;612: 167–192.
  40. 40. Phillips PC, Perez-Cruet MJ. A comparative survey of reef fishes in Caribbean and Pacific Costa Rica | Revista de Biología Tropical. Revista de Biología Tropical. 1984;32: 95–102.
  41. 41. Dominici-Arosemena A, Brugnoli-Olivera E, Núñez JC, Molina-Ureña H, Quesada-Alpizar M. Community structure of eastern Pacific reef fishes (Gulf of Papagayo, Costa Rica). Tecnociencia. 2005;7: 19–41.
  42. 42. Dominici-Arosemena A, Wolff M. Reef fish community structure in the Tropical Eastern Pacific (Panamá): living on a relatively stable rocky reef environment. Helgolander Marine Research. 2006;60: 287–305.
  43. 43. Aburto-Oropeza O, Balart EF. Community Structure of Reef Fish in Several Habitats of a Rocky Reef in the Gulf of California. Marine Ecology. 2001;22: 283–305.
  44. 44. Cortés J. The coral reefs of Golfo Dulce, Costa Rica: distribution and community structure. Atoll Research Bulletin. 1990;344: 1–37.
  45. 45. Chicas FA. Peces juveniles en una poza de marea, Reserva Forestal Térraba-Sierpe, Puntarenas, Costa Rica. Revista de Biología Tropical. 2001; 307–314. pmid:15264544
  46. 46. Sibaja-Cordero JA, Echeverría-Sáenz S. Polychaetes associated with decaying wood in Térraba mangrove, South Pacific, Costa Rica. Revista de Biología Tropical. 2015;63: 61–74.
  47. 47. Cortés J, Fonseca AC, Hebbeln D. Bottom topography and sediments around Isla del Caño, Pacific of Costa Rica. Revista de Biología Tropical. 1996;44(Suppl 3): 11–17.
  48. 48. Bainbridge Z, Lewis S, Bartley R, Fabricius K, Collier C, Waterhouse J, et al. Fine sediment and particulate organic matter: A review and case study on ridge-to-reef transport, transformations, fates, and impacts on marine ecosystems. Marine Pollution Bulletin. 2018;135: 1205–1220. pmid:30301020
  49. 49. Fonseca AC. Bioerosión y bioacreción en arrecifes coralinos del Pacífico sur de Costa Rica. MSc Thesis, Universidad de Costa Rica. 1999.
  50. 50. Fonseca AC, Cortés J. Coral borers of the eastern Pacific: the sipunculan Aspidosiphon (A.) elegans and the crustacean Pomatogebia rugosa. Pacific Science. 1998;52: 170–175.
  51. 51. Oviedo L, Solís M. Underwater topography determines critical breeding habitat for humpback whales near Osa Peninsula, Costa Rica: implications for Marine Protected Areas. Revista de Biología Tropical. 2008;56: 591–602. pmid:19256430
  52. 52. Venegas-Li R, Palacios Alfaro JD, Martínez Fernández D. Coastal aggregation of marine turtles in waters of the South Pacific of Costa Rica. Revista de Biología Tropical. 2015;63(Suppl 1): 1–8.
  53. 53. Hunt CA, Durham WH, Driscoll L, Honey M. Can ecotourism deliver real economic, social, and environmental benefits? A study of the Osa Peninsula, Costa Rica. Journal of Sustainable Tourism. 2015;23: 339–357.
  54. 54. Carlson RR, Foo SA, Asner GP. Land Use Impacts on Coral Reef Health: A Ridge-to-Reef Perspective. Frontiers in Marine Science. 2019;6. Available: https://www.frontiersin.org/article/10.3389/fmars.2019.00562
  55. 55. Delevaux JMS, Whittier R, Stamoulis KA, Bremer LL, Jupiter S, Friedlander AM, et al. A linked land-sea modeling framework to inform ridge-to-reef management in high oceanic islands. PLOS ONE. 2018;13: e0193230. pmid:29538392
  56. 56. Fache E, Pauwels S. The ridge-to-reef approach on Cicia Island, Fiji. Ambio. 2022 [cited 28 Apr 2022]. pmid:35018620
  57. 57. Guzmán-Mora AG, Molina-Ureña H. La Pesquería Artesanal de Peces en Golfo Dulce, Costa Rica. Proceedings of the 60th Gulf and Caribbean Fisheries Institute. 2008;60: 137–142.
  58. 58. Dapp D, Arauz R, Spotila JR, O’Connor MP. Impact of Costa Rican longline fishery on its bycatch of sharks, stingrays, bony fish and olive ridley turtles (Lepidochelys olivacea). Journal of Experimental Marine Biology and Ecology. 2013;448: 228–239.
  59. 59. Campos JA, Burgos B, Gamboa C. Effect of shrimp trawling on the commercial ichthyofauna of the Gulf of Nicoya, Costa Rica. Revista de Biología Tropical. 1984;32: 203–207.
  60. 60. Clarke T, Espinoza M, Ahrens R, Wehrtmann I. Elasmobranch bycatch associated with the shrimp trawl fishery off the Pacific coast of Costa Rica, Central America. Fishery Bulletin- National Oceanic and Atmospheric Administration. 2015;114.
  61. 61. Molina Bustamante C, Crowder L. Human dimension of the marine ecosystems of the Osa and Golfito Region. San José, Costa Rica: Osa and Golfito Initiative, INOGO, Stanford Woods Institute for the Environment, Stanford University; 2014.