观潮 · 模式预报与气象 AI 观察
波罗的海波浪谱 12 小时回报：物理约束只把波高 RMSE 压低 5.01%，且分目标互相打架
> 用浮标实测的一维波浪谱做 12 小时回报，把「监督哪些谱描述符」和「加什么物理约束」逐个拆开比。结果是收益分目标：监督方差加峰值频率（A2）谱误差最低（0.40 m²·Hz⁻¹），再加峰值振幅监督（A3）才是波高误差最低（18.78 cm）；物理侧，局地松弛约束把波高 RMSE 压低 5.01% 到 19.26 cm，把松弛延伸到 13—24 小时（那段没有谱目标）才把谱 RMSE 降到 0.39，而经验风浪上限约束反而让谱误差涨 8.44%、波高涨 4.89%。作者的定性是「modest, objective-dependent benefits」，并且明确说这是回溯性的（hindcast），业务化预报和向独立站点的迁移都还没评估。
缺的是哪一块
有效波高是一个积分量，但方差相同的两段谱可以在峰的位置、峰的密度和形状上完全不同，所以只看波高精度不足以评价谱预报。既有工作里，LSTM 已被用于有效波高，也有方法直接预测谱本身——有的用 ERA5 风场建模单点方向谱，有的把局地风与远端风结合、以数值回报谱为目标，有的做潜压缩重建的全球方向谱预报，还有的把业务方向预报订正到一维浮标谱。
物理知识的引入方式大致分两类：一是对观测量（方差、峰频、形状）加监督，二是对演变过程加约束（如 PINN 把控制方程残差放进损失）。前者约束的是「与观测一致」，后者约束的是「与某个假定物理关系一致」，两者并不等价。作者指出既有研究很少把这两类放在同一个受控框架里、针对「实测浮标谱随后如何演变」这个具体任务做对比。
方法怎么搭的
任务是：用一个小时的浮标谱初始化，输入局地风、方向性风区长度和水深，预测未来 12 小时的谱。风由 ERA5 再分析在整个预报区间内全程提供——这决定了本研究是回溯性的 hindcast，而不是只用到起报时刻可用信息的预报。
数据来自芬兰气象局的三套 Datawell Directional Waverider Mark 4 浮标：
站 · 纬度 · 经度 · 水深 · 训练 · 验证 · 测试
NBP（北波罗的海本部） · 59.25°N · 21.00°E · 107.0 m · 27,685 · 8,553 · 7,695
BS（波的尼亚海） · 61.80°N · 20.23°E · 112.2 m · 23,844 · 8,515 · 8,633
GoF（芬兰湾） · 59.96°N · 25.24°E · 49.7 m · 22,264 · 6,794 · 6,851
谱重采样到小时分辨率并对齐小时风，线性插值到 0.085—0.595 Hz 之间的 52 个频率点（标称间隔 0.01 Hz）。这套插值不显式保持积分方差，因此下文报的方差与波高是按这个截断频率网格算的，会与用原始谱全频段算出的结果不同。训练目标年 2017—2020，验证 2021，测试 2022，样本量 73,793 / 23,862 / 23,179 个起报时刻，三站都参与各时期，相邻起报时刻重叠、不代表独立事件。
对比 variants 共用同一套输入和网络，只在损失项上做加法：
编号 · 名称 · 目标函数
A0 · 谱基线 · 对数谱重建 + 时间增量
A1 · 方差监督 · A0 + 实测谱方差
A2 · 峰频监督 · A1 + 实测软峰频
A3 · 峰幅监督 · A2 + 实测峰值谱密度
P1 · 风浪上限 · A2 + 风浪方差与峰周期的上界
P2 · 松弛 · A2 + 第 1—12 小时的局地松弛残差
P3 · 扩展松弛 · P2 + 第 13—24 小时的残差（该段无谱目标）
结果数字
描述符组（1—12 小时，全部 23,179 个测试起报时刻）：
变体 · E RMSE (m²·Hz⁻¹) · fp RMSE (Hz) · Ap RMSE (m²·Hz⁻¹) · 波高 RMSE (cm) · 偏差 (cm) · r · SI
A0 · 0.41 · 4.95×10⁻² · 1.60 · 19.54 · −6.26 · 0.97 · 0.20
A1 · 0.41 · 4.90×10⁻² · 1.48 · 19.00 · −3.98 · 0.97 · 0.19
**A2** · **0.40** · **4.80×10⁻²** · 1.56 · 19.62 · −5.08 · 0.97 · 0.20
A3 · 0.43 · 4.88×10⁻² · **1.33** · **18.78** · −4.32 · 0.97 · 0.19
物理约束组（相对配对的 A2 参考）：
变体 · E RMSE · fp RMSE · Ap RMSE · 波高 RMSE (cm) · 偏差 (cm) · r · SI
A2（参考） · 0.40 · 4.80×10⁻² · 1.61 · 20.28 · −6.64 · 0.97 · 0.21
P1 · 0.43 · 4.83×10⁻² · 1.80 · 21.27 · −8.02 · 0.96 · 0.22
**P2** · 0.40 · 4.82×10⁻² · 1.57 · **19.26** · −5.71 · 0.97 · 0.20
**P3** · **0.39** · 4.80×10⁻² · 1.46 · 19.66 · −6.08 · 0.97 · 0.20
换算成相对变化：A2 的谱误差与峰频误差比基线低 2.63% 与 3.09%；P1 让谱 RMSE 涨 8.44%、波高 RMSE 涨 4.89%；P2 让波高 RMSE 降 5.01% 到 19.26 cm，谱误差略有上升；P3 谱与波高误差分别降 2.31% 与 3.04%。技巧评分上，P3 的谱技巧为 36.35%，P2 的波高技巧为 51.55%。所有学习模型平均都偏低估，P2 偏差为 −5.71 cm，相关系数 0.97。误差随时效增长，但比持续性预报慢。
分站看（相对持续性预报的技巧评分 S）：
站 · 变体 · 谱 RMSE · 谱 S · 波高 RMSE (cm) · 波高 S
BS · A2 · 0.48 · 33.06% · 22.54 · 49.69%
BS · P3 · 0.47 · 34.55% · 21.47 · 52.10%
GoF · A2 · 0.14 · 37.77% · 13.87 · 50.11%
GoF · P3 · 0.14 · 39.29% · 13.96 · 49.78%
NBP · A2 · 0.45 · 36.65% · 22.30 · 47.74%
NBP · P3 · 0.44 · 38.16% · 21.77 · 48.98%
风况分层的结果：扩展松弛带来的谱收益出现在风速增大、减小和稳定的强风情形下，稳定弱风下几乎无变化。
作者承认的限制
- 这是 hindcast，不是预报。预报区间内的风来自 ERA5 再分析，原文明确「the study therefore concerns retrospective hindcasting」，并说业务化预报仍有待评估。
- 平衡关系的假设与实测状态之间隔了好几层。原文列出：局地水深是对风区上 bathymetry 的近似；瞬时风区长度加上一个很小的端点风速变化，并不能证明强迫已经持续到足以达到平衡；频率截断与固定的 JONSWAP 投影分辨不了交叉的方向性系统。这些在波罗的海尤其要紧，因为风向会改变有效风区，风减弱后局地生成的浪还能维持。
- 基于波龄和风—浪方向的涌浪标签，不等于长周期大洋涌浪。
- 标定允许偏离估计的平衡态，并不独立地验证这个平衡态。容差与站点排除规则是从训练观测里推出来的；12 小时的松弛时间是按残差准则选的，不是实测的物理调整时间。
- 选出的松弛闭合本身跑不过持续性预报。原文说对于符合条件的训练转变，单独的松弛闭合并不优于小时持续性预报，它的价值只在于作为对学习谱的约束。
- 标定规则、站点掩码和损失权重的影响没有被分离，因此原文说结果支持的是这一套具体设定，而不是普适的参数选择。
- 谱是带限的。0.085—0.595 Hz 的截断使得方差与波高与用原始谱全频段计算的结果不同。
- 评估只在训练所用的三个站点上做，向独立站点的迁移尚未评估。
一句话判断
把「监督什么」和「约束什么」放在同一张表里对比，得到的是一句不讨喜但很实在的结论：物理引导的收益只有 5% 上下，而且改善谱和改善波高走的是两条不同的路子——追求谱精度就选 A2 或把松弛外推到无目标时段，追求波高就选 A3 或 12 小时松弛，经验风浪上限约束则两头都不讨好。加上风用的是再分析、站点又是训练站，这份结果目前只能算方法学对照，离业务可用还有明确距离。