基于RS 的于桥水库藻华时空分布特征及影响因子分析岳㊀昂1㊀张㊀赞1㊀魏㊀巍1∗㊀韩静敏1㊀刘红磊2(1.天津市生态环境监测中心,天津300191;2.天津市环境保护科学研究院,天津300191)摘要:于桥水库蓝藻水华的出现造成水体富营养化,影响饮用水安全㊂遥感技术为监测和预防藻华提供了高效快捷的手段,以2008 2017年的Landsat 卫星影像为数据源,通过NDVI 指数反演,分析了于桥水库历年藻华分布特征和变化趋势㊂从空间分布特征来看,于桥水库藻华主要分布在东岸和北岸,严重时会发展到水库中心;从历年藻华面积变化趋势看,2008 2013年,于桥水库藻华爆发面积较小,2013 2015年藻华现象逐年严重,2016 2017年水库藻华面积显著减少并趋于平稳㊂根据遥感监测结果,结合实地监测数据,分析得出:磷是于桥水库藻华爆发的限制性因子,气温的阶段性上升可能导致藻华大面积爆发㊂基于此研究结果,提出以磷输入控制为核心目标的于桥水库藻华防治措施,以减少水体富营养化和藻华爆发㊂关键词:于桥水库;遥感;藻华;影响因子;磷SPATIOTEMPORAL DISTRIBUTION CHARACTERISTICS OF ALGAL BLOOM IN YUQIAORESERVOIR BASED ON RS AND ITS INFLUENCE FACTORS ANALYSISYue Ang 1㊀Zhang Zan 1㊀Wei Wei 1∗㊀Han Jingmin 1㊀Liu Honglei 2(1.Tianjin Eco-Environmental Monitoring Center,Tianjin 300191,China;2.Tianjin Academy of Environmental Sciences,Tianjin 300191,China)Abstract :The appearance of cyanobacteria blooms in Yuqiao Reservoir caused eutrophication of the water body and affectedthe safety of drinking water.Remote sensing technology provides an efficient and fast way to monitor and prevent algal blooms.Landsat satellite remote sensing images from 2008to 2017were used to analyze the distribution characteristics and changing trends of algal blooms in Yuqiao Reservoir through NDVI index inversion.From the perspective of spatial distributioncharacteristics,the algal blooms are mainly distributed on the east and north shores while in severe cases will develop to the center of the reservoir.From the perspective of the algal bloom area changes over the years,the area is relatively small,but more serious year by year from 2013to 2015,and decreased significantly and stabilized from 2016to 2017.Based on theresults of remote sensing monitoring and field monitoring data,the research showed that phosphorus is the limiting factor for algal blooms in the Yuqiao Reservoir,and the stepwise increase in temperature may cause large-scale blooms of algal blooms.Based on this result,algae bloom prevention measures in Yuqiao Reservoir were proposed with the core objective of phosphorus input control to reduce water eutrophication and algal blooms.Keywords :Yuqiao Reservoir;Remote sensing;Algal bloom;Impact factor;Phosphorus㊀㊀㊀㊀㊀㊀㊀㊀收稿日期:2020-03-21基金项目:天津市科技计划项目(16YFXTSF00380)㊂第一作者:岳昂,工程师,研究方向为水环境监测㊁生态环境遥感监测研究㊂∗通信作者:魏巍(1970-),男,高级工程师,主要从事生态环境监测相关工作㊂1003813938@0㊀引㊀言水是人类生存和进步的生命之源,是社会经济发展不可替代的基础自然资源㊂近年来,随着社会经济的迅猛发展,城镇化程度加速扩张,工农业生产和生活污水的排放对内陆湖泊河流造成了不同程度的污染,导致水体富营养化问题出现,水华蓝藻现象爆发,成为水生态系统最重要的问题之一,给内陆湖泊生态环境治理带来了严峻的考验㊂传统的藻华测量方法是通过现场采样后在实验室进行测量,这种方法不仅费时费力,而且无法监测大范围水体的藻华分布情况㊂遥感技术的出现和快速发展,为监测和预防水华蓝藻提供了高效快捷的手段㊂遥感监测具有监测范围广㊁速度快㊁成本低㊁便于进行长期动态监测的优势,利用遥感技术进行蓝藻水华监测能够较好的反映藻华的时空分布和演化规律,可以帮助环境监测相关部门快速掌握藻华的变化趋势,为藻华预警提供科学依据[1]㊂在内陆水体水质遥感监测中,对蓝藻水华的监测主要是针对藻华覆盖的分布范围和动态趋势的监测㊂利用遥感技术监测蓝藻水华的方法包括单波段阈值法㊁光谱指数法㊁监督分类法和水质参数反演法等㊂单波段阈值法是利用藻华在近红外波段的高反射特性,通过选取合适的阈值区分藻华和水体,从而获取藻华分布范围㊂段洪涛等[2]利用MODIS数据,选取0.1为阈值,认为band>0.1时为藻华,提取了海岸带水域的蓝藻水华㊂光谱指数法是利用不同的波段组合方式,通过计算波段比值的方法实现藻华的有效监测㊂李旭文[3]等利用Landsat TM数据通过构建差分植被指数DVI来计算藻类生物量,并对藻类生物量进行分级,对太湖梅梁湖区藻类分布特征进行了监测和分析㊂Hu等[4]通过建立浮游藻类光谱指数提取了青岛海岸的藻类,并将改方法应用于太湖,对太湖2000 2009年藻华分布情况进行了监测㊂监督分类法是利用蓝藻和水体的光谱特征和纹理信息的差异,通过训练样本选择利用遥感图像监督分类的方法提取蓝藻水华㊂林怡[5]利用Landsat ETM+数据构建的归一化蓝藻指数,采用支持向量机(Support Vector Machine,SVM)方法进行监督分类提取淀山湖蓝藻水华,有效地识别了各种密度的蓝藻空间分布信息㊂水质参数反演法是以蓝藻素或叶绿素a为目标,通过遥感定量反演,分析蓝藻的分布状态㊂例如,梅长青[6]㊁怀红燕[7]等利用MODIS数据对我国巢湖㊁淀山湖进行叶绿素反演,通过叶绿素浓度监测藻类水华的暴发㊂于桥水库是天津市唯一的城市集中式应用水水源地,是天津市市区及滨海新区千万人口的生命线,对天津市经济社会发展具有极其重要的作用㊂但是近年来,由于农业面源污染和上游城市污水的排放[8-9],污染加剧,给于桥水库水质状况带来了严重影响,造成了水体富营养化,从1997年就已经出现了蓝藻水华爆发[10],给天津市饮用水安全带来了极大隐患,因此于桥水库藻华监测及爆发因素分析尤为重要㊂本文利用Landsat8和landsat7卫星遥感影像,通过光谱指数计算,提取于桥水库近10年来的藻华空间分布情况,分析其时空分布特征,并结合实地水质监测数据和区域温度数据,分析于桥水库总磷㊁叶绿素a与藻华面积的关系,探讨藻华爆发的影响因子,为于桥水库藻华治理和水质监管提供支撑㊂1㊀藻华遥感解译1.1㊀研究区概况于桥水库位于天津市蓟县城东4km处的州河上游,是一座山谷型水库㊂于桥水库始建于1959年, 1983年引滦入津工程建成后,于桥水库正式纳入引滦入津工程管理,成为天津唯一的水源地,其主要功能以防洪㊁城市供水为主,兼顾灌溉㊁发电等㊂于桥水库上游主要入库河流为淋河㊁沙河和黎河,其中沙河和黎河汇聚而成果河,黎河为引滦输水通道,上游连接输水隧洞,由于输水影响,果河入库的氮㊁磷质量浓度持续超标[11],导致库区水质无法达到规划要求的Ⅲ类标准,直接影响下游的蓟县㊁宝坻㊁宁河㊁玉田㊁汉沽等各县(区)的低洼地区近百万人口和300余万亩耕地㊂1.2㊀数据介绍遥感数据采用了Landsat-8OLI和Landsat-7ETM+卫星遥感图像,数据源来自地理空间数据云(http:// /)㊂Landsat-7是美国的陆地卫星计划(Landsat)中的第七颗,卫星携带增强型专题制图仪(Enhanced Thematic Mapper,ETM+)传感器,与Landsat-5相比增加了分辨率为15m的全色波段(PAN波段)㊂Landsat-8是太阳同步轨道卫星,携带陆地成像仪(Operational Land Imager,OLI)和热红外传感器(Thermal Infrared Sensor,TIRS),OLI陆地成像仪包括9个波段,空间分辨率为30m,其中包括一个15m的全色波段㊂与Landsat5/Landsat7相比, Landsat8涵盖的波段更多㊁波段划分更精细,数据量也提高到以往的3倍[12,13]㊂本文选取2008 2017年8 9月的Landsat-8OLI和Landsat-7ETM+影像为主要数据源进行于桥水库藻华解译分析㊂气象数据来自中国气象科学共享服务系统(http://图1㊀于桥水库遥感影像示意/),获取了流域内气象站(站点编号54439,坐标为117ʎ53ᶄE㊁40ʎ12ᶄN)的日平均和最大气温,分析气温与于桥水库叶绿素a质量浓度和藻华爆发面积的关系㊂1.3㊀藻华提取方法研究1.3.1㊀Landsat数据预处理利用ENVI软件对Landsat数据进行区域裁剪㊁辐射定标㊁大气校正等预处理[14],消除传感器误差㊁大气误差等外界干扰因素的影响,并采用单波段阈值法制作水体区域掩膜[15],获得准确的水体区域,在此基础上再进行蓝藻水华的提取㊂1.3.2㊀蓝藻水华提取方法蓝藻水华具有与植被相似的光谱特征[16],在550nm(绿波段)附近有一个反射峰,在675~710nm 波段具有类似于植被的陡坡效应,而在710~900nm (近红外波段)由于藻类细胞结构特点强烈反射太阳光谱能量而具有一个极高的反射率平台[17],这是蓝藻水华区别于正常水体的最明显的特征,因此利用植被指数可以有效的进行蓝藻水华的提取㊂归一化植被指数(Normalized Difference Vegetation Index,NDVI)是最常用且提取效果较好的植被指数之一,由于蓝藻光谱特征具有与植被相似的陡坡效应,所以利用NDVI可以有效的识别蓝藻水华㊂NDVI的计算公式如下:NDVI=NIR-RNIR+R(1)式中:NIR为近红外波段反射值;R为红波段反射值㊂1.4㊀藻华空间分布及演变趋势分析1.4.1㊀藻华面积解译获得于桥水库的水体区域后,利用ENVI的波段运算工具进行NDVI反演,得到于桥水库的NDVI图㊂在NDVI图像中,NDVI>0的区域为藻华区域,且NDVI值越大藻华浓度越高,NDVI值越小藻华浓度越低;NDVI<0的区域为正常水体㊂对得到的NDVI 图像进行阈值分割,得到于桥水库藻华分布图㊂为了消除其他易混淆物体的干扰,得到精准的藻华分布,通过对象矩形度进一步判断藻华的面积分布㊂对象矩形度即为表明形状能被矩形描述的程度,由多边形面积比多边形外接矩形的面积表示㊂矩形的矩形度为1,非矩形的矩形度小于1㊂对象矩形度是识别藻华的核心因子,根据对目标对象的判断将对象矩形度设定为0.7~0.8,可以有效识别出农田大体特征,消除对藻华提取的影响,但应注意前一步分割参数应设定相对细微㊂此外,由于水流㊁风向等外部环境都会影响到藻华的聚集,因此在卫星影像上含有藻华的水体与正常水体表现了不同的表面特征,即藻华在水体表面一般表现为条带状或者絮状形态,容易与周围湖水相区分,有利于目标物的目视解译[18],利用这一特征对获得藻华分布图进行后处理,最终得到准确的于桥水库藻华面积分布图㊂1.4.2㊀于桥水库近10年藻华空间分布特征及演变趋势分析利用2008 2017年8 9月于桥水库Landsat数据,采用NDVI指数反演结合人工目视解译的遥感方法,提取于桥水库2008 2017年的藻华空间分布图,分析10年间于桥水库藻华变化趋势,结果见图2㊂从空间分布特征来看,于桥水库藻华主要分布在东岸和北岸,严重时会发展到水库中心,西南沿岸基本没有藻华出现,水质状况较好,造成这一分布特征的原因可能与盛行风向㊁入库支流营养盐数据[19]以图2㊀于桥水库2008 2017年夏秋季藻华时空分布和藻华面积变化趋势及水库周边地貌状况有关[20]㊂从历年藻华面积变化趋势看,2008 2013年,于桥水库藻华爆发面积较小,水质状况较好,其中2009年藻华面积最小为1.87km2,2010 2013年藻华面积相对平稳,基本维持在3km2左右;2013 2015年藻华现象逐年严重,到2015年达到一个峰值,藻华面积达到5.65km2; 2016 2017年水库藻华面积显著减少并趋于平稳,藻华面积分别为1.5km2和1.95km2㊂2㊀藻华影响因子分析2.1㊀于桥水库近10年水质情况为研究水质环境与藻华面积变化的关系,本文分析了2008 2017年于桥坝下㊁水库中心㊁东马坊㊁九百户㊁三岔口不同监测点位的总磷(Total Phosphorus, TP)和叶绿素a(Chl-a)的质量浓度变化情况,如图3所示㊂整体来看,在同一监测时间,不同监测点位的TP和Chl-a质量浓度没有明显差异,且五个监测点位上的TP质量浓度与Chl-a质量浓度的历年变化趋势基本一致㊂2008 2012年TP和Chl-a质量浓度均较为平稳,没有明显变化,从2013年开始,TP和Chl-a质量浓度大幅上升,并到2016年达到最大,2017年开始有回落趋势㊂2008年于桥水库TP和Chl-a平均质量浓度分别是0.027mg/L和6.7mg/m3,到2016年分别上升至0.110mg/L和24.3mg/m3,2017年又回落至0.035mg/L和13.3mg/m3㊂ʏ TP; ʻ Chl-a㊂图3㊀于桥水库5个监测站点的总磷和叶绿素a质量浓度变化(2008 2017年)2.2㊀于桥水库磷质量浓度、叶绿素a与藻华面积关系分析分析于桥水库2008 2017年实测的TP和Chl-a 质量浓度可以发现,二者呈现明显的线性关系(R2= 0.944),如图4a所示㊂由此可见,TP质量浓度的增加会促进藻类生长,进而造成水体Chl-a质量浓度的上升,这一结果与Smith等[21]的研究结果是一致的,也进一步证明了磷是于桥水库藻华爆发的限制性因子㊂如图4b是藻华面积与叶绿素a的关系图,从图中可以看出二者没有明显关联㊂结合磷质量浓度与藻华面积的历年变化情况可以看出,磷质量浓度与藻华面积变化也不完全一致㊂从藻华面积变化趋势图中可知,于桥水库藻华面积在2015年达到最大值5.65km2,而磷质量浓度和叶绿素a浓度的最大值均是在2016年,其浓度分别为0.109mg /L 和52.44mg /m 3,而此时藻华面积相较2015年已大幅减少,甚至为历年最少㊂由此可知,虽然TP 浓度增加会促进藻华爆发,然而藻华爆发面积与磷质量浓度没有必然联系㊂这是因为磷并非是引起藻华爆发的唯一因素,藻华发生会受到营养盐浓度㊁温度㊁湖泊地形㊁湖流等多种环境因素的影响㊂于桥水库的输入方式主要是河流输入[22],河流入库口由于流速降低,吸附于颗粒上的磷大量沉积,并水土界面发生交换作用重新进入上覆水体,导致水中磷质量浓度升高[23],造成有支流输入的河口和近岸区域藻华滋生和大面积聚集(见图2)㊂因此,在于桥水库藻华防治工作中,应重点关注水流较为缓慢的岸边区域和有支流输入的河口缓流区域㊂图4㊀叶绿素质量浓度对总磷质量浓度的响应关系及藻华面积与叶绿素质量浓度2.3㊀区域气温与于桥水库藻华爆发关系分析水温对藻类生长具有重要影响,其中25ħ左右是较适宜藻类生长的温度[24],因此每年的8 9月是于桥水库藻华爆发最为严重的时间㊂利用国家气象站(遵化站)2008 2017年8 9月的日平均和最大气温数据,分析不同年份气温变化与于桥水库藻华面积变化的关系㊂结果显示,2015年8 9月日平均气温分布与其他年份有明显不同㊂2015年8 9月,日平均气温在22~28ħ的累计频率为65%左右,明显高于其他年份,说明该部分时段温度高于同期温度㊂进一步分析发现在藻华面积分布较高的年份均有类似的温度分布特征,充分证明了在营养盐充足的条件下,气温的升高会促进藻类大量繁殖,导致藻华大面积爆发㊂ʻ 2008; ʏ 2009; ▶ 2010; һ 2011; ▽ 2012;⎔ 2013; ʻˑ 2014; Җ 2015; ◁ 2016; Ѳˑ 2017㊂图5㊀于桥水库2008 2017年8 9月份日平均气温累积频度分布3㊀结㊀论从上述分析可以看出,于桥水库藻华主要出现库区东岸和北岸的河流入口处,西南及水库中心水体状况较好,藻华出现较少㊂引起藻华现象的原因包括营养盐输入㊁磷质量浓度的增加㊁温度等因子,其中营养盐是藻华发生的物质基础,总磷与水中叶绿素a 浓度有着直接影响,是于桥水库藻华爆发的限制性因子;温度也与藻华爆发密切相关,对藻华爆发面积有着重要影响㊂基于此,在于桥水库蓝藻水华防治方面,建议以磷输入控制为核心目标,具体对策如下:1)控制面源磷输入和废污水排放,实施基于入库磷总量达标的流域水质目标管理方案[25];2)建造人工湿地,开展入库河流岸边带修复和重建,减低入河磷污染和净化水质;3)在东岸和北岸污染较重区域开展生态清淤,减少沉积物内源磷污染㊂参考文献[1]㊀赵丹.典型内陆湖库蓝藻水华遥感监测[D].西安科技大学,2018.[2]㊀张渊智,段洪涛.芬兰海岸带水域蓝藻水华遥感监测[J].湖泊科学,2008,20(2):167-172.[3]㊀李旭文,季耿善,杨静.太湖藻类的卫星遥感监测[J].湖泊科学,1995(7):65-68.[4]㊀HU C,LEE Z,MA R,et al.Moderate resolution imagingspectroradiometer (MODIS)observations of cyanobacteria bloomsin Taihu Lake,China [J ].Journal of Geophysical Research:Oceans,2010,115(C4).(下转第83页)。