Responses of the ocean carbon cycle to climate change:Results from an earth system climate model simulation

WANG Shuang-Jing,CAO Long,LI Na*

Department of Earth Sciences,Zhejiang University,Hangzhou 310027,China

Responses of the ocean carbon cycle to climate change:Results from an earth system climate model simulation

WANG Shuang-Jing,CAO Long,LI Na*

Department of Earth Sciences,Zhejiang University,Hangzhou 310027,China

Based on simulations using the University of Victoria's Earth System Climate Model,we analyzed the responses of the ocean carbon cycle to increasing atmospheric CO2levels and climate change from 1800 to 2500 following the RCP 8.5 scenario and its extension.Compared to simulations without climate change,the simulation with a climate sensitivity of 3.0 K shows that in 2100,due to increased atmospheric CO2concentrations,the simulated sea surface temperature increases by 2.7 K,the intensity of the North Atlantic deep water formation reduces by 4.5 Sv,and the oceanic uptake of anthropogenic CO2decreases by 0.8 Pg C.Climate change is also found to have a large effect on the North Atlantic's ocean column inventory of anthropogenic CO2.Between the years 1800 and 2500,compared with the simulation with no climate change,the simulation with climate change causes a reduction in the total anthropogenic CO2column inventory over the entire ocean and in North Atlantic by 23.1%and 32.0%,respectively.A set of simulations with climate sensitivity variations from 0.5 K to 4.5 K show that with greater climate sensitivity climate change would have a greater effect in reducing the ocean's ability to absorb CO2from the atmosphere.

Climate change;Ocean carbon cycle;Carbon cycle modeling

1.Introduction

According to the Intergovernmental Panel on Climate Change(IPCC)Fifth Assessment Report(IPCC,2013),many currently observed climate changes are unprecedented and the existence of global warming since the 1950s is unequivocal. Due to human activities such as fossil fuel use and land use changes since pre-industrial times,CO2concentrations in the atmosphere have increased by about 40%.Also,only about half of the total anthropogenic CO2emissions remain in the atmosphere((240±10)Pg C);the rest is absorbed by land((160±90)Pg C)and the ocean((155±90)Pg C)(IPCC 2013).Based on simulations using 11 coupled climate-carbon cycle models,Friedlingstein et al.(2006)concluded tha climate change will reduce the efficiency of the land and the ocean in absorbing atmospheric CO2.Furthermore,the IPCC reported that CO2increases and climate change have differen effects on the land's and ocean's carbon storage.CO2increases will lead to the increased carbon storage by the land and ocean Climate change will reduce the land and ocean's capacity to absorb atmospheric CO2due to the increasing temperatures of both the land and ocean and increasing oceanic stratification The ocean,which has absorbed 27.9%of the anthropogenic CO2in the past 200 years(IPCC,2013),plays a crucial role in the global carbon cycle.

The ocean carbon cycle is determined by a series of complex interactions involving the air-sea exchange,inorganic carbon chemistry,ocean general circulation,and marine biological processes(Plattner et al.,2001).Recently,manystudies have discussed the ocean carbon cycle impact of increased atmospheric CO2and climate change.These studies indicate that rising atmospheric CO2levels lead to increased radiative forcing produced by greenhouse gases(IPCC,2013), which causes sea surface warming.Warming seawater will decrease the CO2solubility of seawater(Plattner et al.,2001), resulting in a reduction in oceanic CO2uptake.Furthermore, global warming will cause increased ocean stratification (Sarmiento et al.,1998),which then will cause a reorganization of the thermohaline circulation and the reduction or even collapse of the North Atlantic deep water formation(Manabe and Stouffer,1994;Stocker and Schmittner,1997).These changes will reduce the total transport of oceanic CO2from the ocean surface to the deep ocean(Maier-Reimer et al., 1996;Matebr and Hirst,1999).The reduction in oceanic carbon uptake,caused by the processes mentioned above,will in turn accelerate the increase in atmospheric CO2(Joos et al., 1999).Ocean carbon cycle studies in China began later than other studies worldwide,but have developed rapidly.The early studies mainly used two-dimensional carbon cycle models.For example,Xu et al.(1997)discussed the inorganic carbon cycle,and the sinks and sources of oceanic carbon uptake.Pu and Wang(2000,2001)analyzed the distribution of chemical components related to carbon and the critical factors influencing the distribution of carbon in the Indian Ocean.At present,three-dimensional ocean carbon cycle models are being widely used in climate change research.Xing(2000) and Jin and Shi(2001)pointed out that the biological pump plays an important role in the ocean carbon cycle.Li and Xu (2012)used a perturbation approach to compare the different effects on oceanic CO2uptake in the Pacific Ocean and also modeled its biological processes.In addition to discussing the factors that influence oceanic CO2uptake,other studies(Xu and Li,2009;Bao et al.,2012)have analyzed oceanic CO2uptake and distribution based on simulations.Moreover,Wei et al.(2012,2014)employed an earth system model(BNUESM)to investigate the different climate change and ocean warming responses,including oceanic CO2uptake,carbon sequestration,and ocean acidification,to various anthropogenic CO2emission scenarios.

Based on the previous work noted above,this study uses an earth climate model to simulate the effect of increased atmospheric CO2and associated climate changes on oceanic CO2uptake and the ocean carbon cycle.

2.Model and methods

2.1.Model description

The University of Victoria Earth System Climate Model (UVic)was used,which consists of an energy-moisture balance atmospheric model(Fanning and Weaver,1996),a dynamic-thermodynamic sea-ice model(Bitz et al.,2001; Hibler,1979;Hunke and Dukowicz,1997),and a primitive equation ocean general circulation model(Pacanowski,1995). This model has a horizontal resolution of 3.6°longitude and 1.8°latitude,and divides the ocean into 19 vertical layers.

In order to guarantee the model's computational efficiency, we used a simplified atmospheric model with only one layer, and included only CO2forcing,with no other greenhouse gas forcing such as CH4or aerosol forcing.The vertically integrated thermodynamic energy balance equations assume that the energy and specific humidity decrease vertically with specified scale heights.Momentum conservation equations are replaced by specified wind data,and the atmospheric heat and moisture transport by diffusion are parameterized for simplification(Weaver et al.,2001).The model's wind stress data, 1958-1998 daily reanalysis data(Kalnay et al.,1996),are used to force the ocean and ice components and to calculate the latent heat and sensible heat fluxes between the atmosphere and ocean or ice components.The coupled model's ocean component is Modular Ocean Model(MOM)version 2.2,which is based on Navier Stokes equations that are conditional on Boussinesq and hydrostatic approximations.In the sea-ice component,an elastic-viscous-plastic rheology represents the sea-ice dynamics,and various options for the sea-ice thermodynamics and thickness distribution are included.More information and equations may be found in Weaver et al. (2001).

In addition,our model simulations of the ocean carbon cycle are based not only on the inorganic carbon cycle, following the protocol of the Ocean Carbon Cycle Model Intercomparison Project(Orr et al.,1999),but also take into consideration the ocean ecosystem,including the interactions between nutrients,phytoplankton,zooplankton,and detritus (Schmittner et al.,2008).Model simulations of the terrestrial carbon cycle and vegetation are based on the Met Office-HadleyCentreTRIFFIDdynamicvegetationmodel (Meissner et al.,2003).The model we used in this paper simulates well the large-scale distribution of key climate variables(Weaver et al.,2001)and the ocean carbon cycle (Schmittner et al.,2008),and has been widely used in various global climate and carbon cycle studies(Friedlingstein et al., 2006;Schmittner et al.,2008).

2.2.Simulation experiments

Before conducting the sensitivity simulations,we first set the atmospheric CO2concentration at the pre-industrial level of 280×10-6,and set the mean state of the climate from 1960 to 1990 as the initial condition(Weaver et al.,2001).Then we spun up the model for 10,000 years to reach a quasiequilibrium state between the climate and carbon cycle.At this time,the air-sea carbon flux is close to zero.Then we conducted transient simulations based on specified atmospheric CO2levels.Starting from the end of the spin-up simulation,the atmospheric CO2concentration is based on observations made from 1800 to 2005,and under the RCP 8.5 scenario after 2005.According to the IPCC(2013),the Representative Concentration Pathways(RCPs)usually refer to the pathways extending to 2100 and their corresponding emission scenarios.Extended Concentration Pathways(ECPs) describe the RCPs'extensions from 2100 to 2500.Specifically, RCP 8.5 is a relatively high CO2concentration pathway;theatmospheric CO2concentration reaches to 935.9×10-6by year 2100,then continues to rise over time,and the corresponding ECP maintains a constant concentration after 2250 (Fig.1a).

Based on the speci fied atmospheric CO2concentration,two simulations were performed.One simulation considers climate change due to atmospheric CO2increase and the impact on the ocean carbon cycle due to climate change.The other simulation does not consider climate change due to the atmospheric CO2increase,but assumes that while the change in atmospheric CO2has no direct radiation effect on global climate, the atmospheric CO2increase has a small effect on global warming due to its in fluence on terrestrial vegetation.In the first simulation,we speci fied different climate sensitivity values.The equilibrium climate sensitivity refers to the change in the annual and global mean surface temperature when the climate system reaches a steady state,following a doubling of the atmospheric CO2concentration.Empirical equations are adopted in the model to calculate the radiative forcing by the change in the atmospheric CO2concentration,and therefore radiative forcing can be modi fied by changing the empirical equations'coef ficient,thus changing the climate sensitivity.In our simulations,the climate sensitivity ΔT2xvaries from 0.5 K to 4.5 K at intervals of 0.5 K(ΔT2x=0.5 K,1.0 K,1.5 K, 2.0 K,2.5 K,3.0 K,3.5 K,4.0 K,and 4.5 K).Current estimates of climate sensitivity are 1.5-4.5 K,according to IPCC Fifth Assessment Report(IPCC,2013).In this study,we focus on the difference between the simulation with a climate sensitivity of 3.0 K and the simulation with no climate change with respect to the ocean carbon cycle and the ocean's ability to uptake atmospheric CO2.We also analyze the possible responses of the ocean carbon cycle to different climate sensitivities.

3.Results

Given a specified CO2concentration(RCP 8.5,Fig.1a) the responses of the oceanic CO2uptake to increasing atmospheric CO2concentrations and/or climate change are as follows.If increased atmospheric CO2has no influence on climate change through radiation,the changes in the globa mean land surface air and sea surface temperatures are smal enough to be ignored(Fig.1b and c),and can be regarded as resulting in no climate change(climate sensitivity is 0).When the climate sensitivity is set to 3.0 K,global mean land surface air and sea surface temperatures show a rising trend Even if the CO2concentration remains constant after the year 2250,these temperatures still continue to rise at a relatively slower rate(Fig.1b and c).Relative to the year 1800,the global mean land surface air and sea surface temperatures increase by nearly 8.6 K and 6.6 K by the year 2500 respectively.Research has indicated that,under the same salinity,alkalinity,and dissolved inorganic carbon concentrations,a sea surface warming of 1 K will increase the sea surface CO2partial pressure about 4%(Takahashi et al. 1993).Meanwhile,seasurfacewarmingwillreduce seawaterCO2solubility(Plattneretal.,2001),thus decreasing the oceanic CO2uptake.As shown in Fig.1e-g compared with the simulations with no climate change,decreases in the global mean surface dissolved inorganic carbon(DIC)concentrations,alkalinity,and salinity by year 2500 are 83.4 μmol kg-1(3.7%),57.2 μmol kg-1(2.5%),and 0.2 psu (0.4%)respectively,all of which reach maximum values.In addition,the changes in DIC concentrations and alkalinity have opposite effects on sea surface CO2partial pressure whereby one can nearly offset the other(Plattner et al.,2001). Moreover,Zhou et al.(2005)reported that global warming leads to a rise in sea surface temperature and a dilution in sea surface water at the high latitudes of the North Atlantic Ocean.This will reduce the density of seawater and weaken the meridional density gradient between the high and low latitudes,which will result in the weakening of the thermohaline circulation.This weakening will mainly be limited to the North Atlantic rather than the whole conveyor belt.As shown in Fig.1d,in our simulation,compared with the no climate change case,the intensity of the North Atlantic deep water formation,represented by the maximum meridional overturning stream function,decreases by about 7 Sv at most (1 Sv=106m3s-1,year 2185).The weakening of the North Atlantic deep water formation affects the transport of carbon from the sea surface to deep water,thereby influencing the oceanic CO2uptake(Friedlingstein et al.,2006).Since the rate of oceanic CO2uptake depends on CO2solubility,the differences between atmospheric and sea surface CO2partial pressure,and the ocean dynamic transport of DIC,climate change reduces the oceanic uptake of anthropogenic CO2emissions,and the specific tendencies over time are shown in Fig.1h.In our simulations,the oceanic uptake of anthropogenic CO2(this year's total oceanic CO2uptake minus that in 1800)reaches a maximum of 6.0 Pg C in 2093 with no climate change,while it reaches a maximum of 5.3 Pg C in 2083 with climate change.And the maximum difference of 1.2 Pg C between the two simulations occurs in 2180.

In order to investigate the distribution of anthropogenic CO2in the ocean,we vertically integrated the DIC concentration change(DIC concentration in a given year minus that in 1800)in all 19 layers to produce a global distribution map of the anthropogenic CO2column inventory with and without climate change,as shown in Fig.2.Anthropogenic CO2is not uniformly distributed throughout the ocean.Whether or not there is climate change,the maximum vertically integrated anthropogenic CO2concentration is always found in the North Atlantic Ocean.Due to the formation of North Atlantic deep water,this basin is a major area where seawater sinks(Stocker et al.,1994;Sabine et al.,2004).As it sinks,this seawater can transport the anthropogenic CO2absorbed by the ocean surface to depth more efficiently than in other basins,and thus accumulates a large concentration of anthropogenic CO2. According to our model,the North Atlantic(north to the equator)and the Arctic Ocean account for 10%of the global ocean area and 11.9%of the global ocean volume,but they stored 24.0%(without climate change)and 23.5%(with climate change)of the global oceanic uptake of anthropogenic CO2in 2010,respectively,and 22.0%and 19.5%in 2500, respectively.In 2010,there were no obvious differences in the vertically integrated anthropogenic CO2concentrations between the simulations with and without climate change;the difference in the column inventory of anthropogenic CO2is just 4.1 Pg C(Table 1).However,as atmospheric CO2is absorbed continuously by the ocean and is then transported to the deep ocean,the difference becomes more significant (Fig.2b and c).By 2500,compared to the simulation with no climate change,climate change results in a decrease of 23.1% in the column inventory of anthropogenic CO2in the entire ocean,and of 32.0%in the North Atlantic-Arctic ocean (Table 1).

Next,we divided the ocean into the Atlantic-Arctic and the Pacific-Indian basins and analyzed the latitude-depth distribution of anthropogenic CO2to gain a better understanding of the distribution and transport of anthropogenic CO2.

Anthropogenic CO2enters the ocean via the air-sea exchange and the contours of the highest anthropogenic CO2concentrations appear in basins shallower than 500 m under each scenario(Fig.3).This means that anthropogenic CO2concentrates in the near surface water,and the thermocline is a crucial area for the oceanic uptake of anthropogenic CO2. As shown in Fig.3,the highest concentration of anthropogenic CO2is found in the subtropics where near surface seawater shows convergence and subsidence,thus allowing DIC to be transported to the deep ocean,and benefitting further oceanic uptake of anthropogenic CO2(Bao et al., 2012).

The penetration depth of anthropogenic CO2in the ocean depends on the rate that the near sea surface anthropogenic CO2is transported to the interior(Sabine et al.,2004).And the transport rate mainly depends on ocean circulation factors, such as thermocline ventilation,deep water and intermediate water formation,and so on(Friedlingstein et al.,2006).As it is affected by these factors,the penetration of anthropogenic CO2absorbed by the Pacific-Indian ocean is shallower and slower,and the penetration of anthropogenic CO2absorbed by the basins with deep water formation,such as the North Atlantic Ocean,is deeper and faster(Fig.3).One study reports that regardless of whether basins are characterized by deep water formation or not,the time scales of the near surface water being mixed to the deep ocean may be centuries(Sabine et al.,2004).Climate change increases sea surface temperature (Fig.1b),which causes a decrease in the oceanic CO2uptake by reducing the seawater's CO2solubility(Plattner et al., 2001).This also leads to increased stratification(Sarmiento et al.,1998).The increased stratification leads to a decrease in vertical mixing along the isopycnals and overturning circulation,which decreases the carbon vertical transport.The increased stratification also causes a gradual collapse of thermohaline circulation throughout the entire deep oceanthe thermohaline circulation plays the major role in the carbon balance on a century time scale(Sarmiento et al.,1998),and its collapse has been verified in related early simulations(Joos et al.,1999).As Fig.3 shows,in the year 2500 the 240 μmol kg-1contour of the DIC concentration can extend to a depth of 4,000 m in the Atlantic-Arctic ocean without climatechange(Fig.3a),butinthesimulationwith ΔT2x=3.0 K,it is confined to the sea surface in most of the world ocean(Fig.3b).

To further analyze the effects of climate change on the ocean carbon cycle,we calculated changes in the global oceanic anthropogenic CO2uptake for the years 2100,2500, and 1800-2100,as well as the 1800-2500 cumulative uptake for different climate sensitivities.The climate sensitivities range from 0.5 K to 4.5 K at intervals of 0.5 K,and with no climate change the climate sensitivity is 0,as shown in Fig.4. In 2100,the oceanic uptake of anthropogenic CO2shows a linear change along with the climate sensitivity.When climate sensitivity increases by 0.5 K,the annual uptake decreases by about 0.1 Pg C,and cumulative uptake decreases by about 6.5 Pg C.In 2500,the change does not display a linear trend. Without climate change,the oceanic anthropogenic CO2uptake is 1.5 Pg C and the cumulative uptake is 1,996.0 Pg C. When the climate sensitivity reaches 4.5 K,the annual uptake is 0.9 Pg C and the cumulative uptake is 1,502.8 Pg C,representing decreases of 39.9%and 24.7%,respectively.As shown in Table 2,both the global mean surface air warming and global mean ocean warming values become larger as climate sensitivity increases.Regarding the global mean ocean warming,simulated under different climate sensitivities(horizontal axis),the oceanic carbon uptake and cumulative uptake of anthropogenic CO2show a similar trend in the Figures whose horizontal axis is climate sensitivity.This series of simulations illustrate that climate warming decreases oceanic carbon uptake,and the higher the climate sensitivity or the greater the ocean warming,the larger is the reduction in the oceanic uptake of anthropogenic CO2.

4.Conclusions and discussion

(1)For the same CO2concentrations,relative to simulations without climate change(climate sensitivity is 0),when the climate sensitivity is set as 3.0 K,the changes in global mean surface air temperature and global mean sea surface temperature become larger over time,showing increases of 8.2 K and 6.4 K,respectively,by the year 2500.This warming has direct impacts on the sea surface partial pressure and the solubility of CO2.Owing to climate change,the sea surface DIC concentration, alkalinity,and salinity also decrease by 3.7%,2.5%,and 0.4%,respectively.Moreover,the North Atlantic thermohaline circulation weakens as a result of ocean warming and freshening,and the intensity of the North Atlantic deep water formation decreases by 7 Sv by the year 2185.Due to the net effect of climate change,the oceanic uptake of anthropogenic CO2decreases by 1.2 Pg C by the year 2180.

(2)Anthropogenic CO2is not evenly distributed in the ocean,and the maximum vertically integrated values are found in the North Atlantic basin because of the characteristic North Atlantic deep water formation.The largest vertically integrated differences induced by climate change(ΔT2x=3.0 K)are also found in the North Atlantic basin.By the year 2500,the total anthropogenic CO2in the entire ocean and in the North Atlantic-Arctic basins decrease by 23.1%and 32.0%,respectively,as a result of climate change,and the reduction in the North Atlantic-Arctic basin accounts for 30.4%of the total reduction,despite the fact that the area accounts for just 16.1%of the total ocean area,the volume accounts for just 11.9%of the total ocean volume.

(3)Anthropogenic CO2is also unevenly distributed with increasing depth in the ocean.Anthropogenic CO2mainly concentrates in the thermocline and its highest concentrations are found in subtropical surface waters. In the simulation with a climate sensitivity of 3.0 K,the concentrations of anthropogenic CO2in subtropical surface water decreased by about 20%in 2500 relative to the case without climate change.Compared with the Pacific and Indian oceans,the transport rate is faster and the concentration of anthropogenic CO2is higher in deep water formation regions such as the North Atlantic basin.

(4)Different climate sensitivities have different effects on the ocean carbon cycle.The greater the climate sensitivity,the less the annual oceanic and cumulative uptakes.

In this study,we chose a specified atmospheric CO2concentration to simulate resultant changes in the oceanic CO2uptake and the ocean carbon cycle affected by climate change The ocean is an important carbon sink.In reality,if the CO2uptake by the ocean decreases,the anthropogenic CO2accumulating in the atmosphere will increase,which will in turn affect climate change.Therefore,it is necessary to study the impact of climate change on the ocean carbon cycle and the resulting feedback of the ocean carbon cycle to climate change.In future studies,we will use a specified CO2emission to investigate the feedback interactions between climate change,atmospheric CO2concentrations,and the ocean carbon cycle.

Acknowledgements

This work is supported by National Natural Science FoundationofChina(41276073)andtheFundamenta Research Funds for the Central Universities(2102XZZX012)

Bao,Y.,Qiao,F.-L.,Song,Z.-Y.,2012.The 3-dimensional numerical simulation of global ocean carbon cycle.Acta Oceanol.Sin.34(3),19-26(in Chinese).

Bitz,C.M.,Holland,M.M.,Weaver,A.J.,et al.,2001.Simulating the icethickness distribution in a coupled climate model.J.Geophys.Res. Oceans 106(C2),2441-2463.

Fanning,A.F.,Weaver,A.J.,1996.An atmospheric energy-moisture balance model:climatology,interpentadal climate change,and coupling to an ocean general circulation model.J.Geophys.Res.Atmospheres 101(D10), 15111-15128.

Friedlingstein,P.,Cox,P.,Betts,R.,et al.,2006.Climate-carbon cycle feedback analysis:results from the C4MIP model intercomparison.J.Clim.19 (14),3337-3353.

Hibler III,W.D.,1979.A dynamic thermodynamic sea ice model.J.Phys. Oceanography 9(4),815-846.

Hunke,E.C.,Dukowicz,J.K.,1997.An elastic-viscous-plastic model for sea ice dynamics.J.Phys.Oceanography 27(9),1849-1867.

IPCC,2013.Climate Change 2013:The Physical Science Basis.Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change.Cambridge University Press,Cambridge. Accessed.http://www.climatechange2013.org/images/report/WG1AR5_ ALL_FINAL.pdf.

Jin,X.,Shi,G.-Y.,2001.The role of biological pump in ocean carbon cycle. Chin.J.Atmos.Sci.25(5),683-688(in Chinese).

Joos,F.,Plattner,G.K.,Stocker,T.F.,et al.,1999.Global warming and marine carbon cycle feedbacks on future atmospheric CO2.Science 284(5413), 464-467.

Kalnay,E.,Kanamitsu,M.,Kistler,R.,et al.,1996.The NCEP/NCAR 40-year reanalysis project.Bull.Am.Meteorol.Soc.77(3),437-471.

Li,Y.-C.,Xu,Y.-F.,2012.Uptake and storage of anthropogenic CO2in the Pacific Ocean estimated using two modeling approaches.Adv.Atmos.Sci. 29(4),795-809.

Maier-Reimer,E.,Mikolajewicz,U.,Winguth,A.,1996.Future ocean uptake of CO2:interaction between ocean circulation and biology.Clim.Dyn.12 (10),711-722.

Manabe,S.,Stouffer,R.J.,1994.Multiple-century response of a coupled ocean-atmosphere model to an increase of atmospheric carbon dioxide.J. Clim.7(1),5-23.

Matebr,R.J.,Hirst,A.C.,1999.Climate change feedback on the future oceanic CO2uptake.Tellus B 51(3),722-733.

Meissner,K.J.,Weaver,A.J.,Matthews,H.-D.,et al.,2003.The role of land surface dynamics in glacial inception:a study with the UVic Earth System Model.Clim.Dyn.21(7-8),515-537.

Orr,J.C.,Najjar,R.,Sabine,C.L.,et al.,1999.Abiotic-howto,Internal OCMIP Report.LSCE/CEA Saclay,p.25.

Pacanowski,R.C.,1995.MOM 2 documentation,user's guide and reference manual.GFDL Ocean Group Tech.Rep.3(3),232.

Plattner,G.K.,Joos,F.,Stocker,T.F.,et al.,2001.Feedback mechanisms and sensitivities of ocean carbon uptake under global warming.Tellus B 53(5), 564-592.

Pu,Y.-F.,Wang,M.-X.,2000.An ocean carbon cycle model part I:establishing of carbon model including an oceanic dynamic general circulation field,chemical,physical and biological processes occurred in the ocean. Climatic Environ.Res.5(2),129-140(in Chinese).

Pu,Y.-F.,Wang,M.-X.,2001.An ocean carbon cycle model part II:simulation analysis on the Indian Ocean.Climatic Environ.Res.6(1),67-76(in Chinese).

Sabine,C.L.,Feely,R.A.,Gruber,N.,et al.,2004.The oceanic sink for anthropogenic CO2.Science 305(5682),367-371.

Sarmiento,J.L.,Hughes,T.M.C.,Stouffer,R.J.,et al.,1998.Simulated response of the ocean carbon cycle to anthropogenic climate warming. Nature 393(6682),245-249.

Schmittner,A.,Oschlies,A.,Matthews,H.D.,et al.,2008.Future changes in climate,ocean circulation,ecosystems,and biogeochemical cycling simulated for a business-as-usual CO2emission scenario until year 4000 AD.Glob.Biogeochem.Cycles 22,1013.

Stocker,T.F.,Schmittner,A.,1997.In fluence of CO2emission rates on the stability of the thermohaline circulation.Nature 388(6645), 862-865.

Stocker,T.F.,Broecker,W.S.,Wright,D.G.,1994.Carbon uptake experiments with a zonally-averaged global ocean circulation model.Tellus B 46(2), 103-122.

Takahashi,T.,Olafsson,J.,Goddard,J.G.,et al.,1993.Seasonal variation of CO2and nutrients in the high-latitude surface oceans:a comparative study. Glob.Biogeochem.Cycles 7(4),843-878.

Weaver,A.J.,Eby,M.,Wiebe,E.C.,et al.,2001.The UVic Earth System Climate Model:model description,climatology,and applications to past, present and future climates.Atmos.Ocean 39(4),361-428.

Wei,T.,Yang,S.,Moore,J.C.,et al.,2012.Developed and developing world responsibilities for historical climate change and CO2mitigation.Proc. Natl.Acad.Sci.109(32),12911-12915.

Wei,T.,Dong,W.,Yuan,W.,et al.,2014.In fluence of the carbon cycle on the attribution of responsibility for climate change.Chin.Sci.Bull.59(19), 2356-2362.

Xing,R.-N.,2000.A three-dimensional world ocean carbon cycle model with ocean biota.Chin.J.Atmos.Sci.24(3),333-340(in Chinese).

Xu,Y.-F.,Li,Y.-C.,2009.Estimates of anthropogenic CO2uptake in a global ocean model.Adv.Atmos.Sci.26(2),265-274.

Xu,Y.-F.,Wang,M.-X.,Jin,X.-Z.,1997.A two-dimensional ocean thermohaline circulation carbon cycle model.Sci.Atmos.Sin.21(5),573-580 (in Chinese).

Zhou,T.-R.,Yu,X.,Liu,Y.,et al.,2005.Weak response of the Atlantic thermohaline circulation to an increase of atmospheric carbon dioxide in IAP/LASG Climate System Model.Chin.Sci.Bull.50(3),269-275(in Chinese).

Received 26 May 2014;revised 24 July 2014;accepted 1 August 2014

Available online 30 December 2014

*Corresponding author.

E-mail address:lina2012@zju.edu.cn(LI N.).

Peer review under responsibility of National Climate Center(China Meteorological Administration).

http://dx.doi.org/10.1016/j.accre.2014.11.004

1674-9278/Copyright?2014,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(http://creativecommons.org/licenses/by-nc-nd/3.0/).