Skip to main content

Marine protected areas do not buffer corals from bleaching under global warming

A Correction to this article was published on 21 June 2022

This article has been updated



The rising temperature of the oceans has been identified as the primary driver of mass coral reef declines via coral bleaching (expulsion of photosynthetic endosymbionts). Marine protected areas (MPAs) have been implemented throughout the oceans with the aim of mitigating the impact of local stressors, enhancing fish biomass, and sustaining biodiversity overall. In coral reef regions specifically, protection from local stressors and the enhanced ecosystem function contributed by MPAs are expected to increase coral resistance to global-scale stressors such as marine heatwaves. However, MPAs still suffer from limitations in design, or fail to be adequately enforced, potentially reducing their intended efficacy. Here, we address the hypothesis that the local-scale benefits resulting from MPAs moderate coral bleaching under global warming related stress.


Bayesian analyses reveal that bleaching is expected to occur in both larger and older MPAs when corals are under thermal stress from marine heatwaves (quantified as Degree Heating Weeks, DHW), but this is partially moderated in comparison to the effects of DHW alone. Further analyses failed to identify differences in bleaching prevalence in MPAs relative to non-MPAs for coral reefs experiencing different levels of thermal stress. Finally, no difference in temperatures where bleaching occurs between MPA and non-MPA sites was found.


Our findings suggest that bleaching is likely to occur under global warming regardless of protected status. Thus, while protected areas have key roles for maintaining ecosystem function and local livelihoods, combatting the source of global warming remains the best way to prevent the decline of coral reefs via coral bleaching.

Peer Review reports


Rising ocean temperatures and increased frequency and duration of marine heatwaves [1] are causing the decline of coral reefs at alarming rates via coral bleaching [2, 3]—the process whereby photosynthetic endosymbionts are expelled, revealing the coral skeleton [4,5,6,7,8]. Sustained ocean heat stress can lead to mass bleaching events which may result in mortality of entire coral colonies [3, 9, 10]. If lethal bleaching occurs, the loss of coral cover results in habitat homogenisation and consequently reduced biodiversity [10, 11]. Ultimately, such reductions in coral reef biodiversity inhibits the ecosystem function of coral reefs, critical for supporting > 25% of marine species [12] and for providing ecosystem services to over 100 million people circumtropically [13].

While global warming is unequivocally the predominant driver of mass coral bleaching, a myriad of local scale factors can also induce bleaching of corals. Factors such as turbidity [14, 15], eutrophication [16, 17], hypoxia [16, 18], and sedimentation [14] have been documented to independently induce coral bleaching. However, pioneering studies identified reduced mortality from coral bleaching under higher levels of sedimentation when exposed to heat stress, likely as a result of reduced solar irradiance [14]. Yet, bleaching is known to occur under high temperature regimes regardless of irradiance [19]. Despite convoluted evidence of these interactions influencing bleaching, a widespread expectation that these local stressors interact either additively or synergistically with global warming to exacerbate coral bleaching exists, with field evidence from Mesoamerican reefs [20, 21], and French Polynesia [17].

To mitigate the additive effects of local scale stressors on marine biodiversity overall, Marine Protected Areas (MPAs) have been implemented across different regions of the world, which when effectively designed, are established to enhance regional biodiversity and general ecosystem health [22]. Large and long-established MPAs are often especially effective for enhancing multiple metrics used for monitoring ecosystem health [23, 24], such as coral cover [25], fish biomass [26] and biodiversity [22]. For corals specifically, the role MPAs perform for reducing local stressors intend to enhance coral health through a variety of physiological mechanisms [27], thereby promoting resistance of reef building corals to disturbance. Furthermore, MPAs have the potential to promote resilience to disturbance events, such as marine heatwaves, disease outbreaks, and hurricanes, via ecological processes [28]. This enhanced resilience intrinsically promotes resistance to future bleaching by facilitating full recovery from bleaching before the next disturbance event [11]. Given that marine heatwaves and bleaching events are increasing in frequency and intensity through time [1, 2, 11, 29], the benefits of MPAs for promoting resilience in reef building corals are subsequently crucial for also enhancing the resistance of corals to future bleaching—i.e. managed resilience [22].

However, the effects of MPAs for mitigating coral reef decline remain contested. For example, decline in coral cover attributed to thermal stress is not mitigated by MPAs [30], suggesting that the preservation of coral reefs does not depend significantly on MPAs, but on actions that mitigate the degree of climate warming [22, 31]. Furthermore, multiple stressors on coral reefs tend to be antagonistic rather than synergistic, especially interactions between local stressors, and global warming [22, 32,33,34,35]. This is likely owing to co-sensitivity and co-tolerance of coral species exposed to stressors, along with the frightening prospect that climate warming eclipses the potential advantages that could be expected to result from the mitigation of local stressors [31].

Given both the convoluted relationship between global and local stressors exerted on coral reefs, along with the diversity of primary objectives different MPA’s aim to achieve, it is crucial to discern whether MPA’s have any moderation effect on bleaching under global warming. However, an explicit global scale test to examine the prevalence of bleaching in relation to their protected status remains lacking. To test this hypothesis for coral bleaching specifically, we examine (1) the probability of coral bleaching under thermal stress (quantified as Degree Heating Weeks, DHW) for key MPA attributes—the size and age of MPAs, using a Bayesian Generalised Linear Mixed Model; (2) implement quantitative comparisons of bleaching prevalence on coral reefs within and outside protected regions under different levels of thermal stress; and (3) compare thermal thresholds where the onset of bleaching occurs between protected and non-protected coral reefs based on the gamma distributions of DHWs. To address this hypothesis we utilise a global scale data set containing 8,766 coral bleaching surveys (Fig. 1) over a 16 year period.

Fig. 1
figure 1

The code employed to produce figure 1, and to access CoRTAD and MODIS data used in this study is available from Robert van Woesik and Shannon Sully (

The richness and global scale distribution of Reef Check surveys used to examine the effects of MPAs on coral bleaching. a Represents 5393 surveys which do not fall within an MPA. b Are the 3391 surveys which do fall within the jurisdiction of an MPA.


Bleaching within MPAs

Our Bayesian model from 3391 reef check surveys reveals that older MPAs reduce the probability of coral bleaching, while the size of a MPA shows little evidence for predicting bleaching (Fig. 2).

However, the interaction between DHW and MPA age shows evidence for predicting an increased probability in bleaching based on the 80% credible interval (CI), which becomes weaker when interpreting the 95% CI. Meanwhile, the interaction between DHW and MPA size show little evidence for predicting bleaching as both the 80% and 95% credible intervals cross zero (Fig. 2). These interactions suggest bleaching is likely to occur when under heat stress, but do show a partial moderation in comparison to the sole predictor of DHW.

Fig. 2
figure 2

Model coefficients of predictions from the Bayesian Generalized Mixed Effect Model. DHW are degree heating weeks, MPA attributes are size and age. Blue dot represents mean coefficient, thick dark blue bars are 80% credible intervals (CI), while the whiskers (thin grey bars) are 95% CIs. Data are from 3391 coral bleaching surveys which fell within the jurisdiction of an MPA as categorized by the IUCNs World Database of Protected Areas

Bleaching comparisons between MPA and non-MPA coral reefs

When comparing the average bleaching prevalence for reef check surveys across different DHW categories (none 0, low < 1.5, medium 1.5–4, and high > 4) with a likelihood ratio test there were no statistically significant differences between MPA and non-MPA sites (Fig. 3; Table 1), suggesting that the sites protected status had negligible effect on bleaching prevalence.

Fig. 3
figure 3

Density plots of average bleaching for each degree heating week (DHW) category; none (zero DHW), low (0–1.5 DHW), medium (1.5–4 DHW) and high (> 4 DHW). Average bleaching were log transformed for visual display purposes only. The top panel represents the distribution density of bleaching between MPA and non-MPA sites at high, medium, low, and none DHW categories. The bottom panel show the stacked proportional densities for MPA (left) and non-MPA sites (right) for each DHW category

Table 1 Summarised average bleaching prevalence for coral reef sites and their likelihood of being statistically different between MPA and non-MPA sites based on their thermal stress (DHW category)

Finally, when comparing the DHW values where the onset of coral bleaching has occurred (i.e. the mean DHW value associated with coral bleaching) with a likelihood ratio test, there is no significant difference between MPA and non-MPA sites (Fig. 4) (Likelihood ratio test, Pr(>χ2) = 0.998). Thus, protected status appears to have little, if any, moderation effect on coral bleaching prevalence under thermal stress from climate change when comparing protected reefs to non-protected reefs.

Fig. 4
figure 4

Probability densities for the degree heating week (DHW) values where the onset of coral bleaching occurred based on gamma distributions between MPA and non-MPA Reef Check survey sites. A constant of 1 was added to DHW values to meet gamma distribution assumptions where values are required to be greater than zero [36]


Our analyses reveal that MPAs play a negligible role in mitigating the onset of coral bleaching when under climate change related thermal stress. However, weak evidence for moderation under the interaction of DHW with MPA size, and MPA age, in comparison to the sole predictor of DHW does exist (Fig. 2). Crucially, however, bleaching is still predicted to occur under these interactions. Given our findings, with no discernible difference in temperatures where the onset of bleaching occurs, we expect that coral resistance to thermally induced bleaching is unlikely enhanced by the implementation of protected status. These findings are in accordance with previous studies investigating loss of coral cover within MPAs under climate change [30], and further challenge the assumed benefits of managed resilience for promoting resistance in corals to guide reefs through the gauntlet of climate change [22, 31].

Given the similar bleaching responses between coral reefs residing within and outside MPAs, our findings add to the growing complexity of the relationship between local and global scale stressors for degrading coral reefs. Comparable levels of bleaching where no temperature stress is present (0 DHW category) indicates no difference in bleaching under ambient conditions, which is likely a result of localised conditions [29], survey error [37], and perhaps lack of recovery from a previous disturbance event - i.e. non-branching corals which are able to survive longer while bleached [38]. Furthermore, an identical thermal threshold where the onset of bleaching occurs suggests there is a not a synergistic relationship between local and global factors which exert stress onto coral reefs. Non-synergistic and antagonistic relationships have been widely reported on coral reefs over the last 10 years [32,33,34] challenging the previous supposition that stressors exerted onto coral reefs act synergistically [39]. Our findings could indicate that bleaching is unlikely to be synergistically exacerbated by local stressors [35], which are assumed to be moderated by MPAs, given the identical bleaching responses between MPA and non-MPA environments (Fig. 4). It is likely that the effects of climate change are far eclipsing the role of localised stressors and thus localised mitigation [11, 22, 40]. It should be noted, however, that many MPAs are not adequately managed in many marine regions [41], and have aims focused on social-economic and biodiversity benefits [42], thus our findings may also reflect this. Furthermore, owing to the spatial variability of coral bleaching [29], which is often specific to a wide range of factors such as turbidity [43, 44], internal waves [45, 46], evolutionary history [47], and ecological memory [48], exceptions to the global scale pattern will exist [20, 21].

Our findings also identify insufficient evidence to support the managed resilience hypothesis for reef corals, because bleaching responses are similar between MPA and non-MPA sites (Figs. 3, 4). Consequently, our findings suggests that the resistance of reef building corals will not be enhanced through the implementation of MPAs, which aim to mitigate local stressors and ameliorate physiological performance of corals [27]. However, it is critical to note there are many other benefits of MPAs which ensure food provision, vital for human livelihood, and maintain biodiversity which is critical for ecosystem function around the globe [23, 24, 42]. Yet, the assumption that MPAs will help support coral reefs by preventing impacts of warming (i.e. bleaching) on corals is likely incorrect based on these findings. Rather, continued warming linked to anthropogenic activity will incessantly bleach corals more often through the Anthropocene [2] regardless of protected status.


Collectively, our findings add to the growing evidence that protected status will have little impact for alleviating the effects global stressors such as marine heatwaves which will continue to be exerted on coral reefs. While the implementation of effectively designed MPAs can be beneficial for coral cover and maintaining functional species [22, 25], and most critically support communities dependent on coral reefs [49], they will not mitigate the effects of coral bleaching induced by global warming [22, 30]. Consequently, actions targeting the source of rising global temperatures (i.e. greenhouse gas emissions) remains the most effective way to moderate future coral bleaching caused by global warming, and thus mitigating continued global coral reef decline [11, 22, 39, 50].

Materials and methods

Bleaching data

Bleaching data were collated from, combining 8798 surveys from 3,067 sites across 73 countries (Fig. 1), from the years 2002–2018. Validity of the Reef Check data are well established with less than 7% sampling error for identifying components of benthic cover [37]. These data have also been used in previous macroecological coral bleaching studies [29, 44, 51]. From Reef Check data, the percentage of bleached coral populations within each reef check survey (i.e. % reef bleached) were extracted (i.e. site wide bleaching), along with the date of the survey and their geographic coordinates. This allowed for each survey to be spatially and temporally associated with environmental data. The global scale of these data, along with the long term time frame and high sample size provide a robust basis for elucidating relationships between environmental drivers associated with bleaching.

MPA data

Each Reef Check survey was temporally and spatially associated as being inside or outside an MPA as designated by the IUCNs world data base of protected areas [52] at the time of survey. The reef check survey data were spatially associated by overlaying the MPA shape file with the reef check coordinates using the ‘sp’ package [53] in R studio 4.0 [54]. From the MPA shape file, all data attributes were extracted for sites which fell within the MPA polygons, summarised in Additional file 1: Table S1. A survey was considered to be in a protected area if it fell within a designated protected region as specified by the IUCN protected area categories [55], shown in Table 2. These protected status categories represent the global definition of a protected area, however all will vary with their level of enforcement, and type of protection incurred. For example, this study has 877 surveys which fell within a protected region that is also a no-take zone, while 1876 surveys within MPAs were not reported (i.e., are unknown) to have a no-take zone. Therefore, for simplicity, we are using the simple definition of whether a survey fell within a protected area to determine how any level of protected status interacts with global warming for influencing coral bleaching.

Table 2 IUCN Protected Area Management Category as defined in the IUCN Manual of World Database on Protected Areas User Manual [55]

Environmental data

To determine whether bleaching differs between MPAs and non-MPAs under thermal stress, we used the degree heating weeks (DHWs) metric to quantify thermal anomalies. DHWs represent the global standard for predicting bleaching from thermal stress and are highly robust in their predictions [56]. These data were extracted on a weekly time series from the Coral Reef Temperature Anomaly database (CoRTAD version 6), supplied from the National Oceanic and Atmospheric Administration [57] at a resolution of ~ 4.6 km at the equator. The DHW values were temporally and spatially associated with each coral bleaching survey on a weekly time series using the ‘ncdf4’ [58] package in R. DHW are the global standard for determining bleaching likelihood under thermal stress, where one DHW represents 1 °C increase in the local mean climatic temperature for one week over the last 12 weeks.

We included turbidity within our model to account for other environmental drivers which influence coral physiology, and thus bleaching responses [15, 44]. The diffuse attenuation coefficient of light at the 490 nm wavelength is positively related to turbidity, and has been ubiquitously utilised for deriving turbidity measurements in coral reef studies [44, 59, 60]. The kd490 values were extracted from the Modis-Aqua satellite database ( maintained by NASA’s Earth Observation System Data and Information System (EOSDIS). These data were spatially and temporally matched up with each coral bleaching survey on a weekly time series, at a 4 km resolution.

Statistical analyses

We developed a collinearity matrix using a Pearson’s correlation coefficient on the log transformed MPA attributes (Additional file 1: Fig. S1). We excluded highly collinear data using a conservative 0.65 cut off to prevent collinearity inhibiting convergence of the Bayesian model. If MPA attributes for a survey were missing, the survey was also excluded from the analysis. In total, 3391 Reef Check surveys fell within MPAs at the time of survey which were used in the Bayesian model.

We used a Generalised Linear Model with group specific terms in ‘rstanarm’ [61], which uses the STAN language [62] in R 4.0 [54]. The response variable, average bleaching prevalence, was modelled with a Negative Binomial distribution as data were dominated by zeros, and is a distribution which has been successfully used on these bleaching data in the past [29]. Coral ecoregions were run as a random effect to account for spatial variation in bleaching patterns (Additional file 1: Fig S4). Coral ecoregions were extracted from Veron et al. [63], Corals Ecoregions of The World (COTW). Any survey which fell outside the ecoregion polygons was excluded from analysis. These coral ecoregions represent consistent patterns of taxonomic configuration, dispersal and isolation processes, and patterns of evolutionary history [63, 64], utilised in multiple global scale coral reef studies to account for spatial variation [29, 44, 59]. The default weakly informative priors which internally adjust in scale based on regression coefficients were used in the model [61]. Predictors were standardised to aid convergence. Posterior predictive checks from the ‘Bayesplot’ package [65] in R were used to assess for model fit (Additional file 1: Fig S2). Convergence was visually assessed through trace plots (Additional file 1: Fig. S3) and achieved when the rhat value (Gelman-Rubin statistic) reached 1.0 [61]. The model ran with 4 chains, for 5,000 iterations with 2500 warmups.

Differences in bleaching prevalence between MPA and non-MPA sites were statistically analysed by transformation of bleaching prevalence data into negative binomial distributions using the ‘fitdist’ function in the ‘MASS’ package [36]. The average bleaching prevalence data were split into 4 groupings of thermal stress thresholds of none (0), low (0–1.49), medium (1.5–3.9), and high (> 4) DHW. Bleaching prevalence distributions for each DHW category were then analysed with a likelihood ratio test with 2 degrees of freedom and a true lower tail [29, 44]. Finally, differences in the DHW values associated with bleaching at MPA and non-MPA sites were also statistically analysed with a Likelihood ratio test. The DHW values were fit into a gamma distribution from the ‘MASS’ package. The DHW values > 12 were removed from analysis here, with a constant of 1 also added so data could fit the gamma distribution [36]. The degrees of freedom of the likelihood ratio test were again 2, with the lower tail specified as true.

Availability of data and materials

All data and code used in this study are available on our Github repository ( All data used in this study are also publically available from cited sources in the methods.

Change history


  1. Oliver ECJ, Donat MG, Burrows MT, Moore PJ, Smale DA, Alexander LV, et al. Longer and more frequent marine heatwaves over the past century. Nat Commun. 2018;9:1324.

    PubMed  PubMed Central  Google Scholar 

  2. Hughes TP, Anderson KD, Connolly SR, Heron SF, Kerry JT, Lough JM, et al. Spatial and temporal patterns of mass bleaching of corals in the Anthropocene. Science. 2018;359:80–3.

    CAS  PubMed  Google Scholar 

  3. Hughes TP, Kerry JT, Baird AH, Connolly SR, Dietzel A, Eakin CM, et al. Global warming transforms coral reef assemblages. Nature. 2018;556:492–6.

    CAS  PubMed  Google Scholar 

  4. Warner ME, Fitt WK, Schmidt GW. Damage to photosystem II in symbiotic dinoflagellates: a determinant of coral bleaching. PNAS. 1999;96:8007–12.

    CAS  PubMed  PubMed Central  Google Scholar 

  5. Fitt WK, Brown BE, Warner ME, Dunne RP. Coral bleaching: interpretation of thermal tolerance limits and thermal thresholds in tropical corals. Coral Reefs. 2001;20:51–65.

    Google Scholar 

  6. Douglas AE. Coral bleaching––how and why? Mar Pollut Bull. 2003;46:385–92.

    CAS  PubMed  Google Scholar 

  7. Suggett DJ, Smith DJ. Interpreting the sign of coral bleaching as friend vs. foe. Glob Chang Biol. 2011;17:45–55.

    Google Scholar 

  8. Boilard A, Dubé CE, Gruet C, Mercière A, Hernandez-Agreda A, Derome N. Defining coral bleaching as a microbial dysbiosis within the coral holobiont. Microorganisms. 2020;8:1682.

    CAS  PubMed Central  Google Scholar 

  9. Hughes TP, Kerry JT, Baird AH, Connolly SR, Chase TJ, Dietzel A, et al. Global warming impairs stock–recruitment dynamics of corals. Nature. 2019;568:387–90.

    CAS  PubMed  Google Scholar 

  10. Magel JMT, Burns JHR, Gates RD, Baum JK. Effects of bleaching-associated mass coral mortality on reef structural complexity across a gradient of local disturbance. Sci Rep. 2019;9:2512.

    PubMed  PubMed Central  Google Scholar 

  11. Hughes TP, Barnes ML, Bellwood DR, Cinner JE, Cumming GS, Jackson JBC, et al. Coral reefs in the Anthropocene. Nature. 2017;546:82–90.

    CAS  PubMed  Google Scholar 

  12. Fisher R, O’Leary RA, Low-Choy S, Mengersen K, Knowlton N, Brainard RE, et al. Species richness on coral reefs and the pursuit of convergent global estimates. Curr Biol. 2015;25:500–5.

    CAS  PubMed  Google Scholar 

  13. Moberg F, Folke C. Ecological goods and services of coral reef ecosystems. Ecol Econ. 1999;29:215–33.

    Google Scholar 

  14. Anthony KRN, Connolly SR, Hoegh-Guldberg O. Bleaching, energetics, and coral mortality risk: effects of temperature, light, and sediment regime. Limnol Oceanogr. 2007;52:716–26.

    Google Scholar 

  15. Bessell-Browne P, Negri AP, Fisher R, Clode PL, Duckworth A, Jones R. Impacts of turbidity on corals: the relative importance of light limitation and suspended sediments. Mar Pollut Bull. 2017;117:161–70.

    CAS  PubMed  Google Scholar 

  16. Zhu B, Wang G, Huang B, Tseng CK. Effects of temperature, hypoxia, ammonia and nitrate on the bleaching among three coral species. Chin Sci Bull. 2004;49:1923–8.

    CAS  Google Scholar 

  17. Donovan MK, Adam TC, Shantz AA, Speare KE, Munsterman KS, Rice MM, et al. Nitrogen pollution interacts with heat stress to increase coral bleaching across the seascape. PNAS. 2020;117:5351–7.

    CAS  PubMed  PubMed Central  Google Scholar 

  18. Altieri AH, Harrison SB, Seemann J, Collin R, Diaz RJ, Knowlton N. Tropical dead zones and mass mortalities on coral reefs. PNAS. 2017;114:3660–5.

    CAS  PubMed  PubMed Central  Google Scholar 

  19. Tolleter D, Seneca FO, DeNofrio JC, Krediet CJ, Palumbi SR, Pringle JR, et al. Coral bleaching independent of photosynthetic activity. Curr Biol. 2013;23:1782–6.

    CAS  PubMed  Google Scholar 

  20. Carilli JE, Norris RD, Black BA, Walsh SM, McField M. Local stressors reduce coral resilience to bleaching. PLoS ONE. 2009;4:e6324.

    PubMed  PubMed Central  Google Scholar 

  21. Carilli JE, Norris RD, Black B, Walsh SM, McField M. Century-scale records of coral growth rates indicate that local stressors reduce coral thermal tolerance threshold. Glob Chang Biol. 2010;16:1247–57.

    Google Scholar 

  22. Bruno JF, Côté IM, Toth LT, Climate Change. Coral loss, and the curious case of the parrotfish paradigm: why don’t marine protected areas improve reef resilience? Ann Rev Mar Sci. 2019;11:307–34.

    PubMed  Google Scholar 

  23. Edgar GJ, Stuart-Smith RD, Willis TJ, Kininmonth S, Baker SC, Banks S, et al. Global conservation outcomes depend on marine protected areas with five key features. Nature. 2014;506:216–20.

    CAS  PubMed  Google Scholar 

  24. Strain EMA, Edgar GJ, Ceccarelli D, Stuart-Smith RD, Hosack GR, Thomson RJ. A global assessment of the direct and indirect benefits of marine protected areas for coral reef conservation. Divers Distrib. 2019;25:9–20.

    Google Scholar 

  25. Selig ER, Bruno JF. A global analysis of the effectiveness of marine protected areas in preventing coral loss. PLoS ONE. 2010;5:e9278.

    PubMed  PubMed Central  Google Scholar 

  26. Graham NAJ, Robinson JPW, Smith SE, Govinden R, Gendron G, Wilson SK. Changing role of coral reef marine reserves in a warming climate. Nat Commun. 2020;11:2000.

    CAS  PubMed  PubMed Central  Google Scholar 

  27. Lamb JB, Wenger AS, Devlin MJ, Ceccarelli DM, Williamson DH, Willis BL. Reserves as tools for alleviating impacts of marine disease. Philos Trans R Soc B Biol Sci. 2016;371:20150210.

    Google Scholar 

  28. Mellin C, MacNeil MA, Cheal AJ, Emslie MJ, Caley MJ. Marine protected areas increase resilience among coral reef communities. Ecol Lett. 2016;19:629–37.

    PubMed  Google Scholar 

  29. Sully S, Burkepile DE, Donovan MK, Hodgson G, van Woesik R. A global analysis of coral bleaching over the past two decades. Nat Commun. 2019;10:1264.

    CAS  PubMed  PubMed Central  Google Scholar 

  30. Selig ER, Casey KS, Bruno JF. Temperature-driven coral decline: the role of marine protected areas. Glob Chang Biol. 2012;18:1561–70.

    Google Scholar 

  31. Bruno JF, Bates AE, Cacciapaglia C, Pike EP, Amstrup SC, van Hooidonk R, et al. Climate change threatens the world’s marine protected areas. Nat Clim Chang. 2018;8:499–503.

    Google Scholar 

  32. Darling ES, McClanahan TR, Côté IM. Combined effects of two stressors on Kenyan coral reefs are additive or antagonistic, not synergistic. Conserv Lett. 2010;3:122–30.

    Google Scholar 

  33. Bruno JF, Valdivia A. Coral reef degradation is not correlated with local human population density. Sci Rep. 2016;6:29778.

    CAS  PubMed  PubMed Central  Google Scholar 

  34. Côté IM, Darling ES, Brown CJ. Interactions among ecosystem stressors and their importance in conservation. Proc R Soc B Biol Sci. 2016;283:20152592.

    Google Scholar 

  35. Johnson JV, Dick JTA, Pincheira-Donoso D. Local anthropogenic stress does not exacerbate coral bleaching under global climate change. Glob Ecol Biogeogr. 

  36. Ripley B, Venables B, Bates DM, (ca 1998) KH (partial port, ca 1998) AG (partial port, Firth D. MASS: support functions and datasets for venables and Ripley’s MASS. 2022.

  37. Done T, Roelfsema C, Harvey A, Schuller L, Hill J, Schläppy M-L, et al. Reliability and utility of citizen science reef monitoring data collected by Reef Check Australia, 2002–2015. Mar Pollut Bull. 2017;117:148–55.

    CAS  PubMed  Google Scholar 

  38. van Woesik R, Sakai K, Ganase A, Loya Y. Revisiting the winners and the losers a decade after coral bleaching. Mar Ecol Prog Ser. 2011;434:67–76.

    Google Scholar 

  39. Hughes TP, Graham NAJ, Jackson JBC, Mumby PJ, Steneck RS. Rising to the challenge of sustaining coral reef resilience. Trends Ecol Evol. 2010;25:633–42.

    PubMed  Google Scholar 

  40. Weijerman M, Veazey L, Yee S, Vaché K, Delevaux JMS, Donovan MK, et al. Managing local stressors for coral reef condition and ecosystem services delivery under climate scenarios. Front Mar Sci. 2018;5.

  41. Gill DA, Mascia MB, Ahmadia GN, Glew L, Lester SE, Barnes M, et al. Capacity shortfalls hinder the performance of marine protected areas globally. Nature. 2017;543:665–9.

    CAS  PubMed  Google Scholar 

  42. Belgrano A, Novaglio C, Svedäng H, Villasante S, Melián CJ, Blenckner T, et al. Mapping and Evaluating marine protected areas and ecosystem services: a transdisciplinary delphi forecasting process framework. Front Ecol Evol. 2021;9.

  43. Cacciapaglia C, van Woesik R. Climate-change refugia: shading reef corals by turbidity. Glob Chang Biol. 2016;22:1145–54.

    PubMed  Google Scholar 

  44. Sully S, van Woesik R. Turbid reefs moderate coral bleaching under climate-related temperature stress. Glob Chang Biol. 2020;26:1367–73.

    PubMed  PubMed Central  Google Scholar 

  45. Wall M, Putchim L, Schmidt GM, Jantzen C, Khokiattiwong S, Richter C. Large-amplitude internal waves benefit corals during thermal stress. Proc R Soc B Biol Sci. 2015;282:20140650.

    CAS  Google Scholar 

  46. Storlazzi CD, Cheriton OM, van Hooidonk R, Zhao Z, Brainard R. Internal tides can provide thermal refugia that will buffer some coral reefs from future global warming. Sci Rep. 2020;10:13435.

    CAS  PubMed  PubMed Central  Google Scholar 

  47. Osman EO, Smith DJ, Ziegler M, Kürten B, Conrad C, El-Haddad KM, et al. Thermal refugia against coral bleaching throughout the northern Red Sea. Glob Chang Biol. 2018;24:e474-84.

    PubMed  Google Scholar 

  48. Hackerott S, Martell HA, Eirin-Lopez JM. Coral environmental memory: causes, mechanisms, and consequences for future reefs. Trends Ecol Evol. 2021.

    Article  PubMed  Google Scholar 

  49. Mcleod E, Anthony KRN, Mumby PJ, Maynard J, Beeden R, Graham NAJ, et al. The future of resilience-based management in coral reef ecosystems. J Environ Manage. 2019;233:291–301.

    PubMed  Google Scholar 

  50. Bates AE, Cooke RSC, Duncan MI, Edgar GJ, Bruno JF, Benedetti-Cecchi L, et al. Climate resilience in marine protected areas and the ‘Protection Paradox.’ Biol Conserv. 2019;236:305–14.

    Google Scholar 

  51. Johnson J, Dick J, Pincheira-Donoso D. Mangrove forests mitigate coral bleaching under thermal stress from climate change. bioRxiv. 2021;2021.06.04.447049.

  52. UNEP-WCMC and IUCN. The world database on protected areas (WDPA). 2021.

  53. Pebesma E, Bivand R, Rowlingson B, Gomez-Rubio V, Hijmans R, Sumner M, et al. sp: Classes and Methods for Spatial Data. 2021.

  54. R Core Team. R: a language and environment for statistical computing. Vienna, Austria: R Foundation for Statistical Computing; 2021.

    Google Scholar 

  55. UNEP-WCMC. World Database on Protected Areas User Manual 1.5. UNEP-WCMC. Cambridge, UK; 2017.

  56. Kayanne H. Validation of degree heating weeks as a coral bleaching index in the northwestern Pacific. Coral Reefs. 2017;36:63–70.

    Google Scholar 

  57. US Department of Commerce NNC for EI. NODC-UNC CoRTAD. Accessed 21 Oct 2020.

  58. Pierce D, Pierce MD. Package ‘ncdf4.’ 2019.

  59. Donovan MK, Burkepile DE, Kratochwill C, Shlesinger T, Sully S, Oliver TA, et al. Local conditions magnify coral loss after marine heatwaves. Science. 2021;372:977–80.

    CAS  PubMed  Google Scholar 

  60. van Woesik R, Kratochwill C. A global coral-bleaching database, 1980–2020. Sci Data. 2022;9:20.

    PubMed  PubMed Central  Google Scholar 

  61. Goodrich B, Gabry J, Ali I, Brilleman S. rstanarm: Bayesian applied regression modeling via Stan. R package version. 2018;2:1758.

    Google Scholar 

  62. Carpenter B, Gelman A, Hoffman MD, Lee D, Goodrich B, Betancourt M, et al. Stan: a probabilistic programming language. Grantee Submission. 2017;76:1–32.

    Google Scholar 

  63. Veron J, Stafford-Smith M, DeVantier L, Turak E. Overview of distribution patterns of zooxanthellate Scleractinia. Front Mar Sci. 2015;1.

  64. Spalding MD, Fox HE, Allen GR, Davidson N, Ferdaña ZA, Finlayson M, et al. Marine ecoregions of the world: a bioregionalization of coastal and shelf areas. Bioscience. 2007;57:573–83.

    Google Scholar 

  65. Gabry J, Mahr T, Bürkner P-C, Modrák M, Barrett M, Weber F, et al. Bayesplot: plotting for Bayesian Models. 2021.

Download references


JVJ is funded by the Department for Economy (DfE) Northern Ireland. We thank the 1,000s of citizen scientist that have contributed to Reef check over the last 20 years, and the founders of reef check for availability of their global dataset. We also thank NOAA for maintenance of the CoRTAD, and the public accessibility of their data.


JVJ is funded by the Department for Economy, Northern Ireland.

Author information

Authors and Affiliations



JVJ and DPD designed the study. JVJ collated data and analysed results. JVJ, DPD, and JTAD wrote the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Daniel Pincheira-Donoso.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

The original online version of this article was revised: Following the publication of the original article, we were notified that the authors had overlooked (and thus, failed to explicitly acknowledge) the source for code availability which facilitated the production of figure 1 and data collection in the study

The code employed to produce figure 1, and to access CoRTAD and MODIS data used in this study is available from Robert van Woesik and Shannon Sully (

Supplementary Information

Additional file 1. Table S1.

Environmental covariates and MPAattributes used in this study, with a definition and reference to their source. Fig. S1 Collinearity matrixof explanatory variables for coral bleaching from the Coral Reef TemperatureAnomaly Database (CoRTAD V6) and the World Database of Protected Areas (WDPA).Attributes with 0.65 collinearity score or higher were excluded from the Bayesianmodel. Fig. S2 Posterior predictivechecks of Bayesian model fit with our bleaching data. Figure S3. Trace plots after burnins discardedfrom the Bayesian Generalized linear model with group specific terms ranthrough STAN for each covariate analyzed in the model. Figure S4. Intercept variance for each ecoregionran as a random effect within the STAN Bayesian model.

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 The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Johnson, J.V., Dick, J.T.A. & Pincheira-Donoso, D. Marine protected areas do not buffer corals from bleaching under global warming. BMC Ecol Evo 22, 58 (2022).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: