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

初始误差对双环流变异可预报性影响张坤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层浅水模式基金项目: 国家自然科学基金项目（41230420）；中国科学院知识创新工程重要方向项目（KZCX2-EW-201）；青岛市基础研究计划项目（11-1-4-95-jch）作者简介：张坤（1988-），男，山东济宁人，硕士，主要从事黑潮及其可预报性研究，电话E-mail: ; 王强，通信作者，助理研究员，E-mail: 。1前言双环流是副热带海盆风生环流的典型模态。黑潮和湾流延伸体均为双环流的重要组成部分，其流经海域是全球海气相互作用最强烈的区域之一1,2。研究发现黑潮和湾流延伸体存在拉伸（elongation）模态和收缩(contraction)模态间的年际变化3,4，这种年际变化对副热带海表温度和气候产生重要影响3-5。研究表明黑潮和湾流延伸体两种模态间的年际变化可以很好的被双环流射流变化表征6，因此，研究双环流变异对于理解黑潮和湾流延伸体低频变化有重要的指导意义。目前多数工作是从参数角度出发研究双环流的：Nauw等7研究了侧向摩擦对风生双环流射流不对称的影响；Moore8，Sura等9通过风应力研究了双环流变化；Primeau6,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。参数名称参数值参数名称参数值上层厚度尺度地球半径流速尺度地球自转角速度上层密度重力加速度风应力大小约化重力加速度侧向摩擦系数界面摩擦系数海盆东西方向长 海盆南北方向长表1 浅水模式中主要参数设置Table .1 The standard values of parameters in the shallow-water model研究表明动能可以较好区分双环流射流状态18,19：能量高时，双环流射流处于收缩模态；能量低时，射流处于拉伸模态。区域（0Kmx2500Km, 900Kmy1700Km）包含了双环流西边界以及射流与涡旋强烈作用的区域，因此将的动能作为区分双环流射流路径状态的指数，公式为： (6)图1 区域动能随时间变化序列Fig. 1 Time series of the kinetic energy in region 完成以上设置后，首先积分模式10年作为模式启动阶段，然后继续积分30年。图1为模式10年到40年区域内动能随时间变化序列，可以看出动能大小存在明显的年际变化，这段时间内出现了9次高低能量振荡。在约化重力模式中，海洋上层厚度等值线可以近似的表征表层流场15。因此为了直观了解高低动能对应的射流状态，图2（a）画出模式10年到40年内高动能对应的上层厚度平均场，（b）画出了低动能对应的上层厚度平均场，其中加粗黑线为海洋上层厚度600m等值线。可以看出，当动能高时，双环流射流处于收缩模态，当动能低时，双环流射流处于拉伸模态，这与Qiu等20的关于黑潮延伸体年际变化的研究结果是一致的。 图2 模式10年到40年内高低动能对应的平均上层厚度 单位：米 (a)图为高能量对应的上层厚度平均场，(b)图为低能量对应的上层厚度平均场, 加粗黑线为海洋上层厚度600m等值线，图中红色框为区域Fig. 2 Mean upper-layer thickness field under high(a)and low(b) kinetic energy Units: m3 CNOP方法简介及计算CNOP的基本设置条件非线性最优扰动（CNOP）方法是由Mu21等首次提出，主要用于寻找满足给定约束条件并在优化终止时刻得到最大发展的一类初始扰动。方程（1）离散后可以表示为： (7)和是零时刻和t时刻满足浅水模式（方程1）的状态向量。于是条件非线性最优扰动问题可描述为求解，使其满足以下方程： (8)其中，为初始误差，为目标函数。基于动能可以很好表征双环流射流状态的原因，定义目标函数为预报时刻区域内扰动动能，即： (9)方程（8）中，表示t时刻的误差。为初始扰动约束半径,为了估计所有可能的扰动对预报的影响，本文定义初始扰动约束范数为积分区域内初始扰动的总能量，即： (10)其中，为模式参数（表1），为约化重力加速度，本文中。背景场选取：本文主要是研究双环流变异过程的可预报性，选模式24.33年到25.25 25.08年的这次双环流变异过程为背景场，优化时间为240天。以24.3324.42年为第零天，图3给出了0天、90天、180天和240天（模式25.25年）时背景场的上层厚度分布。在t=0天时，双环流射流处于拉伸模态；当t=90天时，南北两支环流交汇处流向南偏，同时弯曲开始形成；当t=180天时，弯曲已经很明显；当t=240天时，弯曲程度达到最大。优化算法：采用Spectral Projected Gradient Version 2 (SPG2)22方法计算CNOP。图3 t时刻背景场的上层厚度 单位：米Fig .3 The upper-layer thickness of background field at time t Units: m4 初始误差对双环流变异可预报性影响4.1 CNOP和LCNOP在第3部分设置基础上，计算条件非线性最优初始扰动，发现除了存在CNOP（全局极大值点）外，目标函数在相空间内还存在一个局部极大值点，即LCNOP。CNOP和LCNOP均在约束条件的边界上，即。研究表明局部极大值也有比较明确的物理意义23，因此本文考虑CNOP和LCNOP对双环流变异可预报性影响。图4中阴影部分为初始误差上层厚度分量，等值线为背景场初始时刻的上层厚度场。图4（a）为CNOP的上层厚度分布，CNOP中有两个范围比较大的反气旋涡分布在海盆西侧射流的南北两侧，处于南侧的强度较大，此外，在两个反气旋中间存在一个范围小强度大的气旋涡；图4（b）为LCNOP的上层厚度分量分布，LCNOP中有两个气旋涡分布在海盆西侧射流南北两侧，处于南侧的强度较大，此外，在两个气旋涡之间以及下游区域存在强度小的反气旋涡。CNOP和LCNOP的上层厚度分量存在负的相似关系，下面计算二者的相似系数15，计算公式为： (11) 其中。经计算，相似系数为-0.85。图4 CNOP和LCNOP上层厚度分布 单位: 米Fig. 4 Upper-layer thickness component of CNOP and LCNOP Units: m.4.2 预报误差本部分着重研究两种类型初始误差对双环流变异预报的影响。我们可以把目标函数值当作衡量预报误差大小的一个标准，CNOP和LCNOP对应的目标函数值分为和，说明CNOP引起的预报误差大。为了考察这两类初始误差所导致的预报误差是否显著，分别将CNOP和LCNOP叠加到初始时刻背景场上(图3)，积分非线性模式240天，得到两个预报场。计算两个预报场相对于背景场区域内上层厚度的均方根误差(RMSE，图5)，可知均方根误差大，预报技巧低；均方根误差小，预报技巧高。图5中，0到120天内CNOP和LCNOP导致的预报误差大小基本相同，但从120天以后，CNOP导致的预报误差增长更快。在预报时刻t=240时，CNOP和LCNOP导致的预报误差分别为36.8m和28.6m。图5 区域内预报场与背景场上层厚度均方根误差（RMSE）随预报时间变化序列Fig .5 Time series of root-mean-square error (RMSE) of upper-layer thickness between forecast states and background states in area 图6画出了t=240天时海洋上层厚度600米等值线（在模式中可以近似代表射流路径）其中黑色线代表背景场，绿色线代表叠加CNOP后的预报场，红色线代表叠加LCNOP后的预报场。叠加CNOP后，在预报时刻射流弯曲处出现了涡的脱落，可知CNOP促进了射流弯曲的形成，弯曲在预报时刻前就达到最大，在预报时刻时涡已从弯曲处脱落; 相反，叠加了LCNOP后，在预报时刻射流弯曲程度相比背景场射流弯曲程度变小，说明LCNOP延缓了射流弯曲的形成。在双环流射流从拉伸模态到收缩模态的变化过程中，CNOP和LCNOP对射流弯曲形成的作用是相反的。图6 预报时刻t=240天时背景场的双环流射流路径（黑线），叠加CNOP后积分非线性模式得到的射流路径（绿线）,叠加LCNOP后积分非线性模式得到的射流路径（红线），其中双环流射流路径由海洋上层厚度600米等值线近似表示。Fig. 6 Double -gyre jet path at predicting time (t=240d) for background state(black line), obtained after perturbed by CNOP(green line) and obtained after perturbed by LCNOP(red line). The double-gyre jet is approximately represented by 600-m contour of the upper-layer thickness. 4.3误差发展及分析为清楚了解CNOP和LCNOP是如何影响双环流变化的，下面分析误差在预报时间0-240天的发展情况。图7画出了叠加CNOP后，积分非线性模式0、90、180和240天后误差的上层厚度发展（阴影部分），等值线为预报场的上层厚度分布。积分90天时，t=0时CNOP的两个反气旋涡在平流作用下相应移动，位于射流北侧的反气旋在x=1100Km,y=1400Km附近得到了比较大的发展，该位置靠近弯曲发生的区域，而初始时刻位于两个反气旋中间的气旋涡则移动到x=800Km,y=1200Km附近，强度得到很大发展；积分180天时，射流南侧的反气旋涡被耗散掉，相比积分90天时的情况，x=800Km,y=1200Km附近的气旋涡有较小增强，而x=1100Km,y=1400Km附近的反气旋涡则得到极大发展，该涡使得此处流场上层厚度增加。此外在x=1300Km,y=1400Km附近形成了一个强度大的气旋涡，并且使得此处流场上层厚度减小。此时在射流弯曲附近的反气旋涡和气旋涡的共同作用下，双环流射流弯曲程度进一步变大，随着积分时间增加，这对气旋和反气旋涡强度增大，弯曲不断变大。到积分t=240天时，双环流射流弯曲处已出现涡旋脱落。图7 预报时间t时CNOP的上层厚度发展（阴影部分）以及此时流场的上层厚度场（等值线） 单位：米Fig. 7 The upper-layer thickness component evolution of CNOP(shading) and the upper-layer thickness of the double-gyre(contour) at different predicting time t. Units: m.从误差发展角度了解到CNOP导致射流弯曲区域附近涡脱落，但是扰动为何在x=800Km,y=1200Km；x=1100Km,y=1400Km；x=1300Km,y=1400Km三个位置附近得到较大发展的原因仍不清楚。图8画出了叠加CNOP后240天内预报场的平均流场。图中区域A、B、C所分别对应x=1100Km,y=1400Km；x=800Km,y=1200Km ；x=1300Km,y=1400Km三个位置附近。区域A内因处在北面气旋涡和南面反气旋涡相互作用的区域，东西方向流速非常小，区域内有向北和向南的海流同时存在，流速切变很大；区域B内流速普遍较大，在其左侧海流流向是西北方向，而右侧流向基本为正西方向，甚至有向西南流的趋势，这个区域内流速方向有比较大的改变，流速切变也比较大；区域C内流速相对来说较小，但是同样存在较大的流向改变。A、B、C三个区域内流速切变都比较大，在这些区域很可能存在正压不稳定，使得误差不断地从背景场吸收能量，进而得到发展。图 8 叠加CNOP后预报时间段内平均流场 单位：米/秒Fig. 8 The mean velocity field of predicting period when perturbed by CNOP Units: m/s图9给出了叠加LCNOP后，积分非线性模式0、90、180和240天后误差的上层厚度发展（阴影部分），等值线为预报场的上层厚度场分布。积分90天时，t=0时LCNOP中的射流南北两侧气旋涡在平流作用下移动并得到发展，而且在x=800Km,y=1100Km附近误差发展成一个反气旋涡；积分180天时，位于射流南侧的气旋涡被耗散掉，而北侧气旋涡则在x=1000Km,y=1500Km附近得到发展，该涡使流场此处上层厚度减小，此外在x=1300Km,y=1400Km附近误差发展反气旋涡，该涡使得流场此处上层厚度增加,这对气涡和反气旋涡的共同作用阻碍了射流弯曲变大。随着积分时间增加，这对气旋和反气旋涡强度增大，到预报时刻t=240天时，预报场射流弯曲程度相比背景场（图3）变弱。图10画出了叠加LCNOP后240天内预报场的平均流场，图中用A，B，C三个矩形框分别代表x=1000Km,y=1500Km; x=800Km,y=1100Km和x=1300Km,y=1400Km附近的区域。同CNOP得到较大发展的三个区域相同，这三个区域内流速也存在较大切变，这些区域中很可能存在正压不稳定，使得误差不断吸收背景场能量，从而得到发展。图9 预报时间t时LCNOP的上层厚度发展（阴影部分）以及此时流场的上层厚度场（等值线） 单位：米Fig. 9 The upper-layer thickness component evolution of LCNOP (shading) and the upper-layer thickness of the double-gyre(contour) at different predicting time t. Units: m.图 10 叠加LCNOP后预报时间段内平均流场 单位：米/秒Fig. 10 The mean velocity field of predicting period when perturbed by LCNOP Units: m/s5 结论本文使用球坐标下1.5层浅水模式模拟海洋风生双环流，利用CNOP方法研究初始误差对双环流变异可预报性影响。通过计算条件非线性最优扰动，得到CNOP型初始误差和LCNOP型初始误差，两者相似系数为-0.85。在预报时刻t=240天时，相比LCNOP，CNOP导致的预报误差更大，导致双环流变异的预报技巧变的更低。通过研究误差发展，发现CNOP使得射流弯曲程度变大，并且在预报时刻射流弯曲处已出现涡的脱落；相反，LCNOP则使得射流弯曲程度变小。另外我们发现CNOP和LCNOP中误差得到较大发展区域都存在大的速度切变，很大可能是由于这些区域存在正压不稳定，使得误差能够不断吸收背景场能量，进而发展起来。本文研究给出了两类使得双环流变异预报技巧下降最大的初始误差：CNOP型和LCNOP型，并且探讨了其空间结构及发展，如果能够在数值预报中减少甚至消除这两种类型误差，将对黑潮和湾流延伸体的预报提供很大帮助。本研究仍存在一些不足：一方面本研究使用理想模型，并没有考虑地形、风场、海温等实际因素，并不能与实测情况进行对比验证；另一方面，模式采用1.5层浅水模式，不能够考虑斜压等因素。下一步计划采用多层浅水模式，进一步考虑正压、斜压等因素的影响，并且用物理过程更加详细的海洋模式进行模拟验证。致谢衷心感谢Henk A.Dijkstra and Arjen Terwisscha van Scheltinga等提供的球坐标下1.5层浅水模式以及良好的模式运行环境，这为本文工作提供了极大便利。参考文献：1 S. Dong and K. A. Kelly. Heat budget in the Gulf Stream region: The importance of heat storage and advectionJ.Journal of physical oceanography,2004, 34 (5): 1214-12312 K. A. Kelly. The relationship between oceanic heat transport and surface fluxes in the western North Pacific: 1970-2000J.Journal of climate,2004, 17 (3): 573-5883 K. A. Kelly, M. J. Caruso, S. Singh, et al. Observations of atmosphere-ocean coupling in midlatitude western boundary currentsJ.Journal of Geophysical Research-All Series-,1996, 101 6295-63124 B. Qiu. Interannual variability of the Kuroshio Extension system and its impact on the wintertime SST fieldJ.Journal of physical oceanography,2000, 30 (6): 1486-15025 H. Nakamura, G. Lin and T. Yamagata. Decadal climate variability in the North Pacific during the recent decadesJ.Bulletin of the American Meteorological Society,1997, 78 (10): 5-66 F. Primeau and D. Newman. Elongation and contraction of the western boundary current extension in a shallow-water model: A bifurcation analysisJ.Journal of Physical Oceanography,2008, 38 (7): 1469-14857 J. J. Nauw, H. A. Dijkstra and E. P. Chassignet. Frictionally induced asymmetries in wind-driven flowsJ.Journal of physical oceanography,2004, 34 (9): 2057-20728 A. M. Moore. Wind-induced variability of ocean gyresJ.Dynamics of atmospheres and oceans,1999, 29 (2): 335-3649 P. Sura and C. Penland. Sensitivity of a double-gyre ocean model to details of stochastic forcingJ.Ocean Modelling,2002, 4 (3): 327-34510 F. Primeau. Multiple equilibria and low-frequency variability of the wind-driven ocean circulationJ.Journal of physical oceanography,2002, 32 (8): 2236-225611 A. Mahadevan, J. Lu, S. Meacham, et al. The predictability of large-scale wind-driven flowsJ.Nonlinear Processes in Geophysics,2001, 8 (6): 449-46512 A. T. Van Scheltinga and H. Dijkstra. Conditional nonlinear optimal perturbations of the double-gyre ocean circulationJ.Nonlinear processes in geophysics,2008, 15 (5): 727-73413 R. X. Huang. Numerical simulation of wind-driven circulation in a subtropical/subpolar basinJ.Journal of physical oceanography,1986, 16 (10): 1636-165014 S. Shimokawa and T. Matsuura. The asymmetry of recirculation of a double gyre in a two layer oceanJ.Journal of oceanography,1999, 55 (3): 449-46215 Q. Wang, M. Mu and H. A. Dijkstra. Application of the conditional nonlinear optimal perturbation method to the predictability study of the Kuroshio large meanderJ.Advances in atmospheric sciences,2012, 29 (1): 118-13416 Q. Wang, M. Mu and H. A. Dijkstra. The similarity between optimal precursor and optimally growing initial error in prediction of Kuroshio large meander and its application to yargeted observationJ.Journal of Geophysical Research: Oceans,2013, 17 M. Mu and W. Duan. A new approach to studying ENSO predictability: Conditional nonlinear optimal perturbationJ.Chinese Science Bulletin,2003, 48 (10): 1045-104718 S. Jiang, F.-f. Jin and M. Ghil. Multiple equilibria, periodic, and aperiodic solutions in a wind-driven, double-gyre, shallow-water modelJ.Journal of physical oceanography,1995, 25 (5): 764-78619 S. Pierini. A Kuroshio Extension system model study: Decadal chaotic self-sustained oscillationsJ.Journal of physical oceanography,2006, 36 (8): 1605-162520 B. Qiu and S. Chen. Variability of the Kuroshio Extension jet, recirculation gyre, and mesoscale eddies on decadal time scalesJ.Journal of Physical Oceanography,2005, 35 (11): 2090-210321 M. Mu, W. Duan and B. Wang. Conditional nonlinear optimal perturbation and its applicationsJ.Nonlinear Processes in Geophysics,2003, 10 (6): 493-50122 E. G. Birgin, J. M. Martnez and M. Raydan. Nonmonotone spectral projected gradient methods on convex setsJ.SIAM Journal on Optimization,2000, 10 (4): 1196-121123 W. Duan, M. Mu and B. Wang. Conditional nonlinear optimal perturbations as the optimal precursors for El NinoSouthern Oscillation eventsJ.Journal of geophysical research,2004, 109 (D23): D23105 The impact of initial error on predictability of Double-gyre variabilityZhang Kun1,2, Mu Mu1,3, Wang Qiang1,3 (1. Institute of Oceanology, the Chinese Academy of Sciences, Qingdao 266071, China; 2. Graduate School, university of Chinese Academy of Sciences, Beijing 100049, China; 3. Key Laboratory of Ocean Circulation and Wave, the Chinese Academy of Sciences, Qingdao 266071, China)Key words: double-gyre variability; conditional nonlinear optimal perturbation(CNOP); shallow-water modelAbstract: A 1.5-layer reduced-gravi