1. 引言
降水作为洪水预报的输入条件,直接决定了水文模型预报精度,也会直接影响洪涝灾害防治与水库兴利减灾等决策工作。随着数值天气预报技术的发展,降雨预报产品迅速兴起并广泛应用于洪水预报中 [1] [2] [3] [4]。目前国内外常见的数值预报产品有欧洲ECWMF、美国GFS、日本JMA、中国气象局NMCP等。各数值模式预报产品在不同流域、不同降水阈值、不同的预见期存在差异 [5] - [12]。因此,在选取预报产品用于洪水预报之前,需对其进行精度评估,判断模式在研究区的适用性是有必要的。
目前,国内外学者对降雨预报能力的评估做了大量研究,并取得了显著成果,例如:何超禄等 [13] 综合评估了TIGGE中的多模式数据在中国干旱半干旱区的预报效果,研究表明ECWMF模式总体表现上更稳定;潘留杰等 [14] 采用邻域法分析ECWMF、JMA、T639高分辨率模式的降水预报性能,研究表明增大邻域半径能提高技巧评分;张玉兰等 [15] 认为现有研究主要考虑总量误差,对降雨预报的时空分布研究较少,提出了一种综合总量和时空分布的降雨预报方法。
然而,大多数预报评估研究着眼于同时刻站点对站点、网格对网格和面雨量与面雨量的总量误差,研究尺度较大,对于引发中小流域洪涝的短时强降水研究较少。传统的检验评估方法忽略了降水气候概率在时间和空间上的差异,但这些时空差异会影响洪水预报的洪峰、峰现时间 [16] - [21]。因此,评估高分辨模式建立一套考虑时空误差的评估体系是有必要的。
因此,从总量误差、空间分布和时间分布等多角度出发建立一套评价体系,并采用上述方法对NMCP模式在诺敏河研究区的预报效果进行评价,以期为NMCP模式的应用和多模式集成提供一定的理论支持。
2. 研究区域及数据介绍
2.1. 研究区域
诺敏河发源于大兴安岭支脉伊勒呼里山南麓,在呼伦贝尔盟莫力达瓦达斡尔族自治旗尼尔基镇附近分两支注入嫩江,是嫩江右岸7条一级支流之一。诺敏河流域位于北纬48˚00'~50˚30',东经121˚45'~124˚35',属于温带大陆性气候,夏季炎热多雨,冬季寒冷干燥。流域内降雨量年内、年际分配不均匀,主要集中在7~8月,多年平均年降水量约480 mm。本文研究的古城子水文站距离河口42 km,断面以上河长398 km,控制流域面积25,230 km2,每年7~9月降雨历时长,强度高,极易造成洪涝灾害 [22]。古城子流域概况见图1。
![](//html.hanspub.org/file/3-2950940x8_hanspub.png?20140125001100336)
Figure 1. Overview of Guchengzi watershed
图1. 古城子流域概况图
2.2. 数据及预处理
2.2.1. 预报数据来源
本次拟评价的降水数值预报模式为来自于中央气象局的NMCP降水数据。NMCP是创新性运用机器学习和数值模式的双预报引擎模式,具有高分辨率、高精准度的特征,能促进国内天气预报业务以及数值预报技术的发展,在洪涝灾害防治与水库兴利减灾等有一定的应用前景,有效降低洪涝等自然灾害造成的损失。NMCP数据每日世界时(Universal Time Coordinated, UTC) 00:00和12:00发布,共计2次,对应北京时间8:00和20:00,预报时长为72 h,时间分辨率为1 h,空间分辨率为0.05˚ × 0.05˚。默认数据格式为flt类型,可用Visual Studio 2017(VS)中ReadBinSingle方法进行查看和转换。
2.2.2. 实测数据来源
实测资料来自中国气象科学数据信息网CMORPM (NOAA Climate Prediction Center Morphing Method)卫星与全国3 × 104个自动观测站的逐时降水量融合产品。时间序列选取与预报数据保持一致,空间分辨率为0.05˚ × 0.05˚。有研究表明 [17],CMORPH卫星与自动站降水融合资料在降水量值和空间分布上能够更好的反映实际情况,并取得较高的融合降水精度。不同降水量级下的融合产品质量均优于国际上可接受的误差标准。该融合产品能够准确地抓住典型区强降水过程,对于定量监测降水过程具有优势。因此,该数据集可以作为参照物来评估预报模式的预报性能。
2.2.3. 数据预处理
考虑数据序列的完整性和代表性,选取了2021年07月~09月的预报数据与实测数据。预报数据和实测数据均为网格数据,且时间尺度和空间尺度均一致。依据古城子流域边界利用VS软件进行剪裁并提取出研究区范围数据。实测数据与预报数据的流域面雨量均采用网格平均法,通过将流域内所有网格降雨值取平均求得。
3. 建立评价体系
3.1. 量级评价指标
常见的降雨量级评价指标有流域面雨量、降雨中心雨量、降雨标准差、累积最大1 h,3 h,6 h,12 h降雨量平均值、累积最大1 h,3 h,6 h,12 h降雨量最大值等 [23]。基于洪水预报中重点关注洪水总量以及洪峰峰量等,引入相对误差对流域雨量进行精度评价。为使评分与技巧成正相关,取如下指标:
(1)
式中,
为累积最大i小时降雨的总量相对误差,取值范围在0~1之间,数值越大表明预报效果越好;
为累积最大i小时预报降雨量;
为累积最大i小时实测降雨量。在洪水模拟通常会关注最大1,3,6,12 h特征时段和全时段累积降雨量,因此i取1,3,6,12和all。
最后,需要将全时段、各个特征时段评价指标综合,得到最终的量级评价指标为:
(2)
式中:
为各特征时段量级评价所占权重,可根据流域产汇流特性进行赋值,或按照算术平均取值。本文选用平均取值,故
。
3.2. 空间分布评价指标
常见的降水空间分布指标有不均匀度、标准差与平均值之比、均方根误差(RMSE)等 [23]。由于中小尺度流域降水存在着明显的空间分布不均匀性,高分辨率数值模式的降水预报存在着较大的预报偏差。TS评分、RMSE等传统检验指标对观测与预报空间分布一致性要求过于苛刻,导致高分辨率预报模式的评分低于粗网格模式,无法客观地反映高分辨率模式地优势。利用邻域法邻域窗可变的特点,引入FSS评分来评价模式预报效果的空间尺度变化,进而得到优化的降水预报结果。FSS主要评价预报与观测面积的偏差。具体步骤如下:
a) 按照雨量等级划分标准,将12小时雨量划分为7个等级(见表1),但考虑到流域内发生大暴雨和特大暴雨次数较少,故将其归入暴雨等级内一并讨论,且无雨对于洪水的形成几乎没有直接的影响,故不做讨论。因此,定义降水阈值Gk (k = 1, 2, 3, 4),依次对应0.1 mm,5 mm,15 mm和30 mm。
![](Images/Table_Tmp.jpg)
Table 1. Classification of areal rainfall
表1. 面雨量划分等级
b) 对模式格点进行“0/1”二分化处理,即超过降水阈值Gk的格点记为1,反之为0:
(3)
(4)
式中:
和
分别代表观测和预报数据二分化处理后的值;
和
分别代表观测和预报数据,i为经向格点,j为纬向格点。
c) 邻域窗选取。模式实际应用的核心问题在于选择合适的空间尺度使得预报模式表现最好。对任意邻域半径的空间窗来水,区域内超过阈值的格点与总格点比值为:
(5)
(6)
式中:i,j分别为区域x轴和y轴方向的格点,取值范围为1−Nx和1−Ny;k,l分别为邻域窗内x轴和y轴方向的格点数,n为邻域窗半径,当n = 1时,表示格点本身的预报技巧。邻域半径n的均方根误差FBS (Fraction Brier Score)为:
(7)
FBS本质上是概率预报,且是一个负向评分指数,结果越接近0越好。为此对FBS进行转换,得到FSS。
(8)
其中:FSS区间为[0, 1],FSS评分越高,预报能力越好。FSS为0时表示实测与预报无相关性。理想情况下在保证实测与观测降水事件频率一致时,FSS随着邻域半径的增加趋于1。
3.3. 时间分布评价指标
由于气象预报在预报时刻存在误差,为了准确评价实测和预报降雨时间分布的吻合程度,引入时间域半径对预报时刻进行模糊化,利用TFSS评分来作为时间分布评价指标。
a) 参考中小流域面雨量等级划分的相关标准,定义降水阈值Dk (k = 1, 2, 3, 4),依次对应0.1 mm,0.5 mm,1.5 mm和5 mm。
b) 对各时刻面雨量进行“0/1”二分化处理,即超过降水阈值Dk的格点记为1,反之为0
(9)
式中:
为t时刻的降雨量,
为t时刻转换后的值。
c) 选择时间域半径n。
(10)
式中:
为t时刻的邻域概率,取值区间为[0, 1]。
d) 计算TFSS。
(11)
式中:
为t时刻预报数据的邻域概率,
为t时刻观测数据的邻域概率。TFSS取值范围为0~1,数值越大表明预报效果越好。
4. 结果与分析
4.1. 时空尺度确定
为了更好的反映预报数据的时空分布,使用邻域法将预报降水和实测降水进行时空尺度放大,采用FSS指标评估空间分布,采用TFSS指标评估时间分布。
图2为不同降水阈值和空间尺度下NMCP模式的FSS评分。由图2可知,FSS评分随邻域半径的增大而增大,且在同一邻域半径内,随之降水阈值的增大,FSS评分降低。NMCP模式在小雨、中雨和大雨预报的FSS评分较高,表明预报产品对降雨的空间分布和大小预报准确性较高。邻域半径的增大对小雨的FSS评分影响较小,对暴雨表现为显著增长。如邻域半径为5 km时,小雨FSS评分为0.670,暴雨FSS评分仅0.110;当邻域半径增长至60 km时,小雨FSS评分为0.751,暴雨FSS评分增长至0.296。当邻域半径达到60 km时,各个等级降水的FSS增长均已成明显变缓的趋势。考虑到计算效率和FSS评分,在实际应用中较优的邻域半径为60 km。
![](//html.hanspub.org/file/3-2950940x34_hanspub.png?20140125001100336)
Figure 2. FSS score of NMCP model under different precipitation thresholds and spatial scales
图2. 不同降水阈值和空间尺度下NMCP模式的FSS评分
图3为不同邻域半径下不同降雨量级的逐12小时FSS的时间演变。可以看出,小雨FSS的变化范围在0.67~0.73之间,中雨FSS的变化范围在0.4~0.54之间,说明NMCP模式对中雨及以下的预报能力较好,FSS基本上随着预报时效的增加而下降,其中12 h的FSS评分最高。但随着降水量级的增加,预报能力逐渐减弱,NMCP的FSS评分呈现先升后降的趋势,暴雨的FSS评分基本在0.2以下,尤其在60 h以后表现为无任何技巧。对于不同量级降雨,FSS总是随着邻域半径的增加而增加,尤其对前48 h暴雨来说,增加邻域半径可以获得较为明显的FSS增长。
![](//html.hanspub.org/file/3-2950940x35_hanspub.png?20140125001100336)
Figure 3. Time evolution of 12 h FSS scores under different neighborhood radii
图3. 不同邻域半径下,逐12 h FSS评分的时间演变
图4为不同降水阈值和时间尺度下NMCP产品的TFSS评分。由图4可得,时间域半径越大,TFSS评分越好,且在同一时间域半径,TFSS随降雨阈值的增大而降低。暴雨TFSS评分出现0的现象,分析数据可知,研究时间内从面雨量角度出发发生暴雨的次数极少。以小时降雨量超过1.5 mm为例,对比时间域半径为1 h和2 h的邻域概率场,TFSS评分从0.345上升至0.441,上升0.096,但当时间域半径继续增加时,TFSS评分上升不足0.05。随着时间域半径增加,TFSS上升的趋势更加缓慢。在实际应用中,时间域半径定为2 h为较合理的选择。
![](//html.hanspub.org/file/3-2950940x36_hanspub.png?20140125001100336)
Figure 4. TFSS score of NMCP products under different precipitation thresholds and time scales
图4. 不同降水阈值和时间尺度下NMCP产品的TFSS评分
4.2. 场次降雨指标评价结果
选取诺敏河古城子流域2021年07月17日1时至2021年07月20日0时共72 h降雨过程。NMCP模式每12 h预报1次,按预见期0~12 h,12~24 h,24~36 h,36~48 h,48~60 h,60~72 h将NMCP模式分为6段,分别记作NMCP12,NMCP24,NMCP36,NMCP48,NMCP60,NMCP72。
“7.17”诺敏河极端降雨过程中的强降水落区主要集中在流域的右下角,如图5所示,古城子流域的降雨中心值约250 mm/3d。NMCP12和NMCP24预测的强降水范围大致与实况一致,但位置偏北偏西,部分区域的强降水分布以及强度上与实际存在一定的空间差异,主要表现为低估强降雨。NMCP48、NMCP60及NMCP72在估测强降水落区和强度上与实际存在明显差异,降水中心分散且严重低估。
![](//html.hanspub.org/file/3-2950940x37_hanspub.png?20140125001100336)
Figure 5. “7.17” Guchengzi watershed precipitation process cumulative precipitation distribution (mm/3d)
图5. “7.17”古城子流域降水过程累积降水量分布(mm/3d)
为了量化NMCP模式在“7.17”诺敏河极端暴雨过程中的降水误差,采用本文的评价指标对NMCP模式进行评价与分析。
图6为实测降雨与不同预见期预报降雨数据对比图。从图6中可以得出,NMCP12、NMCP24、NMCP36表现较好。在降水发生的前13 h内,预报数据均表现为明显高估实际的降水量。预报数据均提前捕捉到强降水的峰值。在强降水过程发生的57 h后,预报数据均表现低估实际的降水量。NMCP48、NMCP60、NMCP72表现较差,严重低估强降水的峰值。
![](//html.hanspub.org/file/3-2950940x38_hanspub.png?20140125001100336)
Figure 6. Change of hourly area rainfall of “7.17” Guchengzi station (mm/h)
图6. “7.17”古城子站逐小时面雨量的变化(mm/h)
从表2中可以得出,在全时段累积降水量中,NMCP12,NMCP24,NMCP36明显高估于实际降水量,NMCP60,NMCP72表现为明显的低估,NMCP48表现较好。在累积最大1 h,3 h和6 h降水量中,NMCP12、NMCP24、NMCP36与实际降水吻合度较好,
趋于1;NMCP48,NMCP60,NMCP72明显低估实际降水量,
基本小于0.5;在累积最大12 h降水量中,前36 h的预报数据表现为高估,后36 h表现为低估。综合分析量级评价指标,前36 h预报数据
维持在0.8以上,大致与实况一致,尤其对短时强降雨量预报精度较高。
![](Images/Table_Tmp.jpg)
Table 2. NMCP cumulative precipitation and relative error in characteristic period and whole period
表2. NMCP累积特征时段和全时段降水量以及相对误差
根据上文参数选取结果,设置60 km为邻域半径,2 h为时间域半径,计算出不同预报时段内的FSS评分和TFSS评分,结果如表3所示。由表3可知:
1) NMCP降雨预报精度在相同预见期下,评分总体上随着降雨量级的增大而降低,这一结果符合天气统计学的规律。
2) 预报数据对于小雨和中雨的预报精度较高,对于小雨的FSS和TFSS评分基本维持在0.9以上,对于中雨的FSS和TFSS评分在0.8左右,且预报精度随预报时效的变化波动幅度较小。对于大雨,FSS和TFSS均呈现出先升后降的趋势,在预见期48 h时达到最高,48 h后评分迅速下降,72 h的TFSS评分仅0.154。对于暴雨,FSS仍呈现先升后降的趋势,在48 h达到最高,为0.540,在48 h后的FSS均为0;而TFSS仅在24 h大于0,其他时段均表现为预报和观测完全分离。
![](Images/Table_Tmp.jpg)
Table 3. FSS and TFSS scores of NMCP in different forecast periods and rainfall levels
表3. NMCP在不同预见期、不同降雨等级的FSS与TFSS评分
5. 结论
高分辨率数值预报存在一定的时空不确定性,传统的降雨预报指标如预报偏差、误报率、CSI评分等由于无法反映模式预报的空间结构与尺度变化,导致低估高分辨率模式的预报精度,无法应用于洪水预报等场景中。因此,降水量级、空间分布和时间分布多角度出发建立降雨预报评估方法,选取诺敏河古城子流域作为研究对象,采用上述评价体系对NMCP模式进行了评价与分析:在量级上,综合考虑洪水特征时段和总降雨过程的面雨量,引入相对误差作为评价指标;在空间分布上,采用邻域法对NMCP模式降水概率计算时,综合考虑计算效率以及FSS,选用60 km作为最优邻域半径;在时间分布上,引入时间域半径考虑预报的时间误差,对比TFSS与时间域半径热力图,选用2 h为最优时间域半径。
选取辽宁诺敏河流域“7.17”典型降水事件进行分析:在降水空间分布上,预见期在36 h内的预报数据能较好地反映降水空间分布,但仍表现出大雨低估,小雨高估。在量级上,前36 h预报数据
维持在0.8以上,大致与实况一致,尤其对短时强降雨量预报精度较高。在时空分布上,NMCP对中雨及以下的降水事件预报效果较好,FSS和TFSS评分均在0.8以上;对于大雨,FSS和TFSS随预见期的延长均表现为先升后降的趋势,在48 h达到最高;NMCP对暴雨的检验效果较差,这与暴雨样本次数少有着一定的关系。
综上可见,采用该评价方法能考虑预报数据的时空误差从而能得到更真实地评价指标,评价的结果可以为今后降雨数值预报的校正和应用以及多源预报模式的集成提供参考。由于本文只评价了NMCP在半干旱小流域的适用性,对湿润区流域以及大流域未进行研究,今后可进行对比分析其适用性。另外,本文主要针对于防汛需求,该评价体系对于其他需求的适用性在今后需进一步分析。
基金项目
国家研发计划重点项目山区暴雨山洪水沙灾害预报预警关键技术研究与示范(2019YFC1510700)。
NOTES
*通讯作者。