于怡秋,李忠贤,杨宝钢,魏麟骁,陈晨

① 中国气象局 沈阳大气环境研究所,辽宁 沈阳 110166;② 南京信息工程大学 气象灾害教育部重点实验室/气候与环境变化国际联合实验室/气象灾害预报预警与评估协同创新中心,江苏 南京 210044;③ 沈阳区域气候中心,辽宁 沈阳 110166;④ 重庆市气候中心,重庆 401147;⑤ 齐齐哈尔市气象局,黑龙江 齐齐哈尔 161006

北大西洋海表温度存在明显的多时间尺度变化(battisti et al.,1995;molinari et al.,1997;sutton and allen,1997;czaja and marshall,2001;wu and liu,2005)。在季节和年际尺度上,北大西洋三极型海温异常是北大西洋变率的主要模态。北大西洋三极型海温异常的形成很大程度上受北大西洋涛动(north atlantic oscillation,nao)影响(deser and timlin,1997;李忠贤等,2019)。三极型海温异常反过来可以对nao形成一个正反馈过程,在观测数据分析和海气耦合模式模拟研究中均得到了验证(deser and blackmon,1993;kushnir,1994;czaja and frankignoul,1999,2002;watanabe and kimoto,2000;drévillon et al.,2001;pan,2005)。

不同季节的北大西洋海温对大气环流和气候异常均有重要影响。冬季北大西洋海温异常激发欧亚遥相关型波列,从而影响冬季西伯利亚高压、乌拉尔山阻塞高压和东亚大槽,造成我国冬季气温异常(宗海锋等,2008;李崇银和顾薇,2010;liu et al.,2014;李栋梁和蓝柳茹,2017)。冬季北大西洋副极地的海温上升会激发欧亚大陆遥相关波列,造成乌拉尔山北部地区的冬春雪盖深度减小,引起贝加尔湖北侧高度场的升高,从而导致蒙古气旋的减少(黄鑫等,2017)。冬季北大西洋海温异常可通过激发下游异常波列来影响中东急流的变化(李忠贤等,2013)。冬季北大西洋海温异常可导致春季中国东北地区温度异常(张霏燕和徐海明,2011)。前期冬春北大西洋海温异常可以激发下游波列,影响乌拉尔山高压和鄂霍次克海阻塞高压强弱,从而影响江淮入梅时间(徐海明等,2001;gu et al.,2009)。春季北大西洋三极型海温异常对欧亚地面气温具有显著的影响(chen et al.,2016,2017)。cui et al.(2015)研究指出,春季北大西洋三极型海温异常激发出来的波列能影响春季青藏高原上空副热带西风急流的强度。wu et al.(2009)研究认为,从春季持续到夏季的北大西洋三极型海温有利于欧亚遥相关向下游发展,所激发的波列呈现出沿副热带西风急流向极一侧分布,认为北大西洋三极型海温是联系春季北大西洋涛动和东亚夏季风的重要桥梁。春季北大西洋三极型海温异常与我国东北夏季气温关系密切,是我国东北夏季气温的前兆信号(wu et al.,2011)。zuo et al.(2013)通过观测分析和数值试验发现,夏季北大西洋三极型海温异常可以引起东亚夏季风强弱变化,并影响长江中下游降水(wu et al.,2009;zuo et al.,2013)。此外,夏季北大西洋海温还可以通过影响热带纬向环流,从而对西北太平洋对流层低层异常反气旋、西太平洋副热带高压和亚洲副热带西风急流起调制作用(容新尧等,2010;hong et al.,2014;何超等,2015;霍利微等,2016;任宏昌等,2017)。

综上所述,前人的研究表明北大西洋海温异常对欧亚大气环流和东亚气候具有重要的影响。我国气温变化复杂,具有明显的季节性和区域性,高温干旱和低温冷冻等灾害天气可能会造成不同地区的经济损失和较大的社会影响,因此提高我国气温预测的准确性十分必要(陈丽华等,2010)。近年来,不少学者研究指出,北大西洋海温异常对我国气温变化具有重要的影响(wu et al.,2011;张霏燕和徐海明,2011;孙建奇,2014)。然而,以往的研究侧重于分析北大西洋某一特定海区对我国气温的影响,针对夏季北大西洋三极型海温异常与同期我国气温联系的研究相对较少。因此,本论文初步探究夏季北大西洋三极型海温与我国同期气温在年际尺度上的联系,旨在进一步提升我国短期气候预测能力。

1)英国气象局hadley中心提供的逐月海温资料(rayner et al.,2003),水平分辨率为1°×1°,时间段为1951—2016年;2)美国ncep/ncar(national centers for environmental prediction/national center for atmospheric research)提供的再分析资料(kalnay et al.,1996),包括高度场、风场等,水平分辨率为2.5°×2.5°,垂直方向为17层,时间段为1951—2016年;3)中国气象局国家气候中心提供的全国160站逐月气温资料,时间段为1951—2016年。

所用到的方法主要包括:经验正交函数分解(eof)、相关分析和回归分析等。夏季为6—8月平均(jja)。本文的研究主要关注年际变化,因此所有资料均去除9 a以上的年代际变化和线性趋势。文中所涉及的海温场、气温场和大气环流场均为年际变化分量,本文在对年际分量进行相关和回归分析时,应用施能(2002)给出的方法计算有效自由度,然后根据有效自由度进行置信度检验。

对1951—2016年共66 a的夏季北大西洋海温距平场进行eof分解,第一模态在经向上呈现出“-+-”的三极结构(图略),其方差贡献率是19.5%。为方便讨论,本文将北大西洋海温距平从热带至副极地呈现出“-+-”空间型结构称为三极型海温异常正位相结构,“+-+”空间型结构称为三极型海温异常负位相结构。将eof分解得到的第一模态时间系数标准化值定义为夏季北大西洋三极型海温异常(nat,north atlantic sea surface temperature tripole)指数。

研究表明,北大西洋海温异常和中国夏季气候异常受到enso的影响(klein et al.,1999;朱益民等,2007;wu and zhang,2010;wu et al.,2011)。因此,在分析夏季北大西洋海温年际异常与我国气温年际变化的联系时,去除夏季北大西洋三极型海温异常指数、气温和环流场中的前期冬季enso信号。

图1为1951—2016年夏季nat指数与我国同期气温年际异常的回归系数分布。如图1所示,当夏季nat正(负)位相年,我国长江中上游气温异常偏低(高),尤其是西南地区,这一异常关系通过了95%置信度的显著性水平检验。当夏季nat变化一个标准差时,我国西南地区气温异常变化0.15 ℃左右。由此可知,夏季nat与我国同期西南地区气温年际异常关系密切。

图1 夏季nat指数与中国同期气温年际异常的回归系数(单位：℃；阴影区表示通过了95%置信度的显著性水平检验)fig.1 simultaneous regression coefficients of interannual temperature anomalies in china onto nat index in summer (units:℃;shaded areas indicate the coefficients statistically significant at 95% confidence level)

图2 夏季swct指数与同期北大西洋海温年际异常的相关系数(阴影区表示通过了95%置信度的显著性水平检验)fig.2 simultaneous correlation coefficients between swct index and interannual north atlantic sst anomalies in summer (shaded areas indicate the coefficients statistically significant at 95% confidence level)

图3 夏季nat指数(实线)和swct指数(虚线)的标准化曲线(均去除前期enso信号)fig.3 normalized curves of nat index (solid line) and swct index (dashed line) in summer (the previous enso signal are both removed)

将我国西南地区(98°～108°e,22°～30°n)平均的气温年际异常的标准化值定义为年际尺度的西南地区气温(southwest china temperature,swct)指数。由图2可知,我国夏季西南地区气温年际异常与同期北大西洋海温年际异常相关系数大致呈现出由北至南“+-+”的三极结构,与夏季nat负位相的结构类似,相关系数高值区主要分布在北大西洋中部。图3给出夏季nat指数与swct指数的标准化曲线,可以发现二者表现出较好的反向变化特征,计算表明,夏季nat指数与swct指数的相关系数为-0.42(表1),通过95%置信度的显著性水平检验。

表1 nat、swct、beat指数之间的相关系数

为了揭示夏季nat与我国同期气温年际异常联系的可能物理过程,将夏季nat指数与同期位势高度和波活动通量年际异常进行回归分析(图4)。由图4可知,在200 hpa位势高度上,从北大西洋至东亚中高纬地区呈现出一个“-+-+-”波列结构,负异常中心分别位于冰岛、乌拉尔山和东亚上空,正异常中心分别位于西欧和贝加尔湖上空,在500 hpa和700 hpa位势高度上呈现出与200 hpa类似的波列结构,北半球中高纬对流层各层大气环流场呈现出相当正压结构。这意味着夏季nat正(负)位相可激发出从北大西洋至东亚中高纬的遥相关波列,东亚中纬度地区上空为高度场负(正)距平。波活动通量以50°w附近的北大西洋为起点,沿北半球中高纬向下游传播至东亚地区,传播路径与正负位势高度中心相符。能量传播至东亚地区,进而引起东亚中纬度地区位势高度降低,产生异常气旋性环流。rossby的传播路径与波列的路径一致,进一步证明了欧亚中高纬的遥相关波列是由夏季nat异常所激发出的。

图4 夏季nat指数对同期位势高度(阴影区；单位：gpm)和波活动通量(箭矢；单位：m2/s2)年际异常的回归系数(点阴影区表示通过了95%置信度的显著性水平检验；仅给出大于等于0.01 m2/s2的波活动通量)：(a)200 hpa；(b)500 hpa；(c)700 hpafig.4 simultaneous regression coefficients of interannual geopotential height anomalies (shaded areas;units:gpm) and interannual wave activity flux anomalies (arrows;units:m2/s2) onto nat index in summer (stippled shadings indicate the coefficients statistically significant at 95% confidence level,and only wave activity fluxes greater than or equal to 0.01 m2/s2 are given):(a)200 hpa;(b)500 hpa;(c)700 hpa

图5为夏季nat指数与500 hpa、700 hpa旋转风场年际异常的回归系数。由图5可知,夏季nat激发的波列在贝加尔湖地区表现为反气旋性环流异常,东亚地区表现为气旋性环流异常。在500 hpa和700h pa位势高度上,夏季nat正位相年,西南地区在对流层中下层位于贝加尔湖异常反气旋性环流的东部和东亚异常气旋性环流的西部,受异常偏北气流的影响,偏北风将北方的冷空气向南输送,有利于我国西南地区气温降低;反之亦然。

图5 夏季nat指数对同期旋转风年际异常的回归系数(箭矢，单位：m/s；阴影区表示通过了95%置信度的显著性水平检验)：(a)500 hpa；(b)700 hpafig.5 simultaneous regression coefficients of interannual rotational wind anomalies onto nat index in summer (arrows,units:m/s;shaded areas indicate the coefficients statistically significant at 95% confidence level):(a)500 hpa;(b)700 hpa

上述分析表明,与夏季nat所激发的北半球欧亚中高纬遥相关波列相关的贝加尔湖反气旋性环流异常和东亚气旋性环流异常在联系北大西洋海温与我国西南地区气温年际变化中起到了重要作用。为了进一步反映贝加尔湖反气旋性环流异常和东亚气旋性环流异常与北大西洋三极型海温及我国西南地区夏季气温年际变化的联系,本文定义了一个贝加尔湖-东亚遥相关(baikal-east asia teleconnection pattern,beat)指数,具体为贝加尔湖附近(100°～120°e,50°～60°n)与东亚地区(110°～135°e,30°～40°n)的500 hpa位势高度区域平均值之差,并将平均值之差再进行一次标准化处理。

图6为夏季beat指数和swct指数与500 hpa位势高度场年际异常的回归系数。由6a可知,在北半球中高纬地区,从北大西洋到东亚上空基本上呈现出“-+-+-”的遥相关波列,与图4b类似。由此可见,贝加尔湖-东亚遥相关指数与北大西洋三极型海温异常激发的遥相关波列具有很好的联系。由图6b可知,夏季swct指数回归的500 hpa位势高度年际异常场在贝加尔湖上空为位势高度负异常,东亚上空为位势高度正异常,与贝加尔湖-东亚遥相关指数回归的位势高度场基本上呈反向结构。夏季beat指数与我国同期气温年际异常的相关分布场(图7)可知,夏季beat指数与长江流域大部分地区的气温异常存在负相关关系。

图6 夏季beat指数(a)和夏季swct指数(b)对同期500 hpa位势高度年际异常的回归系数(单位：gpm；阴影区表示通过了95%置信度的显著性水平检验)fig.6 simultaneous regression coefficients of interannual geopotential height anomalies onto (a)beat index and (b)swct index in summer (units:gpm;shaded areas indicate the coefficients statistically significant at 95% confidence level)

图7 夏季beat指数与中国同期气温年际异常的相关系数(阴影区表示通过了95%置信度的显著性水平检验)fig.7 simultaneous correlation coefficients between beat index and interannual temperature anomalies in china in summer (shaded areas indicate the coefficients statistically significant at 95% confidence level)

计算nat、swct和beat指数之间的相关系数(表1),进一步分析了它们之间的相关关系。beat指数与nat和swct指数的相关系数分别为0.55和-0.49,均通过了99%置信度的显著性水平检验,表明三者是密切相关的。这进一步证明了夏季nat可能通过贝加尔湖-东亚遥相关型来影响我国西南地区气温年际变化。

夏季东亚副热带西风急流是影响我国天气和气候的重要环流系统(况雪源和张耀存,2006)。图8为夏季nat指数和swct指数与200 hpa纬向风场年际异常的回归系数。由图8a可以看出,当夏季nat正位相时,在东亚地区,40°～55°n之间纬向风负异常,15°～40°n之间纬向风正异常,呈现出南北反向变化的结构,对应东亚副热带西风急流位置偏南,反映出西风急流南北位置的移动。在850 hpa到150 hpa之间,长江中上游地区有辐合上升运动(图略),有利于西南地区的低温的发生;反之亦然。由图8b可知,我国西南气温偏高时,在东亚地区,40°～50°n纬向风正异常,25°～40°n纬向风负异常,也表现为东亚副热带西风急流南北位置的变化,对应西风急流位置偏北。可见,夏季nat也可以通过影响夏季东亚副热带西风急流的南北位置变化,进而对我国西南地区气温变化产生影响。

图8 夏季nat指数(a)和swct指数(b)对同期200 hpa纬向风年际异常的回归系数(单位：m/s；阴影区表示通过了95%置信度的显著性水平检验)fig.8 simultaneous regression coefficients of interannual zonal wind anomalies at 200 hpa onto (a)nat index and (b)swct index in summer (units:m/s;shaded areas indicate the coefficients statistically significant at 95% confidence level)

本文利用hadley中心的海温资料、ncep/ncar再分析资料及全国160站逐月气温资料,探讨了夏季北大西洋三极型海温异常与我国气温年际变化的联系及可能机制,得到以下主要结论:

1)夏季北大西洋三极型海温异常与我国夏季西南地区气温年际变化为显著的负相关关系,两者之间的相关系数为-0.42,通过了99%置信度的显著性水平检验。当夏季北大西洋三极型海温异常处于正(负)位相,我国夏季西南地区气温异常偏低(高)。

2)夏季北大西洋三极型海温正位相异常激发的北半球中高纬欧亚遥相关波列引起贝加尔湖地区位势高度升高,产生反气旋性环流异常;东亚地区位势高度降低,产生气旋性环流异常。西南地区在对流层中下层位于贝加尔湖异常反气旋性环流的东部和东亚异常气旋性环流的西部,在两者的共同作用下,受异常偏北气流的影响,有利于冷空气南下和堆积。此外,北大西洋三极型海温异常也可以通过影响东亚副热带西风急流来影响我国西南地区气温。当夏季北大西洋三极型海温异常正位相时,对流层高层的东亚副热带西风急流位置偏南,西南地区位于急流南侧的异常上升气流中,此种环流形势配置将导致该地区气温降低;反之亦然。