中国气象学会主办。
文章信息
- 徐士琦, 刘刚, 廉毅, 曲绍君, 杨雪艳, 梁洪海. 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 改回
2. 吉林省气候中心,长春,130062;
3. 吉林省气象信息网络中心,长春,130062;
4. 吉林省气象科学研究所,长春,130062;
5. 吉林省气象探测保障中心,长春,130062
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
影响中国天气、气候异常的因子除来自低纬度外强迫外,中高纬度大气环流尤其是阻塞高压(阻高)亦有着举足轻重的作用。由于阻高在大气中的持续时间较长且位置稳定少变,它的出现往往通过上下游效应给周边地区带来严重的天气、气候灾害(胡钰玲等,2017;顾伟宗等,2016;杜小玲等,2014;刘刚等,2012)。阻塞高压活动集中区域的观测研究虽早于阻塞理论研究,但至今尚未形成较为接近的共识。有关北半球阻塞形势活动的地理分布和频次的统计结果不尽相同,造成差异的主要原因可能是确定阻塞过程的定义及持续天数的标准不同(丁一汇,2005)。由此可见,较为客观的检测和划分北半球阻高活动较为频繁的区域仍是需要研究的科学问题。
北半球阻高中,对中国天气、气候有最直接影响的当属东北亚阻高。因此,有关东北亚阻高频繁活动区变异特征的机制及其前期征兆等科学问题的研究历来受到重视。潘婕等(2008)发现乌拉尔山阻高崩溃后,罗斯贝波从乌拉尔山活动中心向东传播,最终形成贝加尔湖和鄂霍次克海地区的持续异常环流;罗斯贝波的传播、青藏高原及其附近地区暖湿空气的传输和赤道西太平洋海表温度异常激发的遥相关波列等对东亚夏季阻塞过程的形成有重要影响(张存杰等,2004)。另外,欧亚阻高间的配置也一直被众多学者所关注。李峰等(2004)统计表明,6月阻塞活动多发生在乌拉尔山和鄂霍次克海地区,以双阻高为主要形势;1998年双阻高形势下稳定在鄂霍次克海附近的阻高,其前部的东北气流引导高纬度冷空气长驱南下,是副热带高压短期变异的关键(任荣彩等,2004);乌拉尔山—鄂霍次克海双阻高活动在1977年前后发生显著突变,突变后双阻高活动天数呈显著减少,突变前(后)乌拉尔山与鄂霍次克海阻高活动呈现同(反)位向变化(Liu,et al,2017);20世纪90年代以来,夏季乌拉尔山地区阻高处于明显偏弱期,鄂霍次克海阻高处于偏强期,但21世纪初,夏季乌拉尔山地区阻高活动频次开始呈上升趋势,而鄂霍次克海阻高频次却呈明显下降趋势(刘刚等,2012)。在对阻高持续活动的大气内动力过程研究的同时,学者们(Wallace,et al,1981;Ding,et al,1987;施能,1996;谭本馗等,2014)还发现了其与气候低频变率的遥相关型联系。李子仪等(2013)研究的个例分析强调苏门答腊西侧的热带印度洋区域存在明显的对流增强,可能激发了类似东亚-太平洋遥相关型的波列向中高纬度传播,从而调制行星波活动再影响乌拉尔山地区位势高度异常的形成和持续。近些年来,有些学者(Luo,et al,2016,2017;Yao,2016;Gong,et al,2017)利用改进的非线性多尺度交互模式,研究了北大西洋急流通过北大西洋涛动正位相影响欧洲阻高的物理机制,乌拉尔山阻高对冬季北极暖—欧亚冷异常的影响,以及北极变暖对冬季乌拉尔山阻高增强的准静止和持续性及欧亚极寒事件的影响,还揭示乌拉尔山阻高对冬季北极海冰减少的促进作用等。以往关于阻高的工作多是基于全年或冬、夏季统计特征或机理的研究(Hartmann,et al,1980;Dole,1982;Shukla,et al,1983;陆日宇等,1998;黄菲等,2002;鹿翠华等,2010),涉及春夏之交的研究较少。5—6月是北半球大气环流春夏转换的关键期,主要特征是北半球西风带行星波发生一次大调整,东亚大槽东移减弱、副热带西风急流突然北跳,这一阶段的东北亚阻高不仅与东亚夏季风的向北推进有密切联系(刘长征等,2004),还与东北冷涡共同对中国天气、气候产生重要影响(孙力,1997;何金海等,2006;Lian,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。对于线性系统X1和X2,从后者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) |
式中,Cij为Xi与Xj的样本协方差,
${{\overset{\bullet }{\mathop{\rm{X}}}\,}_{j,n}} = \frac{{{X_{j,n + k}} - {X_{j,n}}}}{{k\Delta t}}$ | (2) |
式中,
信息流是Liang(2014)在转移熵和Granger因果检验基础上改进的以一种定量方式来表征变量(或事件)之间因果关系的物理量,通过一个变量序列到另一个变量序列的信息事件速率来衡量因果关系,并且这种因果关系是单向的,两个变量之间的信息交换量不仅说明了因果关系的强弱,而且指明了方向。由于北半球阻高的主要活动区都与各自活动区的海陆分布和地形影响有关,很难判断谁对谁的影响更重要,因而本研究引进Liang-Kleeman信息流方法,尝试揭示北半球初夏季节各阻高频繁活动区间的相互联系究竟谁对谁的影响更重要。
3 北半球阻高频繁活动区域的划分及变化特征 3.1 活动频繁区域的划分图1给出了由NCEP/NCAR-Ⅰ再分析资料得到的1979—2018年5月和6月北半球中高纬度地区阻高逐日活动累计天数的分布。不难看出,北半球阻高在5月和6月均存在明显的地域差异:欧洲范围内,以乌拉尔山阻塞为主,纬度偏低;东北亚地区以鄂霍次克海阻塞为主,纬度偏高;阿拉斯加和加拿大西部各有一大值中心,纬度位置与东北亚阻塞同样偏高。5月,自欧洲开始经乌拉尔山—青藏高原北区—鄂霍次克海—北太平洋—北美洲—北大西洋Ⅱ区等阻高活动中心纬度分布,具有“正弦波”特征(北大西洋Ⅰ区阻高除外),即较高纬度—中纬度—高纬度—中纬度的半球尺度分布形状。6月阻高活动中心纬度分布亦具有“正弦波”特征(北大西洋Ⅰ区阻高除外),除欧洲大陆上空阻高较5月有所减少外,其他大部分地区阻高明显增多。
根据阻高出现的地理位置和活动集中程度占北半球活动总天数的百分比,将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高度场上出现了清晰的阻塞环流,文中称为“青藏高原北阻高”。
表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,是(除阿拉斯加阻高外)北半球最高的,这点很值得关注!
区域 | 持续天数(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平均位势高度场进行相关分析(图4、5)发现,各类阻高年活动累计天数与500 hPa高度场均存在相应的显著相关区,其中心区显著性水平均超过95%的信度t检验,这不仅可以进一步验证北半球阻高活动检测和划分结果的科学性,更主要可以揭示阻高各主要活动区引起北半球500 hPa高度场异常的不同特征。
5月,首先,只有乌拉尔山阻高频繁活动能引起整个北半球500 hPa高度场异常,欧洲至亚洲中西部地区和北美洲地区呈现相似的鞍型显著相关场(由北向南为“+”,由西向东为“−”);其次,北大西洋Ⅰ区阻高频繁活动引起北极区与东亚至北大西洋的中高纬度地区为自北向南“+−”异常的半球尺度相关波列;第三,北太平洋阻高频繁活动亦引起北美洲自西向东的纬向“+−”相关波列,以及与其南侧自亚洲西部至太平洋东部的纬向“−”相关显著波状带;其他阻高活动区只产生局地显著相关区。
6月,第一,仍然是乌拉尔山阻高频繁活动能引起整个北半球500 hPa高度场异常,分别引起欧洲、鄂霍次克海至西北太平洋和北美洲中高纬度地区且具有半球尺度的“+”遥相关(图5b);第二,青藏高原北区阻高频繁活动引起斯堪纳维亚半岛附近北冰洋至青藏高原北麓显著的西北—东南向“−+”相关波列,还引起青藏高原下游副热带地区的纬向波状“−”相关带(图5d),而贝加尔湖阻高频繁活动区则相反,形成由地中海—欧洲—贝加尔湖呈西南—东北向“+−+”相关波列(图5c);第三,东北亚阻高频繁活动引起与青藏高原类似的东北亚至中太平洋副热带地区自西北向东南的“+−+”显著相关波列,还引起北美洲与其南面副热带区域的经向较显著的“+−”遥相关波列;其他阻高频繁活动区大都引起局域显著相关区,其中北大西洋Ⅰ和Ⅱ区,以及贝加尔湖阻高频繁活动区所形成的“+”与“−”相间经向相关波列较为引人关注。
3.2 主要变化特征从5—6月北半球阻高天数的标准化时间序列(图6、7)可以看出,各区阻高40 a总体线性趋势均不显著,但阶段性变化明显。从谐波分析(周晓霞等,2007;Li,et al,2004)分离出的年际和年代际分量结果的方差贡献看,各区阻高活动天数的年际分量占显著优势(均在55%以上)。这表明,初夏北半球各区阻高活动天数的年际变化分量起主要作用,预测5—6月北半球阻高活动的关键在年际尺度上。
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月青藏高原北阻高天数的不稳定影响源。
相关系数 | 北美 | 北大西洋Ⅰ | 北大西洋Ⅱ | 欧洲 | 乌拉尔山 | 贝加尔湖 | 青藏高原北 | 东北亚 | 北太平洋 | |
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检验。 |
信息流结果 | 北美 | 北大西洋Ⅰ | 北大西洋Ⅱ | 欧洲 | 乌拉尔山 | 贝加尔湖 | 青藏高原北 | 东北亚 | 北太平洋 | ||
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月东北亚阻高天数的稳定影响源。
相关系数 | 北美 | 北大西洋Ⅰ | 北大西洋Ⅱ | 欧洲 | 乌拉尔山 | 贝加尔湖 | 青藏高原北 | 东北亚 | 北太平洋 | |
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月 | 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月东北亚阻高偏多的形势。
5 2019年6月北半球阻高与极端高温事件2019年6月26—30日,高温热浪席卷欧洲,多地刷新最高温度纪录,还引发多起火灾、加剧空气污染,并造成人员伤亡。从2019年6月北半球地表温度距平分布(图10a)上不难发现,2019年6月欧洲中纬度地区、贝加尔湖—东北亚—北太平洋的高纬度地区、以及北美东部—北大西洋的高纬度地区的地表温度距平值异常偏高。2019年6月500 hPa高度场距平分布(图10b)与1979—2018年6月欧洲阻高活动累计天数和500 hPa高度场的相关系数分布(图5a)相似,中高纬度异常的高压系统稳定维持在欧洲、贝加尔湖—东北亚—北太平洋的高纬度地区、以及北美东部—北大西洋的高纬度地区上空,也就是说上述高温地区均被异常高压系统所占据,持续的下沉气流造成连续多日晴热高温天气。
2019年6月,欧亚中高纬度地区阻高活动频繁,其中欧洲阻高和东亚阻高活动天数异常偏多,分别为19和21 d,较常年平均值分别多12和5 d,反映了阻塞高压作为中高纬度地区典型的持续性环流异常是导致天气、气候异常的重要大尺度环流条件。对2019年6月阻高分布情况(图11)进行统计发现,北半球阻高整体有偏南的趋势,唯独东北亚阻高偏北;具体而言,欧洲阻高异常偏多偏南,且分布集中;东北亚阻高异常偏多偏东偏北,但主要分布在中国东北北部和俄罗斯东北部楚科奇山脉附近,这与以往认知的东北亚阻高有所不同。此外,青藏高原阻高偏少,贝加尔湖阻高异常偏南。以往研究(Black,et al,2004;Stefanon,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,相比之下,最后一次阻高活动强度明显高于前两次。
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 |