中国气象学会主办。
文章信息
- 李巧, 戚友存, 朱自伟, 杨毅, 闵锦忠, 师春香, 张哲, 李东欢, 王楠, 胡启元. 2021.
- LI Qiao, QI Youcun, ZHU Ziwei, YANG Yi, MIN Jinzhong, SHI Chunxiang, ZHANG Zhe, LI Donghuan, WANG Nan, HU Qiyuan. 2021.
- 复杂地形下C波段雷达定量降水估计算法
- Quantitative precipitation estimation algorithm for C-band radar situated in complex topographical regions
- 气象学报, 79(4): 689-702.
- Acta Meteorologica Sinica, 79(4): 689-702.
- http://dx.doi.org/10.11676/qxxb2021.038
文章历史
-
2020-10-10 收稿
2021-04-10 改回
2. 中国科学院地理科学与资源研究所陆地水循环实验室,北京,100101;
3. 中国科学院大学,北京,100864;
4. 南京信息工程大学大气科学学院,南京,210044;
5. 国家气象信息中心气象数据研究室,北京,100081;
6. 陕西省气象局,西安,710014
2. The Key Laboratory of Water Cycle and Related Land Surface Processes,Institute of Geographic Sciences and Natural Resources Research,Chinese Academy of Sciences,Beijing 100101,China;
3. The University of Chinese Academy of Sciences,Beijing 100864,China;
4. School of Atmospheric Sciences,Nanjing University of Information Science and Technology,Nanjing 210044,China;
5. Meteorological Data Research Office,National Meteorological Information Center,Beijing,100081,China;
6. Shaanxi Provincial Meteorological Service,Xi’an 710014,China
降水对气象、水文、地质灾害、农业畜牧业生产以及其他诸多科学和社会应用起着至关重要的作用,但是研发高精度、高时空分辨率降水预报产品具有极高的挑战性。众所周知,地面雨量站观测可以准确获取雨量计所在位置的降水量,但是基于地面雨量站观测的降水插值产品精度严重受到地面雨量站点的疏密和空间分布影响(沈艳等,2012;肖艳姣等,2011;Morrissey,et al,1995;Villarini,et al,2008)。在全球范围内,目前大部分地区地面雨量站的布设密度不足以捕捉到完整的中尺度系统引起的降水分布。对于经常诱发山洪地质灾害的对流尺度强降水,即便在发达国家,地面雨量计分布密度仍无法捕捉对流尺度强降水空间分布。静止气象卫星观测可以覆盖全球,时间分辨率也较高,但由于距离遥远,观测精度有一定的局限(高晓荣等,2013)。
为了弥补地面雨量计和卫星观测的不足,中国近20年来布设了220部CINRAD天气雷达。目前中国雷达站网布设密集,多部雷达组网拼图可有效扩大雷达观测范围,因此,利用雷达资料进行定量降水估计(QPE)获取地面降水信息具有极大的优势(王红艳等,2009)。基于雷达扫描观测特征,天气雷达可以精确地捕捉地面降水的时、空分布。天气雷达最大观测半径可达450 km,并且,单雷达的空间分辨率可以达到1 km × 1 km,时间分辨率可达5 min(Zhang,et al,2011,2016)。然而,雷达定量降水估计存在雷达固有的观测特征、反演技术和空间地理等方面的问题,包括地形、建筑物和地物的遮挡,降水粒子垂直方向的不均匀,降水粒子滴谱空间分布不确定等(Zhang ,et al,2010,2011,2016;Qi,et al,2013a,2013b)。本研究主要从雷达探测资料质量控制、降水类型分类、雷达波波束遮挡,Z-R关系和雷达拼图5个方面来提高雷达定量估计降水的精度。
天气雷达接收的回波信号除了气象回波,还包括地物回波、电磁干扰回波等非气象回波,这些非气象回波会严重影响雷达定量估计降水的精度。因此,在应用雷达资料前,需对其进行质量控制,剔除非气象回波,以提高定量降水估计精度。早期很多学者对雷达探测资料质量控制算法开展了大量研究工作,根据降水的时、空连续性,提出通过检查雷达回波的空间和时间连续性(Smith,1990),分析反射率因子的水平梯度和垂直梯度(Mueller,1977),对雷达探测资料进行质量控制,很大程度上去除了雷达非气象回波。近些年,很多学者利用概率分析方法进行雷达探测资料质量控制,以去除雷达非气象回波。Kessinger等(2003)提出了采用模糊逻辑算法的雷达回波分类技术,从雷达资料(回波强度、径向速度和速度谱宽)中提取特征场,用于区分降水回波、地物和海浪回波等不同类型的雷达回波。刘黎平等(2007)基于Kessinger等(2003)的模糊逻辑算法,提出了分步式识别地物回波的方法,对SA波段雷达进行了质量控制,改进了地物回波识别的效果。江源(2013)在已有的模糊逻辑算法识别地物回波基础上,引入人工智能识别方法,有效剔除了电磁干扰回波。本研究主要利用已有的模糊逻辑质量控制算法(江源,2013),对C波段雷达进行非气象回波识别及去除,以提高雷达定量降水估计精度。
在雷达扫描平面内,若地形复杂,雷达电磁波波束会受到遮挡,这就会给雷达定量降水估计引入误差。为了去除地形遮挡导致的雷达定量降水估计的低估误差,学者(Fulton,et al,1998;Smith,1998)提出了根据雷达不同方位角波束遮挡情况,选择雷达复合平面扫描方法进行定量降水估计。雷达复合平面扫描方法已经业务应用于美国雷达定量降水估计系统(NMQ,Zhang,et al,2011),有效提高了雷达监测降水的能力,提高了雷达定量估测降水的精度。肖艳姣等(2008)将雷达复合平面扫描方法引入到中国S波段雷达定量降水估计中,较大程度上提高了雷达定量降水估测的精度。本研究在已有研究基础上,基于陕西C波段雷达研发雷达复合平面扫描技术,以解决低仰角雷达电磁波受到遮挡的问题。
基于计算得到的混合仰角反射率因子(以下简称混合仰角反射率),如何反演得到准确的降水率是雷达定量降水估计要解决的关键问题。目前雷达定量降水估计最常用的反演方法是基于单一雷达气象方程Z-R关系法(张培昌等,2001)来计算降水率。但本质上,降水率是由降水粒子滴谱决定的。Rosenfeld等(2003)研究了不同区域、不同类型的降水,通过对降水粒子滴谱数据统计分析得到200多种不同的Z-R关系,这说明单一的Z-R关系难以准确描述由于降水过程、降水类型和地区不同导致的降水滴谱特征差异,针对降水粒子滴谱特征不同使用相应的Z-R关系对于提高雷达定量降水估计的准确性具有重要意义。实际情况是,同一个雷达扫描平面内观测到的降水信息可能来自不同的降水类型,使用单一的Z-R关系会给雷达定量降水估计引入偏高或偏低的误差(Balakrishnan,et al,1989;Lee,et al,2006)。因此,需要针对不同类型降水采用不同的Z-R关系将雷达反射率转换成降水率,以提高雷达定量估计降水精度。文中基于雷达观测提供的体扫数据,结合不同类型降水的水平和垂直结构特征,研发降水类型分类模块,识别出不同类型降水。针对不同的降水类型选用具有代表性的Z-R关系将反射率转换成降水率,以提高雷达定量降水估计精度。
由于单部雷达定量降水估计受到雷达本身探测能力的限制,并且雷达波束随径向距离增大而升高,导致径向距离越远,观测信息对于地面降水的代表性往往越低。利用多部雷达降水观测信息进行相互补充,可以提供比单部雷达降水观测范围更广、精度更高的观测信息。因此,需要通过拼图的方式把多部雷达的降水观测信息组合到一起。目前,针对雷达拼图算法已经开展了很多研究,可以总结为以下两种方式:(1)先用单部雷达估测降水,再进行雷达组网拼图(Tabary,2007;高晓荣等,2012);(2)基于雷达反射率场拼图,然后进行降水估测(肖艳姣等,2008;Zhang,et al,2011)。文中将先用单部雷达观测资料进行降水估测,然后再对单雷达降水估测场进行拼图。相较于高晓荣等(2012)的方法,文中基于降水结构及降水粒子垂直相态变化特征,确定单雷达在拼图时的权重系数,最终将单雷达观测的降水信息高效地组合到一起。
文中将从雷达探测资料质量控制、降水类型分类、雷达混合仰角反射率计算、Z-R转换关系和雷达拼图5个方面介绍C波段雷达定量降水估计方面的研究工作,最后通过对比地面雨量站降水观测产品,定量评估新开发C波段雷达定量降水估计算法的性能。
2 雷达定量降水估计算法选取陕西省为研究区域(图1)。陕西省地势南北高、中间低。北部为黄土高原,平均海拔1000 m左右,占全省面积的40%。中部为关中平原,平均海拔500 m左右。南部为秦巴山区,包括秦岭-大巴山脉以及汉江谷地,平均海拔1500 m以上。高原、山地、平原和盆地等多种地形错综复杂,降水导致的滑坡、泥石流等地质灾害经常发生。因此,提高定量降水估计产品的精度对于陕西省气象地质灾害的预警、防灾和减灾都具有重要意义。
所用降水量数据为陕西地面雨量站1 h降水观测数据,陕西省西安、延安、榆林、安康、商洛、宝鸡、汉中7部C波段天气雷达原始体扫数据。图1a给出了陕西省的地形空间分布以及7部C波段天气雷达和地面雨量站的分布情况。基于地形空间分布,结合雷达扫描特征,可以计算出每部天气雷达的最大有效探测降水距离,图1b给出了7部雷达探测降水的最大有效距离。7部雷达从北到南依次为榆林、延安、宝鸡、西安、商洛、汉中、安康。陕西省C波段雷达波束径向探测距离为400 km,径向分辨率为500 m,方位角分辨率及波束宽度均为1°。体扫模式为VCP21,即一个体扫包括9个仰角的扫描,仰角分别为0.5°、1.45°、2.4°、3.35°、4.3°、6°、9.9°、14.6°、19.5°,完成一次体扫的时间大约是6 min。7部C波段天气雷达和地面雨量站基本可覆盖全省。
雷达定量降水估计算法包括以下5个相对独立的模块,图2给出C波段天气雷达定量降水估计算法流程,主要包括:(1)雷达探测资料质量控制,(2)雷达降水类型划分,(3)雷达混合仰角反射率计算,(4)Z-R关系统计分析,(5)雷达拼图。
2.1 雷达探测资料质量控制地物、晴空回波和海浪回波会影响雷达反射率因子数据质量,从而影响雷达定量降水估计精度,并导致风暴监测或跟踪失败。因此,雷达反射率因子数据在应用之前,有必要对其进行质量控制。对雷达基数据采用刘黎平等(2007)提出的模糊逻辑分布式地物回波识别方法对地物进行识别和剔除。模糊逻辑方法是区分降水回波、地物回波、晴空回波和海浪回波等物理量,然后根据这些回波的特征建立隶属函数(刘黎平等,2007),对这些物理量进行特征值计算,得到所有物理量对于不同类型回波的0—1取值范围的判据。回波点物理量对应的判据值越大,那么该回波点属于该类型回波的可能性越大。对这些判据值进行加权累加,如果回波点的地物回波判据超过事先设定的阈值,那么,该回波点的回波将被识别成地物回波。具体计算方法与Kessinger等(2003)的方法相同,首先,计算反映地物和降水回波差异的4个物理量:(1)回波强度纹理(TdBZ),(2)回波垂直变化(GdBZ),(3)回波沿径向变号(SIGN),(4)回波沿径向变化程度(SPIN)。其次,计算反映地物和降水差异的径向速度和谱宽的3个物理量:(1)径向速度的区域平均值(MDVE),(2)径向速度的方差(SDVE),(3)速度谱宽的区域平均值(MDSW)。这些物理量的具体定义见刘黎平等(2007)。
图3给出了2019年4月28日03时24分(世界时,下同)质量控制前(图3a)、后(图3b)西安雷达0.5°仰角反射率。从图3a可以看到,在雷达附近有明显的地物回波,在40°左右有明显的尖峰(spike),质量控制之后(图3b),这些地物回波和尖峰基本被准确识别出来,降水区域的信息基本保留。结果表明:该方法对地物回波和降水回波区分效果较好,多数地物回波能够被有效滤除,甚至某些距离模糊区,速度不可用的位置,地物也被识别出来。需要注意的是,在降水系统边缘,部分气象回波较弱,被误识别为非气象回波从而被滤除。
2.2 雷达降水类型分类雷达降水类型分类算法基于单雷达体扫降水观测资料,自动识别降水类型,本研究中将降水类型分为对流性降水(简称CV),不受0°C层亮带影响的层状云降水(简称ST)和受0°C层亮带影响的层状云降水(简称BB)。降水类型识别是提高雷达定量降水估计精度的重要前提之一,正确的降水分类使得定量降水估计时可以针对不同降水类型引入更具代表性的Z-R关系,以有效提高定量降水估计的精度。
降水类型分类模块主要是利用质量控制后的雷达反射率数据,对不同降水类型进行识别,得到降水类型产品。
2.2.1 0°C层高度检查在进行对流核识别前,先检查0°C层高度,该信息可从模式数据或探空资料获取。文中应用的二维0°C层高度数据从美国国家环境预报中心的全球预报系统(GFS)GFS-ANL资料中获取,该资料空间分辨率为1°×1°,时间分辨率为6 h。当0°C层高度过低(低于雷达高度2 km),则认为没有对流产生,这样的设定是基于如果地面温度较低,对流不容易产生的前提,避免了把低高度的0°C层亮带误识别为对流而造成定量降水估计高估的现象。
2.2.2 对流核识别对流核识别算法是基于单雷达体扫数据,在极坐标下实现。算法考虑了地面雷达探测的特性,在不同区域采用不同的判据来识别对流核。识别对流核基于多个仰角的反射率和垂直液态水含量等,避免了由于仅判断特定层次的反射率或者组合反射率(用体积扫描获取的回波强度数据,在每个方位-斜距库上,对所有仰角的反射率因子在垂直方向进行比较,挑选出最大的反射率因子)而把强0°C层亮带误识别为对流核的问题,具体算法详见Qi等(2013b)和Zhang等(2021)。
2.2.3 对流区识别在识别出所有对流核之后,采用区域增长方法来识别整个对流区。区域增长方法通过检查对流核格点四周的8个格点,如果这8个格点中的任何一个格点满足下述特定条件,则将其标记为对流区,并以其为新的对流核继续检查其邻近的8个格点,直到所有的对流核都被检查完毕。在检查完所有格点后,满足条件的格点即为算法识别出的对流区。
降水类型分类算法中综合组合反射率、最大反射率高度、反射率垂直梯度这3个物理量进一步判断识别出的对流区是否正确。如果该格点组合反射率超过35 dBz且满足下面条件之一,则判断该格点为对流:(1)组合反射率大于45 dBz(0℃层亮带反射率一般不超过45 dBz),(2)最大反射率高度不在0℃层亮带影响的高度(说明高反射率值不是由于融化效应所导致的),(3)反射率垂直梯度小于4 dBz/km(对流在垂直方向上相对均匀,反射率垂直梯度较小)。
2.2.4 0℃层亮带区识别在识别出整个对流区之后,其余非对流降水区域即为层状云降水区,层状云降水区存在0℃层亮带型和非0℃层亮带型降水。0℃层亮带识别算法将首先识别0℃层亮带核,即在层状云区中雷达组合反射率超过35 dBz的格点。然后,通过区域增长法,基于0℃层亮带核识别出整个0℃层亮带区,该区域定义为从最高仰角接触到0℃层高度的位置到最低仰角接触到0℃层高度的位置,区域增长法的判定条件为组合反射率超过30 dBz。
图4给出了雷达组合反射率(图4a)和对应的降水类型分类结果(图4b)。可以看到,该算法很好地将飑线和飑线后方的强0℃层亮带区识别出来,而在雷达北部的小范围对流区也被识别出来,即降水类型分类算法能够较好地识别对流性降水,0℃层亮带型层云降水和非0℃层亮带层状云降水。
2.3 雷达混合仰角反射率计算在利用雷达反射率进行定量降水估计时,当雷达探测到的信息离地面越近,对于地面真实的降水强度描述就越具备代表性,故使用最贴近地面的雷达仰角观测反射率信息进行降水估计更为准确。然而在实际操作中,低仰角在某些方位角上因受到地形遮挡而无法获得降水信息。为得到完整有效的降水信息,需要利用更高仰角的信息进行补充。
在雷达复合平面扫描中,高精度地形数据是美国航空航天局(NASA)和国防部国家测绘局(NIMA)联合测量的机载雷达地形探测任务(SRTM)数字高程模型(DEM),空间分辨率为30 m×30 m。雷达复合平面扫描方法主要是利用高精度地形数据,结合雷达探测原理,计算每部雷达电磁波遮挡百分比,从而得到每部雷达的有效最低混合仰角。计算雷达电磁波能量遮挡百分比,需要确定雷达电磁波在传播过程中,雷达波束对应位置的地形高度空间分布。假设雷达电磁波在标准大气中传播,基于等效地球半径模型(Langston,et al,2007;Zhang,et al,2005)可以得到
${h_{\rm{c}}} = {\left({R_{\rm{e}}^2 + {r^2} + 2{R_{\rm{e}}}r\sin t} \right)^{1/2}} - {R_{\rm{e}}}$ | (1) |
${R_{\rm{e}}} = \frac{\,4\,}{\,3\,}{R_{\rm{0}}}$ | (2) |
式中,
此外,由于雷达波束随着传播距离的增大而变宽,设定波束宽度为
$d = 2r\sin \frac{\,\theta\, }{\,2\,}$ | (3) |
根据雷达电磁波能量在垂直截面上满足高斯分布,并结合SRTM数字高程模型计算
$P = \int_0^d {\int_{{h_{\rm{b}}} - {h_{\rm{c}}}}^{y\left(x \right) - {h_{\rm{c}}}} {f\left({x,y} \right)} {\rm{d}}x{\rm{d}}y} $ | (4) |
$f\left({x,y} \right) = \frac{1}{{2{{\text π }}{{\left| \boldsymbol\varSigma \right|}^{{1 / 2}}}}}\exp \left({ - \frac{1}{2\boldsymbol\varSigma}{{\left({{ X} - \overline { X} } \right)}^{\rm T}}\left({{ X} - \overline { X} } \right)} \right)$ | (5) |
${\boldsymbol X} = \left[ \begin{gathered} x \\ y \\ \end{gathered} \right],\quad\overline {\boldsymbol X} = \left[ \begin{gathered} \overline x \\ \overline y \\ \end{gathered} \right],\quad{{\boldsymbol\varSigma}} = \left[ {\begin{array}{*{20}{c}} {\sigma _{{x}}^2}&0 \\ 0&{\sigma _{{y}}^2} \end{array}} \right]$ | (6) |
式中,
当地形对雷达电磁波能量遮挡超过50%时,认为雷达电磁波受地形遮挡严重,雷达观测的降水信息缺失严重。此时,依次查看高层仰角在相同方位角、相同径向距离上的能量遮挡百分比,找到地形对雷达电磁波能量遮挡小于50%的最低仰角,即为雷达在该方位角和径向距离上的最低有效仰角。同时,对于最低有效仰角上的反射率,按照能量遮挡百分比进行能量补偿。由于雷达平面内地形的遮挡情况分布可能不一致,不同位置对应的最低有效仰角可能不同,因此将最低有效仰角在雷达平面上的分布称为最低有效混合仰角,在最低有效混合仰角上得到的反射率数据称为混合仰角反射率。
图5给出了延安雷达的地形空间分布(图5a),雷达0.5°仰角(图5b)、1.45°仰角(图5c)的电磁波遮挡百分比空间分布以及得到的雷达混合仰角产品(图5d)。从图中可以看出,延安雷达周围地形复杂,在0.5°仰角,雷达西南、东南、西北偏北方向电磁波均受到严重遮挡,而在1.45°仰角,雷达仅在西北方向受到部分遮挡,根据雷达在不同仰角不同方位角的地形遮挡情况生成的混合仰角则综合了各个方位角的最优反射率情况,可以较全面地获得降水信息。
2.4 基于GPM星载双频雷达(DPR)观测统计分析Z-R关系Zhu等(2020a,2020b)利用全球降水测量任务(GPM)星载双频雷达(DPR)对陕西省降水时、空变化特征的分析得到了以下结论:(1)不同降水类型的垂直结构(微物理过程)具有明显的差异,(2)相同降水类型和强度的垂直反射率廓线在陕西省不同季节和不同地区的差异并不显著。因此,利用GPM卫星雷达反演的滴谱数据对不同降水类型的Z-R关系进行了统计分析,文中所用GPM卫星数据从地球数据网站(https://storm.pps.eosdis.nasa.gov/storm/)获取。
在利用卫星雷达的近地面反射率和降雨率进行Z-R关系拟合之前,需要将卫星雷达Ku波段的反射率转换为与地面雷达相同频段(C波段)的反射率(Zhu,et al,2020b),这样形成的Z-R关系才能直接应用于地面雷达的定量降水估计。雷达的反射率可以用以下公式表示(Zhang,et al,2001)
$Z{\left(\lambda \right)_{\rm{h}}} = \frac{{4{\lambda ^4}}}{{{{{\text π }}^4}{{\left| {K_{\rm w}} \right|}^2}}}\int_0^{ + \infty } {{{\left| {{F_{{\rm{hh}}}}} \right|}^2}N\left(D \right){\rm{d}}D} $ | (7) |
${K_{\rm{w}}} = ({\varepsilon _{\rm{w}}} - 1)/({\varepsilon _{\rm{w}}} + 2)$ | (8) |
$\begin{split}\int_0^{ + \infty } {{{\left| {{F_{{\rm{hh}}}}} \right|}^2}N\left(D \right){\rm{d}}D} =& \left({1 - 2{\sigma ^2}} \right)\int_0^{ + \infty } {{{\left| {{f_{\rm{a}}}} \right|}^2}N\left(D \right){\rm{d}}D} +\\& 2{\sigma ^2}\int_0^{ + \infty } {{f_{\rm{a}}}f_{\rm{b}}^ * N\left(D \right){\rm{d}}D} \end{split}$ | (9) |
$N\left(D \right) = {N_{\rm{w}}}f\left(\mu \right){\left({\frac{D}{{{D_{\rm{0}}}}}} \right)^\mu }\exp \left({ - \varLambda D} \right)$ | (10) |
$f\left(\mu \right) = \frac{6}{{{{3.67}^4}}}\frac{{{{\left({3.67 + \mu } \right)}^{\mu + 4}}}}{{\varGamma (\mu + 4)}},\quad\varLambda = \frac{{3.67 + \mu }}{{{D_{\rm{0}}}}}$ | (11) |
式中,
通过改变式(10)中滴谱参数
通过图6b形成的转换关系,选取2014年3月—2019年10月GPM卫星雷达在陕西省内的观测样本,将GPM卫星雷达Ku波段在不同降水类型中观测到的近地面反射率转换为C波段对应的反射率,结合由GPM卫星雷达双波段Ku和Ka得出的雨滴谱估计的地面降雨率,得到C波段雷达定量降水估计的Z-R关系如图7所示。可见,用Z-R关系进行雷达定量降水估计总体上较为准确,尽管对于中等强度降水估计的准确度会出现一定的下降。由于降水率本质上由雨滴谱决定,这里也隐含着中等强度降水雨滴谱的不确定性更大。
2.5 基于雷达观测原理的雷达拼图基于上述模块生成的雷达小时降水率产品和降水类型产品均为单部雷达产品,为便于雷达产品的应用,通常需要将雷达单站产品进行拼图,最终以规则的经纬度网格表示。雷达拼图算法主要是对陕西省7部雷达的雷达反射率、小时降水率和降水类型分类场进行处理,从而得到空间范围(31.01°—40.00° N,104.00°—113.00° E),空间分辨率为0.01°×0.01°,时间分辨率为6 min的拼图产品。
对各网格单元而言,拼图需要确定天气雷达网中对其有贡献的雷达及其在雷达平面上对应的位置,即根据网格和雷达站点经纬度计算网格单元在雷达平面上的方位角和径向距离。两点经纬度坐标与两点间距离与方位角的转换关系通常用大圆模型表达(Langston,et al,2007)。
雷达反射率场和小时降水率产品的拼图方法是对各网格单元有贡献的单站雷达反射率和小时降水率进行加权平均,n部雷达各自的权重
${k_i} = \frac{{1000}}{{{h_{{\rm c},i}}}} + \frac{{1000}}{{{d_i}}}$ | (12) |
${w_i} = \frac{{{k_i}}}{{\displaystyle\sum\limits_{i = 1}^n {{k_i}} }}$ | (13) |
式中,
图8是2018年8月21日11时陕西省7部天气雷达混合仰角反射率(图8a),降水类型分类场(图8b),1 h雷达降水率场(图8c)的拼图结果。可以看出,通过赋给雷达不同权重,保证了近距离和低高度的雷达做出更多贡献,保留了反射率场和1 h降水率场的更多细节,使得拼图没有明显的不连续情况,并且能够显示明显的降水强中心,降水类型分类场中对流降水、0℃层亮带降水和非 0℃层亮带层状云降水与反射率场、1 h降水率场基本对应。表明雷达拼图算法能够较好地结合多部雷达,提供较为完整准确的降水信息。
3 结果分析从2019年5月26日14时—29日02时地面雨量站和雷达定量降水估计算法得到的累计降水结果(图9)可以看到,同雨量站累计降水(图9a)相比,定量降水估计算法(图9b)基于最低混合有效仰角对不同的降水类型使用不同Z-R关系进行反演,能够较好地体现降水过程的降水中心、降水范围,对于降水强度的把握比较接近雨量站观测。需要注意的是,陕西雷达仍存在一些非地形遮挡现象,这些遮挡情况无法通过地形数据计算得到,导致部分雷达(如汉中)仍存在部分方位角降水信息不完整的问题。
为了更系统地评估文中雷达定量降水估计算法的性能,利用下面3个统计评分指标进行评估:
(1)均方根误差(RMSE)
${\rm{RMSE}} = {\left[\frac{1}{N}\sum\limits_{k = 1}^N {{{({r_{\rm{k}}} - {g_{\rm{k}}})}^2}} \right]^{1/2}}$ | (14) |
(2)相对误差(RMAE)
${\rm{RMAE}} = \tfrac{{\displaystyle\sum\limits_{k = 1}^N {\left| {{r_{\rm{k}}} - {g_{\rm{k}}}} \right|} }}{{\displaystyle\sum\limits_{k = 1}^N {{g_{\rm{k}}}} }}$ | (15) |
(3)相对偏差(RMB)
${\rm{RMB}} = \tfrac{{\displaystyle\sum\limits_{k = 1}^N {\left({{r_{\rm{k}}} - {g_{\rm{k}}}} \right)} }}{{\displaystyle\sum\limits_{k = 1}^N {{g_{\rm{k}}}} }}$ | (16) |
式中,
图10为2019年5月26日14时—29日02时累计降水与地面雨量站观测降水散点对比。可以看到,雷达定量降水估计算法基于混合最低有效仰角进行了降水类型分类,并根据降水类型采用相应的Z-R关系,得到的雷达定量估计降水结果其整体趋势基本符合理想曲线,从统计评分指标(RMSE、RMB)也能看出,雷达定量降水估计算法同地面雨量站观测降水较为接近。
以地面小时雨量站观测为基准,在陕西省选取2019年4—7月持续时间较长、影响范围较大的7次降水过程对文中算法生成的雷达QPE产品进行评估。表1给出了这7次个例的具体信息。图11是这7次个例两种降水产品的统计评分指标。图11a中,红点为该次过程中雨量计1 h平均降水量,蓝点为该次过程中雷达QPE产品1 h平均降水量。可以看出,文中算法生成的雷达QPE产品同地面雨量站观测降水较为接近,特别是过程3和4,其平均偏差接近于0。总体上,文中算法生成的雷达QPE产品有一定的低估(RMB为负值),但基本控制在−30%以内,过程5和6达到−40%(图11c)。这说明该算法生成的雷达QPE产品在不同过程中的表现比较稳定,这一点从RMAE指标(RMAE稳定在50%左右,图11b)也可以反映出来。
个例 | 日期-时间 | 个例介绍 |
1 | 04-07 14:00—04-09 04:00 | 降水主要集中在陕西北部及南部地区,其中南部地区主要为对流性降水 |
2 | 04-26 01:00—04-28 09:00 | 飑线过程,自西北向东南移动,逐渐减弱成团状对流系统,降水主要集中在陕西中部及南部地区 |
3 | 05-05 15:00—05-08 22:00 | 降水主要集中在陕西中部及南部地区,其中中部地区有对流性降水过程 |
4 | 05-26 14:00—05-29 02:00 | 降水主要集中在陕西中部及南部地区,有对流性降水过程 |
5 | 06-26 14:00—06-28 14:00 | 全省范围内均有降水,其中南部地区主要为对流性降水 |
6 | 07-20 14:00—07-23 01:00 | 降水主要集中在陕西北部、中部和西南部地区,为对流性降水 |
7 | 07-28 14:00—07-31 08:00 | 全省范围内均有降水,其中北部和南部地区主要为对流性降水 |
总计 | 406 h | 7次个例 |
从均方根误差上看,文中算法生成的雷达QPE产品同地面雨量站观测降水较为接近,7次过程的均方根误差都在3 mm以下,过程3和4的均方根误差小于1 mm(图11a),雨量计1 h平均降水量同雷达QPE产品1 h平均降水量也比较接近。从相对误差看,雷达QPE产品较雨量计观测降水误差基本在60%左右。不同过程之间的差异主要是降水平均强度差异导致的。7次降水过程中,过程6和7的降水强度较大,而对应的均方根误差相对较大(3 mm左右),过程6的相对偏差也在−30%左右,表明对于中等及以上强度的降水,文中定量降水估计算法估计的准确度会有一定下降,但总体而言,文中算法仍较为稳定,生成的雷达QPE产品与地面雨量站降水观测相比较为接近。
4 结论与讨论基于陕西省7部C波段天气雷达观测数据,结合雷达探测资料质量控制、降水类型分类、地形遮挡、Z-R关系和雷达拼图等提出了C波段天气雷达定量降水估测(QPE)算法,生成空间分辨率为0.01°(约1 km),时间分辨率为6 min的雷达QPE拼图产品。通过地面小时雨量站的观测数据,对新的雷达QPE产品进行了评估和比较,主要得到以下结论:
(1)雷达探测资料质量控制算法对地物回波、尖峰等非降水回波和降水回波识别效果较好,能够较好地保留降水信息,但对于降水系统边缘较弱的气象回波有误识别现象,算法需要进一步改进。
(2)通过雷达降水类型分类算法,能够较准确地识别对流性降水、非0℃层亮带层状云降水和0℃层亮带层状云降水。
(3)根据陕西地区雷达遮挡情况,设计了混合最低有效仰角反射率算法,可有效降低地形遮挡导致的降水低估或缺失现象。
(4)利用GPM卫星雷达反演的滴谱数据对不同类型降水的Z-R关系进行了统计分析,针对不同降水类型采用对应的Z-R关系有利于提高QPE精度。
(5)文中选取了7次降水过程,对新算法生成的雷达QPE产品进行了评估,雷达算法相对稳定,相对误差、相对偏差指标也均表明该算法结果与地面雨量站降水观测结果接近。
与其他地区相比,目前陕西地区雷达QPE结果和雨量计降水观测之间的偏差依然较大,主要可能有以下几个原因:(1)雨量计质量问题;(2)雷达受地形阻挡探测不到;(3)雷达电磁波部分受地形遮挡导致能量削弱;(4)陕西雷达存在一定的非地形遮挡现象;(5)部分地区可能存在地形对降水的增强,导致在近地面降水具有较大增幅,尤其是强的对流性降水,使得雷达QPE偏低。另外,虽然C波段电磁波衰减远小于Ku波段,但C波段天气雷达的水平观测距离为近百千米,而GPM星载双频雷达(DPR)Ku波段在垂直方向只有几千米雨区,因此C波段天气雷达的衰减可能远大于DPR Ku波段的衰减,从而导致雷达QPE偏低,尤其是有较大范围强降水和强冰雹时,C波段雷达衰减问题会更加严重。同时,由于DPR数据本身也存在误差,例如由于地物杂波、非均匀充塞导致的Z值观测误差,虽然在算法中对这些误差进行了订正,但不能完全消除;此外,根据滴谱参数反演的R值也存在一定误差。由于陕西地形复杂,雷达电磁波遮挡问题较为严重,因此,仅依靠雷达观测数据,在地形遮挡严重的区域难以获取较为准确的降水量信息。未来应结合卫星数据和质量控制后的地面雨量站数据对雷达受遮挡范围进行补充,进一步提高QPE精度。
目前,中国正在进行双偏振雷达的升级工作,在未来几年中,陕西省雷达也会逐渐升级为双偏振雷达,双偏振雷达提供的多个偏振量能够反映降水的滴谱信息,从而更好地描述降水特征,引入偏振量信息,将促进定量降水估计精度的提高。
高晓荣, 梁建茵, 李春晖. 2012. 雷达定量降水估计技术及效果评估. 热带气象学报, 28(1): 77-88. Gao X R, Liang J Y, Li C H. 2012. Radar quantitative precipitation estimation techniques and effect evaluation. J Trop Meteor, 28(1): 77-88. DOI:10.3969/j.issn.1004-4965.2012.01.009 (in Chinese) |
高晓荣, 梁建茵, 李春晖等. 2013. 多平台(雷达、卫星、雨量计)降水信息的融合技术初探. 高原气象, 32(2): 549-555. Gao X R, Liang J Y, Li C H, et al. 2013. Preliminary studies on merged techniques based on precipitation information from multiplatform (radar, satellite and rain gauge). Plateau Meteor, 32(2): 549-555. (in Chinese) |
江源. 2013. 天气雷达观测资料质量控制方法研究及其应用[D]. 北京: 中国气象科学研究院. Jiang Y. 2013. Meteorological radar data quality control study and application[D]. Beijing: Chinese Academy of Meteorological Sciences(in Chinese)
|
刘黎平, 吴林林, 杨引明. 2007. 基于模糊逻辑的分步式超折射地物回波识别方法的建立和效果分析. 气象学报, 65(2): 252-260. Liu L P, Wu L L, Yang Y M. 2007. Development of fuzzy-logical two-step ground clutter detection algorithm. Acta Meteor Sinica, 65(2): 252-260. DOI:10.3321/j.issn:0577-6619.2007.02.011 (in Chinese) |
沈艳, 潘旸, 徐宾等. 2012. 最优插值法在对中国自动站降水量空间分析中的参数优化. 成都信息工程学院学报, 27(2): 219-224. Shen Y, Pan Y, Xu B, et al. 2012. Parameter improvements of hourly automatic weather stations precipitation analysis by optimal interpolation over China. J Chengdu Univ Inf Technol, 27(2): 219-224. (in Chinese) |
王红艳, 刘黎平, 王改利等. 2009. 多普勒天气雷达三维数字组网系统开发及应用. 应用气象学报, 20(2): 214-224. Wang H Y, Liu L P, Wang G L, et al. 2009. Development and application of the doppler weather radar 3-D digital mosaic system. J Appl Meteor Sci, 20(2): 214-224. DOI:10.3969/j.issn.1001-7313.2009.02.011 (in Chinese) |
肖艳姣, 刘黎平, 杨洪平. 2008. 基于天气雷达网三维拼图的混合反射率因子生成技术. 气象学报, 66(3): 470-473. Xiao Y J, Liu L P, Yang H P. 2008. Technique for generating hybrid reflectivity field based on 3-D mosaicked reflectivity of weather radar network. Acta Meteor Sinica, 66(3): 470-473. (in Chinese) |
肖艳姣, 万玉发, 王珏等. 2011. 改进的C波段天气雷达定量降水估计算法研究. 暴雨灾害, 30(4): 321-327. Xiao Y J, Wan Y F, Wang J, et al. 2011. Study of an improved Modified C-band radar quantitative precipitation estimation algorithm. Torrential Rain Disaster, 30(4): 321-327. (in Chinese) |
张培昌, 杜秉玉, 戴铁丕. 2001. 雷达气象学. 北京: 气象出版社, 177-187. Zhang P C, Du B Y, Dai T P. 2001. Radar Meteorology. Beijing: China Meteorological Press, 177-187(in Chinese)
|
Balakrishnan N, Zrnić D S, Goldhirsh J, et al. 1989. Comparison of simulated rain rates from disdrometer data employing polarimetric radar algorithms. J Atmos Ocean Technol, 6(3): 476-486. DOI:10.1175/1520-0426(1989)006<0476:COSRRF>2.0.CO;2 |
Fulton R A, Breidenbach J P, Seo D J, et al. 1998. The WSR-88D rainfall algorithm. Wea Forecasting, 13(2): 377-395. DOI:10.1175/1520-0434(1998)013<0377:TWRA>2.0.CO;2 |
Kessinger C, Ellis S, Van Andel J, et al. 2003. The AP clutter mitigation scheme for the WSR-88D∥Preprints, 31st Conf. on Radar Meteorology. Seattle Washington: Amer Meteor Soc, 526-529
|
Langston C, Zhang J, Howard K. 2007. Four-dimensional dynamic radar mosaic. J Atmos Ocean Technol, 24(5): 776-790. DOI:10.1175/JTECH2001.1 |
Lee G W, Zawadzki I. 2006. Radar calibration by gage, disdrometer, and polarimetry: Theoretical limit caused by the variability of drop size distribution and application to fast scanning operational radar data. J Hydrol, 328(1-2): 83-97. DOI:10.1016/j.jhydrol.2005.11.046 |
Mishchenko M I, Travis L D, Mackowski D W. 1996. T-matrix computations of light scattering by nonspherical particles: A review. J Quant Spectrosc Radiat Transf, 55(5): 535-575. DOI:10.1016/0022-4073(96)00002-7 |
Mishchenko M I. 2000. Calculation of the amplitude matrix for a nonspherical particle in a fixed orientation. Appl Opt, 39(6): 1026-1031. DOI:10.1364/AO.39.001026 |
Morrissey M L, Maliekal J A, Greene J S, et al. 1995. The uncertainty of simple spatial averages using rain gauge networks. Water Resour Res, 31(8): 2011-2017. DOI:10.1029/95WR01232 |
Mueller E A. 1977. Statistics of high radar reflectivity gradients. J Appl Meteor, 16(5): 511-513. DOI:10.1175/1520-0450(1977)016<0511:SOHRRG>2.0.CO;2 |
Qi Y C, Zhang J, Cao Q, et al. 2013a. Correction of radar QPE errors for nonuniform VPRs in mesoscale convective systems using TRMM observations. J Hydrometeorol, 14(5): 1672-1682. DOI:10.1175/JHM-D-12-0165.1 |
Qi Y C, Zhang J, Zhang P F. 2013b. A real-time automated convective and stratiform precipitation segregation algorithm in native radar coordinates. Quart J Roy Meteor Soc, 139(677): 2233-2240. DOI:10.1002/qj.2095 |
Qi Y C, Zhang J. 2017. A physically based two-dimensional seamless reflectivity mosaic for radar QPE in the MRMS system. J Hydrometeorol, 18(5): 1327-1340. DOI:10.1175/JHM-D-16-0197.1 |
Rosenfeld D, Ulbrich C W. 2003. Cloud microphysical properties, processes, and rainfall estimation opportunities∥Wakimoto R W, Srivastava R. Radar and Atmospheric Science: A Collection of Essays in Honor of David Atlas. Boston, MA: American Meteorological Society, 237-258
|
Smith P L. 1990. Precipitation measurement and hydrology: Panel Report∥Atlas D. Radar in Meteorology. Boston, MA: American Meteorological Society, 607-618
|
Smith P L. 1998. On the minimum useful elevation angle for weather surveillance radar scans. J Atmos Ocean Technol, 15(3): 841-843. DOI:10.1175/1520-0426(1998)015<0841:OTMUEA>2.0.CO;2 |
Tabary P. 2007. The New French operational radar rainfall product. Part Ⅰ: Methodology. Wea Forecasting, 22(3): 393-408. DOI:10.1175/WAF1004.1 |
Villarini G, Krajewski W F. 2008. Empirically-based modeling of spatial sampling uncertainties associated with rainfall measurements by rain gauges. Adv Water Resour, 31(7): 1015-1023. DOI:10.1016/j.advwatres.2008.04.007 |
Zhang G, Vivekanandan J, Brandes E. 2001. A method for estimating rain rate and drop size distribution from polarimetric radar measurements. IEEE Trans Geosci Remote Sens, 39(4): 830-841. DOI:10.1109/36.917906 |
Zhang J, Howard K, Gourley J J. 2005. Constructing three-dimensional multiple-radar reflectivity mosaics: Examples of convective storms and stratiform rain echoes. J Atmos Ocean Technol, 22(1): 30-42. DOI:10.1175/JTECH-1689.1 |
Zhang J, Qi Y C. 2010. A real-time algorithm for the correction of brightband effects in radar-derived QPE. J Hydrometeorol, 11(5): 1157-1171. DOI:10.1175/2010JHM1201.1 |
Zhang J, Howard K, 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. DOI:10.1175/2011BAMS-D-11-00047.1 |
Zhang J, Howard K, Langston C, et al. 2016. Multi-radar multi-sensor(MRMS)quantitative precipitation estimation: Initial operating capabilities. Bull Amer Meteor Soc, 97(4): 621-638. DOI:10.1175/BAMS-D-14-00174.1 |
Zhang Z, Qi Y C, Li D H, et al. 2021. A real-time algorithm to identify convective precipitation adjacent to or within bright band in the radar scan domain. J Hydrometeorol, 22(5): 1139-1151. DOI:10.1175/JHM-D-20-0005.1 |
Zhu Z W, Qi Y C, Cao Q, et al. 2020a. Particle size distribution characteristics within different regions of mature squall-line based on the analysis of global precipitation measurement dual-frequency precipitation radar retrieval. IEEE Geosci Remote Sens Lett. DOI:10.1109/LGRS.2020.3019384 |
Zhu Z W, Qi Y C, Cao Q, et al. 2020b. Conversion of the vertical profile of reflectivity from Ku-band to C-band based on the drop size distribution measurements of the global precipitation measurement mission dual-frequency precipitation radar. IEEE Trans Geosci Remote Sens. DOI:10.1109/TGRS.2020.3025803 |