调水前后南四湖藻密度及水环境要素变化特征*

谭 浪，王宗志，白 莹，舒博宁，王卫光

(1：南京水利科学研究院水文水资源与水利工程科学国家重点实验室，南京 210029) (2：河海大学水文水资源学院，南京 210098) (3：长江保护与绿色发展研究院，南京 210098) (4：山东省济宁市水文中心，济宁 272004)

南四湖是我国华北地区最大的淡水湖泊，是鲁南、苏北地区的重要水源地，也是南水北调东线工程的重要调蓄场所，其水生态稳定性对沿线水质安全及调水效益起决定性作用[1]. 南水北调东线工程在2013年12月正式开始运行后，增加了湖泊水量补给和物质输入，也改变了水动力条件、水文循环过程与水生动植物栖息环境[2]. 在2002年南水北调东线工程建设初期，由于围湖造田、过度捕捞、工业排放等活动影响，南四湖入湖河流和湖区污染严重[3]. 徐好等[4]在南四湖水质时空变化特征研究中，强调了“治”“用”“保”并举等流域综合治理措施，发现2006-2016年间南四湖水质有所改善，基本稳定在地表Ⅲ类水标准. 在浅水湖泊富营养化波动与高强度治理发展互动的态势下[5]，湖泊水生态环境的演变是亟需关注的重大问题.

受流域内外源污染、外调水、水文、气象多重影响，南四湖水环境水生态安全面临严峻挑战. 在调水影响方面，调水工程的建设及其常态化运行，对庞大而复杂的湖泊生态系统有何影响是目前研究中还未解决的难点. 对南水北调东线6个调蓄湖泊水质状况的研究表明[6]，从上游至下游湖泊水质呈改善趋势，调水期存在TN浓度较高，非调水期TP浓度较高；李庆等[7]对比分析调水期与非调水期，湖泊浮游植物与水质状况，发现调水引起的水文脉冲对浮游植物有明显影响；张春梅等[8]认为应优化调水方式，调控水动力条件，从而有效防控藻类突增. 在水文气象影响因素方面，气温升高、极端天气增多、日照时间延长等均有利于藻类生长[9]；邓文文等[10]讨论了气候变化和人类活动导致的生物扰动、水位波动增强等变化，会造成突变阈值减小，并对湖泊稳态造成影响. 藻密度作为湖泊富营养化水平的关键指标之一[11]，其长序列年内、年际变化趋势及水环境要素演变特征对识别调水的生态效应具有重要指示作用. 董贯仓等[12]在南水北调东线工程运行3年后对东平湖浮游植物及水质的调查，发现调水工程可能通过物质带入影响水质, 进而引起藻密度较通水初期稍有升高；2015-2019对南水北调中线干渠的研究表明[13]，藻密度每年都有一致性季节变化，在春季增长，随后从夏季到冬季持续下降，藻密度最大值的夏季和最小值的冬季分别为572.95×104和157.09×104cells/L，氨氮是影响藻类生长的最关键水质参数，其次是COD. 大规模调水背景下，浅水湖泊营养水平发生剧烈波动，藻类生长与相关水环境要素的互馈关系不确定性增加. 目前，少见周期性调水行为前后湖泊藻类生长趋势对比及其影响因素识别研究.

本文依据2010-2020年长系列监测资料，综合分析南水北调东线一期工程正式运行前后，南四湖藻密度及相关环境因子年内、年际变化特征，在此基础上识别南四湖水生态水环境的关键影响因素，以期为南四湖科学管理及南水北调东线工程良性运行提供科学依据.

1 材料与方法

1.1 研究区概况

南四湖位于山东西南部，南北狭长约125 km，东西宽5～25 km，湖内面积1266 km2. 湖泊中部建有二级坝枢纽工程，分成上级湖和下级湖，其中上级湖由南阳湖、独山湖和昭阳湖(大部分)组成，下级湖由昭阳湖(小部分)和微山湖组成. 南四湖入湖河流合计53条，河流同时用于行洪、农灌、纳污，主要入流河道宽浅，水流较缓，并且南四湖各入湖河流上修建了较多的水库及闸坝工程，起到拦蓄调洪作用，湖滨区沟渠串通严重，引进、引出水量不能全面控制，大部分河流无记录年径流量，无法精确控制入湖水量.

南水北调东线工程将长江水通过各级泵站提水向北输送，经南四湖等多个调蓄湖泊后，流至胶东地区和天津. 南水北调东线一期工程调水期在上半年，调水量逐年增加，工程拟利用南四湖进行输、蓄水. 南四湖地势北高南低，天然水流方向是从上级湖至下级湖，在调水工程中，上级湖建长沟泵站抽水入梁济运河，二级坝建泵站抽水入上级湖，下级湖建韩庄闸泵站，上级湖起到只输水不调蓄的功能，下级湖入流较少，作为调蓄湖泊. 《南水北调总体规划》对输水干线的水质要求达到地表水III类标准，由于输水干线穿过城市人口密集区、工农业区、航运区等，目前输水水质保障工作仍面临着挑战.

1.2 采样与监测

在南四湖设置10个采样点，空间分布见图1，包括3个上级湖点位(S1、S2、S3)、3个下级湖点位(S4、S5、S6)、4个距离工业区和生活区较近的主要入湖河流断面(S7、S8、S9、S10). 总氮(TN)、总磷(TP)、氨氮、高锰酸盐指数(CODMn)、五日生化需氧量(BOD5)采样点为S1～S10. 藻密度、叶绿素a(Chl.a)采样点为S1、S2、S3、S5. 水文气象数据由济宁市水文中心提供，水文气象站点位置与水质监测站点位置一致. 监测时段为2010年1月-2020年12月，频率为月度监测，每月上旬采样及送样检测.

图1 采样点分布Fig.1 Distribution of sampling sites

柱状采水器采集水面下0.5 m处水样4 L，2 L聚乙烯采样瓶低温避光保存，用于水质测定；2 L加入鲁哥试剂固定后在实验室内静置24～48 h，用虹吸管缓慢吸走上清液，保留50 mL浓缩样品. TN浓度采用硫酸钾氧化-紫外分光光度法测定；TP浓度采用钼锑抗分光光度法测定；氨氮浓度采用纳氏试剂光度法测定；CODMn采用酸性高锰酸钾法测定；BOD5采用稀释接种法测定；Chl.a浓度参照《水质叶绿素a的测定分光光度法》(SL 88-2012)进行测定；藻密度参照《内陆水域浮游植物监测技术规程》(SL 733-2016)移液样品至浮游植物计数框中，用40～600倍光学显微镜计数. 取上级湖各样点检测平均值代表上级湖指标；取下级湖各样点检测平均值代表下级湖指标.

1.3 分析方法

1.3.1 污染指数计算 综合污染指数是一种常用的水质评价方法，是利用实测数据和评价标准的比值来综合评价水质状态，评价标准采用《地表水环境质量标准》(GB 3838-2002) Ⅲ类水质标准，计算公式：

Pi=Ci/Si

(1)

(2)

式中,Pi为i污染物的污染指数；Ci为i污染物实测平均值，mg/L；Si为《地表水环境质量标准》(GB 3838-2002)中水质因子i的Ⅲ类水质标准值，mg/L；n是选取污染物的项数；P为综合污染指数.

1.3.2 时间序列分析方法 采用Mann-Kendall(M-K)检验法分析时间序列趋势[14]. 可以利用构造检验统计量Z，判断序列在置信水平上，时间序列数据是否存在明显的上升或下降趋势，当Z的绝对值大于等于1.28、1.64、2.32时，分别表示通过了信度90%、95%、99%的显著性检验. 也可以利用M-K检验对各项水质参数进行突变检验，建立统计序列UFk和UBk，进一步分析水质序列的变化，明确突变的区域，UFk为M-K统计值正向序列，UBk为逆向序列，(-1.96，1.96)为显著水平置信区间，若UFk超过显著水平置信区间(临界线)，表明上升或下降趋势显著，若UFk和UBk序列出现交点且交点在临界直线之间，那么交点即为突变点. 在Mann-Whitney-Pettitt(MWP)检验中[15]，确定显著水平α=0.05后，寻找具有显著性的变异点并可求得整个变异区间. 本文使用M-K检验法与MWP检验法找出逐月时间序列所有变化节点，并应用累积距平曲线法[16]对变化趋势和突变点进行识别和综合诊断，确定水质指标序列中最合理的突变年份. 时间序列自相关性对趋势检验结果有重要影响，所以在本次趋势检验之前，先进行了时间序列的自相关性分析，各序列具有良好的独立性. 采用Spearman相关系数分析南四湖与主要入湖河流之间水质指标的相关性、藻密度与相关环境因子的相关性. 采用独立样本t检验分析指标浓度不同时段的差异性，其中统计量P值为原假设为真时所得到的样本观察结果或更极端结果出现的概率. 数据统计及分析在Excel 2010、MATLAB、SPSS 20.0中完成.

2 结果

2.1 水文气象要素变化特征

2010-2020年南四湖湖区平均水温为15.7℃，夏季最高月平均水温为29.5℃，出现在2013年8月. 2010-2020年南四湖湖区年降水量均值为686 mm，流域内75%的降水和入湖径流集中在6-9月份的汛期，其中2013年与2014年降水量分别为597、610 mm，属于偏枯水年份(图2). 南四湖径流分布格局基本符合降水的分布规律[17]，而且入湖径流量与湖区水位涨幅基本对应，湖泊水位是反映水体水情变化较为直观的指标[18]. 由于2013年与2014年沂沭泗流域降水量较同期偏少，入湖径流量减少，湖区蓄水量持续减少，2014年6月12日下级湖水位降至死水位(33.0 m)以下，7月28日南四湖下级湖降至生态水位(31.05 m)以下，部分湖区裸露干裂. 在持续的干旱后有较强降雨，2016年南四湖湖区年降水量为794 mm，出现旱涝交替的现象.

图2 2010-2020年南四湖水文气象因素变化Fig.2 Changes in hydrometeorological factors in Lake Nansi from 2010 to 2020

2.2 藻密度及水环境要素年际变化特征

根据2010-2020年南四湖藻密度及其水环境要素的统计结果(表1)，南四湖藻密度变幅最大，其中上级湖藻密度变异系数为2.046，藻密度最大值是最小值的120倍. 对2010-2020年南四湖上、下级湖藻密度进行M-K趋势性检验，诊断结果中统计量Z分别为-2.02、-1.04，其中上级湖藻密度序列通过了信度95%的显著性检验，序列存在明显下降趋势. 但在M-K检验方法中，趋势性只受到样本秩与样本个数n的影响，从而忽略了样本中排序靠后的较大数据值[19]，从图3可见，南四湖藻密度具有显著的季节特征，上级湖藻密度在2016年后部分月份出现极大值，且远远高于其他月份，故不能得出上级湖藻密度在2010-2020年呈显著减少趋势. 其中上级湖2016-2020年10月藻密度分别为520.0×104、515.0×104、500.0×104、480.0×104、470.0×104cells/L.

表1 2010-2020年南四湖藻密度及其水环境要素统计量*

图3 2010-2020年南四湖藻密度及水环境要素变化(灰色区域表示上、下级湖UFk同时超越临界线的范围)Fig.3 Changes of algae density and water environment factors in Lake Nansi from 2010 to 2020(The gray area indicates the range where the UFk exceeds the critical line)

在突变点分析过程中，全部样本数据序列两端的部分样本点就不能做突变点的判断，故靠近2010年与2020年两端的突变点不予采信. 结合M-K检验、MWP检验、累积距平曲线法(图略)分析各站点水质序列的突变点，其中在M-K检验中，若UFk超越临界直线，表示上升或下降趋势显著. 藻密度及水环境要素序列在2010-2020年间存在多个突变点(表2)，局部年份变化明显，突变点主要集中在2013-2015年. 南四湖上级湖藻密度在2013年12月出现突变点，下级湖在2014年4月出现突变点，均在突变发生后开始增长至2016年，但不显著；上、下级湖TP、氨氮浓度分别在2014年、2013年有突变点，且均在突变点后一直保持显著上升趋势；Chl.a浓度在2015年存在突变点，随后浓度有所上升；TN、氮磷质量比(TN/TP)、BOD5、CODMn在上、下级湖中没有识别出突变点.

表2 南四湖藻密度及其水环境要素序列趋势突变点*

南四湖水质总体良好(图4)，2010-2020年南四湖上、下级湖水质综合污染指数年均值范围分别为0.759～0.945、0.719～0.926，基本达到了东线工程输水水质地表Ⅲ类水要求，但仍有个别污染物的污染指数出现较大值，造成总体水质不能稳定达标，其中上、下级湖水质综合污染指数最大值均出现在2013年. 上、下级湖BOD5污染指数部分年份超过1，TN污染指数一直大于1，TN单项指标污染情况最严重. 对2010-2020年南四湖上、下级湖水质综合污染指数进行M-K趋势性检验，诊断结果中统计量Z分别为-4.11、-3.36，通过了信度99%的显著性检验，序列存在明显下降趋势. 2020年南四湖上、下级湖水质综合污染指数年均值分别为0.796、0.767，2010年分别为0.886、0.907，与2010年相比分别下降10.2%、15.4%，湖泊水环境质量总体出现一定程度的好转.

图4 2010-2020年南四湖水质综合污染指数变化Fig.4 Changes of comprehensive pollution index of water quality in Lake Nansi from 2010 to 2020

2.3 藻密度及水环境要素年内变化特征

根据南水北调东线一期工程的运行规则和南四湖流域的水文气象特征，将南四湖藻密度及水环境要素从年际划分为南水北调东线一期工程运行前后两个部分(2010年1月-2013年12月、2014年1月-2020年12月)，从年内划分为汛前、汛中、汛后(2-5月、6-9月、10-1月)，从而明晰南四湖湖区2010-2020年藻密度及水环境要素年内分布情况(图5).

图5 调水前后南四湖藻密度及其水环境要素年内分布特征对比(图中横括弧显示统计量P值，**表示在0.01级别呈显著差异，*表示在0.05级别呈显著差异)Fig.5 Comparison of the distribution characteristics of algae density and water environment factors in Lake Nansi before and after water transfer

结果表明，调水工程使湖区藻密度年内分布发生了明显变化，汛前藻密度平均值由33.96×104cells/L(工程运行前)下降至6.42×104cells/L(工程运行后)，并具有显著性差异. 调水工程引起的汛中、汛后藻密度平均值变化并没有通过差异显著性检验，但南四湖藻密度每年峰值主要集中在汛中(工程运行前)，转为出现在汛后(工程运行后). TN、TP的年内各时期浓度在调水工程前后变化规律并不相同，但从氮磷比的角度分析，调水对南四湖湖区氮磷比造成了较大影响，氮磷比汛前、汛中、汛后在工程运行后相对工程运行前显著下降，分别降低了30.7%、45.3%、39.5%，没有改变氮磷比汛前高、汛后低的分布情况. BOD5年内各时期在调水工程前后也显著下降，分别降低了32.1%、42.5%、29.7%. Chl.a、氨氮、CODMn在工程运行前后，年内分布没有明显变化，均在汛中出现较高值.

3 讨论

3.1 南四湖营养元素变化原因分析

对比2006-2010年[20]，2010年以来南四湖整体水质明显改善，但较高浓度的氮使东线一期工程水质保护问题凸显. 林雪原等[21]指出流域TN排放中，农田化肥和畜禽养殖贡献率达到85.39%，而农业面源污染型治理难度大，短期内改善程度较低；二级坝的建立，使南四湖分为上、下级湖，湖中水体流速减缓、污染物滞留时间增加，使水质进一步恶化. 南四湖湖面狭长，受纳苏、鲁、豫、皖4省32县(市)来水， TN污染发生随机、机理复杂、时空差异性大、来源多和来量大让南四湖流域的TN污染治理问题难以突破.

水文气象通过影响水体中污染源的来源、迁移、生态效益等过程直接或间接对水生态环境造成影响. 南四湖在2013-2015年水文气象条件有所变化(图2)，在人为干预下，2014年与2015年二级湖闸站全年未开闸放水，下级湖总出湖年径流量为“0”. 上、下级湖与运河之间阻隔，减少湖泊的连通性与流动性，大大增加了流域入湖污染物在水体滞留时间，水体自净能力和纳污能力都会减弱，加之高温与水浅，水质持续恶化，生态环境受到严重威胁[22]. 为应对严峻的旱情，紧急制定了2014年南四湖生态应急调水方案，通过南水北调东线从长江向南四湖实施调水. 调水工程引起的外源输入变化，同时南四湖的严重旱情、水体隔绝、营养富集，可能是2013年与2015年期间南四湖水环境指标的时序发生突变和突变点后上升趋势的原因.

3.1.1 入湖河流的影响 入湖河流的水质状况是南四湖水质状况优劣的决定性因素[23]. 分析2010-2020年南四湖上、下级湖水质与对应入湖河流水质相关程度(表3)，南四湖主要入湖河流水质对南四湖整体水质影响显著，其中上级湖与入湖河流的TN、BOD5呈显著相关，下级湖与入湖河流的各水质指标都呈显著相关. 并结合图3和表4，上、下级湖TN处于较高水平，这与入湖河流长期的高浓度氮输入有一定关系. 解文静等[24]也认为南四湖水体接纳入湖河流的污染，河流中有机物污染相对严重，含氮量高.

表3 2010-2020年南四湖与主要入湖河流水质指标Spearman相关系数

随着南四湖流域城市建设用地和农村建设用地整体分布从点到面的扩大，下垫面条件的改变[25]，导致不透水面积增加和农业面源污染加重. 将含暴雨等级以上降雨强度的月份从总体时序中筛选出来，对比降雨时期河流污染程度与多年平均值，见表4. 南四湖流域存在强降水时，河流中各营养元素都存在不同程度的增加，一方面可能是南四湖作为流域内污染源排放的唯一归宿，在受汛期暴雨冲刷的影响，大量的营养元素随排水或雨水淋溶进入水体，湖泊水体污染负荷势必会增加，导致汛中整体水质较差，另一方面可能是河流中径流量与流速增大，使河道中氮磷污染物从底泥中释放再悬浮，进一步增大河流中污染负荷.

表4 主要河流降雨时期水质与多年均值对比

3.1.2 调水工程的影响 调水是通过改变水文条件和营养盐浓度间接影响藻类的生长与消落，而且湖泊本身存在点面源污染，湖区水生态环境不会因为调水立刻发生变化，增加了影响因子量化的困难性，导致难以具体分析受水区水环境与调水工程响应关系，但从水环境指标的长序列变化来看，上半年的规律性调水营造出了年内与年际的周期性水文脉冲，影响了水环境中各要素或各要素之间的转化和传输，改变了湖泊生态系统的自然规律[26].

对比调水工程建设运行前后(图5)，汛前调水期氮和磷浓度有所增加，可能是调水泄流产生的水动力扰动会让底泥中大量悬浮物进入水体，而底泥作为污染物的天然蓄积池，沉积物营养盐含量远远高于背景值[27]，释放的加快造成水体二次污染. 从外源污染输入方面分析，南四湖调水期水质的改善可能是经过骆马湖与京杭运河的来水量大且水质较好. 从内源污染方面分析，由于南四湖窄长的湖泊形态和泵站抽排水，湖泊水位波动较大(图2)，在自然变化和人为调水的干预下，水位是受到影响最显著的水文过程[28]，有研究表明[29-31]，水位波动强烈时，会导致湖泊的自净能力和环境承载能力降低. 湖泊过高水位对入湖河道的顶托和水漫坡地面积增加，导致地面径流和地质带来的营养盐增多，导致湖泊富营养化；湖泊过低水位会对污染输入响应更剧烈，水环境指标易出现峰值，增加突发环境污染事件的风险.

3.2 藻密度的影响因素

通过南四湖藻密度与各指标Spearman相关性分析结果表明，水温、降水对南四湖藻密度的变化起着至关重要的作用(表5). 由上述分析可知，降雨导致大量外源污染物进入南四湖水体，通过改变藻类所需营养盐浓度而间接影响藻类生长. 南四湖藻密度四季变化明显，有大量研究表明[32-34]，当温度较低时，水体藻类生长受到限制，藻在春季生物量并不高，夏季日照增强、水温较高，随着温度的升高，绿藻、蓝藻等喜温浮游植物大量生长. 温度升高与藻类生长消落过程会增强水体中氮的吸收利用及脱氮能力[35]，这可能是造成藻的生物量与氮呈负相关，湖体氮浓度峰值一般出现在初春的原因(图5).

表5 南四湖藻密度与主要水环境要素Spearman相关系数

对不同营养盐的相关关系对藻密度的影响程度并不一致，其中与氮磷比呈显著负相关，与TP呈显著正相关. 氮和磷作为藻类生长的必需营养元素，两者之比对藻密度的变化有着重要的影响[36]，Redfield等[37]认为浮游植物生长最适氮磷比为16∶1，大于此值时为磷素限制状态，也有研究发现[38]最适氮磷比根据不同水体而发生变化. 南四湖氮磷比远远高于最适氮磷比，结合图3、图5分析，氮磷比在调水工程建设运行后明显下降，而且藻密度在调水后的增长趋势与TP浓度上升基本上保持一致，表明南四湖可能为磷限制性湖泊. 藻密度细胞与Chl.a浓度都为表征藻类生物量的2种重要方法，两者之间也具有显著的相关性. 南四湖藻密度与CODMn呈显著正相关，与BOD5相关性较低，与田时弥[39]、Qu等[6]的研究结果相似，水中有机物经过氧化分解后也是藻类可以利用营养盐的直接来源.

4 结论

1)2010-2020年南四湖上、下级湖水质综合污染指数年均值范围分别为0.759～0.945、0.719～0.926，2020年南四湖上、下级湖水质综合污染指数年均值分别为0.796、0.767，2010年为0.886、0.907，与2010年相比分别下降10.2%、15.4%，南四湖整体水质逐渐改善，但TN单项污染指数仍超标，与入湖河流长期的高浓度TN输入有关. 南四湖水污染防范与治理可将重点落脚于入湖河流关键性水质指标污染的治理.

2)同时受南水北调东线工程运行初期与严峻旱情的叠加效应影响，南四湖藻密度、TP、氨氮、Chl.a监测序列在2013-2015年有突变点，且在突变点后指标浓度均呈增加趋势.

3)南水北调东线工程的运行，改变了南四湖藻密度年内分布情况，改善了调水期总体水质情况. 藻密度年内波动较为强烈，在调水期结束后的10月份出现峰值，其中2016年上级湖藻密度高达520.0×104cells/L.

4)南四湖藻密度的变化不仅仅与TP、氮磷比、CODMn呈显著相关，还与水温和降水关系密切.