利用一种自动识别算法移除天气雷达反射率因子中的亮带
An Automatic Identification Algorithm for the Removal of Bright Band from Reflectivity of CINRAD/SA
-
摘要: 该文介绍了一种自动识别和移除雷达反射率因子资料中亮带的算法, 并对该算法进行了初步测试。该算法利用的是插值到直角坐标系中的雷达反射率因子资料, 其配置和运行也相对简单, 但却对移除亮带比较有效。首先, 设定一套雷达反射率因子垂直廓线的理想模板, 这些理想的模板能够在最大程度上反映不同亮带存在区域的雷达实际反射率因子的垂直廓线特征。然后, 在水平方向每个点上, 进行理想模板和实际反射率因子垂直廓线在垂直和水平两个方向上的拟合和差异计算, 来自动识别雷达反射率因子中存在的连续亮带区域。最后, 利用亮带之上和亮带之下的反射率因子值对亮带中的反射率因子值进行插值纠正, 就可以移除亮带。利用位于天津塘沽的我国新一代天气雷达 (CINRAD/SA) 的反射率因子资料, 通过个例分析和准业务运行试验, 均表明这个简单算法可以识别和移除绝大多数影响雷达定量降水估计的反射率因子亮带区域, 但是实际雷暴区域的反射率因子特征受到该算法的影响比较小。计算分析还表明, 在京津地区的初夏, 上述亮带区域一般容易出现在2.5 km左右的高度处。Abstract: An automatic identification algorithm for the removal of bright band from radar reflectivity data is introduced and tested by case analysis and pre-operational run. The algorithm has good results of removal of bright band based on the radar reflectivity data which have been transformed into a regular 3-dimensional Cartesian grid while configuration and running setup are relative simple. Firstly, the algorithm performs the recognition of vertical bright band profile. The algorithm considers each point independently, and attempts to find the evidence of a bright band by using a pattern recognition approach. Based on the typical reflectivity of the vertical profile attributes in the bright band deducted from statistics, a set of ideal templates of radar reflectivity vertical profiles is established, which can resemble the actual radar reflectivity characteristics in different bright band areas at the utmost. For each horizontal grid the algorithm attempts to couple the templates to the measured vertical profile of radar reflectivity. The couple function is obtained by computing the slope and the correlation of the least squares regression between the template and the measured profile. Secondly, the algorithm finds continuous regions of bright band. Regions of bright band may be further recognized by the horizontal uniformity of the signature. The algorithm attempts to use the goodness-of-fit and height-of-best-fit field computed above to gain further clues of a bright band. A grid kernel-based approach is used, in which a kernel is composed of 5×5 grids in the field. For the grids in the kernel, the algorithm counts the number of difference between contiguous goodness-of-fit values exceeding a given threshold or the difference in the height-of-best-fit values for contiguous points exceeding a given threshold. A high interested value is assigned to points with low numbers, since this indicates uniformity. So bright band areas are identified as those fields with high interested values. Finally, the algorithm removes the bright band signature. It removes the continuous regions of bright band by correcting the reflectivity within the bright band, which employs linear or non-linear interpolation methodology based on the measured reflectivity values in the profile above the band and the measured values below the band. Some thunderstorm cases analysis and pre-operational running tests suggest that the algorithm is able to successfully identify and remove most bright band regions, which influences quantitative precipitation estimation, based on CINRAD/SA radar reflectivity located at Tanggu of Tianjin. However, the algorithm has less effect on actual storms reflectivity data. Results of the analysis also demonstrate that the bright band regions usually appear at 2.5 km over Beijing and Tianjin in early summer. Because the algorithm is relatively easy to configure and run, it is appropriate to use in real-time or operational mode. Of course the algorithm has limitations and needs to be continuously improved for actual operational applications.
-
Key words:
- radar reflectivity;
- brightband;
- identification;
- removal
-
图 4 2003年6月27日14:48天津雷达2.5 km高度的反射率因子 (a) 和移除亮带后的反射率因子 (b)
(从1到2是图 5垂直剖面的位置)
-
[1] Fabry F, Zawadzki I. Long-term radar observations of the melting layer of precipitation and their interpretation. J Atmos Sci, 1995, 52: 838-851. doi: 10.1175/1520-0469(1995)052<0838:LTROOT>2.0.CO;2 [2] 俞小鼎, 姚秀萍, 熊廷南, 等.新一代天气雷达原理和应用讲义 (修订本).北京:中国气象局培训中心, 2004. [3] Smith C J. The reduction of errors caused by bright bands in quantitative rainfall measurements made using radar. J Atmos Oceanic Technol, 1986, 3: 129-141. doi: 10.1175/1520-0426%281986%29003%3C0129%3ATROECB%3E2.0.CO%3B2 [4] Kitchen M, Brown R, Davies A G. Real-time correction of weather radar data for the effects of bright band, range and orographic growth in wide-spread precipitation. Q J R Meteor Soc, 1994, 120: 1231-1254. doi: 10.1002/(ISSN)1477-870X [5] Hardaker P J, Holt A R, Collier C G. Amelting layer model and its use in correcting for the bright band in single polarisation radar echoes. Quart J Roy Meteor Soc, 1995, 121: 495-525. doi: 10.1002/qj.49712152303/full [6] Kitchen M. Towards improved radar estimation of surface precipitation at long range. Q J R Meteor Soc, 1997, 123: 145-163. doi: 10.1002/(ISSN)1477-870X [7] Sanches-Diezma R, Zawadzki I, Sempere-Torres D. Identification of the bright band through the analysis of volumetric radar data. J Geophys Res, 2000, 105 (D2): 2225-2236. http://adsabs.harvard.edu/abs/2000JGR...105.2225S [8] Gourley J, Calvert C. Automated detection of the bright band using WSR-88D data. Wea Forecasting, 2003, 18: 585-599. doi: 10.1175/1520-0434(2003)018%3C0585%3AADOTBB%3E2.0.CO%3B2 [9] 李柏. 中国新一代天气雷达计划. 中国气象局北京城市气象研究所/中国气象学会城市气象学委员会/美国国家大气研究中心 (NCAR) 编. 中美强对流天气临近预报技术国际研讨会文集. 北京: 气象出版社, 2004: 68-82.