Abstract
Thermal-stress events are changing the composition of many coral reefs worldwide. Yet, determining the rates of coral recovery and their long-term responses to increasing sea-surface temperatures is challenging. To do so, we first estimated coral recovery rates following past disturbances on reefs in southern Japan and Western Australia. Recovery rates varied between regions, with the reefs in southern Japan showing more rapid recovery rates (intrinsic rate of increase, r = 0.38 year−1) than reefs in Western Australia (r = 0.17 year−1). Second, we input these recovery rates into a novel, nonlinear hybrid-stochastic-dynamical system to predict the responses of Indo-Pacific coral populations to complex inter-annual temperature cycles into the year 2100. The coral recovery rates were overlaid on background increases in global sea-surface temperatures, under three different climate-change scenarios. The models predicted rapid recovery at both localities with the infrequent and low-magnitude temperature anomalies expected under a conservative climate-change scenario, Representative Concentration Pathway (RCP) 4.5. With moderate increases in ocean temperatures (RCP 6.0) the coral populations showed a bimodal response, with model runs showing either recovery or collapse. Under a business-as-usual climate-change scenario (RCP 8.5), with frequent and intense temperature anomalies, coral recovery was unlikely.
Similar content being viewed by others
Introduction
The rapid rate of contemporary climate change is seriously affecting terrestrial and marine populations1,2, particularly on tropical coral reefs where corals have been living close to their upper thermal limits for millennia3,4,5. Although corals on tropical reefs are adapted to warm waters, they rarely experience an annual temperature range greater than a few degrees Celsius. Warmer than average temperatures combined with high seasonal irradiance can cause a dysfunction in the coral-dinoflagellate symbiosis that leads to coral bleaching, and under extreme conditions leads to mortality6,7,8,9,10. There is, however, considerable variability in the range of thermal tolerances among the eight hundred or more extant coral species11,12,13, and there are substantial differences in the rates of population and community recovery from thermal disturbances14,15.
The rate of recovery from a disturbance is a function of a system’s stability and resilience16,17. Stable, resilient systems recover rapidly; whereas unstable, non-resilient systems recover slowly, or even collapse after disturbances18. Moreover, the slowing of recovery usually indicates that a system is deteriorating and may be approaching a critical threshold, beyond which the system switches to an alternate and often undesirable state19,20. Indeed, understanding the rates of change and the resilience of systems has become central to our understanding of modern ecology21,22.
On modern coral reefs, we frequently ask the question: how quickly will a reef recover from a given disturbance? It is challenging to answer this question, particularly given the complexity of coral reefs and the multitude of nuances that influence recovery23,24. Recovery depends on many interacting factors, including species composition and environmental and geographic vicissitudes. Understanding these factors and accurately predicting the recovery rates of coral populations is critical in a rapidly changing climate25, which is characterized by ocean warming and an increase in the frequency of acute thermal-stress events5. There is, therefore, a need for models that accurately predict the dynamics of coral populations and determine the likelihood of recovery under future climate-change scenarios RCP 4.5, 6.0, and 8.526.
Many models explore the dynamics of coral populations25,26,27,28, although most models are built on the traditional foundation of Leslie and Lefkovitch matrices29,30,31. The matrix approach is convenient, but does not capture the instantaneous dynamics of a population. A system of differential equations captures the dynamics of change32,33, but rarely includes environmental stochasticity and ecological uncertainty. Nonlinear stochastic models, however, do allow for nuances involving uncertain return-periods of thermal events. Here we constructed a nonlinear hybrid-stochastic-dynamical system to model the responses of coral populations to the complexities of sea-surface temperature cycles, which vary stochastically in both frequency and intensity. We estimated recovery rates from disturbances to the coral reefs in southern Japan and Western Australia, which have similar coral-species composition. We then input these recovery rates into a generic model that predicted disturbance and recovery of Indo-Pacific coral populations into the year 2100, superimposed on increasing global sea-surface temperatures.
Methods
Following a disturbance on a coral reef, and in the absence of further extreme environmental conditions, the change in average precentage coral cover (P) through time (Fig. 1) often can be approximated by the logistic growth equation. To add environmental complexity and stochastic extremes we considered that the physiological response of corals depends on temperature stress, which in turn depends on irradiance32. A newly constructed hybrid stochastic model, given below, governs the dynamics of that response and the occurrences of extreme-temperature events:
where r is the intrinsic rate of increase in coral cover per year, K is the steady-state coral cover equilibrium point (%), T is the sea-surface temperature (°C), and where tk’s are the years when extreme thermal-stress events occur. N is an integer that indicates the number of years we wish to run the simulations, and in our case it was through to the year 2100. In addition, h is the duration of a thermal-stress anomaly (weeks), which was extracted from a truncated gamma distribution with a minimum of 2 weeks and a maximum of 8 weeks, and γ is a coefficient that impacts coral populations through temperature changes, which affects percentage coral cover. During an extreme-temperature event, coral cover is modeled by Equation 1b, where ε, which we call an ‘extreme-event coefficient’, is the sudden increase in temperature, herein governed by short-term climate oscillations such as El Niño events, which result in mass coral bleaching and mortality. The function f(ε) controls the dynamics, specific to each coral species or coral assemblage, depending on sudden temperature increases (measured by ε). Since the coral populations decrease during temperature anomalies, f(ε) is a negative function for all times. Using the present biological data the function was estimated as −ε3. These values would most likely vary geographically and for different coral assemblages, and could vary into the future, with anomalous temperature events not necessarily coupled with El Niño cycles5, although the generic construct of the model will remain the same.
Sea-surface temperature (in Celsius), T(t), was obtained by curve fitting satellite data (see below), using the following equation:
where Iave is the annual average of irradiance (in photosynthetic available radiation, PAR), and the parameters a1 and a2 are the rates at which irradiance affects the change in temperature (or the rate of change in Celsius per PAR) (in Celsius/PAR). We assume that seasonal irradiance does not vary among years. Here, λ is defined as the ‘climate-change coefficient’ (Celsius time−1), where λt controls the temperature increase over time, in years, and a3 is the annual average of sea surface temperature (in Celsius) for λ = 0 under normal, non-anomalous years (Table 1). To estimate these parameters we used nonlinear optimization using the program Mathematica®.
Biological data
To estimate the parameters in Equation 1, we first estimated recovery rate (r) and carrying capacity (K). We obtained data sets of coral cover from two geographic regions: Sesoko Island, Okinawa, in southern Japan11,14, and Scott Reef in Western Australia15. Total coral cover estimates were obtained from Sesoko Island from 1997–201011,14. The corals at Sesoko Island experienced an anomalous thermal-stress event in 1998, which was driven by the El Niño Southern Oscillation, which reduced coral populations significantly11. The Sesoko Island reefs also experienced a milder thermal-stress event in 2001 (Fig. 1). In Western Australia, Gilmour et al., see ref.15, collected data on total cover, and the cover of dominant reef-building genera, including Acropora and Porites, from Scott Reef from 1994–2010 (Fig. 1). The reefs also suffered coral losses during the 1998 anomalous thermal-stress event (Fig. 1).
We could not, however, with any degree of certainty, simultaneously fit all the unknown parameters in Equation 1a and also predict coral cover into the future. Therefore, we used a two-step process. First, we estimated r, K and γ as inverse problems using a Bayesian platform. The values from the Bayesian posterior distributions were then fitted to parametric distributions. Second, we used these best-fit distributions to predict the response of coral populations to likely thermal-stress scenarios (see below) through to the year 2100. We estimated r and K values (Equation 1a) at each study locality as an inverse problem using the logistic growth function within a Bayesian framework using a lag-1 temporal auto-regressive error structure on the residuals (after ref.34). The successive r-values for the different localities were derived from a sample size of 50,000 from the posterior distribution, estimated using OpenBUGS. These r-values, from the posterior distribution, were low, <0.4, and were therefore fitted to Beta distributions using the R package ‘fitdistrplus’35. Random samples were taken from these estimated Beta distributions for every time step in the predictive model. We used a similar process to estimate γ-values using Equation 3 below. Although temperatures change through the seasons, in the tropics they only range about 3 °C, which we assumed was small enough to fix T, and thereby use mean temperature and the solution of Equation 1a (i.e., Equation 3) to estimate γ-values. Therefore, to derive a meaningful estimate of γ, we treated temperature,T, and K as constants, and treated r as a stochastic parameter, using the following:
Temperatures, Thermal stress, El Niño, and Climate Change
Sea-surface temperature (SST) data were obtained for each reef from the Moderate-resolution Imaging Spectroradiometer (MODIS) on board the Aqua satellite platform that was launched in 200236. Average monthly MODIS Aqua Global Level 3 mapped mid-infrared nighttime SSTs from January 2003 (the earliest available date) through December 2010 were obtained at a 4.6 km spatial resolution for each locality (Sesoko Island 26.65854°N 127.8611°E; Scott Reef 14.167959°S, 121.87968°E). Monthly SST records were extracted from each study site and mean monthly SST temperature was calculated.
To estimate the distribution of the frequency of thermal-stress events, we obtained Oceanic Niño Index (ONI) data from 1950–2016. The ONI is the standard metric used to identify El Niño (warm) and La Niña (cool) events in the tropical Pacific, and is a running 3-month average sea-surface temperature anomaly for the Niño 3.4 region37 (5°N–5°S and 120°W–170°W), based on the Extended Reconstructed Sea Surface Temperature (ERSST v4) dataset from the National Oceanic and Atmospheric Administration’s National Climate Data Center. These data were derived from the International Comprehensive Ocean-Atmosphere Dataset at a 2° by 2° spatial resolution (https://www.ncdc.noaa.gov). El Niño events were defined as five consecutive three-month periods, at or above the 0.5 °C anomaly. We categorized El Niño events into one of four strengths based on the magnitude of the anomaly: (1) weak events were characterized as 0.5–0.9 °C anomalies, (2) moderate events were defined as 1.0–1.4 °C anomalies, (3) strong events were defined as 1.5–1.9 °C anomalies, and (4) very strong events were greater than or equal to 2.0 °C anomalies (Table 2). The frequency of thermal-stress events was then estimated for each of four categories: (1) all events; (2) only moderate, strong, and very strong events; (3) only strong and very strong events; and (4) only very strong events (Fig. 2).
To estimate the distribution of the frequency of thermal-stress events of each category of El Niño, the number of years in between the start of each event was calculated. The mean, median, and standard deviation of the period in between the anomalous-thermal events were also calculated (Table 2). The frequency of extreme events in the model simulations was drawn from a Poisson distribution where the rate parameter was estimated from the mean of the frequency of thermal-stress events of categories 1–3 (i.e., between 3 and 9 years). The historical record shows that very strong thermal-stress events occurred every 21.7 ± 9 years (median = 18) (Fig. 2); therefore very strong events were excluded from this study because the infrequent return period lowers the confidence in these events for our relatively short timeframe (to the year 2100).
It is expected that, by 2100, sea-surface temperatures will increase by approximately 1.8 °C, 2.2 °C and 3.7 °C for the Representative Concentration Pathways (RCPs) 4.5, 6 and 8.5 respectively38. Therefore, the maximum temperature anomaly was set to 3.7 °C bove the mean background temperatures in the models38. Based on the RCP predictions, the corresponding climate-change coefficient values, λ, in Equation 2 were re-calculated using curve fitting (Table 1).
Model implementation
The intervals between thermal-stress events were extracted randomly from a Poisson distribution, with values selected between 3- and 9-year intervals (Table 2), and the duration of each thermal-stress event was randomly selected from a gamma distribution lasting anywhere from 2 to 8 weeks over any given summer period. Because seasonal irradiance is unlikely to vary among years39, Iave was set at 45 in Equation 2 for all iterations. There was a minimal effect of the λ parameter (in Equation 2) on temperature increases between the years 2000 and 2010, therefore λ (Equation 2) was set to zero for the initial iteration (Table 1). The algorithm used for the model is available in Appendix 1; all the R code and the OpenBugs code is available in the supplementary document.
Results
The estimated carrying capacities (K) were similar between reefs, at ~53% in southern Japan and 64% in Western Australia, although uncertainty, that is the 95% credible intervals, was greater for the reefs in southern Japan than for the reefs in Western Australia (Table 3). The recovery rates, or the intrinsic rates of increase, of the coral populations were higher in southern Japan (r = 0.39) than in Western Australia (r = 0.17), again with considerably greater uncertainty for southern Japan than for Western Australia (Table 3). The recovery rates of Acropora (r = 0.19) were higher than the recovery rates of Porites (r = 0.09) in Western Australia, with similar uncertainty intervals for each genus (Table 4). Under modern temperature conditions with intermittent temperature anomalies, the model predicted that both geographic localities maintained relatively high coral cover, near their carrying capacities (Fig. 3). The predictions for Acropora and Porites varied, with Acropora populations showing greater fluctuations through time than Porites (Fig. 4) in Western Australia.
Considering a negative effect of high sea-surface temperature on coral populations, the λ-values (Equation 2), or the ‘climate-change coefficients’, were estimated at 0.015, 0.021, and 0.039 (Table 1), for RCP 4.5, RCP 6.0, and RCP 8.5, respectively. Increases to the thermal susceptibility, γ, in Equation 1a, resulted in reductions in overall coral cover (compare Fig 3a to 5). Similarly, increasing the ‘extreme-event coefficient’, resulted in bimodal responses, with some model runs predicting population collapse, whereas other runs predicted recovery (Fig. 6). When the frequency of thermal-stress events was increased from 3–9 years to 3–6 years, the corals did not recover, particularly when thermal anomalies also increased in intensity (Fig. 7a), as is expected under a business-as-usual climate-change scenario (RCP 8.5). It is notable that when the two localities, with different coral recovery rates, were modeled, they showed different responses; Western Australia, with lower coral recovery rates (r = 0.17), completely collapsed under high frequency and intensity temperature anomalies (Fig. 7b), whereas reefs in southern Japan, with considerably higher coral recovery rates (r = 0.39), maintained low coral cover (Fig. 7a).
Discussion
Forecasting when and where thermal stresses will occur is particularly challenging because of stochastic fluctuations in short-term inter-annual climatic cycles, and because of regional differences in both the ocean temperature and the rates of change in ocean temperature caused by global warming38,39. Still, while such short-term climate events are unpredictable in any given year, they are predictable over the course of a decade or more. Our model captured that predictability by using Poisson distributions that randomly selected years with anomalous thermal stress. Although the actual years when thermal-anomaly events will take place in the coming century are unknown, the intensity and frequency of the events used in our models are realistic enough to predict the responses of the coral populations through time.
Our nonlinear hybrid-stochastic-dynamical system allowed us to predict and compare a suite of coral-population responses under various climate-change scenarios. The resultant model predicted that coral populations responded more strongly to intensifying stochastic temperature anomalies than gradual increases in background global temperatures. Three parameters consistently affected the coral-population estimates in our model: (1) ε, which is the ‘extreme-event coefficient’, or the magnitude of the temperature anomaly; (2) r, which determines the rate of recovery of each coral population; and (3) γ, which is the impact rate of background temperature on the coral population.
Firstly, the coral-population models were particularly sensitive to the intensity of temperature anomalies (ε). Coral bleaching occurs during the most severe thermal-stress events10 when summer temperatures are exacerbated by climatic phenomenon, particularly El Niño events; although recently, thermal-stress events have been also occurring during non-El Niño years5. Recovery rates dampened through time when the intensity (ε) and the frequency of thermal-stress events increased, which are predicted to occur under the business-as-usual climate-change scenario RCP 8.538. Our results showed that after intense and frequent thermal-stress events, the coral populations often hovered around <10% coral cover, from which the populations were unlikely to recover. Extremely low coral cover may reduce larval supply, which would further reduce recruitment through the Allee effect, with fitness declining in step with declines in population size.
Secondly, the data used to parameterize the model showed that coral populations in the two locations recovered at different rates. Coral recovery is dependent on numerous conditions, including recent history and what remains on the reef, and the period of recovery between thermal-stress events25. Indeed, the capacity of a population to regrow from remnant fragments has a considerable impact on the recovery rates of coral populations40. For example, small isolated remnants of massive Porites species regrow rapidly after disturbances41, but isolated remnants of corymbose Acropora species rarely regrow27. Consequently, corymbose Acropora populations are mostly dependent on larval recruitment for population recovery, whereas massive Porites can regrow from local remnants. These differences were reflected in the high variance in the Acropora response, compared with the lower variance of Porites response (Fig. 4). Recovery on a particular reef also depends upon the location of the nearest larval source. Scott Reef is considerably more isolated than Sesoko Island, which is adjacent to Okinawa where remnant coral populations survived through the 1998-thermal stress event. Isolation may have played a role in the differences in recovery rates between Scott Reef and Sesoko Island.
Thirdly, our model showed that populations with high γ values were more sensitive to high temperatures than populations with low γ values (Fig. 5). These results are not surprising since reef ecologists have long known that coral species differ in their tolerance to temperature anomalies11. Still, γ-values were strongly coupled with r-values — the two are not independent of each other. In fact, if we assume that T is constant, Equation 1a has two equilibrium points, at P1 = 0 and at P2=\(\,\frac{K}{r}\,(r-\gamma T).\) On the one hand, P2 must be positive for the population to survive, but when \(\frac{\gamma T}{r}=1\), then P2=0, which effectively becomes P1, and is therefore a bifurcation point. On the other hand if \(\frac{\gamma T}{r} < 1\), then P2 > 0, the population is positive and at a stable equilibrium. Therefore, the long-term response of a population to temperature is dependent on the dynamics between the population’s sensitivity to temperature and its capacity to recover.
The model also showed considerable fluctuations in Acropora populations compared with the more stable Porites populations. Further studies on thermal tolerance may lead to more realistic and dynamic distributions of thermal tolerance γ values, which may even vary geographically. For example, modeling coral populations on isolated reefs with low genetic diversity42 will most likely predict leptokurtic distributions for γ values, with little variance, which may project considerably worse outcomes for those populations than for populations with high genotypic diversity.
Although carrying capacities, or equilibrium points, on coral reefs have rarely been quantified and discussed, coral-reef habitats do vary in their carrying capacity. For example, Gouezo et al., ref.43, recently showed long-term stability on the reefs of Palau, with the greatest differences in coral cover being apparent among habitats. Over a 15-year period, the nearshore reefs of Palau at 3 m depth supported ~50% coral cover, whereas the outer reefs supported on average 40% coral cover, and the patch reefs supported on average 20% coral cover43. In the present study, the carrying capacities (K) were similar, but were higher in Western Australia (64%) than in southern Japan (53%), although there was considerable variance around the means. This similarity should not be entirely surprising given the considerable overlap of the coral species composition, and the similarity of habitats between regions (van Woesik, pers. obs.). Yet, carrying capacities depend on numerous interacting processes, including exposure to waves and sedimentation rates. Although these processes have been rarely explored, we need further studies that more accurately quantify carrying capacities of coral reefs across habitats, and that determine the extent to which climate change is influencing those carrying capacities.
Notably, the model found characteristic bifurcations when temperature anomalies increased in both intensity and frequency. These bifurcations indicate that there is some uncertainty in reef responses to thermal-stress events, which may lead to either reef recovery or collapse. That collapse may be, in turn, dependent on specific thresholds, beyond which recovery is unlikely. The thresholds and bifurcations may also offer an estimate of just how little atmospheric CO2 it takes before ocean temperatures drastically increase a coral population’s probability of collapsing.
Thermal anomalies are becoming increasingly common. Especially prominent were the recent back-to-back bleaching events on the northern and central Great Barrier Reef in 2016 and 20175. Estimates of recovery rates on the Great Barrier Reef are still preliminary and mostly remain unknown, although new evidence suggests that recovery rates on the Great Barrier Reef are becoming suppressed with chronic disturbances44. Our study provides a guide as to what recovery rates may be expected under future emissions scenarios. A rate of less than r = 0.38, or its discrete-equation equivalent of lambda 1.46, may indicate that the system is deteriorating and that the corals have already passed a critical threshold and lost their capacity to recover20. Understanding rates of recovery has become the cornerstone of resilience studies21,22. On contemporary reefs, increases in ocean temperature, and potential changes in inter-annual temperature cycles, are spatially variable across ocean regions45,46. Our model shows that coral populations are most sensitive to both the intensity and frequency of thermal-stress events, rather than to incremental ocean warming. However, ocean warming is likely to increase both the intensity and frequency of thermal-stress events. These thermal anomalies also may become more frequent outside of El Niño events5, considerably shortening recovery periods. These temperature thresholds, and the recovery rates of coral populations, are solely dependent on the rates of future greenhouse-gas emissions. Therefore, to ensure the best coral - recovery rates, greenhouse-gas emission rates need to track RCP 6.0 or lower.
References
Walther, G. R. et al. Ecological responses to recent climate change. Nature 416, 389–395 (2002).
Parmesan, C. Ecological and evolutionary responses to recent climate change. Annual Review of Ecology, Evolution, and Systematics 37, 637–669 (2006).
Hughes, T. P. et al. Climate change, human impacts, and the resilience of coral reefs. Science 301, 929–933 (2003).
Hoegh-Guldberg, O. et al. Coral reefs under rapidclimate change and ocean acidification. Science 318, 1737–1742 (2007).
Hughes, T. P. et al. Spatial and temporal patterns of mass bleaching of corals in the Anthropocene. Science 359, 80–83 (2018).
Glynn, P. W. Coral reef bleaching in the 1980s and possible connections with global warming. Trends in Ecology and Evolution 6, 175–1797 (1991).
Brown, B. E. Coral bleaching: causes and consequences. Coral Reefs 16, S129–S138 (1997).
Warner, M. E., Fitt, W. K. & Schmidt, G. W. Damage to photosystem II in symbiotic dinoflagellates: a determinant of coral bleaching. Proceedings of the National Academy of Sciences 96, 8007–8012 (1999).
Takahashi, S., Makahura, T., Sakamizu, M., van Woesik, R. & Yamasaki, H. Repair machinery of symbiotic photosynthesis as the primary target of heat stress for reef-building corals. Plant and Cell Physiology 45, 251–255 (2004).
Baker, A. C., Glynn, P. W. & Riegl, B. Climate change and coral reef bleaching: An ecological assessment of the long-term impacts, recovery trends and future outlook. Estuarine, Coastal and Shelf Science 80, 435–471 (2008).
Loya, Y. et al. Coral bleaching: the winners and losers. Ecology Letters 4, 122–131 (2001).
Baird, A. H. & Marshall, P. A. Mortality, growth and reproduction in scleractinian corals following bleaching on the Great Barrier Reef. Marine Ecology Progress Series 237, 133–141 (2002).
McClanahan, T. Coral bleaching, disease and mortality in the western Indian Ocean. Pages 157–176 in Rosenberg, R. & Loya, Y. editors. Coral Health and Disease, Springer, Berlin, Heidelberg, Germany (2004).
van Woesik, R., Sakai, K., Ganase, A. & Loya, Y. Revisiting the winners and the losers a decade after coral bleaching. Marine Ecology Progress Series 434, 67–76 (2011).
Gilmour, J. P., Smith, L. D., Heyward, A. J., Baird, A. H. & Pratchett, M. S. Recovery of an isolated coral reef system following severe disturbance. Science 340, 69–71 (2013).
Holling, C. S. Resilience and stability of ecological systems. Annual Review of Ecology and Systematics 4, 1–23 (1973).
Folke, C. et al. Regime shifts, resilience, and biodiversity in ecosystem management. Annual Review Of Ecology, Evolution, and Systematics 35, 557–581 (2004).
Scheffer, M. et al. Early-warning signals for critical transitions. Nature 461, 53–59 (2009).
Veraart, A. J. et al. Recovery rates reflect distance to a tipping point in a living system. Nature 481, 357–359 (2012).
Scheffer, M. et al. Anticipating critical transitions. Science 338, 344–348 (2012).
Scheffer, M. & Carpenter, S. R. Catastrophic regime shifts in ecosystems: linking theory to observation. Trends in Ecology and Evolution 18, 648–656 (2003).
Mumby, P. J., Hastings, A. & Edwards, H. J. Thresholds and the resilience of Caribbean coral reefs. Nature 450, 98–101 (2007).
Mumby, P. & van Woesik, R. Consequences of ecological, evolutionary and biogeochemical uncertainty for coral reef responses to climatic stress. Current Biology 24, 413–423 (2014).
Mumby, P. J. Embracing a world of subtly and nuance on coral reefs. Coral Reefs 36, 1003–1011 (2017).
Osborne, K. et al. Delayed coral recovery in a warming ocean. Global Change Biology 23, 3869–3881 (2017).
Crabbe, M. J. C. Global warming and coral reefs: modelling the effects of temperature on Acropora palmata colony growth. Computational Biology and Chemistry 31, 294–297 (2007).
Roth, L., Koksal, S. & van Woesik, R. Effects of thermal stress on key processes driving coral-population dynamics. Marine Ecology Progress Series 411, 73–87 (2010).
Li, X., Wang, H., Zhang, Z. & Hastings, A. Mathematical analysis of coral reef models. J. Math. Anal. Appl. 416, 352–373 (2014).
Lefkovitch, L. P. The study of population growth in organisms grouped by stages. Biometrics 21, 1–18 (1965).
Hughes, T. P. Population dynamics based on individual size rather than age: a general model with a reef coral example. American Naturalist 123, 778–795 (1984).
Done, T. J. Transition matrix models, crown-of-thorns and corals. Pages 291–297 in R. H. Bradbury, editor. Acanthaster and the Coral Reef: A Theoretical Perspective: Proceedings of a Workshop held at the Australian Institute of Marine Science, Townsville, Aug. 6–7, 1998. Vol 88. Springer, Berlin, Heidelberg, Germany (1990).
van Woesik, R. & Koksal, S. A coral population response (CPR) model for thermal stress. Coral Reefs and Climate Change: Science and Management 61, 129–144 (2006).
Baskett, M. L., Gaines, S. D. & Nisbet, R. M. Symbiont diversity may help coral reefs survive moderateclimate change. Ecological Applications 19, 3–17 (2009).
van Woesik, R. Quantifying uncertainty and resilience on coral reefs using a Bayesian approach. Environmental Research Letters 8, 044051 (2013).
Delignette-Muller, M. L. & Dutang, C. fitdistrplus: An R Package for Fitting Distributions. Journal of Statistical Software 64, 1–34 (2015).
Minnett, P. J. et al. Sea-surface temperature measurements from the Moderate-Resolution Imaging Spectroradiometer (MODIS) on Aqua and Terra. In Geoscience and Remote Sensing Symposium, 2004. IGARSS'04. Proceedings. 2004 IEEE International Vol. 7, 4576–4579 (2004).
Trenberth, K. E. The definition of El Nino. Bull. American Meteorological Society 78, 2771–2777 (1997).
IPCC (Intergovernmental Panel on Climate Change). Climate Change 2013: The Physical Science Basis. Contribution of Working Group I to the FifthAssessment Report of the Intergovernmental Panel on Climate Change [Stocker, T. F. et al. (eds)]. Cambridge University Press, Cambridge, United Kingdom and New York, NY, USA, 1535, https://doi.org/10.1017/CBO9781107415324 (2013).
Stocker, T. F. et al. IPCC 2013: climate change 2013: the physical science basis. Contribution of working group I to the fifth assessment report of the intergovernmental panel on climate change. Cambridge, United Kingdom, and New York, USA, pp 203 (2013).
Done, T. J. et al. Coral growth on three reefs: development of recovery benchmarks using a space for time approach. Coral Reefs 29, 815–833 (2010).
Brown, B. E., Le Tissier, M. D. A., Scoffin, T. P. & Tudhop, A. W. Evaluation of the environmental impact of dredging on intertidal coral reefs at Ko Phuket, Thailand, using ecological and physiological parameters. Marine Ecology Progress Series 65, 273–281 (1990).
Baums, I. B., Miller, M. W. & Hellberg, M. E. Regionally isolated populations of an imperiled Caribbean coral, Acropora palmata. Molecular Ecology 14, 1377–1390 (2005).
Gouezo, M. et al. Fifteen years of coral reef monitoring demonstrates the resilience of palau’s coral reefs. Palau International Coral Reef Research Center, Technical Report 17-01, 20 (2017).
Ortiz, J.-C. et al. Impaired recovery of the Great Barrier Reef under cumulative stress. Science Advances 4(7), eaar6127 (2018).
Thompson, D. M. & van Woesik, R. Corals escape bleaching in regions that recently and historically experienced frequent thermal stress. Proceedings of the Royal Society B: Biological Sciences 276, 2893–2901 (2009).
Randall, C. J. & van Woesik, R. Contemporary white-band disease in Caribbean corals driven by climate change. Nature Climate Change 5, 375–379 (2015).
Acknowledgements
We would like to thank Y. Loya and K. Sakai for sharing long-term datasets and J. Gilmour and colleagues for providing all data in their supplementary document. Thanks also to S. J. van Woesik for editorial comments. A. Unal was supported by The Scientific and Technological Research Council of Turkey (TUBITAK), which provided academic leave to A. Unal to visit the Florida Institute of Technology. CWC and RvW would like to thank the National Science Foundation, award NSF OCE-1657633 and OCE-1829393 to Robert van Woesik. CJR was supported by a National Science Foundation (NSF) OCE grant 1535007 to Richard B. Aronson. This is Contribution No. 135 from the Institute for Global Ecology at the Florida Institute of Technology.
Author information
Authors and Affiliations
Contributions
R.v.W., S.K., C.J.R., and A.U. conceptualized the model, A.U. complied the climate change information, C.J.R. compiled the El nino information, C.W.C. wrote the R code, R.v.W. wrote the Open BUGS code, R.v.W. wrote the first draft, and R.v.W., C.J.R., C.W.C., S.K. and A.U. revised and edited the manuscript.
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.
Electronic supplementary material
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 license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license 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 license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Woesik, R.v., Köksal, S., Ünal, A. et al. Predicting coral dynamics through climate change. Sci Rep 8, 17997 (2018). https://doi.org/10.1038/s41598-018-36169-7
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41598-018-36169-7
This article is cited by
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.