MJO与北半球高纬地区冬季地表气温异常的联系*

崔 静 杨双艳 LI Tim,2

1 南京信息工程大学气象灾害教育部重点实验室/气候与环境变化国际合作联合实验室/气象灾害预报预警与评估协同创新中心，南京 210044 2 夏威夷大学国际太平洋研究中心，夏威夷 美国

提 要： 基于1979—2016年NCEP-NCAR逐日再分析资料研究了热带季节内振荡(MJO)和北半球冬季高纬地区地表气温(SAT)之间的联系。利用实时多变量MJO(RMM)指数，将MJO分为8个位相，其中位相2(位相6)对应于位于印度洋地区的正(负)对流。不同MJO位相下的SAT合成结果显示MJO第二位相后的5～15 d，北半球高纬地区(60°～90°N、180°～60°W)有明显的负SAT异常；由于热带异常加热信号的改变，在MJO第六位相后的5～15 d该地区则对应于显著的正SAT异常。该地区温度的垂直结构在各个位相下也表现出类似的分布特征。合成的500 hPa位势高度异常场显示，在温度负(正)异常的位相对应有明显的位势高度负(正)异常，这种环流异常主要是由与热带对流异常相联系的向东北方向传播的罗斯贝波列所引起的。通过对波活动通量的计算，推断该东北方向传播的罗斯贝波列很可能是罗斯贝波能量频散的结果。合成的700 hPa比湿异常场和SAT之间在存在着较好的对应关系，考虑到对流层中层的比湿与向下长波辐射之间存在着正相关关系，说明该温度异常也可能与辐射过程相关。上述分析表明与MJO对流相关的大尺度环流异常对高纬地区季节内SAT变率有重要影响，该异常SAT信号可能来自平流输送和辐射过程等。准确把握MJO位相与SAT异常信号的联系也可以为北半球高纬地区SAT的延伸期预报提供一些可靠线索。

引 言

热带大气季节内振荡(Madden-Julian oscillation，MJO)是热带地区大气环流异常的一个最重要的模态，主要的时间尺度为30～60 d，主要特征表现为沿全球热带地区东传的对流活动，具有纬向1波为主的行星尺度空间结构(Madden and Julian，1971；1972)。MJO活动表现出明显的季节演变特征。在北半球冬季，MJO以沿赤道向东的传播为主；在北半球夏季，向东传播会减弱，在印度夏季风区以及东亚季风区以北传的对流活动为主导，在西北太平洋区以西北向的传播为主导。Li(2014)指出北半球夏季热赤道的北移是造成北印度洋和中国南海地区MJO对流北向传播的重要原因之一，同时背景场中的东风切变和水汽的经向分布也起到了重要贡献(Jiang et al，2004)。关于MJO动力机制(包括其触发机制和东传机制)的研究工作已取得了显著的进展(Wang and Rui，1990；Wang and Li，1994；Hsu and Li，2012；Sobel and Maloney，2013；Zhao et al，2013)。

MJO作为热带最主要的环流演变模态可以直接调节热带地区的天气气候(Madden and Julian，1994)。比如，它会对热带气旋和飓风的生成和加强(Maloney and Hartmann，2000；Hall et al，2001；祝从文等，2004；李崇银等，2012)产生影响。MJO还会影响印度洋偶极子事件(Rao and Yamagata，2004)、北美环流异常(Higgins and Mo，1997)、南海夏季风(林爱兰等，2016；李春晖等，2016)、孟加拉湾西南夏季风季节内振荡(李汀和琚建华，2013)以及ENSO事件的发展(Zhang and Gottschalck，2002)。

前人的许多研究表明，MJO对环流异常的影响不仅仅局限于热带，由于对流加热异常可以激发罗斯贝波列(Ferranti et al，1990；Kiladis and Weickmann，1992)，这种影响可以延伸至中高纬地区(Lau and Phillips,1986)。比如，Bond and Vecchi(2003)指出了MJO对于俄勒冈州和华盛顿早冬和晚冬降水的重要影响；Lin and Brunet(2009)表明MJO位相与加拿大冬季地表气温(SAT)之间存在显著的相关关系。除热带季节内振荡以外，中高纬的季节内振荡(ISO)同样会对热带和中高纬的天气气候产生重要影响(Wen et al，2011；许乐心等，2017；Yang and Li,2016;2017;Yang et al，2019；李健颖和毛江玉，2019)。近年来，针对MJO与我国气候异常间的联系也开展了一系列研究，研究发现MJO对我国降水有重要的调制作用(吴捷等，2018)，对中国东部地区夏季(白旭旭等，2011)和冬季(刘冬晴和杨修群，2010)、华南前汛期(李文铠等，2014；章丽娜等，2011)、夏季西南(李永华等，2016)、云南地区(牛法宝等，2013；李汀等，2012)以及夏季广东(林爱兰等，2013)的降水异常均存在一定程度的影响。同时，MJO对我国气温异常也存在显著影响，比如2008年初南方持续性冰冻雨雪天气(吴俊杰等，2009；马宁等，2011)以及2009年11月我国东部大范围雨雪天气(贾小龙和梁潇云，2011)等极端天气都与MJO存在着密切的联系。因此，MJO的发展和演变，对于把握这种气候异常有很大的帮助。MJO作为次季节时间尺度上最具可预报性的模态(Gottschalck et al，2010)，是延伸期(10～30 d)预报主要的可预报来源(丁一汇和梁萍，2010；贾小龙等，2012)。MJO和大气低频季节内振荡的准确预报对于填补短期天气预报和长期气候预测之间的空隙至关重要。

在Vecchi and Bond(2004)的工作中分析了MJO对高纬地区冬季SAT的影响，他们的合成分析均为基于MJO各个位相同时次的合成。前人的研究表明中高纬地区对于热带强迫的响应在时间上存在一定的滞后，完全发展建立需要两周左右的时间(Jin and Hoskins，1995)。为使这种异常SAT信号更加显著，我们的工作主要关注了MJO与北半球冬季高纬地区SAT之间的滞后相关关系，从环流分析等角度进一步揭示两者之间的联系。旨在将MJO和高纬地区的大气变率相联系，从而为该地区季节内SAT变率的延伸期预报提供参考。

1 资料和方法

1.1 资 料

本文使用了美国国家海洋与大气中心(NOAA)的逐日向外长波辐射(OLR)资料以表征热带大尺度对流特征；同时使用了美国国家环境预测中心-国家大气研究中心(NCEP-NCAR)发布的全球逐日再分析资料(水平分辨率为2.5°×2.5°)。本文的分析时段为1979/1980年到2015/2016年共37个冬季，冬季定义为11月到次年3月共5个月。再分析资料中使用到的气象要素包括2 m处的地表气温(T2 m)以及在各个高度层次的纬向风(u)、经向风(v)、温度(T)、p坐标系下的垂直速度(ω)等。

采用由Wheeler and Hendon(2004)所定义的实时多变量MJO指数(real-time multivariate MJO indices，简称RMM指数)来描述MJO的活动特征，RMM指数最常用于MJO的监测和预报。本文使用到的逐日RMM指数可以在澳大利亚气象局官网下载(http:∥www.bom.gov.au/climate/mjo/)。该指数是将逐日资料投影到热带地区(15°S～15°N)纬向平均后的逐日的850 hPa纬向风(u850)、200 hPa 纬向风(u200)和OLR的前两个多变量EOF模态上，得到的前两个主成分时间序列，将该时间序列标准化即可得到RMM指数，分别记为RMM1和RMM2。上述资料均去除了年循环和年际变率。根据RMM1和RMM2所确定的二维位相空间，可以将强的MJO 沿全球热带地区的纬向传播分为8个空间位相。MJO的地理位置可以由一个以RMM1和RMM2指数所定义的二维位相空间中的一个点来表示，而该点与整个空间中心的距离则表示MJO的振幅，代表MJO事件的强度，这8个位相分别对应于MJO生命史的特定阶段，表示MJO对流从印度洋到太平洋的向东传播。

1.2 方 法

为提取与MJO相关的季节内信号，对原始数据进行了相关处理。首先去除了1979—2016年多年气候态的0～3波以去除季节循环，然后通过使用Lanczos带通滤波方法(Duchon，1979)进行30～60 d 的带通滤波，滤波后选择了11月至次年3月的资料用于分析，但1979年1月至次年3月以及2016年11—12月的资料未包括在内。

为检验每个MJO位相合成结果是否显著，使用到的t检验公式为：

式中：ρ1为滞后一天的相关系数，m2是yi的样本数。

2 MJO位相合成的高纬地表气温异常

2.1 MJO位相

为了证实赤道对流活动是和MJO的位相相联系的，图1给出了各个位相下合成的OLR场及850 hPa的经向风和纬向风异常场。在位相1，MJO对流中心位于非洲和西印度洋，此时太平洋以西风异常为主，而印度洋存在东风异常；在位相2和位相3，对流中心主要位于印度洋，继续发展并逐渐东移至海洋性大陆(位相4、位相5)；位相6～位相8，对流中心由西太平洋继续东移至中太平洋。

使用到的RMM指数时间段包括从1979/1980—2015/2016年的37个冬季，因此可以得到每个位相所对应的天数。图2给出了每个MJO位相所对应的天数及平均振幅的大小。可以看到位相2～位相4以及位相6～位相8 MJO事件的发生频率较高，其中平均振幅的最大值出现位相3。

图1 MJO 8个位相(a～h)下合成的30～60 d滤波后的冬季OLR异常场(填色，单位：W·m-2)和850 hPa异常风场(箭头，单位：m·s-1)(阴影区为通过0.05显著性水平检验区域，t检验；右下角数字为每个位相对应的天数)

图2 二维MJO指数(RMM1和RMM2)定义的MJO空间位相以及各位相对应的天数和平均振幅

2.2 与MJO相联系的高纬地区地表气温异常

为探究北半球冬季SAT和MJO之间的联系，对不同MJO位相下的SAT进行了合成分析。在北半球高纬地区这种SAT异常信号较为显著。Jinand Hoskins(1995)研究表明：给定的热带加热响应下的罗斯贝波列会在15 d以内完全建立。因此在下文的合成分析中将分别采取0，5，10，15 d的滞后。

MJO的位相1～位相4代表其一半的生命周期，而位相5～位相8则代表相反的另一半周期，这种生命史的对称特征也表现在SAT的合成场上，随着滞后时间的变化，每个位相下高纬SAT异常的强度也会有一定的变化，总体来说在滞后5～15 d时位相2～位相7下高纬地区SAT信号比较显著，第0天时(即同时次的合成图中)这种异常SAT信号相对较弱，随着滞后时间的增加，该异常信号得到加强。这里给出了滞后5～15 d的MJO的8个位相下的异常SAT合成图(图3)，可以看到位相2和位相6以及位相3和位相7下这种SAT信号比较显著，并且该异常SAT信号存在着明显的反向对应关系(即位相2、6表现为明显的负异常而位相3、位相7则表现为明显的正异常)，后续为了简洁起见选择位相2和位相6进行研究(对位相3和位相7也做了同样分析，得到相似结论)。从图4中可以看到在MJO的第2(第6)位相后5 d左右北半球高纬地区开始出现负(正)的SAT异常，在滞后10～15 d时这种异常信号表现得更加明显。

此外，值得关注的是在高纬地区(60°N以北)，其SAT的异常情况并不是完全一致的，部分区域会呈现出相反的变化，在滞后5 d以后其SAT异常信号均较为显著和一致的主要为60°～90°N、180°～60°W这部分区域。为了探究该区域温度的垂直结构是否也有类似的变化特征，图5给出了MJO位相2和位相6滞后0，5，10，15 d后的温度异常及位势高度异常沿纬向平均(180°～60°W)后的经向-垂直分布。和地表气温类似，整个对流层的温度异常也呈现出相似的分布特征，即在滞后5～10 d时，位相2(位相6)呈现出明显的负(正)温度异常。另一个显著的特征是：例如在位相2后的5 d，一个位于300 hPa左右的负位势高度异常中心出现在70°～80°N，在该中心以上(以下)温度异常表现为暖(冷)异常，较好地满足了静力平衡关系。在第6位相也能看到这种较好的对应关系，静力平衡关系使得低层的与MJO相关的温度异常信号与对流层高层的位势高度异常之间建立了密切的关系。最后，由于在北半球高纬地区冬季有相似幅度的增暖(冷却)发生在MJO第1，3，4(5，7，8)位相后的1～2周，所以该地区冬季SAT与MJO位相2和位相6之间的联系同样应该可适用于MJO的其他位相。

图3 MJO 8个位相(a～h)下冬季SAT滞后5～15 d的合成(打点代表通过0.05的显著性水平检验)

图4 MJO位相2(a)和位相6(b)下冬季SAT滞后0 d(a1,b1)，5 d(a2,b2)，10 d(a3,b3)，15 d(a4,b4)的合成(打点代表通过0.05的显著性水平检验)

图5 MJO位相2(a)和位相6(b)下滞后0 d(a1,b1)，5 d(a2,b2)，10 d(a3,b3)，15 d(a4,b4)的冬季温度异常场(填色)以及位势高度异常场(等值线，单位：gpm)沿180°～60°W平均后的经向-垂直分布(打点代表通过0.05显著性水平检验)

3 大尺度环流异常

为了进一步理解上述北半球高纬地区与MJO相联系的SAT信号，探究了与MJO位相相联系的大尺度异常环流场的分布特征。Ferranti et al(1990)，Matthews et al(2004)研究表明：与MJO相关的热带非绝热加热异常会激发罗斯贝波列并影响中高纬的大尺度环流异常。

图6给出了在冬季MJO位相2和位相6下合成的北半球500 hPa位势高度异常场。对于MJO位相2，在滞后5 d时，一个显著的特征是有一个位于北太平洋，中心位于40°N、180°附近的正位势高度异常中心。在阿拉斯加和加拿大西岸有一个负位势高度异常中心，加拿大东部则为正位势高度异常，可以看到波的模态主要表现出向下游频散的特征，与PNA的模态很类似，但是相对于PNA有略微的西移；到滞后10 d时，北大西洋分支有明显的环流异常在发展，形成了近似于正的NAO的模态。位相6的合成图表现出和位相2相似的演变特征，但由于此时的热带强迫发生了变化，因此呈现出的波列模态也相反。

由图6可以看到北半球高纬表现出负(正)SAT异常的地区主要对应于负(正)位势高度异常，这表明MJO对北半球高纬地区SAT信号的影响可能与水平温度平流的输送相关。这种位势高度的异常在对流层整层基本是一致的，呈现出正压的结构。参考前人的工作，可以推测与MJO相关的热带对流所引起的加热异常可能会触发罗斯贝波列，由此持续性地影响了北半球高纬地区的环流异常，从而对季节内SAT信号产生重要影响。

图6 同图5，但为合成的冬季500 hPa位势高度异常场(等值线，间隔为2 gpm)(阴影为通过0.05的显著性水平检验)

为探究热力因素与异常SAT信号之间是否也存在着一定的对应关系，对MJO位相下合成的700 hPa比湿场异常场(图8)进行了简单分析。可以看到MJO位相2(位相6)下较冷(暖)的SAT分别对应于700 hPa较干(湿)的空气，考虑到对流层中层比湿和向下长波辐射之间又存在正相关关系(Overland et al，1997; Adams et al，2000)，可以推测不同MJO位相下的比湿正(负)异常会导致向下长波辐射的正(负)异常，从而会对该地区的正(负)SAT信号有一定的正贡献。

由上述分析发现与MJO对流相联系热带加热异常会触发东北向传播的罗斯贝波列，该波列会影响到研究所关注的高纬地区的大尺度环流异常，从而对该地区的季节内SAT异常信号产生重要影响。此外，除异常环流场以外，700 hPa比湿场和高纬SAT信号间也有较好的对应关系，表明与辐射过程相关的热力因素可能对该地区的SAT异常有一定的正贡献。而具体哪一过程在温度异常的变化中起主导作用还需借助于温度方程进行诊断分析，在之后的工作中我们会作进一步的研究。

图7 MJO位相2(a)和位相6(b)合成的200 hPa波活动通量(箭头，单位： m2·s-2)及罗斯贝波源(填色，单位：10-10 s-2)

图8 同图5，但为合成的冬季SAT(填色,仅显示通过0.05显著性水平检验的部分)和700 hPa比湿异常场(等值线，间隔为2×10-5 kg·kg-1)

4 结论和讨论

本文研究了MJO和北半球冬季高纬地区SAT之间的联系，基于38年(1979—2016年)的逐日2 m地表气温资料以及RMM指数，得到了MJO热带对流和北半球冬季高纬地区SAT异常之间显著的相关关系。全文主要关注了MJO生命周期中的两个位相：位相2和位相6，分别对应于位于东印度洋和西太平洋地区的MJO正对流。考虑到中高纬对于热带对流的响应存在1～2周的延迟，为使与MJO相联系的异常信号更显著，在进行MJO位相的合成分析时考虑了5～15 d的滞后。可以明显看到:在北半球高纬地区，尤其是60°～90°N、180°～60°W这部分区域，在MJO位于位相2后的5～15 d有明显的负SAT异常，此时的正对流主要位于印度洋和海洋性大陆地区，负对流位于太平洋西部及中部；同样的，在MJO位于位相6后的5～15 d，该区域有显著的正SAT异常，对应位于印度洋和海洋性大陆的负对流以及位于太平洋西部及中部的正对流。该地区的温度异常的垂直分布特征与SAT类似，即在MJO位于位相2后的5～15 d表现为负温度异常，而在MJO位于位相6后的5～15 d则表现为正温度异常。滞后合成的大尺度异常环流场表明热带对流活动会对中高纬的SAT产生影响，通常是以罗斯贝波列的形式。500 hPa位势高度异常的合成场上可以明显看到在太平洋以及北美地区形成了类似于PNA模态的波列。罗斯贝波源及波活动通量的分析表明该波列应该是由与MJO对流相联系的热带加热异常所触发的，同时在加热场东北侧的正波源区有东北指向的波活动通量，表明观测的罗斯贝波列可能是东北方向频散的罗斯贝能量的结果。热带对流通过这种遥相关的方式对中高纬的环流形势产生影响，从而影响季节内SAT变率。此外700 hPa的比湿异常场的合成分析图也呈现出与SAT异常信号间较好的对应关系，在高纬地区温度负异常时有明显的比湿负异常，相反的，在温度正异常时则有明显的比湿正异常，而对流层中层的比湿和向下长波辐射之间存在正相关关系，表明辐射过程等热力因素也可能会对该地区的SAT异常有一定的贡献。

本文中MJO位相和SAT异常间显著的相关关系可能会给北半球高纬地区SAT的预报提供一些帮助，因为热带地区缓慢变化的对流活动可以为预报提供季节内时间尺度上的信息。但是也认识到了本次工作中的几个局限，首先，上文所讨论的滞后合成分析仅仅是多个MJO事件的平均结果，而这些MJO事件通常具有不同的传播特征和强度。此外，北极的SAT异常更多是受到高纬地区自身的变化的影响，比如西伯利亚高压和北极涛动(AO)等，所以MJO的影响只能被视作是一个附加部分并且它的影响也是有限的。最后，还应该注意到，研究不可避免地会受到和厄尔尼诺等相关的更长时间尺度变率的影响。热带MJO是进行10～30 d延伸期预报的重要可预报来源，如果能准确预测MJO的发展和演变，对于气候异常的预测可以提供一些依据。