易文军,石兴琼,何炳文,梁 可
(1.张家界市气象局,湖南张家界 427000;2.石棉县气象局,四川石棉 625400;3.常德市气象局,湖南常德 415000)
湘西北山区是我国南方多雹灾地区,也是湖南省春季冰雹的易发区,由于其特殊的地形和气候特征,春季冰雹多发,雹灾较重。冰雹是一种具有突发性、范围小、持续时间短、破坏性强的灾害性天气,短时间内可以对农业、交通、电力、建筑和人民生命财产造成重大损失。近年来,国内外有关冰雹天气的研究已有很多重要成果。廖玉芳等[1]详细地讨论了我国S波段多普勒天气雷达三体散射的统计特征。朱敏华等[2]分析了强烈雹暴三体散射的多普勒天气雷达基本反射率因子的形成机制。汤兴芝等[3]归纳总结得出利用45 dBz雷达回波高度可较好地识别冰雹云。刘嘉慧敏等[4]通过商南地区一次强雹暴天气的环境条件和雷达回波结构特征分析,得出垂直积分液态含水量密度(VILD)跃增能更好地提前指示大冰雹。张昆等[5]对怀化地区近15年降雹时空分布以及环境参量特征分析总结,得出高架对流类是怀化地区降雹最主要的天气形势,且多发生在2—3月的午后或傍晚,冰雹直径较小。鲍向东等[6]对河南省冰雹天气环境参量特征分析统计,得出一些适合河南地区的降雹参量阈值。刘兵等[7]对张家界多个冰雹个例进行了分析,给出了冰雹预警的几个实用指标。本文利用常规观测资料和张家界多普勒雷达监测资料,对2020年3月21日晚—22日凌晨(北京时,下同)湘西北冰雹天气个例开展研究,着重分析三体散射(TBSS)、旁瓣回波、垂直积分液态含水量(VIL)和VIL密度等特征,以期为湘西北冰雹预警提供参考依据。
2020年3月21日20时—22日02时湘西北地区出现了强冰雹天气过程,并伴有短时强降雨、雷暴大风天气,其具有移速快、局地性强、冰雹颗粒大小分布不均的特点(图1a)。桃源县沙坪镇附近冰雹直径45 mm,永顺县郎溪乡、沅陵县五强溪镇和桃源县竹园村冰雹直径30 mm,持续时间9~18 min。此外,永定区城区和慈利县高桥镇均出现直径小于20 mm的小冰雹,持续时间约5~8 min。
从3月21日20时环流背景形势场的综合分析(图1b)可以看出,500 hPa上有一短波槽位于四川中部一带,在24 h内将快速东移影响湘北地区,850 hPa上湘北至贵北一带存在南北风向的切变线,该切变线以南西南急流发展旺盛,急流核最大风速达18 m/s。地面场上整个湘北地区受暖低压控制,鄂南一带有冷空气侵入暖低压,冷暖交汇加剧了不稳定能量的释放,此外,在850 hPa切变线附近有地面辐合线,有利于触发强对流天气。综合来看,湘西北地区在21日20时处于高温高湿场中,大气层结不稳定,在地面辐合线和冷空气的触发下,引发了此次强对流天气。
图1 2020-03-21T20:00—22T08:00强天气实况分布图(a)和2020-03-21T20:00环流背景综合分析(b)
图2为怀化、恩施3月21日20时的探空曲线。湘西北地区整层都是“上干下湿”的喇叭口层结特征,非常有利于强对流天气的发生。怀化站探空显示CAPE值1 180 J/kg,K指数38.4 ℃,SI指数-2.76 ℃,0 ℃层高度4.6 km,-20 ℃层高度7.0 km,0~6 km垂直风切达18 m/s,700 hPa以下为暖平流,中高层为冷平流,低层暖湿高层干冷,大气层结不稳定。恩施探空显示CAPE值611.3 J/kg,K指数36.4 ℃,SI指数-0.25 ℃,0 ℃层高度3.7 km,-20 ℃层高度6.8 km,0~6 km垂直风切达30 m/s。综合来看,0~6 km垂直风切变属于强垂直风切变,可以增强中层干冷空气的吸入,加强风暴中下沉气流和低层冷空气外流[8];“上干下湿”的层结不稳定,0 ℃层高度约4.2 km,-20 ℃层高度约6.9 km,适宜湘西北地区冰雹天气发生。环境条件有利于超级单体风暴的形成和发展以及大冰雹和雷暴大风的产生。
图2 2020-03-21T20温度对数压力图(a 怀化;b 恩施)
21日21:34永顺县大坝乡附近对流风暴开始发展,21:52最强反射率因子达63 dBz,回波顶高达11 km,从风暴开始发展到回波达到60 dBz用时仅18 min,可见该风暴发展非常迅速且强烈,在极短时间内发展成为超级单体风暴。21:52在2.4°、3.4°、4.3°仰角上均出现了三体散射“长钉”(TBSS)(图3a),风暴发展到22:16,在3.4°和4.3°仰角上出现了旁瓣回波(图3b)。23:16,该超级单体风暴分裂为两个超级单体风暴A和B,A风暴位于永定区后坪镇附近,B风暴位于永顺县和沅陵县交界处。A、B风暴最大回波强度均超过60 dBz,且B风暴比A风暴更强盛。A、B风暴移动路径的区域内均出现了不同程度的冰雹,并伴有雷雨大风天气。
图3 2020-03-21张家界雷达反射率因子(单位为dBz;a 2.4°,b 3.4°)
A风暴东移过程中略有发展后便迅速减弱,23:34在永定城区南部发展到最强,地面在23:56左右出现小冰雹,持续时间约5~8 min,22日00:04该风暴减弱。A风暴减弱消散以后,一个体扫内,在其消散地附近又有多个小尺度强风暴对流发展,最强回波超过65 dBz,22日01:13左右在慈利县高桥镇附近出现小冰雹。降雹后,风暴减弱消散。分析A风暴东移减弱消散的原因,一是地形因素,A风暴轨迹区域多为峡谷平原,而非利于降雹的丘陵山地;其次是降雹天气系统偏南,中低层切变线和地面辐合线更靠近B风暴东移路径。
B风暴从主体风暴分裂以后继续东移发展加强。该风暴沿着沅陵县北部东移,22日00:28,B风暴迅速加强,在2.4°、3.4°、4.3°和6.0°仰角上均出现了三体散射“长钉”和旁瓣回波(图4a~图4d),TBSS和旁瓣回波呈现加强的趋势,且维持时间达47 min。从风暴径向速度图上看,在桃源县西南部还出现了强中气旋(图4d),核区直径约6 km,旋转速度≥27 m/s。沅陵县五强溪附近冰雹直径达30 mm,桃源县沙坪镇附近冰雹直径达45 mm,并伴有7~9级大风。B风暴在22日01::33左右移出桃源县。
2.2.1 超级单体风暴分裂前 21:52,2.4°仰角上首次出现TBSS(图3a),而在6.0°仰角上并未看到TBSS(图略),说明此时冰雹云正在对流发展。随着超级单体风暴的不断发展,22:16,3.4°仰角上出现了明显的三体散射和旁瓣回波(图3b)。以≤15 dBz反射率因子计算,“长钉”长度为14 km,最大反射率因子超过了60 dBz。23:10,超级单体风暴发展到最强盛阶段,6.0°仰角出现了强的三体散射,最大反射率因子达72 dBz,“长钉”长度达23 km。TBSS的长度越长,表明反射率因子核心越大,高反射率因子的区域面积也越大[1]。23:18左右永顺县郎溪乡出现了直径30 mm的大冰雹。由于冰雹落地较快,从发现TBSS到地面降雹时间间隔较短,这给冰雹预警增加了难度。
23:16,超级单体风暴分裂为两个超级单体风暴A、B,三体散射迅速消失,从首次出现TBSS特征到结束,三体散射持续时间达84 min。A风暴在东移过程中并未加强,在永定区南部出现直径小于20 mm冰雹后,其迅速减弱消散,而B风暴继续保持超级单体结构东移。
2.2.2 超级单体风暴分裂后 B风暴沿着沅陵县北部山区维持70 dBz的最强回波东移,22日00:28,2.4°~6.0°仰角均出现了三体散射,在风暴移动的右后侧还出现了旁瓣回波(图略),而在低仰角0.5°和1.5°上,由于风暴径向外侧存在大于25 dBz的真实回波,无法清晰监测到三体散射。B风暴发展到此时,其超级单体风暴结构特征较分裂前更加明显和突出。
00:34,6.0°仰角上TBSS特征最清楚(图4a)。以≤15 dBz反射率因子计算,“长钉”长度达26 km,60 dBz的回波高度达9.1 km,在-30 ℃层高度以上,当≥50 dBz的回波高度超过-20 ℃层高度以上,有利于大冰雹的产生。00:52,2.4°~4.3°仰角上TBSS特征明显,6.0°仰角上“长钉”已经消失,只能看到旁瓣回波,此时TBSS特征最强盛,出现在4.3°仰角上(图4b),“长钉”长度达32 km,60 dBz的回波高度达8.5 km。00:57,最强反射率因子达74 dBz,在3.4°仰角上TBSS特征最明显(图4c),60 dBz的回波高度达8.7 km,“长钉”长度达47 km,这是本次降雹过程出现的最长“长钉”,在湘西北地区十分罕见,表明强反射率因子面积和反射率因子核心均达到了最强盛时期。01:03,在6.0°仰角径向速度图上出现了风暴顶强辐散(图略),正负速度差值达63 m/s。对比分析来看,“长钉”长度越长,冰雹直径越大。“长钉”特征随时间由高层向低层呈现,60 dBz的回波高度也随时间降低,反映了一个典型强烈超级单体风暴降雹的过程。
在速度图上监测到一个深厚持久的强中气旋,该中气旋持续了约58 min,00:52在3.4 °仰角上中气旋显示最强(图4d)。中气旋距离雷达约90 km,核区直径约6 km,旋转速度≥27 m/s,属于强中气旋。与此相应,00:50—01:30在桃源县西南部多个乡镇出现了7~9级大风。
从00:58反射率因子剖面图(图4e)可以看出,风暴呈现出典型的强烈超级单体风暴结构特征。反射率因子核强度超过了70 dBz,回波顶高达到16 km;强回波区结构密实紧凑,大于60 dBz的反射率因子伸展高度达到9 km,存在明显的回波穹窿和其上强大的悬垂回波;弱回波区右侧出现了回波墙以及狭长的反射率因子高值区,反射率高值区从9 km高度一直伸展到低层,表明地面即将出现冰雹。22日01:15左右,桃源县西南部出现了直径45 mm的大冰雹。此外,三体散射造成的回波假象在图中强回波区的上部也清晰可见。综合以上分析来看,B风暴呈现了一个典型的强烈超级单体风暴结构,受其影响多个乡镇造成了冰雹、雷雨大风和短时强降雨等灾害性天气。
图4 2020-03-21张家界雷达不同时刻不同仰角反射率因子(单位为dBz;a 6.0°;b 4.3°;c 3.4°)、00:52径向速度(单位为m/s;d 3.4 °;黑圆圈为中气旋)和00:58反射率因子垂直剖面(单位为dBz;e)
旁瓣回波也是大冰雹的一个预警指标[9]。从表1可以看出:不同仰角上首次出现旁瓣回波的时间基本一致,高仰角上旁瓣回波持续时间比低仰角更短;旁瓣回波最大切向长度达46 km,最大径向宽度达17 km,其出现时最强反射率因子≥73 dBz,表明风暴核反射率因子面积大,强雹暴对流发展旺盛;旁瓣回波出现时的高度为1.9~11.5 km,表明旁瓣回波在每个层次均可出现。由于低仰角上强反射率因子径向外侧有真实回波存在,TBSS特征并不清晰,而旁瓣回波却很明显,这就为冰雹预警提供了另一个可供参考的特征指标,即在强反射率因子径向外侧有真实回波存在时,可以依据旁瓣回波发布冰雹预警。
表1 湘西北2020-03-22凌晨B风暴旁瓣回波特征统计
垂直积分液态含水量(VIL)和VIL密度均是较好的冰雹预警指标。刘兵等[7]总结出湘西北大冰雹的VIL和VIL密度阈值分别为VIL≥43 kg/m2,VIL密度≥4.0 g/m3。图5给出了此次湘西北降雹过程中当VIL≥43 kg/m2时的最大VIL变化曲线及对应的VIL密度变化曲线。可以看出,21日22:46 VIL达54 kg/m2,VIL密度为4.3 g/m3,22日00:57 VIL达68 kg/m2,VIL密度为5.2 g/m3,均满足湘西北大冰雹的预警指标,在之后18~32 min地面出现30~45 mm大冰雹。值得注意的是,21日23:58 VIL达56 kg/m2,而VIL密度却仅为3.0 g/m3。超级单体风暴此时已减弱分裂为A、B两个风暴,地面亦无相应降雹报告,三体散射、旁瓣回波均未出现,说明冰雹预警不能仅看VIL≥43 kg/m2这个指标,而要更加注重VIL密度≥4.0 g/m3的判定指标。分析其原因,可能与雹胚成长有关。23:58的风暴顶高度达到了18.6 km,远超-40 ℃的高度,此时雹胚在充分增长之前就很快的被强上升气流抬升到-40 ℃高度以上,而-40 ℃高度以上没有过冷水滴,不存在有利于冰雹增长的环境[10],故而不会降冰雹。因此在实际业务工作中,尤其要关注VIL密度,它对大冰雹有更好的预警作用。
图5 2020-03-21—22湘西北最大VIL变化曲线(a)、VIL密度变化曲线(b)(黑圆点处表示VIL值虽很大但未降雹;黑三角处表示VIL和VIL密度最大值)
(1)“上干下湿”的喇叭口层结特征,低层高温高湿和不稳定能量输送,配合适宜的0 ℃层高度4.2 km,-20≥层高度6.9 km,在地面冷空气和辐合线的触发下,有利于冰雹、雷暴大风、短时强降雨等强对流天气的发生。此次降雹过程是典型的强超级单体风暴造成的,其具有明显的三体散射“长钉”、旁瓣回波等特征,以及较大的VIL和VIL密度,这些特征对冰雹天气有很好的预警作用。
(2)三体散射在不同仰角上均有出现,且持续时间长,TBSS最强反射率因子≥63 dBz,“长钉”长度达47 km,在其后的12~41 min内地面有大冰雹降落,TBSS特征可作为判定大冰雹的预警指标。三体散射特征越明显,表明强反射率因子面积越大、反射率因子核心越大。
(3)旁瓣回波在不同仰角上比TBSS出现的时间更长、更清晰,旁瓣回波出现时的最强反射率因子≥60 dBz,最大切向长度达46 km,最大径向宽度达17 km。旁瓣回波也是判定冰雹天气的预警指标,特别是当强反射率因子外侧有真实回波而TBSS特征不清晰时,旁瓣回波的预警作用更好。
(4)VIL≥43 kg/m2和VILD≥4.0 g/m3是冰雹天气的有效预警指标,当风暴对流厚度超过-40 ℃高度时,VIL密度能更好地判定冰雹。