1. 引言

北极极涡是北半球重要的环流系统，极涡的活动变化会对全球的大气环流、气候变化产生重要影响。极涡一般是指位于对流层中上层且中心在极区的大尺度气旋性环流系统，也称绕极涡 [1] 。贯穿平流层和对流层，是大规模极冷空气的源地和象征 [2] 。

对于北极极涡的位置和强度等特征量的表述，目前有许多不同的定义方法。Laseur首次定义并计算了描述极涡的面积大小、位置的物理量 [3] ，对于极涡特征量的研究主要基于天气图进行 [4] 。到了20世纪90年代，除Angell以外 [5] ，其他研究者多借助位势高度场对极涡的大小进行定义和相应的计算 [6] [7] 。从20世纪80年代开始，我国极涡与气温长期预报课题协作组开展了对描述极涡特征的物理量及其与我国气温变化关系的研究，对描述极涡的各个物理量用不同的方法进行定义和分析，最终确定了一套能够反映极涡的环流变化特征的参数方案 [8] 。

在过去几十年里，北极地区的地表气温正在以几乎是全球平均速度两倍以上的速度明显变暖，这种现象被称为北极放大(Arctic Amplification, AA) [9] [10] 。北极放大可以引发一系列气候效应，导致北极海冰减少，永久冻土融化等 [11] [12] [13] [14] ，引起了人们极大的关注。在全球不断变暖的大背景下，北极地区的海冰融化加剧、不断减少，北极、中纬度地区环流特征也有许多新的变化 [15] 。对极涡的研究不仅有利于进一步探究极区与全球大气环流的作用机理，而且对我国极端天气气候事件的研究与预报也很重要。

由于500 hPa等压面大体上代表了地球整层大气的平均层，又接近对流层大气的无辐散层，该层上的大气环流对地面天气系统有很好的引导气流作用，因此本文选取500 hPa极涡南界的特征等值线来表征对流层北极极涡边界，进而计算极涡的面积指数、强度指数和综合指数，分析极涡面积和强度的时空变化规律，为深入理解北极极涡的变化规律并提高其预报水平提供了重要基础。

2. 研究区域与观测数据

2.1. 数据来源

本研究使用的大气资料来自美国国家环境预报中心和大气研究中心(National Centers for Environmental Prediction/National Center for Atmospheric Research, NCEP/NCAR)的逐月再分析数据集，分辨率为2. 5˚ × 2. 5˚，包含位势高度、温度等数据。

2.2. 极涡环流指数定义

2.2.1. 极涡范围的确定

对长周期的预报来讲，可以把北半球(大部分地区)的天气气候变化归因于北极极涡和副热带高压两个环流系统的相互作用 [16] 。因此，可以把北半球极涡外围边界定义在极涡和副热带高压的交界区域上，也就是中纬度(30˚N~60˚N之间)西风急流的急流轴上(或者可以看成位势高度场上经向梯度最大的区域)，近似选取最接近位势高度最大梯度的特征等高线作为极涡的南界。从气候态的500 hPa位势高度场出发，计算位势高度场的经向梯度，选取急流轴在各个经度上对应的点，该点确定为极涡在该经度上的最南界，以此类推，求出所有经线上急流轴对应的点处的位势高度均值，即北极极涡的特征等高线，特征等值线包围的区域我们将其定义为极涡。通过对各月、各季节的极涡特征线的分布变化，结合功率谱分析等方法，我们可以得到极涡范围的季节变化等特征。

根据上述定义，计算500 hPa描述北极极涡南界的特征等值线，即5560 gpm等高线(图1)。

2.2.2. 极涡环流指数定义

北极极涡的范围就是500 hPa特征等值线以北所包围的区域，本文将特征等值线5560 gpm以北包围的格点数定义为极涡的范围指数(Polar Vortex extent, PVE)。

从逐日位势高度数据出发，北半球极涡综合指数(Polar Vortex Comprehensive Index, PVC)的计算公式如下：

$PVC=-{\displaystyle \underset{i=1}{\overset{n}{\sum }}\left(hg{t}_i-\overline{hgt}\right)}$，

其中 $hg{t}_i$为某格点500 hPa的位势高度， $\overline{hgt}$为气候态的极涡范围内位势高度均值，n为极涡范围内格点总数，也即PVE指数。该指数综合考虑了极涡范围和强度对极涡变化的贡献，极涡面积越大，极涡区域位势高度越低，则PVC越大，代表极涡越强。

北半球极涡强度指数(Polar Vortex Intensity, PVI)的计算公式如下：

$PVI=PVC/n=PVC/PVE$，

代表极涡区域单位面积(格点)上的低涡强度。从逐月的极涡强度指数出发，我们可以从多个尺度分析北极极涡强度的时间、空间变化规律。

3. 北极极涡的季节变化

极涡综合指数是一个既可以表示极涡强度，又可以反映极涡面积的综合性指数。从逐月的极涡综合指数出发，我们可以得到极涡综合指数的逐月分布(图2)。结合北极海冰范围的季节变化特征和极涡综合指数在不同月份的差异，12月~3月这四个月的极涡综合指数最大，我们将其定义为冬季(DJFM)；6~9月这四个月极涡综合指数最小，我们将其定义为夏季(JJAS)；4~5月、10~11月为两个过渡季节，我们将其分别定义为春季(AM)和秋季(ON)。后文中的季节变化特征均以此作为季节划分的标准。

3.1. 极涡南界和空间形态的季节变化

分析5560 gpm这条特征等高线逐月的分布，我们可以得到北极极涡南界的季节变化(图3)。可以看出，极涡范围随季节变化显著，7月份极涡范围最小，极涡南界最偏北，大约在75˚N的位置；1月份极涡范围达到最大，极涡南界达到最南端的32˚N。从极涡南界在四个季节的分布特征来看(图4)，春、夏、秋、冬四个季节极涡南界的平均位置分别位于50˚N、62.5˚N、47.5˚N、42.5˚N，极涡南界包围的极涡范围也分别从冬、秋、春、夏季依次递减。

从不同季节极涡范围内等值线的形态分布特征可以发现(图4)，夏季(图4(b))极涡范围最小，极涡呈现出绕极型的分布特征，极涡南界特征等高线表现出四波结构，最强的波槽出现在北太平洋白令海区域，可能与6月阿留申上空的极涡南伸有关(图3(f))。冬季(图4(d))极涡南界最偏南，主要呈现出绕极的三波结构，极涡中心则主要表现为双涡的非对称结构，极涡分别偏向北美和东亚一侧，分别对应东亚大槽和北美大槽这两支定常波，其中偏向北美的极涡中心强度更强(图3(a)~(c)，图3(l))。春季(图4(a))和秋季(图4(c))北极极涡主要表现为绕极的非纬向对称结构，极涡中心呈椭圆型分布，长轴分别在北美和欧亚中部地区。

3.2. 极涡范围的季节变化

北极极涡范围指数的季节变化(图5)可以近似地代表极涡面积的季节变化。极涡范围指数的逐月变化(图5(a))表明，从夏到冬的7~12月，极涡面积逐渐增大，极涡向南扩张；其中9月极涡面积增幅最大，扩张速度最快。从冬到夏的1~7月，极涡面积逐渐减小，极涡收缩；其中6~7月极涡收缩幅度最大。极涡面积极小值出现在7月份，极涡面积收缩到一年中最小；极涡面积在2月达到最大，极涡面积扩张到一年中最大。四个季节的极涡范围变化(图5(b))表明，极涡的范围指数与综合指数(图2(b))呈现相同的位相变化，总体表现为冬季极涡面积最大，夏季极涡面积最小，春季和秋季范围居中，秋季的面积略大于春季的。

3.3. 极涡强度的季节变化

极涡强度指数的逐月变化(图6(a))表明，平均每年1~7月，为极涡的减弱阶段，极涡强度的极小值出现在7月份；其中4~7月即春季和初夏极涡强度减弱的幅度最大，与极涡面积收缩幅度最大的时期相对应。7~12月进入极涡的增强阶段，7~10月即夏末和初秋极涡增强幅度较为明显，极涡强度的极大值出现在1月。极涡在冬季强度最强，极涡外扩；夏季其强度最弱，极涡收缩；秋季比春季略强(图6(b))。

4. 北极极涡的年代际减弱及原因分析

4.1. 极涡环流指数的年代际变化特征

1960~2020年北极极涡的环流指数(范围、强度)在不同季节均表现出显著的代际减弱的特征(图略)，即近几十年来北极极涡范围向北收缩，强度减弱，这在极涡综合指数的时间序列图上清楚可见(图7)。极涡综合指数既能反映极涡面积的变化，也能反映极涡强度的变化，从表1可以发现，PVC在四个季节都与PVE和PVI具有显著的正相关关系，都表征了北极极涡范围退缩、强度减弱的一致变化特征，PVC分别与PVE和PVI在夏季和秋季达到最大正相关，相关系数分别为0.805和0.999。极涡综合指数主要反映了极涡强度的变化，其与PVI之间的相关除了夏季的0.768以外，其余季节都超过0.99。而PVC与PVE的相关只有在夏季比较高，春、秋季节二者之间的相关系数接近0.5，在冬季二者之间的相关性较弱，相关系数只有0.25。这些结果表明夏季极涡范围北缩到极区附近62.5˚N (图4(b))，极涡强度与极涡面积变化呈高度正相关，极涡越强，向中纬度扩张得范围越大；而冬季由于极涡南边界太偏南包含了中纬度罗斯贝波活动区域，长波脊的入侵会减弱北半球极涡的低压环流，使得极涡强度与极涡面积之间的正相关关系大大降低。

总体来看，各个季节的极涡均呈现出振荡下降的特征，从20世纪90年代中期之前的强极涡环流指数变为其后的弱环流指数，突变检验结果(图略)表明北极极涡的年代际转折大体发生在1994/1995前后。这种极涡的年代际减弱特征在秋季最强(图7(c))，夏季最弱(图7(b))，与北极放大在秋季最强、夏季不显著的特征相一致，表明北极极涡近二十多年来的年代际减弱可能与全球变暖背景下北极海冰加速融化，从而造成北极放大增暖的正反馈机制有关 [11] 。空间上相比于1995年之前，近二十年来500 hPa位势高度场在各个季节均有显著的异常升高(图8)，表明北极极涡在四个季节均有显著的年代际减弱特征，特别是北半球秋季，北极区域出现最强的位势高度正异常(图8(c))，极涡的减弱信号最显著，这可能由于夏秋季北极海冰范围最小，开阔水域增加，冰雪反照率正反馈使得上层海洋吸收储存的热量增加，在秋季释放加热大气进而导致气柱受热膨胀位势高度异常升高，极涡减弱。

4.2. 北极极涡的年代际减弱与北极放大

前面的结果表明，北极极涡的年代际减弱可能与全球变暖背景下北极放大相关的海–冰–气耦合的正反馈作用有关。为了研究北极极涡变化与全球变暖、北极放大效应的关系，分别计算了全球、北极(66.5˚N以北的区域)的区域平均的表面温度异常，代表全球变暖和北极放大指数(图7中的绿线和橙色线)。可以发现，除夏季外，伴随着北极极涡在20世纪90年代中期以后的年代际转折，北极极区表面气温的

变暖速率是全球的两倍以上(图7)。北极放大指数与极涡综合指数在各个季节均存在显著的负相关(表2)，总体表现为北极地区变暖，极涡强度减弱，极涡面积收缩。秋季北极异常增暖与极涡综合指数的相关性最高，相关系数可达−0.85，且极涡综合指数和北极气温均在1995年前后存在显著的年代际转折，二者的变化趋势呈现负相关，极涡综合指数的变化主要由极涡强度的年代际减弱所主导(表1)，与图8的结果相一致。

4.3. 北极极涡空间形态的年代际变化和季节差异

北极极涡指数在1990s中期发生的年代际减弱不仅主要体现在北极极涡强度的年代际减弱，极涡的形态和强度变化在空间上也同样存在显著的年代际差异和季节差异。从图9中可以发现，1960~1994年期间，四个季节的500 hPa位势高度距平场均为负异常，表明相对于气候平均态而言极涡加强；而1995~2020年期间则正好相反，北半球各季节几乎均为位势高度正异常，北极极涡整体减弱。但极涡的这种年代际减弱在不同季节具有空间差异性。在极涡较强的秋季和冬季，北极极涡主要呈现出偏心结构特征，极涡中心偏离北极点主要集中在北美加拿大群岛附近，秋季由于北极放大的增暖异常最强，北极极涡中心整体呈年代际减弱变化特征(图9(c)，图9(g))；冬季由于海冰融化北极增暖主要在大西洋扇区一侧，因此北极极涡的年代际减弱也主要体现在极涡偏心在北美东部的极涡中心显著地减弱范围变小(图9(d)，图9(h))。夏季位势高度异常的显著变化主要发生在中纬度四个大值区，位于极涡南边界附近(图9(b)，图9(f))，因此夏季极涡主要呈现出绕极的结构，但极涡边界范围在后一年代时段略向北收缩，沿极涡边界依然呈四波结构特征。春季北极极涡形态的年代际变化最显著，绕极非对称的极涡中心分裂成东西半球的欧亚和北美两个中心，使极涡从单涡绕极型转变为双涡偶极型结构，这可能是因为北极放大增暖改变了大气经向位涡梯度分布，进而影响定常波临界波数和罗斯贝波破碎过程 [17] 所导致的。

5. 结论

本文主要通过定义500 hPa描述北极极涡边界的特征等值线(5560 gpm)，计算了对流层描述北极极涡物理特征的范围指数、强度指数和综合指数，分析了1960~2020年在全球变暖背景下对流层北极极涡强度和形态的年代际变化及其季节性特征，初步探究了其与北极放大的关系。

从季节变化特征来看，北极极涡的范围和强度具有一致的季节变化特征，夏季7月份极涡范围最小，强度最弱；冬季1月份极涡范围最大，强度最强。夏季极涡呈现出绕极型的分布特征，极涡南界特征等高线表现出四波结构；冬季，极涡南界特征等高线主要呈现出绕极的三波结构，极涡中心则主要表现为双涡的非对称结构，极涡分别偏向北美和东亚一侧；春季和秋季主要表现为绕极的非纬向对称结构，极涡中心呈椭圆型分布，长轴分别在北美和欧亚中部地区。

从年际年代际变化特征来看，北极极涡的综合指数与强度指数具有高度一致的年际年代际变化特征，PVC与PVI之间的相关除了夏季的0.768以外，其余季节都超过了0.99。而PVC与PVE的相关只有在夏季比较高(0.805)，春、秋季节二者之间的相关系数接近0.5，在冬季二者之间的相关性较弱，相关系数只有0.25。北极极涡的年际年代际变化主要体现为极涡强度的年代际减弱，年代际转折时间发生在1995年前后，其中秋季的年代际减弱最显著，夏季最弱，与北极放大在秋季最强、夏季不显著的特征相一致。

北极极涡的空间形态在1995年前后也发生了显著的年代际变化，且季节性差异明显。秋、冬季节北极极涡的空间形态均表现为极涡中心偏向北美加拿大群岛附近的偏心结构，北极极涡强度年代际减弱后偏心结构的极涡中心范围显著地减小。夏季，绕极型的北极极涡形态年代际变化特征主要表现为极涡边界的向北收缩以及沿极涡边界四波结构长波槽的减弱。春季，绕极非对称的北极极涡空间形态的年代际变化主要表现为极涡中心分裂成东西半球的欧亚和北美两个中心，使极涡从单涡绕极型转变为双涡偶极型结构。

基金项目

国家自然科学基金(42075024)；山东省自然科学基金重大基础研究项目(ZR2019ZD12)。

NOTES

^*通讯作者。

参考文献