大气科学  2015, Vol. 39 Issue (1): 210-220   PDF    
中国大范围持续性极端低温事件的一类平流层前兆信号
施宁1, 布和朝鲁2,3    
1 南京信息工程大学气象灾害预报预警与评估协同创新中心, 气象灾害教育部重点实验室, 南京210044;
2 中国科学院大气物理研究所国际气候与环境科学中心, 北京100029;
3 中国科学院大气物理研究所大气科学和地球流体力学数值模拟国家重点实验室, 北京100029
摘要:本文利用NCEP/NCAR逐日再分析资料,研究了中国大范围持续性极端低温事件(EPECE)的平流层前兆信号及其对对流层环流异常的影响。结果表明,在52个EPECE中,有17个EPECE具有一类共同的平流层前兆信号:(1)在EPECE发生前10天左右,在巴伦支海一带的100 hPa位势高度场呈现较强的正距平特征;(2)随后,该正异常环流逐渐向东移动,并在EPECE发生的前5天左右,使贝加尔湖西北侧200 hPa纬向风显著减弱。本文据此提出了关于EPECE平流层前兆信号的两个判别条件,并以这两个判别条件对1949~2009年冬季(11月至次年3月)所有逐日低频场进行了后查和检验。结果表明,上述两个判别条件对EPECE的发生日期具有一定的预测能力。本文进一步通过位涡(PV)反演探讨了平流层前期异常环流如何影响EPECE发生的物理机制。分析表明,巴伦支海附近的平流层中低层PV异常有利于该地区对流层中上层正高度异常的维持,后者强度的25%来源于前者的作用。
关键词平流层信号     大范围持续性极端低温事件     位涡反演    
A Specific Stratospheric Precursory Signal for the Extensive and Persistent Extreme Cold Events in China
SHI Ning1, BUEH Cholaw2,3    
1 Collaborative Innovation Center on Forecast and Evaluation of Meteorological Disasters, Key Laboratory of Meteorological Disaster of Ministry of Education, Nanjing University of Information Science and Technology, Nanjing 210044;
2 International Center for Climate and Environment Sciences, Institute of Atmospheric Physics, Chinese Academy of Sciences, Beijing 100029;
3 State Key Laboratory of Numerical Modeling for Atmospheric Sciences and Geophysical Fluid Dynamics, Institute of Atmospheric Physics, Chinese Academy of Sciences, Beijing 100029
Abstract: On the basis of daily NCEP/NCAR data sets, this study investigates a stratospheric circulation precursor for the extensive and persistent extreme cold events (EPECEs) in China. Results show that 17 out of 52 EPECEs are characterized by a specific stratospheric circulation precursor. The precursor exhibits itself clearly in two aspects. (1) Strong positive geopotential height anomalies occur around the Barents Sea in the lower stratosphere (100 hPa) 10 days prior to the occurrence of the EPECE. (2) In relation to the eastward movement of the positive height anomalies, the zonal wind over the northwest of Lake Baikal is weakened significantly at 200 hPa about 5 days before the occurrence of an EPECE. According to these stratospheric precursor circulation features, two initial assessment conditions are proposed. To test prediction accuracy, the two assessment conditions are applied to the daily low-frequency data of all winter seasons (November to March) from 1949 to 2009. Results show that our two proposed assessment conditions, at least to some extent, can enable the prediction of the occurrence date of an EPECE. Regarding the physical mechanism of the influence of the stratospheric circulation precursor on the EPECE, potential vorticity (PV) inversion analysis shows that the mid and lower stratospheric PV anomalies around the Barents Sea can induce the geopotential height anomalies at the upper and midtroposphere. In particular, 25% of the real amplitude of the latter originates from the influence of former.
Key words: Stratospheric precursory signal     Extensive and persistent extreme cold events     PV inversion    
1 引言

长期以来,人们认为平流层大气是受对流层影响的被动系统。但自20世纪末期以来,大气科学界逐渐认识到平流层大气和对流层大气是相互耦合在一起。许多观测研究、诊断分析和数值模拟均表明平流层的环流异常也会影响到对流层环流异常(Baldwin and Dunkerton,1999Colucci,2010Hinssen et al.,2010Kodera et al.,1990)。由于平流层异常持续时间较长,其影响不仅具有天气意义,而且有气候意义。因此,充分考虑平流层环流异常的影响应当可以提高地面预测水平(Baldwin et al.,2003Christiansen,2005)。目前,模式的发展也存在着这样一种趋势和共识:能够较好地模拟平流层环流的数值模式能更好地模拟对流层天气系统,如阻塞高压的发展和冷空气的爆发。由此可见,充分利用平流层信号是提高天气气候预报和预测的一个有效途径。通常而言,平流层影响对流层环流异常的动力机制主要是(1)行星波的异常传播、(2)位涡(PV)的影响(Ambaum and Hoskins,2002Hartley et al.,1998Hinssen et al.,2010)以及(3)Rossby波的折射(Nishii and Nakamura,2005)。关于第一个机制,可参考有关综述文献(陈文和魏科,2009胡永云,2006)。

2008年1月(简记为“0801”事件)我国南方地区出现了50年一遇的大范围冰冻雨雪天气(陶诗言和卫捷,2008)。一些研究发现,在此次大范围持续性低温事件(简记为EPECE)前期,平流层极涡出现了明显的加强和偏移(刘毅等,2008王东海等,2008向纯怡等,2009)。但如何提取和利用这些平流层信号以改进我们对EPECE的预测仍是一个难题。因为仅依据个例研究(如“0801”事件),难以提取普遍意义的平流层信号,这需对更多的样本进行深入分析。依据Peng and Bueh(2011)的界定标准,我国自1951年至2009年共发生过52次EPECE。这52个事件均满足极端低温(至少低于10个百分位)、范围大(至少占我国面积的20%)、持续时间长(至少8天)这三个条件。本研究将基于这52个EPECE,重点研究以下两个方面:(1)如何提取最具普遍性的前期平流层信号;(2)讨论平流层信号与对流层异常环流间的动力学链接特征。

2 资料和分析方法

本文采用美国国家环境预报中心/美国国家大气研究中心(NCEP/NCAR)的逐日再分析资料(Kalnay et al.,1996),气象要素场为位势高度(h,后文简称为高度)、水平风速(uv)、气温(T)以及对流层顶高度场。时间段为1948年至2009年的每年冬季(11月至次年3月),这可以覆盖Peng and Bueh(2011)所挑出的52次EPECE发生时 间。除对流层顶高度场外,其他要素场的高度层次均为从1000 hPa至10 hPa(共17层)。为突出环流的缓变演变特征,所有变量场资料均进行了低通滤波,分离出8天以上的低频环流。选用的滤波器为递归六极正切Butterworth滤波器(Kaylor,1977)。图 1中给出了62年冬季平均对流层顶高度随纬度变化特征。可见,对流层顶高度随纬度的增加而降低。在20°N地区,对流层顶位于100 hPa,而在 40°N以北地区则位于300 hPa附近。据此,本研究主要在200 hPa至10 hPa共7层上寻找中高纬度平流层信号。

图 1 1948年至2009年共61年冬季(11月至次年3月)气候平均的对流层顶高度Fig.1 61-year climatological winter-mean(November to March)tropopause level from 1948 to 2009

本文的逐日低频异常场为全场减去其对应日期的气候平均场。为消除年代际变化的影响,这里的逐日气候平均场取为其对应年份的前后5年共11年滑动平均后再31天滑动平均。对于前5年(1948年至1952年)和后5年(2005年至2009年)而言,其逐日气候平均场分别取为前11年(1948年至1958年)和后11年(1999年至2009年)的平均再31天滑动平均。这里取31天滑动平均的目的是在保留月际变率的同时,尽可能消除或减弱逐日气候场中较大的日际变率。在求取距平后,环流异常的逐日变化即可更好地被刻画。

通常而言,“信号”是相对于“噪音”而言。本文将利用局地环流1个标准差(1σ)作标准来区分环流“信号”和环流“噪音”。具体地,我们将强度大于1σ(小于―1σ)的正(负)异常环流称为“正信号”(“负信号”),其余则为“噪音”。

准地转位涡(后文简称位涡,PV)扰动通常定义为

其中, 为准地转扰动流函数;线性算子为静力稳定度,它仅是气压p的函数,这里取为40°N以北的气候平均值。α 为比容。为简洁起见,文中所提到的区域平均是指以格点面积为权重的区域平均值。PV一个重要的性质是其可反演性(Charney and Stern,1962Hoskins et al.,1985)。具体而言,在给定q’分布以及边界条件下,通过反算Lg可得到Ψ'的分布,即。本文采用分段位涡反演方法Davis(1992)分析平流层信号与对流层环流异常之间的联系。上、下边界条件取为Neumann条件:

3 平流层信号的提取

本文以每个EPECE的开始日期为参考日期,记为第0天。开始日期前(后)n天,记为第―n(n)天。需指出的是,若直接对52个EPECE进行合成分析,易导致正、负信号相消或被削弱。为解决该问题,本文首先定义了平流层信号一致性指数CISS(it),其中,i为北半球在平流层中(200 hPa至10 hPa)任意一空间格点的位置,t为时间。CISS的数值为某气象要素场在第t天第i格点上出现局地“信号”的EPECE个数,最后再将CISS转换成百分率(除以52,再乘以100)。CISS极大值区即是52个事件中较多的事件在该时间、该空间位置上出现的一致“信号”。按此方法,可分别得到正、负信号的时空分布特征。

经统计,正值信号的CISS值大值区(大于30%)最早出现在h场上第―10天左右100 hPa巴伦支海附近(图 2)。值得注意的是,在更早或其他层次上也能出现CISS大于30%的区域,但由于其空间范围较小或在随后的演变过程中不连续,这里暂不作进一步讨论。根据图 2,利用标准化的位势高度场()定义了区域信号指数HI,即在(70°~80°N,0°~60°E)范围内(简记为A区)的区域平均值。经统计,52个EPECE中在第―10天100 hPa上HI大于1的事件共有17个(表 1),这对应着图 2中局地30%的CISS指数。这里HI的阈值仍取为1σ,以保证其强度达到本文所提到的“信号”的标准。图 3a给出了这17个事件的HI的高度—时间剖面图。可以看出,信号(大于1σ)最早出现在第―13天100 hPa上,并在第―10天左右达到极值(1.4σ)。此后,HI有所减弱。至第―5天,HI在各层上的值已小于1,“信号”消失。因此,虽然大于1σ的信号能在1000~20 hPa各层出现,但最早和最强的信号均出现在100 hPa上。根据这一前期平流层信号,本文提出EPECE开始日期的判别条件一(简称条件一):100 hPa上HI>1σ,且至少能持续5天,则这5天后的第8天左右可能会发生EPECE。

图 2 第―10天100 hPa CISS指数(%),即局地高度异常大于+ 的事件个数百分率Fig.2 Distribution of CISS(consistency index of stratospheric signal)defined at 100 hPa at day ―10,which represent the percentage of events with the local height anomalies exceed +

图 4给出了17个EPECE合成的100 hPa上及其随时间的演变特征。由图 4可见,巴伦支海附近的正高度异常随时间缓慢向东移动。这也对应着HI在第―10天后的局地减弱(图 3)。至第0天时,其中心已位于(70°N,85°E)附近(图 4c)。值得注意的是,这与“0801”事件中斯堪的纳维亚型环流极罕见的向东移动过程类似(Bueh et al.,2011b)。需注意的是,在该异常环流东移过程中,巴尔喀什湖地区出现了显著负值,它逐渐向东伸展至贝加尔湖及以东地区。与的演变相对应,贝加尔湖西北侧地区的标准化的纬向风场()出现了强烈的信号(图 3b)。这里采用与定义100 hPa上HI类似的方法,把在区域(55°~65°N,80°~100°E)(简记为B区)内的平均值定义为UI。从合成结果来看(图 3b),显著UI负值首先出现在第―18天100 hPa上,只是其强度较小(―0.4σ)。随着巴伦支海附近的正高度异常缓慢向东移动,UI逐渐增大。第―5天,UI<―1σ几乎出现在850~100 hPa各等压面上,其中心位于对流层中上层及平流层底层。这种形势一直维持至第0天左右。因此,就平流层信号而言,已由第―10天左右100 hPa上的HI信号转变为第―5天平流层底层的UI信号。据此提出判别条件二(简称条件二):若第―5天至第0天中至少有一天出现200 hPa上UI<―0.5σ,则进一步确认由条件一预测的事件可能是EPECE。实际上,条件二也间接表征了平流层上巴伦支海地区正高度异常的东移特征。比较条件一和条件二后可以看出,无论从持续时间(1天)还是从强度(0.5σ)上来看,条件二的要求均低于条件一。这样做的目的是突出在更早时间上用条件一作出的预测结果,而条件二对应临近EPECE的发生,它更多地起着定性修正条件一预测结果的作用。

图 3 (a)17个EPECE合成的(70°~80°N,0°~60°E)区域平均信号指数HI的时间—高度剖面,阴影为通过99%显著性检验的区域;(b)同(a),但为标准化纬向风异常场的(55°~65°N,80°~100°E)区域平均值。(a)和(b)中粗线分别为+1 、―,代表强信号Fig.3 (a)Time-pressure cross section of the composited HI averaged over(70°-80°N,0°-60°E)of 17 EPECE events. Shading indicates the area with 99% confidence level;(b)same as(a),but for the normalized zonal wind anomalies averaged over(55°-65°N,80°-100°E). Thick lines in(a) and (b)represent positive and negative 1s respectively,both of which indicate the strong signals

图 4 17个EPECE合成的100 hPa标准化高度异常场:(a)第―10天;(b)第―5天;(c)第0天。等值线间隔为0.2 ;实(虚)线代表正(负)异常,零线已略去;阴影为通过99%显著性检验的区域Fig. 4 Evolution of the composited normalized height anomalies at 100 hPa for 17 EPECE events.(a)day ―10;(b)day ―5;(c)day 0. Contour interval is 0.2s. Solid(dashed)contours indicate the positive(negative)anomalies,zero lines are omitted. Shading indicates the area with 99% confidence level
4 平流层信号的检验

为检验我们的判别方法(条件一及条件二)在预测EPECE上的准确性,首先将条件一应用于1948年至2009年冬季(11月至次年3月)的逐日低频场上,对所有可能的EPECE的开始日期做出预测,并允许前后3天的误差。需指出的是,考虑到100 hPa上出现HI>1σ的连续天数可能远多于5天,我们仅在第一次连续5天HI>1σ时预测其后8天可能出现一次事件,而不是每5天就预测一次事件。最终预测出97个事件(这里不称其为EPECE,仅称其为事件)。表 1中17个EPECE有14个事件(“*”标注)被挑出,这其中包含着“0801”事件。剩余3个事件(表 1中斜粗体)则由于前期HI>1σ的持续时间较长,预测的EPECE开始日期要早于实际日期6天以上,因而没被记入。因此,如不考虑除这17次EPECE外的其他35次EPECE(Peng and Bueh,2011),仅依靠条件一做出的预测准确率大约为14/97×100%≈14.4%。

表 1 满足第―10100 hPa上HI>1σ的17EPECE的开始日期和用判别条件一预测的开始日期以及两者之差 Table 1 Start dates of every 17 EPECE events with HI>1σ at 100 hPa at day―10,start dates of events predicted by the judging condition one and the differences between the real dates and the predicted ones

利用条件二对97个事件进一步分类后发现,除上述14个EPECE满足条件二外,还有50个事件满足条件二,这里称其为潜在的EPECE(简记为P_EPECE)。不满足条件二的剩余事件称为一般事件(简记为CE),共有33个。因此,若加入条件二,可将预测准确率提高至14/(97-33)×100%≈21.9%。此外需要注意的是,14个EPECE加上50个P_EPECE共64个事件,几乎是CE个数(33个)的两倍。这表明用条件一挑出的所有事件中有2/3事件会在5天后左右出现纬向风场在B区显著减弱的现象。

图 5给出了这14个EPECE、50个P_EPECE以及33个CE在第0候(即连续5天HI>1σ后的第6至第10天平均)的异常环流合成图。可以看出,EPECE与P_EPECE无论在300 hPa上高度异常场(图 5a,5b)还是在地表气温异常(Ts图 5d、e)的分布上均十分类似。正、负高度异常分别位于新地岛附近和咸海经贝加尔湖至我国东北北侧,即呈现EPECE中典型的大型斜脊/斜槽式的环流异常分布(Bueh et al.,2011a)。显著Ts负异常则在里海经贝加尔湖至我国东北北侧出现,并有向我国东部南伸的趋势,其中心位于贝加尔湖附近。在下一候(图略),随着环流系统的东移南下,显著的Ts负异常在两类事件中均向南伸至我国东南地区。而在CE事件中,第0候异常环流型(图 5c)几乎与前两者反相(图 5a,5b)。在巴伦支海的下游乌拉尔山以东地区出现了显著负异常(图 5c),该负异常的稳定维持可能和前期在上游巴伦支海及斯堪的纳维亚地区维持的正异常(图略)向下游频散波能量有关,呈现出驻波的特征,是一种较为典型的斯堪的纳维亚环流型(Bueh and Nakamura,2007)。Ts负异常主要位于新地岛附近(图 5f)。以上结果表明,条件二有助于定性修正条件一的预测结果。

图 5 合成的14个EPECE(a,d)、50个P_EPECE(b,e)以及33个CE(c,f)在第0候(以第0天为中心日期的前后2天共5天平均)的异常环流场。左列为300 hPa高度异常场,单位为m;右列为地表气温异常场,单位为°C。等值线间隔在(a)中为50 m,(b)、(c)中为30 m,(d)中为3 °C,(e)、(f)为2 °C;实(虚)线代表正(负)异常,零线已略去;阴影为通过99%显著性检验的区域Fig. 5 Composited anomalies at pentad 0(5-day mean centered at day 0)for(a) and (d)14 EPECE events,(b) and (e)50 P_EPECE events and (c) and (f)33 CE events. Left column is height anomalies at 300 hPa,units are m; right one is the surface air temperature anomalies,units are °C. Contour interval is 50 m in(a),30 m in(b) and (c),3 °C in(d) and 2 °C in(e) and (f). Solid(dashed)lines represent the positive(negative)anomalies,zero lines are omitted. Shading indicates the area with 99% confidence level

总之,(1)条件二的加入有助于我们进一步判断由条件一做出的预测是否定性成立,进而提高 预测准确率;(2)虽然由两个条件预测出的事件个数多于实际EPECE个数,但这些事件在开始发生时仍表现出与实际EPECE十分类似的异常环流特征。 5 可能的影响机制

图 6给出了17个EPECE的(60°N~90°N,30°W~120°E)区域平均的合成PV扰动场。可以看出,PV异常首先出现在30~50 hPa的中平流层(第―20天左右,强度约为―1×10−5 s−1),随后逐渐向下伸展并增强。第―15天左右,对流层中上层开始出现较强的PV负异常。第―10天,出现了两个PV异常的极小值区,分别位于对流层400 hPa左右和平流层50~70 hPa左右。从PV异常的水平分布来看,第―10天,50 hPa上主要存在两个异号的PV异常(图 7a),其中心分别位于东半球的巴伦支海和西半球的维多利亚岛附近。由于西半球的平流层PV异常对东半球对流层高度异常的影响较小(图略),后文将重点讨论东半球平流层PV异常的影响。与平流层情形不同的是(图 7a),在400 hPa(图 7b)上,PV异常的分布形势较复杂,出现了多个大范围的正、负PV异常中心。经分析,尽管这些中高纬度的正、负PV异常各自能在平流层产生引起较强的高度异常,但它们的影响存在相互抵消作用。因此,后文将分析对流层(45°N~90°N,60°W~120°E)范围内所有PV异常对平流层高度异常的整体影响。若对整个北半球的对流层PV异常反演,后文结果不会被定性改变。需指出的是,位于巴伦支海及其下游巴尔喀什湖北侧的负、正PV异常(图 7b)基本对应着EPECE典型的大型斜脊、斜槽的异常环流形势。Bueh et al.(2011b)指出,斜脊与斜槽之间存在着强烈的波能量向下游频散 的特征。换而言之,斜脊的出现有利于斜槽的出 现。据此,本节将以第―10天(平流层信号最强)为例,重点分析平流层PV异常如何影响巴伦支海附近的环流异常。

图 6图 3(a),但为合成的PV异常场,区域平均的范围为(60°~90°N,30°W~120°E),单位为10−5 s−1。浅阴影和深阴影分别为小于―1×10−5 s−1和―2×10−5 s−1的区域Fig. 6 Same as Fig. 3(a),but for composited PV anomalies averaged over(60°-90°N,30°W-120°E). Units: 10−5 s−1. Light and heavy shading marks the regions less than―1×10−5 s−1 and ―2×10−5 s−1,respectively

图 7图 6,但为(a)50 hPa和(b)400 hPa上的PV异常。深、浅阴影分别为小于―3×10−5 s−1、大于3×10−5 s−1的区域;等值线间隔在(a)和(b)中分别为1.5×10−5 s−1和3×10−5 s−1;零线均已略去Fig. 7 Same as Fig. 6,but for composited PV anomalies at(a)50 hPa and (b)400 hPa. Heavy and light shading marks the regions less than ―3×10−5 s−1 and larger than 3×10−5 s−1. Contour intervals in(a) and (b)are 1.5×10−5 s−1 and 3×10−5 s−1,respectively. Zero lines are omitted

图 8给出了17个EPECE第―10天300 hPa上的高度异常场(图 8a)及通过反演PV异常场得到的高度异常场(图 8b、c)。可以看出,通过反演整个PV异常(对流层和平流层)后,巴伦支海附近的高度正异常(图 8b)的强度与原高度异常基本一致(图 8a),中心强度均约为+280 m。这为后文进一步利用分部PV反演分析提供了基础。实际上,若为分析局地影响而将公式(1)中的σ 取为巴伦支海附近[如(50°N~90°N,40°W~130°E)]区域平均的气候平均值,本文结果不会被定性的改变。若仅考虑平流层(200 hPa至10 hPa)巴伦支海附近(60°N~90°N,30°W~120°E)的PV异常,反演得到的300 hPa上的高度异常的中心位置与实际位置基本吻合,其强度达+80 m左右,约占整个高度异常(图 8a)强度的1/4~1/3。500 hPa的情况与此类似(图略)。实际上,通过反演平流层各层PV异常后发现,这种影响主要来自于该地区70 hPa及其以下层次上的平流层PV异常(图略),这与前人认为的平流层中低层的环流异常更易引起对流层环流异常的结论相一致(Colucci,2010Hinssen et al.,2010)。

图 8 第―10天300 hPa上的高度异常场:(a)17个EPECE的合成;(b)整个PV异常的反演场;(c)仅反演平流层PV异常;(d)同(a),但为100 hPa上的高度异常场;(e)仅反演对流层PV异常得到的100 hPa高度异常。单位为m;(a)、(b)和(d)中的等值线间隔为40 m,(c)和(e)中为20 m,零线均已略去Fig. 8 Height anomalies at 300 hPa for 17 EPECE events at day ―10.(a)composite results for 17 EPECE events,(b)height anomalies by inversing the whole PV anomalies,(c)height anomalies by inversing the only stratospheric PV anomalies.(d)Same as(a),but at 100 hPa.(e)Height anomalies by inversing the only tropospheric PV anomalies. Units are m. Contour intervals are 40 m in(a),(b),and (d),20 m in(c) and (e). Zero lines are omitted

若仅反演对流层(45°N~90°N,60°W~120°E)范围内所有PV异常,可发现其在平流层中的影响随高度减弱。这里仅讨论100 hPa上的反演高度异常场(图 8e)。与图 8d比较,反演出的100 hPa正高度异常中心的位置处于实际中心位置的偏西北,且其中心强度较弱(+20 m左右),约占实际异常中心强度的1/16。这表明来自对流层PV异常的整体局地影响较弱。实际上,该情况在更早时刻(比如第―15天)也是如此。

综上所述,前期巴伦支海附近的平流层PV异常可影响到该地区的对流层高度异常,进而与EPECE的发生建立联系。

6 结论和讨论

本文利用NCEP/NCAR再分析资料,对52次我国大范围持续性极端低温事件的前期平流层信号进行了研究。结果表明,有17个事件在其开始日期前10天左右巴伦支海附近出现正高度异常信号(>1σ)。该正高度异常信号随后缓慢向东移动,并在开始日期前5天左右引起贝加尔湖西北侧地区纬向风场的异常。至于平流层高度场信号如何联系到10天后EPECE的发生,本文通过PV反演分析发现,巴伦支海地区的平流层PV异常可解释对流层中高度异常强度的约1/4。这同时也表明,平流层前兆信号并不是导致EPECE事件是否发生的唯一或决定因子。

根据上述环流特征,本文提出了EPECE的两个平流层前兆信号的判别条件。条件一为巴伦支海附近出现正信号(>1σ),若其能够持续至少5天,则认为在该第5天后的第8天可能会发生EPECE。条件二,由条件一预期的发生日期前5天内至少有一天在200 hPa上贝加尔湖西北侧地区出现纬向风减弱(<―0.5 σ)。通过对历史资料进行分析后,除实际发生的14个EPECE外,还有50个事件满足这两个条件。进一步对这50个事件进行合成分析发现,对流层上层300 hPa的高度异常呈大型斜脊/斜槽式分布,我国北方地区出现了地面气温负异常。该结果与实际EPECE的合成结果十分类似。因此,至少可以说,本文提出的判别条件一基本可以提前一周判断是否会发生一次低温事件。而条件二可以进一步提高其预测的准确率。但必须指出,预测出的低温事件能否持续或是否极端(即满足Peng and Bueh(2011)提出的EPECE标准),还需加强对对流层信号和其他平流层信号的监测。

从个例研究来看(Bueh et al.,2011b),对流层中的低频Rossby波的传播及高频瞬变涡动的反馈强迫作用对EPECE的发生和维持起着十分重要的作用。此外,Zeng(1983)从理论上指出,扰动的形态与急流位置之间的配置关系能影响扰动的维持和发展。符仙月和布和朝鲁(2013)也发现,降水偏多的和偏少的EPECE在斜脊/斜槽空间尺度、高空的急流异常、低纬度环流系统的异常方面上均表现出较大的差异。因此,若要进一步排除50个P_EPECE,还需进一步分析这50个事件在对流层动力学特征以及对流层环流异常上与14个EPECE的差异。这是我们下一步研究的重点。

此外,依据Peng and Bueh(2011)提出的5类EPECE空间分布型,前期出现共同平流层信号的17个事件中有9个事件属于全国类,4个东部类,3个西北—江南类,1个东北—华北类。可见,低温事件空间分布情况复杂,同一平流层预兆信号并不唯一对应某一空间型事件。这可能与不同类型事件中(1)巴伦支海附近的平流层环流异常在下传过程中出现差异进而引起对流层中高纬度环流异常差异、(2)对流层中低纬度环流异常差异以及(3)辐射异常差异等因素有关。若能进一步细化区分不同类型低温事件前兆信号之间的差异,必将有利于提高我国对EPECE的预报能力。

为检验本文基于NCEP/NCAR再分析资料所提出的平流层信号判别方法的适用性,本文分析比较了欧洲中期天气预报中心(ECMWF)的逐日再分析资料ERA-40(Uppala et al.,2005)。具体做法是将本文第3节提出的判别方法再应用于ERA-40,重点比较两种再分析资料所得结果的异同。由于ERA-40再分析资料仅覆盖1957年9月至2002年8月时段,分析比较仅在该时段中进行。结果表明,在这42年冬季中,ERA-40资料可挑出46个P_EPECE,而NCEP/NCAR资料可挑出47个(表 2)。若允许两种资料所预测的P_EPECE开始日期之间相差3天、2天和1天以内,则分别有44个、42个和39个事件可被两种再分析资料共同预测得到。这表明本文的判别方法应用于两种再分析资料后所得分析结果基本一致,具有较好的适用性。该结果的出现应当与本文在提取平流层信号时并不使用距平场本身,而是用标准化的距平场有关。由于后者考虑进局地的气候变率(标准差),是一种相对量,因而在讨论环流异常时可以较好地减小不同资料中的系统性误差或差异。

表 2 两种再分析资料得到的P_EPECE开始日期及其两者之间的偏差日数。“—”表示该再分析资料没有分析得到P_EPECE或偏差日缺省 Table 2 Start dates of P_EPECE events based on the two reanalysis data sets and the corresponding different days between them. “—” denotes that the particular event cannot be predicted by the reanalysis data set or that the different days is default

此外需要特别指出的是,52个EPECE中还存在一些事件在前期平流层(10~30 hPa)的中低纬度出现位势高度场的负信号(>―1σ)。有研究表明,对于整个纬向平均的全球大气而言,存在着异常环流从平流层低纬向高纬并向下传播的特征,当抵达极区对流层顶附近后,异常环流会在对流层中从高纬向低纬传播,整个周期较长,约为110天(Cai and Ren,20062007Ren and Cai,20062007陈文和魏科,2009)。实际上,这种起始于热带地区的环流异常的经向传播特征存在于“0801”事件中向纯怡等(2009),它加强了2个月后偏向东亚地区的平流层极涡。但对于其他EPECE,其极涡的变化是否能够联系到2个月前的热带地区环流异常?对该问题的深入研究将为我们在更长时间上预测EPECE的发生提供可能。

参考文献
[1] Ambaum M H P, Hoskins B J. 2002. The NAO troposphere-stratosphere connection[J]. J. Climate, 15 (14): 1969-1978.
[2] Baldwin M P, Dunkerton T J. 1999. Propagation of the Arctic Oscillation from the stratosphere to the troposphere[J]. J. Geophys. Res., 104 (D24): 30937-30946.
[3] Baldwin M P, Stephenson D B, Thompson D W J, et al. 2003. Stratospheric memory and skill of extended-range weather forecasts[J]. Science, 301 (5633): 636-640.
[4] Bueh C, Nakamura H. 2007. Scandinavian pattern and its climatic impact[J]. Quart. J. Roy. Meteor. Soc., 133 (629): 2117-2131.
[5] Bueh C, Fu X Y, Xie Z W. 2011a. Large-scale circulation features typical of wintertime extensive and persistent low temperature events in China[J]. Atmos. Oceanic Sci. Lett., 4 (4): 235-241.
[6] Bueh C, Shi N, Xie Z W. 2011b. Large-scale circulation anomalies associated with persistent low temperature over southern China in January 2008[J]. Atmos. Sci. Lett., 12 (3): 273-280.
[7] Cai M, Ren R C. 2006. 40-70 day meridional propagation of global circulation anomalies[J]. Geophys. Res. Lett., 33 (6): L06818, doi:10. 1029/2005GL025024.
[8] Cai M, Ren R C. 2007. Meridional and downward propagation of atmospheric circulation anomalies. Part I: Northern Hemisphere cold season variability[J]. J. Atmos. Sci., 64 (6): 1880-1901.
[9] Charney J G, Stern M E. 1962. On the stability of internal baroclinic jets in a rotating atmosphere[J]. J. Atmos. Sci., 19 (2): 159-172.
[10] 陈文, 魏科. 2009. 大气准定常行星波异常传播及其在平流层影响东亚冬季气候中的作用[J]. 地球科学进展, 24 (3): 272-285. Chen Wen, Wei Ke. 2009. Anomalous propagation of the quasi-stationary planetary waves in the atmosphere and its roles in the impact of the stratosphere on the East Asian winter climate[J]. Advances in Earth Science (in Chinese), 24 (3): 272-285.
[11] Christiansen B. 2005. Downward propagation and statistical forecast of the near-surface weather[J]. J. Geophys. Res., 110 (D14), doi:10.1029/2004JD005431.
[12] Colucci S J. 2010. Stratospheric influences on tropospheric weather systems[J]. J. Atmos. Sci., 67 (2): 324-344.
[13] Davis C A. 1992. Piecewise potential vorticity inversion[J]. J. Atmos. Sci., 49 (16): 1397-1411.
[14] 符仙月, 布和朝鲁. 2013. 中国大范围持续性低温事件与中国南方降水异常[J]. 大气科学, 37 (6): 1247-1260. Fu Xianyue, Bueh Cholaw. 2013. Wintertime extensive and persistent low-temperature events of China and anomalous precipitation over southern China[J]. Chinese J. Atmos. Sci. (in Chinese), 37 (6): 1247-1260.
[15] Hartley D E, Villarin J T, Black R X, et al. 1998. A new perspective on the dynamical link between the stratosphere and troposphere[J]. Nature, 391 (6666): 471-473.
[16] Hinssen Y, van Delden A, Opsteegh T, et al. 2010. Stratospheric impact on tropospheric winds deduced from potential vorticity inversion in relation to the Arctic Oscillation[J]. Quart. J. Roy. Meteor. Soc., 136 (646): 20- 29.
[17] Hoskins B J, McIntyre M E, Robertson A W. 1985. On the use and significance of isentropic potential vorticity maps[J]. Quart. J. Roy. Meteor. Soc., 111 (470): 877-946.
[18] 胡永云. 2006. 关于平流层异常影响对流层天气系统的研究进展[J]. 地球科学进展, 21 (7): 713-720. Hu Yongyun. 2006. On the influence of stratospheric anomalies on tropospheric weather systems[J]. Advances in Earth Science (in Chinese), 21 (7): 713-720.
[19] Kalnay E, Kanamitsu M, Kistler R, et al. 1996. The NCEP/NCAR 40-year reanalysis project[J]. Bull. Amer. Meteor. Soc., 77 (3): 437-470.
[20] Kaylor R E. 1977. Filtering and decimation of digital time series[R]. Tech. Rep. Note BN 850, 14 pp. (Available from Engineering and Physical Science Library, University of Maryland at College Park, College Park, MD 20740.)
[21] Kodera K, Yamazaki K, Chiba M, et al. 1990. Downward propagation of upper stratospheric mean zonal wind perturbation to the troposphere[J]. Geophys. Res. Lett., 17 (9): 1263-1266.
[22] 刘毅, 赵燕华, 管兆勇. 2008. 平流层环流异常对2008年1月雪灾过程的影响[J]. 气候与环境研究, 13 (4): 548-555. Liu Yi, Zhao Yanhua, Guan Zhaoyong. 2008. Influences of stratospheric circulation anomalies on tropospheric weather of the heavy snowfall in January 2008[J]. Climatic and Environmental Research (in Chinese), 13 (4): 548- 555.
[23] Nishii K, Nakamura H. 2005. Upward and downward injection of Rossby wave activity across the tropopause: A new aspect of the troposphere- stratosphere dynamical linkage[J]. Quart. J. Roy. Meteor. Soc., 131 (606): 545-564.
[24] Peng J B, Bueh C. 2011. The definition and classification of extensive and persistent extreme cold events in China[J]. Atmos. Oceanic Sci. Lett., 4 (5): 281-286.
[25] Ren R C, Cai M. 2006. Polar vortex oscillation viewed in an isentropic potential vorticity coordinate[J]. Adv. Atmos. Sci., 23 (6): 884-900.
[26] Ren R C, Cai M. 2007. Meridional and vertical out-of-phase relationships of temperature anomalies associated with the Northern Annular Mode variability[J]. Geophys. Res. Lett., 34 (7): L07704, doi:10.1029/2006GL028729.
[27] 陶诗言, 卫捷. 2008. 2008年1月我国南方严重冰雪灾害过程分析[J]. 气候与环境研究, 13 (4): 337-350. Tao Shiyan, Wei Jie. 2008. Severe snow and freezing-rain in January 2008 in the southern China[J]. Climatic and Environmental Research (in Chinese), 13 (4): 337-350.
[28] Uppala S M, Kållberg P W, Simmons A J, et al. 2005. The ERA-40 re-analysis[J]. Quart. J. Roy. Meteor. Soc., 131 (612): 2961-3012.
[29] 王东海, 柳崇健, 刘英, 等. 2008. 2008年1月中国南方低温雨雪冰冻天气特征及其天气动力学成因的初步分析[J]. 气象学报, 66 (3): 405- 422. Wang Donghai, Liu Chongjian, Liu Ying, et al. 2008. A preliminary analysis of features and causes of the snow storm event over the Southern China in January 2008[J]. Acta Meteorologica Sinica (in Chinese), 66 (3): 405-422.
[30] 向纯怡, 何金海, 任荣彩. 2009. 2007/2008年冬季平流层环流异常及平流层—对流层耦合特征[J]. 地球科学进展, 24 (3): 338-348. Xiang Chunyi, He Jinhai, Ren Rongcai. 2009. Stratospheric oscillation and stratosphere-troposphere coupling in 2007/2008 winter[J]. Advances in Earth Science (in Chinese), 24 (3): 338-348.
[31] Zeng Q C. 1983. The evolution of a Rossby-wave packet in a three- dimensional baroclinic atmosphere[J]. J. Atmos. Sci., 40 (1): 73-84.