地球物理学报  2018, Vol. 61 Issue (1): 106-118   PDF    
BCC_CSM1.1(m)模式对于夏季亚洲-太平洋涛动的模拟
单机坤1, 梁潇云2 , 吴统文2, 刘向文2, 李巧萍2     
1. 南京信息工程大学大气科学学院, 南京 210044;
2. 中国气象局国家气候中心, 北京 100081
摘要:亚洲-太平洋涛动是北半球夏季亚洲大陆和北太平洋副热带地区对流层中高层扰动温度场上大尺度的东西反相的遥相关现象,其异常变化与亚洲-太平洋地区夏季风气候有着密切的联系.基于欧洲中心的ERA-40再分析资料和国家气候中心BCC_CSM1.1(m)气候系统模式多年的数值模拟结果,本文主要评估了BCC_CSM1.1(m)模式对于夏季亚洲-太平洋涛动的空间分布、指数的时间演变及与其变化所对应的亚洲地区夏季环流异常等方面的模拟能力,结果表明:BCC_CSM1.1(m)模式能够较好地模拟出北半球夏季对流层中高层扰动温度在亚-太地区中纬度存在的西高东低"跷跷板"现象;模式能够模拟出夏季亚洲-太平洋涛动指数的年际变率,但是不能模拟出该指数在20世纪60-70年代明显下降的年代际趋势;模式还能较好地模拟出亚洲-太平洋涛动高低指数年亚洲-太平洋地区夏季环流的异常:指数偏高年份,南亚高压增强,高空西风急流带和热带东风急流均加强,索马里越赤道气流增强,南亚热带季风和东亚副热带季风均增强,东亚季风低压槽加强,西北太平洋副热带高压增强,南亚和东亚北部降水增加,菲律宾地区、中国长江流域-朝鲜半岛-日本一带地区降水减少,反之亦然.
关键词: BCC_CSM1.1(m)模式      亚洲-太平洋涛动      数值模拟      亚洲夏季风     
Simulation of the Asian-Pacific Oscillation using the BCC_CSM1.1(m) model
SHAN JiKun1, LIANG XiaoYun2, WU TongWen2, LIU XiangWen2, LI QiaoPing2     
1. College of Atmospheric Science, Nanjing University of Information Science and Technology, Nanjing 210044, China;
2. National Climate Center, China Meteorological Administration, Beijing 100081, China
Abstract: The Asian-Pacific Oscillation (APO) is a zonal teleconnection pattern of the eddy temperature field in the middle and high troposphere over the extra-tropical Asian-Pacific region. Its anomalous change is closely related to the Asian and Pacific summer monsoon. Based on the monthly European Center for Medium-Range Weather Forecast reanalysis data ERA-40 and numerical simulation of the BCC_CSM1.1(m) model in Beijing Climate Center, this paper evaluates the ability of BCC_CSM1.1(m) model in simulating the summer APO spatial distribution, evolution of APO index and the climate differences between the higher and lower APO index years. Results show that BCC_CSM1.1(m) can well simulate the out-of-phase relationship in the upper-tropospheric eddy temperature between Asia and North Pacific mid-latitudes in Northern Hemisphere summer, and the model is also able to successfully simulate the interannual variability of APO index, but fails to catch the significantly decreasing trend in 1960s-1970s.Furthermore, BCC_CSM1.1(m) can successfully reproduce the Asian and Pacific summer monsoon climate changes corresponding to the APO index variation. With higher APO index conditions, the summer South Asian high and the North Pacific trough are stronger, while the westerly jet stream over Asia and the easterly jet stream over South Asia both strengthen in the upper troposphere. Meanwhile, the Asian low and the North Pacific subtropical high are stronger in the lower troposphere. The anomalous southerlies prevail at 40°E-50°E longitudes and the mid-latitudes of East Asia, and the anomalous westerlies prevail over South Asia. Summer rainfall increases in South Asia and northern of East Asia, while it decreases from the Yangtze River Basin to Japan and the Philippines region, and vice versa.
Key words: BCC_CSM1.1(m) model    Asian-Pacific Oscillation    Numerical simulation    Asian summer monsoon    
0 引言

认识北半球大气环流的内在联系以及理解其中的物理过程是国际大气科学研究的热点之一, 其中夏季北半球的遥相关及其对区域气候的影响长期受到关注.20世纪80年代,Nitta(1987, 1989)首先通过研究菲律宾和日本附近云量反映的对流状况,指出二者存在相反的振荡,并称之为太平洋—日本(PJ)型.与此同时,黄荣辉等(黄荣辉和李维京, 1988, 黄荣辉, 1990, 黄荣辉和孙凤英, 1992)发现东亚夏季风异常不但受到印度季风的影响,更多与西北太平洋副热带高压(简称:西太副高)的季节内和年际变化有关.当西太平洋暖池位于暖位相时,暖水堆积,菲律宾附近对流异常加强,西太副高偏北,长江流域和日本中南部地区偏旱,在位势高度场上表现为菲律宾、日本中部和鄂霍次克海附近的“-+-”波列,并称为东亚—太平洋(EAP)遥相关型.之后,Lau等(Lau, 1992, Lau and Weng, 2002, Lau et al., 2004)的研究指出, 夏季东亚—北美遥相关型使得东亚降水与北美之间存在着密切联系, 并且欧亚和北美夏季气候异常受到北太平洋的调节.Rodwell和Hoskins(2001)认为, 北美中部的气候异常可能是对亚洲季风区加热的一种响应.Ding和Wang(2005)揭示了一个环半球的夏季遥相关型, 这种遥相关在纬向表现为5波的特征, 并与西欧、俄罗斯的欧洲部分、印度、东亚和北美的降水和温度异常有密切关系.最近,Zhao等(2007)发现夏季亚洲和太平洋中纬度对流层扰动温度之间存在着“跷跷板”现象,当亚洲大陆中纬度对流层偏冷时,中、东太平洋中纬度对流层偏暖,反之亦然.与以前的遥相关相比,此遥相关在空间尺度上更大,由于遥相关中心主要出现在亚洲与太平洋区域内,把这种遥相关变率称作亚洲—太平洋涛动(Asian-Pacific Oscillation,简称APO).

已有的研究表明,夏季亚洲—太平洋区域大气环流(包括南亚高压,西太副高,副热带西风和热带东风急流)、亚洲季风降水及西太平洋(包括中国近海)热带气旋均与APO有着密切的联系(Zhao et al., 2007; 赵平等,2008周波涛等,2008; 邹燕和赵平,2009Zhou and Zhao, 2010Zhao et al., 2011a, 2011b, 2012Liu et al., 2011Chen et al., 2013).因此,研究夏季APO的变化规律对于理解东亚夏季风的变化机理具有重要意义.

气候模式是理解气候演变规律,探讨气候现象机理的重要工具.由于气候模式的本身尚存在系统性的误差,这在很大程度上增加了气候模拟的不确定性,而模式评估本身也就显得尤为重要.由于APO这种遥相关现象的主要特征和与之相联系的气候异常可以被部分气候模式所描述(Zhao et al., 2010, 2011b; Chen et al., 2013; 陈晓龙等, 2013; Huang et al., 2013, 2014),且近年来,中国气象局国家气候中心发展了新一代的海-陆-冰-气多圈层相互作用的气候系统模式(BCC_CSM1.1(m)),因此本文作者想评估一下BCC_CSM1.1(m)模式对于夏季APO的模拟效果.BCC_CSM1.1(m)模式不仅参与了CMIP5比较计划(Xin et al., 2012),为气候变化研究提供了大量的数据,目前还被用于短期气候预测业务.研究表明,BCC_CSM1.1(m)模式对20世纪气温的平均气候及其年际变化趋势具有较强的模拟能力,对全球百年气候变化也有一定的模拟能力(Xin et al., 2012, 2013a, 2013b),能够模拟出1971—2000年全球地表气温的空间分布特征(Wu et al., 2014);能合理再现全球年平均降水的基本分布特征,也能较合理再现热带降水年循环模态的基本分布特征,尤其季风模态中降水与环流关于赤道反对称的特征(张莉等,2013);能够模拟出中亚地区显著增温以及感热通量、长/短波净辐射等要素由南向北递减的总体趋势(姜燕敏等, 2015);能模拟出ENSO现象的一些主要特征,但是存在ENSO周期偏短问题(张芳等,2013)等.在预测方面,BCC_CSM1.1(m)模式可以预测出亚洲—太平洋地区夏季风的季节内和年际变率(Liu X W et al., 2014, 2015);对夏季中低纬海温具有一定的预测能力,且在低纬地区的预测技巧尤为出色(汪栩加等,2015).虽然已有一些关于BCC_CSM1.1(m)的模拟和预测效果的评估工作,但是,该模式对影响季风变化的大尺度海-陆热力对比和APO的模拟能力如何,此前缺乏相关研究.

我们知道,目前的气候模式对于东亚夏季降水的模拟和预测,技巧非常低(Wang et al., 2009; Huang et al., 2013; Jiang et al., 2016).如果能够找到一个与东亚夏季降水具有较好的相关关系,并在气候模式中具有较高预测技巧的一个预测因子,就可以利用间接的方式做好东亚降水预测.APO指数是温度场的一种表现形式,在气候模式中具有较好的模拟和预测技巧(Zhao et al., 2010, 2011b; Chen et al., 2013Huang et al., 2013, 2014Liu G et al., 2015),已有的研究也表明该指数和东亚夏季降水具有较好的相关(Zhao et al., 2008, 2010, 2012).基于以上考虑,本文从对流层中上层温度场和APO指数的角度,研究BCC_CSM1.1(m)模式对纬向海-陆热力差异的描述能力,分析模式结果与观测产生偏差的原因,旨在回答如下的科学问题:BCC_CSM1.1(m)模式对海-陆热力差异特别是APO活动中心的模拟能力如何?模式是否可以模拟出APO指数的年际变化?与APO异常对应的区域气候异常在模式中是否有体现?

文章结构如下:第1节是所用观测资料和模式介绍,第2节是对BCC_CSM1.1(m)模式模拟结果分析与讨论,第3节是结论.

1 模式和资料

本文研究所用的观测资料为欧洲中期数值预报中心(European Center for Medium-Range Weather Forecast, ECMWF)提供的月平均的1957年9月—2002年8月的ERA-40再分析资料(Uppala et al., 2005),水平分辨率为2.5°×2.5°.

所用的模式为BCC_CSM1.1(m)模式.该模式是由中国气象局国家气候中心研发的第二代气候系统模式,最初研发的是BCC_CSM1.1版本,BCC_CSM1.1(m)为升级版本.BCC_CSM1.1(m)版本模式的大气模式分量为BCC_AGCM2.2,水平分辨率为T106(网格距约为1°),垂直方向为26层,关于BCC_AGCM2.2动力框架和物理过程的具体介绍可参见文献(Wu et al., 2013, 2014).它的陆面模式分量为BCC_AVIM1.0,是一个大气、植被、土壤相互作用模式,是在NCAR陆面模式CLM3物理模块的基础上,引入了中国发展的动态植被和土壤碳循环模型AVIM2(Ji et al., 2008; Wu et al., 2013).其海洋模式分量是在美国地球物理流体力学实验室(The Geophysical Fluid Dynamics Laboratory, GFDL)发展的MOM4_L40基础上改进的MOM4_L40v2,水平分辨率为1°×1°,经向在热带地区加密到(1/3)°,垂直方向为40层(Wu et al., 2013).其海冰模式分量为GFDL发展的海冰模式SIS,水平分辨率与MOM_L40相同(Winton, 2000).

文中所用模式数据为BCC_CSM1.1(m)模式参与CMIP5比较计划中的历史试验数据,历史试验的时间长度为:1850—2012年,这里为了和ERA-40资料进行对比,截取了相同的时间跨度1958—2002年.文中的夏季是指6、7、8三个月的平均,气候平均是指1958—2002年45年的平均.

2 结果分析

赵平等(2007)将亚洲和太平洋地区上空对流层中高层(500~200 hPa)的扰动温度(T′)之差定义为APO指数,即

其中,T′=T-T, T是空气温度,TT的纬向平均.模式对于500~200 hPa扰动温度场的模拟能力关系到APO的模拟.因此,本文首先考察BCC_CSM1.1(m)模式对500~200 hPa(本文中提到对流层中高层时,特指这个层次,以下不再赘述)扰动温度的模拟能力.在此基础上再分析对APO指数的年际变化和年代际变化趋势的模拟水平,以及与APO变率相联系的亚洲地区气候的异常.

2.1 扰动温度气候特征的模拟

图 1a是ERA-40再分析资料得到的气候平均的夏季对流层中高层T′的分布,亚洲大陆地区上空为正的T′值区,中心位于青藏高原上空,中心值达6 ℃以上;太平洋中东部上空为负的T′值区,中心位于东太平洋的副热带地区,中心值达-3 ℃以下.两个区的T′之间存在着显著的反位相关系.BCC_CSM1.1(m)模式对对流层中高层T′气候特征的模拟与ERA-40再分析资料得到的基本一致(图 1b),只是模式模拟的负值中心位置较观测要偏东,正、负值中心强度均较观测的要略偏弱.由于APO关注的是大尺度陆地和海洋上空温度的纬向对比,这些局地差异并不影响模式对大尺度纬向温度平均态的模拟能力.模式较好地再现了欧亚大陆和北太平洋大尺度海-陆分布所造成的纬向热力差异.

图 1 夏季对流层中高层T′的气候分布(单位:℃) (a) ERA-40再分析资料;(b) BCC_CSM1.1(m)模式结果. Fig. 1 Climatology of summer mean upper-tropospheric T′ (unit: ℃) during 1958—2002 for (a) ERA-40 reanalysis and (b) the BCC_CSM1.1(m) model

图 1a显示北半球夏季对流层中高层扰动温度的正负值中心均出现在30°N纬度附近,因此沿30°N的经度-高度剖面再分析扰动温度的垂直分布.由于对流层的正、负T′中心均出现在500~200 hPa高度上(图 2a),因此用500~200 hPa平均代表对流层中高层的状况是合理的.图 2b显示,模式能够再现对流层中亚洲大陆上空上下一致的暖和北太平洋上空上下一致的冷.

图 2 夏季T′沿30°N的经向-高度剖面(单位:℃) (a) ERA-40再分析资料;(b) BCC_CSM1.1(m)模式结果.图中阴影为地形高度. Fig. 2 Longitude-height cross section of summer mean T′ (unit:℃) along latitude 30°N for (a) ERA-40 reanalysis and (b) BCC_CSM1.1(m) model Shaded areas denote terrain height.

将北半球夏季对流层中高层T′做EOF分析,得到的第一分量(EOF1)模态见图 3.在ERA-40中第一模态解释了总方差的26%(图 3a),亚洲大陆和北太平洋之间的“跷跷板”现象非常明显,亚洲大陆上空为正值区,中心强度超过0.02,北太平洋上空为负值区,中心强度低于-0.015.表现为北半球夏季中纬度大尺度海陆热力差异.BCC_CSM1.1(m)模式对流层中高层T′的EOF1模态与观测的EOF1模态基本一致(图 3b),这说明模式能较好地模拟出北半球夏季对流层高层温度场的主要特征.

图 3 夏季对流层中高层T′的EOF1模态(×0.01) (a) ERA-40再分析资料;(b) BCC_CSM1.1(m)模式结果. Fig. 3 Unrotated EOF1 of the normalized summer mean upper-tropospheric T′ (×0.01) for (a) ERA-40 reanalysis and (b) the BCC_CSM1.1(m) model

以上分析显示,BCC_CSM1.1(m)模式能够较好再现北半球中纬度大尺度纬向温度梯度的平均状况,也能较好地刻画出整个对流层高低层温度一致的分布,这是利用该模式研究APO变化机理的基础.

2.2 APO指数的模拟

图 4a4b分别为ERA-40再分析资料和BCC_CSM1.1(m)模式模拟的1958—2002年APO指数序列,代表东亚—北太平洋上空的纬向海陆热力差异.模式模拟的APO指数年际变率偏弱,其标准差(0.28)比ERA-40再分析资料(0.48)的偏小.这是由于APO指数定义中的区域陆地和海洋上空温度指数序列在模式的模拟结果中均较ERA-40的要偏弱的原因(图 4c4d).虽然,模式模拟的陆地和海洋上空的温度指数均较ERA-40要偏弱,但是在模式中两者之间仍是显著的反相关系,相关系数为-0.75,这和用ERA-40再分析资料计算得到的两者相关系数-0.77非常接近.再来分析模式对APO指数年代际变化趋势的模拟.ERA-40再分析资料结果显示,1975年之前,APO基本处于正位相状态,1975年之后APO维持一个负位相状态(图 4a).20世纪60和70年代,APO指数明显下降,主要有两个原因,一是由于陆地区域对流层中高层在此时段有个显著的降温(降温趋势为0.65 ℃/10a),二是由于海洋区域对流层高层在此时段有个明显的升温(升温趋势为0.37 ℃/10a),两者共同作用造成的(图 4c).这与Yu等(2004)Zhou和Zhang(2009)的研究结论一致.20世纪80和90年代陆地和海洋上空的温度都趋于维持一个稳定状态,因此,APO指数没有明显的变化趋势.与ERA-40再分析资料相比,模式未能模拟出APO指数60和70年代的明显下降趋势(图 4b),主要原因是由于陆地上空60和70年代的明显减弱,模式没能模拟出来(图 4d).模式能够模拟出60和70年代海洋上空对流层中高层温度的升温趋势,但是模式模拟的升温率(升温趋势为0.16 ℃/10a)较观测要偏小.Huang等(2014)的研究显示,不仅BCC_CSM1.1 (m)模式不能模拟出APO指数20世纪60年代之后明显的减弱趋势,欧洲中心、英国气象局和法国国家气象局的耦合模式也都存在这个问题.模式模拟的APO指数与观测的APO指数的相关仅为0.16,低于90%信度检验.而将观测的APO指数去掉线性趋势后,两者的相关系数为0.30,通过了95%信度检验.这说明,模式还是可以模拟出APO指数的年际变率的.

图 4 夏季APO(a,b)指数和亚洲、北太平洋(c,d)上空T′指数(单位:℃) (a, c) ERA-40再分析资料;(b, d) BCC_CSM1.1(m)模式结果. Fig. 4 Summer APO index (a, b) and Asian and Pacific T′ indices (c, d) (unit:℃) for (a, c) ERA-40 reanalysis and (b, d) BCC_CSM1.1(m) model
2.3 APO与T′关系的模拟

用APO指数和对流层中高层T′计算相关(图 5a),得到相关系数的空间分布与T′做EOF分析得到的第一模态的空间分布非常相似(图 3a).亚洲大陆上空为正的相关区,北太平洋中东部上空为负的相关区,且大部分区域都通过了95%的信度检验.说明在北半球夏季中纬度对流层中高层确实存在一个扰动温度场的遥相关型.这种大尺度的遥相关在BCC_CSM1.1(m)模式中也可以得到较好的模拟(图 5b).

图 5 夏季APO指数和对流层中高层T′的相关(×0.1) (a) ERA-40再分析资料;(b) BCC_CSM1.1 (m)模式结果.阴影区为通过95%信度检验区. Fig. 5 One-point correlation coefficient (×0.1) between APO index and upper-tropospheric T′ during summer for (a) ERA-40 reanalysis and (b) the BCC_CSM1.1(m) model Shaded areas are at 95% confidence level.

再用APO指数回归得到的T′,检验遥相关型的垂直分布.图 6a显示了沿30°N回归的T′,显著的正异常区从亚洲大陆上空近地面延伸至对流层高层200 hPa,中心从青藏高原中西部上空对流层中层,随着高度的上升向东延伸到对流层高层青藏高原中部上空,中心值达0.8以上.显著的负异常区域位于北太平洋上空800~200 hPa之间的对流层,中心值也达-0.8以下,中心位置在对流层高层200~300 hPa之间.在垂直分布上,APO是一个倾斜的“跷跷板”.

图 6 夏季APO指数回归的T (a, c, e) ERA-40再分析资料;(b, d, f) BCC_CSM1.1 (m)模式结果.(a, b)沿30°N的高度-经度剖面;(c, d) 250 hPa层;(e, f) 1000 hPa层.阴影区为通过95%信度检验区. Fig. 6 Regressed T′ against APO index during summer for (a, c, e) ERA-40 reanalysis and (b, d, f) BCC_CSM1.1(m) model (a, b) Longitude-height cross section along latitude 30°N. (c, d) At 250 hPa. (e, f) At 1000 hPa. Black shaded areas denote terrain height in (a, b). Grey shaded areas are at 95% confidence level.

在250 hPa层(图 6c)上,亚洲地区上空为显著的正异常区域,中心位于亚洲的中纬度地区,中心值大于1.5.显著的负异常区域出现在北太平洋中纬度地区上空,中心值超过了0.9.与对流层高层回归的T′相比,地面回归的T′显著的正异常区域仅仅出现在东亚的东北部,中心强度也可以达0.9,但是北太平洋地区仅仅在东太平洋沿岸为负异常区域,其余大部为正异常区域,北太平洋中纬度地区还存在一个显著的正异常区域(图 6e),这说明APO在对流层中上层比较显著.

与ERA-40再分析资料的结果对比,BCC_CSM1.1(m)模式可以模拟出APO在垂直方向上的倾斜“跷跷板”结构(图 6b),也可以模拟出对流层高层T′的亚洲陆地上空大范围的显著正异常和北太平洋上空显著的大范围的负异常(图 6d),以及近地面大范围的反向现象(图 6f).

2.4 APO异常年亚-太季风区环流和降水的模拟

为了考察与APO相关的亚-太季风区夏季气候异常特征,挑选了图 3a中APO指数最强的7年(1961,1962,1973,1960,1963,1967和1978),称为高指数(HI)年,最弱的7年(1992,1987,1983,1993,1976,1997和1980),称为低指数(LI)年.通过高低指数年合成的差异,分析与APO指数相联系的亚-太季风区夏季环流及降水的异常.

夏季,亚洲对流层高层高度场上最明显的特征是南亚高压.图 7a为气候平均的对流层高层200 hPa扰动高度场H′,这里与扰动温度的定义相同,扰动高度H′=H-HH为位势高度场,H为位势高度场H的纬向平均.与对流层中高层T′的气候分布对应比较一致(图 1a),对流层高层H′的正值亦是主要位于亚洲大陆的青藏高原地区上空,在北太平洋上空的对流层高层为负的H′,对应于对流层高层北太平洋地区上空的洋中槽.而在对流层低层(图 7b),亚洲大陆上空为高度场的负值区,对应夏季的大陆热低压,北太平洋上空为高度场的正值区,对应洋面上的冷高压,即副热带高压.

图 7 夏季H′的气候分布(a, b)以及APO高低指数年H′的差异(c—f)(单位:dagpm) (a, b, c, d) ERA-40再分析资料;(e, f) BCC_CSM1.1 (m)模式结果. (a, c, e) 200 hPa层;(b, d, f) 850 hPa层.阴影区为通过95%信度检验区. Fig. 7 Climatology of the summer mean H′ (unit: dagpm) during 1958—2002 at (a) 200 hPa and (b) 850 hPa, composite differences of H′ between HI and LI cases at (c, e) 200 hPa and (d, f) 850 hPa for (a, b, c, d) ERA-40 reanalysis and (e, f) the BCC_CSM1.1(m) model Shaded areas are at 95% confidence level.

图 7c为APO高低指数年对流层高层200 hPa H′的差异.显著的正异常出现在亚洲地区上空,显著的负异常出现在北太平洋副热带上空.表明APO指数偏高的年份,对流层高层的南亚高压增强,北太平洋上空的洋中槽加强.图 7d为APO高低指数年对流层低层850 hPa H′的差异.显著的负异常出现在亚洲大陆上空,显著的正异常出现在北太平洋中东部.表明APO指数偏高的年份,亚洲低压偏强,副热带高压偏强、偏东、偏北.图 7e图 7f分别是BCC_CSM1.1(m)模式模拟的APO高低指数年对流层高层200 hPa和对流层低层850 hPa H′的差异.可以看到,BCC_CSM1.1 (m)模式可以模拟出APO高低指数年对应的对流层高层高度场上东西“跷跷板”现象,也可以模拟出对流层低层高度场上的反相“跷跷板”现象.

西风急流是亚-太地区夏季季风系统的重要成员之一.气候上夏季西风急流位于青藏高原的北侧,热带东风急流位于青藏高原的南侧.急流的位置与青藏高原上空温度的分布关系非常密切.图 8为夏季高低APO指数年纬向风u沿30°N的差异,当APO指数偏高的年份,青藏高原北侧为正的纬向风异常,表明西风急流有所加强,青藏高原南侧为负的纬向风异常,表明热带东风急流也有所加强(图 8a).对于高空急流在APO指数高低年的变化,BCC_CSM1.1(m)模式可以再现出来,只是与ERA-40再分析资料结果相比模式模拟的变化强度偏弱(图 8b).

图 8 沿90°E夏季APO高低指数年纬向风(单位:m·s-1)的差异 (a) ERA-40再分析资料;(b) BCC_CSM1.1(m)模式结果.阴影区为通过95%信度检验区. Fig. 8 Latitude-height cross section of composite difference of the summer mean zonal wind (unit: m·s-1) between HI and LI cases along longitude 90°E for (a) ERA-40 reanalysis and (b) the BCC_CSM1.1(m) model Shaded areas are at 95% confidence level.

我们知道,亚洲季风区夏季对流层低层盛行偏南风,南亚季风区盛行偏西南风,东亚季风区盛行偏东南风,来自海洋上的偏南暖湿气流给亚洲大陆带来丰沛的降水和雨季.图 9a为高低APO指数年夏季对流层低层850 hPa流场的差异.异常的南风出现在40°E—50°E附近非洲大陆东部的索马里地区和东亚的北部地区,异常的偏西风出现在阿拉伯海北部—印度半岛—孟加拉湾—中南半岛一带地区,这表明夏季APO指数偏高的年份,索马里越赤道气流将加强,南亚夏季风和东亚夏季风均加强.与低层风场相对应,在APO指数偏高年份,印度半岛的中南部、中南半岛的中部、中国的华南地区和华北地区降水将增加,而菲律宾及其附近的西太平洋地区和中国长江流域地区—朝鲜半岛南部—日本岛地区一带的降水将减少(图 10a).对于APO高低指数年对流层低层风场的变化(图 9b)和对应降水量的变化(图 10b),BCC_CSM1.1(m)模式都可以较好地模拟出来.

图 9 夏季APO高低指数年850 hPa风场的差异(单位:m·s-1) (a) ERA-40再分析资料;(b) BCC_CSM1.1 (m)模式结果.阴影区为通过95%信度检验区. Fig. 9 Composite difference of the summer mean 850 hPa wind (unit: m·s-1) between HI and LI cases for (a) ERA-40 reanalysis and (b) the BCC_CSM1.1(m) model Shaded areas are at 95% confidence level.
图 10图 9,但为降水量(单位:mm·day-1) Fig. 10 Same as Fig. 9 but for precipitation (unit: mm·day-1)
2.5 讨论

APO实质是夏季亚洲大陆和北太平洋海陆东西热力差异的结果,Zhao等(2010, 2011)的研究表明,热带太平洋海温热源仅能在对流层低层形成较小尺度的类APO形态,只有青藏高原热源的存在,才可以形成与观测一致的APO形态.为什么BCC_CSM1.1(m)模式能够较好地模拟出夏季APO遥相关模态以及与之异常相关的亚-太地区环流与降水异常特征,其原因在哪里?为了回答这个问题,也参考以往的研究成果,作者考察了模式对于视热源Q1的模拟情况.参考Yanai等(1992)的工作计算了Q1,这里Q1=SH+R+LHSH为感热通量,R为气柱吸收的净辐射,LH为潜热通量.Q1为正值表示为热源,Q1为负值则为冷源.ERA-40再分析资料的计算结果显示(图 11a),夏季欧亚大陆大部Q1为正值,是热源,中心位于中纬度40°N附近,强度在60 W·m-2以上,北太平洋的热带地区大部Q1为正值,也是热源,而北太平洋的中纬度地区大部Q1为负值,是冷源.在北半球夏季的中纬度地区,亚洲大陆和北太平洋东西热力差异最大.热的亚洲大陆加热其上空大气,造成温度的正偏差,位于亚洲大陆南部的青藏高原由于地形较高,其热源的影响可以延伸至对流层中高层,北太平洋东部冷源造成其上空的气温降低,产生温度的负偏差.BCC_CSM1.1(m)模式可以较好地再现北半球夏季Q1的空间分布(图 11b).图 11c为夏季亚洲大陆区域平均的视热源Q1距平的时间演变.它和图 4c的亚洲地区对流层中高层T′指数之间具有较高的相关,相关系数为0.34,超过了95%的信度检验,同时它还与图 4a的APO指数之间也具有较高的相关(相关系数为0.3,达到了95%的信度检验).从图 11c可以看到Q1也在20世纪60、70年代存在一个明显减弱的趋势,与APO指数的年代际趋势十分一致.这进一步证明,在夏季APO的形成中,亚洲大陆地区热源是起主要作用的.图 11d为BCC_CSM1.1(m)模式模拟的亚洲大陆地区Q1的情况,模式没能模拟出亚洲大陆夏季Q1减弱的年代际趋势,但是模式模拟的亚洲大陆的Q1与去除线性趋势后的ERA-40资料的亚洲大陆Q1两者之间的相关系数为0.34,超过了95%信度检验.这说明BCC_CSM1.1(m)模式既可以模拟出Q1的空间分布,又可以较好地模拟出对APO形成中起主要作用的亚洲大陆地区Q1的年际变化,所以它可以较好地模拟出夏季APO遥相关模态以及与之异常相关的亚-太地区环流与降水异常特征.

图 11 视热源Q1的气候分布(a,b)和亚洲区域平均Q1的时间演变(c,d)(单位:W·m-2) (a, c) ERA-40再分析资料;(b, d) BCC_CSM1.1 (m)模式结果. Fig. 11 Climatology of summer mean apparent heat source Q1 (unit: W·m-2) and temporal evolution of Asian regional mean Q1 during 1958—2002 for (a) ERA-40 reanalysis and (b) the BCC_CSM1.1(m) model

对比图 1a图 1b,模式对于位于北太平洋地区上空的对流层负的扰动温度,模拟的中心略偏东,这可能是模式对于北太平洋东部的冷源模拟较观测偏弱的原因所造成.亚洲大陆和青藏高原上空Q1值在模式的模拟结果中较观测要偏小,这正好解释了模式对于对流层中高层正的扰动温度强度模拟的偏弱.

模式对于夏季APO指数在20世纪60—70年代的年代际减弱趋势没能模拟出来,虽然这是目前气候系统模式普遍存在的问题,但是作为身处亚洲季风区的国家之一,中国的气候系统模式研发者应该在这个问题上进行深入分析,查找原因,努力改善动力模式系统对于东亚地区气候特征的模拟效果.模式对于亚洲大陆及青藏高原地区视热源Q1的模拟存在偏低,也需要进行深入分析,也许对这个问题的改进是提高海气耦合模式对东亚地区环流和气候特征模拟效果的关键所在,这对于提高模式在东亚地区气候预测也具有重要意义.

3 结论

APO是夏季亚洲大陆和北太平洋地区热力差异的一个主要模态,与亚洲—太平洋地区夏季风气候有着密切的联系(Zhao et al., 2012),因此评估BCC_CSM1.1(m)模式对于夏季APO的模拟效果,可以更好地了解该模式对东亚季风区夏季气候特征的模拟能力.

基于欧洲中心的ERA-40再分析月平均资料和国家气候中心BCC_CSM1.1(m)气候系统模式45年的数值模拟结果,分析了该模式对于夏季亚—太地区对流层高层温度场存在的大尺度遥相关APO的空间分布、指数的时间演变和与其变化所对应的亚洲地区的夏季环流异常等方面的模拟能力,结果表明:

(1) BCC_CSM1.1(m)模式能够较好地模拟出北半球夏季对流层中高层扰动温度在亚—太地区中纬度存在的西高东低“跷跷板”现象,也可以较好地模拟出扰动温度第一主分量模态的空间分布.

(2) BCC_CSM1.1(m)模式能够模拟出夏季APO指数的年际变率,但是不能模拟出APO指数在20世纪60—70年代明显下降的年代际趋势.原因是模式不能模拟出亚洲大陆上空扰动温度在60—70年代的变冷趋势.

(3) APO指数和扰动温度之间的关系在BCC_CSM1.1(m)模式中也能够得到较好的再现.

(4) BCC_CSM1.1(m)模式还能较好地模拟出APO高低指数年亚洲—太平洋地区夏季环流的异常.APO高指数年:对流层高层,南亚高压增强,高压北侧的西风急流带加强,高压南侧的热带东风急流也加强,北太平洋上空的洋中槽加深,东西热力差异增大;对流层低层,索马里越赤道气流增强,南亚季风区西风增强,东亚季风低压槽加强,西北太平洋副热带高压增强,东亚季风也加强;南亚和东亚北部降水增加,菲律宾地区、中国长江流域—朝鲜半岛—日本一带地区降水减少.反之亦然.

致谢

感谢中国气象局国家气候中心李维京研究员和何文平研究员对本研究提出的宝贵建议.

参考文献
Chen J M, Zhao P, Yang S, et al. 2013. Simulation and dynamical prediction of the summer Asian-Pacific oscillation and associated climate anomalies by the NCEP CFSv2. Journal of Climate, 26(11): 3644-3656. DOI:10.1175/JCLI-D-12-00368.1
Chen X L, Zhou T J, Zou L W. 2013. Variation characteristics of the Asian-Pacific Oscillation in boreal summer as simulated by the LASG/IAP climate system model FGOALS_gl. Acta Meteorologica Sinica, 71(1): 23-37.
Ding Q H, Wang B. 2005. Circumglobal teleconnection in the Northern Hemisphere summer. Journal of Climate, 18(17): 3484-3505.
Huang R H, Li W J. 1988. Influence of heat source anomaly over the western tropical Pacific on the subtropical high over East Asia and its physical mechanism. Chinese Journal of Atmospheric Sciences, 12(S1): 107-116.
Huang R H. 1990. Studies on the teleconnections of the general circulation anomalies of East Asia causing the summer drought and flood in China and their physical mechanism. Chinese Journal of Atmospheric Sciences, 14(1): 108-117.
Huang R H, Sun F Y. 1992. Interannual variation of the summer teleconnection pattern over the northern hemisphere and its numerical simulation. Chinese Journal of Atmospheric Sciences, 16(1): 52-61.
Huang Y Y, Wang H J, Zhao P. 2013. Is the interannual variability of the summer Asian-Pacific Oscillation predictable?. Journal of Climate, 26(11): 3865-3876. DOI:10.1175/JCLI-D-12-00450.1
Huang Y Y, Wang H J, Fan K. 2014. Improving the prediction of the summer Asian-Pacific Oscillation using the interannual increment approach. Journal of Climate, 27(21): 8126-8134. DOI:10.1175/JCLI-D-14-00209.1
Ji J J, Huang M, Li K R. 2008. Prediction of carbon exchange between China terrestrial ecosystem and atmosphere in 21st century. Science in China Series D:Earth Sciences, 51(6): 885-898. DOI:10.1007/s11430-008-0039-y
Jiang D B, Tian Z P, Lang X M. 2016. Reliability of climate models for China through the IPCC Third to Fifth Assessment Reports. International Journal of Climatology, 36(3): 1114-1133. DOI:10.1002/joc.2016.36.issue-3
Jiang Y M, Huang A N, Wu H M. 2015. Evaluation of the performance of Beijing Climate Center Climate System Model with different horizontal resolution in simulating the annual surface temperature over Central Asia. Chinese Journal of Atmospheric Sciences, 39(3): 535-547.
Lau K M. 1992. East Asian summer monsoon rainfall variability and climate teleconnection. Journal of the Meteorological Society of Japan, 70(1B): 211-241. DOI:10.2151/jmsj1965.70.1B_211
Lau K M, Weng H Y. 2002. Recurrent teleconnection patterns linking summertime precipitation variability over East Asia and North America. Journal of the Meteorological Society of Japan, 80(6): 1309-1324. DOI:10.2151/jmsj.80.1309
Lau K M, Lee J Y, Kim K M, et al. 2004. The North Pacific as a regulator of summertime climate over Eurasia and North America. Journal of Climate, 17(4): 819-833. DOI:10.1175/1520-0442(2004)017<0819:TNPAAR>2.0.CO;2
Liu G, Zhao P, Chen J M. 2011. A 150-year reconstructed summer Asian-Pacific Oscillation index and its association with precipitation over eastern China. Theoretical and Applied Climatology, 103(1-2): 239-248. DOI:10.1007/s00704-010-0294-7
Liu G, Zhao P, Chen J M, et al. 2015. Preceding factors of summer Asian-Pacific Oscillation and the physical mechanism for their potential influences. Journal of Climate, 28(7): 2531-2543. DOI:10.1175/JCLI-D-14-00327.1
Liu X W, Wu T W, Yang S, et al. 2014. Relationships between interannual and intraseasonal variations of the Asian-western Pacific summer monsoon hindcasted by BCC_CSM1. 1(m). Advances in Atmospheric Sciences, 31(5): 1051-1064. DOI:10.1007/s00376-014-3192-6
Liu X W, Wu T W, Yang S, et al. 2015. Performance of the seasonal forecasting of the Asian Summer Monsoon by BCC_CSM1. 1(m). Advances in Atmospheric Sciences, 32(8): 1156-1172. DOI:10.1007/s00376-015-4194-8
Nitta T. 1987. Convective activities in the tropical western Pacific and their impacts on the Northern Hemisphere summer circulation. Journal of the Meteorological Society of Japan, 65(3): 373-390. DOI:10.2151/jmsj1965.65.3_373
Nitta T. 1989. Global features of the Pacific-Japan Oscillation. Meteorology and Atmospheric Physics, 41(1): 5-12.
Rodwell M J, Hoskins B J. 2001. Subtropical anticyclones and summer monsoons. Journal of Climate, 14(15): 3192-3211. DOI:10.1175/1520-0442(2001)014<3192:SAASM>2.0.CO;2
Uppala S M, Kållberg P W, Simmons A J, et al. 2005. The ERA-40 re-analysis. Quarterly Journal of the Royal Meteorological Society, 131(621): 2961-3012.
Wang B, Lee J Y, Kang I S, et al. 2009. Advance and prospectus of seasonal prediction, Assessment of the APCC/CliPAS 14-model ensemble retrospective seasonal prediction (1980-2004). Climate Dynamics, 33(1): 93-117. DOI:10.1007/s00382-008-0460-0
Wang X J, Zheng Z H, Feng G L, et al. 2015. Summer prediction of sea surface temperature in key areas in BCC_CSM model. Chinese Journal of Atmospheric Sciences, 39(2): 271-288.
Winton M. 2000. A reformulated three-layer sea ice model. Journal of Atmospheric and Oceanic Technology, 17(4): 525-531. DOI:10.1175/1520-0426(2000)017<0525:ARTLSI>2.0.CO;2
Wu T W, Li W P, Ji J J, et al. 2013. Global carbon budgets simulated by the Beijing Climate Center Climate System Model for the last century. Journal of Geophysical Research Atmospheres, 118(10): 4326-4347. DOI:10.1002/jgrd.50320
Wu T W, Song L C, Li W P, et al. 2014. An overview on progress in Beijing Climate Center Climate System Model-Its development and application to climate change studies. Acta Meteorologica Sinica, 72(1): 12-29.
Xin X G, Wu T W, Zhang J. 2012. Introduction of CMIP5 experiments carried out by BCC climate system model. Progressus Inquisitiones De Mutatione Climatis, 8(5): 378-382.
Xin X G, Cheng Y J, Wang F, et al. 2012. Asymmetry of surface climate change under RCP2. 6 projections from the CMIP5 models. Advances in Atmospheric Sciences, 30(3): 796-805.
Xin X G, Wu T W, Zhang J. 2013b. Introduction of CMIP5 experiments carried out with the climate system models of Beijing Climate Center. Advances in Climate Change Research, 4(1): 41-49.
Yanai M, Li C F, Song Z S. 1992. Seasonal heating of the Tibetan Plateau and its effects on the evolution of the Asian summer monsoon. Journal of the Meteorological Society of Japan, 70(1B): 319-351. DOI:10.2151/jmsj1965.70.1B_319
Yu R C, Wang B, Zhou T J. 2004. Tropospheric cooling and summer monsoon weakening trend over East Asia. Geophysical Research Letters, 31(22): L22212, 217-244.
Zhang F, Dong M, Wu T W. 2014. Evaluation of the ENSO features simulations as done by the CMIP5 models. Acta Meteorologica Sinica, 72(1): 30-48.
Zhang L, Wu T W, Xin X G, et al. 2013. The annual modes of tropical precipitation simulated by the Beijing Climate Center Climate System Model (BCC_CSM). Chinese Journal of Atmospheric Sciences, 37(5): 994-1012.
Zhao P, Zhu Y N, Zhang R H. 2007. An Asian-Pacific teleconnection in summer tropospheric temperature and associated Asian climate variability. Climate Dynamics, 29(2-3): 293-303. DOI:10.1007/s00382-007-0236-y
Zhao P, Chen J M, Xiao D, et al. 2008. Summer Asian-Pacific oscillation and its relationship with atmospheric circulation and monsoon rainfall. Acta Meteorologica Sinica, 66(5): 716-729.
Zhao P, Cao Z H, Chen J M. 2010. A summer teleconnection pattern over the extratropical Northern Hemisphere and associated mechanisms. Climate Dynamics, 35(2): 523-534.
Zhao P, Yang S, Wang H J, et al. 2011a. Interdecadal relationships between the Asian-Pacific Oscillation and summer climate anomalies over Asian, North Pacific and North America during recent 100 years. Journal of Climate, 24(18): 4793-4799. DOI:10.1175/JCLI-D-11-00054.1
Zhao P, Yang S, Jian M Q, et al. 2011b. Relative controls of Asian-Pacific summer climate by Asian land and tropical-North Pacific sea surface temperature. Journal of Climate, 24(15): 4165-4188. DOI:10.1175/2011JCLI3915.1
Zhao P, Wang B, Zhou X J. 2012. Boreal summer continental monsoon rainfall and hydroclimate anomalies associated with the Asian-Pacific Oscillation. Climate Dynamics, 39(5): 1197-1207. DOI:10.1007/s00382-012-1348-6
Zhou B T, Cui X, Zhao P. 2008. Relationship between the Asian-Pacific Oscillation and the tropical cyclone frequency in the Western North Pacific. Science in China Series D:Earth Sciences, 51(3): 380-385. DOI:10.1007/s11430-008-0014-7
Zhou B T, Zhao P. 2010. Influence of the Asian-Pacific Oscillation on spring precipitation over central eastern China. Advances in Atmospheric Sciences, 27(3): 575-582. DOI:10.1007/s00376-009-9058-7
Zhou T J, Zhang J. 2009. Harmonious inter-decadal changes of July-August upper tropospheric temperature across the North Atlantic, Eurasian continent and North Pacific. Advances in Atmospheric Sciences, 26(4): 656-665. DOI:10.1007/s00376-009-9020-8
Zou Y, Zhao P. 2009. Relation of summer Asian-Pacific oscillation to tropical cyclone activities over the coastal waters of China. Acta Meteorologica Sinica, 67(5): 708-715.
陈晓龙, 周天军, 邹立维. 2013. 夏季亚洲-太平洋涛动的耦合模式模拟. 气象学报, 71(1): 23–37.
黄荣辉, 李维京. 1988. 夏季热带西太平洋上空的热源异常对东亚上空副热带高压的影响及其物理机制. 大气科学, 12(S1): 107–116. DOI:10.3878/j.issn.1006-9895.1988.t1.08
黄荣辉. 1990. 引起我国夏季旱涝的东亚大气环流异常遥相关及其物理机制的研究. 大气科学, 14(1): 108–117.
黄荣辉, 孙凤英. 1992. 北半球夏季遥相关型的年际变化及其数值模拟. 大气科学, 16(1): 52–61.
姜燕敏, 黄安宁, 吴昊旻. 2015. 不同水平分辨率BCC_CSM模式对中亚地面气温模拟能力评估. 大气科学, 39(3): 535–547. DOI:10.3878/j.issn.1006-9895.1408.14133
汪栩加, 郑志海, 封国林, 等. 2015. BCC_CSM模式夏季关键区海温回报评估. 大气科学, 39(2): 271–288. DOI:10.3878/j.issn.1006-9895.1408.13329
吴统文, 宋连春, 李伟平, 等. 2014. 北京气候中心气候系统模式研发进展——在气候变化研究中的应用. 气象学报, 72(1): 12–29. DOI:10.11676/qxxb2013.084
辛晓歌, 吴统文, 张洁. 2012. BCC气候系统模式开展的CMIP5试验介绍. 气候变化研究进展, 8(5): 378–382.
张芳, 董敏, 吴统文. 2014. CMIP5模式对ENSO现象的模拟能力评估. 气象学报, 72(1): 30–48. DOI:10.11676/qxxb2014.011
张莉, 吴统文, 辛晓歌, 等. 2013. BCC_CSM模式对热带降水年循环模态的模拟. 大气科学, 37(5): 994–1012. DOI:10.3878/j.issn.1006-9895.2012.12126
赵平, 陈军明, 肖栋, 等. 2008. 夏季亚洲-太平洋涛动与大气环流和季风降水. 气象学报, 66(5): 716–729. DOI:10.11676/qxxb2008.066
周波涛, 崔绚, 赵平. 2008. 亚洲-太平洋涛动与西北太平洋热带气旋频数的关系. 中国科学D辑:地球科学, 38(1): 118–123.
邹燕, 赵平. 2009. 夏季亚洲-太平洋涛动与中国近海热带气旋活动的关系. 气象学报, 67(5): 708–715. DOI:10.11676/qxxb2009.071