世界气象组织(WMO)为了促进气象科学的发展,将天气和气候预报更加广泛地应用于世界各国,尤其是发展中国家,从而提升全球人民福祉,近年来特别加强和突出了对世界气象中心的建设。在1967年认定了美国、俄罗斯和澳大利亚等3个世界气象中心之后,2017年,欧洲中期天气预报中心(ECMWF)、英国、加拿大、日本和中国成为第二批被认定的世界气象中心。
WMO作为联合国旗下各国气象部门的协调机构,一直关注作为现代天气预报能力核心工具的各国业务数值天气预报模式能力的提升。在2011年第16次WMO大会上,通过了在基本系统委员会(CBS)建立一个评估确定性NWP预报的验证中心(WMO LC-DNV)的决议,ECMWF被指定为该中心的责任机构,实时和分阶段针对各国NWP模式指标给出汇总和评价。本文所指的“主要气象业务中心”,包含了截至2018年7月,在LC-DNV网页上针对各国业务预报模式能力进行实时验证的约10个国家气象业务中心。目前,进入实时验证结果监测的国家气象业务中心的数量还在增加。本文借助近年来针对主要气象业务中心模式预报能力的对比分析结果,用核心预报指标梳理了这些中心预报能力的进步,并在此基础上结合“领先水平”的发展痕迹,特别是过去10年里预报预测水平提升的幅度,在给出当前水平(2016—2017年)的同时,对未来10年(2025年)和20年(2035年)国际上预报预测的先进水平将可能达到的高度进行讨论。经典天气预报理论中,洛伦茨基于大气混沌的本质,认为天气的可预报性为16.8 d。随着模式模拟技术的发展、学界对非线性过程认识的进步以及计算、模拟技术改进带来的手段更新等,可预报性问题得到了更加清晰的解析和拓展。Buizza等重新综合审视可预报性时,从模式与预报对象的时空尺度出发,给出了新的可预报性2D分布(图1)。
图1 洛伦茨经典可预报理论(a)和最新研究给出的可预报性2D分布(b)
全球模式是气象业务中心开展确定性预报的核心工具,基于以往研究,表1梳理了目前部分气象中心业务数值预报模式的情况。预报系统性能领先的ECMWF,其模式水平分辨率达9km,垂直层数为137层,预报时效在10~15d。2017年,ECMWF将其综合预报系统(IFS)升级为IFS Cycle 43r3:改进了对流模拟,采用了新的辐射方案和气溶胶气候学,并且更好地利用了下投式探空仪数据和其他观测数据。未来5~10年,随着各国气象部门业务主模式的更新换代以及高性能计算机的发展,全球NWP业务模式分辨率将达到千米尺度。其中ECMWF的战略目标是,到2025年以5 km的分辨率进行全球集合预报。2.2 ECMWF 为代表的全球领先模式核心指标现状ECMWF给出的截至2017年底各种核心业务指标统计中,影响最大的预报产品——北半球中纬度24h降水集合预报的进步明显,从21世纪初的2.5d增加到7d左右。尽管如此,对比21世纪前10年和最近10年的预报进步幅度,可以发现预报的改进速度明显“减缓”(图2),且在更多变量和时次预报中也是如此,但变缓的进步仍使10d左右的预报具有价值。
图2 ECMWF北半球中纬度24 h降水总量集合预报有效天数变化(2002—2017年)(图中给出预报时效(CRPSS≥0.1)的12个月滑动平均值)从全球主要气象中心的预报指标演变和相互比较中看出,21世纪以来ECMWF的发展趋势也代表了世界预报业务整体上的进步。从图3中多个气象中心指标进步速度的比较可以看出,21世纪以来,第二个10年与第一个10年相比,各主要气象业务中心的预报系统性能提升速率都明显放缓。

图3 世界上6个主要气象中心2002—2017年业务系统北半球中纬度500 hPa位势高度预报有效(距平相关≥80%)预报的天数变化,图中还给出2个再分析结果用于参考
3 全球NWP模式2025和2035年性能进步预测对世界上主要气象业务中心在未来10年甚至更长时间内预报水平提升幅度的估计,有多种方法和渠道。首先,各国气象部门在其战略中,定性或定量提出未来核心指标的改进目标,例如ECMWF制定的2016—2025年战略,在预报方面明确提出了“2+4+1”战略,即到2025年高影响天气的有效集合预报提前2周,大尺度形势和机制转化预报提前4周,全球尺度异常预测提前1年。NOAA天气气候预报机构NCEP提出,未来5年部分预报系统升级的主要特征为:水平分辨率进入10 km,区域模式和预警服务进入1~3 km,集合成员数稳定在20~30个,最长预报时段超过1年,达到15个月。此外,在最新发表的一些评述文章中,基于过去模式性能改进的技术和科学推动以及未来新技术在获取新资料和改进模式分辨率等方面发力的分析,给出了未来预报模式改进的乐观估计,以及面临的挑战(图4)。

图4 未来数值天气预报的关键挑战领域。
预测技术的进步来自计算科技的革新、
物理过程的参数化表达、
地球系统构件的耦合、
先进观测资料同化算法的使用、
通过集成方法对不确定性的一致描述
以及它们在不同尺度的相互作用等。
椭圆部分表示数值模式中解析的
10-2~104 km尺度范围内的主要现象
并展示出从小尺度范围流动
到完全耦合地球系统的模拟复杂性;
框表示未来预测技术改进遇到的最重要挑战;
箭头表示跨分辨率尺度和地球系统显现的误差传播重要性
3.1 基于全球预报过去 20 年的改进 :2025 和2035 年预报时效将分别提高到 8.5 和 10.5 d图5分别给出过去22和18年ECMWF与NCEP全球预报系统的预报时效的改进进程。ECMWF当前5 d预报与22年前的2 d预报相当,而目前的15 d预报与过去的10 d预报相当;NCEP在截止2017年的过去18年里,预报时效提高了4 d。综合两个机构的预报时效提升速率,即如果平均每10年提高近2d的进程在未来20年里持续,那么目前主要预报指标大约6.7d的时效,在2025和2035年可分别提高到8.5和10.5d。

图5 ECMWF(a)和NCEP(b)预报系统对500 hPa位势高度预报指标的进步历程
3.2 基于高分辨率模式过去 22 年的改进 :2025和 2035 年时效分别提高到 7.6 和 8.4 d高分辨率模式是目前各国气象部门、研究机构和企业重点研发的预报工具,其意义不仅在于大幅度提升观测大数据和区域预报两端的数据吸纳和服务能力,更在于未来人工智能(AI)更多地介入NWP后,高分辨率模式可成为支持全球模式通过机器学习解读更多模拟中“灰色”地带的有效手段。图6给出过去20年,ECMWF高分辨率模式系统预报时效的改进:最近10年预报时效增加了0.8 d,从大约6.0 d提高到6.8 d。未来10~20年,预报时效持续进步将使高分辨率预报系统的时效分别提升到7.6和8.4 d,接近同期集合预报水平。
图6 ECMWF高分辨预报系统HRES对500 hPa位势高度预报的进步历程
3.3 多种模式要素预报 :过去 10 年进步在 1~2 d图7给出过去20年,ECMWF在温度和降水预报上的改进:其中最近10年预报时效增加了0.8~2 d,这些要素在未来10~20年持续改进,将让这些关键预报分别得到可期待的1~2 d和2~4 d的改进。对比前面针对可预报性的讨论,到2035年前后,各类集合和确定性模式预报的水平,预报有效性将全面接近和进入10 d,也更加接近可预报性的极限。
图7 ECMWF 1998—2017年850 hPa温度(a)、降水总量(SEEPS验证,b)和降水(分级概率技巧,c)预报进步历程和最近10年的进步幅度
1967年,Manabe等论述一维辐射对流平衡模式的论文,被认为是现代数值气候模式研究的起点。在20世纪随后的时间里,气候模式被加入了更多的要素,模式的结构和计算也走向精细化。20世纪90年代起,世界气候研究计划(WCRP)耦合模拟工作组(WGCM)组织开展耦合模式比较计划(CMIP)。前两个阶段影响范围局限在气候模式界,第三阶段(CMIP3)起,模式结果的公开性大大增加,开始为其他研究团体所用。2008年启动的CMIP5,已成为政府间气候变化专门委员会(IPCC)评估报告的重要参考。当前试验阶段为CMIP6,共44家机构参与。参与CMIP6的气候中心及其耦合模式或气候系统模式的情况,能够给出当前此类模式的全球概览。世界天气研究计划(WWRP)和WCRP S2S项目选择的关键指标,一为MJO预测技巧,二为天气型转变预测。在2013—2017年的第一阶段实施中,10个S2S模式中,有7个MJO提前20天预报的双变量相关技巧(bivariate correlation skill)超过0.5,只有1个模式(ECMWF)预报时效达到30天(图8);在天气型转变预报方面,北大西洋涛动(NAO)正负转型预报可提前3周,其他天气型的转化提前16天预报。
图8 10个S2S模式模型集合均值与ERA Interim之间的MJO双变量相关性相对提前期的演变
根据ECMWF模式对MJO预报时效的发展,2006—2016年,预报时效(双变量相关≥0.5)10年约增加12 d。如保持此进步速度,到2025年将达到49 d。在2004年NCEP召开的纪念全球NWP业务开展50周年的学术会议上,Lord对那时NCEP的预报水平进行了阐述,用Niño3.4海温距平指标,给出从1997年冬季到2003年冬季的历史回报结果(图9a),其中提前3和6个月的相关性分别约为85%和76%(取6个预报指标的中数)。同样的指标,当前预报也有很多可比较的量化结果,例如气候系统历史预报项目(CHFP)集合世界9个中心对1992—2010年历史回报3和6个月技巧分别能达到90%和82%(图9b)。这样的对比结果表明,在Niño 3.4海温距平指标为代表的短期气候信号预测方面,过去15年提前3和6个月的预测水平,相关性分别提高了5%和6%。
图9 2000年前后NCEP(a)和CHFP项目给出的9个主要气候预测中心(b)针对Nino3.4海温距平指标的历史回报结果
综合考虑,上述相差大约15年的气候预测指数的进展,以及目前指标已经非常接近100%的情况,未来3和6个月针对Niño 3.4海温距平预测的指标,分别可能达到93%和86%(2025年)以及96%和90%(2035年)的水平。天气和气候模式的发展,主要以业务预报需求为导向。在当前预报时效已经接近可预报性极限的情况下,借助高性能计算新技术和AI等新手段和思路,以及围绕预报开发新产品支持各种业务服务,成为推动模式预报有效性及其应用的主要手段。其中,新预报量的设计和业务化、下一代预报模式和资料同化技术的研发等,成为目前NWP领域发展的新趋势。ECMWF最新给出了世界领先的5个气象业务中心2018年和2017年3—5月预报指标的比较(图10)。其中,在几乎所有的预报时段,1年里的进步都是“可测量”且是正面的。ECMWF在不断改进和升级其业务预报模式的同时,2018年,在其模式输出端新增加了3个重要变量的模式预报场:闪电、总水汽输送和最大对流有效位能,以支持更多的模式应用和研究。这些新的预报产品提升了服务(如闪电)和天气诊断(如总水汽输送和对流有效位能)利用的可能。
图10 2017和2018年ECMWF(红色)、UKMO(深蓝)、NCEP(绿色)、JMA(橙色)和CMC(浅蓝)对北半球中纬度地区500 hPa位势高度预报技巧得分对比
5.2 2025 和 2035 年两个时段可期的预报指标针对学者最新给出的可预报性2D分布图,基于文献研究成果,本节选取了反映世界气象中心的核心业务指标,给出当前以及未来10和20年可能的发展趋势(表2)。在分析这些指标过去10~20年进步的基础上,同时考虑了指标接近可预报极限,增幅减缓的可能。表2中反映气象业务中心整体能力的所谓核心指标,实际上是气象业务中心一系列研发工作及其综合能力在一个总体上的考量。要有效地提升核心能力指标,WMO天气公开科学会议提出至少以下8个方面高水平的研发工作:2)重大国际外场试验和观测技术(雷达、廓线仪和卫星等)复杂研发互动;表2 世界领先水平核心预报预测指标目前值和未来(2025年和2035年)估计