A Brief Introduction to BNU-HESM1.0 and Its Earth Surface Temperature Simulations

YANG Shili,DONG Wenjie?,CHOU Jieming,FENG Jinming,YAN Xiaodong, WEI Zhigang,YUAN Wenping,GUO Yan,TANG Yanli,and HU Jiacong

1State Key Laboratory of Earth Surface Processes and Resource Ecology,Beijing Normal University,Beijing 100875

2Zhuhai Joint Innovative Center for Climate-Environment-Ecosystem,Future Earth Research Institute, Beijing Normal University,Zhuhai 519087

3Key Laboratory of Regional Climate–Environment for East Asia,Institute of Atmospheric Physics, Chinese Academy of Sciences,Beijing 100029

A Brief Introduction to BNU-HESM1.0 and Its Earth Surface Temperature Simulations

YANG Shili1,2,DONG Wenjie?1,2,CHOU Jieming1,2,FENG Jinming3,YAN Xiaodong1,2, WEI Zhigang1,2,YUAN Wenping1,2,GUO Yan1,2,TANG Yanli1,and HU Jiacong2

1State Key Laboratory of Earth Surface Processes and Resource Ecology,Beijing Normal University,Beijing 100875

2Zhuhai Joint Innovative Center for Climate-Environment-Ecosystem,Future Earth Research Institute, Beijing Normal University,Zhuhai 519087

3Key Laboratory of Regional Climate–Environment for East Asia,Institute of Atmospheric Physics, Chinese Academy of Sciences,Beijing 100029

Integrated assessment models and coupled earth system models both have their limitations in understanding the interactions between human activity and the physical earth system.In this paper,a new human–earth system model,BNUHESM1.0,constructed by combining the economic and climate damage components of the Dynamic Integrated Model of Climate Change and Economy to the BNU-ESM model,is introduced.The ability of BNU-HESM1.0 in simulating the global CO2concentration and surface temperature is also evaluated.We f i nd that,compared to observation,BNU-HESM1.0 underestimates the global CO2concentration and its rising trend during 1965–2005,due to the uncertainty in the economic components.However,the surface temperature simulated by BNU-HESM1.0 is much closer to observation,resulting from the overestimates of surface temperature by the original BNU-ESM model.The uncertainty of BNU-ESM falls within the range of present earth system uncertainty,so it is the economic and climate damage component of BNU-HESM1.0 that needs to be improved through further study.However,the main purpose of this paper is to introduce a new approach to investigate the complex relationship between human activity and the earth system.It is hoped that it will inspire further ideas that prove valuable in guiding human activities appropriate for a sustainable future climate.

economic model component,earth system model,human activity,global change

1.Introduction

Integrated assessment models(IAMs)and coupled atmosphere–land–sea-ice earth system models are two major tools used to investigate issues relates to climate change. However,both IAMs and coupled earth system models have their advantages and limitations in simulating and projecting emissions scenarios and climate changes.

IAMs are widely used to provide information on global and national emissions and the costs of different climate polices.They can depict the cause–effect chain between emissions and rising temperature,but with relatively simple equations(Matsuoka et al.,1995;Brenkert et al.,2003;Sokolov et al.,2005;Bouwman et al.,2006;Moss et al.,2010).For example,the Dynamic Integrated Model of Climate Change and Economy(DICE)uses just one or two simple equations to calculate the global CO2concentration and temperature change,and most IAMs use a linearized representation of ocean carbon uptake(van Vuuren et al.,2011;Hof et al., 2012).This is an obvious oversimplif i cation of the complex physical and biochemical processes of the earth system.Indeed,recent studies have shown that the lack of complexity in the representation of physical processes in IAMs is likely to have signifi cant impacts on the outputs of policy cost,carbon tax,and so on(Schneider and Thompson,1981;Schultz and Kasting,1997;Smith and Edmonds,2006).For example, Smith and Edmonds(2006)demonstrated that the uncertainty in the carbon cycle results in a broader range of the cost of achieving certain CO2concentrations,and is equivalent to a change in the concentration target of up to 100 ppmv.

Coupled earth systems models,which include ocean and terrestrial carbon cycle feedbacks,are designed to capture the biophysical processes of the real climate system(Taylor et al.,2012).They possess unique advantages in simulating, projecting and attributing climate change under given naturaland anthropogenic forces(Wei et al.,2012;Gillett and Fyfe, 2013;Weller and Cai,2013;Bellenger et al.,2014;Wei et al.,2014)and can be used to compensate for the def i ciencies of IAMs.For example,Wei et al.(2012)investigated the responsibilities of developed and developing countries in terms of historical climate change and CO2mitigation using two earth system models,and the results showed contributions to historical climate change of 2/3 and 1/3,respectively, while their contribution to future climate mitigation was predicted to be 1/3 and 2/3.Additionally,coupled earth system models provide larger datasets for the Intergovernmental Panel on Climate Change(IPCC)reports,which deepen our understanding of climate change.However,these models still do not include economic and climate-damage modules; they can only use the carbon emissions provided by separate economic models as external forces,and are thus unable to ref l ect the interactions between human activity and the physical earth system.So,to achieve more realistic simulations,it is essential and necessary to combine the economic and climate-damage components of IAMs with coupled earth system models,which should be possible in this era of increasing supercomputing power.A successful precedent is the integrated global system model framework(IGSM)of the Massachusetts Institute of Technology.The IGSM includes an emissions prediction and policy analysis model(EPPA) and uses an earth system model with intermediate complexity as its earth system component.It has been widely used in addressing scientif i c goals of earth system modeling and helping to inform the policy-making process(Prinn,2012; Reilly et al.,2013).However,in China,very little progress has been made in this area of study.

In this paper,a new coupled human–earth system model (BNU-HESM1.0)is introduced.Designed to simulate the interaction between human activity and the earth system,the model was developed by integrating a simplified version of DICE and the coupled earth system model from Beijing Normal University(BNU-ESM).The construction of BNUHESM1.0 is described in section 2.In section 3,the model’s capability in simulating the global CO2concentration and surface temperature is evaluated.A summary and further discussion are provided in section 4.

2.Methods

2.1.DICE

DICE,which was first developed by Nordhause(1992), an economist from Yale University,formed the basis of IAMs in climate change.DICE was based on the Ramey model of optimal economic growth and was designed to maximize the discounted“utility”under a number of economic and climate constraints.The model consists of three parts:(1)an economic–emissions relationship;(2)an objective function; and(3)an emissions–climate relationship.In the following, we provide the main equations for the first and second parts, but for a more detailed description of DICE readers are referred to Nordhaus(1992).

In the first part,the output is given by the Cobb–Douglas function,Eq.(1),whereQ(t)represents the net outputs andγis the elasticity of output with respect to the capital.The termsK(t)andL(t)represent the capital and technology,respectively.The term ?(t)represents the net output ratio, which is the fraction of the rest of the output that the impact of climate damage and policy interference on the total output is removed—see Eq.(2).In Eq.(2),the termμ(t) represents the fraction of the reduction of the emission relative to uncontrolled emission,andT(t)is the global averaged temperature.The empirical coeff icientsb1,θ1andθ2in Eq. (2)are 0.0686,2.887,0.00144 and 2.0,respectively.

DICE also includes two equations that describe the distribution of net output and the process of capital accumulation. The net output is divided into investment(I)and consumption(C),Eq.(3),and the termδtin Eq.(4)is the depreciation rate.

Equation(5)describes the relationship between the net outputs and emissions,whereσ(t)is the exogenous technology parameter,which describes the trend of the output and emissions ratio.

DICE maximizes a social welfare function that is the discounted sum of the utility of per capital consumption.The mathematical expression is given as follows,

wherec(t)represents the per capital consumption andρis the pure rate of social time preference.The utility function is given in Eq.(6),in whichαrepresents the rate of inequality aversion.

In this paper,we assume theμ(t)=0 and use Eqs.(1–5) to calculate the global emissions amount by prescribing the exogenous variablesA(t),K(t),L(t)andσ(t).This simplified approach is designed to match the integration algorithm of BNU-ESM,and is inconsistent with the original algorithms that calculate the optimal emissions path by maximizing the discounted utility.

2.2.BNU-ESM

The earth system model of Beijing Normal University (BNU-ESM)is an atmosphere–land–sea-ice fully coupled model,and was one of the models that participated in the Coupled Model Intercomparison Project Phase 5(CMIP5). Recently,various aspects of the simulation ability of BNUESM have been evaluated,and it is has also been widely usedin climate attribution and projection studies(Wuebbles et al., 2013;Bellengeretal.,2014;Mehran etal.,2014).Brief l y,the atmosphere component is CAM3.5(Community Atmosphere Model Version 3.5)from the NCAR(National Center for Atmospheric Research),the land component is CoLM(Common Land Model),which includes the LPJ(Lund–Postdam–Jena)DGVM(dynamic global vegetation model),and its ocean component is the idealized biogeochemical cycling (iBGC)module.Wu et al.(2013)provides a detailed description of BNU-ESM.Comparison between BNU-ESM results, observation,and other Chinese coupled climate model outputs indicates reasonable performance(Zhou et al.,2014).

2.3.BNU-HESM1.0 model construction

BNU-HESM1.0 was built by coupling the economic parts of DICE with BNU-ESM.The coupled model includes the following two steps:First,the economic model calculates global emissions using exogenous variables(e.g.,capital, population,technology etc.),and then BNU-ESM uses these emissions as external forcing to calculate the possible climate changes.The second step connects the climate change results and the economic components through the damage function, describing the climate impact on the economic process(e.g., the GDP or outputs).These two processes are similar to the integrated approach of general climate–economic models.

In the firststep,BNU-ESMusually usesthe monthly gridded CO2fossil fuel flux(units:kg m?2s?1)as the external forcing,but DICE only outputs the annual global emission (units:Gt C).So,we have to convert the annual global CO2emissions amount to CO2flux,which has a certain spatial and temporal resolution,and we realize this goal by def i ning theMindex as follows:

whereEmandEannrepresent the observed monthly gridded CO2flux and annual emissions amount,respectively.The termsrepresent the monthly gridded CO2flux and annual emissions amount calculated by DICE,and we can thus obtainas follows:

In the second step,the global annual averaged temperatureT(t)is calculated using parallel computing methods,and is called by the damage function,Eq.(2),in the economic components.Besides,it should be noted that we modif i ed the original time step of DICE(10 yr)to 1 yr,to coordinate it with BNU-ESM.As the exogenous variable of DICE is expressed in exponential form,this modif i cation is reasonable and has no effect on the outputs(Figs.1 and 2).The economic components calculate the emissions each year,and the climate damage also feeds back to the economic components each year.Therefore,BNU-HESM1.0 can describe the complete cycle of the socioeconomic–earth system.

3.Experimental design

We conducted two simulations to evaluate the simulation ability of BNU-HESM1.0.One was the historical simulation carried out by the original model(BNU-ESM),referred to as ESM-Hist,which followed the CMIP5 historical simulation requirements(Taylor et al.,2012)and whose simulation period ran from 1850 to 2005.The other simulation,referred to as HESM-Hist,was carried out by BNU-HESM1.0.In terms ofthe inputdata ofthe economic components,which began in 1965,HESM-Hist was simulated from 1965 to 2005 and used the simulation results of ESM-Hist in 1965 as the initial condition.In the HESM-Hist simulation,BNU-HESM1.0 calculated the CO2flux with its economic component,so there was no need to provide the observed CO2flux as external forcing. Apart from the CO2flux,the other forcing conditions were the same as those used in ESM-Hist during the period 1965–2005.Further details regarding the natural and anthropogenic forcing conditions are provided in Table 1.

4.Results

4.1.C O2emissions and concentration

Figure 3 shows the observed and simulated global CO2concentrations from BNU-ESM and BNU-HESM1.0.It isclear that BNU-ESM overestimates the CO2concentration and its rising trend after the 1970s,while BNU-HESM1.0 produces a much lower CO2concentration during 1965–2005,resulting from its underestimations of global annual emissions(Fig.4).Figure 4 also indicates that BNUHESM1.0 can reproduce the increasing trend of global emissions,but it does not have the ability to capture the interannual variability of global emissions.We consider this inability mainly result from the uncertainty of the economic exogenous variables,such as the population,capital,and so on.For 2005,the CO2concentration simulated by BNU-HESM1.0 model is about 25 ppm less than that observed,which is much greater than the BNU-ESM uncertainty.

Table 1.Forcing data for the experiments.

4.2.Global temperature

Figure 5 shows the observed and simulated global temperature anomalies relative to 1965–99.It indicates that BNU-ESM underestimates the global temperature during 1850–1970,but overestimates it after the 1970s.In 2005, the temperature simulated by BNU-ESM is 0.4?C above the observation.This deviation is within the range of CMIP5 model uncertainty(?0.1?C to 0.7?C).Relatively,the temperature simulated by BNU-HESM1.0 during 1965–2005 is much closer to that observed.The lineartrend ofthe temperature over1965–2005 is0.16?C(10 yr)?1,0.3?C(10 yr)?1and 0.14?C(10 yr)?1for the observation,BNU-ESM and BNUHESM1.0,respectively.However,this result does not mean that BNU-HESM1.0 is better than BNU-ESM at simulating the temperature,because BNU-HESM1.0 underestimates the global emissions amount.Both BNU-ESM and BNUHESM1.0 can reproduce the globaltemperature gradientwith latitude,with little difference between the two models(Fig. 6).This means that BNU-HESM1.0 can maintain BNUEMS’ability in simulating the spatial distribution of climate variables.

4.3.Conclusion and discussion

This paper describes the first attempt to construct a human–earth system model(BNU-HESM1.0)by coupling an economic model to an earth system model,and provides a new approach to investigating the complex interactions between human activity and the earth system.BNU-HESM1.0 is more advanced compared to DICE or BNU-ESM.Specifically,compared to DICE,BNU-HESM1.0 includes a fully coupled earth system model,which is superior in simulating various aspects of climate change.Meanwhile,compared to BNU-ESM,it not only retains the simulation ability of BNU-ESM,but is also able to quantify emissions levels and climate damage by itself,which represents an advancement over BNU-ESM.Thus,we could use BNU-HESM1.0 to project future emissions levels and climate change under given macroeconomic conditions(e.g.,population,capital etc.).However,there are still many unresolved problems and uncertainties with BNU-HESM1.0.For example,the results of the present reported experiments showed it underestimates global CO2emissions and cannot reproduce its interannual variability,due to the uncertainty of the empirical parameters and exogenous variables.As we know,economics comprises complex issues affected by numerous social activities (e.g.,policy,employment etc.).These complex processes are highly parameterized or even ignored in DICE,resulting in its limitations in capturing a complete socioeconomic picture. Besides,economic processes also possess large regional disparities,and these features are not well expressed in DICE. Therefore,in future work,we should try to use,instead of DICE,much more complex economic models,such as RICE (Nordhaus,2010)or EPPA(Paltsev et al.,2005),which feature more detailed economic processes and can re flect the regional disparity of the economy.Alternatively,we could seek to add new parameters to make the economic model much more re flective of reality.In terms of the climate damage component,this is excessively simple,possessing many de fi ciencies.For example,as well as economic output,climate change also has signi ficant impacts on agriculture,ecology,health,and so on,which is not well expressed in this study.So,other methods to develop BNU-HESM1.0 must include attempts to improve the climate-damage function.A fi nal but very important point is that the ability of BNUHESM1.0 in simulating various aspects of the natural system is largely dependent on that of BNU-ESM.So,we also need to improve BNU-ESM using traditional development methods,such as enhancing the resolution,adding physical processes with smaller scales,and so on.

Despite the limitations acknowledged in the above discussion,this paper represents a pioneering attempt to build a human–earth system model,and several ways to improve or complete BNU-HESM1.0 have been identif i ed.We hope that the development of BNU-HESM1.0 in the future produces a comprehensive tool that not only provide guidance for appropriate human activity(Fu and Ye,1995;Ye et al.,2001,2003, 2009)and global and regional sustainable development,but also strengthens the communications and exchanges of the IPCC’s three working groups.

Acknowledgements.This work was funded by the National Key Program for Global Change Research of China(Grant No. 2010CB950504),the National Basic Research Program of China (973 Program)(Grant No.2012CB95570001),and the nationallevel major cultivation project of Guangdong Province entitled:“The construction and application of coupled earth system and social economic models”(Grant No.2014GKXM058).

REFERENCES

Bellenger,H.,E.Guilyardi,J.Leloup,M.Lengaigne,and J. Vialard,2014:ENSO representation in climate models:FromCMIP3 to CMIP5.Climate Dyn.,42(7–8),1999–2018.

Bouwman,A.,T.Kram,and K.K.Goldewijk,2006:Integrated Modelling of Global Environmental Change,an Overview of IMAGE 2.4.The Netherlands Environmental Assessment Agency,Bilthoven.228 pp.

Brenkert,A.L.,A.Smith,S.H.Kim,and H.M.Pitcher,2003:Model Documentation for the MiniCAM.Pacific Northwest National Laboratory,Richland,WA.188 pp.

Fu,C.B.,and D.Z.Ye,1995:Global Change and the Future Trend of Ecological Environment Evolution in China.Scientia Atmospherica Sinica,19(1),116–126.(in Chinese)

Gillett,N.P.,and J.C.Fyfe,2013:Annular mode changes in the CMIP5 simulations.Geophys.Res.Lett.,40(6),1189–1193.

Hof,A.F.,C.W.Hope,J.Lowe,M.D.Mastrandrea,M.Meinshausen,and D.P.van Vuuren,2012:The benefits of climate change mitigation in integrated assessment models:The role of the carbon cycle and climate component.Climatic Change, 113(3-4),897–917.

Matsuoka,Y.,M.Kainuma,and T.Morita,1995:Scenario analysis of global warming using the Asian Pacific Integrated Model (AIM).Energy Policy,23(4-5),357–371.

Mehran,A.,A.AghaKouchak,and T.J.Phillips,2014:Evaluation of CMIP5 continental precipitation simulations relative to satellite-based gauge-adjusted observations.J.Geophys. Res.,119(4),1695–1707.

Moss,R.H.,and Coauthors,2010:The next generation of scenarios for climate change research and assessment.Nature, 463(7282),747–756.

Nordhaus,W.D.,1992:The‘DICE’Model:Background and structure of a dynamic integrated climate-economy model of the economics of global warming.Cowles Foundation for Research in Economics,Yale University.[Available online from http://cowles.econ.yale.edu/P/cd/d10a/d1009.pdf]

Nordhaus,W.D.,2010:Economic aspects of global warming in a post-Copenhagen environment.Proceedings of the National Academy of Sciences of the United States of America, 107(26),11 721–11 726.

Paltsev,S.,J.M.Reilly,H.D.Jacoby,R.S.Eckaus,J.R.Mc-Farland,M.C.Sarof i m,M.O.Asadoorian,and M.H.M. Babiker,2005:The MIT Emissions Prediction and Policy Analysis(EPPA)Model:Version 4.MIT Joint Program on the Science and Policy of Global Change.[Available online from http://hdl.handle.net/1721.1/29790]

Prinn,R.G.,2012:Development and application of earth system models.Proceedings of the National Academy of Sciences of the United States of America,110,3673–3680.

Reilly,J.,and Coauthors,2013:Valuing climate impacts in integrated assessmentmodels:The MIT IGSM.Climatic Change, 117(3),561–573.

Schneider,S.H.,and S.L.Thompson,1981:Atmospheric CO2and climate:Importance of the transient response.J.Geophys.Res.,86,3135–3147.

Schultz,P.A.,and J.F.Kasting,1997:Optimal reductions in CO2emissions.Energy Policy,25,491–500.

Smith,S.J.,and J.A.Edmonds,2006:The economic implications of carbon cycle uncertainty.Tellus B,58,586–590.

Sokolov,A.P.,and Coauthors,2005:Mit Integrated Global System Model(IGSM)Version 2:Model description and baseline evaluation.MIT Joint Program on the Science and Policy of Global Change,Report No.124.[Available online at http://hdl.handle.net/1721.1/29789.]

Taylor,K.E.,R.J.Stouffer,and G.A.Meehl,2012:An overview of CMIP5 and the experiment design.Bull.Amer.Meteor. Soc.,93(4),485–498.

van Vuuren,D.P.,and Coauthors,2011:How well do integrated assessment models simulate climate change?.Climatic Change,104(2),255–285.

Wu,Q.Z.,J.M.Feng,W.J.Dong,L.N.Wang,D.Y.Ji,and H.Q.Cheng,2013:Introduction of the CMIP5 experiments carried out by BNU-ESM.Progressus Inquisitlones de Mutatlone Climatis,9(4),291–294.

Wei,T.,and Coauthors,2012:Developed and developing world responsibilities for historical climate change and CO2mitigation.Proceedings of the National Academy of Sciences of the United States of America,109(32),12 911–12 915.

Wei,T,W.J.Dong,W.P.Yuan,X.D.Yan,and Y.Guo,2014: Inf l uence of the carbon cycle on the attribution of responsibility for climate change.Chinese Science Bulletin,59(15), 1459–1467.

Weller,E.,and W.J.Cai,2013:Realism ofthe Indian Ocean dipole in CMIP5 models:The implications for climate projections.J.Climate,26(17),6649–6659.

Wuebbles,D.,and Coauthors,2013:CMIP5 climate model analyses:Climate extremes in the united states.Bull.Amer.Meteor. Soc.,95(4),571–583.

Ye,D.Z,C.B.Fu,J.J.Ji,W.J.Dong,J.H.Lu,G.Wen,and X.D.Yan,2001:Orderly human activities and subsistence environment.Advances in Earth Science,16(4),453–460,(in Chinese).

Ye,D.Z.,C.B.Fu,W.J.Dong,G.Wen,and X.D.Yan,2003: Some advance in global change science study.Chinese J.Atmos.Sci.,27(4),435–450.(in Chinese)

Ye,D.Z,J.J.Ji,Z.W.Yan,Y.Wang,and X.D.Yan,2009: Anthroposphere—An interactive component in the Earth system.Chinese J.Atmos.Sci.,33(3),409–415.(in Chinese)

Zhou,T.J.,and Coauthors,2014:Chinese contribution to CMIP5: An overview of f i ve Chinese models’performances.Journal of Meteorology Research,28(4),481–509.

:Yang,S.L.,and Coauthors,2015:A brief introduction to BNU-HESM1.0 and its earth surface temperature simulations.Adv.Atmos.Sci.,32(12),1683–1688,

10.1007/s00376-015-5050-6.

11 February 2015;revised 10 June 2015;accepted 15 June 2015)?

DONG Wenjie Email:dongwj@bnu.edu.cn