反射率质量控制.docVIP

WSR-88D 反射率质量控制 作者：Jian Zhang, Shunxin Wang, and Beth Clarke 翻译：张深寿 福建省龙岩市气象局 364000 1.介绍 WSR-88D（天气警戒雷达-1988多普勒）能探测气象和非气象目标的散射，非气象目标的散射能污染雷达观测数据和由它产生的下游数据。而且，受污染的雷达数据会使下游应用和天气雷达产品（比如降水估测）产生错误。这篇文章介绍一个反射率质量控制（QC）算法，以能够在反射率中自动分辨和去除非降水回波。 在下一节将展现反射率质量控制（QC）算法。第三节展现不同非降水回波通过反射率质量控制后的对比结果。第四节为总结概要。 2.方法 2.1不同回波的特征 雷达反射率回波包含两个总要的内容：降水回波和非降水回波。非降水回波在反射率主要有三种类型：第一种类型为由昆虫、鸟类、尘土、同频干扰和其他在大气中的微粒足够大时散射产生的一些杂波；第二种为山脉和高大建筑物的地物回波，或海浪回波；第三种为由于大气折射产生的超折射地物或海浪回波，它是由于雷达波束在特殊的大气折射条件时向地面弯曲传播，它们就是在雷达气象中所谓的“异常传播（AP）”，雷达波束正常的传播路径通常弯曲向上传播。 第一种非降水回波（除了同频干扰）的产生，通常低于离地面3公里高度。而且，这些生物或灰尘目标产生噪声强度很弱（低于35dBZ）和水平尺度很小。地物和海浪回波，包括超折射地物海浪回波通常和地形或海洋有关系，并且出现在雷达探测的最低仰角切面。 降水回波可以分为两大类：对流的和层状云的。对流降水系统常和强降水和强的水平反射率梯度有关，它们水平尺度通常比层状云的降水系统小。对流降水回波一般有远大于5公里高度的垂直尺度。对流性降水回波大的垂直尺度使得我们很容易把它和非降水回波分辨开来。 层状云降水回波的特征表现为轻度到中等的降水、弱的水平反射率梯度、大范围覆盖（通常远大于500 km2）。层状云降水的云顶在4.5公里以内或甚至3公里（离地面）的情况也是有的，这种浅高度使它很难和非降水回波分辨开来。然而，层状云降水回波的降水强度比较均匀，它的水平强度梯度很小。对于主要的层状云降水包含小尺度的噪声，提出一个标准的水平平滑度有助于分辨层状云降水回波和非降水回波。这个值得关注，虽然，有时候非正常传播（AP）产生的回波比如在严重的超折射条件下一个平原的地形，也产生水平梯度小、平滑的反射率回波分布。在这个情形，非正常传播（AP）产生的回波很难和层状云降水回波分辨开来，正确的分辨它们需要使用另外的标准，其中一个标准就是回波宽广度。前面有提到，层状云降水回波的面积通常在100到1000平方公里。非正常传播（AP）产生的回波，除了非常少的情况，通常局部化、在径向上很小的连续。因此，径向上相邻连续的检查，能分辨和去除一些非正常传播（AP）产生的回波。 2.2质量控制方法的步骤 图1 是一个质量控制的流程图。下面详细描述这些步骤。 图1 质量控制的流程图 2.2.1 数据输入和噪声过滤 质量控制需要的数据是3维球面坐标（方位、距离、仰角）的体扫资料。先对原始的反射率数据进行一个简单的过滤，去除孤立的点和线。对于给定的反射率库X，以它为中心点，取5×5的数据，得到非空反射率探测值个数N，那么非空值的百分比（Px）的计算公式为： Px＝N/Ntotal 这里Ntotal等于25代表5×5数据的个数。如果Px小于一个给定的阈值（缺省为75％），则认为这个点X是孤立的值，把它删除设为空。这个过滤对于去除斑点的噪声回波是非常有效的（如图2）。 图2 0.5度仰角反射率图噪声过滤前后，a过滤前，b过滤后 2.2.2 水平反射率结构 反射率水平平滑度是一个测量参数，叫做反射率水平结构（TDBZ，Kessinger er. al.，2003）。它表示顺着径向相邻的反射率值的方差。给定库的TDBZ计算公式是： 其中i，j分别是反射率库在方位和距离方向上的序号值，Z代表反射率值，单位dBZ，nrays和ngates是以给定库为中心在方位和距离方向上的库的个数。N＝nrays×ngates代表用于计算的数据矩形的数据个数。nrays和ngates的缺省值都为7。 图3显示了一个飑线的反射率图和TDBZ图。反射率的仰角是0.5°。大的TDBZ值和对流单体有关，也一样和在飑线后面的非正常传播（AP）产生的回波一样有关。所以，单单用结构参数还不能正确分辨降水和非降水回波。 图3 飑线在0.5°的基础反射率（a）和计算后的TDBZ（b）图 2.2.3 更高仰角的参考切面 象2.1节描述，对流性降水和许多层状云降水的云顶高度比3公里高，而非降水回波（除了同频干扰）通常在3公里以下。因而，离地面3公里和以上高度的垂直方向的反射率连续性就成为分辨降水和非降水回波的非常