Spatial variation of reef fishes and the relative influence of biotic and abiotic habitat traits

Abstract

Patterns of distribution of reef fishes were examined across three spatial scales and related to habitat traits along 25 km of the northern Portuguese coast. Response variables included the multivariate assemblage structure, the total number of taxa and individuals, and the abundance of single groups categorized according to their preference for the benthic, proximo-benthic or pelagic environment, feeding and reproductive behaviour. Habitat traits included topographic elements (small and large ‘drops’ like cracks and crevices) and the extent of dominant morpho-functional types of macroalgae (kelp, large foliose, small erect, turf-forming filamentous, and encrusting). All fish responses were characterized by the largest variance at the smallest scale (among transects tens m apart), followed by that among reefs (hundreds m to 1 km apart) and almost null variance among sites (some km apart). Small and large ‘drops’ of the substratum explained, respectively, considerable variation of assemblage structure and the total abundance of individuals, while the extent of bare rock influenced the richness of taxa and that of benthic fishes, fishes feeding on sessile invertebrates and fishes laying benthic eggs or having nesting behaviour. Combinations of abiotic and biotic structural attributes of reefs influenced proximo-benthic fishes, the predators of mobile animals and fishes releasing pelagic eggs. The here reported associations between patterns of distribution of reef fishes and habitat traits have implications for the design of future protection schemes suitable to guarantee the conservation of reef fish communities and of the processes responsible for their variation. Within the SLOSS (single-large vs. several-small) debate in the design of marine reserves, for example, effective protection to the studied reef fishes would be provided by a set of small reserves, rather than a single large which might be appropriate for fishes having wider home ranges.

Background

There is evidence that patterns of variation of natural populations and assemblages are more evident at some scales than others [1,2,3,4]. This characteristic is crucial for a number of basic and applied ecological issues, including the notion that, since different processes are likely to operate differently in space and time, the identification of relevant scales of variation is essential to propose and then test their potential explanatory models [5, 6]. For example, such an identification may be important to assess and predict the effects of anthropogenic disturbances operating at a hierarchy of spatial scales [7] and to support decisions about the design and management of effective protected areas [8,9,10].

The spatial patterns of biodiversity in coastal marine, especially rocky intertidal and shallow subtidal, systems have been widely examined, although across a range of scales usually smaller compared to that of terrestrial studies [but see 3, 11,12,13,14]. However, most analyses have identified a common pattern of relevant small- (tens to hundreds of centimetres) to middle-scale (tens to hundreds of metres) variation across several organisms and habitats [2,3,4, 6, 15, 16], highlighting the general principle that small-scale processes may be as, if not more, important as large-scale processes in driving patterns of distribution, abundance and diversity of coastal assemblages. As a consequence, small-scale natural variability should not be considered a statistical inconvenience, but a key issue to address in order to increase ecological understanding [17].

Shallow subtidal rocky reefs provide essential and heterogeneous habitat for a large diversity of organisms, supporting some of the most productive communities in the world [18]. Physical and biological variables have been indicated as major direct or indirect drivers of the small-scale spatial variation described for various reef populations and assemblages [19,20,21,22]. Different rock types, changes in the biological coverage, including the dominance of habitat-forming algae and invertebrates, of the bottom, and topographic complexity can drastically affect patterns of distribution of organisms living in close association to the substratum [19, 23,24,25]. Larger-scale variation in the structure and dynamics of reef communities, instead, may depend on factors such as differences in wave exposure [26] or food and propagule supply [27, 28]. Reef-associated fish assemblages, in particular, vary in abundance and diversity in response to variability in environmental conditions over a range of scales [29,30,31,32,33]. These may include the variability of oceanographic conditions as driver of fish variation at spatial scales of hundreds of kilometres or more [34, 35], and the physical structure of the reef as driver of fish variation at scales from metres to some kilometres [36].

Several studies have examined the relationships between reef habitat traits and associated fish assemblages in tropical [37,38,39,40,41] and temperate [31, 33, 42,43,44] systems. The physical nature of the substratum [19], the relative availability of vertical, horizontal or gently sloping substratum [45] and the cover and physical structure of algal assemblages [30, 44, 46, 47], possibly affecting the intensity of local biological processes such as predation [48,49,50], are among the main habitat traits documented to determine the structure of reef fish assemblages. This knowledge made the direction and strength of the relationships between fish assemblages and the reef habitat a criterion for the implementation of management strategies of fisheries and the design and assessment of the effectiveness of marine protected areas [30, 31]. None of these objectives can be achieved by focusing on fish assemblages without taking into account their positive versus negative, direct versus indirect, more or less specialized association with relevant habitat traits. Moreover, most studies aimed at assessing the effectiveness of marine reserves are based on comparing one protected area with one or more adjacent non-protected areas, but such an approach can provide unreliable results if the variability of variables intended to indicate the success of the implemented protection is not properly taken into account [51, 52]. Reserve effects and their responsible processes can be correctly identified and predicted only through sampling and monitoring designs including adequate replication in space and/or time [53, 54]. As a direct consequence, the possibility to generalise such issues to not yet implemented circumstances critically depends on identifying patterns of distribution of fishes and the relative influence of a range of biological and abiotic traits in relatively unstudied and ecologically peculiar systems.

In this context, the Portuguese continental coast provides an ideal case-study due to the occurrence of several species at their southern or northern distributional range edges and the documentation of clear latitudinal shifts in patterns of distribution and abundance of a number of macroalgal, invertebrate and fish organisms [33, 55,56,57,58]. The cooler northern region, in particular, appears as where kelps, once dominating the entire Portuguese coast, are now restricted [33, 59]. This is likely due to the combination of the recent reduction of the intensity of upwelling events and warming, which could have cause a northward retraction and reduction of the environmental conditions suitable for kelp species, and the reduced grazing intensity in the northern compared to the southern regions [56, 57, 60].

This study adopted a hierarchical sampling design to partition the spatial variability at three spatial scales (tens of metres, hundreds to 1000 metres, and kilometres) of reef fish assemblages from the northern Portuguese coast. Such variability was correlated to a range of potentially responsible abiotic and biological habitat traits. In addition to indentifying and comparing the relevant scales of variation in the fish assemblage structure, richness of identified taxa and total number of individuals, the study specifically tested the following hypotheses regarding the spatial patterns of different fish groups categorized according to several life-history traits [61]: (1) the distribution of fishes belonging to the benthic, proximo-benthic and pelagic spatial guild [62] would be characterized, respectively, by increasingly larger relevant scales of variation; (2) smaller scales of variation would be displayed by fishes with feeding habits focusing mostly on sessile invertebrates compared to omnivores (those feeding on available food, including organic debris) and active predators of animals with larger movement abilities [63]; (3) the relative greater importance of smaller versus larger spatial scales would vary between fishes having benthic eggs or nesting behaviour and fishes having pelagic eggs and larvae which disperse over long distances [64]; (4) according to such functional differences, variation in local abiotic (type of substrate and topographic heterogeneity) and biotic (type and extent of dominant algal beds) habitat traits would explain a relatively larger proportion of variation of fishes having living, feeding or reproductive habits relying on smaller-scale processes. In general, by partitioning the examined fish assemblages into their component functional groups that exploit, for different purposes, different parts of the spatial domain, testing such hypotheses aimed at investigating how local selection is related to spatial patterns of distribution [65, 66]. Fishes living in open habitats would likely rely on defences based on large numbers spread over large areas [67, 68], while fishes associated with soft and rocky bottom habitat traits would likely rely on relatively smaller, but more heterogeneously distributed numbers [62, 69].

Methods

Study area and sampling design

The study area (Fig. 1) spanned about 25 km along the northern Portuguese coast, where numerous rocky reefs occur interspersed within the predominantly sandy bottom [60]. Within the whole study area, the sampling was conducted at each of three haphazardly selected sites (about 5 km long, some km apart), each one including five randomly chosen rocky reefs (hundreds of metres up to about 1 km apart, 5–12 m depth). The almost rectilinear coastline, oriented from north to south, made all reefs comparable in terms of exposure to prevailing waves and wind (Fig. 1). These reefs are not subjected to any implemented protection status and, to our knowledge, are equally exposed to local fishing, typically artisanal, activities.

Fig. 1

Map of the study area in northern Portugal, showing the location of the five reefs (from r 1 to r 5, interspersed within extensive areas of sandy bottom) sampled in each of three sites

At each reef, the number of all fishes was estimated by a SCUBA diver by means of visual census performed along five, randomly chosen, transects (25 m × 5 m belts, 20–50 m apart) [36, 70]. Most fishes were identified to the species level, in some cases lumping different individuals to the same genus when poor water transparency did not allow to distinguish species in a fully consistent way. This was the case, for instance, of labrids of the genus Symphodus, especially when represented by juvenile individuals. The level of identification of each fish taxon, however, was the same for all sampled transects and reefs. The number of large (> 1 m) and small (< 1 m) topographic elements (hereafter indicated as large and small ‘drops’, including cracks, crevices, caves and holes) of the substratum was also counted within each 125 m² transect used for the fish counts [71]. Finally, a second diver measured, along each transect, the relative extent of the portions which were characterised by different types of substratum (continuous rock, large boulders, cobbles/pebbles, sand) and dominated by different morpho-functional groups of algae, including kelp (e.g., Laminaria hyperborea and Sacchoriza polyschides), small erect (e.g., Corallina officinalis), large foliose (e.g., Dyctiota dichotoma), turf-forming filamentous (e.g., Polysiphonia spp.) and encrusting algae (e.g., Lithophyllum incrustans), plus bare rock.

The study was carried out, by the same two divers, between September 2010 and July 2012. During this period, safe and effective diving opportunities were critically constrained by the harsh weather and sea conditions characterizing the study area [72], making the sampling occasions rather opportunistic. Such a logistic difficulty prevented to formally test any time-related hypothesis. However, each set of five reefs within each site was sampled over a period of about 1.5 years and guaranteeing that at least two reefs were sampled in each of the ‘cool’ and the ‘warm’ period (December–May and June–November, respectively) of sea water in the region [73]. This allowed to keep the temporal variability of the sampling of each set of five reefs as much as possible comparable among sites, thus reducing, though not completely removing, the potential confounding between the specifically examined spatial variation and the not fully controlled temporal variation. In addition, the sampling was always done during the day and avoiding crepuscular periods.

Data analyses

Spatial patterns in the structure of whole fish assemblages were examined by a two-way permutational multivariate analysis of variance (PERMANOVA) [74], with ‘Reef’ nested in ‘Site’ and the five transects sampled in each reef providing the replicates. This analysis was based on Bray–Curtis dissimilarities calculated from square root-transformed data. The same model of analysis was used in analysis of variance (ANOVA) to examine variation in the total number of identified taxa (a surrogate measure for species richness), the total number of fish individuals and the abundance of each functional group of fishes. Before each ANOVA, homogeneity of variances was assessed with Cochran’s C test and data were transformed when necessary to remove heterogeneity and avoid increased probability of type I error.

To quantify and compare the spatial variability of univariate and multivariate response variables at each sampled scale (among transects, reefs and sites), variance components were calculated from ANOVA (as well as pseudo-variance components from PERMANOVA) [30] by equating empirical and expected mean squares [75, 76]. Obtained negative estimates of variance were interpreted as sample underestimates of very small to null variances and set to zero [76, 77]. All variance components were calculated from untransformed data [e.g., 78] and the contribution of each scale was finally expressed as the percentage to the total variance of each multivariate and univariate response variable.

The distance-based linear modelling (DISTLM) [79] routine was used to run multiple regression models to assess the contribution of each abiotic (small vs. large ‘drops’ and four types of substratum) and biotic (five groups of dominant algae and bare rock) habitat traits to the total variation of each response variable. Model selection was performed by first calculating second order Akaike information criterion values (AIC_C), which are more appropriate than raw AIC values when, such as in the present case, the number of observations is relatively small compared to the number of parameters in the model [80,81,82]. Then, for each model, the differences in AIC_C with respect to the AIC_C of the best candidate model were computed. Finally, second order Akaike weights were calculated as described in [80,81,82].

Results

Structure, richness and overall abundance of reef fish assemblage

A total of 22 fish taxa were identified during the study, with about 90% of the total number of individuals provided by eight most common taxa (Additional file 1: Supplementary material and Fig. 2). The structure of the whole assemblage varied significantly among reefs, although the smallest scale (among transects) provided the greatest contribution to the total (pseudo) variance. This was about six times and seventy-seven times greater than that of the reef and the site scale, respectively (Table 1 and Fig. 3a). Both the richness of taxa and the total number of individuals showed patterns of spatial variation very similar to those of the whole assemblage, with the former variable, in particular, being characterized by a null variance at the scale of sites (Table 1, Fig. 3b, c).

Fig. 2

Occurrence (proportion of transects where a fish taxon was found, out of the total number of sampled transects) of the eight most numerous fish species or genera observed during the study at each sampled site. The reported fish taxa are those collectively exceeding 90% of the total number of individuals sampled throughout the study

Table 1 Results of multivariate (PERMANOVA on assemblage structure) and univariate (ANOVA on all other variables) analyses testing for the effects of ‘sites’ and ‘reefs’ on several response variables of fish assemblages

Fig. 3

Variance components of each of ten response variables of fish assemblages at three spatial scales

Functional groups of reef fishes

Analogously, spatial patterns of each functional group of fishes showed the greatest variance at the transect scale, followed by the reef and the site scales (Fig. 3d–j). The largest examined scale provided a null contribution to the variance of the benthic guild (Fig. 3d), the omnivores (Fig. 3f), the browsers of benthic invertebrates (Fig. 3g), and both the fishes laying benthic eggs, or showing a nesting behaviour (Fig. 3i), and pelagic eggs (Fig. 3j). The browsers of benthic invertebrates and fishes deploying benthic eggs were also the only response variables for which the contribution of the among-reefs variance reached at least 50% of that provided by the transect scale (Fig. 3g, i, respectively). A significant variation among reefs, however, was found only for the browsers of benthic invertebrates and the omnivorous fishes (Table 1), although such effect might have been, to some extent, actually masked by the larger smaller-scale variability of all examined responses.

The planktivore and the pelagic group included, respectively, only the ‘big-scale sand smelt’, Atherina boyeri (15 individuals), and the ‘horse mackerel’, Trachurus sp. (one individual), which were all found along a single transect of one reef at site 3. Therefore, these two groups were not analysed individually.

Influence of habitat traits on reef fishes

The ‘best’ model explaining patterns of variation of the whole fish assemblage structure included only the small ‘drops’ as a relevant habitat trait (Table 2). This habitat variable provided a small, but significant, contribution to variation of the assemblage structure. The availability of bare rock and the number of large ‘drops’, instead, were included into the best model regarding the richness of fish taxa and the total number of fish individuals, respectively (Table 2). Both habitat traits were significantly related to the corresponding fish variable (Table 3).

Table 2 Results of AIC_C analysis for competing models assessing the abiotic (sd small ‘drops’, ld large ‘drops’, br bare rock, cr continuous rock, pc pebbles/cobbles, s sand) and biotic (en encrusting algae, er erect algae, f foliose algae, k kelp, t turf) habitat variables (no. var. number of habitat variables for the model) to the variation of different response variables of fish assemblages

Table 3 Results of multiple regression models testing the relationship between habitat traits and different variables of fish assemblages

Concerning individual fish groups, bare rock was also the only habitat variable included in the best model explaining variation of the benthic fishes, the browsers of benthic invertebrates and the fishes laying benthic eggs/having nesting behaviour (Table 2), although its influence was significant for the last two groups only (Table 3). The same set of four variables (large and small ‘drops’, extent of turf-formers as dominant algae, and bare rock) was selected by the best model explaining the variation of the proximo-benthic fishes and of the predators of mobile animals (Table 2), with a cumulative contribution of 20.6 and 23.5%, respectively (Table 3). A significant influence, however, was found only for the number of large drops on both such response variables and for the number of small topographic elements on active predators (Table 3). The variation of the abundance of omnivorous fishes was best explained by the model including only the large foliose algae (Table 2), though the influence of this habitat variable was not significant (Table 3). Finally, the variation of the abundance of fishes releasing pelagic eggs was best explained by the model including foliose algae and sandy bottom as the main available habitats (Table 2). Such habitat variables were both significant and provided a cumulative contribution of 25.5% to the total variation of pelagic spawners (Table 3).

Discussion

The most evident finding of this study is the far larger variance of reef fish distribution at the transect compared to the larger scales. Such pattern was observed consistently across all examined response variables, although a significant among-reefs variability was also documented for the structure and aggregate measures of richness and abundance (total number of taxa and total number of individuals) of the whole assemblage and for the abundance of both fish groups having omnivorous or browser benthic invertebrates feeding habit.

Heterogeneity at small scale (tens of metres) is increasingly emerging as an ubiquitous feature of reef fishes. It was documented in the last decades by studies, in most cases even encompassing a larger extent and more numerous scales than the present one, in New Zealand [30, see also 83 for a review], Australia [29, 37, 39, 71], Aleutian Archipelago [84], Mexican Caribbean coast [85] and the Mediterranean Sea [31, 86]. This can be explained considering that several tropical [87,88,89] and temperate [45, 90,91,92] reef fish species show home ranges that, though changing over their life-cycles, are within a few thousand m² or even confined to less than 10 m².

Spatial variation at the scale of tens of metres can be driven by variation in habitat structure [31, 42, 71, 84] and its possible interaction with recruitment [93] and post-recruitment [94] processes. In this respect, the present finding identifying small or large drops as explaining spatial variation of the whole structure and overall abundance of fish assemblages is consistent with their expected role in increasing the number and diversity of microhabitats. These may be associated to a greater small-scale heterogeneity in the distribution of reef fishes [31, 45, 71, 84, 85]. In principle, a larger availability of microhabitats might also allow relatively more species to coexist in the same area [e.g., 86], but both small and large drops, as well as each dominant algal group, did not explain variation of the total number of taxa. This variable, however, was influenced by the availability of bare rock, potentially providing indirect support to analogous considerations. In fact, the whole set of sampled transects was almost evenly divided into those showing bare rock as the dominant habitat along more than 5% of the 25 m extent, and those where bare rock was less available or completely absent (36 and 39 transects, respectively). On average, a single taxon characterized the first circumstance, while such number was almost doubled in the second one. Of course, such small absolute numbers might just depend on chance rather than any specific ecological mechanism, thus suggesting caution in attempting to generalise ecological mechanisms possibly responsible for spatial patterns of richness of reef fishes. However, they could provide further support to the general importance of habitat heterogeneity in the present system.

Adding to the small-scale variation, previous investigations have detected a considerable proportion of variance in the distribution of reef fishes at scales comparable to that of the present sites. This was often explained with spatial changes of hydrodynamics [95], the relative availability of algae playing different functional roles [30, 84, 96, 97], or a combination of biotic and abiotic processes [32, 98, 99]. Such differences, instead, were unlikely to occur in the present system, where all sampled reefs and sites were chosen as being comparable in terms of wave exposure and general occurrence of macroalgal beds and might have critically contributed to the smaller variation at these scales compared to the transect scale. The only response variables showing a considerable proportion of variance at the reef scale were the abundance of fishes feeding mostly on sessile invertebrates and that of fishes generally having a nesting behaviour. Both life-traits may imply patterns of variation consistent with those of factors and processes affecting the availability of suitable food or protection from predators and harsh environmental conditions. For example, Tuya et al. [44] described relevant variation in the abundance of labrid fishes (the numerically dominant taxon in the present study, collectively providing about 42% of all sampled individuals) as a positive function of among-reefs variation of physical (rock topography) and biotic (cover of red algae) habitat traits. Heterogeneous habitats may facilitate the access to invertebrate preys, at the same time providing enough shelter and protection from predators. Here, the same variables were not found to explain any amount of variation of such fish groups, but the influence of bare rock might suggest, once again, analogous underlying mechanisms. In fact, the three sampled reefs where the browsers of sessile invertebrates and the nesting fishes were more abundant (respectively, 17 and 18 individuals averaged over reef 4 at site 1, reef 4 at site 2 and reef 2 at site 3) were characterized by a proportion of bare rock as dominant habitat about three times lower (10.9 vs. 31.3%) than that of the three reefs where they were less abundant (in both cases, < 1 individual averaged across reef 3 at site 1 and reefs 3 and 4 at site 3).

Contrarily to the initial hypothesis predicting different relevant scales of variation for fish groups with different spatial use of the habitat, an increasingly smaller variance at increasing scale was quantified for all fish groups. Our expectation was based on the role of factors, such as high mobility and long-distance feeding and reproductive behaviour, as direct drivers of patterns of fish distribution. Although this type of relationship can not be unambiguously ruled out by the present findings, the influence of factors potentially determining relevant variation of some fish groups at relatively large scale might be secondary compared to that of complex interactions with within-reef processes. This could occur, in particular, in the present system, where isolated rocky reefs occur interspersed within extensive portions of sandy bottom. Under such circumstances, even species whose life-traits would make minimal small-scale variation expectable might find conditions that are critical for aspects of their life history near or within rocky reefs and not on the surrounding sandy bottom. This could be the case, for example, of the here common ‘European sea bass’, Dicentrarchus labrax, which has high mobility and releases pelagic eggs [61], but may depend on within-reef occurring spawning grounds [100]. In a system, such as the present one, lacking marked among-sites differences in potentially relevant factors, including wave exposure [95], position relatively to the coast [101] and fishing pressure [102], smaller-scale variation might logically predominate irrespectively of fish life-histories relying on obvious local or larger-scale processes. Nonetheless, the general lack of among-sites significant variation should be interpreted with caution since the present impossibility to tease apart spatial and temporal variation likely increased the among-reefs variation, potentially masking, at least in part, among-sites differences. It is worth noting, however, that the Site test remained non-significant also when the non-significant Reef source of variation was eliminated from the analyzed model (e.g., ANOVA on fishes deploying pelagic eggs), suggesting that such spatial homogeneity was due to an actual ecological feature rather than biased by sampling constraints.

Finally, a note of caution is needed when interpreting present findings only in terms of spatial patterns of variation and their possible underlying processes, since logistic difficulties eventually led to a sampling design where the estimated patterns of spatial variance were, at least in part, likely affected by temporal variance. The availability, within each site, of replicate reefs sampled in both the ‘cool’ and the ‘warm’ period reported for the study region [73] allowed to reduce, but not to fully remove, such a confounding effect. Indeed, reef-associated fish assemblages may respond to environmental conditions with changes in abundance not just at different spatial scales, but also at a range of temporal scales [29,30,31,32]. For instance, oceanographic differences can drive seasonal patterns of fish recruitment differing between adjacent islands within the same archipelago [102]. There is evidence, however, that temporal variation of adult reef fishes is much smaller than spatial variation, at least on scales ranging from days and weeks to months, and as long as crepuscular periods are avoided [103, 104]. In spite of the mentioned limitations, the present discussion on spatial patterns of the studied fish assemblages is built on the notion that variability in space (from metres to tens up to 100 km) is often more important than variability in time (weeks to months and even years), especially when response variables are examined at the assemblage level [29, 30, 32, 37, 71].

Conclusion

Although this study did not address any issue specifically related to marine protected areas, it provides basic ecological knowledge that is relevant for the management and conservation of reef fishes. The identification of relevant scales of variation of natural populations and assemblages is essential to design protection schemes suitable to guarantee the protection of their responsible processes [10, 31]. The implementation of local strategies intended to conserve reef fish assemblages, for instance, is likely to fail if their patterns of distribution, abundance and diversity are mainly shaped by larger-scale (e.g., regional) drivers. This kind of problem was highlighted by Benedetti-Cecchi et al. [8] who did not find clear evidence of realized protection on algal and invertebrate assemblages of shallow Mediterranean rocky coasts, likely due to the fact that the relevant scales of variation were not properly represented within the established protected area. At the same time, multi-scale sampling designs able to quantify and compare the scales of variation of ecological variables of interest are required to assess the positive or negative impacts of existing human activities, including already implemented protected areas and fishing [99].

The northern Portuguese coast, in particular, is subjected to small-scale fishing using artisanal gears [105] that may negatively impact reef fishes and habitats more than it is normally assumed [106]. Therefore, the present findings are crucial firstly to assess whether such anthropogenic activities have the potential to alter the relevant scales of variation of reef fishes through interfering with their likely responsible processes, secondly to optimise the design of planned protection schemes. In the context of the well-known SLOSS dichotomy (a single-large vs. several-small options debated in reserve design since the 1970s) [see 107], effective protection to the studied reef fishes might be provided by a set of small reserves, likely preserving the important processes, rather than a single large one, which might be more appropriate for fishes having wide home ranges [e.g., 108]. Interestingly, the same conclusion was drawn by Hattori and Shibuno [109] examining reefs lacking, analogously to the present ones, territorial herbivore fishes. Moreover, present data may help to design monitoring programmes specifically focusing on the relevant spatial scales, thus improving their cost–benefit analysis in logistically challenging systems in terms of sampling conditions.

In conclusion, although further investigations, including experimental manipulations of habitat traits and specific observations of fish behaviour, competition and predation, are needed to unambiguously reveal cause-effect relationships, this study contributes to disclose the multi-scale drivers of the structure of reef fishes in a relatively unstudied region. This is a key goal of fish ecologists and provides fundamental information for the management and conservation of a coastal system hosting fishery- and ecologically important species and habitats. Moreover, present data may provide a reference for future comparisons aimed at assessing the possibility to generalize the same patterns to larger scales and systems.