最优插值方法在西北太平洋海温同化中的应用研究.pdfVIP

第25 卷 第2 期 海海海海 洋洋洋洋 预预预预 报报报报 V o l . 25 , N o . 2 2 0 0 8 年 5 月 MMMMAAAARRRRIIIINNNNEEEE FFFFOOOORRRREEEECCCCAAAASSSSTTTTSSSS May. 2 0 0 8 最优插值方法在西北太平洋海温同化中的应用研究 李 云，刘钦政，张建华，王 旭，李 燕 （国家海洋环境预报中心，北京100081） 摘 要：本文开发出一套与三维温盐流数值模式POM 相匹配的最优插值数值同化模块，具 有将船舶报资料和 Argo 海温数据加入数值模拟的能力。利用该同化系统，本文对 2001 ～ 2002 年进行了数值模拟试验，通过模拟结果和未加入同化的模拟结果的对比表明，该同化 系统能有效地改进三维海温模拟结果，使之更为接近于观测值。 关键词：资料同化；最优插值；船舶报资料；Argo 资料；三维海洋数值模式 中图分类号：P731 文献标识码：A 文章编号：1003 －0239 （2008）2 －0025 －08 引言 随着海洋观测数据的不断增加，在物理海洋的科学研究中，作为把空间上不规则分 布的离散观测资料融合到数值模式中的资料同化技术受到越来越多的重视，在海洋状态 数值模拟中也得到越来越多的应用[ 1 ～5] ，是海洋环境数值预报初试场建立的核心和关键。 最优插值方法是国际上应用较多也较为成熟的一种资料同化方法，它能考虑背景场 和观测误差的统计特征，即包含了观测、预报和分析之间的内在关系，决定出最小二乘 意义上的最佳权重，从而使得分析误差能在统计意义上达到最小[ 2 ] 。现今日本气象厅的 天气预报系统和海洋预报系统、美国NCEP 的ENSO 预报系统中的海洋资料同化系统中 都使用最优插值方法作为观测数据的同化方案之一。 本文结合国家海洋环境预报中心西北太平洋三维温盐流预报模式，利用船舶报资料 和Argo 海温资料进行温度同化，以检验最优插值方法在三维温盐流预报中所起的作用。 海洋模式 本文采用国家海洋环境预报中心在 POM 模式基础上所发展的西北太平洋温盐流数 值预报模式[6 ] ，模拟区域的东西范围从99°～150°E ，南北范围从2°～45°N ，包括了中国 海全海域、日本海的大部和西北太平洋的大部分；模式采用经纬度网格，水平分辨率为 15’×15’，垂直方向采用 坐标，分为不等距的15 层，最大水深限定为4000m ，使模式 最上层厚度不大于 10m ；时间积分采用内、外模分解算法，外模时间步长为 15s，内模 时间步长为600s；模式的温盐背景场采用SODA 全球0.5°×0.5°×40 层的三维温度盐度月 平均资料，利用双线性插值方法计算得到模式网格上的月平均气候温盐场；上边界强迫 收稿日期：2007-07-06 ；修订日期：2008-03-11 基金资助：国家基础研究重大项目前期研究专项（编号2001CCB00200） 作者简介：李 云（1981-），男，研究实习员，主要从事三维温盐流预报研究。 26 海 洋 预 报 25 卷 采用NCEP 一日4 次的再分析资料，采用块体公式计算海气界面通量和海面感热潜热通 量；侧边界采用辐射边界条件。 同化方案 温度同化方案采用最优插值同化方法，该方案基本算法为： 1 = + = + 其中 和 分别为分析场和背景场， 为最优权重， 为观测资料， 为背景场值 插值到观测位置的线性算子，上标 表示矩阵的转置， 为背景场误差协方差矩阵， 为观测场误差协方差矩阵。 为了充分利用有限的海洋观测资料，采用四维资料同化技术来研制和建立海洋资料 [7 ] 同化系统 。四维海洋资料同化技术的基本思想是利用海洋动力模式将观测到的海洋资 料进行时空外推，在给定的某一时刻和某一区域，即便是没有海洋观测资料，可以根据 其他区域以及过去时刻的观测资料，利用动