Combining Methods to Better Estimate Total Fish Richness on Temperate Reefs: The Case of a Mediterranean Coralligenous Cliff

: Coralligenous habitat is considered as one of the most important special habitat types in the Mediterranean; however, due to its inaccessibility, little is known about it, although it is considered as one of the Mediterranean’s richest habitats in terms of species. Due to a low number of studies, it was presumed that the richness of coralligenous ﬁsh assemblages is underestimated using traditional visual census methods which are not applicable to the deep, steep, and vertical slopes of coralligenous cliffs and do not capture exhaustively cryptobenthic species commonly found in this habitat. This paper aims at producing a more complete assessment of ﬁsh assemblages on a coralligenous cliff by combining different methods, particularly the deep vertical transect visual census and square with anesthetics method. A total of 76 ﬁsh species were recorded on a single coralligenous cliff, supporting the opinion that coralligenous cliffs are important Mediterranean biodiversity hotspots. The analysis of species traits between species recorded by the different methods showed how complementary they are to better describe species compositions. Hence, the result of this study demonstrates that the combined use of methods is essential for a more exhaustive description of the whole ﬁsh community structure and for accurate estimates of the abundance and diversity patterns, particularly in complex habitats such as coralligenous cliffs. investigation, A.S., I.G. and M.K.; resources, A.S., I.G. and M.K.; data curation, A.S., I.G. and M.K.; writing—original draft preparation, A.S.; writing—review and editing, A.S., I.G. and M.K.; visualization, A.S.; supervision, A.S. and M.K.; project administration, A.S. authors


Introduction
The marine environment is extraordinarily diverse and made up of a large number of different marine habitats. They are identified by physical structures or specific environmental conditions. However, a large portion of the ocean remains unexplored, even the coastal habitats, particularly benthic assemblages thriving from 20 m depth down to the end of the continental shelf. Therefore, their distribution, ecological roles, threats, and status remain poorly known. In that zone, a typical Mediterranean habitat includes various coralligenous formations. Due to their inaccessibility, little is known about coralligenous habitats; however, they are considered as one of the Mediterranean's richest habitats in terms of species [1]. The ecological importance of coralligenous habitats and their scientific and biodiversity interest is recognized by international conventions (e.g., Barcelona Convention); hence, they can be considered one of the most important "special habitat types" that should be assessed under the Marine Strategy Framework Directive [2,3] through accurate monitoring plans [1].
Coralligenous formations are indicative of a circalittoral biocoenosis consisting of a three-dimensional biogenic build-up that forms a solid substrate primarily dominated by coralline algae [4]. Due to the variety of formations, there is no generally accepted definition; however, one of the most used defines coralligenous formations as a hard substrate of biogenic origin mainly produced by the accumulation of calcareous encrusting algae growing in dim light conditions [5]. They develop on both hard and soft bottoms, at 4 to 160 m of water depth, and adapt to low levels of nutrients and temperature from 10 to 23 • C in waters with moderate hydrodynamics [4,5]. Although the morphology and inner structure of coralligenous formations depend on depth, topography, and the nature of prevailing algal builders, generally, the two geomorphotypes of coralligenous formations can be distinct: cliffs that are vertical or near-vertical walls from a steep littoral rock face, and banks that can be defined as flat frameworks with thickness ranging from 0.5 to 4 m mainly built over more or less horizontal substrate [4][5][6]. Coralligenous formations are composed of a wide number of species that create a complex habitat attracting variability of mobile species. Although sessile coralligenous assemblages have been studied, little is known about coralligenous fish assemblages [7]. The species richness of coralligenous fish assemblages was considered to be lower than that of shallow rocky habitats [8], but later studies showed that the species richness number per single coralligenous site is similar to fish assemblages of shallow rocky habitats [7]. Moreover, even those few studies presumed that the richness of coralligenous fish assemblages is underestimated using conventional visual census methods [9] and because of cryptobenthic species that are not recorded by visual census methods [7,10]. Such presumption was supported by the recent study of Soldo and Glavičić [11], who researched coralligenous habitat and suggested that underwater steep rocky coralligenous reefs have the highest fish species richness of any rocky habitat in the Mediterranean.
Monitoring methods change according to the objectives of each study. Although traditional methods using fishing gear have been used for the assessment of coralligenous fish assemblages [12], they are not applicable for steep cliffs. Currently, the most common observational method for reefs, in general, is an underwater visual census (UVC). A variety of UVC methods have been developed, but those focusing on reef fishes are mainly based on the use of SCUBA (Self-Contained Underwater Breathing Apparatus) divers conducting a survey, due to their better ability to search within such complex habitat compared to a stationary camera used by different remote underwater video methods [11,13]. However, as steep and vertical coralligenous reefs usually extend to depths deeper than 30 m, the main disadvantage of conventional SCUBA diving is depth restriction because of safe diving considerations. The build-up of residual nitrogen in the diver's blood dictates short maximum dive times. Furthermore, as the dive depth increases, the allowable dive time decreases and requires a sufficient no-dive recovery period between dives [11]. The additional reason for depth restrictions in conventional SCUBA diving derives from limitations in commonly used existing technology [11]. Most divers worldwide use opencircuit diving gear with compressed air as breathing gas. Recently, an increasing number of divers have been experimenting with mixed-gas diving technology, including closed-circuit rebreather gear, to extend the depth limits of conventional SCUBA [14,15], but the high cost of equipment and requirement for extensive advanced training restricts such application to a small group of highly trained scientific divers [11].
In addition to diving related problems in applying UVC methods on coralligenous cliffs, an issue which impacts the quality of the assessment is the high small-scale heterogeneity of coralligenous outcrops that are very complex structures composed of several microhabitats [5]. Such habitat complexity usually impacts the quality of the standard UVC methods that can miss a large number of fish species, particularly those belonging to a cryptobenthic group [16][17][18], as defined by Kovačić et al. [18]. Soldo and Glavičić [11] recently reported six gobiid species within 51 fish species recorded from the surface to 50 m depth on the coralligenous cliff using the deep vertical transect (DVT) method. However, as the visibility of the species is a prerequisite of any visual census, bearing in mind the definition of cryptobenthic species by Kovačić et al. [18], it can be presumed that some species cannot be detected by the visual census, regardless of the method and technology used, as they are simply not visible at any moment.
Thus, it can be hypothesized that, for an improved assessment of fish assemblages on highly complex habitats, a combination of methods should be used. Considering the high complexity of the coralligenous cliffs and the recent results that suggest those habitats have the highest fish species richness of any rocky habitat in the Mediterranean [11], this paper aims to combine and compare methods for producing a complete assessment of fish species composition on a coralligenous cliff.

Sampling Location
The study was carried out at the underwater coralligenous cliff in the area of the island of Hvar channel (Figure 1), central eastern Adriatic (43 • 14 02 N, 16 • 34 29 E). The cliff Smočiguzica is elliptically shaped with a truncated and oblique cone located 300 m from the coastline. gardless of the method and technology used, as they are simply not visible at any moment.
Thus, it can be hypothesized that, for an improved assessment of fish assemblages on highly complex habitats, a combination of methods should be used. Considering the high complexity of the coralligenous cliffs and the recent results that suggest those habitats have the highest fish species richness of any rocky habitat in the Mediterranean [11], this paper aims to combine and compare methods for producing a complete assessment of fish species composition on a coralligenous cliff.

Sampling Location
The study was carried out at the underwater coralligenous cliff in the area of the island of Hvar channel (Figure 1), central eastern Adriatic (43°14′02″ N, 16°34′29″ E). The cliff Smočiguzica is elliptically shaped with a truncated and oblique cone located 300 m from the coastline. A large part of the cliff has a high slope with the highest point of 9.1 m and the deepest at 60.7 m. The cover is coralligenous and it is surrounded by a sandy bottom. Apart from the depth changes, this single habitat on a single locality is surprisingly homogenous, since the steep inclination prevents the forming of the mixed bottom by sand sediment or the accumulation of multilayers of boulders and pebbles across the bedrock surface, except at the foot of the cliff at 60+ m depth. Sessile species covering each transect varied by depth and included both algae and animals, but the dominant species, apart from red algae, were green algae (Chlorophyta) Codium bursa (Olivi) C. Agardh 1817 (from 9 to 15 m) and  A large part of the cliff has a high slope with the highest point of 9.1 m and the deepest at 60.7 m. The cover is coralligenous and it is surrounded by a sandy bottom. Apart from the depth changes, this single habitat on a single locality is surprisingly homogenous, since the steep inclination prevents the forming of the mixed bottom by sand sediment or the accumulation of multilayers of boulders and pebbles across the bedrock surface, except at the foot of the cliff at 60+ m depth. Sessile species covering each transect varied by depth and included both algae and animals, but the dominant species, apart from red algae, were green algae (Chlorophyta) Codium bursa (Olivi) C. Agardh 1817 (from 9 to 15 m) and

Data Collection
Two methods were used for sampling: the standard (DVT) method [11] and the square with anesthetics method (SAM) [18]. In addition, every diving expedition had some dive time not spent on the methods but on preparation and swimming from the boat to the studied surfaces and back, as well as the dive time and swim between two squares for 4 of 15 SAM. The species recorded during that time and not recorded by any of the methods on any of the studied surfaces are separately listed in the results.
The diving expeditions with DVT method were performed on sunny days from mid-June until the end of August 2020. The total number of performed transects was 216. The DVT method is based on the use of 50 m long transects with 5 m height (being a vertical transect) as described in Soldo and Glavičić [11]. Two observers dive together, during which the first diver uses a high-definition (full HD) video camera and swims through the transect by browsing within it (S-type transect), varying in depth and angle and zooming when needed or notified by another diver of the presence of small or hidden fish individuals [11]. Due to the low light level at the dive depths, custom-made LED lights 2 × 25 W, producing a light temperature of 6400 K (sunlight color temperature), were used during video recording for easier post-identification of the species. Considering the high efficiency of the DVT method for the assessment of fish assemblages on vertical coralligenous cliffs [11], that method was used for collecting data on most of the fishes inhabiting or associating with the coralligenous cliff. Additionally, due to the wider base of the cliff, which was partly outside of the boundaries of the standard 50 m transect, a modified DVT method was used in a way that the length of the transect was not strictly 50 m but usually longer in order to sample the whole cliff area. Hence, the modified DVT method was applied mostly on areas down to 62 m depth, while transects shallower than 10 m were not sampled, as those areas were part of the standard transects. All DVT data were pooled together.
The SAM using diving and anesthetics (or alternatively ichthyocides) is the only method in general that can ensure efficient collection of truly cryptobenthic fishes [18]. According to Kovačić et al. [18], a cryptobenthic fish is a species (or a life history stage of a fish species) whose individuals exclusively or predominantly spend their lifetime in cryptobenthic microhabitats, i.e., in the restricted living spaces underneath the bottom surface of the substrate or biocover, with a physical barrier to open spaces. Hence, the same authors defined an epibenthic fish as a species or life history stage of a species that lives exclusively or predominantly on top of the bottom surface. Under this definition, cryptobenthic fishes should be distinguished from epibenthic fishes that use camouflage colors and shapes or shallowly bury themselves in the soft sediment surface. Moreover, cryptobenthic fishes should also be distinguished from any other predominantly epibenthic species or life history stages of species, which only seek shelter underneath the surface for a limited period of their total time, for example, after escape from predators, for resting, or for breeding. Otherwise, for example, all or almost all benthic gobiid species would be assigned to cryptobenthic fishes. Cryptobenthic fishes should also be distinguished from marine cave-dwelling fishes if the latter live epibenthically on top of the surface of a substrate or they swim free inside large cave spaces [18]. The SAM data were collected in the period 2012-2017 and were partially published for the data in 2013 [15] and for the data from autumn 2014 [19]. Considering the qualitative characteristic of the used data, the sample time interval length and the sample time difference from the DVT method did not influence the results. The total number of performed squares was 54, with each square having 1 m 2 surface size. On all SAM diving expeditions, the method followed the protocol on squares, as described by Kovačić et al. [18], including the visual census of squares before anesthetic application. The anesthetic used was quinaldine diluted 1:15 with 96% ethanol and then mixed 1:5 in 750 mL bottles with seawater. Two bottles were usually used for a square, i.e., the total volume of the deployed quinaldine-ethanol solution was about 300 mL/m 2 . Therefore, the SAM results consisted of the specimens observed by visual census on the square and of the specimens found only after deploying anesthetic, summarized separately.
In all methods, data and specimens were always collected by two SCUBA divers. Each dive started around 1:00 p.m., in order to have the best possible light conditions. DVT sampling ranged from the surface to 62 m depth, while SAM was performed from 15 m to 62 m, since the coralligenous outcrops as a complex structure composed of several microhabitats were not present between the surface and 15 m depth. Depending on the depth of sampling, different breathing gases were used. For depths down to 50 m, compressed air was used as a breathing gas, while, for deeper depths, trimix 20/30 was used as a bottom gas. During all dives, nitrox 50 was used as a decompression gas. Trimix enables safe diving at greater depths and nitrox enables a significant reduction in decompression time (a 50% reduction compared to decompression with compressed air), as well as safer diving [11,19]. The change of breathing gases was performed at a depth of 21 m, corresponding to oxygen partial pressure in nitrox of 1.6 bar.
Post-analysis of the video recorded by the DVT method was carried out on a 55 screen, using standard software that enables pausing, slow-motion viewing, and zooming of the images. Video analysis was conducted within 24 h of the dives (to use information still fresh in the memory of divers if necessary) with all fish identified and counted per species. The specimens collected by the square method were put to death after reaching the surface by a quinaldine overdose and stored in 70% ethanol solution. The preliminary field species identification of cryptobenthic specimens was later rechecked in the lab.

Data Analysis
The species compositions of the DVT method and SAM were incomparable quantitatively due to differences in surface size, time, and number of samples between methods. Therefore, the species compositions were compared qualitatively. The similarity of DVT and SAM qualitative species compositions was calculated by the Jaccard/Tanimoto coefficient, a binary similarity coefficient which is the ratio of the samples' intersection to their union. For each species observed by the DVT method, the frequency of occurrence was calculated according to the number of dives when the species was recorded divided by the total number of dives. Similarly, the species frequency of occurrence on SAM squares was calculated as the number of squares with species present divided by the total number of performed squares. The values of frequencies of occurrence were also not comparable as absolute values between methods, considering the large difference in the surface size of more than two orders of magnitude between DVT transects and SAM squares, which causes more frequent species occurrences on DVT transects. Therefore, the ranks were compared, and the significance of association between the ranked orders of species frequency of occurrence from DVT transects and SAM squares was tested by the nonparametric Kendall coefficient of rank correlation, τ [20]. The traits of species occurring in the DVT and SAM samples were compared. The species traits were divided into categories as follows: maximum size (very small fish > 6 cm total length (TL), small fish > 6-10 cm TL, medium sized fish > 10-50 cm TL, large fishes > 50 cm (TL); relationship to the bottom (pelagic fish, benthopelagic fish, hyperbenthic, epibenthic, cryptobenthic); relationship to the reef (reef-associated, occasional), trophic level (2-2.99, 3-3.99, 4-4.7), catchability by commercial and recreational fishing tools (catchable, non-catchable), tangible species vs. elusive species, measured indirectly by the century of description (18th, 19th, 20th, 21st), and taxonomy (by the taxon level of orders). Maximum size and trophic level followed Froese and Pauly [21]. Trophic levels were estimated as described in Pauly et al. [22]. Trophic levels in aquatic environments range from 2 for herbivores and detritivores to 5.5 for specialized predators of marine mammals; for fishes in particular, the trophic levels generally range from 2 for the detritus feeding species to 4.7 for piscivorous species [22]. The relationship to the bottom, the relationship to the reef, and catchability were estimated from the species data in Dulčić and Kovačić [23]. The taxonomy and the century of description for each species followed Van der Laan et al. [24]. The significance of differences in each trait was tested by the Fisher exact test [25]. The Fisher exact test, despite being conservative, is more appropriate than other tests of independence for the contingency tables with some cell counts lower than 5 [25]. Statistical analyses were performed in PAST version 2.17c [26].

Results
Total observed fish richness on this single reef with a homogenous habitat of steep bedrock was 76 fish species, with 56 of them recorded by DVT, 28 recorded by SAM, and five randomly noticed during DVT and SAM divings, but outside of the studied area and applied methods ( Table 1). The species composition overlap between DVT and SAM was just 13 species. Therefore, visual census on large transects, as the main method, was able to capture 73.7% of total fish species richness on locality, while SAM recorded 36.8% of the total fish species richness on locality. Only 17.1% of all species overlapped between DVT and SAM ( Table 1). The DVT and SAM low overlap in species composition represented about 1/4 of DVT species and about 1/2 of SAM species. The Jaccard/Tanimoto coefficient of similarity between DVT and SAM species composition was low, 0.13, showing low similarity in species composition. The DVT method recorded 45 reef associated species and 11 occasional species, the latter being defined as mainly active and more mobile predators that use the underwater cliff as an important feeding area (Table 1). Within the total observed fish richness, 19.7% of fishes were recorded exclusively by SAM. From 28 SAM species, 13 were observed visually on squares, and 22 were driven out by anesthetic; therefore, seven species were ambivalent (i.e., epibenthic/cryptobenthic) in relation to the bottom and 15 were exclusively cryptobenthic (Table 1). Only five species of total fish species richness were recorded outside of the transects and squares of both methods and were randomly recorded during diving, three of them being tidal blenniid species recorded at the start or end of diving expeditions. The remaining two species are real inhabitants of infralittoral reefs, while the lack of regular occurrence of L. candolii could be explained by the absence of multilayers of boulders and pebbles on a vertical cliff, and the lack of regular occurrence of S. ocellatus could be explained by the absence of high vegetation or highly structured shallow rocky bottom. The ranked orders of species frequency of occurrence of each method showed a general mismatch between the two methods ( Table 2). The Kendall coefficient of rank correlation showed an absence of a significant association between the ranked orders of species frequency of occurrence between DVT and SAM (τ = −0.119, p = 0.127). The species recorded by the DVT method and SAM show a significant difference in all species trait compositions ( Table 3). The DVT recorded species significantly larger and less associated with the bottom than SAM, with the DVT species size modal value at 10-50 cm total length vs. the SAM species size modal values at ≤6 cm total length (Table 3). Hyperbenthic species were the most numerous DVT species vs. cryptobenthic as the most numerous SAM species (Table 3). The DVT species composition was significantly different from the SAM method in relation to the reef, with 11 occasional species recorded, while none of the occasional species were found by SAM ( Table 3). The DVT recorded species had significantly different trophic levels than the SAM recorded species, having more omnivorous and herbivorous species than SAM, as well as species with a trophic level above 4, suggesting a wider trophic range, although the modal value for both groups was at 3-3.99 (Table 3). The DVT species composition was significantly different from the SAM species composition in exposure to fishing tools, having prevailingly species catchable by fishing tools, while the majority of the SAM species were not possible to collect by any conventional fishing method (Table 3). SAM also mostly recorded species described later (19th, 20th, 21st), which was a significantly different composition by century of description to DVT species which were mostly described in the 18th century (Table 3). SAM had significantly different composition by taxonomy on the taxon level of order, compared to DVT species composition, with the species richest order being Gobiiformes, with 14 species, versus the species richest order Perciformes in the DVT species composition, with 31 species. Table 1. Species recorded by DVT method: + reef-associated species or species which permanently occupy the coralligenous surface or habitats around the reef base (sandy), (+) occasional species or species which are mainly active and more mobile predators that use reef as feeding areas. Species recorded by SAM: + visually observed specimens on 1 m 2 squares, * specimens collected with anesthetics on 1 m 2 squares. Species from DVT, modified DVT, and SAM dives visually noticed during those dives, but not recorded by the methods, marked with + in a separate column.

Family
Method

Discussion
La Mesa et al. [27] compared fish richness in different locations and habitats of the Mediterranean including sandy habitat to Posidonia oceanica beds and rocky reefs. Overall, they recorded 59 species with rocky habitat being the richest, having a total of 43 species. Higher species richness on rocky reefs can be explained by habitat structure, which is recognized as one of the factors likely to explain the variability of Mediterranean fish assemblages [28]. García-Charton et al. [28] presumed that higher heterogeneity and complexity of habitat will increase the number of species and individuals. Nevertheless, despite coralligenous habitat complexity, the species richness of coralligenous fish assemblages was previously considered to be lower than that of shallow rocky habitats [8]. However, Piazzi [7] recorded 61 fish species in a coralligenous habitat, but at 14 different locations, where the number of species per location varied between 35 and 46. Thus, Piazzi [7] concluded that the mean value of species number per single coralligenous site is similar to fish assemblages of shallow rocky habitats. Hence, it was presumed that the richness of coralligenous fish assemblages was underestimated using visual census methods which are not applicable to the deeper, steep, and vertical slopes of Mediterranean coralligenous cliffs [9] and do not exhaustively record cryptobenthic species that are commonly found in this habitat [7,10]. Consequently, Soldo and Glavičić [11] developed a new transect method for coralligenous habitat to an extended depth of 50 m and recorded a total of 51 species, thus suggesting that underwater steep rocky coralligenous reefs have the highest fish biodiversity and population density of any rocky habitat in the Mediterranean. The same conclusion can be derived from the results of this study, as the total of 76 fish species recorded on a single coralligenous cliff certainly supports the presumption about coralligenous cliffs as Mediterranean biodiversity hotspots.
Although this study was conducted up to 62 m of depth, the use of new diving techniques to extend the depth limits of conventional SCUBA worked well with sampling methodologies, which demonstrates that usual depth restrictions in conventional SCUBA diving, derived from the use of the open-circuit diving gear with compressed air as breathing gas, could be surpassed by mixed-gas diving technology.
Despite discussed limitations, visual census on large transects remains the crucial method for the assessment of fish assemblages, especially on simple and/or relatively large habitats, being able to capture the majority of species. However, results from this study show that, for more complex habitats, such as coralligenous cliffs, a significant portion of total species richness, about one-quarter, is not recorded by the transect method. Although the DVT method, both standard and modified, recorded a significant portion of small body size species, owing to advances in video high-resolution technology and appropriate sampling technique, it is obvious that some species cannot be detected by any method of visual census, because they entirely or predominantly spend their lifetime hidden and, thus, not visible. Therefore, in such cases, the SAM method was shown to be a necessary component in total species richness estimation, since nearly one-fifth of total species richness was recorded only by SAM. The DVT and SAM poor overlap in species composition showed that these methods are truly complementary and should be applied together in studying total reef fish diversity. The ranked orders of species frequency of occurrence of each method showed a general mismatch and no significant association between the two methods ( Table 2). Until now, combined methods for estimating total fish species richness of a single habitat have not been used, but some studies have compared the effectiveness of the visual census and the anesthetic/ichthyocide census techniques for determining richness of cryptobenthic species [29][30][31]. The reported portion of cryptobenthic species underestimated by the visual census method from those studies is in accordance with the results of this study.
While the higher effectiveness of combined methods versus the visual census alone was expected [29][30][31], the complementary characteristic of these methods and the degree of effectiveness when they were combined were revealed by the present data. The species recorded by the DVT method and SAM showed a significant difference in all species trait compositions (Table 3). A similar comparison of species trait compositions between these methods was not done before. The DVT species were, on average, significantly larger and more distant to the bottom or surface of the cliff. The DVT method also recorded occasional species, not found by SAM. The significantly different trophic levels recorded between methods were caused by the wider trophic distribution of DVT-recorded species compared to the SAM-recorded species. The DVT method revealed more omnivorous and herbivorous species at the lower trophic levels, as well as a number of species with a trophic level above 4, not present in any of the SAM-recorded species. The DVT method had prevailingly species catchable by fishing tools, while the majority of SAM species were not possible to collect by any conventional fishing method. SAM also covered significantly more poorly known and elusive species compared to the DVT method, revealed indirectly by the later century of species description. SAM, compared to DVT species composition, showed different species composition by taxonomy on the taxon level of order. Therefore, the analysis of species traits showed in which way and to what degree SAM has complementary species composition to the DVT method; SAM targeted smaller, hidden species not collectible by any conventional method and, therefore, elusive, mostly belonging to different taxonomical orders than DVT species. With just two real inhabitants of infralittoral reefs found independently of the two methods, the combination of DVT and SAM was shown to be an efficient tool in estimating total fish species richness on coralligenous cliffs.
The number of species found on a single coralligenous cliff during this study by combining methods is clearly higher than that found during previous single-method studies [7,8,12] when 35-46, 43, and 45 species were identified, respectively. Moreover, the number is higher even compared to studies from the Adriatic Sea [11,32,33] where 51, 39, and 33 species were reported, respectively. The most significant difference regarding species composition between this and previous studies is in the recorded number of species belonging to Gobiidae, Blennidae, and Gobiesocidae families (17,5, and 2 species, respectively). Thus, this study shows that, when the aim is to determine total fish species richness of the marine complex habitat, relying on only one method results in more or less a significant underestimation of species richness, depending on the complexity of the investigated habitat. The incomplete and inaccurate registration of the fish community composition affects the efficiency of conservation planning and the implementation of appropriate conservation plans, which can be, consequently, inadequate and deficient.
The results of this study demonstrate that the combined use of the methods is essential for a better description of the whole fish community structure and for more accurate estimates of the abundance and diversity patterns, particularly in complex habitats such as coralligenous cliffs. Only when a community is completely known, further actions can be conducted to fully understand its trophic structure and ecosystem functioning. Furthermore, having precise and complete data on community structure is essential to define monitoring schemes required for the assessment of the conservation status of these highly sensitive and threatened Mediterranean marine habitats, especially related to potential changes in species structure and composition in reaction to multiple threats.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author.

Conflicts of Interest:
The authors declare no conflict of interest.