Letter The following article is Open access

Risks to coral reefs from ocean carbonate chemistry changes in recent earth system model projections

, , and

Published 3 July 2013 © 2013 IOP Publishing Ltd
, , Citation K L Ricke et al 2013 Environ. Res. Lett. 8 034003 DOI 10.1088/1748-9326/8/3/034003

1748-9326/8/3/034003

Abstract

Coral reefs are among the most biodiverse ecosystems in the world. Today they are threatened by numerous stressors, including warming ocean waters and coastal pollution. Here we focus on the implications of ocean acidification for the open ocean chemistry surrounding coral reefs, as estimated from earth system models participating in the Coupled Model Intercomparison Project, Phase 5 (CMIP5). We project risks to reefs in the context of three potential aragonite saturation (Ωa) thresholds. We find that in preindustrial times, 99.9% of reefs adjacent to open ocean in the CMIP5 ensemble were located in regions with Ωa > 3.5. Under a business-as-usual scenario (RCP 8.5), every coral reef considered will be surrounded by water with Ωa < 3 by the end of the 21st century and the reefs' long-term fate is independent of their specific saturation threshold. However, under scenarios with significant CO2 emissions abatement, the Ωa threshold for reefs is critical to projecting their fate. Our results indicate that to maintain a majority of reefs surrounded by waters with Ωa > 3.5 to the end of the century, very aggressive reductions in emissions are required. The spread of Ωa projections across models in the CMIP5 ensemble is narrow, justifying a high level of confidence in these results.

Export citation and abstract BibTeX RIS

Content from this work may be used under the terms of the Creative Commons Attribution 3.0 licence. Any further distribution of this work must maintain attribution to the author(s) and the title of the work, journal citation and DOI.

1. Introduction

When carbon dioxide dissolves into seawater it forms carbonic acid, decreasing the ocean's pH, in a process generally referred to as ocean acidification [1]. Rising concentrations of carbon dioxide (CO2) in the atmosphere have already increased the concentrations of hydrogen ions in ocean surface waters by 26%, i.e., a 0.1 unit reduction in global mean ocean surface pH [2]. These hydrogen ions react with carbonate ions to form bicarbonate ions, reducing the saturation of all mineral forms of calcium carbonate in seawater, including aragonite, which corals precipitate to build their skeletons. As the Ωa decreases, so too do calcification rates of corals [3, 4]. The IPCC's Fourth Assessment Report stated that doubling of atmospheric CO2 would reduce calcification of corals by 20–60% [5]. As acidification increases, damage will occur to the coral reefs that serve as protective coastal barriers, breeding and spawning grounds for fish and other marine life, sources of tourism income and hotspots of ocean biodiversity [68].

Thirteen Earth System Models (ESMs) participating in the Coupled Model Intercomparison Project 5 (CMIP5) each have interactive ocean biogeochemistry models that output state variables including total dissolved inorganic carbon (DIC), alkalinity, sea surface temperature (SST), and salinity. Several recent studies have used the CMIP5 ensemble to investigate the implications of temperature projections on coral bleaching [9, 10]. Our study focuses on implications of geochemical projections of chemical stress on corals. We use results from these ESMs to examine the changing saturation state of aragonite (Ωa), the mineral form of calcium carbonate that corals use to make their skeletons, of open ocean waters surrounding approximately 6000 coral reefs. In the remainder of this letter, when we refer to the chemistry of water surrounding coral reefs we mean the chemistry predicted by the Earth system models for the location of the reefs. For locations adjacent to but not in areas resolved by a model, we use the chemistry of water in the nearest location represented in the model.

2. Data and methods

We analyze Ωa surrounding reefs in the CMIP5 simulations under several representative concentration pathway (RCP) scenarios. In the main text, we focus on three scenarios: RCP8.5, a business-as-usual scenario in which very little mitigation is undertaken this century; RCP4.5, under which greenhouse gas emissions peak in mid-century and are reduced and stabilized by late this century; and RCP2.6, a scenario in which emissions are dramatically curtailed in the near future and are negative (CO2 is captured from the air) by the end of the century [11, 12]. Results for another scenario—RCP6.0—which fewer of the ensemble's models simulated, are shown in the supplementary information (available at stacks.iop.org/ERL/8/034003/mmedia).

For each model, Ωa was computed from simulated DIC, alkalinity, temperature, and salinity. But before computation each simulated input field was corrected for its decadal mean bias relative to modern observations, using a slight variant of the approach used in [2, 13]. That is, Ωa was computed after first removing from each model field, the average difference between the model mean during 1989–1998 and the observational reference. For observed fields, we used the GLODAP gridded data product [14] for DIC and alkalinity along with the 2009 World Ocean Atlas climatology for temperature [15], salinity [16], and concentrations of phosphate and silica [17]. Changes to the concentrations of phosphate and silica were assumed to be zero, because all models did not provide those variables4. We computed Ωa using routines based on the standard OCMIP carbonate chemistry adapted for earlier studies [2, 13] to compute all carbonate system variables and use recommended constants from the Guide to Best Practices for Ocean CO2 Measurements [18]. All equilibrium constants are on the total pH scale, except for KS, which is referenced the free pH scale. Total boron is computed as a function of salinity following [19]. To calculate saturation, we used the apparent solubility product for aragonite KsP from [20] and the calcium concentration as a function of salinity [21].

For analysis, model data were regridded to a common 1° by 1°  grid using bilinear interpolation. Reef locations were taken from the Reef Base Global Database [22] and non-reef coral colonies were excluded from the analysis. If a reef was located in a water cell in a particular model, the aragonite saturation value for that cell was assigned to the reef. If a reef was in a masked-value cell (land or regions of the ocean not simulated), but was within 1° of a water cell (i.e., chemistry of one or more of the eight surrounding cells was simulated), the reef was assigned the inverse-distance weighted mean value of all surrounding cells with data. Only reefs with an Ωa value in all ensemble models were included in our analysis. Including interpolated values, 5901 reefs had Ωa data values in all 13 ESMs. Of the ReefBase coral reefs considered, 3821 were excluded, mostly the zone in where there is no GLODAP gridded data (near Indonesia between the Indian and Pacific Oceans). Summary statistics were calculated by aggregating results for a given ESM first, then calculating quartiles and ranges across the ensemble.

See supplementary table S1 (available at stacks.iop.org/ERL/8/034003/mmedia) for a list of models and modeling groups, a summary of the scenarios simulated using each ESM and the years of data available for each model scenario.

3. Results

The aragonite saturation threshold (Ωacrit) for viable reef ecosystems is uncertain [23, 24], but the preindustrial distribution of water chemistry surrounding coral reefs may nevertheless provide some indication of viability. While Ωa levels are different within coral reefs than in the surrounding ocean water, reef chemistry and calcification rates are strongly influenced by the open ocean aragonite saturation [2527]. Thus we use open ocean Ωa relative to preindustrial at reef locations as a proxy for a similar ratio within the reefs.

Historical simulations of ocean chemistry indicate that in preindustrial times, seawater surrounding coral reefs were typically oversaturated with respect to aragonite by a factor of at least 3.5. Given evidence that corals have some biological resilience to acidification [28] and in interest of presenting a conservative analysis, we use this lower bound of preindustrial chemistry conditions as an upper bound viability threshold and also test two lower values of Ωacrit, 3.25 and 3, which are similar to thresholds used in some past studies [2931]. Cao and Caldeira [32] conducted a similar analysis using simulations from a single intermediate-complexity model and found that with atmospheric CO2 concentrations of about 550 ppm or more, all reefs will be in waters with Ωa values of 3.5 or less.

Figure 1 shows global maps of Ωa data and histograms of Ωa at the included coral reef locations, designated on the maps, for four concentrations of atmospheric CO2. The values presented are ten-year means around the year where atmospheric CO2 reaches the designated values of 290, 390, 550 and 900 ppm. Each cell displayed in the maps and each bar in the histograms represent the multimodel ensemble median value. The Ωacrit thresholds that we focus on in our analysis are shown with contour lines on the maps. Error bars on the histograms show the ensemble quartiles and range.

Figure 1.

Figure 1. Aragonite saturations surrounding coral reefs as projected for RCP 8.5. Maps show coral reef locations (in gray) and Ωa in areas where it can be computed from the CMIP5 ensemble and bias corrected with the GLODAP dataset. Histograms show those reefs binned by aragonite saturation of surrounding waters. Panels show ensemble median results for ensemble-member 10-year mean Ωa at: (a) the beginning of the ensemble simulations in the 1860s (287 ppm xCO2); (b) around present day (400 ppm); and for future conditions at (c) 550 ppm (mid-century for RCP 8.5 and end-of-century for RCP 4.5), and (d) 900 ppm (end-of-century under RCP 8.5). Contour lines show Ωacrit thresholds of 3, 3.25 and 3.5; panel (d) has no contours because all ocean waters have Ωa < 3. Histogram error bars show ensemble quartile values and the maximum and minimum range.

Standard image High-resolution image

While the data in figure 1 are taken from the RCP8.5 simulations, the proportion of reefs that lie above a threshold depends primarily on atmospheric CO2 concentration (xCO2) and is not highly path dependent, as discussed below. As such, the time slices shown in the figure can also be considered approximately representative of the corresponding xCO2, regardless of scenario. In the 1860s (287 ppm CO2), all the coral reefs analyzed lie above the three Ωacrit values under preindustrial conditions. Also illustrated are simulated conditions in the approximate present day (390 ppm); at 550 ppm, which corresponds to mid-century for RCP8.5 and end-of-century for RCP4.5; and 900 ppm, which corresponds to the end of the century under RCP8.5.

The importance of the correct value of Ωacrit varies substantially by the scenario analyzed. Figure 2 shows time series of percentage of reefs analyzed above Ωacrit values of 3.5, 3.25 and 3 for the three RCP emissions scenarios in ensemble summary statistics (median, quartiles and range). For RCP8.5, the timing of reefs crossing the thresholds varies with the value of Ωacrit, but by the end of the 21st century all the seawater surrounding reefs analyzed have crossed all three thresholds, exiting the range of ocean chemistry conditions that surrounded coral reefs in preindustrial times. For RCP4.5, the fate of coral reefs strongly depends on the value of Ωacrit; by the end of the 21st century aragonite saturation states have stabilized and approximately 1–5% of reefs are located in ocean waters with Ωa > 3.5, 21–26% are located in ocean waters with Ωa > 3.25, and 42–47% are located in ocean waters with Ωa > 3. Under an aggressive emissions reduction scenario like RCP2.6, both the thresholds and the timescales for reef resilience are pertinent. Saturation states reach minimum values at mid-century with 32–37% of reefs located in ocean waters with Ωa > 3.5, 47–57% in ocean waters with Ωa > 3.25, and 77–87% in ocean waters with Ωa > 3; afterwards, Ωa begins to increase again.

Figure 2.

Figure 2. Percentage of coral reefs surrounded by open ocean waters above aragonite saturation thresholds as a function of time in three RCP scenarios. The ensemble median (lines), quartiles (darker shading) and range (lighter shading) are displayed for aragonite saturation thresholds (Ωacrit) of 3.5, 3.25 and 3. Panels show RCP emissions scenarios (a) 8.5, (b) 4.5 and (c) 2.6. The future of coral reefs depends on both future emissions trajectories and the sensitivity of reefs to changes in Ωa.

Standard image High-resolution image

Figure 3 shows the same results from all scenario simulations as a function of atmospheric CO2 concentration instead of time. The ensemble median and quartile Ωa values align closely between emissions scenarios, though the range varies considerably between the low and high emissions scenarios. Variability in the outlier models appears to depend on the rate of change of xCO2 as well as its absolute concentration. Between the point at which the population of reefs begins crossing the Ωacrit threshold and that at which about 10% reefs are below the threshold, the decrease in the percentage of reefs above Ωacrit is roughly linear with xCO2 increase.

Figure 3.

Figure 3. Percentage of coral reefs surrounded by open ocean waters above aragonite saturation thresholds as a function of atmospheric CO2 concentration. The ensemble median (lines), quartiles (darker shading) and range (lighter shading) are displayed for aragonite saturation thresholds (Ωacrit) of 3.5, 3.25 and 3. The narrow spread across all models and emissions scenarios considered indicates that xCO2 is a good predictor of the distribution of aragonite saturation around coral reefs.

Standard image High-resolution image

Results for individual models are plotted in the supplementary materials (see supplementary figures S1 and S2 available at stacks.iop.org/ERL/8/034003/mmedia). The 13 models' ocean aragonite saturation responses to rising atmospheric CO2 track closely across the ensemble and display similar interannual variabilities. This agreement suggests that these ocean chemistry projections, unlike some other simulated climatological indicators such as precipitation, are fairly robust in the current generation of climate models.

4. Discussion

Many factors including light, temperature, salinity, nutrients, and sediment loading affect the distribution of coral reefs, and some of these factors are correlated with aragonite saturation state (Ωa). Therefore, the absence (or near-absence) of coral reefs in waters with Ωa levels less than 3, 3.25 or 3.5, does not necessarily mean that a threshold aragonite saturation limits the biogeographic distribution of coral reefs. However, while coral reefs have sustained existence under higher temperatures in the Earth's past, such as during the Cretaceous, they have not thrived far outside of the tropics [33] where high aragonite saturations have prevailed. Such regions of high aragonite saturation are projected to disappear by the end of the century under RCP8.5. In the light of experimental data showing reduced calcification with decreasing aragonite saturation, the lack of overlap between the aragonite saturation levels surrounding coral reefs before the industrial revolution and those in a future ocean under continued CO2 emissions provide reasons to doubt future reef sustainability. In the absence of substantial reductions in CO2 emissions, surrounding waters will go outside the bounds of the open ocean saturation states that surrounded all coral reefs before the industrial revolution.

Reefs that did not fall within 1° of a water cell in all 13 CMIP5 GLODAP-corrected datasets are excluded from our analysis. Coral reefs in certain regions are therefore not represented in our results. In particular, the reef-rich Indo-Pacific between the Korean peninsula and the northern coast of Australia was masked in the ensemble. Such lower latitude regions tend to have higher Ωa values on average, which could indicate that our results are slightly biased.

A number of factors in addition to aragonite saturation will influence the fate of coral reefs in a changing climate, perhaps even more so than ocean acidification. Several recent analyses have examined potential future heat stress as predicted by the CMIP5 ensemble [9, 10]. Other human stresses on the environment in and around coral reefs also constitute significant threats to reef viability in coming decades, including pollution, overdevelopment and overfishing [3437]. The correlation (or lack thereof) between changing carbonate chemistry and these other projected stresses is not clear, although the analysis in Meissner et al [31, 37] using one intermediate-complexity model suggests that the reefs that are chemically stressed first under the RCP scenarios are often different than those that are thermally stressed first [38].

Our results indicate that if civilization continues along its current carbon dioxide emissions trajectory (which is currently steeper than even RCP8.5), by the end of the century all existing coral reefs will be surrounded by ocean chemistry conditions well outside of preindustrial (or even present day) conditions. This result is consistent across all of the Earth System Models simulating ocean biogeochemistry in the CMIP5 ensemble. On the timescales and within the range of scenarios considered here, carbonate chemistry of waters surrounding reefs is closely related to atmospheric CO2 concentrations and is not sensitive to the pathway by which those concentrations are reached. With deep emission cuts, the fraction of reefs that remains sustainable from a chemical point of view depends sensitively upon the critical value of Ωa. Only with large emissions abatement efforts will a significant fraction of reefs remain in ocean waters with aragonite saturation states of three or more, and only under the most aggressive policies, possibly including direct capture of CO2 from the atmosphere, will reef-amenable carbonate chemistry conditions be preserved around a majority of reefs.

Acknowledgments

We acknowledge the World Climate Research Programme's Working Group on Coupled Modelling, which is responsible for CMIP, and we thank the climate modeling groups (listed in table S1 of the supplementary materials available at stacks.iop.org/ERL/8/034003/mmedia) for producing and making available their model output. For CMIP the US Department of Energy's Program for Climate Model Diagnosis and Intercomparison provides coordinating support and led development of software infrastructure in partnership with the Global Organization for Earth System Science Portals. We thank Long Cao and two anonymous reviewers for helpful comments on the manuscript. To produce the derived CMIP5 data for saturation states, JCO acknowledges support from two EU FP7 projects: EPOCA (grant 211384) and CARBOCHANGE (grant 264879).

Footnotes

  • In the low latitude oceans, phosphate and silica generally have low concentrations. Together their combined alkalinity affects the computed saturation state for aragonite by less than 0.005 in terms of the zonal mean. However, in some areas with high nutrients such as the equatorial Pacific, the contribution is larger. We account for phosphate and silicate contributions in our calculations. As our baseline for those nutrients, we use the gridded data product from the World Ocean Atlas [17], and those are combined with other input, namely for DIC and Alk from GLODAP [14] as well as temperature and salinity from [15]. However, when adding the perturbations from the model simulations we neglect changes in phosphate and nutrients. We consider that those changes are relatively small and will have very little impact on computed saturation states.

Please wait… references are loading.