《Frontiers in Marine Science》：Assessing tidal sensitivity to bottom friction: analytical insights for model optimization
编辑推荐：
底摩擦系数（bottom friction coefficients，BFCs）是潮汐模拟中的关键参数，对水动力条件和沉积物输运具有显著影响。全面理解BFC的变异性及其对潮汐模型的影响对于模型优化至关重要。本研究通过Dirac delta函数将BFC扰动引入一
底摩擦系数（bottom friction coefficients，BFCs）是潮汐模拟中的关键参数，对水动力条件和沉积物输运具有显著影响。全面理解BFC的变异性及其对潮汐模型的影响对于模型优化至关重要。本研究通过Dirac delta函数将BFC扰动引入一个修正的、深度均匀的理想化浅水模型中。通过将控制方程变换到极坐标系，非齐次系统被重构为具有动力学一致偶极子源的齐次形式。由Bessel函数组合构建的解析解揭示了一系列旋转Poincaré波，这些波以基态潮汐频率的奇数倍频率围绕扰动点传播。Bessel函数的权重非线性地依赖于局部水深、潮汐流振幅和相位滞后，决定了旋转方向（顺时针、逆时针或静止），并为潮汐系统对BFC扰动的敏感性提供了理论指标。解析解捕捉了测试案例中的主要空间模态和相对响应特征。理想化M2实验与解析解表现出最佳的定量一致性。渤海模拟保持了与理论的定性一致性，尽管近岸区域的定量一致性有所下降。借助该框架，研究人员从理论上评估了渤海、黄海和东海区域的BFC敏感性。分析表明，浅水强流区域（包括部分潮波入口和出口区域）的敏感性相对较高。这一发现可为针对性的模型校准提供基于物理的先验信息。
论文解读：潮汐对底摩擦系数扰动的敏感性——解析框架与模型优化启示
1. 研究背景与问题
在海洋与海岸带数值模拟中，底摩擦系数（BFCs）是控制能量耗散及动量、物质垂向输运的关键参数。然而，BFC受地形、底床粗糙度和流场等多因素影响，直接观测几乎不可行，通常需借助参数化方案或数据同化方法进行估计。数据同化虽能在大区域优化BFC，但高维优化问题计算代价高昂且易陷入局部极小。为提升效率，研究者常选取若干关键点进行BFC优化，其余区域通过插值获得。关键点位置的选择对优化成败至关重要，但现有研究多采用均匀网格布点，缺乏对潮汐系统响应BFC扰动内在机理的理解。因此，从正向敏感性角度出发，解析地刻画局部BFC扰动所引发的潮汐波场响应，可为关键点选取和模型校准提供物理基础。该研究正是针对这一需求展开。
2. 研究内容与主要结论
研究人员基于线性浅水方程，在极坐标系中推导了描述局部底摩擦系数扰动响应的解析解。通过引入Dirac delta函数将局地扰动数学化，并将非齐次方程转化为带偶极子源的齐次问题，利用分离变量法和Bessel函数组合构建了通解。该解被识别为一族旋转Poincaré波，其频率为基频潮汐频率的奇数倍，权重系数非线性依赖于局部水深、潮流振幅及相位。该解析框架随后在理想化矩形海域和真实渤海环境中通过数值实验验证，并进一步应用于渤海、黄海和东海的大范围敏感性评估。研究表明：理想化M2实验定量吻合最佳；渤海区域定性一致但近岸定量偏差增大；浅水强流区域（包括潮波入口和出口）往往表现出更高的BFC敏感性。该框架可作为物理驱动的先验信息，辅助潮汐模型的目标化校准与优化。
3. 主要技术方法
研究人员采用了以下关键技术方法：一是线性化扰动方法，将BFC扰动视为点源并利用Dirac delta函数表示，使控制方程转化为可解析处理的非齐次系统；二是极坐标变换与分离变量法，将原始方程化为极坐标下齐次波动方程，求得Bessel函数型通解；三是Fourier-Bessel级数展开与系数匹配，通过比较偶极子源表达式确定奇异N1分量的解析权重；四是数值验证与最小二乘拟合，使用Princeton Ocean Model（POM）外部模态开展理想化实验和渤海实验，并通过最小二乘方法估计正则J1分量；五是敏感性指标构建，利用解析N1系数的平均绝对值（MAV）作为BFC敏感性的对比指标，数据来源包括FES全球潮汐模型数据集和GEBCO水深数据集。
4. 研究结果
4.1 理想化矩形域实验
研究人员在均匀水深1000 m、范围100°E–180°E、20°N–50°N的理想化矩形域中，以0.5°分辨率运行POM外部模态，仅受M2引潮力驱动。在域中心（140°E, 35°N）施加ΔCd=1.0×10??的BFC扰动后，数值模拟的高频扰动场呈现顺时针旋转的Poincaré波，周期约13小时，振幅量级5×10?? m。基于解析N1系数的纯理论解以及拟合J1分量的重构波场均能捕捉主要空间模态。定量评估表明：仅用纯N1特定解时相关系数R=0.461，拟合J1后R升至0.880，同时拟合N1与J1时R=0.931，说明理论预定的N1系数已接近数值模型的最优系数。
4.2 网格分辨率敏感性
研究人员在0.5°、0.25°、0.125°三种分辨率下重复理想化实验，以检验有限网格面积与点源假设的偏差。结果显示，纯N1解的相关系数从0.514逐步升至0.546，相对L2误差从0.859降至0.838，改善有限但呈单调趋势。理论约束J1拟合后，相关系数（0.937、0.932、0.925）未随分辨率提高而增加，说明有限面积效应只是解析-数值偏差的一个次要来源，有限域调整和其他数值效应仍起重要作用。
4.3 渤海真实环境实验
研究人员使用1/25°分辨率POM外部模态，水深取自GEBCO数据，开边界由FES全球潮汐模型提供，仅考虑M2分潮，背景BFC设为0.0013。当扰动位于渤海中部（38.56°N, 119.86°E）时，旋转Poincaré波周期约13小时，振幅达5×10?? m，纯N1解与数值结果的相关性R=0.537，Willmott指数d=0.721；同时拟合N1和J1后R=0.920，RMSE降至4.50×10?? m。当扰动位于近岸（39.6°N, 119.1°E）时，振幅仅约1×10?? m，波形出现更短波长成分，被解释为非线性生成的M4倍潮；纯N1解失效（R=0.181），但拟合N1和J1后R=0.852。近岸实验说明理论应用边界受到水深梯度、边界约束和高阶非线性效应的限制。
4.4 渤黄东海敏感性分布
研究人员利用FES数据提供的M2潮流调和常数和GEBCO水深，计算了各网格点上解析N1系数的平均绝对值（MAV），作为BFC敏感性指标。结果显示：较高指标值主要出现在浅水强流海岸区域，包括黄海部分海岸及强潮交换区；水深较深的开阔海域敏感性较低。指标近似满足U2H?3/2的标度规律，并受科氏参数调制。该分布特征为潮汐模型的定向校准提供了物理依据。
5. 讨论与结论
讨论部分指出，扰动点附近水流减速导致上游壅水和下游亏水，形成偶极子状压力异常；科氏力使该结构旋转，产生向外辐射的Poincaré波，这是旋转波结构形成的物理解释。该解析框架与海山内潮生成、岛屿近惯性运动研究在数学方法上相似，但研究动机不同——本研究聚焦于BFC正向敏感性分析而非内波激发机制。框架的局限性包括：假设局地均匀水深、常值科氏参数、小扰动、线性正压动力学及指定背景潮，在缓变地形和弱非线性区域表现最佳，近岸强非线性环境不适用。未来可扩展至有限面积源、变水深、多分潮和非线性相互作用情形。结论强调：该理论框架推导了BFC扰动引起的旋转Poincaré波解析解；数值实验验证了解析解对主导扰动模态的捕捉能力；在渤海、黄海、东海M2应用中，浅水强流区如潮波入口和出口区域敏感性更高；该框架可为潮汐与海洋模型的敏感性分析和校准设计提供基于物理的先验信息。
