(C) PLOS One This story was originally published by PLOS One and is unaltered. . . . . . . . . . . Cloudiness delays projected impact of climate change on coral reefs [1] ['Pedro C. González-Espinosa', 'Department Of Geography', 'University Of British Columbia', 'Vancouver', 'Bc', 'Institute For The Oceans', 'Fisheries', 'Simon D. Donner'] Date: 2023-03 Abstract The increasing frequency of mass coral bleaching and associated coral mortality threaten the future of warmwater coral reefs. Although thermal stress is widely recognized as the main driver of coral bleaching, exposure to light also plays a central role. Future projections of the impacts of climate change on coral reefs have to date focused on temperature change and not considered the role of clouds in attenuating the bleaching response of corals. In this study, we develop temperature- and light-based bleaching prediction algorithms using historical sea surface temperature, cloud cover fraction and downwelling shortwave radiation data together with a global-scale observational bleaching dataset observations. The model is applied to CMIP6 output from the GFDL-ESM4 Earth System Model under four different future scenarios to estimate the effect of incorporating cloudiness on future bleaching frequency, with and without thermal adaptation or acclimation by corals. The results show that in the low emission scenario SSP1-2.6 incorporating clouds into the model delays the bleaching frequency conditions by multiple decades in some regions, yet the majority (>70%) of coral reef cells still experience dangerously frequent bleaching conditions by the end of the century. In the moderate scenario SSP2-4.5, however, the increase in thermal stress is sufficient to overwhelm the mitigating effect of clouds by mid-century. Thermal adaptation or acclimation by corals could further shift the bleaching projections by up to 40 years, yet coral reefs would still experience dangerously frequent bleaching conditions by the end of century in SPP2-4.5. The findings show that multivariate models incorporating factors like light may improve the near-term outlook for coral reefs and help identify future climate refugia. Nonetheless, the long-term future of coral reefs remains questionable if the world stays on a moderate or higher emissions path. Citation: González-Espinosa PC, Donner SD (2023) Cloudiness delays projected impact of climate change on coral reefs. PLOS Clim 2(2): e0000090. https://doi.org/10.1371/journal.pclm.0000090 Editor: Lin Liu, First Institute of Oceanography Ministry of Natural Resources, CHINA Received: July 4, 2022; Accepted: December 1, 2022; Published: February 8, 2023 Copyright: © 2023 González-Espinosa, Donner. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. Data Availability: The coral bleaching database is available at https://zenodo.org/record/6780843, the python and R codes, and CDO script is available in the repository https://doi.org/10.5281/zenodo.6795835. The NOAA Coral Reef Watch CoralTemp Version 3.1 Daily Global 5 km is publicly available at https://www.star.nesdis.noaa.gov/pub/sod/mecb/crw/data/5km/v3.1/nc/v1.0/daily/sst/. The Downward short-wave radiation at the surface (RSDS) and total cloud area (CLT) data from the Geophysical Fluid Dynamics Laboratory’s ESM4 model (GFDL-ESM4) are publicly available at https://esgf-node.llnl.gov/search/cmip6/ (source ID = GFDL-ESM4). Funding: This work is supported by the Mexican Council for Science and Technology (CONACYT) doctoral scholarship (number 438530) to P.C.G.E. and by the Natural Sciences and Engineering Research Council of Canada (NSERC) Discovery Grant to S.D.D. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. Competing interests: The authors have declared that no competing interests exist. 1 Introduction Climate change is a key risk for coral reefs [1,2]. Tropical corals are composed of coral polyps and microalgal symbionts living within the corals’ tissue which provide up to 95% of their photosynthetic products to the coral [3]. This symbiotic relationship can be disrupted by increases in ocean temperatures of as little as 1–2°C above those usually experienced in summer, resulting in the phenomenon known as coral bleaching [2,4]. Climate change has led to more frequent and more extreme mass coral bleaching events worldwide over the last four decades, and widespread loss of living corals [2,5,6]. Moreover, coral reefs have been projected to face dangerously frequent bleaching conditions by mid-century, even in a 1.5°C warming scenario [7,8]. The heterogeneity of coral responses to environmental stress represents a challenge for predicting mass coral bleaching events in real-time and for projecting the response of corals to climate warming [9,10]. In addition to thermal or heat stress, other stressors such as solar radiation, turbidity, salinity, air exposure, or pollution can effect the coral-algal relationship [11–13]. Although elevated temperature is often reported as the primary cause of mass coral bleaching, the interaction between temperature and other environmental variables modulates the bleaching response [14]. Light exposure in particular plays a central role for corals [15,16]. Light is critical to coral growth but excess exposure can inhibit or damage the photosynthetic processes of the microalgal symbionts [17–19]. High irradiance, when coupled with high temperatures, can disrupt the balance between light absorption and electron transport rates, promoting the production of reactive oxygen species (ROS) that can trigger photoinhibition, impair the photosynthetic apparatus of the microalgal symbionts [20] and, in turn, exacerbate thermal stress [15,16,21]. Moreover, the reductions in zooxanthallae pigment density and increased light absorption can cause a positive feedback which worsens the effects of elevated temperatures [22]. To minimize the potential photodamage due to excess light, corals and their symbionts have acclimation strategies including regulating the density of symbionts and/or the ratio of photosynthetic pigments [23], restructuring their skeletal morphologies [24], hosting different Symbiodinium species [25,26] or increasing antioxidant contents [27]. While light reaching the surface of the corals is modulated by properties of seawater itself, field observations indicate that clouds can reduce the amount of incoming solar radiation and limit the effects of thermal stress on corals [28–31]. A recent global analysis found that incorporating cloud cover can improve the accuracy of mass coral bleaching prediction [32]. While models are beginning to include evolutionary [33] and ecological processes [34,35] in projecting the future of coral reefs, the forcing variables in those models are generally still restricted to temperature. There has yet to be an analysis of the integrated effect of thermal stress, solar radiation and cloudiness on future bleaching conditions. In this study, we test for the first time whether including cloudiness in coral bleaching projections models affects the projected frequency of mass coral bleaching conditions at regional and global scales under climate change. First, building on previous work [32], we use a historical mass bleaching dataset to develop a new bleaching prediction algorithm incorporating thermal stress and cloudiness that are designed for use with climate model output. Second, we contrast the projected frequency of mass coral bleaching through the year 2100 with and without the role of clouds under four future scenarios, using CMIP6 output from the GFDL-ESM4 Earth System Model. To test for robustness, the analysis is repeating using an algorithm incorporating incoming shortwave radiation rather than cloudiness. Finally, we assessed whether thermal adaptation or acclimation by corals and their symbionts could influence the role of clouds in future bleaching projections. The results indicate that incorporating light, in the form of either cloudiness or incoming shortwave radiation, can delay bleaching projections in low-to-moderate emissions scenarios and help identify future climate refugia. 2 Methods 2.1 Historical bleaching reports The coral bleaching data used in this study were derived from version 2.0 [36] of a high-resolution global mass coral bleaching database [6] that includes reports through the year 2017. Each report includes a categorized bleaching severity from levels 0 to 3 based on the protocol of ReefBase, where category 0 refers to no bleaching (<1% of observed living coral area is bleached), category 1 refers to mild bleaching (1–10% bleached), category 2 refers to moderate bleaching (11–50% bleached) and category 3 refers to severe bleaching (>50% bleached). To align with the available historical high-resolution climate data [see sections 2.2 & 2.3], we focused on all bleaching reports from 2001 through 2017 (n = 33,768) which reported the month or date of observation. 2.2 Historical sea surface temperature and thermal stress To describe historical coral thermal exposure, sea surface temperature (SST) data from 0.05° x 0.05° lat-long resolution CoralTemp V3.1 [37] was used to compute the Degree Heating Month (DHM) value corresponding to each coral bleaching report in the database. The DHM is a monthly-scale version of the Degree Heating Week index used in real-time bleaching prediction; the monthly rather than weekly index from [32] is employed in order to develop an algorithim that aligns with the available resolution of most archived climate model output [38–40]. Here, DHM is computed by summing the SST anomaly in excess of the maximum monthly mean (MMM), a maximum value from a climatology, over a 3-month rolling window [7,40,41] using the following formula: The MMM Max [42], represents the mean of the warmest monthly SST from each year during the period 1985–2004 calculated for each grid cell. The MMMmax is a variation of the MMM threshold that is used in recent models [33,40] because it better characterizes the bleaching threshold in regions where the timing of the seasonal peak in SST varies from year to year. 2.3 Historical cloudiness and incoming radiation To represent the amount of light energy that corals receive, 1° x 1° lat-long total cloud area (CLT) and incoming shortwave radiation at the surface data (represented by the variable downwelling shortwave radiation or RSDS) for the 2001–2017 period were retrieved from CERES project (CERES-EBAF Edition 4.0; [43,44]). Since corals acclimate to background light conditions, monthly anomalies for each variable (CLTa, RSDSa) were computed by subtracting monthly values from their respective 2001–2017 monthly climatologies [32]. This data can be used to test the hypothesis that exposure to excessive light (i.e. positive RSDSa or incoming shortwave radiation anomaly) for a given location and time of year due to anomalously low cloud cover (i.e. negative CLTa or total cloud area anomaly) correlates with higher bleaching severity, and vice versa. Using the RSDSa and the CLTa data provided two different methods for characterizing incoming light in a bleaching prediction algorithm. The effect of each variable was tested separately as a check on the robustness of the underlying hypothesis that incorporating light influences the accuracy of bleaching prediction. 2.4 Bleaching prediction model development We used the historical bleaching observations together with the climate variables to develop three different bleaching prediction models by performing mixed-effect models with a random effect structure (lmerTest package in R [45]). In the first or control model, only thermal stress (DHM) was used as a fixed effect, with geographical position set as a random intercept to control the spatial structure in the data. In the other two models, thermal stress (DHM) and either cloud anomaly (CLTa) or incoming shortwave radiation anomaly (RSDSa) were used as fixed effects, with geographical position again set as a random intercept. The models also include an interactive term, as follows: Where: severity = is the predicted value of the projected severity of bleaching β 0 = is the severity intercept β 1 X 1 = the regression coefficient (β 1 ) of the first independent variable (e.g., DHM) β 2 X 2 = the regression coefficient of the second independent variable (e.g., CLTa, the total cloud cover anomaly) β 12 X 12 = the regression coefficient (β 12 ) of the interaction of both variables (e.g., DHM:CLTa) Latitude was selected to represent geographical position, because sensitivity analyses performed with latitude, longitude, and both latitude and longitude as random intercepts produced similar results (Table A in S1 Text). To assess the performance of the bleaching prediction models we computed the accuracy (the percentage of correctly classified instances out of all instances) but also the Cohen´s Kappa (the difference between the observed overall accuracy of the model and the overall accuracy that can be obtained by chance) for which higher values represent better model performance ([46]). 2.5 Future projections We applied the two multiple variable bleaching severity models (Section 2.4) to 1985–2100 output from the Geophysical Fluid Dynamics Laboratory’s ESM4 model (GFDL-ESM4) for four of the socio-economic pathways (SSPs) employed by the climate modelling community. Analysis was conducted for all the 1° x 1° lat-long model cells worldwide containing coral reefs, according to the Millennium Coral Reef Mapping Project [47]. The GFDL-ESM4 model was selected because it has moderate climate sensitivity (1.6°C for transient climate sensitivity and 2.6°C for equilibrium climate sensitivity) among ESMs and coupled ocean-atmosphere general circulation models (GCMs) from the Coupled Model Intercomparison Project 6 (CMIP6; [48]). The workflow was as follows: independent variables (DHM, CLTa, RSDSa) required by the bleaching severity models (Section 2.4) were computed from GFDL-ESM4 output. Monthly values for the SST, total cloud area (CLT) and incoming shortwave radiation at the surface (RSDS) at 1° x 1° lat-long grid resolution were retrieved from the https://esgf-node.llnl.gov/search/cmip6 node for the four SSPs described above and from a historical all-forcing simulation. Historical (1985–2014) and projected (2015–2100) output data were concatenated to create a complete-time series. As is common in climate change impacts research, we corrected for biases in model output (Fig A in S1 Text) using the delta method [49]; raw output for SST, CLT and RSDS were adjusted by adding the model anomalies (simulated projected values–simulated historical climatology) to the observed climatologies for respective variables. The adjusted SST values were then used to compute a DHM time series, using the method described in Section 2.2, and the adjusted CLT and RSDS values were used to compute the total cloud area anomaly (CLTa) and the incoming shortwave radiation anomaly (RSDSa) time series, as described in Section 2.3. For each grid cell containing coral reefs, we used the coefficients from the bleaching projection models (Section 2.4) to compute the projected bleaching severity based on the computed DHM, CLTa or RSDSa time series (1985–2100). Given the documented uncertainty in model cloud projections [50], a sensitivity analysis was performed in which the bleaching severity calculations were repeated using projected DHMs but observed CLTa and RSDSa from the climatology rather than the projected values. To assess the effects of incorporating cloudiness on future projections of coral bleaching frequency and severity, the results from the bleaching severity models are compared to that of the common DHM-only method used in past climate change projection studies. Past studies used a threshold of DHM ≥ 2°C·month to indicate the likelihood of occurrence of severe bleaching conditions [33,38–40]. Here, a threshold of DHM ≥ 2.5°C·month was employed, since the DHM ≥ 2°C·month threshold tends to overpredict the bleaching frequency when interacting with another variable (Fig B in S1 Text), i.e. even though cloudiness anomalies reduce bleaching severity, it depends on its magnitude and the interaction effect with DHM [32]. Following previous studies, we express results in terms of year that the probability of bleaching conditions in each grid cell exceeds 0.2 (bleaching occuring more than twice per decade) [33,39,40] and the year that bleaching conditions in each grid cell become an annual occurrence (p = 1) van Hooindonk [38]. In those results, the year refers to the midpoint of a 10-year running analysis, e.g., for bleaching probability for the year 2020 is based on occurrences from the period 2015–2024. To examine if thermal adaptation or acclimation by corals,their algal symbionts or human interventions [51–54] changes the effect of incorporating cloudiness on bleaching projections, a series of simple additional sensitivity analyses were conducted. The bleaching projections were recomputed three times, assuming a +0.5°C, +1.0°C and +1.5°C adjustment in the threshold above which thermal stress accumulates (MMM max ), following estimates of maximum likely adaptation employed in previous projection studies [33,39]. For example, if the thermal threshold is estimated to be 29.5°C (for a certain moment), it is expected that some corals would bleach under these conditions, while some others would resist because their threshold is 0.5°C higher (30°C) due to natural physiological acclimation, assisted evolution or other mechanisms. Finally, to explore the temporal and spatial change in cloudiness but also to visualize the magnitude of the change between scenarios we compared the total cloud trend (CLT trend ), for the period 2015–2100. All data retrievals, processing and statistical tests were conducted using Python (version 3.7.4), R project (version 3.6.3), Climate Data Operators (CDO; version 1.9.8) and ArcGIS Pro (version 2.6). 4 Discussion and conclusions Ocean warming has been clearly documented as the key driver of the observed increase in mass coral bleaching events over the last four decades [2,5,6]. However, geographical patterns in mass coral bleaching are increasingly being linked to multiple environmental variables including, but not restricted, to temperature [9,55–57]. Global coral bleaching projections, including those driving summary statements in the recent Intergovernmental Panel on Climate Change reports [58,59], have nevertheless relied on model-based estimates of bleaching-level thermal stress [7,8,38,39,60]. Although there have been recent advances in projecting coral ecological and evolutionary responses to climate change [33–35],the location of potential future climate refuges [61], such large-scale projections are still largely reliant on temperature as the independent variable. Those hindcast or smaller scale projection studies which included multiple environmental variables have generally excluded light. Here, we find that incorporating cloud cover into global projections can delay the onset of high frequency bleaching conditions, based on the documented effect that clouds have had on mitigating against thermal stress. For example, in the most optimistic scenario analyzed in this study, SSP1-2.6, clouds could delay the onset of high-frequency bleaching by multiple decades, although the majority of the world’s reefs still experience dangerously frequent bleaching conditions by the end of the century (~90% of reef cells with the DHM-only model vs. ~70% with the cloud model). However, the results from higher emissions scenarios point to a limit to the mitigating effect of clouds, as the capacity for DHM to increase exceeds that of clouds (i.e. the model incorporates cloud cover, which has a maximum value of 100%). Even in the moderate scenario SSP2-4.5, more than 95% of reef cells would experience high frequency bleaching conditions by 2050 (Fig 1). Bleaching conditions are projected to become dangerously frequent across time and space in the absence of adaptation, physiological acclimation, shifts in coral community structure or human interventions like active restoration or assisted evolution [58,62]. More advanced models have found that incorporating evolutionary processes [33,40,63] or population dynamics [34,35] improves the forecast for the world’s coral reefs. The coarse adaptation sensitivity test applied here finds that the mitigating effect of clouds increases assuming some positive adjustment of the thermal threshold for corals, but that this mitigating effect of clouds is still limited in moderate or high emissions scenarios. For example, in the moderate emission scenario SSP2-4.5, the cloud effect combined with a +1.5° C increase in the bleaching thresholds still leads to ~75% of reef cells experiencing bleaching conditions at least once every five years by the end of the century (Fig 3B). Although this research does not incorporate ecological or evolutionary processes, the difference in the projected bleaching frequencies between simulations, including cloud effects and both cloud effects and adaptation or acclimation, highlight the importance of considering additional driving variables in climate change impact projections. Moreover, results from studies like this could help to reinforce the accomplishment of the objective of curbing greenhouse gas emissions to tackle climate change. Consideration of additional driving variables, like incoming light and clouds, is particularly critical for identifying potential climate change refugia. These projections agree with past work [32] finding that the central and south Pacific, including parts of French Polynesia, and the central equatorial Pacific, including the Republic of Kiribati’s Gilbert Islands, the Phoenix Islands and Line Islands, are regions where coral reefs have the highest likelihood of some cloud protection during future thermal stress. Although these regions are still at risk, and have experienced extensive bleaching events, field evidence does confirm that cloudiness may there have reduced coral response to past thermal stress events [28,31]. Coral bleaching is often associated with extreme ENSO conditions, the enhanced convective activity in these parts of the Pacific, associated with anomalously warm SSTs, also drive cumulus cloud development [64]. An additional consideration, not included in this study, is that light reaching the coral surface can also be influence by turbidity which increases light attentuation in the water column [65]. A study based on present-day turbidity concluded that parts the Coral Triangle, Northern Great Barrier Reef, or some Indian Ocean reefs in Kenya and Tanzania may receive protection from thermal stress [66]. Notably, the potential refugia identified in these models incorporating light differ from those identified in studies only employing thermal stress [39,61]. This analysis provides the first assessment on how the interaction of thermal stress and solar radiation could influence the response of coral reefs at regional and global scales to climate change. Although output from only one ESM was employed to test the hypothesis that including cloudiness would alter future bleaching projections, the similarity between the results using observed vs. model projected cloudiness suggests that the simulated influence of clouds on bleaching projections would be broadly similar with other models. In reality, at the local level, the effect of clouds on bleaching likelihood will depend on factors like cloud height, cloud composition, and light attenuation in the water column, which vary at finer spatial resolution than models can simulate. The negative feedback whereby high SST can enhance convection and low-level cloudiness, which then reduces downward short-wave radiation, may protect some low latitude reefs, as has been observed in Kiribati [28]. Further analyses incorporating cloud types, temperatures, height, and radiative properties for different regions would unveil fine-tuned responses of corals due to thermal stress and cloudiness. Performing an analysis using fine-scale variability of cloud type and other properties could be done as a case study in a location with high resolution bleaching severity data (e.g., the Mesoamerican Barrier Reef System or the Great Barrier Reef). Notwithstanding, while a limitation of this research is that cloud properties (like mentioned above) were not considered our findings broadly show that clouds may spare some coral reefs from dangerously frequent bleaching conditions in the near term. Even with the mitigating effect of clouds, this study finds that coral reefs remain under existential threat in the long-term without sharp reductions in greenhouse gas emissions. Acknowledgments The authors thank Alex Tso and Alejandra Virgen-Urcelay for their assistance in developing the coral bleaching database. We also specifically thank Mark Eakin and Denise Devotta of NOAA Coral Reef Watch for contributing data for the 2014–2016 global coral bleaching event to the database, and William Skirving of NOAA Coral Reef Watch for assistance with the CoralTemp data. We thank John Dunne and Peter Kalmus for helpful insights in developing the research. We also thank Andrés Cisneros-Montemayor and Xinru Li for their comments and suggestions. [END] --- [1] Url: https://journals.plos.org/climate/article?id=10.1371/journal.pclm.0000090 Published and (C) by PLOS One Content appears here under this condition or license: Creative Commons - Attribution BY 4.0. via Magical.Fish Gopher News Feeds: gopher://magical.fish/1/feeds/news/plosone/