由于人类大量使用矿物燃料并过度排放二氧化碳和其他污染物所产生的温室效应使得全球变暖，而在近百年全球气候显著变暖的特征之中，冬季变暖最为明显 [1] - [4] 。针对我国冬季平均气温的主要模态所展开的研究已有许多十分有意义的成果 [5] - [10] 。大多数研究发现我国冬季平均气温主要呈全国一致性变化，即一致偏冷或偏暖。

然而，对于我国冬季气温的研究，大多数是针对我国冬季三个月或五个月平均气温的年际与年代际变化或长期趋势展开分析 [11] [12] [13] 。然而从该角度分析我国冬季气温可能会平滑掉冬季气温在季节内尺度的变化 [14] 。已有一些观测事实和研究指出，我国冬季气温在季节内尺度具有显著变化的特征 [15] [16] 。一些观测事实表明，前冬和后冬频繁出现冷暖交替的现象：如2014/2015年冬季气温呈前冷后暖现象，而2015/2016年冬季上旬和中旬为前暖转为后冷 [17] [18] 等。2020/2021年我国冬季为前冬冷转为后冬暖 [19] ，而之后2021/2022的冬季却为前冬暖转为后冬冷。如果对于气温的研究仅仅考虑传统的冬季季节平均，则会掩盖掉其中的次季节变化信号。为了更准确地认识我国冬季气温的变化特征，对其季节内变化的研究是非常必要的。

前人对于前、后冬气温反转的成因及物理机制也有一些研究且尚未达成共识。韦炜等发现，前冬和后冬温度反向演变受到了大气定常波和瞬变波的影响 [20] ，即当我国处于前冬偏暖、后冬偏冷的反向演变模态时，环流异常及定常波的传播、瞬变波的反馈作用也与该情形相对应。黄嘉佑等指出中国冬季的环流系统及近地面气温年际和年代际变化在前冬和后冬的不同致使前、后冬的气温反向变化 [16] 。有研究认为，中国及欧亚大陆的气温存在显著的季节内转变，并且中国气温的主要年际变化与欧亚大陆更大范围的高纬–中低纬气温的反号变化有关 [14] 。Lü等人研究结果表明，秋季北极海冰在前冬和后冬通过不同机制影响西伯利亚高压使得其强度在冬季季节内发生反向变化，而西伯利亚高压的强度与我国冬季气温呈显著负相关，从而使得我国冬季气温在季节内出现明显反向变化 [21] 。Geng等人北大西洋涛动在强厄尔尼诺年的冬季前期出现明显的正负位相反转，并且激发大气瞬变波异常活动，通过大气遥相关影响我国冬季气温在冬季前期和后期出现冷暖交替 [22] 。由此可见，冬季气温存在显著的前、后冬差异，且造成反转的成因也各不相同。

因此，为了更加深入地了解我国冬季气温的变化特征及机理，本文将重点分析我国冬季气温的前、后冬冷暖交替的变化特征并探究季节内冷暖交替的成因。

文使用了中国气象局国家信息中心提供的 1648站逐日平均温度和降水资料，为保证研究期间气温和降水数据的连续性，剔除有缺测及观测站点迁移的站点，最后得到426个台站。以35˚N为标准区分我国南、北方，得到187个北方站点，239个南方站点。大气环流资料为欧洲中期天气预报中心(European Centre for Medium-Range Weather Forecasts，ECMWF)提供的第五代大气逐日平均再分析资料ERA5，水平分辨率为0.25˚ × 0.25˚。本文所用数据的时间段为1980~2020年冬季。上述变量的异常值均去除其对应时间段的气候平均和长期趋势。本文中的冬季指的是上一年的12月到此年的2月，即1980年到2020年的冬季为1979年12月到2020年2月。

本文主要采用了多变量经验正交分解(MV-EOF)，回归分析，合成分析等分析方法。

由于中纬度天气变化与北极冷空气活动关系密切，而北极冷空气的向南爆发往往与北极极涡的变化相联系，因此本文将从北极极涡变化的角度探究中国冬季气温前后冬冷暖交替的成因。分别对PC2高、低值年的前冬和后冬100 hPa位势高度及其异常场的合成( 图5 )分析，可以发现，对于中国冬季发生冷暖交替季节内振荡的模态，北极极涡均表现为纬圈非对称的“纺锤型”形态特征，极涡长轴分别向北美和欧亚大陆两侧伸展，短轴两侧分别对应北欧和北太平洋区域。在前暖后冷年( 图5(a)~(c) )，无论前冬还是后冬，北极极涡总体上异常加强，后冬比前冬极涡更强( 图5(c) )，极区为位势高度负异常，周围中纬度为位势高度正异常，体现了北极和中纬度之间跷跷板式反位相振荡的类北极涛动(Arctic Oscillation, AO)正位相的特征。极涡边界以4波结构为主，且主要受到北欧高压脊和阿拉斯加高压脊的影响。前冬( 图5(a) )，北欧高压脊向北入侵北极，使得极涡形态相对狭长，而后冬( 图5(b) )北欧高压脊位置偏南，入侵北极程度减少，极涡“纺锤型”的短轴相对于前冬更加向中纬度扩张。从 图5(c) 的100 hPa位势高度差值场也能清楚地看到北欧偏北极处有显著的正位势高度异常，即前冬的北欧高压脊强于后冬，且位置偏北。在前冷后暖年( 图5(d)~(f) )，北极极涡在前、后冬均异常减弱，呈现出AO负位相的异常空间分布特征，且后冬AO负位相异常分布更强，极涡减弱也更强( 图5(f) )。中纬度极涡边界主要呈2波结构，纬圈非对称的极涡“纺锤型”形态结构更加明显，同样是受到北欧高压脊和阿拉斯加高压脊向北入侵的影响，特别是在后冬阿拉斯加高压脊的明显增强且向北入侵使得极涡受到挤压后形态内缩且呈狭长偏心形状。

前文的分析可以发现，中国冬季气温季节内冷暖交替的振荡变化不仅与北极极涡的强度异常(AO的正负位相)有关，还与极涡在前后冬的形态变化有关。为了探究冬季北极极涡的形态变化是否存在前后冬反转的季节内振荡，分别对前暖后冷年和前冷后暖年北方和南方区域平均气温距平序列( 图4 )与100 hPa位势高度异常场进行回归分析。结果表明( 图6 )，北极极涡区域主要出现三个位势高度异常显著的区域，分别位于北欧、白令海–阿拉斯加地区和北美地区，呈“三叶型”异常分布，其中北欧和白令海–阿拉斯加地区的位势高度异常与北美区域符号相反。因此采用白令海–阿拉斯加(130˚E~130˚W, 50˚~90˚N) (区域1)、北欧(0˚~60˚E, 50˚~90˚N) (区域2)和北美(110˚W~50˚W, 40˚~75˚N) (区域3)三个区域的100 hPa位势高度异常定义“三叶型”极涡振荡指数：Index = 0.25 h100_1 + 0.25 h100_2 – 0.5 h100_3，其中h100_1、h100_2和h100_3分别代表区域1、区域2和区域3平均的100 hPa位势高度异常。

分别合成PC2高值年( 图7(a) )和低值年( 图7(b) )的“三叶型”极涡振荡指数，可以发现，无论是PC2的高值年还是低值年份，该极涡指数均存在季节内振荡和在一月中旬前后的反转特征。在PC2的高值年份，白令海–阿拉斯加区域的振荡大约从1月12日开始转折，其后依次为北欧区域(1月15日)和北美区域(1月19日)；在PC2的低值年份，“三叶型”极涡振荡中仍为白令海–阿拉斯加区域最先开始振荡，而北欧区域和北美区域近乎同时发生振荡。这些结果表明，在中国发生前后冬冷暖交替现象时，无论是前冬暖转变为后冬冷，还是前冬冷转变为后冬暖，极涡总是呈现出纬向不对称的“纺锤”状形态。当冬季前暖后冷时，极涡形态由前冬的短轴狭窄型转变为后冬的短轴变宽型( 图5 )，反映了极涡形态从前冬的短轴向北收缩到后冬的短轴向南扩张的前后冬振荡特征，这种极涡形态的前后冬振荡变化主要由“三叶型”极涡振荡指数所表征，且总是在前冬北欧和白令海–阿拉斯加地区先后有中纬度长波脊的向北极入侵，进而挤压北极极涡中心向北美一侧加强，形成两个正异常中心和一个北美负异常中心的“三叶型”极涡异常结构；在后冬该振荡指数为负位相( 图7(a) )，表明北欧和白令海–阿拉斯加脊异常减弱，“纺锤”形极涡短轴同时向南扩张，这样就构成了北极极涡从前冬向极收缩到后冬向南扩张的极涡季节内振荡特征。而前冷后暖的冬季里，“三叶型”极涡振荡指数由前冬整体为负位相转为后冬的正位相( 图7(b) )表征北欧和白令海–阿拉斯加脊异常在后冬强度增强，使得极涡从前冬的短轴向南扩张到后冬的短轴向北收缩。

北极极涡形态在前、后冬的这种收缩–扩张的振荡与中国冬季前后冬冷暖交替的变化特征密切相关。在前暖后冷年，“三叶型”极涡振荡指数与中国南、北方区域平均的气温距平序列之间的相关系数分别为0.62和0.52，与MV-EOF2的南、北方特征向量序列的相关性分别高达0.78和0.76 ( 表1 )，均远远超过95%的显著性水平。“三叶型”极涡振荡的三个关键区中，北欧区域指数与中国冬季气温冷暖交替振荡的相关性最高，其与北方和南方MV-EOF2特征向量序列的相关性分别为0.64和0.69 ( 表1 )，表明北极极涡形态的前后冬振荡可能通过北极极涡强度和形态的变化影响极区冷空气向南爆发的强度和范围，进而影响中国气温在前后冬的冷暖交替振荡。其中北欧地区中纬度阻塞高压异常向极入侵和南撤造成的极涡振荡起着主要贡献。对于前冷后暖年，“三叶型”极涡振荡指数与中国南北方气温的冷暖交替振荡依然存在显著的相关关系，但相关系数均较前暖后冷年有所降低。极涡振荡指数与中国北方和南方区域平均的气温距平序列之间的相关系数分别为0.4和0.38，与北方和南方MV-EOF2的特征向量序列之间的相关性分别为−0.54和−0.52。与前暖后冷年不同的是，在前冷后暖年，白令海–阿拉斯加区域的季节内振荡与我国冬季气温冷暖交替变化的关系最密切，其分别与南、北方区域平均气温距平序列的相关系数为0.36和0.32，与南、北方MV-EOF2特征向量序列的相关性分别为−0.42，−0.44 ( 表1 )，表明造成中国冬季前冷后暖年的极涡振荡主要是由白令海–阿拉斯加区域的异常高压脊的南退–北进的前后冬振荡影响的。

本文利用站点观测资料和再分析资料，从次季节时间尺度上探究了1980~2020年我国冬季气温变化的时空特征，并从北极极涡振荡的角度解释了我国冬季气温冷暖交替的成因，现得结论如下：

(1) 通过改进的MV-EOF时时矩阵分解的第二模态结果，发现1980~2020年中国冬季气温存在季节内的冷暖交替现象，前冬和后冬季节内冷暖交替的时间转折点为1月13日，即前冬为12月1日至1月13日；后冬为1月14日至2月28日。

(2) 本文从北极极涡振荡的角度解释了1980~2020年中国冬季季节内冷暖交替的成因。通过回归分析等方法定义了“三叶型”极涡振荡指数，从而发现“三叶型”极涡振荡与我国冬季气温前后冬反转存在密切联系。冬季极涡主要受到北欧区域、白令海–阿拉斯加区域以及北美区域的振荡影响，使得中国前、后冬气温出现季节内振荡的反转特征。前暖后冷年，北极涛动处于正位相，北极极涡总体上异常加强，前冬比后冬极涡更弱。北极极涡呈“三叶型”极涡振荡正位相，即前冬向极收缩后冬向南扩张，前冬北欧高压脊向北入侵北极，使得极涡形态相对狭长，从而减少了极地冷空气对中国的影响，前冬气温升温。而后冬北欧高压脊的强度减弱且南退，极涡短轴向南扩张，极地冷空气南下入侵中国，导致后冬降温。而前冷后暖年则是相反情况，北极涛动呈负位相，北极极涡总体上异常减弱，后冬比前冬极涡更弱。“三叶型”极涡振荡为负位相，即前冬极涡向南扩张后冬极涡向极收缩，极涡主要受到后冬白令海–阿拉斯加脊异常增强的影响，使得极涡横轴由前冬的向南扩张转为后冬向北缩小，对应中国的前冬降温转为后冬升温。

本文改进的MV-EOF时时分解方法可以客观地区分出中国冬季季节内冷暖交替的主模态年际变化特征，这种冷暖冬交替的发生是否有前期预报因子还值得继续深入探究。本文做过初步尝试，用PC2分别回归了冬季之前的春、夏、秋季节的海表温度和北极海冰密集的变化特征，并没有发现非常显著的前期外强迫异常信号，这对寻找前期稳定的预报因子带来了困难，未来可能需要从大气内部的环流异常变化寻找前期预报因子。注意到PC2高、低值年合成的北极极涡形态均为非纬向对称的“纺锤形”特征，可能意味着极涡边缘中纬度大气长波呈现为2波结构时，易于产生极涡振荡，进而引起中国冬季前后冬冷暖交替的振荡特征。在未来全球变暖北极放大背景下，极涡形态的变化和中纬度Rossby波动力特征的变化机理，可能是未来研究中国冬季季节内冷暖交替机制的可能途径。

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

^*通讯作者。