基于增长模繁殖法的集合预报初扰饱和分析S5.pdfVIP

* 基于增长模繁殖法的集合预报初扰饱和分析 (S5) 于永锋1 张立凤1 2 1 解放军理工大学气象学院 南京 211101 2 中国科学院大气物理研究所大气科学和地球流体力学数值模拟国家重点实验室 北京 100029 摘 要 利用T63L9全球谱模式讨论了用增长模繁殖法求取集合预报初始扰动的过程中 增长模 扰动 的增长和饱和情况 通过分析增长率和形态相似系数随繁殖时间的演变 得到以下结论 扰动增长模在繁殖3 4天后可呈现明显的饱和特征 不同变量扰动增长模的饱和时间有所差别 同一变量在不同等压面上饱和时间亦不同步 不同的初始模大小也会对扰动 饱和过程产生影响 关键词 增长模繁殖 初始扰动 饱和 集合预报 1 引 言 [1] 经典的集合预报是进行初值集合 集合预报的理论和试验均表明 如何产生更有意义的初始扰动 场对集合预报的效果是至关重要的 上世纪 90年代 基于快速增长模 的思 [2],[3] [4] [3],[5] 想 形成了两种最具代表性的方法 奇异矢量法(SVs)和增长模繁殖法(BGM)而增长模繁殖法 在效果并不逊于奇异矢量法的前提下 又以其独特的构思 明确的动力学意义以及相对低廉的计算代价 受到越来也多的关注[6] 用BGM方法求取集合预报初始扰动的基本思想是 模仿分析循环中误差的增长过程 由一任意的初 始随机扰动出发 通过模式的循环积分使其得到 繁殖 我们可将该扰动在某种意义上达到的饱和态 作为对实际分析误差快速增长模的一个估计 并由此产生集合预报的初始扰动 具体繁殖步骤可参考文 献[3][5]初始扰动的生成是集合预报的关键技术 当用BGM方法生成初始扰动时 繁殖循环何时结 束也是此方法的关键技术 从增长模繁殖法的基本思想可知 这与饱和测度的定义有关 NCEP的试验 [5] 结果 表明 初始扰动经过3-4天的繁殖即可达到饱和 他们对 饱和 的定义是扰动增长率达到1.5/ 天 这一增长率接近实际短期预报误差的增长 BGM方法求取初始扰动是通过数值预报模式积分来实现 [7] 的 扰动达到饱和的时间 即繁殖长度 和饱和增长率也因模式而异 Mitchell等 的试验就表明好模 式中的扰动增长会更快一些 所以在利用BGM方法制作集合预报时必须针对集合预报所用模式 对繁殖 循环中扰动饱和的情况进行分析 以便确定繁殖循环时间 求取最佳的集合预报初始扰动 初始扰动的饱和不仅仅是扰动增长率的饱和 即繁殖模大小方面的饱和 为了寻找更能代表初始误 差分布的集合预报初始扰动 每一步循环中初始扰动的分布也是必须要考虑的一个因子 即扰动形态上 也必须达到饱和 所以扰动的饱和应包含两方面 大小的饱和和分布的饱和 关于误差分布的饱和的研 究工作还很少 是一个很值得研究的问题 本文针对所用模式首先对扰动增长率的饱和情况作了分析 然后定义了形态饱和的测度 对形态饱 和的情况进行了考察 这些研究是制作集合预报的基础工作 也是集合预报中的关键技术 2 模式选取和试验设计 本文所用模式为全球中期数值预报谱模式T63L9这是由解放军理工大学气象学院开发的微机版本 关于模式的详细介绍参见文献[8]试验所用资料为 1998年NCEP逐6小时资料 * 中国科学院大气物理研究所大气科学和地球流体力学数值模拟国家重点实验室开放课题资助 7 繁殖循环初始扰动的选取是实施增长模繁殖循环的第一步 也是用增长模繁殖法制作集合预报的第 一步 为反映初始场中的随机误差 参照文献[9]的做法 将该模式 24小时预报的均方根误差乘随机数 作为繁殖循环初始扰动,该随机数服从[-1,+1]区间上的均匀分布 其由下式确定