中国气象学会主办。
文章信息
- 万玉发, 王钰, 金鸿祥 . 2013.
- WAN Yufa, WANG Yu, JIN Hongxiong . 2013.
- 雷达与雨量计同步结合区域型估算降水方程的误差分析
- Error analysis of the regional radar-raingauge synchronously integrated precipitation estimation equation
- 气象学报, (2): 332-343
- Acta Meteorologica Sinica, (2): 332-343.
- http://dx.doi.org/10.11676/qxxb2013.023
-
文章历史
- 收稿日期:2012-04-20
- 改回日期:2012-10-25
2. 武汉中心气象台, 武汉, 430074
2. Wuhan Central Meteorological Observatory, Wuhan 430074, China
目前主要存在两个可业务应用的降水测量系统:一是以雨量计为代表的地面测量系统,二是以天气雷达为代表的空中测量系统。前者用雨量筒测量地面降水量,直接而准确,可以获取降水的地面真值;但雨量筒测雨面积只约3.14 dm2,只能视作降水的点测量且在地面上呈现离散分布,因此,仅用雨量站网资料难以得到可靠的细网格分辨率地面降水量分布。后者是对空中雨滴大小及密度的一种电磁学间接测量,难以直接应用;但雷达能在半径约200 km内快速、实时得到分辨率为1 km×1 km反射率因子Z(mm6/m3)的平面分布,为其最为突出的优点。雷达定量降水估算(QPE,Quantitative Precipitation Estimation)的研究任务(Zhang et al,2011)是将这两类系统结合起来,即以雷达探测为基础、以地面雨量站网为标准化系统,将雷达所测反射率因子Z的某种空中分布,按照一定的“空-地”点对应方式转换成为地面雨强R(mm/h)分布,并要求在雨量站处通过转换和调整后的估算值,尽量和雨量计的实测值相等。由此可以研究和设计出多种雷达降水的估算和评估理念、技术和算法,从而形成多种雷达降水的估算和评估体系。
无论何种雷达定量降水估算技术,以雨量站网测量数据为标准去衡量评估雷达估算,必然会表现出差异(误差),并且,会有相当大的不确定性。许多学者(Joss et al,1990;Zawadzki,1984;Collier,1984)对此做过研究,认为影响因素多而复杂。仅就雷达和雨量计结合估算降水而论,主要可概括为5个方面的因素。其一,空-地对应点位置确定的合理性问题。联系雷达反射率因子Z和地面降水强度R确定幂律参数已是一个长时间的问题(Krajewski et al,1991),其核心是雷达和雨量计的数据采样问题。一种是时-空差异采样方法,即按照雨团下落并随风飘移的轨迹,找出每个雨量站在雷达探测层中雨团对应的初始点位置,这是相当困难的,但Mittermaier等(2004)和Lack等(2004)提出了雷达测雨因风漂移引起的误差订正方法,这样的空-地点对应具有直接的物理联系,可称为有理对应;另一种是瞬时垂直同步采样,显然这是一种不考虑雨滴下落时间和风对雨滴移动影响而存在时-空差异的采样方式,这种空-地对应点的采样不可能在同一雨团中进行,这和经典Z-R关系(Atlas et al,1990)所体现的同雨团样本是不相符的,本文将这种没有直接物理联系的空-地对应称为无理对应。其二,空-地对应点的面积不对称问题。雷达和雨量计空-地对应点采样面积之比约为亿分之三,降水量的时-空变化非常大,这给具有不同时-空分辨率的仪器观测做比较带来了重大影响(Steiner et al,2001);Zawadzki(1975)也明确指出了雷达测量的平滑误差,并断言试图从雷达和雨量计比较中建立Z-R关系仅在均匀的降水中才是成功的。其三,雨量站网布设引起的统计抽样随机性问题。雷达结合雨量计估算降水,实质上是由N个雨量站组成的随机抽取子样本所求到的统计规律,来推测到整个雷达探测面积这个统计总体的每个点值,其结果必定因这个子样本的选择差异(雨量站密度及站址分布)而发生变化,其实用雨量计校准面临着同样的问题(Br and es,1975; Ahnert et al,1986; Cheze,1999;田付友等,2010)。其四,降水天气条件及其变化产生的随机问题。实践中发现,即便在同一雨量站网条件下,采用同一个估算技术也会得到不同的估算结果,Joss等(1990)指出,雨滴谱的自然可变性是雷达测量降水不确定性的一个重要原因,本文认为这可能反映了降水天气条件的变化。其五,区域型雷达估算降水技术带来的误差问题。虽然Doviak等(2006)认为,选择合适的Z-R关系可以补偿某些系统误差,但在一个覆盖区域各个格点采用一个公用的转换关系,即估算参数不随地理位置不同而变化,这显然会产生误差。
由万玉发等(2008)提出的“基于准同雨团样本概念雷达和雨量计的实时同步结合方法”,即RASIM(RAdar-raingauge Synchronously Integrated Method in real-time)(Wan et al,2010),可与大多数估算技术一起归并于区域型雷达降水估算技术,这些估算技术概括有:最优化法、A值平均法、概率配对法以及Z-R关系结合校准的方法,如平均校准法、距离加权法、卡尔曼滤波法、变分法等(张培昌等,1992,2001;Rosenfeld et al,1994; 尹忠海等,2005)。但RASIM采用了时间积分垂直同步采样,在无风和弱风条件下也可得到大部分的同雨团样本,可称为雷达和雨量计的积分对应或准有理对应。RASIM在固定指数bf(Smith et al,1997)基础上,将“小时等效反射率因子ZB(mm6/m3)”和相应的雨量计采样QG(mm)直接对应并构成具有幂律形式的ZB-QG关系,做到了转换和调整的合二为一。本文进一步综合提出4种RASIM式的降水估算方程,经过比较,ABS区域云雨转换公共系数的估算方程最优;其次,文中提出了单站降水误差估算因子μi是造成估算误差的基本元素,并分析了其产生的原因;在此基础上又提出了几种误差表示和评估的公式,研究了各自的特点及相互关联;最后,建立一种可通过设定μi合理区间来判别雷达和雨量计错误数据对的质量控制新方法。2 RASIM的区域型降水累积量估算方程
本节将对RASIM研究成果(万玉发等,2008;Wan et al,2010)与本文相关的内容做出概括的提要和阐述,并在概念和方法上做进一步的扩充。2.1 单站云雨转换方程及转换系数
首先,将雷达和雨量计的空-地对应关系做一个概念代换。若将雷达在空中探测的降水粒子看作云的话,即可将空-地对应ZB-QG关系称为云雨的对应关系或云雨转换方程。虽然这种概念代换对于大气物理中云和雨滴的区分是不严格的,但本文的目的仅仅是为了便于表述而已。需要强调的是本文所研究的对象是指云雨1 h累积量的转换方程和转换系数。2.1.1 以AB为云雨转换系数的单站云雨转换方程
若雷达覆盖区域内有N个雨量站点,对于第i个站点(x,y)处Tj时次内存在有雷达和雨量计的垂直空-地积分对应关系

在式(1)中求出ABi,并且,标以时-空坐标,则

在Tj时次的1 h内,取第i个雨量站点处QGi直接与其垂直对应的反射率因子Z平均值ZMi构成幂律关系

在雷达有效探测区域S内与Tj时次形成一个{ABi}S或{AMi}S的数据集合体,一般各个ABi(AMi)不尽相同,与自然界中降水具有很大的时-空变率相关,为各种降水天气系统的尺度、种类和部位以及所处地理环境等综合影响所致。从{ABi}S和{AMi}S的集合体中,需要研究的是代表区域天气总体共性的区域公共系数;区域内{ABi}S和{AMi}S数据的离散和聚合状态。2.2.1 AB类型的区域云雨转换公共系数
{ABi}S集合的第1个统计代表值是ABi的算术平均值AB,为AB类型的第1个区域云雨转换公共系数




同样,{AMi}S集合的第1个统计代表值就是AMi的算术平均值AM,为AM类型的第1个区域云雨转换公共系数


在区域S内各个雨量站点处其空中的ZBi或ZMi和地面的QGi都是实测的,根据式(1)、(3)即可计算各个单站云雨转换系数ABi或AMi。而对于S区域内非雨量站点(x,y)处,即使已拥有经雷达实测的ZB(x,y)或ZM(x,y)的面分布也不能得到AB(x,y)或AM(x,y)的面分布。为了使得这个问题的近似解决,可利用固定指数bf和区域云雨转换公共系数以及由雷达实测的ZB(x,y)或ZM(x,y)的面分布,按照幂律形式来构建区域型降水累积量估算方程,在这类方程中需要待定的只有地面降水量的估算值QR(x,y)了。本文提出了4个区域型云雨转换公共系数,可相应建立如下4个区域型降水估算方程,分别称为ABS型、AB型、AMS型和AM型




这4个方程之所以称为区域型降水估算方程,是因为采取了区域型云雨转换公共系数,根据这些方程能估算出区域S内1 km×1 km每个点的1 h降水累积量。必须着重指出,由于ABS型方程在区域S内所有站点的降水估算总量与地面实测总量具有相等的特点(详见第4节),所以在5 a多来的业务实践中被优先采用,本文也对此类降水估算方程做重点分析和研究。3 雷达降水估算误差的基本因子
在第2节中所介绍的区域型降水估算方程是用一个公共系数(如ABS)来代替区域S中各点真实的AB(x,y)进行地面降水估算的,因此,难免或多或少存在估算误差。雷达降水估算误差应该发生在整个区域S内各个格点上,但实际上只能在雨量站点上进行检验或评估。3.1 单站误差估算因子μi
现设在Tj时次内QRi为第i个雨量站点的雷达估算雨量,而QGi为相应实测雨量,ΔQRi为该点降水估算的误差量或称单站降水误差估算量,并根据式(1)、(8),经过适当代换则有


第一,式(13)的第1个等式是一个通用的表达式,表示μi是单站误差估算量与单站降水量的比值,其物理意义为单站单位实测雨量所含有的误差估算量,故μi也可称为单站误差估算率。
第二,式(13)的第2个等式(包括第3个等式)则是RASIM方法所特有的表达,表示μi是该站云雨转换系数在为指数的幂值上相对于公共转换系数ABS的相对离差,反映了μi产生的原因在于ABi对于ABS的偏离度。显然,当ABi偏离ABS越大,则μi的绝对值越大;反之越小;当ABi=ABS时μi=0(图 1a)。图中等AB线即为等μi线,等ABS线将图像区域分成两部分,左上部分区域为μi<0,右下部分区域为μi>0。
![]() |
图 1 ABS型降水估算方程误差分析(a.单站降水估算误差ΔQRi=μi×QGi示意图,曲线为等AB(或等μi)的ZB-QG关系线,b.单站误差估算因子μi的函数曲线,每根曲线描绘了在不同bf指数条件下μi随ABi/ABS变化) Fig. 1 Analysis chart of the error from the ABS regional radar precipitation estimation equation(a. Sketch map of ΔQRi=μi×QGi,the misestimate of single station,in which each curve is the ZB-QG relationship line with equal AB(or equal μi); and b. function curves of μi,the misestimate factor of single station,in which each curve shows the change of μi with ABi/ABS under the different bf factor) |
第三,式(13)的第3个等式表示μi是ABi/ABS和bf的函数(图 1b)。在此图中横坐标为ABi/ABS,纵坐标为μi,原点设在(1,0)上,在图中第一和第三象限的每根曲线描绘了在不同bf指数条件下μi随ABi/ABS变化。图中可见ABi<ABS时μi<0,表示降水被ABS方程低估;而ABi>ABS时μi>0,则为高估。
第四,式(12)表示了单站降水误差估算量ΔQRi是单站误差估算因子μi和单站实测降水QGi的乘积。在图 1a中存在两种特殊情况,其一是当μi为常数时,即沿某根等AB线上各点ΔQRi随QG的变化成正比;其二,当QG为常数即均匀降水时,沿等QG线上各点的ΔQRi随μi变化成正比,于是离ABS线更远的ABi点就具有较大的降水估算误差量。3.2 区域平均绝对误差估算因子μA
现定义|μi|的算术平均为区域平均误差估算因子μA,并将式(13)代入则为

在第3节中,主要研究了单个雨量站位置上降水估算误差的表达式及其意义,同时指出了区域平均绝对误差估算因子μA及其特性。本节将研究区域S内N个雨量站上降水误差估算的几种整体表达式及其意义,它比单站误差估算的表达式更为重要。4.1 区域降水误差估算率μS
首先定义区域降水误差估算率μS为带有正负的各单站降水误差估算量在区域内的总和与区域对应各站点降水实测总量的比值,用于衡量区域降水总量的估算误差。区域S在Tj时次内,令QRS和QGS分别表示区域各站点降水估算总量和实测总量,因此,区域降水估算的总误差简称区域降水误差估算量ΔQRS为

在第2节中已经介绍了RASIM体系下的4个区域型降水估算式(8)、(9)、(10)、(11),因而式(15)中QRi在同一时次和同一区域内会有4种不同的估算值,相应产生4种不同的区域降水误差估算量,分别用ΔQRS(ABS)、ΔQRS(AMS)、ΔQRS(AAB)、ΔQRS(AAM)表示。4.1.1 ABS型区域站点降水总量估算的特点
用ABS型降水估算方程时产生的区域降水误差估算通过式(15)、(12)、(13)可表示为:



应该充分注意到,4个不同类型降水估算方程(式(8)、(9)、(10)、(11))都有相同的幂律数学形式,区域云雨转换系数分别是ABS、AB、AMS、AM。而小时等效反射率因子前两者用ZB,后两者则用ZM,其算法有所不同。因此,可以参照4.1.1 节的原理,分别求出另外3种区域降水估算误差表达式



AB型优于AM的根本原因在于AB是基于反射率因子Z积分而得,而AM则是基于Z的简单平均。
时间 | ABS | AMS | AB | AM | ||||
μS | μ|S| | μS | μ|S| | μS | μ|S| | μS | μ|S| | |
01时 | 0.00 | 0.30 | -0.13 | 0.34 | -0.58 | 0.60 | -0.62 | 0.64 |
02时 | 0.00 | 0.33 | -0.15 | 0.37 | -0.64 | 0.65 | -0.65 | 0.67 |
03时 | 0.00 | 0.31 | -0.20 | 0.37 | -0.51 | 0.56 | -0.55 | 0.59 |
04时 | 0.00 | 0.36 | -0.20 | 0.45 | -0.67 | 0.71 | -0.67 | 0.72 |
05时 | 0.00 | 0.29 | -0.20 | 0.38 | -0.43 | 0.49 | -0.47 | 0.54 |
06时 | 0.00 | 0.25 | -0.24 | 0.37 | -0.44 | 0.48 | -0.52 | 0.56 |
07时 | 0.00 | 0.26 | -0.22 | 0.37 | -0.92 | 0.93 | -0.92 | 0.92 |
08时 | 0.00 | 0.29 | -0.20 | 0.37 | -0.50 | 0.53 | -0.60 | 0.63 |
09时 | 0.00 | 0.32 | -0.21 | 0.41 | -0.62 | 0.65 | -0.66 | 0.68 |
10时 | 0.00 | 0.33 | -0.20 | 0.39 | -0.69 | 0.71 | -0.71 | 0.73 |
11时 | 0.00 | 0.33 | -0.24 | 0.38 | -0.62 | 0.64 | -0.66 | 0.68 |
12时 | 0.00 | 0.33 | -0.27 | 0.44 | -0.73 | 0.75 | -0.76 | 0.79 |
![]() |
图 2 图 2 武汉雷达200 km半径范围的雨量站网分布 (·为常规站,·为加密站) Fig. 2 Distribution of the rain gauge network within a radius of 200 km from the Wuhan radar (·: conventional observation station and ·: intensive observation station) |
4.1节中已经证明唯有在用ABS型降水估算方程的条件下,在S区域内Tj时次必有区域误差估算率μS=0,即ΔQRi=0,但是具体到每个站点误差估算量一般不为0,应有
|ΔQRi|≠0,由此给出区域降水绝对误差估算率μ|S|的定义:各单站降水误差估算量绝对值在区域S内的总和与区域对应各站点降水实测总量的比值,用于衡量区域降水总的绝对误差估算率,是应用最为广泛的降水误差估算的评估判据。结合式(12)、(13)可有

现引入区域平均绝对误差估算量EN,根据式(22)并令区域S中各站实测平均降水量为,QGi,则


这两个量常用于降水估算评估中,但少有分析其关系的。在式(22)的第2式的分子和分母中引入μA和两个量,整理后有:



第一,在区域N个站点中,当|μi|和QGi两者均为常数或者分别单独为常数,此时η=1,显然μ|S|=μA。其中,当QGi=QG时,表明S区域内发生均匀降雨。有意义的是|μi|与QG存在一个统计关系(图 3a),即QG较小时一般|μi|较大,而QG较大时|μi|值较小。在此给出物理解释:万玉发等(2008)、Wan等(2010)论证了雷达和雨量计采
样时-空差异取决于雨滴的直径、雨量站离雷达的距离和环境场的风速,其中雨滴的直径是一个最主要因素。当雨滴较小时雨滴下落速度较慢,易受到雨量站距离和风速的影响,在这种情况下,时间积分垂直同步采样难以取得准同雨团样本,于是各站点的|μi|就较大,所以,在区域S内出现较小的均匀降雨时必然出现较大的μA或μ|S|。反之,当区域S内出现较大的均匀降雨时,雨滴下落速度较快,受到雨量站距离和风速的影响较小,时间积分垂直同步采样近似于准同雨团样本采样,各站点的|μi|相应较小,μA或μ|S|必然变小。在实际大气降水条件中,理想的均匀降雨极为罕见,但μ|S|与μA相差较小的情况还是时有发生,此时可以认为在S区域Tj时次产生了一场比较均匀的降水。
![]() |
图 3 2008年7月2日01—15时降水过程μ|S|和μA的关系(a. QG-|μi|点聚分布,黑色粗实线为散点趋势线; b. 同一个过程的μ|S|-μA点聚分布,黑色粗实线为μ|S|=μA线) Fig. 3 Analysis chart of the relation between μ|S| and μA for the rain event during 01:00-15:00 BT 2 July 2008(a. QG-|μi| scatter diagram,black thick real line is scatter trendline; and b. μ|S|-μA scatter diagram,black thick real line is line of μ|S|=μA) |
第二,在区域内N个站点中,当|μi|和QGi两者均不为常数,此时η≠1,一般μ|S|≠μA。此时在区域S中出现了大、中、小雨混合纷呈不均匀的现象。统计表明,小雨的站点频次一般高于大雨的频次(图 3a)。由式(14)可见,用μA作评估判据会过于偏向于小雨站点的贡献(Klazura et al,1995);而μ|S|是各|μi|按QGi分配的加权平均值,因此,一般μ|S|<μA(图 3b)。由此可见,有时人们将μA作为降水评估判据是不太合理的,从本文的分析来看μA,只是各站点降水误差估算因子绝对值的平均,可诠释为多个单站的ABi相对于ABS的聚合与离散程度的整体表示,但μA与雨量大小无关,无力判定降水估算的误差量。4.5 降水过程的绝对误差估算率
为科学评估降水过程的误差估算,区域S内N个站点在一次降水过程T(有M个Tj时次)中雷达降水误差估算的统计特征量可有下列3种表示方式或3种检验评估判据。
(1)区域平均单站降水过程绝对误差估算率(μT|S|)

(2)降水过程平均的区域绝对误差估算率(μT|S|)

(3)过程和区域降水总体绝对误差估算率μT|S|

计算2008年7月2日01—15时降水过程估算评估得出μT|S|、μT|S|、μT|S|分别为0.25、0.27、0.41。由此可见,对于同样一个降水过程其评估的3个量是有差别的,相比而言,μT|S|的算法比较严格,在开展降水过程评估时应优先选用。而μT|S|值偏大的原因正是4.5节所指出的,先单站时间平均就等于将那些小量级降水站的高误差估算率保留下来,加之这些站数相对偏多,由此再进行区域平均时就自然偏大了小量级而高误差估算率的贡献。5 雷达和雨量计数据对的质量控制5.1 基本概念
雷达降水探测的质量控制分为两类,其中,第一类可简称为降水单变量的质量控制,如雷达反射率因子Z场的质量控制。其主要是利用Z场本身空间结构自然特性,去鉴别和剔除非降水回波(Steiner et al,2002;王佑兵等,2006; 庄薇等,2012),这对于除去Z场中的地物回波和超折射回波是相当有效的。第二类可简称为降水双变量的质量控制,利用空中Z场和地面降水QG场互相在对应点的结合关系来判断数据对(ZB,QG)的质量,这其中包括雷达和雨量计各自数据或者两者均存在问题时都将认为数据对是不可靠的,或者是错误的,经判断后应予以舍去,这就是所谓雷达和雨量计数据对的双变量质量控制。
Gabella等(2000)指出执行雨量计-雷达测量数据对质量控制的重要性。其原因首先在于几乎所有的降水估算技术方法都是从研究和开发这些数据对的对应关系脱胎而成,其次雷达降水误差估算的检验和评估也是从这些对应的关系计算出来的。一件有意义的事在于一组数据对(ZBi,QGi)通过第2节中式(1)的云雨转换方程就唯一地变为云雨转换系数ABi,于是ABi就含有这组数据对的空-地对应关系特性,于是雷达和雨量计数据对的双变量质量控制问题就变成为云雨转换系数ABi单变量的质量控制问题。5.2 误差估算因子控制法
万玉发等(2008)、Wan等(2010)中已经初步论述了(ZBi,QGi)数据对或ABi质量控制方法,并分为二级。第一级方法很简单,即将凡是1 h累积雨量 小于1.0 mm的数据对全部删去,其理由在原文中已有详述,简单可归结为在此条件下ABi的离散性非常大,几乎没有出现明显的聚合性,因此,将这样的数据对用作降水估算计算评估的基础明显是不可靠的,所以,应该舍去。第二级方法简单可叙述为将QG大于1.0 mm以上的雨量分为5个雨量子区间,根据一些降水估算实例,分别计算出各个雨量区间较为合适使用的AB可允许范围,将此暂时固定下来,用作二级质量控制的取舍标准。由此可见,这种方法不随天气过程和时次而变化显得较为机械和呆板。
在本节中将介绍一种与原有二级质量控制有所不同的一种新方法,称为单站误差估算因子μi控制法。前已提及雷达和雨量计数据对的质量控制可归结为云雨转换系数ABi的质量控制,而从第3节式(13)可知,当在固定指数bf的条件下以及当某时次ABS计算确定后,ABi与单站误差估算因子μi是唯一对应的,于是雷达和雨量计数据对的质量控制可归结为误差估算因子μi的质量控制。其基本操作方法是:必须设立一个较为合理的μi的允许区间(如μi∈[-0.8,1.5]),在μi值超出此允许区间的所有(ZBi,QGi)将被舍掉。余下的数据对重新融合得到新的ABS,建立新的区域降水估算方程。这样一来,ABi相对于ABS离散性将会变小而聚合性将会增大,区域平均误差估算因子μA和区域降水绝对误差估算率μ|S|等一般也随之变小。
不论是原有二级质量控制或现在的误差估算因子控制,从原理上都是将|μi|的大值舍掉,这就等价于将相对ABS太小的或太大的ABi舍去,这样做的理由在于当很小的ABi出现时这些雨量站点上必将出现淡云大雨,而反之就会出现浓云小雨的奇怪现象。出现这两种现象的站点一般不会是多数而是属于一种存在有疑点和不稳定的少数。其原因有多种可能,一种是由于雷达和雨量计存在时、空采样差异和分辨率差异而造成的,另一种可能是雷达或雨量计的观测错误或者是强度和位置的标定错误,当然也不能排除确实是大气中的自然现象,即使这样也会是极为罕见的。
综上所述,雷达和雨量计数据对的质量控制可形成一级控制、误差估算因子控制及双重控制3种控制。误差估算因子μi的允许区间可根据多个个例试验总结规律来设定。一旦设定了μi的允许区间,相应ABi的允许区间可随天气过程和时次而自动变化,这就是误差估算因子控制能适应降水天气变化的优越性。图 4a为未经质量控制数据ZB-QG点聚分布,图 4b为未经质量控制和在3种质量控制情况下的μi频率分布。从μi的频率分布可以看出,总体分布还是比较合理的,大都集中在μi=0值两旁;只是在μi=1.5后仍出现一定频率分布,特别在μi=2.0后出现一个较高频率状况,这显然是不合理的,因此,将μi的合理范围确定为-0.8—1.5。经过质量控制之后,μi=2.0后的高频率极端值去掉。
![]() |
图 4 2010年7月11日06—20时质量控制分析(a.未经质量控制数据ZB-QG点聚分布;曲线为等AB(或等μi)的ZB-QG关系线;b.未经质量控制和在3种质量控制情况下的μi频率分布) Fig. 4 Analysis chart of the quality control for the period from 06:00 BT to 20:00 BT 11 July 2010(a. ZB-QG scatter diagram without the quality control,each curve is the ZB-QG relationship line with equal AB(or equal μi); and b. μi the frequency distribution without the quality control and with the three kinds of the quality control) |
比较上述4种情况可以看出,对数据进行质量控制之后,μ|S|有所下降(表 2)。比较3种质量控制方式的μ|S|、μA可以看出,双重控制与单独一级(将凡是1 h累积雨量(QG)小于0.5 mm的数据对全部删去)或误差估算因子控制比较,μ|S|并不表现出明显的改善;而单独一级与误差估算因子比较,虽然μ|S|值后者改善不多,但μA值则从0.51下降到0.35。由此可见,在实际操作过程只需进行单独误差估算因子控制即可。
控制方式 | 总站数 | 平均雨量 | μ|S| | μA |
无控制 | 155 | 10.1 mm | 0.26 | 3.12 |
一级控制 | 136 | 11.5 mm | 0.23 | 0.51 |
误差估算因子 | 135 | 11.5 mm | 0.21 | 0.35 |
双重控制 | 129 | 12.1 mm | 0.20 | 0.33 |
本文对武汉雷达探测到的4次降水过程进行数据对质量控制前、后的降水估算结果对比评估。这4次降水过程是:(1)2007年7月1日18时—2日13时,湖北省东部有13个气象站出现了暴雨,8个气象站出现了大暴雨。(2)2008年6月8日12—21时,湖北省中东部出现大到暴雨天气,16个气象站降暴雨。(3)2008年7月2日01—15时,湖北省东北部出现暴雨天气,3个气象站降大暴雨。(4)2010年7月10日07时—11日04时,湖北省中东部地区降大到暴雨,3个气象站降大暴雨。使用武汉雷达探测半径200 km范围内雷达反射率因子和599个区域雨量站资料进行数据质量控制对比分析。评估结果见表 3,在没有开展数据对质量控制前,过程和区域降水总体绝对误差率μT|S|均高于30%,经过本研究提出的误差估算因子控制法处理后,过程和区域降水总体绝对误差率μT|S|平均下降7%。
个例1 | 个例2 | 个例3 | 个例4 | 平均 | |
无控制 | 0.36 | 0.34 | 0.34 | 0.39 | 0.36 |
控制后 | 0.30 | 0.27 | 0.28 | 0.30 | 0.29 |
(1)RASIM 构建了具有幂律形式的雷达和雨量计直接积分对应的单站关系式ZBi=ABi×QbfGi,这是一个准同雨团样本概念下的准有理对应关系,随时-空可变。其中系数ABi特称为单站云雨转换系数,它表明了在某时某站空中雨团ZB下落到地面雨量计成为QG的转换能力,其大小显然与天气系统和地理环境密切相关,值得今后深入研究。由单站云雨转换系数的数据集{ABi}S可以产生区域云雨转换公共系数ABS,从而构成区域型降水估算方程。这种用共性来代表个性的方法不可避免会带来一定的估算误差,而对于各个站点估算误差的表示和分析实质上是基于ABi的个性方程和基于ABS的共性方程的比较,使得误差分析具有清晰的物理意义和明确的数学表达式。
(2)单站误差估算因子μi是RASIM方法估算误差表示的最基本元素,以它为基础可构造出一系列误差表示和评估判据,如区域平均绝对误差估算因子μA、区域降水误差估算率μS、区域降水绝对误差估算率μ|S|以及区域平均绝对误差估算量EN等。在RASIM中μi正好等于ABi与ABS以bf-1为指数幂值的相对离差,从而揭示了引起降水估算误差的根本原因是单站ABi相对偏离于公共系数ABS所致。μA是区域中各μi绝对值的平均,其值大小体现了区域内各站降水条件差异的总体程度,与ABi相对ABS的聚合和离散密切有关,但μA与降水量大小无关。文中证明了只有在区域内出现均匀降雨的条件下,μA与μ|S|才相等。而μ|S|和EN既联系着μA,而且也与雨量有关,是较为合理而实用的降水估算误差判据。有意义的是,ABS区域型降水估算方程在区域内所有站点的高估总量必定与低估总量相等,这是一个十分突出的优点。
(3)RASIM区域型雷达降水估算方程中的估算参数ABS是经过对单站空-地对应关系式的面积分而得到的,在雷达覆盖区域内是一个公用常数,它可随时间而变化但不随站点而变化,换言之ABS只对时间变率具有一定的敏感度而对于空间变率即天气的地理分布差异是无力反映的。同时,雷达和雨量计是积分对应,模糊了不同性质降水的差异,会影响降水估算的精度。目前RASIM方法的估算精度在天气条件差异较小,即μA较小时,μ|S|可达到20%左右,而在天气条件差异较大即μA较大时,μ|S|在30%以上。因业务需要对武汉雷达8次和广州雷达4次降水过程进行了分雨量的评估,小时降水量0.5—5 mm段的μ|S|比5—20 mm段和20 mm以上段要大近5%,由此说明RASIM方法对小雨估算的准确性要差一些。为了对雷达估算降水作进一步的改进,必须研究区域云雨转换分布函数AB(x,y)S来替换公共系数ABS,从而建立RASIM体系下的局地型降水估算方程。
田付友, 程明虎, 张亚萍等. 2010. 校准雨量计密度对雷达联合雨量计估测流域平均面雨量的影响. 气象学报, 68(5): 717-730 |
万玉发, 吴翠红, 金鸿祥. 2008. 基于准同雨团样本概念雷达和雨量计的实时同步结合方法. 气象学报, 66(2): 262-273 |
王佑兵, 万玉发. 2006. 雷达体扫反射率场的自动质量控制. 气象科技, 34(5): 615-619 |
肖艳姣, 刘黎平, 杨洪平. 2008. 基于天气雷达网三维拼图的混合反射率因子生成技术. 气象学报, 66(3): 470-473 |
尹忠海, 张沛源. 2005. 利用卡尔曼滤波校准方法估算区域降水量. 应用气象学报, 16(2): 213-219 |
张培昌, 戴铁丕, 傅德胜等. 1992. 用变分方法校准数字化天气雷达测定区域降水量基本原理和精度. 大气科学, 16(2): 248-256 |
张培昌, 杜秉玉, 戴铁丕. 2001. 雷达气象学. 北京: 气象出版社, 511pp |
庄薇, 刘黎平, 余燕群等. 2012. 雷达地物回波模糊逻辑识别法的改进及效果检验. 气象学报, 70(3): 576-584 |
Ahnert P R, Krajewski W F, Johnson E R. 1986. Kalman filter estimation of radar-rainfall field bias//Preprints, 23rd Radar Meteorology Conference and Cloud Physics Conference. Amer Meteor Soc, JP33-JP37 |
Atlas D, Ulbrich C W. 1990. Early foundations of the measurement of rainfall by radar//Atlas D. Radar in Meteorology. Amer Meteor Soc, 86-97 |
Brandes E A. 1975. Optimizing rainfall estimates with the aid of radar. J Appl Meteor, 14: 1339-1345 |
Cheze J. 1999. The HYDRAM project of Meteo-France: A quality-based approach for quantative radar measurements//Preprints, 29th Conference on Radar Meteorology. Amer Meteor Soc, 321-324 |
Collier C G. 1984. The operational performance in estimating surface rainfall of a raingauge-calibrated radar system//22nd Conference on Radar Meteorology. Amer Meteor Soc, 257-262 |
Doviak R J, Zrnic D S. 2006. Doppler Radar and Weather Observations. New York: Dover Publications, INC, 223-235 |
Gabella M, Amitai E. 2000. Radar rainfall estimates in an Alpine environment using different gage-adjustment techniques. Phys Chem Earth (B), 25(10-12): 927-931 |
Joss J, Waldvogel A. 1990. Precipitation measurement and hydrology//Atlas D. Radar in Meteorology. Amer Meteor Soc, 577-606 |
Klazura G E, Kelly D S. 1995. A comparison of high resolution rainfall accumulation estimates from the WSR-88D precipitation algorithm with rain gage data//Preprints, 27th Conference on Radar Meteorology. Amer Meteor Soc, 31-34 |
Krajewski W F, Smith J A. 1991. On the estimation of climatological Z-R relationship. J Appl Meteor, 30(10): 1436-1445 |
Lack S A, Fox N I. 2004. Errors in surface rainfall rates retrieved from radar due to wind drift//Sixth International Symposium on Hydrological Applications of Weather Radar. 2004 |
Mittermaier M P, Hogan R J, Illingworth A J. 2004. Using mesoscale model winds for correcting wind-drift errors in radar estimates of surface rainfall. Quart J Roy Meteor Soc, 130(610): 2105-2123 |
Rosenfeld D, Wolff D B, Amitai E. 1994. The window probability matching method for rainfall measurement with radar. J Appl Meteor, 33(6): 682-693 |
Smith P L, Joss J. 1997. Use of a fixed exponent in “adjustable” Z-R relationships//Preprints, 28th Conference on Radar Meteorology. Amer Meteor Soc, 254-255 |
Steiner M, Smith J A, Uijlenhoet R, et al. 2001. Rainfall variability and its effect on rainfall measurement//Preprints, 30th Conference on Radar Meteorology. Amer Meteor Soc, 562-569 |
Steiner M, Smith J A. 2002. Use of three-dimensional reflectivity structure for automated detection and removel of non-precipitation echoes in radar data. J Atmos Ocean Technol, 19(5): 673-686 |
Wan Y F, Wu C H, Jin H X. 2010. Real-time synchronous integration of radar and raingauge measurements based on the quasi same-rain-volume sampling. Acta Meteor Sinina, 24(3): 340-353 |
Zawadzki I. 1975. On radar-raingage comparison. J Appl Meteor, 14: 1430-1436 |
Zawadzki I. 1984. Factors affecting the precision of radar measurement of rain//22nd Conference on Radar Meteorology. Amer Meteor Soc, 251-256 |
Zhang J, Howard D, Langston C, et al. 2011. National mosaic and multi-sensor QPE(NMQ) system: Description, results, and future plans. Bull Amer Meteor Soc, 92(10): 1321-1338 |