张东方,成青燕,2*,何慧根,杨 蓉,韦 玮,汪 玲
(1.成都市气象局,成都 611133;2.高原与盆地暴雨旱涝灾害四川省重点实验室,成都 611133;3.重庆市气候中心,重庆 401147)
长期以来,公众所熟知的气象服务是3d 以内的短期天气预报,随着社会经济的发展,中期预报甚至月、季、年尺度的短期气候预测逐渐成为公众气象服务所迫切需求的内容,目前中期数值预报产品最长预报时效是6~10d,而短期气候预测只考虑气象要素的变化趋势,并未细化到具体的时间分布或天气过程。10~30d 的延伸期预报填补了中期天气预报和短期气候预测之间的缝隙,在防灾减灾决策服务中起着重要作用,一直是科学研究和业务预报关注的重点和难点[1-2]。近10 多年来,诸如2008 年南方低温冰冻雨雪天气、2009 年西南五省持续干旱、2013年初雾霾笼罩全国、2017 年北方大范围极端高温天气等阶段性、持续性高频发高影响的异常天气事件越来越多,而延伸期预报工作的薄弱及产品的匮乏也使得气象部门在相关决策服务工作中略显被动。中国气象局印发的《全面推进气象现代化行动计划(2018—2020 年)》和《智能网格预报行动计划(2018—2020 年)》,已明确提出将气象要素和重要天气过程的延伸期预报列为重点任务和着力攻关的关键核心技术之一。国内学者和业务单位在延伸期预报方法的探索和业务产品的研发领域持续努力,取得可喜成效,其中包括应用最广泛的上海市气候中心研发的低频天气图方法[3-5],国家气候中心的月动力延伸模式DERF2.0 形势场预测及降尺度释用[6-8],河北省气候中心异常相似(释用)方法[9-10],西北区域气候中心的准150d 韵律方法[11]以及美国环境研究中心以及欧洲数值预报中心的延伸期预测。前4 种方法作为延伸期预测的重要手段,前期工作繁琐,需建立统计模型,且有一定的地域局限性;欧洲数值预报中心的延伸期预测产品属于非公开,不能通过一般途径得到。美国环境预报中心延伸期模式CFSv2 提供逐日向后延伸45d 的要素预报,本文利用其公开预测产品,直接进行双线性插值对成都地区未来45d的逐日气温和降水进行预报应用,对近5a 应用效果进行全面检验和评估,从而为该模式在成都市延伸期预报中有效利用提供理论依据。通过检验和评估提出订正方案,可以提高该模式在成都地区延伸期预测的准确率。
本文所用资料主要包括:(1)2016—2020 年逐日CFSv2 地面降水、850hPa 气温未来45d 预报结果,以及相应气候平均值(1991—2010 年)。(2)成都地区14 个国家气象站2016—2020 年逐日的降水观测资料及对应的气候平均值。
本文采用双线性插值法[12]将格点数据插值到站点。采用距平相关系数ACC、趋势异常综合检验Ps对模式预测降水的总体趋势进行评估,采用Zs、Cs评分对降水过程的预测进行定量评估。其中空间距平相关系数ACC、Ps 方法在文献[13]中已进行了详细描述,Cs/Zs 评分方法详见文献[14],本文不予赘述。
由于模式每日输出未来45d 预测,为配合本地短期气候预测业务使用,本文仅以每月28 日起报预测结果作为检验对象。
距平相关系数ACC 常用来描述某一气象变量距平的预测值与观测值的年际变化或空间分布的相似程度,广泛应用于动力模式和统计方法对不同区域范围不同变量的预测评分,是世界气象组织(WMO)确定并建议使用的指标。
图1 为CFSv2 模式在2016 年1 月—2020 年12月间逐月28 日起报的未来45d 降水与实况的距平相关系数ACC。60 个样本中,模式预测与实况数据的距平相关系数ACC 在-0.6~0.6 之间,其中有21个样本为负值,20 个样本ACC>0.30,最大值为2019年3 月28 日起报结果,预报与实况ACC 达到0.63。全部样本的ACC 平均值为0.13,虽然评分较低,但高于2016—2020 年四川省气候中心发布的月降水预测的ACC 年平均评分0.06,且高于同时段国家气候中心指导发布的月降水预测的平均ACC 评分0.05(评分来自中国气象业务内网气候预测产品检验频道http://10.1.64.154/portal/web-link_ycpg.index#)。由此表明,CFSv2 模式对成都地区的未来45d 降水预报效果略高于目前短期气候预测的业务水平,且由图1 可知,2016—2020 年ACC 正值频次逐渐增加,一定程度表明CFSv2 模式的预报准确率近5a 逐步提高。
图1 逐月延伸期降水预报与实况的距平相关系数ACC
距平相关系数在成都地区不同季节的空间分布如图2 所示,在不同季节,CFSv2 模式对成都地区降水的预测性能都存在一定的地区差异,春季西南部和东南部ACC 高于其它地区,夏季ACC 高分区主要在东南区,秋季西部和东南部为高分区,冬季除东部外ACC 基本都大于0.2。总体上,成都西南部和东南部的ACC 评分高于其它地区。
图2 不同季节距平相关系数的空间分布图
趋势异常综合检验Ps 评分法是中国气象局为促进气候预测业务由平均态趋势预测向异常气候趋势预测转变制定的,同时考虑到预报的趋势项、异常项和漏报项,于2014 年开始应用于业务评分。如图3,CFSv2 模式对成都地区延伸期降水的趋势异常综合检验Ps 评分的月际变化明显,季节差异不大,没有明显的持续较差或持续较好。60 个样本中评分最低为30 分,有20 个样本Ps 低于60 分,最高为94.7分,平均为65.7 分,与四川省气候中心发布的2016—2020 年月降水预测的Ps 平均评分65.4 基本一致,略高于同时期国家气候中心发布的指导预报Ps 平均评分61.7(评分来源于中国气象业务内网气候预测产品检验频道)。2016—2020 年,Ps 评分总体趋于稳定并逐渐升高。
图3 逐月降水的趋势异常综合检验Ps 评分
从Ps 评分的空间分布图来看,同样存在显著空间差异,但并无明显变化规律,春季西南部PS 评分显著高于其它地区,夏季与春季基本呈相反的变化,秋季PS 高分区位于成都市南部,冬季Ps 高分区分布在西部及偏北地区。总体上西南部PS 评分略高于其它地区,这与距平相关系数ACC 评分得到的结论具有相似性。
利用850hPa 变温为负、降水量非0 作为有降水的判断依据,认为同时有5 站或以上站点产生降水即为1 次明显的降水过程。表1 为对2016 年1 月—2020 年12 月期间降水过程Cs、Zs 评估,由表可知,结合负变温和降水,CFSv2 模式对成都地区的延伸期降水过程纯客观预报结果具有一定的参考价值。Cs 评分重在考查对延伸期降水过程的日期预报效果。60 个样本中,评分最高出现于2016 年4 月28 日的预报,Cs 为0.80;最低评分0.23 为2018 年2 月28日起报。其中有25 个样本Cs 值大于0.50,14 个样本Cs 值小于0.3,成都地区平均的Cs 评分为0.48。Zs 评分重在考查对延伸期内主要降水过程的把握,因此明显高于Cs 评分,60 个样本平均的Zs 评分为0.64,其中有22 个样本的Zs 值大于0.7,2017 年2月28 日起报的Zs 评分高达0.95。无论是Cs 评分还是Zs 评分,从高评分的分布来看,模式对延伸期内降水过程的预报并没有明显的季节分布特征。
表1 CFSv2 模式对月内主要降水过程的Cs/Zs 评估
图4 不同季节PS 评分的空间分布
空报率和漏报率也是反映降水过程预报效果的重要指标,参照中国气象局《中短期天气预报质量检验办法》(气发〔2005〕109 号)计算月内主要降水过程的空报率和漏报率,CFSv2 模式对成都地区延伸期主要降水过程的空报率为38%,漏报率为34%。从高空报和高漏报的分布来看(表略),冬半年空报较多,夏半年漏报较多。总体评分逐年提高。
综合以上应用和检验可知,CFSv2 模式对成都地区的延伸期降水预测具有一定程度的预报能力,但是作为延伸期预报的主要参考依据,需要对模式结果进行订正。不同于气温的连续变化,降水是大气环流系统之间共同协调作用产生,数据更加离散化,单纯依靠模式对于降水的预测,准确率是不确定的,想要建立一个参数恒定的订正模型比较困难,也不利于准确率的提高,在消除偏差订正方案[14]的基础上,将观测值平均值改为模式偏差,然后采用一元线性回归的方法进行误差订正,具体用公式表示为:
误差为预报值与实况值之间的偏差:
式中,i 为起报时间,j 为第j 个预报时次,bij 为i 起报时间预报的j 时次的预报误差,fij 为i 起报时间预报的j 时次预报值,aij 为j 时次的实况值。
对于离散数据而言,误差可表示为:
式中,w 为权重系数,也称为回归系数。利用最小二乘法原理,通过一元线性回归法,计算权重系数。
经过订正后的预报可表示为:
即模式的订正预报为模式的初始预报-模式误差。
根据公式(1),对2016—2020 年逐月28 日起报的逐6h 降水数据进行误差计算。再利用公式(2),基于最小二乘法,求出各个时次的误差权重系数。利用公式(3)计算订正后的预报值。
订正后模式预报降水与实况的距平相关系数ACC 与订正前相比,以提高为主,其中2017 年5 月28 日起报的降水与实况距平相关系数较订正前提高了0.37,其余月份提高不明显,基本小于0.05,个别月份ACC 较订正前有所下降。订正后平均的ACC较订正前升高0.04。订正后的ACC 空间分布与订正前一致,西南部和东南部略高于其它地区(图略)。
订正后模式预报降水与实况的Ps 评分与订正前相比,总体提高较明显,其中2020 年5 月28 日起报的模式降水与实况的Ps 评分较订正前提高了51分,另有4 个样本的订正后Ps 评分较订正前提高幅度在20~35 分之间,有20 个样本为负订正,订正后的Ps 评分较订正前降低1~10 分,其余各月均以正订正为主,Ps 评分的提高幅度在1~20 分之间。Ps 评分的空间分布订正前后并无显著差异(图略)。
表2 为订正前后模式预报的降水过程的Cs/Zs评分差异,无论是对主要降水过程的把握还是对降水过程发生日期的把握,订正后的结果大多好于订正前,个别月份存在负订正。订正后Cs 评分平均提高0.05 分,Zs 评分平均提高0.14 分。订正后空报率明显下降,较订正前下降22%,漏报率下降5%。
表2 订正前后CFSv2 模式对月内主要降水过程的评估差异
综上,利用滚动式一元线性回归的方法来对CFSv2 模式的延伸期降水预报进行订正,可以有效提高降水预报的准确率,对于模式资料在成都地区的直接应用具有重要的贡献。
利用2016 年1 月至2020 年12 月CFSv2 模式的延伸期降水预报结果进行成都市的延伸期降水预测试验,经过与实况对比发现:
(1)CFSv2 模式预报的延伸期降水与实况的距平相关系数ACC 平均为0.13,总体较低,但略高于省级、国家气象局短期气候预测业务水平。ACC 的空间分布反映出CFSv2 模式对成都西南部和东南部的预测效果好于其它地区。
(2)CFSv2 模式预报的延伸期降水与实况的Ps评分平均为65.7 分,与省级发布预测评分相当,高于国家级对四川地区的指导预报评分。
(3)CFSv2 模式预报的延伸期降水过程在成都地区平均Cs 为0.48,Zs 为0.64,空报率38%,漏报率34%。Cs 评分明显低于Zs 评分,模式对主要过程的把握好于对过程发生日期的把握。冬半年空报多,夏半年漏报多。
(4)改进的消除偏差一元线性回归对2016 年1月至2020 年12 月期间的模式预报进行订正计算和重新评估,订正后ACC 平均提高0.04,Ps 提高4.5分。订正后模式延伸期降水过程Cs 评分平均提高0.05 分,Zs 评分平均提高0.14 分。空报率较订正前降低22%,漏报率降低5%。
通过以上分析可以看到,经过订正后的CFSv2模式的延伸期预报结果可以作为成都地区延伸期降水预报的重要决策参考依据,但由于资料序列时间较短,本文仅采用双线性差值的方法将模式结果进行直接应用,还需要今后持续检验并结合动力与统计的降尺度方法作出改进,进一步提升延伸期预报服务能力。