XIE Bing,ZHANG Hua,*,YANG Dong-Dong,WANG Zhi-Li

aLaboratory for Climate Studies of China Meteorological Administration,National Climate Center,China Meteorological Administration,Beijing 100081,China

bCollaborative Innovation Center on Forecast and Evaluation of Meteorological Disasters,Nanjing University of Information Science&Technology,Nanjing 210044,China

cCollege of Atmospheric Science,Nanjing University of Information Science&Technology,Nanjing 210044,China

dState Key Laboratory of Severe Weather&Key Laboratory of Atmospheric Chemistry of CMA,Chinese Academy of Meteorological Sciences,Beijing 100081,

China

A modeling study of effective radiative forcing and climate response due to increased methane concentration

XIE Binga,b,ZHANG Huaa,b,*,YANG Dong-Dongc,WANG Zhi-Lid

aLaboratory for Climate Studies of China Meteorological Administration,National Climate Center,China Meteorological Administration,Beijing 100081,China

bCollaborative Innovation Center on Forecast and Evaluation of Meteorological Disasters,Nanjing University of Information Science&Technology,Nanjing 210044,China

cCollege of Atmospheric Science,Nanjing University of Information Science&Technology,Nanjing 210044,China

dState Key Laboratory of Severe Weather&Key Laboratory of Atmospheric Chemistry of CMA,Chinese Academy of Meteorological Sciences,Beijing 100081,

China

An atmospheric general circulation model BCC_AGCM2.0 and observation data from ARIS were used to calculate the effective radiative forcing(ERF)due to increased methane concentration since pre-industrial times and its impacts on climate.The ERF of methane from 1750 to 2011 was 0.46 W m-2by taking it as a well-mixed greenhouse gas,and the inhomogeneity of methane increased its ERF by about 0.02 W m-2. The change of methane concentration since pre-industrial led to an increase of 0.31°C in global mean surface air temperature and 0.02 mm d-1in global mean precipitation.The warming was prominent over the middle and high latitudes of the Northern Hemisphere(with a maximum increase exceeding 1.4°C).The precipitation notably increased(maximum increase of 1.8 mm d-1)over the ocean between 10°N and 20°N and signifcantly decreased(maximum decrease＞-0.6 mm d-1)between 10°S and 10°N.These changes caused a northward movement of precipitation cell in the Intertropical Convergence Zone(ITCZ).Cloud cover signifcantly increased(by approximately 4%)in the high latitudes in both hemispheres,and sharply decreased(by approximately 3%)in tropical areas.

Methane;Effective radiative forcing;Climate change

1.Introduction

Global surface air temperatures have increased due to increases in emissions of anthropogenic greenhouse gases(GHGs)since the start of the industrial era(IPCC,2013; Zhang et al.,2014b).Methane(CH4),with a relatively short lifetime(about 12 years),is the most important anthropogenic GHG besides CO2(IPCC,2013).Although the burden of CH4in the atmosphere is signifcantly smaller than that of CO2,the radiative effciency(RE)of CH4is 26.5 times as much as that of CO2(Yashiro et al.,2008;Renaud and Caillol,2011).The absorption of CH4on radiative fux infuence the temperature, especially near the Earth's surface.

The atmospheric CH4concentration almost doubled from 1750 to 2011,the volume mixing ratio from 722×10-9to 1803×10-9(IPCC,2013).CH4increased rapidly until 2000. Then,after a decade of stabilization or slightly decreasing concentrations,the global CH4concentration showed a well-defned increase again in 2007(Dlugokencky et al.,2003; Rigby et al.,2008),measured using ground-based observations(Cunnold et al.,2002;Langenfelds et al.,2002; Dlugokencky et al.,2009)and aircraft profles(Wecht et al., 2012;Worden et al.,2012).The IPCC Fifth Assessment Report(IPCC,2013)showed that the radiative forcing(RF)of CH4was 0.48±0.05 W m-2from 1750 to 2011 based on Myhre et al.(1998).CH4is well mixed in the atmosphere,but its concentrations vary with latitude and altitude,contributing to 2%of the uncertainty in its RF(Freckleton et al.,1998).

In this study,we estimated the effective radiative forcing (ERF)and climate response due to changes in atmospheric methane concentration from 1750 to 2011 using the general circulation model BCC_AGCM2.0 from the National Climate Center of China.

2.Descriptions of model and method

2.1.Model

We used the general circulation model BCC_AGCM2.0 developed by the National Climate Center of China.The horizontal resolution of the model is approximately 2.8°×2.8°,and the vertical direction includes 26 layers,with a rigid lid at 2.9 hPa.This model was based on the Community Atmosphere Model Version 3(CAM3.0)from the National Center for Atmospheric Research,U.S.BCC_AGCM2.0 contains several enhancements in the physics: BCC_AGCM2.0 uses the radiation scheme of BCC_RAD (Beijing Climate Center Radiative Transfer Model),developed by Zhang et al.(2003,2006a,2006b),and the cloud overlap scheme of the Monte Carlo independent column approximation(Zhang et al.,2014a).These schemes increase the accuracy of the sub-grid cloud structure and its radiative transfer process.Further details on BCC_AGCM2.0 can be found in Wu et al.(2010).The model has been used to study the RFs and the subsequent effects on climate due to aerosols(e.g., Zhang et al.,2012;Wang et al.,2013a,2013b,2014,2015; Zhao et al.,2015)and tropospheric ozone(Xie et al.,2016).

2.2.Satellite data

CH4profle data observed by the Atmospheric Infrared Sounder(AIRS)(URL:http://www.nasa.gov/mission_pages/ aqua)was used.AIRS was onboard the Aqua spacecraft and launched by the National Aeronautic and Space Administration(NASA)in May 2002.An overview of the AIRS instrument is given by Aumann et al.(2003).AIRS has 2378 channels,which cover from 649 to 1136,1217 to 1613,and 2169 to 2674 cm-1with high spectral resolution(λ/ Δλ=1200),and the absorption band of CH4is included.The AIRS dataset used in this study is on a global spatial resolution of 1°×1°,and on vertical pressure levels from 1000 h Pa to 1 hPa(divided into 24 levels).More information about characterization and validation of methane products from AIRS can be found in Xiong et al.(2008)and Zhang et al.(2014c).

2.3.Experimental design

Our aim was to calculate the ERF and climate response due to changes in atmospheric CH4concentration.To this end,fve simulations(EXP1,EXP2,EXP3,EXP4,and EXP5)were conducted.EXP1,EXP2 and EXP3 were used to calculate the ERF of CH4at fxed sea surface temperature(SST)(Hurrell et al.,2008).The same model settings were used in these three simulations,only the CH4concentrations were different (Table 1 shows the details).The differences between EXP1 and EXP2(EXP2 minus EXP1)were regarded as the ERF of well mixed CH4,and differences between EXP1 and EXP3 (EXP3 minus EXP1)were the ERF of CH4with spatial variation.Each simulation was run for 15 years.Kristj?ansson et al. (2005)reported that,after a period of adjustment(generally 5 years for the model with prescribed SST and 30 years for the model with a coupled slab ocean model),the global mean surface air temperature reached equilibrium.Therefore,the results from the last 10 years of the 15-year simulations of EXP1,EXP2 and EXP3 were used to calculate the ERF,as follows:

where ΔFis the net radiation fux(the difference between incoming and outgoing shortwave and longwave radiative fux)at the top of the atmosphere(TOA).

EXP4 and EXP5 were used to calculate the climate response of CH4.We used the same CH4volume mixture ratios in EXP4 and EXP5 as in EXP1 and EXP2,respectively. However,to consider the feedback of the oceans to CH4forcing,a slab ocean model was coupled with BCC_AGCM2.0 in the simulations instead of using the fxed SST.The two simulations were run for 70 years,and we used the results from the last 40 years to discuss the climate response due to changes in CH4concentration.

3.Simulation results and analysis

3.1.ERF

Human activities have increased anthropogenic emissions of CH4since pre-industrial times,leading to increases in its atmospheric concentration and,subsequently,changes in itsERF.Fig.1 shows the simulated ERF(the difference between EXP2 and EXP1)due to change in CH4as well mixed greenhouse gas(WMGHG).A well-defned positive ERF was generally observed near 60°in both hemispheres,with the maximum value of 7 W m-2.However,the ERF was negative in western Siberia,southern Africa,Greenland,and most parts of South America,with a value of approximately-6 W m-2over southern Africa.The negative ERF in those areas might be explained by an increase in low cloud cover.The simulated global mean CH4ERF was 0.46 W m-2,which is consistent with the value reported in IPCC AR5.As Fig.2 shown,CH4concentrations vary with latitude and sharply decrease above the tropopause.In lower troposphere,the concentration of methane was mainly in zonal division,and asymmetry on the Northern and Southern Hemispheres.The volume mixing ratio of CH4was reduced from north to south.The asymmetry of CH4concentration was becoming less distinct with the increase of altitude.In the stratosphere,the volume mixing ratio was symmetrically distributed in both hemispheres,and it got less with the increasing latitude.These spatial variations of CH4had little impact on ERF(less than 0.02 W m-2).

Table 1Experimental design.

Fig.1.Distribution of the effective radiative forcing(ERF)of well mixed atmospheric methane from 1750 to 2011(units:W m-2).Shaded area represents the value at 0.05 signifcance level.

Fig.2.Zonal distribution of the volume mixing ratio of CH4in 2011(×10-6), observed by AIRS.

3.2.Surface air temperature and cloud cover

CH4is a key long-lived GHG that strongly absorbs longwave radiation.The ERF of CH4is generally positive,leading to a warming effect on the Earth's climate system and thus the surface.The difference between EXP5 and EXP4 showed that an increase in the atmospheric CH4concentration since preindustrial times caused an increase of 0.31°C in global mean surface air temperature.As shown in Fig.3a,the surface air temperature increased over the globe except for the small decreases in several high-latitude areas in both hemispheres. Warming over the middle latitudes of the Northern Hemisphere was prominent,with the maximum temperature increase exceeding 1.4°C.There was also signifcant warming (approximately 1.0°C)in the Antarctic area.Fig.3b shows the response of the surface net radiation fux(SNRF)due to the change in CH4.The distribution of change in SNRF was consistent with that of surface temperature over the ocean. There were signifcant increases in SNRF over the high latitudes in both hemispheres.For example,the SNRF over the North Pacifc Ocean increased by more than 6.0 W m-2,and the surface air temperature also increased signifcantly. Changes in cloud cover and heat transportation can also affect surface air temperature.Although the SNRF showed a welldefned decrease over the Indian Ocean,South Pacifc,and the high latitudes of the Southern Hemisphere,the surface air temperature in the same regions did not change accordingly.

Fig.3c and d shows the changes in low-level(below 680 hPa)and high-level(above 440 hPa)cloud cover.Changes in cloud cover directly affect SNRF,thereby infuencing surface air temperature.Increases in low-level cloud result in decreases in SNRF and a cooling effect at the surface,whereas increases in high-level cloud cause increases in surface air temperature due to high-level cloud's warming effect on the Earth's climate system.As shown in Fig.1,the ERF was clearly negative in the western and southern regions of South America,and low-level cloud cover in these areas increased by about 20%(Fig.3c),resulting in marked decreases in surface temperature due to the scattering effect of low-level clouds to solar radiation.The increase in temperature observed over the eastern Japan Sea and Mediterranean regions might be due to increase in high-level cloud cover(Fig.3d).

Fig.3e and f shows the zonally averaged distributions of the changes in cloud cover and relative humidity.There is a high level of correlation between the two variables.The relative humidity showed signifcant increases in most of the troposphere near 70°N and between 10°N and 20°N in the lower troposphere in the Southern Hemisphere,in the higher troposphere in tropical areas,and in most of the troposphere over the Antarctic,and the cloud cover increased by 0.2%-1% in these regions.These increases led to decreases in the SNRF (Fig.3b).In contrast,the relative humidity and cloud cover clearly decreased in the most of troposphere near 60°S and between 30°N and 40°N in the middle to upper troposphere near the equator and in most of stratosphere,resulting in increases in the SNRF(Fig.3b)in some areas.

Fig.3.Climate responses due to changes in atmospheric CH4concentration since pre-industrial times.Distribution of(a)surface air temperature,(b)surface net radiation fux,(c)low cloud,and(d)high cloud.Zonal average distribution of(e)cloud and(f)relative humidity.Shaded area represents the values at 0.05 signifcance level.

3.3.Precipitation and surface water fux

The increase in CH4concentration resulted in a warming effect in the atmosphere and at the surface due to positive ERF at the TOA,which caused an increase of 0.02 kg m-2d-1in global mean surface water fux(SWF)(Fig.4b).The spatial distributions of the changes in SWF and SNRF were similar (Figs.2b and 3b).The SWF dramatically increased(by＞0.12 kg m-2d-1)over most areas of the ocean,especially in the northern Pacifc,western Atlantic,and equatorial Pacifc. In contrast,the SWF showed well-defned decreases due to the decreased SNRF in most areas.In particular,the SWF decreased by approximately 0.14 kg m-2d-1in eastern South America and central Africa.

Fig.4a shows the changes in precipitation due to CH4, which were notable in the Intertropical Convergence Zone. Precipitation signifcantly increased(by＞0.5 mm d-1,with a maximum increase of 1.8 mm d-1)over the ocean between 10°N and 20°N.However,precipitation signifcantly decreased(maximum decrease＞0.6 mm d-1)over the ocean between 10°S and 10°N.Hence,there was a negative correlation between changes in precipitation over the tropics in each hemisphere,with precipitation increased in the Northern Hemisphere and decreased in the Southern Hemisphere.

Fig.4.Climate responses due to changes in atmospheric CH4concentration since pre-industrial times.Distribution of(a)precipitation and(b)surface water fux. Shaded area represents the values at 0.05 signifcance level.

4.Conclusions

The ERF and climate responses due to the change in atmospheric CH4concentration from pre-industrial times(1750) to 2011 were investigated using the atmospheric general circulation model BCC AGCM2.0,in combination with CH4volume mixture ratios from IPCC AR5.The global mean ERF for CH4as WMGHG was 0.46 W m-2,and the spatial variation of methane infuenced the ERF by 0.02 W m-2.The increase in atmospheric CH4led to an increase of 0.31°C and 0.02 mm d-1in global mean surface air temperature and precipitation,respectively.Warming was signifcant in the middle and high latitudes,especially in the Northern Hemisphere,with the maximum warming exceeding 1.4°C.The global distribution of change in precipitation was in line with that of changes in cloud cover,especially near the equator.The precipitation notably increased(maximum increase of 1.8 mm d-1)over the tropical regions of the Northern Hemisphere and sharply decreased(maximum decrease＞-0.6 mm d-1)between 10°S and 10°N,and these changes led the precipitation cell in ITCZ to move northward.In the most of high latitudes in both hemispheres,cloud cover was signifcantly increased(by approximately 4%)and decreased (by approximately 3%)in tropical areas.

Acknowledgments

This work was supported by the National Natural Science Foundation of China(41575002,91644211).

Aumann,H.,Chahine,M.T.,Gautier,C.,et al.,2003.AIRS/AMSU/HSB on the Aqua mission:design,science objectives,data products,and processing systems.IEEE Trans.Geosci.Remote Sens.41,253-264.

Cunnold,D.M.,Steele,L.P.,Fraser,P.J.,et al.,2002.In situ measurements of atmospheric methane at GAGE/AGAGE sites during 1985-2000 and resulting source inferences.J.Geophys.Res.Atmos.107.ACH20-1,Cite ID 4225.

Dlugokencky,E.,Houweling,S.,Bruhwiler,L.,et al.,2003.Atmospheric methane levels off:temporary pause or a new steady-state.Geophys.Res. Lett.30(19),1992.http://dx.doi.org/10.1029/2003GL018126.

Dlugokencky,E.,Bruhwiler,L.,White,J.W.C.,et al.,2009.Observational constraints on recent increases in the atmospheric CH4burden.J.Geophys. Res.Lett.36,L18803.http://dx.doi.org/10.1029/2009GL039780.

Freckleton,R.,Highwood,E.,Shine,K.,et al.,1998.Greenhouse gas radiative forcing:effects of averaging and inhomogeneities in trace gas distribution. Q.J.R.Meteorol.Soc.124,2099-2127.

Hurrell,J.W.,Hack,J.J.,Shea,D.,et al.,2008.A new sea surface temperature and sea ice boundary dataset for the community atmosphere modal.J. Clim.21,5145-5153.

IPCC.Climate Change,2013.The Physical Science Basis.Contribution of Working Group I to the Fifth Assessment Report of the IPCC.Cambridge University Press,Cambridge and New York.

Kristj?ansson,J.E.,Iversen,T.,Kirkev?g,A.,et al.,2005.Response of the climate system to aerosol direct and indirect forcing:role of cloud feedbacks.J.Geophys.Res.110,D24206.http://dx.doi.org/10.1029/ 2005JD006299.

Langenfelds,R.L.,Francey,R.J.,Pak,B.C.,et al.,2002.Interannual growth rate variations of atmospheric CO2and its δ13C,H2,CH4,and CO between 1992 and 1999 linked to biomass burning.Glob.Biogeochem.Cycles 16 (3),21-22.

Myhre,G.,Highwood,E.J.,Shine,K.P.,et al.,1998.New estimates of radiative forcing due to well mixed greenhouse gases.Geophys.Res.Lett.25, 2715-2718.

Renaud,de R.,Caillol,S.,2011.Fighting global warming:the potential of photocatalysis against CO2,CH4,N2O,CFCs,tropospheric O3,BC and other major contributors to climate change.J.Photochem.Photobiol.C 12, 1-19.

Rigby,M.,Prinn,R.G.,Fraser,P.J.,et al.,2008.Renewed growth of atmospheric methane.Geophys.Res.Lett.35,L22805.http://dx.doi.org/ 10.1029/2008GL036037.

Wang,Z.L.,Zhang,H.,Lu,P.,2014.Improvement of cloud microphysics in the aerosol-climate model BCC_AGCM2.0.1_CUACE/Aero,evaluation against observations,and updated aerosol indirect effect.J.Geophys.Res. 119(13),8400-8417.

Wang,Z.L.,Zhang,H.,Zhang,X.Y.,2015.Simultaneous reductions in emissions of black carbon and co-emitted species will weaken the aerosol net cooling effect.Atmos.Chem.Phys.15(7),3671-3685.

Wang,Z.L.,Zhang,H.,Jing,X.W.,et al.,2013a.Effect of non-spherical dust aerosol on its direct radiative forcing.Atmos.Res.120,112-126.

Wang,Z.L.,Zhang,H.,Li,J.,et al.,2013b.Radiative forcing and climate response due to the presence of black carbon in cloud droplets.J.Geophys. Res.Atmos.118,3662-3675.

Wecht,K.J.,Jacob,D.J.,Wofsy,S.C.,et al.,2012.Validation of TES methane with HIPPO aircraft observations:implications for inverse modeling of methane sources.Atmos.Chem.Phys.12,1823-1832.

Worden,J.,Kulawik,S.,Frankenberg,C.,et al.,2012.Profles of CH4,HDO, H2O,and N2O with improved lower tropospheric vertical resolution from Aura TES radiances.Atmos.Meas.Techn.5,397-411.

Wu,T.,Yu,R.C.,Zhang,F.,et al.,2010.The Beijing Climate Center atmospheric general circulation model:description and its performance for the present-day.Clim.Dyn.34,123-147.http://dx.doi.org/10.1007/s00382-009-0594-8.

Xie,B.,Zhang,H.,Wang,Z.,et al.,2016.A modeling study of effective radiative forcing and climate response due to tropospheric ozone.Adv. Atmos.Sci.33(7),819-828.http://dx.doi.org/10.1007/s00376-016-5193-0.

Xiong,X.,Barnet,C.,Maddy,E.,et al.,2008.Characterization and validation of methane products from the Atmospheric Infrared Sounder(AIRS).J. Geophys.Res.113,G00A01.http://dx.doi.org/10.1029/2007JG000500.

Yashiro,Y.,Kadir,W.R.,Okuda,T.,et al.,2008.The effects of logging on soil greenhouse gas(CO2,CH4,and N2O)fux in a tropical rain forest, Peninsular Malaysia.Agric.For.Meteorol.148,799-806.

Zhang,H.,Jing,X.W.,Li,J.,2014a.Application and evaluation of a new radiation code under McICA scheme in BCC_AGCM_2.0.1.Geosci. Model Dev.7,737-754.

Zhang,H.,Nakajima,T.,Shi,G.-Y.,et al.,2003.An optimal approach to overlapping bands with correlated k distribution method and its application to radiative calculations.J.Geophys.Res.108(D20),4641.http:// dx.doi.org/10.1029/2002JD003358.

Zhang,H.,Shi,G.,Nakajima,T.,et al.,2006a.The effects of the choice of the‘K’-interval number on radiative calculations.J.Quant.Spectrosc.Radiat. Transf.98,31-43.

Zhang,H.,Suzuki,T.,Nakajima,T.,et al.,2006b.Effects of band division on radiative calculations.Opt.Eng.45,016002,1-10.

Zhang,H.,Wang,Z.L.,Wang,Z.,et al.,2012.Simulation of direct radiative forcing of aerosols and their effects on East Asian climate using an interactive AGCM-aerosol coupled system.Clim.Dyn.38, 1675-1693.

Zhang,H.,Xie,B.,Zhao,S.Y.,et al.,2014b.PM2.5 and tropospheric O3in China and an analysis of the impact of pollutant emission control.Adv. Clim.Change Res.5(3),136-141.

Zhang,Y.,Xiong,X.Z.,Tao,J.H.,et al.,2014c.Methane retrieval from Atmospheric Infrared Sounder using EOF-based regression algorithm and its validation.Chin.Sci.Bull.59(14),1508-1518.

Zhao,S.Y.,Zhang,H.,Feng,S.,et al.,2015.Simulating direct effects of dust aerosol on arid and semi-arid regions using an aerosol-climate model system.Int.J.Climatol.35(8),1858-1866.http://dx.doi.org/10.1002/ joc.4093.

*Corresponding author.Laboratory for Climate Studies of China Meteorological Administration,National Climate Center,China Meteorological Administration,Beijing 100081,China.

E-mail address:huazhang@cma.gov.cn(ZHANG H.).

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

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

1674-9278/Copyright?2016,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/4.0/).

Received 27 July 2016;accepted 8 December 2016 Available online 14 December 2016