Introduction

Coral reefs are currently facing multiple stressors arising from a combination of local and global threats that have contributed to coral cover loss1,2,3. Increasing sea surface temperature (SST) has resulted in more frequent and intense climate-induced bleaching events, impacting coral reefs globally4,5,6,7. Thermal anomalies and coral bleaching have increased in frequency, from once every 25 to 30 years at the beginning of the 80’s, to once every 5.9 years by 20166, affecting even regions considered coral refuges8,9,10,11. Regardless of the projected future climate scenarios, coral reefs are expected to be degraded, resulting in changes in species composition12, including a general loss of biodiversity and ecosystem services13,14. Indeed, by 2050, it is predicted that coral cover will decrease globally by 5–14%, while by the end of the century more than 40% of coral may be lost from the world’s oceans15.

Coral bleaching occurs when, under sustained stress, corals expel their algal symbionts from their transparent polyp tissue, thereby exposing the white, calcareous coral skeleton16. In a zooxanthellate coral, these symbionts may provide over 90% of the energetic requirements for the host17. If the duration is long enough, bleaching can cause coral death18 and consequently reef degradation16,19, thus directly affecting reef functioning20,21. Reef degradation often manifests as a shift towards dominance by macroalgae or other non-calcifying organisms (e.g. soft corals, sponges) which are able to outcompete calcifiers (e.g. hard corals) after a disturbance3,22,23. This shift in community structure and the potential loss in coral reef-associated biodiversity has a direct effect on reef functioning, affecting diverse ecological processes such as biogeochemical fluxes, primary production, herbivory, and predator-prey interactions24,25,26. Ultimately, these ecological changes will be translated into negative economic effects27.

The assessment of long-term changes to coral reef benthic community composition as a response to stressors associated with climate change requires the investigation across large temporal scales, supported wherever possible by well-established and standardized monitoring programs28. Nevertheless, in many tropical countries long-term monitoring programs are scanty and rarely standardized; not only are temporal data scarce or missing, but often baseline data themselves is absent, making evaluation of the ecological effects of coral bleaching challenging29. Such lack of data hampers the ability to make projections based on future scenarios. This is particularly true coral reef communities in the Red Sea, which remain largely understudied compared to other major reef systems (e.g. the Great Barrier Reef)30. Red Sea coral reef community responses (and recovery trajectories) to bleaching are not always documented as national monitoring programs have been limited in scope and duration. Monitoring by local stakeholders often faces the same obstacles (but see PERSGA31).

Over the last decade, the Red Sea’s coral reefs have seen several heat-related coral bleaching events32,33,34, and more recently, seasonal bleaching due to cold stress35. Though heat-related bleaching events in the Red Sea are beginning to garner attention, Genevier et al.36 reported that thermal anomalies in the Red Sea have increased in frequency and intensity since 1998, mirroring the trend seen worldwide37,38. In the aforementioned study, the authors determined historical marine heatwave events in the region by analyzing 35 years of SST data, identifying areas at risk, and highlighting potentially undetected bleaching events. According to the study, marine heatwaves encompassed more extensive areas and occurred more frequently than previously reported, even in the non-El Niño years, suggesting that coral bleaching may have been underestimated in the past.

The Red Sea is characterized by well-known north-south gradients of salinity, SST, and primary productivity39, making the basin an ideal natural laboratory to study coral reef dynamics. Communities in the three regions of the Red Sea (northern, central, and southern) are exposed to different environmental conditions and, therefore, may respond to SST anomalies in contrasting ways due to their different environmental history.

Based on a network of reefs that have been monitored since 2014 along the Saudi Arabian Red Sea coast (revisited nearly every two years) using standardized photo-transects, we investigated spatial and temporal changes in the community composition of coral reef benthic organisms over a six-year period. This period (2014–2019) was comprised of surveys conducted before and after the third and longest global-scale coral bleaching event caused by the El Niño-Southern Oscillation (ENSO) in 2015/201619,40. This study provides insights into the dynamics of benthic community structure across regional and temporal scales in the hypersaline and naturally warm Red Sea, and lays the foundation to track future trajectories of recovery or changes in these coral reefs. The relationship between patterns in Degree Heating Weeks (DHW) and SST data was also investigated (following findings by DeCarlo et al.41), as were the reported regional differences in bleaching prevalence. We hypothesize that coral communities of the northern Red Sea will be the most resistant to marine heat waves due to their higher thermal tolerance; indeed, this region is often considered a coral refuge that could potentially help to repopulate other regions42,43.

Materials and methods

A total of 59 reefs along the Saudi Arabian Red Sea coast (28° N to 16.68° N), were sampled every two years (2014–2019) for three time points (Fig. 1, Supplementary Table S1). Sampled reefs were grouped into three regions (North, Central, South) based on observations in previous studies on spatial patterns of the coral reef cryptobiome44, spatiotemporal patterns of salinity and SST, and the distribution and abundance of surface phytoplankton biomass45. The regions and sites are as follows: Northern Red Sea (NRS) region: 20 reefs among two locations, Central Red Sea (CRS) region: 17 reefs among two locations, and the Southern Red Sea (SRS) region: 22 reefs among three locations (Fig. 1). Not all the reefs were sampled at every period due to limited human resources and complications in obtaining permits for some localities.

Fig. 1
figure 1

Location of sampling sites along the Saudi Arabian coast of the Red Sea. The 59 reefs are distributed across seven locations within three regions: North Red Sea (NRS) represented in shades of blue (Duba, Al Wajh), Central Red Sea (CRS) in shades of pink (Thuwal, Jeddah), and South Red Sea (SRS) in shades of green (Al Lith, Farasan Banks, Farasan Islands). Horizontal lines delimit the different regions. Maps throughout this paper were created using ArcGIS Pro software by Esri (software version 2.5.22; www.esri.com). ArcGIS and ArcMap are the intellectual property of Esri and are used herein under license. Copyright Esri. All rights reserved.

Sea surface temperature

To characterize thermal conditions at the study sites, we retrieved SST remote sensing data at a resolution of 4 km from the MODIS system. For this, a running 8-day mean of longwave (11 μm) night-time Level-3 data spanning seven years (June 2013 to May 2020) was obtained for each site from the Aqua satellite46 via NASA’s Ocean color group website (http://oceancolor.gsfc.nasa.gov/). Occurrence of heat stress during the same period was assessed using NOAA’s Coral Reef Watch Degree Heating Week (DHW) product47. Daily DHW data at a resolution of 5 km was downloaded from the Pacific Islands Ocean Observing System (PacIOOS) via the Environmental Research Division’s Data Access Program (ERDDAP) data server (https://pae-aha.pacioos.hawaii.edu/erddap/griddap/dhw_5km.html). According to Coral Reef Watch48, coral bleaching is likely to occur in an area when DHW reaches values of 4 and above, and severe bleaching (with subsequent mortality) is likely for DHW values of 8 and above. Given the resolution of the remote sensing products used here (i.e., 4 km for SST and 5 km for DHW data, respectively), some sites located near each other were associated with the same remote sensing data tile.

Benthic assessments

To record benthic cover at the selected reefs, semipermanent transects were surveyed at a depth of around 10 m. At each site, GPS coordinates were taken and used to deploy a transect tape along which surveys were conducted. Local observations of the reef structure were also taken to guide future assessments. In some cases, these were related to the presence of Autonomous Reef Monitoring Structures (ARMS) deployed at many of the reefs surveyed that marked the starting point of the transects. From this starting point, a series of three replicate 20 m long and 1 m wide benthic transects were conducted taking sequential photos (1 m x 1 m) for future analysis. A consistent distance to the bottom was the benthic composition (i.e. percent cover of benthic groups) for each photo was analyzed using Coral Point Count software with Microsoft Excel extensions (CPCe)49. In each 1 m2 frame, 48 points were randomly generated. This point count method has previously been applied to infer percent cover and has been shown to provide comparable results to line intercept methods50, with the added benefit of allowing for an archive of photos that can be used for future assessments and taxonomic validation. Here, we averaged the calculated percent cover of benthic categories across all photos of each transect. Subsequently, the percent cover at the reef level was calculated as the average of the three replicate transects for further analysis.

The benthic composition of each randomly generated point was assigned one of the following categories: hard corals (including calcifying hydrozoans of the genera Distichopora and Millepora), soft corals (including gorgonians), other anthozoans (zoanthids, anemones, and corallimorpharians), fleshy macroalgae, turf algae, Crustose Coralline Algae (CCA), other invertebrates (e.g., sponges, bivalves), pavement/rock, debris (plastic debris, litter, nylon nets), rubble (small hard fragments, rocks or lose coral fragments), and sand. Hard corals, soft corals, and macroalgae were recorded to the genus level whenever possible. Acropora and Porites specimens were further separated based on their growth forms into the following categories: table (Acropora), branching (Acropora, Porites), columnar (Porites), massive (Porites), and encrusting (Porites). The full list of benthic categories is shown in Supplementary Table S2.

Data analysis

The study period was divided into three different periods: period 1 (2014/2015), period 2 (2016/2017), and period 3 (2018/2019). This division was based on the occurrence of a major bleaching event in the summer of 2015 that differentially affected coral reefs along the Red Sea34,41. All surveys conducted in 2015 were carried out before the summer, therefore they were combined with the 2014 surveys and considered the reference or pre-bleaching period. For the assessment of changes in benthic community composition over time and across locations as well as regions, we performed a Non-Metric-Multidimensional Scaling (NMDS) analysis of all benthic categories. Percent cover data was square root transformed for normality, and a resemblance matrix was created using Bray-Curtis distance. A cluster analysis was conducted using the same resemblance matrix to confirm the associations of samples. We generated a heat map with the relative percent cover of the main benthic categories per site and period (i.e. any category with percent cover greater than 2%): hard corals, soft corals, algae, turf algae, CCA, rubble, pavement/rock, and sand. To visualize differences in the abundance of hard and soft corals pre- and post-bleaching, a bar graph was created using the ggplot251 package in R v4.3.052 via RStudio v2023.09.153. This graph showcases the relative cover of both hard and soft corals on the reef that were surveyed before and after the anomalous event. To investigate the coral genera (including in some cases morphotypes) that contributed to the differences among the regions within each period, a SIMPER analysis was conducted using a Bray-Curtis dissimilarity index based on square root transformed data. All analyses were conducted in PRIMER-v754 with the PERMANOVA + extension55 unless otherwise noted.

Differences in the benthic community composition were further analyzed using a four-way Permutational Multivariate Analysis of Variance (PERMANOVA)56 on the following factors: region (three levels: NRS, CRS, SRS, ), period (three levels: 2014/2015, 2016/2017, 2018/2019), location (seven levels: Duba, Al Wajh, Thuwal, Jeddah, Al Lith, Farasan Banks, Farasan Islands) and reef (59 levels), based on the same dissimilarity matrix described above. The factor ‘region’ was fixed, and the rest of the factors were considered as random; ‘location’ was nested within ‘region’, and ‘reef’ was nested within ‘location’. To investigate comparisons within the different factors, pair-wise PERMANOVA comparisons were conducted in PRIMER-v7. To assess changes in the percentage of coral cover over time, the difference in coral cover was calculated for reefs sampled in more than one period and subsequently binned in the following categories indicating percent gain or loss in coral cover: +10 to 15%, + 6 to 10%, 0 to 5%, -1 to -5%, -6 to 10%, -11 to 15%, and ≤-16%. These reefs were plotted on a map according to the geographical location of the reef using the ArcGIS Pro software (Version 2.5.22) by Esri. To analyze the changes in coral cover per genus in the three regions, the species that contributed to the dissimilarity in the regional coral assemblages were selected according to the SIMPER analysis mentioned above.

Differences in thermal regimes (i.e. SST and DHW) were determined using a non-parametric Kruskal-Wallis rank sum tests with R’s base function kruskal.test and post-hoc Dunn’s tests using rank sums with the function dunn.test from the package of the same name (v1.3.557). Tests were carried out for the following factors: ‘year’ (seven levels: 2013/14, 2014/15, 2015/16, 2016/17, 2017/18, 2018/19, 2019/2020), ‘location’ (seven levels, see above) and year-location (49 levels; no post-hoc test was carried out for this interaction factor). To reduce zero-inflation, DHW data were subset to the months June through December for statistical tests (no heat stress occurred in any year at any of the reefs during January-May). Mean and standard deviation of SST and DHW for each location and year were calculated using R base functions. Data were formatted and visualized using the R packages reshape2 v1.4.458, ggpubr v0.4.059, gridExtra v2.360 and packages from the tidyverse collection v1.3.161. All analyses and visualizations of SST and DHW data were performed in R via RStudio and the code is found in Supplementary File 3.

Results

Thermal regimes

SST exhibited seasonal variability during the study period, with minima of ~ 21.5 °C in higher latitudes (i.e. in the NRS region in winter) and maxima of ~ 33.0 °C in lower latitudes (i.e. the SRS region in fall) (Fig. 2A). Annual SST variability within locations decreased from north to south, which was expected given generally more pronounced seasonality in higher latitudes62. In most years, locations towards the south showed a second minor SST peak in early summer before temperatures dropped and then increased again to the annual maximum (Fig. 2A).

Timing of northern-hemisphere summer fall temperature peaks showed a latitudinal trend, with northern locations exhibiting maxima in August/September and central and southern locations in October / November. Highest annual SST maxima were observed in the south (i.e., Farasan Islands and Farasan Banks). See Supplementary File 1 for the full SST data set. The year 2015/16 displayed the highest annual SST during the study period, with an average of 28.46 ± 2.72 °C. This year also represented the warmest period on average within each location, except for Thuwal and Jeddah (for which 2018/19 and 2016/17 were the warmest years, respectively). The mean SST during 2015/16 at Farasan Islands represents the sole incidence of our study locations exceeding an annual mean SST of 30 °C (30.23 ± 1.93 °C). The three SRS locations (i.e., Al Lith, Farasan Banks, Farasan Islands) exhibit a clear plateau of elevated SST in late summer/fall of 2015 (Fig. 2A). Though peaks of similar magnitude occurred in other years, they were ephemeral. The distinguishing factor in 2015 is the clear plateau of elevated temperatures lasting between one to two months occurring at all locations in the SRS (and to some extent in the CRS) (Fig. 2A; Supplementary File 1).

The SST profiles differed significantly among years (Kruskal-Wallis χ­2 = 103.29, p < 0.001), and SST was significantly higher in 2015/16 than during the preceding two years and the following year (pairwise adjusted p < 0.001). Thermal conditions at the seven locations mirrored the separation of study sites into three regions. Sites at Duba and Al Wajh (i.e., NRS), Thuwal and Jeddah (i.e., CRS) and Al Lith, Farasan Banks and Farasan Islands (i.e., SRS), respectively, exhibited similar temporal trends in SST among them (Fig. 2A). This was supported by results of the statistical tests, where SST regimes displayed significant differences among locations during the entire study period (Kruskal-Wallis χ­2 = 4371.60, p < 0.001). Pairwise comparisons supported the regional grouping of the study locations. Intra-regional differences in location SST were all non-significant, while inter-regional differences were significant in every case (all pairwise adjusted p < 0.001). SST profiles were significantly shaped by the interaction of the factors ‘location’ and ’year’ (p-value < 0.001) which was expected given the large sample size. See Supplementary Table S4 for detailed information on mean SST per location and year and Supplementary Material S2 for results of statistical tests.

The occurrence of heat stress measured as DHW showed marked differences among locations and years, although the general trend of DHW peaking in summer/fall at all locations was consistent (Fig. 2B). As with SST, locations in the north showed a DHW peak in August/September, while this peak shifted towards October/November with decreasing latitude (Fig. 2B). Duba exhibited the highest DHW across all locations in each summer/fall season during the study period (maxima ranged between 2.57 °C-weeks between August 31 and September 26, 2017 and 12.58 °C-weeks between September 15 and 18, 2019), with the exception of 2015/16 when Al Lith experienced the highest DHW (maximum of 8.82 °C-weeks between October 23 and December 1, 2015). Duba also showed highest overall mean DHW throughout the study period (3.07 ± 3.79 °C-weeks). Overall mean DHW consistently decreased towards the south, with the exception of the Farasan Banks which showed the lowest mean DHW (0.32 ± 1.05 °C-weeks). See Supplementary Table S5 for data on means per year and location. Total mean DHW per year across the entire Red Sea was highest in 2015/16 (2.33 ± 2.60 °C-weeks) and lowest in 2017/18 (0.22 ± 0.53 °C-weeks). All locations southwards from (and including) Jeddah experienced their highest heat stress during the study period in 2015/16. Remarkably, 2015/16 represents the only period of significant heat stress in the SRS (i.e. Al Lith, Farasan Banks and Farasan Islands) (Fig. 2B and Supplementary Table S5). See supplementary File 2 for the full DHW data set.

Fig. 2
figure 2

Time series of Sea Surface Temperature (SST, A) and Degree Heating Weeks (DHW, B) between June 2013 and May 2020 for the seven locations examined in this study. Profiles represent the averages at each point in time for all surveyed reefs within each location.

Benthic cover

A total of 106 living benthic categories were identified on the reefs surveyed, including 57 hard corals, 18 macroalgae, 13 soft corals, 11 other invertebrates, 5 taxa of other cnidarians, CCA and turf algae (Supplementary material Table S2). The main benthic categories with a relevant contribution in benthic percentage cover (i.e., ≥ 2%) in the reefs along the Saudi coast of the Red Sea were hard corals (17.5%±12.5), soft corals (7.58%±8.6), CCA (14.43%±13.1), turf algae (21.94%±12), macroalgae (6.61%±12.7), sand (21.01%±19.4), rubble (2.79%±5.9) and, pavement/rock (6.32%±9.9).

NMDS analysis showed that the benthic communities relatively homogeneous in composition and did not differ among locations before the bleaching event (Fig. 3). NRS and CRS locations generally clustered together, whereas those from the SRS showed the largest variability, with some reefs presenting high affinity to both central and northern regions. After the 2015 mass bleaching event (period 2), communities from the south had a different composition from those in the north and central regions (Table 1) and were clearly segregated as a unique region in the multivariate space based on our data (Fig. 3). However, three-four years post-bleaching (period 3), most of the reefs in the south were still generally distinct from reefs in the CRS and the NRS regions, although some southern reefs again showed affinity in composition to the northern and central locations (Fig. 3).

Fig. 3
figure 3

Non-metric multidimensional Scaling Plot (NMDS) of the coral reef benthic communities, for the three geographical regions: NRS (blue), CRS (pink) and SRS (green); (A) Period 1: 2014/2015, (B) Period 2: 2016/2017, (C) Period 3: 2018/2019. Hierarchical cluster analysis using the Bray-Curtis dissimilarity index for seven different locations (DU: Duba, AW: Al Wajh, TH: Thuwal, JD: Jeddah, AL: Al Lith, FI: Farasan Islands, FB: Farasan Banks) in three periods (Period 1: 2014/2015, Period 2: 2016/2017, Period 3: 2018/2019.). Heatmaps show the percentage cover of the benthic categories with a cover greater than 2% (hard corals, soft corals, CCA turf algae, algae sand, rubble, pavement/rock); the darker the shade, the higher the percentage cover.

In period 1, the benthic communities in the south were dominated by hard and soft corals (Fig. 4). In period 2 (post-bleaching) the abundance and richness in corals decreased, and 15 hard coral categories disappeared completely, including Porites branching, Porites encrusting, Psammocora, Euphyllia, and Plerogyra (Fig. 3). In contrast, macroalgal percent cover increased, particularly the genera Dyctiota, Lobophora, and Peyssonnelia; turf algae also increased in percent cover, though to a lesser degree. When analyzing the magnitude of change among periods, the SRS reefs were the most affected, showing limited recovery four years after the bleaching (Fig. 5).

Fig. 4
figure 4

Soft and hard coral cover before and after the 2015–2016 bleaching event among the three regions. The bars show relative cover, and codes to the right represent the reef surveyed before and after the bleaching order latitudinally.

Significant statistical differences in benthic composition were found among regions (Pseudo-F = 2.1216, P = 0.007), periods (Pseudo-F = 13.224, P = 0.001), and zones (Pseudo-F = 1.5249, P = 0.032) (Table 1). Interaction effects between region and period (Pseudo-F = 2.8754, P = 0.009) and between reef and period (Pseudo-F = 4.6073, P = 0.001) became apparent with the PERMANOVA analysis.

Table 1 Results of four-way permutational multivariate analysis of variance (PERMANOVA) on factors region (R), period (P), locations (Lo) and reef (re), based on square-root transformed percentage benthic cover data and Bray Curtis similarity index. Asterisks indicate statistically significant values. Df = degrees of freedom, MF = mean of squares, Pseudo-F = variance ratio, P = significance level.
Fig. 5
figure 5

Changes in hard coral cover at the study locations along the Saudi Arabian Red Sea coast. Colored points indicate the locations where an increase (shades of green) or decrease (from yellow to dark red) in coral cover was recorded across the different periods. (a) From period 1 (2014/2015) to period 2 (2016/2017), (b) Form period 2 (2016/2017) to period 3 (2018/2019), and (c) from period 1(2014/2015) to period 3 (2018/2019). Created using the ArcGIS Pro software (Version 2.5.22).

In the NRS, the most affected (i.e. biggest decrease in percent cover) genera were Pavona, Acropora branching, and Montipora (Fig. 6). In the CRS region, the genera with the highest decrease in coral cover recorded were Pavona, and Acropora table, although there was an increase in the relative proportion of Porites massive. In the SRS, Pavona, Acropora branching, and Montipora showed the largest reductions, but here we also assisted to a decrease in cover of Porites massive; they persisted in period 3, while the rest of the genera on the list disappeared in this period (Fig. 6).

Fig. 6
figure 6

Differences in percentage of cover of the coral genera per region in the three study periods represented by different colors (period 1: 2014/2015; period 2: 2016/2017; period 3: 2018/2019). Coral genera were selected based on SIMPER analysis (Supplementary Table S3). Colored dots represent the average and the lines the standard deviation. Black dots on the right represent the overall changes in percentage cover between periods 3 and 1.

Discussion

This study documented latitudinal and temporal changes in the composition and structure of coral reef benthic communities over a 5-year period encompassing the 2015–2016 mass bleaching event along the Saudi Arabian coast of the Red Sea. The results show that Red Sea coral reef benthic communities responded differently to a mass bleaching event dependent on the environmental context they are set in. This study also highlights the value of in situ biological data regards to assessing the effects and recovery trajectory of bleaching events. Indeed, we show here a strong mismatch between the predictions of bleaching using DHW metrics and observations, reinforcing the synergistic effects between SST and other variables, like nutrients, triggering bleaching as previously reported by DeCarlo et al.41 in the Red Sea.

Interestingly, responses of coral genera appear to be region-specific; although some genera seem to generally be more sensitive (e.g., Pavona, Montipora, and Acropora branching), or more resistant (e.g., Porites massive), this does not hold true across all sites or regions. The nearshore reefs of the most affected region (SRS) did not show clear signs of recovery four years after the mass bleaching event. The variability in coral bleaching across different sites, where some corals are affected while others remain unscathed, can be attributed to multiple factors. This diversity in bleaching patterns suggests that the phenomenon is influenced by a complex interplay of elements, probably including water temperature, water quality, light intensity, local adaptation of coral organisms and acclimatization, individual thermal history, and the type of symbiotic zooxanthellae present63,64.

The increased CCA cover reported in this study may signify an encouraging possibility for coral reef regeneration, as CCA plays a pivotal role in fortifying the reef substrate through the secretion of calcium carbonate, and is linked to the settlement and attachment of coral larvae65,66. This calcified deposit provides an anchor for coral development and acts as a protective shield, mitigating the effects of physical erosion67. Even so, any potential for positive outcomes are counterbalanced by the influence of human-induced pressures. Factors such as the overfishing of reef herbivores and nutrient enhancement can impede these beneficial effects, thus hindering coral reef recovery68.

Region-specific responses of Red Sea benthic communities to the 2015/2016 mass bleaching event

Investigating the patterns of variability in coral reef benthic communities through space and time, we found a region-specific response that aligned with the prevailing SST trends and the known productivity patterns in the region45. Specifically, benthic communities from the southern Red Sea were characterized by high coral cover (25.9%±19.4 in period 1), a low cover of fleshy macroalgae (3.1%± 6.0), and a substantial cover of soft corals (8.0%±8.4) before the bleaching event, but were severely affected by the 2015/2016 mass bleaching event. This result based on in situ photo transects supports the findings by DeCarlo et al.41 The authors detected a higher intensity of the 2015/2016 bleaching event in the southern Red Sea compared to previous bleaching events affecting the Red Sea region (i.e., 1982, 1995, and 2002), hence a region-specific signature of these events. Although the SRS region is known for its comparatively high primary productivity due to Gulf of Aden Intermediate Water (GAIW) intrusion69,70 and monsoon-driven upwelling71,72, corals are still very sensitive to the synergistic effects of both temperature and nutrient stressors when above certain thresholds. Indeed, it seems that the combination of thermal anomalies and an increase in nutrients from an early upwelling caused by the effects of El Niño may determine the intensity of bleaching events in the region41,72. In other words, corals in the Red Sea can resist an increase in SST if not accompanied by an increase in nutrients, as observed in 1998 or 2010 in the southern Red Sea41. Another study reported that environmental pressures like heat stress and higher levels of dissolved organic carbon (DOC) could potentially accelerate the breakdown of the coral-algae symbiotic by boosting the rates of nitrogen fixation linked to corals73.

In this study, a relatively high turf cover was also recorded in period 1 (16.6%), and it increased after the bleaching event (24.1%); this trend was strongest in the SRS which increased from 16.5% in period 1 to 28.4% in period 2). Although there are several factors related to the decline and loss of coral species, including coral disease, invasive species, and climate change, PERSGA31 reports that major threats in this region are primarily linked to human pressures such as overfishing, coastal development, population expansion, and industrial pollution. Thus, the observed increase in turf algae cover may be a part of the broader pattern of ecological shifts driven by these anthropogenic factors and not an outcome of the bleaching event. It is important to note that these factors can contribute to decrease the tolerance of coral species to additional stressors like those related to increased SST usually associated with bleaching events.

The most remarkable, yet expected, response to the 2015/2016 bleaching event was a significant decrease in coral cover. Particularly for hard corals (both branching and massive forms demonstrated clear reductions, showing a 20-fold drop in cover (from 25.8 to 1.2%). A similar, though less drastic, decrease in soft coral cover was also detected (from 8 to 0.3%). Such sharp decreases in coral cover have been reported in response to bleaching events worldwide6,74,75, followed by either an increase in turf or macroalgae (an increase from 15.3 to 25.9%, and 3.1–34.8%, respectively for the present study)76,77.

Even though Red Sea coral reef communities have widely been considered resilient, namely due to their high recruitment rates78,79, more frequent bleaching events in the region may hinder their potential for recovery, as observed elsewhere16,80. This can be particularly true for nearshore reefs that are more affected by the bleaching events due to, among other stressors, more drastic variation in SST than offshore reefs81. Indeed, considering that four years after bleaching the percent cover of hard corals was about 12% of the pre-bleaching levels ( pre-bleaching: 29.8%±18.8%; four years post-bleaching: 3.6%±3.2%), signs of solid recovery were not detected. Corals in the southern Red Sea may face a serious threat as mass bleaching events are expected to occur on average every 5.9 years globally6. Nevertheless, though slow, signs of recovery were noticed. Average hard coral cover almost doubled in two years (period 2: 1.2 ± 0.8%; period 3: 3.6%±3.2%). It is important to note that such signs of recovery ere not homogenous, as clarified by the high standard deviation.

Taxa susceptibility to bleaching: toward a prevalence of massive forms with a loss of rapid growing branching corals

During a heat stress event, coral species will respond differently to changes in environmental conditions. In some cases, high plasticity may be observed, and response patterns may not be consistent across reef habitats or across a shelf gradient82. In general, branching corals are considered more prone to bleaching and subsequent mortality83,84,85,86. A similar pattern has been detected in previous studies in the Red Sea32,34 and is supported by our findings. A previous study on the Saudi Arabian Red Sea coast showed that the coral genera more sensitive to changes in environmental conditions were indeed mainly genera with branching growth forms, such as Acropora, Stylophora, Seriatopora, and some species of Porites87. In our study, Porites, Pocillopora, Pavona, and tabular and branching Acropora were the most abundant corals in the south in the first sampling period contributing to ~ 12% total cover. After the bleaching event, these genera specifically (and all branching corals in general) were severely affected and almost disappeared, confirming the growth form’s sensitivity to bleaching. In our study, the massive Porites category was the main driver of the community composition. Consistent levels (and a general dominance) of coverage highlights its relative tolerance to global warming effects. This is supported by a previous study in the Great Barrier Reef, which found that acroporids bleach faster and present a higher whole-colony mortality rate, while species such as Platygyra daedalea and Porites lobata, both massive species, were characterized by a later onset of bleaching and partial mortality88.

Changes in benthic community composition does not necessarily correspond with heat stress metrics

It has previously been shown that the occurrence of widespread coral bleaching and heat stress in the northern Red Sea have so far been disconnected from each other43. The present study also showed a discrepancy between heat stress and observed changes in benthic reef communities. Although showing constantly high heat stress each summer during the study period, reefs in the northern Red Sea experienced changes in benthic cover on a much lower scale compared to reefs in the central and southern regions (with the exceptions of some reefs at Al Wajh). These observations are in line with the proposal of the northern Red Sea as a thermal refuge for scleractinian corals42,43. However, this region is not without its own unique set of environmental stressors. Intense coastal development in the region associated with two major development projects (NEOM, and the Red Sea Global89,90) which can affect turbidity, nutrient availability, and local hydrodynamics, all of which have been shown to influence the response of corals to bleaching91,92,93. Temperature variability patterns experienced by corals can, in some cases, be more relevant to the susceptibility of corals than the maximum temperature populations experience over a certain period of time92,93,94,95. Significant coral cover decline in some of the CRS and SRS locations observed here may indicate that corals in these areas are already close to their thermal tolerance thresholds, while this may not be the case yet for corals in the NRS43. However, it remains unclear for how long reefs in the NRS will be able to withstand widespread bleaching and coral cover loss in the future given the already high rates of ocean warming in the north96.

The DHW peaks particularly in the southern locations combined with the early summer upwelling / GAIW-intrusion most likely played a major part in the bleaching pattern observed in the southern Red Sea, as reported by DeCarlo et al.41,72. Upwelling in the southern Red Sea in June to August 2015 was the strongest recorded since the 1980s but ceased earlier than usual, leading to relatively low SST during the upwelling but a longer than usual subsequent warming period41,72. This strong upwelling of GAIW also led to comparatively high nutrient levels in otherwise oligotrophic waters where (sub-) tropical coral reefs thrive, another major stressor for corals41. The combination of a long and intense warming period and a strong upwelling event leading to high ambient nutrient levels are most probably what caused bleaching in southern-central and southern Saudi Arabian reefs in 2015 and potentially the subsequent coral cover decline observed at these sites.

Conclusions

This study investigated the response of benthic communities along a latitudinal gradient on the Saudi coast of the Red Sea. We observed a significant regional variation in the effects of the mass bleaching event in 2015–2016. The southern region was the most affected, showing a major decline in hard coral cover, a major loss of branching corals, and a shift in the benthic composition towards an algal (macro- and turf) community, accompanied by a loss in coral diversity. In contrast, the northern and central reefs registered a small decline in coral cover (1–5%), and in some cases coral cover remained stable or even slightly increased.