请留言
slide1

全球视野 媒体聚焦

接受央视报道共计15次 党媒、央媒、地方官媒总传播量达0.5亿人次

0101
新闻分类

洪泽湖水体叶绿素a浓度反演

2026-04-18

概述

叶绿素 a 是浮游植物的主要光合色素,常被用作浮游植物生物量的指标,是水环境管理的重要组成部分。叶绿素 a 浓度(Chla)是评估水体富营养化程度和研究水生环境中碳循环的关键指标,其快速、准确获取对湖泊生态管理至关重要。

传统的湖泊营养状态评价方法主要依赖于对水体的实地取样分析,然而该方法易受局部天气与环境影响,取样与测试过程也比较耗时耗力、成本较高,较难实现对湖泊富营养状态在精细时空尺度上的监测。与之相比,遥感技术具有覆盖范围广、获取资料快、周期短等优点,可弥补传统水质采样的诸多不足,因此,已被广泛应用于湖泊水环境和水生态等方面的监测。

数据

数据源:

(1)采样点数据:在本次研究中,共收集并分析了19个有效采样点数据,这些数据点经过严格筛选和质量控制,确保后续分析的准确性和可靠性。

(2)遥感数据:采用了Sentinel-2多光谱仪(MSI)提供的Level-2A级标准产品作为主要的遥感数据源。该数据产品涵盖了10个关键的光谱波段,包括B1至B9波段以及B8A波段,能够全面捕捉地表水体的光谱特征。所有遥感数据均通过Google Earth Engine平台高效获取和处理。

数据处理:

(1)实测数据清洗:首先剔除叶绿素a浓度(Chl-a)小于或等于0以及数值缺失的数据点,以确保数据的有效性和准确性。

(2)水体提取:应用归一化水体指数(NDWI),计算公式为(B3 - B8) / (B3 + B8),通过设定阈值大于0.1生成水体掩膜,从而精确提取出水体区域,排除非水体像元的干扰。

(3)时空匹配处理:在空间上,提取每个采样点周围3x3像素窗口的平均反射率值,以减小像元级噪声和空间异质性的影响;在时间上,依据实际采样日期进行精确匹配,确保遥感数据与实地观测数据在时间上高度一致,提升分析的时空准确性。

技术路线

图1. 技术路线图

方法

在水质参数的反演工作中,获取与水质参数关系密切的敏感波段,将其作为模型的输入因子,这样建立的模型具有更高的预测精度。本次通过 Pearson 相关分析来获取与叶绿素 a 质量浓度关系密切的敏感波段,Pearson 相关系数是一种线性相关系数,通过 Pearson 相关分析可以得出不同波段或者波段组合与叶绿素 a 质量浓度的相关性强弱,可进一步剔除弱相关的、可能干扰反演模型建立的波段信息。由于水体的反射特性主要位于可见光和近红外波段,且有研究明,可见光和近红外的反射率可以成功地用以反演水体的叶绿素 a 质量浓度。

本次研究基于已有的实测数据和卫星影像数据,通过对比不同单波段以及波段组合的相关系数,分别以相关性最高的波段组合值作为自变量,水质参数叶绿素a浓度作为因变量,构建统计回归模型进行水质反演,并以决定系数R2为评估标准检验模型是否准确,最终得到反演模型。

结果与分析

基于上述反演模型,得到采样点当日洪泽湖水体叶绿素a浓度分布图:

图2.洪泽湖水体叶绿素a浓度反演(mg/L)

如图所示,基于Sentinel-2反演的洪泽水体叶绿素a(Chl-a)浓度全湖整体处于较低营养水平,Chl-a浓度最大值仅为0.0145 mg/L,符合冬季水体浮游植物生物量受低温与弱光照抑制的季节性特征。空间上呈现明显的非均质格局,表现为“湖心低、边缘高”的分布模式:湖心开阔水域Chl-a浓度普遍偏低,水体较为清洁;而沿岸区域,尤其是入湖河口及邻近城镇的浅水湾汊,出现局部高值斑块,该处人类活动频繁,存在对水体富营养化的潜在影响。总体来看,此次反演时间为1 月(冬季),洪泽湖水温较低,藻类生长速率整体受抑制,因此主湖区 Chla 浓度处于全年较低水平;但湖岸带受陆源输入的 “点源式” 影响,仍形成局部高值区,这与洪泽湖 Chla 浓度 “冬春季近岸高、核心区低” 的季节分布规律一致。