四维变分同化技术在风暴潮数值模拟中的应用.pdfVIP

第32卷 第 1期 海 洋 预 报 vo1．32．NO．1 2015年2月 M ARINE FORECASTS Feb．20l5 DOI：10．11737~．issn．1003—0239．2015．O1．001 四维变分同化技术在风暴潮数值模拟中的应用 王宗辰 ，于福 江 ，原野 (1．国家海洋环境预报中心，北京 100081；2．国家海洋局海洋灾害预报技术研究重点实验室，北京 100081) 摘 要：利用区域海洋模式ROMS(RegionalOceanModellingSystem)及其四维变分同化模块，建 立了具有资料同化能力的东中国海风暴潮数值模式，通过将海洋站水位观测资料同化到风暴潮模 式中，提高了模式对风暴潮的模拟精度。四维变分同化技术能够在整个同化时间窗口保持动力协 调 ，使模拟结果在该时间窗 口内最大程度 的靠近观测，同时，得到了最优预报初始场。利用该模 式，对两次温带风暴潮过程进行 了数值模拟，结果表明：在同化窗口内，同化对模拟精度有明显的 提高；结束同化之后，得到的最优预报初始场对临近预报精度也有一定提高。 关键词：区域海洋模式ROMS；四维变分同化；风暴潮数值模拟；最优预报初始场 中图分类号：P731．23文献标识码：A 文章编号：1003．0239(2015)01—0001—09 1 引言 资料序列时问间隔不等的情况下 ，模式仍可以顺利 运行得到结果。Sasaki[“第一次提出了满足控制方 目前风暴潮数值预报结果多是由确定性数值 程强约束条件的伴随方法，完成了四维变分同化的 模式计算得到，但是采用确定性风暴潮数值模式，有 基本理论。20世纪末 ，伴随方法被引入到气象数据 时模拟结果和实际情况相比仍然有较大的误差 。这 分析和同化领域n，之后 ，又被引入到全球海洋状态 些误差的来源主要包括3个方面：一是动力模式对 分析和海洋环流模式中。Zhang等 。通过将模式得 真实物理过程的各种近似处理，例如模式物理参数 到的 “伪观测数据”(包括余水位和近岸潮汐 ，余水 化过程；二是大气风应力和模式开边界等输入资料 位指总水位减去潮汐)同化到二维POM (Princeton 的不确定性；三是由数值计算方法引起的误差。解 OceanMode1)中，估计了风拖曳系数和侧边界的潮 决这一问题比较有效的方法之一是将风暴潮数值 汐等。Peng等 在台风风暴潮研究中应用四维变 模式和宝贵的观测资料通过最优的方法融合起来， 分同化技术 ，将 “伪观测数据”(由高分辨率模式产 发展具有同化能力的数值模式。近年来 ，各种同化 生)同化到模式中，对台风风暴潮进行了后报试验 ， 方法在风暴潮预报领域开展了较为广泛的应用n ， 结果表明只同化水位就可以有效 的改进预报结 其 中比较有代表性的是Kalman滤波在西北欧国家 果 。但遗憾 的是 ，超过同化 的3h窗 口，模式结果的 风暴潮业务化 中的应用[2-4]，和变分 同化在参数估 改进大概只能持续 5h左右 ，这也是 由台风风暴潮 计、状态估计中的应用 ]。 的特点决定的，其风场变化剧烈，增水速度快，持续 Kalman滤波方法的缺点是对观测资料序列有 时间短 ，一般在6h到 12h之问。本文则选取了温 严格限制，即观测资料不能出现缺测，并且要求资 带风暴潮作为试验对象，主要考虑其风场变化相对 料序列必须是等时间间隔 。本文所采用的四维变 缓慢 ，增减水的时间尺度更长 ，可 以持续 12h到 分同化算法突破了这些限制