华南“4.21”弓状回波的γ⁃中尺度涡旋雷达观测分析
1.
2.
Radar observation analysis of meso⁃γ⁃scale vortices in bow echo on 21 April in South China
1.
2.
通讯作者:
收稿日期: 2022-07-29
| 基金资助: |
|
Received: 2022-07-29
关键词:
Keywords:
本文引用格式
居圆圆, 徐昕, 唐滢.
Ju Yuanyuan, Xu Xin, Tang Ying.
弓状回波是一类典型的准线状中尺度对流系统(Quasi⁃Linear Convective Systems,QLCS),能够在地表产生破坏性的直线大风[1](Straight⁃Line Wind),最高风速可超过40 m·s-1,导致严重的人员伤亡和经济财产损失.早期Fujita[2]的研究认为,弓状回波顶点后侧下沉的后向入流急流(Rear⁃Inflow Jet,RIJ)是地表灾害性大风的主要成因.近年来,观测和数值模拟结果表明,在弓状回波内部常观测到γ⁃中尺度(2~20 km)[3]涡旋环流(在北半球以气旋式环流为主),其出现的位置多在弓状回波前端对流线(Leading Convective Line)的低层,被称为中涡旋(Mesovortex,MV).弓状回波近地面大风的形成同样与中涡旋密切相关,强度较强的中涡旋甚至能够造成龙卷风[4].
Wakimoto et al[5]基于机载多普勒雷达结合地面灾害调查结果发现中涡旋对弓形回波强灾害性大风发生的位置有重要作用,中涡旋的涡旋流与强的环境背景流叠加能够产生F0~F1[6]级强度的地表灾害性大风.Wheatley et al[7]也借助观测证实了中涡旋与地表灾害性直线大风存在密切关联.此外,关于中涡旋生成的位置,Weisman and Trapp[8]开展了中涡旋的理想数值模拟,将对流有效位能(Convective Available Potential Energy,CAPE)固定在2200 J·kg-1,而不断改变低层(0~2.5 km)和深层(0~5 km)的垂直风切变大小.当低层风切变较小时,生成的中涡旋沿着系统前端对流线,与大多数情况相一致,而采用较大的深层风切变时产生的中涡旋则沿着或者位于弓状回波顶端的北侧.Xu et al[9]将观测和数值模拟对比也得到了类似的结果,弓形回波中的中涡旋多位于其顶端的北侧,且该区域的中涡旋与其他位置相比生命史较长.该模拟研究同样发现中涡旋流对此次灾害性大风的产生有显著贡献,占比约为30%~50%.
关于中涡旋的生成存在多种解释,由于斜压性是导致许多旋转性结构产生的至关重要的原因,水平涡度的斜压项能够使其扭转为垂直,因此早期的研究将中涡旋的生成归因为斜压性的作用[10].Trapp and Weisman[11]基于理想数值模拟认为降水下沉拖曳导致沿着冷池出流边界的斜压水平涡度向下扭转从而产生了一对反向旋转的涡旋对,其中气旋式涡旋由于行星涡度的辐合拉伸增强生成气旋式中涡旋;Wheatley and Trapp[12]通过对2001年10月24日发生在美国密歇根州南部的弓形回波进行研究发现切变不稳定可以用来解释中涡旋的生成.Rotunno et al[13]提出中尺度冷池和低层环境垂直风切变的相互作用是飑线发展维持的重要机制,即“冷池⁃风切平衡理论(RKW Theory)”.Atkins and St. Laurent[14]利用该理论解释了强的低层环境风切变有利于产生强度更强的中涡旋,并且还认为强冷池以及科氏力强迫对中涡旋的生成具有重要作用.Xu et al[9]利用WRF模式对2009年5月8日发生在美国中部的弓形回波展开数值模拟研究,采用拉格朗日粒子轨迹追踪和三维涡度收支对其中两个较为典型的中涡旋进行分析,得出拖曳摩擦对低层中涡旋的生成也有重要贡献,这一结论在Schenkman et al[15]有关龙卷风生成的研究中也同样被提及.
过去的研究主要关注美国地区的弓状回波中涡旋及其形成机制.郑媛媛等[16]利用多普勒天气雷达对2002年5月27日发生在皖北的一次超级单体风暴过程进行分析,发现超级单体低层径向速度图上有明显的中气旋特征,并在之后加强,旋转速度可达22 m·s-1,与钩状回波的缺口区相对应.Tang et al[17]将对流单体中产生的γ⁃中尺度涡旋,如中气旋,称为单体类型的中涡旋(Cell⁃Type MV),并利用多年雷达观测对我国江淮地区暖季的中涡旋进行统计分析,揭示了该地区中涡旋的时空分布特征(如月变化、日变化)和大小、强度、生命史等特征.但是Tang et al[17]研究主要关注江淮地区的中涡旋,对于华南地区中涡旋的特征仍缺乏深入了解.对此,本研究利用华南地区的多普勒天气雷达观测对一次弓状回波产生的γ⁃中尺度涡旋进行分析.
1 数据与方法
1.1 数据
图1
图1
广州雷达站点位置(GZRD)及观测范围图
The black circle represents the coverage with a radius of 230 km and the radar composite reflectivity of the bow echo (shading) at 06:00 (unit: dBZ) and the red box indicates the domain of 2°×2° centered on the radar site
Fig.1
The location and range of Guangzhou radar site (GZRD)
1.2 MV识别和追踪
本文主要通过广州雷达来研究本次天气过程中产生的中涡旋,该识别系统使用的水平网格分辨率为500 m,共计953×953个格点,由于中涡旋为γ⁃中尺度对流系统,因此能够满足对中涡旋识别精度的要求,同时也能够更精细化地再现此次天气过程.
中涡旋的识别方法采用的是线性最小二乘导数法(Linear Least Square Derivative,LLSD)[22],且要求方位切变强度大于10-3 s-1时才可被识别.由于中涡旋一般生成于对流层低层,因此只计算0.5°仰角面的中涡旋的切变,并将识别区域限制在距雷达中心半径150 km的范围内,大于该距离雷达辐角方向的分辨率较低且延伸的高度较高,识别的可信度较差.同样是针对径向速度的计算,但该方法与过去常用的通过径向速度场出现的近距离正负速度偶判断中涡旋的方法相比更为先进,能够缓解径向速度噪音误差带来的干扰.
中涡旋的追踪使用的是客观的追踪方法,通过前向与后向追踪相结合得到中涡旋持续的时间.本研究将前后两个时次中涡旋的相关系数设为0.8,大于等于该值则将其看作同一个中涡旋展开追踪,若出现一对多或者多对一的情况,则选择相关系数最大的作为其在该时刻的延续.由于雷达扫描时次为6 min,根据前向和后向的追踪要求,本文在识别到的众多中涡旋中选取持续时间至少为18 min作为研究对象,由此得到表1.
表1 系统识别的中涡旋特征
Table 1
| 序号 | 经度(°E) | 纬度(°N) | 方位切变(10-3 s-1) | 持续时长(min) | 初次识别时间 (UTC) | 相对弓状回波 顶点的位置 |
|---|---|---|---|---|---|---|
| 1 | 111.89 | 22.98 | 5.46 | 18 | 3:24 | M |
| 2 | 111.95 | 23.23 | 2.89 | 24 | 3:36 | M |
| 3 | 112.59 | 23.69 | 3.06 | 30 | 3:48 | N |
| 4 | 112.89 | 23.68 | 3.23 | 18 | 4:24 | N |
| 5 | 112.88 | 23.57 | 2.64 | 18 | 4:24 | N |
| 6 | 112.98 | 23.60 | 3.02 | 18 | 4:36 | N |
| 7 | 113.21 | 23.61 | 5.46 | 24 | 4:54 | N |
| 8 | 112.72 | 23.02 | 3.1 | 24 | 4:54 | M |
| 9 | 113.34 | 23.59 | 3.53 | 30 | 5:18 | N |
| 10 | 113.34 | 23.56 | 3.0 | 18 | 5:18 | N |
| 11 | 113.09 | 23.04 | 2.92 | 18 | 5:30 | M |
| 12 | 113.06 | 22.54 | 2.48 | 18 | 5:54 | S |
| 13 | 113.67 | 23.04 | 3.83 | 18 | 6:18 | M |
| 14 | 113.58 | 23.05 | 2.62 | 24 | 6:24 | M |
| 15 | 113.06 | 22.34 | 2.5 | 18 | 6:30 | S |
| 16 | 113.23 | 22.25 | 2.61 | 18 | 6:48 | S |
| 17 | 113.43 | 22.29 | 2.69 | 48 | 7:00 | S |
| 18 | 113.48 | 22.36 | 2.52 | 18 | 7:00 | S |
| 19 | 113.1 | 22.11 | 2.51 | 18 | 7:00 | S |
| 20 | 113.54 | 22.54 | 3.43 | 30 | 7:06 | S |
“N”表示北侧,“M”表示中间,“S”表示南侧
2 中尺度天气系统概况及过程分析
2.1 个例概述
该弓状回波在广东省境内历时约7 h,图2给出的雷达组合反射率是该过程不同时刻的形态特征及位置.2017年4月21日02时(世界时,以下同),该弓状回波从广西自西向东移动,03:30时整体位于广东省境内,处于发展阶段,弓状回波特征明显.06:00时,弓状回波与周围的分散性对流单体合并成一个整体,处于成熟阶段,回波强度最大可达55 dBZ,从雷达图像上来看回波主体位于珠三角附近(图1).07:30时,系统弓状回波特征逐渐消失,南端位于海岸线附近.09:00时,该弓状回波东移入海,对广东的影响逐渐消失.广东省大部分地区都遭受了此次强对流天气的影响,有多地出现大到暴雨,并伴有短时雷暴大风(≥17.2 m·s-1).05:00~07:00,地面大风集中出现在具有明显弓状结构的区域附近,其中最大阵风位于广州市海珠区华洲街,可达32.4 m·s-1 (11级).此外,部分站点观测到有冰雹产生[23].
图2
图2
多时刻雷达组合反射率拼图(阴影,2017年4月21日03:30~07:30,间隔1 h,单位:dBZ),中涡旋位置及移动路径
The blue dotted lines represent the location of MVs and their
moving path
Fig.2
Hourly reflectivity from 0330 UTC to 0730 UTC 21 April 2017 (shading,unit:dBZ),the location of MVs and their moving path
2.2 中涡旋的时空分布特征
从时间上来看,符合要求的中涡旋最先在03:24时被识别到,最后一个被观测到的时刻为07:06,分布上相对均匀,最密集出现的时刻为07:00,有三个.识别到的中涡旋持续时间多维持在约18 min,占60%,持续时间最长的中涡旋出现在07:00,共计48 min,出现在系统南侧,位于海陆交界处,这与Xu et al[9]对美国发生的弓状回波事件中产生的中涡旋多位于顶点北侧的研究结论有差别.
2.3 中涡旋的雷达径向速度特征
针对不同类型中涡旋移动路径的空间分布特征(图2),选取典型时刻结合雷达径向速度展开进一步的分析.广州多普勒雷达站点位置位于113°21′E,23°N(图1),其中北侧代表性中涡旋选取的时刻为03:48和04:24(图3a~b),中间选取的时刻为03:24和06:18(图3c~d).03:48时(图3a),中涡旋位于系统北侧,径向速度大小在5~10 m·s-1的区域,该区域径向速度等值线分布密集.04:24时(图3b),系统北侧识别的中涡旋有两个.与图3a相比,它们出现的位置偏向大小为5 m·s-1的径向速度的等值线之间,处于速度辐合区,等值线较为稀疏,但该时刻中涡旋西侧的地面大风最大可达21.8 m·s-1.03:24时(图3c),系统弓状回波特征明显,后侧边缘径向速度最大可达20 m·s-1,有明显的后向入流急流缺口,中涡旋位于弓状回波顶点附近,方位切变强度最大为5.46×10-3 s-1并伴有风速超过20 m·s-1的地面大风,由此推测该区域的局地大风可能是中涡旋和后向入流急流共同作用的结果.仅从雷达分析无法准确了解中涡旋的生成机制,后续将利用高分辨率数值模拟作进一步的研究.相比之下,虽然图3d径向速度为20 m·s-1,地面大风实际维持在22 m·s-1左右,但方位切变强度弱,后向入流急流缺口不明显,仅在反射率回波上表现出局地化特征.
图3
图3
雷达组合反射率(阴影,单位:dBZ),径向速度(等值线,间隔5 m·s-1)及中涡旋(绿色实心三角形)位置分布图:
(a) 03:48,(b) 04:24代表北侧;(c) 03:24,(d) 06:18代表中间
Fig.3
Radar composite reflectivity (shading,unit:dBZ),radial velocity (contours,interval 5 m·s-1) and MV (green solid triangle) locations at (a) 0348 UTC and (b) 0424 UTC in the north of the QLCS, (c) 0324 UTC and (d) 0618 UTC in the middle
The solid lines of the contour represent positive values (away from the radar),and the dashed lines represents negative values
(towards the radar). The value of radial velocity equal to 5 m·s-1 is in blue,and greater than 5 m·s-1 in red
对于位于弓状回波系统南侧的中涡旋,07:00时不仅同时观测到三个中涡旋,且生命史最长的中涡旋也位列其中.因此选取该时刻的两个中涡旋作为研究对象,按照时间变化选择雷达连续六个扫描时次(图4).07:00时(图4a),初次观测的两个中涡旋均位于5 m·s-1径向速度的等值线之间,地面大风区在南侧中涡旋的西南方向,最大风速可达23.2 m·s-1.北侧的中涡旋生命史较短,为18 min.在该时段内,朝向雷达方向的径向速度等值线(蓝色虚线)覆盖的范围逐渐缩小至偏南的生命史最长的中涡旋一侧,因此认为径向速度辐合逐渐减弱直至消失是图4a中北侧中涡旋消失的重要原因.持续时间较长的中涡旋伴随系统逐渐东移,周围的径向速度逐渐变大,位于5 m·s-1与10 m·s-1之间,且10 m·s-1的等值线逐渐东移靠近(图4e~f).
图4
图4
雷达组合反射率(阴影,单位:dBZ),径向速度(等值线,间隔5 m·s-1)及中涡旋(绿色实心三角形)位置分布图:
(a) 07:00,(b) 07:06,(c) 07:12,(d) 07:18,(e) 07:24,(f) 07:30代表南侧
Fig.4
Radar composite reflectivity (shading,unit:dBZ),radial velocity (contours,interval 5 m·s-1) and MV (green solid triangle) locations at (a) 0700 UTC,(b) 0706 UTC,(c) 0712 UTC,(d) 0718 UTC,(e) 0724 UTC and (f) 0730 UTC in the south of the QLCS
2.4 中涡旋的环境特征
Weisman and Trapp[8]通过理想数值模拟研究指出,较强的环境垂直风切和较大的CAPE均有利于弓状回波中涡旋的生成.图5分别给出不同时刻南北两侧的垂直风切变(0.1~3 km)和CAPE分布情况.两个时刻的垂直风切变方向都是由西南指向东北,且低层风主导方向为西南风有利于低层来自海洋的水汽输送(图6),从而导致CAPE值自南向北逐渐减少.05:00时(图5a),弓状回波自中间向北位于垂直风切变大于17 m·s-1的区域,该范围内0.1~1 km(图略)和0.1~3 km的CAPE值均小于1500 J·kg-1,这表明在弓状回波发展阶段强的环境垂直风切变对中涡旋的生成和发展具有更大的贡献.07:00时(图5b),系统南侧中涡旋生成区域的风切小于15 m·s-1,但此时CAPE值最大,可达2000 J·kg-1.低层0.1~1 km具有相似的环境特征,垂直风切变为7 m·s-1,CAPE同样位于大值区(图略),说明该弓状回波发展到后期高CAPE值对中涡旋的生成有更明显的作用.虽然与北侧相比,该环境条件中产生的中涡旋有较长的生命史,但平均强度较弱.
图5
图5
(a) 05:00,(b) 07:00的对流有效位能(灰色阴影,单位:J·kg-1),0.1~3 km垂直风切变(箭头)及其大小(等值线,单位:m·s-1)叠加对应的雷达组合反射率(彩色阴影,单位:dBZ):(a)中蓝色实心三角形代表03:48,04:24,04:36,04:54时刻的中涡旋位置,(b)中蓝色实心三角形代表06:30,06:48,07:00时刻
(b) 0700 UTC: the blue triangles in (a) represent the MV positions at 0348,0424,0436 and 0454 UTC,
and the blue triangles in (b) at 0630,0648 and 0700 UTC
Fig.5
Convective available potential energy (gray shading,unit:J·kg-1),0.1~3 km vertical wind shear (vector) and values (contour line,unit:m·s-1),and composite reflectivity (colorful shading,unit:dBZ ) at (a) 0500 UTC and
图6
3 结论与讨论
本文利用广州业务多普勒天气雷达观测对2017年4月20日16时至21日16时华南一次弓状回波过程中的中涡旋进行研究,结合ERA5再分析资料对识别到的中涡旋的环境特征进行初步分析.主要结果如下:
(1)弓状回波在2017年4月20日02时从广西省移入,03:30时对流系统整体位于广东省境内,并在此后一段时间内具有明显的局地化弓形回波特征.09:00时,系统处于消散阶段并东移入海.该过程导致地面大风及冰雹等灾害性强对流天气,对广东省大范围区域都有影响.
(2)本次天气过程中识别到符合条件的中涡旋共计20个,根据相对弓状回波顶点的位置分为三类(南侧、北侧和中间),数量上并无明显差异.北侧中涡旋多发生在系统的发展阶段,平均方位切变强度为3.42×10-3 s-1,而南侧出现在弓状回波消散阶段,平均方位切变强度偏弱,为2.68×10-3 s-1.中间的中涡旋常伴有明显的局地化弓形回波结构,平均方位切变与北侧相当,雷达径向速度显示其后侧有RIJ缺口,且周围的径向速度最大.南侧的中涡旋位于径向速度辐合区,且最长生命史的中涡旋也出现在该侧,说明风场辐合有利于中涡旋的长时间维持.
(3)强垂直风切变以及高CAPE值均有利于中涡旋的生成.在弓状回波发展阶段,北侧环境特征为高风切、低CAPE,受高风切影响产生强度较强的中涡旋;消散阶段,南侧环境特征表现为低风切、高CAPE,其中高CAPE对中涡旋的生成有重要贡献,但生成的中涡旋强度相对较弱.
(4)区域平均的风场廓线特征在925~600 hPa变化最为明显,深厚的低空西南急流水汽输送有利于弓状回波的维持与发展.
本研究主要基于雷达观测初步分析了中涡旋位置、强度等基本特征.未来将开展高分辨率数值模拟,详细分析本次弓状回波中涡旋的形成机制及其对地表大风的影响.
参考文献
Vortex structure and evolution within bow echoes
Part Ⅰ:Single⁃Doppler and damage analysis of the 29 June 1998 Derecho.
Manual of downburst identification for project NIMROD
A rational subdivision of scales for atmospheric processes
A concentrated outbreak of tornadoes,downbursts and microbursts,and implications regarding vortex classification
High winds generated by bow echoes
Part Ⅰ:Overview of the Omaha bow echo 5 July 2003 storm during BAMEX.
Tornadoes and downbursts in the context of generalized planetary scales
Radar and damage analysis of severe bow echoes observed during BAMEX
Mesovortices within the 8 May 2009 bow echo over the Central United States:Analyses of the characteristics and evolution based on Doppler radar observations and a high⁃resolution model simulation
On the Rotation and propagation of simulated supercell thunderstorms
The effect of mesoscale heterogeneity on the genesis and structure of mesovortices within quasi⁃linear convective systems
A theory for strong,long⁃lived squall lines
Tornadogenesis in a simulated mesovortex within a mesoscale convective system
一次典型超级单体风暴的多普勒天气雷达观测分析
Analysis of a strong classic supercell storm with Doppler weather radar data
Characteristics of low⁃level meso⁃γ⁃scale vortices in the warm season over East China
The ERA5 global reanalysis
Efficient assimilation of radar data at high resolution for short⁃range numerical weather prediction
∥Extend Abstract.
The use of radial velocity derivatives to diagnose rotation and divergence
∥
/
| 〈 |
|
〉 |

