快速检索
  气象   2024, Vol. 50 Issue (11): 1397-1408.  DOI: 10.7519/j.issn.1000-0526.2024.050601

技术交流

引用本文 [复制中英文]

李妮娜, 刘凑华, 林建, 等, 2024. 2019—2021年北京春季温度预报的精细化检验评估[J]. 气象, 50(11): 1397-1408. DOI: 10.7519/j.issn.1000-0526.2024.050601.
[复制中文]
LI Nina, LIU Couhua, LIN Jian, et al, 2024. Refined Evaluation of Spring Temperature Forecast in Beijing During 2019-2021[J]. Meteorological Monthly, 50(11): 1397-1408. DOI: 10.7519/j.issn.1000-0526.2024.050601.
[复制英文]

资助项目

国家重点研发计划(2021YFC3000904)资助

第一作者

李妮娜, 主要从事精细化预报检验评估工作.E-mail: linn@cma.gov.cn

文章历史

2023年6月25日收稿
2024年5月6日收修定稿
2019—2021年北京春季温度预报的精细化检验评估
李妮娜 , 刘凑华 , 林建 , 代刊 , 韦青 , 赵声蓉     
国家气象中心,北京 100081
摘要:利用2019—2021年2—4月北京地区国家级气象台站逐3 h观测资料,基于统计分析方法对欧洲中期天气预报中心全球数值模式(ECMWF)、中国气象局全球区域同化预报系统(CMA-GFS)以及智能网格预报国家级指导报(SCMOC)、国省融合预报(SMERGE)的温度预报开展了精细化检验评估。结果表明,ECMWF和CMA-GFS对北京春季的温度预报多出现负偏差,且山区和平原地区无明显差异,但在夜间更加突出。SCMOC和SMERGE对温度预报有较好订正能力,温度预报平均误差集中在-1~1℃,温度预报准确率较高且平均绝对误差也较小。四种产品对于24 h变温以及日内昼夜温差的预报存在问题为:所有产品对强变温的变化幅度预报不足,且智能网格预报产品也未体现明显的订正能力;此外,所有产品对昼夜温差幅度的预报存在1~3℃的正偏差,SCMOC有一定订正能力,而SMERGE高估昼夜温差的情况更突出。模式预报昼夜温差偏大与低温(05时温度)预报偏低密切相关,而智能网格预报的高温(14时)预报偏高也不容忽视。通过精细化检验分析,说明智能网格预报在关注整体准确率提升(平均绝对误差,减小)的同时,还应关注天气过程发展演变特征。
关键词精细化检验    偏差特征    日变化    过程演变    温度预报    
Refined Evaluation of Spring Temperature Forecast in Beijing During 2019-2021
LI Nina, LIU Couhua, LIN Jian, DAI Kan, WEI Qing, ZHAO Shengrong    
National Meteorological Centre, Beijing 100081
Abstract: Using the 3 h observation data from national-level stations in Beijing Region in spring (February to April) during 2019-2021, we evaluate the temperature forecasts made by numerical models including European Centre for Medium-Range Weather Forecasting (ECMWF) and Global and Regional Assimilation and Prediction System (CMA-GFS), and by national gridded guidance forecast product (SCMOC) as well as the revised feedback gridded forecast product at provincial level (SMERGE). The results show that the spring temperature forecasts in Beijing by the models of ECMWF and CMA-GFS often show negative biases. The biases have no significant difference in mountainous and plain areas, but more prominent during nighttime periods. The gridded forecast products (SCMOC and SMERGE) have a good ability to correct the temperature forecasted by the models (ECMWF and CMA-GFS). The temperature forecast biases of the gridded products are concentrated in the range of -1-1℃, meanwhile the forecast accuracy is higher and the mean absolute error is lower than that of model forecast. Some problems are found in the forecasts of 24 h temperature change and diurnal temperature range from the four products. The amplitude of intense 24 h temperature change forecasted by all products is relatively smaller than that of observation, and the gridded forecast products fail to demonstrate significant correction ability. In addition, the diurnal temperature range forecasted by all products has a positive deviation of 1-3℃ compared to the observation. SCMOC has a better correction ability for the diurnal temperature range by the models, while SMERGE overestimates the difference more prominently than that by models. The positive deviation of diurnal temperature range forecast is closely related to the underestimation of low temperature at 05:00 BT in the models, while the overestimation of high temperature at 14:00 BT can not be ignored in the grid forecasts. The refined analysis suggests that gridded forecast products should not only focus on improving overall accuracy (reducing biases), but also on the development and evolution of synoptic processes.
Key words: refined evaluation    bias characteristic    diurnal cycle    process evolution    temperaure forecast    
引言

华北地区位于我国第二、三级地形阶梯,西倚太行山,北靠燕山山脉,东临渤海。多变的地形及特殊的地理环境使得该地区的天气状况复杂。特别地,在全球变化背景下,极端降水、高温热浪、低温寒潮等灾害性天气频发,这类恶劣天气过程严重影响了社会生产和人民生活(Zhai et al, 2005; Zhang and Zhai, 2011; Liao et al, 2021; Chen et al, 2021)。对于发展程度高,人口密集的北京更是如此。气象灾害的防御亟需提高天气预报的准确性和及时性。为应对人们对气象服务更高的需求,在天气预报业务的现代化进程中,无缝隙精细化智能网格天气预报业务得以迅速发展(赵声蓉等, 2012; Hamill et al, 2015; 毕宝贵等, 2016; 代刊等, 2016; 金荣花等, 2019),即便如此,智能网格天气预报也很难满足极端性天气的预报需求。为了进一步提高极端灾害性天气的预报准确率和防灾减灾能力,通过精细化的检验评估来促进数值模式、智能网格预报技术的发展势在必行(毛冬艳等, 2014; 陈昊明等, 2021)。

早期的检验评估多聚焦于单模式或单变量(王雨, 2003; 管成功等, 2006; 2008),且多关注逐日或逐12 h时间尺度的预报(韩洁等, 2015; 季晓东和漆梁波, 2018)。随着精细化预报产品的增加,为应对传统检验方法所带来的问题,出现了一系列空间检验方法(Davis et al, 2009; Gilleland et al, 2008; 刘凑华和牛若芸, 2013; 潘留杰等, 2014; 符娇兰和代刊, 2016; 赵滨和张博, 2018; 陈法敬等, 2019)。但是这些研究多注重空间上的对比,对于天气过程在发生、持续和结束等特征的关注较少(金荣花等,2019)。因此,越来越多的研究开展了各产品间逐3 h或逐1 h的对比检验,以便了解数值模式或客观产品对日内尺度的精细化特征的预报能力(卢冰等, 2017; 许晨璐等, 2017; 陈龙等, 2024)。韦青等(2020)通过对国家级智能网格指导报、省市级订正反馈网格预报以及模式输出的降水和温度预报进行全面检验,发现智能网格预报可以改善午后降水峰值预报显著提前的偏差,且最高/最低温预报较数值模式改进明显,并指出智能网格预报需进一步改进强降水预报空报率较高、最低温预报模型构建等问题。Yuan et al(2020)指出ECMWF预报的西南地区降水事件开始以及达到峰值的时间明显偏早且午后降水峰值被高估,这些偏差在西南地区西部更加突出。值得注意的是,这些精细化的检验多针对降水,针对温度的检验相对较少。

北京春季受西风槽频繁活动影响,地面冷锋活跃,天气复杂多变。气温变化快,昼夜温差大且大风天气频发,加之受北京周边复杂地形影响,给精细化预报带来了极大挑战。因此,本研究聚焦北京春季的温度预报,针对当前业务中常用的数值模式以及智能网格预报产品,开展精细化对比检验,分析不同产品的预报偏差特征、预报难点及业务参考价值。一方面为客观产品在业务预报中的合理应用提供参考,另一方面为智能网格预报产品后期研发与改进提供可靠依据。

1 检验数据和方法

考虑到近年来北京的入春时间整体呈现偏早趋势,本文选取2—4月为研究时段进行北京春季温度预报检验分析。实况数据为2019—2021年2—4月北京地区20个国家级气象台站的温度及风场观测资料。预报数据分为两类:两套数值模式[欧洲中期天气预报中心全球数值模式(ECMWF)、中国气象局全球区域同化预报系统(CMA-GFS)]和两套智能网格预报[国家级指导报(SCMOC)以及国省融合预报(SMERGE)]的2019—2021年2—4月逐日08时(全文均为北京时)和20时起报的逐3 h(预报时效, 108 h后为逐6 h)温度预报产品。其中,ECMWF预报产品的空间分辨率为12.5 km,CMA-GFS于2016年6月正式进入业务运行并向全国下发相关产品,其空间分辨率为25 km(沈学顺等,2017)。两套智能网格预报产品空间分辨率均为5 km(金荣花等,2019韦青等,2020),其中,SCMOC是基于CMA-GFS采用GMOS方法的精细化订正预报产品,由国家气象中心制作下发;SMERGE为省级业务单位在SCMOC基础上修改订正后上传得到的。将对这两类(四套)产品进行对比检验,以期深入认识不同产品对北京春季温度预报的偏差特征。

采用双线性插值方法将逐3 h温度预报产品插值到北京地区20个国家级气象台站后开展检验分析。检验内容包括平均误差、平均绝对误差、预报准确率、分类误差频率统计以及日内尺度的变化特征等。

其中,温度预报准确率:参照中国气象局中短期天气预报质量检验相关办法,定义为观测与预报温度之差的绝对值≤2℃的占比,即:$R_c=\frac{N_{r 2}}{N} $;其中,N为总站次,Nr2为温度绝对误差≤2℃的站次。

24 h变温:利用逐3 h温度数据,计算每个t时刻及t-24时刻的温度差,获得逐3 h滑动的24 h变温,即:ΔT24=Tt-Tt-24,其中t为研究时段内逐3 h时刻。定义ΔT24≥6℃为强升温时刻,ΔT24≤-6℃为强降温时刻。检验针对ΔT24开展。

昼夜温差:利用14时和05时的温差表征昼夜温差,以检验不同产品对温度日内变化幅度的预报能力。

2 检验结果分析 2.1 温度预报基本性能检验

从2019—2021年2—4月全球模式ECMWF、CMA-GFS以及智能网格预报产品SCMOC、SMERGE对北京地区的温度预报来看(图 1),大部分产品的温度预报准确率(绝对误差≤2℃)随预报时效增加而逐渐降低。CMA-GFS温度预报准确率随预报时效变化较小,集中在0.30~0.45,较其他预报产品的准确率偏低。ECMWF预报准确率明显高于CMA-GFS,8 d以上的预报准确率高于0.35。SCMOC与ECMWF预报准确率相当,在短期预报时效(<72 h)优于ECMWF。SMERGE订正作用明显,在绝大部分预报时效下准确率最高。平均绝对误差随着预报时效增加而逐渐增大。CMA-GFS在各预报时效下绝对误差均最大,大部分预报时效下误差均超过3℃,10 d内平均绝对误差达3.40℃。除CMA-GFS外,其余产品5 d内的平均绝对误差小于3℃,5 d以上平均绝对误差不超过4℃;SMERGE误差最小,10 d内平均绝对误差为2.45℃。对于平均误差,ECMWF和CMA-GFS在不同时效下均呈现温度预报偏低的特征,分别偏低1℃和2℃左右;对于SCMOC和SMERGE这两种预报产品来说,大部分预报时效下温度预报平均误差均略低于0℃。

图 1 2019—2021年2—4月北京地区不同产品温度预报(a)准确率(绝对误差≤2℃),(b)平均绝对误差和(c)平均误差随预报时效的变化 Fig. 1 Variation of (a) accuracy (absolute error ≤2℃), (b) mean absolute error and (c) mean error of temperature by different forecast products in Beijing with valid time from February to April during 2019-2021

图 2a给出北京春季温度预报准确率的空间分布,整体来看,各预报产品均呈现出低海拔地区的预报准确率略高于高海拔地区的特征。从平均误差的空间分布来看(图 2b),四种预报产品对春季北京地区温度预报整体偏低。CMA-GFS偏低尤为明显,中部地区/城区大多在3℃以上,ECMWF温度预报平均误差集中在-2~-1℃。SCMOC以及SMERGE对温度预报有明显订正能力,平均误差集中在-1~1℃。

图 2 2019—2021年2—4月北京地区不同产品温度预报(a)准确率和(b)平均误差的空间分布 注:灰色阴影为地形高度。 Fig. 2 Spatial distribution of (a) accuracy and (b) mean error of temperature by different forecast products in Beijing from February to April during 2019-2021

为了解北京山区(站点海拔≥200 m)和平原地区(站点海拔<200 m)的昼夜温度预报特征差异,图 3a, 3c分别给出两类地形区不同时段温度预报的平均误差。ECMWF和CMA-GFS预报的温度存在系统性偏差,无论是山区还是平原地区,两个模式的温度预报平均误差都呈现为负偏差,这一特征在白天(08—17时)和夜间(20—05时)均存在。值得注意的是,夜间时段温度预报的平均误差幅度大于白天时段。智能网格预报SCMOC和SMERGE平均误差幅度明显低于模式预报:夜间时段平均误差由模式的-3~-2℃减小至-0.5℃左右;白天时段的预报则表现为较小的正偏差。分析不同大小的误差出现频率(图 3b),发现ECMWF、CMA-GFS温度预报误差为负的出现频率明显偏高,分别为67.35%和73.29%。SCMOC和SMERGE的温度预报误差较小,-1~1℃频率最高,但负偏差出现频率仍略高于正偏差。考虑到温度预报以负偏差为主,图 3d进一步分析了各预报产品负偏差出现频率的日变化特征。ECMWF与CMA-GFS,预报负偏差在23时出现频率最高。SCMOC以及SMERGE负偏差出现频率的峰值时刻较模式预报偏晚,分别出现在02时和05时。

图 3 2019—2021年2—4月北京(a)山区和(c)平原地区不同产品(a, c)不同时段温度预报平均误差, (b)预报误差频率统计和(d)预报负偏差出现频率日变化 注:图c为分布在各温度区间的频率;如,横坐标零值两侧的柱状,分别表示预报误差在[-1, 0)、[0, 1)的频率。 Fig. 3 (a, c) Mean error of temperature in different time periods in (a) mountains and (c) plain areas, (b) frequency statistics of different error bins, (d) diurnal variation of occurrence frequency of negative temperature forecast deviation by different forecast products in Beijing from February to April during 2019-2021

总体来说,CMA-GFS对北京春季温度预报效果较差,但是SCMOC产品预报效果明显提升,与ECMWF效果相当,甚至在短期时效优于ECMWF,SMERGE产品对北京春季温度预报效果最优。这四种预报产品均倾向于低估北京地区的春季温度,以CMA-GFS模式预报的温度偏低最为明显。北京的山区和平原地区均存在温度预报的负偏差,这种偏差在夜间时段更大,且夜间时段出现负偏差的频率更高。

2.2 24 h变温

考虑到北京春季经常受西风槽频繁活动影响,转折性天气过程较多,气温升降幅度大,24 h变温是温度预报关注的重点问题之一,这部分将聚焦24 h变温开展检验评估。

基于逐3 h滑动的24 h变温进行统计,图 4给出观测与预报中不同变温幅度的出现频率以及北京不同地区强升温及强降温的频率。从图 4a可以看到,观测中北京春季出现升温的频率(58.21%)高于降温频率(41.79%),其中以升温0~2℃的频率最高,可达22.49%。但关注强变温过程可以发现,北京地区强降温频率(7.88%)略高于强升温频率(6.21%),且山区发生强变温的频率高于平原地区(图 4b, 4c)。四种产品对变温幅度<2℃的过程预报频率偏高,而对强变温过程的预报频率明显不足。这种特征在SMERGE中表现尤为明显,强变温预报频率甚至不足观测的1/2。尽管从预报准确率、整体平均误差来看SMERGE温度预报较好,但SMERGE并没有很好地把握升降温过程。

图 4 2019—2021年2—4月北京地区观测和各预报产品(a)不同变温幅度的频率,及(b,c)不同区域强变温的发生频率(b)强升温(ΔT24≥6℃),(c)强降温(ΔT24≤-6℃) Fig. 4 (a) Frequency statistics in different 24 h temperature change ranges, and (b, c) the occurrence frequency of (b) intense temperature rising (ΔT24≥ 6℃) and (c) intense temperature cooling (ΔT24≤-6℃) by observation and different forecast products in Beijing from February to April during 2019 to 2021

聚焦强变温过程,表 1给出发生强变温的个例情况,正如前文所提到的,强降温的发生频次高于强升温频次。基于挑选的2132/2642个强升温/强降温站次,分析强变温发生前及发生时的预报特征。图 5给出强变温过程发生前一天及发生时不同温度预报误差的频率分布。对于强升温过程(图 5a5c),升温前ECMWF与CMA-GFS温度预报的正、负偏差基本相当,SCMOC和SMERGE温度预报偏高的频率略高于温度预报偏低的频率(图 5a)。而升温日,温度预报出现了明显的负偏差。四种产品出现预报负偏差的频率均在80%以上,CMA-GFS中预报误差超过-7℃的频率高达35%,其余产品温度预报误差在该区间的频率也超过了15%(图 5c)。对于强降温过程(图 5b5d),降温发生前温度预报多呈现负偏差(图 5b),四种产品预报出现负偏差的频率在78.56%~92.89%。降温发生时,ECMWF和CMA-GFS温度预报正、负偏差出现的频率较为接近,SCMOC及SMERGE倾向于高估温度;四种产品中温度预报偏高超过7℃的频率均在5%以上。

表 1 强变温主要发生日期、发生日数、时次数以及总站次数 Table 1 Main occurrence date, number of day, frequency of intense temperature change, as well as the total number of times and stations

图 5 2019—2021年2—4月北京地区不同产品在强变温(a,b)发生前及(c,d)发生日的温度预报误差频率分布(a,c)强升温(ΔT24≥6℃),(b,d)强降温(ΔT24≤-6℃) 注:图中柱状为分布在各温度区间的频率;如,横坐标零值两侧的柱状,分别表示预报误差在[-1, 0)、[0, 1)的频率。 Fig. 5 Frequency statistics at different temperature forecast error bins by different forecast products in Beijing from February to April during 2019-2021 for (a, b) one-day before and (c, d) on the day of occurrence of intense temperature change (a, c) intense temperature rising (ΔT24≥6℃), (b, d) intense temperature cooling (ΔT24≤-6℃)

综上所述,当出现强变温过程时,四种产品的变温幅度预报不足。对于强升温过程来说,主要是升温日温度预报偏低造成升温幅度预报不足。对于强降温过程来说,降温前温度预报偏低、降温日温度预报偏高共同影响了降温幅度的预报。

2.3 昼夜温差

北京春季温度另一突出特征是昼夜温差较大。本文利用14时(高温)和05时(低温)的温度差考察不同产品对逐3 h定时温度日内变化幅度的预报情况。从预报平均误差图中可以看出四种预报产品均高估了北京春季的昼夜温差(图 6)。ECMWF和CMA-GFS中,大部分站点的昼夜温差预报偏高2℃左右,CMA-GFS预报平均误差略小于ECMWF。SCMOC对于昼夜温差的预报有明显订正能力,其预报平均误差减小至1℃以内,仅为0.70℃。但对于SMERGE来说,北京地区昼夜温差预报平均误差甚至高于模式预报结果,最大昼夜温差的平均误差可达3.85℃,最小也超过2℃。区分08时和20时起报的预报发现,在两种不同起报时间下昼夜温差的预报平均误差基本一致(图略),表明昼夜温差幅度预报偏大属于系统性误差。

图 6 2019—2021年2—4月北京地区不同产品昼夜温差预报的平均误差空间分布 注:灰色阴影为地形高度。 Fig. 6 Spatial distribution of mean error of diurnal temperature range by different forecast products in Beijing from February to April during 2019-2021

不同昼夜温差出现频率的统计结果显示,北京地区春季昼夜温差集中在6~14℃(出现频率为48.83%),8~10℃出现频率最高(达12.76%)(图 7)。对于昼夜温差小于8℃的情况,四种预报均明显低估其出现频次,SMERGE对昼夜温差在6℃以内的频率预报甚至不足观测的1/2。对于昼夜温差较大的情况,四种预报产品均会高估其频率。另外,预报中最频繁出现的昼夜温差区间也较观测偏大:ECMWF(CMA-GFS)预报的昼夜温差更多集中在12~14℃(10~12℃),SMERGE预报的昼夜温差主要集中在12~16℃,仅SCMOC能够做出较为合理的预报,其预报的昼夜温差主要集中在8~12℃,最接近观测结果。

图 7 2019—2021年2—4月北京地区观测和各预报产品不同昼夜温差的出现频率 Fig. 7 Frequency statistics of different diurnal temperature ranges by observation and different forecast products in Beijing from February to April during 2019-2021

进一步分析昼夜温差预报偏大的情况,挑选各产品中预报误差超过4℃的日期考察该日14时和05时温度预报平均误差。如图 8所示,昼夜温差幅度预报误差较大时,模式产品的05时温度预报误差幅度明显大于14时温度预报误差,低温预报误差是高温预报误差幅度的3倍甚至更高;对于SCMOC和SMERGE来说,高低温预报误差幅度基本相当。综合高低温的预报特征,将昼夜温差预报呈现正偏差的情况分为5种类型,分别是:(1)低温预报较好(误差在±2℃),但高温预报偏高(误差>2℃);(2)高温预报较好,低温预报偏低(误差<-2℃);(3)高温预报偏高且低温预报偏低;(4)高低温预报均偏低(误差均<-2℃), 但低温预报偏低幅度更大;(5)高低温预报均偏高(误差均>2℃), 但高温预报偏高幅度更大。

图 8 2019—2021年2—4月北京地区不同产品昼夜温差幅度预报偏大日的(a)高低温的预报平均误差, (b)5种情况概率分布 注:图b柱状图上方数值为昼夜温差幅度预报偏高4℃以上的总站次数(单位:站·h)。 Fig. 8 (a) Mean error of high and low temperatures in days and (b) frequency statistics of five cases with large positive biases of diurnal temperature range by different forecast products in Beijing from February to April during 2019 to 2021

图 8b给出各类型预报误差的出现频率以及不同产品昼夜温差预报过大的站次。从图中可以看出四种产品昼夜温差预报偏大的站次存在一定差异。其中,CMA-GFS昼夜温差预报偏大的站次数高于ECMWF预报。SCMOC对昼夜温差的预报略有改善,正偏差较大的站次数较模式预报结果有所减少。与图 6中较高的昼夜温差预报平均误差对应,SMERGE产品中昼夜温差偏大的站次数也明显高于其余三种产品。关注不同类型偏差出现情况可以发现,模式预报与智能网格预报产品昼夜温差偏大的情况有所区别。对于ECMWF和CMA-GFS而言,昼夜温差幅度预报偏大很大程度上是由于低温预报偏低(绿色柱)造成的,这种情况出现频率分别为53.02%和51.09%。第二类主要的偏差为高温预报偏高与低温预报偏低,分别为31.86%和27.89%。值得注意的是,相较于智能网格预报产品,模式产品高低温预报均偏低, 但低温偏低幅度更大(紫色柱), 这种情况的出现频率也较高。整体来说,低温预报偏低是模式产品昼夜温差幅度预报偏大的最主要原因。

对于SCMOC和SMERGE来说,低温预报偏低、高低温预报均偏低的情况则明显减少。昼夜温差幅度预报偏大,在很大程度上是由于高温预报偏高和低温预报偏低(红色柱)共同导致的。此外,仅高温预报偏高(橙色柱)的频率明显增加;对应于图 8a,智能网格预报产品中低温预报偏低幅度略小于模式产品,但高温预报正偏差明显增大。结合图 2图 3分析可以推断,智能网格预报产品在订正温度时倾向于系统性增加温度来减小温度预报误差,但聚焦昼夜温差的预报时,就会过多地出现高温预报明显偏高的情况。后续需要结合模式预报偏差的主要问题,针对性做出温度预报的订正,使其在整体温度预报准确率提升的同时,也能对温度的日变化演变特征做出更精准预报。

2.4 预报误差的可能原因

整体来看,北京春季温度预报难点为:(1)变温过程中预报变温幅度偏小; (2)温度的日内变化(昼夜温差)幅度预报偏大。针对上述两方面的预报难点,进一步探讨其可能原因。

通过分析研究时段内北京地区所有出现强变温日的风场预报情况(图 9)可以发现,对于强降温过程来说,无论是ECMWF还是CMA-GFS,其预报的西北—偏西方向的风速仅为2~3 m·s-1,明显小于观测风速(3~4 m·s-1)。北京地区通常是在以西北气流为主的冷空气影响下发生降温,而模式预报的西北—偏西风风速偏弱,说明北京地区的降温幅度预报不足可能与冷空气的平流输送不足有关。对于强升温过程,模式预报的风场在各方向的风速、出现频率基本相当,与观测的风场特征较为相近,说明北京地区的升温幅度预报不足可能更多与辐射过程相关。

图 9 2019—2021年2—4月北京地区观测和模式预报的(a)发生强降温日(ΔT24≤-6℃)以及(b)发生强升温日(ΔT24≥6℃)风矢量分布 注:红色、蓝色分别表示观测、预报风速的一致性。 Fig. 9 Observed and model-forecasted wind vector distribution in the occurrence days of (a) intense temperature cooling (ΔT24≤-6℃) and (b) intense temperature rising (ΔT24≥6℃) in Beijing from February to April during 2019-2021

考虑到云辐射强迫对气温有重要影响,针对昼夜温差幅度预报偏大这一问题,进一步分析研究时段内所有误差偏大日,14时和05时的云量预报平均误差。从图 10可以看出,除ECMWF预报的白天云量较观测略偏高外,模式预报的云量通常较观测偏少,平均误差约-10%,甚至超过-20%。白天的少云会带来更强的辐射增温,而夜间云量的偏少会使向外长波辐射增强,日内的降温幅度更加明显。同时,图 10也显示了夜间的云量偏少较白天更加突出,这与前文提到的模式产品低温预报偏低是最主要误差原因相对应。也就是说,昼夜温差幅度预报偏高与当日云量的预报误差具有一定关联。以2021年3月1—5日过程为例,通过检验北京站当时总云量的预报可以发现,观测总云量可以达到90%左右,而对应的模式预报结果显示基本无云(总云量约为0%)(图略),造成这段时间白天(夜间)温度预报偏高(偏低)。正如秦庆昌等(2022)指出,此次春季温度预报失误的原因主要与观测中持续性低云天气过程有关,而模式无法再现这类低云特征。

图 10 2019—2021年2—4月北京地区昼夜温差幅度预报偏高日14时和05时云量预报平均误差 Fig. 10 Mean error of cloud cover forecast at 14:00 BT and 05:00 BT in days with large positive biases of diurnal temperature range forecast in Beijing from February to April during 2019-2021
3 结论和讨论

聚焦北京春季的温度,通过开展精细化检验评估,分析了业务常用数值模式ECMWF、CMA-GFS以及智能网格预报产品SCMOC和SMERGE的预报能力以及预报偏差特征,主要结论如下:

(1) CMA-GFS对北京地区春季温度预报准确率低,平均绝对误差大;但基于CMA-GFS采用GMOS方法的SCMOC产品预报效果明显提升,与ECMWF效果相当;SMERGE产品对北京春季温度预报准确率最高,平均绝对误差最小。所有预报产品均倾向于低估北京地区的春季温度。温度预报偏低这一问题在山区和平原地区无明显差异,但在夜间时段更加突出。SCMOC和SMERGE对温度预报的负偏差具有较好的订正能力,温度预报平均误差集中在-1~1℃。

(2) 北京春季转折性天气较多,24 h变温剧烈。当前四种产品对转折性天气预报能力较差,强变温过程中24 h变温幅度预报偏小;6℃及以上的强升温过程中幅度预报不足主要是受升温日温度预报偏低影响,-6℃及以下强降温过程幅度预报不足是降温前温度预报偏低和降温日温度预报偏高两者共同影响;强升温(强降温)幅度预报不足可能与辐射过程(平流输送)预报偏差相关。

(3) 北京春季日内昼夜温差较大,以8~10℃的温差最常出现。模式预报产品均过高估计了温度的日内变化幅度,昼夜温差预报偏高1~3℃。SCMOC对昼夜温差幅度高估的问题有所改善,但SMERGE对昼夜温差的预报甚至不如模式预报产品。对模式产品来说,低温预报偏低是昼夜温差幅度预报偏大最主要的原因,这可能与云量的预报偏差相关;智能网格预报产品中,高温预报偏高且低温预报偏低的出现频率有所增加,且仅高温预报偏高的情况也不容忽视。

本文通过细致分类,从统计检验和个例分析的角度全面给出了温度预报的偏差特征,并通过建立24 h变温、昼夜温差幅度这两个具有物理意义的评估指标进行检验,指出所有的预报产品对北京春季温度具有较好预报能力,但面对重要天气过程时,依旧存在预报误差大、预报难度高的情况。针对这些问题,所有预报产品应该在关注整体预报准确率提升的同时,重点聚焦天气过程演变特征,做出精准预报。另外,本文对偏差的成因进行了初步的探讨,发现强降温幅度预报偏低可能与模式预报的西北气流强度偏弱有关。但影响强降温的环流型众多,除了主导的西北气流型外,还可能存在偏北气流型、偏西气流型以及西南气流型等环流形势的影响(杜晓辉等, 2016杨晓玲等, 2016)。未来将拓展研究时段,对影响强变温的环流型进行分类,探究不同类型下的预报误差特征,明晰偏差出现时的主导环流系统以及关键影响因子。在此基础上开展针对性客观预报技术研究,完善当前技术体系中的预报因子和建模样本方案、偏差订正方案、预报融合集成方案等,减小针对转折性、高影响天气预报误差偏大的问题,为最终提升预报服务能力提供支撑。

参考文献
毕宝贵, 代刊, 王毅, 等, 2016. 定量降水预报技术进展[J]. 应用气象学报, 27(5): 534-549. Bi B G, Dai K, Wang Y, et al, 2016. Advances in techniques of quantitative precipitation forecast[J]. J Appl Meteor Sci, 27(5): 534-549 (in Chinese).
陈法敬, 陈静, 韦青, 等, 2019. 一种基于可预报性的暴雨预报评分新方法Ⅱ: 暴雨检验评分模型及评估试验[J]. 气象学报, 77(1): 28-42. Chen F J, Chen J, Wei Q, et al, 2019. A new verification method for heavy rainfall forecast based on predictability Ⅱ: verification method and test[J]. Acta Meteor Sin, 77(1): 28-42 (in Chinese).
陈昊明, 李普曦, 赵妍, 2021. 千米尺度模式降水的检验评估进展及展望[J]. 气象科技进展, 11(3): 155-164. Chen H M, Li P X, Zhao Y, 2021. A review and outlook of verification and evaluation of precipitation forecast at convection-permitting resolution[J]. Adv Meteor Sci Technol, 11(3): 155-164 (in Chinese).
陈龙, 陈静静, 胡媚, 等, 2024. 华南快速循环同化模式在湖南不同环流型下的小时降水预报性能检验[J]. 气象, 50(4): 434-448. Chen L, Chen J J, Hu M, et al, 2024. Evaluation of hourly quantitative preci-pitation forecast of the rapid updating cycle assimilation and forecasting system in south china based on circulation classification in Hunan Province[J]. Meteor Mon, 50(4): 434-448 (in Chinese). DOI:10.7519/j.issn.1000-0526.2023.061401
代刊, 曹勇, 钱奇峰, 等, 2016. 中短期数字化天气预报技术现状及趋势[J]. 气象, 42(12): 1445-1455. Dai K, Cao Y, Qian Q F, et al, 2016. Situation and tendency of operational technologies in short- and medium-range weather forecast[J]. Meteor Mon, 42(12): 1445-1455 (in Chinese). DOI:10.7519/j.issn.1000-0526.2016.12.002
杜晓辉, 代华光, 德庆, 2016. 西藏强降温时空分布特征与天气概念模型[J]. 西藏科技, (8): 62-64. Du X H, Dai H G, De Q, 2016. Spatial and temporal distribution characteristics and synoptic conceptual type of strong cooling in Tibet[J]. Tibet Sci Technol, (8): 62-64 (in Chinese).
符娇兰, 代刊, 2016. 基于CRA空间检验技术的西南地区东部强降水EC模式预报误差分析[J]. 气象, 42(12): 1456-1464. Fu J L, Dai K, 2016. The ECMWF model precipitation systematic error in the east of Southwest China based on the contiguous rain area method for spatial forecast verification[J]. Meteor Mon, 42(12): 1456-1464 (in Chinese). DOI:10.7519/j.issn.1000-0526.2016.12.003
管成功, 陈起英, 佟华, 等, 2008. T639L60全球中期预报系统预报试验和性能评估[J]. 气象, 34(6): 11-16. Guan C G, Chen Q Y, Tong H, et al, 2008. Experiments and evaluations of global medium range forecast system of T639L60[J]. Meteor Mon, 34(6): 11-16 (in Chinese). DOI:10.7519/j.issn.1000-0526.2008.06.002
管成功, 王克敏, 陈晓红, 2006. 2002—2005年T213数值降水预报产品分析检验[J]. 气象, 32(8): 70-76. Guan C G, Wang K M, Chen X H, 2006. Analyses and verifications of T213 rainfall forecast from 2002 to 2005[J]. Meteor Mon, 32(8): 70-76 (in Chinese). DOI:10.3969/j.issn.1000-0526.2006.08.012
韩洁, 庞翻, 陈卫东, 2015. WRF模式对宝鸡地区2013年夏季温度和降水预报的评估检验[J]. 陕西气象, (1): 1-6. Han J, Pang F, Chen W D, 2015. Evaluation and verification of summer temperature and precipitation forecast in Baoji, Shaanxi Province in 2013 by WRF model[J]. J Shaanxi Meteor, (1): 1-6 (in Chinese).
季晓东, 漆梁波, 2018. ECMWF模式降水预报与极端天气预报指数在暴雨预报中的评估与应用[J]. 暴雨灾害, 37(6): 566-573. Ji X D, Qi L B, 2018. Evaluation and application of ECMWF model precipitation and extreme weather forecast index of precipitation on heavy rainfall forecast[J]. Torr Rain Dis, 37(6): 566-573 (in Chinese).
金荣花, 代刊, 赵瑞霞, 等, 2019. 我国无缝隙精细化网格天气预报技术进展与挑战[J]. 气象, 45(4): 445-457. Jin R H, Dai K, Zhao R X, et al, 2019. Progress and challenge of seamless fine gridded weather forecasting technology in China[J]. Meteor Mon, 45(4): 445-457 (in Chinese). DOI:10.7519/j.issn.1000-0526.2019.04.001
刘凑华, 牛若芸, 2013. 基于目标的降水检验方法及应用[J]. 气象, 39(6): 681-690. Liu C H, Niu R Y, 2013. Object-based precipitation verification method and its application[J]. Meteor Mon, 39(6): 681-690 (in Chinese). DOI:10.7519/j.issn.1000-0526.2013.06.003
卢冰, 孙继松, 仲跻芹, 等, 2017. 区域数值预报系统在北京地区的降水日变化预报偏差特征及成因分析[J]. 气象学报, 75(2): 248-259. Lu B, Sun J S, Zhong J Q, et al, 2017. Analysis of characteristic bias in diurnal precipitation variation forecasts and possible reasons in a regional forecast system over Beijing Area[J]. Acta Meteor Sin, 75(2): 248-259 (in Chinese).
毛冬艳, 朱文剑, 樊利强, 等, 2014. GRAPES-MESO V3.3模式强天气预报性能的初步检验[J]. 气象, 40(12): 1429-1438. Mao D Y, Zhu W J, Fan L Q, et al, 2014. Preliminary synoptic verification on the performance of GRAPES-MESO V3.3 in severe weather forecasting[J]. Meteor Mon, 40(12): 1429-1438 (in Chinese). DOI:10.7519/j.issn.1000-0526.2014.12.001
潘留杰, 张宏芳, 王建鹏, 2014. 数值天气预报检验方法研究进展[J]. 地球科学进展, 29(3): 327-335. Pan L J, Zhang H F, Wang J P, 2014. Progress on verification methods of numerical weather prediction[J]. Adv Earth Sci, 29(3): 327-335 (in Chinese).
秦庆昌, 于波, 徐路扬, 等, 2022. 造成北京春季气温预报失误的持续性低云过程分析[J]. 气象, 48(6): 705-718. Qin Q C, Yu B, Xu L Y, et al, 2022. Analysis of a continuous low-cloud process causing spring temperature forecast errors in Beijing[J]. Meteor Mon, 48(6): 705-718 (in Chinese). DOI:10.7519/j.issn.1000-0526.2022.030702
沈学顺, 苏勇, 胡江林, 等, 2017. GRAPES-GFS全球中期预报系统的研发和业务化[J]. 应用气象学报, 28(1): 1-10. Shen X S, Su Y, Hu J L, et al, 2017. Development and operation transformation of GRAPES global middle-range forecast system[J]. J Appl Meteor Sci, 28(1): 1-10 (in Chinese).
王雨, 2003. 2002年主汛期国家气象中心主客观降水预报对比检验[J]. 气象, 29(5): 21-25. Wang Y, 2003. Verification of NMC subjective and objective precipitation prediction during the main flood season in 2002[J]. Meteor Mon, 29(5): 21-25 (in Chinese). DOI:10.3969/j.issn.1000-0526.2003.05.005
韦青, 代刊, 林建, 等, 2020. 2016—2018年全国智能网格降水及温度预报检验评估[J]. 气象, 46(10): 1272-1285. Wei Q, Dai K, Lin J, et al, 2020. Evaluation on the 2016-2018 fine gridded precipitation and temperature forecasting[J]. Meteor Mon, 46(10): 1272-1285 (in Chinese). DOI:10.7519/j.issn.1000-0526.2020.10.002
许晨璐, 王建捷, 黄丽萍, 2017. 千米尺度分辨率下GRAPES-Meso4.0模式定量降水预报性能评估[J]. 气象学报, 75(6): 851-876. Xu C L, Wang J J, Huang L P, 2017. Evaluation on QPF of GRAPES-Meso4.0 model at convection-permitting resolution[J]. Acta Meteor Sin, 75(6): 851-876 (in Chinese).
杨晓玲, 丁文魁, 马中华, 等, 2016. 河西走廊东部强降温变化特征和典型环流型[J]. 气象, 42(6): 756-763. Yang X L, Ding W K, Ma Z H, et al, 2016. Change characteristic and typical circulation type of strong cooling in Eastern Hexi Corridor[J]. Meteor Mon, 42(6): 756-763 (in Chinese).
赵滨, 张博, 2018. 邻域空间检验方法在降水评估中的应用[J]. 暴雨灾害, 37(1): 1-7. Zhao B, Zhang B, 2018. Application of neighborhood spatial verification method on precipitation evaluation[J]. Torr Rain Dis, 37(1): 1-7 (in Chinese).
赵声蓉, 赵翠光, 赵瑞霞, 等, 2012. 我国精细化客观气象要素预报进展[J]. 气象科技进展, 2(5): 12-21. Zhao S R, Zhao C G, Zhao R X, et al, 2012. The Development of objective meteorological element forecast in China[J]. Adv Meteor Sci Technol, 2(5): 12-21 (in Chinese).
Chen Y, Liao Z, Shi Y, et al, 2021. Detectable increases in sequential flood-heatwave events across China during 1961-2018[J]. Geophys Res Lett, 48(6): e2021GL092549.
Davis C A, Brown B G, Bullock R, et al, 2009. The Method for Object-Based Diagnostic Evaluation (MODE) applied to numerical forecasts from the 2005 NSSL/SPC spring program[J]. Wea Forecasting, 24(5): 1252-1267.
Gilleland E, Lee T C M, Gotway J H, et al, 2008. Computationally efficient spatial forecast verification using Baddeley's delta image metric[J]. Mon Wea Rev, 136(5): 1747-1757.
Hamill T M, Scheuerer M, Bates G T, 2015. Analog probabilistic precipitation forecasts using GEFS reforecasts and climatology-calibrated precipitation analyses[J]. Mon Wea Rev, 143(8): 3300-3309.
Liao Z, Chen Y, Li W, et al, 2021. Growing threats from unprecedented sequential flood-hot extremes across China[J]. Geophys Res Lett, 48(18): e2021GL094505.
Yuan W H, Hu X L, Li Y S, 2020. Evaluation of the hourly rainfall in the ECMWF forecasting over southwestern China[J]. Meteor Appl, 27(4): e1936.
Zhai P M, Zhang X B, Wan H, et al, 2005. Trends in total precipitation and frequency of daily precipitation extremes over China[J]. J Clim, 18(7): 1096-1108.
Zhang H, Zhai P M, 2011. Temporal and spatial characteristics of extreme hourly precipitation over eastern China in the warm season[J]. Adv Atmos Sci, 28(5): 1177-1183.