气象学报  2018, Vol. 76 Issue (3): 361-378   PDF    
http://dx.doi.org/10.11676/qxxb2018.005
中国气象学会主办。
0

文章信息

罗雄, 李国平. 2018.
LUO Xiong, LI Guoping. 2018.
高空急流对青藏高原切变线影响的数值试验与动力诊断
Numerical experiments and dynamic diagnosis of the upper-level jet stream impact on the Tibetan Plateau shear line
气象学报, 76(3): 361-378.
Acta Meteorologica Sinica, 76(3): 361-378.
http://dx.doi.org/10.11676/qxxb2018.005

文章历史

2017-08-31 收稿
2017-12-08 改回
高空急流对青藏高原切变线影响的数值试验与动力诊断
罗雄1, 李国平1,2     
1. 成都信息工程大学大气科学学院, 成都, 610225;
2. 气象灾害预报预警与评估协同创新中心, 南京, 210044
摘要: 利用NCEP 1°×1° FNL分析资料和中尺度数值模式WRF对一次青藏高原(简称高原)切变线过程进行了数值试验,主要研究高空急流强度对高原切变线的影响,并结合ω方程分析了影响高原切变线上垂直上升运动的若干因子。研究得出高空急流的强度对低层风场有重要影响,急流增强会使高原切变线上的风切变增大,切变线变长,同时高空急流强度的增强也有利于高原切变线上水汽的辐合。高空急流可通过影响高层辐散、低层辐合的散度场垂直配置对高原切变线上的正涡度柱与辐合上升运动产生作用。ω方程的诊断分析表明,温度平流的拉普拉斯项对高原切变线上的垂直上升运动起主导作用,低层暖平流有利于切变线上产生上升运动。高空急流强度的变化对差动涡度平流项的影响要大于温度平流拉普拉斯项,高空急流强度的增强会放大差动涡度平流项和温度平流项的正贡献,从而更加有利于上升运动及高原切变线的维持。
关键词: 高空急流     青藏高原切变线     强度     数值试验     ω方程    
Numerical experiments and dynamic diagnosis of the upper-level jet stream impact on the Tibetan Plateau shear line
LUO Xiong1, LI Guoping1,2     
1. School of Atmospheric Sciences, Chengdu University of Information Technology, Chengdu 610225, China;
2. Collaborative Innovation Center on Forecast and Evaluation of Meteorological Disasters, Nanjing 210044, China
Abstract: Using the NCEP 1°×1° FNL analysis data and the mesoscale model WRF, numerical experiments have been conducted to simulate a Tibetan Plateau shear line process. The study focuses on the impact of the upper-level jet stream intensity on the Tibetan Plateau shear line. In combination with the ω equation, factors affecting the ascending motion in the Tibetan Plateau shear line are also analyzed. The results show that the intensity of the upper-level jet stream has an important influence on low-level wind fields. The increases in jet stream will intensify the wind shear along the Tibetan Plateau shear line, and the shear line will become longer. Meanwhile, the enhancement of the upper-level jet intensity is also favorable for the convergence of water vapor along the Tibetan Plateau shear line. The upper jet can affect positive vorticity, convergence and ascending motion along the Tibetan Plateau shear line by influencing the vertical configuration of lower-level convergence and upper-level divergence. The diagnostic analysis of the ω equation shows that the Laplace term of temperature advection plays a dominant role in the ascending motion along the Tibetan Plateau shear line, and the low-level warm advection is favorable for the development of ascending motion along the shear line. The influence of the variation of upper-level jet stream intensity on the differential term of the vorticity advection is greater than that of the Laplace term of temperature advection, and the enhancement of the upper-level jet intensity will enlarge the positive contributions of the differential vorticity advection term and the temperature advection term, which is more favorable for the development of ascending motion and the maintenance of the Tibetan Plateau shear line.
Key words: Upper-level jet stream     Tibetan Plateau shear line     Intensity     Numerical experiments     ω equation    
1 引言

青藏高原(简称高原)切变线通常指500 hPa等压面上温度梯度小,3个测站风向对吹或2个测站风向对吹且长度大于5个经(纬)距的辐合线(李山山等,2017a),是在特殊地形条件下形成的一类高原天气系统,其活动不仅会影响高原地区,而且在东移、发展过程中会对高原以东地区的天气产生重要影响(郁淑华等, 1997, 2007徐建芬等,2000高文良等,2007屠妮妮等,2010李山山等,2017b)。关于高空急流的研究,Tao等(1987)指出,200 hPa高空西风急流是影响东亚天气、气候的重要系统。高空急流及其附近的散度场和与高空急流相伴随的次级环流对上升运动、低层气旋、暴雨以及强对流等均有重要影响(Brill,et al,1985Prezerakos,et al,2006Whitney,2009廖移山等,2011程正泉等,2012)。林毅等(2006)通过模拟高空急流在一次大暴雨过程中的作用,指出高、低空急流在适宜的配置下,高空急流南侧的辐散场和低空急流右侧的辐合场耦合形成强烈的上升运动,并且高层的强辐散场有利于低层的辐合,有助于低层切变线的维持。周兵等(2003)开展了高空西风急流强度对长江中下游暴雨影响的数值试验,研究表明,高空急流的强度减弱对低空急流的变化也有重要影响,高空急流的减弱有利于低层偏西气流的减弱。师锐等(2011)指出,在对流层200 hPa层上,相对较强的高空急流有利于高原切变线的东移,无论是在移出高原切变线的生成时,还是在移出高原主体时,高原切变线都位于高空急流入口区的南侧,而急流南侧所产生的高空辐散有利于500 hPa高原切变线的发展和移出高原。由此可见,高空急流对低层系统有着重要的影响。

目前关于高原切变线及其伴随的高原涡、西南涡等高原低值系统相互作用的研究已取得一定成果(姚秀萍等,2014),但对高空急流和高原切变线这两个高、低层系统关系的研究还不多见。所以,本研究以高原切变线为主要对象,对一次发生在高空急流入口区南侧的高原切变线个例进行数值模拟,然后开展高空西风急流强度变化对高原切变线影响的数值试验,通过对比分析敏感性试验与对照试验的差异,并结合ω方程的诊断计算,以探寻高空西风急流对低层高原切变线的影响和作用。

2 资料与天气个例 2.1 资料

所用资料有NCEP/FNL(National Centers for Environment Prediction Final Analyses)分析资料(逐6 h,水平分辨率1°×1°),中国自动气象站与CMORPH(Climate Prediction Center Morphing Technique)融合的水平分辨率为0.1°×0.1°逐时降水数据。

2.2 天气个例

2014年6月29日至7月1日出现了一次高原切变线天气过程,从500 hPa天气图(图略)可知,此次高原切变线过程发生于稳定的大尺度环流背景下,欧亚大陆高纬度地区呈“两槽一脊”的环流形势,中亚为脊区,乌拉尔山以东的西西伯利亚和亚洲东岸为低槽,西太平洋副热带高压(副高)588 dagpm线逐步西伸,呈东西向的带状分布;200 hPa上高空西风急流稳定维持,并且此次高原切变线位于急流入口区南侧。6月29日00时(世界时,下同)在高原东南部玛多至索县有一横切变线生成,随后切变线向东北方向移动,在移动过程中切变线南侧风速先增大后减小,北侧风速逐渐增大。30日00时高原切变线北侧风速增大,切变线开始向西南方向移动;7月1日00时副高的西脊点达到105°E,副高脊线位于25°N附近,西伸的副高与高原上的高压逐渐合并,使得高原切变线趋于减弱消失。受此次高原切变线过程的影响,西藏、四川部分地区发生大于25 mm的降水,降水中心位于西藏林芝,48 h累计降水量达75.5 mm。

3 数值模拟、结果验证及试验方案 3.1 模式方案

利用WRF模式V3.7.1对上述高原切变线过程进行数值模拟。模拟区域为两重双向嵌套的网格,采用兰勃特投影,模拟区域中心为(31°N,95°E),两个模拟区域的格距分别为30、10 km,格点数分别为221×152、313×196,模式垂直方向为不等距28层,顶层气压为50 hPa,初始场和侧边界条件均采用NCEP/FNL 1°×1°逐6 h的分析资料。两重网格均采用WSM6微物理方案、Noah陆面方案、YSU边界层方案、KF积云对流参数化方案、RRTM长波辐射方案和Dudhia短波辐射方案。模拟时间步长为180 s,逐时输出一次数据,积分时间从2014年6月29日00时到7月1日06时,共54 h,涵盖了此次高原切变线生成、发展、减弱和消失的全过程。

3.2 模拟结果验证

图 1为模式积分18、24、36、54 h与FNL分析资料的500 hPa风场的对比,可知模拟的风场与FNL分析资料基本一致,高原切变线的位置与FNL分析资料也对应较好,只是风速偏大。需要指出的是,模式积分18和54 h切变线西段的“对吹”特征与FNL分析资料相比不明显,其表现为风向的切变,但是整体上看模式较成功地模拟出了此次切变线过程。由图 2可知,模式积分18(图 2ab)和24 h(图 2cd),模拟的200 hPa风场与FNL分析资料基本吻合,但是急流中心的风速偏大,随着积分时间的延长(图 2e-h),模拟的200 hPa风场与FNL分析资料基本一致,并且高空急流位置和强度对应也较好。

图 1 FNL分析资料(a、c、e、g)与模拟(b、d、f、h)的500 hPa风场 (a、b.6月29日18时,c、d. 30日00时,e、f.30日12时,g、h.7月1日06时;其中,黑色实线为高原切变线) Figure 1 Wind fields at 500 hPa from the FNL analysis data (a, c, e, g) and simulation (b, d, f, h) (a, b. 18:00 UTC 29 June; c, d. 00:00 UTC 30 June; e, f. 12:00 UTC 30 June; g, h. 06:00 UTC 1 July; The black solid line represents the Tibetan Plateau shear line)
图 2 FNL分析资料(a、c、e、g)与模拟(b、d、f、h)的200 hPa风场(矢线)和高空急流(色阶,单位:m/s) (a、b.6月29日18时,c、d.30日00时,e、f.30日12时,g、h.7月1日06时) Figure 2 Wind fields (vectors) and upper level jet stream (shaded, unit: m/s) at 200 hPa from the FNL analysis data (a, c, e, g) and simulation (b, d, f, h) (a, b. 18:00 UTC 29 June; c, d. 00:00 UTC 30 June; e, f. 12:00 UTC 30 June; g, h. 06:00 UTC 1 July)

通过对比模式细网格和CMORPH融合降水资料的24 h累计降水量可知,6月29日(图 3a)在高原中东部降雨带基本呈东西走向,降水带的中心雨量为48.68 mm,在林芝附近存在一个降水的大值中心44.81 mm。从雨量模拟结果可见(图 3c),在高原中东部也有一条准东西走向的雨带,中心雨量为41.45 mm,比融合降水资料略小;而与林芝地区相对应的高原南部的中心降水值为52.57 mm,比融合降水量偏大。6月30日CMORPH融合降水中心在(30°N,94.3°E),中心降水值为49.89 mm(图 3b),模拟结果的中心位置在(30.1°N,94.32°E),与融合降水中心的位置基本一致,中心雨量为67.05 mm(图 3d),模拟结果偏大。总体而言,WRF模拟的降水范围和强度与CMORPH融合降水量基本一致;但积分24 h后模拟的降水量略偏大。江志红等(2013)研究表明,中国区域平均的融合降水资料的时、空精度远高于CMORPH和TRMM卫星降水资料,但是在青藏高原地区相对较差,这是由于西部地区观测资料缺测较多,难以捕捉到小尺度强降水过程的结果,并且网格内站点稀疏,导致该地区融合降水资料主要取决于卫星降水资料。这可能是降水模拟结果比融合资料偏大的一个原因。此外,模拟的风速偏大也有可能是引起降水偏大的一个原因,考虑到高原降雨的复杂性,可认为WRF模式较为成功地模拟了此次高原切变线引起的降雨过程。总体来看,WRF模式较为成功地再现了此次高原切变线过程,可利用开展下一步敏感性试验。

图 3 2014年6月29-30日CMORPH融合降水(a、b)与模拟(c、d)的高原区域降水量分布(单位:mm) (a、c.6月29日00时-30日00时,b、d.6月30日00时-7月1日00时) Figure 3 Distributions of precipitation in the Tibetan Plateau region based on the CMORPH (a, b) and simulation (c, d) from 29 June to 30 June in 2014 (unit: mm) (a, c. from 00:00 UTC 29 June to 00:00 UTC 30 June; b, d. from 00:00 UTC 30 June to 00:00 UTC 1 July)
3.3 试验方案

为了进一步探讨高空急流的强度变化对高原切变线的影响,以对照试验(Ctrl)为基础设计了高空急流增强试验(S-Jet)和高空急流减弱试验(W-Jet),具体方案为:将粗、细网格中模式大气上层(19-23 eta层,相当于265-143 hPa)高空急流所在区域(32°-40°N,90°-112°E)的风速强迫增大(减小)0.5倍,但不改变风向,其余同对照试验。需要指出的是,在大尺度条件下,大气风场是向气压场适应的,若试验中只改变模式风场的初值,模式大气积分会较快失去对初始风场扰动的记忆,高空急流强度只在模式积分开始后较短时间内有变化,并且在试验过程中随着模式积分时间的延长,由于模式的调整和适应作用,风场会重新调整到与对照试验几乎相同的状态。因此,为了保证试验结果的合理性,先修改前处理系统中由metgrid程序得到的各时刻插值数据的风场,然后再执行模式的主程序过程,这样不仅使风场具有初始扰动,而且又能通过更新侧边界条件使模式上层的风场每隔6 h强迫增强或减弱。

4 数值试验结果分析 4.1 高空急流强度对500 hPa风场的影响

高守亭等(1984)认为,叠加在由高、低空质量调整引起的地转风上的变压风能引起低空急流的发展。斯公望(1989)指出,高层急流入口右侧的辐散引起的变压风与低空急流的形成有关。由此可见,高空急流和低层风场存在一定联系,结合图 45可知,高空急流强度的改变会直接导致高、低空的风场发生改变。模式积分18 h(图 4a),S-Jet试验的200 hPa急流中心的风速以增强为主,并且最大增强了10 m/s;在W-Jet试验中(图 4b),高空急流中心风速最大减弱了12 m/s,除了试验区域风速有变化外,其他区域的风速也会有一定变化,并且风速的变化方向不是单一的。模式积分54 h后(图 4cd),200 hPa的风场也发生了类似的改变,但是随着积分时间的延长,由于模式的调整和适应作用,高空急流区风速的变化幅度有所减小。林毅等(2006)指出,高空急流在一定程度上制约着低层切变线的位置和辐合强度。从6月29日18时各组试验500 hPa的风场(图 5a-c)可以看出,高原切变线的位置和走向基本一致,但其长度和两侧的风速有明显差异,增强高空急流后,切变线西段的风切变区相对对照试验向南略有偏移,使得切变线长度变长,并且切变线南侧的西南风也有一定增强;而W-Jet试验与对照试验相比,切变线东段的风切变明显减弱,导致切变线的长度变短。在切变线减弱阶段(图 5d-f),S-Jet试验中切变线的长度也要长于对照试验,虽然在W-Jet试验中切变线长度与对照试验相差不大,但是风切变和西南风均明显减弱。

图 4 各试验与对照试验200 hPa的风速差值分布(单位:m/s) (a、b.6月29日18时,c、d.7月1日06时) Figure 4 Distributions of wind speed differences at 200 hPa between control experiment and other experiments (unit: m/s) (a, b. 18:00 UTC 29 June; c, d. 06:00 UTC 1 July)
图 5 各试验500 hPa风场 (a、b、c.6月29日18时,d、e、f.7月1日06时;黑色实线为高原切变线) Figure 5 Wind fields at 500 hPa from all experiments (a, b, c. 18:00 UTC 29 June; d, e, f. 06:00 UTC 1 July; The black solid line represents the Tibetan Plateau shear line)
4.2 高空急流强度对切变线上水汽特征的影响

郁淑华等(2003)研究表明,切变线水汽带的变化对高原切变线活动有指示意义。切变线水汽带范围变宽时,高原切变线稳定维持;切变线水汽带减弱时,则高原切变线会向东南移。500 hPa切变线是水汽聚集带,切变线附近上空的水汽和不稳定能量聚集是切变线引发强降水的重要机制(何光碧等,2014);可见水汽辐合是高原切变线的一个重要特征。从500 hPa的水汽通量演变(图略)可知,维持此次高原切变线的水汽来源主要有两支,一支是来自孟加拉湾的西南水汽输送,另一支是来自中国南海沿副高西侧的东南水汽输送,这两支水汽在高原东侧汇合后形成较强的水汽区。结合图 6可见,高空急流的强度变化对低层水汽含量有一定影响,急流强度增强,低层的水汽含量也相应增加,这可能是由于切变线南侧的西南风增强从而更有利于水汽的输送。图 7为各试验高原切变线上水汽通量散度的垂直剖面。由图可见,6月29日18时(图 7a-c),从地面到400 hPa左右沿切变线有明显的水汽辐合,其中S-Jet试验切变线上水汽通量散度中心最强为-1.4×10-6 g/(s·cm2· hPa),其强度和范围要大于对照试验。而在W-Jet试验中,水汽通量散度中心最强为-0.2×10-6 g/(s·cm2·hPa),相比对照试验明显减小。在切变线减弱阶段,水汽的辐合明显减弱,从图 7d-f可以看到,S-Jet试验中在对流层低层高原切变线东段水汽辐合的强度要明显强于对照试验,而W-Jet试验中,高原切变线上水汽的辐合相对较弱。

图 6 在高原切变线活动范围内(29°-32°N,92°-104.5°E)各试验500 hPa比湿(单位:g/kg)区域平均的时间廓线 Figure 6 Temporal evolutions of regionally averaged specific humidity (unit: g/kg) at 500 hPa from all experiments over the area (29°-32°N, 92°- 104.5°E) covered by the Tibetan Plateau shear line
图 7 各试验沿高原切变线位置(30°-32°N)平均的水汽通量散度(单位:10-6 g/(s·cm2·hPa))的纬向-垂直剖面 (a、b、c.6月29日18时,d、e、f.7月1日06时) Figure 7 Zonal-vertical cross sections of water vapor flux divergence (unit: 10-6 g/(s·cm2·hPa) from all experiments averaged between 30° and 32°N along the Tibetan Plateau shear line (a, b, c. 18:00 UTC 29 June; d, e, f. 06:00 UTC 1 July)

上述分析表明,高空急流强度增强有利于高原切变线上水汽的辐合,反之急流减弱会使切变线上的水汽辐合减弱。李山山等(2017b)研究指出,水平风辐合带也是水汽辐合带,切变线上风场的分布对水汽的辐合、辐散以及切变线的维持有重要作用。因此,高空急流强度的变化会引起低层高原切变线风场发生相应变化,从而改变高原切变线附近的水汽特征。

4.3 高空急流强度对切变线动力结构的影响

何光碧等(2011, 2014)指出,平均的500 hPa横切变线上为正涡度与辐合上升运动;在切变线活动期间有正涡度、辐合上升运动与切变线配合。由图 8a-c可见,6月29日18时,在高空急流右侧(此次高原切变线南北活动的主要范围为29°-32°N)为明显的低层辐合、高层辐散。在S-Jet试验中,低层辐合中心强度达-7×10-5 s-1,对照试验中低层辐合中心强度为-5×10-5 s-1;而在W-Jet试验中,低层辐合中心强度减小为-4×10-5 s-1。从涡度的经向垂直剖面(图 9a-c)可知,在辐合区有正涡度柱与之对应,最强可伸展到300 hPa以上。并且当高空急流增强时,正涡度中心强度也相应增大,反之则会减弱。结合图 10a-c可以看出,高空急流的右侧有明显的次级环流圈,下沉支较弱,上升气流较强并与切变线上的辐合上升区重叠,造成在切变线上出现强的上升运动。在S-Jet试验中,上升区的垂直速度最大可达-0.6 Pa/s,强于对照试验的-0.4 Pa/s;而在W-Jet试验中,上升区中心的垂直速度减小为-0.2 Pa/s。由于风切变增强有利于垂直速度的增大和维持(郑淋淋等,2016),而前述的研究结果表明高空急流增强会使切变线上的风切变增强,这可能是切变线上垂直上升运动增强的一个原因。模式积分54 h后,即在高原切变线的减弱消亡阶段,对照试验中低层的辐合中心(图 8e)、正涡度柱(图 9e)以及所对应的上升运动(图 10e)均明显减弱,并且从S-Jet试验(图 8d9d10d)和W-Jet试验(图 8f9f10f)的结果可见,随着高空急流强度的改变,切变线上动力特征量的强度也随之改变。

图 8 各试验沿高原切变线位置(90°-105°E)平均的散度(单位:10-5 s-1)经向-垂直剖面 (a、b、c.6月29日18时,d、e、f.7月1日06时) Figure 8 Meridional-vertical cross sections of divergence (unit: 10-5 s-1) from all experiments averaged between 90° and 105°E along the Tibetan Plateau shear line (a, b, c. 18:00 UTC 29 June; d, e, f. 06:00 UTC 1 July)
图 9 各试验沿高原切变线位置(90°-105°E)平均的涡度(单位:10-5 s-1)经向-垂直剖面 (a、b、c.6月29日18时,d、e、f.7月1日06时) Figure 9 Meridional-vertical cross sections of vorticity from all experiments averaged between 90° and 105°E along the Tibetan Plateau shear line (unit: 10-5 s-1) (a, b, c. 18:00 UTC 29 June; d, e, f. 06:00 UTC 1 July)
图 10 各试验沿高原切变线位置(90°-105°E)平均的垂直速度经向-垂直剖面 (a、b、c.6月29日18时,d、e、f.7月1日06时;矢量表示v×(-100ω), 水平风速单位为m/s, 垂直速度ω单位为Pa/s) Figure 10 Meridional-vertical cross sections of vertical velocity from all experiments averaged between 90° and 105°E along the Tibetan Plateau shear line (a, b, c. 18:00 UTC 29 June; d, e, f. 06:00 UTC 1 July; vectors represent v×(-100ω), the unit of v is m/s, and the unit of ω is Pa/s)
5 基于ω方程的动力诊断

为了进一步探讨影响前述垂直上升运动的因子及贡献,采用不考虑非绝热加热项的准地转ω方程(朱乾根等,2000李国平,2014)

(1)

式中,σ为层结稳定参数,f为地转参数,Vg为地转风,ζg为相对涡度,ϕ为位势高度。方程右边第1项为差动涡度平流项(简称涡度平流项),第2项为温度平流的拉普拉斯项(简称温度平流项)。当涡度平流随气压降低(高度升高)而增强,或者在暖平流的中心区,均有利于形成上升运动。除了上述简单的定性分析以外,还可以通过对比式(1)涡度平流强迫项和温度平流强迫项的大小,分析两者对上升运动的贡献,从而分析不同影响因子在天气过程中的作用,揭示天气过程的发生、发展机制。

图 1112分别给出了涡度平流项和温度平流项的经向-垂直剖面,6月29日18时,从对照试验中(图 11b)可见,在上升运动的主要区域(29°-32°N)由地面至200 hPa涡度平流项基本上随高度增强,但是在32°N左侧从地面至400 hPa温度平流项却随高度减弱。S-Jet试验中上升区的涡度平流项随高度的变化和对照试验基本一致,在30°N附近从低层至高层负涡度平流项转换为正涡度平流项,有利于上升运动的产生。在W-Jet试验中29°-32°N涡度平流项随高度增强,但低层的正涡度平流项却明显小于对照试验。结合图 12a-c可以看出,各组试验在上升运动区伴有暖平流,并且随着高空急流强度的变化,低层暖平流中心的范围和强度也相应改变。在S-Jet试验中暖平流中心值可达12×10-5 K/s,虽然对照试验和W-Jet试验暖平流中心的最大值均为9×10-5 K/s,但对照试验中暖平流中心的范围明显大于W-Jet试验。温度平流从高原地面向上随高度减弱反映了温度平流在垂直方向上的变化,而这种“上冷下暖”的温度平流垂直差异会导致气温直减率增大,使大气层结趋于不稳定(朱禾等,2016)。积分54 h后,由于高原切变线进入减弱消亡阶段,上升运动开始减弱,从图 11d-f可见,各组试验涡度平流项随高度增强的趋势变得不明显,甚至开始出现随高度减弱,但暖平流项仍然存在(图 12d-f)。

图 11 各试验沿高原切变线位置(90°-105°E)平均的涡度平流(单位:10-9 s-2)的经向-垂直剖面 (a、b、c.6月29日18时,d、e、f.7月1日06时) Figure 11 Meridional-vertical cross sections of vorticity advection (unit: 10-9 s-2) from all experiments averaged between 90° and 105°E along the Tibetan Plateau shear line (a, b, c. 18:00 UTC 29 June; d, e, f. 06:00 UTC 1 July)
图 12 各试验沿高原切变线位置(90°-105°E)平均的温度平流(单位:10-5 K/s)的经向-垂直剖面 (a、b、c.6月29日18时,d、e、f.7月1日06时) Figure 12 Meridional-vertical cross sections of temperature advection (unit: 10-5 K/s) from all experiments averaged between 90° and 105°E along the Tibetan Plateau shear line (a, b, c. 18:00 UTC 29 June; d, e, f. 06:00 UTC 1 July)

为了进一步了解二者对于上升运动的定量贡献,图 13给出了此次高原切变线活动范围内(29.5°-32°N,92°-104.5°E)各组试验ω方程涡度平流强迫项和温度平流强迫项区域平均的垂直廓线。从图 13b可见,6月29日18时对照试验中500-250 hPa温度平流项廓线与两项之和的廓线基本重合,最大值出现在400 hPa(15.2×10-15 hPa-1·s-3),而涡度平流项基本处于0线附近,说明温度平流项在上升运动中起主导作用。在S-Jet试验中(图 13a),400 hPa温度平流项的值为20.0×10-15 hPa-1·s-3,其贡献比对照试验明显增大,并且500-250 hPa附近的涡度平流项廓线基本上为正值,由此可见涡度平流项的贡献要强于对照试验。在W-Jet试验中(图 13c),500-400 hPa温度平流项要强于对照试验,而此时涡度平流项却为负值,但由于温度平流项的正值较大,所以二者之和仍为正值,仍然会有上升运动的产生;400-250 hPa各项垂直廓线与对照试验基本一致,只是两项之和要比对照试验小。

图 13 各试验高原切变线活动范围内(29°-32°N,92°-104.5°E) ω方程各项区域平均的垂直廓线 (a、b、c.6月29日18时,d、e、f.7月1日06时;单位:10-15 hPa-1·s-3;点线为涡度平流项,虚线为温度平流项,实线为两项之和) Figure 13 Vertical profiles of the terms in the ω equation for all experiments regionally averaged over the area (29°-32°N, 92°-104.5°E) covered by the Tibetan Plateau shear line (a, b, c. 18:00 UTC 29 June; d, e, f. 06:00 UTC 1 July; The unit is 10-15 hPa-1·s-3, the dotted line represents the term of vorticity advection, the dashed line represents the term of temperature advection, and the solid line represents the sum of the two components)

7月1日06时,高原切变线处于减弱阶段,从图 13d可知,对照试验中上升运动依然是温度平流项起主导作用,500-400 hPa涡度平流项与温度平流项之和变为负值,说明此时不利于低层的上升运动。在S-Jet试验中(图 13e),500-400 hPa两项之和的负值减小;而在W-Jet试验中(图 13f),两项之和的负值增大,说明高空急流强度的增强可在一定程度上减缓高原切变线减弱阶段垂直上升运动的衰减。

综上所述,ω方程中温度平流项对上升运动的贡献最大,暖平流中心有利于产生上升运动;高空急流强度的变化对涡度平流项的影响要大于温度平流项,高空急流强度增强会放大涡度平流项和温度平流项的贡献,有利于低层高原切变线对应的上升运动的加强、维持(图 14)。

图 14 高空急流强度对高原切变线影响机理示意 Figure 14 The schematic diagram of the impact of the upper-level jet stream intensity on the Tibetan Plateau shear line
6 结论与讨论

利用WRF模式V3.7.1对一次发生在高空急流入口区右侧的高原横切变线过程进行了数值模拟,在模拟成功的基础上通过敏感性试验和ω方程的诊断计算研究了高空急流强度对高原切变线的影响,得到以下结果:

(1) 高空急流的强度对低层风场有重要影响,急流加强会使高原切变线上的风切变增强,切变线长度延长。

(2) 高空急流强度的增强有利于高原切变线上水汽的辐合,即高空急流可通过改变高原切变线周边水汽的分布对切变线产生影响。

(3) 高空急流可通过影响高层辐散、低层辐合的散度场垂直配置,对高原切变线上的正涡度柱与辐合上升运动产生作用。

(4) ω方程的动力诊断表明温度平流项对高原切变线上的垂直上升运动起主导作用,暖平流中心区有利于产生上升运动。高空急流强度增强会放大涡度平流项和温度平流项的贡献,从而更有利于上升运动的加强及高原切变线的维持。

本研究通过数值试验研究了高空急流强度变化对高原切变线的影响,这一工作对于认识二者的关系具有一定的意义,但所做研究仍是初步的,需要更多个例的数值试验来支持文中所揭示的现象及影响机理。此外,非绝热加热对于高原切变线的影响也不可忽视,由于高原上地面感热通量和潜热通量以及降水过程中释放的凝结潜热对于高原切变线的影响较为复杂,加之篇幅有限,文中对此暂不做讨论;还有诸如高空急流强度实际的气候变化差异与数值试验设计的差异对高原切变线的影响有何异同,以及高、低空急流耦合对高原切变线的作用等深层次问题,都值得在后续的工作中加以考虑和研究。

参考文献
程正泉, 陈联寿, 李英. 2012. 登陆热带气旋海马(0421)变性加强的诊断研究. 气象学报, 70(4): 628–641. Cheng Z Q, Chen L S, Li Y. 2012. Diagnosis of the extra-tropical transition of landfalling tropical cyclone Haima (2004). Acta Meteor Sinica, 70(4): 628–641. DOI:10.11676/qxxb2012.051 (in Chinese)
高守亭, 孙淑清. 1984. 次天气尺度低空急流的形成. 大气科学, 8(2): 178–188. Gao S T, Sun S Q. 1984. The forming of subsynoptic scale low-level jet stream. Chinese J Atmos Sci, 8(2): 178–188. (in Chinese)
高文良, 郁淑华. 2007. 高原低涡东移出高原的平均环流场分析. 高原气象, 26(1): 206–212. Gao W L, Yu S H. 2007. Analyses on mean circulation field of the plateau low vortex moving out of Tibetan Plateau. Plateau Meteor, 26(1): 206–212. (in Chinese)
何光碧, 师锐. 2011. 夏季青藏高原不同类型切变线的动力、热力特征分析. 高原气象, 30(3): 568–575. He G B, Shi R. 2011. Studies on dynamic and thermal characteristics of different shear lines over Tibetan Plateau in summer. Plateau Meteor, 30(3): 568–575. (in Chinese)
何光碧, 师锐. 2014. 三次高原切变线过程演变特征及其对降水的影响. 高原气象, 33(3): 615–625. He G B, Shi R. 2014. Analysis on evolution characteristics of three plateau shear lines and their effect on precipitation. Plateau Meteor, 33(3): 615–625. DOI:10.7522/j.issn.1000-0534.2013.00023 (in Chinese)
江志红, 卢尧, 丁裕国. 2013. 基于时空结构指标的中国融合降水资料质量评估. 气象学报, 71(5): 891–900. Jiang Z H, Lu Y, Ding Y G. 2013. Analysis of the high-resolution merged precipitation products over China based on the temporal and spatial structure score indices. Acta Meteor Sinica, 71(5): 891–900. DOI:10.11676/qxxb2013.076 (in Chinese)
李国平. 2014. 新编动力气象学. 2版. 北京: 气象出版社: 131-132. Li G P. 2014. The Dynamic Meteorology Newly Organized. 2nd ed. Beijing: China Meteorological Press: 131-132.
李山山, 李国平. 2017a. 一次高原低涡与高原切变线演变过程与机理分析. 大气科学, 41(4): 713–726. Li S S, Li G P. 2017a. Evolution and mechanism analysis of a plateau vortex and plateau shear line. Chinese J Atmos Sci, 41(4): 713–726. (in Chinese)
李山山, 李国平. 2017b. 一次鞍型场环流背景下高原东部切变线降水的湿Q矢量诊断分析. 高原气象, 36(2): 317–329. Li S S, Li G P. 2017b. Diagnostic analysis based on wet Q-vector of a shear line with rain on the east side of Qinghai-Tibetan Plateau under the saddle pattern circulation background Field. Plateau Meteor, 36(2): 317–329. (in Chinese)
廖移山, 冯新, 石燕, 等. 2011. 2008年"7.22"襄樊特大暴雨的天气学机理分析及地形的影响. 气象学报, 69(6): 945–955. Liao Y S, Feng X, Shi Y, et al. 2011. Analysis of the mechanism for "2008.7.22" excessive rain event in Xiangfan with a focus on the terrain effect. Acta Meteor Sinica, 69(6): 945–955. DOI:10.11676/qxxb2011.082 (in Chinese)
林毅, 刘铭, 刘爱鸣. 2006. 高低空急流在闽西北大暴雨过程中的作用及数值模拟. 气象科学, 26(4): 449–455. Lin Y, Liu M, Liu A M. 2006. Numreical study and action of high and low level jet stream in "2002.6" torrential rainstorm of northwest Fujian. Scientia Meteor Sinica, 26(4): 449–455. (in Chinese)
师锐, 何光碧. 2011. 移出与未移出高原的高原切变线背景环流对比分析. 高原气象, 30(6): 1453–1461. Shi R, He G B. 2011. Contrast analysis on background circulation of plateau shear line moving out and not moving out of the Tibetan Plateau. Plateau Meteor, 30(6): 1453–1461. (in Chinese)
斯公望. 1989. 论东亚梅雨锋的大尺度环流及其次天气尺度扰动. 气象学报, 47(3): 312–323. Si G W. 1989. On the large-scale circulation of Mei-Yu system over East Asia. Acta Meteor Sinica, 47(3): 312–323. DOI:10.11676/qxxb1989.041 (in Chinese)
屠妮妮, 何光碧. 2010. 两次高原切变线诱发低涡活动的个例分析. 高原气象, 29(1): 90–98. Tu N N, He G B. 2010. Case analysis on two low vortexes induced by Tibetan Plateau shear line. Plateau Meteor, 29(1): 90–98. (in Chinese)
徐建芬, 陶健红, 夏建平. 2000. 青藏高原切变线暴雪中尺度分析及其涡源研究. 高原气象, 19(2): 187–197. Xu J F, Tao J H, Xia J P. 2000. A meso-scale analysis of Qinghai-Xizang Plateau snow storm and its vortex source study. Plateau Meteor, 19(2): 187–197. (in Chinese)
姚秀萍, 孙建元, 康岚, 等. 2014. 高原切变线研究的若干进展. 高原气象, 33(1): 294–300. Yao X P, Sun J Y, Kang L, et al. 2014. Advances on research of shear convergence line over Qinghai-Xizang Plateau. Plateau Meteor, 33(1): 294–300. DOI:10.7522/j.issn.1000-0534.2013.00164 (in Chinese)
郁淑华, 高文良, 顾清源. 2007. 近年来影响我国东部洪涝的高原东移涡环流场特征分析. 高原气象, 26(3): 466–475. Yu S H, Gao W L, Gu Q Y. 2007. The middle-upper circulation analyses of the plateau vortex moving out of Plateau and influencing flood in east China in recent years. Plateau Meteor, 26(3): 466–475. (in Chinese)
郁淑华, 何光碧. 2003. 水汽图像在高原天气预报中应用的初步分析. 高原气象, 22(S1): 75–82. Yu S H, He G B. 2003. Analysis on application of statellite vapour image to weather prediction over Plateau. Plateau Meteor, 22(S1): 75–82. (in Chinese)
郁淑华, 何光碧, 滕家谟. 1997. 青藏高原切变线对四川盆地西部突发性暴雨影响的数值试验. 高原气象, 16(3): 306–311. Yu S H, He G B, Teng J M. 1997. The numerical experiments of influence of shear line over Qinhai-Xizang Plateau on a sharp heavy rain in west Sichuan basin. Plateau Meteor, 16(3): 306–311. (in Chinese)
郑淋淋, 孙建华. 2016. 风切变对中尺度对流系统强度和组织结构影响的数值试验. 大气科学, 40(2): 324–340. Zheng L L, Sun J H. 2016. The impact of vertical wind shear on the intensity and organizational mode of mesoscale convective systems using numerical experiments. Chinese J Atmos Sci, 40(2): 324–340. (in Chinese)
周兵, 韩桂荣, 何金海. 2003. 高空西风急流对长江中下游暴雨影响的数值试验. 南京气象学院学报, 26(5): 595–604. Zhou B, Han G R, He J H. 2003. Numerical experiment of effects of the upper level westerly jet on a torrential rain over the middle-lower reaches of the Yangtze river. J Nanjing Inst Meteor, 26(5): 595–604. (in Chinese)
朱禾, 王秀明. 2016. 影响气温直减率变化的物理过程讨论. 气象科技进展, 6(5): 51–54. Zhu H, Wang X M. 2016. A discussion about the physical process of environmental temperature lapse rate tendency. Adv Meteor Sci Technol, 6(5): 51–54. (in Chinese)
朱乾根, 林锦瑞, 寿绍文, 等. 2000. 天气学原理和方法. 3版. 北京: 气象出版社: 120-122. Zhu Q G, Lin J R, Shou S W, et al. 2000. The Principle and Method of Weather. 3rd ed. Beijing: China Meteorological Press: 120-122.
Brill K F, Uccellini L W, Burkhart R P, et al. 1985. Numerical simulations of a transverse indirect circulation and low-level jet in the exit region of an upper-level jet. J Atmos Sci, 42(12): 1306–1320. DOI:10.1175/1520-0469(1985)042<1306:NSOATI>2.0.CO;2
Prezerakos N G, Flocas H A, Brikas D. 2006. The role of the interaction between polar and subtropical jet in a case of depression rejuvenation over the eastern Mediterranean. Meteor Atmos Phys, 92(1-2): 139–151. DOI:10.1007/s00703-005-0142-y
Tao S Y, Chen L. 1987. A review of recent research on the East Asian summer monsoon in China//Chang C P, Krishnamurti T N. Monsoon Meteorology. Oxford: Oxford University Press, 60-92
Whitney L F Jr. 2009. Relationship of the subtropical jet stream to severe local storms. Mon Wea Rev, 105(4): 398–412.