于桥水库蓝藻水华遥感长时序监测研究
1.
2.
3.
Remote Sensing Long-term Monitoring of Cyanobacterial Blooms in Yuqiao Reservoir
1.
2.
3.
通讯作者:
收稿日期: 2019-03-27 修回日期: 2020-04-22 网络出版日期: 2020-07-03
基金资助: |
|
Received: 2019-03-27 Revised: 2020-04-22 Online: 2020-07-03
作者简介 About authors
岳昂(1986-),男,天津人,硕士,工程师,主要从事生态环境遥感监测与评估方面的研究E⁃mail:
关键词:
Keywords:
本文引用格式
岳昂, 曾庆伟, 王怀警.
Yue Ang, Zeng Qingwei, Wang Huaijing.
1 引 言
于桥水库是引滦入津和天津本地产水最重要的调蓄与汇流湖泊,保证着天津市及滨海新区千万人口的饮用水供给,是天津市赖以生存和发展的生命线[1],保障饮用水安全对天津经济建设和社会发展意义重大。由于水库上游及周边地区社会经济的发展,大量污染物直接或间接排入输水河道,导致水体富营养化趋势明显,水库水质状况不容乐观,特别是汛期存在藻类爆发的风险,致使全市饮用水存在安全隐患。由于缺乏先进的水环境监测预警体系,使监测长期滞后于环境管理发展的需求。目前蓝藻水华遥感监测方法较为成熟,且应用广泛,但多数方法是围绕大型浅水湖泊(如太湖、巢湖等),多利用MODIS、MERIS等高光谱数据和水质实测数据展开研究,建立特征波段与遥感数据间的经验、半经验模型或利用辐射传输模型,从而实现水质参数的定量化反演[2,3,4,5,6,7,8,9];构建的反演模型涵盖了常见的统计模型、机器学习和深度学习等算法。而针对小型浅水湖库(如于桥水库),MODIS、MERIS等低空间分辨率数据已不再适用,基于该类卫星数据的蓝藻水华遥感监测算法也不再适用。与此同时,中小型水库蓝藻水华监测的需求却日益迫切,缺乏水库水华现象快速遥感监测方法,为生态工程效应评估提供参考和辅助决策依据。本文利用30 m分辨率Landsat ETM/TM/OLI遥感时序数据,构建蓝藻水华的快速监测模型,从而实现中小型水库蓝藻水华的动态监测,并结合气象因子和生态工程建设的协同分析,揭示水华爆发的驱动因素,为及时发现并跟踪监测突发性水环境污染事件、水环境治理工程生态效益评估、制定和修改应对措施等方面提供有力的决策依据[10,11,12,13,14]。
2 研究区与数据
2.1 研究区概况
于桥水库位于天津北部(39°56′~40°23′N,117°26′~118°12′E),是国家重点中型水库之一,北临河北省兴隆县,西接天津蓟县,东南部与河北省遵化市、迁西县毗邻,水库流域地形差异显著,以流水剥蚀地貌为主,实际控制流域面积约2 060 km²,占整个流域面积的96%。流域内年均温为10.4~11.5 ℃、多年平均降雨750 mm,最热月均温为25.4 ℃(7月),极端最高温度为41.2 ℃,最冷月均温为-6.88 ℃(1月)、历史最低温为-28.6 ℃,无霜期约170~195 d,温度昼夜差别随地区分布略有不同,总的温度趋势呈自北向南递升。水库库区位于盆地,最大回水长东西约30 km,南北宽约8 km,最大淹没面积250 km²(正常蓄水位时淹没面积86.8 km²),具体地理位置如图1所示。该水库属富营养型水体[15],夏季处于营养化初级阶段,秋季为蓝藻水华频发期,冬季稍好。水库库区浮游植物种类繁多,藻类种类多达近百种,以蓝藻为主。
图1
2.2 数据源
卫星数据为美国陆地卫星系列(Landsat),该数据已应用于包括土地利用变化监测、城市扩展监测、水资源监测和农林监测等诸多领域,对该数据的详细介绍已多见报道,此处不再加以赘述。鉴于水库蓝藻爆发的关键时期为6~9月,本研究获取了2008~2017年7~9月共10景无云的Landsat系列数据(如表1),开展于桥水库的蓝藻水华监测和时空变化分析。
表1 数据列表
Table 1
数据源 | 时相 | 分辨率 |
---|---|---|
Landsat 5 TM | 2008年8月2日 | 30 m |
2009年7月20日 | ||
2010年8月1日 | ||
2011年9月5日 | ||
Landsat 7 ETM | 2012年7月4日 | |
Landsat 8 OLI | 2013年7月31日 | |
2014年7月11日 | ||
2015年7月14日 | ||
2016年7月16日 | ||
2017年7月10日 |
2 研究方法
2.1 数据预处理
因为目前Landsat数据均为正射校正后数据,所以对于Landsat系列遥感影像的预处理流程较为简单,主要包括影像辐射校正和大气校正,即将原始影像DN值转化为反射率值的过程。采用6S辐射传输模型对数据进行大气校正和瑞利散射校正,消除成像过程中因大气散射、水体折射、水面瑞利散射等现象带来的误差。分别对2008~2017年共10 a间的TM/Landsat系列影像进行上述预处理操作。
2.2 水华监测方法
水体中氮、磷、钾等含量增加,适宜浮游植物生长,浮游藻类颗粒短时间内大量聚集形成水华现象,最终导致水体富营养化。当蓝藻覆盖水体表面时,不仅会提高水体中叶绿素a含量,而且会导致水体透明度的降低、水体温度的升高以及水体色度的变化等,最终导致水体的光谱反射特性发生变化。对于光学遥感卫星数据而言,蓝光波段以及红光波段的反射率会明显降低,因藻类细胞内含有丰富的叶绿素,故绿光波段和近红外波段反射率升高,使得水体具有与植被类似的光谱特征,形成的反射率光谱曲线既具有水体的光谱特征,又具有植被的光谱特征[17]。
针对水华等级划分方法,有学者依据水体中蓝藻密度来区分,此方法虽然较为精准,但该法时效性差、成本高、实际应用难度较大,且难以大范围应用[18];亦有研究采用遥感影像假彩色合成法识别水华分布,具有直观明了的优点,但其影像直方图拉伸后的色调显示效果可能因人而异,给对比分析工作带来不利影响。归一化差值植被指数可以消除传感器定标、大气条件等因素引起的有关辐照度变化,可以增强对植被的响应能力,得到广泛应用。蓝藻有着与植被相似的光谱反射特征,绿色植物能吸收大部分的可见光,在近红外波段却具有强烈的反射,因此可以利用归一化差值植被指数(Normalized Difference Vegetation Index, NDVI)监测蓝藻和水华的分布和程度,其计算公式如
其中:
表2 实测样地水质参数
Table 2
样点号 | 总磷(mg/L) | 总氮(mg/L) | 叶绿素a(mg/L) | NDVI |
---|---|---|---|---|
1 | 0.023 | 0.243 | 1.93 | 0.351 |
2 | 0.010 | 0.151 | 1.14 | 0.110 |
3 | 0.004 | 0.050 | 0.22 | -0.12 |
4 | 0.052 | 0.500 | 9.13 | 0.212 |
5 | 0.234 | 2.780 | 12.0 | 0.402 |
6 | 0.010 | 0.049 | 0.19 | -0.111 |
7 | 0.232 | 4.740 | 12.0 | 0.419 |
8 | 0.223 | 4.550 | 10.3 | 0.399 |
9 | 0.224 | 4.620 | 10.0 | 0.398 |
10 | 0.035 | 0.331 | 9.23 | 0.302 |
表3 水质参数与NDVI相关性分析
Table 3
参数 | 相关系数(R) |
---|---|
总磷 | 0.748 |
总氮 | 0.702 |
叶绿素a | 0.817 |
图2
图2
水质参数与NDVI变化关系
Fig.2
Relation between water guality parameters and NDVI change
NDVI的取值范围在[-1,1]之间,对于植被而言,由于近红外反射率值远大于红光波段,故植被覆盖区域NDVI值一般大于0;对于纯净水体来说,由于水体的红光波段反射率大于近红外反射率,故 NDVI一般小于0。当水体覆盖区域的NDVI大于0时,则表明水体被植被覆盖或者水体中叶绿素含量较高。
4 结果和分析
图3
图3
2008~2017年于桥水库库区水华等级分布图
Fig.3
Grade distribution of water blooms in Yuqiao reservoir area from 2008 to 2017
表4 2008~2017年于桥水库水华等级及面积统计表
Table 4
年份 | 2008 | 2009 | 2010 | 2011 | 2012 | 2013 | 2014 | 2015 | 2016 | 2017 |
---|---|---|---|---|---|---|---|---|---|---|
轻度水华区/km² | 39.11 | 52.66 | 60.99 | 64.31 | 67.09 | 37.63 | 5.36 | 63.93 | 44.22 | 4.81 |
中度水华区/km² | 3.72 | 1.71 | 5.61 | 2.53 | 0.88 | 2.01 | 1.67 | 2.67 | 2.97 | 1.59 |
重度水华区/km² | 3.53 | 0.34 | 2.29 | 1.18 | 0.07 | 1.39 | 2.92 | 2.32 | 5.40 | 4.77 |
为了直观表达2008~2017年10 a间水华分布及不同等级水华面积的变化情况,以年份为横轴,以各水华等级面积所占百分比为纵轴,绘制如图4所示的堆叠柱状图,从图中可以发现:①2008~2012年的5 a间,于桥水库库区发生水华的总面积总体呈现上升趋势,其中总面积和中度级别面积的峰值均出现在2010年表明2008~2010年间水库库区水生环境在不断发生恶化,2010~2012年间总面积基本持平,但中度以上水华区面积逐年减少,以轻度水华为主,到2012年中度、重度级别水华区为自2008年以来最小;②2013~2014年,轻度水华面积显著减少,但重度水华区逐渐增加,这个主要是2014年天津市水务局斥资105亿元大力推动清水河道行动、引滦水源保护工程、南水北调市内配套工程、防洪排水工程及农村水利基础设施工程等五大重点工程建设。由此表明生态治理工程对水华治理效果显著,2014年水华总面积达到最小;③2015年引滦水量少于往年[25],且上游来水水质较差,导致于桥水库水深变浅,水库中多年积累的污染物质,受风浪和水温影响,底泥被搅动,悬浮物被快速释放,含氮、磷的营养物质进入水体,为藻类爆发性生长提供营养物质,造成于桥水库水体富营养化再次暴发;④2016~2017年,水库库区水华面积整体呈逐年减少趋势,至2017年降低到最小值。进一步分析可以发现,水华面积减少主要为轻度水华区域,中和重度水华在轻度水华显著减少时先增加后略减少态势,表明轻度水华在较短时间内可以恢复,而中、重度水华则难以在短期内得到有效恢复和治理,水华等级与治理恢复难度成正相关关系;此外,重度水华集中分布在距河岸较近区域,表明人为活动是内陆水华暴发的重要诱因之一。
图4
图5
图5
2008~2017年天津市年均降雨量及年均温变化图
Fig.5
Change map of annual mean rainfall and annual mean temperature in Tianjin from 2008 to 2017
5 结 语
通过分析NDVI值与实测水质参数之间的相关关系,选择水华等级的NDVI分割阈值,开展了于桥水库库区蓝藻水华现象年际变化遥感快速监测,掌握了库区10 a间水华发生面积和程度的动态变化情况,并结合当地气象资料协同分析水华爆发的影响因素。结果:①水体的叶绿素a含量与NDVI相关性较强,表明利用卫星遥感技术监测中小型水面水体蓝藻水华具有可行性;②2008~2017年间于桥水库水华爆发面积呈先上升后下降的趋势;③水华爆发受自然(气温和降水等)和社会活动等诸多因子共同驱动和影响;④外界人为扰动造成的影响尤为显著,生态修复工程实施后水华爆发趋于减缓,水质状况得以有效改善。
需要说明的是分析用时间序列卫星遥感数据为3个传感器获取,时相差异最大达57 d。实际上年际间夏秋季气温、降水的平均值和极端值等的差异都会对水华形成有重要影响,而且库区周边人类活动(如排放等)的干扰也有待深入探究。随着国产高分卫星技术发展,后续研究可利用多星组网联合观测建立高空间、高频度卫星监测体系,为业务管理部门提供准实时监测预警服务,全力保障饮用水水源地生态环境的安全。
参考文献
Eutrophication Trend and Causes of Yuqiao Reservoir
[J]
于桥水库富营养化趋势及成因
[J].
Research on the Spatial Variation of Water Quality Parameters in East Dongting Lake based on Hyperion
[J].
基于Hyperion的东洞庭湖水质参数空间分异规律
[J].
Muti-model Collaborative Retrieval of Chlorophyll a in Taihu Lake based on Data Assimilation
[J].李渊
基于数据同化的太湖叶绿素多模型协同反演
[J].
An Operational Algorithm to Estimate Chlorophyll-a Concentrations in Lake Chaohu from MODIS Imagery
[J].
一种基于MODIS影像可业务化运行的巢湖水体叶绿素a估算算法
[J].
A Novel MODIS Algorithm to Estimate Chlorophyll a Concentration in Eutrophic Turbid Lakes
[J].
Deriving Concentration of TP, TN based on AISA Hyper Spectral Reflectivity
[J].崔爱红
基于AISA高光谱数据的总磷总氮反演
[J].
Water Quality Retrievals from MERIS Satellite Data in Yangtze Estuary
[J].
长江口水质MERIS卫星数据遥感反演研究
[J].
A MODIS-based Novel Method to Distinguish Surface Cyanobacterial Scums and Aquatic Macrophytes in Lake Taihu
[J].
Ensemble Modeling Methods for Remote Sensing Retrieval of Water Quality Parameters in Inland Water
[J].
内陆水体水质参数遥感反演集合建模方法
[J].
Environment Satellite HJ-1 CCD Data Dynamically Monitoring Cyanobacteria Bloom in Lake Dianchi, China
[C]∥
Investigation and Analysis of Water Pollution Factors in Yuqiao Reservoir
[J].
于桥水库水污染因素调查分析
[J].
Development of Water Quality Parameter Retrieval Algorithms for Estimating Total Suspended Solids and Chlorophyll-a Concentration Using Landsat 8 Imagery at Poteran Island Water
[C]∥
Preliminary Study on Water Stratification Characteristics of Yuqiao Reservoir based on Water Quality Field Monitoring Parameters
[J].
基于水质现场监测参数的于桥水库水体分层特征初步研究
[J].
. Remote Monitoring of Cyanobacterial Blooms Using Multi-source Satellite Data: A Case of Yuqiao Reservoir, Tianjin
[J].
基于多源卫星数据的小型水体蓝藻水华联合监测--以天津于桥水库为例
[J].
Investigation of Summer Plankton Community Structure in Yuqiao Reservoir,Tianjin
[J].武丹,王海英,张震
天津于桥水库夏季浮游生物调查及群落结构变化
[J].
Progress and Prospects on Cyanobacteria Bloom-forming Mechanism in Lakes
[J].
湖泊蓝藻水华发生机理研究进展
[J].
Evaluation Method and Application on Cyanobacteria Bloom Degree Classification With Algal Density
[J].
关于用藻密度对蓝藻水华程度进行分级评价的方法和运用
[J].
Extraction Methods of Cyanobacteria Bloom in Lake Taihu based on RS Data
[J].
基于遥感数据的太湖蓝藻水华信息识别方法
[J].
CyanObacteria Bloom Monitoring with Remote Sensing in Lake Taihu
[J].
太湖蓝藻水华遥感监测方法
[J].
Spectral Features, Remote Sensing Idengtification and Breaking-out Meteorological Conditions of Algal Bloom in Lake Dianchi
[J].
滇池蓝藻水华光谱特征、遥感识别及暴发气象条件
[J].
Recognition of Cyanobacteria Bloom based on Spectral Analysis of Remote Sensing Imagery
[J].
基于遥感影像光谱分析的蓝藻水华识别方法
[J].
Monitoring and Annual Statistical Analysis of Algal Blooms in Chaohu based on Remote Sensing
[J].
巢湖水华遥感监测与年度统计分析研究
[J].
Research on Eutrophic Evaluation and Control Measures for Yuqiao Reservoir
[J].
于桥水库水质富营养化评价及防治对策研究
[J].
. Impact of A "Source-Sink" Landscape Pattern in An Urbanized Watershed on Nitrogen and Phosphorus Spatial Variations in Rivers: A Case Study of Yuqiao Reservoir Watershed, Tianjin, China
[J].
城镇化流域“源-汇”景观格局对河流氮磷空间分异的影响—以天津于桥水库流域为例
[J].
Research on the Meteorological and Hydrological Factors of Cyanobacteria Bloom in Lake Dianshan
[J].
淀山湖蓝藻水华暴发的气象水文因素探讨
[J].
Remote Sensing and Analysis on Meteorological Factors of Blue Algal Bloom in Lake Tai
[J].
卫星遥感太湖蓝藻水华分布及其气象影响要素分析
[J].
Effect of Meteorological Conditions on Blue Algal Bloom Distribution in Dianchi Lake
[J].
气象条件对滇池水华分布的影响
[J].
Analysis on the Distributions and Meteorological Conditions of Cyan Obacteria Bloom in Chaohu Lake
[J].
巢湖蓝藻水华分布特征和气象条件分析
[J].
/
〈 |
|
〉 |
