Quantifying Mesophotic Fish Assemblages of Hawai‘i’s Au‘au Channel: Associations With Benthic Habitats and Depth

Boland, Raymond C.; Hyrenbach, K. David; DeMartini, Edward E.; Parrish, Frank A.; Rooney, John J.

doi:10.3389/fmars.2021.785308

ORIGINAL RESEARCH article

Front. Mar. Sci., 05 January 2022

Sec. Coral Reef Research

Volume 8 - 2021 | https://doi.org/10.3389/fmars.2021.785308

Quantifying Mesophotic Fish Assemblages of Hawai‘i’s Au‘au Channel: Associations With Benthic Habitats and Depth

$\r\nRaymond C. Boland*$ Raymond C. Boland^1*

K. David Hyrenbach²

Edward E. DeMartini¹

Frank A. Parrish¹ $John J. Rooney&#x;\r\n$ John J. Rooney^1†

¹NOAA NMFS Pacific Islands Fisheries Science Center, Honolulu, HI, United States
²Department of Natural Sciences, Hawai’i Pacific University, Kaneohe, HI, United States

Mesophotic reefs (30–150 m) occur in the tropics and subtropics at depths beyond most scientific diving, thereby making conventional surveys challenging. Towed cameras, submersibles, and mixed-gas divers were used to survey the mesophotic reef fish assemblages and benthic substrates of the Au‘au Channel, between the Hawaiian Islands of Maui and Lāna‘i. Non-parametric multivariate analysis: Non-metric Multidimensional Scaling (NMDS), Hierarchical Cluster Analysis (HCA), Multi-Response Permutation Procedure (MRPP), and Indicator Species Analysis (ISA) were used to determine the association of mesophotic reef fish species with benthic substrates and depth. Between 53 and 115-m depths, 82 species and 10 genera of fish were observed together with 10 types of benthic substrate. Eight species of fish (Apolemichthys arcuatus, Centropyge potteri, Chaetodon kleinii, Chromis leucura, Chromis verater, Forcipiger sp., Naso hexacanthus, and Parupeneus multifasciatus) were positively associated with increasing depth, Leptoseris sp. coral cover, and hard-bottom cover, and one species (Oxycheilinus bimaculatus) of fish was positively associated with increasing Halimeda sp. algae cover. Fish assemblages associated with rubble were not significantly different from those associated with sand, Montipora coral beds and Leptoseris coral beds, but were distinct from fish assemblages associated with hard bottom. The patterns in the data suggested two depth assemblages, one “upper mesophotic” between 53 and 95 m and the other deeper, possibly part of a “lower mesophotic” assemblage between 96 and 115 m at the edge of the rariphotic and bottomfish complex.

Introduction

Mesophotic Coral Ecosystems (MCEs) are the most understudied coral reef ecosystems in the world. MCEs are found in the tropics and subtropics between 30 and 150-m depths, a range defined by the limit of scientific scuba diving (30 m) and corals’ ability to photosynthesize (150 m) (Hinderstein et al., 2010). The mapping of mesophotic benthic habitats has documented distinct coral and algal substrates (Rooney et al., 2010; Blythe-Skyrme et al., 2013), but little data exist for the fish assemblages that inhabit these substrates, particularly in Hawai‘i. MCEs have been described as transition zones between shallow and deep-reef fish assemblages (Brokovich et al., 2008; Weijerman et al., 2019), and represent a frontier for research as demonstrated by the continuing discovery of new and conspicuous reef fish species (Pyle, 1996; Pyle et al., 2008; Rocha et al., 2017; Arango et al., 2019). Despite increasing interest in MCEs during the last decade (Lindfield et al., 2016; Abesamis et al., 2020), the effective management of this habitat remains problematic due to limited knowledge about the community structure and habitat associations of mesophotic fish assemblages (Dennis and Bright, 1988; Hinderstein et al., 2010; Bejarano et al., 2014; Pinheiro et al., 2016).

In Hawai‘i, MCEs have recently been described within select geographic regions in association with specific bathymetric habitat features (e.g., banks, ledges) or benthic cover (e.g., algal meadows, high coral cover) (Rooney et al., 2010; Blythe-Skyrme et al., 2013; Pyle et al., 2016). Many studies have investigated various aspects of MCE benthic habitats in the main Hawaiian Islands including coral beds, algal meadows, coral physiology, fish life histories, and fish trophic ecology (Boland and Parrish, 2005; Kahng and Kelley, 2007; Rooney et al., 2010; Blythe-Skyrme et al., 2013; Costa et al., 2015; Pochon et al., 2015; Bradley et al., 2016; Pyle et al., 2016; Spalding et al., 2016; Asher et al., 2017; Weijerman et al., 2019), with the number and spatial extent of surveys increasing in recent years (2009–2018) (Spalding et al., 2019). Much less is known of MCEs in the Northwestern Hawaiian Islands (NWHI) as a result of the additional logistical constraints of research at these distant locations. Seven studies conducted in the NWHI have been published to date, focusing on benthic substrate cover, fish assemblages on banks, a comparison of shallow reef vs. mesophotic fish, foraging habitat of the Hawaiian Monk Seal and/or levels of fish endemism (Parrish et al., 2000; Parrish and Boland, 2004; Rooney et al., 2010; Kane et al., 2014; Fukunaga et al., 2016, 2017; Kosaki et al., 2016).

Studies in Hawai‘i and elsewhere have described zonation in fish assemblages across mesophotic depths. Weijerman et al. (2019) divided fish assemblages of the mesophotic zone in Hawai‘i into an upper mesophotic zone (30–59 m) and a lower mesophotic zone (60–129 m) with deeper depths described as “rariphotic” and “bathyal” zones.

Mesophotic fish assemblages have been associated with various benthic habitats (Dennis and Bright, 1988; Bejarano et al., 2010; Garcia-Sais, 2010; Asher et al., 2017; Stefanoudis et al., 2019; Williams et al., 2019) and with depth (Brokovich et al., 2008; Pearson and Stevens, 2015; Semmler et al., 2016; Asher et al., 2017; Sih et al., 2017; Coleman et al., 2018; Stefanoudis et al., 2019; Williams et al., 2019) and observed as early as 1981 in the Cayman Islands (Lukens, 1981). Thresher and Colin (1986) perceived multiple “fish communities” within the mesophotic assemblage of fishes on the outer wall at Enewetak Atoll in the Marshall Islands; however, their categorization was not rigorously quantitative. Dennis and Bright (1988) recognized multiple benthic-fish associations that were related to depth at the Flower Garden banks in the northwestern Gulf of Mexico.

This study investigates the MCE reef fish assemblages of the Au‘au Channel between Lāna‘i and Maui in the main Hawaiian Islands using archived MCE survey data collected between 2007 and 2011 by the National Oceanic and Atmospheric Administration (NOAA) National Marine Fisheries Service (NMFS) These surveys utilized three distinct platforms: Towed Optical Assessment Device (TOAD), submersibles Pisces IV and V (SUB), and mixed-gas scuba diving (SCUBA).

Materials and Methods

Study Site

Seafloor depth of Au‘au Channel (Figure 1) ranges from 0 to approximately 250 m, with much of the area located in the mesophotic zone (Costa et al., 2015). During the last glacial maximum 21,000 years ago, this channel was a dry land bridge between Maui and Lāna‘i. Lakes and lagoons on this land bridge were submerged by rising sea level and became ideal habitat for hard corals and precious corals in the genera Antipathes and Myripathes (Grigg et al., 2002; Price and Elliott-Fisk, 2004). The Au‘au Channel bottom is heterogeneous and features a variety of substrate types such as large beds of Halimeda sp. algae, sand, Montipora sp. coral and Microdictyon sp. algae, hard pavement, Leptoseris sp. coral, and rubble (Costa et al., 2015).

FIGURE 1

Figure 1. Bathymetry of the Au‘au Channel with survey areas. Lāna‘i Island is on the left and Maui Island is on the right. TOAD surveys denoted by purple, Submersible surveys in blue, Diver surveys in green and Submersible and Diver surveys in black.

Survey Methods

Visual surveys of fish and substrate cover were conducted using strip transects of various widths and distances, dependent on the survey platform. Fish were identified to the finest taxonomic level possible (ideally to species or minimally to genus) and their abundance recorded. All fish were sized (standard length, SL) to the nearest centimeter using platform-specific techniques. The TOAD and submersibles had two or more lasers spaced a specific distance apart to aid in length estimation during observations (Kelley et al., 2005; Rooney et al., 2010). Both the TOAD and submersible maintained video recordings of the survey with the laser scale within the field of view. These video recordings were matched against the observer’s observations to control or correct for bias. To minimize possible biases and variability in fish identification and sizing, one observer (the senior author) was primarily used, with two additional trained observers. Observations from the additional observers were compared to the video log and to the primary observer’s observations to control for bias.

Divers conducted transects that were 25 m long by 2 m wide (Brokovich et al., 2008). One diver counted and sized fish; the second diver followed and video-recorded the substrate along the transect tape.

Submersible transects followed a specific isobath and habitat and finished either after 20 min or when the end of the habitat was reached (Kelley et al., 2005). The position of the submersible was continuously recorded via TrackLink 5000HA (LinkQuest Inc.) and plotted on a bathymetric chart. Area surveyed was determined by the transect distance multiplied by a strip width of 6 m, which is the width of area of a single viewport on the submersible (Kelley et al., 2005). Due to the constraints of surveying fish out of a submersible viewport, only fish ≥10 cm were identified and sized. Observations were verbally recorded onto the mission videotape via microphone within the command sphere and later transcribed and scored.

TOAD transects followed an approximately straight track, regardless of depth or substrate. The TOAD was deployed at a specific geographic location, towed continuously, tracked using GPS, and recovered at transect end. Thus, transects varied in both duration and length, with the distance of the track determined by mapping the GPS locations and calculating the linear distances between all consecutive points. Transect area surveyed was determined using transect length multiplied by the camera’s view width of 22 m (Rooney et al., 2010). Due to the limitations of the camera only fish ≥10 cm could be identified and sized.

Transect substrate was recorded by video. TOAD and SUB transects used the point intersect method to determine percent cover. Every 30 s of video was scored by an evenly spaced 5-point screen overlay to calculate percent cover (Blythe-Skyrme et al., 2013). Because divers surveyed small specific areas limited to one type of substrate, percent cover for SCUBA surveys was visually estimated along the transect. These substrate observations were classified into 10 categories: sand, hard bottom, rubble, Leptoseris sp. Coral, Montipora sp. Coral, Porites sp. Coral, Black Coral, Halimeda sp. Algae, Microdictyon sp. Algae, and Ulva sp. Algae. Percent cover data for the 10 categories were estimated for every survey.

Data Analysis

Survey methods varied in transect strip width and length, fish size detectability, and spanned different depths and substrates, which produced unstandardized data that were not normally distributed and contained numerous zeros and some high counts. Data were transformed into presence/absence and non-parametric multivariate analyses were employed to counter the sampling limitations and data deficiencies.

Fish species were divided into two life stages (juvenile, adult) using 50% of the maximum published size (Randall, 2005, 2007; Froese and Pauly, 2014)¹ as the threshold value between juvenile and adult (Binohlan and Froese, 2009; Froese and Pauly, 2014; see text footnote 1). This produced a primary matrix of 180 species-life stages and 52 surveys.

Mean depth and mean percent cover by substrate type were estimated for each survey, resulting in a secondary matrix of 11 environmental variables (depth plus 10 substrates) and 52 surveys. Using the PC-Ord program, these matrices were used to perform four complementary multivariate analyses: Non-metric Multidimensional Scaling (NMDS), Hierarchical Cluster Analysis (HCA), Multi-response Permutation Procedure (MRPP), and Indicator Species Analysis (ISA) (McCune and Grace, 2002). To make the NMDS and the HCA results comparable, both relied on the relative Sorensen distance measure. The significance for association in NMDS was estimated by 999 randomized iterations and for ISA it was 4999 permutations (Legendre and Legendre, 1998; McCune and Grace, 2002).

Two Mantel tests (McCune and Grace, 2002) were conducted to assess the null hypothesis of no relationship between the similarity of the fish community (matrix 1) and the substrate data (matrix 2) for all pair-wise surveys. NMDS was chosen over Principal Component Analysis (PCA) to compare matrix 1 with matrix 2 because of the large number of zero values and occasional large counts in data matrices (Kenkel and Orloci, 1986; McCune and Grace, 2002). HCA was used to organize the surveys into discrete groups based on similarities among fish species and utilized the overall percent of variance explained by NMDS (50%) as the threshold for these groups. These groups also indicated the five dominant benthic habitats and was the basis for the ranges of the six depth bins. MRPP tested the null hypothesis of no difference between all pair-wise group comparisons of the five benthic habitat types, and six binned depth categories. Depth bin categories were defined by 10 m depth increments, except for the first bin which was a 15 m depth increment. ISA was used to identify which fish species were indicators of the distinct MRPP groups identified on the basis of depth categories, and benthic habitat types.

Results

Fifty-two surveys (9 SCUBA diver, 29 SUB and 14 TOAD) were conducted between 53 and 115 meters utilizing the three survey platforms (Table 1). Eighty-two fish species and an additional ten fish genera were identified (92 total taxa). Ten different benthic substrates were observed and classified as: sand, hard bottom, rubble, Leptoseris sp. coral, Montipora sp. coral, Porites sp. coral, black coral, Halimeda sp. algae, Microdictyon sp. algae, and Ulva sp. algae. There was a significant but weak, positive relationship between fish assemblage and depth or substrate (Mantel test, r = 0.200, P < 0.001, 999 randomizations).

TABLE 1

Table 1. Survey methods used in this study, showing effort and number of fish recorded.

NMDS yielded a significant (P < 0.001) 3-dimensional solution that altogether explained 50.4% of the variance (21.7% Axis 1: Depth, Leptoseris and Hard Bottom, 14.8% Axis 2: Montipora and 13.9% Axis 3: Halimeda) with a stress of 11.567. Using Clarke’s rule of thumb where 5–10 is good ordination and 10–20 is a usable portrayal, a stress of 11.567 indicates a low risk of misinterpretation (McCune and Grace, 2002; Figure 2). Pair-wise orthogonality of the 3 axes were 99.3% (axis 1 vs. axis 2), 99.3% (axis 1 vs. axis 3), and 99.7% (axis 2 vs. axis 3), indicating that all axes were aligned at nearly 90 degrees to each other. Mean depth, 2 substrate types, and 8 fish species were significantly correlated with Axis 1 (Table 2). Both age classes of Cirrhilabrus jordani but no substrates or depth, were significantly correlated with Axis 2. Finally, Halimeda substrate and Adult Oxycheilinus bimaculatus were significantly correlated with Axis 3 (Table 2).

FIGURE 2

Figure 2. Non-metric Multidimensional Scaling (NMDS) plot of fish species abundance and substrate categories. Samples are color-coded according to their hierarchical cluster membership. Environmental variables are shown using gradients in depth and four substrate types (Hard Bottom, Halimeda, Leptoseris, and Montipora).

TABLE 2

Table 2. Significant Kendall Tau correlations (P < 0.001) of environmental variables and fish species within the axes-space of Non-metric Multidimensional Scaling.

NMDS ordination was used to identify discrete groupings of samples with similar fish community structure, using Hierarchical Clustering with a threshold of 50% for defining groupings (Figure 3). This produced ten clusters and three ungrouped surveys (Table 3 and Figure 3). These clusters, which were grouped primarily by fish species and secondarily by substrate types had descriptors assigned based upon this hierarchy of definition (Table 3). These different groups were then plotted on the NMDS ordination graph (Figure 2) and illustrated the relation of the groups with the substrate-defined axes.

FIGURE 3

Figure 3. Hierarchical clustering analysis of fish abundance by survey. Red line indicates the 50% explained variance cut-off from the NMDS, red numbers indicate the 10 clusters and blue numbers indicate surveys that were not part of a cluster. Survey by method (T, TOAD; S, Submersible; DD, Deep Diver) and date.

TABLE 3

Table 3. Results of the hierarchical clustering analysis identifying and quantifying the components of clusters.

The grouping pattern appeared to correlate well with the plots of substrate (Figure 4) and fish species (Figure 5) indicating fish assemblages were affected by substrate. There are two clusters of 7 surveys each that are defined by a damsel fish species (Chromis verater or Chromis leucura) and each species association with Leptoseris coral beds (Table 3). Another cluster of six surveys defines shallow Montipora coral beds as habitats with juvenile Cirrhilabrus jordani. A cluster of six surveys shows the relationship of adult Oxychelinus bimaculatus with Halimeda algal beds. A cluster of five surveys defined a transition zone that had adult Cirrhilabrus jordani as the dominant fish species. Clusters defined transition zones, sand and rubble zones and even hard bottom zones (Table 3). The clusters are grouped consistently in regards to their member surveys when plotted (Figure 2) and suggests a strong positive association with substrate. These clusters defined five primary benthic habitats: Sand, Hard Bottom, Leptoseris coral beds, Montipora coral beds and Rubble. While Rubble was never a dominant substrate, it accounted for a large percentage of benthic habitat especially in the transition zone. Three singular surveys had unusual fish abundance that impacted their compatibility with other groups. The two surveys that were dominated by large schools were submersible surveys. In one case at a depth of 99 m, 85% of the observed fish were Aphraeus rutilans. In the other case at a depth 93 m, 95% of the observed fish were Luzonichthys earlei. A TOAD survey recorded a single Oxychelinus bimaculatus at a depth of 73 m (Table 3).

FIGURE 4

Figure 4. Non-metric Multidimensional Scaling plot of environmental variables, showing their Kendall Tau correlations with the three axes (vectors). Axis one is correlated with depth, Leptoseris and hard bottom. Axis two did not correlate with any substrate and axis three correlated with Halimeda.

FIGURE 5

Figure 5. Non-metric Multidimensional Scaling plot of fish species, showing their Kendall Tau correlations with the three axes (vectors). Axis one is correlated with depth, Leptoseris and hard bottom. Axis three is correlated with Halimeda. Axis two is not correlated with any environmental variable.

MRPP showed that fish assemblages differed (P < 0.001) among the five dominant benthic habitat categories (Table 4). Pairwise comparisons of fish assemblages indicated similarities among the benthic habitat types, with Rubble being central, sharing similarities with Sand, Montipora sp. coral bed, and Leptoseris sp. coral bed fish assemblages, whereas the fish assemblages of Hard Bottom only related to those of Leptoseris sp. coral beds (Figure 6). Fish assemblages also significantly differed (P < 0.001) between the six depth bins (Table 5). However, pairwise comparisons indicated relationships among the four shallowest depth bins (50–95 m), with the fish assemblages of the 2 deepest depth bins (96–105 m, 106–115 m) showing no similarities to the shallower depth bins or to each other (Figure 7).

TABLE 4

Table 4. MRPP results of benthic habitat comparison.

FIGURE 6

Figure 6. Multi-response permutation procedure (MRPP) relationships for benthic habitat categories. Connected benthic habitat types are not significantly different from one another in terms of their fish assemblages. P-values between the benthic habitat categories are next to the connecting lines.

TABLE 5

Table 5. MRPP results of depth categories comparison.

FIGURE 7

Figure 7. Multi-response permutation procedure relationships for the depth domains. Connected depth types are not significantly different from one another in terms of their fish assemblages. P-values between the substrate categories are next to the connecting lines.

Four significant Indicator Species, with values greater than 50%, were documented (Table 6). Two species were indicators of the deepest depth category (106–115 m): Pseudanthias hawaiiensis and juvenile Parapercis schauinslandii. Another two species were indicators of a hard-bottom substrate: Pseudanthias hawaiiensis and Odontanthias fuscipinnis. Centropyge potteri, was an indicator of Leptoseris sp. coral beds.

TABLE 6

Table 6. Groups and their indicator species.

Discussion

NMDS and HCA discerned fish assemblage associations with benthic substrates and depth in our Au‘au Channel data. NMDS Axis 1 showed the strongest association with increasing depth, hard bottom, and Leptoseris coral cover, whereas Halimeda algal fields were associated with Axis 3. Surveys conducted throughout the Hawaiian Islands by Rooney et al. (2010), Blythe-Skyrme et al. (2013), Costa et al. (2015), and Pyle et al. (2016) corroborate these findings, with Leptoseris coral sp. coverage increasing between the depths of 75 and 110 m and then diminishing at greater depths. Halimeda algal fields are prevalent throughout the Au‘au Channel at 90–100 m depths (Rooney et al., 2010). Leptoseris coral sp. and hard-bottom substrates increased with depth and supported a diversity of fishes including highly abundant planktivores such as Chromis verater and C. leucura. Increases in specific fish trophic groups associated with increasing depth and hard bottom have been described by Thresher and Colin (1986) who saw an increase in planktivores and Asher et al. (2017) who observed greater abundance and species in mobile invertivores.

Oxycheilinus bimaculatus is a wrasse that has a preference for algal beds and rubble fields (Randall, 2007). Halimeda fields comprise much bottom area in the Au‘au Channel (Rooney et al., 2010) and are bordered by transition zones of rubble or sand. Both Oxycheilinus bimaculatus and Halimeda are associated with axis three as indicators of an algal field habitat.

Montipora coral beds, another prevalent substrate in Au‘au Channel (Rooney et al., 2010), occurred only between 50 and 70-m depths; this habitat was populated by small-bodied and juvenile fish, primarily by juvenile and adult Cirrhilabrus jordani. Montipora coral beds were not strongly associated with a specific NMDS axis; however, they were most positively associated (tau value = 0.127) with Axis 2 and both life stages of Cirrhilabrus jordani. Axis 2 thus represents the Montipora coral bed habitat found at 50–70 m.

Large algal beds became obscured in our comparison between fish assemblages and substrate because the surveys that had algal beds contained larger areas of sand, rubble, or Montipora coral. The dominant cover was represented by five groups: Sand, Rubble, Hard Bottom, and Leptoseris and Montipora corals. Rubble was the central substrate with Sand, Montipora, and Leptoseris all associated with Rubble, but with no inter-relationship between each other. This result indicates that rubble is a transitional substrate that borders all of these habitats. Leptoseris coral also was associated with Hard Bottom, most likely because Leptoseris usually occurred in areas that had exposed Hard Bottom (Rooney et al., 2010; Blythe-Skyrme et al., 2013). This is reflected by the axes on our NMDS plot: Axis 3 represents the Halimeda sp. algae (Sand) and transition zones (Rubble), whereas Axis 1 represents the separate Leptoseris coral and Hard Bottom habitats.

Fish assemblages thus were inter-related in a manner that created two groups. The first set of inter-relationships indicates a clear depth range of associated fish assemblages: 50–95 m. This depth pattern seems to form an “upper” mesophotic fish assemblage in the Au‘au Channel. Deeper than 95 m the assemblages are unrelated to each other or the “upper” group which may mean that these depths are either a transition zone to deep ocean fish assemblages or that they form a “lower” mesophotic reef fish assemblage. Pinheiro et al. (2016), Semmler et al. (2016), Baldwin et al. (2018), Rocha et al. (2018), and Weijerman et al. (2019) describe similar breaks in “upper,” “lower,” and “rariphotic” fish assemblages at similar depths.

Four fish species showed relationships with either depth or substrate, with one appearing as an indicator for both. Adult Pseudanthias hawaiiensis and juvenile Parapercis schauinslandi were indicators of the deepest depth range of 106–115 m. Parapercis schauinslandi primarily occupies the sand and rubble that is found to greater extents at this depth. The P. schauinslandi that occur here appear to be mostly juveniles or small adults, which might indicate that the deep sand fields at the edge of Leptoseris beds are locations for juvenile recruitment, analogous to what Brokovich et al. (2007) found for a pomancanthid angelfish in the Gulf of Aqaba.

Centropyge potteri, Pseudanthias hawaiiensis, and Odontanthias fuscipinnis were indicator species for Hard Bottom. Pseudanthias hawaiiensis and O. fuscipinnis are usually found in areas of high relief, such as steep slopes or walls (Randall, 2007), corresponding to our Hard Bottom category. Pseudanthias hawaiiensis is a known deep-water species (Randall, 2007) and Hard Bottom increased as depth increased, further linking this species to the substrate and depth strata. NMDS results associating Centropyge potteri with rugose Leptoseris coral beds are consistent with the species preference for coral reef habitat (Randall, 2007).

ISA and MRPP results are further consistent with the original NMDS analysis, where Pseudanthias hawaiiensis and Centropyge potteri positively associated with Axis 1, the axis of Depth (Pseudanthias hawaiiensis), Hard Bottom (Pseudanthias hawaiiensis), and Leptoseris sp. coral (Centropyge potteri). This result indicates that mesophotic fish assemblages are associated primarily with specific benthic habitats similar to what was observed in the Gulf of Mexico by Dennis and Bright (1988), and in Puerto Rico by Bejarano et al. (2010) and Garcia-Sais (2010).

In the Au‘au Channel of Hawai‘I we observed what appeared to be a distinct “upper” mesophotic fish assemblage and either a “lower” or transitional assemblage as observed by Pinheiro et al. (2016), Semmler et al. (2016), Baldwin et al. (2018), and Weijerman et al. (2019).

Our analysis (NMDS, HCA, MRPP, and ISA) indicate that the mesophotic fish assemblage of the Au‘au Channel can be divided into an “upper” (50–95 m) and “lower” (96–115 m) group with specific species associated with specific benthic habitats. This is similar to what Stefanoudis et al. (2019) observed with mesophotic fish assemblage variance associated primarily with depth followed by benthic habitat. Further studies of fish (and other) assemblage components of mesophotic reef ecosystems are needed to further determine zonation of assemblages with depth and/or benthic habitat.

Data Availability Statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Ethics Statement

No animals were sampled by this project and thus there was no requirement for ethical review and approval for an animal study.

Author Contributions

RB, FP, and JR conducted the surveys. KH, ED, and RB conducted the data analysis. RB wrote the manuscript. All authors contributed to the article and approved the submitted version.

Funding

Funding for remote camera surveys, submersible dives and technical scuba diving were awarded under NOAA grants: NA10OAR4300132, NA07NOS4780187, NA07NOS4780188, NA07NOS4780189, NA07NOS4780190, NA05OAR4301108, and NA09OAR4300219.

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Acknowledgments

We thank C. Kelley for translating the submersible tracks and plots and C. Bradley for her advice and for being a great sub diving buddy. Funding and support was provided through the NOAA NCCOS Deep-CRES Program.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmars.2021.785308/full#supplementary-material

Footnotes

^ FishBase.org

References

Abesamis, R. A., Utzurrum, J. A., Raterta, L. J., and Russ, G. P. (2020). Shore-fish assemblage structure in the central Philippines from shallow coral reefs to the mesophotic zone. Mar. Biol. 167, 1–15.