Abstract
Demographic analyses offer insight into the state of a population. Here, we surveyed different reef flat zones (exposed, midreef and sheltered) of six reefs over a cross-shelf gradient to characterize the population structure of Stylophora pistillata, a coral species which dominates reef flats in the central Red Sea. Phototransects were conducted at each reef flat zone, and the density of S. pistillata, the planar area of each colony, and the occurrence of partial mortality were calculated using the program ImageJ. Each colony was also assigned a color morph (yellow, purple or mixed colors). Density and mean size were extremely variable, both among reef flat zones and reefs, but overall, both metrics were lower on the midshelf reefs. The yellow color morph accounted for nearly 90% of colonies surveyed and dominated most reef flats assessed, with the exception of one site where 81% of colonies were purple morphs. There were no spatial trends in the percentage of colonies suffering partial mortality, but overall there is a positive correlation with size class and proportion of colonies with partial mortality. Despite few trends emerging from assessing individual parameters, a PERMANOVA analysis revealed differences among reef flat zones in most of the reefs, highlighting the importance of multivariate analysis. The data presented here serve as a baseline for monitoring and may identify possible future demographic changes to this important coral species in a region increasingly affected by bleaching events.
Similar content being viewed by others
Introduction
Coral reef ecosystems are underpinned by hermatypic coral species, which construct the reef structure providing habitat for approximately one quarter of marine life1. At a fundamental level the growth and maintenance of the reef environment are dependent on the demographic rates of individual coral species (e.g., reproductive output, recruitment and settlement rates, colony growth rates, population density, and mortality)2. In recent decades, coral species have been severely affected by anthropogenic disturbances, compromising the stability of their populations and the resilience of coral reefs worldwide.
Long-term monitoring is ideal for gathering information on trends of population growth and decline, which can shed light on how species respond to changes in environmental conditions or disturbances. However, long-term data often are not available for many regions, especially remote coral reefs or historically understudied systems such as the Red Sea3. Historically, coral percent cover has been used as the standard monitoring variable in benthic reef communities, and is a useful metric that encapsulates very general information about the state of a coral reef. However, this traditional approach obscures details that can be useful in understanding population dynamics. This method ignores the actual abundance and sizes of individual corals, since “percent coral cover” only measures the total space occupied by a coral category. Assessing size structure (e.g., the abundance of individuals within different size classes in a population) is a better tool for understanding the population trajectories of a species because it captures demographic information not contained in other methods like percent cover4,5,6. Therefore, simple metrics such as size structure can offer a more in-depth “snapshot” of the state of a population. Often, size structure is used to infer the state of coral populations under anthropogenic stressors, such as pollution7, climate change8,9,10,11, and fishing pressure12,13. Therefore, in an era of intense degradation as a result of local and global disturbances, a more comprehensive evaluation of the status of coral reefs, potential causes of decline, and projected trajectories is essential to guide mitigation and conservation measures.
Size structure comprises two key metrics: abundance and the size of all individuals7,14. In corals, abundance is the product of several demographic processes, including asexual reproduction, settlement and mortality rates, while the size of individual coral colonies is dependent on growth rates and partial mortality15,16,17. The resulting size frequency distribution provides a way to compare demographic information among different species, populations or sites14. For example, a size-frequency distribution that is negatively skewed indicates an overabundance of large individuals in a population, while positive skew denotes a population dominated by small individuals. Previous studies have suggested that negatively skewed size-frequency distributions are potentially indicative of a degraded reef14,18. The implication is that there are few young individuals to replenish older individuals in a population, and may be due to a lack of (or unsuccessful) settlement/recruitment or high mortality rates for smaller size classes. While skewness alone cannot indicate the nature of degradation, it can be relevant for future investigations to understand the cause of demographic shifts.
Coral size structure varies naturally in time and space. Abiotic factors such as depth19,20, temperature21 and water quality22 can influence coral size structure. Environmental conditions can vary drastically at the scale of the reef, particularly as a function of wave exposure and water depth; sites near the exposed side of a reef usually experience less variability than sheltered sites in a variety of parameters (e.g., temperature, oxygen concentration, etc.)23. Perhaps the most extreme gradient that exists over this small scale is over the shallow reef flat from exposed to sheltered sites. The reef flat is typically the shallowest part of the reef and its height is limited by sea level24. The reef flat habitat is subjected to frequent disturbances like low tides, high temperature variability, and extreme wave action25. Though the benthic community on the reef flat is less diverse than the reef slope, the reef flat habitat is important for trophodynamics of the entire reef system26.
Typically reef flats can be defined as either coral- or rubble-dominated, and within each reef flat clear zonation in the biological community can be observed from the high wave energy windward side of the reef to the low wave energy leeward side24. In the central Red Sea, the branching coral Stylophora pistillata is common in shallow areas and can be found across all reef flat zones27,28,29. Stylophora pistillata has been called an r-strategist because it quickly colonizes unstable environments, reproduces early in its lifetime, grows quickly, and has high population turnover30. These characteristics allow S. pistillata to be dominant in the disturbance-prone reef flat habitat, while in more stable environments it tends to be outcompeted by other coral species30. The reef flat habitat is therefore an ideal system to study S. pistillata size structure, as this coral’s high population turnover is likely to reflect recent natural or anthropogenic disturbance history.
Previous studies in the Red Sea have examined the abundance and size class structure of juvenile pocilliporid corals31 and common coral species over cross-shelf32 and large-scale gradients28. However, despite being the most abundant species on reef flats in the region27, no studies assessing size-class structure have been conducted on S. pistillata populations, and intra-reef comparisons are lacking. Such information is critical for understanding population dynamics in a region that is increasingly affected by coral bleaching events33,34. Here, we conduct a population structure assessment of S. pistillata by evaluating size structure, color morph distribution, and incidence of partial mortality to investigate large and small-scale patterns of variability across the shelf and within each reef across the reef flat.
Materials and methods
Study sites
This study was carried out between 16–20 December 2020, on six reef flats in the central Red Sea over a cross-shelf gradient, with the closest reefs situated ~ 2 km from shore and the farthest ~ 20 km offshore (Fig. 1). We selected two nearshore (Abu Shousha and Tahala North), two midshelf (Al Fahal and Shark Reef), and two offshore (Cement Wreck and Shib Nazar) reefs for this study. The reef flats are characterized by consolidated pavement, rubble, or sand, and are dominated by the branching coral S. pistillata (Fig. 2a). The width of each reef flat (i.e., axis that is perpendicular to shore) varies from 125 to 300 m, and all are at depths between 50 and 150 cm. Within each reef, three sites were chosen to represent the exposed, midreef, and sheltered zones of the reef. The exposed sites receive the most wave action, and are typically on the western edge of the reef flat (within ~ 5 m in front of the breaking waves). The midreef sites are roughly halfway across the reef flat, while the sheltered sites are on the leeward reef edge, generally on the eastern side of the reef (Fig. 1 inset).
Temperature profiles
As water flows across the reef flat, it warms during the day and can result in a temperature gradient from the exposed to sheltered sides35, which may influence the ecology of S. pistillata in different reef zones. As part of a separate ongoing study for monitoring fine-scale differences in temperature across reef flat zones, Onset Hobo pendant temperature loggers were deployed at each reef flat zone of all six reefs, starting in June 2019. The calibrated loggers were programed to log temperature every 10 min, and were secured to lead weights and deployed on the substrate at 50–100 cm depth. Both the weights and loggers were wrapped in white electrical tape to minimize the effects of intense sunlight in shallow water on the logger’s temperature readings. Loggers were swapped out every two to twelve months depending on logistical constraints (e.g., national COVID-19 lockdowns). Temperature was logged for some reefs starting in June 2019 and near-continuous profiles were obtained until the time of the study. Unfortunately, many sites are missing data for long periods of time due to logistical constraints. We therefore chose to focus on the two reefs with the most complete datasets: the offshore reef Shib Nazar and the nearshore reef Abu Shousha. We highlighted a six-week period of 2020 when temperatures were highest and calculated the degree-heating weeks for each reef zone. Degree-heating weeks is a measure of accumulated heat stress commonly used in coral biology36. We used a baseline of 31.8 °C, one degree above the mean maximum for this region, which has been used in previous studies for this area37.
Phototransect surveys and image processing
At each of the reef sites, we conducted three replicate phototransect surveys of 10 m length. Transect locations were haphazardly selected to begin near the temperature logger deployed at each of the sites, and were laid parallel to the reef edge with 3–5 m separating the start and end of each replicate. Photos were taken along the entire length of the transect via snorkel using a Canon G7 X camera (Tokyo, Japan). Photos were later processed in the software ImageJ to count the number and size of S. pistillata colonies38. Size was measured by calibrating each image to a known length (referencing the transect tape in the photo) and using the “area” measurement tool in ImageJ to trace the contours of each colony, resulting in the planar area39 (Fig. 2b). All colonies of S. pistillata that fell completely within the area of the photo were counted and measured; colonies that were partially outside of the photo were not considered in the analysis. Many colonies suffered partial mortality; in these cases, a separate measurement was taken to quantify the amount of living tissue relative to the size of the colony. For measuring abundance, we only considered colonies with greater than 50% living tissue (sensu19). Since the water depth on reef flats varied slightly between reefs and with tidal oscillations, the distance between the camera from the substrate was not equal for all surveys, resulting in differences in the area covered for each photo. To standardize abundance counts between reefs, we measured the area of each photo to obtain the total area surveyed for each site, and expressed abundance in terms of density (individual colonies m−2). For all colonies considered in the analysis (> 50% living tissue), we also quantified partial mortality between sites and reefs. For each colony we also assigned a color morph category: yellow, purple, or mixed colors, if the colony was partially both colors. We also calculated the proportion of juvenile colonies based on previous studies that consider colonies of S. pistillata and related species less than 10 cm22 planar area likely to be juveniles11,40. Finally, summing the total area of all colonies allowed us to calculate the percent cover of S. pistillata for each reef flat zone to compare with colony density.
Data analysis
We plotted the size class structure of each reef site, the total for each reef, and the entire sampled population of S. pistillata in our study system. Colony size was first log transformed to normalize the distributions. We separated size classes into bins of 0.5 on the log scale for visualizing size structure (a dashed line was added at the juvenile adult cut-off, 10 cm241). The descriptive statistics (skewness, kurtosis, mean size) were calculated for each site individually, as well as for each reef and the total population sampled; these analyses were performed in R (R Development Core Team 2018) using the package “moments”42. We compared log-transformed size-frequency distributions using two-sided Kolmogorov–Smirnov tests for reef zones within each reef. Due to extreme variability among reefs, we did not perform these tests for reef flat zones between different reefs. For each of the assessed parameters (e.g., density, mean size, etc.), we performed a non-parametric Kruskal–Wallis analysis to test for differences among reef flat zones, reefs, and shelf position. When significant differences were found, pairwise Mann–Whitney U tests were performed.
A Permutational Multivariate Analysis of Variance (PERMANOVA) was carried out in Primer v7 (Permanova + for Primer43). The model included three factors: shelf (fixed, orthogonal, three levels – nearshore, midshelf, offshore); reef (random, nested in shelf, two levels per shelf position); and reef flat zone (fixed, nested in reef, three levels – exposed, midreef, sheltered). The analysis was conducted on the distance matrix generated using the Euclidean distance and considering mean colony size, colony density, percentage of purple colonies, and percentage of colonies with partial mortality as variables. Percent juveniles were not included in the analysis as this parameter correlated strongly with mean colony size (r2 = 0.846, p < 0.05). Data was normalized and transformed (square root) before the analysis.
Results
Temperature profiles
Overall, temperatures on the offshore reef Shib Nazar were lower than for the nearshore reef Abu Shousha. Comparing 18 months of temperature profiles of just the exposed reef flat zones of each reef shows that, generally, Abu Shousha experienced higher temperatures during peak summer months, but also lower temperatures during the winter months (Fig. 3a). To assess fine-scale variability across the reef flat zones, we focused on the warmest period of 2020, from August to mid October (Fig. 3b,c). Across reef zones within each reef, generally the exposed zone is less variable and has lower maximum temperatures than the midreef and sheltered reef flat zones, but this pattern is less pronounced for Abu Shousha (Fig. 3b) than for Shib Nazar (Fig. 3c). For this period, the mean daily temperature ranges of Abu Shousha are 1.88 ± 0.60°, 1.93 ± 0.80°, and 2.27 ± 0.44 °C (mean ± se) for the exposed, midreef and sheltered reef flat zones, respectively. For Shib Nazar, the values are 1.36 ± 0.96°, 2.46 ± 0.86°, and 2.68 ± 0.71 °C for the exposed, midreef and sheltered reef flat zones, respectively. Despite the differences in temperature variability, the mean temperatures during this period are similar across all reef flat zones of each reef (Abu Shousha: 33.10 ± 0.73°, 32.99 ± 0.73°, 33.11 ± 0.81 °C ; Shib Nazar: 32.61 ± 0.57°, 32.82 ± 0.87°, 32.88 ± 0.92 °C, exposed, midreef, and sheltered, respectively). The accumulated heat stress of Abu Shousha was higher than for Shib Nazar for all reef flat zones (Abu Shousha: 14.09, 12.92, 14.15 °C-weeks ; Shib Nazar: 8.66, 10.79, 11.58 °C-weeks, exposed, midreef, and sheltered, respectively).
Coral size structure
Overall, an area of 279.16 m2 was surveyed in this study (Table 1). A total of 1,146 colonies of S. pistillata were counted across the study sites, ranging in size from 0.14 to 537.39 cm2 planar area (Table 1). The size structure is extremely variable, both across reef flat zones and between reefs, and there are no clear patterns across reef flat zones or reefs in terms of mean size. The highest mean colony size was 153.83 ± 34.32 cm2 at the exposed reef flat zone on Shark Reef, while the smallest mean colony size was found at the exposed reef flat zone of Al Fahal at just 6.37 ± 1.35 cm2 (Table 1), which are both midshelf reefs. Furthermore, within-reef mean colony sizes are variable and there is no clear trend between exposed, midreef or sheltered reef flat zones among reefs. For each individual reef, there were no differences in mean size across the different reef flat zones (Kruskal–Wallis, p > 0.05 for each). Reefs had significantly different mean colony sizes (Kruskal–Wallis, p < 0.05). Abu Shousha has the smallest mean colony size, Tahala North has the largest, and the other reefs had intermediate mean colony sizes. Considering shelf position, the midshelf reefs have significantly smaller mean colony sizes (39.99 ± 4.60 cm2) compared to the offshore reefs ( 60.27 ± 3.87 cm2) and nearshore reefs (50.28 ± 4.99 cm2) (Mann–Whitney U test, p < 0.05 for each). Offshore and nearshore reefs did not have significantly different colony sizes (Mann–Whitney U test, p = 0.583).
Density of S. pistillata was similarly variable over multiple spatial scales, and ranged from 0.44 colonies m−2 at the exposed reef flat zone of Tahala North to 11.93 colonies m−2 at the midreef reef flat zone of Cement Wreck (Table 1). For each individual reef, there were no differences in density across the different reef flat zones (Kruskal–Wallis, p > 0.05 for each). However, reefs did have statistically different densities, with Tahala North and Shark Reef having the lowest densities, Al Fahal, Cement Wreck and Shib Nazar having intermediate densities, and Abu Shousha having the highest density (Kruskal–Wallis, p < 0.05). The offshore and nearshore reefs have similar overall density of S. pistillata (5.85 and 3.93 individuals m−2, respectively; Mann–Whitney U test, p = 0.606), which is higher than the midshelf reefs (2.57 individuals m−2; Mann–Whitney U test, p < 0.05 for each).
Considering size-frequency distributions for each site (Fig. 4) and the entire population sampled in this study (Fig. 5), we found that skewness (g1) ranged from − 0.66 at the Shark Reef exposed zone to 1.46 at the Tahala North exposed zone (Table 1). Abu Shousha showed a negatively skewed distribution for the exposed reef flat zone, but positively distributed for the midreef and exposed reef flat zones. Tahala North showed the opposite trend, with its exposed reef flat zone being the only negatively skewed zone. Al Fahal had positively skewed distributions for all reef flat zones, while the only positively skewed site for Shark Reef was the sheltered reef flat zone. Cement Wreck had negatively skewed distributions for all reef flat zones, while Shib Nazar only had a negative skew for the midreef reef flat zone. Considering the size frequency of entire reefs, two are negatively skewed (Tahala North and Shib Nazar) and the rest are positively skewed. Kurtosis (g2), which is a measure of the effect of outliers on a distribution, was variable among sites but mostly platykurtic (few outliers in the population; g2 < 3) or mesokurtic (roughly approximating a normal distribution; g2 ~ 3). The only cases where sites’ kurtosis was leptokurtic (many outliers in the population; g2 > 3) were the exposed reef flat zones of Tahala North and Shark Reef (Table 1). The Kolmogorov–Smirnov tests revealed that within every reef, there were differences in size-frequency distributions between one reef flat zone with the other two; however, there was no consistency of which reef flat zone was unique (exposed for Abu Shousha, Tahala North and Shark Reef; midreef for Cement Wreck; and sheltered for Al Fahal and Shib Nazar) (Fig. 4).
The proportion of juvenile colonies varied between reefs and reef flat zones. The highest percentage was observed at the exposed zone of Al Fahal, with 86.54% of colonies under 10 cm2 planar area, while the lowest was 5.45% at the midreef zone of Tahala North (p < 0.05, Mann–Whitney U test) (Table 1, Fig. 4). For each individual reef, there were no differences in proportion of juveniles across the different reef flat zones (Kruskal–Wallis, p > 0.05 for each). Reefs did differ in the proportion of juvenile colonies, with Tahala North and Al Fahal having higher proportions than the other reefs (Kruskal–Wallis, p < 0.05). The proportion of juveniles was not significantly different at the level of shelf position (Kruskal–Wallis, p = 0.198).
Size structure compared to coral cover
Summing the planar area of S. pistillata for each transect allows us to compare the percent cover to some aspects of size structure. Colony density is positively correlated with percent cover, but this correlation is weak (r2 = 0.373, p < 0.05) (Fig. 6). Since percent cover is a broad measurement that obscures more detailed information, sites with similar values of percent cover can have substantially different size class distributions. For example, the exposed and midreef reef flat zones of Cement Wreck have almost identical percent cover of S. pistillata (4.19 and 4.18%, respectively). However, the colony density of the exposed zone is nearly three times lower than the midreef zone (4.03 and 11.93 ind. m−2, respectively), and the mean colony sizes differ significantly (104.08 ± 12.29 and 35.04 ± 4.44 cm2, respectively; p < 0.05, Mann–Whitney U test). In addition, their size-frequency distributions are significantly different (Kolmogorov–Smirnov, p < 0.05).
Frequency of colonies with partial mortality
Incidence of partial mortality ranged from 0 to 47% of colonies between reef flat zones (Table 1). There were no differences in frequency of partial mortality among reef flat zones, reefs, or shelf position (Kruskal–Wallis, p > 0.05 for each). Across the entire study there is a trend of increasing proportion of colonies with partial mortality with increasing size class (Fig. 5), reaching 50% for S. pistillata found in the largest size class.
Color morph abundance and distribution
The proportion of purple colonies per site varied widely, from 0% at several sites to 81% in the exposed reef flat zone of Al Fahal. Mixed color morphs are uncommon, only accounting for 0.7% of colonies from three reefs. Overall, yellow colonies dominated reef flats in this study, accounting for nearly 90% of all colonies surveyed (Table 1). For each individual reef, there were no differences in proportion of purple colonies across the different reef flat zones (Kruskal–Wallis, p > 0.05 for each). However, there were differences at the level of reef (Kruskal–Wallis, p < 0.05), with Al Fahal having a higher of purple colonies than all other reefs (Mann Whitney U test, p < 0.05 for each). In additon, no purple colonies were found in Shark Reef, resulting in a lower proportion than all other reefs (Mann Whitney U test, p < 0.05 for each). Though not statistically significant (Kruskal–Wallis, p = 0.058), there is an overall trend of a higher percentage of purple colonies in the exposed zones (18.33%) than the other reef zones (11.54 and 3.80% for sheltered and midreef zones, respectively). There is no difference in the proportion of purple colonies at the level of shelf position (Kruskal–Wallis, p = 0.327).
Population structure
As no broad trends were found in the analysis of coral size structure, a PERMANOVA analysis was conducted including data on partial mortality and color morph assignment. The PERMANOVA analysis showed no differences for the factor shelf, but a significant interaction between the factors reef and reef flat zone, showing a high level of variability within the sampling area. Pair-wise tests of reefs within each shelf position revealed that the offshore reefs Cement Wreck and Shib Nazar are not significantly different, but the nearshore (Abu Shousha and Tahala North) and midshelf (Al Fahal and Shark Reef) reefs are different from each other within groups. Considering reef flat zones within each reef, in Abu Shousha the reef flat zones were not significantly different from each other, while Tahala North, Al Fahal and Cement Wreck showed differences between all reef flat zones. Shark Reef only had significant differences between the exposed and sheltered reef flat zones, while Shib Nazar had significant differences between sheltered and both midreef and exposed reef flat zones (Table 2).
Discussion
Assessments of size class structure can yield insights into coral population status5,6. Using size frequency distributions instead of percent coral cover is more informative for understanding demographic trajectories. In this study, we found large variability of the size class structure of S. pistillata both within and among reefs. We did not observe a clear pattern for either cross-shelf or within-reef scales for any of the parameters analyzed separately (colony densities, size, partial mortality, and color morphs). However, multivariate analysis conducted on these parameters showed that most reefs had significant differences between some or all of their reef flat zones. Our results for two reefs with differing thermal profiles suggest that environmental factors such as temperature may be partly responsible for driving spatial differences in S. pistillata population structure. In the Red Sea, few historical datasets are available, especially outside of the northern Gulf of Aqaba region. Although S. pistillata from this region tend to be resilient to high temperatures, further warming due to climate change may result in temperatures that exceed their thermal tolerance, especially when coupled with local stressors like pollutants44. This demographic “snapshot” will therefore serve as a baseline assessment of this common species, which is essential for future biomonitoring. Such baselines have revealed changes in size structure of Stylophora corals in other regions. For example, in the Great Barrier Reef, the genus Stylophora has experienced major declines in overall abundance over the past two decades, with concomitant shifts in size class distribution10. In Kenya, Stylophora underwent a significant reduction in median colony size following the mass bleaching event of 199812.
Temperature variability in space and time
The within-reef patterns of temperature were broadly similar between the two reefs assessed. The midreef and sheltered reef flat zones tended to have larger daily ranges and higher maximum temperatures than the exposed reef flat zones of both reefs (Fig. 3b,c). In addition, for both reefs, the accumulated heat stress was highest at the sheltered reef flat zone. However, the offshore reef Shib Nazar exhibited a more pronounced difference between the exposed and other reef flat zones, which is reflected in the accumulated heat stress (the midreef and sheltered reef flat zones are 25 and 34% higher, respectively, than the exposed reef flat zone) (Fig. 3c). Interestingly, for Abu Shousha, the midreef reef flat zone had the lowest accumulated heat stress (9% lower than the exposed reef flat zone), while the sheltered reef flat zone was similar to the exposed reef flat zone (0.5% higher). The difference in variability between the two reefs may be due to two phenomena: the shelf position of each reef, and the distance across each reef flat (Shib Nazar, ~ 300 m; Abu Shousha, ~ 140 m). As an offshore reef, Shib Nazar receives wave-driven water flow from deeper, more thermally stable waters, particularly on the exposed side. Abu Shousha, on the other hand, is surrounded by shallow coastal waters that are subjected to greater temperature variability, even on the exposed side of the reef. This is supported by comparing the long-term temperature profiles of each reef’s exposed reef flat zone, where Shib Nazar tends to have less extreme temperatures in both summer and winter (Fig. 3a). In addition, the greater distance across Shib Nazar results in a longer residence time of seawater, allowing for greater diurnal heating potential, while the reef flat of Abu Shousha may not be large enough to induce such variation (i.e., the residence time is too short to cause substantial heating). This may explain why the midreef and sheltered reef flat zones of Shib Nazar tend to have much higher variability compared to the exposed reef flat zone, while all reef flat zones of Abu Shousha exhibit similar variability (Fig. 3b,c).
Size structure
The results presented here offer a snapshot of S. pistillata populations on several reef flats in the central Red Sea. The information contained in the size frequency distributions reveal extreme variability on both inter-and intra-reef scales. We chose to use belt transects as a survey method, which may have a small bias towards counting smaller colonies. However, based on our personal observations, we are confident our results accurately reflect the true size structure. Previous studies of coral populations have also found variable size class structure across sites. Adjeroud et al. (2007) note that although they found significant variation in the size structure of two coral species in Moorea, they found no clear patterns in depth or location that could explain the differences, and suggest that other stochastic processes may be at play.
There are few patterns in any of the measured parameters at the level of reef flat zones, and over the entire population there is a fairly even distribution of all size classes among reef flat zones (Fig. 5). The high overall variability may be indicative of stochastic processes occurring in the area, as has been suggested in previous studies19,20. Considering shelf position, the offshore reefs have the highest density and largest colonies of S. pistillata, perhaps indicating these are more suitable habitats for sustaining larger populations of S. pistillata than the reefs closer to shore (Table 1). A paucity of smaller colonies suggests either a lack of recruitment or high mortality of younger colonies. Unlike many coral species which spawn, S. pistillata is a brooding coral that retains its larvae in the parent colony until planulation, and the majority of planula tend to settle within hours45,46. This suggests S. pistillata may have limited dispersal ability compared to broadcast spawning coral species. Indeed, some studies conducted in the same area as this study have suggested that there is genetic structure of S. pistillata among reefs47. The percentage of juveniles in different reef flat zones seems to change over the cross-shelf gradient, from more juveniles on the exposed reef flat zones of Abu Shousha, Tahala North and Al Fahal, to a lower proportion of juveniles on the exposed reef flat zones of Shark Reef, Cement Wreck and Shib Nazar (Table 1). It is unclear what causes this pattern, or if it simply is a stochastic process. Though not directly measured here, in our experience visiting these sites, the exposed reef flat zones of Shark Reef, Cement Wreck and Shib Nazar tend to experience higher wave action, potentially affecting settlement rates of S. pistillata on these reefs48.
Evaluating the shape (i.e., skewness) of the size class distribution can reveal important details about population structure. Considering the sites individually, half are negatively skewed, indicating a lack of smaller colonies and over-abundance of large colonies (Table 1). However, considering the reef as one entity, only two of six are negatively skewed, and the entire population is slightly positively skewed, showing that the scale of the survey is important for interpreting trends. Negatively skewed populations, with a lower proportion of small individuals, have been suggested to be indicative of a degraded environment14. However, some studies suggest that coral populations in degraded environments can also be positively skewed due to high adult mortality12,28,49. Without repeated monitoring, it is difficult to determine the drivers of skewness observed here.
Two of the reefs in this study are of note and deserve further discussion. Tahala North exhibits different substrate characteristics than the other reefs, with a layer of sand covering the pavement across most of the reef flat. Although S. pistillata is still common across much of the reef, the sandy substrate may influence growth and juvenile mortality, which affects population density. In addition, the exposed zone of Tahala North experiences seasonal Sargassum blooms which can completely cover any corals growing there for several months per year. We speculate the extremely low density of S. pistillata colonies there is due to periodic Sargassum blooms which may negatively affect S. pistillata growth50 or recruitment51. The second reef which deserves more attention is Shark Reef, which also had low overall density of S. pistillata. During the late summer of 2020, a bleaching event affected the reef system in the study region, but Shark Reef was the most impacted reef (pers. obvs.). During an October 2020 research expedition, we noted parts of the reef crest on Shark Reef experienced 100% bleaching (Fig. S2), whereas the other reefs in this study were minimally impacted. While we were not able to assess the reef flat on this expedition, it’s possible that habitat was also heavily impacted and the population of S. pistillata could have experienced high mortality rates. Interestingly, during a mass bleaching event in the central Red Sea in 2015/2016, Shark Reef (refered to as Qita Al-Kirsh) was also heavily impacted by bleaching, moreso than other midshelf and offshore reefs assessed34. This may be an indication that Shark Reef is more susceptible to disturbance, although it is not clear what factors may be responsible.
Frequency of partial mortality
As modular organisms, corals frequently suffer partial mortality52, and the percentage of colonies with partial mortality tends to increase with increasing colony size14,53,54,55. Our results support this pattern, with the smallest size classes containing no colonies suffering partial mortality and an increasing proportion of partial mortality with increasing size class (Fig. 5). Furthermore, the site with the highest percentage of colonies with partial mortality is the exposed reef flat zone of Shark Reef, which also exhibits the largest mean size of any site in the study (Table 1). The largest size class in this study (S. pistillata larger than ~ 400 cm2) reached the highest incidence of partial mortality of 50% (Fig. 5), which is comparable to previous studies on other coral species14. Partial mortality has been shown to have little effect on per-polyp reproductive output56, meaning overall reproduction will decline proportional to tissue loss. As protoanderous simultaneous hermaphrodites, S. pistillata tends to produce only male gametes at smaller sizes (diameter < 4 cm) and produce both male and female gametes at larger sizes (diameter > 4 cm)41. Approximating the colony’s planar area as a circle, this equates to a cutoff of roughly 12.57 cm2, meaning only the smallest adults will be producing male-only gametes. Given that our results show the frequency of partial mortality is very rare in colonies of small size classes, partial mortality is unlikely to drastically affect the ratio of male/female gametes produced by S. pistillata colonies.
Color morphs
Purple colonies accounted for just 10% of S. pistillata colonies over the entire study are. The dominance of yellow colonies on these reef flats is notable, as previous work has shown that purple morphs of S. pistillata typically outcompete the yellow morphs57. However, a related species that also exhibits yellow and purple color morphs (Pocillipora damnicornis) shows a similar pattern of yellow morphs dominating shallow areas and purple colonies more common below 3 m depth58. The spatial distribution of color morphs of other coral species are often related to the light environment59, and can influence disease or bleaching susceptibility60,61. However, the biological significance of color morphs in S. pistillata is not clear. Studies on the photochemistry of the pocilliporin pigments, which give the purple morphs their characteristic color, showed no evidence of photoprotection, function as an accessory pigment, or UV protection62. It may be that the yellow morphs of S. pistillata are better adapted to warmer temperatures, higher irradiances, or some other factor unique to reef flat environments, and the purple morphs are competitively dominant only in deeper areas. Color morphs may also harbor different symbiont communities, which may have different tolerances to temperature and irradiance59. Although S. pistillata from the Red Sea tend to exhibit a depth zonation in their Symbiodiniaceae communities (Symbiodinium in shallow habitats and Cladacopium in deeper habitats), a recent study conducted in the same area as our study found that both genera can be present in shallow-water colonies63. Regardless, the high percentage of purple morphs on Al Fahal is notable. The exposed reef flat zone is dominated by juvenile purple S. pistillata, which could indicate a recent recruitment pulse for that site. However, there is also a high proportion of adult purple colonies across all reef flat zones of Al Fahal. The higher abundance of purple morphs here may be the result of other stochastic processes.
Multivariate analysis
While each of the individual factors examined here do not show any clear spatial patterns, the PERMANOVA analysis reveals most reefs have significantly different population structures across their reef flat zones. The different pattern of within-reef temperature variability of Abu Shosuha and Shib Nazar may help explain some patterns observed in the populations of S. pistillata at these sites. The multivariate analysis revealed that Abu Shousha was the only reef assessed which did not have significant differences between its reef flat zones. While other factors are likely also responsible, the similarity in temperature profiles between Abu Shousha’s reef flat zones may indicate that the environmental conditions are similar across the entire reef flat, which is reflected in the similarity of populations of S. pistillata. On the other hand, the sheltered reef flat zone of Shib Nazar was different from both the midreef and exposed reef flat zones (Table 2). Unfortunately we lack contemporaneous temperature data for the other reefs in this study, making it difficult to ascertain the importance of temperature on population structure at each site. However, these two contrasting patterns in terms of temperature variability within a reef flat indicates that environmental conditions can vary at small spatial scales, and may partly explain the lack of biological patterns observed in the asssessment of the population structure of S. pistillata in this study.
Conclusion
As the dominant branching coral species inhabitating reef flats in the region, S. pistillata plays an important role in the ecology and geology of this sometimes overlooked habitat. Here we show that there is high variability over small spatial scales of size structure of this species, perhaps in large part driven by stochastic processes. It is also important to highlight that a multivariate analysis can reveal differences in populations which are not apparent in analyzing each parameter seperately. Our comparison of percent cover to colony density shows that the former metric obscures detailed information about the state of a population that can be gleaned from size structure data. Though our surveys reveal the offshore reefs may have “healthier” (higher colony density and mean colony size) populations of S. pistillata, future monitoring is needed to track population changes and relate them to potential disturbances (e.g., bleaching events, pollution, etc.). We hope the data presented here will foster future studies addressing which other abiotic processes and variables may play a role in this small-scale variability of size structure. Such information may help explain why some reefs experience mass bleaching of corals while other nearby reefs are largely spared. Finally, evaluating the role of biotic interactions (e.g., symbiodiniaceae community, microbiome, predators/herbivores, competition for space, etc.) at different spatial scales may also shed light on differences in population structure of S. pistillata in the region.
Data availability
The datasets used during the current study are available from the corresponding author on reasonable request.
References
Reaka-Kudla, M. L. The global biodiversity of coral reefs: a comparison with rain forests. Biodivers. II. Underst. Prot. Our Biol. Resour. 2, 551 (1997).
Connell, J. H. Population ecology of reef-building corals. in Biology and Geology of Coral Reefs (eds. Jones, O. A. & Endean, R.) 205–245 (Academic Press, 1973). doi:https://doi.org/10.1016/B978-0-12-395526-5.50015-8.
Berumen, M. L. et al. The status of coral reef ecology research in the Red Sea. Coral Reefs 32, 737–748 (2013).
Hughes, T. P., Graham, N. A., Jackson, J. B., Mumby, P. J. & Steneck, R. S. Rising to the challenge of sustaining coral reef resilience. Trends Ecol. Evol. 25, 633–642 (2010).
Edmunds, P. J. & Riegl, B. Urgent need for coral demography in a world where corals are disappearing. Mar. Ecol. Prog. Ser. 635, 233–242 (2020).
Pisapia, C. et al. Projected shifts in coral size structure in the Anthropocene. Adv Mar Biol 87, 31–60 (2020).
Meesters, E. et al. Colony size-frequency distributions of scleractinian coral populations: spatial and interspecific variation. Mar. Ecol. Prog. Ser. 209, 43–54 (2001).
Riegl, B. et al. Demographic mechanisms of reef coral species winnowing from communities under increased environmental stress. Front. Mar. Sci. 4, 344 (2017).
Pisapia, C., Burn, D. & Pratchett, M. Changes in the population and community structure of corals during recent disturbances (February 2016-October 2017) on Maldivian coral reefs. Sci. Rep. 9, 1–12 (2019).
Dietzel, A., Bode, M., Connolly, S. R. & Hughes, T. P. Long-term shifts in the colony size structure of coral populations along the Great Barrier Reef. Proc. R. Soc. B 287, 20201432 (2020).
Lachs, L. et al. Linking population size structure, heat stress and bleaching responses in a subtropical endemic coral. Coral Reefs 40, 777–790 (2021).
McClanahan, T., Ateweberhan, M. & Omukoto, J. Long-term changes in coral colony size distributions on Kenyan reefs under different management regimes and across the 1998 bleaching event. Mar. Biol. 153, 755–768 (2008).
Grimsditch, G. et al. Variation in size frequency distribution of coral populations under different fishing pressures in two contrasting locations in the Indian Ocean. Mar. Environ. Res. 131, 146–155 (2017).
Bak, R. P. & Meesters, E. H. Coral population structure: the hidden information of colony size-frequency distributions. Mar. Ecol. Prog. Ser. 162, 301–306 (1998).
Hughes, T. & Jackson, J. Do corals lie about their age? Some demographic consequences of partial mortality, fission, and fusion. Science 209, 713–715 (1980).
Hughes, T. P. & Jackson, J. Population dynamics and life histories of foliaceous corals. Ecol. Monogr. 55, 141–166 (1985).
Soong, K. Colony size as a species character in massive reef corals. Coral Reefs 12, 77–83 (1993).
Bak, R. P. & Meesters, E. H. Population structure as a response of coral communities to global change. Am. Zool. 39, 56–65 (1999).
Adjeroud, M., Pratchett, M. S., Kospartov, M. C., Lejeusne, C. & Penin, L. Small-scale variability in the size structure of scleractinian corals around Moorea, French Polynesia: patterns across depths and locations. Hydrobiologia 589, 117–126 (2007).
Adjeroud, M., Mauguit, Q. & Penin, L. The size-structure of corals with contrasting life-histories: A multi-scale analysis across environmental conditions. Mar. Environ. Res. 112, 131–139 (2015).
Bauman, A. G. et al. Variation in the size structure of corals is related to environmental extremes in the Persian Gulf. Mar. Environ. Res. 84, 43–50 (2013).
Smith, L., Devlin, M., Haynes, D. & Gilmour, J. A demographic approach to monitoring the health of coral reefs. Mar. Pollut. Bull. 51, 399–407 (2005).
Lowe, R. J. & Falter, J. L. Oceanic forcing of coral reefs. Annu. Rev. Mar. Sci. 7, 43–66 (2015).
Thornborough, K., Davies, P. Reef flats. Encycl. Mod. Coral Reefs 869–876 (2011).
Camp, E. F. et al. The future of coral reefs subject to rapid climate change: lessons from natural extreme environments. Front. Mar. Sci. 5, 4 (2018).
Bellwood, D. R. et al. The role of the reef flat in coral reef trophodynamics: Past, present, and future. Ecol. Evol. 8, 4108–4119 (2018).
Pineda, J. et al. Two spatial scales in a bleaching event: Corals from the mildest and the most extreme thermal environments escape mortality. Limnol. Oceanogr. https://doi.org/10.4319/lo.2013.58.5.1531 (2013).
Riegl, B. M., Bruckner, A. W., Rowlands, G. P., Purkis, S. J. & Renaud, P. Red Sea coral reef trajectories over 2 decades suggest increasing community homogenization and decline in coral size. PLoS ONE 7, e38396 (2012).
Riegl, B., Berumen, M. & Bruckner, A. Coral population trajectories, increased disturbance and management intervention: A sensitivity analysis. Ecol. Evol. https://doi.org/10.1002/ece3.519 (2013).
Loya, Y. The red sea coral Stylophora pistillata is an r strategist. Nature https://doi.org/10.1038/259478a0 (1976).
Lozano-Cortés, D. F. & Berumen, M. L. Colony size-frequency distribution of pocilloporid juvenile corals along a natural environmental gradient in the Red Sea. Mar. Pollut. Bull. 105, 546–552 (2016).
Ellis, J. et al. Cross shelf benthic biodiversity patterns in the Southern Red Sea. Sci. Rep. 7, 1–14 (2017).
Furby, K. A., Bouwmeester, J. & Berumen, M. L. Susceptibility of central Red Sea corals during a major bleaching event. Coral Reefs 32, 505–513 (2013).
Monroe, A. A. et al. In situ observations of coral bleaching in the central Saudi Arabian Red Sea during the 2015/2016 global coral bleaching event. PLoS ONE 13, e0195814 (2018).
Davis, K. et al. Observations of the thermal environment on Red Sea platform reefs: A heat budget analysis. Coral Reefs 30, 25–36 (2011).
Liu, G. et al. Reef-scale thermal stress monitoring of coral ecosystems: new 5-km global products from NOAA Coral Reef Watch. Remote Sens. 6, 11579–11606 (2014).
Voolstra, C. R. et al. Standardized short-term acute heat stress assays resolve historical differences in coral thermotolerance across microhabitat reef sites. Glob. Change Biol. https://doi.org/10.1111/gcb.15148 (2020).
Abràmoff, M. D., Magalhães, P. J. & Ram, S. J. Image processing with ImageJ. Biophotonics Int. 11, 36–42 (2004).
Morais, J., Morais, R. A., Tebbett, S. B., Pratchett, M. S. & Bellwood, D. R. Dangerous demographics in post-bleach corals reveal boom-bust versus protracted declines. Sci. Rep. 11, 1–7 (2021).
Hall, V. R. & Hughes, T. P. Reproductive strategies of modular organisms: Comparative studies of reef-building corals. Ecology https://doi.org/10.2307/2265514 (1996).
Rinkevich, B. & Loya, Y. Reproduction of the Red Sea coral Stylophora pistillata. 2. Synchronization in breeding and seasonality of planulae shedding. Mar. Ecol. Prog. Ser. 1, 145–152 (1979).
Komsta, L. & Novomestky, F. Moments, cumulants, skewness, kurtosis and related tests. R Package Version 14, (2015).
Anderson, M., Gorley, R. & Clarke, K. PERMANOVA+ for PRIMER: guide to software and statistical methods. Primer-E Plymouth UK (2008).
Meziere, Z. et al. Stylophora under stress: A review of research trends and impacts of stressors on a model coral species. Sci. Total Environ. 816, 151639 (2022).
Rinkevich, B. & Loya, Y. Reproduction of the Red Sea coral Stylophora pistillata 1. Gonads and planulae. Mar. Ecol. Prog. Ser. 1, 133–144 (1979).
Nishikawa, A., Katoh, M. & Sakai, K. Larval settlement rates and gene flow of broadcast-spawning (Acropora tenuis) and planula-brooding (Stylophora pistillata) corals. Mar. Ecol. Progress Ser. https://doi.org/10.3354/meps256087 (2003).
Monroe, A. Genetic differentiation across multiple spatial scales of the Red Sea of the corals Stylophora pistillata and Pocillopora verrucosa. M.S. thesis, King Abdullah University of Science and Technology, Thuwal, Saudi Arabia (2015).
Gouezo, M. et al. Relative roles of biological and physical processes influencing coral recruitment during the lag phase of reef community recovery. Sci. Rep. 10, 1–12 (2020).
Boco, S. R., Cabansag, J. B. P., Jamodiong, E. A. & Ticzon, V. S. Size-frequency distributions of scleractinian coral (Porites spp.) colonies inside and outside a marine reserve in Leyte Gulf, central Philippines. Reg. Stud. Mar. Sci. 35, 101147 (2020).
River, G. F. & Edmunds, P. J. Mechanisms of interaction between macroalgae and scleractinians on a coral reef in Jamaica. J. Exp. Mar. Biol. Ecol. 261, 159–172 (2001).
Kuffner, I. B. et al. Inhibition of coral recruitment by macroalgae and cyanobacteria. Mar. Ecol. Prog. Ser. 323, 107–117 (2006).
Hughes, T. & Jackson, J. Do corals lie about their age? Some demographic consequences of partial mortality, fission, and fusion. Science 209, 713–715 (1980).
Lewis, J. B. Abundance, distribution and partial mortality of the massive coral Siderastrea siderea on degrading coral reefs at Barbados West Indies. Mar. Pollut. Bull. 34, 622–627 (1997).
Meesters, E. H., Wesseling, I. & Bak, R. P. Coral colony tissue damage in six species of reef-building corals: partial mortality in relation with depth and surface area. J. Sea Res. 37, 131–144 (1997).
Meesters, E. H., Wesseling, I. & Bak, R. P. Partial mortality in three species of reef-building corals and the relation with colony morphology. Bull. Mar. Sci. 58, 838–852 (1996).
Graham, J. & Van Woesik, R. The effects of partial mortality on the fecundity of three common Caribbean corals. Mar. Biol. 160, 2561–2565 (2013).
Rinkevich, B. & Loya, Y. Intraspecific competitive networks in the Red Sea coral Stylophora pistillata. Coral Reefs https://doi.org/10.1007/BF00571193 (1983).
Takabayashi, M. & Hoegh-Guldberg, O. Ecological and physiological differences between two colour morphs of the coral Pocillopora damicornis. Mar. Biol. 123, 705–714 (1995).
Innis, T., Cunning, R., Ritson-Williams, R., Wall, C. & Gates, R. Coral color and depth drive symbiosis ecology of Montipora capitata in Kāne ‘ohe Bay, O ‘ahu, Hawai ‘i. Coral Reefs 37, 423–430 (2018).
Gochfeld, D., Ansley, M., Ankisetty, S. & Aeby, G. Antibacterial chemical resistance to disease in the Hawaiian coral Montipora capitata. Planta Med. 80, CL31 (2014).
Shore-Maggio, A., Callahan, S. M. & Aeby, G. S. Trade-offs in disease and bleaching susceptibility among two color morphs of the Hawaiian reef coral Montipora capitata. Coral Reefs 37, 507–517 (2018).
Dove, S. G., Takabayashi, M. & Hoegh-Guldberg, O. Isolation and partial characterization of the pink and blue pigments of pocilloporid and acroporid corals. Biol. Bull. 189, 288–297 (1995).
Hume, B. C. C., Mejia-Restrepo, A., Voolstra, C. R. & Berumen, M. L. Fine-scale delineation of Symbiodiniaceae genotypes on a previously bleached central Red Sea reef system demonstrates a prevalence of coral host-specific associations. Coral Reefs https://doi.org/10.1007/s00338-020-01917-7 (2020).
Acknowledgements
This study was performed during a research cruise operated by Dream Divers aboard the Dream Master vessel. We thank Dominador Osido, Fernando Canadilla, Sammy Tayo, Benjamin Mascarenhas, Hilario Abidan, and Romeo Lagrason of Dream Divers for their assistance during the cruise. We also thank Ute Langer for creating the map and two anonymous reviewers whose comments improved the quality of the manuscript.
Author information
Authors and Affiliations
Contributions
W.A.R. conceived and designed the study, oversaw field data acquisition, assisted in image analysis, analyzed data, made figures 2-5, and wrote the main manuscript text. S.C. contributed to study design, assisted writing the first drafts of the manuscript, performed the PERMANOVA analysis, and supervised the project. R.C., G.G., and K.G. collected data in the field and assisted in image analysis. M.L.B. contributed to study design, supervised the project, and provided funding. All authors reviewed the manuscript and contributed to the final version.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Publisher's note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Rich, W.A., Carvalho, S., Cadiz, R. et al. Size structure of the coral Stylophora pistillata across reef flat zones in the central Red Sea. Sci Rep 12, 13979 (2022). https://doi.org/10.1038/s41598-022-17908-3
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41598-022-17908-3
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.