- Open Access
Marine protected areas do not buffer corals from bleaching under global warming
BMC Ecology and Evolution volume 22, Article number: 58 (2022)
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.
Rising ocean temperatures and increased frequency and duration of marine heatwaves  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  and for providing ecosystem services to over 100 million people circumtropically .
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  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 . Yet, bleaching is known to occur under high temperature regimes regardless of irradiance . 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 .
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 . Large and long-established MPAs are often especially effective for enhancing multiple metrics used for monitoring ecosystem health [23, 24], such as coral cover , fish biomass  and biodiversity . For corals specifically, the role MPAs perform for reducing local stressors intend to enhance coral health through a variety of physiological mechanisms , 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 . This enhanced resilience intrinsically promotes resistance to future bleaching by facilitating full recovery from bleaching before the next disturbance event . 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 .
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 , 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 .
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.
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.
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.
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.
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 , 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 , survey error , and perhaps lack of recovery from a previous disturbance event - i.e. non-branching corals which are able to survive longer while bleached . 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 . Our findings could indicate that bleaching is unlikely to be synergistically exacerbated by local stressors , 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 , and have aims focused on social-economic and biodiversity benefits , thus our findings may also reflect this. Furthermore, owing to the spatial variability of coral bleaching , which is often specific to a wide range of factors such as turbidity [43, 44], internal waves [45, 46], evolutionary history , and ecological memory , 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 . 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  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 , 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 were collated from reefcheck.org, 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 . 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.
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  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  in R studio 4.0 . 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 , 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.
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 . 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  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’  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 (https://oceandata.sci.gsfc.nasa.gov/MODIS-Aqua/Mapped/Monthly/4km/Kd_490/) 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.
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’ , which uses the STAN language  in R 4.0 . 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 . 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. , 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 . Predictors were standardised to aid convergence. Posterior predictive checks from the ‘Bayesplot’ package  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 . 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 . 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 . 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 (https://github.com/JackVJohnson/Bleaching-in-MPAs). All data used in this study are also publically available from cited sources in the methods.
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.
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.
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.
Warner ME, Fitt WK, Schmidt GW. Damage to photosystem II in symbiotic dinoflagellates: a determinant of coral bleaching. PNAS. 1999;96:8007–12.
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.
Douglas AE. Coral bleaching––how and why? Mar Pollut Bull. 2003;46:385–92.
Suggett DJ, Smith DJ. Interpreting the sign of coral bleaching as friend vs. foe. Glob Chang Biol. 2011;17:45–55.
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.
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.
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.
Hughes TP, Barnes ML, Bellwood DR, Cinner JE, Cumming GS, Jackson JBC, et al. Coral reefs in the Anthropocene. Nature. 2017;546:82–90.
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.
Moberg F, Folke C. Ecological goods and services of coral reef ecosystems. Ecol Econ. 1999;29:215–33.
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.
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.
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.
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.
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.
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.
Carilli JE, Norris RD, Black BA, Walsh SM, McField M. Local stressors reduce coral resilience to bleaching. PLoS ONE. 2009;4:e6324.
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.
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.
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.
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.
Selig ER, Bruno JF. A global analysis of the effectiveness of marine protected areas in preventing coral loss. PLoS ONE. 2010;5:e9278.
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.
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.
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.
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.
Selig ER, Casey KS, Bruno JF. Temperature-driven coral decline: the role of marine protected areas. Glob Chang Biol. 2012;18:1561–70.
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.
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.
Bruno JF, Valdivia A. Coral reef degradation is not correlated with local human population density. Sci Rep. 2016;6:29778.
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.
Johnson JV, Dick JTA, Pincheira-Donoso D. Local anthropogenic stress does not exacerbate coral bleaching under global climate change. Glob Ecol Biogeogr.
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.
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.
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.
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.
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.
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.
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.
Cacciapaglia C, van Woesik R. Climate-change refugia: shading reef corals by turbidity. Glob Chang Biol. 2016;22:1145–54.
Sully S, van Woesik R. Turbid reefs moderate coral bleaching under climate-related temperature stress. Glob Chang Biol. 2020;26:1367–73.
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.
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.
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.
Hackerott S, Martell HA, Eirin-Lopez JM. Coral environmental memory: causes, mechanisms, and consequences for future reefs. Trends Ecol Evol. 2021. https://doi.org/10.1016/j.tree.2021.06.014.
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.
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.
Johnson J, Dick J, Pincheira-Donoso D. Mangrove forests mitigate coral bleaching under thermal stress from climate change. bioRxiv. 2021;2021.06.04.447049.
UNEP-WCMC and IUCN. The world database on protected areas (WDPA). 2021. www.protectedplanet.net.
Pebesma E, Bivand R, Rowlingson B, Gomez-Rubio V, Hijmans R, Sumner M, et al. sp: Classes and Methods for Spatial Data. 2021.
R Core Team. R: a language and environment for statistical computing. Vienna, Austria: R Foundation for Statistical Computing; 2021.
UNEP-WCMC. World Database on Protected Areas User Manual 1.5. UNEP-WCMC. Cambridge, UK; 2017.
Kayanne H. Validation of degree heating weeks as a coral bleaching index in the northwestern Pacific. Coral Reefs. 2017;36:63–70.
US Department of Commerce NNC for EI. NODC-UNC CoRTAD. https://www.nodc.noaa.gov/sog/cortad/. Accessed 21 Oct 2020.
Pierce D, Pierce MD. Package ‘ncdf4.’ 2019.
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.
van Woesik R, Kratochwill C. A global coral-bleaching database, 1980–2020. Sci Data. 2022;9:20.
Goodrich B, Gabry J, Ali I, Brilleman S. rstanarm: Bayesian applied regression modeling via Stan. R package version. 2018;2:1758.
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.
Veron J, Stafford-Smith M, DeVantier L, Turak E. Overview of distribution patterns of zooxanthellate Scleractinia. Front Mar Sci. 2015;1.
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.
Gabry J, Mahr T, Bürkner P-C, Modrák M, Barrett M, Weber F, et al. Bayesplot: plotting for Bayesian Models. 2021.
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.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
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.
About this article
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). https://doi.org/10.1186/s12862-022-02011-y
- Climate change
- Coral reefs
- Local stressors