1. 引言
雷暴大风是一种由大气对流活动导致的灾害性天气，具有突发性强、持续时间短、致灾性强等特点，经常给农业、工业、养殖业等造成大量经济损失，对其的准确预报一直是短临预报工作的重难点。它有多种产生机制，主要因对流风暴的强下沉气流而造成，或者因为冷池密度流或者高空水平动量下传的作用而产生。研究发现，当对流层中层有干空气侵入时，下沉气流具有较大的不稳定性，同时对流层低层环境大气的温度直减率较大、中低层水平风垂直切变大均有利于雷暴大风的出现 [1] [2] [3] ；梁建宇等 [4] 、孙虎林 [5] 等认为雷暴大风的产生，高层冷空气侵入、低层增暖是大尺度条件，低层和地面辐合、干线是触发条件。随着我国多普勒雷达组网的建成，从中小尺度的角度分析雷暴大风逐渐成为研究重点，多从大风的雷达回波特征及其演变等角度入手。经概括，对于雷暴大风的雷达识别指标主要包括：对流组织类型、中层径向辐合、风暴移动速度、回波质心下降、低仰角风速大值区和辐散等 [6] [7] [8] 。
近年来重庆铜梁也多次出现雷暴大风天气。本文将利用2012~2021年的雷达资料配合大风的实况观测资料，对铜梁地区发生的雷暴大风天气的时空特征进行统计分析。通过对其天气现象发生时雷达反射率、回波顶高、累积液态水含量等指标进行分析，提炼出适应于本地化的雷暴大风天气发生发展的雷达特征产品的指标范围，有效提高大风天气监测预报预警服务水平和气象防灾减灾能力。
2. 资料与方法
2.1. 雷暴大风个例的筛选
雷暴大风是指，某天(即当日20时至翌日20时)既观测到瞬时风力 ≥ 17.2 m/s的大风又观测到雷暴，则计为一个雷雨大风日 [9] 。根据对2012~2021年间的雷暴大风个例进行统计分析，近10年铜梁本站出现17.2 m/s以上的雷暴大风个例有18个，均发生在汛期(4~9月) (如表1)。
Table 1. Screening results of thunderstorm gale processes in Tongliang from 2012 to 2021
表1. 2012~2021年铜梁雷暴大风过程筛选结果
2.2. 雷达特征产品的选取
考虑强对流生成的三要素(层结不稳定、水汽、抬升触发机制)，因而在选取雷达特征产品的时候也应围绕着对流单体发展强度(回波顶高ET、回波组合反射率CR)、水汽条件(垂直液态水含量VIL)来重点分析与总结。
3. 雷暴大风的回波特征分析
Figure 1. Monthly statistics of thunderstorm gale in Tongliang from 2012 to 2021
图1. 2012~2021年铜梁雷暴大风逐月统计图
雷暴大风在8月份出现次数最多，占比36%，其次为5月，占比分别为27%，从逐月的雷暴大风统计结果也可以看到(图1)，近十年里8月是汛期发生雷暴大风频率最高的，其次是5月、4月，6、7月发生雷暴大风次数较少，9月份几乎不发生。4月、5月作为春夏转换的过渡时期，较易发生雷暴大风比较正常。反而在副高稳定控制重庆地区的盛夏(8月)，是雷暴大风的高发期。实际业务工作中又很难及时准确的把握住完整清晰的出流边界，因而这类型的雷暴大风天气存在着一定的预警难度。
3.1. 回波反射率的统计分析
由于每次雷暴大风过程对应的雷暴单体可能是多个同时发展或单一强雷暴单体长时间维持，因而对于雷暴大风过程最强反射率中心的分析也应从单一雷暴单体出发，首先对从历史实时资料中筛选得到的18个雷暴大风个例进行相应的组合反射率产品特征统计分析。
通过分析所有过程发生时段及前后的雷达回波的组合反射率图像，对18个雷暴大风天气过程的反射率形态进行了分类，类别主要包括孤立块状、零散块状、带状、片状及多形态五种(图2)，孤立块状指反射率图像呈现单一的块状单体分布，零散块状指图像呈多个独立的块状单体相对独立的分布，带状指图像呈现多个单体且生成与移动方向一致的线状分布，片状指图像由大片的低反射率回波伴随较强反射率中心的回波分布，多形态指图像由以上两种或多种的组合或者是存在有一种到另一种的演变而成。
Figure 2. Characteristic morphology statistics of the combined reflectivity of thunderstorm gale in Tongliang from 2012 to 2021
图2. 2012~2021年铜梁雷暴大风组合反射率CR的特征形态统计
分析发现近十年雷暴大风个例的回波形态分布较为平均，几种回波形态的个例均有发生，且发生频率也较为接近。大风发生持续时间相对更长的过程往往伴随着多个单体组合而成的多形态特征分布，以零散块状–带状组合更为常见。同时，从组合反射率最强反射率因子的强度统计结果来看(表2所示)，强度表现的更强的大多在零散块状、孤立块状这两类形态中，因为有超级单体的发生发展使得相应的最强回波反射率也更强，能发展到65 dbz以上的强度，而带状、片状及多形态的多单体回波形态则表现出中等强度特征，大多发展到55 dbz左右的强度。从所有统计结果来看，CR ≥ 50 dbz的最大组合反射率因子强度是雷暴大风(8级)产生的预警阈值，通常产生8级大风的最大组合反射率因子范围在50~55 dBz之间，而产生极端大风(9级以上大风)的最大组合反射率因子阈值为≥60 dbz。
Table 2. Characteristic values of combined reflectivity during the thunderstorm gale processes in Tongliang from 2012 to 2021
表2. 2012~2021年铜梁雷暴大风过程组合反射率因子特征值
3.2. 回波顶高ET的统计分析
回波顶高是在≥18 dbz反射率因子被探测到时，显示以最高仰角为基础的回波顶高度。因而常用来对最高顶定位，判识风暴体。
从表3给出的各个过程发生时段的最高回波顶高可见，对应雷暴大风发生的回波单体的发展高度整体还是较高，ET值 ≥ 12 km可以视为一个理论阈值，且对流发展最强阶段很多单体能发展到14 km以上的高度。经对每种回波形态对应的过程中最高回波顶高，发现以下特征：片状回波的回波顶相对来说较低，ET值多在12 km以内，块状、带状回波ET值偏高，多在14 km以上，个例中最大发展到了21 km。
同时考虑对雷暴大风过程时刻前后的回波顶高度也进行了统计，以期能够对雷暴大风开始及结束的一个相应特征值进行总结，以提升雷暴大风预报的提前量及准确率。经过统计分析，将雷暴大风过程回波顶高发生前、中、后的时序变化分为五个类别：1、逐渐增高；2、逐渐降低；3、先减再增；4、先增再减；5、稳定变化。经分析发现，从雷暴大风发生前、中、后期保持同一变化趋势特征的回波顶高，即逐渐增加及逐渐减少的回波顶高对应的大风过程表现出的最显著特征为过程持续时间较短，大风只出现并维持在两到五个体扫(12~30分钟左右)的时间。而带有趋势变化的先增后减及先减后增的类型则对应了较长的大风过程持续时间(维持六个体扫以上)。
Table 3. Characteristic values of ET during the thunderstorm gale processes in Tongliang from 2012 to 2021
表3. 2012~2021年铜梁雷暴大风过程回波顶高特征值
3.3. 垂直积分液态水含量VIL的统计分析
东高红等 [10] 利用天津雷达垂直积分液态水含量(VIL)资料的分析结果中指出VIL值达到30 kg/m2是地面灾害性大风出现的阈值，而当VIL值达到40 kg/m2可视作大风的预警指标，地面灾害性大风的出现往往伴随VIL值达到最大后的快速减小，当再次跃增时就是大风开始的标志。为了验证这一结论在铜梁的适应性，重点对雷暴大风过程开始前及发生初期的几个时次进行统计分析。
Table 4. Characteristic values of VIL during the thunderstorm gale processes in Tongliang from 2012 to 2021 (unit: g/m2)
表4. 2012~2021年铜梁雷暴大风过程垂直液态水含量特征值(单位：g/m2)
表4给出了重庆铜梁近十年的典型雷暴大风过程发生前的4个体扫以及大风发生时段的前两个体扫对应的垂直液态水含量值统计情况，从现有个例的VIL值来看，大多数个例表现出随大风时刻临近VIL值逐渐增加的趋势(4/18)，部分个例伴随着VIL值逐渐减小(2/18)，部分个例伴随着VIL值的稳定维持(10/18)，能满足文献中提到的先达到最大值再快速减小，并再次跃增的个例只有2次。
同时，文献中提到的灾害性大风出现的阈值为30 kg/m2，而结合我们现有的个例资料来看，铜梁雷暴大风的阈值应该低于该标准，VIL值15 kg/m2是本地区雷暴大风过程可能发生的预警阈值(参考所有雷暴大风过程VIL平均值与众数)。
4. 结论与分析
本文利用铜梁地区2012~2021年的雷达资料配合大风的实况观测资料，根据极端强天气的定义，筛选出18例雷暴大风天气个例，并详细统计分析各个个例的雷达回波特征，得到以下结论：
1) 铜梁的雷暴大风过程基本发生在4~9月。并且8月是汛期发生雷暴大风频率最高的，其次是4月、5月，6、7月发生雷暴大风次数较少，9月份几乎不发生。
2) 大风发生持续时间更长的过程伴随着多个单体组合而成的多形态特征分布，强度表现得更强的大多在片状、孤立块状这两类形态中。CR ≥ 50 dbz的最大组合反射率因子强度是雷暴大风(8级)产生的阈值，CR ≥ 60 dbz是产生9级以上雷暴大风的阈值。
3) ET值 ≥ 12 km可以视为雷暴大风发生的理论阈值，且对流发展最强阶段很多单体能发展到14 km以上的高度。逐渐增加及逐渐减少的回波顶高对应的大风过程表现出的最显著特征为过程持续时间较短；先增后减及先减后增的类型则对应了较长的大风过程持续时间。
4) VIL值15 kg/m2是本地区雷暴大风过程可能发生的预警阈值。
基金项目
重庆市铜梁区科技计划项目(2021015)。
NOTES
*第一作者。
#通讯作者。