基于切线伴随技术计算GRAPES-Meso 模式的奇异向量.PDF

26 卷 4 期 热 带 气 象 学 报 Vol.26, No.4 2010 年08 月 JOURNAL OF TROPICAL METEOROLOGY Aug., 2010 文章编号 1004-4965(2010)04-0421-08 基于切线伴随技术计算GRAPES-Meso 模式的奇异向量 1, 2 2, 3 1 刘永柱 ， 杨学胜 ， 王洪庆 (1. 北京大学物理学院大气科学系, 北京 100871; 2. 中国气象科学研究院灾害天气国家重点实验室, 北京 100081; 3. 国家气象中心, 北京 100081) 摘 要: 集合数值天气预报的关键问题就是如何生成有效的初始扰动。奇异向量反映了初始扰动在大 气系统相空间中演变发展的最不稳定方向，基于奇异向量产生的集合样本是模拟概率密度函数的最合理 方法。以非静力、半隐式半拉格朗日 GRAPES-Meso 中尺度数值预报模式为基础，采用 Lanczos 迭代算 法，利用GRAPES-Meso 的切线伴随模式计算GRAPES 奇异向量。为了检验求得的奇异向量的正确性， 提出了两种检验奇异向量正确性的方法：一是比较计算的奇异值的一致性；二是依据特征向量在矩阵变 换中的方向不变性意义，验证 GRAPES 奇异向量空间结构的正确性。最后研究了不同的时间间隔对 GRAPES 奇异向量的影响，结果表明 GRAPES 奇异向量在 36 小时的最优时间间隔误差增长速度最快， 这表明在非静力、半隐式半拉格朗日格点模式中利用切线伴随技术计算奇异向量是可行的。 关 键 词 数值天气预报；集合预报；奇异向量；初始扰动；GRAPES 模式；切线伴随 中图分类号 P456.7 文献标识码 A Doi ：10.3969/j.issn.1004-4965.2010.04.005 方法总是以分析资料作为模式的初始条件，所以当 1 言 资料处理的方法改变时，集合预报中初始扰动场的 设计方法也应有相应的变化。另一种方法是以欧洲 集合数值天气预报是用不同的初始条件，对数 中期天气预报中心(European Center for Medium 值模式积分得到一组数值预报结果的方法，其最大 -Range Weather Forecast ，ECWMF)为代表的奇异向 - 的难点在于初始扰动初值的产生方式。而不管采用 量法(Singular Vector, SV)[5 9] 。SV 方法有着完善的 何种方法，求得的初始扰动都应与基本气流的动力 数学理论基础，求得的一些主奇异向量能代表初始