环境科学  2020, Vol. 41 Issue (11): 4914-4923   PDF    
2009~2018年太湖湖泛强度变化及其影响因素
辛华荣1, 朱广伟2, 王雪松3, 李骏3, 陆隽1     
1. 江苏省水文水资源勘测局, 南京 210029;
2. 中国科学院南京地理与湖泊研究所湖泊与环境国家重点实验室, 太湖湖泊生态系统研究站, 南京 210008;
3. 江苏省水文水资源勘测局常州分局, 常州 213001
摘要: 基于2009~2018年的每年4~10月对太湖湖泛易发区的逐日巡查成果,以及相关太湖气象、水文、水质和水华监测资料的综合分析,揭示了太湖湖泛发生的环境条件及年际差异性特征,分析了影响太湖湖泛强度年际变化的因素及防控途径.结果表明,2009~2018年10 a中太湖共发生大小湖泛事件75次,平均湖泛面积1.35 km2,最大面积9.20 km2,平均持续时间3 d,最长持续16 d;湖泛发生引起了水体有机质、氮和磷等水质指标明显增高;湖泛的发生均在水温20℃以上,发生首日均处在5~9月;年际湖泛强度变化较大,2017年湖泛强度最大,2018年次之,而2014年湖泛强度最小;年际湖泛强度与蓝藻水华强度及5~9月的湖水平均温度显著正相关,而与水体营养盐等指标关系不密切,表明年际气候条件的波动对湖泛情势影响很大;湖泛的形成与河口区蓝藻水华的堆积关系密切:除了5次草源性湖泛外,70次藻源性湖泛均发生在入湖河口附近水域,河口区的底泥污染状况及底泥悬浮状况可能对湖泛的发生具有促进作用.根据形成湖泛的腐烂物质、来源、发生地点和持续时间的不同,可将湖泛分为藻源性湖泛和草源性湖泛这2种类型,藻源性湖泛又可分为港源型、迁移型和原发型这3种.结果表明,控制蓝藻水华强度是降低湖泛风险的根本途径,而避免河口区蓝藻水华过度堆积及厌氧分解的措施,如对蓝藻水华的及时打捞、工程措施形成的局部流场改变、蓝藻水华的离岸打捞工程和滨岸水华堆积区的应急曝气工程,以及河口区的污染底泥疏浚等工程措施,是降低湖泛发生风险的工程措施选项.
关键词: 太湖      湖泛      蓝藻水华      气象条件      环境效应      防控途径     
Variation and Driving Factors of Black Water Event Intensity in Lake Taihu During 2009 to 2018
XIN Hua-rong1 , ZHU Guang-wei2 , WANG Xue-song3 , LI Jun3 , LU Jun1     
1. Jiangsu Hydrology & Water Resources Investigation Bureau, Nanjing 210029, China;
2. Taihu Laboratory for Lake Ecosystem Research, State Key Laboratory of Lake Science and Environment, Nanjing Institute of Geography and Limnology, Chinese Academy of Sciences, Nanjing 210008, China;
3. Jiangsu Hydrology & Water Resources Investigation Bureau Branch Office in Changzhou, Changzhou 213001, China
Abstract: Using daily survey and monitoring "black water event" (BWE) results in the sensitive area of Lake Taihu from April to October, 2009-2018, as well as the BWE analyzed data for relative meteorological, hydrological, chemical, and algal bloom conditions, the characteristics and yearly differences of BWEs were summarized. A BWE control strategy was suggested. There were 75 BWE occurrences detected in the past 10 years. The average area of a BWE was 1.35 km2, with a maximum area of 9.20 km2. The BWEs lasted for an average of three days, while the longest lasted 16 days. The BWEs significantly increased organic matter, total nitrogen, total phosphorus, ammonia, and sulfate, among others. All the BWEs occurred at water temperatures over 20℃. All the BWE occurrences started between May and September. The yearly BWE intensity (BWEI) varied significantly among years, with the strongest intensity in 2017 and the second strongest in 2018. The BWEI was significantly positively related to yearly algal bloom intensity (ABI) and average daily water temperature from May to September, while there was no significant relationship with major nutrient indicators. This suggests that climatic variation among years will significantly influence the risk of a BWE in Lake Taihu. The occurrence of a BWE was significantly influenced by a polluted river mouth. Almost all the BWEs occurred near river mouths, except for five macrophyte-related BWEs. This suggests that sediment pollution and its resuspension could be strongly related to the occurrence of a BWE. When considering degradation factors, a BWE could be classified as an algal-related BWE and a macrophyte-related BWE. The algal-related BWE could be further classified into three types:river-related BWE, bloom transport BWE, and local origin BWE. This research suggests that algal bloom control will be the fundamental countermeasure to decrease the risk of a BWE. Bloom abatement treatments, including mechanical bloom cleaning, water current adjustment engineering, coast-away bloom cleaning engineering, mechanical aeration treatment, and sediment dredging near river mouths would be effective methods to abate the risk of a BWE.
Key words: Lake Taihu      black water event(BWF)      cyanobacterial bloom      meteorological condition      environmental effect      control countermeasure     

湖泛, 亦称黑水团或污水团, 是指富营养化湖泊水体中, 在藻类或水草大量暴发、积聚和死亡后, 在适宜的气象和水文条件下, 伴随底泥中的有机物在缺氧和厌氧条件下产生的生化反应, 释放硫化物、甲烷和二甲基三硫等硫醚类物质, 形成褐黑色伴有恶臭的“黑水团”, 从而导致水体水质迅速恶化、生态系统受到严重破坏的现象[1].湖泛现象对湖泊的供水功能、旅游功能、渔业功能等危害极大, 是引发2007年5月太湖无锡供水危机的元凶[2], 也成为太湖等蓝藻水华肆虐的湖库中水质安全保障与治理的重点对象[3]. 2008年太湖湖西区又发生持续近16 d、面积近20 km2的湖泛, 引进了社会与学术界的高度重视[3].自2009年起, 江苏省水文水资源勘测局启动了每年4~10月的太湖湖泛易发区逐日巡查监测和日报工作[4].

认知湖泛发生的关键条件, 摸清湖泛发展变化规律, 是湖泛灾害防控的基础.陆桂华等[1]通过2007~2008年2次湖泛发生条件的分析, 揭示了水体污染中、蓝藻水华物质堆积和湖底流泥存在等是湖泛发生的物质条件, 而水位偏低、气温偏高和稳定的向岸风是湖泛形成的促发条件.此后, 陆桂华等[5]又在2009年“湖泛”巡查成果的基础上, 进一步分析了藻情、底泥、水质及水位与湖泛发生的关系, 提出了水位偏高有利于湖泛防控的观点.范成新等[6]从湖泛发生的化学机制分析, 梳理了湖泛生消过程原理, 并试图定量化回答湖泛形成的物质基础和微生物条件.王成林等[7]则重点从风场和温度陡变等角度, 提出湖泛事件的暴发可能与风场及温度的突变有关.上述工作基本都是基于2007~2009年湖泛监测结果进行的分析, 当时获得的具体湖泛案例比较少, 对湖泛形成条件的认知有一定的局限性.

近年来太湖湖泛状况又出现了新变化.刘俊杰等[4]基于9 a的湖泛巡查成果分析表明, 2017年太湖湖泛的暴发强度和次数均超过了有湖泛记录以来的所有年份, 该研究还发现, 水位波动强度等水文条件与湖泛强度之间存在显著正相关.因此, 影响湖泛环境因素仍需进一步研究.在进行湖泛日巡查10 a之后, 有必要针对所获得的大量湖泛案例, 进一步分析湖泛发生的环境特征, 探讨影响湖泛强度的机制与关键条件.

本文基于江苏省水文水资源勘测局太湖10 a(2009~2018年)的湖泛巡查资料, 结合中国科学院太湖湖泊生态系统研究站太湖蓝藻水华遥感监测数据, 以及对太湖逐日水温、水位、降雨量和逐月水质观测资料, 分析了太湖湖泛发生频次及其区域特征, 探讨了影响年尺度上湖泛强度的环境条件, 以期为湖泛灾害的精准防控提供科学支撑.

1 材料与方法 1.1 太湖湖泛巡测方案及湖泛强度评估

本文采用2009~2018年江苏省水文水资源勘测局湖泛巡测数据.基于2007年5月无锡饮用水危机的教训, 以及对历史上湖泛发生规律的初步分析[2, 5], 江苏省太湖水污染防治办公室每年发布“江苏省太湖应急防控工作实施方案”, 将4月10日至10月20日定为“太湖安全度夏期”, 组织环太湖各城市的环境、水利、气象和城建等多部门, 联合应对太湖蓝藻水华带来的水环境风险和饮用水安全保障.按照该方案的布置, 江苏水文部门自2009年4月开始对太湖开展湖泛巡查, 至10月20日结束, 此后每年在“太湖安全度夏期”实施湖泛巡查, 巡查区域和固定监测点位见图 1.

图 1 太湖湖泛巡查线路及固定水质测定点位示意 Fig. 1 Black patch event survey area and the water monitor sites in Lake Taihu

根据巡查规范, 巡查水域包括了太湖的湖泛易发区及重要饮用水源地, 具体巡查区域分为3个:①太湖西岸带八房港至竺山湖区域, 由常州市水文分局执行巡查任务; ②马山月亮湾、梅梁湾及贡湖湾区域, 由无锡市水文分局执行巡查任务; ③东部湖区, 包括苏州市在贡湖、上山、渔洋山及胥口湾的4个水厂附近水域, 由苏州市水文分局执行巡查任务[5].湖泛巡查规范要求, 巡查船只在巡查区域内来回往复穿梭、蛇行曲折向前巡查及监测.

巡查过程中还设置固定水质测定点位.太湖西岸带八房港至竺山湖区域:设置八房港东、陈东港东、官渎港东、新渎港东、社渎港东、师渎港东、符渎港东、沙塘港东、雅浦港南、竺山湖中及竺山湖北这11个巡查监测点, 巡查水域面积为190 km2; 马山月亮湾、梅梁湾及贡湖湾区域:设置月亮湾、拖山、南泉水厂、马山水厂、直湖港东、小湾里水厂、南龚(外)、张桥(外)、锡东水厂和白洋湾水厂这10个巡查监测点, 巡查水域面积为272 km2; 东部湖区:设置白洋湾水厂、镇湖水厂、渔洋山水厂和园区水厂这4个巡查监测点, 巡查水域面积为150 km2.在固定巡测点监测时, 用手持式风速风向仪测定瞬时的风速风向, 用便携式水质仪测定水体温度(WT)、溶解氧(DO), 并采集表层水样, 带回实验室测定水体总氮(TN)、总磷(TP)、氨氮(NH4+-N)、化学需氧量(COD)、硫化物(S2-)和溶解性铁(Fe2+)等湖泛相关水质指标.

巡查中对于湖泛的判别, 主要靠目测水色及嗅觉气味等初步判别, 并用便携式水质仪的DO指标进行确认.目测表观有4种情况:①湖泛水体浑浊, 颜色为浆褐色, 与未发生湖泛水体有清晰的界限; ②空气中弥漫着类似下水道的恶臭和硫化氢的气味; ③湖泛发生湖水面时有气泡冒出; ④水面发现湖鳅等死鱼现象.当初步判断某水域有湖泛发生时, 立即进一步测定水体溶解氧加以确认.并进一步查看湖泛水域的边界, 根据区域的大致长、宽距离, 估算湖泛面积, 详细测量水体DO的空间和垂向分布, 测定湖泛区的WT等现场特征指标, 并立即上报太湖应急防控部门, 此后每天巡测时, 再次调查该湖泛区的面积和水质状况等, 观测其变化, 直至湖泛消失.

为定量表征每次湖泛的强度(black water event index, BWEI), 采用每场湖泛的平均面积与持续天数的乘积来表示:

式中, ABWE为单次湖泛的平均面积, km2, TBWE为该场湖泛持续天数, d, BWEI的单位为km2·d.

1.2 太湖水质与水文气象观测

湖泛的出现是偶然的.因此湖泛现场水质、水文、气象监测记录可能会存在不连续等问题.为此, 调取了中国科学院太湖湖泊生态系统观测站(TLLER)的逐日水温、水位(WL)及降雨量(RF)等观测场记录数据.此外, TLLER每月在蓝藻水华易发的北太湖布设14个固定采样点, 每月采集水样分析浮游植物叶绿素a(Chla)、TN、TP、NH4+-N、生物需氧量(BOD5)、高锰酸盐指数、电导率(EC)及硫酸根离子(SO42-)等相关水质指标, 样点布设及指标分析方法见文献[8].

1.3 太湖蓝藻水华的年度强度

TLLER的遥感组采用每月多幅遥感影像, 获得全年太湖蓝藻水华面积平均值, 记作ABL.这是一个定量反映全湖水华分布及其动态变化的指标.蓝藻水华面积及其叶绿素a(ChlaRS)的遥感反演方法见文献[9], 其中, 遥感影像数据采用MODIS aqua L-0数据, 2009年1月1日至2018年12月31日的3 000余景的逐日影像从NASA的Goddard Space Flight Center Web网站下载(http://oceancolor.gsfc.nasa.gov), 经目视遴选, 共获得太湖湖面无云可用的1 000余景影像, 实施水体叶绿素a浓度反演, 估算ABL, 具体方法见文献[10].

TLLER在每年2、5、8和11月4次在全太湖32个观测点采样(分上、中和下这3层采样后混合水样)分析获得的水体Chla浓度, 按照泰森多边形的方式获得各个点位的水量权重, 加权计算出全年全湖水体Chla平均值, 计作ChlaLB, 具体布点及加权平均值计算方法见文献[11].

太湖年度蓝藻水华强度(algal bloom index, ABI)的表征, 采用遥感获得的蓝藻水华面积参数ABL及现场调查获得的水体分层浮游植物叶绿素a的ChlaLB年度均值综合获得.

式中, ABL为当年蓝藻水华解译的蓝藻水华平均面积, km2, ChlaLB为太湖站32个点位4个季度3层水体中浮游植物叶绿素a的年度加权平均值, μg·L-1.

2 结果与分析 2.1 2009~2018年太湖湖泛发生的基本特征

2009~2018年的10 a中, 共发现湖泛事件75起, 年发生次数3~17次不等.其中2017年发生次数最多, 为17次; 2012年发生次数为第2多, 共发生12次; 而2014年发生次数最少, 只有3次(图 2).就湖泛强度而言, 2017年和2018年显著强于其他年份(图 2), 两年平均值为185.45 km2·d, 而其余8 a的平均值仅有16.75 km2·d, 相差超过10倍.单次湖泛平均持续时间为3~4 d, 最长的达16 d之久, 最短的只有1 d.单次湖泛发生面积在0.01~9.2 km2之间, 其中湖西区湖泛平均面积最大, 达2.1 km2, 贡湖湾湖泛平均面积相对最小, 为0.2 km2.

图 2 宜兴水域(湖西区)及无锡水域(梅梁湾及贡湖湾之和)湖泛发生次数及强度 Fig. 2 Annual times and intensity of black water events in Lake Taihu from 2009 to 2018

就发生时间而言, 湖泛集中发生在每年的5~9月, 尤其以6~7月最多. 75次湖泛的开始日期均在5~9月, 其中5月发生12次、6月25次、7月23次、8月9次和9月6次.除了2012年9月30日的湖泛延迟到了10月7日才结束之外, 其余湖泛都在9月就已经结束.发生湖泛的最早时间为5月11日.

空间上, 隶属宜兴市的湖西区沿岸带是湖泛的高发水域, 共计41次, 超过总次数的一半, 主要发生在北至沙塘港和南至八房港的25 km长岸线的各港入湖水域.马山半岛南部的月亮湾暴发频次也很高, 10年发生了16次.梅梁湾的发生频次也较高, 10年共发生11次, 集中发生在北部的渔港套闸及西北角的闾江口附近, 偶尔有湖泛发生在三山岛北的杨湾及马山港入梅梁湾的河口区.贡湖湾北岸带水域10年共发生7次湖泛, 主要发生在小溪港、黄泥田港、庙港、许仙港和壬子港等主要入湖河口区.几乎所有的湖泛都发生在靠近入湖河道的入口附近水域.

根据湖泛形成的生源物质的不同, 可对75场湖泛事件细分为草源性湖泛和藻源性湖泛两种类型.顾名思义, 草源性湖泛是因大量水草腐烂而引发的, 腐烂物质主要为死亡水草残体; 而藻源性湖泛是因大量蓝藻聚积、死亡、腐烂而造成的, 腐烂物质主要为死亡藻体. 10 a来的75次湖泛中绝大多数为藻源性湖泛, 仅有5次为草源性湖泛, 其中4次发生在贡湖湾北岸的小溪港、许仙港、大溪港一带, 另有1次发生在梅梁湾北岸沿岸的杨湾至三山岛一带水域, 与这些属于春末夏初大量菹草的季节性死亡有关.草源性湖泛发生的强度都不大, 平均面积为0.25 km2, 最大面积0.80 km2, 持续时间一般3~5 d.

根据对湖泛事件的连续跟踪分析, 在藻源性湖泛中, 针对藻的来源及其发生的地点和持续时间, 又可细分为港源型、迁移型和原发型这3种类型.港源型湖泛是指由于降雨形成的径流将港口内堆积的大量死亡蓝藻冲至湖泊, 在港口附近水域腐烂并诱发湖泛, 湖泛区湖-河联系紧密, 呈现由河至湖扩散的特征.这种湖泛面积较小, 持续时间较短.迁移型湖泛是指在持续高温作用下, 离岸的大湖面水域蓝藻水华先期大暴发, 然后在连续多日东南风或东风的驱动下, 湖面上的蓝藻向岸边或湖湾浅水区聚积, 进而大量死亡、下沉并腐烂而形成湖泛, 多发生在太湖西岸或北部湖湾.这种面积较大, 持续时间较长.原发型湖泛是指发生湖泛区前期未见外围蓝藻水华的明显移动, 湖泛区本地藻类水华大量持续暴发、堆积, 进而在持续高温和弱风的气象条件下, 蓝藻大量死亡所致.这种湖泛的面积一般也不大, 但持续时间长.

2.2 湖泛发生的水温与水质特征

2009~2018年, 湖泛发生首日, 湖泛发生核心区的水温变化范围为20.2~34.2℃, 这表明, 水温20℃以上可能是湖泛发生的临界温度.湖泛发生首日, 湖泛水域的DO均值为1.95 mg·L-1, NH4+-N含量均值为2.70 mg·L-1, COD的均值为68.0 mg·L-1, TP的均值为0.478 mg·L-1, TN的均值为5.99 mg·L-1, Fe2+离子的浓度均值为0.53 mg·L-1.与TLLER的5~9月多年监测均值相比, 氨氮高出湖区平均值的6倍, 而总氮和总磷高出背景2~3倍(表 1).鉴于湖泛区本来就处于沿岸带河口区, 营养盐本就高于湖心区, 这样的营养盐浓度差别并不是特别高.

表 1 湖泛发生首日核心区水温及水质特征 Table 1 Water temperatures and water parameters of the first day of black water events

各湖区首日湖泛核心区的水质情况存在差异.其中湖泛首日的湖西区水质最差, 而贡湖相对情况较好; 这与太湖水体营养盐浓度本就存在的空间差异相符.湖泛发生首日除水温与溶解氧外的其它水质指标浓度最大值均出现在湖西, 而浓度最小值大部分出现在贡湖.据金颖薇等的调查, 无论是氮、磷含量, 还是浮游植物叶绿素a浓度, 2013年夏季太湖西北区水体TN、TP和Chla浓度分别为3.96 mg·L-1、0.259 mg·L-1和47.39 μg·L-1, 而梅梁湾区域则分别为2.34 mg·L-1、0.089 mg·L-1和21.49 μg·L-1 [12], 西北沿岸带的3项指标本就比梅梁湾高出不少.

随着湖泛的发展, 整个湖泛过程期间水域营养盐等水质指标明显高于未发生湖泛的一般水体.以湖西区的3种湖泛类型整个湖泛期间水质监测结果为例, 湖泛核心区水域在整个湖泛期间多日平均的NH4+-N、COD、TP和TN分别为《地表水环境质量标准》(GB 3838-2002)Ⅲ类标准上限值的35.5、56.4、207和119倍(表 2).这说明, 随着湖泛的持续, 对周边水域的水质产生了极大地影响.其中, 迁移性湖泛对水体有机质、氮、磷的影响最大, 但原发性湖泛的水体硫化物浓度最高, 恶臭问题最严重.此外, 从表 2中也可以看出, 各指标的变幅较大, 显示出湖泛区水质的空间异质性较大.

表 2 湖西区3种类型湖泛核心区水质指标/mg·L-1 Table 2 Water quality in different types of black water event zones in the western coast of Lake Taihu/mg·L-1

2.3 湖泛发生的气象水文条件

刘俊杰等[4]对太湖湖泛的前期10 d内气象条件进行了分析表明, 气象要素是湖泛形成的关键促发因素, 往往决定了湖泛发生的时间、地点及强度.本研究在更大样本量的基础上, 进一步梳理了太湖湖泛形成前5 d的综合气象水文条件特征.其中气象水文要素主要考虑的是气温、水温、风向风速、日照、气压和水位.根据统计学原理, 考虑到不同次湖泛发生时的观测数据具有很大的差异性和不确定性, 数据复杂且多处于偏态分布, 因此采用SPSS 22.0软件, 采用四分位法确定各气象水文条件阈值, 即以下四分位数至上四分位数区间值作为湖泛发生各气象水文条件的预警阈值.不同湖泛类型的环境条件统计如表 3.

表 3 不同湖泛类型的发生前5 d的环境因素 Table 3 Meteorological and hydrological conditions of in five days ahead of the black water events

2.4 北太湖近10 a来的水质和藻情变化

2009~2018年太湖藻情总体呈加重趋势, 但是年际之间差别甚大.其中2017年蓝藻水华大暴发, 异常严重, ABI达到10 a中的最大(图 3).就蓝藻水华强度指标ABI而言, 10 a平均值为3 644 km2·μg·L-1.而2017年则达到11 858 km2·μg·L-1, 是多年均值的3倍多, 是2009~2014年均值(2 309 km2·μg·L-1)的5倍多.同时, 夏季北太湖14个点位的浮游植物叶绿素a月均值(ChlaN)及全湖加权平均值(ChlaLB)也呈类似的趋势, 显示出2017年的异常高值.ChlaN与ChlaLB的变化趋势有所不同.比如5~9月ChlaN在2015、2016和2018年均表现出明显的高值, 但是全湖全年平均值ChlaLB的增高不明显.

图 3 2009~2018年太湖水华强度及浮游植物叶绿素a年均值 Fig. 3 Annual algal bloom intensity and chlorophyll a concentrations of Lake Taihu from 2009 to 2018

根据TLLER在每年5~9月北太湖14个监测点获得的水体主要水质指标如表 4.总体上TN呈下降的趋势, 然而2017年北太湖TN出现了明显的反弹.北太湖TP在湖泛季均较高, 但是2017年也高于其余9年均值52%.总体上, 5~9月北太湖的水体氨氮浓度不高, 10 a的均值为0.42 mg·L-1.而BOD5则在2016年达到最大值, 2017年也较高.同样反映水体有机污染状况的指标高锰酸盐指数与BOD5的变化非常一致.而反映水体电解质状况的指标EC变化也颇大, 最高值出现在2014年, 为590 μS·cm-1, 最低值出现在2016年, 为386μS·cm-1. 2009~2018年EC、SO42-受降雨量影响明显, 均与RF呈反相关.EC在2016年出现最低值, 与降雨量最大值吻合.SO42-是太湖水体电解质的主要贡献者之一[13].北太湖5~9月水体SO42-年均值变化也较大, 与EC呈显著的正相关.

表 4 2009~2018年北太湖5~9月水体主要水质指标及水文气象指标均值 Table 4 Average values of water quality parameters and meteorological and hydrological conditions in north Taihu from May to September

基于太湖站逐日观测, 北太湖的水温值5~9月的多年均值为27.4℃.其中最大值出现在2018年, 为28.5℃.第二高值出现在2017年, 也达到了28.0℃.而最低值出现在2014年, 为26.5℃.北太湖每年5~9月的水位日均值介于3.214~3.729 m之间, 年际之间相差最大为0.515 m. 2016年的洪水年显著高于其他年份.而2015年也较高, 与当年的总体降雨量偏大一致.

年与年之间5~9月的总降雨量差别比较大.发生大洪水的2016年5~9月降雨量达1 059.4 mm, 而较枯的年份2010年仅有341.0 mm, 相差3倍以上, 显示出较大的年际波动性.

3 讨论 3.1 影响太湖湖泛强度年际变化的主要因素

由于湖泛的形成中, 犹如有机肥的堆肥过程, 需要一定的化学反应时间[14], 因此, 湖泛的发生看似具有较大的偶然性, 关于湖泛形成条件的研究, 大都关注在短期气象、水文和物质基础方面[7, 15, 16].刘俊杰等[4]在分析之前9 a的湖泛发生规律时, 初步发现了水位的年变幅与年际尺度的湖泛强度之间具有一定的关联性.然而, 其结论主要受2017年洪水的影响, 具有一定的偶然性.在年际尺度上, 何种条件下湖泛事件更容易发生, 仍未有较为清晰地认识.

近10年太湖湖泛的发生情况存在较大的变动性.从湖泛的成因上, 温度、风浪扰动、水华强度、水质污染程度和底泥均可能对湖泛的形成和强度产生影响.从图 2可以看出, 早期开展湖泛发生机制、条件研究的2009~2016年, 太湖的湖泛强度还是相对较低的.而2017~2018年湖泛强度反而特别强.为探讨影响年际湖泛强度的因子, 利用近10年的数据, 对可能影响湖泛强度的多种因子(5~9月均值)与湖泛年际强度进行了相关性分析, 如表 5.

表 5 太湖年度湖泛强度与相关环境因子之间的相关性1) Table 5 Person correlation index between the black water event index and environmental factors

表 5中可以看出, 环境因子中, 与湖泛强度BWEI显著相关的只有蓝藻水华强度和水温这2个指标.这表明, 水华情势是影响湖泛情势的首要因素.这一结果, 与Yu等在梅梁湾做的模拟实验结果的认知[17], 以及Zhou等[18]在太湖的现场观测认知基本一致, 即水体叶绿素a浓度, 或者说蓝藻水华强度, 是影响湖泛的第一要素. 5~9月的平均水温与湖泛强度显著相关(图 4), 表明高温条件是诱发湖泛的最主要气象条件, 这与湖泛形成中的微生物作用有关(图 4).蓝藻水华的野外分解模拟实验表明, 蓝藻水华从细胞活体到死亡分解, 大量的微生物参与其中[19], 而微生物的活性与温度密切关联.除了统计上显著相关的2个因子之外, 反映藻类生物量的Chla、反映水体易腐烂有机物状况的高锰酸盐指数, 以及水体TP, 也与BWEI具有较高的关联性.之所以它们也与湖泛有较明显的统计相关, 可能是这些指标与水体中藻类生物量密切相关, 高的氮磷浓度、有机物浓度, 往往是藻类生物量升高的一个结果, 与蓝藻水华严重程度密切关联.

图 4 太湖年度湖泛强度指数与水华指数及水温的关系 Fig. 4 Relationship between BWEI and ABI or water temperature in Lake Taihu from 2009 to 2018

表 5可以看出, 太湖的湖泛强度与各种营养盐的关系大都不好.这说明, 太湖水体营养盐浓度高低对湖泛强度的影响是十分有限.Xu等[20]的观测与模拟实验分析表明, 目前太湖水体营养盐的浓度大都明显高于能够限制蓝藻水华生长的临界阈值, 也即太湖蓝藻水华强度基本不受水体氮和磷等营养盐浓度的限制.因此, 太湖蓝藻水华的年际变化受营养盐浓度影响弱, 水体氮、磷浓度常常是蓝藻水华变化的一个结果, 而非原因[8].对于蓝藻水华在特定环境条件下发生腐烂分解等化学和生物作用之后产生的湖泛而言, 受营养盐的影响肯定就更低了.湖泛伴随着水体营养盐浓度的剧增, 然而, 也应该是湖泛的结果, 而非原因.北太湖5~9月水体TP和BOD5等反映营养状况的水质参数受藻情影响很大, 也说明了这一点.

在湖泛的形成中, 污染底泥的参与也极为重要[5, 6].据Chen等[21]的模拟实验研究, 实施底泥疏浚后藻-泥培养系统中, 挥发性有机硫化合物(VOSCs), 如二甲基二硫(DMDS)和二甲基三硫(DMTS)等的释放强度明显低于对照.本研究发现几乎所有的湖泛都发生在靠近入湖河口区域, 一个重要的原因是这些区域易于高有机质含量的污染底泥蓄积, 加上岸线的变化和来自河流和湖泊中污染物质的交汇, 这些区域既容易产生底泥再悬浮, 也容易形成水质突变, 打破原有生物地球化学平衡, 发生极端化学反应.由于本文关于底泥方面的调查资料不足, 并未能将底泥污染状况的内容置于其中, 底泥的污染程度及其再悬浮过程对湖泛的影响还有待进一步考虑.

由于太湖的蓝藻水华情势较大程度上受年际气温、风速等气候波动的影响[9], 而5~9月的高温状况又对当年湖泛强度产生直接的影响, 因此, 太湖的年度湖泛强度在较大程度上受到气候波动的影响.近年来, 随着高温热浪等极端气象事件的增多[22], 高温热浪等极端天气事件大大增加了湖体有机质暴发式分解的几率[23], 太湖的湖泛风险有增高的可能.

3.2 降低太湖湖泛风险的途径分析

控制蓝藻水华强度是消除湖泛灾害风险的根本途径.从前面的分析看, 蓝藻水华强度与湖泛强度密切相关. 2017年与2018年湖泛强度的增加, 其根本原因是太湖蓝藻水华强度增加的必然结果.根据对湖泛发生条件分析, 气象水文要素是湖泛发生的诱因, 物质要素则为湖泛的发生提供了基础[1], 因此, 湖泛防控措施主要是以控制蓝藻水华强度为核心的各种控制工程, 包括滨岸带的生态修复及外源营养盐的控制.可行的措施包括:①努力大幅度减少营养盐的外源入湖量, 使之降至能显著影响蓝藻水华强度的浓度区间.入湖河道是湖泊污染负荷的重要来源[24], 要继续推进重点区域治理、小流域综合整治及水系畅通工程, 推动并强化入湖河流两岸生态保护红线划分与管理, 全面深化河湖长制, 使河湖长有名有实在责, 加强入湖河道的水质管控.②控制内源以降低湖体营养盐浓度.湖体营养主要来自外源与内源, 外源主要是通湖河港输入, 而比较可行的控制内源的工程措施主要是蓝藻打捞和生态清淤.构建专业化打捞队伍, 持续开展蓝藻打捞, 固定打捞与机动打捞相结合, 深化藻泥的无害化处置、资源化利用, 强化河荡入湖处及湖泛易发区的生态清淤工作.③加强生态修复, 包括滨岸带的生态修复及浅水区的沉水植被恢复等.作为一个大型浅水湖泊, 沉水植物存在与否, 底泥再悬浮的营养盐通量可差数十倍[25].没有足够的沉水植被等吸收水体营养盐、抑制沉积物再悬浮等湖体生态系统的自我净化能力提升, 想要把水体营养盐控制到抑制蓝藻水华发生的水平是极其困难的.因此, 应当大力推进环湖滨岸带、东部湖湾等传统上水草茂盛的水域尽快实施生态修复, 提升湖体生态控藻的能力, 从根本上降低湖泛发生的物质基础.

降低太湖湖泛风险, 还需要在严密监控的基础上及时开展预防蓝藻水华长时间聚集腐烂的工程措施.加强对入湖河口区蓝藻水华堆积的机械打捞、曝气和动力驱散等人工干预是降低湖泛发生的有效途径.从前期的研究成果看, 无论是蓝藻水华物质基础, 还是风速和气温等气象条件, 以及湖底流泥等湖泛形成的环境条件, 本质上是有利于蓝藻水华在特定区域稳定、长期堆积腐烂.从过去10a的75场湖泛调查发现, 几乎所有的湖泛都发生在太湖西北及北部易于堆积蓝藻水华的河口区.而这些河口区大都设置了蓝藻水华打捞点.因此, 加强打捞点和河口区的藻情监控和及时处置, 是降低湖泛发生风险的有效途径.此外, 从底泥污染与湖泛形成的关系看, 实施入湖河口区域的污染底泥清淤, 对降低湖泛风险应该有用.从太湖多年的湖泛防控工作实践看, 蓝藻水华的离岸打捞、调水引流、湖底藻体有机质清扫、高效增氧曝气船干预等河口区湖泛干预技术的实施, 已经在太湖的湖泛灾害防控中发挥了相当大的作用.

4 结论

(1) 从太湖10 a的湖泛巡查调查看, 太湖湖泛均发生在水温20℃以上的5~9月, 基本都发生在污染较重的入湖河口区附近, 绝大多数湖泛是蓝藻水华堆积导致的藻源性湖泛, 只有5次(占总数的7%)为草源性湖泛.而藻源性湖泛因藻的来源、发生地点和持续时间的不同, 又可分为港源型、迁移型和原发型3种.

(2) 藻情和温度是影响湖泛强度年际差别的主要原因.年度蓝藻水华越严重, 湖泛发生的风险就越大, 总体年度湖泛的强度就越大. 5~9月的平均水温越高, 湖泛发生的强度也越大.由于太湖蓝藻水华情势本身也受到气候条件变化的明显影响, 所以, 太湖的年度湖泛强度在很大程度上受到气候波动的影响.

(3) 控制太湖蓝藻水华是降低湖泛强度的根本途径, 而加强河口区蓝藻水华堆积的监控与工程干预, 加强引流调水、离岸打捞、应急打捞与曝气等河口区蓝藻水华堆积干预, 是降低湖泛灾害的有效途径.

致谢: 感谢江苏省水文水资源勘测局常州分局、无锡分局和苏州分局十年来的艰苦湖泛巡测; 感谢施坤和李娜提供十年的蓝藻水华遥感监测数据.

参考文献
[1] 陆桂华, 马倩. 太湖水域"湖泛"及其成因研究[J]. 水科学进展, 2009, 20(3): 438-442.
Lu G H, Ma Q. Analysis on the causes of forming black water cluster in Taihu Lake[J]. Advances in Water Science, 2009, 20(3): 43-442.
[2] Qin B Q, Zhu G W, Gao G, et al. A drinking water crisis in Lake Taihu, China: linkage to climatic variability and lake management[J]. Environmental Management, 2010, 45(1): 105-112. DOI:10.1007/s00267-009-9393-6
[3] 陆桂华, 张建华. 太湖蓝藻监测处置与湖泛成因[M]. 北京: 科学出版社, 2011: 219-261.
[4] 刘俊杰, 陆隽, 朱广伟, 等. 2009-2017年太湖湖泛发生特征及其影响因素[J]. 湖泊科学, 2018, 30(5): 1196-1205.
Liu J J, Lu J, Zhu G W, et al. Occurrence characteristics of black patch events and their influencing factors in Lake Taihu during 2009 and 2017[J]. Journal of Lake Sciences, 2018, 30(5): 1196-1205.
[5] 陆桂华, 马倩. 2009年太湖水域"湖泛"监测与分析[J]. 湖泊科学, 2010, 22(4): 481-487.
Lu G H, Ma Q. Monitoring and analysis on "Black Water Aggregation" in Lake Taihu, 2009[J]. Journal of Lake Sciences, 2010, 22(4): 481-487.
[6] 范成新. 太湖湖泛形成研究进展与展望[J]. 湖泊科学, 2015, 27(4): 553-566.
Fan C X. Progress and prospect in formation of black bloom in Lake Taihu: a review[J]. Journal of Lake Sciences, 2015, 27(4): 553-566.
[7] 王成林, 张咏, 张宁红, 等. 太湖藻源性"湖泛"形成机制的气象因素分析[J]. 环境科学, 2011, 32(2): 401-408.
Wang C L, Zhang Y, Zhang N H, et al. Analysis of meteorological factors of forming feculent and anaerobic water aggregation induced by algal bloom in Taihu Lake[J]. Environmental Science, 2011, 32(2): 401-408.
[8] 朱广伟, 秦伯强, 张运林, 等. 2005-2017年北部太湖水体叶绿素a和营养盐变化及影响因素[J]. 湖泊科学, 2018, 30(2): 279-295.
Zhu G W, Qin B Q, Zhang Y L, et al. Variation and driving factors of nutrients and Chlorophyll-a concentrations in northern region of Lake Taihu, China, 2005-2017[J]. Journal of Lake Sciences, 2018, 30(2): 279-295.
[9] Shi K, Zhang Y L, Zhang Y B, et al. Phenology of phytoplankton blooms in a trophic lake observed from long-term MODIS data[J]. Environmental Science & Technology, 2019, 53(5): 2324-2331.
[10] 朱广伟, 施坤, 李未, 等. 太湖蓝藻水华的年度情势预测方法探讨[J]. 湖泊科学, 2020, 32(5): 1421-1431.
Zhu G W, Shi K, Li W, et al. Seasonal forecast method of cyanobacterial bloom intensity in eutrophic Lake Taihu, China[J]. Journal of Lake Sciences, 2020, 32(5): 1421-1431.
[11] 朱广伟, 邹伟, 国超旋, 等. 太湖水体磷浓度与赋存量长期变化(2005-2018年)及其对未来磷控制目标管理的启示[J]. 湖泊科学, 2020, 32(1): 21-35.
Zhu G W, Zou W, Guo C X, et al. Long-term variations of phosphorus concentration and capacity in Lake Taihu, 2005-2018: implications for future phosphorus reduction target management[J]. Journal of Lake Sciences, 2020, 32(1): 21-35.
[12] 金颖薇, 朱广伟, 许海, 等. 太湖水华期营养盐空间分异特征与赋存量估算[J]. 环境科学, 2015, 36(3): 936-945.
Jin Y W, Zhu G W, Xu H, et al. Spatial distribution pattern and stock estimation of nutrients during bloom season in Lake Taihu[J]. Environmental Science, 2015, 36(3): 936-945.
[13] Wu T H, Zhu G W, Zhu M Y, et al. Use of conductivity to indicate long-term changes in pollution processes in Lake Taihu, a large shallow lake[J]. Environmental Science and Pollution Research, 2020, 27(17): 21376-21385. DOI:10.1007/s11356-020-08590-x
[14] 申秋实, 邵世光, 王兆德, 等. 太湖月亮湾湖泛发生过程模拟及水土物化性质的响应[J]. 水科学进展, 2011, 22(5): 710-719.
Shen Q S, Shao S G, Wang Z D, et al. Simulation of black bloom in Moon Bay of Lake Taihu and physical and chemical responses of water and sediment[J]. Advances in Water Science, 2011, 22(5): 710-719.
[15] 王成林, 张宁红, 张咏, 等. 基于气象条件的太湖湖泛预警研究[J]. 环境监控与预警, 2010, 2(5): 1-4.
Wang C L, Zhang N H, Zhang Y, et al. The research on warning of black water cluster in Taihu Lake based on meteorological factors[J]. Environmental Monitoring and Forewarning, 2010, 2(5): 1-4.
[16] 李未, 秦伯强, 张运林, 等. 富营养化浅水湖泊藻源性湖泛的短期数值预报方法——以太湖为例[J]. 湖泊科学, 2016, 28(4): 701-709.
Li W, Qin B Q, Zhang Y L, et al. Numerical forecasting of short-term algae-induced black bloom in eutrophic shallow lake: a case study of Lake Taihu[J]. Journal of Lake Sciences, 2016, 28(4): 701-709.
[17] Yu D Z, Xie P, Zeng C, et al. In situ enclosure experiments on the occurrence, development and decline of black bloom and the dynamics of its associated taste and odor compounds[J]. Ecological Engineering, 2016, 87: 246-253. DOI:10.1016/j.ecoleng.2015.11.039
[18] Zhou Y Q, Jeppesen E, Zhang Y L, et al. Chromophoric dissolved organic matter of black waters in a highly eutrophic Chinese lake: freshly produced from algal scums?[J]. Journal of Hazardous Materials, 2015, 299: 222-230. DOI:10.1016/j.jhazmat.2015.06.024
[19] Shao K Q, Zhang L, Wang Y P, et al. The responses of the taxa composition of particle-attached bacterial community to the decomposition of Microcystis blooms[J]. Science of the Total Environment, 2014, 488-489: 236-242. DOI:10.1016/j.scitotenv.2014.04.101
[20] Xu H, Paerl H, Qin B, et al. Determining critical nutrient thresholds needed to control harmful cyanobacterial blooms in eutrophic Lake Taihu, China[J]. Environmental Science & Technology, 2015, 49(2): 1051-1059.
[21] Chen C, Zhong J C, Yu J H, et al. Optimum dredging time for inhibition and prevention of algae-induced black blooms in Lake Taihu, China[J]. Environmental Science and Pollution Research, 2016, 23(14): 14636-14645. DOI:10.1007/s11356-016-6627-z
[22] Lopez H, West R, Dong S F, et al. Early emergence of anthropogenically forced heat waves in the western United States and Great Lakes[J]. Nature Climate Change, 2018, 8(5): 414-420. DOI:10.1038/s41558-018-0116-y
[23] Bartosiewicz M, Laurion I, Clayer F, et al. Heat-wave effects on oxygen, nutrients, and phytoplankton can alter global warming potential of gases emitted from a small shallow lake[J]. Environmental Science & Technology, 2016, 50(12): 6267-6275.
[24] 陈洁, 朱广伟, 许海, 等. 不同雨强对太湖河网区河道入湖营养盐负荷影响[J]. 环境科学, 2019, 40(11): 4924-4931.
Chen J, Zhu G W, Xu H, et al. Influence of rainfall intensity on the nutrient loading from an inflowing river in the plain river network of the Taihu catchment[J]. Environmental Science, 2019, 40(11): 4924-4931.
[25] Zhu M Y, Zhu G W, Nurminen L, et al. The influence of macrophytes on sediment resuspension and the effect of associated nutrients in a shallow and large lake (Lake Taihu, China)[J]. PLoS One, 2015, 10(6). DOI:10.1371/journal.pone.0127915