Argo资料的顺次三维变分同化方法研究.pdfVIP

Ａｒｇｏ资料的顺次三维变分同化方法研究 １，３ ２ １，３ ３ 何忠杰 　谢元富　李威 　韩桂军 ３ ３ ３ ３ 刘克修 　马继瑞 　李冬 　张学峰 １．中国海洋大学海洋环境学院，青岛 山东２６６００３ ２．ＮＯＡＡＲｅｓｅａｒｃｈ－ＥａｒｔｈＳｙｓｔｅｍＲｅｓｅａｒｃｈＬａｂｏｒａｔｏｒｙ，ＢｏｕｌｔｅｒＣＯ，ＵＳＡ ３．国家海洋局 海洋信息中心，天津３００１７１ 摘要：传统三维变分采用预先设定相关长度分布或递归滤波方法来设定背景场的 误差协方差矩阵，具有较大的任意性，且不能对多种尺度的海洋现象同时得到很 好的同化效果。本文引入顺次三维变分（Ｓ３ＤＶＡＲ）的方法，克服了以上不足之处。 顺次三维变分方法通过依次改变背景场的相关长度，逐次同化不同尺度的观测信 息，可以构造出非均匀的、各向异性的背景场误差协方差矩阵，具有更大的灵活性 和适应性。本文将该方法应用到全球Ａｒｇｏ资料的同化系统中，使同化效果得到明 显提高。 关键词：三维变分；同化方法；Ａｒｇｏ资料 １　引言 在三维变分同化过程中，背景场误差协方差矩阵Ｂ起着重要作用，它决定了分析场相 对于背景场的修正程度。而矩阵Ｂ的确定较为复杂，起初通常根据经验方法确定，即确定 ［１］ 某一海区的水平相关长度，再假定任意两点的误差协方差为距离的ｅ指数函数 ，从而可 以直接设定Ｂ矩阵的形式，这被称之为相关长度法。该方法具有很大的任意性，同时由于 Ｂ矩阵的存储和求逆需要很大的内存，相关长度法在应用中遇到了很多困难。 递归滤波（ＲｅｃｕｒｉｖｅＦｉｌｔｅ）方法克服了相关长度法的主要缺陷，可以保证Ｂ矩阵的对称 和正定，同时它具有计算速度快，节省内存等优点，因而自２０世纪８０年代以来，递归滤波逐  基金项目：国家自然科学基金项目４０４７６００６资助。 作者简介：何忠杰（１９７７～），男，助理研究员，从事海洋数据同化、卫星遥感等研究。 Ｅ－ｍａｉｌ：ｈｅｚｈｏｎｇｊｉｅ＿７７＠ｙａｈｏｏ．ｃｏｍ．ｃｎ １８３ 渐成为三维变分数据同化过程中构造背景场误差协方差矩阵的一种重要手段［２～４］。 Ａｒｇｏ浮标资料以其分布范围广、数据量大、精度高的优势，逐渐成为海洋温度观测信息 的重要来源，尤其在海温的数值再分析方面已得到广泛应用。由于Ａｒｇｏ资料的空间分布密 度较大，其中包含有多种尺度的温度变化信息，而递归滤波法和相关长度法只对特定波长 的信息敏感，无法同时对多种尺度的观测信息得到较好的同化效果，因此采用这些同化方 法对Ａｒｇｏ资料进行同化时会丢失许多信息，相关长度法在这方面的问题则更突出。如何更 加灵活地构造误差协方差矩阵，使之能够准确地反映背景场的误差分布，成为近来三维变 分数据同化方法应用研究的一个重要问题［５～７］。 ［８］ Ｘｉｅ等 首次提出用顺次三维变分法依次对不同尺度的观测信息进行同化，解决了相 关长度法和递归滤波法等传统三维变分方法无法同时提取多尺度信息的问题，本文将这种 方法应用于Ａｒｇｏ温度数据的三维同化系统中。 ２　顺次三维变分的基本原理 在三维变分同化系统中，要求背景场的误差协方差矩阵Ｂ是正定、对称的，但直接确定 －１ ［９］ Ｂ比较困难，而直接求解Ｂ 更加困难。Ｈｕａｎｇ 证明了递归滤波器可以近似地起到算符 １／２ Ｂ 的作用，于是利用递归滤波方法构造背景场误差协方差矩阵，可以方便的计算出目标泛 函。 递归滤波器是一种低通滤波器，以一维的情形为例，递归滤波操作可表示为： ∞