初始误差对双环流变异可预报性影响张坤.docVIP

初始误差对双环流变异可预报性影响 张坤1,2，穆穆1,3，王强1,3 (1. 中国科学院 海洋研究所, 山东 青岛 266071; 2. 中国科学院大学 研究生院, 北京 100049; 3. 中国科学院海洋环流与波动重点实验室, 山东 青岛 266071) 摘要：使用球坐标下1.5层约化重力浅水模式模拟海洋风生双环流（double-gyre），模拟出的双环流射流存在拉伸和收缩两种模态间的年际变化，这与观测中黑潮和湾流延伸体的年际变化是类似的。取模式24.33年到25.25年的模拟结果为背景场，利用条件非线性最优扰动（CNOP）方法研究初始误差对双环流变异可预报性的影响，得到两种类型初始误差：全局CNOP和局部CNOP（LCNOP），两者相似系数为-0.85。CNOP和LCNOP不仅在结构上存在负的相似，它们对双环流变异的影响也是相反的。分别将CNOP和LCNOP叠加到初始时刻背景场上，考察误差随时间的发展，发现在射流从拉伸模态向收缩模态转变过程中，CNOP型初始误差使得射流弯曲程度变大，并在预报时刻已有涡的脱落；而LCNOP型初始误差则使得射流弯曲程度变小。相比LCNOP，CNOP引起的预报误差更大，导致双环流变异的预报技巧下降更多。两种类型误差得到较大发展的区域均有较大的流速切变，极有可能是因为存在正压不稳定，使误差能够不断从背景场吸收能量进而发展起来。 关键字：双环流变异；条件非线性最优扰动（CNOP）方法；1.5层浅水模式 1前言 双环流是副热带海盆风生环流的典型模态。黑潮和湾流延伸体均为双环流的重要组成部分，其流经海域是全球海气相互作用最强烈的区域之一[1,2]。研究发现黑潮和湾流延伸体存在拉伸（elongation）模态和收缩(contraction)模态间的年际变化[3,4]，这种年际变化对副热带海表温度和气候产生重要影响[3-5]。研究表明黑潮和湾流延伸体两种模态间的年际变化可以很好的被双环流射流变化表征[6]，因此，研究双环流变异对于理解黑潮和湾流延伸体低频变化有重要的指导意义。 目前多数工作是从参数角度出发研究双环流的：Nauw等[7]研究了侧向摩擦对风生双环流射流不对称的影响；Moore[8]，Sura等[9]通过风应力研究了双环流变化；Primeau[6,10]分别在准地转和浅水模式下通过粘性系数和风应力等参数考察了双环流的多平衡态和低频变化。目前从初始条件角度研究双环流的工作还比较少：在准地转模式下，Mahadevan等[11]通过线性奇异向量研究了双环流可预报性，Van Scheltinga等[12]通过条件非线性最优扰动（CNOP）方法研究了双环流稳定性。但是这些研究存在一些不足：准地转模式忽略了上层厚度变化及其伴随的非线性作用[13,14]；线性奇异向量采用的是线性近似的方法，在考察非线性作用方面有一定局限性。因此，本文利用浅水模式模拟双环流，采用CNOP方法研究初始误差对双环流变异可预报性的影响。一方面，浅水模式能够考虑上层厚度变化的影响；另一方面，CNOP方法能够克服线性近似方法的不足，能更好地考虑非线性过程对误差发展的作用，并且已经在黑潮大弯曲路径可预报性及最优前期征兆[15,16]、ENSO可预报性[17]、风生环流稳定性[12]等领域得到成功应用。文章主要安排如下：第2部分为模式介绍和模拟结果；第3部分为CNOP方法简介及计算CNOP的基本设置；第4部分为初始误差对双环流变异可预报性影响；第5部分为结论。 2 模式介绍及模拟结果 2.1 模式介绍 本文使用的数值模式是球坐标下1.5-层浅水方程模式，首先给出模式无量纲化后的控制方程： (1) 在无量纲化过程中，上层厚度平均深度、长度尺度、水平速度尺度、平流时间尺度以及风应力大小尺度分别记为、(地球半径)、、，。方程(1)中，表示纬向流速，表示经向流速，是海洋上层的厚度，和分别表示纬向和经向坐标，和分别为： ， (2) 其中和分别表示风应力的纬向分量和经向分量。物质导数可表示为： (3) 无量纲的模式参数的表达式为： ，, , , (4) 其中，，，，分别表示地球自转角速度、约化重力加速度、侧向摩擦系数和界面摩擦系数。 2.2参数设置及模拟结果 本文研究采用理想模型：模拟区域为矩形海盆，格点数为，边界为无滑动边界，主要参数设置参见表1。模式采用隐式差分方法，时间步长设为10天，风应力方程[6]为： (5) 其中,，为风应力大小控制参数，为风应力结构控制参数，本文中取0.1。 参数名称 参数