期刊检索:
  暴雨灾害   2021, Vol. 40 Issue (5): 474-483.  DOI: 10.3969/j.issn.1004-9045.2021.05.004

论文

DOI

10.3969/j.issn.1004-9045.2021.05.004

资助项目

江苏省气象局面上项目(KM201910)

第一作者

乔娜, 主要从事中小尺度天气与数值模拟研究。E-mail: qixiang0820@163.com.

文章历史

收稿日期:2020-02-27
定稿日期:2020-10-27
东北冷涡背景下中尺度低涡的演变成因及其对MCS的影响分析
乔娜 , 钱鹏 , 周勍 , 张晨昕 , 吴昕悦     
江苏省镇江市气象局, 镇江 212000
摘要:利用WRF中尺度数值模式,NCEP/NCAR分析资料、多普勒雷达观测资料等,对2015年8月3日发生在山东地区附近的一次东北冷涡背景下中尺度低涡伴随MCS的过程进行数值模拟和分析,研究了东北冷涡影响下中尺度低涡发生发展的原因及其对MCS发展的影响。结果表明:(1)在天气尺度的东北冷涡槽前,高层高空急流入口区右侧存在强辐散场,低层有不稳定能量释放为中尺度低涡以及对流发生提供有利背景场,低层中尺度低涡的生成促进了MCS的形成;(2)中尺度低涡形成和发展时期,低层正涡度主要受到水平散度项和垂直输送项的影响,水平散度项作用更大,中层正涡度主要受到垂直输送项的影响,由于两项的正贡献作用,中低层正涡度增加,低涡发展;中尺度低涡减弱时期,中低层水平散度项减弱,垂直输送项变为负值,正涡度制造减弱,负涡度制造增加,低涡消亡;(3)在MCS发展过程中,扰动有效位能的增加主要由于潜热释放作用和中尺度低涡纬向平均有效位能向扰动有效位能转化,之后扰动有效位能通过垂直运动向扰动动能转化,扰动动能的增加促进强对流发展。中尺度低涡不仅为强对流的发展提供能量,其中风场变化也是影响强对流演变的重要原因。
关键词中尺度低涡    MCS    东北冷涡    涡度收支    能量方程    
Analysis of the evolution reason of the mesoscale vortex and its influence on MCS under the background of the northeast cold vortex
QIAO Na , QIAN Peng , ZHOU Qing , ZHANG Chenxin , WU Xinyue     
Zhenjiang Meteorological Bureau of Jiangsu Province, Zhenjiang 212000
Abstract: An event of mesoscale vortex with MCS occurred in Shandong region on 3 August 2015 was investigated and analyzed using the WRF mesoscale numerical model based on NCEP/NCAR analysis data and Doppler radar observation data. The study revealed the reason for the occurrence and development of the mesoscale vortex and its effect on the mesoscale convection system under the environment of the northeast cold vortex. The results are as follows. (1) The upper-level jet stream' s strong divergence in front of the trough of the synoptic-scale northeast cold vortex and the low-level instability energy release provided favorable condition for the lower mesoscale vortex and convection occurrence. The generation of the low-level mesoscale vortex promoted the occurrence of the MCS. (2) When the mesoscale vortex forms and develops, the positive vorticity of the lower layer is mainly affected by the horizontal divergence term and the vertical transport term, with the former having a greater effect. The positive vorticity of the middle layer is mainly affected by the vertical transport term. Due to positive contributions of the two terms, the positive vorticity increases, and the vortex develops.When the mesoscale vortex weakens, the horizontal divergence term in the middle and lower layers weakens and the vertical transport term becomes negative, leading to the decrease of positive vorticity generation and the increase of negative vorticity generation, and therefore the vortex decreases. (3) When the MCS develops, the increase of perturbation effective potential energy mainly comes from the latent heat release and the conversion from zonal mean effective potential energy to disturbance effective potential energy, and then the perturbation effective potential energy is transformed into the perturbation kinetic energy through the vertical motion. The increase of the perturbation kinetic energy promotes the development of convections. The mesoscale vortex not only provides energy for the severe convection, but also influences the evolution of the strong convection.
Key words: mesoscale vortex    mesoscale convection system    northeast cold vortex    vorticity budget    energy equation    
引言

暴雨是中国主要的灾害性天气之一,国内外学者早在20世纪80年代就已经开始关注引发强降水的中小尺度系统。研究表明,中尺度对流系统(MCS)在同一个地区重复出现往往是产生持续性暴雨的直接原因(谌伟等,2017),MCS的发生发展不仅受到天气尺度系统的影响,同时受到中尺度系统的制约(Orlanski,1975贝耐芳等,2002),因此在研究强对流的同时,加强对不同系统之间相互影响的分析十分必要。

东北冷涡是指我国东北地区具有一定强度,能维持一段时间,且有深厚冷空气的高空冷性的气旋性涡旋。气象研究中常用的定义为500 hPa天气图上,115°—45°E、35°—60°N范围内,具有一定强度(至少有一根闭合等高线),能维持3 d或以上,且有冷中心或者冷槽配合的高空气旋性涡旋(郑秀雅,1992Zhao et al., 2007Fu and Sun, 2012)。东北冷涡的发生、发展、滞留、衰亡过程均对大气环流起着至关重要的作用,在其演变过程中,它与高空急流、锋面等系统相配合形成温压不对称结构,进而为强对流发生提供良好的环境场,并且它经常与中尺度系统相互作用,为强对流的发展提供有利的环流背景(白人海等,1998孙力等,2000张立祥等, 2008, 2009王宇欣等,2014)。大尺度流场为强对流的发生提供了有利的环流背景,而中尺度系统同样也是影响强对流发生发展的关键,研究表明,中尺度低涡的活动会直接影响强对流的演变及发展,进而引发暴雨的出现(东高红等,2013何光碧等,2014高守亭等,2019)。程麟生等(2001)利用MM5模式对“98.7”特大暴雨过程进行模拟指出,低层环境场中存在一个明显的中尺度低涡,低涡区域强烈的辐合作用是这次大暴雨出现的主要原因。不少学者指出对流层中低层正垂直涡度扰动有利于中尺度低涡的生成,而低涡辐合作用又为对流增强提供了有利的环境场,进而引发暴雨出现(徐亚梅等,2002王智等,2003孙建华等,2004王芬等,2015)。张庆红等(2000)对一次MCS过程进行数值模拟指出,当中尺度对流发生时,潜热释放会引发中尺度低压出现,该低压会直接影响低空急流的强度,进而影响对流的发展。王海东等(2008)指出低涡东移会引发大暴雨的生成。

以上研究成果表明,关于东北冷涡和中尺度低涡分别从环流场上影响强对流的分析较多,但是对于不同系统相互作用及其如何影响强对流发生发展的研究较少,因此本文利用WRF数值模式对东北冷涡背景下的一次中尺度低涡伴随MCS的过程进行模拟分析与诊断,深入了解东北冷涡影响下中尺度低涡发生发展的原因及其对MCS的物理影响机制,以期提高低涡影响下强对流预报的准确率。

1 资料和方法 1.1 选取资料

使用的资料包括:美国国家大气研究中心(NCAR)和美国国家环境预报中心(NCEP)提供的FNL (Final Operational Global Analysis)全球分析资料,时间间隔为每6 h一次,空间分辨率为1°×1°,属于GRIB2码全球资料;全国常规地面观测和探空资料;SWAN系统处理后的多普勒雷达资料,包括滨州、沧州等11个雷达站的观测资料,时间间隔为每6 min一次,空间分辨率为0.01°×0.01°;中尺度数值模式WRFV3.6 (ARW)输出的高分辨率模拟资料。

1.2 模式方案设计

本文采用非静力中尺度数值模式WRFV3.6进行数值模拟研究。模拟时间从2015年8月3日00时—3日18时(世界时,下同)共18 h,初始场及侧边界场采用NCEP 1°×1°逐6 h的FNL全球分析资料。模拟采用三重双向嵌套,中心经纬度为(117.037°E、36.316°N),三层嵌套网格数分别是349×349,454×496,598×595,水平分辨率分别是18.0 km、6.0 km、2.0 km,垂直方向共分为35层,第一、二层每1 h输出一次模拟结果,第三层每10 min输出一次模拟结果。采用的参数化方案包括: Thompson aerosol-aware微物理方案(胡向军等,2008),RRTM长波辐射方案(Mlawer et al., 1997),Dud? hia短波辐射方案(Lacis,1974),K-F积云对流方案(第三层不使用) (Kain et al., 1990, 1993),Monin-Obukhov近地面层方案(章国材,2004),YSU边界层方案(Hong et al., 2006),Noah陆面过程方案(赵向军,2017)等,文中主要分析第三区域D03的模拟结果。

2 天气过程实况 2.1 环流背景

从天气形势来看,2015年8月3日06时(世界时,下同) (图 1a),500 hPa呈现两脊一槽的形势,高压脊分别位于贝加尔湖和鄂霍次克海西部,槽线位于内蒙古高原-黄土高原-四川盆地一带。东北冷涡位于黑龙江以西,中心在(120°E、53°N)附近,闭合环流尺度达1 500 km以上,强度在556 dagpm以下,是深厚的冷性低压系统。低纬西北太平洋副热带高压(以下简称副高)脊线西伸至110°E附近并稳定维持,副高西北侧的西南暖湿气流为强对流天气发生提供良好水汽条件。200 hPa高空急流位于东北冷涡中心位置以南,中心最大风速达50 m·s-1以上,强度较强。中尺度低涡和对流发生在500 hPa东北冷涡槽前的偏西气流中,200 hPa高空急流入口区右侧,高层对应强辐散区。850 hPa (图 1c)由于对流区和高空急流右后方强烈的辐散加强低空急流左前方辐合,从而导致中尺度低涡(400 km)生成,随后10时(图略)低涡辐合作用进一步促进对流单体组织成带状中尺度对流系统(MCS)。

图 1 2015年8月3日06时(a)、12时(b)的200 hPa全风速场(阴影, 单位: m·s-1)和500 hPa高度场(黑色实线, 单位: dagpm)、温度场(红色实线, 单位: ℃)和风场(风向标, 单位: m·s-1)与06时(c)、12时(d)的850 hPa全风速场(阴影, 单位: m·s-1)、高度场(蓝色实线, 单位: dagpm)和流场(黑色流线, 单位: m·s-1) (红色三角表示强对流区, 棕色实线表示500 hPa槽线) Fig. 1 Wind speed field (shadow, unit: m·s-1) at 200 hPa and geopotential height field (black solid line, unit: dagpm), temperature field (red solid line, unit: ℃) and wind field (wind vane, unit: m·s-1) at 500 hPa at (a) 06∶00 UTC and (b) 12∶00 UTC on 3August 2015, and wind speed field (shadow, unit: m·s-1), geopotential height field (blue solid line, unit: dagpm) and flow field (black line, unit: m·s-1) at 850 hPa at (c) 06∶00 UTC and (d) 12∶00 UTC on 3 August 2015 (The red triangle represents severe convection zone and the solid brown line represents the trough-line at 500 hPa).

12时(图 1b),500 hPa东北冷涡缓慢东移,中心位置向东偏移,槽线移至东北平原—华北平原一带,200 hPa高空急流中心也略微东移,强度有小幅度减弱,中尺度低涡和MCS仍位于东北冷涡槽前、高空急流的右后方。850 hPa低空急流向东北方向移动(图 1d),风速增强,低涡中心略微东北移,伴随143 dagpm等高线范围的扩大而强度增加,高度延伸至800 hPa,MCS处于低空急流与中尺度低涡之间辐合区内。

2.2 中尺度低涡和MCS实况分析

低涡是低压在流场上的表现,雷达回波分布可以体现MCS的演变特征。通过上节分析可知,东北冷涡为中尺度低涡和对流的发生提供了良好的天气尺度环流背景,同时中尺度低涡为MCS的组织生成有重要贡献。由图 2可知,06时(图 2a)东北冷涡与高低空急流等系统共同作用使得对流单体在山东西部、河南北部生成,同时由于对流区上方对应假相当位温的密集带,锋生使得对流单体快速组织发展;高空急流右后方的辐散加强了低空急流左前方气流的辐合,同时对流区辐合作用推进了偏北风与西南风交汇,进而导致低层中尺度低涡在此处生成。12时(图 2b)中尺度低涡范围扩大,强度增加,它与低空急流间的切变辐合加强,风场的辐合作用使得散乱的对流单体组织成带状中尺度对流系统,即MCS生成发展。

图 2 2015年8月3日06时(a)、12时(b)的850 hPa实况雷达回波反射率(阴影, 单位: dBz)、流场(黑色流线, 单位: m·s-1)和假相当位温分布(红色实线, 单位: K, 黑色实心点代表沧州站) Fig. 2 Observed radar reflectivity (shadow, unit: dBz), flow field (black line, unit: m·s-1) and distribution of pseudo-potential temperature (red solid line, unit: K) at 850 hPa at (a) 06∶00 UTC and (b) 12∶00 UTC on 3August 2015 (The black solid point represents Cangzhou station).

由于14时左右为中尺度低涡和MCS发展强盛的阶段,因此图 3给出了该时刻沧州站上空的雷达探测资料,由径向风速(图 3a)可以看出,山东西北部存在完整的中尺度低涡形态,此时低涡正强盛发展。同时,由沧州站上空雷达反射率(图 3b)可见,紧邻低涡东南部的带状对流系统清晰,可见低涡与MCS的发展存在密切联系。

图 3 2015年8月3日14时沧州站仰角1.5°的径向风(a, 单位: m·s-1, 黄圈表示中尺度低涡)和仰角0.5°的雷达反射率(b, 单位: dBz) Fig. 3 (a) Radial wind (unit: m·s-1) at elevation angle 1.5° and (b) radar reflectance (unit: dBz) at elevation angle 0.5° from Cangzhou station at 14∶00 UTC 3 August 2015 (The yellow circle represents the mesoscale vortex).

总体而言,东北冷涡等天气尺度系统为中尺度低涡和强对流的发生提供了有利的环流背景,同时中尺度低涡的发展又是促进MCS生成的重要因素,因此此次过程中尺度低涡的演变成因和它对MCS发展的影响将是本文讨论的主要问题。

3 模拟结果检验

图 4为本次模拟得到的环流场,从中可见,06时(图 4ac) 500 hPa呈现两脊一槽形态,东北冷涡位于黑龙江西部,低纬副高脊线明显西伸,200 hPa高空急流位于东北冷涡南部,850 hPa存在明显的低空急流,高低层系统共同作用促使对流和中尺度低涡生成,与实况类似。12时(图 4bd) 各层系统东移,强度增加,MCS处于低空急流和中尺度低涡之间的切变辐合区内,与实况一致。

图 4 2015年8月3日06时(a)、12时(b)模拟的200 hPa全风速场(阴影, 单位: m·s-1)和500 hPa高度场(黑色实线, 单位: dagpm)、温度场(红色实线, 单位: ℃)和风场(风向标, 单位: m·s-1)与06时(c)、12时(d)模拟的850 hPa全风速场(阴影, 单位: m·s-1)、高度场(蓝色实线, 单位: dagpm)和流场(黑色流线, 单位: m·s-1) Fig. 4 Simulated wind speed field (shadow, unit: m·s-1) at 200 hPa and geopotential height field (black solid line, unit: dagpm), temperature field (red solid line, unit: ℃) and wind field (wind vane, unit: m·s-1) at 500 hPa at (a) 06∶00 UTC and (b) 12∶00 UTC on 3 August 2015, and simulated wind speed field (shadow, unit: m·s-1), geopotential height field (blue solid line, unit: dagpm) and flow field (black line, unit: m·s-1) at 850 hPa at (c) 06∶00 UTC and (d) 12∶00 UTC on 3August 2015.

通过模拟的中尺度低涡(图 5)可以看出:一个中尺度低涡于06时(图 5a)在113°E、38°N附近生成,此处正是对流发生区域,随后低涡中心东移,强度增加,范围扩大;12时(图 5b)低涡范围扩大到400 km,散乱的对流单体已组织成东北-西南向的带状中尺度对流系统(MCS);14时(图 5c)低涡和MCS向东北方向移动,强度和范围不断增加;15时30分(图 5d)低涡继续向东北方向移动,并逐渐减弱消失,MCS也逐渐东移减弱。通过与实况(图 23)对比可知,本次模拟较好地再现了中尺度低涡发生的位置、时间和演变过程,以及它与强对流的对应关系,因此可以利用本次模拟结果进行进一步的中尺度分析。

图 5 2015年8月3日06时(a)、12时(b)、14时(c)、15时30分(d)的850 hPa模拟雷达回波反射率(阴影, 单位: dBz)和流场(黑色流线, 单位: m·s-1, 紫色实心圆代表沧州站) Fig. 5 Simulated radar reflectivity (shadow, unit: dBz) and flow field (black line, unit: m·s-1) at 850 hPa at (a) 06∶00 UTC, (b) 12∶00 UTC, (c) 14∶00 UTC, and (d) 15∶30 UTC on 3 August 2015 (The purple solid point represents Cangzhou station).
4 中尺度低涡和MCS的涡度与能量分析

由上述分析可知,在东北冷涡的背景下,中尺度低涡的发生发展和MCS的演变关系密切,因此下面将从涡度方程出发探究中尺度低涡发生发展的原因,并利用能量转换方程分析低涡是如何影响MCS发展的。

4.1 中尺度低涡涡度收支分析

为了研究中尺度低涡的发生发展机制,下面利用涡度方程(式(1))对低涡区域进行涡度收支分析。p坐标系下的涡度方程(不考虑摩擦项) (朱乾根等,2007)

$ \begin{array}{l} \frac{{\frac{{\partial \varsigma }}{{\partial t}}}}{{\rm{A}}} = \frac{{\left[ {u\frac{{\partial \varsigma }}{{\partial x}} + v\left({\beta + \frac{{\partial \varsigma }}{{\partial y}}} \right)} \right]}}{{\rm{B}}}\frac{{ - \omega \frac{{\partial \varsigma }}{{\partial p}}}}{{\rm{C}}}\frac{{ - \left({f + \beta } \right)\nabla \cdot \vec V}}{{\rm{D}}}\\ \;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\frac{{ + \left({\frac{{\partial \omega }}{{\partial y}}\frac{{\partial u}}{{\partial p}} - \frac{{\partial \omega }}{{\partial x}}\frac{{\partial v}}{{\partial p}}} \right)}}{{\rm{E}}} \end{array} $ (1)

其中,A为相对涡度的局地变化项,B为水平平流项,C为垂直输送项,D为水平散度项,E为倾侧项,uv为水平风速,ω为垂直速度,ζ为相对涡度,f为地转涡度,β = $ {}^{\partial f}\!\!\diagup\!\!{}_{\partial y}\;$

4.1.1 局地相对涡度变化

为了更好的观察中尺度低涡发生发展过程中垂直方向上涡度的变化,选取低涡初生到消亡阶段每半小时进行时间平均,低涡区域5°×4°的范围进行区域平均,进而对相对涡度进行计算分析。从图 6可见,在中尺度低涡演变过程中,500 hPa以下始终以正涡度为主,最大正涡度集中在700—800 hPa,500 hPa以上以负涡度为主;低涡最高延伸至800 hPa,06—15时低涡不断加强,800 hPa以下相对涡度不断增大,850 hPa相对涡度最大增至8×10-5 s-1以上,可以推测在此期间有中层正涡度持续的垂直下传;15时之后,低涡开始减弱,800 hPa以下相对涡度有所减小,850 hPa相对涡度减小明显。

图 6 区域(5°×4°)时间(每0.5 h)平均的中尺度低涡相对涡度的垂直廓线(单位: 10-5s-1) (不同颜色的线条表示不同时次, 横虚线表示850 hPa高度, 横实线表示950 hPa高度) Fig. 6 The vertical profile of relative vorticity of the mesoscale vortex with region (5°×4°) and time (0.5 h) mean (unit: 10-5 s-1) (Different color lines indicate different times, and the horizontal dotted line represents 850 hPa, and the horizontal solid line represents 950 hPa).
4.1.2 涡度方程各项贡献

为了探究在中尺度低涡发生发展过程中涡度变化的机制,根据涡度方程,分别选取了中尺度低涡发生(06时)、强盛(14时30分)、减弱(15时)阶段的三个时次,计算垂直速度、散度和涡度方程中各项的变化来分析涡度变化的原因。

在中尺度低涡形成时的8月3日06时,由图 7a可知,低涡区域整层的垂直速度为正值,最大值在600 hPa附近,600 hPa以下为水平辐合,600 hPa以上为水平辐散,850 hPa的垂直涡度为3×10-5s-1;由图 7d可知,低涡区域中低层正的垂直涡度主要受到垂直输送项和水平散度项的影响,其中700 hPa以下以水平散度项的正贡献为主,700—500 hPa以垂直输送项的正贡献为主,中高层负的垂直涡度主要受到水平平流项和倾侧项影响。结合图 7ad可知,低层水平辐合引起垂直涡度的增加,垂直速度的峰值和垂直输送项的峰值均集中在600 hPa附近,这表明低层由于水平散度的辐合作用增加的正涡度,将通过垂直平流向上输送,从而有利于增强气旋性涡度的局地变化。

图 7 2015年8月3日06时(a)、14时30分(b)、15时(c)的中尺度低涡区域平均(5°×4°)和时间平均(每0.5 h)的垂直速度(W, 单位: 10-2m·s-1)、涡度(Vor, 单位: 10-5s-1)和散度(Div, 单位: 10-5s-1)的垂直廓线和06时(d)、14时30分(e)、15时(f)的涡度方程各项(单位: 10-9s-2)的垂直廓线(A为相对涡度的局地变化项, B为水平平流项, C为垂直输送项, D为水平散度项, E为倾侧项) Fig. 7 Vertical profiles on the mesoscale vortex region (5°×4°) and time (0.5 h) mean of vertical velocity (W, unit: 10-2 m·s-1), vertical vorticity (Vor, unit: 10-5 s-1) and divergence (Div, unit: 10-5 s-1) at (a) 06∶00 UTC, (b) 14∶30 UTC, and (c) 15∶00 UTC on 3 August 2015, and vertical profiles of terms of the vorticity equation (unit: 10-9 s-2) at (d) 06∶00 UTC, (e) 14∶30 UTC, and (f) 15∶00 UTC on 3 August 2015 (A represents local variation term of relative vorticity, and B represents horizontal advection term, and Crepresents vertical transport term, and D represents horizontal divergence term, and E represents heeling term).

在中尺度低涡强盛时的8月3日14时30分,由图 7b可知,低涡区域整层的垂直速度为正值,峰值有所增加,且位于300 hPa附近,300 hPa以下为弱辐合,300 hPa以上为强辐散,这可能与200 hPa高空急流的东移增强有关,850 hPa的垂直涡度增大至8×10-5 s-1;由图 7e可知,低涡区域中低层正涡度主要受到垂直输送项和水平散度项的积极作用,其中800 hPa以下以水平散度项的积极作用为主,且强度比06时有所增加,800—400 hPa以垂直输送项的积极作用为主,高层负的垂直涡度主要受到倾侧项影响。结合图 7be可知,低层水平散度项的正作用加强,使得低层低涡区域有正的相对涡度局地变化,正涡度加大,由于800 hPa以上垂直速度和垂直输送项均为正值且有所加强,因此将低层正涡度不断向中高层输送,使得800—400 hPa相对涡度局地变化加强,气旋性涡度增加。

在中尺度低涡减弱时的8月3日15时,由图 7c可知,低涡区域的上升运动有所减弱,峰值仍位于300 hPa附近,300 hPa以下有弱辐合,300 hPa以上有弱辐散,这与200 hPa高空急流的减弱相关联,850 hPa的垂直涡度减小至6×10-5 s-1;由图 7f可知,800 hPa以下水平散度项的正值减小,800—400 hPa垂直输送项的正值减小,850 hPa以下垂直输送项变为负值,中高层负垂直涡度仍主要受到倾侧项作用。结合图 7cf可知,800 hPa以下水平散度项的正贡献减小,垂直输送项变为负贡献,相对涡度的局地变化由正值变为负值,使得低层气旋性涡度开始减弱,并且中高层垂直运动和垂直输送项有所减弱,中尺度低涡开始减弱。

4.1.3 涡度方程各项变化分析

为进一步探究涡度方程各项变化,作各项的时间- 高度剖面如图 8。06—15时为中尺度低涡发生发展的阶段,水平散度项和垂直输送项正贡献最大,水平平流项和倾侧项负贡献最大,800 hPa以下水平散度项正贡献更明显,即辐合作用强,水平平流项负贡献更明显,800—400 hPa垂直输送项正贡献更明显,即向上输送强,倾侧项负贡献更明显。水平散度项在800 hPa以下对局地涡度变化有明显的积极作用,强辐合引起正的涡度收支,使得低层气旋性涡度增大,中层有较弱的消极作用,引起负的涡度收支,使得中层反气旋性涡度略微增大,这种低层气旋性涡度,中层反气旋性涡度的配置有利于涡旋和上升运动的发展。垂直输送项在800—400 hPa始终表现为很强的积极作用,使得局地相对涡度稳健增长,并将低层正涡度向上输送,有利于涡旋的发展。水平平流项在800 hPa以下对局地相对涡度为明显的负贡献,有正涡度向周围流失,800 hPa以上存在正贡献,周围有正涡度向涡旋中心输送,这种涡度平流下负上正的垂直分布形态,有利于上升运动的发展,从而进一步加强低层辐合。倾侧项在800— 400 hPa表现为负贡献,多有垂直涡度向水平涡度的转换,产生负涡度收支,不利于气旋性涡旋向上伸展。

图 8 中尺度低涡区域平均(5°×4°)和时间平均(0.5 h)的涡度方程各项中水平平流项(a)、垂直输送项(b)、水平散度项(c)、倾侧项(d)和相对涡度的局地变化项(e)的时间-高度剖面图(单位: 10-9s-2) Fig. 8 Time-height cross-sections for (a) horizontal advection term, (b) vertical transport term, (c) horizontal divergence term, (d) heeling term and (e) local variation term of relative vorticity(unit: 10-9 s-2) of the vorticity equation on the average mesoscale vortex.

15时之后中尺度低涡开始减弱,800 hPa以下水平散度项逐渐减弱,由正值转为负值,垂直输送项负值区有所伸展,水平平流项仍为负值,使得相对涡度的局地变化由正值变为负值,正的涡度收支有所减弱,负的涡度收支逐渐出现,从而引起低层的气旋性涡度环流减弱;800—400 hPa垂直输送项有所减弱,使得正的涡度收支减弱,倾侧项有所减弱,但仍表现为负值,引起局地相对涡度负变化。

综合上述分析,中尺度低涡形成和发展时期,低层的水平辐合和低层向高层正涡度的垂直输送,使得正的垂直涡度增加,气旋性环流加强,低涡发展;中尺度低涡减弱时期,低层辐合减弱、垂直输送变为负值、垂直输送项有所减小,使得正涡度制造减弱,负涡度制造增加,气旋性环流减弱,低涡消亡。

4.2 强对流带能量收支分析

上面利用涡度方程分析了中尺度低涡发生发展的原因,而伴随着低涡的演变,MCS生成并发展,那么中尺度低涡是如何影响强对流的,下面将利用Jiang等(1995)的能量收支方法,详细探讨强对流带与中尺度低涡间的能量转换。

将相关变量分为纬向平均场和扰动场,进而将大气动能和有效位能分成纬向平均部分和扰动部分来进行研究,其中有效位能可用来衡量在绝热过程中位能转化为动能的最大量。本节主要分析式(2)和(3)中涉及的三种转化过程和一种生成过程,其中忽略S(KE) 的作用,相应转化过程和生成过程的计算方法如式(4)—(8)所示。

$ \frac{{\partial {K_E}}}{{\partial t}} = C\left({{A_E} \to {K_E}} \right) + C\left({{K_Z} \to {K_E}} \right) + S\left({{K_E}} \right) $ (2)
$ \frac{{\partial {A_E}}}{{\partial t}} = C\left({{A_Z} \to {A_E}} \right) - C\left({{A_E} \to {K_E}} \right) + G\left({{A_E}} \right) $ (3)
$ C\left({{A_E} \to {K_E}} \right) = \frac{g}{{[\theta ]}}\left[ {{w^\prime }{\theta ^\prime }} \right] $ (4)
$ C\left({{A_Z} \to {A_E}} \right) = - \frac{g}{\sigma }\left[ {{v^\prime }{\theta ^\prime }} \right]\frac{{\partial [\theta ]}}{{\partial y}} $ (5)
$ \begin{array}{l} C\left({{K_Z} \to {K_E}} \right) = - \left[ {{u^\prime }{v^\prime }} \right]\frac{{\partial [u]}}{{\partial y}} - \left[ {{u^\prime }{w^\prime }} \right]\frac{{\partial [u]}}{{\partial z}} - \left[ {{v^\prime }{v^\prime }} \right]\frac{{\partial [v]}}{{\partial y}}\\ \;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\; - \left[ {{v^\prime }{w^\prime }} \right]\frac{{\partial [v]}}{{\partial z}} \end{array} $ (6)
$ G\left({{A_E}} \right) = \frac{g}{\sigma }\left[ {{H^\prime }{\theta ^\prime }} \right] $ (7)
$ \sigma = [\theta ]\frac{{d[\theta ]}}{{dz}} $ (8)

式(2)和(3)分别表示扰动动能和扰动有效位能随时间的变化,其中AE表示扰动有效位能,AZ表示纬向平均有效位能,KE表示扰动动能,KZ表示纬向平均动能,G(AE) 表示扰动有效位能在非绝热加热过程中的生成率,S(KE) 表示非保守外力对能量的变化率。C(XY)表示能量从X形式转化到Y形式。在以上方程中,方括号表示物理量的纬向平均,物理量右上角的撇号表示相对于该平均值的偏差,其中u'、v'、w ' 分别表示xyz方向上的扰动速度,θ' 表示扰动位温,H ' 表示绝热加热率的扰动。

图 9是根据上式计算出的模式第三层D03 (111°— 123°E、31°—41°N)区域时间平均的各种能量转换随时间的变化,具体算法如下:计算各变量纬向平均值、扰动值和扰动乘积的纬向平均值等;依据式(4)—(8)计算各转换项;计算各转换项的垂直积分;计算各转换项垂直积分后的区域时间平均;进而分析扰动动能和扰动有效位能的变化。其中由于中尺度低涡范围很大,基本覆盖D03区域,因此纬向平均能量可以代表中尺度低涡的能量,而扰动能量可以表示强对流的能量。由图可知,06—12时,扰动有效位能的生成项G(AE)和纬向平均有效位能向扰动有效位能的转化项C(AZAE)均为正值,且G(AE)在10时达到最大,最大值为20×10-4W,C (AZAE)在12时达到最大,最大值约为34×10-4W,扰动有效位能向扰动动能的转化项C(AEKE)较小,说明在此期间强对流获得并积累了较高的扰动有效位能,而非绝热加热项和纬向平均有效位能向扰动有效位能的转化项均是扰动有效位能增加的主要来源。同时纬向平均动能向扰动动能的转化项C(KZKE)也为正值,且不断增长,说明此段时间中尺度低涡不断向强对流区域输送动能,以助其发生发展。12—15时,扰动有效位能的生成项G(AE)和纬向平均有效位能向扰动有效位能的转化项C(AZAE)呈下降趋势,纬向平均动能向扰动动能的转化项C(KZKE)也开始减弱,但仍表现为正值,同时,扰动有效位能向扰动动能的转化项C(AEKE)呈现出明显的上升趋势,说明在此期间,扰动有效位能转化为扰动动能明显,扰动动能的增加促进MCS加强。15时之后中尺度低涡和强对流带开始减弱,不做具体分析。整个过程中尺度低涡和MCS之间的主要能量转化如图 10所示,这样的能量转化有利于中尺度对流系统的发展增强。

图 9 D03区域时间平均(每1 h)的垂直积分的不同形式能量转化随时间的变化(单位: 10-4W) (黑色实线表示C(AEKE), 红色长虚线表示C(AZAE), 深紫色短虚线表示C(KZKE), 青色点划线表示G(AE)) Fig. 9 The energy transformation chart of vertical integral different forms of D03 region and time (1 h) mean over time (unit: 10-4W), and the solid black line represents C (AEKE), and the long dashed red line represents C (AZAE), and the short dashed dark purple line represents C(KZKE), and the dotted blue line represents G(AE)

图 10 能量转化示意图 Fig. 10 Diagram of energy transformation.

依据上述分析,进一步分析风场对于强对流带的影响。图 11给出了D03区域时间平均的垂直积分的纬向平均风场和扰动风场随时间的变化,纬向平均风场代表中尺度低涡,扰动风场代表强对流。分析可知,11时之前,纬向平均风场主要以西风为主,伴有较弱的南风,垂直上升运动较强,且不断加大,扰动风场主要以西风为主,伴有较弱的北风分量,存在一定的下沉运动,由此可知,中尺度低涡输送的南方暖空气与对流区域的北方冷空气交汇,产生垂直运动,进而促进能量转化,强对流发生。11—15时,纬向平均风场西风减弱,南风加强,12时之后主要以西南风为主,垂直运动减弱,但仍为正值,扰动风场西风先增强后减弱,北风扰动逐渐变为南风扰动,14时之后主要以西南风为主,下沉运动减弱,出现上升运动,由此可知,由于中尺度低涡带来的强烈西南风作用,使得强对流带处的北风逐渐消失,随后主要受到西南气流的控制,水汽条件较好,且扰动下沉转为上升运动,进而促进强对流发展。

图 11 D03区域时间平均(每1 h)的垂直积分的纬向平均风场(a, 单位: m·s-1)和扰动风场(b, 单位: 10-2m·s-1)随时间的变化 Fig. 11 Vertical integral charts of (a) zonal mean wind field (unit: m·s-1) and (b) disturbed wind field (unit: 10-2 m·s-1) of D03 region and time (1 h) mean over time.
5 结论

由于中尺度低涡和MCS常引发暴雨等灾害,本文对2015年8月3日东北冷涡背景下的一次中尺度低涡伴随MCS过程进行数值模拟与分析,研究了中尺度低涡演变成因和它对强对流发生发展的影响,得到以下结论:

(1) 此次过程发生于200 hPa高空急流右后方辐散区,500 hPa东北冷涡南侧,槽前偏西气流中,850 hPa低空急流西北侧,高低层系统相配合为中尺度低涡和强对流的生成提供了有利的环流背景,同时,中尺度低涡的辐合作用促使散乱的对流单体组织成带状中尺度对流系统,即低涡的生成有利于MCS的组织形成。

(2) 利用涡度方程分析中尺度低涡发生发展的原因,低层正涡度的变化主要受到水平散度项和垂直输送项的影响,其中水平散度项作用更明显,而中层正涡度的变化主要受到垂直输送项的影响,高层负涡度的变化主要受到水平平流项和倾侧项的影响,其中倾侧项作用更显著。中尺度低涡形成和发展时期,低层的水平辐合和低层向上正涡度的垂直输送,使得正涡度增加,低涡发展;中尺度低涡减弱时期,低层辐合减弱、垂直输送变为负值,使得正涡度制造减弱,负涡度制造增加,低涡消亡。

(3) 利用能量收支方程分析中尺度低涡和MCS之间的能量转换,纬向平均场代表中尺度低涡,扰动场代表强对流,可知在MCS发展过程中,潜热释放作用和纬向平均有效位能向扰动有效位能的转化是扰动有效位能增加的主要来源,而风场变化所引起的垂直上升下沉运动是扰动有效位能转化为扰动动能的主要机制,扰动动能增加,强对流发展。中尺度低涡不仅为强对流的发展提供能量,低涡附近风场的变化也是影响强对流演变的重要原因。

参考文献
白人海, 谢安. 1998. 东北冷涡过程中的飑线分析[J]. 气象, 24(4): 37-40. DOI:10.3969/j.issn.1673-8411.1998.04.012
贝耐芳, 赵思雄. 2002. 1998二度梅期间突发强暴雨系统的中尺度分析[J]. 大气科学, 26(4): 526-540. DOI:10.3878/j.issn.1006-9895.2002.04.10
程麟生, 冯伍虎. 2001. "98.7" 突发大暴雨及中尺度低涡结构的分析和数值模拟[J]. 大气科学, 25(4): 465-478. DOI:10.3878/j.issn.1006-9895.2001.04.04
东高红, 韩素芹, 刘一玮, 等. 2013. 一次大暴雨过程中尺度涡旋系统特征分析[J]. 暴雨灾害, 32(2): 97-104. DOI:10.3969/j.issn.1004-9045.2013.02.001
高守亭, 周玉淑. 2019. 近年来中尺度涡动力学研究进展[J]. 暴雨灾害, 38(5): 431-439. DOI:10.3969/j.issn.1004-9045.2019.05.005
何光碧, 屠妮妮, 张利红. 2014. 一次低涡暴雨过程发生机制及其模式预报分析[J]. 暴雨灾害, 33(3): 239-246. DOI:10.3969/j.issn.1004-9045.2014.03.006
胡向军, 陶健红, 郑飞, 等. 2008. WRF模式物理过程参数化方案简介[J]. 甘肃科技, 24(20): 73-75. DOI:10.3969/j.issn.1000-0952.2008.20.028
谌伟, 岳阳, 刘佩廷, 等. 2017. 鄂东北一次特大暴雨过程的两个中尺度对流系统分析[J]. 暴雨灾害, 36(4): 357-364. DOI:10.3969/j.issn.1004-9045.2017.04.008
孙建华, 张小玲, 齐琳琳, 等. 2004. 2002年中国暴雨试验期间一次低涡切变上发生发展的中尺度对流系统研究[J]. 大气科学, 28(5): 675-691. DOI:10.3878/j.issn.1006-9895.2004.05.03
孙力, 安刚, 廉毅, 等. 2000. 夏季东北冷涡持续性活动及其大气环流异常特征的分析[J]. 气象学报, 58(6): 704-714. DOI:10.3321/j.issn:0577-6619.2000.06.006
王芬, 唐浩鹏, 陈晓燕. 2015. 黔西南一次低涡切变型暴雨的中尺度分析[J]. 沙漠与绿洲气象, 9(5): 41-46.
王海东, 李怀川, 吴正可. 2008. 一次低涡东移引发的大暴雨过程诊断分析[J]. 暴雨灾害, 27(4): 341-345. DOI:10.3969/j.issn.1004-9045.2008.04.010
王宇欣, 宋瑶. 2014. 东北冷涡引发的强雷暴个例分析[J]. 暴雨灾害, 33(3): 264-272. DOI:10.3969/j.issn.1004-9045.2014.03.009
王智, 翟国庆, 高坤. 2003. 长江中游一次β中尺度低涡的数值模拟[J]. 气象学报, 60(1): 66-77. DOI:10.3321/j.issn:0577-6619.2003.01.007
徐亚梅, 高坤. 2002. 1998年7月22日长江中游中β尺度低涡的数值模拟及分析[J]. 气象学报, 60(1): 85-95. DOI:10.3321/j.issn:0577-6619.2002.01.010
章国材. 2004. 美国WRF模式的进展和应用前景[J]. 气象, 30(12): 27-31.
张立祥. 2008. 东北冷涡中尺度对流系统研究[D]. 南京: 南京信息工程大学
张立祥, 李泽椿. 2009. 东北冷涡研究概述[J]. 气候与环境研究, 14(2): 218-228.
张庆红, 刘启汉, 陈受钧. 2000. 台湾海峡中尺度对流系统的数值研究Ⅳ: 动量收支[G]//海峡两岸及邻近地区暴雨试验研究, 北京: 气象出版社
赵向军. 2017. 飑线发展过程中水平涡度与垂直速度变化的特征分析及成因研究[D]. 南京: 南京信息工程大学
郑秀雅, 张廷治, 白人海. 1992. 东北暴雨[M]. 北京: 气象出版社, 129.
朱乾根, 林锦瑞, 寿绍文, 等. 2007. 天气学原理和方法[M]. 北京: 气象出版社, 108-116.
Fu S M, Sun J H. 2012. Circulation and eddy einetic energy budget analyses on the evolution of a northeast China cold vortex (NCCV) in May 2010[J]. J Meteor Soc Japan, 90(4): 553-573. DOI:10.2151/jmsj.2012-408
Hong S Y, Noh Y, Dudhia J. 2006. A new vertical diffusion package with an explicit treatment of entrainment processes[J]. Monthly Weather Review, 134: 2318-2341. DOI:10.1175/MWR3199.1
Jiang H, Raymond D J. 1995. Simulation of a mature mesoscale convective system using a nonlinear balance mode[J]. Journal of the Atmospheric Sciences, 52: 161-175. DOI:10.1175/1520-0469(1995)052<0161:SOAMMC>2.0.CO;2
Kain J S, Fritsch J M. 1990. A one-dimensional entraining/detraining plume model and its application in convective parameterization[J]. Journal of the Atmospheric Sciences, 47(23): 2784-2802. DOI:10.1175/1520-0469(1990)047<2784:AODEPM>2.0.CO;2
Kain J S, Fritsch J M. 1993. Convective parameterization for mesoscale models: the Kain-Fritsch scheme. The representation of cumulus convection in numerical models[J]. American Meteorological Society, 246: 165-170.
Lacis A A. 1974. A parametrization for the absorption of solar radiation in the Earth's atmosphere[J]. Journal of the Atmospheric Sciences, 31(1): 118-133. DOI:10.1175/1520-0469(1974)031<0118:APFTAO>2.0.CO;2
Mlawer E J, Taubman S J, Brown P D, et al. 1997. Radiative transfer for inhomogeneous atmospheres: RRTM, a validated correlated-k model for the longwave[J]. Journal of Geophysical Research Atmospheres, 102(14): 16663-16682.
Orlanski I. 1975. A rational subdivision of scales for atmospheric processes[J]. Bulletin of the American Meteorological Society, 56(5): 527-530. DOI:10.1175/1520-0477-56.5.527
Zhao S X, Sun J H. 2007. Study on cut-off low-pressure systems with floods over northeast Asia[J]. Meteorology and Atmospheric Physics, 96: 159-180. DOI:10.1007/s00703-006-0226-3