伊犁河谷冰雹及其雷达回波特征分析
2022-02-03孙鸣婧刘俊哲陈春艳施俊杰郑育琳
孙鸣婧,刘俊哲,陈春艳,施俊杰,郑育琳
(1.中国气象局乌鲁木齐沙漠气象研究所,新疆 乌鲁木齐 830002;2.新疆气象台,新疆 乌鲁木齐 830002;3.巴州气象局,新疆 库尔勒 841000)
冰雹是新疆的主要灾害性天气之一,常给农牧业生产造成惨重损失,严重制约着新疆粮棉、特色林果业的发展[1]。由于冰雹空间尺度小、生命史短、突发性强,准确预报冰雹一直是短时临近预报中的难点和重点。目前,天气雷达是国内外冰雹预警及地面人工防雹作业最有效和最快速的探测手段[2-3],预报员如何根据雷达探测资料识别雹云回波、确定落区及强度对冰雹的预报和预警十分重要。
我国常用回波强度、回波顶高、累积垂直液态含水量来做冰雹预警,然而在不同地区,总结出的冰雹预警指标也不尽相同,有些指标甚至差异很大[4-8]。福建龙岩将回波强度≥60 dBZ、回波顶高≥8 km作为冰雹预警指标[9];安徽地区冰雹云的最大回波强度基本在55 dBZ以上,回波顶高平均在12.1 km[10];湘西北山区夏季冰雹云的回波顶高均在9 km以上,垂直液态水含量一般为30~35 kg/m2[11]。以往的研究表明,在新疆的不同地区,冰雹的雷达回波表现特征并不一致,如在新疆喀什地区冰雹预警指标为组合反射率≥50 dBZ、回波顶高≥9 km、垂直液体水含量≥20 kg/m2等[12];在阿克苏地区判别冰雹的初步指标为反射率因子为50.7 dBZ、回波顶高为9.4 km、垂直积分液态水含量为12.8 kg/m2[13]。由于冰雹的强局地性,建立合理的分区预报、预警指标才能有效提高冰雹预报预警能力。
新疆天山西部的伊犁河和特克斯河流域是当地冰雹重要源地,以往对此地冰雹的研究多为昭苏县冰雹个例分析[14],对于河谷其他县市的冰雹分析或雷达回波指标凝练很少,日常业务中科学支撑薄弱。本文将统计近12 a发生在伊犁河谷地区冰雹的时空分布、雹云雷达回波产品各项指标,并分析典型雹云个例的雷达回波特征,为伊犁河谷地区冰雹的临近预警提供参考。
1 资料和方法
冰雹常由中小尺度天气系统造成,其尺度小于现行的气象观测网,因此使用气象台站观测到的地面降雹资料存在一定局限性。为最大限度减少冰雹观测数据遗漏,本文收集了中国气象局灾情直报信息作为冰雹实况的补充。2010—2021年共统计出冰雹152次(其中2021年只有灾情资料,未收集到地面观测资料,因此在分析时间变化时所用资料时间为2010—2020年),根据其发生时间和地点分析其时空分布特征。剔除雷达资料缺失和回波上无明显对流风暴的个例,最终得到23次有较完整冰雹实况和雷达资料的过程(表1)。一次冰雹过程可能受多个对流单体影响,23次冰雹过程共32个雹暴单体,通过追踪这32个雹暴的演变,归纳总结伊犁河谷地区冰雹的雷达回波特征。
表1 有完整雷达探测数据的冰雹个例集
2 伊犁河谷冰雹的时空分布特征
2.1 年、月变化
伊犁河谷每年都有冰雹发生,但次数分布不均匀,逐年变化呈多波动型,可预测性低。2010—2020年平均每年14次,其线性趋势呈略减少态势。2013年出现冰雹最多,2014年最少(图1)。
图1 2010—2020年冰雹次数的年(a)、月(b)变化
伊犁河谷降雹都发生在4—10月,集中出现在4—9月,与新疆汛期时间一致,与内地沿海城市有较大差异[15]。其中5—8月较多,占比86%,6月出现冰雹频次最高,占比30%;其次是7月,占比29%。记录中,最早的冰雹出现在4月5日(2010年新源县),最晚出现在10月1日(2011年特克斯县、昭苏县)。
2.2 日变化
伊犁河谷冰雹日变化特征显著,主要出现在午后到夜间(13—23时),占比95%,其中16—19时为冰雹高发期,占比59%;最易出现冰雹的时次是17时,占比20%(图2)。与我国东部城市相比,伊犁河谷出现冰雹的时间跨度大(前半夜也易发生),其他省市多发生在傍晚,如安徽的冰雹主要集中在15—18时[10]、唐山地区冰雹高发期在16—18时[16]等。伊犁河谷冰雹易发时间与新疆其他地市也不完全相同,阿克苏地区冰雹易发时段为15—20时,高发时段为18—20时[13],这对新疆强对流天气的监测及预警提出了更大的挑战。
图2 冰雹次数的日变化
2.3 空间分布
昭苏县出现冰雹次数最多,11 a共有84次,占伊犁河谷地区冰雹总数(152次)的一半以上,平均每年发生冰雹8次;其次是特克斯县,共19次,1 a约8次;霍城县次之,共14次;尼勒克县共出现9次、霍尔果斯市出现7次、察布查尔县和伊宁市出现6次、新源县出现4次、巩留县出现2次。
图3为多普勒雷达监测的32个雹暴单体的移动路径图。冰雹多发在南部沿山的昭苏县和特克斯县及北部沿山的霍城县、尼勒克县和霍尔果斯市,这与河谷的山脉地貌息息相关,山脉的地形强迫抬升更利于强对流天气的发生。雹云的移动路径基本与山脉走向一致,在伊犁河谷北部雹云自西北向东南移动,在南部雹云自西南向东北移动。
图3 伊犁河谷冰雹云移动路径
3 冰雹云的雷达回波特征
3.1 回波强度
冰雹云是强对流风暴的产物,根据微波散射理论,当产生的冰雹尺寸越大时,其回波也会越强,因此冰雹云的最强回波常被作为冰雹预报预测指标。在伊犁河谷地区,发展旺盛的冰雹云回波强度可达70 dBZ以上。表2统计了32次冰雹云的最大反射率因子值,其范围为49~72 dBZ,与其他地区相比,最强回波反射率因子跨度大。总体来看99%的最大反射率因子在50 dBZ以上(美国国家气象局环境预报中心给出的参考值也是50 dBZ[17]),冰雹云最强回波常分布在55~60 dBZ,占比47%,与喀什地区一致[12]。因此在日常业务监测伊犁河谷地区雷达回波时,如反射率因子发展到50 dBZ以上时,就应警惕出现冰雹的可能性。
表2 反射率因子强度分级统计
3.2 回波高度
回波顶高是用来分析对流云强弱和云体内上升气流伸展高度的指标。伊犁河谷地区冰雹云的最大回波顶高为13.3 km,最小为6.8 km,平均回波顶高为11 km,与甘肃、湘西北地区、喀什地区相当。99%的冰雹云回波顶高在8 km以上,冰雹云出现在11~12 km最多,占比31%。几乎所有冰雹云的回波顶高都在8 km以上,因此当回波顶高发展到8 km时,就应警惕出现冰雹的可能,及时做好预警。
回波顶高的逐月分布呈单峰型,峰值在6月,即伊犁河谷地区的冰雹云在6月发展最高(图4)。4—6月,回波顶高随时间推移增大,冰雹云的强度越来越强,而7—8月回波顶高则呈下降趋势,冰雹云的强度逐渐减弱。产生冰雹的强对流风暴最显著的特征是反射率因子高值区向上扩展到较高的高度,伊犁河谷地区50 dBZ的强回波核也同样在6月发展的最高,2021年6月6日50 dBZ强回波核伸至10.6 km。
图4 伊犁河谷冰雹云(强)回波顶高
3.3 冰雹云垂直累积液体水含量
垂直累积液体水含量(VIL)是将反射率因子数据转化为等价的液体水值,它对于判别带有大冰雹的风暴有重要参考意义,如果有VIL值异常大,则发生冰雹的可能性很大[17]。将冰雹云单体发展过程中VIL的最大值记为VILmax,表3统计了伊犁河谷31个单体的VILmax(剔除2017年4月14日伊宁个例,雷暴处于静椎区,VIL值不可用),其平均值为19.6 kg/m2。不同雹暴VILmax差异很大,最小值为5.02 kg/m2,出现在2020年8月11日的昭苏县,最大冰雹直径0.5~1 cm;最大值出现在2021年6月6日的特克斯县,为47.7 kg/m2(无冰雹直径记录)。有记录的最大冰雹直径3 cm个例中(2016年5月28日的霍尔果斯市)对应单体VILmax仅27.1 kg/m2,可见单体的VILmax与降雹直径没有明显的对应关系。
表3 单体最大VIL值分段统计
多项研究指出VIL的变化,尤其是跃增的特性,可作为判断降雹的指标之一。图5给出雹暴单体VIL的时间序列变化(单体VILmax出现时记为0,时间单位为雷达体扫间隔,约5~6 min),多数单体VIL的时间变化序列呈波动状态,基本都出现了1~2次跃增,单体VIL出现多次跃增,地面降雹往往也呈间歇性。遗憾的是新疆降雹实况时间记录多为10 min级甚至小时级,无法与5 min级的体扫对应比较,得出更具参考性的指标,但在伊犁河谷降雹过程中,都伴随着VIL的跃增。
图5 冰雹云单体VILmax时间序列变化单体(VILmax出现时间记为0)
4 冰雹个例分析
根据传统的分类方法[2,17]将对流风暴分为4类:普通单体风暴、多单体风暴、线风暴和超级单体风暴。
伊犁河谷冰雹云中最常见的为多单体风暴(表4),9 a共出现11次,其次是普通单体风暴,线风暴和超级单体风暴出现次数较少。普通单体风暴和多单体风暴均在昭苏县出现次数最多,值得注意的是霍城县是线风暴和超级单体风暴出现最多的地方,这可能是霍城县北部山脉对西南气流起辐合抬升作用,为对流风暴的发展提供了有利地形条件。
表4 4类对流风暴参数统计
对比4类对流风暴的回波强度发现,普通单体、多单体、线风暴和超级单体风暴的最强回波范围及平均值依次逐渐增强,这指示了伊犁河谷地区对流风暴的强度排序。但回波顶高却与之相反,超级单体回波顶高在4类中最低,普通单体和多单体的回波顶高相对较高,这与对流风暴越强回波顶高越高[18]相悖,以超级单体风暴为例,风暴成熟时0.5°仰角上观察到30 m/s左右的雷达径向速度,较强的水平风场可能抑制了对流在垂直方向的发展。3次超级单体风暴均伴随着低层钩状回波和有界弱回波区,其中2次出现旁瓣回波,但未发现三体散射,只有普通单体风暴出现过1次三体散射现象。线风暴低层易出现钩状回波和垂直回波悬垂。
4.1 普通单体风暴
2018年6月30日17—18时,伊犁州昭苏县昭苏镇、昭苏马场、洪纳海乡出现暴雨、冰雹等恶劣天气,降雹持续约20 min。造成此次冰雹的对流风暴为普通单体风暴。
17:10,>45 dBZ的回波生成于昭苏以西边境外,随后向西南方向发展。17:15,回波强度>50 dBZ,回波顶高达10 km,GR2软件开始识别出冰雹。随后回波继续加强发展,最强回波达59.5 dBZ、回波顶高达13.2 km,50 dBZ的强回波核高达9.4 km,且在垂直方向上反射率因子大值区基本重叠,没有倾斜,剖面图上表现为高悬的回波质心,无弱回波区或有界弱回波区。对应的径向速度图上存在明显的正负速度对,低层有明显辐合,但高层辐散不明显。19:21雹暴消亡,生命史约2 h。
此次雹暴还伴随着三体散射现象,17:41在2.4°和3.4°仰角雷达径向上均观测到三体散射长钉(TBSS),共维持了5个体扫约20 min,由此判断降雹最晚17:41开始。而VIL在17:31出现跃增,由13.1 kg/m2升至30.1 kg/m2,从VIL出现跃增到降雹,时间10 min。
4.2 多单体风暴
2021年6月6日15—17时,特克斯县阔克铁热克乡、乔拉克铁热克镇、喀拉托海镇出现暴雨冰雹等灾害天气。造成阔克铁热克乡、乔拉克铁热克镇降雹的雹云为普通单体风暴,其在东移北上过程中分裂发展,形成多单体风暴,继而造成喀拉托海镇降雨降雹。
14:20,特克斯西边界生成一对流单体A0(图6 a黄色箭头始端),其在北上东移的过程中逐渐增强,在14:38—15:58的15个体扫中,GR2软件均预警了大冰雹,回波强度均在50 dBZ以上,回波顶高在11 km以上,VIL出现了3次跃增,跃增量分别是13.6、8.6、9 kg/m2,预警指标也适用于此个例。约15:40前后,对流单体A0造成降雹,降雹时间与灾情记录时间一致。
15:46(图6b),A0的西侧出现了新单体,普通单体风暴开始发展为多单体风暴。16:03,A0开始减弱,分裂成A和B两个单体(图6c),并在西侧新生了单体C、D;随后多单体风暴整体东移,在风暴的西侧、西北侧不断有新的单体生成(图6d),这些单体分别增长、成熟、衰减。以G—F—C一条线上的3个单体为例做剖面(图6e),可清晰的看出多单体风暴的传播过程,G在发展并趋于成熟(通常最强烈的天气都由这个阶段的对流单体产生),F在减弱趋于衰亡,C已处于衰亡阶段。
图6 2.4°仰角反射率因子(a、b、c、d)及垂直剖面(e)
普通单体风暴在15:23前后发展为多单体风暴,此后不断有新的单体在多单体固定一侧生成,说明有利的天气条件一直在持续。到20时前后多单体风暴消亡,共历时4.5 h,生命史较长。
4.3 线风暴
2016年6月19日19—20时,霍城县出现强对流天气,其中尺度系统为线风暴,此线风暴由多单体风暴演变而来。18:55霍城县东部有3个孤立对流单体A、B、C(图7),A缓慢东南移并不断发展加强,B和C两侧不断有新的对流单体生成并互相靠拢加强(图7a);19:12这些分散的对流单体结合在一起成,呈断线型排列(图7b),初具线风暴特征,最大回波强度达56 dBZ。随后线风暴继续南下,D单体移速快,19:39与线风暴合并,在线风暴的东侧仍有新的单体生成。19:45线风暴发展到成熟阶段(图7d),长度约40 km,显著特征是回波整体呈弱“弓形”结构。
图7 2.4°仰角反射率因子(a,b,c)、垂直剖面(d)及径向速度(e)
从雷达不同仰角径向速度的演变来看,无论是高层还是低层,径向速度均以线风暴为界,南北两侧分别为正速度区和负速度区(以2.4°仰角为例,图e),即线风暴北侧为强西北风,南侧为弱东南风,线风暴生成于南北气流汇合处。但此线风暴的组织性较差,多个单体分布较零落,强回波顶高不超过5 km。20:00,飑线开始断裂成三段,后逐渐减弱消亡,生命史约1.5 h。
4.4 超级单体风暴
2015年7月5日霍城县多个乡镇出现短时强降水和冰雹等强对流天气(图8),此次风暴几乎具有超级单体风暴的所有特征:反射率因子图上有典型的钩状回波、前侧入流缺口,垂直剖面上有悬垂结构和有界弱回波区、速度图上有中气旋。
16:42,霍城县西北边境线附近出现46 dBZ回波,此后回波向东南传播发展,并不断与其南部的弱回波合并,17:25最强回波超过50 dBZ(图8a),GR2软件提示冰雹预警。对流风暴继续加强,形态由块状回波发展为弓形回波,并在低层反射率因子图上出现典型的钩状回波和弱回波区WER(图8b),风暴已经出现超级单体风暴特征。17:47最大回波中心达到61 dBZ,回波顶高达9.5 km,50 dBZ的强回波顶高为4.8 km。在反射率因子垂直剖面图上可以看到回波悬垂、有界弱回波区BWER以及BWER左侧的回波墙(图8d),速度图上明显的正负速度对为气旋式辐合(图8e),指示着有中气旋存在,此时风暴已发展为成熟的超级单体风暴。随后反射率因子继续加强,最大至72 dBZ,50 dBZ强回波核心增高至5.8 km,17:58,50 dBZ强回波底高已由17:47的802 m降至0,说明强回波核开始接地,正在产生大降水或冰雹。18:16风暴组织结构变得松散,逐渐消亡,整个超级单体生命史约2 h。
在17:42后的4个体扫中,低层均有旁瓣回波(图8b),这是目标物回波强度特别强的表现,伊犁河谷3次超级单体个例2次都出现过旁瓣回波。
图8 2.4°仰角反射率因子(a、b、c)、垂直剖面(d)及径向速度(e)
5 结论与讨论
通过分析2010—2021年伊犁河谷地区152次冰雹的时空分布特征及32个雹云单体的雷达回波特征得出以下结论:
(1)5—8月是伊犁河谷冰雹高发期,其中6月出现冰雹最多,且冰雹云在6月发展最为旺盛;16—19时是冰雹日变化高峰期,其中17时发生冰雹次数最多。
(2)冰雹云最强回波在49~72 dBZ,大部分冰雹的最强回波都在50 dBZ以上;回波顶高在6.8~13.3 km,99%的雹云顶高都在8 km以上;VILmax值跨度大,为5.02~47.7 kg/m2,VIL的跃增对冰雹预警有较好的指示意义。反射率因子达50 dBZ、回波顶高超过8 km、VIL出现明显跃增可作为伊犁河谷冰雹的预警指标。
(3)冰雹的空间分布与河谷喇叭口地形息息相关,南部沿山的昭苏县和特克斯县及北部沿山的霍城县、尼勒克县和霍尔果斯市是冰雹的多发地。昭苏县的冰雹占伊犁河谷地区冰雹总数的一半以上,它也是最常出现普通单体风暴和多单体风暴的地方,而线风暴和超级单体风暴多发在霍城县。伊犁河谷雹暴多为多单体,其次是普通单体,线风暴和超级单体出现较少,其平均回波强度由大到小依次是超级单体、线风暴、多单体、普通单体风暴。
(4)普通单体反射率因子大值区在垂直方向上基本重叠,没有倾斜;多单体一般出现在较为有利的天气条件下,不断有单体在同一侧新生、增长、成熟、衰减,会对同一地区反复影响,且生命史较长;伊犁河谷已出现的线风暴组织性都较差,结构松散,回波顶高不高,低层常出现钩状回波;已出现的超级单体基本具备所有典型雷达回波特征,但有些个例中气旋特征不显著。
伊犁河谷因其特殊地理地形,雹暴也具有其特殊性,在此观测到的超级单体风暴的回波顶高是4类雹暴中最低的,低层较大水平风速,对风暴在垂直方向的发展可能有抑制作用,这还需要更精细的探测资料来论证。本文所能收集到的降雹时间分辨率多为小时级或10 min级,对应5~6 min的雷达体扫时间过粗,雷达回波指标的时间精度有待进一步提高。