T639数值预报产品订正方案.docxVIP

Ｔ６３９数值预报产品订正方案任文斌１杨新２孙潇棵３陈磊４管习权２（１广东省揭西县气象局，揭西５１５４００；２广东省揭阳市气象局，揭阳５２２０００；３湖南省宁乡县气象局，宁乡４１０６００；４浙江省慈溪市气象局，慈溪３１５３００）摘要利用２０１０年１—１２月的Ｔ６３９数值预报产品作为预报形势场，以１９８０—２００９年共３０年的ＮＣＥＰ资料作为历史样本场。利用动态相似统计方法，对Ｔ６３９的预报产品，从３０年的历史样本场中找相似个例。选取因子，采用多元回归法，对Ｔ６３９的预报产品进行订正方案设计。用２０１１年１—６月的Ｔ６３９预报产品对方案的订正效果进行检验。对比分析发现，方案对东北半球３天以上的高度场和４天以上的温度场订正效果明显，高度场的订正效果优于温度场。关键词数值预报Ｔ６３９Ｌ６０模式动态相似统计误差订正引言Ｔ６３９Ｌ６０全球中期数值预报模式是对中期谱模式Ｔ２１３Ｌ３１升级而来的，成功引入ＥＣＭＷＦ的 稳定外插半拉格朗日方案，将归约高斯格点转换为线性高斯格点，解决了模式分辨率提高的关键技术，于２００７年１２月１４日通过准业务运行评审［１］。目前Ｔ６３９数值预报产品通过９２１０系统向全国业务单位广播，在天气预报业务中得到广泛的应用。 然而模式描述的并非真实大气，输出产品总会有误差［２－７］。Ｔ６３９数值模式自２００８年投入业务运行以来，很多学者对其预报效果进行了检验评估［８－１１］，而 对其产品的误差订正研究较少。以往的中期数值模 式误差订正方法有经验修正、ｋａｌｍａｎ滤波、集合ｋａｌｍａｎ滤波、人工神经网络ＢＰ模型、支持向量机、 变分等方法［１２－１６］，张兰慧等人也提出了根据模式前几个时次的误差计算出后面预报时次误差的新方法［１７］。数值预报误差的订正一般都在误差统计分析的 基础上进行，因为统计分析法是改进数值模拟、数值预报所不可缺少的分析研究方法。用统计分析法从 大量的历史气象资料中寻找大气现象间的统计规律或相互关系，能从统计的意义上找出影响某种天气或天气过程的主要因素，甚至可发现一些理论分析中难以推断出的事实［１８－２５］。基于这一思想，本文利用动态相似统计法，在历史资料中找相似个例，利用相似个例构造因子建立回归方程，对Ｔ６３９模式产品进行订正方案设计，从统计角度对Ｔ６３９中期数值预报产品的误差订正提供一种尝试。本文动态相似统计法的动态有两层含义：①相 似个例的选取是动态的，它们是根据逐日Ｔ６３９产 品而定的；②用于挑选相似个例和组建回归方程的历史样本资料是动态的［２６－２８］。研究中取Ｔ６３９产品的５００ｈＰａ高度场和８５０ｈＰａ温度场为预报形势场，取Ｔ６３９产品对应前后各１５天同期３０年的ＮＣＥＰ资料作为历史样本场。在进行动态相似挑取过程中，对于高度场形势相似大于要素相似，而对于温度场要素相似和形势相似对总体的相似性贡献一样。１资料介绍１．１２０１０年１月至２０１１年６月的Ｔ６３９资料Ｔ６３９数值预报产品按每天２次（世界时， ００：００、１２：００）发布，提供未来１０天的预报。以 ｇｒｉｂ１版本编码形式存储在１．０°×１．０°的网格点上。 目前通过９２１０发布业务上使用的Ｔ６３９产品的网格范围为０°～９０°Ｎ，０°～１８０°Ｅ，包含东北半球共１６４７１（１８１×９１）个格点［２９－３０］。每一个文件包含一ｈｔｔｐ：／／ｗｗｗ．ｑｘｋｊ．ｎｅｔ．ｃｎ气象科技作者简介：任文斌，男，１９８４年生，学士，工程师，主要从事地面观测及短期预报工作，Ｅｍａｉｌ：ｒｅｎｗｅｎｂｉｎ＿２００８＠１６３．ｃｏｍ收稿日期：２０１２年１１月２日；定稿日期：２０１３年４月１０日１４６气象科技第４２卷个预报时次的全部高度层所有要素的资料。文件名形式如下：Ｚ＿ＮＡＦＰ＿Ｃ＿ＢＡＢＪ＿２０１１０５０４００＿Ｐ＿ｇｍｆ－６３９－２０１１０５０４０００２４．ｇｒｉｂ１。该文件名表示是以２０１１年５月４日００：００为初始场起报的２０１１年５ 月４日２４ｈ的预报产品，即为５月５日的预报产品。Ｚ＿ＮＡＦＰ＿Ｃ＿ＢＡＢＪ＿２００８０７０１００＿Ｐ＿是通信标识符，２０１１０５０４００表示以２０１１年０５月０４日００：００为起报时刻，后缀表示为ｇｒｉｂ１格式。１．２１９８０—２００９年ＮＣＥＰ资料ＮＣＥＰ／ＮＣＡＲ资料为美国国家环境预报中心和国家大气研究中心联合发布的再分析资料，它能为中尺度数值模式提供初始场和侧边界条件。该资料每天有４个时次的值（即世界时００：００、０６：００、 １２：００、１８：００的资料），用ｎｅｔＣＤＦ格式以压缩二进边界上的点用周围４个最近的点插出。如果Ｔ６３９产品在插值格点上有值，则直接代替。计算距离时， 因为地球上不同纬度两经线之间距离不相等，所以需给ｒ乘纬度的余弦做为修正。３订正方案设计３．１动态相似法选取相似个例用２０１０年１月至２