气象学报  2020, Vol. 78 Issue (2): 221-236   PDF    
http://dx.doi.org/10.11676/qxxb2020.020
中国气象学会主办。
0

文章信息

徐士琦, 刘刚, 廉毅, 曲绍君, 杨雪艳, 梁洪海. 2020.
XU Shiqi, LIU Gang, LIAN Yi, QU Shaojun, YANY Xueyan, LIANG Honghai. 2020.
近40年5—6月东北亚和北半球其他区域阻塞高压活动特征及其联系
Characteristics of blocking high in Northeast Asia and other areas of the northern Hemisphere and their relationships during May to June in recent 40 years
气象学报, 78(2): 221-236.
Acta Meteorologica Sinica, 78(2): 221-236.
http://dx.doi.org/10.11676/qxxb2020.020

文章历史

2019-08-13 收稿
2019-12-08 改回
近40年5—6月东北亚和北半球其他区域阻塞高压活动特征及其联系
徐士琦1,2 , 刘刚1,3 , 廉毅1,4 , 曲绍君5 , 杨雪艳2 , 梁洪海2     
1. 中高纬度环流系统与东亚季风研究开放实验室,长春,130062;
2. 吉林省气候中心,长春,130062;
3. 吉林省气象信息网络中心,长春,130062;
4. 吉林省气象科学研究所,长春,130062;
5. 吉林省气象探测保障中心,长春,130062
摘要: 基于美国环境预报中心/国家大气研究中心(NCEP/NCAR-Ⅰ)再分析资料,按照国际通行的阻塞高压(阻高)定义,采用计算机自动识别和天气图人工主、客观相结合的检测方法,获得了1979—2018年5月和6月逐日北半球中高纬度阻高活动的无缝隙分布;并利用线性相关和信息流因果关系等统计方法,研究了东北亚(鄂霍次克海)阻高与北半球其他阻高活动区特征及其联系。根据阻高活动集中程度与占北半球活动总天数的百分比,选定8个5月阻高主要活动区和10个6月阻高主要活动区。从统计结果来看,5月北太平洋、北美和乌拉尔山阻高活动占前3位,6月东北亚阻高活动天数跃升第1位。值得注意的是,初夏在贝加尔湖以南、青藏高原以北地区出现了明显的阻塞形势,称为青藏高原北阻高。5月和6月北半球各阻高活动中心集中的纬度分布,自西向东具有“正弦波”特征,其主要活动区引起的同期北半球500 hPa高度场的异常虽不尽相同,但均具有一定的遥相关特征;5月和6月东北亚阻高活动天数的年际变化分量具有明显的相对独立性。5月东北亚阻高是6月贝加尔湖阻高活动年际变化的稳定影响源,是6月青藏高原北阻高活动年际变化的不稳定影响源;唯有5月北太平洋阻高是6月东北亚阻高年际变化的稳定影响源。2019年6月北半球阻高活动异常的区域位于欧洲,达19 d,为气候平均值(7 d)的2.7倍,也是导致欧洲大陆初夏高温的主要影响系统。
关键词: 5—6月    北半球阻高    东北亚阻高    青藏高原北阻高    阻高活动特征    信息流    
Characteristics of blocking high in Northeast Asia and other areas of the northern Hemisphere and their relationships during May to June in recent 40 years
XU Shiqi1,2 , LIU Gang1,3 , LIAN Yi1,4 , QU Shaojun5 , YANY Xueyan2 , LIANG Honghai2     
1. Laboratory of Research for Middle-High Latitude Circulation System and East Asian Monsoon,Changchun 130062,China;
2. Climate Center of Jilin Province,Changchun 130062,China;
3. Jilin Meteorological Information Network Center,Changchun 130062,China;
4. Research Institute of Meteorological Science of Jilin Province,Changchun 130062,China;
5. Jilin Meteorological Observation Support Center,Changchun 130062,China
Abstract: Based on NCEP/NCAR-Ⅰ reanalysis and the internationally accepted definition of blocking high, the northern Hemisphere blocking high distribution maps were obtained from May to June over the 40-year period of 1979—2018 using subjective and objective detection methods. Linear correlation and information flow statistical methods are used to study interannual and interdecadal relationships between anomalous variations of blocking highs over Northeast Asia and other areas in the northern Hemisphere. According to the concentration of blocking high activities and their percentages of total days in the northern Hemisphere, 8 major areas in May and 10 major areas in June are selected. In May, the North Pacific blocking, the North America blocking, and the Ural blocking are the top three. In June, the number of days for the Northeast Asia blocking jumps to the top. It is worth noting that in the early summer, there is a significant increase in blocking to the south of the Baikal Lake and north of the Tibetan Plateau, which is called the North Tibetan Plateau blocking high. The latitude distribution of the northern Hemisphere blocking high activity centers in May and June shows a "sine wave" characteristic. Although the 500 hPa height field anomalies caused by blockings in their main active areas in the northern Hemisphere show different characteristics in May and June, they demonstrate teleconnection characteristics. The interannual variation of the Northeast Asia blocking high in May and June shows an obvious independence. The Northeast Asia blocking high in May is a source of the interannual variation of Lake Baikal blocking high in June; it is also an unstable source of the interannual variation of the North Tibetan Plateau blocking high in June. The North Pacific blocking in May is the only stable source of the interannual variation of the Northeast Asia blocking in June. In June 2019, the area of most frequent and abnormal blocking high activities in the northern Hemisphere remains in Europe for 19 d, which is 2.7 times the climatological average of 7 d. The blockings are also the main influence system leading to the early summer high temperature in Europe.
Key words: May to June    Northern Hemisphere blocking high    Northeast Asia blocking high    North Tibetan Plateau blocking high    Characteristics    Information flow    
1 引 言

影响中国天气、气候异常的因子除来自低纬度外强迫外,中高纬度大气环流尤其是阻塞高压(阻高)亦有着举足轻重的作用。由于阻高在大气中的持续时间较长且位置稳定少变,它的出现往往通过上下游效应给周边地区带来严重的天气、气候灾害(胡钰玲等,2017顾伟宗等,2016杜小玲等,2014刘刚等,2012)。阻塞高压活动集中区域的观测研究虽早于阻塞理论研究,但至今尚未形成较为接近的共识。有关北半球阻塞形势活动的地理分布和频次的统计结果不尽相同,造成差异的主要原因可能是确定阻塞过程的定义及持续天数的标准不同(丁一汇,2005)。由此可见,较为客观的检测和划分北半球阻高活动较为频繁的区域仍是需要研究的科学问题。

北半球阻高中,对中国天气、气候有最直接影响的当属东北亚阻高。因此,有关东北亚阻高频繁活动区变异特征的机制及其前期征兆等科学问题的研究历来受到重视。潘婕等(2008)发现乌拉尔山阻高崩溃后,罗斯贝波从乌拉尔山活动中心向东传播,最终形成贝加尔湖和鄂霍次克海地区的持续异常环流;罗斯贝波的传播、青藏高原及其附近地区暖湿空气的传输和赤道西太平洋海表温度异常激发的遥相关波列等对东亚夏季阻塞过程的形成有重要影响(张存杰等,2004)。另外,欧亚阻高间的配置也一直被众多学者所关注。李峰等(2004)统计表明,6月阻塞活动多发生在乌拉尔山和鄂霍次克海地区,以双阻高为主要形势;1998年双阻高形势下稳定在鄂霍次克海附近的阻高,其前部的东北气流引导高纬度冷空气长驱南下,是副热带高压短期变异的关键(任荣彩等,2004);乌拉尔山—鄂霍次克海双阻高活动在1977年前后发生显著突变,突变后双阻高活动天数呈显著减少,突变前(后)乌拉尔山与鄂霍次克海阻高活动呈现同(反)位向变化(Liu,et al,2017);20世纪90年代以来,夏季乌拉尔山地区阻高处于明显偏弱期,鄂霍次克海阻高处于偏强期,但21世纪初,夏季乌拉尔山地区阻高活动频次开始呈上升趋势,而鄂霍次克海阻高频次却呈明显下降趋势(刘刚等,2012)。在对阻高持续活动的大气内动力过程研究的同时,学者们(Wallace,et al,1981Ding,et al,1987施能,1996谭本馗等,2014)还发现了其与气候低频变率的遥相关型联系。李子仪等(2013)研究的个例分析强调苏门答腊西侧的热带印度洋区域存在明显的对流增强,可能激发了类似东亚-太平洋遥相关型的波列向中高纬度传播,从而调制行星波活动再影响乌拉尔山地区位势高度异常的形成和持续。近些年来,有些学者(Luo,et al,20162017Yao,2016Gong,et al,2017)利用改进的非线性多尺度交互模式,研究了北大西洋急流通过北大西洋涛动正位相影响欧洲阻高的物理机制,乌拉尔山阻高对冬季北极暖—欧亚冷异常的影响,以及北极变暖对冬季乌拉尔山阻高增强的准静止和持续性及欧亚极寒事件的影响,还揭示乌拉尔山阻高对冬季北极海冰减少的促进作用等。以往关于阻高的工作多是基于全年或冬、夏季统计特征或机理的研究(Hartmann,et al,1980Dole,1982Shukla,et al,1983陆日宇等,1998黄菲等,2002鹿翠华等,2010),涉及春夏之交的研究较少。5—6月是北半球大气环流春夏转换的关键期,主要特征是北半球西风带行星波发生一次大调整,东亚大槽东移减弱、副热带西风急流突然北跳,这一阶段的东北亚阻高不仅与东亚夏季风的向北推进有密切联系(刘长征等,2004),还与东北冷涡共同对中国天气、气候产生重要影响(孙力,1997何金海等,2006Lian,et al,2013刘刚等,2015)。

根据以往研究(刘刚等,2019),选用美国环境预报中心/ 国家大气研究中心(NCEP/NCAR-Ⅰ)再分析资料和相同的计算机客观识别方法,尝试划分1979—2018年5—6月北半球阻高频繁活动区域,讨论其活动天数年际与年代际变化的主要特征;并以东北亚阻高为例,采用信息流因果关系等方法,进一步研究其与北半球阻高活动频繁区的可能联系。此外,2019年6月末欧洲出现极端高温热浪,多地刷新最高温纪录,造成火灾和人员伤亡。本研究在1979—2018年气候态分析的基础上,用类似方法对2019年6月北半球阻高活动频繁区进行检测,为进一步揭示阻高与极端天气、气候事件的联系奠定科学基础。

2 资料和主要方法

采用1979—2019年5—6月NCEP/NCAR-Ⅰ逐日500 hPa高度场、温度场等再分析资料,资料水平分辨率为2.5°×2.5°。

2.1 阻高定义及识别

依据Rex(1950)提出的阻高天气学定义和朱乾根等(2000)总结出的阻高存在条件,对中高纬度阻高过程给出如下标准(刘刚等,2012):(1)40°—75°N范围内,500 hPa位势高度场上有高压中心;(2)高压中心日移动不超过12.5个经距;(3)高压中心维持不少于3 d。按照上述标准,应用Fortran语言编程对NCEP逐日500 hPa高度场资料进行检索、筛选,通过确定组中心点的位置和移动,可以直接找出阻高形成、移动和消亡的轨迹,确定其存在的始末位置和维持时间,并统计阻高活动的频次和出现的累计天数,实现计算机的客观自动识别。

2.2 Liang-Kleeman信息流简介

因果关系定量地由两时间序列间单位时间内传递的信息来表征。一个事件到另一个事件存在信息传递,逻辑上表明前者与后者因果联系的强弱,若不存在因果联系则信息流为0。对于线性系统X1X2,从后者X2到前者X1的信息流可近似表示为(Liang,2014

${\rm{T}}{}_{2 \to 1} = \frac{{{C_{11}}{C_{12}}{C_{2,d1}} - C_{12}^2{C_{1,d1}}}}{{C_{11}^2{C_{22}} - {C_{11}}C_{12}^2}}$ (1)

式中,CijXiXj的样本协方差, ${C_{i,dj}}$ Xi ${{\overset{\bullet }{\mathop{\rm{X}}}\,}_{j}}$ 的协方差, ${{\overset{\bullet }{\mathop{\rm{X}}}\,}_{j}}$ ${\rm d}{X_j}/{\rm d}t$ 的差分近似,为欧拉前差格式

${{\overset{\bullet }{\mathop{\rm{X}}}\,}_{j,n}} = \frac{{{X_{j,n + k}} - {X_{j,n}}}}{{k\Delta t}}$ (2)

式中, $\Delta t$ 为时间步长。若 ${{\rm{T}}_{2 \to 1}} = 0$ ,则X2不会引起X1变化,即X2不是X1变化的因;若 ${{\rm{T}}_{2 \to 1}} {\text{>}} 0$ ,则X2X1变化的因,且X2使X1变化的不确定性增强,可预报性降低;若 ${{\rm{T}}_{2 \to 1}} {\text{<}} 0$ ,则X2X1变化的因,且X2使X1变化趋于稳定,可预报性增强。相应地,Liang(20142016)给出了检验信息流显著性的方法,对于大量数据而言,参数的最大似然估计近似服从正态分布,该方法已在研究厄尔尼诺与印度洋偶极子(IOD)因果关系时得到验证。

信息流是Liang(2014)在转移熵和Granger因果检验基础上改进的以一种定量方式来表征变量(或事件)之间因果关系的物理量,通过一个变量序列到另一个变量序列的信息事件速率来衡量因果关系,并且这种因果关系是单向的,两个变量之间的信息交换量不仅说明了因果关系的强弱,而且指明了方向。由于北半球阻高的主要活动区都与各自活动区的海陆分布和地形影响有关,很难判断谁对谁的影响更重要,因而本研究引进Liang-Kleeman信息流方法,尝试揭示北半球初夏季节各阻高频繁活动区间的相互联系究竟谁对谁的影响更重要。

3 北半球阻高频繁活动区域的划分及变化特征 3.1 活动频繁区域的划分

图1给出了由NCEP/NCAR-Ⅰ再分析资料得到的1979—2018年5月和6月北半球中高纬度地区阻高逐日活动累计天数的分布。不难看出,北半球阻高在5月和6月均存在明显的地域差异:欧洲范围内,以乌拉尔山阻塞为主,纬度偏低;东北亚地区以鄂霍次克海阻塞为主,纬度偏高;阿拉斯加和加拿大西部各有一大值中心,纬度位置与东北亚阻塞同样偏高。5月,自欧洲开始经乌拉尔山—青藏高原北区—鄂霍次克海—北太平洋—北美洲—北大西洋Ⅱ区等阻高活动中心纬度分布,具有“正弦波”特征(北大西洋Ⅰ区阻高除外),即较高纬度—中纬度—高纬度—中纬度的半球尺度分布形状。6月阻高活动中心纬度分布亦具有“正弦波”特征(北大西洋Ⅰ区阻高除外),除欧洲大陆上空阻高较5月有所减少外,其他大部分地区阻高明显增多。

图 1  1979—2018年5月(a)和6月(b)北半球阻高活动累计天数分布 (黑点为阻高活动累计天数,蓝线为阻高天数随纬度分布呈“正弦波”特征) Fig. 1  Distribution of cumulative days of the northern hemisphere blocking highs in May (a) and June (b) from 1979 to 2018 (Black spots represent the cumulative days of blocking highs,and blue lines represent a "sine wave" characteristic of the cumulative days of blocking highs)

根据阻高出现的地理位置和活动集中程度占北半球活动总天数的百分比,将5月和6月北半球阻高活动分别划分为8和10个区域。其中,5月阻高频繁活动区包括:①欧洲区(45°—67.5°N,0°—35°E,下文称A区),②乌拉尔山区(45°—62.5°N,37.5°—77.5°E,下文称B区),③青藏高原北区(40°—47.5°N,80°—117.5°E,下文称C区),④东北亚区,也称“鄂霍次克海区”(52.5°—72.5°N,120°—165°E,下文称D区),⑤北太平洋区(50°—75°N,167.5°E—135°W,下文称E区),⑥北美区(47.5°—72.5°N,135°—60°W,下文称F区),⑦北大西洋Ⅰ区(60°—75°N,57.5°—2.5°W,下文称G区),⑧北大西洋Ⅱ区(40°—57.5°N,50°—2.5°W,下文称H区);6月阻高活动频繁区包括:①欧洲区(45°—70°N,0°—42.5°E,下文称A区),②乌拉尔山区(42.5°—62.5°N,45°—75°E,下文称B区),③贝加尔湖区(50°—75°N,77.5°—117.5°E,下文称C区),④青藏高原北区(40°—47.5°N,80°—117.5°E,下文称D区),⑤东北亚区,也称鄂霍次克海区(42.5°—75°N,120°—162.5°E,下文称E区),⑥中太平洋区(57.5°—72.5°N,165°E—167.5°W,下文称F区),⑦阿拉斯加区(60°—75°N,165°—132.5°W,下文称G区),⑧北美区(52.5°—72.5°N,130°—60°W,下文称H区),⑨北大西洋Ⅰ区(60°—75°N,57.5°—10°W,下文称I区),⑩北大西洋Ⅱ区(40°—55°N,42.5°—2.5°W,下文称J区)。

选取1980年6月18—20日(图2)和2009年6月5—7日(图3)两个个例,采用主观识别方法对上述NCEP/NCAR-Ⅰ再分析资料计算得到的青藏高原北阻高和东北亚阻高结果进行验证。可以清晰地在两组图中分别看到,(40°—47.5°N,90°—110°E)和(45°—70°N,120°—160°E)范围内500 hPa高度场上出现了闭合环流,由此证明了计算机客观识别结果的科学合理。特别地,1980年个例显示,在贝加尔湖以南、青藏高原以北地区500 hPa高度场上出现了清晰的阻塞环流,文中称为“青藏高原北阻高”。

图 2  1980年6月18—20日 (a—c) 500 hPa高度场 (单位:gpm) Fig. 2  500 hPa height field on 18—20 (a—c) June 1980 (unit:gpm)
图 3  2009年6月5—7日(a—c)500 hPa高度场 (单位:gpm) Fig. 3  500 hPa height field on 5—7 (a—c) June 2009 (unit:gpm)

表1统计了1979—2018年5月和6月北半球各区域阻高在不同持续时间下的活动频数、40 a总活动天数、以每个经度×纬度为一间隔的平均天数,以及所占北半球总阻高天数的比例。5月和6月各区阻高合计天数分别占北半球总阻高天数的89.96%和89.88%,证明所划分的阻高区域基本代表整个北半球阻高的分布。北半球阻高在不同区域的活跃程度差异较大,各区阻高的持续天数多在7 d以内。5—6月北半球阻高活动最明显的变化是,5月北太平洋、北美洲和乌拉尔山阻高活动占前3位,6月东北亚阻高活动天数跃升第1位;欧洲和乌拉尔山区阻高减少较为显著。另外,尽管6月青藏高原北阻高活动天数与北半球阻高总体天数相比并不占优势,但在每个经度×纬度的平均天数为0.55 d,是(除阿拉斯加阻高外)北半球最高的,这点很值得关注!

表 1  1979—2018年5月和6月北半球各区阻高活动在不同天数下发生的频次 Table 1  Frequency of the northern Hemisphere blocking high activities in different days in May and June from 1979 to 2018
区域 持续天数(d) 合计天数(d) 平均天数
(d/(经度×纬度))
所占比重(%)
3 4 5 6 7 8 9 10 11 12 13 14 15 16 17
5月 A 22 18 12 8 3 3 2 309 0.37 10.61
B 22 15 13 12 8 1 3 2 1 392 0.52 13.46
C 16 5 4 1 1 1 119 0.36 4.09
D 48 19 13 8 3 1 362 0.37 12.43
E 41 33 20 6 6 2 1 1 472 0.31 16.21
F 50 35 7 4 5 6 2 1 1 471 0.24 16.17
G 27 15 6 3 3 3 1 243 0.27 8.34
H 11 20 9 4 3 2 2 1 252 0.28 8.65
6月 A 23 13 9 9 4 1 1 1 1 290 0.26 7.85
B 17 11 9 7 8 2 5 1 2 331 0.51 8.96
C 26 13 9 7 7 3 1 1 2 2 1 374 0.35 10.12
D 19 8 8 3 2 1 1 181 0.55 4.90
E 48 32 17 14 12 6 6 4 2 691 0.47 18.70
F 23 17 6 2 4 2 1 234 0.51 6.33
G 28 13 9 6 2 3 2 1 1 1 306 0.57 8.28
H 38 19 15 7 6 3 4 409 0.27 11.07
I 25 14 13 2 1 2 1 1 259 0.33 7.01
J 22 7 7 4 4 3 1 2 1 246 0.38 6.66

对1979—2018年逐年5月和6月各频繁区阻高活动天数与500 hPa平均位势高度场进行相关分析(图45)发现,各类阻高年活动累计天数与500 hPa高度场均存在相应的显著相关区,其中心区显著性水平均超过95%的信度t检验,这不仅可以进一步验证北半球阻高活动检测和划分结果的科学性,更主要可以揭示阻高各主要活动区引起北半球500 hPa高度场异常的不同特征。

图 4  1979—2018年5月北半球阻高活动累计天数与500 hPa高度场的相关系数场 (a—h分别表示A、B、C、D、E、F、G、H区域阻高;阴影区通过95%信度t检验) Fig. 4  Correlation coefficient maps between cumulative days of various types of blocking high and average 500 hPa height in May from 1979 to 2018 (a—h represent types A,B,C,D,E,F,G and H respectively;shaded areas are for values at the 95% confidence level)
图 5  1979—2018年6月北半球阻高活动累计天数与500 hPa高度场的相关系数场 (a—j分别表示A、B、C、D、E、F、G、H、I、J区域阻高;阴影区通过95%信度t检验) Fig. 5  Correlation coefficient maps between cumulative days of various types of blocking high and average 500 hPa height in June from 1979 to 2018 (a—j represent types A,B,C,D,E,F,G,H,I and J respectively;shaded areas are for values at the 95% confidence level)

5月,首先,只有乌拉尔山阻高频繁活动能引起整个北半球500 hPa高度场异常,欧洲至亚洲中西部地区和北美洲地区呈现相似的鞍型显著相关场(由北向南为“+”,由西向东为“−”);其次,北大西洋Ⅰ区阻高频繁活动引起北极区与东亚至北大西洋的中高纬度地区为自北向南“+−”异常的半球尺度相关波列;第三,北太平洋阻高频繁活动亦引起北美洲自西向东的纬向“+−”相关波列,以及与其南侧自亚洲西部至太平洋东部的纬向“−”相关显著波状带;其他阻高活动区只产生局地显著相关区。

6月,第一,仍然是乌拉尔山阻高频繁活动能引起整个北半球500 hPa高度场异常,分别引起欧洲、鄂霍次克海至西北太平洋和北美洲中高纬度地区且具有半球尺度的“+”遥相关(图5b);第二,青藏高原北区阻高频繁活动引起斯堪纳维亚半岛附近北冰洋至青藏高原北麓显著的西北—东南向“−+”相关波列,还引起青藏高原下游副热带地区的纬向波状“−”相关带(图5d),而贝加尔湖阻高频繁活动区则相反,形成由地中海—欧洲—贝加尔湖呈西南—东北向“+−+”相关波列(图5c);第三,东北亚阻高频繁活动引起与青藏高原类似的东北亚至中太平洋副热带地区自西北向东南的“+−+”显著相关波列,还引起北美洲与其南面副热带区域的经向较显著的“+−”遥相关波列;其他阻高频繁活动区大都引起局域显著相关区,其中北大西洋Ⅰ和Ⅱ区,以及贝加尔湖阻高频繁活动区所形成的“+”与“−”相间经向相关波列较为引人关注。

3.2 主要变化特征

从5—6月北半球阻高天数的标准化时间序列(图67)可以看出,各区阻高40 a总体线性趋势均不显著,但阶段性变化明显。从谐波分析(周晓霞等,2007Li,et al,2004)分离出的年际和年代际分量结果的方差贡献看,各区阻高活动天数的年际分量占显著优势(均在55%以上)。这表明,初夏北半球各区阻高活动天数的年际变化分量起主要作用,预测5—6月北半球阻高活动的关键在年际尺度上。

图 6  1979—2018年5月北半球阻高天数的标准化时间序列 (a—h分别表示A、B、C、D、E、F、G、H区域阻高;柱状为原序列,虚线为年际分量,实线为年代际分量) Fig. 6  Normalized time series of the northern hemisphere blocking high days in May from 1979 to 2018 (a-h represent types A,B,C,D,E,F,G and H, respectively;the columns are for the original sequence,the dotted line denotes the interannual component,and the solid line represents the interdecadal component)
图 7  1979—2018年6月北半球阻高天数的标准化时间序列 (a—j分别表示A、B、C、D、E、F、G、H、I、J区域阻高;柱状为原序列,虚线为年际分量,实线为年代际分量) Fig. 7  Normalized time series of the northern hemisphere blocking high days in June from 1979 to 2018 (a-j represent types A,B,C,D,E,F,G,H,I and J, respectively;the columns are for the original sequence,the dotted line shows the interannual component,and the solid line shows the interdecadal component)

5月和6月北半球各阻高频繁活动区的年代际变化和周期性各不相同。需要特别指出的年代际变化特点是:进入21世纪以来,欧洲区阻高活动5月偏多而6月却呈明显偏少,青藏高原区阻高活动变化正相反,5月明显偏少而6月偏多;21世纪最初10年以来,6月北半球10个阻高频繁活动区有6个(欧洲、乌拉尔山、鄂海、阿拉斯加、北大西洋Ⅰ和Ⅱ区)阻高活动天数年代际分量明显偏少。

4 东北亚阻高活动天数同北半球其他活动区阻高间的联系

进一步统计5—6月东北亚阻高活动天数与北半球其他区域阻高的关系,由于资料时间长度所限,主要分析东北亚阻高与北半球其他区域阻高的年际分量关系。值得注意的是,文中涉及到的时间序列长度较短,故后文统计的信息流结果普遍偏小。

从1979—2018年5月东北亚阻高与5—6月北半球其他区域阻高活动天数的相关关系(表2)上看,5月东北亚阻高天数仅与6月欧洲阻高天数的年际分量存在显著负相关。由信息流结果(表3)可见,5月只有青藏高原北阻高天数是同期东北亚阻高天数的稳定影响源;年际分量上,5月东北亚阻高天数是同期北大西洋Ⅱ区阻高和6月北大西洋Ⅰ区阻高及贝加尔湖阻高天数的稳定影响源,是6月北美、北大西洋Ⅱ区和青藏高原北阻高的不稳定影响源。由此可以看到,年际分量上,5月东北亚阻高天数不受同期北半球其他区域阻高活动的影响,是明显相对独立的;同时,5月东北亚阻高的年际分量又是5—6月北半球其他区域阻高的影响源。特别地,5月青藏高原北阻高是5月东北亚阻高天数的稳定影响源,而5月东北亚阻高又是6月青藏高原北阻高天数的不稳定影响源。

表 2  1979—2018年5月东北亚阻高与北半球其他区域阻高活动天数的相关系数 Table 2  Correlation coefficient between days of the Northeast Asia blocking high and blocking high in other regions of the northern Hemisphere in May from 1979 to 2018
相关系数 北美 北大西洋Ⅰ 北大西洋Ⅱ 欧洲 乌拉尔山 贝加尔湖 青藏高原北 东北亚 北太平洋
5月 原序列 0.04 −0.28 0.15 0.16 0.19 0.12 1 −0.28
年际 0 −0.25 0.20 0.02 0.13 0.18 1 −0.19
6月 原序列 0.19 −0.02 0.05 −0.28 0.07 −0.04 −0.08 −0.05 0.17
年际 0.23 −0.01 0.09 −0.44* 0.23 −0.17 −0.06 0.05 0.24
 *表示通过99%信度t检验。
表 3  1979—2018年5月东北亚阻高(O)与北半球其他区域阻高(X)活动天数的信息流 Table 3  Information flow of the Northeast Asia blocking high(O)days and blocking high(X)days in other areas of the northern hemisphere in May from 1979 to 2018
信息流结果 北美 北大西洋Ⅰ 北大西洋Ⅱ 欧洲 乌拉尔山 贝加尔湖 青藏高原北 东北亚 北太平洋
5月 TO→X 原序列 0.0106 −0.0023 −0.0460 0.0088 −0.0123 −0.0207 0.0422
年际 0.0059 0.0098 −0.0531* −0.0098 −0.0135 8.15×10-4 0.0175
TX→O 原序列 −0.0011 −0.0380 −0.0200 −0.0065 0.0181 −0.0342* 0.0428
年际 −0.0052 −0.0398 0.0203 0.0039 0.0022 4.49×10-4 0.0243
6月 TO→X 原序列 0.0301 −0.0080 0.0206* −0.0600 −0.0046 −0.0026** 0.0421** −0.0038 −0.0129
年际 0.0428* −0.0111* 0.0314* −0.0608 −0.0130 −0.0678** 0.0186** 0.0079 −0.0270
 *表示通过95%信度t检验,**表示通过99%信度t检验。

线性相关分析(表4)显示,6月东北亚阻高天数与5月乌拉尔山阻高天数呈显著负相关;年际分量上,6月东北亚阻高天数与5月北大西洋Ⅱ区、6月青藏高原北和北太平洋阻高天数呈显著负相关。注意到,6月东北亚阻高与同期乌拉尔山和贝加尔湖阻高年际分量关系不显著,这可能与北半球哈得来环流和费雷尔环流等活动有关。信息流结果(表5)表明,5月北大西洋Ⅰ区阻高天数是6月东北亚阻高天数的不稳定影响源;年际分量上,只有5月北太平洋阻高天数是6月东北亚阻高天数的稳定影响源。

表 4  1979—2018年6月东北亚阻高与北半球其他区域阻高活动天数的相关系数 Table 4  Correlation coefficient between days of the Northeast Asia blocking high and day of blocking high in other regions of the northern Hemisphere in June from 1979 to 2018
相关系数 北美 北大西洋Ⅰ 北大西洋Ⅱ 欧洲 乌拉尔山 贝加尔湖 青藏高原北 东北亚 北太平洋
5月 原序列 0 0.12 −0.27 0.02 −0.32* 0.27 −0.05 0.08
年际 −0.02 −0.09 −0.44** 0.24 −0.28 0.29 0.05 0.03
6月 原序列 0.23 −0.19 0.08 0.20 0.22 −0.20 −0.30 1 −0.31
年际 0.25 −0.26 −0.01 0.29 −0.14 0.10 −0.36* 1 −0.33*
 *表示通过95%信度t检验,**表示通过99%信度t检验。
表 5  1979—2018年6月东北亚阻高(O)与北半球其他区域阻高(X)活动天数的信息流 Table 5  Information flow of the Northeast Asia blocking high(O)days and the blocking high(X)days in other areas of the northern hemisphere in June from 1979 to 2018
信息流结果 北美 北大西洋Ⅰ 北大西洋Ⅱ 欧洲 乌拉尔山 贝加尔湖 青藏高原北 东北亚 北太平洋
5月 TX→O 原序列 0.0011 0.0386* −0.0439 −0.0032 −0.0288 0.0182 −0.0038 −0.0151
年际 4.21×10-5 −0.0154 −0.0870 −0.0194 −0.0101 0.0367 0.0079 −0.0129*
6月 TO→X 原序列 −0.0206 2.59×10-4 0.0054 0.0147 0.0284 0.0487 0.0189 0.0590
年际 −0.0567 −0.0227 −0.0011 0.0263 0.0186 0.0014 0.0587 0.0292
TX→O 原序列 0.0116 −0.0287 −0.0020 −0.0501 0.0590 0.0599 −0.0290 −0.0034
年际 0.0701 −0.0414 0.0032 −0.0747 −0.0200 −0.0014 −0.0127 −0.0294
 *表示通过95%信度t检验,**表示通过99%信度t检验。

从以上的讨论中不难发现,6月东北亚阻高天数的年际分量很大程度上受5月北太平洋阻高活动影响。初夏东北亚阻高和北太平洋阻高活动均会产生一定的大气遥相关,由上游至下游传来的遥相关异常很可能为气候预测提供有用的信号。为此,分别统计5月北太平洋阻高活动天数和6月东北亚阻高活动天数与同期北半球10种大气遥相关指数年际分量的关系。其中,所使用的大气遥相关指数为每月北半球500 hPa高度异常场旋转主成分分析的前10个主要模态,指数序列源自美国气候预测中心(http://www.cpc.ncep.noaa.gov/data)。

由图8可以看到,5月北太平洋阻高天数与同期东太平洋/北太平洋(EP/NP)型遥相关相似,二者年际分量的相关系数达0.50(通过0.01显著性t检验);图9可以看到,6月东北亚阻高天数与同期西太平洋(WP)型遥相关相似,二者年际分量相关系数达−0.48(通过0.01显著性t检验)。近40年中,共有10 a 5月东太平洋/北太平洋型正位相转为6月西太平洋型负位相,其中有7 a 6月东北亚阻高天数偏多;共有9 a 5月东太平洋/北太平洋型负位相转为6月西太平洋型正位相,其中有8 a 6月东北亚阻高天数偏少。当5月北太平洋频繁活动区阻高天数偏多时,在500 hPa高度场上,形成以北太平洋正距平为中心,在其南侧经亚洲西侧至北美太平洋西海岸负距平带状区,横跨亚洲至北美的副热带的显著相关区域,这种在东北亚地区由北向南的“+−”显著相关波列(图8a),恰为6月东北亚阻高频繁活动区阻高天数偏多与同期500 hPa高度场相关形势(图9a)创造有利条件,若5月北太平洋阻高所引起的波列显著西退,很容易形成6月东北亚阻高偏多的形势。

图 8  1979—2018年5月北太平洋阻高天数 (a) 和东太平洋/北太平洋型遥相关指数 (b) 与同期500 hPa高度场年际分量的相关系数场(阴影表示通过95%信度t检验) Fig. 8  Interannual component correlation coefficient maps between cumulative days of North Pacific blocking high (a),EP/NP index (b) and average 500 hPa height in May from 1979 to 2018 (shaded areas are for values at the 95% confidence level)
图 9  1979—2018年6月东北亚阻高天数 (a) 和西太平洋型遥相关指数 (b) 与同期500 hPa高度场年际分量的相关系数场(阴影表示通过95%信度t检验) Fig. 9  Interannual component correlation coefficient maps between cumulative days of Northeast Asia blocking high (a),WP index (b) and average 500 hPa height in June from 1979 to 2018 (shaded areas are for values at the 95% confidence level)
5 2019年6月北半球阻高与极端高温事件

2019年6月26—30日,高温热浪席卷欧洲,多地刷新最高温度纪录,还引发多起火灾、加剧空气污染,并造成人员伤亡。从2019年6月北半球地表温度距平分布(图10a)上不难发现,2019年6月欧洲中纬度地区、贝加尔湖—东北亚—北太平洋的高纬度地区、以及北美东部—北大西洋的高纬度地区的地表温度距平值异常偏高。2019年6月500 hPa高度场距平分布(图10b)与1979—2018年6月欧洲阻高活动累计天数和500 hPa高度场的相关系数分布(图5a)相似,中高纬度异常的高压系统稳定维持在欧洲、贝加尔湖—东北亚—北太平洋的高纬度地区、以及北美东部—北大西洋的高纬度地区上空,也就是说上述高温地区均被异常高压系统所占据,持续的下沉气流造成连续多日晴热高温天气。

图 10  2019年6月北半球地表温度 (a) 及500 hPa高度场距平 (b) Fig. 10  Surface temperature in northern Hemisphere in June 2019 (a) and 500 hPa height field anomaly (b)

2019年6月,欧亚中高纬度地区阻高活动频繁,其中欧洲阻高和东亚阻高活动天数异常偏多,分别为19和21 d,较常年平均值分别多12和5 d,反映了阻塞高压作为中高纬度地区典型的持续性环流异常是导致天气、气候异常的重要大尺度环流条件。对2019年6月阻高分布情况(图11)进行统计发现,北半球阻高整体有偏南的趋势,唯独东北亚阻高偏北;具体而言,欧洲阻高异常偏多偏南,且分布集中;东北亚阻高异常偏多偏东偏北,但主要分布在中国东北北部和俄罗斯东北部楚科奇山脉附近,这与以往认知的东北亚阻高有所不同。此外,青藏高原阻高偏少,贝加尔湖阻高异常偏南。以往研究(Black,et al,2004Stefanon,et al,2012王倩云等,2016)表明,欧洲夏季热浪的发生与同期欧洲大陆的持续反气旋环流(阻高)有关,阻高的位置决定了热浪发生的区域。从阻高活动过程来看,欧洲阻高共出现了5次,分别在6月5—9日、13—15日、16—24日、23—25日、25—27日;从阻高的演变趋势来看,可以将上述阻高过程归纳为3次阻高活动过程,第一次为6月5—12日,源自乌拉尔山阻高西退、发展、东移、逐渐减弱和崩溃过程;第二次为13—25日,源自欧洲东部阻高东移、发展、稳定、减弱、崩溃;第三次为23—30日,欧洲西部产生阻高、发展、东移,而后在更西部又产生新的阻高,继续发展、东移、后期减弱、崩溃。3次阻高活动过程中,最后一次阻高中心平均强度为5928 gpm,前两次阻高中心平均强度分别为5843和5831 gpm,相比之下,最后一次阻高活动强度明显高于前两次。

图 11  2019年6月北半球阻高活动累计天数 (黑点) 分布 Fig. 11  Distribution of cumulative days (black spot) of the northern Hemisphere blocking highs in June 2019
6 结论与讨论

为深入研究初夏东北亚阻高频繁活动区的变异特征及其前期征兆,有必要考虑整个北半球阻高活动区的变异及其之间的联系。文中基于美国NCEP/NCAR-Ⅰ再分析资料,按照国际通行的阻高定义,采用计算机自动识别和天气图人工主、客观相结合的检测方法,获得了40 年来5月和6月逐日北半球中高纬度阻高活动的无缝隙分布;并利用线性相关和信息流因果关系等统计方法,研究了东北亚阻高与北半球其他频繁活动区阻高年际变化的联系。主要结论如下:

(1)5月阻高主要活动区按活动集中程度与占北半球活动总天数百分比选定为8个,排序如下:北太平洋区(占16.21%)、北美区(占16.17%)、乌拉尔山区(占13.46%)、鄂霍次克海区(占12.43%)、欧洲区(占10.61%)、北大西洋Ⅱ区(占8.65%)、北大西洋Ⅰ区(占8.34%);此外,将出现在贝加尔湖南侧(40°—47.5°N,80°—117.5°E)区域的阻高,称为“青藏高原北”阻高,占5月北半球总阻高天数的4.09%。

6月阻高主要活动区按活动集中程度与占北半球活动总天数百分比选定为10个:鄂霍次克海区(占18.70%)、北美区(占11.07%)、贝加尔湖区(占10.21%)、乌拉尔山区(占8.96%)、阿拉斯加区(占8.28%)、欧洲区(占7.85%)、北大西洋Ⅰ区(占7.01%)、北大西洋Ⅱ区(占6.66%)、中太平洋区(占6.33%)、青藏高原北区(占4.90%)。尽管青藏高原北阻高只占北半球阻高活动天数的4.90%,但是其活动中心集聚程度很高。

将文中检测结果与20世纪70—90年代有关北半球阻塞高压活动频繁区的观测进行对比:Treidl等(1981)发现北半球阻塞主要产生在北太平洋和北大西洋上,最高频率发生在欧洲和北美洲西海岸附近地区,即大西洋与太平洋上天气尺度扰动活跃区的下游(丑纪范等,1994)。White等(1975)统计1950—1970年北太平洋中部阻塞高压在秋冬季(10—3月)更活跃,春夏季(4—9月)出现机会少。统计发现,北太平洋中部阻高仅5月就达到472 d,出现频率占北半球总阻高的16.21%,远比White等(1975)所给出的北太平洋春季阻高活动频率百分之几高得多;Sumner(1954)用完全主观的方法指出欧洲阻塞活动5月最频繁,这与本研究检测发现5月比6月多是一致的。

(2)初夏北半球阻高活动分布特征与月尺度变异特征表明,5月,自欧洲开始经乌拉尔山—青藏高原北区—鄂霍次克海—北太平洋—北美洲—北大西洋Ⅱ区等阻高活动中心纬度分布(北大西洋Ⅰ区阻高除外),具有“正弦波”特征,即较高纬度—中纬度—高纬度—中纬度的半球尺度分布形状。6月阻高活动中心纬度分布亦具有“正弦波”特征,可能揭示阻高的形成与维持的大气低频准定常波路径高、低纬“反射”有关(徐祥德等,2002)。另外,5—6月北半球阻高活动最明显的变化是,东北亚阻高活动天数跃升北半球第1位;欧洲和乌拉尔山区阻高减少较为显著;青藏高原北区阻高活动尽管天数与北半球总体天数相比并不占优势,但是在每1经度×纬度的平均天数为0.61 d,是北半球(除阿拉斯加阻高外)最高的,很值得关注。

(3)5月和6月只有乌拉尔山阻高频繁活动能引起整个北半球500 hPa高度场异常。青藏高原北区阻高频繁活动引起斯堪纳维亚半岛附近北冰洋至青藏高原北麓的显著“−+”相关波列,还引起青藏高原下游副热带地区的纬向波状“−”相关带;东北亚阻高频繁活动引起与青藏高原类似的东北亚至中太平洋副热带地区自西北向东南的“+−+”显著相关波列,还引起北美洲与其南面副热带区域的经向较显著的“+−”遥相关波列;其他阻高频繁活动区大都引起局域高度场的显著变化。魏东等(2006)定义了一个鄂霍次克海阻高指数,给出了活动强、弱年以及和同期500 hPa位势高度场的相关系数分布,本研究结果与他们的结论类似。

(4)本研究从气候态统计出发,认为初夏东北亚(鄂霍次克海)阻高年际变化明显,且具有相对独立性。5月东北亚阻高活动天数的年际变化仅与6月欧洲阻高存在显著的负相关,与其他北半球阻高频繁活动区天数年际变化没有任何同期相互影响的显著相关和信息流因果关系。值得注意的是,5月东北亚阻高活动天数的年际变化不仅是5月北大西洋Ⅱ区阻高、6月北大西洋Ⅰ区和贝加尔湖阻高天数年际变化的稳定影响源,还是6月北美、北大西洋Ⅱ区和青藏高原北阻高天数年际变化的不稳定影响源。6月东北亚阻高活动天数的年际变化仅与同期青藏高原北阻高存在显著负相关,还只与5月北大西洋Ⅱ区存在滞后的负相关;除5月北太平洋区阻高活动是6月东北亚阻高活动年际分量的稳定影响源外,6月东北亚阻高活动天数的年际分量并不受初夏北半球其他阻高活动影响。

文中对初夏北半球阻高活动天数的主要变化特征及东北亚阻高与北半球其他区域阻高活动天数的联系进行了初步分析和讨论。此外,回顾了2019年6月欧洲地区出现的极端高温事件,初步诊断分析认为欧洲阻高异常偏多且位置偏南可能是影响高温热浪的主要因子之一。值得思考的是,6月东北亚阻高天数的年际分量受5月北太平洋阻高活动的影响,但这种影响通过什么途径产生?中高纬度遥相关和外强迫对初夏东北亚阻高变异影响如何协同作用?研究中发现,在贝加尔湖以南、青藏高原以北地区500 hPa上空也存在着阻高,这种阻高如何形成?它与东北亚阻高有何联系,对中国天气、气候又有何影响?这些问题在今后的工作中将深入探讨。另外,文中主要探讨东北亚阻高与北半球其他阻高的关系,而北半球其他各阻高间的关系尚待进一步研究。

参考文献
丑纪范, 刘式达, 刘式适. 1994. 非线性动力学. 北京: 气象出版社, 200pp. Chou J F, Liu S D, Liu S S. 1994. Nonlinear Dynamics. Beijing: China Meteorological Press, 200pp (in Chinese)
丁一汇. 2005. 高等天气学. 北京: 气象出版社, 523pp. Ding Y H. 2005. Advanced Synoptic Meteorology. Beijing: China Meteorological Press, 523pp (in Chinese)
杜小玲, 高守亭, 彭芳. 2014. 2011年初贵州持续低温雨雪冰冻天气成因研究. 大气科学, 38(1): 61-72. DOI:10.3878/j.issn.1006-9895.2013.12119
顾伟宗, 陈丽娟, 左金清等. 2016. 多因子协同作用对1992年和1998年黄淮地区夏季降水异常的影响. 大气科学, 40(4): 743-755.
何金海, 吴志伟, 江志红等. 2006. 东北冷涡的“气候效应”及其对梅雨的影响. 科学通报, 51(23): 2803-2809. He J H, Wu Z W, Jiang Z H, et al. 2006. “Climate effect” of the northeast cold vortex and its influences on Meiyu. Chinese Sci Bull, 52(5): 671-679.
胡钰玲, 赵中军, 康彩燕等. 2017. 中国南方2008年与2016年两次低温雨雪冰冻天气过程对比分析. 冰川冻土, 39(6): 1180-1191.
黄菲, 姜治娜. 2002. 欧亚大陆阻塞高压的统计特征及其与中国东部夏季降水的关系. 青岛海洋大学学报, 32(2): 186-192. DOI:10.3969/j.issn.1672-5174.2002.02.004
李峰, 丁一汇. 2004. 近30年夏季亚欧大陆中高纬度阻塞高压的统计特征. 气象学报, 62(3): 347-354. DOI:10.3321/j.issn:0577-6619.2004.03.009
李子仪, 李双林, 李国平. 2013. 夏季乌拉尔长阻塞过程及其与热带热力异常关系个例分析. 大气科学, 37(3): 731-744. DOI:10.3878/j.issn.1006-9895.2012.12006
刘长征, 王会军, 姜大膀. 2004. 东亚季风区夏季风强度和降水的配置关系. 大气科学, 28(5): 700-712. DOI:10.3878/j.issn.1006-9895.2004.05.05
刘刚, 沈柏竹, 廉毅等. 2012. 亚洲阻塞高压分类及其与东北冷涡活动和东北夏季低温的联系. 地理科学, 32(10): 1269-1274.
刘刚, 廉毅, 颜鹏程等. 2015. 5~8月东北冷涡客观识别分类及北半球大气环流特征分析. 地理科学, 35(8): 1042-1050.
刘刚, 徐士琦, 廉毅. 2019. 夏季亚洲阻塞高压识别及其对中国东北气候异常的可能影响: 不同再分析资料对比. 气象学报, 77(2): 303-314.
鹿翠华, 罗德海. 2010. 近几十年北半球孤立子阻塞和偶极子阻塞活动的对比研究. 中国海洋大学学报, 40(8): 11-17.
陆日宇, 黄荣辉. 1998. 东亚-太平洋遥相关型波列对夏季东北亚阻塞高压年际变化的影响. 大气科学, 22(5): 727-734. DOI:10.3878/j.issn.1006-9895.1998.05.07
潘婕, 布和朝鲁, 纪立人等. 2008. 夏季欧亚中高纬环流持续异常事件的Rossby波传播特征. 大气科学, 32(3): 615-628. DOI:10.3878/j.issn.1006-9895.2008.03.17
任荣彩, 刘屹岷, 吴国雄. 2004. 中高纬环流对1998年7月西太平洋副热带高压短期变化的影响机制. 大气科学, 28(4): 571-578. DOI:10.3878/j.issn.1006-9895.2004.04.09
施能. 1996. 近40年东亚冬季风强度的多时间尺度变化特征及其与气候的关系. 应用气象学报, 7(2): 175-182.
孙力. 1997. 东北冷涡持续活动的分析研究. 大气科学, 21(3): 297-307. DOI:10.3878/j.issn.1006-9895.1997.03.06
谭本馗, 陈文. 2014. 中高纬度大气遥相关动力学及其对东亚冬季气候影响的研究进展. 气象学报, 72(5): 908-925.
王倩云, 罗德海, 王佳妮. 2016. 欧洲地区夏季热浪的特征及其与阻塞环流的联系. 气候与环境研究, 21(4): 367-379.
魏东, 王亚非, 董敏. 2006. 日本东部附近海域海温异常对鄂霍次克海高压的影响. 气象学报, 64(4): 518-526. DOI:10.3321/j.issn:0577-6619.2006.04.012
徐祥德, 高守亭. 2002. 外源强迫与波流作用动力学原理. 2版. 北京: 海洋出版社, 242-243. Xu X D, Gao S T. 2002. Exogenous Forcing and Wave-current Interaction Dynamics. 2nd ed. Beijing: China Ocean Press, 242-243 (in Chinese)
张存杰, 宋连春, 李耀辉. 2004. 东亚地区夏季阻塞过程的研究进展. 气象学报, 62(1): 119-127.
周晓霞, 王盘兴, 段明铿等. 2007. 我国季平均气温和降水局地同时相关的时空特征. 应用气象学报, 18(5): 601-609. DOI:10.3969/j.issn.1001-7313.2007.05.004
朱乾根, 林锦瑞, 寿绍文等. 2000. 天气学原理和方法. 3版. 北京: 气象出版社, 187pp. Zhu Q G, Lin J R, Shou S W, et al. 2000. Weather Principle and Method. 3rd ed. Beijing: China Meteorological Press, 187pp (in Chinese)
Black E, Blackburn M, Harrison G, et al. 2004. Factors contributing to the summer 2003 European heatwave. Weather, 59(8): 217-223. DOI:10.1256/wea.74.04
Ding Y H, Krishnamurti T N. 1987. Heat budget of the Siberian High and the winter monsoon. Mon Wea Rev, 115(10): 2428-2449. DOI:10.1175/1520-0493(1987)115<2428:HBOTSH>2.0.CO;2
Dole R M. 1982. Persistent anomalies of the extratropical northern hemisphere wintertime circulation[D]. Massachusetts: Massachusetts Institute of Technology
Gong T T, Luo D H. 2017. Ural blocking as an amplifier of the Arctic sea ice decline in winter. J Climate, 30(7): 2639-2654. DOI:10.1175/JCLI-D-16-0548.1
Hartmann D L, Ghan S J. 1980. A statistical study of the dynamics of blocking. Mon Wea Rev, 108(8): 1144-1159. DOI:10.1175/1520-0493(1980)108<1144:ASSOTD>2.0.CO;2
Li L P, Wang P X, Li H. 2004. Interdecadal and interannual variabilities of air and sea and their relations over the Pacific. Acta Meteor Sinica, 18(2): 227-244.
Lian Y, Shen B Z, Li S F, et al. 2013. Impacts of polar vortex, NPO, and SST configurations on unusually cool summers in Northeast China. Part Ⅰ: Analysis and diagnosis. Adv Atmos Sci, 30: 193-209. DOI:10.1007/s00376-012-1258-x
Liang X S. 2014. Unraveling the cause-effect relation between time series. Phys Rev E, 90(5): 052150. DOI:10.1103/PhysRevE.90.052150
Liang X S. 2016. Information flow and causality as rigorous notions ab initio . Phys Rev E, 94(5): 052201. DOI:10.1103/PhysRevE.94.052201
Liu G, Wang T, Yang X Y, et al. 2017. Climate characteristics of abnormal double-blocking activities over the Ural mountains and sea of Okhotsk. J Meteor Res, 31(4): 694-707. DOI:10.1007/s13351-017-6048-z
Luo D H, Xiao Y Q, Yao Y, et al. 2016. Impact of Ural blocking on winter warm Arctic-cold Eurasian anomalies. Part Ⅰ: Blocking-induced amplification. J Climate, 29(11): 3925-3947. DOI:10.1175/JCLI-D-15-0611.1
Luo D H, Yao Y, Dai A, et al. 2017. Increased quasi stationarity and persistence of winter Ural blocking and Eurasian extreme cold events in response to arctic warming. Part Ⅱ: A theoretical explanation. J Climate, 30(10): 3569-3587. DOI:10.1175/JCLI-D-16-0262.1
Rex D F. 1950. Blocking action in the middle troposphere and its effect upon regional climate: Ⅰ. An aerological study of blocking action. Tellus, 2(3): 196-211.
Shukla J, Mo K C. 1983. Seasonal and geographical variation of blocking. Mon Wea Rev, 111(2): 388-402. DOI:10.1175/1520-0493(1983)111<0388:SAGVOB>2.0.CO;2
Stefanon M, D’Andrea F, Drobinski P. 2012. Heatwave classification over Europe and the Mediterranean region. Environ Res Lett, 7(1): 014023. DOI:10.1088/1748-9326/7/1/014023
Sumner E J. 1954. A study of blocking in the Atlantic-European sector of the northern hemisphere. Quart J Roy Meteor Soc, 80(345): 402-416. DOI:10.1002/qj.49708034510
Treidl R A, Birch E C, Sajecki P. 1981. Blocking action in the northern hemisphere: A climatological study. Atmos Ocean, 19(1): 1-23. DOI:10.1080/07055900.1981.9649096
Wallace J M, Gutzler D S. 1981. Teleconnections in the geopotential height field during the Northern Hemisphere winter. Mon Wea Rev, 109(4): 784-812. DOI:10.1175/1520-0493(1981)109<0784:TITGHF>2.0.CO;2
White W B, Clark N E. 1975. On the development of blocking ridge activity over the central North Pacific. J Atmos Sci, 32(3): 489-502. DOI:10.1175/1520-0469(1975)032<0489:OTDOBR>2.0.CO;2
Yao Y, Luo D H, Dai A G, et al. 2016. The positive North Atlantic oscillation with downstream blocking and middle east snowstorms: Impacts of the North Atlantic Jet. J Climate, 29(5): 1853-1876. DOI:10.1175/JCLI-D-15-0350.1