引言

政府间气候变化专门委员会第五次评估报告(IPCC，2014)指出，自上世纪中期以来，全球年平均气温上升了0.8℃左右，尤以冬季北半球中高纬地区增暖最明显。因而暖冬成为普遍关注的一个热点问题，但在暖冬里，寒冷的音符并未休止，21世纪以来伴随东亚冬季风再次转强(梁苏洁等，2014)，大规模寒潮、暴雪频繁侵袭欧亚和北美大陆，极端低温事件频发，全球气候似乎又进入一个新的波动阶段。

值得注意的是，近年来一些观测事实表明，中国冬季气温季节内差异明显，前冬和后冬冷暖异常急转的现象频繁出现：如2001年前冬(2001年12月)全国一致偏冷，后冬(2002年1—2月)转为偏暖；2007/2008年的情况则与2001/2002年相反：前冬整体偏暖，而后冬异常偏冷，并发生了历史罕见的低温雨雪冰冻灾害(丁一汇等，2008；王东海等，2008；布和朝鲁等，2018)；2014/2015年前、后冬出现冷-暖转折以及2015/2016年前、隆冬出现暖-冷急转现象(司东等，2016；聂羽等，2016)等。据统计，前冬出现全国性偏暖(冷)或北冷(暖)南暖(冷)的气温异常时，而后冬出现类似和相反气温异常的概率分别在50%左右(韦玮等，2014)，从预测角度来讲具有较大不确定性，给实际业务会带来一定困扰和难度。与此同时，北半球前、后冬气温变化趋势也呈现明显非对称性(董文杰等，2003；Cohen et al，2014)：前冬呈增温趋势，而后冬则为中性或降温的反转趋势；黄嘉佑和胡永云(2006)也曾指出，中国内陆前冬和后冬的气温变化也存在不同的趋向性，且这种趋势与整个冬季平均的趋势是不同的。Wei et al(2011)利用月平均冬季风强度指数的隔月差定义季风变化的倾向指数(MTI)，描述冬季风演变过程，1月前为正而进入1月后为负，对应中国气温倾向在前冬一致变冷，后冬一致变暖，这同北半球前后冬趋势正好相反。同时，前、后冬气温的这些非对称性变化对整个冬季增暖趋势具有一定调控作用。

以上情况表明，我国冬季气温变化存在季节内转变的可能性，且这种转折极易造成极端气候事件的发生。然而，通常科研和业务上主要是针对整个冬季的平均气温，这存在局限性并会掩盖一些问题，自然季节的这种季内转折特征在整个冬季平均场上无法体现。随着社会不断发展，针对整个冬季平均气温的预测已逐渐不能满足各行业的实际需求，季节内尺度如次季节预测正越来越受到广泛关注(Osman and Alvarez, 2018；Liang et al，2018；蓝柳茹和李栋梁，2019；李思和武炳义，2019)。中国冬季气温变化在不同时间尺度受不同因子影响(康丽华等，2006)，此前已得到广泛关注和深入研究，而关于季节内尺度方面的研究较少，一方面是对气温季节内尺度变化的认识还存在一定局限，对其相关机制及可预报性来源的研究尚无定论，另一方面已有研究中发现的影响因子间复杂的非线性相互叠加且在季节内尺度上存在较大的不稳定性(封国林等，2008；何溪澄等，2008)。

气温的变化不仅包括均值的演变，还包括均方差、极值等统计参数的变化。以往研究多针对冬季平均气温的年际、年代际变率特征展开研究，为冷暖冬预测提供极其重要的科学参考价值，但对于冬季气温季节内变率方面研究较少，如季节内气温波动和变化稳定性等。均方差作为描述变量变化离散程度的重要统计量，代表均值变化的稳定性及强弱程度，均方差越大表明均值稳定性越弱，极端性越强，因此均方差可用来衡量要素变化的稳定性情况，进而反映变率特征。本文利用气温均方差这一统计量来刻画气温变率，分析在气候变暖背景下中国冬季气温季节内变率特征，并进一步探讨其影响因子及可能机制，突破仅从季节平均角度认识冷暖冬异常的局限性，尝试从另一角度认识冬季气温的变率特征及其机制，以期为我国冬季气温季节内预测提供新的视角和参考依据，进一步提升冬季气候预测的精细化水平。

1 资料和方法 1.1 资料说明

站点资料来自国家气象信息中心的台站观测数据集，选用序列完整、可靠性较高、分布均匀的中国364个站冬季逐日平均气温；格点资料选用美国国家环境预测中心(National Centers for Environmental Prediction, NCEP)提供的风场、海平面气压场和位势高度场再分析资料(Kalnay et al, 1996)，水平格距为2.5°×2.5°。资料起止时间均为1961—2018年，冬季指12月至次年2月。

1.2 研究方法

本文利用冬季平均气温均方差作为度量冬季气温季节内变率的指标，均方差值越大，表明冬季平均气温稳定性越弱，其高频波动越强；均方差值的年际变率越大，冬季平均气温变化趋势的稳定性也越弱，因而气温均方差可反映出冬季气温冷暖波动及其变化趋势的平稳性。

针对1961—2018年冬季全国364个站逐日平均气温序列，利用公式(1)计算历年各站平均气温均方差序列σ_jt(j=1，2，…，s；t=1，2，…，n)：

$ \sigma_{j t}=\sqrt{\frac{1}{m} \sum\limits_{i=1}^{m}\left(x_{i j t}-\bar{x}_{j t}\right)^{2}} $

(1)

式中：j为测站数，t为年份，i为冬季日数(i=1, 2, …, m)，x_ijt表示逐日平均气温，x_jt为整个冬季的平均气温。

文中需进一步分析气温变率和环流变率的联系，冬季环流场的均方差也参考上述方法计算，用来表示大气环流季节内变率的强弱。不同尺度环流分量信号采用Butterworth带通滤波器对逐日环流场进行滤波，分离出冬季不同尺度的环流分量，本文主要提取了2~10、10~30、30~60和60~90 d四种尺度分量，代表天气尺度、延伸期尺度、次季节尺度和季节尺度信号，再并根据式(1)计算出不同尺度环流场的均方差，代表不同尺度环流分量的季节内变率。

2 中国冬季平均气温季节内变率特征及环流分析 2.1 冬季气温季节内变率特征

根据式(1)计算全国364个站1961—2018年冬季气温均方差场，首先对该场在空间维平均，即对全国范围进行空间平均后得到时间序列，该序列的变化可反映中国冬季气温季节内变率的时间演变特征(图 1a)，全国平均气温均方差总体呈线性减弱趋势，表明近几十年来气温季节内变率逐渐趋弱。同时也表现出明显的年代际特征，20世纪80年代中期以前，处于变率偏强时期，意味着中国冬季气温冷暖波动频繁，处在稳定性偏弱的阶段；20世纪80年代后期至21世纪初，明显转弱，表明此阶段冬季气温季节内稳定性增强；2004年以后，再次增强。以上年代际变化同东亚冬季风年代际转折时间相匹配，特别是20世纪80年代后期东亚冬季风偏弱至2005年再次转强(梁苏洁等, 2014; Wang and Chen, 2014)，这与气温季节内变率的转折期相吻合。由此表明，东亚冬季风年代际变化也可能会影响中国冬季气温的季节内变率，对气温波动及强弱程度产生影响，进而影响气候变暖的稳定性。此外，气温季节内变率的年际振荡较明显，尤其在近十年，这与近年冷暖冬交替现象频繁也具有一定联系。与此同时，对气温均方差场在时间维平均，得到全国各站的气温均方差的多年平均值，分析其空间特征(图 1b)，发现中国北方大部和长江以南地区冬季气温季节内变率明显大于其他地区，尤其在东北和新疆北部地区(变率大值区)，这是由于我国北方易受冷空气频繁侵袭所致。值得注意的是，长江以南存在另一变率高值区，这可能同东亚冬季风“北弱南强”模态(彭京备和孙淑清，2017)存在一定联系。变率低值区主要位于中国西南地区，表明该地区冬季气温稳定性较高，高频波动相对较弱。

图 2a为冬季平均气温及气温均方差变化趋势，冬季我国大范围地区平均气温呈增加趋势，尤其青藏高原地区最明显，增温速率达0.4 ℃·(10 a)^-1以上，而华南和西南东部则呈变冷趋势。在此增暖背景下，均方差趋势系数分布显示全国为一致减弱趋势，尤其东北地区的减弱趋势最明显。根据各站平均气温和均方差的相关图来看，全国大范围地区为负相关区，显著区域位于我国北方地区，在新疆和东北地区表现最突出，说明北方冬季平均气温异常对气温季节内变率具有较好的指示性，气候增暖(变冷)则意味着季节内变率减弱(增强)。

2.2 影响气温季节内变率的环流异常特征 2.2.1 冬季平均环流场

大范围或局地地面气温异常与大尺度平均环流异常具有直接联系，季节平均环流异常是季内气温变率的背景场。图 3给出了1961—2018年中国冬季气温均方差对同期冬季地表气温、位势高度场及海平面气压场的相关分布，当气温均方差偏强时，我国北方及中东部地区气温显著偏低，而我国西南、华南等区域气温略微偏高，意味着冬季西伯利亚地区气温异常对我国北方冬季气温高频波动具有明确的指示意义，这是由于冷空气易在西伯利亚地区聚集加强后南下沿偏东路径影响北方地区(图 3a)。高度场上自高层至低层呈类似准正压结构，平流层表现为极涡明显偏弱，北大西洋、北美地区和欧亚中纬度上空为显著负相关(图 3b)，研究表明，当平流层异常强的极涡向下传播时，引起低层北极涛动(Arctic Oscillation，AO)正异常发展，西伯利亚高压增强，导致东亚沿岸东西气压梯度增加，有异常偏强的季风(Chen and Zhou, 2012；Li and Yang, 2010)，我国南部地区温度易偏低，而东北和亚洲中高纬度地区温度增加(Wang and Chen, 2014)。反之，当异常弱极涡向下传播时，会引起大致相反的过程，这与上述均方差偏强时平流层极涡偏弱特征一致，因而平流层弱极涡信号下传时，可能会导致我国北方地区冬季气温季节内变率增强，易出现极端气温事件；对流层中层(图 3c)，自西欧—东亚呈现“-+-”型波列分布，其中大西洋东北部和贝加尔湖附近为负异常中心，斯堪的纳维亚上空为明显正异常，这种水平环流异常同SCA(Barnston and Livezey, 1987)正位相特征非常一致，这种分布型可为极地冷空气向南爆发影响中纬度地区提供有利的大尺度背景；海平面气压场上，极地同北大西洋、欧洲、北美地区反向，表明当气温均方差偏强时，对应北大西洋涛动(North Atlantic Oscillation, NAO)负位相特征，同时西伯利亚地区为显著正相关，表明西伯利亚高压偏强(图 3d)。

通过大气遥相关型和异常行星波活动等途径，冬季NAO的影响可扩展到东亚地区，NAO负(正)位相时对应乌拉尔山阻塞高压偏强(偏弱)，以及东亚大槽加深(变浅)(Chen and Zhou, 2012)。值得注意的是，冬季期间NAO与东亚气候之间的联系存在明显的季节内差异，Watanabe and Kimoto(2010)指出前冬NAO的信号局限在欧洲—大西洋区域，但后冬NAO的信号可扩展到东亚和北太平洋。中国位于NAO相关环流系统下游，谭桂容等(2010)研究指出，2008年1月中国气温异常与同期北大西洋西风急流和NAO遥相关环流异常存在紧密联系。2009/2010年冬季，中国北方气温异常与NAO负位相密切相关(孙诚和李建平，2011)，那么NAO是否与中国冬季气温季节内变率同样具有密切联系？通过分析，冬季NAO与中国冬季气温均方差的相关系数为-0.496，通过0.01显著性水平检验。其中显著负相关区域主要位于我国东部地区，表明冬季NAO为负位相时上述地区气温季节内变率增强，气温稳定性偏弱，易发生极端气温事件，由此说明冬季NAO的位相对我国冬季气温季节内变率也具有较好的指示性。

相关研究还发现，NAO/AO可通过影响东亚大槽、西伯利亚高压、东亚冬季风等系统，使其变得更加活跃，使中国冬季大部分地区气温偏低，进而导致极端低温频发(武炳义等，1999)。气温高频波动、稳定性减弱的直接因素是环流系统波动异常产生的天气过程，因而NAO可能通过相关环流的高频变化来影响气温季节变率。通过计算冬季逐日海平面气压的方差，代表近地面大气环流的季节内变率，反映其高频变化和稳定性。根据中国冬季气温均方差同海平面气压(SLP)方差的相关场(图 5a)显示，亚速尔群岛、赤道大西洋地区、北极—西伯利亚—贝加尔湖及以南大范围地区呈明显正相关区，且一直向南延伸至中国北方及东部地区，其显著中心区与西伯利亚高压活动中心区一致，冰岛附近为显著负相关区，意味着NAO负位相偏强时，有利于西伯利亚高压、冷涡系统活跃，易引起我国气温波动异常。提取位于西伯利亚中心关键区(45°~65°N、65°~120°E)范围内各格点SLP均方差的区域平均值，代表西伯利亚高压中心区的季节内变率，中国冬季气温均方差与该中心区变率相关系数高达0.578，显著通过0.001显著性水平检验，表明冬季西伯利亚高压的高频变化与我国冬季气温季节内变率具有紧密联系。此外，NAO同SLP变率场的关系也很好反映出西伯利亚地区的高频波动对NAO异常的响应(图 5b)，二者呈显著负相关关系，即NAO负位相增强时，西伯利亚至贝加尔湖及以南地区季节内变率偏强。可见，中国冬季气温季节内变率异常，很大程度上与NAO相关环流系统的季节内变率有关，可能是NAO通过影响西伯利亚高压天气尺度上的高频波动进而引起中国冬季气温的季节内变率异常。

2.2.2 不同时间尺度环流分量的影响

为了揭示大气内部变率对气温季节内变率的影响，利用逐日再分析资料，选取海平面气压场和位势高度场两个物理量，采用Butterworth带通滤波方法提取冬季内包括2~10、10~30、30~60和60~90 d四种尺度环流分量，接着计算各尺度环流分量的均方差，用来反映冬季不同时间尺度上大气内部环流分量的季节内变率。首先，1961—2018年中国平均冬季气温均方差同海平面气压场的不同尺度分量变率关系显示，西伯利亚—贝加尔湖—蒙古—我国北方均存在显著正相关区，西伯利亚高压在不同时间尺度上的波动异常均会对中国冬季气温季节内变率产生影响，因而近地层西伯利亚附近环流变率异常尤其重要。在此提取海平面气压场关键区(45°~65°N、65°~120°E)平均变率，分析气温均方差与该关键区四种尺度平均变率的关系(图 6)，均显著通过0.01显著性水平检验，同次季节尺度环流变率的相关性最高，延伸期尺度次之，而天气尺度则最弱，因而冬季内次季节尺度上的西伯利亚高压变率异常对我国气温变率起到决定性作用。

与此同时，在天气尺度分量场上(图 6a)，鄂霍次克海—白令海地区、中南半岛及南海地区也出现显著正异常区，而显著负异常区位于北美东部至西欧地区，表明阿留申低压在天气尺度上的频繁波动及上游欧洲脊的稳定少动与气温变率联系紧密；延伸期和次季节尺度分量场(图 6b，6c)均显示为显著正异常特征，异常关键区主要集中在贝加尔湖及以南延伸至中国南海及赤道地区，北大西洋地区也存在显著异常区，其中31~60 d尺度分量变率的影响范围更集中、空间尺度亦更大。季节尺度分量的显著异常区主要集中在新地岛和西伯利亚—巴尔喀什湖—贝加尔湖及以南的蒙古高原—中国北方(图 6d)，体现出冷空气源地和西伯利亚高压在较长时间尺度上的波动对气温变率的重要影响。

图 7为冬季500 hPa高度场上各尺度分量变率与中国冬季气温变率的相关场，在2~10 d尺度上，北大西洋—欧洲—亚洲中西部地区为显著负相关，北大西洋和北太平洋均呈“+-+”三极型分布，表明中国冬季气温季节变率偏大时，对应欧亚中高纬阻塞系统变率减弱或稳定维持，而鄂霍次克海阻塞系统变率增强或易波动；11~30 d尺度上，格陵兰岛、西伯利亚及以北、西北太平洋地区和中南半岛及以南地区呈显著正相关，反映中高纬地区波列异常及东亚大槽的高频波动所带来的影响；其余两种尺度的影响较相似，主要体现为欧亚阻塞异常的影响，如乌拉尔山地区为正相关，表明乌拉尔山阻塞高压在月—季节尺度上的变率对冬季气温高频波动具有重要影响。

3 结论与讨论

本文从冬季气温季节内变化的角度，分析了我国冬季气温季节内变率和稳定性特征，在时间上总体呈线性减弱趋势，表明近几十年来气温季节内变率逐渐趋弱，从而季节内气温变化稳定性增强，同时气温季节内变率还具有明显的年代际特征，20世纪80年代中期以前，变率偏强，冬季冷暖波动频繁，稳定性较弱，80年代后期至21世纪初，变率明显转弱，2004年以后，变率再次增强，这同东亚冬季风年代际转折时间相吻合，因而东亚冬季风年代际变化可能会对中国冬季气温的季节内波动的强弱程度产生影响，进而影响气候变暖的稳定性。就空间上来看，中国北方大部和长江以南地区变率偏大，是由于我国北方易受冷空气频繁侵袭所致，长江以南高值区可能同东亚冬季风“北弱南强”模态存在一定联系。中国西南地区变率较小，表明该地区冬季气温稳定性较高，高频波动相对较弱。当气温季节内变率异常偏强时，冬季平均环流场上呈类似准正压结构，平流层极涡偏弱，对流层中高纬呈类似斯堪的纳维亚遥相关型分布，中低纬地区为大范围负异常区，近地层上西伯利亚高压偏强，NAO为负位相；NAO同我国东部气温变率联系密切，进一步分析揭示出NAO是通过影响西伯利亚高压的高频变化来作用于气温季节内变率。最后，通过提取天气—季节四种时间尺度上的大气环流内部变率，发现在各个尺度上，气温季节内变率均受西伯利亚高压和东亚冷涡的调控作用；此外，在天气尺度上，阿留申低压频繁波动及上游欧洲脊的稳定少动与气温变率有密切联系，季节尺度上包括欧亚阻塞和鄂霍次克海阻塞异常对气温变率有显著影响。

关于气温季节内变率，主要从中高纬分析了可能的影响，那么热带地区是否也会对我国冬季气温季节内变率产生一定影响，比如典型ENSO背景下，是否会通过热带环流异常进而影响到气温季节内变率，中高纬和热带的影响作用是否会存在一定的年代际变化，以上问题值得继续探讨。