QUALITY EVALUATION OF X-BAND POLARIMETRIC PHASED ARRAY RADAR QPE PRODUCT IN ZHUHAI
-
摘要: 以地面雨量站观测资料为基准,选取2021年6—8月的6个强降水过程,利用珠海4部X波段双偏振相控阵雷达QPE组网产品和对应区域雨量站逐小时降水资料,使用相关系数(CC,Correlation Coefficient)、相对偏差(RB,Relative Bias)、相对平均绝对误差(RMAE,Mean Relative Absolute Error)、均方根误差(RMSE,Root Mean Squared Error)、探测率(POD,Probability of Detection)、误报率(FAR,False Alarm Ratio)、临界成功指数(CSI,Critical Success Index),对相控阵雷达QPE产品质量进行评估。结果表明:QPE产品总体CC、RB、RMAE、RMSE分别为0.83、1.62%、47.59%、6.68 mm。有无降水探测率为88.66%,FAR为12.03%,CSI为83.99%;能够准确捕捉强降水落区、把握小时雨强变化趋势,雨强峰值谷值出现的时间匹配的较好。小时降水分级检验结果显示,QPE产品对于大雨强的评估效果相对较好,小时雨量≥16.1 mm时CC、POD、CSI均最大;雨强较小时,QPE产品明显高估降水强度,当小时雨强超过16.0 mm时,低估降水强度;RMAE随着雨强的增大逐渐减小,当雨强≥8.1 mm时,RMAE趋于稳定,为34% 左右。Abstract: This study evaluates the QPE product of X-band Phased Array Radar in Zhuhai with hourly precipitation data covering 73 rain gauges. Six heavy rainfall cases were selected and seven statistical analysis indexes were used in this study. These indexes include Root Mean Squared Error(RMSE), Correlation Coefficient (CC), Relative Bias (RB), Mean Relative Absolute Error(RMAE), Probability of Detection(POD), False Alarm Ratio (FAR) and Critical Success Index (CSI). Results show that QPE product can well capture the spatial pattern of heavy rain and the trend of rainfall intensity in terms of time with CC about 0.83, RB about 1.62%, RMAE about 47.59%, and RMSE about 6.68mm; the accuracy rate of clear or rain is 88.66%, while the FAR and CSI is 12.03% and 83.99%. Additionally, the QPE product performs best when hourly rainfall rate exceeds16.0 mm / h, with bigger CC, POD and CSI. The deviation of QPE product is positive for the lower rainfall rate, and it turn negative as the hourly rainfall rate over 16.0 mm / h. RMAE decrease as rainfall rate increase, but it stabilizes at around 34% as hourly rainfall rate over 8.0 mm/h, which means that the quality of QPE product is better in heavy rain.
-
Key words:
- Phased Array Radar /
- QPE product /
- evaluation /
- heavy rainfall
-
表 1 珠海X波段相控阵天气雷达QPE产品评估个例情况
过程编号 时间(北京时) 数据对 最大小时雨量/mm 1 6月1日02—24时 1 270 148.6 2 6月24日05—14时 719 56.8 3 6月28日02—13时 849 89.0 4 7月18日22时—19日15时 1 314 39.5 5 8月10日02—8时 509 97.2 6 8月11日01—9时 646 51.8 表 2 雨量站降水与QPE估测产品列联表
QPE估测产品 雨量站在阈值内 雨量站不在阈值内 QPE在阈值内 A B QPE不在阈值内 C D 表 3 QPE产品在不同雨强下的误差估计
小时雨量/mm CC RMSE/mm RB RMAE 0.1~1.0 0.28 1.15 257.74% 280.50 1.2~2.4 0.26 1.70 102.70% 118.49 2.5~8.0 0.34 1.88 42.72% 60.86 8.1~16.0 0.31 0.87 7.60% 34.93 ≥16.1 0.54 9.02 -26.53% 34.91 -
[1] 汪瑛, 冯业荣, 蔡锦辉, 等. 雷达定量降水动态分级Z-Ⅰ关系估算方法[J]. 热带气象学报, 2011, 27(4): 601-608. doi: 10.3969/j.issn.1004-4965.2011.04.018 [2] 高晓荣, 梁建茵, 李春晖. 雷达定量降水估计技术及效果评估[J]. 热带气象学报, 2012, 28(1): 77-88. doi: 10.3969/j.issn.1004-4965.2012.01.009 [3] 卓健, 陈少斌, 周冬静, 等. 快速动态分级法进行雷达定量估测降水方法研究[J]. 热带气象学报, 2018, 34(6): 856-864. doi: 10.16032/j.issn.1004-4965.2018.06.013 [4] 王红艳, 刘黎平, 王改利, 等. 多普勒天气雷达三维数字组网系统开发及应用[J]. 应用气象学报, 2009, 20(2): 214-224. doi: 10.3969/j.issn.1001-7313.2009.02.011 [5] 王红艳, 王改利, 刘黎平, 等. 利用雷达资料对自动雨量计实时质量控制的方法研究[J]. 大气科学, 2015, 39(1): 59-67. https://www.cnki.com.cn/Article/CJFDTOTAL-DQXK201501006.htm [6] 李建通, 郭林, 杨洪平. 雷达-雨量计联合估测降水初值场形成方法探讨[J]. 大气科学, 2005, 29(6): 168-178. https://www.cnki.com.cn/Article/CJFDTOTAL-DQXK200506015.htm [7] 邵月红, 张万昌, 刘永和, 等. BP神经网络在多普勒雷达降水量的估测中的应用[J]. 高原气象, 2009, 28(4): 846-853. https://www.cnki.com.cn/Article/CJFDTOTAL-GYQX200904016.htm [8] 张哲, 戚友存, 朱自伟, 等. 深圳S波段与X波段双偏振雷达在定量降水估计中的应用[J]. 气象学报, 2021, 79(5): 786-803. https://www.cnki.com.cn/Article/CJFDTOTAL-QXXB202105006.htm [9] 吕晓娜, 牛淑贞, 袁春风, 等. SWAN中定量降水估测和预报产品的检验与误差分析[J]. 暴雨灾害, 2013, 32(2): 142-150. https://www.cnki.com.cn/Article/CJFDTOTAL-HBQX201302006.htm [10] 何彩芬, 朱龙彪, 董杏燕, 等. 宁波夏季强对流和台风短时暴雨雷达回波特征对比分析[J]. 气象, 2006, 32(11): 67-72. https://www.cnki.com.cn/Article/CJFDTOTAL-QXXX200611010.htm [11] 宋林烨, 陈明轩, 程丛兰, 等. 京津冀夏季雷达定量降水估测的误差统计及定量气候校准[J]. 气象学报, 2019, 77(3): 497-515. https://www.cnki.com.cn/Article/CJFDTOTAL-QXXB201903009.htm [12] 李晶, 包云轩, 蒋义芳, 等. 多普勒天气雷达1h降水产品的质量评估[J]. 气象科学, 2015, 35(3): 340-345. https://www.cnki.com.cn/Article/CJFDTOTAL-QXKX201503014.htm [13] 程元慧, 傅佩玲, 胡东明, 等. 广州相控阵天气雷达组网方案设计及其观测试验[J]. 气象, 2020, 46(6): 823-836. https://www.cnki.com.cn/Article/CJFDTOTAL-QXXX202006009.htm [14] 吴罛, 刘黎平, 仰美霖, 等. X波段双偏振雷达相态识别与拼图的关键技术[J]. 应用气象学报, 2021, 32(2): 200-216. https://www.cnki.com.cn/Article/CJFDTOTAL-YYQX202102006.htm [15] 林良勋, 冯业荣, 黄忠, 等. 广东省天气预报技术手册[M]. 北京: 气象出版社, 2009: 104-105. [16] CHANDRASEKAR V, BRINGI V N, BALAKRISHNAN N, et al. Error structure of multiparameter radar and surface measurements of rainfall. Part Ⅲ: Specific differential phase[J]. J Atmos Oceanic Technol, 1990, 7(5): 621-629. [17] MATROSOV S Y, KROPFLI R A, REINKING R F, et al. Prospects for measuring rainfall using propagation differential phase in X- and Ka-radar bands[J]. J Appl Meteor, 1999, 38(6): 766-776. [18] MAKI M, IWANAMI K, MISUMI R, et al. Semi-operational rainfall observations with X-band multi-parameter radar[J]. Atmos Sci Lett, 2005, 6(1): 12-18. [19] 张羽, 田聪聪, 傅佩玲, 等. 广州X波段双偏振相控阵天气雷达观测试验进展[J]. 气象科技进展, 2020, 10(6): 80-85. https://www.cnki.com.cn/Article/CJFDTOTAL-QXKZ202006017.htm [20] CHEN S, HONG Y, CAO Q, et al. Similarity and difference of the two successive V6 and V7 TRMM multisatellite precipitation analysis performance over China[J]. J Geophys Res: Atmos, 2013, 118(23): 13 060-13 074. [21] CHEN S, HONG Y, CAO Q, et al. Performance evaluation of radar and satellite rainfalls for Typhoon Morakot over Taiwan: Are remote- sensing products ready for gauge denial scenario of extreme events?[J]. J Hydrol, 2013, 506(25): 4-13. [22] 刘俊, 黄兴友, 何雨芩, 等. X波段相控阵气象雷达回波数据的对比分析[J]. 高原气象, 2015, 34(4) : 1 167-1 176. https://www.cnki.com.cn/Article/CJFDTOTAL-GYQX201504028.htm [23] 胡志群, 刘黎平, 楚荣忠, 等. X波段双线偏振雷达不同衰减订正方法对比及其对降水估测影响研究[J]. 气象学报, 2008, 66(2): 251-261. https://www.cnki.com.cn/Article/CJFDTOTAL-QXXB200802010.htm [24] 张蔚然, 吴翀, 刘黎平, 等. 双偏振相控阵雷达与业务雷达的定量对比及观测精度研究[J]. 高原气象, 2021, 40(2): 424-435. https://www.cnki.com.cn/Article/CJFDTOTAL-GYQX202102019.htm -