王婧 呂俊梅

中國(guó)氣象科學(xué)研究院災(zāi)害天氣國(guó)家重點(diǎn)實(shí)驗(yàn)室，北京 100081

1 引言

觀測(cè)表明，在全球變暖背景下，北極地區(qū)顯著增暖，近二十年歐亞大陸卻頻繁遭受冷空氣侵襲，出現(xiàn)大范圍降溫趨勢(shì)，引起廣泛關(guān)注（Zhang et al.,2008;Petoukhov and Semenov,2010;Overland et al.,2011;Inoue et al.,2012;Mori et al.,2014;Cohen et al.,2020）。在這種氣候背景下，2008年1月中國(guó)南方發(fā)生大范圍罕見(jiàn)低溫雨雪天氣，造成近一億人口受災(zāi)，幾億畝農(nóng)田凍害及幾十萬(wàn)間房屋倒塌，交通運(yùn)輸受阻，造成直接經(jīng)濟(jì)損失近千億（丁一匯等,2008;顧雷等,2008;李崇銀等,2008）。2009/2010年冬季歐洲多地遭受暴雪天氣，新疆、內(nèi)蒙古等地同樣發(fā)生暴雪天氣，致使幾十人喪生（李如琦等,2015;王遵婭和周波濤,2018）。2012年1～2月，歐亞大陸遭受極寒天氣侵害，俄羅斯等地出現(xiàn)近百年最低氣溫，許多地區(qū)發(fā)生暴雪天氣，嚴(yán)重影響交通及社會(huì)秩序（WMO,2012;韓哲等,2014）。歐亞大陸冷冬頻繁出現(xiàn)及寒潮、暴雪災(zāi)害頻發(fā)可能造成農(nóng)作物凍害、電線積冰、路面積雪以及能見(jiàn)度降低，嚴(yán)重影響通信與交通運(yùn)輸，危害公共安全，阻礙經(jīng)濟(jì)發(fā)展并造成嚴(yán)重經(jīng)濟(jì)損失。

近年來(lái)科學(xué)家們將北極增暖而歐亞大陸變冷的地表氣溫模態(tài)稱為“暖北極—冷歐亞”模態(tài)（Warm Arctic–Cold Eurasian，簡(jiǎn) 稱WACE;Mori et al.,2014;Luo et al.,2016;Sung et al.,2018），也有科學(xué)家將其稱為“暖北極—冷大陸”或“暖北極—冷西伯利亞”模態(tài)（Overland et al.,2011;Inoue et al.,2012;Sorokina et al.,2016），并且科學(xué)家們?cè)谀觌H時(shí)間尺度上廣泛地探討了WACE模態(tài)形成的原因。有的研究認(rèn)為，受到全球變暖影響，北極海冰減少導(dǎo)致WACE形成（Inoue et al.,2012;Tang et al.,2013;Mori et al.,2014;Kug et al.,2015）。北極海冰影響WACE模態(tài)的物理機(jī)制為：北極海冰融化，有利于北大西洋濤動(dòng)處于負(fù)位相且東亞大槽加深，烏拉爾阻塞發(fā)生頻率增加，造成歐亞大陸溫度偏低（Inoue et al.,2012;Liu et al.,2012;Mori et al.,2014;Kug et al.,2015;Luo et al.,2016;Yao et al.,2017;Wegmann et al.,2018）。Zhang et al.（2016）從行星波上傳的角度解釋北極海冰融化影響歐亞大陸溫度的過(guò)程：巴倫支—喀拉海海冰減少，緯向一波增強(qiáng)，導(dǎo)致行星波上傳至平流層，引起平均流異常。平流層將這種信號(hào)下傳至對(duì)流層，使極渦向歐亞大陸偏移并減弱，烏拉爾山高壓脊異常發(fā)展并伴隨東亞大槽加深，導(dǎo)致冷空氣侵襲西伯利亞，進(jìn)而歐亞與北美大陸溫度在近幾十年降低。

然而，另外一些工作指出，大氣內(nèi)部變率對(duì)WACE模態(tài)的影響不可忽視。一些科學(xué)家進(jìn)行了數(shù)值試驗(yàn)，結(jié)果表明北極海冰減少及人為強(qiáng)迫沒(méi)有顯著影響歐亞大陸冬季溫度，WACE的形成源于大氣內(nèi)部自身變率（McCusker et al.,2016;Sun et al.,2016）。西風(fēng)急流減弱、西伯利亞高壓加強(qiáng)、阿留申低壓加強(qiáng)可能造成暖空氣進(jìn)入極區(qū)而極地增暖，冷空氣南侵，歐亞大陸溫度降低，最終WACE形成（Bengtsson et al.,2004;Zhang et al.,2008;Sorokina et al., 2016）。另外，Graversen et al.（2008）認(rèn)為，雖然海冰的融化對(duì)表層大氣的增暖有貢獻(xiàn)，但是大氣能量傳輸是北極上空大氣增暖的主要原因。Woollings et al.（2014）發(fā)現(xiàn)去線性趨勢(shì)后的巴倫支—喀拉海增暖與WACE并無(wú)顯著關(guān)聯(lián)。Sorokina et al.（2016）認(rèn)為冬季巴倫支海湍流熱通量的EOF（Empirical Orthogonal Function）第一模態(tài)與海冰減少關(guān)聯(lián)較弱，主要受到大氣內(nèi)部變率影響：西風(fēng)急流使海冰邊緣北退，巴倫支—喀拉海出現(xiàn)暖平流，湍流熱通量減弱，WACE形成。

另外，WACE同樣受到海洋內(nèi)部變率影響（Sato et al.,2014;Nakanowatari et al.,2014;Park et al.,2015;Woods and Caballero,2016;Jung et al.,2017）。北大西洋海溫正異常，激發(fā)羅斯貝波，導(dǎo)致水汽與熱量向極輸送，進(jìn)而北極增暖，進(jìn)一步改變風(fēng)暴軸與急流位置、行星波與能量傳播進(jìn)而影響歐亞大陸溫度，這對(duì)WACE的形成十分重要（Sato et al.,2014;Woods and Caballero,2016;Jung et al.,2017）。此外，海溫可能通過(guò)影響北極海冰融化進(jìn)而影響WACE（Nakanowatari et al.,2014;Park et al.,2015），即北極海冰融化對(duì)WACE的影響是氣候系統(tǒng)內(nèi)部變率的體現(xiàn)（Luo et al.,2016;Wegmann et al.,2018）。北大西洋海溫異?？赡芡ㄟ^(guò)洋流輸送至極區(qū)造成海冰融化（Nakanowatari et al.,2014）。印度洋及西太平洋海溫正異常，可能導(dǎo)致洋面上空對(duì)流活動(dòng)加強(qiáng)，進(jìn)而增加極區(qū)向下長(zhǎng)波輻射，導(dǎo)致北極海冰減少（Park et al.,2015）。

目前對(duì)于WACE年代際變化特征及成因研究較少。何金海等（2015）發(fā)現(xiàn)北極與歐亞大陸冬季溫度的反向變化在21世紀(jì)初由北極冷—大陸暖轉(zhuǎn)為北極暖—大陸冷。Sung et al.（2018）認(rèn)為羅斯貝波與西伯利亞高壓相互作用導(dǎo)致了WACE模態(tài)年代際尺度上的加強(qiáng)。另外，大西洋多年代際振蕩（Atlantic Multidecadal Oscillation，簡(jiǎn)稱AMO）與太平洋年代際振蕩（Pacific Decadal Oscillation，簡(jiǎn)稱PDO）可能導(dǎo)致了20世紀(jì)初北極增暖（Miles et al.,2014;Tokinaga et al.,2017）。

綜上所述，過(guò)去的研究主要集中于WACE年際變化物理原因的探討。然而，WACE還具有顯著的年代際變化特征（Sung et al.,2018）。目前對(duì)于WACE年代際變化特征的研究，研究時(shí)段多集中于20世紀(jì)中后期至今，并且未去除外強(qiáng)迫影響，研究結(jié)果可能包含外強(qiáng)迫貢獻(xiàn)（Sung et al.,2018）。20世紀(jì)以來(lái)WACE模態(tài)年代際變化特征及與其年代際變化相聯(lián)系的大氣環(huán)流特征仍不清楚，我們對(duì)導(dǎo)致其年代際變化的物理原因的了解還不夠深入。然而，掌握WACE年代際變化規(guī)律和形成機(jī)制，有助于我們研究歐亞冬季溫度演變規(guī)律，為極端冷事件、寒潮天氣、暴雪等事件的預(yù)測(cè)提供年代際背景。因此，本文利用第六次國(guó)際耦合模式比較計(jì)劃的多模式Historical試驗(yàn)，首先去除人為活動(dòng)等外強(qiáng)迫的影響作用，分析20世紀(jì)以來(lái)WACE模態(tài)年代際內(nèi)部變率，以及WACE不同位相對(duì)應(yīng)的大氣環(huán)流場(chǎng)差異，并進(jìn)一步探討北大西洋海溫影響WACE年代際變化的物理機(jī)制。

2 資料和方法

2.1 資料

本文使用的資料為1910～2020年美國(guó)航空航天局戈達(dá)德空間研究所（GISS）2°×2°的月平均地表氣溫?cái)?shù)據(jù)集（GISTEMP v4；Lenssen et al.,2019），數(shù)據(jù)主要由氣象站點(diǎn)及船舶、浮標(biāo)、其他傳感器得到并對(duì)城市熱島效應(yīng)進(jìn)行訂正，是目前時(shí)間長(zhǎng)度較長(zhǎng)，對(duì)地表氣溫代表性較好的觀測(cè)資料（唐國(guó)利等,2011；秦大河,2018）。海溫為1910～2014年英國(guó)氣象局哈德萊中心1°×1°的月平均資料（HadISST；Rayner et al.,2003）。此外，還使用1910～2014年美國(guó)國(guó)家海洋和大氣局—環(huán)境科學(xué)協(xié)作研究所（National Oceanic and Atmospheric Administration–Cooperative Institute for Research in Environmental Sciences，簡(jiǎn)稱NOAA-CIRES）20世紀(jì)月平均再分析資料（20th Century Reanalysis，簡(jiǎn)稱20CR；Compo et al.,2011），2°×2°網(wǎng)格要素包括：海平面氣壓、850 hPa水平風(fēng)場(chǎng)和比濕、500 hPa和200 hPa位勢(shì)高度場(chǎng)及水平風(fēng)場(chǎng)；高斯網(wǎng)格（192×94）要素包括：地表氣溫、地表凈長(zhǎng)波輻射通量、地表凈短波輻射通量、地表感熱通量、地表潛熱通量和降水。日平均要素使用500 hPa位勢(shì)高度場(chǎng)，以研究阻塞發(fā)生頻率。由于1910年以前80°N以北地區(qū)GISS地表氣溫存在較多缺測(cè)，因此研究WACE年代際變化的時(shí)段為1910～2020年，其中較小范圍的缺測(cè)采用自然鄰點(diǎn)插值法進(jìn)行插值。由于20CR大氣環(huán)流資料僅更新至2014年，因此大氣環(huán)流及海溫的時(shí)段取1910～2014年。

許多科學(xué)家利用模式數(shù)據(jù)集合平均對(duì)地表氣溫去除趨勢(shì)（Dai et al.,2015;Dong and Dai,2015;Yang et al.,2020;王岱等,2020），以期去除外部強(qiáng)迫（例如，人類活動(dòng)、太陽(yáng)輻射、火山爆發(fā)引起的氣溶膠），對(duì)氣候系統(tǒng)內(nèi)部自然變率的影響作用。為了獲得更為可信的20世紀(jì)以來(lái)WACE的年代際內(nèi)部變率及其影響因素，我們參照Dai et al.（2015）的方法，選用參與第六次國(guó)際耦合模式比較計(jì)劃（Coupled Model Intercomparison Project Phase 6，簡(jiǎn)稱CMIP6）的15個(gè)模式：AWI-CM-1-1-MR、BCC-ESM1、CAMS-CSM1-0、CanESM5、CASESM2-0、CESM2、CESM-WACCM、CMCC-CM2-SR5、E3SM1-1、EC-Earth3-Veg-LR、FGOALS-f3-L、GFDL-ESM4、IPSL-CM6A-LR、MPI-ESM1-2-HR、MRI-ESM2-0。對(duì)于地表氣溫，將1910～2014年Historical試驗(yàn)結(jié)果與2015～2020年的SSP2-4.5排放情景下的試驗(yàn)結(jié)果進(jìn)行拼接；對(duì)于海溫，選用1910～2014年Historical試驗(yàn)結(jié)果。模式數(shù)據(jù)集合平均反映外強(qiáng)迫對(duì)溫度變量的影響，在外強(qiáng)迫影響下，20世紀(jì)60年代前地表氣溫及海溫變暖趨勢(shì)較平緩，60年代后變暖趨勢(shì)較明顯（圖1）。相較于傳統(tǒng)去除線性趨勢(shì)方法，模式對(duì)全球變暖趨勢(shì)擬合較合理，這與Dai et al.（2015）研究結(jié)果相似。因此，我們利用CMIP6的15個(gè)模式集合平均后全球平均冬季地表氣溫與海溫序列，分別對(duì)觀測(cè)得到的冬季地表氣溫與海溫去除外部強(qiáng)迫造成的趨勢(shì)。

2.2 方法

我們參照Daiet al.（2015）的方法，對(duì)地表氣溫及海溫去除趨勢(shì)，以盡可能消除外強(qiáng)迫影響：將1910/1911～2019/2020年冬季平均地表氣溫或海溫觀測(cè)場(chǎng)視為變量T(n,i)，其中n為年份，i為格點(diǎn)。將T(n,i)視為外部強(qiáng)迫TF(n,i)加上非外部強(qiáng)迫TI(n,i)的結(jié)果，即：

其中，TF(n,i)包含了大量的外強(qiáng)迫信息，TI(n,i)包含了少量外強(qiáng)迫信息。我們的目的是盡可能去除T(n,i)中包含的外強(qiáng)迫信息，因此，我們將TF(n,i)視作人類活動(dòng)等外強(qiáng)迫造成的溫度變化，將TI(n,i)視作由內(nèi)部變率造成的溫度變化，TF(n,i)可表示為

其中，Tm(n)表示15個(gè)模式集合平均后全球平均地表氣溫或海溫序列（圖1）。bm(i)可以通過(guò)回歸方程T(n,i)=bm(i)Tm(n)來(lái) 估計(jì)。bm(i)Tm(n)已包含大部分外強(qiáng)迫信息，對(duì)于年代際變率可忽略局地變化項(xiàng)εF(n,i)（王岱等,2020）。因此，溫度內(nèi)部變率TI(n,i)為

圖1 CMIP6模式Historical試驗(yàn)集合平均（藍(lán)色實(shí)線）與觀測(cè)（黑色實(shí)線）的（a）1910/1911～2019/2020年冬季全球平均地表氣溫距平（單位：°C）與（b）1910/1911～2013/2014年冬季全球平均海溫距平（單位：°C）的3年滑動(dòng)平均。距平參考時(shí)段為1960/1961～1989/1990年冬季Fig.1 Time series of CMIP6 models’ ensemble mean(blue solid lines)and the observed anomalies(black solid lines)of(a)three-year moving average of global mean surface air temperature(SAT,units:°C)during boreal winter of 1910/1911–2019/2020,and(b)three-year moving average of sea surface temperature(SST, units:°C)during boreal winter of 1910/1911–2013/2014.The anomalies reference period is 1960/1961–1989/1990 winter

其余變量由于全球平均變暖趨勢(shì)不顯著（Dong and Dai,2015），采用線性回歸方法去除線性趨勢(shì)。

根據(jù)過(guò)去的研究（Moriet al.,2014;Sung et al.,2018），我們對(duì)去除趨勢(shì)并緯度加權(quán)后的歐亞中高緯地區(qū)冬季地表氣溫求取距平，然后進(jìn)行經(jīng)驗(yàn)正交函數(shù)分解，并將第二模態(tài)相應(yīng)的主成分（PC2）定義為WACE指數(shù)。由于GISS提供的地表氣溫資料在80°N以北缺測(cè)，此套數(shù)據(jù)計(jì)算出的WACE指數(shù)相較傳統(tǒng)WACE指數(shù)的計(jì)算范圍偏小。因此，我們分別采用GISS與20CR兩套數(shù)據(jù)的地表氣溫計(jì)算WACE指數(shù)，比較兩套資料描述WACE模態(tài)的能力。其中，GISS地表氣溫進(jìn)行EOF分解的范圍是（20°～80°N，0°～150°E）；20CR地表氣溫進(jìn)行EOF分解的范圍是（20°～90°N，0°～150°E）。為判斷EOF分析所得模態(tài)是否相互獨(dú)立，我們使用North方法（North et al.,1982）進(jìn)行顯著性檢驗(yàn)。

對(duì)WACE指數(shù)進(jìn)行滑動(dòng)平均、滑動(dòng)t檢驗(yàn)、墨西哥小帽小波分析獲得WACE年代際變化特征。還使用差值分析以及21年滑動(dòng)平均以后的相關(guān)分析獲得與WACE的年代際變化相聯(lián)系的大氣環(huán)流以及北大西洋海溫異常。對(duì)于21年滑動(dòng)平均以后的序列A與序列B的相關(guān)分析，計(jì)算有效自由度n′：

其中，n為樣本數(shù)，ai為序列A自相關(guān)系數(shù)，bi為序列B自相關(guān)系數(shù)，i為滯后時(shí)長(zhǎng)（Quenouille,1952）。根據(jù)有效自由度來(lái)計(jì)算相關(guān)系數(shù)的顯著性檢驗(yàn)的閾值。差值分析與相關(guān)分析均使用Studentt檢驗(yàn)方法進(jìn)行顯著性檢驗(yàn)。

為分析WACE年代際變化的主要影響因素，根據(jù)溫度傾向方程：

為了研究阻塞發(fā)生頻率，參照Tibaldi and Molteni（1990）及Davini et al.（2012）定義阻塞指數(shù)：

其中，經(jīng)度λ取值范圍為60°W～100°E、緯度φ取值范 圍為30°N～75°N，φs=φ?15°，φn=φ+15°，Z為位勢(shì)高度。Gn、Gs分別代表(λ,φ)與北側(cè)、南側(cè)的位勢(shì)高度差梯度，當(dāng)Gn0 gpm deg?1（deg表示緯距或經(jīng)距）時(shí)，判斷當(dāng)日(λ,φ)發(fā)生阻塞，冬季（12月1日至次年2月28日）發(fā)生阻塞的次數(shù)，占整個(gè)冬季的比例為冬季阻塞發(fā)生頻率。

2.3 模式試驗(yàn)

本文利用美國(guó)國(guó)家大氣研究中心（National Center for Atmosphere Research，簡(jiǎn)稱NCAR）研制的NCAR CAM3.0模式探究北大西洋海溫異常對(duì)大氣環(huán)流及WACE年代際變化的影響。該模式為全球大氣環(huán)流模式，水平分辨率為T42，即緯向均勻分布128個(gè)格點(diǎn)，經(jīng)向分布64個(gè)高斯格點(diǎn)，垂直方向?yàn)棣舀Cp混合坐標(biāo)，共26層，模式層頂為2.917 hPa（Collins et al.,2004）。我們?cè)O(shè)計(jì)了控制試驗(yàn)及敏感性試驗(yàn)，分析冬季北大西洋關(guān)鍵區(qū)[區(qū)域1（50°～60°N，60°～20°W）與區(qū)域2（25°～32°N，55°～45°W）]海溫正異常影響下的大氣環(huán)流及地表氣溫變化。

控制試驗(yàn)：采用模式提供的多年平均氣候態(tài)海溫場(chǎng)作為外強(qiáng)迫，積分50年，將模式輸出的氣象要素的后30年平均作為控制試驗(yàn)的結(jié)果。

敏感性試驗(yàn)：在模式提供的多年平均氣候態(tài)海溫的基礎(chǔ)上，將冬季北大西洋關(guān)鍵區(qū)海溫疊加正異常，其他區(qū)域仍使用氣候態(tài)海溫。將此海溫異常場(chǎng)作為外強(qiáng)迫，積分50年，將模式輸出的氣象要素的后30年平均作為敏感性試驗(yàn)的結(jié)果。

敏感性試驗(yàn)與控制試驗(yàn)結(jié)果的差值場(chǎng)可以反映冬季北大西洋關(guān)鍵區(qū)海溫正異常影響下的大氣環(huán)流及地表氣溫的異常。

3 WACE的年代際變化和相應(yīng)的大氣環(huán)流異常

3.1 WACE的年代際變化特征

觀測(cè)資料顯示，1910/1911～2019/2020年冬季全球平均地表氣溫具有變暖趨勢(shì)（圖1a）。CMIP6模式集合平均序列反映外強(qiáng)迫的影響（Dai et al.,2015）。因此，計(jì)算WACE模態(tài)及指數(shù)前，利用模式集合平均序列對(duì)地表氣溫去除趨勢(shì)，即去除外強(qiáng)迫對(duì)地表氣溫內(nèi)部自然變率的影響。

首先，比較GISS與20CR兩套資料對(duì)WACE模態(tài)的描述能力。對(duì)來(lái)源于GISS觀測(cè)資料的1910/1911～2019/2020年歐亞大陸中高緯地區(qū)（20°～80°N，0°～150°E）冬季地表氣溫距平進(jìn)行EOF分析，第一模態(tài)為歐洲大陸溫度一致變化型，解釋方差為39.0%（圖2a）；第二模態(tài)為歐亞大陸和北極溫度反向變化型，解釋方差為20.6%（圖2c）。這兩個(gè)EOF模態(tài)通過(guò)了North檢驗(yàn)，說(shuō)明它們是相互獨(dú)立和正交的（North et al.,1982）。利用20CR再分析資料，將1910/1911～2013/2014年歐亞中高緯地區(qū)（20°～90°N，0°～150°E）冬季地表氣溫距平進(jìn)行EOF，得到與GISS資料相似模態(tài)（圖2e、g），第一模態(tài)解釋方差為28.7%；第二模態(tài)解釋方差為13.6%，兩個(gè)模態(tài)同樣通過(guò)North檢驗(yàn)。GISS資料EOF第一模態(tài)主成分（GISS PC1）（圖2b）與20CR資料EOF第一模態(tài)主成分（20CR PC1）（圖2f）相關(guān)系數(shù)為0.91，GISS PC2（圖2d）與20CR PC2（圖2h）相關(guān)系數(shù)為0.82，均通過(guò)0.01顯著性水平的顯著性t檢驗(yàn)?？偟膩?lái)說(shuō)，兩套資料的結(jié)果與以前的研究結(jié)果相符合（Moriet al.,2014;Sung et al.,2018），但是GISS資料的第一模態(tài)和第二模態(tài)所占的解釋方差更高。由于GISS為觀測(cè)資料，可信度更高，我們之后的分析將使用根據(jù)GISS資料計(jì)算的WACE指數(shù)。

圖2 （a–d）1910/1911～2019/2020年GISS觀測(cè)資料、（e–h）1910/1911～2013/2014年20CR再分析資料的冬季地表氣溫距平場(chǎng)經(jīng)驗(yàn)正交分解（EOF）的（a、e）第一模態(tài)及（b、f）主成分（PC1）、（c、g）第二模態(tài)及（d、h）主成分（PC2）Fig.2(a,e)The first mode and(b,f)time coefficient(PC1),and(c,g)the second mode and(d, h) time coefficient(PC2)of EOF for surface air temperature anomalies in boreal winter obtained from(a–d)GISS(Goddard Institute for Space Studies)observation data for period of 1910/1911–2019/2020,(e–h)20CR (Twentieth Century Reanalysis)reanalysisdata for the period of 1910/1911–2013/2014

經(jīng)小波分析發(fā)現(xiàn)，WACE指數(shù)具有20年、70年左右的顯著振蕩周期，20年振蕩周期在1935、1965、1975、2000年左右顯著（圖3）。通過(guò)11年滑動(dòng)平均與11年滑動(dòng)t檢驗(yàn)，得到WACE的年代際躍變點(diǎn)為：1929年、1940年、1958年、1967年、1978年、2005年（圖4）。其中，1940年雖為顯著躍變點(diǎn)，但其11年滑動(dòng)平均序列在1929～1957年處于正值，并且1958年這個(gè)躍變點(diǎn)仍然處于1940年躍變的下方。因此，我們將WACE的年代際位相劃分為：P1（1929～1957年），P2（1958～1966年），P3（1967～1977年），P4（1978～2004年），P5（2005年以后）5個(gè)時(shí)期，其中P1、P3、P5為正位相，P2、P4為負(fù)位相，這里的1929年指1928/1929年冬季，其余年份以此類推。

圖3 1910/1911～2019/2020年GISS資料計(jì)算的（a）冬季W(wǎng)ACE指數(shù)墨西哥帽小波變換系數(shù)（藍(lán)色實(shí)線）和（b）小波全譜（藍(lán)色實(shí)線）。水平黑色點(diǎn)劃線表示顯著振蕩周期。圖a中黑色粗實(shí)線表示小波變換系數(shù)通過(guò)0.1顯著性水平的白噪音檢驗(yàn)，黑色細(xì)實(shí)線兩側(cè)區(qū)域表示邊界效應(yīng)影響域。圖b中橙色虛線表示白噪音檢驗(yàn)的0.1顯著性水平線Fig.3(a) Wavelet analysis of Mexican cap WACE(Warm Arctic–Cold Eurasia)index(blue solid lines)and(b) the global wavelet power spectrum(blue solid line)calculated from GISSdatasets during boreal winter for the period of 1910/1911–2019/2020.The black horizontal dotted-dashed lines indicate the significant oscillation periods.In Fig.a,the black thick lines indicate passing the white noise test at 0.1 significance level, the regions on both sides of the black thin solid lines represent the area with boundary effect.In Fig. b,the orange dotted line indicates the 0.1 significance level for white noise test

圖4 1910/1911～2019/2020年冬季（a）WACE指數(shù)（黑色曲線）及11年滑動(dòng)平均（橙色曲線），（b）WACE指數(shù)的11年滑動(dòng)t檢驗(yàn)（藍(lán)色曲線），黑色水平實(shí)線為0.1顯著性水平線Fig.4 (a)The WACEindex(black line)and 11-year moving average(orange line)and (b)11-yr moving t-test of WACEindex during boreal winter from 1910/1911 to 2019/2020.In Fig. b,the black solid lineis the significancelevel at 0.1

3.2 WACE不同年代際位相時(shí)期大氣環(huán)流異常

為了獲得導(dǎo)致WACE年代際變化的異常大氣環(huán)流型，我們對(duì)WACE正、負(fù)位相的冬季大氣環(huán)流場(chǎng)進(jìn)行差值分析。結(jié)果表明，WACE正位相時(shí)，500 hPa及200 hPa歐亞中高緯地區(qū)存在顯著的位勢(shì)高度正異常，也即歐亞地區(qū)極渦強(qiáng)度在WACE正位相時(shí)期異常減弱。200 hPa風(fēng)場(chǎng)中歐亞中高緯地區(qū)存在異常反氣旋（圖5a）；對(duì)于低層，海平面氣壓場(chǎng)上表現(xiàn)為歐洲大陸高壓異常增強(qiáng)，伴隨850 hPa風(fēng)場(chǎng)在歐亞中高緯地區(qū)存在異常反氣旋（圖5b）。因此，WACE正位相時(shí)期，大氣環(huán)流最顯著的特征是歐亞中高緯地區(qū)由低層至高層的異常高壓，并且異常高壓中心隨高度向北極傾斜。

圖5 冬季W(wǎng)ACE正、負(fù)位相的差值：（a）200 hPa位勢(shì)高度場(chǎng)（陰影，單位：gpm）、200 hPa風(fēng)場(chǎng)（矢量箭頭，單位：m s?1）和500 hPa位勢(shì)高度場(chǎng)（等值線，單位：gpm）；（b）海平面氣壓（陰影，單位：Pa）、850 hPa風(fēng)場(chǎng)（矢量箭頭，單位：m s?1）。圖a、b中的打點(diǎn)區(qū)域表示200 hPa位勢(shì)高度場(chǎng)、海平面氣壓通過(guò)0.1顯著性水平的t檢驗(yàn)；紫色等值線表示500 hPa位勢(shì)高度通過(guò)0.1顯著性水平的t檢驗(yàn)；風(fēng)場(chǎng)只繪制出緯向風(fēng)通過(guò)0.1顯著性水平的t檢驗(yàn)部分Fig.5 Differences between positive and negative phases of WACE in winter:(a)200-hPa geopotential height field(shadings,units:gpm),200-hPa wind field(vector arrows,units:m s?1),and 500-hPa geopotential height field(contours, units:gpm);(b)sea level pressure(shadings,units:Pa)and 850-hPa wind field (vector arrows,units: m s?1).In Figs.a,b,the black dotsareasindicate 200-hPa geopotential height field,sea level pressure passed the t-test at 0.1 significance level,and the purple contours areas indicate 500-hPa geopotential height field passed the t-test at 0.1 significance level,only thezonal wind field passed the t-test at 0.1 significancelevel aredrew

為分析影響WACE年代際變化的主要因素，首先診斷了冬季950 hPa溫度平流在WACE正、負(fù)位相的差值，發(fā)現(xiàn)在WACE年代際正位相時(shí)，北極存在顯著異常暖平流，歐亞大陸存在顯著異常冷平流（圖6a），因此，溫度平流可能是導(dǎo)致WACE發(fā)生年代際變化的原因之一。另外，診斷950 hPa垂直運(yùn)動(dòng)導(dǎo)致的溫度變化在WACE正、負(fù)位相的差值，發(fā)現(xiàn)WACE年代際正位相時(shí)，西伯利亞地區(qū)存在由垂直運(yùn)動(dòng)導(dǎo)致的溫度正負(fù)異常交替分布的特征（圖6b），但這與WACE模態(tài)不一致，說(shuō)明WACE的年代際變化可能與垂直運(yùn)動(dòng)導(dǎo)致的加熱關(guān)聯(lián)較小。

進(jìn)一步對(duì)WACE正、負(fù)位相時(shí)期地表輻射通量及感熱和潛熱通量進(jìn)行差值分析，取向上為正，結(jié)果表明WACE年代際正位相時(shí)，北極大部分地區(qū)的地表感熱通量為負(fù)異常，新地島以東的小部分極區(qū)地表感熱通量為正異常但未通過(guò)顯著性檢驗(yàn)，歐亞大陸地表感熱通量為正異常（圖6c）。這說(shuō)明在北極地區(qū)，大氣向地表傳輸感熱，而在歐亞大陸，地表向大氣傳輸感熱，即北極大氣并未得到地表傳輸?shù)母袩?，歐亞大陸的大氣得到地表傳輸?shù)母袩?。這意味著感熱不會(huì)造成北極地表氣溫增加以及歐亞大陸地表氣溫降低，感熱可能不是導(dǎo)致WACE發(fā)生年代際變化的主要原因。另外，在WACE年代際正位相時(shí)，新地島以東的地表潛熱通量為顯著正異常，新地島以西為負(fù)異常，歐亞大陸的地表潛熱通量為顯著負(fù)異常（圖6d），說(shuō)明潛熱可能是影響WACE年代際變化的原因之一。因此，WACE的年代際變化可能受到潛熱的影響，而感熱可能不是導(dǎo)致其發(fā)生年代際變化的主要因素。

地表凈長(zhǎng)波輻射通量的差值分析表明，在WACE年代際正位相時(shí)，北極存在地表凈長(zhǎng)波輻射顯著負(fù)異常，而歐亞大陸存在顯著正異常（圖6e），這說(shuō)明可能極區(qū)大氣吸收并向下放射長(zhǎng)波輻射，對(duì)北極大氣起到保溫作用，而歐亞大陸反之。另外，在WACE年代際正位相時(shí)，地表凈短波輻射在北極地區(qū)無(wú)顯著異常，在歐亞大陸存在負(fù)異常（圖6f），這可能是歐亞大陸地表氣溫負(fù)異常的原因之一。

通過(guò)上面的分析，我們發(fā)現(xiàn)溫度平流、長(zhǎng)波輻射和潛熱造成的非絕熱加熱可能是導(dǎo)致WACE年代際變化的主要影響因素，短波輻射可能主要影響歐亞大陸氣溫，而對(duì)北極地表氣溫異常的影響較小，感熱可能不是使WACE發(fā)生年代際變化的主要因素。長(zhǎng)波輻射、潛熱與水汽的關(guān)聯(lián)十分密切，水汽吸收長(zhǎng)波輻射，使大氣對(duì)于長(zhǎng)波輻射的透過(guò)率降低，可能導(dǎo)致向下長(zhǎng)波輻射的異常，對(duì)大氣起到保溫作用（Park et al.,2015;Woods and Caballero,2016），同時(shí)，水汽的異常還可能影響降水的異常，降水進(jìn)一步影響云的形成以及潛熱釋放（秦大河,2018；鄭彬等,2019）。

通過(guò)診斷整層水汽通量及散度在WACE正、負(fù)位相的差值發(fā)現(xiàn)，WACE正位相時(shí)，歐亞大陸有向西的異常水汽輸送，歐亞大陸西北部的水汽向極區(qū)輸送，極區(qū)存在顯著的水汽輻合。同時(shí)，大西洋向歐亞大陸平流的水汽減少，歐亞中緯度地區(qū)出現(xiàn)顯著的水汽輻散（圖6g）。另外，WACE正位相時(shí)，極區(qū)降水正異常，有利于潛熱釋放；歐亞大陸降水負(fù)異常，不利于潛熱釋放（圖6h）。綜上，WACE正位相時(shí)，北極地區(qū)水汽輻合且降水正異常，可能造成向下長(zhǎng)波輻射增加以及潛熱的釋放導(dǎo)致極區(qū)增溫；與此同時(shí)，歐亞大陸水汽輻散并出現(xiàn)降水負(fù)異常，可能導(dǎo)致向下長(zhǎng)波輻射減少且不利于潛熱釋放，進(jìn)而歐亞大陸溫度偏低。

圖6 冬季W(wǎng)ACE正、負(fù)位相時(shí)期的差值：（a）950 hPa溫度平流（單位：10?5 K s?1）、（b）950 hPa垂直運(yùn)動(dòng)導(dǎo)致的溫度變化（單位：10?5 K s?1）、（c）地表感熱通量（單位：W s?1）、（d）地表潛熱通量（單位：W s?1）、（e）地表凈長(zhǎng)波輻射通量（單位：W s?1）、（f）地表凈短波輻射通量（單位：W s?1）、（g）整層（1000～300 hPa）水汽通量（箭頭，單位：kg m?1 s?1）及散度（陰影，單位：10?6 kg m?2 s?1）、（h）地表降水率（單位：mm s?1）。打點(diǎn)區(qū)域通過(guò)0.1顯著性水平的t檢驗(yàn)Fig.6 Differences between positive and negative phases of WACE in winter:(a)950-hPa temperature advection(units:10?5 K s?1);(b)950-hPa temperaturechange(units:10?5 K s?1)dueto vertical motion;(c)surface sensibleheat fluxes(units:W s?1);(d)surfacelatent heat fluxes(units:W s?1);(e)net surface long-wave radiation fluxes(units:W s?1);(f)net surfaceshort-wave radiation fluxes(units:W s?1);(g)water vapor fluxes(arrows,units:kg m?1 s?1)in thewhole layer and their divergences(shadings,units:10?6 kg m?2 s?1);(h)surface precipitation rate(units: mm s?1).The black dots area passed the t-test at 0.1 significancelevel

為了探討阻塞發(fā)生的頻率和WACE年代際變化之間的關(guān)系，本文將WACE正位相與負(fù)位相時(shí)期冬季阻塞發(fā)生頻率進(jìn)行差值（圖7）。在WACE正位相時(shí)，40°～80°E的中高緯地區(qū)冬季阻塞發(fā)生頻率為正異常，顯著的正異常中心位于60°E附近，即烏拉爾阻塞發(fā)生頻率最高的經(jīng)度（Diao et al.,2006）。烏拉爾阻塞發(fā)生頻率偏高，一方面意味著阻塞的維持使暖空氣向極區(qū)輸送而導(dǎo)致極區(qū)增暖，另一方面表示冷空氣侵襲歐亞大陸更加頻繁，可能導(dǎo)致了歐亞大陸降溫（Luo et al.,2016）。

圖7 冬季W(wǎng)ACE正、負(fù)位相時(shí)期阻塞發(fā)生頻率（單位：d?1）差值，陰影區(qū)域通過(guò)0.1顯著性水平的t檢驗(yàn)Fig.7 Differences between positiveand negative phasesof WACE in blocking frequency (units:d?1)during boreal winter, the shaded area passed the t-test at 0.1 significancelevel

4 北大西洋海溫與WACE年代際變化的聯(lián)系

根據(jù)前面的分析，WACE正位相時(shí)期，水汽向極區(qū)輸送，影響向下長(zhǎng)波輻射與潛熱釋放，導(dǎo)致極區(qū)增暖，歐亞大陸反之。北大西洋位于歐亞大陸上游地區(qū)，向大氣提供水汽，且北大西洋海溫具有顯著的多年代際變化特征，其年代際變化影響著北半球大部分地區(qū)的氣候（秦大河,2018）。這意味著北大西洋海溫有可能通過(guò)中高緯的海氣相互作用在WACE的年代際變化中起作用。1910/1911～2013/2014年冬季全球平均海溫與地表氣溫序列具有相似變暖趨勢(shì)，亦即20世紀(jì)60年代前海溫變暖趨勢(shì)較平緩，60年代后變暖趨勢(shì)較陡峭（圖1b）。因此，分析年代際尺度上WACE與北大西洋海溫的聯(lián)系前，我們利用CMIP6模式數(shù)據(jù)對(duì)北大西洋海溫去除趨勢(shì)，即去除外強(qiáng)迫對(duì)海溫的影響，著重研究海溫的內(nèi)部自然變率對(duì)WACE年代際變化的可能影響作用。

為了獲得影響WACE年代際變化的海溫關(guān)鍵區(qū)，我們將1910/1911～2013/2014年去除外強(qiáng)迫并進(jìn)行21年滑動(dòng)平均后的冬季W(wǎng)ACE指數(shù)及北大西洋海溫進(jìn)行相關(guān)分析，結(jié)果表明北大西洋海溫存在兩個(gè)顯著正相關(guān)區(qū)，分別位于北大西洋中緯度和副熱帶地區(qū)（圖8a）。根據(jù)海溫的兩個(gè)顯著正相關(guān)區(qū)，定義北大西洋區(qū)域1（50°～60°N，60°～20°W）與區(qū)域2（25°～32°N，55°～45°W）的平均冬季海溫為北大西洋海溫指數(shù)（NAS）（圖8b）。21年滑動(dòng)平均后的WACE與NAS指數(shù)相關(guān)系數(shù)為0.75，通過(guò)0.1顯著性水平的顯著性t檢驗(yàn)，表明兩者在年代際時(shí)間尺度上有著密切的聯(lián)系。另外，將21年滑動(dòng)平均后的NAS指數(shù)與大氣環(huán)流場(chǎng)進(jìn)行相關(guān)分析，發(fā)現(xiàn)結(jié)果與影響WACE年代際變化的大氣環(huán)流特征相似（圖略）。

圖8 1910/1911～2013/2014年（a）21年滑動(dòng)平均的WACE指數(shù)與冬季海溫相關(guān)系數(shù)，（b）21年滑動(dòng)平均的WACE指數(shù)（橙色實(shí)線）、北大西洋區(qū)域1（圖a上方藍(lán)色矩形區(qū)域，50°～60°N，60°～20°W）與區(qū)域2（圖a下方藍(lán)色矩形區(qū)域，25°～32°N，55°～45°W）平均的去趨勢(shì)后的北大西洋海溫指數(shù)（NAS，灰色實(shí)線）、21年滑動(dòng)平均的NAS指數(shù)（紅色實(shí)線）、11年滑動(dòng)平均的NAS指數(shù)（藍(lán)色實(shí)線）。圖a中，打點(diǎn)區(qū)域通過(guò)0.1顯著性水平的t檢驗(yàn)Fig.8(a)Correlation coefficients between WACEindex and winter sea surface temperature after the 21-yr moving average,(b)WACEindex(solid orange line)after the 21-yr moving average,sea surface temperature index(solid gray line)of the North Atlantic(NAS)averaged in region 1(50°–60°N,60°–20°W)and region 2(25°–32°N,55°–45°W)after the trend was removed,NAS(solid red line)after the 21-yr moving average, NAS(solid red line)after the11-yr moving average from 1910/1911 to 2013/2014.In Fig.a,the dotted area passed the t-test at 0.1 significance level

為了驗(yàn)證冬季北大西洋關(guān)鍵區(qū)的海溫異常通過(guò)影響大氣環(huán)流，進(jìn)而導(dǎo)致WACE的年代際變化，利用NCAR CAM3.0模式設(shè)計(jì)了控制試驗(yàn)與敏感性試驗(yàn)。圖9a表示敏感性試驗(yàn)與控制試驗(yàn)的海溫外強(qiáng)迫差值。由于NCAR CAM3.0模式是大氣環(huán)流模式，因此敏感性試驗(yàn)與控制試驗(yàn)結(jié)果的差值表示冬季北大西洋關(guān)鍵區(qū)海溫的正異常對(duì)大氣環(huán)流的影響。

冬季地表氣溫的模擬結(jié)果差值場(chǎng)顯示，歐亞的北極地區(qū)出現(xiàn)溫度正異常，貝加爾湖西側(cè)大陸出現(xiàn)溫度的負(fù)異常，呈現(xiàn)北極暖而歐亞冷的特征（圖9b），溫度異常的中心位置較WACE整體偏南（圖2c）。這說(shuō)明在冬季北大西洋關(guān)鍵區(qū)海溫正異常的影響下，下游的中高緯地表氣溫出現(xiàn)WACE模態(tài)，數(shù)值模擬的結(jié)果較好。

圖9 （a）疊加在氣候態(tài)（1949～2001）海溫上的冬季海溫異常強(qiáng)迫（單位：K）及（b）模式模擬的冬季地表氣溫的差值（敏感性試驗(yàn)結(jié)果減控制試驗(yàn)結(jié)果，單位：K）Fig.9(a)Winter SST anomalies(units:K)forcing superimposed on SST of climatic states(1949–2001),(b)the differences(sensitivity experiment results minus control experiment results)in simulated surfaceair temperature(units:K)

前面已經(jīng)獲得導(dǎo)致WACE年代際變化的大氣環(huán)流異常型，為了驗(yàn)證冬季北大西洋關(guān)鍵區(qū)海溫的年代際變化是否導(dǎo)致了這種大氣環(huán)流異常，將200 hPa風(fēng)場(chǎng)和地表氣壓進(jìn)行敏感性試驗(yàn)與控制試驗(yàn)的差值，發(fā)現(xiàn)關(guān)鍵區(qū)海溫正異常導(dǎo)致北大西洋東側(cè)的西歐大陸200 hPa產(chǎn)生異常氣旋性環(huán)流（圖10a）、地表氣壓負(fù)異常（圖10b），歐亞大陸中高緯地區(qū)出現(xiàn)異常反氣旋性環(huán)流，貝加爾湖西側(cè)出現(xiàn)東風(fēng)異常（圖10a），地表歐亞大陸高壓增強(qiáng)（圖10b），導(dǎo)致冷空氣容易侵襲歐亞大陸導(dǎo)致大陸降溫。這與WACE正位相時(shí)期大氣環(huán)流場(chǎng)的特征相似。

圖10 模式模擬的冬季（a）200 hPa水平風(fēng)場(chǎng)（單位：m s?1）、（b）地表氣壓（單位：Pa）差值（敏感性試驗(yàn)結(jié)果減控制試驗(yàn)結(jié)果）Fig.10 Differences(sensitivity experiment results minus control experiment results)in simulated(a)200-hPa horizontal wind(units:m s?1),(b) surface pressure (units:Pa)in winter

在前面的統(tǒng)計(jì)分析中，發(fā)現(xiàn)水汽可能是影響WACE年代際變化的重要原因。因此，利用數(shù)值試驗(yàn)來(lái)驗(yàn)證冬季北大西洋關(guān)鍵區(qū)海溫正異常是否能影響水汽輸送。將敏感性試驗(yàn)與控制試驗(yàn)的整層水汽通量及散度進(jìn)行差值，發(fā)現(xiàn)關(guān)鍵區(qū)海溫正異常導(dǎo)致非洲北側(cè)的西歐大陸出現(xiàn)水汽的氣旋式輸送，而歐亞大陸出現(xiàn)反氣旋式水汽輸送，二者相互配合，導(dǎo)致水汽向極區(qū)輸送造成極區(qū)水汽輻合，而阻擋水汽向歐亞大陸平流進(jìn)而導(dǎo)致歐亞大陸水汽輻散（圖11），這與WACE正位相時(shí)期水汽輸送特征相似。這表明了冬季北大西洋關(guān)鍵區(qū)海溫正異常通過(guò)影響水汽輸送，進(jìn)而導(dǎo)致了極區(qū)增溫，歐亞大陸降溫。綜上所述，統(tǒng)計(jì)分析與數(shù)值試驗(yàn)結(jié)果充分說(shuō)明了，WACE的年代際變化與北大西洋關(guān)鍵區(qū)海溫存在密切聯(lián)系。

圖11 模式模擬的冬季整層（1000～300 hPa）水汽通量（箭頭，單位：kg m?1 s?1）及水汽通量散度差值（敏感性試驗(yàn)結(jié)果減控制試驗(yàn)結(jié)果，陰影，單位：10?6 kg m?2 s?1）Fig.11 Differences(sensitivity experiment results minuscontrol experiment results)in simulated water vapor fluxes(arrows, units:kg m?1 s?1)in the wholelayer (1000–300 hPa)and their divergence (shadings,units:10?6 kg m?2 s?1)

5 結(jié)果與討論

本文利用CMIP6模式數(shù)據(jù)集合平均得到的全球平均地表氣溫與海溫序列，去除人類活動(dòng)等外部強(qiáng)迫影響作用后，研究年代際時(shí)間尺度上WACE模態(tài)的內(nèi)部變率及其成因。在獲得WACE不同年代際位相時(shí)期大氣環(huán)流場(chǎng)異常型的基礎(chǔ)上，進(jìn)一步揭示了北大西洋海溫異常影響WACE年代際變化的物理機(jī)制并用數(shù)值模式加以驗(yàn)證。主要得到以下結(jié)論：

（1）GISS資料與20CR資料計(jì)算的冬季W(wǎng)ACE模態(tài)與指數(shù)大致相同，兩套資料得到的WACE指數(shù)相關(guān)系數(shù)為0.82。1910/1911～2019/2020年WACE指數(shù)具有顯著的年代際變化特征，

其顯著振蕩周期為20年、70年左右，其在1929年、1940年、1958年、1967年、1978年、2005年發(fā)生年代際躍變。因此，結(jié)合11年滑動(dòng)平均以及滑動(dòng)t檢驗(yàn)結(jié)果，在年代際尺度上將WACE劃分為五個(gè)位相：P1（1929～1957年），P2（1958～1966年），P3（1967～1977年），P4（1978～2004年），P5（2005年以后），其中P1、P3、P5為正位相，P2、P4為負(fù)位相。

（2）WACE正位相時(shí)期，歐亞中高緯地區(qū)存在一個(gè)異常的高壓，高壓中心隨高度向北極傾斜，

呈現(xiàn)正壓結(jié)構(gòu)。具體而言，在對(duì)流層低層，歐亞大陸海平面氣壓異常增強(qiáng)，伴隨著850 hPa的異常反氣旋；在對(duì)流層中高層，歐亞極渦異常減弱，相應(yīng)的200 hPa風(fēng)場(chǎng)存在異常反氣旋。另外，在WACE正位相時(shí)期，40°～80°E阻塞發(fā)生頻率偏高，有利于熱量向極區(qū)輸送，極區(qū)存在異常暖平流，且歐亞大陸中高緯反氣旋環(huán)流向極區(qū)輸送水汽，極區(qū)水汽輻合，向下長(zhǎng)波輻射增加，對(duì)流活動(dòng)增加使得潛熱釋放產(chǎn)生非絕熱加熱，有利于極區(qū)增暖。同時(shí)，阻塞頻繁發(fā)生使得極區(qū)冷空氣更易入侵歐亞大陸，歐亞大陸存在異常冷平流，且水汽由大陸西北側(cè)向極區(qū)輸送，導(dǎo)致歐亞大陸水汽輻散，進(jìn)而向下長(zhǎng)波輻射減小，對(duì)流活動(dòng)減少不利于潛熱釋放，導(dǎo)致歐亞大陸變冷。

（3）數(shù)值模擬結(jié)果表明，冬季北大西洋中緯度及副熱帶地區(qū)海溫的年代際變化可以通過(guò)強(qiáng)迫低層與高層異常大氣環(huán)流型，進(jìn)而導(dǎo)致WACE的年代際變化。冬季北大西洋關(guān)鍵區(qū)海溫處于正異常時(shí)，海溫異常通過(guò)海氣相互作用強(qiáng)迫出下游西歐地區(qū)的氣旋性環(huán)流以及歐亞大陸反氣旋性環(huán)流，有利于水汽向極區(qū)輸送，使得極區(qū)增暖；此外，歐亞大陸50°～60°N出現(xiàn)異常東北風(fēng)，冷空氣更容易侵襲歐亞大陸，同時(shí)大西洋向東的水汽輸送受到阻礙，歐亞大陸水汽輻散，導(dǎo)致歐亞大陸溫度負(fù)異常。

研究表明，極渦在近幾十年向歐亞大陸偏移并減弱，對(duì)歐亞大陸溫度降低具有貢獻(xiàn)（Zhang et al.,2016）。作為北半球主要環(huán)流系統(tǒng)之一，極渦配合北極濤動(dòng)（Arctic Oscillation，簡(jiǎn)稱AO）、北大西洋濤動(dòng)（North Atlantic Oscillation，簡(jiǎn)稱NAO）及阻塞高壓等系統(tǒng)影響著歐亞地區(qū)氣候（張恒德等,2008；藍(lán)柳茹和李棟梁,2016）。根據(jù)本文研究，近幾十年WACE處于年代際正位相，正好與極渦的年代際減弱相對(duì)應(yīng)。

值得注意的是，AMO與格陵蘭海冰在周期為60～90年的時(shí)間尺度上存在協(xié)同變化，且AMO與PDO的位相共變可能是20世紀(jì)初北極變暖的原因（Miles et al.,2014;Tokinaga et al.,2017），這說(shuō)明AMO可能與WACE存在聯(lián)系。為了對(duì)比AMO以及NAS指數(shù)對(duì)WACE年代際變化的不同影響作用，我們計(jì)算了去除外強(qiáng)迫后21年滑動(dòng)平均后的冬季AMO指數(shù)與同期WACE指數(shù)的相關(guān)，相關(guān)系數(shù)為0.61，未通過(guò)0.1顯著性水平的顯著性t檢驗(yàn)。這可能由于AMO的周期為65～80年并且AMO指數(shù)的計(jì)算方式為北大西洋海溫的整體區(qū)域平均（李雙林等,2009）。這可能導(dǎo)致北大西洋海溫中20年左右較短周期的年代際變化信息被掩蓋。因此，本文選擇與WACE年代際尺度上相關(guān)性最顯著的區(qū)域，探究北大西洋海溫與WACE年代際變化的物理聯(lián)系。

通過(guò)本文的研究，我們對(duì)WACE年代際變化特征及物理機(jī)制有了更深刻的認(rèn)識(shí)，并從氣候系統(tǒng)內(nèi)部變率的角度，利用大氣環(huán)流模式探討了北大西洋海溫對(duì)WACE年代際變化的影響，未來(lái)需要深入研究北大西洋海溫影響WACE年代際變化的機(jī)制。另外，WACE年代際變化的成因仍有許多問(wèn)題值得討論，例如，洋流輸送是否影響了WACE的年代際變化？北大西洋海溫是否與太平洋地區(qū)海溫相配合，進(jìn)而影響WACE年代際變化？WACE的年代際變化對(duì)中國(guó)寒潮、暴雪等極端事件的發(fā)生有什么影響？我們將利用更豐富的資料，對(duì)WACE年代際變化成因及影響進(jìn)行探究。

致謝 真誠(chéng)地感謝陳軍明副研究員對(duì)數(shù)值模式試驗(yàn)的指導(dǎo)與幫助,使本研究工作得以順利開(kāi)展。同時(shí)真誠(chéng)感謝審稿專家對(duì)工作提供的寶貴建議與思路，使本研究工作得到很大改進(jìn)。