利用陆基高光谱遥感捕捉太湖蓝藻水华日内快速变化过程*

张运林，张毅博，李 娜,2，孙 晓,2，王玮佳,2，秦伯强,3，朱广伟,3

(1:中国科学院南京地理与湖泊研究所湖泊与环境国家重点实验室，南京 210008)(2:中国科学院大学，北京 100049)(3:南京中科深瞳科技研究院有限公司，南京 210031)

过去40年，我国经济社会快速发展的同时也遭遇持续不断的污染困扰，人类活动导致的大量氮磷污染物排放进入河流和湖库，引发长江中下游及东部沿海地区许多湖库都面临富营养化与蓝藻水华问题，严重威胁湖库生态系统健康和饮用水安全[1-3]. 而全球气候变化与营养盐富集的协同放大效应进一步加剧了湖库富营养化及蓝藻水华暴发程度，导致近年来通过大规模流域控源截污以防止湖库藻类水华也难奏效[4-5]，使得湖库富营养化治理与藻类水华控制成为一项长期而又艰巨的任务. 2007年太湖蓝藻水华事件曾造成无锡自来水被污染，数百万居民自来水供应中断一周[2]. 之后的十多年高强度治理显示，太湖总氮浓度持续下降，但作为浮游植物生长和蓝藻水华发生关键限制因子的总磷浓度则先降后升，并且在2017-2019年居高不下、高位波动，蓝藻水华频繁暴发甚至发生小规模“湖泛”事件的威胁依然存在[6-8]. 据中国科学院太湖湖泊生态系统研究站(简称太湖站)全湖32个站点季度长期定位观测数据显示，2019年太湖水体总氮浓度已由2007年的3.37 mg/L降到2.03 mg/L，但总磷浓度只由0.122 mg/L变为0.118 mg/L. 相伴随的是2017年发生了有记录以来最严重的蓝藻水华，最大水华面积达1582 km2[4]，2018年虽略有减轻，但2019年的蓝藻水华强度依然较高，强度仅次于2017年. 因此，太湖生态环境保护与治理仍面临较长时期内的控源截污、蓝藻水华防控和“两个确保”压力(确保饮用水安全、确保不发生大面积湖泛).

太湖蓝藻水华在风场和水动力作用下呈现高度的时空异质性和快速动态变化，在小时尺度上会发生明显水平迁移和垂直混合[9-12]，传统的湖面定位和断面监测总感觉蓝藻水华有时“来无影、去无踪”[13]. 蓝藻水华在湖面上的漂浮不定给藻华和湖泛等异味水体的防控和治理带来了很大挑战，也威胁到饮用水源地安全. 多源卫星遥感可以实现蓝藻水华空间分布同步观测，但受制于卫星传感器时间、空间和光谱分辨率，加之云雨天气的影响，很难捕捉蓝藻水华快速动态变化过程. 水下探头高频监测为水质快速动态变化过程观测提供了新的途径[14-16]，但太湖湖面开阔、风浪扰动强烈、水体浑浊，水下探头很容易被附着和污染，也容易被风浪损坏，受水况环境干扰造成监测精度不稳定，可监测指标有限，使用受到极大限制，并且仪器价格高昂、后期管理和维护成本非常高[17-19]. 与此同时，2015年国务院发布的《生态环境监测网络建设方案》提出：“强化高新技术、先进装备与系统的应用，提高生态环境监测立体化、自动化、智能化水平”. 因此，迫切需要发展新的自动化和智能化的水环境监测设备和技术方法，实现蓝藻水华及关键水质参数快速动态变化连续高频监测，应对水环境突然变化，补充和完善现有监测体系不足，服务蓝藻水华防控和水环境监管.

本文利用前期自主研发的陆基高光谱多参数水质遥感监测仪，通过秒-分钟频次浮游植物叶绿素a、透明度、总氮、总磷和高锰酸盐指数的连续观测，试图精细刻画蓝藻水华及其关键水质参数快速变化过程，结合逐时风速风向等气象数据分析其快速变化原因，明晰蓝藻水华漂浮的临界风速，量化外部漂浮输入和本地生长对藻类生物量的贡献，有助于深入理解太湖蓝藻水华形成过程与驱动机制，破解蓝藻水华“来无影、去无踪”奥秘，支撑蓝藻水华精准打捞和快速清除等主动防控策略，实现更高水平“两个确保”.

1 数据与方法

1.1 现场观测

2018年我们率先提出近地面离水高光谱多参数水质陆基(地基、岸基)遥感监测理念和实施路径，南京中科深瞳科技研究院有限公司联合杭州海康威视数字技术股份有限公司首创近水面非接触式陆基(地基、岸基或者平台、船舶、桩基等能固定安装的均可)高光谱多参数水质遥感监测仪(ZKDI-R2021型)，仪器集成高光谱采集、视频监控、毫米波雷达水位测定、水质参数反演和深度学习等技术，可以实现复杂天况下总氮、总磷、高锰酸盐指数、悬浮物、透明度、浊度、消光吸收、有色可溶性有机物吸收、叶绿素、藻蓝素、藻密度共11个关键水质参数实时高频监测. 陆基高光谱多参数水质遥感监测仪既克服了以往水下探头的易污染、难维护、参数少和精度低的问题，又解决了卫星遥感监测中时间和光谱分辨率低以及受天气和大气尘埃干扰的影响，实现遥感和断面监测的深度融合.陆基高光谱多参数水质遥感监测仪架设在中国科学院太湖湖泊生态系统研究站(简称太湖站)水上观测场栈桥顶端，离岸边约400 m，观测站点东面紧靠岸边，北、西、南三面为开敞水域，特别是西面没有任何岸边阻断(图1). 仪器距水面约4.0 m(架设高度3～10 m均可)(图1)，光谱波段为400～1000 nm，光谱分辨率为1 nm，观测频次为20 s/次. 仪器通过上光谱相机和下光谱相机同步测量入射光和水面反射光进而计算得到辐照度比. 仪器内已植入前期开发的11个关键水质参数深度学习反演算法，算法基于太湖、富春江水库、千岛湖、梁溪河和南京城市河流等近700个样品构建，可同步输出水体高光谱辐照度比和水质参数实时数据[20]①. 对于400～1000 nm可见光和近红外光学遥感来说，悬浮物、透明度、叶绿素和有色可溶性有机物等光学物质有明显的光学信号，高光谱遥感反演起来比较容易，而总氮、总磷的吸收光谱在紫外波段，其在可见光和近红外确实没有光学信号，但总氮、总磷等非光学物质与悬浮物、透明度、叶绿素和有色可溶性有机物等光学物质具有复杂的耦合关系，因此在反射光谱上也会有响应，只是这种响应不是线性或者一对一的. 目前仪器植入的是神经网络模型、高斯回归过程等机器学习算法，从400～900 nm以每5 nm或者10 nm光谱辐照度比作为输入数据进行模型训练和检验，由于算法包含的光谱通道比较多，可以很好反映总氮、总磷等非光学物质与光谱反射率之间的复杂耦合关系，同时算法训练数据集和检验数据集覆盖不同类型水体，样本量也比较大，相关参数精度都在80%以上[20-22](1)陆基高光谱多参数水质遥感算法2021.05版..

图1 观测点位置(a)、陆基高光谱多参数水质遥感监测仪现场架设照片(b)和藻华现场图片(c)Fig.1 Location of observation site (a), photo of in situ land-based hyperspectral water quality remote sensing monitor instrument (b), and algal bloom photo (c)

1.2 气象水文数据

太湖站架设有气象观测站，实时记录逐时气象观测数据，主要包括气温(日最低、日最高、日平均)、风速风向(平均风速、最大风速)和降水等. 水温数据来自于太湖站栈桥顶端水下探头连续高频监测，监测频率是0.5 h/次.

1.3 蓝藻水华界定

水华界定和对应的生物量阈值判定一直都是个难题[23]，有研究用叶绿素a浓度大于10 μg/L或藻细胞超过1.5×107cells/L，并在水面形成一层肉眼可见的漂浮藻类作为水华判别阈值[24-25]；美国佛罗里达州将湖泊水体叶绿素a浓度大于40 μg/L定义为藻类水华[26]. 根据我们在太湖开展的蓝藻水华卫星遥感监测和同步叶绿素a浓度测定发现，当叶绿素a浓度在30～40 μg/L之间时卫星影像就能观测到明显表面水华，因此以叶绿素a浓度30～40 μg/L作为太湖蓝藻水华判别阈值.

1.4 统计分析

采用SPSS 20软件对数据进行统计分析，包括计算其平均值、线性拟合等，当显著性水平值P≤0.05时表明显著相关，P≤0.01时表明极显著相关. 利用OriginPro 2018绘制文中数据图.

2 结果与讨论

2.1 气象水文条件

从观测期间气象水文来看，8月7-9日没有明显降水过程，7号有零星小雨，当天累积降水量3.6 mm，天气整体晴朗少云，气温较高，经历高温过程，逐时平均气温在24.4～36.6℃. 逐时平均水温在29.6～31.3℃，变化幅度显著小于气温日内变化，受气温日变化影响较小，逐时气温与水温不存在显著性相关(R2=0.04，P=0.14)，这与以往逐日观测结果明显不一样[27]. 整体而言，风速不大，7日逐时平均风速为2.66 m/s，8日更是降至1.60 m/s，属于小风和微风环境；7日风向以东风偏北为主，由太湖站岸边往湖面开敞水域吹，8日风向则以北风和西北风为主，由湖面开敞水域往太湖站岸边吹(图1、图2).

图2 观测站点2021年8月7-9日气温、水温(a)以及风速、风向(b)逐时变化Fig.2 Hourly variations of air temperature and water temperature (a), wind speed and wind direction (b) of observation site from August 7-9, 2021

2.2 叶绿素a快速变化过程

8月8日8:30-18:30期间总共观测了1800组叶绿素a数据，变化范围在8.0～115.0 μg/L，平均值为(59.5±22.63) μg/L. 以叶绿素a浓度30～40 μg/L作为太湖蓝藻水华判别阈值，从图3a叶绿素a高频动态变化过程来看，全天经历了明显的水华快速变化过程. 早上8:30-9:00期间叶绿素a浓度非常低，在10 μg/L波动，没有观测到明显水华；在9:20左右出现一个短暂峰值，叶绿素a浓度在30～35 μg/L之间，出现零星水华；9:30-10:30维持在低值，叶绿素a浓度略高于10 μg/L. 从10:30之后叶绿素a浓度快速拉升，在短短半小时内迅速增加到100 μg/L以上，增加了10倍左右，水面漂浮大量藻颗粒，形成明显的蓝藻水华(图1)，清晰彰显了短期内蓝藻水华“来无影”快速变化过程. 之后尽管叶绿素a浓度有所降低，但到18:30时一直都维持在高值，整体在60～80 μg/L波动，期间13:30也有一个叶绿素a浓度在100 μg/L左右的明显峰值. 此外，在15:40、16:20和17:50还存在多个叶绿素a浓度低峰，大约在80 μg/L左右.

从叶绿素a的高频变化过程来看，在短短的10小时之内，观测站点经历一次从无到有的水华快速变化过程，特别是10:30-11:00期间叶绿素a浓度迅速从10 μg/L左右上升到100 μg/L以上. 在这么短的时间内藻类自身的生长繁殖不可能支撑叶绿素a浓度如此快速的变化过程，只有可能是周边水域藻类水华漂浮集聚到观测站点. 综合分析观测站点位置和观测期间风速风向变化，很容易解释叶绿素a浓度短时间内的快速变化过程. 观测站点东面紧靠岸边，北、西、南三面为开敞水域，特别是西面没有任何岸边阻断(图1). 8月8日8:00-19:00期间观测站点经历微风和小风环境，风速非常低，在0.3～2.3 m/s之间，平均风速为1.70 m/s，这个风速非常有利于蓝藻水华在水面的漂浮集聚，以往研究表明当风速低于3～4 m/s时蓝藻很容易在水面漂浮集聚[9,28]，8:00-16:00盛行西北风，平均风向为325°，恰好将梅梁湾湖面开敞水域表面蓝藻水华吹向观测站点，从而形成明显的水面藻华过程(图1、图3a).

过去关于太湖蓝藻水华形成和快速变化究竟是本区域藻类迅速生长累积形成还是由于风速水动力扰动堆积而成存在很多的争议，本文从秒-分钟高频观测很好地证明了小时尺度上表面蓝藻水华快速形成主要受制于风力形成的水动力作用，水动力引起的蓝藻水华短期迁移堆积可能要远远大于生物的生长累积，半小时内叶绿素a浓度由10 μg/L左右上升到100 μg/L以上，说明90%以上藻类生物量累积是由于风浪对蓝藻水华的迁移堆积所致. 这些研究结果进一步深化了以往基于MODIS逐日和GOCI逐时遥感监测得到的蓝藻水华动态变化过程认识[9-11,29]，量化外部漂浮输入和本地生长对藻类生物量的贡献. 除了遥感监测外，自2007年太湖饮用水危机事件后，为快速捕捉蓝藻水华动态变化过程，国家和地方监测部门以及高校科研院所使用美国YSI多参数水质仪在太湖布设了许多高频自动监测系统，但由于太湖水体浑浊、水下探头很容易被附着和污染，造成YSI多参数水质仪监测的叶绿素a浓度存在很大出入和不确定性，因此鲜少见到基于水下探头分析太湖蓝藻水华和水质动态变化的研究结果[10]. 而基于陆基遥感高频监测能非常真实、客观和准确地反映叶绿素a、藻类水华和关键水质参数的快速动态变化，精细刻画传统的季度、逐月和逐周观测无法捕捉的水环境快速高频变化过程.

2.3 其他水质参数变化过程

由于藻类颗粒物对光的吸收、散射和衰减，随着藻类生物量增加透明度呈现相反的变化趋势. 整体来看，从8:30-10:30当叶绿素a浓度较低时，水体透明度较高，在0.4 m左右，当叶绿素a浓度增加到60 μg/L左右时，透明度大致维持在0.3 m左右(图3b). 由于太湖水体浑浊、除藻类颗粒物外还有大量的非色素无机颗粒物，因此水体透明度都不高. 受藻类生物量影响，总氮、总磷浓度等呈现出与叶绿素a浓度较为一致的变化趋势(图3c、3d、3e)，10:30-11:00期间也呈现快速增加趋势，分别由1.75和0.05 mg/L增加到3.5和0.175 mg/L，增加了100%和250%，之后维持高值动态波动. 相比于叶绿素a浓度短期内增加了900%，总氮、总磷的峰值和增加没有那么明显，反映除了藻类生物量引起的氮磷变化外，总氮、总磷中还包含许多其他颗粒态和溶解态氮磷成分. 高锰酸盐指数随着叶绿素a浓度变化会发生一些变化，但不如氮磷营养盐那么明显，当叶绿素a浓度较低时高锰酸盐指数在4.5 mg/L左右，叶绿素a浓度明显增加后其值维持在6.0 mg/L左右，增加30%左右，说明高锰酸盐指数受藻类生物量的影响相对较小.

图3 2021年8月8日8:30—18:30叶绿素a浓度(a)、透明度(b)、总氮浓度(c)、总磷浓度(d)和高锰酸盐指数(e)高频变化Fig.3 High-frequency observations of chlorophyll-a concentration (a), secchi disc depth (b), total nitrogen concentration (c), total phosphorus concentration (d)and chemical oxygen demand (e) from 8:30 to 18:30 on August 8, 2021

2.4 叶绿素a与其它水质参数间的关系

为定量分析透明度、总氮、总磷和高锰酸盐指数随叶绿素a浓度的变化情况，进行简单线性相关分析. 从图4结果可知，透明度与叶绿素a浓度存在极显著负相关，总氮、总磷、高锰酸盐指数则均与叶绿素a浓度存在显著正相关. 从线性相关的决定系数来看，叶绿素a浓度可以解释总氮、总磷80%以上变化，而只能解释透明度、高锰酸盐指数35%左右的变化，说明无机悬浮颗粒物、溶解性有机物和其他耗氧有机物控制了透明度和高锰酸盐指数的动态变化等. 尽管氮磷营养盐是水体浮游植物生长及湖泊营养状态的重要因素和关键驱动因子，但藻类的大量生长反过来也会造成颗粒态有机氮磷和总氮、总磷的显著增加. 太湖站2005-2020年长期定位观测数据显示，太湖水体叶绿素a浓度与总磷浓度呈极显著幂函数关系，水体叶绿素a浓度越高，水体总磷浓度则越高，表明水体浮游植物生物量高低是水体总磷浓度高低的重要影响因素[8,10,30]. 由于硝化和反硝化作用影响，野外全年逐月监测往往很难观测到叶绿素a与总氮存在显著正相关，有时反而为负相关，但春夏季或者短期内藻类的快速聚集无疑会显著增加总氮浓度[10,30-31]. 而藻类生物量聚集也会明显增加颗粒态有机物，其死亡降解会消耗大量氧气，造成耗氧有机物指标高锰酸盐指数增加. 以往在太湖的研究也发现夏季蓝藻水华死亡腐烂降解释放造成高锰酸盐指数与叶绿素a浓度存在一致的空间分布，在叶绿素a非常高的站点同样会出现高锰酸盐指数异常高值，叶绿素a与高锰酸盐指数存在显著正相关[32-33].

图4 叶绿素a浓度与透明度(a)、总氮浓度(b)、总磷浓度(c)和高锰酸盐指数(d)间的线性关系Fig.4 Linear relationships between chlorophyll-a concentration and secchi disc depth (a), total nitrogen concentration (b), total phosphorus concentration (c), and chemical oxygen demand (d)

3 结论与展望

本文基于自主研发的陆基高光谱多参数水质遥感监测仪，通过秒-分钟频次高频监测成功捕捉了蓝藻水华及其关键水质参数在一天内的短期快速变化过程，并结合气象水文数据分析其形成原因. 研究发现在微风和小风条件下蓝藻容易在表层水体漂浮，盛行西北风驱动湖面开敞水域蓝藻水华在太湖站栈桥顶端观测点集聚，造成叶绿素a浓度在短短半小时内由10 μg/L快速攀升至100 μg/L左右，之后整体在60～80 μg/L波动；相伴随的是透明度明显降低，总氮、总磷和高锰酸盐指数显著增加，说明藻类水华可以引发关键水质参数快速变化.

与此同时，我们研究结果也显示陆基遥感高频观测可以揭示水环境和生态系统精细化过程和突变情况，深化湖泊生态系统动力学认知，为藻类水华和水质水环境监测、诊断提供强大工具，也可以提升未来水环境演化预测预警能力[19]. 实际上除了藻类快速变化外[34]，受极端降水、洪水、台风等自然和人类活动强烈作用，湖泊生态系统物理、化学和生物过程还遭受许多其他短期快速冲击[16,35-37]，传统的逐季、逐月和逐周等低频观测往往会错过许多快速变化的关键过程[23,36]，而陆基高光谱遥感观测除了本文呈现的叶绿素a、透明度、总氮、总磷和高锰酸盐指数外，还可以开展水温、浊度、有色可溶性有机物、悬浮物等连续高频监测，因此能很好捕捉其他快速变化过程，深入揭示极端气象水文等快速事件对河流、湖库和湿地水环境与水生态的影响过程与机制.

相比于航空/航天/近地无人机载的水环境遥感，陆基高光谱水环境遥感由于离水面比较近(5 m左右)，离水辐射信号强，大气的影响基本上可以忽略，无需进行大气校正，因此可以实现阴天、多云和晴朗等复杂天况下更高精度的水环境遥感监测. 同时，陆基高光谱遥感监测仪覆盖400～1000 nm可见光和近红外波段，光谱分辨率为1 nm，可以应用于包括水色参数在内的更多水质参数，实现更精准的遥感反演和监测. 当前，仪器内置了总氮、总磷、高锰酸盐指数、叶绿素、透明度和有色可溶性有机物吸收系数等11个关键水质参数算法. 未来还需继续开展大范围陆基高光谱遥感同步实验，积累覆盖范围更宽的反射率和水质数据集，优化和提升现有算法精度和适用性；构建更多水质参数深度学习算法(如颗粒和溶解性有机碳等)，支撑河流湖库碳输移和迁移转化过程监测，服务于国家“双碳”计划；对接生态环境管理部门，形成技术规程与标准，支撑断面考核和环境管理. 此外，陆基高光谱遥感不单可应用于河流、湖库、湿地、池塘等近岸水体开展监测，事实上在开敞水域的平台、船舶、桩基等固定和移动位置均可以安装陆基高光谱多参数水质遥感监测仪，实现对关键或敏感水域连续高频监测. 集成的雷达水位测定和视频监控还可以对暴雨洪水过程以及人为干扰和破坏进行监控报警.

陆基高光谱多参数水质遥感监测仪适合生态环境、水利、市政、自然资源监测等行业部门，广泛应用于全国地表水监控断面、集中式饮用水源地及其他敏感水域开展连续高频水环境监测，弥补现有的人工和自动监控断面监测在观测频次、观测参数、观测精度和观测成本上的不足. 同时，也可以用于科研院所、高校等开展监测科研研究，安装简单方便、可移动. 此外，目前的陆基高光谱多参数水质遥感监测仪还可以开发手持式移动观测以及船载和无人机载观测系列产品及成像系统，支撑全国地表水体水环境监测与应用研究，为水环境治理与修复提供监测与诊断的系统化解决方案，服务水环境管理.

致谢：太湖湖泊生态系统研究站提供相关气象水文和长期营养盐数据，在此表示感谢. 感谢杭州海康威视数字技术股份有限公司高阳辉、高晶、田野和许毅等在陆基高光谱多参数水质遥感监测仪安装和调试等方面提供的帮助.