文献阅读(45)
题目
Sizing the largest ocean waves using the SWOT mission
第一作者与通讯作者信息
- 第一作者：Fabrice Ardhuin（法国）
- 通讯作者：Fabrice Ardhuin（邮箱：fabrice.ardhuin@gmail.com）
- 单位：Université de Brest, CNRS, Ifremer, IRD, LOPS, IUEM, Plouzané, France
第一作者其他三篇代表性著作
Ardhuin, F., et al. (2024).
Phase-resolved swells across ocean basins in SWOT altimetry data: Revealing centimeter-scale wave heights including coastal reflection.
Geophys. Res. Lett., 51, e2024GL109658.
→ 首次利用 SWOT 二维海面高度数据提取全球涌浪相位，奠定本文方法基础。Ardhuin, F., & De Carlo, M. (2025).
Storm tracks based on wave heights from LOPS WAVEWATCH III hindcast and ERA5 reanalysis, 1991–2024.
Dataset, SEANOE. https://doi.org/10.17882/105148
→ 构建 34 年全球极端海浪事件目录，为本文“Eddie”等风暴对比提供基准。Accensi, M., Ardhuin, F., et al. (2025).
WAVEWATCH III model simulations GLOB-30M_LOPS_2025.
Ifremer Sextant repository.
→ 提供本文校验与谱形更新的高分辨率数值 wave 模式输出。
摘要
> 利用 SWOT 卫星二维海面高度数据，我们追踪到远离风暴的涌浪，其波长与高度随传播距离呈幂律变化：波长 ∝ d²，涌浪高 ∝ d⁻⁹。据此反推出风暴内谱形在峰频以下应为 f⁻¹⁷，而非传统 JONSWAP 的 f⁻⁵。更新后的谱形使最长波能量降低 20 倍，并首次给出卫星测得的最大有效波高 19.7±0.3 m 对应峰周期 20.2±0.6 s。该成果可用于海岸工程、地震噪声解释及古风暴重建。
相关研究的重要性
| 序号 | 重要性维度 | 具体阐述 |
|---|---|---|
| 1 | 极端海况安全 | 现有设计标准对 >16 m 海浪缺乏实测校验；明确峰周期-波高关系可降低海洋结构失效风险。 |
| 2 | 气候变迁指标 | 极端海浪能量过去 30 年是否增长存在争议；长周期涌浪记录可作独立气候指标。 |
| 3 | 地震噪声源 | 18–26 s 微震源长期不明；本研究把海浪峰周期与微震频段直接挂钩，助力“地震海洋学”。 |
| 4 | 古风暴重建 | 沉积物波痕、漂砾尺寸需同时知道波高与周期；新谱形可提高古风暴强度反演精度。 |
前人相关研究与不足
| 研究 | 主要结果 | 不足 |
|---|---|---|
| Pierson & Moskowitz (1964) | 提出成熟风海谱 PM 形，低频 exp(-1.25(f/fp)⁻⁴) | 无 >20 s 数据，低估低频能量 |
| JONSWAP (Hasselmann 1973) | 在 PM 谱上加峰增强 γ=3.3，描述成长风海 | 仅调峰区，f<0.8fp 沿用 PM，与四波相互作用理论不符 |
| Toba (1973) | 给出 Hs ∝ Tp^(3/2) 经验律 | 验证数据 Hs<10 m，极端条件下偏差未知 |
| Munk (1947) | 提出用先导涌浪反推风暴信息 | 仅有 3 个浮标，无法全球应用 |
| Cavaleri (2009) | 指出卫星高度计“错过峰” | 未提供解决办法，无周期信息 |
| Aster et al. (2023) | 发现 20 世纪晚期以来微震能量上升 | 缺乏同期海浪周期数据，无法归因 |
→ 共同不足：
- 无 >16 m 海浪的同步波高-周期实测；
- 低频端谱形经验外推，未用四波作用理论约束；
- 卫星一维高度计无法分辨长周期涌浪；
- 微震与海浪周期未建立定量联系。
本文数据与方法
| 类别 | 内容 |
|---|---|
| 核心数据 | 1. SWOT L2 250 m 二维海面高度（2023-04—2024-12）<br>2. 8 颗卫星高度计（Jason-3, Sentinel-3/6, Cryosat-2, CFOSAT…）<br>3. LOPS WAVEWATCH III 30 km 全球 hindcast（1991-2024） |
| 预处理 | 40 km×40 km 框 FFT → 5 km 瓦片平均 → 128 自由度谱；<br>方向-频率掩膜隔离单风暴涌浪分区 |
| 反演方法 | 1. 线性弥散关系拟合 LE(d) 得风暴中心时空校正；<br>2. 用 Hss(d) ∝ d⁻⁹ 与理论 E∝f⁻¹⁷ 联合最小二乘反推风暴峰周期 SPP 与等效半径 r；<br>3. 更新谱形：保持 f>fp 为 PM，f<fp 改为 f⁻¹⁷，平滑过渡 |
| 验证手段 | 对比同一时刻/地点模式 Hs 与高度计 Hs（r=0.98）；<br>用 Toba 律 Hs-Tp 关系独立校验 SPP |
主要结果
- 全球 100 个最强风暴中，SWOT 观测到 10 个涌浪场，反演 SPP 成功 6 个；
- 2024-12-21“Eddie”风暴：SPP=20.2±0.6 s，同期 SWOT 高度计测得 Hs=19.7±0.3 m，为 1991 以来卫星最高记录；
- 更新谱形使低频（0.5–0.8 fp）能量降为原 JONSWAP 的 1/20，与 SWOT 涌浪衰减一致；
- 发现涌浪高随距离衰减指数 n≈9，对应能量 ∝d⁻¹⁸，与四波作用逆级联理论吻合；
- 提供 2023-2024 年 6 个风暴 SPP 数据集（公开于 SEANOE）。
创新点与贡献
| 创新/贡献 | 说明 |
|---|---|
| 首次用二维卫星高度数据反推风暴峰周期 | 解决“卫星错过峰”难题，无需现场浮标 |
| 提出 f⁻¹⁷ 低频谱形更新 | 把四波作用理论嵌入工程标准谱，填补 PM/JONSWAP 空白 |
| 给出卫星时代最大 Hs-Tp 同步实测 | 19.7 m & 20.2 s，可作为极端设计新标准 |
| 建立涌浪衰减幂律模型 | 只需 d 与 SPP 即可预测远场长波高，服务海啸/微震研究 |
| 公开 34 年极端风暴目录 + SWOT SPP 数据集 | 方便后续气候趋势与工程统计分析 |
本文不足
- SWOT 带宽限制：250 m 空间分辨率→可解最短 18 s 周期，对 TC 更短周期风暴误差大；
- 未考虑海流-地形折射：部分北太平洋路径受岛链影响，衰减指数可能偏差；
- 未做现场浮标同步：缺乏 Hs>16 m 的锚系/漂流浮标“真值”校验；
- 等效风暴半径 r 物理意义模糊：未与风场半径或气压梯度挂钩；
- 微震-海浪转换函数未定量：仅给出周期一致，幅度关系待后续建立。
后续改进与跟进
| 方向 | 具体措施 |
|---|---|
| 提高分辨率 | 2026 年 SWOT+ Ka 波段 100 m 产品，可解 12–15 s 涌浪，扩展至 TC 短周期风暴 |
| 现场加密观测 | 在北大西洋/北太平洋布放极端事件快速浮标阵列（如“extreme-wave glider”），同步记录 Hs>18 m 与 Tp |
| 耦合海流模式 | 将 WAVEWATCH III 与全球 1/12° 海流模式耦合，重新计算折射-绕射，校正衰减指数 |
| 矩张量反演 | 把更新谱作为输入，计算 18–30 s 微震源区辐射强度，与全球地震噪声记录做联合反演 |
| 古风暴重建 | 将新 Hs-Tp 关系代入沉积物波痕模型，重新估算全新世飓风强度，验证气候变暖趋势 |
| 工程标准更新 | 向 ISO 19901-1、DNV-RP-C205 推荐采用“LOPS-2025”谱形，作为极端设计基准场景之一 |
| 长序列趋势 | 利用 2025-2030 年 SWOT+CFOSAT 数据，把 SPP 序列延长至 10 年，分析极端海浪能量是否继续攀升 |