短时暴雨会引起山洪、地质灾害、城市内涝等次生灾害,一直是强对流天气预报中的关注重点。短时暴雨具有突发性强、降水强度大、局地性明显的特点,其发生发展的机理复杂,常规监测手段难以捕捉。随着探测技术的迅速发展,应用多源新型探测资料可对强对流天气发生的潜势提前做出估计和判断,有助于短临预报水平的提高。有学者已将多源新型探测资料应用于强对流天气监测预警中
[1-3]。
微波辐射计可实时、连续遥测地面到10 km高空的大气温度、水汽、湿度、液态水分布等,提供高频温湿廓线演变,对暴雨天气分析和业务预报具有重要的指示意义。目前国内许多学者先后用多通道的微波辐射计研究大气中气态、液态水的分布及其演变规律,取得了成果
[4-6]。高宇星等
[7]利用微波辐射计反演的液态水含量分析,发现在降水时段液态水含量呈上升趋势,与降水强度有很好的对应关系,并指出液态水含量>0.14 g·m
-3可作为降水临近的参考指标。冼星河等
[8]使用微波辐射计资料分析大暴雨过程,指出微波辐射计反演的水汽含量对降雨有指示作用,提前20 min的数值增长具有参考性。
毫米波云雷达通过电磁波探测大气中的云雨粒子,提供精细的云雨廓线结构,有助于分析降水过程中的云物理及演变特征。傅新姝等
[9]利用云雷达分析,发现云雷达能清晰地反映强降水过程中降水云系统的云层高度、厚度以及融化层高度等云结构及其变化特征,对认识强降水过程中的降水云系统演变具有重要参考意义。肖婷等
[10]利用云雷达分析,发现在降水发生前1 h云雷达已经观测到高层云层的发展,相比于天气雷达有更多提前量,对于强降水预报有指示意义。
S波段双偏振天气雷达与传统单偏振雷达相比,增加了多个新型的雷达观测变量,包括差分反射率因子(
ZDR)、差分相移率(
KDP)和相关系数(
CC )等双偏振量,有利于进一步理解降水系统的微物理特征。申高航等
[11]利用多源资料研究了台风雨带强降水过程中降水粒子结构变化的微物理机制,利用双偏振特征参量解释了降水较强时段的冷云降水机制。袁媛等
[12]基于双偏振雷达观测资料对短时强降水监测预警开展了研究工作,指出
KDP与小时雨强有显著正相关。
选取2023年8月6日金华地区一次局地强对流天气过程作为研究对象,对新型探测资料在短临预报预警中的作用进行分析,尝试提炼新型探测设备资料在短时暴雨过程中的应用指标,为短临预报预警提供参考。
1 多源新型探测设备介绍
利用多源新型探测资料对一次局地短时暴雨进行分析,所涉及的探测设备信息如
表1所示。
2 降水实况
2023年8月6日21:00—22:00金华出现一次局地短时暴雨天气过程(
图1a),暴雨主要集中在金华市区,有5站降水量达到50 mm以上,最大累计降水量为金华高铁站(77.7 mm);金华国家气象站位于降水中心区域,累计降水量达55.9 mm。分钟降水资料(
图1b)显示,金华国家气象站21:25开始降水,降水集中期(21:25—22:05)降水效率高,分钟降水量在1 mm以上,最大达2.6 mm。
3 天气背景
2023年8月6日20时500 hPa湘鄂交界附近为大陆高压,东海东部130°E附近为台风“卡努”,金华地区处于大陆高压东侧和“卡努”西侧的偏北气流中(
图2a)。500 hPa存在<0 ℃的冷舌和温度露点差>7 ℃的干舌,850 hPa存在>23 ℃的暖舌和露点温度差<4 ℃的湿舌,上冷下暖的温度层结和上干下湿的湿度层结,为对流性降水的出现提供了有利的不稳定条件;对流层低层湿空气接近饱和,为对流性降水提供了水汽条件;925 hPa金衢盆地有一个显著的中-α尺度气旋性环流中心,有利于次天气尺度上升运动的产生和维持,为对流性强降水出现提供了必要的垂直上升运动条件。
从21时地面图(
图2b)可知,上游地区(金华西北部)由于前期降水而形成冷池,北侧相对冷的空气与金华相对暖的空气对峙,等温线密集,形成中-β尺度冷锋,其南下过程中有利于触发中小尺度对流;金华站附近风场中存在一个中-β尺度的气旋性闭合环流,有助于中小尺度上升运动的建立,为短时暴雨的形成提供有利的辐合抬升条件。地面中-β尺度冷锋和气旋性辐合中心为局地短时暴雨的出现提供了有利的动力和热力条件。
4 多源新型探测资料分析
4.1 <bold>S</bold>波段双偏振天气雷达
从0.5°仰角反射率因子演变可知,8月6日21时对流单体生成并发展,22时以后减弱,22:30消亡。21—22时>45 dBZ强回波区自北向南移动,且强回波区中始终镶嵌有>55 dBZ的回波块。受强回波影响,位于雷达偏南方向10~15 km的5个气象站出现>50 mm·h
-1的短时暴雨(
图1a)。从0.5°仰角差分反射率因子
ZDR和差分相移率
KDP演变可知,强降水区对应
ZDR、
KDP和
CC大值区,说明以浓度较高、直径较大的液态降水粒子为主。
由
图3a可知,21:43(处于降水集中期)0.5°仰角的反射率因子PPI中>45 dBZ的强回波呈挂件状形态,0.5°仰角的径向速度PPI图显示存在明显的辐合区(
图3b)。由反射率因子RHI(
图3c)可知,风暴发展较为旺盛,>45 dBZ的强回波伸展至0 ℃层附近(5.8 km左右),在0 ℃层以下存在弱回波区(WER)和弱悬垂结构,强回波中心值达到57.5 dBZ,根据
Z-R关系,反射率因子与雨强呈非线性正相关,高的反射率因子可侧面反映大的雨强。由径向速度RHI(
图3d)可知,向着雷达的入流支和沿着雷达径向向外的出流支,在距离雷达10~15 km内交汇,产生辐合上升运动,上升气流将对流层低层的大量水汽源源不断输送到高空并凝结,造成此次短时暴雨事件。分析双偏振产品可以发现,
ZDR剖面(
图3e)在弱回波区上方存在
ZDR大值区(
ZDR≥1 dB)即
ZDR柱,
ZDR柱顶部延伸至6 km左右,说明上升气流强盛
[13-14],
ZDR最大值为4.18 dB;由
KDP剖面(
图3f)可知,弱回波区上方存在
KDP大值区(
KDP≥1 (°)/km)即
KDP柱,发展高度接近0 ℃层,说明降水强度强
[13-14],
KDP最大值达4.5(°)/km。
综上,此次短时暴雨天气过程由一个生消速度快、维持时间短的对流单体造成。WER以上存在较为宽广的ZDR柱和KDP柱,上升气流强盛,水汽供应充足,水汽凝结形成深厚降水云,云体以暖云层为主导,大量的云滴经历碰撞合并等微观过程,最终形成浓度较高、直径较大的液态降水粒子并降落,产生短时暴雨。在短临预报预警中,ZDR和KDP迅速增大并形成伸展高度接近0 ℃层的ZDR柱和KDP柱,表征强的上升气流和降水强度,可为短时暴雨的预警提供重要参考。
为检验S波段雷达定量降水估测的效果,对比21:00—22:00基于HSR混合扫描
Z-R算法的OHP产品(单偏振产品)、基于CSU-HIDRO算法的OHP产品(双偏振产品)定量降水估测效果。前者是基于反射率因子得到单偏振降水产品,即78号产品。后者CSU-HIDRO算法是基于偏振量进行相态识别,分液态、固态和混合降水3个类别,联合使用
ZH、
ZDR、
KDP变量进行降水估测,即169号产品。由
图4可知,169号产品存在>50 mm的暴雨中心,中心的估测量达75 mm,接近金华高铁站实测小时雨量(77.7 mm);78号产品降水估测量基本<50 mm,存在明显的低估问题。原因是169号产品通过双偏振参量获取丰富的微观粒子特征,用不同的雷达测雨公式计算出更合理的雨量;78号产品使用的
Z-R关系式为固定的
Z=300
R1.4,对于不同的降水类型(层状云降水或对流云降水),未使用不同的参数值。
4.2 毫米波云雷达
由
图5a可知,21:05云体回波顶高升至12 km附近,积雨云发展旺盛,柱状回波明显,此时高层反射率因子大于低层;21:05—21:25由于降水粒子下落,反射率因子大值中心迅速下传,云体中下部回波强度逐渐增强,反射率因子中心值达32 dBZ,该时段还未出现降水;21:25以后,30 dBZ左右的反射率因子高值区集中在1 km以下,持续时间约40 min,该时段出现较强的降水,基本维持在50 mm·h
-1以上的降水强度。此外,21:25—22:05距地1 km高度以上的反射率因子出现强衰减现象,是因为强降水粒子对毫米波云雷达信号造成非常严重的衰减,导致反射率因子在1 km以上区域表现为空值。
21:05—21:25毫米波云雷达径向速度(
图5b)显示,6 km以上的高层云大部分区域的径向速度为正,说明云体内气流有显著的上升运动,速度最大为10 m·s
-1,金华站上空的水汽被输送到高层大气之后发生凝结,凝结过程中向周围环境释放热量,空气吸收热量形成上升运动;云体中下层垂直速度以负值为主,为-5~0 m·s
-1,表明云体内气流以下沉运动为主,这是因为云体内已经发生降水,粒子降落产生拖曳作用而形成下沉气流。降水集中期(21:25—22:05)距地面0~1 km为一致的下沉运动区,最大下沉速度为-10 m·s
-1。22:15以后降雨渐止,径向速度图谱逐渐由负值演变为正值。
速度谱宽是可以表征云系中粒子径向速度间存在的差异以及大气湍流差别的物理量。由
图5c可知,降水集中期(21:25—22:05),云系中因为上升气流和降水拖曳作用下沉气流的联合作用,导致降水粒子的下降速度差异较大,使云系内粒子的速度谱宽出现大值区,最大值为11 m·s
-1。22:15以后云体内速度谱宽分布较为均匀,都在0.5 m·s
-1以内,说明云体内降水粒子径向速度差异不大、粒子大小差异也不大,并且没有明显的风切变或者湍流,云体内已无明显对流。
可见,毫米波云雷达在对流云的垂直结构观测中具有明显优势,不仅能够获得精细化的云层结构和动态变化信息,还可以获取云系内部粒子微观的垂直运动特征。本次过程毫米波云雷达探测到云体的回波顶高迅速伸展至12 km、反射率因子迅速增强至30 dBZ以上、高层径向速度迅速增加至10 m·s-1,对强降水有一定的指示意义,借助指标可提前约20 min做出短时暴雨的预报预警。
4.3 微波辐射计
综合水汽含量(IWV)和液态水路径(LWP)是微波辐射计常用于分析大气水汽条件的产品。IWV表示单位面积上的水汽质量(单位kg·m-2),可直观反映大气中水汽含量水平。大气中的水以水汽为主,水汽本身无法形成降水,只有云中的液态水才有可能通过自然过程或人工影响产生降水,液态水含量的高低直接决定降水强度;液态水路径也称为积分液态水含量(单位g·m-2),表示单位面积垂直气柱内积分的液态水含量,可以直观反映大气柱中液态水含量。
由
图6a可知,对流云未发展时,IWV稳定在60 kg·m
-2附近。降水发生前20 min(21:05—21:25)IWV增幅较大,由60 kg·m
-2迅速增长至90 kg·m
-2,表明大气水汽含量发生明显的波动,对降水发生有较好的指示意义。降水集中期(21:25—22:05),IWV一直稳定维持高值状态,在100 kg·m
-2附近。22:05以后降水减弱,IWV迅速降回60 kg·m
-2附近。
LWP与本次降水过程对应较好(
图6b),降水发生前(21:05—21:25) LWP增幅较大,由0 g·m
-2迅速增至10 g·m
-2,说明云体中的液态水不断增多;LWP在强降水集中期(21:25—22:05)稳定维持在12 g·m
-2附近的高值状态;22:05以后降水趋于减弱,LWP迅速降低。
微波辐射计IWV和LWP产品的迅速增大,预示着探测站点上空水汽迅速聚集并且通过微物理作用形成大量雨滴,出现降水的可能性非常大,可为监测预报预警工作提供有利的参考指标。从本次降水过程的IWV和LWP产品演变规律来看,IWV骤增至90 kg·m-2、LWP骤增至10 g·m-2,对出现短时暴雨具有较好的指示意义。
4.4 激光气溶胶雷达
激光气溶胶雷达的消光系数表示气溶胶粒子通过散射和吸收作用,对某一特定波长光辐射在单位面积、单位长度的衰减过程。消光系数接近0时,没有污染或者污染很轻;消光系数介于1~2时,通常是水汽与污染的混合造成;当消光系数超过2时,主要是水汽造成,因为水汽的消光系数比PM2.5等污染物大得多。
激光气溶胶雷达的退偏振比对判断气溶胶粒子非球形程度有一定的指示意义,其数值越大表明非球形程度越高,是区分气溶胶粒子类型的一个有效指标。水汽一般是较为规则的球形粒子,退偏振比数值很小;雨滴在降落过程中会受到空气阻力、碰撞等外力作用,呈现为不规则的非球形粒子,出现降水时激光气溶胶雷达会呈现出高退偏的特征。
21:25激光气溶胶雷达垂直消光廓线在2 km以下消光系数迅速增大,尤其是0.5 km以下的近地面,消光系数一般在2以上,在图谱上表现为近地面红色高消光区,22时开始减弱(
图7a)。结合退偏振比图谱(
图7b)可知,21:25—22:05的退偏振比较大,分布在0.3~0.5,这是由于大量的水汽聚集之后,空气中的水汽饱和度明显增加且出现了大量的液态水,雨滴在降落中,激光气溶胶雷达退偏振比呈现出一个大值区。实况降水(
图7c)与激光气溶胶雷达的图谱特征对应较好。
本次降水天气过程对应“高消光、高退偏”的图谱特征,消光系数迅速增大到2以上,指示有大量水汽聚集,同时结合退偏振比迅速增大到0.3以上,表明出现了大量的不规则雨滴。
5 结论
运用高时空分辨率的多源新型探测资料和常规气象观测资料,从大气热动力结构、降水云垂直结构等方面对金华地区一次局地短时暴雨过程进行分析,得出以下结论:
(1)上干冷、下暖湿的温湿层结为暴雨提供了不稳定条件,边界层中-α尺度辐合系统、地面中-β尺度冷锋和气旋性涡旋起到动力和热力触发作用。
(2)S波段天气雷达双偏振参量能较好地反映降水微物理特征,ZDR和KDP迅速增大并形成伸展高度接近0 ℃层的ZDR柱和KDP柱,表征强的上升气流和降水强度,可为短时暴雨的预警提供重要参考。
(3)毫米波云雷达能较好地观测云体的垂直廓线信息,降水发生前20 min探测到高层云开始发展,对短时暴雨有一定指示意义;微波辐射计IWV迅速增大表明探测站点上空的水汽迅速聚集,LWP迅速增大预示着水汽通过微物理作用形成大量雨滴,依据IWV和LWP产品演变特征,提前约20 min做出短时暴雨预警;激光气溶胶雷达消光系数>2、退偏振比>0.3的“高消光、高退偏”图谱特征能较好地指示短时暴雨。
本文只是针对一次短时暴雨天气过程的分析,今后需要收集更多典型的短时暴雨个例进行对比研究,归纳总结出多源新型探测资料有利于形成短时暴雨的特征阈值。