Marine Hydrology

Spatial-temporal variations of the dynamic parameters of internal solitary waves in the Sulu-Celebes Sea

  • XIE Jieshuo , 1, 2 ,
  • GONG Yankun 1, 2 ,
  • NIU Jianwei 1, 2 ,
  • HE Yinghui 1, 2 ,
  • CHEN Zhiwu 1, 2 ,
  • XU Jiexin 1, 2 ,
  • CAI Shuqu , 1, 2
Expand
  • 1. State Key Laboratory of Tropical Oceanography, South China Sea Institute of Oceanology, Chinese Academy of Sciences, Guangzhou 510301, China
  • 2. Southern Marine Science and Engineering Guangdong Laboratory (Guangzhou), Guangzhou, 511458, China
CAI Shuqun. email:

Copy editor: LIN Qiang

Received date: 2021-12-24

  Revised date: 2022-03-17

  Online published: 2022-04-19

Supported by

Project of National Natural Science Foundation of China(42130404)

Project of National Natural Science Foundation of China(41776008)

Key Research Program of Frontier Sciences, CAS(QYZDJ-SSW-DQC034)

Southern Marine Science and Engineering Guangdong Laboratory (Guangzhou)(GML2019ZD0304)

Special project by Department of natural resources of Guangdong Province([2020]017)

Youth Science and Technology Innovation Talent of Guangdong TeZhi Plan(2019TQ05H519)

Natural Science Foundation of Guangdong Province(2020A1515010495)

Natural Science Foundation of Guangdong Province(2021A1515012538)

Natural Science Foundation of Guangdong Province(2021A1515011613)

Abstract

On the basis of the weakly nonlinear theory of internal solitary waves (ISWs) and the reanalyzed assimilation data in the Sulu-Celebes Sea in both summer and winter, we calculate the dynamic parameters of the ISWs, including the linear wave speed, the first- and second-order nonlinear parameters and the linear dispersion parameter, and investigate the spatial-temporal variations of these parameters. We find that, although the intrusion from the north Pacific Ocean into the Celebes Sea is very obvious, the spatial-temporal variations of the dynamic parameters of the ISWs in the Sulu Sea is more significant than that in the Celebes Sea. In the Sulu sea, the linear speed of ISW in summer is ~0.1 m·s-1 larger than that in winter; however, in the Celebes Sea, the linear speed of ISW in summer is ~0.05 m·s-1 smaller than that in winter. Regardless of the first- and second-order nonlinear parameters, their variations in the Celebes Sea are relatively weak, but in the Sulu Sea, their variations are significant. In the Sulu Sea, the first-order nonlinear parameter in summer is ~3×10-3 s-1 higher than that in winter, but the second-order nonlinear parameter in summer is ~3×10-5 m-1·s-1 lower than that in winter. In addition, in both the Sulu and Celebes Sea, the linear dispersion parameter in summer tends to reduce, but the reduction magnitude in Sulu Sea is relatively large. This paper shows that the variation in the depth of the maximum buoyancy frequency in the Sulu-Celebes Sea is the dominant factor that in turn results in the associated spatial-temporal variations of the dynamic parameters of ISWs.

Cite this article

XIE Jieshuo , GONG Yankun , NIU Jianwei , HE Yinghui , CHEN Zhiwu , XU Jiexin , CAI Shuqu . Spatial-temporal variations of the dynamic parameters of internal solitary waves in the Sulu-Celebes Sea[J]. Journal of Tropical Oceanography, 2022 , 41(6) : 132 -142 . DOI: 10.11978/2021183

海洋内孤立波是发生在密度跃层附近的一种强非线性高频波动(Apel et al, 2006; Zhang et al, 2011; Cai et al, 2012; Guo et al, 2014)。内孤立波致近表层辐聚和辐散可在遥感图像上留下清晰可辨的印迹(Wang et al, 2013, 2019; Li et al, 2013; Jia et al, 2018; 孙丽娜 等, 2018; 李子木 等, 2014), 基于该印迹, 海洋学家汇总了全球海洋内孤立波的图集(http://www.internalwaveatlas.com/)。大量遥感相关观测研究表明, 苏禄—苏拉威西海附近海域是海洋内孤立波的活跃区域之一(Jackson et al, 2011; Zhang et al, 2020; 张涛 等, 2020; Gong et al, 2022)。
苏禄海和苏拉威西海位于印度尼西亚与菲律宾之间, 一方面二者被苏禄群岛等诸岛所间隔, 另一方面二者又经锡布图通道产生联系。图1a显示了苏禄海一个典型的在苏禄群岛附近激发的北—西北向传播的内孤立波, 其传播范围可波及大部分苏禄海海盆(张涛 等, 2020)。图1b显示了苏拉威西海一个典型的在锡布图岛附近激发的南—东南向传播的内孤立波。此外, 苏拉威西海还存在大量受苏拉威西海东南部海脊激发的西向传播内孤立波(Zhang et al, 2020)。两个源地激发的内孤立波亦可波及广阔的苏拉威西海盆。苏禄海和苏拉威西海是沟通南海与印尼海的两大主要通道之一(Han et al, 2009; Qu et al, 2009; Du et al, 2010; Gordon et al, 2012; Cheng et al, 2016), 而内波致水体混合可改变该海域局地水团性质(Hurlburt et al, 2011; Nagai et al, 2015), 进而对南海环流格局及印尼贯穿流水体输运等产生潜在影响(Cai et al, 2010; Li et al, 2019)。因此, 深入开展该海域内孤立波研究, 不仅对保障该海域航行安全具有实际意义, 也对南海环流及印尼贯穿流的研究具有重要科学意义。
图1 苏禄—苏拉威西海内孤立波传播影像

红色虚线箭头表示内孤立波传播方向
a. 苏禄海2020年4月11日MODIS影像; b. 苏拉威西海2003年3月21日MODIS影像。审图号为GS(2016)1609

Fig. 1 Images of internal solitary waves in the Sulu-Celebes Sea

相比于南海等热点研究区域(Alford et al, 2015; Xie et al, 2015; Li et al, 2016; Lai et al, 2019; Bai et al, 2019; Xu et al, 2021), 该海域内孤立波研究无论是在观测(Tessler et al, 2012; Liu et al, 2019; Zhao et al, 2021)或模式(Hurlburt et al, 2011; Nagai et al, 2015)方面均相对偏少。尽管如此, Apel等(1985)于1980年开展的现场观测仍在全球海洋内孤立波的早期研究中占有一席之地, 这次观测研究也使得Apel等推测出苏禄海内孤立波的生成机制及其主要源地。2009年, 美国学者Arnold Gordon等在开展菲律宾群岛区域过海槛的溢流水输运及水动力过程现场调查时, 获得了经苏禄群岛激发的内孤立波初期演化特征(Tessler et al, 2012), 这是目前关于苏禄海内孤立波较为近期的一项现场调查研究。基于此次观测, Tessler等比较了有限深理论(Joseph, 1977)与浅水KdV理论(Korteweg et al, 1895)的适应性, 指出二者在刻画苏禄海内孤立波特征方面各有优劣, 并且总体而言, KdV理论模型刻画的内孤立波特征及能量甚至更为接近观测结果。观测研究还显示, 即使在我国南海北部水深约3000多米的海盆处, 弱非线性的浅水KdV理论模型也被发现能够较好地刻画海洋内孤立波的主要特征及相关过程(Klymak et al, 2006; Grimshaw et al, 2014; Huang et al, 2016; Chen et al, 2019)。此外, 大量实验室室内试验也一致性地显示, 即使在受到其模型特征标度律(model scaling regimes)限制的情形下, 浅水KdV理论依然比有限深理论模型具有更好的拟合特性(Koop et al, 1981; Segur et al, 1982; Kao et al, 1985)。实际上, 基于KdV理论的动力参数是描述内孤立波特征的基本运动学参量。这些参数被广泛应用于世界大洋(Grimshaw et al, 2007)及诸多区域海, 包括我国南海(Cai et al, 2014; Liao et al, 2014)及东海(Cho et al, 2016)等海域内孤立波时空变化特征的研究。
研究还表明, 苏禄—苏拉威西海域因受季风及南海与太平洋水的共同调制作用, 环流格局表现出明显的时空差异性(Cai et al, 2009; He et al, 2017; Hao et al, 2021)。该海域环流结构的时空差异调制着该海域内孤立波传播演变的背景环境, 从而可能对该海域内孤立波动力参数的时空变化产生调制作用。然而, 现阶段依然缺乏该海域内孤立波动力参数时空变化特征的研究。鉴于KdV理论仍然在该海域具有较好的适应性, 我们在本文中拟采用eKdV理论模型, 基于HYCOM(hybrid coordinate ocean model)再分析同化产品, 计算苏禄—苏拉威西海的内孤立波动力参数, 并研究其时空变化特征。

1 模型、数据及背景介绍

1.1 理论模型及其动力参数

KdV理论的基本思想是对控制内波特征的小参数进行渐进展开, 忽略高阶非线性项, 并保留低阶非线性项(如Choi et al, 1999; Grimshaw et al, 2004; Yuan et al, 2018)。这里, 我们采用保留至二阶非线性项的eKdV理论模型描述该海域内孤立波动力参数的时空变化特征(Apel et al, 2006; Xie et al, 2016), 模型控制方程写为
$\frac{\partial \eta }{\partial T}+C\frac{\partial \eta }{\partial X}+\alpha \eta \frac{\partial \eta }{\partial X}+\beta \frac{{{\partial }^{3}}\eta }{\partial {{X}^{3}}}+{{\alpha }_{1}}\eta \frac{\partial \eta }{\partial X}=0$
其中, TX分别为时、空变量; η为内孤立波扰动位移; c为内孤立波的线性长波速度, 可通过求解如下特征值问题获得。
${{c}^{2}}\frac{{{\operatorname{d}}^{2}}\varphi }{\operatorname{d}{{z}^{2}}}+N_{\text{B}}^{2}\varphi =0\text{ }$
该式在z=0, -H时, 满足φ=0; H代表水深。式(2)中NB代表背景层化, φ为内孤立波线性长波扰动特征分布函数。尽管除线性速度外, 内孤立波速度还受到其非线性部分的影响, 但线性部分仍然对内孤立波传播速度的影响起主导作用。实际上, 即使对于大振幅强非线性内孤立波, 其线性部分依然可达80%~90%左右。因此, 详细评估内孤立波线性速度的时空变化对于内孤立波准确到达时间的研究仍然具有明显的实际意义(Xie et al, 2021)。
此外, 式(1)中的a为内孤立波的一阶非线性参数, 可表述为
$\alpha =\frac{3c}{2}\frac{\int_{-H}^{0}{{{({\operatorname{d}\varphi }/{\operatorname{d}z}\;)}^{3}}\operatorname{d}z}}{\int_{-H}^{0}{{{({\operatorname{d}\varphi }/{\operatorname{d}z}\;)}^{2}}\operatorname{d}z}}$
α<0时, 内孤立波为凹陷波; 反之, 当α>0时, 内孤立波为抬升波; 在α=0时, 内孤立波发生由凹陷波到抬升波的极性转变。一阶非线性参数α绝对数值的大小代表内孤立波的非线性强度, 其所在的一阶非线性项与色散项的平衡关系对海洋内孤立波传播演化特征起主要的调制作用。式(1)中的β即为色散参数, 可写为
$\beta =\frac{c}{2}\frac{\int_{-H}^{0}{{{\varphi }^{2}}\operatorname{d}z}}{\int_{-H}^{0}{{{({\operatorname{d}\varphi }/{dz}\;)}^{2}}\operatorname{d}z}}$
色散系数β数值的大小代表着对内孤立波色散分裂抑制的强度。当线性色散项占优时(也即$\frac{\alpha }{\beta }{{\eta }_{s}}L_{s}^{2}<1$时), 内孤立波色散分裂会被抑制, 反之, 当一阶非线性项占优时(也即$\frac{\alpha }{\beta }{{\eta }_{s}}L_{s}^{2}<1$时), 内孤立波更容易形成并发生分裂(注: 此处ηsLs分别代表特征扰动的垂向振幅及水平尺度; Xie et al, 2015)。然而, 仅考虑一阶非线性项依然难以满足对海洋大振幅内孤立波强非线性演化特征的刻画, 诸多研究指出考虑二阶非线性的eKdV理论在描述大振幅强非线性内孤立波方面具有更好的适应性(Grimshaw et al, 2004; Helfrich et al, 2006), 式(1)中的α1即为二阶非线性参数, 可写为
${{\alpha }_{1}}=\frac{\int_{-H}^{0}{\begin{align} & \left\{ 3{{c}^{2}} \right.\left[ 3\operatorname{d}\xi /\operatorname{d}z-2{{({\operatorname{d}\varphi }/{\operatorname{d}z}\;)}^{2}} \right]{{({\operatorname{d}\varphi }/{\operatorname{d}z}\;)}^{2}}-{{\alpha }^{2}}{{({\operatorname{d}\varphi }/{\operatorname{d}z}\;)}^{2}} \\ & +\alpha c({\operatorname{d}\varphi }/{\operatorname{d}z}\;)\left. \left[ \text{5}{{({\operatorname{d}\varphi }/{\operatorname{d}z}\;)}^{2}}-4\operatorname{d}\xi /\operatorname{d}z \right] \right\}\operatorname{d}z \\ \end{align}}}{2c\int_{-H}^{0}{{{({\operatorname{d}\varphi }/{\operatorname{d}z}\;)}^{2}}\operatorname{d}z}}$
式中$\xi $表示对内孤立波模态特征函数的非线性修正, 其可通过如下非齐次特征值问题获得
$\frac{{{\operatorname{d}}^{2}}\xi }{\operatorname{d}{{z}^{2}}}+\frac{N_{B}^{2}}{{{c}^{2}}}\xi =\left( \text{3}\frac{\operatorname{d}\varphi }{\operatorname{d}z}-\frac{\alpha }{c} \right)\frac{{{\operatorname{d}}^{2}}\varphi }{\operatorname{d}{{z}^{2}}}$
该式在z=0, -H, zmax时, 满足ξ=0, 其中zmax使得φ (zmax)=1。

1.2 苏禄—苏拉威西海表层环流

为研究苏禄—苏拉威西海内孤立波动力参数的时空变化, 我们采用HYCOM最新的3.1版1/12°高分辨率同化数据(global ocean forecasting system 3.1)。前人研究表明, 苏禄—苏拉威西海背景环流表现出明显的季节性差异。例如, 苏禄海在冬季表现为近似气旋式的海盆尺度环流, 而在夏季其表现为明显的反气旋式海盆尺度环流(Cai et al, 2009; Han et al, 2009; He et al, 2017); 并且, 经验正交函数(EOF)分析结果表明(Cai et al, 2009), 冬季1月份的气旋式环流特征最为显著, 而夏季8月份的反气旋式环流特征最为显著。两种不同的环流格局, 会引起背景水体层化等的时空差异, 从而改变内孤立波相关的动力参数。因此, 我们分别研究冬季及夏季苏禄—苏拉威西海内孤立波动力参数的变化, 对HYCOM数据进行了月平均处理。图2a展示苏禄—苏拉威西海冬季1月份海表环流特征, 图2b展示该海域夏季8月份海表环流特征。可见, 与苏禄海环流格局相反, 冬季1月份苏拉威西海主要近似表现为反气旋式环流, 夏季8月份则主要表现为气旋式环流。而且, 由于苏拉威西海同时受到较强的西北太平洋水入侵影响, 其海表环流强度总体比苏禄海大, 特别在8月份较为明显(图2b)。
图2 苏禄—苏拉威西海海表环流特征

a. 1月; b. 8月。审图号为GS(2016)1609

Fig. 2 Characteristics of surface circulation in the Sulu-Celebes Sea

2 苏禄—苏拉威西海水体层化时空变化特征

浮力频率是反映海洋水体层化及内孤立波动力参数变化的主要物理量。采用上述月平均的1月及8月份HYCOM同化数据, 我们可以计算得到高分辨率的苏禄—苏拉威西海的浮力频率。通过展示最大浮力频率(图3a、b)及其所在深度(图3c、d), 可以获得苏禄—苏拉威西海水体层化特征的时空差异性。图3a、b表明, 总体而言, 无论夏冬, 苏禄海海盆中水体层化的最大浮力频率(约0.015~0.017s-1)均稍大于苏拉威西海海盆的最大浮力频率(约0.13~0.15s-1)。图3c表明, 冬季1月份苏禄海盆的最大浮力频率所在深度(约70~90m)总体上比苏拉威西海海盆的最大浮力频率所在深度(约80~90 m)稍浅。然而, 图3d则表明, 夏季8月份苏禄海盆的最大浮力频率所在深度(约100~120 m)总体上却远深于苏拉威西海海盆的最大浮力频率所在深度(约90~100m)。
图3 苏禄—苏拉威西海水体层化特征

a. 1月最大浮力频率分布; b. 8月最大浮力频率分布; c. 1月最大浮力频率所在深度分布; d. 8月最大浮力频率所在深度分布。审图号为GS(2016)1609

Fig. 3 Characteristics of sea water stratification in the Sulu-Celebes Sea

进一步, 通过用夏季8月份的最大浮力频率及其所在深度分别减去冬季1月份的最大浮力频率及其所在深度, 给出苏禄—苏拉威西海水体层化的季节性差异(图4a、b)。可以看到, 总体而言, 无论是苏禄海还是苏拉威西海, 最大浮力频率的季节性变化幅度(基本小于0.001s-1)在两个海盆中的变化均不明显(见图4a等值线)。尽管如此, 图4b显示, 在苏禄海, 最大浮力频率所在深度的季节性差异较为明显, 最大可达40~50m, 但是在苏拉威西海, 除西北角内孤立波活动相对较为稀疏的小部分海域外, 最大浮力频率所在深度的季节性差异总体较小, 仅为10m左右。
图4 苏禄—苏拉威西海水体层化季节差异

a. 最大浮力频率差异; b. 最大浮力频率所在深度差异。审图号为GS(2016)1609

Fig. 4 Seasonal difference of sea water stratification in the Sulu-Celebes Sea

3 苏禄—苏拉威西海内孤立波动力参数时空变化

图5—8分别展示了基于公式(2)—(5)计算的1月和8月份苏禄—苏拉威西海内波线性速度c、一阶非线性动力参数α、色散动力参数β及二阶非线性参数α1的空间分布。在图9中, 进一步通过用夏季8月份减去冬季1月份的动力参数, 分别获得了这些动力参数在冬季和夏季的季节性差异。
图5 苏禄—苏拉威西海内孤立波线性速度c分布

a. 1月; b. 8月。审图号为GS(2016)1609

Fig. 5 Linear speed of internal solitary waves in the Sulu-Celebes Sea

图6 苏禄—苏拉威西海内孤立波一阶非线性参数α分布

a. 1月; b. 8月。审图号为GS(2016)1609

Fig. 6 First-order nonlinear parameter α of internal solitary waves in the Sulu-Celebes Sea

图7 苏禄—苏拉威西海内孤立波色散参数β分布

a. 1月; b. 8月。审图号为GS(2016)1609

Fig. 7 Linear dispersion parameter β of internal solitary waves in the Sulu-Celebes Sea

图8 苏禄—苏拉威西海内孤立波二阶非线性参数α1分布

a. 1月; b. 8月。审图号为GS(2016)1609

Fig. 8 Second-order nonlinear parameter α1 of internal solitary waves in the Sulu-Celebes Sea

图9 苏禄—苏拉威西海内孤立波动力参数季节差异

a. 线性速度c差异; b. 一阶非线性参数α差异; c. 色散参数β差异; d. 二阶非线性参数β差异。审图号为GS(2016)1609

Fig. 9 Seasonal difference of the dynamic parameters of internal solitary waves in the Sulu-Celebes Sea

图5表明, 无论是冬季1月或夏季8月, 由于苏拉威西海海盆总体水深相对较深, 其线性速度总体约达2.7~2.8m·s-1, 然而, 由于苏禄海海盆总体水深相对较浅, 其线性速度仅约1.7~2.1m·s-1。因此, 如图1a显示, 由于苏禄海海盆总体相对较小的线性传播速度, 相邻内孤立波波包之间的距离仅约81~108km, 但图1b显示, 由于苏拉威西海总体相对较大的线性速度, 相邻内孤立波波包之间的距离则可达122~137km。尽管两个海域内孤立波线性速度的空间差异主要是受到海水深度影响, 但根据特征方程(2)可知, 两个海域内孤立波线性速度的季节性差异则主要是由背景环流格局变化造成的水体层化差异所导致(图34)。例如, 图9a表明, 夏季8月份的苏禄海内波线性速度总体上比冬季1月份约大0.1m·s-1; 与此相反, 在苏拉威西海, 夏季8月份内波线性速度总体上却比冬季1月份约小0.05m·s-1
此外, 空间上, 图5表明, 苏禄海海盆南侧(~8°N以南)线性速度总体比北侧(~8°N以北)偏大约0.2~0.3m·s-1, 因此, 从图1a的内孤立波遥感图像可以看出, 第一、二个内孤立波波包之间的距离相较于后续波包之间的距离偏大约20km。然而, 整个苏拉威西海盆线性速度的变化则不大, 但图1b苏拉威西海盆遥感图像则显示第一、二个波包之间的距离相较于第二、三个之间的距离偏小约15km, 通过观察各波包中孤立子个数可以推测, 这主要是因为后续波包中内孤立波非线性部分速度的显著增加所导致。
图6结果表明, 总体而言, 苏禄海海盆区域内孤立波一阶非线性动力参数α可达1.0×10-2~ 1.3×10-2s-1, 然而, 在苏拉威西海, 其数值相对较小, 仅约0.8×10-2s-1。因此, 相对而言, 苏禄海区域内波更容易获得较强的非线性, 因此, 该区域的内波相对更容易发生分裂并产生非线性内孤立波。图9b中关于一阶非线性动力参数α的季节性差异则表明, 苏拉威西海的一阶非线性动力参数α基本不存在季节性的差异, 但是在苏禄海, 夏季8月份的一阶非线性动力参数α比冬季1月份高出约3×10-3s-1量级。该季节变化显示, 8月份的苏禄海区域相对可能更容易有非线性内孤立波的裂变发生。
图7结果表明, 无论是冬季1月或夏季8月, 苏拉威西海海域的内孤立波色散动力参数β均可达1.0×106~1.4×106m3·s-1, 但在苏禄海, 内孤立波色散动力参数β仅约2.0×105~6.0×105m3·s-1, 远小于苏拉威西海海域。这表明, 在苏拉威西海, 内孤立波的色散分裂会被显著抑制, 不利于内孤立波的形成, 但在苏禄海, 色散动力对内孤立波分裂及形成的抑制性则较低, 因此更有利于内孤立波的形成。例如, 图1中的MODIS (moderate resolution imaging spectroradiometer)遥感图像显示, 苏禄海海域一个内孤立波波包中分裂形成的内孤立子个数可达5~10个左右, 但苏拉威西海海域的一个内孤立波波包中内孤立子的个数至多仅约2~4个。图6图7中一阶非线性与色散动力的平衡还表明, 苏禄海内波总体而言相对更容易分裂并形成高频非线性内孤立波。图9c关于色散动力参数β的季节性差异进一步表明, 夏季8月份苏禄海和苏拉威西海的色散动力参数均减少, 但苏禄海减少幅度更大, 可达1×105m3·s-1, 而苏拉威西海的减弱幅度则相对较小, 约0.6×105m3·s-1。因此, 我们发现一个有意思的现象是, 虽然相比于苏禄海, 苏拉威西海域受到西北太平洋水的入侵更为显著, 但苏禄海域内孤立波动力参数却比苏拉威西海具有更显著的季节性变化。
类似于图6中一阶非线性动力参数α, 图8中关于二阶非线性动力参数α1的结果表明, 苏禄海海盆区域内孤立波二阶非线性动力参数强度(0.8×10-4~0.9×10-4m-1·s-1)同样也表现为强于其在苏拉威西海的变化(0.3×10-4~0.4×10-4m-1·s-1)。图9d的季节性差异则进一步表明, 苏拉威西海的二阶非线性动力参数基本无季节性差异, 但是在苏禄海其季节性变化则比较明显。具体而言, 苏禄海夏季8月份的二阶非线性动力参数可比1月份低约0.3×10-4m-1·s-1。并且, 对比图68表明, 无论是一阶或二阶非线性动力参数, 在苏拉威西海海盆的空间变化均非常微弱, 但在苏禄海具有较为显著的空间变化。例如, 1月份苏禄海的一阶非线性动力参数自南向北逐渐增强, 通过与图3c对比发现, 这主要是由于1月份最大浮力频率所在深度的空间变化导致。
空间上, 苏禄海盆的两个非线性动力参数均比苏拉威西海盆具有更大的变化幅度。苏禄海盆的非线性参数在其东南部及南部海盆区域较大。苏禄海盆的色散参数也具有类似的空间变化特点。但在苏拉威西海盆, 其两个非线性参数的空间变化均较为微弱, 仅色散参数的空间变化较为显著, 主要表现为在其南部的部分区域显著增大。
苏禄海各动力参数较大的季节性差异主要也受控于苏禄海最大浮力频率所在深度的季节性差异。例如图4b表明, 苏禄海最大浮力频率所在深度的季节性变化可达40m以上, 也即苏禄海最大浮力频率深度季节性变化幅度约达40%, 但苏拉威西海最大浮力频率所在深度的季节性差异则仅10m左右, 也即变化幅度仅约10%。此外, 图4a显示, 苏禄—苏拉威西海域最大浮力频率的季节及空间变化总体上均不超过0.001s-1, 即时空变化幅度低于8%~10%。因此, 苏禄—苏拉威西海环流引起的水体层化最大浮力频率所在深度的时空变化可能是造成上述内孤立波动力参数时空变化的根本原因。
总体而言, 上述动力参数时空变化的结果表明, 相比于苏拉威西海, 苏禄海内孤立波动力参数较容易受到海盆尺度环流变化的影响。因此, 我们在开展苏禄—苏拉威西海海域内孤立波特征及动力参数变化研究时, 有必要对局地海盆尺度环流变化的调制作用进行细致分析, 特别是在苏禄海域, 海盆尺度环流的调制作用更加显著且不容忽视。

4 总结与讨论

苏禄—苏拉威西海是海洋内孤立波活跃区域之一。大量前人研究表明, 苏禄—苏拉威西海域因受季风及南海与太平洋水的共同调制作用, 环流格局表现出明显的时空差异性。然而, 关于该海域环流格局时空差异性对局地内孤立波动力参数时空变化影响的研究却较少。文章基于弱非线性假设的eKdV理论模型, 采用HYCOM再分析同化数据, 计算了苏禄—苏拉威西海冬季1月及夏季8月的内孤立波动力参数, 包括内波线性速度c、一阶非线性动力参数α、色散动力参数β及二阶非线性参数α1, 并研究了其时空变化特征。
基于HYCOM再分析数据, 我们对比了苏禄海和苏拉威西海的环流格局, 可以发现两个海域环流格局正好相反。冬季1月份, 苏禄海主要表现为气旋式环流, 而苏拉威西海则主要表现为反气旋式环流, 但是在夏季8月份, 苏禄海表现为反气旋式环流, 而苏拉威西海则表现为气旋式环流。由于环流格局的差异, 苏禄—苏拉威西海水体层化也存在着时空差异性。苏禄海海盆中水体层化的最大浮力频率及其所在深度均大于苏拉威西海海盆水体的最大浮力频率及其所在深度, 并且, 苏禄海水体最大浮力频率所在深度存在极为显著的季节性差异, 其夏季最大浮力频率所在深度可达120m, 比冬季1月份的最大浮力频率深度深约40~50m。
虽然相比于苏禄海, 苏拉威西海域受到西北太平洋水的入侵更为显著, 但我们发现, 苏禄海内孤立波动力参数比苏拉威西海具有更为显著的时空变化特征。夏季苏禄海内孤立波线性速度总体上比冬季约大0.1m·s-1; 与此相反, 夏季苏拉威西海内孤立波线性速度却比冬季1月份小约0.05m·s-1。无论是一阶或二阶非线性动力参数, 其在苏拉威西海海盆的时空变化均非常微弱, 但在苏禄海则具有较为显著的时空变化。例如, 苏禄海夏季一阶非线性动力参数比冬季高出约3×10-3s-1, 但是夏季二阶非线性动力参数却比1月份低约3×10-5m-1·s-1。虽然夏季苏禄海和苏拉威西海的色散动力参数均减少, 但苏禄海减少幅度更大。此外, 除季节差异外, 相比于苏拉威西海, 苏禄海内孤立波动力参数也表现出更为显著的空间差异。这些结果表明, 在开展苏禄—苏拉威西海内孤立波特征及动力参数变化研究时, 有必要对局地海盆尺度环流变化的调制作用进行细致分析, 特别是在变化显著的苏禄海, 海盆尺度环流的调制作用尤为不容忽视。
[1]
李子木, 蔡树群, 陈举, 等, 2014. 2010-2011年吕宋海峡西侧潜标观测的初步分析[J]. 热带海洋学报, 33(1): 10-16.

LI ZIMU, CAI SHUQUN, CHEN JU, et al, 2014. Preliminary analysis of observations by deep submersible mooring in west Luzon Strait during 2010 to 2011[J]. Journal of Tropical Oceanography, 33(1): 10-16. (in Chinese with English abstract)

DOI

[2]
孙丽娜, 张杰, 孟俊敏, 2018. 基于遥感与现场观测数据的南海北部内波传播速度[J]. 海洋与湖沼, 49(3): 471-480.

SUN LINA, ZHANG JIE, MENG JUNMIN, 2018. On propagation velocity of internal solitary waves in the northern South China Sea with remote sensing and in-situ observations data[J]. Oceanologia et Limnologia Sinica, 49(3): 471-480. (in Chinese with English abstract)

[3]
张涛, 张旭东, 2020. 基于MODIS和VIIRS遥感图像的苏禄-苏拉威西海内孤立波特征研究[J]. 海洋与湖沼, 51(5): 991-1000.

ZHANG TAO, ZHANG XUDONG, 2020. Characteristics on internal solitary waves in the sulu-celebes sea based on MODIS and VIIRS remote sensing images[J]. Oceanologia et Limnologia Sinica, 51(5): 991-1000. (in Chinese with English abstract)

[4]
ALFORD M H, PEACOCK T, MACKINNON J A, et al, 2015. The formation and fate of internal waves in the South China Sea[J]. Nature, 521(7550): 65-69.

DOI

[5]
APEL J, OSTROSKY L, STEPANYANTS Y, et al, 2006. Internal solitons in the ocean[R]. Technical Report, Woods Hole Oceanographic Institution. WHOI-2006-04.

[6]
APEL J R, HOLBROOK J R, LIU A K, et al, 1985. The Sulu Sea internal soliton experiment[J]. Journal of Physical Oceanography, 15(12): 1625-1651.

DOI

[7]
BAI XIAOLIN, LIU ZHIYU, ZHENG QUANAN, et al, 2019. Fission of shoaling internal waves on the northeastern shelf of the South China Sea[J]. Journal of Geophysical Research: Oceans, 124(7): 4529-4545.

DOI

[8]
CAI SHUQUN, HE YINGHUI, WANG SHENGAN, et al, 2009. Seasonal upper circulation in the Sulu Sea from satellite altimetry data and a numerical model[J]. Journal of Geophysical Research: Oceans, 114(C3): C03026.

[9]
CAI SHUQUN, HE YINGHUI, 2010. Association of the Sulu Sea surface circulation with the South China Sea[J]. Journal of Marine Systems, 81(4): 335-340.

DOI

[10]
CAI SHUQUN, XIE JIESHUO, HE JIANLING, 2012. An overview of internal solitary waves in the South China Sea[J]. Surveys in Geophysics, 33(5): 927-943.

DOI

[11]
CAI SHUQUN, XIE JIESHUO, XU JIEXIN, et al, 2014. Monthly variation of some parameters about internal solitary waves in the South China Sea[J]. Deep Sea Research Part I: Oceanographic Research Papers, 84: 73-85.

DOI

[12]
CHEN LIANG, ZHENG QUANAN, XIONG XUEJUN, et al, 2019. Dynamic and statistical features of internal solitary waves on the continental slope in the northern South China Sea derived from mooring observations[J]. Journal of Geophysical Research: Oceans, 124(6): 4078-4097.

DOI

[13]
CHENG XUHUA, XIE SHANGPING, DU YAN, et al, 2016. Interannual-to-decadal variability and trends of sea level in the South China Sea[J]. Climate Dynamics, 46(9-10): 3113-3126.

DOI

[14]
CHO C, NAM S H, SONG H, 2016. Seasonal variation of speed and width from kinematic parameters of mode-1 nonlinear internal waves in the northeastern East China Sea[J]. Journal of Geophysical Research: Oceans, 121(8): 5942-5958.

DOI

[15]
CHOI W, CAMASSA R, 1999. Fully nonlinear internal waves in a two-fluid system[J]. Journal of Fluid Mechanics, 396: 1-36.

DOI

[16]
DU YAN, QU TANGDONG, 2010. Three inflow pathways of the Indonesian throughflow as seen from the simple ocean data assimilation[J]. Dynamics of Atmospheres and Oceans, 50(2): 233-256.

DOI

[17]
GONG YANKUN, XIE JIESHUO, XU JIEXIN, et al, 2022. Oceanic internal solitary waves at the Indonesian submarine wreckage site[J]. Acta Oceanologica Sinica, 41(3): 109-113.

[18]
GORDON A L, HUBER B A, METZGER E J, et al, 2012. South China Sea throughflow impact on the Indonesian throughflow[J]. Geophysical Research Letters, 39(11): L11602.

[19]
GRIMSHAW R, PELINOVSKY E, TALIPOVA T, et al, 2004. Simulation of the transformation of internal solitary waves on oceanic shelves[J]. Journal of Physical Oceanography, 34(12): 2774-2791.

DOI

[20]
GRIMSHAW R, PELINOVSKY E, TALIPOVA T, 2007. Modelling internal solitary waves in the coastal ocean[J]. Surveys in Geophysics, 28(4): 273-298

DOI

[21]
GRIMSHAW R, GUO CHUNCHENG, HELFRICH K, et al, 2014. Combined effect of rotation and topography on shoaling oceanic internal solitary waves[J]. Journal of Physical Oceanography, 44(4): 1116-1132.

DOI

[22]
GUO C, CHEN X, 2014. A review of internal solitary wave dynamics in the northern South China Sea[J]. Progress in Oceanography, 121: 7-23.

DOI

[23]
HAN WEIQING, MOORE A M, LEVIN J, et al, 2009. Seasonal surface ocean circulation and dynamics in the Philippine archipelago region during 2004-2008[J]. Dynamics of Atmospheres and Oceans, 47(1-3): 114-137.

DOI

[24]
HAO ZHANJIU, XU ZHENHUA, FENG MING, et al, 2021. Spatiotemporal variability of mesoscale eddies in the Indonesian Seas[J]. Remote Sensing, 13(5): 1017.

DOI

[25]
HE YINGHUI, FENG MING, XIE JIESHUO, et al, 2017. Spatiotemporal variations of mesoscale eddies in the Sulu Sea[J]. Journal of Geophysical Research: Oceans, 122(10): 7867-7879.

DOI

[26]
HELFRICH K R, MELVILLE W K, 2006. Long nonlinear internal waves[J]. Annual Review of Fluid Mechanics, 38(1): 395-425.

DOI

[27]
HUANG XIAODONG, CHEN ZHAOHUI, ZHAO WEI, et al, 2016. An extreme internal solitary wave event observed in the northern South China Sea[J]. Scientific Reports, 6: 30041.

DOI PMID

[28]
HURLBURT H E, METZGER E J, SPRINTALL J, et al, 2011. Circulation in the Philippine archipelago simulated by 1/12° and 1/25° global HYCOM and EAS NCOM[J]. Oceanography, 24(1): 28-47.

DOI

[29]
JACKSON C R, ARVELYNA Y, ASANUMA I, 2011. High-frequency nonlinear internal waves around the Philippines[J]. Oceanography 24(1): 90-99.

DOI

[30]
JIA T, LIANG J J, LI X M, et al, 2018. SAR Observation and numerical simulation of internal solitary wave refraction and reconnectionbehind the Dongsha Atoll[J]. Journal of Geophysical Research: Oceans, 123(1): 74-89.

DOI

[31]
JOSEPH R I, 1977. Solitary waves in a finite depth fluid[J]. Journal of Physics A: Mathematical and General, 10(12): L225-L227.

DOI

[32]
KAO T W, PAN F S, RENOUARD D, 1985. Internal solitons on the pycnocline: Generation, propagation, and shoaling and breaking over a slope[J]. Journal of Fluid Mechanics, 159: 19-53.

DOI

[33]
KLYMAK J M, PINKEL R, LIU C T, et al, 2006. Prototypical solitons in the South China Sea[J]. Geophysical Research Letters, 33(11): L11607.

[34]
KOOP C G, BUTLER G, 1981. An investigation of internal solitary waves in a two-fluid system[J]. Journal of Fluid Mechanics, 112: 225-251.

DOI

[35]
KORTEWEG D J, DE VRIES G, 1895. XLI. On the change of form of long waves advancing in a rectangular canal, and on a new type of long stationary waves[J]. Philosophical Magazine, 39(240): 422-443.

[36]
LAI ZHIGANG, JIN GUANGZHEN, HUANG YONGMAO, et al, 2019. The generation of nonlinear internal waves in the South China Sea: A three-dimensional, nonhydrostatic numerical study[J]. Journal of Geophysical Research: Oceans, 124(12): 8949-8968.

DOI

[37]
LI MINGTING, WEI JUN, WANG DONGXIAO, et al, 2019. Exploring the importance of the Mindoro-Sibutu Pathway to the upper-layer circulation of the South China Sea and the Indonesian Throughflow[J]. Journal of Geophysical Research: Oceans, 124(7): 5054-5066.

DOI

[38]
LI QIANG, WANG BING, CHEN XU, et al, 2016. Variability of nonlinear internal waves in the South China Sea affected by the Kuroshio and mesoscale eddies[J]. Journal of Geophysical Research: Oceans, 121(4): 2098-2118.

DOI

[39]
LI XIAOFENG, JACKSON C R, PICHEL W G, 2013. Internal solitary wave refraction at Dongsha Atoll, South China Sea[J]. Geophysical Research Letters, 40(12): 3128-3132.

DOI

[40]
LIAO GUANGHONG, XU XIAOHUA, LIANG CHUJIN, et al, 2014. Analysis of kinematic parameters of internal solitary waves in the northern South China Sea[J]. Deep Sea Research Part I: Oceanographic Research Papers, 94: 159-172.

DOI

[41]
LIU BINGQING, D’SA E J, 2019. Oceanic internal waves in the Sulu-Celebes Sea under Sunglint and Moonglint[J]. IEEE Transactions on Geoscience and Remote Sensing, 57(8): 6119-6129.

DOI

[42]
NAGAI T, HIBIYA T, 2015. Internal tides and associated vertical mixing in the Indonesian Archipelago[J]. Journal of Geophysical Research: Oceans, 120(5): 3373-3390.

DOI

[43]
QU TANGDONG, SONG Y T, 2009. Mindoro Strait and Sibutu Passage transports estimated from satellite data[J]. Geophysical Research Letters, 36(9): L09601.

[44]
SEGUR H, HAMMACK J L, 1982. Soliton models of long internal waves[J]. Journal of Fluid Mechanics, 118: 285-304.

DOI

[45]
TESSLER Z D, GORDON A L, JACKSON C R, 2012. Early stage soliton observations in the Sulu Sea[J]. Journal of Physical Oceanography, 42(8): 1327-1336.

DOI

[46]
WANG CAIXIA, WANG XIN, DA SILVA J C B, 2019. Studies of internal waves in the strait of Georgia based on remote sensing images[J]. Remote Sensing, 11(1): 96.

DOI

[47]
WANG JUAN, HUANG WEIGEN, YANG JINGSONG, et al, 2013. Study of the propagation direction of the internal waves in the South China Sea using satellite images[J]. Acta Oceanologica Sinica, 32(5): 42-50.

[48]
XIE JIESHUO, HE YINGHUI, CHEN ZHIWU, et al, 2015. Simulations of internal solitary wave interactions with mesoscale eddies in the northeastern South China Sea[J]. Journal of Physical Oceanography, 45(12): 2959-2978.

DOI

[49]
XIE JIESHUO, HE YINGHUI, HAIBIN, et al, 2016. Distortion and broadening of internal solitary wavefront in the northeastern South China Sea deep basin[J]. Geophysical Research Letters, 43(14): 7617-7624.

DOI

[50]
XIE JIESHUO, FANG WENDONG, HE YINGHUI, et al, 2021. Variation of internal solitary wave propagation induced by the typical oceanic circulation patterns in the northern South China Sea deep basin[J]. Geophysical Research Letters, 48(15): e2021GL093969.

[51]
XU ZHENHUA, WANG YANG, LIU ZHIQIANG, et al, 2021. Insight into the dynamics of the radiating internal tide associated with the Kuroshio Current[J]. Journal of Geophysical Research: Oceans, 126(6): e2020JC017018.

[52]
YUAN C, GRIMSHAW R, JOHNSON E, et al, 2018. Topographic effect on oblique internal wave-wave interactions[J]. Journal of Fluid Mechanics, 856: 36-60.

DOI

[53]
ZHANG XUDONG, ZHANG TAO, LI XIAOFENG, 2020. Satellite observation of tansmeridional propagating internal waves in the Celebes Sea[C]// IGARSS 2020-2020 IEEE international geoscience and remote sensing symposium. Waikoloa, HI, USA: IEEE: 6961-6964.

[54]
ZHANG Z, FRINGER O B, RAMP S R, 2011. Three-dimensional, nonhydrostatic numerical simulation of nonlinear internal wave generation and propagation in the South China Sea[J]. Journal of Geophysical Research: Oceans, 116(C5): C05022.

[55]
ZHAO XIAOYU, XU ZHENHUA, FENG MING, et al, 2021. Satellite investigation of semidiurnal internal tides in the Sulu-Sulawesi seas[J]. Remote Sensing, 13(13): 2530.

DOI

Outlines

/