引言

差分反射率(Z_DR)柱是雷达观测中一种垂直柱状回波结构，表现为Z_DR高值区在对流云内呈柱状延伸至0℃层以上(Illingworth et al, 1987)。暖层中连续的Z_DR高值区通常与大雨滴的存在相关(Kumjian and Ryzhkov, 2008; Snyder et al, 2015)，并可用于区分雨滴和冰相粒子(如霰或冰雹)(Bringi et al, 1991)，而Z_DR高值区为何能延伸到0℃层以上，成为自应用双偏振天气雷达后较早被关注的现象(Hall et al, 1980; 1984;Conway and Zrnić, 1993)。飞机观测、地面雷达探测(Brandes et al, 1995)及数值模拟(Tuttle et al, 1989; Kumjian et al, 2014; Ilotoviz et al, 2016; van Lier-Walqui et al, 2016)等和Z_DR柱相关的研究已揭示，强对流上升气流的输送作用和过冷层的液滴碰并促进了过冷大雨滴和对应的Z_DR柱的形成。

Z_DR柱可作为识别强对流天气、预警冰雹和大风的重要雷达观测特征(Snyder et al, 2015)。其形态、高度和面积等参数与强上升气流及冰雹生成密切相关(Ilotoviz et al, 2018; Snyder et al, 2017)。Z_DR柱的延伸高度和面积越大，通常预示更强的对流活动和冰雹潜力(潘佳文等, 2023)。Z_DR柱不仅在强对流单体早期随主要的上升气流出现，在对流中期也可随地面降水出流强迫而次生(Ren et al, 2024)，有助于揭示对流单体内部动力结构的复杂性。近年来的研究还发现高炮防雹作业后Z_DR柱会在短时间内消失(孙跃等, 2023)。这些研究表明Z_DR柱在强对流天气的监测、预警和人工防雹等方面均具有重要的应用价值。

传统机械扫描式天气雷达扫描速度慢，且不同时刻扫描的仰角不同，可能使观测到的对流单体垂直结构出现错位，从而导致无法正确捕捉到Z_DR柱。相比之下，相控阵天气雷达采用方位机械扫描、俯仰相扫方式，具有高时空分辨率扫描能力，可精确捕捉云体内部结构与演变(Wu and Liu, 2014; Palmer et al, 2022; Kollias et al, 2022)，开展强对流天气监测预警研究(Kuster et al, 2020; McKeown et al, 2020)。例如Kuster et al(2019)使用快速更新的雷达数据(体扫时间≤2.3 min)，基于美国俄克拉何马州42次强对流个例分析了Z_DR柱的面积和深度等特征，指出这些Z_DR柱特征相对于-20℃层反射率因子核心出现时间均有不同程度的提前，可为冰雹、大风等强对流预警争取更多时间。

尽管我国近年来陆续部署一批相控阵天气雷达(Zhao et al, 2024)，但是，利用Z_DR柱特征进行冰雹预警方面还存在许多挑战。一方面，类似经验(Kuster et al, 2019)不容易直接移植，因为不同地区的天气、气候背景存在诸多差异(潘佳文等，2020；高丽等，2021；阮悦等，2022；汤兴芝等，2022；王晓君和郑媛媛，2022；龚佃利等，2024；王福侠等，2024；夏凡等，2024；杨湘婧等，2025；余文梅等，2025)，加之不同波段和型号雷达的数据本身存在差异，观测到的Z_DR柱特征也各不相同(王清川等, 2022；刘瑾等，2022；沈雨等，2023；李芳等，2023；郭飞燕等，2023)。另一方面，目前我国不同地区相控阵天气雷达探测冰雹云的个例大都尚在积累阶段，有必要着眼于部署时间相对长、范围广、个例积累相对较多的传统机械扫描天气雷达开展统计研究。

为充分利用强对流云的雷达探测信息，推动冰雹云早期识别和预警研究，本文基于山东近年来的三次局地性强、发展迅速的降雹个例，提出能够诊断Z_DR柱特征的分析方法，并结合不同扫描体制和型号的天气雷达，开展雷达数据特征和冰雹预警研究。

1 数据与方法 1.1 降雹个例与雷达数据

三次降雹个例和使用的雷达见表 1和图 1。三次过程500 hPa上山东均位于西北气流中，低层均受西太平洋副热带高压西侧偏南风控制。文中所用时间均为北京时。

对于2022年6月13日山东滨州的降雹个例(以下简称个例1)，使用滨州X波段双偏振相控阵雷达(X-PAR)观测数据进行分析。雷达站点海拔高度22 m，采用方位机械扫描、俯仰相扫的体制，每约123 s生成一个体扫数据，含有24个仰角(0.9°~42.3°，采样分辨率1.8°)，每个仰角下含400个方位(采样分辨率0.9°)，径向距离分辨率为30 m，最大探测距离为57.96 km。

对于2024年5月18日发生在山东临沂费县的降雹个例(以下简称个例2)，使用济宁S波段业务化天气雷达数据进行分析。该雷达为CINRAD-SA型双偏振雷达，海拔高度168 m，监测降水云时采用VCP21模式进行体扫，在约340 s内于9个仰角(0.5°~19.5°)下进行方位周扫，径向距离分辨率为250 m，最大探测距离为525.75 km。

对于2024年5月18日发生在山东济南的降雹个例(以下简称个例3)，除使用济宁S波段雷达外，还使用距降雹点更近的泰山C波段业务化天气雷达数据，两部雷达直线距离约为120 km。该雷达为CINRAD-CD型双偏振雷达，海拔高度1546 m，在约354 s内于14个仰角(0.5°~19.5°)下进行方位周扫，径向距离分辨率为150 m，最大探测距离为457.05 km。利用上述S波段和C波段雷达分别对同一降雹个例进行观测，也有助于对这两个常见的业务化天气雷达在冰雹预警能力方面的一致性和适用性进行对比。

1.2 数据处理和分析方法 1.2.1 三维直角坐标插值

将雷达数据(基于球坐标)插值到三维直角坐标格点。采用三线性插值方法估算目标格点上的水平反射率(Z_H)和差分反射率(Z_DR)(王超，2023；Ren et al, 2024)。对不同型号雷达的插值格点垂直分辨率均设为0.1 km。为了兼顾X-PAR的精细化观测能力，将其插值出的水平格点分辨率设为0.1 km。而对于S波段和C波段雷达，将水平分辨率设为1 km。这主要是因为本文使用的S波段和C波段雷达的波束宽度和方位角采样宽度均为1°，在100 km距离上已对应超过1.74 km的雷达波束切向距离，设定比1 km更精细的插值水平分辨率无实际意义。在插值前，考虑到X-PAR的统计波动较大，对本文所使用的Z_H和Z_DR沿雷达径向进行500 m窗口的有限脉冲响应(FIR)低通滤波(Ren et al，2024)，并剔除最低两层仰角层中交叉极化相关系数小于0.7的点，屏蔽人工确定为经常性遮挡的径向。对C波段进行1.2 km窗口的FIR低通滤波。涉及数据处理的其他问题见1.2.4节。

1.2.2 冷层组合Z_DR

冷层组合Z_DR(C_ZDRC)是每个水平格点上空0℃层以上高度Z_DR的最大值。其计算方法类似组合反射率(每个水平格点上空Z_H的最大值，以下简称Z_e)。本文三个降雹个例中的0℃层高度由相近时次的FNL分析数据和ERA5再分析数据综合获得。

C_ZDRC可直观反映Z_DR柱的水平位置和强度(Ren et al，2024)。在传统观点中，求不同高度上Z_DR的最大值并无意义，因为云顶存在水平冰晶、雪粒子的区域和云体边缘弱回波中水平和垂直偏振功率不平衡的区域都常存在一些Z_DR的大值，会使得组合Z_DR图中布满无诊断意义的大值点。但X-PAR由于功率小、灵敏度低，其能观测到的回波主要为Z_H>20 dBz的回波，恰好很难采集到云体边缘和云顶等易产生Z_DR大值的弱回波区的数据，由此计算的C_ZDRC上的大值中心便可直观反映Z_DR柱特征。利用C_ZDRC分析Z_DR柱的位置和强度特征，一方面可以避免对不同仰角层数据进行逐层查看、找到疑似Z_DR柱特征又由于高度不同而不易确定温度的局限性，另一方面对由于回波移动而使回波垂直结构在机械扫描天气雷达体扫数据中存在错位的情况也具有一定的兼容性，可避免人工垂直剖面查看时出现遗漏。

基于上述在X-PAR上利用C_ZDRC诊断Z_DR柱的经验，对本文所使用的X-PAR、S波段和C波段雷达，均屏蔽Z_H < 20 dBz的回波后再计算C_ZDRC。对这种屏蔽做法的必要性和影响的讨论见2.4节。

1.2.3 回波面积的时间-高度演变

对Z_H和Z_DR大于一定阈值的格点面积进行统计，形成时间-高度演变，可以弥补C_ZDRC图上只有Z_DR柱水平位置特征而没有垂直结构诊断特征的局限性。本文对Z_H分别取35、45、55 dBz，对Z_DR分别取1、2、3 dB作为阈值进行统计。其中Z_H>35 dBz仅用于表征明显降水回波的出现。屏蔽某个时间和高度上只有1个计数的点，以减少偶发性的杂波或非Z_DR柱内水平取向粒子的干扰。选取回波面积在时间和高度上相对连续出现的最早时刻在图中进行标示，用于讨论Z_H强回波和Z_DR柱特征出现的先后顺序。

对于回波面积时间-高度演变的空间统计范围，为了针对性地研究降雹云的变化特征，以每个个例中雷达监测范围内的首个降雹报告地点为准，屏蔽周边其他单体的回波，选取对应的降雹云的活动区域进行统计。具体做法为：个例1屏蔽37.9°N以北、118°E以东的数据；个例2屏蔽35.5°N以北、117°E以西的数据；个例3屏蔽36.5°N以南、117.5°E以东的数据。关于选定统计区域所造成影响的讨论见2.4节。

1.2.4 质量控制的其他问题

除1.2.1节所述数据处理外，本文不再进行其他杂波滤除操作。这主要是为了避免过度滤波或数据剔除影响到原本可见的Z_DR柱特征。按照1.2.1节的插值方法，只有在上下相邻两个仰角同时具有回波时，才能获得目标点的插值结果，故大部分偶发性杂波和主要集中在第一层仰角中的边界层杂波等非降水性回波将在插值过程中自动被剔除。

本文对所使用的X波段和C波段雷达数据不进行衰减订正操作，主要有以下三个原因：一是本文着眼于冰雹的早期预警，即尚未形成明显的强对流云图时的冰雹预警，此时雷达信号衰减十分有限；二是衰减订正操作难免会引入一些偏差，例如过度订正或波束间额外的条纹状偏差，这些会干扰到本文对Z_DR柱特征的识别；三是滨州X-PAR基数据中储存了衰减订正后的Z_H和Z_DR，但没有储存原始的Z_DR，这使得重新开展Z_DR的衰减订正工作存在一定困难。

2 结果分析 2.1 2022年6月13日滨州降雹个例

在滨州X-PAR监测范围内首个降雹时刻前约31 min(11:43)、降雹点以西可见20 dBz以上的单体云初生(图 2)。此时C_ZDRC图显示有许多分散的1 dB左右的数据点，很难直观反映是否存在Z_DR柱。在降雹前约25 min(11:49)，降雹点西北侧的Z_e核心刚超过40 dBz，但C_ZDRC图上已能看到明显的Z_DR柱特征。在回波核心的偏后侧(西北侧)存在3 dB左右的C_ZDRC大值中心，表明该处存在强上升气流和过冷雨滴。在首个降雹时刻(12:14)，降雹点附近的云团已发展为多个Z_e核心在55 dBz以上的单体。同时，C_ZDRC图显示此时存在多个Z_DR柱，其中一部分是围绕降雹点的，但降雹点上空不具备Z_DR柱特征。这种Z_DR柱水平分布形态和超级单体是相反的，但是与Ren et al(2024)归纳的模型类似，即在强上升气流和过冷雨滴区的对流层中高层下风方向可产生霰和雹，相关差异将在2.6节展开讨论。综上，本个例显示，C_ZDRC图中出现的Z_DR柱特征虽然不能直接提示降雹点，但在回波发展的早期能够提示对流的增强。

回波面积的时间-高度演变(图 3)显示，Z_H大于35 dBz、45 dBz的回波分别在降雹前27 min和23 min出现(图 3a, 3b)。虽然45 dBz是一个监测强对流的阈值，但是，以往山东观测个例中某些融化层亮带也会出现≥45 dBz的回波(Sun et al，2023)，且在降雹前20 min左右45 dBz的强回波面积只有10⁰ km²量级，前期如此小的强回波面积很难明确指示强对流云，而待其扩大到10¹~10² km²量级，距降雹已不足10 min。此外，＞55 dBz的回波在降雹前15 min首先出现在暖层，可能是由于霰粒子融化和雨滴碰并增长的共同作用，但并不能起到提示降雹的作用，而当冷层出现垂直方向上连续大于55 dBz的回波时(12:08)，距降雹时刻(12:14)只有6 min (图 3c)，可提前预警的时间很短。

另一方面，Z_DR大于1 dB、2 dB的回波在11:49前主要集中在暖层，对应暖层雨滴，而11:49先出现在0℃附近冷层的Z_DR回波与暖层不相连，不易直接判断其对应的是Z_DR柱还是其他水平取向的冰相粒子或杂波。11:51(降雹前23 min)，Z_DR大于1 dB、2 dB的回波从暖层延伸到冷层(图 3d, 3e)，符合Z_DR柱特征，且该特征持续出现，故可以指示Z_DR柱的存在。Z_DR柱出现时间与图 2中C_ZDRC诊断结果也基本一致(图 2)。Z_DR＞3 dB时，上述Z_DR柱特征出现的时间为11:57，较降雹时刻提前17 min(图 3f)。此外，Z_DR＞1 dB的回波面积在高空始终有零星的“噪点”，这可能是受水平冰晶、雪或其他杂波影响产生的，会干扰对Z_DR柱是否出现的判断。总体来说，个例1中采用Z_DR大于2 dB、3 dB阈值，可在降雹前20 min左右从回波面积的时间-高度演变中识别出Z_DR柱，从而对冰雹云进行早期识别。

2.2 2024年5月18日费县降雹个例

个例2中，在济宁S波段雷达上可明确看到降雹点附近有对流云形成时(13:57)，单体的Z_e已超过35 dBz(图 4a)。此时C_ZDRC虽然整片超过1 dB，但未出现强中心(图 4d)，Z_DR柱特征不明显。在降雹前35 min(14:03)，Z_e核心超过45 dBz，但回波范围仍然很小(图 4b)，此时C_ZDRC图上已有明显的3 dB左右的中心(图 4e)，表明了Z_DR柱的存在。在降雹时(14:37)，Z_e最大值在60 dBz左右，而C_ZDRC图中显示的Z_DR柱特征与2.1节个例类似，Z_DR柱围绕降雹点且位于降雹点的对流层中层上风方向(图 4c, 4f)。

回波面积的时间-高度演变(图 5)显示，Z_DR大于1、2、3 dB且从暖层延伸到冷层的回波分别比Z_H大于35、45、55 dBz的回波提前5~11 min出现。但结合图 4来看，至少在15:57 C_ZDRC图上大于1、2 dB的回波强中心不明显，不容易确认Z_DR柱的存在，此时根据图 5中1 dB、2 dB阈值的Z_DR诊断Z_DR柱不一定可靠。而Z_DR>3 dB时，出现Z_DR柱时间(14:03)较降雹时刻提前了35 min(图 5f)，图 4中该时刻Z_DR柱也较为明显，更适合作为提前预警冰雹云的指标。这一特征虽然只比Z_H>55 dBz出现提前5 min(图 5c)，且和Z_H>45 dBz出现时间相同(图 5b)，但考虑到初期强Z_H的面积很小，在降雹前约10~20 min才开始具有可观的面积，因此，使用Z_DR>3 dB作为预警指标较Z_H>45 dBz可使预警提前约10 min。

2.3 2024年5月18日济南降雹个例 2.3.1 济宁S波段雷达

在降雹前约1 h，降雹点上空无回波，但其东南侧邻近范围已出现具有Z_DR柱的对流单体(图 6a, 6b)。在降雹前39 min(18:01)，紧邻降雹点的东南侧出现回波，其Z_e虽然刚超过20 dBz，但C_ZDRC图已显示出超过3 dB的中心(图 6b, 6e)，表明此处有强上升气流和过冷雨滴区形成。在降雹时(18:41)，降雹点上空为Z_e超过60 dBz、沿西北—东南方向分布的强对流单体(图 6c)。此时C_ZDRC图上的Z_DR柱特征与前两个个例类似，位于降雹点的偏上风方向(图 6f)。

图 7显示，Z_DR大于1、2、3 dB的回波出现时间相同，且分别比Z_H大于35、45、55 dBz的回波早出现5、17、23 min。与以Z_H>45 dBz或更高阈值的回波面积为特征指标的预警方案相比，基于图 7中Z_DR柱可使预警提前约20 min。

2.3.2 泰山C波段雷达

泰山C波段雷达结果(图 8)和济宁S波段雷达结果(图 6)总体一致，在降雹前40 min(18:00)紧邻降雹点的东南侧出现Z_DR柱。不同之处表现在两方面：一是可能由于C波段雷达本身灵敏度高于S波段，且泰山C波段雷达距离降雹点更近，灵敏度更高，在Z_e图上出现更多0 dBz左右的非降水回波，相应位置C_ZDRC图上也有极个别大值噪点，但并不影响对Z_DR柱的判断，也很容易设置Z_H强度阈值进行数据剔除；二是在降雹时(18:42)，在降雹点西北侧(对流层中层上风方向)有更完整的环绕降雹点的Z_DR柱。

泰山C波段雷达结果(图 9)与济宁S波段雷达(图 7)类似。Z_DR大于1、2、3 dB的回波出现时刻相同，均较降雹时刻提前40 min。对于Z_H，由于C波段灵敏度高于S波段，Z_H大于45 dBz、55 dBz的回波出现时间比S波段早6~12 min。与Z_H>45 dBz或更高阈值的回波面积为特征指标的预警方案比较，基于图 9中的Z_DR柱可使预警提前10~20 min。

2.4 统计区域和屏蔽阈值对冰雹云早期识别的影响

就统计区域而言，前文个例中的回波面积时间-高度演变均是针对降雹点附近范围含有指定冰雹云的空间范围统计而来的。如果不做区域限定，而以单部雷达所能探测到的范围进行统计，以2.2节和2.3节中济宁雷达数据为例，结果如图 10a，10b所示，可见Z_DR大于2 dB、3 dB的回波从暖层延至冷层且连续特征出现时间(14：00左右，图中标记略)均早于当天首次降雹时刻(14:38)，但从图中很难诊断出当天第2、第3次降雹。因此，基于这种单部雷达所有探测范围的统计只对当天首次降雹有提前提示作用，但也不失为一种粗略的提前预警方案。

就屏蔽阈值而言，如果在计算Z_DR回波面积时不屏蔽Z_H < 20 dBz的回波，结果如图 10c，10d所示，可见噪点更多，更难分辨出Z_DR柱。这是因为介电常数较雨滴小、Z_H回波相对较弱的冰晶、雪等粒子也可能导致较大的Z_DR值。此外，图 10中暖层Z_DR面积大值中心经过雷达原始数据分区域溯源分析后被确认为边界层杂波，例如边界层中较弱的非降水性层状云、干扰和部分移动地物等目标会导致近地面出现起伏较大、噪声样分布的Z_DR，进而导致不同Z_DR阈值统计的回波面积在近地面出现大值中心，故在图 10中不应将大范围统计的Z_DR面积大值中心视为降水特征看待。但是，由于上述杂波并不会导致增加虚假的、由暖层延至冷层且连续出现的Z_DR回波，因而其对本文的分析方法没有实质影响。

2.5 Z_DR初生位置、冰雹落区和风廓线的关系

在本文讨论的三个降雹个例中，2024年5月18日的两个个例的首个降雹点均位于Z_DR柱初生位置附近，而2022年6月13日的首个降雹点则与Z_DR柱初生位置相距约10 km。为分析上述差异成因，图 11给出了降雹个例当天降雹点位置的风廓线，可见滨州上空的风向由地面的南风顺转为中层(5 km左右)的西偏北风，且中层最大风速超过10 m·s^-1，有利于形成沿中低层偏西气流移动的降水云系。而费县和济南上空的风向均从南风顺转为北风，中层最大风速不足6 m·s^-1。这种低层和中层风向相反的配置更易促成局地的对流环流，加之风速偏小，可能是导致当地对流云系移动特征不明显的关键因素。这种特征在进行冰雹预警时值得注意，但也需要对更多个例进行机理研究以确认。

2.6 Z_DR柱与Z_H核心的配置

本文三个个例的Z_DR柱特征与以往典型模型不同。典型超级单体0℃层高度的Z_DR柱形态可接近半环状并水平分布于中低层入流一侧(Kumjian et al, 2010; Ryzhkov and Zrnic, 2019)，故也被称为“Z_DR环”，且方向通常与云体延伸方向相近。而本文三个个例降雹时的C_ZDRC图所反映的Z_DR柱位置既与云体延伸方向相反，也不在低层入流一侧，即在云体的后部而非前部。这与Ren et al(2024)讨论过的孤立对流单体中情形更加类似。进一步做垂直剖面分析时，考虑到X波段和C波段中的衰减及其订正可能引入不确定性，下文仅使用S波段雷达数据，选取过Z_DR柱和降雹点的剖面进行分析。

在与图 6个例相同但比降雹更早一个观测时刻(18:35)，过降雹云前部较弱的Z_DR柱和降雹点(图 12a, 12b中AA′)的垂直剖面(图 12c)显示出强Z_H核心和Z_DR柱的配置。综合考虑对流层中层风向和回波的移动、发展方向，将南偏东方向称为云体的前部、北偏西方向称为云体的后部。云体前部Z_H偏弱的区域存在Z_DR柱(图 12c中偏左侧的填色大值区)，通常对应于上升气流区；后部是Z_H超过55 dBz但Z_DR较小的区域，与该个例的降雹点位置对应。过降雹点和降雹云后部Z_DR柱(图 12a, 12b中BB’)的垂直剖面(图 12d)显示，云体后部的Z_DR柱高于前部。Ren et al(2024)指出，在对流单体中期可由降水出流强迫抬升云体后部的中低层气流，从而在后部形成区别于主上升气流区Z_DR柱的、高度较低且强度较弱的“次级”Z_DR柱，但与本个例后部Z_DR柱比前部更高、更强的现象(图 12c、12d中填色大值区)正好相反，原因之一可能是低层和中低风向相反的配置(图 11)。低层入流恰好在中层风的下风方向，会导致降水出现后在一定程度上抑制下风方向近处的上升运动。但是，后部的Z_DR柱则不受这种抑制，从而可导致更强的上升运动。至于后部更强的Z_DR柱中的过冷雨滴是否为冰雹提供了更主要的形成和增长条件，有待在后续工作中配合数值模拟进一步研究。

3 结论与讨论

本文利用可表征Z_DR柱特征的C_ZDRC图和Z_DR回波面积时间-高度演变，结合山东省近年来的三次局地性强、发展迅速的降雹个例，开展了雷达数据特征和冰雹预警研究，主要结论如下：

(1) 在雷达回波面积较小的阶段，C_ZDRC图上已可见降雹点附近云中出现3 dB左右的Z_DR柱特征。在三个个例的回波面积时间-高度演变中，Z_DR大于2 dB、3 dB的Z_DR柱出现时间比降雹时刻分别提前约20~40 min，且比Z_H>45 dBz的面积超过100 km²或55 dBz回波出现时间提前10~20 min。表明C_ZDRC图和Z_DR回波面积时间-高度演变具有将降雹预警时间进一步提前的潜力。

(2) 本文使用的S波段和C波段业务雷达虽然在时间和空间分辨率上不如X-PAR，但C_ZDRC图和Z_DR回波面积时间-高度演变仍可显示出用于强对流预警的Z_DR柱特征，两种业务雷达观测和分析同一降雹个例的结果一致性较好。表明这两种用于冰雹预警的分析图形不仅可应用具有快速扫描能力的相控阵天气雷达，也具有应用于现有传统机械扫描式天气雷达的能力，能够较为直观、便捷地显示Z_DR柱出现的位置和时间。

(3) C_ZDRC图上Z_DR柱在三个个例中均出现了紧邻且环绕降雹点、位于降雹点的西北侧即对流层中层上风方向的特点，对指示降雹地点也具有一定的参考。

需要指出的是，就建立完备的目标识别方法而言，本文研究尚未能涉及虚警率的探讨，主要原因在于很难准确确定有强对流特征但未能产生降雹的观测样本。与可用遥感定位的闪电和可用大量自动雨量站捕获的极端强降雨相比，降雹尚缺乏空间连续性较高的原位观测手段。有限的地面观测站难以全面覆盖潜在的降雹区域，可能存在漏报的情况，导致对虚警率的评估存在很大局限。近年来，山东省人工影响天气办公室开发的“追雹者”微信小程序是一种很好的尝试，其允许群众上传拍摄到的降雹情况，对冰雹监测和相关研究具有很大的补充潜力。可在后续工作中针对性地收集和整编一批既具有强对流云特征又确定没有降雹和防雹作业的观测个例，结合降雹个例一并开展预警准确率的统计研究。

对冰雹云早期识别和预警的主要目标在于更科学、可靠地确定人工防雹作业位置和作业时机，以便在冰雹云形成的早期进行及时干预，而本文所展示的C_ZDRC图上Z_DR柱恰是一种比较明确的冰雹云早期观测特征，下一步可探索将其进一步应用到防雹作业指挥中。对于Z_DR回波面积时间-高度演变的业务化应用，可以选择按特定区域或炮点覆盖范围进行针对性统计，或者设计软件使预报员或观测员通过交互的方式针对性地查看一个手动划定区域中的统计结果。此外，就相控阵雷达数据的深入应用而言，对于一些雷达厂家没有存储订正前的数据和订正导致原始观测可能被污染的问题，也需要进行针对性研究。