《Remote Sensing of Environment》：Fusing SAR and SWIM observations to improve the global stokes drift estimates
编辑推荐：
斯托克斯漂移（Stokes drift）是上层海洋示踪物输运的关键参数，也是波流耦合中的重要物理量。本研究提出一种融合框架，将合成孔径雷达（Synthetic Aperture Radar, SAR）与海浪探测与监测仪（Surface Waves Invest
斯托克斯漂移（Stokes drift）是上层海洋示踪物输运的关键参数，也是波流耦合中的重要物理量。本研究提出一种融合框架，将合成孔径雷达（Synthetic Aperture Radar, SAR）与海浪探测与监测仪（Surface Waves Investigation and Monitoring, SWIM）的观测数据进行协同融合，以改进全球斯托克斯漂移的估算精度。研究利用SAR提供的波向信息解决SWIM波谱存在的方向模糊问题，同时保留SWIM在波幅测量上的准确性。通过浮标观测数据的验证表明，该融合方法在不同海况条件下较单一传感器具有更高的估算精度。研究进一步讨论了观测误差与方法误差的主要来源，包括SAR方位截止效应、SWIM方向模糊以及风浪与涌浪方向不一致等因素。基于融合框架估算的全球斯托克斯漂移季节变化特征显示出良好的应用前景，为海浪谱高阶矩参数的反演提供了新的技术途径。
论文解读
一、研究背景与科学问题
斯托克斯漂移是海洋学中的基础物理现象，源于海浪轨道运动随深度衰减的特性，表现为沿波浪传播方向的净漂移。该现象在海洋上层示踪物输运中扮演核心角色，约占开阔大洋风驱漂移的三分之二，并可通过与流场相互作用产生朗缪尔环流（Langmuir circulations），显著影响上层海洋的垂直混合及生物地球化学示踪物的垂向输运。因此，斯托克斯漂移的准确估算对海洋生态系统和气候研究具有重要意义。
然而，斯托克斯漂移的直接测量极为困难。现有估算途径主要包括：高频雷达（HF radar）反演、基于实测海浪谱的谱矩计算、经验参数化方法以及海浪数值模式模拟。高频雷达方法在从多普勒频移中分离斯托克斯漂移仍存在争议；基于浮标实测谱的方法受限于浮标站点稀疏且集中于近岸区域，难以实现全球覆盖；经验方法虽能较好估算风浪贡献，但对涌浪（swell）成分的估算存在显著困难，因为涌浪反映的是远处历史风况的影响；海浪模式如WAVEWATCH III（WW3）虽可提供全球估算，但计算成本高昂，且缺乏足够的观测数据进行验证和参数校准。
遥感技术为全球海浪谱的连续观测提供了可能。SWIM搭载于中法海洋卫星（CFOSAT），具有较宽的波数测量范围和较小的幅值低估；SAR广泛搭载于Sentinel、ERS等卫星，在波向分辨方面具有优势。然而，两种传感器均存在固有局限：SAR的方位截止效应（azimuth cutoff effect）限制了其可测波数范围并造成波谱畸变；SWIM则存在相反方向的波谱模糊问题。此外，沿飞行轨迹的斑点噪声和调制传递函数效应也会扭曲波谱形状。这些因素对有效波高等低阶矩参数影响不大，但对高阶谱矩及斯托克斯漂移的估算可能引入显著偏差。
现有研究曾通过假设波浪沿风向传播来消除SWIM的方向模糊，但该方法忽略了涌浪导致的风浪方向不一致问题，可能同时影响估算的幅值和方向精度。因此，亟需探索一种融合方法，结合SAR的波向优势和SWIM的波幅优势，以改进斯托克斯漂移的全球估算。
二、研究内容与主要结论
本研究针对SAR与SWIM各自的局限性，提出协同融合框架。研究通过将SWIM波谱划分为不同波浪分量，利用SAR观测进行涌浪分量匹配以确定其传播方向，同时采用波龄方法（wave age method）结合SWIM L2产品中集成的ECMWF风场数据确定风浪方向，从而解决SWIM的方向模糊问题。研究在三个代表不同海况特征的浮标站点进行验证：海洋站Papa（Ocean Station Papa）具有显著季节变化，NDBC 51002站（夏威夷东南部）以涌浪为主，南方海洋通量站（Southern Ocean Flux Station, SOFS）以风浪为主。验证结果表明，融合方法在不同海况条件下均较单一传感器具有更高的斯托克斯漂移估算精度。研究还讨论了观测误差和方法误差的三类主要来源，并将融合框架应用于全球海洋斯托克斯漂移估算，初步结果显示其季节变化特征具有良好的应用前景。该研究发表在《Remote Sensing of Environment》。
三、关键技术方法
研究人员采用的主要技术方法包括：其一，SAR交叉谱分析技术，基于Hasselmann和Hasselmann（1991）建立的SAR图像谱与海浪谱之间的闭合映射关系，利用Sentinel-1 L2 OCN产品提取波浪传播方向信息；其二，SWIM波谱反演技术，利用CFOSAT搭载的SWIM仪器获取二维海浪谱，其波数覆盖范围优于SAR；其三，波谱分量划分与方向解模糊算法，将SWIM波谱按峰值波数升序划分为不同波浪分量，通过SAR分量匹配确定涌浪方向，通过波龄方法结合ECMWF风场数据确定风浪方向；其四，统一的波数域截断计算方案，为保证不同波谱产品间的可比性，研究采用SWIM的一致波数域计算斯托克斯漂移，即"可解析斯托克斯漂移"（resolved Stokes drift），以区别于全波数域积分得到的完整斯托克斯漂移。
四、研究结果
4.1 融合框架的构建与方向解模糊
研究人员通过将SWIM波谱划分为多个波浪分量，利用SAR观测对涌浪分量进行匹配以确定其传播方向，同时通过波龄方法利用ECMWF风场数据确定风浪方向。在假设条件下，第一和第三波浪分量的方向模糊得到解决，而第二分量因超出SAR可测域而未能解模糊。这一结果表明融合框架能够有效利用SAR的波向信息和SWIM的波幅信息，获得优于任一单一传感器的海浪谱和斯托克斯漂移估算。
4.2 浮标验证结果
研究人员在三个浮标站点对融合估算的斯托克斯漂移进行了系统验证。在海洋站Papa，融合方法在季节变化显著的条件下表现出良好的追踪能力；在NDBC 51002站，融合方法在涌浪主导条件下显著改善了方向估算精度；在SOFS站，融合方法在风浪主导条件下保持了较高的幅值精度。验证结果表明，融合方法在不同海况条件下均较单一传感器具有更高的估算精度。
4.3 误差来源分析
研究人员通过典型实例分析了三类误差来源：SAR方位截止效应导致的高波数分量缺失及其对高阶谱矩的影响；SWIM方向模糊未完全消除时对波向估算的偏差；风浪与涌浪方向不一致时对斯托克斯漂移方向和幅值的综合影响。这些误差分析为后续算法改进提供了明确方向。
4.4 全球斯托克斯漂移的季节变化
研究人员将融合框架应用于全球海洋斯托克斯漂移估算，初步结果显示其季节变化特征与海浪气候学的已知规律一致，表明该方法具有全球应用的潜力。
五、讨论与结论
讨论部分指出，融合框架不仅解决了斯托克斯漂移估算的即时问题，还为全球二维波谱重建提供了新框架。现有研究多基于低阶谱矩参数（如有效波高、峰值波向）验证波谱质量，而斯托克斯漂移依赖于高阶谱矩，对波能量、方向和整体谱形高度敏感，因此本研究为波谱测量评估提供了更严格和全面的基准。
研究结论可归纳为：SAR与SWIM的协同融合能够有效克服单一传感器的局限性，显著提升斯托克斯漂移的估算精度；融合方法在多种海况条件下均表现良好，具有全球应用的可行性；研究识别的三类误差来源为后续改进提供了方向；全球斯托克斯漂移的季节变化估算显示出有前景的结果。未来研究可进一步拓展融合框架的应用范围，并针对二维波谱的其他方面开展补充验证。
