
==== Front
PeerJ
PeerJ
PeerJ
PeerJ
2167-8359
PeerJ Inc. San Diego, USA

17763
10.7717/peerj.17763
Aquaculture, Fisheries and Fish Science
Ecology
Ecosystem Science
Marine Biology
Benthic and fish community composition on mesophotic reefs in Grand Cayman
Le Gall Lucas
Johnson Jack V.
Chequer Alex
http://orcid.org/0009-0000-6259-3006
Doherty Matthew Louis
Goodbody-Gringley Gretchen ggoodbody@reefresearch.org

Reef Ecology and Evolution Lab, Central Carribean Marine Institute, Little Cayman, Cayman Islands
Aguilar-Perera Alfonso
29 8 2024
2024
12 e1776327 3 2024
26 6 2024
© 2024 Le Gall et al.
2024
Le Gall et al.
https://creativecommons.org/licenses/by/4.0/ This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, reproduction and adaptation in any medium and for any purpose provided that it is properly attributed. For attribution, the original author(s), title, publication source (PeerJ) and either DOI or URL of the article must be cited.

Mesophotic Coral Ecosystems (MCEs) represent unique ecological habitats that range from 30 to 150 m deep, harbouring phylogenetically distinct species and offering refuge for many taxa during times of environmental stress. Yet owing to inaccessibility of ecosystems at these depths, most MCEs remain unexplored, with quantifications of ecological communities in these habitats lacking across many regions. Here, using open- and closed-circuit technical diving, we quantified benthic and fish community composition at four mesophotic reef sites (45 m depth) in Grand Cayman. We show significant differences in benthic community composition over a small spatial scale driven by disparate coverage of sponges, crustose coralline algae, and sand/rubble, yet consistent patterns of macroalgal dominance representing >50% coverage at each site and low hard coral cover at an average of 2.4%. Reef fish species richness, biomass, and density was consistent across sites, however the relative contribution of individual species to community composition differed significantly. Macrocarnivores were found to be the dominant contributors to biomass, with invertivores the most speciose, and omnivores and planktivores at the highest densities, consistent with previous descriptions of mesophotic fish assemblages in other regions. Similarly, the low hard coral cover and high macroalgae and sponge cover of the benthic communities also appear ecologically similar to several described mesophotic reefs yet is not uniform across the Caribbean. The ecological organisation of Grand Cayman’s MCEs may result from a variety of factors such as isolation from other major land masses, geology, local geography, and anthropogenic activity at both the local and global scale and highlight the importance of continued exploration and documentation of MCE communities.

Mesophotic reef
Benthic
Carribean
Fish community
Cayman Islands
Mesophotic
Coral reef
Assemblages
Macroalgae
Benthic cover
Darwin Plus InitiativeRef: DPlus 162 This work was supported by the Darwin Plus Initiative (Ref: DPlus 162) and a matching fund made by an anonymous private donor. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
==== Body
pmcIntroduction

Mesophotic Coral Ecosystems (MCEs) represent globally distributed, deeper water ecosystems, typically found between depths of 30 and 150 m (Pyle & Copus, 2019). Comprised of light-dependent scleractinian corals, these ecosystems also contain a variety of taxa including sponges, macroalgae, encrusting algae, and soft coral colonies (Hinderstein et al., 2010). Despite the ecological importance of MCEs for harbouring marine biodiversity, which contributes to coastal ecosystem services (Holstein et al., 2019), MCEs have historically received less attention than their shallow water counterparts due primarily to limited access using traditional SCUBA. However, as diving and submersible technology has advanced, so has access to the mesophotic, resulting in an increase in scientific investigation and biodiversity assessments (Pyle & Copus, 2019). This recent surge has primarily focused on Caribbean MCEs, which are typically located in regions with high anthropogenic activity, comprising roughly 55% of MCE focused publications as of 2018 (Laverick et al., 2018). These studies have been crucial in documenting the biodiversity of MCEs, highlighting the ecologically and biologically distinct nature of these systems (Kahng et al., 2010) with prominent differences in both benthic (Lesser, Slattery & Leichter, 2009; Stefanoudis et al., 2019) and fish assemblages (García-Sais et al., 2008; Lesser, Slattery & Leichter, 2009; Pinheiro et al., 2016; Stefanoudis et al., 2019) compared to shallow water ecosystems. Within MCEs, observed breaks in biodiversity are documented at approximately 60 m, distinguishing the upper and lower mesophotic zones (Lesser, Slattery & Leichter, 2009). While the defined depth range is influenced by topography, light intensity is also a major driving factor affecting community composition and defining species boundaries (Laverick, Andradi-Brown & Rogers, 2017; Pérez-Castro et al., 2022). Such physical factors of the environment differ among and within locations, however, resulting in strong differences in communities and distributions among various MCEs, warranting continued assessments of community composition.

The benthic composition of MCEs of the Cayman Islands has been previously described by Slattery & Lesser (2019), who focused primarily on sponge communities over limited spatial distribution and document a “sponge belt” (Slattery & Lesser, 2019). Mesophotic sponges in Grand Cayman were also the focus of Macartney et al. (2021), whose work was concentrated on a single species. While Carpenter et al. (2022) assessed the benthic composition of corals and key functional groups on MCEs across Little Cayman Island, their study did not include Grand Cayman. Importantly, none of the previous mesophotic studies in the Cayman Islands included assessments of fish communities, nor the interaction of benthic habitat and fish assemblages, highlighting a critical knowledge gap in our understanding of these ecosystems in Grand Cayman.

Variations in benthic habitat often corresponds to changes in fish community composition. For example, in the Seychelles shallow reefs with high macroalgal cover had higher densities of generalist fish species and lacked specialists, whereas reefs with high coral cover had higher abundances of specialised corallivores species (Chong-Seng et al., 2012). In the Philippines, a positive association was found between hard coral cover and biomass of planktivores, piscivores, and corallivores, while a negative association was found between hard coral cover and detritivores and sand feeders (Russ et al., 2021). In the Caribbean, a study found a positive correlation between macroalgae density and herbivorous fish assemblages on shallow reefs but fail to show an impact of coral cover on fish populations (Sandin, Sampayo & Vermeij, 2008). Mesophotic studies in the Gulf of Mexico have also shown a difference of fish assemblages in relation with benthic composition, with high fish biomass associated with coralline algal reefs (Voss et al., 2014). As these patterns vary across locations and depths, it is important to include assessments of both benthic habitat and fish communities to accurately characterize MCEs, which was previously lacking for MCEs in Grand Cayman.

Here, we characterize the biodiversity and community composition of mesophotic reefs across four sites on Grand Cayman Island. Using benthic photogrammetry and in situ visual surveys we provide a description of the benthic community composition, including scleractinian corals and key functional groups, as well as of fish community assemblages. Given the hypothesized role of MCEs as areas of refuge for preserving biodiversity under future climate change scenarios (e.g., DRRH; Bongaerts et al., 2010), baseline characterizations of these habitats are crucial for increasing our understanding of mesophotic ecosystem function and informing marine spatial planning and management.

Materials and Methods

The present study was carried out off the coast of Grand Cayman (19.329858–81.252361), the biggest and most populated of the three islands that compose the Cayman Islands (Economics and Statistics Office Government of the Cayman Islands, 2022). Four sites situated around the West Bay peninsula were surveyed: Ghost Mountain (19.402300–81.385680), Roundabout (19.384333–81.318300), La Mesa (19.321680–81.393080), and Lighthouse Point (19.372954–81.421855) (Fig. 1). The sites were situated off the wall, with a steep slope except for lighthouse point that is characterized by a more gradual slope (G. Goodbody-Gringley, 2022, personal observation). All four sites were composed of alternating hard-bottom and sandy substrates.

10.7717/peerj.17763/fig-1 Figure 1 Context of the study and sampling techniques.

Maps of the study sites, Lighthouse Point, Ghost Mountain, La Mesa and Roundabout and the position of Grand Cayman amidst the Caribbean Sea (A). Picture of the dive team laying transect and taking pictures of the quadrat during a 15 m practice dive (B). Site map generated with the “rnaturalearthdata” package (https://cran.r-project.org/package=rnaturalearth). Photo credit: Alex Chequer.

Fish surveys were conducted along 30 m by 2 m belt transects at a depth of 45 m (n = 3–5 per site; 60 m2 per transect). All fishes encountered along the transects were identified to species level, counted, and categorised into size classes based on total length: 0–5, 5–10, 10–20, 20–30, 30–40, and >40 cm. Fish surveys were conducted exclusively by divers using closed-circuit rebreathers (Hollis Prism2) to limit any influence of open vs. closed circuit diving on fish behavior. Following the surveys, fish data were classified into one of five food web classes: herbivore, invertivore, omnivore, macrocarnivore, and planktivore based on data available on www.fishbase.org (Welcomme, 1988; Houde & Zastrow, 1993), a data base by Froese & Pauly (2000). Density was calculated as the total number of individuals per transect (no. of fish/60 m2 transect). Biomass was calculated using the formula:

(a×Lb)N

where L is the estimated length using the upper end of the size class category, N is the number of individuals of that size class (Bohnsack & Harper, 1988; Kulbicki et al., 1993), parameter a is a scaling constant varying between species and parameter b is an allometric scaling factor. When b = 1, it indicates isometry and biomass scales linearly with size, when b < 1, it indicates negative allometry, and biomass increase slower than size and when b > 1, it is positive allometry and biomass increases faster than size. Values for a and b were obtained for each species on www.fishbase.org (Welcomme, 1988, Houde & Zastrow, 1993). When these values were not available for a particular species, data for a congeneric species were used. Lastly, species richness is the total number of species per transect and diversity is based on the Shannon index of diversity (Shannon & Weaver, 1949): −∑Pi(logPi)

with Pi being the proportion of each species in the sample.

At each of the four sites, a 10 m × 10 m plot was laid on the benthos using transect tapes, with the lower edge at 45 m depths. Images of the benthos were taken every second using the time lapse setting by two divers that swam in a lawnmower pattern approximately 1 m above the benthos, using two GoPro Hero 10 cameras mounted onto a PVC frame with two lights (Sola Video 2500 Flood). Eight targets were laid haphazardly on the substrate within the plot to assist with building a digital 3D model of the studied area using the software AGIsoft Metashape Pro (version 1.8.3). This was then converted into a 2D orthomosaic and exported as a PNG that was subsequently split into ~1 m2 replicate quadrats (48 in Ghost Mountain, 90 in La Mesa, 99 in Lighthouse Point and 100 in Roundabout) using Photoshop (version 24.3) and scaled on the pictured measuring tape. ImageJ (version 1.53t) software was then used for further analysis on all 337 images. Within each quadrat all scleractinian corals were identified to the highest taxonomic level possible based on image quality. Species within the genera Agaricia, Siderastrea, and Orbicella, were limited to genus. Surface area of each colony was then measured using the freehand tool (in cm2). Total surface area of all corals relative to the size of the quadrat was then used to calculate percent coral cover.

Benthic composition of other key functional groups was then assessed visually to the nearest 1% cover for eight categories: Crustose Coralline Algae (CCA), Soft Coral (SC), Sand and Rubble (S/R), Turf Algae (TA), Recently Dead Coral (RDC), Macroalgae (MA), and Other (soft coral plumes, fish, targets, blurry and indistinguishable area, etc.). Measurements and visual assessment of the functional groups were conducted by the same person throughout the data collection phase to limit bias and errors. Statistical analyses were conducted on R version 4.1.0 (R Core Team, 2023).

Fish density, biomass, diversity, coral diversity index and richness were compared among sites and among food webs with non-parametric Kruskal-Wallis’ tests (Kruskal & Wallis, 1952), as they did not meet the assumption of normality (Shapiro & Wilk, 1965), followed by post-hoc analysis using the Dunn’s test (Dunn, 1964) in the FSA package (Ogle et al., 2023). The fish diversity index met the assumption of normality (Shapiro & Wilk, 1965) and was compared among sites with a one-way analysis of variance (ANOVA), followed by a post hoc analysis using Tukey’s test (Tukey, 1949).

Benthic functional group composition among sites did not meet the assumptions of normality (Shapiro & Wilk, 1965) and was therefore compared using a Kruskal Wallis with a Dunn’s test as a post hoc. “Other”, while being measured as a functional group, is not integrated into any of the mentioned statistical analyses. By being an imprecise category that is compensating for a various array of non-benthic category measures or photomosaic inconsistence, it was judged irrelevant to the aim of the study and the measured percentages have been scaled accordingly.

To assess dissimilarity in the benthic composition and fish community composition among sites, we used non-Metric Multi-Dimensional scaling (nMDS) models from the Vegan package (Oksanen et al., 2022). Models were based on a Bray-Curtis dissimilarity matrix of community composition for both benthic composition, fish species composition, and reef fish trophic community. A Bray-Curtis dissimilarity is less sensitive to the absence of species within the matrix compared to the distance (metric) based counterparts. Thus, a dissimilarity matrix is preferable for our dataset which reflects community composition rather than a distance matrix that explains the distance between points without assessing the difference in community composition (Jaccard, 1901). Subsequently, to statistically compare the dissimilarity between sites for community composition, we used a PERMANOVA test where we set the dissimilarity matrix as the response, and site as the predicting factor.

Results

In total 48 fish species were recorded across all four sites ranging from 30 at Lighthouse Point to 24 at Roundabout, however, median species richness per transects did not differ significantly among sites (chi-squared = 5.803, P-value = 0.1216; Kruskal Wallis) (Fig. 2A). Yet, according to the species accumulation curve, the number of surveys conducted was not enough to account for the full species composition, and thus this is likely an underestimate of diversity (Fig. S1). Fish biomass also did not differ significantly among sites (chi-squared = 4.6559, P = 0.1988; Kruskal Wallis) (Fig. 2B), Lighthouse Point had the highest median biomass (5,936.296 ± 1,026.23 g/transect), while Ghost Mountain had the lowest (3,211.59 ± 797.70 g/transect) (Fig. 2B). The median fish density was significantly different across sites (chi-squared = 9.1522, df = 3, P-value = 0.02733; Kruskal Wallis). Among sites, the median fish density only presented significant differences (X2 = −2.74, P = 0.03; Dunn’s Test) between Lighthouse Point (157 ± 30.39 fish per transect) and La Mesa (63 ± 10.08 fish per transect) (Fig. 2C). Shannon diversity index was significantly different among sites (F-value = 4.655, P-value = 0.022; ANOVA) between Ghost Mountain and La Mesa (X2= 0.726, P = 0.019; Tukey), being respectively the least diverse with a median index value of 1.65 ± 0.14 and the most diverse with a median value of 2.31 ± 0.15 (Fig. 2D).

10.7717/peerj.17763/fig-2 Figure 2 Comparison of species richness (number of unique species per transect) (A), biomass (estimated weight of fish in g per transect) (B), density (number of fish per transect) (C) and diversity (Shannon diversity index) (D) between sites.

Boxplots’ horizontal bars represent the sites median values, the upper and lower sections (box outlines) represent the interquartile range, the whiskers extend to 90% of the data, with dots being outliers. Biomass is calculated using the formula (a * Lb)/N where a is a scaling constant different for each species, b is an allometric scaling factor that is different for each species, L is the length of the individual in cm and N is the number of individuals of the species. The diversity value is obtained via the Shannon-Weaver index of diversity (Shannon & Weaver, 1949). Based on Kruskal Wallis’ test followed by a Dunn’s test as post hoc, median species richness and biomass among sites did not differ significantly (Kruskal & Wallis, 1952). Median fish density between Lighthouse Point and La Mesa was significantly different (X2 = −2.74, P = 0.03), as well as the Shannon-Weaver diversity index between Ghost Mountain and La Mesa (X2 = 0.726, P = 0.019). *Significant difference (P < 0.05).

For trophic guilds specifically, invertivores were the most species rich group of fish with a median of 10 ± 1.31 species recorded and planktivores were the least with a median species richness of 3 ± 0.4 (Fig. 3A). Macrocarnivores had the highest median biomass (6,693.41 ± 1,434.09 g/transect; Fig. 3B) while planktivores had the lowest median biomass (789.6238 ± 394.81 g/transect). Omnivores had the highest median density (200 ± 4.42 of fish per transect; Fig. 3C). Significant differences were found in species richness (X2 = 12.169, P = 0.016; Kruskal Wallis) (Fig. 3A), biomass (X2 = 9.622, P = 0.047; Kruskal Wallis) (Fig. 3B) and fish density (X2 = 17.031, P = 0.001; Kruskal Wallis) (Fig. 3C) among trophic guilds. Higher species richness was observed for invertivores compared to planktivores (X2 = 3.13, P = 0.017; Dunn’s Test), while macrocarnivores had significantly lower density compared to planktivores (X2 = −2.816, P = 0.043; Dunn’s Test) and omnivores (X2 = −3.697, P = 0.002; Dunn’s Test). Shannon diversity index values did not differ significantly (F-value = 1.236, P = 0.327, ANOVA) (Fig. 3D) among trophic groups.

10.7717/peerj.17763/fig-3 Figure 3 Comparison of species richness (number of unique species per transect) (A), biomass (estimated weight of fish in mg per transect) (B), density (number of fish per transect) (C), and diversity (Shannon diversity index) (D) between trophic guilds.

Boxplots’ horizontal bars represent the sites median values, the upper and lower sections (box outlines) represent the interquartile range, the whiskers extend to 90% of the data, with dots being outliers. Biomass is calculated using the formula (a * Lb)/N where a is a scaling constant different for each species, b is an allometric scaling factor that is different for each species, L is the length of the individual in cm and N is the number of individuals of the species. The diversity value is obtained via the Shannon-Weaver index of diversity (Shannon & Weaver, 1949). Based on Kruskal Wallis’ test followed by a Dunn's test as post hoc, species richness (X2 = 12.169, P = 0.016), fish density (X2 = 17.031, P = 0.001) and biomass (X2 = 9.622, P = 0.047) among trophic guilds were significantly different. However, the Shannon diversity index did not differ significantly when comparing trophic groups. Significant difference *P < 0.05, **P < 0.005.

Across sites, 20 families of fish were recorded (Fig. 4). The basslet family (Grammatidae), had the highest density per square meter with a mean value of 19.048 ± 7.554 % (Fig. 4A), but only accounted for 0.59 ± 0.317 % of the total fish biomass (Fig. 4B). Members of the snapper family (Lutjanidae) had the highest average biomass (10.432 ± 3.974 %) but only accounted for 0.071 ± 0.031% of the recorded density (Figs. 4A and 4B).

10.7717/peerj.17763/fig-4 Figure 4 Biomass (A) and density (B) of the common fish families across all sites. Biomass (C) and density (D) of each species across all sites.

Solid lines in plots (A and B) are the average per common family, and the whiskers are the range of the data. Each rectangle in (C and D) shows at least one occurrence of a species on one of the transects. The absence of a rectangle represents the absence of a species on a site.

Across the surveyed sites, Gramma melacara, had the highest frequency by density at three sites, with 53.78% at Ghost Mountain, 35.91% at Roundabout and 24.75% at Lighthouse Point, while Chromis cyanea had the highest frequency at La Mesa of 20.61% (Fig. 4C). Scarus taeniopterus had the highest frequency by biomass at Ghost Mountain at 28.02%, while Lutjanus analis had the highest frequency by biomass at La Mesa at 41.85%, and Roundabout at 25.46% (Fig. 4D). At Lighthouse Point, Lutjanus mahogoni had the highest frequency by biomass at 20.91% (Fig. 4D). Out of the 48 observed species, only 13 were found across all sites and 22 species were endemic to one site.

Macroalgae percent cover was significantly higher than all other functional groups across all four sites (Fig. 5 and Table 1) with a median surface coverage of 63.378 ± 1.125%, while hard coral cover was significantly lower than all other functional groups, with a median of 0.662 ± 0.285% (Fig. 5 and Table 1). Sponge cover did not differ significantly from either sand/rubble or CCA cover, however, sand/rubble cover was significantly higher than CCA (Fig. 5 and Table 1).

10.7717/peerj.17763/fig-5 Figure 5 Comparison of the percentage cover of the different functional groups (macroalgae, sand/rubbles, sponge, CCA and hard coral).

The horizontal middle bars represent the median values, whilst the upper and lower sections (box outlines) represent the interquartile range. The whiskers extend to 90% of the data, with dots being outliers. Each functional group was compared using a Kruskal Wallis test followed by a Dunn’s test, as the values did not meet the assumptions of normality.

10.7717/peerj.17763/table-1 Table 1 Benthic component’s adjusted P-values throughout sites.

	Sand/rubble	Sponge	CCA	Hard coral	
Macroalgae	<0.001	<0.001	<0.001	<0.001	
Sand/rubble		0.067	0.001	<0.001	
Sponge			0.157	<0.001	
CCA				<0.001	

Agaricia was the dominant coral genus at all four sites (Fig. 6) ranging from 77% at Ghost Mountain to 54% at Lighthouse Point. Siderastrea and Orbicella were present in moderate frequency at all sites, Montastraea and Porites were consistently present but in low coverage, and Mycetophyllia was only recorded in La Mesa. Shannon diversity of coral species showed no significant differences among sites (P > 0.05, Kruskal-Wallis).

10.7717/peerj.17763/fig-6 Figure 6 Comparison of the percent contribution of each Scleracctinian coral genus observed per site.

With Agaricia being the most common genus recorded across all sites representing up to 70% of the coral cover at Ghost Mountain.

Multivariate analyses revealed the benthic composition was significantly dissimilar among sites (PERMANOVA, F value = 11.69, P = 0.0001). Lighthouse Point appears to be the most unique site based on benthic composition (Fig. 7A), related to high sand/rubble coverage (Fig. 7B). The trophic composition of fish communities (Fig. 7C) was also significantly different among sites (PERMANOVA, F value = 2.76, P = 0.0058) as was the species composition (Fig. 7E; PERMANOVA, F value = 2.55, P = 0.003).

10.7717/peerj.17763/fig-7 Figure 7 Ordination plots of the dissimilarity between the sites regarding benthic composition (A), fish trophic guild composition (C) and fish species (E).

Each dot represents a quadrat for the benthic ordination plot (A) and a transect for the fish trophic plots (C and E). Roundabout’s ellipse is not present owing to a deficit of replicates which does not allow the projection (C and E). Each ellipse represents the tendency of a site to drift along the two axes following the vectors on the right of each ordination plot (B and D). When compared via a multivariate analysis, it reveals a significant difference of benthic composition among sites. It also showed a significant difference between fish communities (PERMANOVA, F value = 2.76, P = 0.0058) and species composition among sites (PERMANOVA, F value = 2.55, P = 0.003). Significant difference *P < 0.05, **P < 0.005, ***P < 0.0005.

Discussion

Here, we provide the first characterization of mesophotic reefs in the Cayman Islands that combines benthic habitat and fish assemblages. Our results show that fish populations are primarily composed of basslet (Grammidae) and damselfish (Pomacentridae > Stegastinae) species, with species of wrasse (Labridae) and chromis (Pomacentridae > Chrominae) also present in high relative abundance. Although present in lower relative densities, species of snapper (Lutjanidae) and parrotfish (Scaridae) represent the highest contribution to community biomass on the mesophotic reefs surveyed (Fig. 4). The benthos at 45 m is dominated by macroalgae and sand/rubble, which can be generally described as algae-covered hard substrate dispersed among patches of sand and rubble. Sponges, seafans, and crustose coralline algae are also found in high frequency. However, hard coral cover represents the lowest percent cover (Fig. 5). Among the corals present, the dominant contributing species belong to the genera Agaricia, Siderastrea, and Orbicella (Fig. 6).

Upper mesophotic reefs often represent transition zones between shallow water ecosystems and lower mesophotic ecosystems (Pyle & Copus, 2019). Fish populations in these regions frequently harbor both depth-generalist species as well as distinctive mesophotic fauna (Baldwin, Tornabene & Robertson, 2018; Bejarano, Appeldoorn & Nemeth, 2014; Pinheiro et al., 2016). In fact, many of the species found on the MCEs surveyed here overlap with those documented on neighboring shallow reef sites in Grand Cayman (Johnson, Chequer & Goodbody-Gringley, 2023). Yet, among all the mesophotic sites surveyed, Gramma melacara (blackcap basslet) and Chromis insolata (sunshinefish) are the most abundant in terms of number of fish, species that are absent on nearby shallow reefs (Johnson, Chequer & Goodbody-Gringley, 2023). Previous surveys of Caribbean fish communities also note high densities of C. insolata on mesophotic reefs with their presence generally lacking on shallow reefs, suggesting this species is a mesophotic specialist (Pinheiro et al., 2016; Baldwin, Tornabene & Robertson, 2018; García-Hernández et al., 2018; Gress et al., 2018). Likewise, G. melacara is also a mesophotic specialist and is considered common on MCEs in the Western/Central Caribbean (Dustan & Lang, 2019; Chasqui et al., 2020), however it is not present on MCE species lists from Cozumel (Gress et al., 2018), Curacao (Pinheiro et al., 2016; Baldwin, Tornabene & Robertson, 2018), or Bermuda (Pinheiro et al., 2016; Goodbody-Gringley et al., 2019), and thus its geographic distribution may be limited to the Central Caribbean. Interestingly, despite being abundant in other locations at this depth (Andradi-Brown et al., 2017; Goodbody-Gringley et al., 2019, 2023) very few invasive lionfish (Pterois volitans) were found in this study (two total), suggesting they are currently not a major component of the mesophotic fish populations in Grand Cayman. Overall, the upper mesophotic reefs of Grand Cayman can clearly be categorized as transition ecosystems as they house both depth-generalist and mesophotic specialist species, however, species richness, diversity and abundance were all lower on MCEs (this study) compared to shallow reefs (Johnson, Chequer & Goodbody-Gringley, 2023).

Although species of snapper (Lutjanidae) are the most abundant by biomass across all sites, this is not consistent among the four sites surveyed, with these species being generally absent from Ghost Mountain. Rather, parrotfish are the dominant group at this site, where the most abundant species by biomass is Scarus taenopterus. Within the other three sites, variation also exists among the dominant species of snapper with Lutjanus analis being the highest contributor to biomass at La Mesa and Roundabout, while Lutjanus mahogoni, is the highest contributor at Lighthouse Point. Macrocarnivores, such as snappers, are often transient. So, it is likely that the variations in species contributions of these groups and the absence of snappers from surveys at Ghost Mountain (Fig. 4) may be an artifact of the visual survey method, which is limited to a short period of time and may miss transient species. Fish density and richness was found to differ significantly among sites, being higher in Lighthouse Point compared to the other sites. Lighthouse Point is situated at the northwest corner of Grand Cayman and is often characterized by converging currents that lead to strong mixing and potential upwelling. Increased mixing promotes nutrient availability and fuels the food chain, potentially leading to the significantly higher density and more speciose fish population found at this site. Increased plankton availability due to mixing is supported by the significantly higher density of planktivores compared to other trophic guilds found at Lighthouse Point (Bodungen et al., 1982; Rosa et al., 2015). Biomass, however, does not differ among sites, indicating that although density is higher at Lighthouse Point the average fish size is smaller. Smaller fish at this site may be related to the unique benthic composition of Lighthouse Point, which was dominated by areas of sand/rubble with low rugosity, limiting habitat availability for larger-bodied fishes (Kuffner et al., 2007).

Benthic cover of MCEs in Grand Cayman was found to be similar to mesophotic reefs at other Central Caribbean locations (Carpenter et al., 2022; Dustan & Lang, 2019; Reed et al., 2019). For example, sponges are the highest contributing benthic invertebrate to percent cover on MCEs in Grand Cayman, constituting an average of 10.587 ± 0.560% of the benthos, which conforms to previous reports from Lesser, Slattery & Mobley (2018) who described an increase in sponge abundance and diversity with increasing depth. However, while macroalgae cover in Grand Cayman is similar to that on MCEs at Pulley Ridge (Florida, USA) and in Puerto Rico (Reed et al., 2019; Appeldoorn et al., 2019), it is higher than that reported for several other Caribbean MCEs (Dustan & Lang, 2019; Gress et al., 2019), including Little Cayman (Carpenter et al., 2022). Differences in macroalgae cover between these neighbouring islands may be related to differences in levels of marine protection, as roughly 57% of the nearshore environment in Little Cayman is classified as no-take marine protected areas compared to less than 20% in Grand Cayman. While comparable fish data for mesophotic reefs in Little Cayman is not available, protection of ecologically important reef fish populations, such as herbivores, is expected to lead to reduced algal cover (Jackson et al., 2014).

While coral cover at 45 m depth in Grand Cayman is low, at an average of 2.6%, similar percent coral cover is noted on MCEs in Honduras, Mexico, and at Pulley Ridge. Yet, it is lower than the 8% cover documented in Bonaire and substantially lower than in Puerto Rico and the US Virgin Islands, where percent cover reaches nearly 30% on upper mesophotic sites (Appeldoorn et al., 2019; Frade et al., 2019; Gress et al., 2019; Reed et al., 2019; Smith, Holstein & Ennis, 2019). Despite regional differences in percent cover of scleractinian corals, species presence on MCEs appears relatively homogenous across the Caribbean, with species in the genus Agaricia being predominant (Reed et al., 2019; Gress et al., 2019; Appeldoorn et al., 2019; Smith, Holstein & Ennis, 2019; Frade et al., 2019). In Grand Cayman, coral species assemblages closely resemble those of nearby Little Cayman, where species in the genera Agaricia, Siderastraea and Orbicella have the highest contribution to percent cover at 45 m depth (Carpenter et al., 2022). However, Madracis spp. occur in relatively high frequency in Little Cayman yet are absent from the surveys in Grand Cayman. Conversely, the genus Mycetophylia is a top contributor to MCEs in Grand Cayman but not in Little Cayman. Variations in percent coral cover and species contributions among locations are likely driven by the availability of incident light as a function of reef topography. For example, on the Orbicella reefs in the US Virgin Islands, the reef slope is gradual with extensive flat regions in the upper mesophotic zone (Smith, Holstein & Ennis, 2019). However, the Cayman Islands are surrounded by deep trenches, with near vertical walls beyond the shallow fore reef leading to limited light availability beyond 30 m (Lesser, Slattery & Mobley, 2018; Slattery & Lesser, 2019).

Alternatively, differences in benthic composition, such as the dominance of Agaricia spp on Grand Cayman MCEs, may be related to disturbance history. The proximity of these MCEs to the Island of Grand Cayman exposes these sites to heavy levels of tourism and overfishing, impacting both coral and fish communities (Stallings, 2009; Johnson, Chequer & Goodbody-Gringley, 2023). In fact, benthic composition of MCEs is documented to affect the fish community, with high coral cover being linked to high fish density and richness (Garcia-Sais, 2010). Thus, MCEs likely do not provide a refuge from disturbance for coral and fish communities in regions with high disturbance frequency or severity, and will likely be composed of weedy, fast-growing coral species that can tolerate high disturbance conditions, associated with less dense and less diverse fish populations.

Conclusion

In conclusion, this study provides an overview of benthic and fish assemblages of Grand Cayman’s mesophotic reefs at 45 m, which qualitatively appear different to other regions within the Caribbean (Smith et al., 2010; García-Hernández et al., 2018; Scott et al., 2019). Our study shows that benthic composition is highly variable across a small geographical scale, with the exception of macroalgae, being the dominant benthic component across all sites. Fish assemblages were also distinct, with representatives of both depth-generalist and mesophotic-specific species present at all sites, but differences in biomass and density of trophic guilds among sites. Whether the unique composition of benthic and fish communities is a result of geological and biological processes, or the influence of local scale anthropogenic activity remains to be seen. Future research on remote mesophotic reefs within the same ecoregion will provide further insight into why Grand Caymans mesophotic reefs appear as distinct communities and will allow assessments of the relative influence from local anthropogenic activity.

Supplemental Information

10.7717/peerj.17763/supp-1 Supplemental Information 1 Species, Biomass and number of fish.

The dataset and code used to create Figures 2, 3, 4 and a part of Figure 7 as well as for the data analysis.

10.7717/peerj.17763/supp-2 Supplemental Information 2 Functional group percent cover throughout sites.

The dataset and the codes used to create Figures 5 and 7 and Table 1 and for the data analysis.

10.7717/peerj.17763/supp-3 Supplemental Information 3 Species accumulation curve comparing the cumulative count of observed species over the total number of surveys.

The curve did not reach a plateau after 16 surveys, showing the need for more replicates of the study. Light purple zone is the confidence interval around the mean with a level of confidence of 95%.

10.7717/peerj.17763/supp-4 Supplemental Information 4 Map of Grand Cayman and the Carribean.

Code used to create Figure 1.

10.7717/peerj.17763/supp-5 Supplemental Information 5 Size and genus of all the recorded coral colonies throughout sites.

The dataset and the codes used to create Figure 6 and for the data analysis.

We would like to express our gratitude to all those who contributed to this article, including researchers, interns, and support staff. We also extend our sincere thanks to the reviewers for their time, insightful comments, and valuable feedback, which have significantly improved the quality of this manuscript.

Additional Information and Declarations

Competing Interests

Author Contributions

Field Study Permissions

Data Availability

The authors declare that they have no competing interests.

Lucas Le Gall analyzed the data, prepared figures and/or tables, authored or reviewed drafts of the article, and approved the final draft.

Jack V. Johnson analyzed the data, prepared figures and/or tables, authored or reviewed drafts of the article, and approved the final draft.

Alex Chequer conceived and designed the experiments, performed the experiments, authored or reviewed drafts of the article, and approved the final draft.

Matthew Louis Doherty performed the experiments, authored or reviewed drafts of the article, and approved the final draft.

Gretchen Goodbody-Gringley conceived and designed the experiments, performed the experiments, authored or reviewed drafts of the article, and approved the final draft.

The following information was supplied relating to field study approvals (i.e., approving body and any reference numbers):

We did not need permits for this work.

The following information was supplied regarding data availability:

The raw data is available in the Supplemental Files.
==== Refs
References

Andradi-Brown et al. (2017) Andradi-Brown DA Vermeij MJA Slattery M Lesser M Bejarano I Appeldoorn R Goodbody-Gringley G Chequer AD Pitt JM Eddy C Smith SR Brokovich E Pinheiro HT Jessup E Shepherd B Rocha LA Curtis-Quick J Eyal G Noyes T Rogers AD Exton DA Large-scale invasion of western Atlantic mesophotic reefs by lionfish potentially undermines culling-based management Biological Invasions 2017 19 3 939 954 10.1007/s10530-016-1358-0
Appeldoorn et al. (2019) Appeldoorn R Alfaro M Ballantine D Bejarano I Ruiz H Schizas N Schmidt W Sherman C Weil E Puerto rico Mesophotic Coral Ecosystems. Coral Reefs of the World 2019 Cham Springer 111 129 10.1007/978-3-319-92735-0_7
Baldwin, Tornabene & Robertson (2018) Baldwin CC Tornabene L Robertson DR Below the mesophotic Scientific Reports 2018 8 1 4920 10.1038/s41598-018-23067-1 29559694
Bejarano, Appeldoorn & Nemeth (2014) Bejarano I Appeldoorn RS Nemeth M Fishes associated with mesophotic coral ecosystems in La Parguera, Puerto Rico Coral Reefs 2014 33 2 313 328 10.1007/s00338-014-1125-6
Bodungen et al. (1982) Bodungen BV Jickells TD Smith SR Ward JAD Hillier GB The Bermuda Marine Environment—Volume III. The final report of the Bermuda Inshore Waters Investigations 1975–1980 BBSR Special Publication 1982 18 123
Bohnsack & Harper (1988) Bohnsack JA Harper DE Length-weight relationships of selected marine reef fishes from the southeastern United States and the Caribbean 1988 https://repository.library.noaa.gov/view/noaa/3027 9 October 2023
Bongaerts et al. (2010) Bongaerts P Ridgway T Sampayo EM Hoegh-Guldberg O Assessing the ‘deep reef refugia’ hypothesis: focus on Caribbean reefs Coral Reefs 2010 29 2 309 327 10.1007/s00338-009-0581-x
Carpenter et al. (2022) Carpenter GE Chequer AD Weber S Mass T Goodbody-Gringley G Light and photoacclimatization drive distinct differences between shallow and mesophotic coral communities Ecosphere 2022 13 8 e4200 10.1002/ecs2.4200
Chasqui et al. (2020) Chasqui L Rincón-Díaz N Vanegas MJ Chasqui L Rincón-Díaz N Vanegas MJ Abundance of the invasive lionfish Pterois volitans in the coastal coral reefs from Colombian Caribbean Boletín de Investigaciones Marinas y Costeras—INVEMAR 2020 49 157 170 10.25268/bimc.invemar.2020.49.1.779
Chong-Seng et al. (2012) Chong-Seng KM Mannering TD Pratchett MS Bellwood DR Graham NAJ The influence of coral reef benthic condition on associated fish assemblages PLOS ONE 2012 7 8 e42167 10.1371/journal.pone.0042167 22870294
Dunn (1964) Dunn OJ Multiple comparisons using rank sums Technometrics 1964 6 241 252 10.2307/1266041
Dustan & Lang (2019) Dustan P Lang JC Loya Y Puglise KA Bridge TCL Discovery Bay, Jamaica Mesophotic Coral Ecosystems 2019 Cham Springer International Publishing 85 109
Economics and Statistics Office Government of the Cayman Islands (2022) Economics and Statistics Office Government of the Cayman Islands The Cayman Islands’ 2021 Census of Population and Housing Report 2022 https://www.eso.ky/2021-population-and-housing-census-report.html
Frade et al. (2019) Frade PR Bongaerts P Baldwin CC Trembanis AC Bak RPM Vermeij MJA Loya Y Puglise KA Bridge TCL Bonaire and Curaçao Mesophotic Coral Ecosystems 2019 Cham Springer International Publishing 149 162
Froese & Pauly (2000) Froese R Pauly D FishBase 2000: concepts designs and data sources 2000 Penang, Malaysia WorldFish
Garcia-Sais (2010) Garcia-Sais JR Reef habitats and associated sessile-benthic and fish assemblages across a euphotic-mesophotic depth gradient in Isla Desecheo, Puerto Rico Coral Reefs 2010 29 2 277 288 10.1007/s00338-009-0582-9
García-Hernández et al. (2018) García-Hernández JE Sanchez PJ Hammerman NM Schizas NV Fish, coral, and sponge assemblages associated with Altiphotic and mesophotic reefs along the Guánica biosphere reserve continental shelf edge, Southwest Puerto Rico Frontiers in Marine Science 2018 5 303 10.3389/fmars.2018.00303
García-Sais et al. (2008) García-Sais J Appeldoorn R Battista T Bauer L Bruckner A Caldow C Carrubba L Corredor J Diaz E Lilyestrom C García-Moliner G Hernández-Delgado E Menza C Morell J Pait A Sabater J Weil E Williams E Williams S The state of coral reef ecosystems of Puerto Rico 2008 https://www.researchgate.net/publication/265082327_Puerto_Rico_The_State_of_Coral_Reef_Ecosystems_of_Puerto_Rico_The_State_of_Coral_Reef_Ecosystems_of_Puerto_Rico
Goodbody-Gringley et al. (2023) Goodbody-Gringley G Chequer A Grincavitch C Noyes T Dowell R Lundberg A Corbett E Smith A Impacts of recurrent culling of invasive lionfish on mesophotic reefs in Bermuda Coral Reefs 2023 42 2 443 452 10.1007/s00338-023-02354-y
Goodbody-Gringley et al. (2019) Goodbody-Gringley G Eddy C Pitt JM Chequer AD Smith SR Ecological drivers of invasive lionfish (Pterois volitans and Pterois miles) distribution across mesophotic reefs in Bermuda Frontiers in Marine Science 2019 6 716 10.3389/fmars.2019.00258
Gress et al. (2018) Gress E Arroyo-Gerez MJ Wright G Andradi-Brown DA Assessing mesophotic coral ecosystems inside and outside a Caribbean marine protected area Royal Society Open Science 2018 5 10 180835 10.1098/rsos.180835 30473832
Gress et al. (2019) Gress E Voss JD Eckert RJ Rowlands G Andradi-Brown DA Loya Y Puglise KA Bridge TCL The mesoamerican reef Mesophotic Coral Ecosystems 2019 Cham Springer International Publishing 71 84
Hinderstein et al. (2010) Hinderstein LM Marr JCA Martinez FA Dowgiallo MJ Puglise KA Pyle RL Zawada DG Appeldoorn R Theme section on “mesophotic coral ecosystems: characterization, ecology, and management” Coral Reefs 2010 29 2 247 251 10.1007/s00338-010-0614-5
Holstein et al. (2019) Holstein DM Fletcher P Groves SH Smith TB Loya Y Puglise KA Bridge TCL Ecosystem services of mesophotic coral ecosystems and a call for better accounting Mesophotic Coral Ecosystems. Coral Reefs of the World 2019 Cham Springer International Publishing 943 956
Houde & Zastrow (1993) Houde ED Zastrow CE Ecosystem- and taxon-specific dynamic and energetics properties of larval fish assemblages Bulletin of Marine Science 1993 53 2 290 335
Jaccard (1901) Jaccard P Étude comparative de la distribution florale dans une portion des Alpes et des Jura Bulletin de la Société Vaudoise des Sciences Naturelles 1901 37 547 579 10.5169/seals-266450
Jackson et al. (2014) Jackson J Donovan M Cramer K Lam V Status and trends of Caribbean coral reefs: 1970–2012 2014 Gland, Switzerland Global Coral Reef Monitoring Network; International Union for the Conservation of Nature (IUCN)
Johnson, Chequer & Goodbody-Gringley (2023) Johnson JV Chequer AD Goodbody-Gringley G Insights from the 2-year-long human confinement experiment in Grand Cayman reveal the resilience of coral reef fish communities Scientific Reports 2023 13 21806 10.1038/s41598-023-49221-y 38071390
Kahng et al. (2010) Kahng SE Garcia-Sais JR Spalding HL Brokovich E Wagner D Weil E Hinderstein L Toonen RJ Community ecology of mesophotic coral reef ecosystems Coral Reefs 2010 29 2 255 275 10.1007/s00338-010-0593-6
Kruskal & Wallis (1952) Kruskal WH Wallis WA Use of ranks in one-criterion variance analysis Journal of the American Statistical Association 1952 47 260 583 621 10.1080/01621459.1952.10483441
Kuffner et al. (2007) Kuffner IB Brock JC Grober-Dunsmore R Bonito VE Hickey TD Wright CW Relationships between reef fish communities and remotely sensed rugosity measurements in Biscayne National Park, Florida, USA Environmental Biology of Fishes 2007 78 1 71 82 10.1007/s10641-006-9078-4
Kulbicki et al. (1993) Kulbicki M Moutham G Thollot P Wantiez L Length-weight relationships of fish from the lagoon of New Caledonia Naga Iclarm Q 1993 16 2/3 26 30
Laverick, Andradi-Brown & Rogers (2017) Laverick JH Andradi-Brown DA Rogers AD Using light-dependent scleractinia to define the upper boundary of mesophotic coral ecosystems on the reefs of Utila, Honduras PLOS ONE 2017 12 8 e0183075 10.1371/journal.pone.0183075 28809933
Laverick et al. (2018) Laverick JH Piango S Andradi-Brown DA Exton DA Bongaerts P Bridge TCL Lesser MP Pyle RL Slattery M Wagner D Rogers AD To what extent do mesophotic coral ecosystems and shallow reefs share species of conservation interest? A systematic review Environmental Evidence 2018 7 1 15 10.1186/s13750-018-0127-1
Lesser, Slattery & Leichter (2009) Lesser MP Slattery M Leichter JJ Ecology of mesophotic coral reefs Journal of Experimental Marine Biology and Ecology 2009 375 1–2 1 8 10.1016/j.jembe.2009.05.009
Lesser, Slattery & Mobley (2018) Lesser MP Slattery M Mobley CD Biodiversity and functional ecology of mesophotic coral reefs Annual Review of Ecology, Evolution, and Systematics 2018 49 1 49 71 10.1146/annurev-ecolsys-110617-062423
Macartney et al. (2021) Macartney KJ Clayshulte Abraham A Slattery M Lesser MP Growth and feeding in the sponge Agelas tubulata from shallow to mesophotic depths on Grand Cayman Island Ecosphere 2021 12 9 e03764 10.1002/ecs2.3764
Ogle et al. (2023) Ogle DH Doll JC Wheeler AP Dinno A FSA: simple fisheries stock assessment methods 2023 R package version 0.9.5 https://CRAN.R-project.org/package=FSA
Oksanen et al. (2022) Oksanen J Simpson GL Blanchet FG Kindt R Legendre P Minchin PR O’Hara RB Solymos P Stevens MHH Szoecs E Wagner H Barbour M Bedward M Bolker B Borcard D Carvalho G Chirico M Caceres MD Durand S Evangelista HBA FitzJohn R Friendly M Furneaux B Hannigan G Hill MO Lahti L McGlinn D Ouellette MH Cunha ER Smith T Stier A Braak CJFT Weedon J vegan: community ecology package 2022 https://CRAN.R-project.org/package=vegan
Pérez-Castro et al. (2022) Pérez-Castro MÁ Schubert N Ang-Montes de Oca G Leyte-Morales GE Eyal G Hinojosa-Arango G Mesophotic coral ecosystems in the eastern tropical pacific: the current state of knowledge and the spatial variability of their depth boundaries Science of the Total Environment 2022 806 150576 10.1016/j.scitotenv.2021.150576 34582873
Pinheiro et al. (2016) Pinheiro HT Goodbody-Gringley G Jessup ME Shepherd B Chequer AD Rocha LA Upper and lower mesophotic coral reef fish communities evaluated by underwater visual censuses in two Caribbean locations Coral Reefs 2016 35 1 139 151 10.1007/s00338-015-1381-0
Pyle & Copus (2019) Pyle RL Copus JM Loya Y Puglise KA Bridge TCL Mesophotic coral ecosystems: introduction and overview Mesophotic Coral Ecosystems 2019 Cham Springer International Publishing 3 27
R Core Team (2023) R Core Team R: a language and environment for statistical computing 2023 Vienna, Austria R Foundation for Statistical Computing https://www.R-project.org
Reed et al. (2019) Reed JK Farrington S David A Harter S Pomponi SA Cristina Diaz M Voss JD Spring KD Hine AC Kourafalou VH Smith RH Vaz AC Paris CB Dennis Hanisak M Loya Y Puglise KA Bridge TCL Pulley ridge, Gulf of Mexico, USA Mesophotic Coral Ecosystems 2019 Cham Springer International Publishing 57 69
Rosa et al. (2015) Rosa M Alves A Medeiros D Coni E Ferreira C Ferreira B Rosa R Amado-Filho G Pereira-Filho G Moura R Thompson F Sumida P Francini-Filho R Mesophotic reef fish assemblages of the remote St. Peter and St. Paul’s Archipelago, Mid-Atlantic Ridge, Brazil Coral Reefs 2015 35 113 123 10.1007/s00338-015-1368-x
Russ et al. (2021) Russ GR Rizzari JR Abesamis RA Alcala AC Coral cover a stronger driver of reef fish trophic biomass than fishing Ecological Applications 2021 31 e02224 10.1002/eap.2224 32866333
Sandin, Sampayo & Vermeij (2008) Sandin SA Sampayo EM Vermeij MJA Coral reef fish and benthic community structure of Bonaire and Curaçao, Netherlands Antilles Caribbean Journal of Science 2008 44 137 144 10.18475/cjos.v44i2.a2
Scott et al. (2019) Scott AR Battista TA Blum JE Noren LN Pawlik JR Patterns of benthic cover with depth on Caribbean mesophotic reefs Coral Reefs 2019 38 5 961 972 10.1007/s00338-019-01824-6
Shannon & Weaver (1949) Shannon C Weaver W The mathematical theory of communication 1949 Illinois Urban University of Illinois Press 125
Shapiro & Wilk (1965) Shapiro SS Wilk MB An analysis of variance test for normality (complete samples) Biometrika 1965 52 3/4 591 10.2307/2333709
Slattery & Lesser (2019) Slattery M Lesser MP Loya Y Puglise KA Bridge TCL The Bahamas and Cayman Islands Mesophotic Coral Ecosystems. Coral Reefs of the World 2019 Cham Springer International Publishing 47 56
Smith et al. (2010) Smith TB Blondeau J Nemeth RS Pittman SJ Calnan JM Kadison E Gass J Benthic structure and cryptic mortality in a Caribbean mesophotic coral reef bank system, the Hind Bank Marine Conservation District, U.S. Virgin Islands Coral Reefs 2010 29 2 289 308 10.1007/s00338-009-0575-8
Smith, Holstein & Ennis (2019) Smith TB Holstein DM Ennis RS Loya Y Puglise KA Bridge TCL Disturbance in mesophotic coral ecosystems and linkages to conservation and management Mesophotic Coral Ecosystems 2019 Cham Springer International Publishing 911 929
Stallings (2009) Stallings CD Fishery-independent data reveal negative effect of human population density on Caribbean predatory fish communities PLOS ONE 2009 4 5 e5333 10.1371/journal.pone.0005333 19421312
Stefanoudis et al. (2019) Stefanoudis PV Rivers M Smith SR Schneider CW Wagner D Ford H Rogers AD Woodall LC Low connectivity between shallow, mesophotic and rariphotic zone benthos Royal Society Open Science 2019 6 9 190958 10.1098/rsos.190958 31598316
Tukey (1949) Tukey JW Comparing individual means in the analysis of variance Biometrics 1949 5 2 99 10.2307/3001913 18151955
Voss et al. (2014) Voss J Williams M Reed J Clark R Chapter 5 benthic and fish communities in the mid and lower mesophotic zone of the sanctuary Fish and Benthic Communities of the Flower Garden Banks National Marine Sanctuary: Science to Support Sanctuary Management 2014 Silver Spring, MD NOAA Technical Memorandum NOS NCCOS 201 260
Welcomme (1988) Welcomme RL International introductions of inland aquatic species 1988 Rome Food & Agriculture Organization
