Spatial heterogeneity of coral reef benthic communities in Kenya

Spatial patterns of coral reef benthic communities vary across a range of broad-scale biogeographical levels to fine-scale local habitat conditions. This study described spatial patterns of coral reef benthic communities spanning across the 536-km coast of Kenya. Thirty-eight reef sites representing different geographical zones within an array of habitats and management levels were assessed by benthic cover, coral genera and coral colony size classes. Three geographical zones were identified along the latitudinal gradient based on their benthic community composition. Hard coral dominated the three zones with highest cover in the south and Porites being the most abundant genus. Almost all 15 benthic variables differed significantly between geographical zones. The interaction of habitat factors and management levels created a localised pattern within each zone. Four habitats were identified based on their similarity in benthic community composition; 1. Deep-Exposed Patch reef in Reserve areas (DEPR), 2. Deep-Exposed Fringing reefs in Unprotected areas (DEFU), 3. Shallow Fringing and Lagoon reefs in Protected and Reserve areas (SFLPR) and 4. Shallow Patch and Channel reefs (SPC). DEPR was found in the north zone only and its benthic community was predominantly crustose coralline algae. DEFU was found in central and south zones mainly dominated by soft corals, Acropora, Montipora, juvenile corals and small colonies of adult corals. SFLPR was dominated by macroalgae and turf algae and was found in north and central zones. SPC was found across all geographical zones with a benthic community dominated by hard corals of mostly large colonies of Porites and Echinopora. The north zone exhibits habitat types that support resistance properties, the south supports recovery processes and central zone acts as an ecological corridor between zones. Identifying habitats with different roles in reef resilience is useful information for marine spatial planning and supports the process of designing effective marine protected areas.


Introduction
Coral reefs are among the most biodiverse ecosystems in the world [1], providing an array of ecological services which are important for human well-being [2,3]. However, coral reefs are declining at an alarming rate globally with climate change posing as one of the biggest threats [4]. There are still limited global efforts to abate climate change [5] making reefs a highly threatened ecosystem with climate scenarios showing 99% of the reef could disappear in this century [6]. The persistence of coral reefs will rely heavily on their ability to maintain a coraldominated state and avoid shift to algal-dominated or other alternative stable states amid the inevitable effects of climate change [7][8][9][10][11].
The Kenya coast stretches along 536 km, between latitudes 1 0 and 5 0 S with narrow fringing reefs in the southern part and patchy reefs with low reef development in the north [12,13], (Fig 1). The distribution of coral species in Kenyan reefs is influenced by the large-scale current dynamics with the East African Coastal Current (EACC) bringing coral larvae from the southern 'center of diversity' for the Western Indian Ocean (WIO) region [14]. A cold-water system prevails in northern Kenya due to the convergence of the EACC with the seasonal Somali Current (SC) that is characterized by poor water conditions for reef development [15,12]. The interaction of the EACC and the SC in the north creates a marginal, high-latitude and upwelling system with transitioning communities from the East African to Somali-Arabian fauna [16]. This results in high coral species diversity in the southern parts of Kenya and a decrease in diversity towards the north [16]. In addition, the presence of river systems in the central-northern region introduces small-scale influences in species distribution by creating environmental barriers that further limit the transport of larvae to the north [17].
There has been a substantial amount of coral reef studies in Kenya focusing on a range of issues including: the effect of MPAs in protecting habitats and biodiversity [18][19][20], impact of anthropogenic stressors, including the 1998/99 mass coral bleaching event, on coral and fish communities [21][22][23][24][25]; but without reference to spatial patterns of communities along the latitudinal gradient. These studies, which include coral reef monitoring with different levels of management programs (marine reserve vs. marine park), have been running for more than three decades showing long-term trends in coral reef communities [26,2,27] but have mostly focused on shallow lagoons within the south coast, leaving deeper forereefs and the north zone depauperate of studies. In addition, there is limited information on coral reef communities on reefs that are remote and not within conservation areas along the Kenyan coast. Knowing the spatial patterns of the coral reef benthic community can help identify habitats that offer resistance/protection to, and recovery after a disturbance including climate change and could help the persistence of coral reefs.
In this study, we follow a quantitative approach to characterize coral reef benthic communities along the Kenyan coast. Using extensive data collected at different geographical zones, reef habitats and management levels multivariate analysis was utilized to investigate the effects of the latitudinal gradient, local habitat and management in structuring spatial patterns in benthic communities.

Study area
Kenya is located on the east coast of Africa between Somalia and Tanzania with a coastline running approximately 500km long between 1 0 and 5 0 S (Fig 1), [12]. Two monsoons namely, southwest (December to March) and northeast (May to October) highly influence the climate and seasonality of this area. The Kenyan coast also experiences the EACC originating from Tanzania and flows northwards throughout the year. In addition, the SC mainly driven by northeast monsoon (NEM) counteracts the EACC when the wind blows southwards. Interaction of these two currents during the NEM causes cooling in the northern part of Kenya, and through mixing of nutrients makes this zone highly productive [28].
Coral reefs in Kenya are characterized by narrow fringing reefs in the southern part and patchy and discontinuous fringing reefs in the northern part [12]. These reefs are found at different depths ranging from less than 1 m at low tide to over 20 m. There are different geomorphological characteristics of the reefs with channel reefs in highly sedimented areas of river inflows to ocean, lagoon reefs mostly occurring at leeward side of reef flat and fringing reefs on the seaward side and patch reefs occurring at offshore outcrops of reefs.
Sea surface temperature (SST) for two 'normal years' with no influence from extreme climate events (e.g.the Indian Ocean Dipole) in 2013 and 2014 show that the north zone had a low mean SST (27.45 0 C) but with the highest variation with a standard deviation of 1.23 0 C, a maximum of 30.16 0 C and minimum of 25.27 0 C. The south zone had the highest mean SST of 27.63 0 C and the lowest variation (sd = 1.15 0 C, max = 29.71 0 C, min = 25.71 0 C). The central zone had the lowest mean SST (27.26 0 C) and an intermediate variation (sd = 1.21 0 C, max = 29.58 0 C, min = 25.41 0 C). SST data were derived from MODIS satellite source by plotting georeferenced surveyed sites on a GIS system. Using ArcGis 10.8, points were created along vector lines at two pixels apart (one pixel = 4x4km) for each zone (north = 10points, central = 5points and south = 3points) and SST values were averaged over the two years. There are different levels of management for coral reefs in Kenya; 'marine parks' receive legal protection from fishing activities and collection of shells or corals. 'Marine reserves' receive partial protection with only traditional fishing methods allowed. 'Unprotected' areas are only regulated by fishing licensing and recently some areas receive fishing restrictions set by local communities. The southern zone has 2 parks and 2 reserves, and the central zone 1 park and 1 reserve in south zone. The rest of the reef is within unprotected areas. The northern part has only 1 reserve and the rest are unprotected reef areas.

Sampling design
A hierarchical sampling design (S1 Table) was used to collect benthic community data across Kenya coast. There were three pre-defined geographical zones (North, Central, South). In each geographical zone, sampling sites were selected to reflect different habitat factors: two depth levels (deep [6.5-18m], shallow [0-6m]), two levels of exposure to oceanic waves (exposed, sheltered) and four different reef types (channel, lagoon, fringing, patch). Additionally, three different management states were sampled; park (no-extraction), reserve (regulated extraction) and unprotected (extraction allowed). It should be noted that not all habitat types and management levels were found at the three geographical zones, for example the north zone had sites in reserve and unprotected levels but there is no park in this zone (S1 Table). Furthermore, the north zone had all the four levels of reef type while central and south zones did not. In the end, there were twenty sites in the north, ten sites in central and eight sites in the south. Choice of these study sites was based on historical monitoring of reef resilience [29], with an additional fourteen sites selected haphazardly and stratified within each geographical zone in order to cover habitats that had not been studied before.

Data collection
Surveys for benthic communities were conducted during the north-east monsoon season during 2015 and 2016. A modified resilience assessment protocol [29] was used to record benthic cover and coral size class distributions. Benthic cover was recorded using photoquadrats with an underwater camera positioned at approximately 1 m above the seafloor with the aid of a 1 m PVC stick. Within each belt-transect, twenty five 50 x 50 cm photos were taken at an interval of 1 m. Cover was recorded in fifteen categories; hard corals (HC), coralline algae (CA), macroalgae (MA), turf algae (TA), soft corals (SC), Halimeda (HA), bare substrate (BS), dead standing coral (DSC), recently dead coral (RDC), invertebrates (INV), rubble (RB), sand (SD), seagrass (SG), silt (SL) and others (OT). HA was separated from macroalgae group as it is a calcifying algae that contributes to the production of carbonate and formation of reef sediments [30,31]. Living coral cover (LCC) was recorded by coral genera.
The size class distribution of coral colonies was measured for 23 selected genera of adult corals using the belt transects [29], and for all juvenile corals by subsampling within the belt transects in six 1m 2 quadrats that were 5 m apart along the 25m belt. Adult colonies were defined as those larger than 10cm in diameter and grouped into six size classes (11-20 cm, 21-40 cm, 41-80 cm, 81-160 cm, 161-320 cm and >320 cm). Juvenile corals were defined as those smaller than 10 cm in diameter and were classified into three size classes (1-2.5 cm, 2.6-5 cm and 6-10 cm).
Benthic photos were analysed using CPCe v4.0 software [32] by randomly placing 25 points on each photo. This number of points overlaid on the photo are within the optimum number required for accurately estimating coral cover in a coral reef that has about 30% cover [33]. The type of benthic substrate was identified under each point.

Data analysis
Three datasets were analyzed: 1. benthic cover by major categories (S1 Table), 2. benthic cover by coral genera (S2 Table) and 3. coral size classes (S3 Table). A set of criteria was established to select only the key variables important for determining the spatial pattern as follows. First, only biotic variables from the major benthic categories were considered as these are influenced by latitudinal gradient and habitat environmental conditions. Only coral genera that had a > 1% cover were included in the analysis. In the second stage, analysis was done in PRIMER v7 [34] where variables were square root transformed in order to reduce skewness and make the variances more homogenous. Each of the three datasets were examined by principle coordinate analysis (PCoA) in order to select those key variables that only contributed to the clustering of sites. Only those variables whose vector direction demonstrated influencing on clustering of sites were selected for further analyses. Variables whose vectors had similar direction were either pooled together or representative variables were selected to avoid analysing many variables that otherwise have similar effect on clustering sites. A permutational multivariate analysis of variance (PERMANOVA) [35], was conducted for each of the three datasets to determine if geographic zone had an effect on site clustering. The selected key variables from the three datasets were then pooled together and collectively analysed. A One-way ANOVA was done in excel to determine how each key variable differed across the geographic zones. PERMANOVA was done in R using the Adonis function within the vegan package [36] to test the effect of depth, reef type, exposure to oceanic waves, management and their interaction. Since geographic zone was found to have a significant effect on variables, this analysis was treated as a block design where all other factors were randomly assigned within each block (i.e. geographic zone).

Selection of key variables
Three PCoAs performed on the ten biotic benthic categories, 13 coral genera and ten coral size classes showed clustering of sites (Fig 2A-2C). Of the ten benthic categories, six (hard corals, coralline algae, turf algae, macroalgae, soft corals and dead standing corals), were selected for further analyses (Fig 2A). Halimeda was not included because their vectors pointed in the same direction, so only CCA was selected as the representative variable. Seagrass and 'otherinvertebrate' categories were not included as key variables because they had very low mean percentage cover. Only four representative genera, Acropora, Montipora, Echinopora and Porites, were considered as key variables for further analysis based on the different sets of similareffect vectors on the PCoA (Fig 2B) as well as their abundance (S6 Table). Vectors of Acropora and Montipora clearly pointed in opposite direction showing that they had distinct influence on site clustering. Vectors of Porites and Echinopora pointed in similar direction and most orthogonal to Acropora and Montipora, and were among the most abundant genera recorded in this study, so were selected as indicative of other genera. Vectors for 1-10 cm size classes were closely grouped, as well as vectors for > 40 cm size classes (Fig 2C), so a representative size class from each was selected. The other size-classes (11-20cm and 20-40cm) were considered separately based on their independence in clustering sites.
A PERMANOVA test on the effect of geography for each of the datasets showed significant differences between zones (S8 Table) and therefore further analysis was restricted to within geographic zones (Table 1).
A total of 64,850 coral colonies were recorded with an average density of 79 colonies per 100m 2 across 4 size classes. The density of corals (number of colonies per 100m 2 ) in 1-10cm size class differed significantly (F 2,88 = 5.0, p = 0.01; Table 2) with the highest density in central and south zones (510.8 ± 333.0 colonies and 510.9 ± 352.8 colonies, respectively) (Fig 3). The density of 11-20cm size class was significantly higher in the south zone (208.3 ± 134.5 colonies) followed by central and north zones (108.6 ± 55.3 and 107.5 ± 92.8 colonies,

Habitat factors and management status
A PERMANOVA test on the effect of habitat factors and management showed significant effect (p<0.05) of depth and management on benthic variables (Table 3). While exposure and reef type did not show significant effect, their interaction with depth and management was significant.
Embedded on the four quadrants of a PCoA for the key variables is an illustration of different combinations of habitat and management interactions at the three zones (Fig 4). The upper-left quadrant represents a habitat found only in the north zone with the benthic community highly dominated by crustose coralline algae. This habitat is characterized by Deep-Exposed Patch reefs that occur within Reserves (DEPR). The upper-right quadrant contains sites from central and south zones dominated by soft corals, Acropora and Montipora, and with a high density of juvenile corals (1-10cm) and small adult corals (11-20cm). This habitat is characterized by Deep-Exposed Fringing reefs within Unprotected areas (DEFU). The lower-left quadrant contains sites found mainly in the north and central zones. The benthic community is predominantly turf algae and macroalgae, and is characterized as Shallow Fringing and Lagoon reefs within Parks and Reserves (SFLPR). The quadrant to the lower-right contains sites from all geographic zones and has moderate levels of all variables. This represents a habitat that characterized as Shallow Patch and Channel reefs (SPC).
Of the four habitats described from this study, the north and central zones contained three of them and south contained two habitats (S1 Fig). DEPR was found uniquely in the north, SFLPR was found in all geographic zones, DEFU was found in central and south zones and SPC was found in the north and central zones.

Discussion
Two distinct patterns of coral reef benthic community structure were observed along the Kenyan coast. A latitudinal gradient was observed indicating a differentiation in benthic community composition between northern and southern sites. A second pattern was based on localised habitat and management factors prevailing within each geographic zone.

Latitudinal/geographic gradients in benthic community composition
The observed differences in benthic community composition among geographic zones found in this study is consistent with previous authors that have documented the biogeographic distribution of reef-building corals [37,16]. The south zone had a higher LCC with most genera occurring in abundance but mainly dominated by Porites and Acropora. This area is known to have a high diversity of corals mainly due to its proximity to the WIO center of biodiversity which has been documented to be around the northern Mozambique channel (NMC) [37,16]. The south zone is the receiving point of coral larvae into Kenyan waters as currents transporting larvae progress northwards [17] resulting in higher coral diversity compared to the north. In the north, a smaller number of genera was found with Porites being dominant, indicating a marginal coral community as it is farther from sources of larvae. The north zone is also an area of poor reef development due to the unfavorable conditions caused by the convergence of the colder nutrient rich SC and the EACC, which diminishes chances of high coral diversity in this zone. This convergence of currents makes the north an area of biogeographic transition from southern East African reef communities to Somali-Arabian Gulf communities to the north. This may mean that the north zone can act as an environmental corridor for species movements as the climate warms as has been suggested for other marginal communities [38]. Other studies have reported low coral diversity on marginal reefs e.g. upwelling reefs of the Arabian gulf [39], but there is also some isolated cases of high coral cover and diversity in marginal reefs e.g., at the world's southern-most coral reefs in subtropical Islands of the Pacific Ocean [40,41]. Some of these marginal reefs have the potential to act as habitable spaces for corals as climate warms, creating climate refugia or 'seed banks' that can re-populate coral communities once climate stabilizes to favorable conditions [42].
The benthic communities in the north and south represent two alternate aspects of coral reef resilience. The comparably lower abundance of juvenile coral colonies in the north and comparably higher abundance of large colonies suggests poor coral recruitment but high survival of adult coral colonies. Low recruitment could be attributed to isolation of the north zone from the southern sources of larvae, with only local sources of larvae for recruitment. Higher abundance of large colonies in the north may be due to the coral community in the north being protected from bleaching mortality by the cooling effect of upwellings in the SC. On the other hand, higher abundance of juvenile corals in the south indicates high recruitment likely due to its proximity to sources of larvae [17,16]. The low abundance of large colonies in the south implies low survival and only a few coral colonies surviving to large sizes. The south zone is known to face more stress from bleaching as it is closer to the warmest zones at 10-12 o S, and higher anthropogenic pressures due to higher human population densities as compared to the north. There have been several bleaching episodes since the 1998/99 mass coral bleaching which could have slowed down recovery of coral communities in the south.
Other studies in the WIO region have shown that the coral community in the north is very different from that in the south, with a higher diversity in the south comprising species more closely resembling the northern Tanzanian coral community [43,16]. Studies from other regions have also shown that the latitudinal gradient drives differences in benthic community composition on coral reefs. Coral reefs in the Hawaiian Archipelago in the north Pacific Ocean show a latitudinal gradient contributing to three distinct reef regimes each dominated by either hard coral, turf algae or macroalgae [44]. In a Caribbean-wide study, spatial and temporal differences in coral cover and macroalgae were found between subregions [45]. Eastern Australia coral reefs showed a trend of decreasing cover of Acroporidae as latitude increased.
The diversity of benthic communities differed along the latitudinal gradient and across geographic zones in Taiwan [46].

Influence of habitat factors and management on benthic community composition
Habitat factors and management influenced the pattern of benthic communities within each geographic zone. This offers the opportunity to propose ecological controls that structure benthic communities at this scale, and which may be used to guide resource management decisions. Similar findings have been reported in parts of the Pacific Ocean, where benthic communities were structured by depth, exposure to waves and reef types, in the Hawaiian Islands [47][48][49] and in French Polynesia [50][51][52][53][54].
The heterogeneity and uniqueness of habitat types in the north of Kenya can be attributed to a number of factors including the fact that it is composed of a set of small islands in an archipelago, presenting diverse reef types at different depths. Perpendicular to the shore, in a crossshelf line, there are numerous channel reefs between the islands and mainland, then lagoons, fringing reefs and at the outermost edge, a deep ridge with patch reefs on it. In addition, the northern part of Kenya has been described as an 'ecotone' zone due to the presence of marginal/transitional coral reef communities composed of rare or regional endemic species [16]. Other studies have also shown significant differences in benthic communities within archipelagos such as distinct coral species assemblages among different islands in Penghu in Taiwan [55], a cross-shelf difference in benthic community within the Spermonde Archipelago in the Coral Triangle [56] and differences in benthic communities across habitats in the Abrolhos archipelago of eastern Brazil [57].
The habitat type DEPR that was uniquely found in the north zone is dominated by coralline algae, which may be due to a high biomass of fish on the patch reefs [21,58]. Fish herbivores graze away competitors of crustose coralline algae, and fish predators feed on sea urchins which are the greatest eroders of coralline algae [59]. Also, the upwelling currents in this zone could be a source of nutrients stimulating growth of coralline algae. Previous studies have associated coralline algae as a suitable cue for coral settlement and recruitment [60,61] making these habitats potential spaces for corals to grow. Being deep habitats in an area of upwellingcooler waters could offer an ecological refugia to coral communities as the climate warms [62]. However, light limitation would reduce the opportunity of a diverse coral community growing and dominating these spaces [63][64][65] and only a few taxa that can tolerate these depths would colonise this habitat.
A majority of deep-fringing reefs that are exposed to oceanic waves in the central and north zones were dominated by Acropora, soft corals and small-sized coral colonies (DEFU habitat type). This habitat type is characterized by gradual slopes on the fringing reefs creating an environment where Acropora can dominate competitively. Other studies have reported Acropora occupying similar environments of forereef areas, where there is low stress from sedimentation and high rates of development processes such as reproduction [66]. Notably, the low abundance of large colonies within this habitat is worrying considering that the south zone is an area of high coral diversity. This could be indicative of failure to replenish standing stocks of large colonies in the central and south zone, likely due to bleaching mortality which has been recurrent over the past two decades [26,67]. All zones were badly affected by the 1997/98 mass coral bleaching episode with a 50% -80% average loss of coral cover with some individual sites loosing up to 100% cover [68,58]. Recovery has been slow since then perhaps due to the occurrence of smaller-scale bleaching episodes in 2005, 2007, 2010 and 2016 [26,23,67]. In the south zone, this bleaching impact was likely compounded by other anthropogenic stresses originating from its high human settlement. A similar explanation goes for the central zone reef where this study found few large colonies likely due to a slowed down recovery from the mass coral bleaching. Post-bleaching reports indicate that there was an increase in macroalgae and a decrease in fish biomass followed by an explosion in sea urchin density probably as a response to the decrease of fish biomass which constitute their predators. Consequently, sea urchin grazed the reef substate to bare [25]. The presence of poor water conditions due to proximity to river Sabaki effluents, high fishing pressure and recurrent bleaching episodes could have made perseverance of coral colonies to large sizes a great challenge [69,23].
Shallow fringing or lagoon reefs in the north and central zones (SFLPR habitat type) were characterized by a dominance of macroalgae and turf algae. Most of these habitats are proximal to land and are highly influenced by land-use activities such as effluents from rivers as well as fishing. When hard corals die macroalgae and turf algae take over the space due to increased nutrient availability or reduced herbivory [70][71][72]. Turf algae are the main focus of grazing but their form allows for rapid regrowth [73] making them thrive even in protected areas where grazing is high. Other studies have found that macroalgae are not influenced by latitudinal effects, with more local-scale influences such as nutrient levels strongly driving their distribution [74]. In this study, the abundance of macroalgae within these habitats indicates a worrying level of coral reef degradation since most coral cover and taxonomic diversity has always been thought to occur within shallow depths [75][76][77].
Shallow habitats on patch or channel reefs (SPC) were dominated by large colonies with an abundance of Porites and Echinopora. Similar observations were made in the Great Barrier Reef where large colonies of Porites were associated with shallow inshore reefs [78]. The presence of large coral colonies in shallow reefs indicates a habitat type that supports the survival of corals probably due to acclimatization to bleaching [79]. Also, the high fluctuation of water temperatures in shallow reefs creates a variable environment that may support different coral genotypes. A diversity of genotypes offers a chance for coral colonies to survive to larger sizes and persist over time as seen in an inshore reef of the Great Barrier Reef where Porites dominated through very persistent genotypes [80]. In addition, Porites and Echinopora have been considered as resistant and generalist coral genera within the functional groups of corals, based on their response to bleaching [81]. This makes them capable of inhabiting shallow areas where competitive but less tolerant corals cannot survive the fluctuating temperatures, bleaching events and proximity to human disturbance such as sedimentation.
Porites is a widespread genus that occurs at different habitats and thus reinforces its functional role as a resistant species [81][82][83]. The presence of very large colonies of Porites in the north compared to the south could imply that bleaching mortality of Porites is lower in the north due to the presence of habitat types that support resistance properties for that genus. In the south zone, the high percentage cover of Porites is mainly composed of small colonies. This signifies an area which could be facing high bleaching mortality, as evidenced by high percentage cover of dead standing corals, but experiences high recruitment that recovers the standing stock. This recovery feature may not be present in the north because Acropora communities failed to recover after the 1998 coral bleaching episode and recruitment is persistently low [21]).

Implications for conservation
In order to effectively manage and mitigate the degradation of coral reefs, information on community spatial patterns, habitat and management influences is necessary. Bleaching episodes have become more frequent and under the most recent IPCC report, 99% of coral reefs are expected to be lost this century if temperatures increase 2˚C above the pre-industrial baseline [6]. However, if temperatures can be stabilized at 1.5˚C, then 10~30% of the reefs could be saved [84]. These reefs could serve as refugia [62,85] or 'seed banks' to repopulate other reefs once mechanisms of reducing climate change to desirable levels have been achieved [86].
Management strategies should ensure the protection of habitats that exhibit bleaching protection or recovery characteristics in order to enhance the resilience of coral communities. In addition, there should be mechanisms to reduce local impacts such as sedimentation and pollution. High sedimentation levels are experienced on near-shore coral reefs near rivers with high sediment load, particularly in the central zone due to the river Sabaki and river Tana [87]. Most of the reefs in this area are within parks and reserves but the effect of sediment load is likely to still degrade these reefs. In the south zone, sediment loading is mainly experienced from smaller rivers and creeks, with increasing threat from logging of mangroves, thus exposing sediments to erosion. In all these areas, increasing urbanization elevates discharge of sewage and other pollutants, further exacerbating stress to corals.
In conclusion, the occurrence of a geographical and a local habitat pattern in coral benthic communities that is further mediated by the level of management offers a range of environments that support coral reef communities. The north zone presents a unique environment of a marginal reef that provides protection to coral communities from bleaching [88] even as the climate warms, expressing resistance properties of coral reef resilience [89]. The south zone's high diversity of corals [16] and high recruitment potentially offer a habitat that expresses recovery properties after a bleaching disturbance [88]. The central zone is intermediate, providing a corridor for transfer of coral propagules between the north and south. Similar propositions on the importance of ecological corridors in conservation have been mentioned for large scale studies in the Eastern Tropical Pacific [90] and within the Gulf of Mexico [91].
The heterogeneity of coral reef habitat types in Kenya offers an opportunity for designing MPAs that protect a diverse range of functional traits of coral communities [92]. As such, spatial heterogeneity in MPAs may reduce the risk of catastrophic regime shifts [93].
Supporting information S1