1961—2018年横断山区寒潮时空变化及其影响因素①

任钇潼， 程清平，2，3， 金韩宇

1.西南林业大学 地理与生态旅游学院，昆明 650224；2.西南生态文明研究中心，昆明 650224；3.中国科学院 西北生态环境与资源研究院玉龙雪山冰冻圈与可持续发展国家(云南省)科学野外观测研究站，兰州 730000

近几年全球气候加速变暖，极端气候事件风险进一步加剧[1].寒潮作为我国主要的天气过程之一，伴随的急剧降温所触发的极端低温事件严重影响国民经济及人民生命财产安全[2-3].长期以来，国内学者对寒潮进行了大量研究工作，李宪之[4]、陶诗言[5]、丁一汇[6]等研究寒潮源地、发展机理与侵袭路径，形成了比较系统和经典的寒潮理论[7]；王遵娅和丁一汇[8]、魏凤英[9]、朱晨玉等[10]，从全国尺度对寒潮进行统计分析并分别指出大气环流异常(诸如西伯利亚高压减弱、北极涛动位相增强)对国内寒潮减少有显著影响.当前，长江中下游[11]、青藏高原[12]、四川盆地[13]、河西走廊[14]、云南[15]等区域相关寒潮研究均进一步揭示了寒潮的空间异质性以及大气环流异常与寒潮的显著关联.

横断山区气候敏感，容易受到青藏高原影响，是东南亚极为重要的生态廊道——纵向岭谷区的重要组成部分[16-17].在全球气候变暖的背景下，研究横断山区历史寒潮事件的时空特征，识别寒潮与大尺度环流指数的关联具有重要的现实意义，可以为寒潮监测预警等气候业务提供参考[10].

1 数据来源与方法

1.1 研究区域概况

横断山区(24.68°-34.33°N，97.32°-104.45°E)位于青藏高原东南侧(图1)，包括川滇西部和西藏自治区东部，系南北走向山脉，海拔329～6 304 m，平均海拔3 000 m.横跨亚热带和高原温带的地理条件导致区内的水热条件空间异质性大.

审图号：GS(2020)4619号

1.2 数据来源

横断山区内66个气象站点1961—2019年的逐日最低气温数据来源于国家气象科学数据中心的中国地面气候日值数据集(http：//data.cma.cn/)；为研究寒潮对大尺度环流异常的响应所选取的北极涛动指数(AO)、大西洋数十年震荡指数(AMO)、太平洋年代际振荡指数(PDO)、多变量ENSO指数(MEI)、东大西洋-西俄罗斯遥相关指数(EA/WR)、太平洋-北美遥相关型指数(PNA)、西藏高原-1指数(TPR1)、西藏高原-2指数(TPR2)、西太平洋副热带高压面积指数(WPSHI)、西太平洋副热带高压强度指数(WPSHII)、西太平洋副热带高压脊线指数(WPSHLI)及西太平洋副热带高压西伸脊点指数(WPSHPI)来自中国气象局国家气候中心发布的大尺度环流指数集(http：//ncc-cma.net/cn/)，指数定义详见http：//ncc-cma.net/Website/index.php?ChannelID=43&WCHID=5.

1.3 研究方法

1.3.1 寒潮标准

依据《寒潮等级》(GB/T21987-2017)确定横断山区寒潮标准：将某地日最低气温24 h内降温幅度大于或等于8 ℃，或48 h内降温幅度大于或等于10 ℃，或72 h内降温幅度大于或等于12 ℃，并且使该地日最低气温大于或等于4 ℃的冷空气活动定义为寒潮[18].

以研究区1961—2019年逐日最低气温为基础，分析寒潮频次及初终日的时空变化特征.由于寒潮主要发生在春季、秋季及冬季，本研究将该年9月1日到次年5月31日定义为一个寒潮年.寒潮初日定义为一个寒潮年内初次寒潮的日期；寒潮终日定义为一个寒潮年内末次寒潮的日期.采用儒略日换算方法确定寒潮的初日(寒潮年内首次出现寒潮的日期)、终日(寒潮年内最后出现寒潮的日期)及持续日数，如：以1月1日为起点记为1，2月1日记为32，即将寒潮出现的初日、终日日期转换为距离1月1日的实际日数，得到寒潮初日、终日的时间序列.

1.3.2 趋势特征识别

采用Sen斜率[19]结合修订的Mann-Kendall[20]以及Pettitt突变检验[21]辨识寒潮时空变化特征.

1.3.3 大尺度环流因子检测

采用小波相干分析(Wavelet Transform Coherence，WTC)[22]探测寒潮与大尺度环流指数的共振周期，剖析寒潮与大气环流异常的联系.厄尔尼诺/拉尼娜-南方涛动(El Nio/La Nia-Southern Oscillation，ENSO)是造成全球气候年际变化的最重要的海气耦合现象，对中国气候影响显著，因此单独剖析ENSO下寒潮变化特征对于揭示极端事件的驱动因素具有重要意义[23].ENSO事件的判定采用NOAA地球系统实验室发布的多变量ENSO指数(MEI)双月等级排序数据集进行：若MEI等级排序为1～21，视为一次拉尼娜现象；若MEI等级排序为48～68，视为一次厄尔尼诺现象.

2 结果与分析

2.1 寒潮时间特征分析

由图2可知，横断山区全年寒潮呈显著减少态势(每10年减少0.05次)，其中又以冬季最为显著(每10年减少0.04次)，春季次之(每10年减少0.03次)，秋季最低(每10年减少0.02次)，即冬季寒潮的显著减少主导了区内寒潮变化趋势.Pettitt突变检验发现，横断山区于1983年发生显著突变，1983年以前寒潮年均频次为0.95次/a，1983年以后寒潮年均频次0.78次/a.

注：“*”表修订的Mann-Kendall趋势检验统计值绝对值|Z|≥1.96

2.2 寒潮空间变化特征

空间上横断山区全年单站寒潮(图3a)表现为北多南少，以石渠、红原、色达及若尔盖站最多(≥3次/a).81.82%的站点寒潮趋于减少，以德格、稻城、美姑、香格里拉以及弥渡站最为显著，且沿金沙江峡谷以及云岭呈“V”形分布；少数站点寒潮趋于缓慢增加，多分布在横断山南部.

春季(图3b)，石渠、红原、色达以及若尔盖站寒潮多发(≥1 次/a).68.19%的站点寒潮趋于减少，以若尔盖、红原、色达、稻城以及香格里拉站最为显著，多分布在横断山北部和中部；18.19%的站点寒潮趋于缓慢增加，其余站点寒潮无明显变化趋势.秋季(图3c)，51.52%的站点寒潮呈减少态势，以色达、松潘以及香格里拉站最为显著；33.33%的站点寒潮频次无明显趋势，多分布于横断山区东西两侧；15.16%的站点寒潮缓慢增加，多分布在横断山南部.冬季(图3d)，60.61%的站点寒潮趋于减少，以若尔盖、色达、香格里拉站最为显著，此外有36.37%的站点寒潮频次趋于增长，多分布在“若尔盖-兰坪站”以东.

中国分类号：GS(2020)4619号

2.3 寒潮对大尺度环流异常的响应特征

2.3.1 WTC分析

图4可知，1961-2018年，EA/WR与寒潮在11.69～15.34 a的周期上呈稳定的显著负相关，1996-2008年，在3.68～5.84 a的周期上寒潮显著提前于EA/WR；MEI于1970-1877年、1984-1988年分别在5.21～7.80 a、2.19～2.46 a的周期上与寒潮呈显著正相关.PNA，TPR1以及TPR2在1974-1989年分别在9.28～12.38 a、6.95～10.41 a以及6.59～11.03 a的周期上与寒潮呈显著负相关；WPSHLI在1965-1980年在2.19～2.92 a的周期上与寒潮呈显著负相关，1980年以后表现为在更低周期上的显著负相关.

注：黑色粗轮廓表示通过α=5%红噪声检验；从左向右箭头表示寒潮事件与大尺度环流指数正相关，从右向左箭头表示寒潮事件与大尺度环流指数负相关，垂直向下箭头表示寒潮比大尺度环流指数提前90°，垂直向上的箭头表示大尺度环流指数比寒潮事件提前90°

AMO，AO，PDO，WPSHI，WPSHII，WPSHPI与区内寒潮仅出现短暂的一致性，其中AO，PDO在1996-2004年分别在2.32～3.28 a、2.07～2.15 a的周期上显著提前，滞后于寒潮.

2.3.2 厄尔尼诺下寒潮事件变化特征

如表1所示，总体上，厄尔尼诺和拉尼娜下均以冬季寒潮日所占比重最大(25.88%～39.42%).1961-1980年，拉尼娜较厄尔尼诺更为频繁，拉尼娜下寒潮日比重最高(54.5%～66.04%)；1980年以后，厄尔尼诺较拉尼娜更为频繁，且以厄尔尼诺下寒潮日比重最高(52.56%～53.06%).此外，厄尔尼诺下，春、秋、冬三季的寒潮日比重均趋于上升，拉尼娜下春、秋、冬三季的寒潮日比重均趋于减小，且趋势变化程度均以冬季最高.

表1 厄尔尼诺和拉尼娜下季节寒潮累积日数年代际变化

2.4 寒潮持续日数与初终日统计

2.4.1 寒潮持续日数统计特征

图5a表明单站寒潮累积日数由北向南呈“阶梯式”减少，以32°N最为明显.32°N以北，1961-1980年单站寒潮日相对偏多，1980-2018年偏少但波动强烈；32°N以南，情况相反，即1961-1980年寒潮日偏多，1980-2018年偏少，以28°N地区最为明显.图5b表明，横断山区单次单站寒潮事件可持续1～15 d，但以短持续日寒潮(1～5 d)为主.短持续日寒潮、中持续日寒潮(6～10 d)在1961-2000年趋于减少，2000年后趋于增加.长持续日寒潮(11～15 d)属极小概率事件，在1960-2018年未出现明显增减.

图5 1961-2018年横断山区寒潮结构特征

2.4.2 寒潮初终日统计特征

如图6a所示，横断山区寒潮初日(10月4日至次年5月7日)表现为西北提早南部滞后(1月5日至4月6日).全区56.06%的站点寒潮初日整体趋于提早(每10年提早2 d)，仅弥渡站(每10年提早3.33 d)通过0.05显著性检验；43.94%的站点寒潮初日滞后，其中以香格里拉(每10年滞后4.62 d)、宁蒗站(每10年滞后7.96 d)最为显著.如图6b所示，全区寒潮终日区域差异性较小，即各站寒潮终日较于初日均较为同步(1月10日-3月16日).66.67%的站点寒潮终日趋于提早，其中以昭觉(每10年提早3.70 d)、色达(每10年提早5.48 d)、香格里拉(每10年提早13.8 d)、稻城(每10年提早7.8 d)、喜德(每10年提早0.32 d)、弥渡(每10年提早3.75 d)、美姑站(每10年提早2.61 d)最为显著.其余站点寒潮终日推迟，但均未通过0.05显著性检验.

图6 1961-2018年横断山区寒潮初日及终日时空变化特征

3 讨论与结论

3.1 讨论

本研究发现横断山区寒潮多发于冬季，与王遵娅等[8]所描述的南方寒潮多发于春季的情况不符，这可能与对寒潮的判定标准不同以及国内气温趋于升高[24]导致春季寒潮减少有关.研究区内寒潮于1983年发生显著突变，与朱晨玉等[10]研究得出的国内寒潮突变时间点接近；研究区寒潮趋于减少，与长江中下游[11]、青藏高原[12]、云南[15]地区寒潮变化趋势一致，与四川盆地[13]趋势相反.研究区内寒潮对EA/WR的响应最强，对AO的响应较弱，这与王遵娅等[8]的AO对寒潮具有显著影响的结论不一致.

全球气候变暖，极端低温事件发生的频率和强度与空间范围变化显著[25-26]，因此需加强相关领域的监测分析.尽管横断山区寒潮趋于减少，但寒潮初日、终日趋于提早或对春播作物产生严重影响，因此需加强春季寒潮的监控预防.此外，横断山区是重要生态屏障，是生物的重要庇护所[16]，需进一步从频率、强度等方面对寒潮作更深入的研究.

3.2 结论

1)1961-2018年研究区寒潮趋于显著减少(每10年减少0.05 次)，但持续日上仍以1～5 d寒潮为主；受纬度和海拔的影响，研究区内寒潮北多南少，以32°N为重要分界线.

2)弥渡、宁蒗、香格里拉站的寒潮初日，昭觉、色达、香格里拉、稻城、喜德、弥渡以及美姑站的寒潮终日均已发生显著变化.

3)EA/WR，MEI，PNA，TPR1，TPR2，WPSHLI与寒潮的相关性最为显著；厄尔尼诺下寒潮比重趋于增加.