双分辨率集合卡尔曼滤波同化方案与模拟试验.PDF

第２９卷 第５期 气 象 与 环 境 学 报 Ｖｏｌ．２９Ｎｏ．５ ２０１３年１０月 ＪＯＵＲＮＡＬＯＦＭＥＴＥＯＲＯＬＯＧＹＡＮＤＥＮＶＩＲＯＮＭＥＮＴ Ｏｃｔｏｂｅｒ２０１３ 孙龙，乔小，蒋大凯，等．双分辨率集合卡尔曼滤波同化方案与模拟试验［Ｊ］．气象与环境学报，２０１３，２９（５）：４３－４８． ＳＵＮＬｏｎｇｙｕ，ＱＩＡＯＸｉａｏｓｈｉ，ＪＩＡＮＧＤａｋａｉ，ｅｔａｌ．ＴｈｅｄｕａｌｒｅｓｏｌｕｔｉｏｎａｐｐｒｏａｃｈｆｏｒｅｎｓｅｍｂｌｅＫａｌｍａｎｆｉｌｔｅｒａｎｄｓｉｍｕｌａｔｉｏｎｔｅｓｔｓ［Ｊ］． ＪｏｕｒｎａｌｏｆＭｅｔｅｏｒｏｌｏｇｙａｎｄＥｎｖｉｒｏｎｍｅｎｔ，２０１３，２９（５）：４３－４８． 双分辨率集合卡尔曼滤波同化方案与模拟试验 １ ２，３ ２ ２ ２ ２ 孙龙　乔小 　蒋大凯　陈力强　吴曼丽 　梁寒 （１沈阳市气象局，辽宁 沈阳１１０１６８；２沈阳中心气象台，辽宁 沈阳１１０１６６； ３南京信息工程大学，江苏 南京２１００４４） 　　摘　要：对于高分辨率中尺度模式集合卡尔曼滤波实际应用的困难是集合预报会耗费大量的时间，而双分辨率集合卡尔曼 滤波是由一组低分辨率样本提供同化所需的背景误差协方差矩阵，这种方法可以减少集合预报的时间。为了检验其有效性，利 用模拟资料，与标准高分辨率集合卡尔曼滤波方法比较。结果表明：在第一个同化时次，两者对５００ｈＰａ水平风场和扰动位温 场的分析增量场均与真实增量场的高低值中心位置一致，且结构与真实增量场接近，前者（高分辨率集合卡尔曼滤波）的增量 值比后者（双分辨率集合卡尔曼滤波）的增量值更接近真实情况。在连续的预报—同化循环试验中，随着同化次数的增加，两 种方法分析变量的均方根误差总体上都是下降的，均表现了较好的同化能力，但后者与前者相比仍存在一定的差距；而在相同 的运行环境下，后者的运行时间仅是前者的１／６。 关键词：集合卡尔曼滤波；分辨率；预报误差协方差矩阵；耗时 　　中图分类号：Ｐ４５６７　　文献标识码：Ａ　　ｄｏｉ：１０．３９６９／ｊ．ｉｓｓｎ．１６７３５０３Ｘ．２０１３．０５．００７ 辨率模式背景场的同化，所得分析场再作为高分辨 引言 率模式预报的背景场。正是由于增量法可以节省大 １９９４年 Ｅｖｅｎｓｅｎ提出了基于随机动力预测理 量的时间，四维变分才能得到广泛应用。在增量法 ［１］ 论 的集合卡尔曼滤波方法（ＥｎｓｅｍｂｌｅＫａｌｍａｎＦｉｌ ［８］ 的四维变分方法启发下，Ｇａｏ和Ｘｕｅ 提出了双分辨 ［２］ ［３－４］ ｔｅｒ，简称ＥｎＫＦ） ，与三维变分（３ＤＶＡＲ） 和四 率集合卡尔曼滤波同化法，这种方法既可以减少Ｅｎ ［５］ 维变分（４ＤＶＡＲ） 相比，它的主要优点是可以方便 ＫＦ分析的总耗时，又可以继承ＥｎＫＦ同化的优点。 地得到依流型而变的预报误差协方差，并且算法简单， 其主要原理是一组低分辨率样本为它本身和一个高 不需要编写切线性模式和伴随模式，因此该方法在大 分辨率样本提供同化所需的背景误差协方差矩阵。 气、海洋和水文资料同化的研究中受到广泛关注。 其基于ＡＲＰＳ（ＡｄｖａｎｃｅｄＲｅｉｇｉｎａｌＰｒｅｄｉｃｔｉｏｎＳｙｓｔｅｍ） 集合平方根滤波（ ［６］ ＥｎＳＲＦ） 是由ＥｎＫＦ衍生而 采用ＥｎＳＲＦ方案以美国德尔城，俄克拉何马地区一 来的有效的集合同化方法，可以避免ＥｎＫＦ方法中由 次典型的超级单体风暴为例，利用模拟雷达资料进 于扰动观测所造成的采样误差。但是与ＥｎＫＦ一样， 行同化试验，结果表明了这一方法是有效的。 集合数较小时，蒙特卡罗方法产生的采样误差则会