Contrasting distribution and foraging patterns of herbivorous and detritivorous fishes across multiple habitats in a tropical seascape

Understanding drivers behind patterns of functionally important groups of fishes is crucial for successful management and conservation of tropical seascapes. Herbivorous fishes are the most prominent consumers of marine primary production which can have profound effects on reef resilience. We explored environmental variables affecting distribution and foraging patterns of herbivorous and detritivorous fish assemblages (siganids, acanthurids and parrotfish) across distinct shallow-water habitats (coral reefs, macroalgae beds and seagrass meadows) during September–November 2016 at Mafia Island, Tanzania (8°00 ′ S, 39°41 ′ E). We performed underwater visual census to quantify fish assemblages, measured habitat features, deployed macroalgal assays and conducted inventories of grazing scars. Multi-dimensional scaling and mixed-effects linear models were used to evaluate differences in fish assemblages and environmental variables influencing abundance and foraging patterns of fishes. Fish communities of focal functional groups differed among habitats. Abundance of herbivores and detritivores as well as relative browsing and scraping was highest on coral reefs compared to macroalgae and seagrass meadows. Adult fish were more abundant on coral reefs while juveniles were abundant in macroalgal beds. Coral cover and crustose coralline algal cover had a positive effect on the abundance of fish in coral reef areas, while macroalgal cover had a negative effect. Contrastingly, in macroalgae habitats, macroalgal cover had a positive effect on the abundance of parrotfish. These results highlight the importance of considering connectivity between macroalgal beds and coral reefs through ontogenetic shifts in habitat use by primarily microphagous parrotfish and of incorporating a range of habitats within coastal management plans.


Introduction
Disentangling which habitat characteristics are critical in shaping abundance, distribution, and ecological processes within the available range of coastal habitats will improve our understanding of marine community patterns and connectivity in shallow coastal seascapes. Benthic habitats provide marine organisms with food and shelter from abiotic or biotic stressors, or a combination of these (Friedlander and Parrish 1998), which can further influence competition and predation risk (Almany 2004).
In tropical reef environments, spatial heterogeneity is a strong driver of distribution and diversity patterns of reef fish (Chabanet et al. 1997;Sabater and Tofaeono 2007;Messmer et al. 2011;Samoilys et al. 2018), where characteristics such as structural complexity and rugosity are important predictors of fish assemblage structure and abundance (Bell and Galzin 1984;Gratwicke and Speight 2005;Graham and Nash 2013). In submerged macrophyte habitats, canopy height and cover have been shown to structure fish assemblages in both tropical seagrass meadows (Gullström et al. 2008;Alonso Aller et al. 2014) and macroalgal beds van Lier et al. 2018). However, environmental factors as such may impact coral fish communities on multiple spatial and temporal scales as some fish species might have shifting habitat preferences during different ontogenetic stages and hence utilize alternative habitats like seagrass meadows, mangroves and/or macroalgal beds during their juvenile life stage (Nagelkerken et al. 2000;Wilson et al. 2010;Berkström et al. 2012). For example, trade-offs between foraging success and predation risk depending on body size have been shown to explain habitat choice of juvenile groupers (Dahlgren and Eggleston 2000). Microhabitat preferences might therefore not be consistent across habitats, and environmental predictors might change depending on habitat identity.
Certain functional groups of reef fish have been identified as having a major role in maintaining reef resilience and reef health such as nominally herbivorous and detritivorous fishes (Bellwood et al. 2004;Mumby 2006;Hughes et al. 2007;Goatley and Bellwood 2010). The former group can be further classified into three groups: algivores (i.e., fish that target eukaryotic macroalgae and turfing algae), seagrass-feeders, and microphages (fish feeding on microscopic phototrophs that can be both epilithic and endolithic) (Choat et al. 2002;Clements et al. 2017;Clements and Choat 2018). Drivers predicting and influencing abundance, distribution and feeding patterns of these groups are considered imperative within the management of tropical shallow seascapes. Through feeding activities, these groups are acknowledged to be important and influential forces in shaping and structuring ecosystems through top-down control, especially by balancing algal and coral communities in favour of corals.
Although distribution patterns of herbivorous fish and herbivory have been studied extensively in coral reef habitats (see, e.g., Hay 1981;Bellwood et al. 2014;Adam et al. 2015;Russ et al. 2015;Steneck et al. 2017), fish herbivory can also have consequences in other habitats. In seagrass systems, removal of epiphytes through grazing is hypothesized to enhance seagrass production by decreasing competition for light (van Montfrans et al. 1984), and seagrass browsing species have a major role in trophodynamics as they incorporate carbon and nutrients in the food chain (Unsworth et al. 2007). In temperate macroalgal habitats, spatial patterns of algal distribution have been shown to be structured by top-down control induced by fish (Vergés et al. 2009;Taylor and Schiel 2010). However, foraging patterns of nominally herbivorous fish in naturally occurring tropical macroalgal beds have been largely understudied (but see Lim et al. 2016). Canopy-forming tropical macroalgal beds can harbour large numbers of herbivores and have only recently been acknowledged as important nursery grounds for coral reef fishes, including juvenile acanthurids, parrotfish and rabbitfish (Evans et al. 2014;Tano et al. 2017;Eggertsen et al. 2017;Fulton et al. 2019), highlighting the importance of investigating habitat-specific ecological functions and linkages to other habitats.
A prevailing perspective is that nominally herbivorous fishes are driving benthic structure and habitat differentiation through strong top-down control (Mumby 2006;Jackson et al. 2014;Bonaldo et al. 2014;Adam et al. 2015). As a large part of the studies that propose top-down control on the benthos by herbivorous and detritivorous fishes are performed in the Caribbean, similar mechanisms might or might not be applicable to the same extent in other geographical locations (Roff and Mumby 2012). Instead, bottom-up effects such as benthic characteristics might be more important in other places, considering the strong effects benthos can have on reef fish communities as a whole (Friedlander and Parrish 1998;Friedlander et al. 2003;Messmer et al. 2011). Long-term studies from the Philippines have shown strong bottom-up effects on distribution and abundances of detritivorous fish (Russ et al. 2018) and parrotfish assemblages (Russ et al. 2015), and similar patterns have been observed in the Pacific (Tootell and Steele 2016) and the Indian Ocean . Due to the inconsistencies of responses observed, which might be system-specific, further studies are needed.
In the Western Indian Ocean (WIO) region, fishing has been identified as the strongest factor influencing coral reef fish abundance, although benthic variables also had some effects (McClanahan and Arthur 2001), indicating a mix of top-down and bottom-up factors. The majority of herbivoremacroalgae-coral studies have been performed in Kenya (see, e.g., McClanahan et al. 1999McClanahan and Arthur 2001), where top-down control in the form of herbivory has been identified as an important predictor for macroalgae abundance on coral reefs Mörk et al. 2009). Likewise, a study in seagrass meadows from Zanzibar has shown that both top-down (fishing) and bottom-up (seagrass density) control explained fish assemblage structure (Alonso Aller et al. 2014). However, patterns as such might vary across the WIO region, in different habitat types and with different fishing pressure.
The present study therefore sets out to explore distribution and foraging patterns of herbivorous and detritivorous fishes across multiple habitats within a shallow tropical seascape in a marine protected area (MPA) in Tanzania. Additionally, the extent these patterns are driven by benthic characteristics (bottom-up controlled) is studied. In doing so, we also identify habitat variables which might be important in structuring juvenile communities in nursery/recruitment habitats. This is done by (1) quantifying detritivorous and functional groups of herbivorous fish assemblages in distinct and common shallow-water habitats (coral reefs, macroalgal beds and seagrass meadows), (2) identifying environmental variables structuring these fish assemblages and foraging patterns within the seascape and (3) quantifying browsing and scraping of algae in these habitats. We hypothesized that the observed fish assemblages (abundances and species composition) and foraging patterns would differ across habitats on a broad spatial scale and will be related to differences in habitat characteristics, resource availability, and predatory assemblages at fine spatial scales (1-25 m).

Description of study area
The present study was conducted within the boundaries of the Mafia Island Marine Park (MIMP), in Mafia Island, Tanzania (8°00′S, 39°41′E) (Fig. 1). The Mafia Island area is composed of a small archipelago located 20 km from the Tanzanian mainland and 120 km south of Unguja Island (Zanzibar) and has a high degree of marine biodiversity (Horrill et al. 1996;McClanahan et al. 2008). The MIMP was established in 1995 and is located on the southeastern part of the main island, covering a total area of 822 km 2 , of which 75% is below the high water mark (Gaspare et al. 2015). The protected area is a multi-use park, divided into zones subjected to different degrees of protection, from allowance of subsistence/artisanal fishing in some areas to complete closure in other areas. However, all destructive fishing methods and coral mining practices are forbidden (Berkström et al. 2013;Gaspare et al. 2015).
The archipelagic environment is influenced by the East African Coastal Current (EACC), flowing northwards along the East African coast, and the northeast and southeast monsoon winds with one drier and sunnier period (October-March) and one rainier and more cloudy period (March-October) (Berkström et al. 2013). Tides are semidiurnal with a mean amplitude of 3.3 m and causing strong and complex currents with velocities reaching up to 6 knots (Horrill et al. 1996). The eastern coastline of the Mafia Island is exposed to the open Indian Ocean, but protected by fringing reefs towards the southern part of the archipelago (Garpe and Öhman 2003;Gaspare et al. 2015). A large and relatively shallow bay (< 15 m deep), fringed by mangrove forests is located on the eastern side of the main island (Chole Bay), and shallow areas are also located around and towards the islands of Jibondo, Juani and Chole (Berkström et al. 2013). Shallow areas are characterized by a mosaic of macroalgae-dominated areas, sand flats, seagrass meadows and coral patch reefs (Berkström et al. 2013).

Field study
The field survey was conducted during September-November, 2016. The study coincided with the cooler season in A total of 3 different habitats types were surveyed at 9 sites scattered throughout the shallow seascape (< 4 m of depth) (Fig. 1). Seagrass meadows were monospecific and consisted mainly of Thalassodendron ciliatum and macroalgal areas were dominated by the canopy-forming brown algae Sargassum aquifolium and Turbinaria conoides. While seagrass densities and cover can fluctuate with season (Robbins and Bell 2000), they generally do not undergo such dramatic changes as some Sargassum spp. (Fulton et al. 2014). Many Sargassum species from the southern Western Indian Ocean have slower growth rates during this time of the year, with cover and biomass in algal patches becoming less dense (Gillespie and Critchley 1999a). There is an increase in biomass and cover during the warmer season (December-March) (Gillespie and Critchley 1999b). Coral reef sites were comprised of a heterogenous mosaic of scleractinian corals. Habitat types were chosen on the premises where they were commonly occurring, located at approximately the same depths (1-4 m) and present in large and consistent patches (> 100 m 2 ).

Fish and habitat surveys
Fish assemblages were surveyed by performing underwater visual census (UVCs), (n ≈ 5 habitat and site −1 , n total = 48; n coral = 15, n macroalgae = 20, n seagrass = 13) along 25 × 2 m belt transects by snorkeling. All fishes were identified to lowest taxonomical level and total length (TL) was noted to closest cm, according to the method described in Tano et al. (2017). All surveys were performed between 09:30 and 16:00 and between the low and high tide peak. A snorkeler swam along the transect line ~ 0.1 m s −1 , documenting all mobile fish species and then returning along the line examining substrate and plant/algae cover more thoroughly for small, cryptic species. Each snorkeler was equipped with a camera (Canon Powershot G7x Mark II and Canon WP-DC54 underwater housing), for later identification, if necessary. To avoid potential bias of length estimations, trial estimations were performed prior to the study until snorkelers were calibrated with each other and any possible bias consistent. All UVCs were performed by the same observers (M. Eggertsen and D.H. Chacin).
Every 5th meter of the transect, a 0.5 × 0.5 m square was placed and photographed from above, where benthic composition (percent cover of calcareous rubble, soft substrate, live coral cover, macroalgal cover, seagrass cover, number of macrophyte species, epilithic algal matrix (EAM) cover and crustose coralline algal (CCA) cover and rugosity) was estimated and the height of the 3 tallest macrophytes was measured to the closest cm using a ruler (n = 6 transect −1 ). Rugosity was visually estimated according to a 1-5 grade scale, where 1 denoted no rugosity (flat) and 5 the highest, a method suggested by Gratwicke and Speight (2005).
Depth was measured for each transect with a dive computer (Suunto Vyper Novo) and each transect was georeferenced by marking start and end point with a GPS kept at the surface and placed inside a waterproof bag. For all variables measured, a mean value per transect was calculated.

Browser assays
For estimating relative browsing pressure in different habitats, standardized browsing assays were used. Even though browsing assays only provide a relative estimation of the surrounding grazing pressure, it has been widely used within ecological studies for quantifying grazing pressure across areas and for identifying important browsing species (Fox and Bellwood 2008a, b;Ganesan et al. 2006). However, results have to be interpreted with care, as there is a profound risk that opportunistic fish might feed on assays out of curiosity, an effect which will be impossible to disentangle from data on biomass loss only (Wulff 2017). To circumvent this issue to some degree, 6 of the bioassays were documented by remote underwater video cameras (RUVs) (GoPro Hero Session) for ~ 1 h during each experimental setup (n total = 54), to identify which species of herbivores targeted the tethered algae. Assays were constructed by tethering thalli of Eucheuma denticulatum and Sargassum aquifolium to PVC pipes and placing them in three different habitat types (Thalassodendron ciliatum-dominated seagrass meadows, macroalgal beds, and coral reefs) during a full tidal cycle (~ 24 h). The macroalgal species were chosen because they are a common feature of the macroalgal assemblages in coastal Tanzania (Tano et al. 2016(Tano et al. , 2017, and have been found in relatively large quantities in stomachs of main herbivorous fishes such as Siganus sutor (Eggertsen et al. unpublished), a common browser in Tanzanian coastal waters (Lugendo et al. 2007;Kimirei et al. 2011).
Fronds of E. denticulatum and S. aquifolium were collected from the field and kept in tanks with aerated seawater over night. Suitable sizes of thalli were chosen (< 5 g), dead/ necrotic tissue and epiphytes were removed, and algae were spun in a salad spinner for 10 s to remove excess water. Thalli were then weighed and put into marked zip lock bags for transportation to the site. Two thalli of each species were tethered to a PVC pipe. Every browsing assay constituted a total of 4 pieces of thallus (n browsing assays = 18 habitat −1 , n thallus = 72 habitat −1 , n total = 648) (Fig. 2). Care was taken to place the browsing assays in continuous habitats and at a distance away from habitat edges.
Fishes removing biomass from the browsing assays were identified from the videos and herbivorous fishes within the video frame that were not feeding were also noted. RUVs have in some cases been shown to reveal more information than UVCs regarding observations of herbivorous species (Fox and Bellwood 2008a) whereas in some cases results have been equal (Longo and Floeter 2012). Both methods were included in the current study as RUVs have not been used within the study area before for monitoring herbivores. Browsing assays were retrieved after ~ 24 h, demounted from the PVC pipes and put in marked zip lock bags for transportation to land. Prior to weighing, algae samples were spun in a salad spinner for 10 s to remove excess water. To control for handling and daily growth rates of tethered macroalgae, a caged control (n thallus = 4) was placed in the experimental area at each sampling occasion.

Grazing scar inventories
To quantify grazing activity in different habitats by excavating and scraping parrotfish, an inventory of feeding scar densities was conducted. As feeding by excavating and scraping parrotfish resulted in conspicuous grazing scars in the substratum it is possible to visually quantify these Bellwood 2008, 2009;Bonaldo et al. 2011).
Quadrats of 0.15 × 0.15 m were randomly placed on hard, calcareous (dead) surfaces, where all clean grazing scars were counted and measured to the nearest centimeter (n ~ 20 habitat type −1 ) (Fig. 3). Due to lack of suitable substrate for scrapers to feed on in seagrass habitats, this study was only performed in hard substrate habitats, i.e., coral and macroalgal habitats.
Distribution of focal functional groups of fish in different habitats (coral, macroalgae and seagrass) was explored with non-metrical multidimensional scaling (nMDS) with  Bray-Curtis dissimilarity index from the "vegan" package (Oksanen et al. 2017). Main and interactive effects of habitat and sites for species composition of focal fish communities were tested using PERMANOVA, Bray-Curtis dissimilarities (999 permutations) using "adonis" function from the "vegan" package (Oksanen et al. 2017). Species which occurred ≤ 3 in the UVCs were removed from the nMDS analysis. A similarity of percentage (SIMPER) analysis (Clarke 1993) from the "vegan" package was performed to identify which fish species contributed the most to dissimilarities among the three habitat types. Abundance of functional groups (herbivores, detritivores, and microphages) was analyzed with mixed linear models, "lme4" package, with "site" as random factor and habitat as a fixed factor (Bates et al. 2015). Functional group "territorial grazers" was removed from the analyses as roving herbivorous and detritivorous fishes were the main focus of the present study.
Habitat variables affecting total abundances of all focal functional groups were analyzed separately for each habitat using mixed linear models (multiple regressions) with "site" as a random factor and additional predictor variables as fixed factors (percent cover of calcareous rubble, soft substrate, live coral cover, macroalgal cover, seagrass cover, number of macrophyte species, EAM cover, CCA cover, rugosity and predatory fish abundance). Macroalgal habitats were analyzed with multiple linear regression models, as there was no difference on site level. Model selection was performed by stepwise removal of predictor variables, starting with a model containing all predictor variables allowed. Non-significant variables were then removed one by one and using function 'step' from "lme4" package, based on Aikaike information criterion (AIC) values. If ΔAIC ≥ 2, the model with the lowest AIC value was considered as the most parsimonious one.
Browsing assays (biomass loss for the 2 species of macroalgae) were analyzed with ANOVA, and mixed linear models (total biomass loss in different habitats) with "site" as a random factor and "habitat" as a fixed factor. Mean values of biomass loss for each macroalgal species were calculated per bioassay and the two species of algae were analyzed separately. To identify which habitats differed from each other, Tukey's contrast test was used from package 'multcomp' (Hothorn et al. 2008). Total biomass loss in relation to environmental variables (macroalgal cover, macrophyte cover, rugosity, seagrass and algivorous fish abundances, EAM cover, live coral cover) was tested using mixed linear regression models with "site" and "tethered species ID" as random factors and additional predictor variables as "fixed". Where there was no variance added due to "site" or "tethered species ID", linear regression models were used.
All fish of the focal functional groups < 5 cm were classified as juveniles and larger individuals as adults. Size classes of pooled functional groups of fishes in different habitat types (coral, macroalgae and seagrass) were investigated using mixed linear models, and ANOVA models were performed with abundant juveniles/adults of microphagous fishes as "site" did not add any variation.
Predictor variables were checked for multicollinearity by pairwise comparison using the Spearman rank test and by evaluating variation inflation factor (VIF) values (Zuur et al. 2010). Predictor variables with VIF-values ≥ 2 were removed from the same model. Prior to model fitting, normal distributions of predictor variables were visually examined by basic diagnostic plots and if needed, transformation log(x + 1) and rescaling to size range was performed. All statistical analyses were performed in R version 3.3.1. (R Core Team 2017).

Fish assemblages of focal functional groups
A total of 3,672 fishes from 34 families, 97 genera and 165 species were recorded in the UVCs during the study, of which 1141 were identified as seagrass and algivores, detritivores or microphages. Of these, 645 individuals were identified to species level (27 species), the majority of fish identified only to family level being juvenile labrid scarids. The three habitat types held distinct herbivorous fish assemblages and there was also a significant difference depending on site (PERMANOVA, Fig. 4, Table 1). The algivorous acanthurid Naso brevirostris was one of the most influential species in dissimilarities among habitats as it was abundant in coral reef sites but absent from macroalgal and seagrass habitats (SIMPER, Table 2).
Detritivores were most abundant on coral reef sites and absent from seagrass habitats, while seagrass and algivores and microphages were more evenly distributed across the seascape (Fig. 5).
RUVs revealed a higher diversity of seagrass and algivorous species than the UVCs in all habitat types. For example, N. elegans which were observed feeding on macroalgae in the tethering videos were never observed in the UVCs. Furthermore, few Siganus sutor were observed in UVCs but were common in RUVs (Electronic Supplementary Material; Table 6S).

Environmental variables
Environmental variables identified as food resources targeted by the focal functional groups of fishes (macroalgae, seagrass, EAM, and calcareous rubble) were unevenly distributed between the habitat types (Electronic Supplementary Materials; Table 7S). Calcareous rubble was the most abundant food resource in both coral and macroalgal areas (50.1 ± 5.0 and 52.8 ± 4.6%, respectively). Seagrass habitats had low cover of all food resources except for seagrass (94.1%).
Environmental variables explaining fish abundances differed among the different habitat types (Table 3, Fig. 7, Electronic Supplementary Materials; Fig. 10S). In the coral reef sites, live coral cover and CCA cover had a positive influence on the abundance of all fish (herbivores, detritivores, and microphages pooled) ( Table 3, Fig. 7, Electronic Supplementary materials; Fig. 10S). Likewise, EAM cover had a positive influence on fish abundance in coral reef habitats, both independently and in combination with CCA cover. Within the macroalgal beds, there were significant positive relationships between fish abundances and macroalgal height, macroalgal cover and rugosity (  Fig. 10S). There was also a significant interaction between macroalgal height and rugosity (Linear regression, r2 = 0.43, F (3,16)=5.703, P = 0.019). In seagrass habitats, there were no significant relationships between fish abundances and any of the measured environmental variables.
When separating fish into functional groups, slightly different linear models explained abundances in different habitat types (Table 4, Electronic Supplementary Materials; Fig. 11S). In coral reef habitats, no model was found for seagrass and algivores and detritivores, and for microphages the best model included macroalgal cover, live coral cover, predatory fish abundance and CCA cover (Table 4). In the macroalgal areas, the most parsimonious model for  were recorded on Sargassum aquifolium at the termination of the experiment.

Discussion
To our knowledge, this is the first study examining habitat associations, distribution, and foraging patterns of seagrass and algivorous, detritivorous, and microphagous fishes from varying life stages across several shallow-water habitat types (coral reefs, macroalgae beds, and seagrass meadows) within a tropical seascape. As hypothesized, densities of focal functional groups of fishes and foraging patterns were not equally distributed within the seascape, but varied with habitats. Likewise, species composition of nominally herbivorous and detritivorous fishes also changed significantly with habitat, suggesting that fish communities in this study are subjected to, and to a certain extent, driven by bottom-up processes, which was most obvious for the microphagous parrotfish. explaining abundance of seagrass and algivores was predatory fish abundance and for microphages, macroalgal cover and macrophyte height (Table 4). In the seagrass habitats, there were no significant models.

Browser assays
Biomass loss of tethered algae was significantly higher in the coral reef habitats than in the two macrophyte habitats (mixed linear model, F = 25.454, P < 0.0001, Table 5), with biomass loss being slightly higher in the macroalgal habitats compared to the seagrass meadows. There was no difference in biomass loss between the two macroalgal species. Macrophyte cover (both seagrass and fleshy macroalgae) in the surrounding habitat was the single significant factor explaining biomass loss of tethered macroalgae (linear regression, r2 = 0.436, F (1,16)=12.371, P =0.0028) (Fig. 8).
Very few fish were observed feeding on the tethered macroalgae in the RUVs, and all observations were from the coral reef sites (Chawe Kubwa, Kitutia and Milimani). Acanthurids were the most common fish observed; small (~ 25 cm) Naso brevirostris (n = 4), adult N. elegans (n = 2), N. vlamingii (n = 1) and one siganid species; Siganus sutor (n = 1). All observed fish targeted E. denticulatum, but bites  Life stages of microphagous parrotfish also varied with habitat; smaller individuals (< 5 cm) were more abundant in macroalgal habitats and larger individuals more common in coral reef and seagrass habitats. This was mirrored in the inventory of grazing scars, which were significantly smaller in macroalgal habitats compared to coral reef sites, likely reflecting the distribution patterns of ontogenetic stages of this functional group. Lower numbers of grazing scars in the macroalgal habitats might be the result of fishes feeding more intensely on epiphytes than on hard substrate surfaces. However, that would not explain the significant differences in sizes of grazing scars. Whether or not it was the same species using different habitats during distinct life stages was impossible to tell, as a majority of the juveniles were seldom identified to species level. Tano et al. (2017) also observed large numbers of juvenile parrotfish in macroalgal areas in Zanzibar, Tanzania, compared to neighboring habitats, suggesting that shallow areas in the Western Indian Ocean (WIO) dominated by canopy-forming macroalgae may serve as nurseries for this particular group of fishes.
Habitat features benefiting certain functional groups of fishes in one habitat might not be the same in another since different environmental variables explained fish abundances in different habitat types. "Live coral cover", "CCA cover" and "CCA cover" in combination with "EAM cover", had a positive effect on the abundance of fish (all functional groups pooled) in coral reef habitats, which is similar to results from other studies (Friedlander et al. 2003;Osuka et al. 2018). In the macroalgal beds, habitat quality variables such as "macrophyte height", "macrophyte cover" and "rugosity", had a positive effect, patterns that are consistent with studies of microhabitat selection of labrid fishes in both tropical (Lim et al. 2016) and temperate macroalgal areas van Lier et al. 2017). "Calcareous rubble", which has been shown to have a strong positive effect on abundances of acanthurids and parrotfish (Russ et al. 2015(Russ et al. , 2018, did not have any effects on abundance in the present study. This is probably due to the rather high cover of calcareous rubble in coral and macroalgal habitats (50.3 and 52.8%, respectively), implying that this is currently not a food-limiting resource on Mafia Island. The positive relationship between "live coral cover" and nominally herbivorous fish abundances is probably due to the numerous feeding surfaces in the form of dead hard substrata this benthic category is associated with. Note that "live coral cover" was almost never the dominating substrate in UVCs (48.465 ± 5.26%), not even in the coral reef habitats, which all had a high degree of calcareous rubble. Although not a food source for the majority of the fishes in the study, CCA cover might indicate areas with low sedimentation, which are more attractive feeding grounds for many herbivorous fishes (Bellwood and Fulton 2008).
The lack of significant explanatory variables in the seagrass habitats is likely a consequence of the limited sample size and small value range within measured variables. Factors such as spatial arrangement within the seascape and distance to other habitats or deeper waters might be more important for structuring fish assemblages in these habitats (Gullström et al. 2008;Henderson et al. 2017). Fish likely respond differently to habitat characteristics and environmental variables depending on ontogenetic stage (Macpherson 1998;Almany 2004). The differing bottom-up variables explaining fish abundance in different habitat types in the present study are probably reflecting shifts in habitat preferences for different life stages of fishes. Abundance of microphagous parrotfish was negatively associated with macroalgal cover in coral reef habitats, while in macroalgal areas the relationship was reversed and there was instead a positive effect of macroalgal cover and canopy height on parrotfish abundance. This is likely an effect of smaller/younger fish seeking shelter from predators in the form of structural complexity, provided by height and cover of macroalgae, and by hard underlying structure (rugosity) . The negative relationship with macroalgal cover is in line with results from previous studies from the Great Barrier Reef, Australia, where dense macroalgal stands of Sargassum spp. in coral reef environments were found to induce avoidance behavior of roving herbivores (Hoey and Bellwood 2011). There is likely a trade-off between food and shelter changing throughout ontogeny in these fishes, but also between different types of shelter, provided by live coral or by microhabitat topography in the macroalgal areas. Shelter from predators is an important resource in coral reef environments (Kerry and Bellwood 2016), and can be sizespecific, as inter-structural spaces which are smaller than the potential predator width will act as shelter for a fish that fit within the space provided (Bartholomew et al. 2000;Gullström et al. 2011). Macroalgal areas with both soft and hard complexity in the form of fronds and coral rubble might therefore provide sufficient shelter for smaller individuals, while it is advantageous for larger fish to keep to areas with larger structures, as provided by live coral. Consequently, macroalgal areas in shallow areas with tall canopies and high macrophyte cover might be advantageous for the survival and recruitment of parrotfish and hence contribute to the replenishment of important herbivores and microphages on coral reefs. This is strengthened by recent studies showing that tropical macroalgal habitats host high numbers of juvenile and subadult parrotfish, suggesting that they are important nursery grounds for these species and that linkages between habitats should be highlighted in coral reef management (Tano et al. 2017;Eggertsen et al. 2017). However, we would like to stress that deteriorated coral reefs with high cover of fleshy macroalgae is not what we define as "natural macroalgal areas", and that nursery habitat qualities might depend on where in the seascape a habitat is located. This is, to some extent, illustrated in the present study as high macroalgal cover in the coral reef habitat had a negative effect on the abundance of parrotfish while it was positive in another habitat type.
Algal consumption, measured as the removal of macroalgae from assays, was the highest within coral reef sites, followed by macroalgal beds, and seagrass meadows. However, it does not necessarily mean that grazing intensity per se was the highest in coral reef habitats in the present study. Since biomass loss of tethered macroalgae was negatively related to macroalgal cover in the surrounding habitat, it suggests a density-dependent process, where grazing pressure might be "smoothed out" in areas with high macrophyte cover. Macroalgal feeding fish might therefore not be able to eradicate dense macroalgal stands or "blooms" on coral reefs, if they are already established. Solitary fronds may, however, be targeted and limited by browsing fish since the survival and growth of newly settled Sargassum spp. have been shown to be strongly suppressed by top-down control (Diaz-Pulido and McCook 2003).
Understanding how habitats structure fish assemblages in spatially heterogeneous environments is imperative because it will improve habitat-based approaches of management and conservation. As tropical natural macroalgal habitats have been shown to host high numbers of juvenile and subadult scarine labrids, these linkages between reefs and macroalgal beds should be considered in seascape management (Tano et al. 2017;Eggertsen et al. 2017;Fulton et al. 2019), even though these habitats have only been recently acknowledged. Furthermore, the availability of nursery grounds and connectivity to coral reefs by nominally herbivorous fishes have been found to enhance herbivory (Adam et al. 2011;Olds et al. 2012, Harborne et al. 2016. Therefore, the protection of important ecological functions urges the preservation of seascapes comprising multiple habitats with a high degree of connectedness. This is particularly important in the WIO area, where a large part of the population is living within the coastal region and is heavily dependent on the local marine resources.

Conclusion
The current study illustrates that multiple habitat types within the seascape matrix are needed to host a diverse community of herbivorous and detritivorous fishes. Further, the study suggests that in the absence of overfishing, the herbivorous community may be driven by bottom-up mechanisms, with emphasis on microphagous parrotfish. The study also highlights the importance of habitat quality in terms of cover and complexity for maintaining certain ecological functions. These features might vary across habitats, and affect life stages of nominally herbivorous fishes differently. A habitat feature (e.g., macroalgal cover on coral reefs) that is detrimental for a certain age class of a functional group of fishes (e.g., adult microphages) might simultaneously be important in an adjacent habitat (e.g., natural macroalgal beds) for a different life stage of a fish (juvenile) of the same functional group. This highlights the need of a holistic view of the shallow tropical seascape for efficient management and conservation.