卫星辐射率资料资料直接同化研究及初步试验研究.pdfVIP

卫星辐射率资料资料直接同化研究及初步试验 张华 薛纪善吴雪宝 张风英 朱国富庄世宇 马 刚 摘要 1 前言 随着探测技术的发展，探测系统以常规探空观测为主的时代已成为过去，而以卫星、 雷达、飞机、风廓线仪等非常规观测为主的时代已经到来。这一变化使资料的覆盖面与时 空分辨率大大提高。这些新的探测要素可能具有更好的代表性，但也在资料的处理和同化 分析方面带来新的问题。代表着当今国际资料同化技术发展趋势的变分同化正是能够适应 上述要求的一种有效方法。 目前卫星反演资料对NWP有负面影响，使北半球预报质量降低。因此，人们在改进 探测器的垂直分辨率和全天候探测能力的同时，纷纷寻求在NWP中应用卫星资料的新途 径和新方法。辐射率直接同化是目前最具影响和潜力的方法，该方法避免计算复杂的不适 定问题所掺杂的反演计算误差，实现用正演方法求解反演问题，从方法论上避开反演问题 的复杂性。为了实现卫星辐射率资料直接同化及其它非常规观测资料的同化，中国气象科 学研究院数值预报研究中心先后开发研制了GRAPES三维变分同化系统，将ECMWF开发 的第六代快速辐射传输模式移植到该三维变分同化系统并通过严格地正确性检验，成功地 实现了整条轨道卫星辐射率资料直接同化，并将常规探空观测资料和卫星辐射率资料融合 到一起。初步试验结果显示，卫星探测辐射率直接同化方案是合理、有效的。 2 GRAPS三维变分同化系统 2．1目标函数 三维变分同化方法的基本思想是将目标函数极小化。在本方案中，目标函数定义如下： =J^+d。 其中 x 长度为Ⅳ的分析变量向鬣 p K度为^，的背景变量向鬣 49 少 长度为肌的观测向量 。 口。 NXN背景误差协方差矩阵的逆 口’ MXM观测误差协方差矩阵的逆 Ⅳ 观测算子 Ⅳ 分析场的自由度数目 膨 观测场数目 根据Courtier(1994)嘲、本方案采用增量分析法．目的是减少每一个迭代步的计算 量，令 Dx=x一％ 将观测算子H泰勒展开。取其前两项 小工)=Ⅳ(％+盘)角H(x^)+日钕 令 d=H(xb)一y。 (2．2) k=兰蠡r廖‘‘盘+兰(日慑+drD_(Ⅳ撅+d) 其梯度为 V盘J=B’1蠡+Ⅳ’rO’。似钕+d) (2．3’ V三，=B-’+H盯0_‘H’ (2．4) 2．2控制(分析)变量 选择控制变量应遵循以下原则： a、变蟹之间不相关： b、考虑计算开销．即当考虑观察算子y=H)=日(，￡肺)时，变换c，必须相对简单． 根据上述原则．控制变量取为 a) 流函数v．非平衡速度势函数x。 x-=x·x# (2．5) 在本方案第一版本中·暂不考虑x小即令x一 功 非平衡位势咖． 吼=口一吼 (2．6) 咖。是通过某种平衡关系得到的．目前仅考虑地转平衡．哦和风场在统计上是不 相关的。 d 比漫q 3 观测算子一RTTOV6快速辐射传输模式 观测算子的主体采用由ECMWF开发的用于在数值预报中同化卫星资料的第六代快速 度、湿度、03、云中液态水垂直廓线以及若干云和表面参数作为输入变量，计算卫星相应 通道的辐射值。 模式采用大气辐射传输方程的近似形式，不考虑散射效应，即 L(v，0)=(1一Ⅳ)r’扣(v，0)+Ⅳ∥埘(y，臼)