随着经济发展与人类活动加剧，水体富营养化已成为世界性的水环境问题.近年来，世界各大淡水湖和我国滇池、太湖、巢湖等内陆淡水湖多次发生大规模蓝藻集聚现象，严重破坏水生态系统稳定，威胁当地生活生产用水安全.叶绿素a(Chl.a)是藻类进行光合作用的重要色素，其在水中的浓度通常用于表征水中藻类总体含量以及评判水体的营养状态^[1].

国内外已开展大量淡水湖泊Chl.a浓度时空动态变化相关研究. Ghadouani等研究伊利湖的藻类和营养盐分布，表明Chl.a高浓度区主要集中于西部湖湾及北部沿岸地区，且与营养盐分布密切相关^[2]. Jindal等对Prashar湖藻类随时间动态变化过程进行连续观测，发现Chl.a浓度曲线呈双峰型特征，分别在5月和9月达到峰值^[3]. Wu等针对鄱阳湖藻类时空变化进行分析，得出结论：水下光照强度和透明度是限制藻类生长的关键因子，同时南部湖区Chl.a浓度高于北部湖区主要受水温和营养盐影响^[4]. Yang等对巢湖Chl.a浓度和总氮(TN)、总磷(TP)等营养盐浓度进行空间分布和季节性变化特征研究，总结出巢湖藻类生物量存在明显先下降后上升的季节变化特点，同时东、西部湖区存在明显空间差异，西部湖区营养盐及Chl.a浓度高于东部主要受上游汇水影响^[5-6]. Wu等通过对太湖短期连续观测表明，Chl.a浓度水平分布模式与风生流关系紧密^[7].钱昊钟等分析太湖湖心、东南湖区与其他湖区Chl.a浓度差异，得出TP、水温和溶解氧是主要限制性因子^[8].朱晶晶等发现滇池Chl.a浓度总体呈周期性波动缓慢上升趋势，4-11月上升达到峰值；空间分布上，滇池Chl.a浓度高值区多出现在靠近城镇的人口密集区域^[9-10]. Huang等通过分析滇池水体富营养化与人类活动的响应关系，得出含高浓度营养盐的污水排放对滇池富营养化起到关键性作用^[11].郭劲松等研究三峡小江回水区Chl.a浓度季节变化特征为夏、秋季高，春、冬季低^{[1, 12]}. Zhang等研究三峡水库蓄水后大宁河的藻类季节变化特征，得出水温和TN/TP是藻类暴发的关键因素^[13].分析水体Chl.a浓度时空变化特征有助于深化对藻类水华发生过程机理的认识，从而为有效开展藻类水华防控治理工作提供科学依据.

洪泽湖是苏北地区居民生活、农业用水、发电以及维护当地生境的重要水源地.同时作为南水北调东线工程重要输水线路和调蓄湖泊，其水质不仅关系到南水北调水质安全，也关系到沿河、沿湖乃至整个淮河流域经济可持续发展^[14].因此全面认识洪泽湖富营养化水平和藻类时空分布特征十分重要.但是目前研究较多围绕洪泽湖富营养化水平定性分析和水质理化因子时空特征分析^[15-16]，Chl.a浓度时空变化特征研究相对较少.本文基于洪泽湖全年水质监测数据，采用方差分析和典范对应分析等统计手段研究洪泽湖藻类时空变化特征，进一步探究其藻类时空异质性原因，为宏观掌控洪泽湖营养水平、制定水质改善措施提供科学依据^[16].

洪泽湖地处淮河中游末端(33°06′~33°40′N，118°10′~118°52′E)，为中国第四大淡水湖泊.湖区跨洪泽、淮阴、泗阳、泗洪和盱眙五县.东岸平直，其余岸线曲折多湾.湖区集水面积为15.6×10⁴ km²，补给系数为99.常年平均水位12.37 m，长65.0 km，最大宽55.0 km，平均宽为24.26 km，最大水深4.37 m，平均水深1.77 m.湖区属北亚热带与南温带的过渡气候，多年平均气温16.3℃，7月平均气温最高为28.8℃.多年平均降水量925.5 mm，年内7-9月降水量约占年降水量的65.5%.湖水来源除大气降水外，主要依靠地表径流，入湖河流主要有8条(淮河、新汴河、老汴河、新濉河、老濉河、徐洪河、怀洪新河和安东河)，其中淮河流量最大^[17].洪泽湖为吞吐型湖泊，换水周期约为35 d^[18].湖流以吞吐流和风生流形式并存，在淮河入湖口水流以扇形扩散.

根据湖盆形态，生境差异，水文、水力特征以及出入湖河流特性等因素，本研究将洪泽湖湖区分为3个子区：北部湖区、东部湖区和西部湖区.北部湖区即成子湖区域，水体流动性差，并且受宿迁市城市尾水的影响，水质总体略差；西部湖区包括溧河洼区域，地势较高且水深较浅，水生植物分布面积相对较广，生长茂盛^[19]，水质整体较好；东部湖区包括淮河入湖口、过水通道及蒋坝闸湾，湖区水体流动性好，湖流运动较为剧烈.根据上述湖区形态特征，本研究在全湖共布设8个水质监测站点.北部湖区设有颜圩、高湖2个监测点；西部湖区有临淮、洪泽湖区(宿迁南)2个监测点；东部湖区布设洪泽湖区(淮安北)、洪泽湖区(淮安西)、洪泽湖区(淮安东)和洪泽湖区(淮安南)4个监测点(图 1).

图 1(Fig.1)

图 1 洪泽湖地理位置及水质监测点分布 Fig.1 Location of Lake Hongze and its water sampling sites

本研究选取2012年12月-2013年11月采集的洪泽湖月水质监测数据样本，运用统计软件IBM SPSS Statictis 19和Canoco for Windows 4.5对数据进行统计分析.数据是由江苏省水文水资源勘测局提供的常规监测资料.监测项目包括5项：化学需氧量(COD)、溶解氧(DO)、TP、TN和叶绿素a(Chl.a)浓度.其中COD浓度用高锰酸盐指数法测定，DO浓度用电化学探头法测定，TP和TN浓度分别用钼酸铵分光光度法和碱性过硫酸钾消解紫外分光光度法测定，Chl.a浓度用分光光度法测定.此外，日平均气温(T)数据来自国家气象中心的盱眙站(站点编号：58138).根据研究区全年气温变化规律和地处四季分明的温带地区，综合参考传统气候统计法，将研究区季节划分为：春季为3-5月，夏季为6-8月，秋季为9-11月，冬季为12至次年2月.

洪泽湖Chl.a浓度的季节变化规律在不同湖区有所差异(图 2).北、西部湖区Chl.a浓度在冬末(1月)出现明显的上升趋势，平均浓度达22.73±2.58 μg/L.进入春季(3-5月)，Chl.a浓度稍有波动(13.45±4.91 μg/L)后，变化趋于平缓(17.95±4.94 μg/L).春季湖区平均气温达到15.54℃，日光照时长为6.45 h，光照气温条件适宜.同时春耕施肥等面源污染带来营养盐汇入，营养盐水平较高，水体中TN平均浓度为1.35 mg/L，TP平均浓度为0.045 mg/L.藻类生理、生化活性恢复，生长速率加快，种群生物量快速增加后趋于稳定.夏季(6-8月)洪泽湖区平均气温达到26.99℃，光照时长为5.72 h.虽然光照充足、温度适宜，但大规模藻类生长为浮游动物及其他大型鱼类牧食提供便利，同时营养盐也被快速消耗，遏制藻类数目的进一步增加.造成不同站点Chl.a浓度先后回落至低浓度(3.38±1.59 μg/L)后才逐步上升.进入秋季，藻类种间竞争得到缓解，优势种得以快速增殖，使Chl.a浓度在9月左右达到峰值(36.28±45.92 μg/L)，其中颜圩站点浓度高达104.9 μg/L.随后10-11月Chl.a浓度受气温下降、光照减弱影响，呈波动下降状态.冬季湖区气温降至2.78℃，藻类代谢基本停止，进入衰亡、休眠期^[20]，因此Chl.a浓度逐步下降至最低浓度(1.55±1.98 μg/L).

图 2(Fig.2)

图 2 洪泽湖湖区Chl.a浓度的季节变化曲线 Fig.2 Seasonal variation of Chl.a concentrations in northern, western and eastern zones in Lake Hongze

东部湖区Chl.a浓度季节变化曲线呈“双峰型”(图 2)，分别在3月和8月达到浓度峰值.冬末春初(1-3月)洪泽湖的适宜光照和温度以及较高营养盐水平(TN浓度为4.30 mg/L，TP浓度为0.079 mg/L)使得藻类快速从底泥中复苏，进行大量增殖. Chl.a浓度持续上升直至第1个峰值浓度(13.45±4.89 μg/L).随后Chl.a浓度表现为迅速下降状态，且至峰谷低浓度(1.68±0.67 μg/L).研究表明在大型浅水湖泊中，水动力对浮游生物的数量、分布影响十分明显.如果水体滞留时间变长、流速变缓，浮游植物拥有更长时间生长，在适宜的营养条件下数量往往会显著增长^[21].但由于东区包含过水通道，水体流动性大且流速快.藻类受水流冲刷，不容易集聚，造成Chl.a浓度曲线的下降.夏季(6-8月)藻类在充足光照下达到最大光合效率，藻类生物量成倍增加，使得Chl.a浓度再次上升并于8月达到第2个浓度峰值(25.18±13.45 μg/L).秋季(9-11月)光照减弱且温度持续下降使Chl.a浓度上下波动且总体呈下降趋势.冬季(12月)伴随着藻类大量休眠死亡、群落解体，Chl.a浓度下降为最低值(1.15±0.69 μg/L).

洪泽湖全年Chl.a浓度均值表现为北部湖区(18.50±21.26 μg/L)＞西部湖区(9.8±5.94 μg/L)＞东部湖区(7.68±7.96 μg/L).单因素方差分析(ANOVA，表 1)结果显示3个湖区的Chl.a浓度之间在春季存在显著差异(P＜0.05)，其余季节差异不显著.

表 1(Tab.1)

表 1 洪泽湖不同湖区Chl.a浓度单因素方差分析 Tab.1 Results of ANOVA based on Chl.a concentration of different zones in Lake Hongze 季节 ANOVA Chl.a/(μg/L) 营养状态 F P 均值 变化范围 春季 14.32 0.000 11.26±7.73 0.9~24.7 富营养型 夏季 1.174 0.328 9.75±10.00 1.2~39.7 富营养型 秋季 3.296 0.057 14.44±21.21 2.2~104.9 富营养型 冬季 4.689 0.021 8.25±8.29 0.2~24.8 富营养型

表 1 洪泽湖不同湖区Chl.a浓度单因素方差分析 Tab.1 Results of ANOVA based on Chl.a concentration of different zones in Lake Hongze

春季和冬季东部湖区Chl.a浓度显著低于北、西部湖区(图 3)，春季浓度均值为北部湖区(19.95±3.13 μg/L)＞西部湖区(12.95±3.85 μg/L)＞东部湖区(6.06±6.00 μg/L).冬季浓度均值为北部湖区(14.66±9.39 μg/L)＞西部湖区(10.32±8.40 μg/L)＞东部湖区(4.02±3.41 μg/L). 3个湖区的温度和光照等气象条件相近，但是营养盐水平有所不同.东部湖区的TN和TP浓度均高于北、西部湖区(图 4)，表明营养盐不是春季和冬季东部湖区Chl.a浓度的限制因子.东部湖区Chl.a在春、冬季浓度较低很有可能是由于其过水通道特性，流动性较好.研究表明流速可能对着生藻类生物量具有显著影响，较大流速环境不适合着生藻类的生长^[22].当流速过大时，水流的冲刷作用使藻类的生长、繁殖环境受到破坏，藻类的增长和聚集受到有效抑制^[23-24].北、西部湖区相对封闭，湖流运动相对较弱，藻类较容易进行快速繁殖.同时西、北部湖区周围农耕渔业比东部湖区发达，虽然营养盐浓度低于东区，但是足够为藻类提供生长物质基础，加之适宜的光照温度，藻类生物量可以稳定上升.

图 3(Fig.3)

图 3 洪泽湖不同湖区水体Chl.a浓度的季节变化 Fig.3 Seasonal variation of Chl.a concentrations of waters in different zones of Lake Hongze

图 4(Fig.4)

图 4 洪泽湖不同湖区TN、TP及Chl.a浓度的变化 Fig.4 Changes of TN, TP and Chl.a concentrations in different zones of Lake Hongze

夏季和秋季3个湖区之间Chl.a浓度无显著差异(图 3)，夏季浓度均值表现为东部湖区(13.64±11.84 μg/L)＞北部湖区(7.45±6.86 μg/L)＞西部湖区(5.93±2.37 μg/L).根据Chl.a浓度季节变化规律可知，东部湖区由于春季藻类生物量下降，种内竞争大大减弱，造成在夏季藻类快速增长，而北、西部湖区在夏季则是由于种内、种间竞争加剧造成生物量下降.两者比较，东部湖区Chl.a浓度高于其他2个湖区.秋季洪泽湖Chl.a浓度均值表现为北部湖区(29.75±40.85 μg/L)＞西部湖区(10.04±4.72 μg/L)＞东部湖区(7.84±3.60 μg/L).北部湖区Chl.a浓度远高于其他两个湖区，其中颜圩站点在9月和11月分别出现Chl.a浓度极大值(104.9和48 μg/L).北部湖区由于引水条件较好以及水草丰富，水产养殖业迅猛发展，沿岸多为人工开挖鱼塘和围网养蟹池^{[14, 25]}，造成该区水污染严重，同时雨水冲刷引发农业和养殖业面源污染，造成水质富营养化水平较高(图 4). TN和TP浓度分别达到1.76±1.66和0.14±0.12 mg/L.丰富的营养盐条件促使藻类大规模生长暴发，造成北部湖区局部Chl.a浓度较高.

典范对应分析(canonical correspondence analysis, CCA)提供了分析浮游植物群落与环境因子之间对应关系的工具^[26].本研究进行全年月份和湖区站点Chl.a浓度与水质理化因子(化学需氧量、溶解氧、TP、TN、温度)之间的CCA分析，结果表明，两种CCA分析的轴1特征值均大于其余轴.从浮游植物与环境因子之间的相关系数中，可以看到与轴1和轴2的相关性都较高(表 2).

表 2(Tab.2)

表 2 洪泽湖Chl.a浓度与理化因子CCA分析的统计信息 Tab.2 Summary statistics for the first two axes of CCA performed between Chl.a concentration and water quality factors in Lake Hongze 分析对象 轴1 轴2 月份CCA分析 站点CCA分析 月份CCA分析 站点CCA分析 特征值 0.277 0.255 0.093 0.062 浮游植物与环境因子相关性 0.999 0.959 0.991 0.818 浮游植物与环境因子方差累积百分数 69.1 73.9 92.1 91.9

表 2 洪泽湖Chl.a浓度与理化因子CCA分析的统计信息 Tab.2 Summary statistics for the first two axes of CCA performed between Chl.a concentration and water quality factors in Lake Hongze

洪泽湖12个月的Chl.a浓度都分布在主轴1和主轴2构成的4个象限内，这4个象限基本把研究区4个季节的Chl.a浓度分隔开.第1象限内分布春初(3月)和秋季中旬(10月)的Chl.a浓度，与TN浓度表现为正相关；第2象限分布冬季(1、2月)和春季(4、5月)的Chl.a浓度，与TN/TP表现出正相关；第3象限分布秋季(9、11月)的Chl.a浓度，与COD表现出负相关；第4象限分布夏季(6、7、8月)的Chl.a浓度与TN、TP浓度表现出不同程度的正相关(图 5).以上结果证明，除秋季外，其余季节Chl.a浓度均受不同程度的营养盐浓度影响.氮、磷营养盐是藻类生长不可或缺的元素，同时在一定程度上，温度升高能显著促进藻类对营养盐的吸收利用^[27]，因此夏季Chl.a浓度表现为与TN浓度呈正相关.春、冬两季大多数月份Chl.a浓度与TN/TP呈正相关，研究表明TN/TP对藻类的暴发性生长具有重要意义，是水中浮游植物营养结构特点的重要反映^[28].

图 5(Fig.5)

图 5 洪泽湖水质理化因子与月际Chl.a浓度的CCA分析结果(Chl1~Chl12分别代表1-12月的Chl.a浓度) Fig.5 Results of canonical correspondence analysis between water quality parameters and monthly chlorophyll-a concentration in Lake Hongze

洪泽湖3个湖区的站点可以较好地根据象限区分开，东部湖区站点(洪泽湖区(淮安北)、洪泽湖区(淮安西)、洪泽湖区(淮安东)、洪泽湖区(淮安南))位于轴1右侧，与温度、DO分别呈正、负相关关系(图 6).通常情况下温度上升可以促进藻类生长繁殖.但温度是影响水体DO的关键因素，非夏季月份，温度的升高导致DO浓度降低^[29]，因此两者呈负相关.此外东区是过水通道，湖流运动相对剧烈，一定程度上影响DO浓度，因而东部湖区Chl.a浓度受温度和DO影响.北部湖区(颜圩)Chl.a浓度主要受营养盐浓度影响，受DO浓度等因子影响较弱，主要是由于北部湖区相对封闭，且周围养殖、农耕业相对发达，外源营养盐的排入对北区Chl.a浓度产生重大影响.西部湖区同时受TN和DO浓度的影响，西部湖区周边也进行了一定程度开发，人工投饵与施肥等均会对该区营养盐状况产生影响.同时，淮河从西区老子山附近入湖，湖流运动也会对该区藻类生长造成影响.

图 6(Fig.6)

图 6 洪泽湖水质理化因子与各站点Chl.a浓度的CCA分析结果(yw：颜圩；gh：高湖；lh：临淮；sqn：洪泽湖区(宿迁南)；hab：洪泽湖区(淮安北)；hax：洪泽湖区(淮安西)；had：洪泽湖区(淮安东)；han：洪泽湖区(淮安南)) Fig.6 Results of canonical correspondence analysis between water quality parameters and chlorophyll-a concentration in each water sampling site of Lake Hongze

选取洪泽湖2012年12月至2013年11月水质监测数据，通过方差分析和CCA统计分析归纳总结洪泽湖Chl.a浓度时空变化规律，并进一步探究其原因.从时间维度上看，洪泽湖Chl.a浓度季节变化规律在不同湖区有所差异.东部湖区Chl.a浓度随季节变化曲线呈“双峰型”，北、西部湖区Chl.a浓度在春季保持平缓，夏季有所回落后上升，在秋季达到峰值.从空间维度上看，3个湖区的Chl.a浓度变化在春、冬季均存在显著差异(P＜0.05)，其余季节差异不显著.全年月份和湖区站点Chl.a浓度与水质理化因子(COD、DO、TP、TN、温度)之间的CCA分析(图 6)显示，不同月份和不同湖区Chl.a浓度与水质理化因子之间存在不同的响应关系.从以上结论可以看出：

1) 洪泽湖是典型的过水性湖泊，东区包含过水通道，湖流运动较为剧烈，对藻类生长起到一定的抑制作用.

2) 洪泽湖北区营养盐水平和Chl.a浓度较高，多是由于上游城市汇水以及周围水产养殖、农耕施肥带来的过剩营养盐.丰富的营养盐为藻类大量繁殖提供物质基础，因而造成北部湖区较为严重的水质污染.

3) 在洪泽湖藻类防控治理方面，应该充分考虑上游污水的截留减排以及控制周围农业、渔业的投饵施肥量，以降低湖区营养盐浓度，有效控制藻类生物量.