关于内潮非线性演化的学术研究报告：聚焦超谐波级联机制

一、 研究团队与发表信息

本研究的核心作者为 Bruce R. Sutherland 和 Maninderpal S. Dhaliwal，他们均来自加拿大阿尔伯塔大学（University of Alberta），分属物理系以及地球与大气科学系。这项研究成果以论文形式发表在国际流体力学领域的权威期刊《Journal of Fluid Mechanics》上，发表于2022年，具体卷期为第948卷，文章编号为A21。论文标题为“The nonlinear evolution of internal tides. Part 1: the superharmonic cascade”，即《内潮的非线性演化。第一部分：超谐波级联》。

二、 研究背景与目标

本研究隶属于物理海洋学和流体动力学领域，具体聚焦于海洋内部潮汐（internal tides）的非线性演化过程。海洋内部潮汐是由正压潮（barotropic tide）在变化的海底地形上振荡时，将能量转化为斜压运动而产生的，是全球海洋能量输运和湍流混合的重要驱动力。观测发现，内潮在传播过程中，在某些区域（如南海）会演变为大振幅的内孤立波（internal solitary waves）列，而在另一些区域则能长距离传播而不发生显著变形。理解这种差异背后的物理机制，特别是非线性过程如何导致内潮波形变陡乃至孤立波的形成，是海洋动力学中的一个核心科学问题。

先前的研究（如Baker & Sutherland, 2020）已揭示，在非均匀分层（non-uniform stratification）的海洋中，一个具有单一垂直模态（vertical mode-1）的水平传播内波（即“母波”，parent wave）会通过自相互作用（self-interaction）激发其超谐波（superharmonic），即频率为母波整数倍（本研究主要关注二倍频）的波动。当母波振幅较小时，这种激发的超谐波（波数为2k）因其固有频率（ω₂）与强迫频率（2ω）存在微小失谐（off-resonant），其振幅会呈现周期性的增长与衰减，形成一个拍频现象。

本研究的核心目标是拓展上述理论，以处理更大振幅的母波和/或更强的近共振强迫情况，后者在科里奥利力影响较小的低纬度地区尤为显著。研究旨在探究在强非线性条件下，母波是否会通过级联过程激发一系列更高阶的超谐波（如3k, 4k…），即所谓的“超谐波级联”（superharmonic cascade），以及这种级联效应如何最终导致内潮波形演化为内孤立波列。研究试图建立一个比传统浅水理论限制更少、计算更快捷的数学模型，以深入理解内孤立波的生成机制。

三、 详细研究流程与方法

本研究采用了理论推导与数值模拟相结合的研究范式，主要包含以下几个步骤：

第一步：理论框架构建与方程推导 1. 建立物理模型与运动方程：研究考虑在f-平面（f-plane）上，无粘、无扩散、不可压缩的Boussinesq流体。控制方程包括动量方程和内能守恒方程（浮力方程）。通过引入流函数（streamfunction）简化二维（x-z平面）问题描述。 2. 线性波解与垂直结构：首先求解小振幅线性内波的解。母波被设定为具有特定水平波数k和频率ω的垂直模态-1波，其垂直结构函数 ψ̂₁(z) 由背景浮力频率剖面 N²(z) 和边界条件（自由滑移表面和底面）决定的本征值问题给出。研究采用了两种分层背景：纯指数分层和包含表面混合层的指数分层，以验证模型对实际海洋分层的适用性。 3. 非线性相互作用与模态耦合方程：将流函数表示为母波及其可能激发的所有超谐波（波数为 nk，频率近似为 nω，n=1,2,3…）的叠加。将这些展开式代入包含非线性项（如平流项）的完整运动方程中。 4. 推导振幅演化方程：利用多重尺度分析（multiple-scale analysis），引入慢时间尺度 T = εt（其中ε是度量失谐程度的小参数），并假设高阶时间导数可忽略。通过将非线性强迫项投影到各个垂直模态-1的结构上（利用本征函数的正交性），最终推导出一组描述各阶谐波复振幅 a_n(T) 演化的耦合非线性常微分方程（ODE）系统： [ \frac{da_n}{dT} - i(n-1)\omega b_n an = -i\frac{\alpha}{\omega} \sum{m+l=n, m\ge l} e_{ml} a_m a_l ] 其中，α = a₀/d 是母波的无量纲初始振幅（a₀为初始垂直位移振幅，d为分层特征深度），bn 是与频率失谐相关的系数，e{ml} 是三波相互作用系数，其值通过积分垂直结构函数及其导数得到，完全由背景分层N²(z)、波数k和科里奥利参数f决定。这项工作的重要进展在于，方程包含了母波与任意多个超谐波之间的相互作用，而不仅限于先前研究中的母波与二阶超谐波。

第二步：关键参数确定与系数计算 1. 参数设定：研究参数基于对夏威夷西南海域传播的内潮的远场观测数据（Rainville & Pinkel, 2006）进行设定。关键参数包括：海洋深度h≈5200 m，分层e折深度d≈218 m（对应d/h≈0.04），浮力频率尺度N₀≈0.017 s⁻¹，母波水平波数k≈4.2×10⁻⁵ m⁻¹（对应kh≈0.22），母波频率ω（M₂半日潮频率），以及科里奥利参数f（随纬度变化）。母波无量纲振幅α在0.025到0.125之间变化，覆盖了观测到的范围。 2. 系数计算：针对不同的kh、d/h和f/N₀组合，计算了理论模型中的关键系数ε（失谐参数）、bn和e{ml}。计算表明，对于长波近似（nkh << 1），相互作用系数e_{ml}主要贡献来自背景浮力频率的垂直梯度dN²/dz，且其值对分层细节（如有无表面混合层）相对不敏感，但对f的变化（尤其是f→0，即接近赤道时）导致的ε急剧减小非常敏感。

第三步：完全非线性数值模拟验证 1. 模拟设置：为了验证理论模型的预测，研究进行了完全非线性（fully nonlinear）的数值模拟。模拟在一个二维垂向剖面上进行，使用了一个经过验证的非线性、非静力、不可压缩的Boussinesq流求解器。模拟域水平方向采用周期性边界条件，垂直方向为刚盖自由滑移边界。初始条件设置为理论推导的垂直模态-1内潮波形（即母波）。 2. 模拟内容：在不同的母波初始振幅（α）和不同纬度（即不同f值）条件下运行了一系列数值实验。模拟跟踪了内潮波形随时间的演化，并通过谱分析等方法，定量提取了母波及各阶超谐波（波数k, 2k, 3k…）的振幅和相位随时间的变化。

第四步：理论模型求解与对比分析 1. 模型求解：对推导出的耦合ODE系统进行数值积分，初始条件为a₁(0)=1（母波归一化振幅），a_n(0)=0 for n>1（初始无超谐波）。求解得到各阶谐波振幅a_n(T)的演化过程。 2. 与模拟结果对比：将理论模型预测的各阶谐波振幅和相位演化，与完全非线性数值模拟的对应结果进行详细对比。特别关注不同α/ε比值下（即不同非线性强度和失谐程度下），超谐波级联的发展情况以及最终波形特征。

第五步：与经典浅水理论模型的联系 将本模型在强非线性（α/ε 较大）条件下预测的波形演化结果，与基于浅水理论的经典模型（特别是Korteweg–de Vries方程及其包含旋转效应的扩展——Ostrovsky方程）所预测的孤立波列生成进行对比，探讨两者的一致性，从而将新模型置于更广泛的波动理论框架中。

四、 主要研究结果

1. 理论模型的预测能力：研究发现，推导出的耦合非线性ODE系统能够高度准确地预测完全非线性数值模拟中内潮的演化过程。理论模型不仅重现了母波与二阶超谐波（2k）之间的周期性能量交换（当α/ε较小时），更重要的是，它成功地预测了当α/ε较大时发生的超谐波级联现象。
2. 超谐波级联现象：当母波振幅足够大（α较大）和/或频率失谐很小（ε很小，如在低纬度）时，即α/ε ≥ 1时，母波不仅激发二阶超谐波（2k），二阶超谐波与母波的相互作用还会有效地激发三阶超谐波（3k），进而可能激发更高阶的超谐波。这些相继产生的超谐波会增长到不可忽略的振幅。数值模拟和理论计算均清晰地展示了能量从母波（k）向2k、3k乃至更高波数模态的传递过程。
3. 向孤立波列的转变：理论模型和数值模拟均表明，在超谐波级联发生后，各阶超谐波之间会形成一种特定的相位关系。当所有这些谐波分量叠加时，初始平滑的正弦形式的内潮波形会发生显著变形，演变为一个具有陡峭波前和振荡波尾的孤立波列（solitary wave train）。这一发现为内潮如何通过非线性自相互作用演化为观测到的内孤立波提供了一种新的物理解释机制。
4. 参数依赖性：结果明确显示，内潮非线性陡化和孤立波生成的潜力强烈依赖于纬度。在科里奥利参数f较小的低纬度地区（接近赤道），ε值非常小（ω₂ ≈ 2ω），即使母波振幅α不大，α/ε的比值也可能很大，从而导致强烈的非线性相互作用和超谐波级联。这解释了为何在赤道附近区域更容易观测到内潮向孤立波的转化。
5. 与浅水理论的一致性：研究指出，在强非线性条件下，本模型所预测的波形演化与Ostrovsky方程（包含旋转效应的KdV方程）的预测结果一致。这建立了本基于垂向模态展开的“第一原理”模型与经典的、基于特定假设（如弱非线性、长波）的浅水模型之间的桥梁，证明了后者在描述此类现象时的有效性，同时也明确了本模型具有更少的限制性假设（如不要求长波近似kh << 1）。

五、 研究结论与意义

本研究通过发展一个包含母波与多阶超谐波非线性耦合的理论模型，并结合完全非线性数值模拟，系统地阐明了在非均匀分层海洋中，内潮通过“超谐波级联”过程向孤立波列演化的机制。

科学价值： 1. 理论创新：提出了一个描述内潮非线性演化的新理论框架，该框架通过一组耦合ODE刻画了连续超谐波的激发和相互作用，放松了先前理论中要求α/ε << 1的限制，能够处理更广泛的非线性情况。 2. 机制阐明：明确了“超谐波级联”是连接平稳传播的内潮与局地化孤立波列的关键非线性过程，揭示了各谐波分量间特定的相位锁定是波形变陡的核心。 3. 参数洞察：定量地揭示了内潮非线性行为对母波振幅和纬度（通过科里奥利参数f影响ε）的敏感性，为解释不同海域内潮演化差异提供了理论依据。

应用价值： 1. 高效预测工具：相较于耗时的完全非线性数值模拟，本研究推导的耦合ODE系统求解速度极快，为海洋学家和气候学家提供了一个潜在的高效工具，用于快速评估在不同环境参数（振幅、分层、纬度）下，内潮非线性演化的可能结局（是保持波形还是演化为孤立波列）。 2. 改进参数化方案：该模型有助于改进海洋环流和气候模型中关于内潮能量耗散和混合的参数化方案，因为内潮演化为孤立波列通常会加速能量的局地耗散。

六、 研究亮点

1. 重要的新发现：首次在理论和数值上系统演示并命名了内潮演化中的“超谐波级联”现象，并将其确定为内孤立波生成的一种前驱机制。
2. 新颖的理论方法：建立了基于垂直模态展开的、可处理任意数量超谐波相互作用的耦合振幅方程理论，该理论具有坚实的流体力学基础，且形式简洁。
3. 研究对象的特殊性：聚焦于低纬度（f小）情境，揭示了内潮在赤道附近更强的非线性演化倾向，这与观测事实相符，填补了该区域相关理论研究的空白。
4. 完美的验证闭环：理论推导、快速ODE模型求解、完全非线性高精度数值模拟三者相互印证，构成了一个完整、严谨的研究体系，结论可靠性高。

七、 其他有价值的内容

研究还详细讨论了背景分层结构（如加入表面混合层）对相互作用系数e_{ml}的影响，发现其影响较小（%），这表明模型对实际海洋中常见的分层细节具有一定的鲁棒性。同时，研究也指出，若改变深海的分层结构（如减弱浮力频率的指数衰减），系数可能更敏感，这为后续研究指明了方向。附录中给出了在长波近似下相互作用系数的简化解析表达式，增加了理论的实用性。