Simulation of global ocean acidification and chemical habitats of shallow- and cold-water coral reefs

Transcription

1 Available online at ScienceDirect Advances in Climate Change Research 5 (2014) 189e196 Simulation of global ocean acidification and chemical habitats of shallow- and cold-water coral reefs ZHENG Mei-Di, CAO Long* Department of Earth Science, Zhejiang University, Hangzhou , China Received 29 January 2015; revised 13 March 2015; accepted 26 May 2015 Available online 12 June 2015 Abstract Using the UVic Earth System Model, this study simulated the change of seawater chemistry and analyzed the chemical habitat surrounding shallow- and cold-water coral reefs from the year 1800 to 2300 employing RCP2.6, RCP4.5, RCP6.0, and RCP8.5 scenarios. The model results showed that the global ocean will continue to absorb atmospheric CO 2. Global mean surface ocean temperature will rise 1.1e2.8 K at the end of the 21st century across RCP scenarios. Meanwhile, the global mean surface ocean ph will drop 0.14e0.42 and the ocean surface mean concentration of carbonate will decrease 20%e51% across the RCP scenarios. The saturated state of sea water with respect to calcite carbonate minerals (U) will decrease rapidly. During the pre-industrial period, 99% of the shallow-water coral reefs were surrounded by seawater with U > 3.5 and 87% of the deep-sea coral reefs were surrounded by seawater with aragonite supersaturation. Within the 21st century, except for the high mitigation scenario of RCP2.6, almost none shallow-water coral reefs will be surrounded by seawater with U > 3.5. Under the intensive emission scenario of RCP8.5, by the year 2100, the aragonite saturation horizon will rise to 308 m under the sea surface from 1138 m at the preindustrial period, thus 73% of the cold-water coral reefs will be surrounded by seawater with aragonite undersaturation. By the year 2300, only 5% of the cold-water coral reefs will be surrounded by seawater with aragonite supersaturation. Keywords: Simulation research; Aragonite saturation state; Ocean acidification; Shallow-water coral reefs; Cold-water coral reefs 1. Introduction Since the beginning of the industrial revolution, the CO 2 concentration in atmosphere has increased rapidly, because of human activities such as fossil fuel combustion and land use. Anthropogenic CO 2 emissions have reached 545 Pg C from 1750 to 2011 (IPCC, 2013). Not all of CO 2 emissions remain in the atmosphere, 26% of them are absorbed by the ocean and * Corresponding author. address: longcao@zju.edu.cn (CAO L.). Peer review under responsibility of National Climate Center (China Meteorological Administration). Production and Hosting by Elsevier on behalf of KeAi 28% are absorbed in terrestrial soil (Sabine et al., 2004). Recently, both observation and simulation studies have agreed that the ocean absorbs a large amount of CO 2, which mitigates global climate change. However, this absorption is not entirely harmless. The average ph of the sea surface has dropped from 8.2 to 8.1 since pre-industrial time, which means that the ocean hydrogen ion concentration has increased by 26% (Gattuso and Hansson, 2011). Doney et al. (2009) reported the ocean carbon chemical equilibrium formula as: CO 2ðatmosÞ %CO 2ðaqÞ þ H 2 O%H 2 CO 3 %H þ þ HCO 3 %2Hþ þ CO 2 3 : ð1þ According to this equation, the increase of hydrogen ions (H þ ) will result in the decrease of carbonate ions concentration (CO 2 3 ), which will further reduce the calcium carbonate /Copyright 2015, National Climate Center (China Meteorological Administration). Production and hosting by Elsevier B.V. on behalf of KeAi. This is an open access article under the CC BY-NC-ND license ( This is an English translational work of an article originally published in Progressus Inquisitiones De Mutatione Climatis (Chinese). The original article can be found at /j.issn

2 190 ZHENG M.-D., CAO L. / Advances in Climate Change Research 5 (2014) 189e196 (CaCO 3 ) saturation state. These changes in sea water will harm marine calcified organisms, especially the coral reefs (Doney et al., 2009; Fine and Tchernov, 2007; Fabry et al., 2008; Guinotte and Fabry, 2008). Coral reefs are an important component of marine ecosystems and are composed of aragonite (a relatively soluble mineral form of calcium carbonate). Many laboratory experiments showed that the decline of the ocean aragonite saturation state would lead to the decrease in the shallow-water coral calcification rate, because it is difficult for corals to extract calcium ions and bicarbonate ions from the surrounding seawater to form bones and shells (Ca 2þ þ 2HCO 3 %CaCO 3 þ CO 2 þ H 2 O). In addition, coldwater coral reefs also suffer from the effects of this change in the marine environment (Langdon et al., 2003). Guinotte et al. (2006) suggested that the global distribution of coldwater coral in the deep sea is limited by the surrounding seawater aragonite saturation state. Therefore, an increase in the CO 2 concentration in the atmosphere can have important influences on the marine environment and coral reef systems. Anthropogenic CO 2 emissions will decrease the ocean surface ph and the concentration of carbonate ions (Caldeira and Wickett, 2005; Orr et al., 2005; Cao and Caldeira, 2008; Steinacher et al., 2009). This development will directly affect the ocean aragonite saturation state. Aragonite saturation state (U) is defined as (Feely et al., 1988): U ¼ Ca 2þ. CO 2 3 K * sp ; ð2þ which is calculated from the calcium ion concentration ([Ca 2þ ]), carbonate ion concentration (½CO 2 3 Š) and equilibrium thermodynamic solubility product (Ksp ). The [Ca2þ ] and ½CO 2 3 Š are calculated from the sea water salinity, total alkalinity (TAlk), and dissolved inorganic carbon (DIC). The solubility product Ksp is calculated based on seawater temperature, salinity, and pressure (Mucci, 1983), which is defined as the product of [Ca 2þ ] and ½CO 2 3 Š concentrations at the saturation state. When U < 1 the seawater is undersaturated with respect to aragonite, while U > 1 means it is supersaturated. The aragonite saturation horizon is defined as the depth which U ¼ 1. In the ocean today, the marine aragonite saturation state decreases from the surface to the depth. Thus, the water above the saturation horizon is supersaturated and the water under the saturation horizon is undersaturated. There are numerous simulation studies about ocean acidification under various CO 2 emission scenarios. Using the ocean circulation model of the Lawrence Livermore National Laboratory, Caldeira and Wickett (2005) found that the global mean ocean surface ph would drop 0.3e0.5 at the end of the 21st century under the SRES scenario. Orr et al. (2005) commented that surface seawater would be undersaturated with respect to aragonite at the middle of the 21st century under the IS92a scenarios, which was based on the simulation results of 13 models, which would have a negative impact on polar shell plankton. In a recent research, it was found that 20% of the surface area of the Canadian Basin has been undersaturated with respect to aragonite (Robbins et al., 2013). This mismatch between observations and models is probably due to the low resolution of global models and the inability of the global models to accurately simulate sea ice. The actual speed of the decline of Arctic sea ice is faster than the forecast, because of global warming (Stroeve et al., 2007). On the basis of UVic ESCM simulation, Cao and Caldeira (2008) speculated that the ocean mean ph at high latitude would drop by more than 0.2 when the atmospheric CO 2 concentration stabilized at , and in this case, 7% of the area of the south ocean (south of 60 S) would be undersaturated in aragonite. Fine and Tchernov (2007) experimentally determined that Scleractinian coral would be huge, soft, and skeleton-free in acidified seawater (ph ¼ 7.4). Chinese scientists have conducted several studies using the global carbon cycle model simulation. For example, Xu and Li (2009) simulated the global ocean anthropogenic CO 2 uptake. Cao et al. (2014a) and Wang et al. (2014) used UVic ESCM to simulate and analyze how climate sensitivity could affect the uptake of CO 2 by the ocean under the RCP8.5 scenario. Cao et al. (2014b) analyzed the response of ocean acidification to a gradual increase and decrease of atmospheric CO 2 and found that marine ecosystems would not respond well to recover to their natural chemical habitats even if the atmospheric CO 2 content is lower than future predictions. On the basis of the previous researches, our study used an intermediate complex earth system climate model to compare various ocean acidification rates from before the industrial revolution to the year 2300 under four CO 2 Representative Concentration Pathway scenarios (RCPs). These scenarios have been widely employed in the IPCC AR5 report. At the same time, this study also analyzed the effect of ocean acidification on both shallow- and cold-water coral reefs. Besides being a supplement to the CMIP5 multi-models prediction in the AR5 report, the work also represents the first attempt to analyze the chemical habits of cold-water coral reefs under four RCPs. 2. Model and method 2.1. Model description Developed at the University of Victoria at Canada, the UVic Earth System Climate Model is an intermediate complexity climate model, which consists of a 3D ocean general circulation model with a spherical grid resolution of 3.6 by longitude and 1.8 by latitude with 19 vertical layers in the ocean and one-layer in the atmosphere with an energymoisture balance. It is coupled with a sea-ice model, an atmospheric model and a land-ice model (Weaver et al., 2001). The terrestrial carbon cycle model is based on the TRIFFID land surface and dynamic vegetation scheme of Hadley Centre Met Office (Meissner et al., 2003). In addition, the ocean carbon cycle consists of a CO 2 air-sea interaction process, a marine inorganic carbon process (Orr et al., 1999) and a marine organic carbon process (Schmittner et al., 2008). The organic carbon process includes a simple marine ecosystems, which has nutrient (PO 4 3,NO 3 ), interaction of phytoplankton and zooplankton and a feedback for the marine carbon cycle. This model has been widely used in numerous studies of

3 ZHENG M.-D., CAO L. / Advances in Climate Change Research 5 (2014) 189e climate change, ocean biochemistry cycles (Schmittner et al., 2008), climate feedback of ocean carbon cycles (Zickfeld et al., 2013) and ocean acidification (Cao and Caldeira, 2008; Matthews et al., 2009) Simulation experiments We ran the UVic model for 10,000 model years with an atmospheric CO 2 concentration of as being the CO 2 level of the pre-industrial time (Indermühle et al., 1999). Next, we set this state as the initial condition for the nominal year 1800 and simulated the climate change from 1800 to From 1800 to 2005, the model was driven by the historical concentration of atmospheric CO 2. After 2005, it was forced by the CO 2 concentration in the atmosphere as prescribed by four RCP scenarios (RCP2.6, RCP4.5, RCP6.0, RCP8.5, the numbers behind RCP mean that global radiation forcing will reach 2.6, 4.5, 6.0, 8.5 W m 2 at the year 2100, while atmosphere CO 2 concentration will reach , , , , respectively). mean DIC and TAlk were 2254 and 2363 mmol kg 1, while during the same corresponding period, the model global mean DIC and TAlk were 2242 and 2368 mmol kg 1. The model error of the DIC was 0.53% and the TAlk error was 0.21%. Furthermore, the model simulation distributions of the DIC and TAlk are similar to the observations (Fig. 1). This includes the simulation of the high value region and the data change of the upper ocean with latitude distribution Correction and calculation The UVic model was used to simulate ocean temperature, alkalinity (ALK), salinity, DIC and phosphate. On the basis of these data, other variables were calculated using the chemistry routine from the OCMIP-3 project ( OCMIP/phase3). These variables included, ph, the concentration of carbonate ions and saturated state of aragonite. We estimated the aragonite saturation state for the shallowwater and deep-sea coral reef locations by assuming that each reef possessed the simulated seawater chemistry of the model grid cell where the reef was located. In this study, we used different assignment methods such as liner interpolation. We did not test the accuracy of each method, because of the uncertainty of the cold-water coral reefs' location and each method's own error. The information for the longitude and latitude of the shallow-water coral reef locations was obtained from the Reef Base Database ( and the information concerning the longitude, latitude, and depth of the cold-water coral locations was obtained from the Global Distribution of Cold-Water Coral Reefs (Freiwald et al., 2004). We used the Global Ocean Data Analysis Project (GLO- DAP) data to correct the model output data. The GLODAP data was obtained in the 1990s and was subsequently gridded into 1 horizontal resolution with 33 vertical layers by Key et al. (2004). Because the observation data is more accurate than the model data, we linearly interpolated the former into the grid of the UVic model. 3. Results 3.1. Comparison of the model and observed results In this section, we compare the result of UVic model and the GLODAP observation data to test the reliability of the ocean simulation model. For example, the observed global Fig. 1. Latitudeedepth distribution of ocean DIC and alkalinity from GLO- DAP observation and UVic model simulation. (a1) DIC of GLODAP data, (a2) DIC of model data in year 1994, (b1) total alkalinity of GLODAP data, (b2) total alkalinity of model data in year To reduce the error between model results and the observation data, we corrected the model data following a formula used in previous studies (Caldeira and Wickett, 2005; Orr et al., 2005; Cao and Caldeira, 2008): D co ðnþ¼d m ðnþ D m1994 þ D ob : Where D co is the corrected data, D m is the model data, D m1994 is the model data of the year 1994, D ob is the observation data and n is the each year from 1800 to Here we assumed a simple linear error. The true error is more complex, so there is still some uncertainty Simulation of marine chemistry The atmospheric CO 2 concentration increased at various rates under the RCP scenarios except RCP2.6 (Fig. 2a), while the ocean continues to absorb CO 2 from the atmosphere under all the scenarios (Fig. 2b). Based on the model simulation results, the global oceans will absorb 162e411 Pg C from the year 2010 to 2100 across the RCP scenarios. By the year 2300, the global oceans will absorb Pg C. Compared with ð3þ

4 192 ZHENG M.-D., CAO L. / Advances in Climate Change Research 5 (2014) 189e196 Fig. 2. Time series under 4 RCP scenarios (a) atmospheric CO2 concentration, (b) CO2 uptake by ocean, (c) global mean ocean surface ph, (d) global mean ocean 2 surface [CO2 3 ], (e) global mean ocean surface U, (f) global mean ocean ph, (g) global mean ocean [CO3 ], (h) global mean U, (i) global mean ocean saturation horizon of aragonite, (j) percent of ocean volume with U > 1.0. the pre-industrial period, the global mean ocean surface temperature will rise 1.1e2.8 K at the year 2100, which is in broad agreement with the prediction of CMIP5 (0.8e3.1 K) (IPCC, 2013). At the year 2300, the ocean temperature will increase 0.8e6.1 K. The increase of CO2 concentration is also important in the strength of the North Atlantic thermohaline circulation. The intensity of the North Atlantic Deep Water (NADW) was 21.5 Sv (1 Sv ¼ 106 m3 s 1) at the year This will decrease to 16.6e19.6 Sv under RCP8.5 to RCP2.6 scenarios. After the CO2 concentration stabilizes, the NADW intensity will increase. Meanwhile, the ocean surface water warming will increase faster than that of the deep ocean so there will be a greater temperature gap between shallow-water and deep-water. Thus, ocean water stratification will be more stable and the convection currents will weaken. With respect to ocean chemistry, according to Equation (1), the dissolved CO2 causes the seawater to produce more hydrogen ions, causing the ph to decrease. At the end of the 21st century, the global mean ocean surface ph will drop by 0.14e0.42 compared with the pre-industrial period (Fig. 2c). The ph will rise 0.05 from the year 2100 to 2300 under the RCP2.6 scenario. Under the other three scenarios, the global mean ph will continue to decrease at different rates until the atmospheric CO2 concentration stabilizes. Under the intensive emission scenario of the RCP8.5 scenario, the global mean ph will drop 0.73 from the year 1800 to 2300, which means that the hydrogen ion concentration will increase by 437%. Between the period 1986e2005 and 2081e2100, the UVic model predicted a decrease in global mean surface ph by units for RCP2.6, for RCP4.5, for RCP6.0, and for RCP8.5. These projected ph changes are in close agreement with the CMIP5 Earth System Model results. For instance, between 1986e2005 and 2081e2100 the CMIP5 models project a model-mean decrease in the surface ph of 0.06e0.07 for RCP2.6, 0.14e0.15 for RCP4.5, 0.20e0.21 for RCP6.0, and 0.30e0.32 for RCP8.5. Due to the increase of the hydrogen ion concentration, the global mean ocean surface carbonate concentration will be reduced by 20%e51% by 2100 under the RCP scenarios. By 2300, the change will be 13%e71% (Fig. 2d). The entirety of the ocean will change more moderately than the ocean surface. The physical and chemical properties of the deep ocean lag behind the surface water (Cao et al., 2014b). Even if the change of the surface water has stabilized following the atmospheric CO2 concentration stabilization, the deep ocean will still experience continued acidification (Fig. 2f). Averaged over the whole ocean, the decrease in ½CO2 3 under the RCP scenarios ranges from 8% to 31%. Fig. 3 shows the latitudeedepth distribution of ocean temperature, DIC, ph and ½CO2 3 in years 2100 and 2300 relative to These changes extend from the surface to the deep sea. Temperature changes at different latitudes are relatively uniform and the area with the greatest warming appears on the surface ocean near 45 N latitude. The zonal average ph, DIC and ½CO2 3 exhibit greater changes in subtropical surface region than in the tropic and high-latitudes. In the vertical direction, these changes are delivered from subtropical area to mid-latitude regions. This is due to the subtropical water sinking and carrying anthropogenic CO2 into the deep ocean by these convection currents (Sabine et al., 2004) Chemical habits of coral reefs Aragonite is a relatively soluble mineral form of calcium carbonate, which is an important component of coral reefs. In the aragonite supersaturated water, corals can easily extract calcium and carbonate ion from the surrounding sea to build the reef. With a volume weighted average, 19% of the ocean was supersaturated in aragonite during the pre-industrial period (Fig. 2j). This will drop to 5%e11% under RCP8.5

5 ZHENG M.-D., CAO L. / Advances in Climate Change Research 5 (2014) 189e196 Fig. 3. Latitudeedepth distribution of ocean variables in 2100 and 2300 differ from 1800 (Dashed lines are latitude mean aragonite saturation horizon in the year 1800; solid lines are latitude mean aragonite saturation horizon in the year 2100 and 2300). to RCP2.6 scenarios. By 2300, over 98% of the ocean will be undersaturated under the RCP8.5 scenario, which could be lethal to coral reefs. The increasing atmosphere CO2 concentration leads to higher seawater temperature and a lower aragonite saturation state, so that corals will suffer the combined effects of global warming and ocean acidification (Pandolfi et al., 2011; Reynaud et al., 2003). The surface seawater will be undersaturated in the Arctic and the south ocean by the mid-21st century under the RCP8.5 scenario. This undersaturation will be delayed 20 years and 60 years under the RCP6.0 and RCP4.5 scenarios and undersaturation will not occur under the high mitigation scenario of RCP2.6. However, a relative observation shows that the Canadian basin surface appears as an aragonite undersaturated area. However, that is due to the 193 limited resolution of the UVic model. Although UVic model has a coupled sea-ice model, the present sea-ice model cannot simulate and predict the future sea-ice variation trend. The real melt rate of sea-ice is faster than the model estimation (Stroeve et al., 2007). Averaged over the whole ocean, the aragonite saturation horizon will increase from a depth of 1138 m in the pre-industrial period to 308 m at the year 2100 under RCP8.5 scenario. In each part of the ocean, it will lift from 1967 me530 m in AtlanticeArctic and 805 me196 m in PacificeIndia. The aragonite saturation states around both shallow- and cold-water coral reefs will rapidly decline. Maps were generated using UV-CDAT ( In the preindustrial period, the water surrounding shallow-water coral reefs had a mean U value of 4.1. Over 99% of these areas are surrounded by seawater with U > 3.5 (Fig. 4). Under the RCP scenarios, the average U of the seawater surrounding the shallow-water coral reefs will decrease to 2.2e3.0. At the year 2055, there will be less than 1% of these areas surrounded by seawater with U > 3.5. This situation will be delayed for 20 years and 35 years under the RCP6.0 and RCP4.5 scenarios, respectively. Even if CO2 emissions follow the RCP2.6 scenario, there will be only 27% of shallow-water coral reefs surrounded by seawater with U > 3.5. With respect to the cold-water coral reefs, the pre-industrial aragonite saturation state of the seawater surrounding deep-sea coral reefs was 1.8 (Fig. 5). By the year 2100, the average saturation state of the water around cold-water coral reefs will decline to 0.9e1.4 (across RCP8.5 to RCP2.6). Meanwhile, there is only 27%e72% of the cold-water coral reefs that will be surrounded by seawater with U > 1. By the year 2300, more than 96% of these will be surrounded by seawater with aragonite undersaturation according to the RCP8.5 scenario. Based on the experiments of shallow-water coral reefs, the calcification rate of reef-building corals would be significantly reduced when aragonite saturation decreased slightly under the situation of supersaturated (Fabry et al., 2008; Langdon et al., 2003). What's more, recent observations show that the calcification rate of corals has been declining in recent decades (De'ath et al., 2013; Su, 2012). Therefore, there are reasons to believe that ocean acidification would have significant influence on coral reefs before the seawater becomes undersaturated in aragonite. 4. Conclusions and discussion (1)The ocean will continuously absorb CO2 from the atmosphere under four RCP scenarios. Thus, the global mean ocean surface temperature will rise by 1.1e2.8 K while the global mean ocean surface ph and [CO2 3 ] will drop 0.14e0.42 and 20%e51% across the RCPs scenarios. The physical and chemical changes of the deep sea lag behind the ocean surface. The effects will last longer than the time under which the atmospheric CO2 concentration is stabilized. It is credible that a large amount of anthropogenic CO2 emissions will continue to impact the whole ocean.

6 Fig. 4. Model-simulated ocean surface aragonite saturation state with shallow-water coral reefs under RCP8.5 scenario, (a1ea4) surface aragonite saturation state overlaid with shallow-water coral reef locations (black dots) at CO2 levels of (around the year 1800), (around the year 2050), (around the year 2080), and (around the year 2100), (b1eb4) percentage distribution of shallow-water coral reefs surrounded by seawater at each aragonite saturation bin. Fig. 5. Model-simulated ocean aragonite saturation horizon with cold-water coral reefs under RCP8.5 scenario, (a1ea4) aragonite saturation horizon overlaid with deep-sea coral reef locations. Coral reef locations at different depth range are represented by different colored dots. Areas in pink represent regions where aragonite saturation horizon has reached the surface. (b1eb4) percentage distribution of cold-water coral reefs surrounded by seawater at each aragonite saturation bin.

7 ZHENG M.-D., CAO L. / Advances in Climate Change Research 5 (2014) 189e (2)The aragonite saturation state of sea water surrounding coral reefs will decrease rapidly. During the pre-industrial period, over 99% of shallow-water coral reefs were surrounded by seawater with U > 3.5 and 87% of cold-water coral reefs were surrounded by seawater with U > 1. On the basis of the model-simulations, by the end of the 21st century, less than 1% of shallow-water coral reefs will be surrounded by seawater with U > 3.5 with the exception of the RCP2.6 scenario, while 73% of cold-water coral reefs will suffer from undersaturated aragonite seawater under the RCP8.5 scenario. The calcification rate of coral reefs will decrease because of the lowering of the aragonite saturation state. Thus, both shallow- and cold-water coral reefs will suffer from ocean acidification. (3)The increasing CO 2 concentration will lead to higher temperatures and lower aragonite saturation state in the ocean. This means that coral reefs will be synergistically impacted by global warming and ocean acidification. Analyzing the situation using four RCP scenarios, reducing CO 2 emissions can effectively slow the process of ocean physical and chemical changes. The predictions of the deep ocean acidification have greater uncertainty in comparison to those of the ocean surface. This is mainly due to the ocean currents simulation and physical transmission process between models (Cao et al., 2009). The model prediction and simulation of deep ocean acidification and its influence on the cold-water coral reefs need to be quantified using the simulations from other earth system models. This study focuses on the open ocean chemistry of the aragonite saturation state of seawater surrounding coral reefs, but the carbonate chemistry within the coral reef system might be substantially different from that of the surrounding seawater (McCulloch et al., 2012; Andersson et al., 2014). Furthermore, changes in a number of other environmental factors, many of which are associated with human activities, such as heat stress, light, salinity, the abundance of food and nutrients, overfishing, and pollution could all influence the fate of coral reefs (Pandolfi et al., 2011). Acknowledgements This work was supported by National Natural Science Foundation of China ( , ), National Key Basic Research Program of China (2015CB953601), Zhejiang University K.P. Chao's High Technology Development Foundation and the Fundamental Research Funds for the Central Universities. References Andersson, A.J., Yeakel, K.L., Bates, N.R., et al., Partial offsets in ocean acidification from changing coral reef biogeochemistry. Nat. Clim. Change 4, 56e61. Caldeira, K., Wickett, M.E., Ocean model predictions of chemistry changes from carbon dioxide emissions to the atmosphere and ocean. J. Geophys. Res. Oceans 110, C09S04. Cao, L., Caldeira, K., Atmospheric CO 2 stabilization and ocean acidification. Geophys. Res. Lett. 35 (19), L Cao, L., Eby, M., Ridgwell, A., et al., The role of ocean transport in the uptake of anthropogenic CO 2. J. Biogeosci. 6, 375e390. Cao, L., Wang, S.-J., Zheng, M.-D., et al., 2014a. Sensitivity of ocean acidification and oxygen to the uncertainty in climate change. J. Environ. Res. Lett. 9 (6), Cao, L., Zhang, H., Zheng, M.-D., et al., 2014b. Response of ocean acidification to a gradual increase and decrease of atmospheric CO 2. Environ. Res. Lett. 9 (2), 239e246. De'ath, G., Fabricius, K., Lough, J., Yes d Coral calcification rates have decreased in the last twenty-five years! Mar. Geol. 346, 400e402. Doney, S.C., Fabry, V.J., Feely, R.A., et al., Ocean acidification: the other CO 2 problem. Ann. Rev. Mar. Sci. 1, 169e192. Fabry, V.J., Seibel, B.A., Feely, R.A., et al., Impacts of ocean acidification on marine fauna and ecosystem processes. Ices J. Mar. Sci. 65, 414e432. Feely, R.A., Byrne, R.H., Acker, J.G., et al., Winter-summer variations of calcite and aragonite saturation in the Northeast Pacific. Mar. Chem. 25 (88), 227e241. Fine, M., Tchernov, D., Scleractinian coral species survive and recover from decalcification. Science 315, Freiwald, A., Fosså, J.H., Grehan, A., et al., Cold-water Coral Reefs. UNEP World Conservation Monitoring Centre, Cambridge UK. Gattuso, J.P., Hansson, L., Ocean Acidification: Background and History. Oxford University Press, Oxford. Guinotte, J.M., Fabry, V.J., Ocean acidification and its potential effects on marine ecosystems. Ann. N. Y. Acad. Sci. 1134, 320e342. Guinotte, J., Orr, J., Cairns, S., et al., Will human-induced changes in seawater chemistry alter the distribution of deep-sea scleractinian corals? Front. Ecol. Environ. 3, 141e146. Indermühle, A., Stocker, T.F., Joos, F., et al., Holocene carbon-cycle dynamics based on CO 2 trapped in ice at Taylor Dome, Antarctica. Nature 398 (6723), 121e126. IPCC, Climate Change 2013: the Physical Science Basis. Cambridge University Press, Cambridge, UK. Key, M.R., Kozyr, A.L., Sabine, C., et al., A global ocean carbon climatology: results from Global Data Analysis Project (GLODAP). Glob. Biogeochem. Cycles 18 (4), 357e370. Langdon, C., Broecker, W.S., Hammond, D.E., et al., Effect of elevated CO 2 on the community metabolism of an experimental coral reef. Glob. Biogeochem. Cycles 17 (1), 11. Matthews, H.D., Cao, L., Caldeira, K., Sensitivity of ocean acidification to geoengineered climate stabilization. Geophys. Res. Lett. 36 (10), L McCulloch, M.J., Falter, J., Trotter, J., et al., Coral resilience to ocean acidification and global warming through ph upregulation. Nat. Clim. Change 2, 623e627. Meissner, K.J., Weaver, A.J., Matthews, H.D., et al., The role of land surface dynamics in glacial inception: a study with the UVic Earth System model. Clim. Dyn. 21, 515e537. Mucci, A., The solubility of calcite and aragonite in seawater at various salinities, temperatures and one atmosphere total pressure. Am. J. Sci. 283 (7), 780e799. Orr, J.C., Najjar, R., Sabine, C., et al., Abiotic-HOWTO. Internal OCMIP Report. LSCE/CEA, Gifsur-Yvette, Saclay, France. Orr, J.C., Fabry, V.J., Aumont, O., et al., Anthropogenic ocean acidification over the twenty first century and its impact on calcifying organisms. Nature 437, 681e686. Pandolfi, J.M., Connolly, S.R., Marshall, D.J., et al., Projecting coral reef futures under global warming and ocean acidification. Science 333 (6041), 418e422. Reynaud, S., Leclercq, N., Romainelioud, S., et al., Interacting effects of CO 2 partial pressure and temperature on photosynthesis and calcification in a scleractinian coral. Glob. Change Biol. 9, 1660e1668. Robbins, L.L., Wynn, J.G., Lisle, J.T., et al., Baseline monitoring of the western Arctic Ocean estimates 20% of Canadian basin surface waters are undersaturated with respect to aragonite. PLoS One 8 (9), 274. Sabine, C.L., Feely, R.A., Gruber, N., et al., The oceanic sink for anthropogenic CO 2. Science 305 (5682), 367e371.

8 196 ZHENG M.-D., CAO L. / Advances in Climate Change Research 5 (2014) 189e196 Schmittner, A., Oschlies, A., Matthews, H.D., et al., Future changes in climate, ocean circulation, ecosystems, and biogeochemical cycling simulated for a business-as-usual CO 2 emission scenario until year 4000 AD. Glob. Biogeochem. Cycles Steinacher, M., Joos, F., Fr olicher, T.L., et al., Imminent ocean acidification in the Arctic projected with the NCAR global coupled carbon cycle-climate model. Biogeosciences 6 (4), 515e533. Stroeve, J., Holland, M.M., Meier, W., et al., Arctic sea ice decline: faster than forecast. Geophys. Res. Lett. 34 (9), 529e536. Su, R.-X., Coral calcification under increasing atmospheric CO 2 concentration and global warming in the southern South China Sea. Quat. Int. 279e280, 474 (in Chinese). Wang, S.-J., Cao, L., Li, N., Responses of the ocean carbon cycle to climate change: results from an earth system climate model simulation. Adv. Clim. Change Res. 5, 123e130. Weaver, A.J., Eby, M., Wiebe, E.C., et al., The UVic Earth system climate model: model description, climatology, and applications to past, present and future climate. Atmos. Ocean 39 (4), 361e428. Xu, Y.-F., Li, Y.-C., Estimates of anthropogenic CO 2 uptake in a global ocean model. Adv. Atmos. Sci. 26 (2), 265e274. Zickfeld, K., Eby, M., Alexander, K., et al., Long-term climate change commitment and reversibility: an EMIC intercomparison. J. Clim. 26 (16), 5782e5809.