东海有害藻华多源卫星遥感监测能力评估：GOCI-II、Sentinel-2 MSI与Sentinel-3 OLCI跨传感器比较

《Frontiers in Marine Science》：Multi-source satellite remote sensing of harmful algal blooms in the East China Sea: cross-sensor comparison of GOCI, Sentinel-2 and Sentinel-3

　　

引言：有害藻华（HABs）是反复发生的近岸生态灾害，威胁水产养殖、海洋生态系统和近岸环境安全。卫星遥感为HAB监测提供必要观测能力；然而，针对同一藻华事件，地球静止轨道GOCI-II与极轨Sentinel-2 MSI、Sentinel-3 OLCI之间的系统定量互比较仍有限。本研究探讨三种海洋水色传感器在东海HAB观测中的互补优势与局限。方法：研究人员选取2021年东海两次有充分记录的HAB事件（4月30日和9月22日）。GOCI-II、Sentinel-2 MSI和Sentinel-3 OLCI数据集采用传感器匹配的大气校正方案（GOCI-II标准海洋水色处理、Sentinel-2采用ACOLITE-DSF（暗光谱拟合，DSF）、Sentinel-3 OLCI采用C2RCC（Case-2 Regional Coast Colour））。应用NDCI（归一化差异叶绿素指数）、RBD（红边波段差）、NRTI（归一化赤潮指数）、RAB（赤潮藻华指数）、PDI（东海原甲藻指数）和DI（硅藻指数）6个光谱指数进行藻华检测。实施20 m、300 m和500 m空间退化实验，以区分几何尺度相关效应对藻华面积和边界表征的影响。开展近同步过境比较和GOCI-II逐小时时间序列追踪，以量化传感器间差异和日内藻华动态。结果：在2021年4月30日约10 min近重合过境条件下，GOCI-II制图藻华面积为901.52 km²，Sentinel-2B为571.52 km²；差异共同源于空间分辨率、光谱配置、大气校正和检测阈值，而藻华平流影响较小。2021年9月22日，GOCI-II逐小时观测捕获完整日内藻华循环：藻华光学异常面积从37.30 km²变化至峰值760.20 km²，再降至107.67 km²，伴随45.74 km质心位移；单次极轨快照仅部分捕获该动态。空间退化分析表明，粗分辨率下边界长度损失超过76–85%，总面积偏差为?14.8%至+11.0%，突显Sentinel-2在近岸精细边界刻画中的独特价值。即使进行分辨率归一化，传感器藻华产品间仍存在显著空间错配（交并比，IoU低至0.014）。讨论：单一传感器呈现不同权衡：GOCI-II提供高频时间筛查，Sentinel-2提供精细空间细节，Sentinel-3 OLCI提供丰富光谱信息用于区域背景分析。据此，研究人员提出条件式三阶段多源协同工作流：GOCI-II用于逐小时快速筛查和漂移追踪，Sentinel-2 MSI用于优先区精细检查，Sentinel-3 OLCI用于区域光谱验证和跨传感器一致性检查。传感器间差异源于光谱波段配置、大气校正、指数阈值和采样时间的综合影响。局限性包括仅分析两次藻华事件且缺乏密集像元级现场真值。需要开展更多事件验证并协调大气校正，以将该条件式监测模板推广至光学复杂东亚近岸水体的业务化HAB预警。

研究背景方面，东海是中国HAB高发海域，藻华多发生于4—9月，主要致灾种包括Prorocentrum donghaiense、Karenia mikimotoi和Noctiluca scintillans；该海域受高悬浮泥沙、有色溶解有机物（CDOM）和强河口—近岸光学梯度影响，传统现场监测在空间覆盖、时间频率和响应时效上受限，而单传感器难以兼顾时间连续性、空间细节和光谱诊断能力。已有GOCI日变化监测及Sentinel-2/3跨传感器研究，但针对东海同一藻华事件、近同步观测下面积差异量级及粗分辨率对边界刻画相对面积估计的退化程度仍缺乏系统认识，因此研究人员开展GOCI-II、Sentinel-2 MSI与Sentinel-3 OLCI的跨传感器比较，并提出条件式多源协同流程，其意义在于为光学复杂近岸水体HAB快速筛查、精细核查和区域一致性检查提供可检验模板。

关键技术方法不超过250字：研究人员使用2021年东海两次HAB事件：4月30日GOCI-II与Sentinel-2B近同步（约10 min），9月22日8幅GOCI-II小时场景并配Sentinel-2 MSI和Sentinel-3 OLCI；参考信息包括官方HAB通报和5个地理参考现场记录（Supplementary Table S1），9月22日有温州海洋环境监测中心业务记录佐证A. sanguinea藻华。GOCI-II L2海洋水色产品为250 m，Sentinel-2 MSI L1C经ACOLITE-DSF校正并重采样20 m，Sentinel-3 OLCI L1B经C2RCC处理为300 m；采用NDCI、RBD、NRTI、RAB、PDI、DI及动态/固定阈值检测；开展二值掩膜和反射率域300/500 m退化实验、0.003°共同网格IoU比较、质心Haversine位移与bootstrap估计。

结果部分保留原小标题并浓缩如下。3.1 2021年4月30日藻华事件：通过近同步GOCI-II与Sentinel-2B比较，GOCI-II以RAB+PDI检测901.52 km²（有效水体8.43%），Sentinel-2B以NDCI>0.1检测571.52 km²（4.03%）；约1.6倍面积差异不能归因于单一传感器属性，GOCI-II呈连续沿岸带、Sentinel-2B呈碎片化精细结构，该比较为不同测量响应而非精度排序。3.2 2021年9月22日藻华事件：通过8幅GOCI-II小时场景及OLCI、MSI近时过境比较，10:15 GOCI-II为375.69 km²、10:18 OLCI为606.19 km²、10:25 MSI为156.04 km²；GOCI-II日内面积由37.30 km²升至12:15峰值760.20 km²再降至107.67 km²，说明单次极轨快照可能采样非代表性相位，但不能直接解释为生物量增长。3.3 多源传感器比较：通过不同指数、阈值、大气校正和分辨率产品比较，发现其不能视为真实藻华面积分数的直接可比估计；光谱上藻华水在705 nm具红边峰（R_rs 0.027–0.032 sr^-1，背景0.002–0.005 sr^-1），443和665 nm吸收加深；GOCI-II 660/680 nm跨叶绿素吸收谷与红边，Sentinel-2 665/705 nm适配NDCI，Sentinel-3 OLCI 681 nm唯一支持RBD，555/660 nm支持NRTI。3.4 空间分辨率效应：通过Patch A和Patch B的掩膜尺度与反射率域退化实验，Patch A由20 m的1.394 km²、边界46.72 km降至300 m反射率域1.281 km²（?8.1%）、边界11.12 km（?76.2%），500 m为1.188 km²（?14.8%）、6.82 km（?85.4%）；Patch B由2.928 km²、115.84 km变为300 m的3.060 km²（+4.5%）、22.80 km（?80.3%），500 m为3.250 km²（+11.0%）、20.00 km（?82.7%），表明分辨率降低对边界长度影响强于面积；共同网格MSI–OLCI IoU=0.131、面积比0.179，MSI–GOCI-II IoU=0.014、面积比0.316，说明分辨率匹配不能消除差异。3.5 时间序列分析与漂移追踪：通过GOCI-II小时质心轨迹，08:15—15:15 BJT质心移动45.74±3.2 km（bootstrap），支持快速筛查和后续观测调度，但轨迹可能包含生长、平流和潮汐效应，不能单独归因。

讨论部分总结：GOCI-II的核心优势是逐小时时间采样，可提供筛查、漂移解释和后续观测调度的时间背景；Sentinel-2在近岸边界精细刻画方面价值突出，但退化实验不能复现原生粗分辨率传感器的完整点扩散函数和大气校正差异；Sentinel-3 OLCI以21个海洋水色波段、约300 m采样和681 nm波段提供区域光谱背景及跨传感器一致性检查。阈值敏感性方面，9月22日MSI固定NDCI>0.1面积为156.04 km²，动态阈值T_dyn=0.02993时增至508.85 km²（+226.1%），但仍比同期OLCI的606.19 km²小97.34 km²（16.1%），说明阈值是重要不确定性来源而非已验证最优值。误差来源包括悬浮泥沙、CDOM、浅水底反射、物种色素与垂直分布差异、大气校正链差异、MSI L2A缺乏逐像元不确定层及时间错配；仅两次事件、无密集像元级真值，故因素归因保持描述性。协同框架为条件式三阶段：GOCI-II小时筛查与漂移追踪、Sentinel-2 MSI优先区精细检查、Sentinel-3 OLCI区域光谱背景与一致性检查，属可检验模板而非普适最优处方。研究结论翻译：本研究对2021年东海两次HAB事件开展GOCI-II、Sentinel-2 MSI和Sentinel-3 OLCI系统跨传感器比较，主要结论为：三种传感器分别具有时间连续性、空间细节和光谱诊断的互补能力；4月30日近同步比较显示GOCI-II与Sentinel-2B面积显著不同（901.52 vs. 571.52 km²），但反映分辨率、光谱配置、大气校正和阈值选择的综合效应而非单一传感器效应；9月22日时间对齐观测为GOCI-II 10:15 375.69 km²、OLCI 10:18 606.19 km²、MSI 10:25 156.04 km²，MSI动态阈值T_dyn=0.02993使面积增至508.85 km²但未消除与OLCI差异；GOCI-II小时观测解析37.30→760.20→107.67 km²面积变化和45.74 km质心位移，但在未独立控制光照和大气校正效应时不应解释为直接生物量增长；反射率域退化显示边界损失76.2–85.4%（Patch A）和80.3–82.7%（Patch B），面积变化较小（?14.8%至+11.0%）；匹配分辨率比较MSI–OLCI IoU=0.131、面积比0.179，MSI–GOCI-II IoU=0.014、面积比0.316，说明仅匹配分辨率不能消除跨传感器差异。由于仅含两次事件且无MSI不确定层，上述因素归因仍为描述性，需更多事件、协调大气校正、不确定传播和独立验证。条件式三阶段工作流可作为多传感器任务划分的可检验工作假设，而非普适优越性或已验证业务最优性的断言。