ABSTRACT

Based on the ERA-Interim daily reanalysis data (1979-2014) from the European Center for Medium range Weather Forecasting (ECMWF), storm track is defined as monthly mean variance of geopotential height at 300 hPa with harmony filtering from 2 to 7 days. Temporal and spatial variations of the storm track over northern hemisphere are analyzed. Results indicate that the intensity and location of storm track have a remarkable interdecadal shift around 1998 and an opposite behavior for the Pacific storm track (PST) and the Atlantic storm track (AST). After 1998 when global warming hiatus period began, the “midwinter suppression” phenomenon of PST appears frequently and the PST shifts polar ward. Meanwhile, the AST maximum intensity occurs in October-November early than in February-March and the AST moves toward equator. These “seesaw” characteristics of PST and AST show clearly in the leading mode by a Season-reliant Empirical Orthogonal Function (S-EOF) analysis, especially in cold season. The storm track’s interdecadal shift may relate to the variability of mean flow-wave interaction, which is yet to be investigated.

Keywords:Storm Track, Interdecadal Shift, Pacific, Atlantic, S-EOF

全球变暖停滞背景下北半球风暴路径的年代际变化特征

王刚¹，黄菲^1,2*，王宏¹，樊婷婷¹

¹中国海洋大学物理海洋教育部重点实验室和青岛海洋科学与技术国家实验室，山东 青岛

²宁波大学宁波市非线性海洋和大气灾害系统协同创新中心，浙江 宁波

收稿日期：2017年7月6日；录用日期：2017年7月22日；发布日期：2017年7月25日

摘 要

本文基于1979-2014年欧洲中期天气预报中心ERA-Interim的逐日再分析资料，用300 hPa位势高度场天气尺度扰动方差代表风暴路径，分析了北半球风暴路径的时空变化特征。结果表明，北半球风暴路径的强度和位置变化在1998年前后存在显著的年代际差异，且北太平洋风暴路径(PST)和大西洋风暴路径(AST)表现不同。在1998年后全球变暖停滞阶段PST出现“深冬抑制”现象加强的特征，且风暴轴位置向极移动，AST则出现最大强度从晚冬(2~3月)向深秋(10~11月)提前的变化趋势，风暴轴位置向南移动。这些特征均在风暴路径S-EOF的第一主模态时空特征中得以体现，特别是秋冬季节PST和AST“跷跷板”式的反位相年代际变化，可能与波–流相互作用的年代际差异有关，仍需要更深入的物理机制分析。

关键词 :风暴轴，年代际变化，太平洋，大西洋，S-EOF

Copyright © 2017 by authors and Hans Publishers Inc.

This work is licensed under the Creative Commons Attribution International License (CC BY).

http://creativecommons.org/licenses/by/4.0/

1. 引言

风暴路径是斜压性气旋盛行的带状区域，气旋活动高发区和强气旋活动区同风暴路径活动区域相对应 [1] 。风暴路径对大气环流的维持和天气气候变化起着重要的作用，风暴路径强度和位置的变化同中纬度天气和短期气候变率密切相关。所以，风暴路径一直是气候动力学和短期气候预测理论研究的重点之一。北半球天气尺度上瞬变扰动最活跃的区域主要集中在中纬度两个大洋区域，分别为北太平洋风暴路径(PST)和大西洋风暴路径(AST)。

北半球风暴路径活动具有显著的年代际变化，Chang和Fu [2] 通过分析大气再分析数据和探空仪观测数据，指出早在1970s年代初期其就由强度较弱的风暴路径转变为较强的风暴路径。Nakamura等 [3] 发现在1980s晚期和1990s初期，北太平洋冬季风暴路径也存在显著年代际变化。有分析认为风暴路径发生这种年代际变化可能与AO和ENSO等低频变化密切相关 [2] ，也有研究指出可能还与背景环流和SST的低频变化有关 [4] 。但是近三十年来，在全球气候变暖的大背景下，北极海冰迅速融化，中纬度异常天气频繁发生，而风暴路径的强度和位置也发生了变化。南北半球的风暴路径都开始向极地和垂直向上的方向发生移动 [5] 。Zhang等 [6] 认为风暴轴向极的移动可能是由于冰面融化，表面能量通量释放进入大气底层，海冰边缘向极地中心移动，导致气旋生成区向北移动。同时，由于北极海冰融化，北极温室增暖效应被放大 [7] ，这导致极地向赤道的温度梯度进一步减小，从而通过大气环流的作用改变风暴路径的强度和位置变化。

早期的研究受限于资料的时间长度，且以往对风暴路径年代际变化的研究都集中在上个世纪70年代的冬季 [2] [8] ，近三十年里全球不仅经历了快速变暖的过程，同时在1990s年代末发生了全球变暖停滞的现象 [9] ，更是引起了科学家的关注。在这种全球变暖停滞的背景下，北半球两大洋上的风暴路径会发生怎样的变化？其在全球变暖和变暖停滞不同的年代际背景下时空特征如何？这将是本文重点分析的内容。

2. 资料和方法

本文研究采用的大气资料来源于欧洲中期天气预报中心(ECMWF)提供的ERA-Interim全球逐日再分析资料，包含风场、高度场、温度场以及热通量数据等，垂直方向从1000 hPa到10 hPa共17层，水平分辨率为2.5˚ × 2.5˚，资料起止时间为1979年1月1日到2014年12月31日。

考虑到本文研究侧重于风暴路径的统计特征，我们采用天气尺度的瞬变扰动方差来表征风暴路径。对300 hPa位势高度场的逐日资料进行2~7天的谐波滤波，然后用月平均均方差来表征风暴轴活动强度。以滤波方差最大的格点作为风暴轴中心，其中心所在经纬度作为风暴轴的位置。通过36年(1979~2014年)北半球300 hPa位势高度扰动方差的平均分布，将大于等于30 dagpm²区域定为风暴路径主轴区(图1)。

可以看出，北半球主要存在两大风暴路径，分别位于北太平洋和北大西洋上空：北太平洋风暴路径(Pacific Storm track, PST)起始于120˚E跨过太平洋抵达北美西海岸120˚W；北大西洋风暴路径(Atlantic Storm track, AST)则位于90˚W向东跨越整个大西洋到0˚附近。因此，分别取120˚E~120˚W，35˚N~70˚N作为PST的中心活动区域，取90˚W~0˚E，35˚N~70˚N作为AST的中心活动区域(如图1红色扇形区域)。

3. 风暴路径的年代际变化特征

3.1. 两大洋上风暴路径强度的年代际变化

北半球风暴路径强度存在很强的季节变化，冬季最强位置偏南，夏季最弱位置偏北。这种季节变化与大气斜压性和急流的变化基本一致，其强度和位置主要受到太阳辐射的影响。我们分别选取PST和

图1. 北半球风暴路径气候平均(1979~2014)的空间分布，单位：dagpm²

AST在300 hPa位势高度场2~7天滤波方差的中心最大强度表征风暴轴强度，分析其强度随着季节和年际变率的分布情况(图2)。

可以发现，北半球两大洋上的风暴活动均表现为秋冬半年(10月~次年3月)强，春夏半年(4~9月)弱的特征，最小值发生在7月，但秋冬季节的最大值特征两大洋上表现不同。AST在整个秋冬季节保持强盛，在1998年以前的全球变暖阶段，AST最大强度的出现时间从深秋季节(11~12月)逐渐向隆冬季节(1~3月)推迟，随后在全球变暖停滞阶段AST最大强度又逐渐从晚冬(2~3月)提前到深秋时节(10~11月)。AST的强度变化存在明显的年际振荡特征，大约每隔3~5年会有一个较强中心出现，这与胡增臻和黄荣辉 [10] 的结论相一致。

PST在秋冬季节并不像AST一样一直处于较强的状态，而是在11月、12月份强度较大，却在1月、2月强度发生减弱，3月、4月继续增强，呈现出双峰结构的特征。Nakamura [11] 把这种在晚秋和早春活动强盛，而在深冬季节显著减弱的现象称为PST的“深冬抑制”。有研究认为这可能与非绝热加热效应、地形 [12] 以及急流在对流层上层发生的卷吸效应有关 [3] 。值得注意的是，PST的这种“深冬抑制”特征存在显著的年代际变化。特别是1998年之后的全球变暖停滞阶段，PST的“深冬抑制”特征加强，PST强度在晚冬时节显著加强。

图2. 1979-2014年(a)AST和(b)PST风暴路径强度随着季节–年的变化分布图，单位：dagpm²

3.2. 两大洋上风暴轴位置的年代际变化

从以上分析，我们发现，PST和AST的强度变化在二十世纪九十年代后期全球变暖停滞阶段均存在显著的年代际转型特征，这与朱伟军和李莹 [13] 曾指出的北太平洋冬季风暴路径强度增强或减弱的主体一致变化型模态的年代际变化不尽相同，同时他们还指出北太平洋风暴路径年代际变化的第二模态是风暴路径中东部位置南北移动的经向异常型。因此我们就冬季(1~3月)PST和AST中心纬度位置的变化进行分析(图3)。

从图3可以很清楚地看到，PST和AST的中心纬度变化存在显著的年代际变化特征，其转折点同强度类似，都发生在二十世纪九十年代后期的全球变暖停滞开始阶段，但二者呈现出相反的南北移动特征。冬季PST强度由弱变强(图2(b))，位置由南向北移动(图3(b))；AST在强度和南北位置上都与PST呈现相反的年代际变化，即冬季AST强度由强变弱(图2(a))，位置由北向南移动(图3(a))。我们分别将大于70 dagpm²的等值线所在最高和最低纬度值定义为风暴轴的南、北边界，值得注意的是，随着PST强度加强，其南北的边缘范围也有所扩张；而AST除了在2010年中心位置异常偏南导致其北部边缘出现了极度收缩以外，其南北边缘在整个全球变暖停滞期并没有明显的收窄。

风暴路径这种中心位置出现的年代际变化，可能主要由于急流核和风暴轴中心位置的相互配合，在年代际尺度上显示出一致性的南北摆动 [14] [15] 。韩博等 [16] 在研究PST异常变化特征及其与纬向风场的关系时，指出高层纬向风场异常与PST发展存在着相互反馈的作用，是大气上游异常能量向下游地区频散的一种表现。在全球变暖停滞背景下，为什么北半球风暴路径出现PST位置向极加强，而AST位置向极减弱这种“跷跷板”式的反位相变化现象，值得我们进一步深入分析。

图3. 冬季(a)AST(红色实线)和(b)PST(蓝色实线)中心纬度的时间变化序列，灰色阴影区为风暴路径强度的南北边界(大于70 dagpm²的等值线所在最高和最低纬度值) 单位：N˚

3.3. 北半球风暴路径的时空主模态特征

为了能系统地分析风暴路径活动在时间和空间上的主要模态特征，经验正交函数(EOF)方法可以提供很好的分析手段，但它只会提取风暴路径最强发生季节–冬季的时空分布信息。而我们知道风暴路径具有非常强的季节变化特征，为了真实地反映风暴路径的季节演变随年份变化的时空主模态特征，本文首次引入Wang和An [17] 提出的季节经验正交函数(Season-reliant Empirical Orthogonal Function; S-EOF)分析方法对北半球风暴路径进行了时空特征分析(图4)。在S-EOF分解中，分别取7~9月(JAS)代表夏季，10~12月(OND)代表秋季，次年的1~3月(JFM + 1)代表冬季，次年4~6月(AMJ + 1)代表春季，得到从夏季到次年春季随季节演变的北半球风暴路径时空分布的主模态，前8个模态的总方差贡献率为47.75%，其中第一主模态(图4)的方差贡献百分率为10.17%。

可以发现，第一模态主要反映了北半球风暴路径的年代际变化特征，在1998年之前，风暴路径的主成分(PC1)主要以负值为主，而在1998年之后转为以正值为主，表明在1998年之后的全球变暖停滞阶段北半球风暴路径发生了显著的年代际转型。PC1的极大值分别出现在2000、2005、2008、2009、2012和1984年，其中大部分出现在2000年以后；而极小值则分别出现在1982、1992、1995、1997和2004年，其中大部分都出现在上世纪，极大极小值出现的年份也从一个侧面体现了风暴路径第一主模态的这种年代际变化特征。另外，对照ENSO(厄尔尼诺–南方涛动)指数(MEI, http://www.esrl.noaa.gov/psd/enso/mei/#data)，可以发现PC1的极值年份所体现的年际振荡与ENSO有一

图4. 风暴路径SEOF分解的第一模态不同季节的空间分布型(a)~(d)以及对应的标准化时间系数(e)。图(a)~(d)中黑色等值线为多年气候平均的风暴路径主轴区分布，图e中的红色线为时间系数(左轴)，绿色柱为5年滑动平均序列(右轴)

定的相关性，PC1极小值年多为厄尔尼诺年，而极大值年多为拉尼娜年。

风暴路径第一主模态的空间上主要反映出风暴路径冬强夏弱的季节锁相特征，且秋、冬季节PST和AST呈现出一种“跷跷板”式的反位相变化特征，即当PST增强时，AST异常减弱。结合PC1的时间变化可以发现，在全球变暖停滞的秋季，PST活动减弱，AST活动加强。而在冬季，PST异常加强区位于风暴轴心区的东侧和北侧，表明PST向极异常加强；AST则在风暴路径的北侧东侧异常减弱，而在轴心区的东侧南侧异常加强，并向下游的欧洲地区扩展，表明AST有向极异常减弱的特征。这与前述分析中两大洋上冬季风暴路径强度和位置的年代际变化特征相一致，PC1与PST中心位置变化曲线(图3(b))相关系数达到−0.51，超过了99%的显著性水平。这些结果与前人 [18] [19] 分析的全球变暖未来情景下AST和PST一致加强或减弱的结果不一致，可能与未来情景模拟中气候模式不能准确模拟大气内部变率的年代际变化有关。

为了证实SEOF的第一模态主要反映了北半球风暴路径的年代际变化特征，根据PC1在1998年前后出现年代际突变，将1979-2014年分为两个阶段，即P1阶段(1979/80~1997/98)和P2阶段(1998/99~ 2013/2014)，分别对秋、冬季节两个阶段的风暴路径做差值场(图5)。可以看到秋季(OND)时，PST在P2期间比P1显著减弱，负异常从风暴路径入口区向东南一直延伸至北太平洋中部，而AST的西南部减弱，东北部略有抬升并向下游伸展加强。在冬季(JFM)时，PST和AST的反位相变化更为显著。PST在P2阶段显著比P1阶段强，整体有明显向极地方向移动并向下游扩展至北美大陆，南部强度有所减弱。然而，

图5. 风暴路径分别在秋季(OND)和冬季(JFM)P1阶段(1979~1997年)与P2阶段(1998~2013年)的差值场，单位：dagpm²，打点区域为过置信度95%的区域，黑色实线为风暴路径气候平均态(1979~2013年)，绿色等值线为300 hPa纬向风场，等值线间隔为2 m/s

AST在气候平均位置的东北部强度则显著减弱。对比S-EOF的空间分布(图4(b)，图4(c))可以发现P2-P1的年代际差值场与S-EOF空间型分布基本一致，表明第一模态主要反映了北半球风暴路径的年代际变化特征，特别是冬季PST向极加强和AST向极减弱的反位相特征主要出现在全球变暖停滞阶段。

诸多研究表明 [14] [15] [20] ，风暴路径的变化同中纬度高空急流的异常变化具有相当高的相关性，这体现了波–流之间的相互作用结果。从图5中可以发现异常纬向风的中心与风暴路径的正负异常中心位置基本上相配合，表现出全球变暖停滞期PST的向极加强和AST的向极减弱均对应着西风急流的在太平洋上的向极加强和大西洋上的向极减弱。两大洋上这种低频纬向流–高频天气尺度波之间的反位相关系在年代际尺度上为什么会发生反转？特别是在太平洋上秋季风暴活动的年代际减弱主要出现在风暴路径的上游日本海一带和风暴轴的南侧，而冬季风暴路径年代际加强区为上游日本海一带以及风暴轴的北侧太平洋海盆区，这种季节上的年代际差异可能原因是什么？都值得我们关注。限于篇幅原因，有关机制探讨的工作将在后续文章中进行深入研究。

4. 总结和讨论

本文使用1979~2014年欧洲中期天气预报中心ERA-Interim的逐日再分析资料，用300 hPa位势高度场天气尺度扰动方差代表风暴路径，分析了北半球风暴路径的时空变化特征。结果表明，北半球风暴路径的强度和位置变化在1998年前后存在显著的年代际差异，且北太平洋风暴路径(PST)和大西洋风暴路径(AST)表现不同。在1998年后全球变暖停滞阶段PST出现“深冬抑制”现象加强的特征，且风暴轴位置向极移动，AST则出现最大强度从晚冬(2~3月)向深秋(10~11月)提前的变化趋势，风暴轴位置向南移动。

通过引入能够表征季节演变时空特征的S-EOF分析方法，对北半球风暴路径进行S-EOF分析，发现风暴路径的第一主模态时空特征很好地体现了上述不同季节PST和ASP的年代际变化特征，第一模态方差贡献10.47%，时间上主要反映了1998年前后的年代际转变特征，且在年际尺度上与ENSO有一定的关系。空间上反映出PST和AST呈“跷跷板”式的反相关变化，表现为全球变暖停滞阶段的冬季PST向北加强、AST向北减弱的年代际变化特征，特别是PST在秋冬季节在风暴轴的上游和南北两侧存在相反的表现。这些结果显然与前人 [18] [19] 分析的全球变暖未来情景下AST和PST一致加强或减弱的结果非常不一致，其原因可能是未来情景模拟中气候模式不能准确地模拟大气内部变率的年代际变化特征，而观测中全球变暖和变暖停滞可能与北太平洋年代际涛动(PDO)等海洋–大气内部变率有关。

研究结果还表明，两大洋上风暴路径的强度和位置的年代际反位相变化与高空纬向风异常密切相关。风暴路径实际是斜压不稳定波活动的一种表征 [1] ，平流对于斜压增长的贡献，可以调节斜压波动的增长环境，这同风暴路径生成发展机制相一致。Hoskins和Valdes [21] 指出，平均气流中的斜压不稳定和非绝热加热作用对风暴路径的维持起主要作用。在全球变暖停滞背景下，北极海冰的快速消融，热带海温的异常变化，都有可能直接或间接地影响到中纬度的西风急流和风暴活动。因此北半球两大洋上风暴路径的这种反相的年代际变化特征，可能与波–流相互作用的年代际差异有关，但这其中的复杂物理过程，仍需要更深入的机制分析。

基金项目

国家重大科学研究计划项目(2015CB953904)和国家自然科学基金项目(41575067)资助。

文章引用

王 刚,黄 菲,王 宏,樊婷婷. 全球变暖停滞背景下北半球风暴路径的年代际变化特征
Interdecadal Shift of the Northern Hemisphere Storm Track under Global Warming Hiatus[J]. 气候变化研究快报, 2017, 06(03): 186-194. http://dx.doi.org/10.12677/CCRL.2017.63020

参考文献 (References)