马鲁古海卫星可见内孤立波（ISWs）波前特征与两阶段预测框架：基于MODIS和SWOT的2002—2025年观测

《Applied Ocean Research》：Satellite-observed characteristics and prediction of internal solitary waves in the Maluku Sea

　　

内孤立波（ISWs）是层化海洋中常见的高能非线性内波，可诱导强流、垂向等密度面位移和增强混合。尽管其主要动力学过程发生在海表以下，ISWs常可通过海表粗糙度调制被空间观测识别，并在光学和合成孔径雷达（SAR）影像中表现为明暗相间条带；近期宽幅测高亦揭示了与ISWs相关的海表高度异常信号。卫星遥感因而为记录ISWs的发生、空间分布和海表表现提供了有效手段，但此类观测本质上限于海表，可揭示卫星可见波前和海表表现，而波振幅、垂向结构和能量通量的直接信息仍需现场观测或高分辨率数值模拟。印度尼西亚海域是全球最活跃的ISW区域之一，强正压潮与陡坎、复杂地形和强层结的相互作用有利于其活动。以往研究已利用卫星影像、现场观测、理论分析和数值模拟研究了Sulu Sea、Celebes Sea、Lombok Strait、Flores Sea和Banda Sea等区域的ISWs，表明半日潮强迫、地形转化和背景层结是控制ISW生成和传播的关键因素。近期研究已从分布制图转向利用遥感和机器学习方法预测ISW发生、波前演化和传播。马鲁古海位于Sulawesi与Halmahera之间，连接西太平洋、Celebes Sea和Banda Sea，是印度尼西亚海域内动力学上重要的区域，并与Banda Sea共同构成太平洋与印度洋路径之间水团交换和区域环流的关键枢纽。以往研究报道了马鲁古海及邻近印度尼西亚水域中明显的ISW信号，提供了ISW反复活动及区域潮汐和地形影响的证据；在邻近Banda Sea，合成孔径雷达（SAR）观测已识别出长波前和相对较高的传播速度，报道的相速度为2.2–3.0 m/s。然而，与Sulu Sea、Celebes Sea和Lombok Strait相比，马鲁古海受到的体系性关注较少，该区域卫星可见ISW波前的长期卫星统计和面向预测的分析仍有限。近期卫星观测进展为研究此类数据稀疏区域的ISW海表信号提供了新机会。传统光学和SAR传感器主要捕捉海表粗糙度型态，使波前和波包结构可在大范围制图；地表水和海洋地形（SWOT）任务通过观测与内波信号相关的二维海表高度异常提供了互补视角。SWOT采用被星下点间隙分隔的双50 km条带，可分辨约15 km或更长波长的海洋信号；因此，对ISW研究而言，SWOT观测可通过提供海表高度扰动信息补充光学和SAR粗糙度信号。将长期光学影像与SWOT观测整合，可能改善马鲁古海卫星可见ISW海表信号的表征。本研究建立在以往区域研究基础上，并非声称首次发现马鲁古海ISWs；其贡献在于汇编多年卫星影像数据集，绘制卫星可见ISW波前的海盆尺度分布，将观测路径与潮汐和地形诊断比较，并检验历史遥感波前几何能否支持海表波前位置预测。因此，本研究目标是刻画并预测马鲁古海卫星观测ISWs的海表信号；研究人员不尝试直接反演ISW振幅、湍流混合或能量通量。尽管如此，在現場监测稀疏区域预测ISW海表信号具有价值，因为基于遥感的预测框架可为卫星任务规划、水下航行器路由和海洋灾害评估提供有用先验信息。基于固定源区和平均相速度的传统估计可提供近似到达时间，但不能表征二维波前几何或源区变异性。为满足该需求，研究人员开发了卫星可见ISW波前的两阶段预测框架：第一阶段预测源区附近的初始波前，第二阶段使用基于卷积长短期记忆（ConvLSTM）的序列模型预测后续传播速度和方向。MODIS存档显示卫星可见ISW波前的两条主要路径：主导分支起源于Lifamatola Passage附近并向北穿越海盆传播，而较弱分支从Celebes Sea进入。波前长度从源区附近数十公里增至中央海盆数百公里，随后在北部海槽受到约束。Terra/Aqua匹配影像对表明传播主要向北或西北向，表观相速度主要在0.5–2.5 m/s之间。SWOT观测揭示与ISW相关的SSHA信号具有日变化和大潮—小潮调制。这些结果表明卫星可见ISW信号与潮汐强迫之间存在统计联系，但未直接验证次表层振幅、位移或能量通量。研究人员进一步开发了预测卫星可见ISW波前的两阶段框架：STEP1根据时间和潮汐条件输入预测源区波前掩膜，STEP2使用ConvLSTM模型估计后续传播速度和方向。STEP1在测试集取得F1和交并比（IoU）分数分别为0.831和0.711；STEP2取得速度平均绝对误差（MAE）/均方根误差（RMSE）/决定系数（R²）为0.09 m s^?1/0.13 m s^?1/0.97，方向MAE/RMSE为3.0°/5.7°。结果表明，历史卫星反演波前几何为预测海表ISW信号提供了有用先验信息，但该框架不应被解释为完整的次表层动力学模型。

研究背景与问题：内孤立波（internal solitary waves, ISWs）是层化海洋中常见的高能非线性内波，可诱导强流、等密度面位移和增强混合；其主要动力学位于海表以下，但可通过海表粗糙度调制在光学和SAR影像中形成明暗条带，宽幅测高亦可揭示海表高度异常（sea surface height anomaly, SSHA）信号。卫星遥感适合记录ISWs发生、空间分布和海表表现，但本质上限海表，不能直接提供振幅、垂向结构和能量通量。印度尼西亚海域因强正压潮、陡坎、复杂地形和强层结而ISW活跃；马鲁古海（Maluku Sea）连接西太平洋、Celebes Sea和Banda Sea，是水团交换和区域环流关键区，但相较Sulu Sea、Celebes Sea和Lombok Strait，其长期卫星统计和面向预测的分析仍有限。研究人员因此刻画并预测马鲁古海卫星观测ISWs海表信号，不直接反演振幅、混合或能量通量，并构建两阶段框架以提供遥感先验信息，服务卫星任务规划、水下航行器路由和海洋灾害评估。研究人员汇编多年卫星影像，绘制卫星可见波前分布，比较路径与潮汐/地形诊断，并检验历史遥感波前几何能否支持海表波前位置预测；结论识别两条复发路径，主导路径自Lifamatola Passage附近向北延伸，较弱路径从Celebes Sea进入；SWOT SSHA呈现潮汐调制；两阶段模型给出可用的波前位置和传播先验，但并非完整次表层动力学模型。

关键技术方法（≤250字）：研究人员构建2002—2025年MODIS Terra/Aqua长时序数据集并人工标注清晰波前，利用Terra/Aqua约3 h匹配估计表观位移；采用FY-4A/B、GK-2A提高源区时间采样，并以LLC4320小时输出和1/48°分辨率作为辅助波前样本；使用SWOT KaRIn L3低速率250 m产品，经质控、陆雨掩膜和大尺度背景去除分析Pass 021 SSHA。STEP1以月相和半日潮相位编码及波前掩膜/坐标构建1102个样本（163个MODIS、68个静止卫星、871个LLC4320；772/220/110划分），采用编码器—解码器分割；STEP2以8136条MODIS向量训练ConvLSTM（卷积长短期记忆网络），输入经纬度、水深、月均密度代理及滞后速度/方向正余弦，并加入时间平滑损失。

研究结果：**3.1 空间分布与波前长度**：通过2002—2025年MODIS存档人工标注（472幅Aqua、457幅Terra含ISW特征）得出两组卫星可见波前，主导组从Lifamatola Passage附近向北穿越海盆，较弱组从Celebes Sea进入；波前长度多为60—200 km，近源常<120 km，赤道附近约200 km，1°N附近偶超350 km，1.5°N以北受西部海槽约束缩至50—150 km；该分布应解释为可观测波前发生型态，而非全部次表层ISW。**3.2 匹配MODIS影像的传播速度与方向**：通过113对Terra/Aqua约3 h匹配获得8724个位移矢量，主导分支主要向北至西北传播，Celebes Sea分支有少量向东传播；表观速度多0.5—2.5 m/s，少数>3 m/s；3月平均速度最高2.2 m/s、12月最低1.6 m/s；多数发生于>1000 m水深，集中1500—3000 m，南部深区1—2 m/s；这些矢量表示卫星可见波前的表观水平位移，而非直接观测的次表层界面运动。**3.3 SWOT观测的SSHA信号与潮汐强迫**：通过Pass 021校准/验证期98个有效观测日，去除大尺度背景后识别窄而相干SSHA扰动；2023年4月20日至23日量值分别为13.3、14.75、13.05、10.93 cm；日平均SSHA随源区潮高大小潮调制，较本地潮汐有1—2天偏移，被描述为传播一致而非追踪单个波包；谱峰约13.7 d，0.5°N—1°N信号较强；结果支持潮汐调制解释，但不能验证次表层密度或垂向位移。**4.1 STEP1-ISW起始位置预测**：通过将观测时间、月相和半日潮相位与波前几何关联，构建源区初始波前概率掩膜模型；测试F1/IoU为0.831/0.711；结论为在相似时间和潮汐条件下提供源区波前先验估计，而非确定性探测每个ISW事件。**4.2 STEP2-ISW传播预测**：通过8136条MODIS向量序列训练两层ConvLSTM（隐藏维64、核3），输入经纬度、水深、月均上下层密度代理及连续性滞后变量，目标为下一步速度与方向正余弦，损失含时间平滑（α=0.05）；速度MAE/RMSE/R²=0.089 m/s/0.128 m/s/0.968，方向MAE/中位/RMSE=3.03°/1.82°/5.65°；模型再现主要传播变化，方向样本集中和稀疏角度可能影响相关评价。**4.3 评估指标**：通过像素真阳性/假阳性/假阴性计算F1与IoU，以MAE、RMSE、R²及圆形RMSE评价STEP2，角度差包裹至?180°—180°；STEP1训练F1/IoU=0.906/0.828、测试=0.831/0.711，训练测试差距中等，IoU低于F1反映窄波前对像素偏移敏感。**4.4 ISW预测结果**：通过STEP1四个代表性案例比较预测与观测初始波前，距离误差5.5—19.5 km、形状相似度80%—88%；进一步以2 h步长耦合STEP1—STEP2开展海盆尺度预测，与MODIS比较，距离误差8.5—13.7 km，方向误差约13°降至约6—9°，多数形状相似度>90%、最高95%；模型具有稳健预测技能，非线性增强和多周期潮汐强迫下误差增大，比较存在约一个潮周期时间偏移，主要反映传播结构和轨迹一致性。

讨论总结：**5.1 潮汐调制与源区可预测性**：SWOT显示日变化和 fortnightly 调制；LLC4320 36小时（2011年10月1日03:00至10月2日15:00 UTC）三个半日潮周期与TPXO9比较显示，前两较强周期波前更大更连续，第三较弱周期波包更小更短；波前在北向潮流减弱并转南向时出现，呈相位锁定，支持STEP1潮汐描述符，但不是生成机制直接验证。**5.2 模型误差特征与应用价值**：误差图显示STEP1真阳性沿复发源区路径，假阳性和假阴性多在波前带边缘；STEP2速度误差多低于约0.13 m/s，方向误差多在约3.5°内（平滑图），受路径几何和地形影响；框架可提供二维波前几何和方向依赖传播先验，支持卫星任务、海上作业规划和区域ISW风险评估，尤其无实时现场观测区。**5.3 局限与未来工作**：仅基于海表粗糙度和SSHA，不能直接获取等密度面位移、振幅和能量通量；云、太阳耀斑、单传感器日采样和Terra/Aqua约3 h限制完整生命周期；LLC4320仅作源区波前几何辅助，不能解析非线性ISW；STEP2未显式含背景流和时变层结；未来需Lifamatola Passage及北向路径系泊观测、高分辨率非静力模拟、源区分层试验和前瞻性时间分块验证。

研究结论（第6节）翻译：本研究刻画马鲁古海卫星可见ISW信号，并评估其复发几何能否支持海表波前预测。分析结合2002—2025年MODIS存档、匹配Terra/Aqua观测、SWOT SSHA、潮汐诊断和两阶段机器学习框架。识别两条复发海表波前路径：主导路径初现于Lifamatola Passage附近，并在部分场景向北延伸>700 km；较少见路径从Celebes Sea进入。观测波前长度常为60—200 km，偶超300 km。表观传播主要为向北至西北向，速度0.5—2.5 m s^?1。这些结果描述卫星可见波峰几何和运动，不确立次表层振幅演变。SWOT观测显示日际SSHA变化和约13.7 d谱峰，符合大小潮调制；观测相位关系与辅助LLC4320型态支持STEP1使用潮汐描述符，但不构成生成机制直接验证。基于物理认识，研究人员建立混合预报框架，将潮汐相位信息与卫星反演波前耦合，预测ISW起始和传播。框架采用两步设计：STEP1联系潮汐相位、月周期与观测波前形态重建初始波前位置；STEP2用ConvLSTM预测受水深和层结约束的传播速度与方向。STEP1混合源测试F1=0.831、IoU=0.711；STEP2速度MAE/RMSE/R²=0.09 m s^?1/0.13 m s^?1/0.97，方向MAE/RMSE=3.0°/5.7°。案例应用表明框架可在真实卫星观测条件下提供波前位置和传播路径有用先验估计。该框架可为数据稀疏区遥感规划和海洋决策支持提供海表波前位置与方向先验；目前范围不是业务化次表层ISW预报。需现场剖面、系泊流、高分辨率非静力模拟、源区分层试验和前瞻性时间分块验证后，才能声称次表层动力学或业务可靠性。