Our study is the first to use functional approach in a temporal scale for reef benthic communities. Our results show that Rocas Atoll benthic community is temporally dominated by functional entities responsible for providing low reef tridimensional complexity (e.g. turf algae, cyanobacteria and Zoanthus sociatus), high tridimensional complexity and primary production (e.g. macroalgae as Sargassum spp., Dictyopteris spp.), with low contribution of carbonate accretion (e.g. crustose coralline algae and corals as Siderastrea stellata, Porites astreoides and Favia gravida) and nutrient cycling (e.g. sponges and/or ascidians) functional entities. This dominance reflected in low values of functional entropy (FEnt), due to uneven abundances distribution between unique functional entities, that is, those that have unique trait combination and are in functional space vertices. Functional richness (FRic) oscillated over years, but did not show great changes in functional spaces, maintaining an equity in the number of functional entities and indicating stability of reef functions in Rocas Atoll, even with unequal abundances’ distribution.
Diversity patterns in Brazilian oceanic islands are well known for reef fishes assemblages and barely known for benthic organisms (Floeter and Gasparini 2000; Floeter et al. 2001, 2008; Pinheiro et al. 2017, Aued et al. 2018). Marine oceanic communities generally present low species richness when compared to coastal communities, specially due to their isolation and different oceanographic conditions (Pinheiro et al. 2017), and functional richness (FRic) usually responds to species richness pattern (Whittaker et al. 2014; Jacquet et al. 2017; Wong et al. 2018). Our results indicate that organisms’ richness did not affect the number of functional entities and FRic values, which presented small variations through the years. In 2013, when we identified 36 morphospecies, and in 2019 - the lowest morphospecies richness -, FRic was the same in our temporal series. If compared to other regions in the Brazilian province - mainly those located at mid-latitudes, but also with Fernando de Noronha Archipelago, other oceanic island located less than 100 miles from Rocas Atoll -, the number of taxa and FRic are both low (Aued et al. 2018; Aued et al. in prep).
Low tridimensional complexity functional entities are composed by benthic organisms from different groups, such as epilithic algal matrix, articulate coralline algae, cyanobacteria and zoanthids. Different from entities that are bioconstructors or provide high tridimensional complexity, they form widespread patches and occupy larger areas on the substrate, with small vertical growth. Epilithic algal matrix (EAM or turf algae) is the main component of reef benthic communities among Brazilian reefs, including Rocas Atoll (Longo et al. 2015; Aued et al. 2018). Our results showed higher abundances of low tridimensional complexity functional entities every year, indicating that the community is also dominated temporally by turf. On EAM composition, it is possible to find a variety of algal groups (e.g. diatoms, cyanobacteria, Chlorophyta, Rodophyta), associated cryptofauna and sediment (Connell et al. 2014). Therefore, EAM contribute to sediment trap, primary productivity, and provide habitat and food for other marine organisms (Vroom et al. 2006; Dijkstra et al. 2017), but offer smaller tridimensional complexity. Therefore, due to its high percent cover, the functional entities composed by EAM are also the main responsible for primary production in Rocas Atoll.
High tridimensional complexity in marine environments, such as tropical forests in terrestrial environments, provides habitat for many marine species, promotes interactions between species and the environment (Marx and Herrkind 1985; Williams et al. 2002; Graham and Nash 2013; Richardson et al. 2017; Fulton et al. 2019; Lamy et al. 2020), and also influences reef fish recruitment and development (Evans et al. 2014; Fontoura et al. 2019). This function is generally associated to branching corals in the Caribbean and Indo-Pacific sites (Alvarez-Filip et al. 2011; McWilliam et al. 2018; Richardson et al. 2020), but for Rocas Atoll, due to the absence of branching bioconstructors, the ‘high tridimensional complexity´ is associated to macroalgae. Opposed to massive corals, like Siderastrea stellata, or benthic organisms that form widespread patches, like EAM and zoanthids, macroalgae can act as refuge from larger predators and affect the local diversity (Ware et al. 2019). Besides providing local tridimensional complexity, primary production is also a function performed by macroalgae and EAM (Brandl et al. 2019; Fulton et al. 2019). Primary producers usually dominate populated atolls and regions with high human disturbances (Sandin et al. 2008). Conversely, for Rocas Atoll, where anthropic impacts are controlled, primary production is associated to ten of the 17 functional entities found in our study. Entities associated to low and high tridimensional complexity occupy together around 80% of the substrate in every sampling year. Due to their dominance over the community and to their extreme trait combination, FEnt was low at Rocas Atoll every year, but still the small oscillations in FEnt occurred as a function of the oscillations of functional entities percent cover, indicating that these functional entities are substantially contributing to ecosystem functioning (D’Agata et al. 2016).
Both carbonate accretion (bioconstructors) and nutrient cycling (suspention/filder feeders) functional entities, have low percent cover at Brazilian Province when compared to the other groups (Aued et al. 2018), and our results showed that this has been a consistent pattern at Rocas Atoll through time. The distribution and abundance of bioconstructors can create habitat heterogeneity and provide the main framework of reef ecosystems (Perry and Alvarez-Filip 2018). FRic revealed four functional entities responsible by calcium carbonate accretion across time, but due to its unique trait combination and low percent cover, they might be considered as vulnerable functional entities. Crustose coralline algae (CCA) is the main reef builder in Rocas Atoll (Kikuchi and Leão 1997), but during the years of sampling it presented low and oscillating percent cover. This pattern may be related to a potential limitation of the photoquadrat method used (Preskitt et al. 2004). The same occurred to the functional entities composed by corals that act as massive bioconstructors. Still, the maintenance of massive bioconstructors percent cover indicates that calcium carbonate accretion function is stable over the years. Sponges and/or ascidians (suspension/filter feeder functional entities) are the main responsible by nutrient uptake and cycling, and bioerosion in reef ecosystems (Bell 2008; Brandl et al. 2019). Although suspension/filter feeders were persistent, their percent cover was extremely low. Even providing important functional roles in the ecosystem (Bell 2008; De Goeij et al. 2013), due to its low percent cover, suspension/filter feeders might not be considerably influencing the nutrient cycling in Rocas Atoll shallow reefs.
Functional richness and functional entropy are indices that represent the actual state of communities functional structure (Galland et al. 2020). Therefore, applying them in temporal scales might act as a baseline for future studies. Our results indicated a temporal stability in Rocas Atoll benthic communities, in the opposite way of many places around the world, where not just benthic, but also fish communities, are eroding (Jackson 2001; Gardner et al. 2003; Schutte et al. 2010; Alvarez-Filip et al. 2011; Cruz et al. 2015a; Smith et al. 2016; de Bakker et al. 2017; Leggat et al. 2019; Muñiz-Castillo et al. 2019). Changes in FRic over the years occurred when two specific functional entities, a massive bioconstructor and a suspension/filter feeder, were found in our analysis. Higher FRic values were in 2015 and 2017, when the massive bioconstructors functional entity appeared in one photoquadrat in both years, while in 2018 – the second higher FRic -, the suspension/filter feeder was detected. In the other years, when both functional entities were not found in our analysis, FRic in Rocas Atoll presented the lowest values, but no ecosystem function was lost.
FEnt was low every year due to uneven abundances distribution between functional entities with unique trait combination (i.e. those located at functional space vertices). The uneven abundance distribution over the years is conduced manly by the dominance of functional entities responsible for low and high tridimensional complexity, and low percent cover of other groups as bioconstructors and suspension/filter-feeders, a pattern already observed for Brazilian Province and Rocas Atoll (Longo et al. 2015; Aued et al. 2018). Oscillations in FEnt are, therefore, related to these dominant groups and their ephemeral oscillations, given that they are groups that respond faster to different environmental conditions (Cheroske et al. 2000; Vroom et al. 2006). By using both FRic and FEnt, it is possible to have a notion of communities functional redundancy (Galland et al. 2020). Functional redundancy buffers ecosystems functions loss against species loss in a community, ensuring stability and resilience to communities structures and functions (Bender et al. 2017; Biggs et al. 2020). FEnt low values and high FRic over the years might be indicating low functional redundancy, that is, there are many functional entities with unique trait combination with low percent cover in Rocas Atoll benthic community. If somehow these functional entities and morphospecies go extinct locally, considering all impacts that might happen with climate change over benthic communities (i.e. phase shifts, reef flattening) (Alvarez-Filip et al. 2011; de Bakker et al. 2017; McWilliam et al. 2020), there will be no other functional entities able to substitute the same function in the functional space and, consequentially, there will be erosion in ecosystems functions (Wong et al. 2018).
Species occurrence and abundances can vary in different ways in a community, depending on natural changes of environmental conditions such as light penetration, temperature and sediment (Putman and Wratten 1984; Perry and Alvarez-Filip 2018). Consequently, changes on environmental conditions can cause impacts in species traits and abundances and, therefore, in ecosystem functions (Wong et al. 2018; Zawada et al. 2019a). Being an oceanic island, Rocas Atoll is a very dynamic system (Longo et al. 2015), with intense water circulation and wave action that can influence the availability and sediment dynamics (Hearn et al. 2001). For example, depending on the quantity of sediment, it can cause shading and hence, affect primary production (Riul et al. 2008), also it can promote the temporal burrowing and smothering of coral colonies (Brown et al. 2002). This probably explains differences in percent cover of massive bioconstructors functional entities. The intense hydrodynamic at shallow tide pools can also favor the presence of more resistant functional entities (Steneck and Dethier 1994; Cheroske et al. 2000), such as the ones composed by epilithic algal matrix, macroalgae and crustose coralline algae.
Considering the more common and frequent global impacts caused by ocean warming and thermal anomalies, and other emerging impacts of the Anthropocene (Hughes et al. 2017, 2018b; Lough et al. 2018; Perry and Alvarez-Filip 2018; Williams et al. 2019), there is a chance of Rocas Atoll be impacted either. In view of these changes around the world, it is more than necessary to enhance our knowledge of ecosystems actual state to serve as baselines for future studies (Christensen et al. 1996; Sandin et al. 2008; Downs et al. 2011; Gatti et al. 2015; Chollett and Robertson 2020), as well as to keep monitoring these communities structures and functions (Alvarez-Filip et al. 2011; Cruz et al. 2015b; Heron et al. 2016; Smith et al. 2016; de Bakker et al. 2017). Our study, inserted in Long Term Ecological Research of Brazilian Oceanic Islands, therefore acts as a baseline for Brazilian Province, since it provides the actual state of reef benthic communities using a functional approach, in an environment with no direct anthropic impacts.
Wilderness areas as Rocas Atoll support unique ecological values, different from the ones near the coast, due to well-managed efforts and limited human pressure, being the only environment considered as near pristine in Brazil (Longo et al. 2015). In our research, due to the maintenance in number of functional entities, small oscillations in FRic and FEnt, ecosystem functions do not seem to vary substantially, indicating a stable state dominated by primary producers and tridimensional complexity functional entities over the years. However, we speculate that the changes in the Anthropocene may result in gains in some functions (e.g. low complexity) and loss of other functions (e.g. carbonate accretion), and hence cause erosions in FRic and even smaller FEnt values. Indeed, a strong bleaching event has stricken Rocas Atoll bioconstructors recently (Gaspar et al. 2021). Therefore, reef communities are expected to show changes in the next years. Although COVID-19 pandemic and political dismantling of science and institutions in Brazil are impairing our monitoring, the continuity of temporal sampling is needed to understand how local and global variables are affecting benthic ecosystem functions even at pristine places, like Rocas Atoll. This can help to predict the effects on some ecosystem function (e.g. complexity) caused by global changes and it´s consequence for the ecosystem services.