湖泊科学   2023, Vol. 35 Issue (4): 1359-1369.  DOI: 10.18307/2023.0427
0

研究论文——生物地球化学与水环境保护

引用本文 [复制中英文]

任鹏, 范中亚, 杨忠勇, 朱士江, 王文才, 崔玉洁, 和玉芳, 沿江城市中小型闸控通江湖泊总磷变化特征及成因分析——以安庆市石塘湖为例. 湖泊科学, 2023, 35(4): 1359-1369. DOI: 10.18307/2023.0427
[复制中文]
Ren Peng, Fan Zhongya, Yang Zhongyong, Zhu Shijiang, Wang Wencai, Cui Yujie, He Yufang. Variation characteristics and causes of total phosphorus in small and medium-sized gate-controlled lakes along the Yangtze River—A case study of Lake Shitang in Anqing City. Journal of Lake Sciences, 2023, 35(4): 1359-1369. DOI: 10.18307/2023.0427
[复制英文]

基金项目

国家自然科学基金项目(52009066, U2040220)、长江生态环境保护修复联合研究项目(第一期)(2019-LHYJ-01-0212, 2019-LHYJ-01-0212-17)、长江生态保护与修复安全驻点研究项目(2018CJA030301-014)、长江科学院开放研究基金项目(CKWV20221003/KY)和湖北省高校优秀中青年科技创新团队项目(T2021003)联合资助

通信作者

朱士江, E-mail: zhusjiang@gmail.com
王文才, E-mail: wangwencai@scies.org

文章历史

2022-01-16 收稿
2022-09-30 收修改稿

码上扫一扫

沿江城市中小型闸控通江湖泊总磷变化特征及成因分析——以安庆市石塘湖为例
任鹏1,3 , 范中亚2 , 杨忠勇1 , 朱士江1 , 王文才2 , 崔玉洁1 , 和玉芳4     
(1: 三峡大学,三峡库区生态环境教育部工程研究中心,宜昌 443002)
(2: 生态环境部华南环境科学研究所,广州 510655)
(3: 浙江华东工程咨询有限公司,杭州 311122)
(4: 河北省水利规划设计研究院有限公司,石家庄 050021)
摘要:长江中下游沿江城郊闸控湖泊普遍面临总磷浓度偏高的现象,解析其总磷时空变化特征及影响驱动机制成为精准治理与修复此类湖泊前亟待解决的关键性问题。本文基于安庆市沿江城郊中小型闸控湖泊——石塘湖实测气象降雨、水文、河湖水质等数据,采用多因子相关性分析、变异系数法和主成分分析的方法,研究各指标因子与湖泊总磷浓度时空变化的响应关系。结果表明:(1)湖泊水质呈现丰(5-8月)、平(3-4月和9-10月)、枯(11-次年2月)水期聚类效果显著,但湖泊空间差异不明显;湖泊水质的季节性变化受总磷浓度变化控制,其他理化指标影响较弱。(2)高强度降雨和汛期闸站调度下的水动力变化决定污染物迁移速率,是导致湖泊总磷浓度在丰水期更容易受入湖河流输入影响而达到峰值的主要驱动因子。(3)高强度人类活动导致入湖河流季节性输入是石塘湖总磷上升的决定性因素,从单位土地利用类型产生单位总磷负荷来看,农业用地远大于城镇建设用地,林地和草地充当污染物进入湖泊的预前“汇”。因此,从治理策略和途径来看,可在控制外源输入的同时,适当调节闸站以减弱丰水期水动力强度,减少入湖河流总磷输入影响的同时,发挥湖泊营养盐的滞留净化能力。本文主要从人类活动导致的外源污染及闸站抽排调度角度分析沿江城郊闸控湖泊石塘湖总磷变化、影响因素及驱动机制,可以为此类城郊中小型沿江闸控型湖泊污染治理提供理论参考。
关键词总磷    闸站调度    水动力    驱动机制    石塘湖    长江    
Variation characteristics and causes of total phosphorus in small and medium-sized gate-controlled lakes along the Yangtze River—A case study of Lake Shitang in Anqing City
Ren Peng1,3 , Fan Zhongya2 , Yang Zhongyong1 , Zhu Shijiang1 , Wang Wencai2 , Cui Yujie1 , He Yufang4     
(1: Engineering Research Center of Eco-Environment in Three Gorges Reservoir Area, Ministry of Education, China Three Gorges University, Yichang 443002, P.R.China)
(2: South China Institute of Environmental Sciences, Ministry of Ecology and Environment, Guangzhou 510655, P.R.China)
(3: Power China Zhejiang Huadong Engineering Consulting Corporation Limited, Hangzhou 311122, P.R.China)
(4: Hebei Water Resources Planning and Design Institute Co. LTD, Shijiazhuang 050021, P.R.China)
Abstract: The gate-controlled lakes along the middle and lower reaches of theYangze River are generally facing with the phenomenon of high total phosphorus concentration. Analyzing the spatiotemporal variation characteristics and impact driving mechanisms of total phosphorus has become a key issue that needs to be urgently addressed before precise governance and restoration of such lakes. Based on the measured meteorological rainfall, hydrology, water quality and other data of Lake Shitang, a small and medium-sized gate-controlled lake in the suburbs of Anqing City along the Yangtze River, this paper uses the methods of multi factor correlation analysis, coefficient of variation and principal component analysis to study the response relationship between each index factor and the temporal and spatial changes of total phosphorus concentration in the lake. The results show that: (1) the clustering effect of lake water quality is significant during the abundant (May-August), flat (March-April and September-October), and dry (November-February of the following year), but the spatial differences of lakes are not significant. The seasonal changes in lake water quality are controlled by changes in total phosphorus concentration, while other physical and chemical indicators have a weaker impact. (2) The hydrodynamic changes under high intensity rainfall and flood season gate station scheduling determine the migration rate of pollutants, which is the main driving factor that makes the total phosphorus concentration in lakes more susceptible to the impact of river input during the flood season and reaches its peak. (3) The seasonal input of rivers entering the lake caused by high-intensity human activities is the decisive factor for the increase of total phosphorus in Lake Shitang. From the perspective of the unit total phosphorus load generated by unit land use type, agricultural land is far greater than urban construction land, and forest and grassland serve as the pre sink of pollutants entering the lake. Therefore, from the perspective of governance strategies and approaches, it is possible to control external inputs while appropriately adjusting gate stations to reduce the hydrodynamic intensity during the high water season, reduce the impact of total phosphorus input into rivers entering the lake, and at the same time, exert the retention and purification capacity of lake nutrients. This article mainly analyzes the changes, influencing factors, and driving mechanisms of total phosphorus in Lake Shitang, a gate controlled lake in the suburbs of the Yangtze River, from the perspective of external pollution caused by human activities and gate station discharge scheduling. It can provide theoretical reference for the pollution control of small and medium-sized gate controlled lakes in such suburbs.
Keywords: Total phosphorus    regulation station    hydrodynamic    driving mechanism    Lake Shitang    Yangtze River    

我国地形和气候条件复杂多样,湖泊广泛分布,且以中小型湖泊为主,尤其是在长江流域以中小型湖泊居多。根据第二次全国湖泊调查结果显示,我国水面面积大于1 km2的湖泊有2693个,其中中小型湖泊数为2666个,占全国湖泊总数的99%,沿长江流域中小型湖泊数量接近全国湖泊总数的1/4,占流域湖泊总数的99.38%[1]。自1950s以来,长江中下游通江浅水湖泊兴建闸坝、筑堤等水利工程,极大地改变了江湖原有的自然连通状态,导致绝大多数湖泊成为人为阻隔湖泊,目前长江流域建设的水闸数量已超过5×104个,仅剩洞庭湖、鄱阳湖、石臼湖等少数湖泊为自然通江状态[2-3]。闸控性湖泊物质循环及迁移输运过程受人为调度影响极为显著,特别是近年来,随着城市化快速发展和农业密集生产,产生大量污染汇入湖泊水体,城市和城郊湖泊水质恶化、水体富营养化等问题突出[4]

湖、库等水体富营养化过程主要是水体中汇入过量的氮、磷等营养元素,从而导致在竞争中占优势的藻类大量繁殖,致使湖泊水生态系统中初级生产力的异常暴发[5-9]。根据加拿大湖沼研究中心沃伦威德博士研究表明,80%的湖泊富营养化受磷元素制约,10%的湖泊富营养化与氮、磷元素直接相关,其余10%的湖泊是氮与其他元素起作用[10],表明湖泊富营养化问题更容易受到磷元素的影响,而且近5年《中国生态环境状况公报》显示,长江流域湖泊水质改善主要受磷浓度限制,但是现有的湖泊研究多集中在区域性重点大型湖泊,如太湖[11],或水库,如三峡水库[12]、密云水库[13]等,或自然性小型养殖性池塘[14],对于与绝大部分经济和人类生活息息相关的沿江城郊闸控型湖泊富营养化研究相对较少,尤其是富营养化控制因子总磷与闸站调度水动力的相关关系还有待解析。因此,本文以安庆市沿江城市中小型闸控通江湖泊石塘湖为研究区域,分析在闸控调度下石塘湖流域土地利用类型、入湖河流输入、气象气候条件及闸站调度过程对石塘湖总磷时空分布影响及成因解析,结合控水闸门的调度方式,进一步探究石塘湖总磷控制及改善途径,为研究不同水位波动下河湖物质迁移输运过程与污染治理提供参考借鉴,也可以为此类河湖水系连通战略实施,提高水资源统筹配置能力,提高生态环境质量和加强防洪保障能力提供理论基础。

1 研究区域概况

石塘湖位于安庆市城区东北部,是安庆市的备用水源地,水体与破罡湖相连,水流自西南流向东北,并经破罡湖闸站进入长江,属于典型的沿江城郊闸控型湖泊,水动力过程受破罡闸闸站调度影响。石塘湖湖泊面积约14.0 km2,多年平均水深1.5 m,最大水深3.93 m,所在区域属亚热带湿润季风气候,多年平均气温16.5℃,历史上极端最高气温为1934年7月13日的44.7℃,极端最低温度为-12.5℃,最热月(7月)平均气温为28.8℃,最冷月(1月)平均气温为3.5℃,全年主导风向为东北风,夏季为西南风,年平均风速3.1 m/s,最大风速20 m/s[15]。安庆市年内降雨分布极不均匀,春、冬季降雨量较小,夏、秋季进入梅雨季节,尤其是在丰水期(5-8月)降水最为集中,且强度较大,冬季(12月次年2月)是一年内降水最少的季节,12月和1月份降水量均不足50 mm,年均降雨量1389.1 mm,多年平均蒸发量1611.4 mm,多年平均降雨时间139.1 d,多年平均降雪时间12.8 d[15]。其中,4-8月降雨量占全年的65.7%[16],在丰水期极易受洪涝灾害的影响。石塘湖主要污染特征为水体磷浓度季节性偏高,且随着社会发展和历史变迁,石塘湖TP浓度呈现出阶段性变化,1981年以前水体总磷浓度为50~60 μg/L,湖泊一直维持在中等营养水平;1981年以后水体总磷浓度呈现明显升高趋势(>100 μg/L)[16];2017年总磷浓度更是远超地表水Ⅴ类水标准[17],其余指标浓度相对较低,基本优于Ⅲ类水;石塘湖总磷浓度在年内波动较大,丰水期远超Ⅲ类限值,枯水期总磷浓度处于波谷,基本优于Ⅲ类水。石塘湖流域如图 1所示,由于石塘湖属于小型湖泊,监测点位空间上呈全湖平均分布,S1~S4代表石塘湖监测点,R1~R4分别代表跃进大沟、606大沟、界河、余敦排洪河4条主要入湖河流监测点。

图 1 研究区域 Fig.1 Study area
2 材料与方法 2.1 数据采集

采集指标包括水质理化指标:水温(WT)、溶解氧(DO)、电导率(σ)、透明度(SD)、pH值;富营养化指标:总磷(TP)、总氮(TN)、氨氮(NH3-N);有机污染物:高锰酸盐指数(CODMn)、化学需氧量(COD);水文气象数据:破罡闸站降雨、水位、闸站调度数据。气象降雨(主要是降雨量)、水文和闸站调度数据均来源于安庆市水旱情网(安庆市水旱情网:http://aqswj.cn:9000/aqshqw/shqmh/index.html#/home),时间范围均为2017-2019年,分辨率均为逐日;闸站调度主要关注闸站开关时间,可以根据人工记录和水位变化反推数据。湖泊水质指标来源于2017-2019年每月地表水常规监测;入湖河流水质为2018年逐月数据。水温、pH值、电导率和溶解氧采用便携式水质分析仪HQ40d(美国)现场测定,氨氮浓度采用纳式试剂光度法(HJ 535-2009)测定;TP浓度采用钼酸铵分光光度法(GB 11893-1989)测定;TN浓度采用碱性过硫酸钾消解紫外分光光度法(HJ 636-2012)测定;化学需氧量采用快速消解分光光度法(HJ/T 399-2007)测定; 高锰酸盐指数采用酸性法(GB/T 11892-1989)测定;石塘湖周边土地利用数据来源于2017年卫星遥感反演数据。

2.2 研究方法

运用主成分分析的方法提取石塘湖流域水质变化的主导因素,并通过多元统计评价方法[18]和相关性分析来探究石塘湖总磷负荷来源及成因解析,再利用变异系数法判别石塘湖总磷时空分布差异。在进行主成分分析之前,先对原始数据进行KMO和Bartlett球度检验,当KMO>0.6且P < 0.05时可以做主成分分析,为消除各指标之间因不同量纲的影响,将原始数据进行Z-Score标准化处理[19](标准化数据=(原始数据-均值)/标准差),根据特征值λ>1、累计方差贡献率接近甚至高于80%的原则,选取主成分方法分析水质数据时空分布特征,指标因子相关性及荷载;采用相关性研究湖泊总磷与主要入湖河流总磷相关性,解析其来源;变异系数法对石塘湖4个位点时间尺度和每月空间尺度进行柱状图分析,解析总磷因子在时间和空间离散程度。主成分分析主要采用Canoco5分析和作图;相关性分析采用SPSS 26分析和梳理,统计分析则采用Origin 2018绘制完成。

3 结果与讨论 3.1 石塘湖水质主要影响因子

根据统计检验结果显示KMO值为0.602,P值小于0.0001,水质理化指标之间存在较强相关性。按照特征值均大于1的原则,提取了影响石塘湖水质变化的4个主成分(表 1),其累计方差贡献率为85.472%,大于80%的要求[20],基本能充分反映原始数据信息。结果显示,第一主成分方差贡献率为46.296%,关联的主要参数有TP、WT、DO、σ、COD和CODMn,主要反映水体营养盐磷素(TP),有机污染物(CODMn、COD)和水质理化指标(WT、DO、σ)主导着水质变化;第二、三主成分的方差贡献率分别为14.866%和13.386%,关联的主要参数有NH3-N和TN,表示营养盐氮素对石塘湖水质变化的影响,第四主成分的方差贡献率为10.924%,主要表现为pH、SD等理化指标影响着水质变化。

表 1 石塘湖水质理化指标主成分分析结果* Tab. 1 Principal component analysis results of physical and chemical indicators of water quality in Lake Shitang

图 2是PCA分析双序(荷载系数和得分值)图,用采样样品得分和环境因子的荷载系数代替原始数据对水质进行分析[21],由图 2可以看出,采样数据没有呈现点位上的趋势性或聚集性分布,但在月份上呈现较为明显聚集性分布规律。3-4月主要分布于第一象限,5-6月和9-10月处于第二象限,7-8月采样点基本聚集在第三象限,11-2月则位于第四象限,丰、平、枯水期聚类效应显著。图中环境因子的箭头越长表示其因子荷载越大;箭头间的夹角越小,表示相关性越高。图中箭头较长的环境因子分为TP、COD、CODMn,σ、WT和TN、SD、NH3-N三簇排列,主导着石塘湖水质变化,且每簇环境因子夹角较小,相关性较高。现将样点向环境因子箭头方向投影,投影位置在箭头方向越远,表示数值越大,该环境因子对其正贡献越大,反之越小。由图 2可知,TP、COD主导着7-8月水质变化,其中在7月TP浓度最高;9-10月CODMn、σ和WT值最高,对水质有较大影响;TN、SD、NH3-N影响着3-6月石塘湖水质。

图 2 石塘湖水质因子PCA分析双序(荷载系数和得分值)图 (圆圈代表监测位点和月份,如S1-1代表石塘湖S1采样点-1月份;箭头代表环境因子) Fig.2 Double-sequence (load factor and score value) diagram of PCA analysis of water quality factors in Lake Shitang (The circle represents the monitoring site and month, S1-1 represents S1 site in Lake Shitang-January, the arrow represents the environmental factor)
3.2 石塘湖水质理化指标对闸控湖泊藻类的影响特征与总磷的响应关系

湖泊水质恶化过程中,总磷等富营养化控制因子常常成为影响水质变化的主控因子,有机污染物和水质理化指标也对水质正常波动产生干扰。姚焕玫、尹炜等对武汉市城市湖泊东湖[22]和丹江口[23]的深入研究结果表明,TP等是主要污染指标,占据第一、二主成分;水质理化指标和有机污染分别占据第三和四主成分,影响着湖泊水质时空变化。石塘湖有机污染物COD、CODMn浓度基本处于目标水质Ⅲ类范围,对石塘湖水质影响有限;氮素年内呈小范围变化,整体优于Ⅲ类水范畴;与之相比总磷浓度年内季节波动变化较大,丰水期浓度高,枯水期浓度偏低,主导着石塘湖水质时空波动。藻类等水生植物的生长需要充足的营养盐和适宜的环境条件,分析水生植物群落结构和分布及其与总磷的关系,可揭示其对总磷浓度的影响[24-25]。根据2018年夏季石塘湖水生植被现状调查结果显示,石塘湖湖面已基本没有水生植被;北部靠近出口的湖区有少量野菱群丛,周边沿岸有少量芦苇、菰和芡实,植被覆盖度极低;从石塘湖藻类叶绿素a监测结果结果看,叶绿素a浓度呈现平水期(3-4月和9-10月)>丰水期(5-8月)>枯水期(11-次年2月),平水期相比丰水期闸站频繁调度下水动力条件更弱,与之枯水期来说水温和营养盐浓度更高,促进了浮游藻类生长,叶绿素a浓度达到年最高值;丰水期则受水动力条件的制约,叶绿素a浓度次之;枯水期由于低温和营养盐浓度限制,藻类小范围零星分布,叶绿素a浓度处于年最低值,其年变化范围在20~29 μg/L之间,呈季节性小范围波动状态,反映在石塘湖藻类生物量来看,年内变化相对稳定,对石塘湖营养盐尤其是总磷变化影响相对较小。

3.3 石塘湖总磷浓度与水文气象和闸站调度之间的关系

图 3为石塘湖水深、总磷浓度和降雨量变化过程图,水深变化反应石塘湖水位波动特征,受破罡闸调度管控,年内水位呈现丰、平、枯水期规律波动变化,丰水期会因安庆市防洪需要,破罡闸闸门开启,石塘湖水位迅速降至汛限水位以下,处于年内低水位;枯水期水位一般保持在正常蓄水位偏上;平水期则处于丰、枯水期之间变化。石塘湖降雨量在2018年1-3月比较平均且强度较小,水深基本无较大波动,总磷浓度基本稳定在0.05 mg/L,为地表水Ⅲ类水,随着降雨量快速增长到5月峰值94.5 mm,水深随即攀升至第一个波峰2.53 m,在这期间破罡闸关闭处于蓄水状态,水动力变化不大,与之同时年初来流水质较好,入湖含磷污染物主要在重力的作用下沉积,但入湖总磷浓度仍大于石塘湖初始浓度,石塘湖水位对入湖总磷浓度起到一定的稀释作用[26],总磷浓度表现缓慢增长,此时石塘湖表现出湖泊特性,接纳来自外源输入性磷;为满足安庆市丰水期防洪要求,破罡闸闸门开启,石塘湖水深急剧下降到年内最低值1.87 m,水动力达到峰值,驱动含磷颗粒物快速输运,伴随着高强度的降雨导致大量面源冲刷,外源高浓度的含磷微团快速在石塘湖汇集,石塘湖总磷浓度急剧攀升,进而达到峰值0.36 mg/L,远超作为安庆市备用水源地Ⅲ类水标准(湖泊Ⅲ类水标准为0.05 mg/L),此时湖泊主要表现出河流特性,受入湖河流总磷输入的影响,水质较好水体被排出石塘湖;丰水期结束,破罡闸关闭,但来流还未降到最低值,水流持续向石塘湖汇集,水深快速回升到第二个峰值2.57 m,入湖泥沙易于落淤,导致水体内以吸附在悬浮颗粒物(SS)上的TP下降[27],8—12月降雨量普遍偏小,且较为分散,外源汇集效率降低,且在湖泊自净作用下浓度持续下降。我们从石塘湖各时期换水周期也可以看出,相比平水期和枯水期来说,石塘湖在丰水期换水周期明显缩短,2017和2019年丰水期换水周期不到平、枯水期换水周期的一半,2018年丰水期换水周期相当于同年平、枯水期的58.18%和54.08%,即丰水期石塘湖更偏向于河流特征,换水周期较短。

图 3 石塘湖水深、总磷浓度和降雨量变化 Fig.3 Changes of water depth, total phosphorus and rainfall in Lake Shitang

石塘湖总磷浓度变化特征、水文气象条件、闸站调度三者之间存在极大关联性,降雨量的季节性变化驱动大量外源磷在石塘湖内汇集,为石塘湖总磷浓度升高提供必备条件,而闸站调度所引起的水位变动(水动力变化)为石塘湖内总磷迁移和浓度在短时间内快速上升提供驱动力。丰水期降雨较集中,且强度较大,外源磷向石塘湖内汇集速率加快;闸站调度(致破罡闸闸门开启)使水位快速下降到防洪水位,湖泊水动力加强,驱动高浓度含磷污染物在石塘湖内迁移输运,这期间石塘湖内总磷浓度剧增;平水期和枯水期降雨强度较小,破罡闸闸门关闭,水动力条件减弱,主要以沉积作用为主,伴随着稀释作用,致使总磷浓度下降。

3.3.1 闸控湖泊通常兼具湖泊和河流两种特性

石塘湖总磷浓度年内季节性变化过程是由水文气象条件与闸站调度的结果,也与石塘湖湖泊特性和湖泊面积息息相关。与自然性湖泊相比,闸控小型湖泊通常兼具湖泊和河流两种形态特征,不同水期下的总磷变化特征也存在着较大区别,在平水期和枯水期主要表现出阻隔湖泊特征,对入湖污染物进行接纳和净化;而丰水期开闸后更接近于河流特征,石塘湖总磷浓度变化与入湖河流总磷输入同步且差异较小,入湖河流在丰水期带来大量营养盐,对湖区水质产生较大的影响[28],且人为控制闸站湖泊其换水周期也远高于自然连通湖泊[29],但是在季节性换水周期变化上,闸控湖泊在丰水期闸站调度使换水周期大幅度减短,平水期和枯水期由于闸门关闭,换水周期普遍增加。

3.3.2 中小型闸控湖泊丰水期水体磷滞留作用较弱

石塘湖属于中小型湖泊,对于大型闸控湖泊而言,水动力的强弱会影响营养盐在湖中的迁移输运距离,根据朱伟等[30]的研究,高强度降雨带来的大量含磷微团在向出湖口输运过程中,磷会逐步滞留。湖泊水体在丰水期较大的水动力作用下,加大总磷的输运距离及速率,丰水期结束,水动力减弱,相应总磷高浓度等值线出现回缩,因此大型闸控湖泊水体磷影响主要集中在入湖区域,例如太湖富营养化区域主要集中在主要入湖河流梅梁湾、竺山湾乃至整个西部湖区[31-32],东部水域无论是富营养水平还是污染物浓度偏低,中小型闸控湖泊受外源输入及湖体水动力影响较大,全湖总磷浓度空间差异较小,尤其是丰水期闸站调度使整个湖区流域水动力条件处于较大值,高浓度含磷污染物在全湖迁移输运,水体磷的滞留效应与大型湖泊相比更加微弱。

3.3.3 防洪调度对城区湖泊总磷浓度产生影响

非防洪调度形湖泊,受枯水期入湖河流径流减少等影响,其湖泊总磷浓度通常呈现丰水期低,枯水期浓度较高的特点,此外,大型浅水湖泊还会因为枯水期由于低水位致使沉积物再悬浮及内源磷动态释放而总磷浓度更高,如华阳河湖群[33]和菜子湖[34],城郊闸控湖泊承担城市景观和生活用水,常年保持较高水位,从而抑制了沉积物的释放,极大的改变了城郊湖泊原有水文波动过程和水动力特性[35-36],从而影响营养物质的输移、沉降与释放过程,湖泊的营养状态自然随之改变[37]。石塘湖总磷浓度季节性变化就是由破罡闸闸站的季节性调度所致,水动力条件改变,最终影响石塘湖总磷的演变进程。

3.3.4 水文气象条件对中小型闸控湖泊总磷浓度的影响

水利工程仅仅是对湖泊磷的时空分布和演变特征产生影响,而气候变化则会加剧湖泊接受更多的营养物质,促进湖泊富营养化的发生[38],研究显示降雨量与TP呈显著正相关性[39],降雨增强会加大外源营养物质汇集;气象变化还改变入湖径流量,进而影响湖泊水位变化,致使湖泊营养程度改变。湖泊水文节律和换水周期是影响长江中下游湖泊营养盐的重要因素,丰水期破罡闸闸门调度时间提前和减慢水位下降速率,都有利于减小丰水期石塘湖水动力强度,从而降低丰水期对含磷污染物的迁移输运过程。

3.4 石塘湖总磷浓度与外源输入的相关性分析

2018年石塘湖及入湖河流总磷相关性系数均大于0.80,且满足显著性检验要求(P < 0.01),所有采样点总磷浓度均呈显著正相关,考虑到入湖河流周边土地利用类型存在较大区别,跃进大沟、606大沟主要受纳城市建成区来水,界河、余敦排洪河排口则主要受纳农村地区来水,表明石塘湖全湖总磷通量主要来自于外源人为活动的影响,受丰水期降雨农田面源冲刷汇集和城市汇流输入影响为主。

2018年石塘湖及入湖河流监测位点总磷浓度变化曲线基本呈现出一致的变化趋势,年内1-4月总磷浓度最低,4-12月总磷浓度呈现先上升后下降趋势,其中7月总磷浓度达到峰值,石塘湖与入湖河流的总磷平均浓度均为0.36 mg/L(图 4a)。结合表 2,石塘湖总磷负荷主要对石塘湖南部(S1号点位)影响最大,导致S1号点位总磷浓度相对较高,在1-8月最高值小于R1~R4 4条入湖河流总磷浓度最大值;9-11月丰水期结束,外源入湖总磷负荷大幅度减少,破罡闸闸门关闭,石塘湖水动力减弱,石塘湖总磷浓度因重力作用下含磷颗粒浓度沉降作用而降低,但是在相当长时间内,其浓度仍处于高值[40],S1~S4号位点稍微略高于外源入湖河流总磷浓度,12月份石塘湖绝大部分颗粒态磷基本已经沉淀,入湖含磷微团主要沉积在入湖河口,石塘湖和入湖河流总磷浓度都呈下降趋势。图 4b与c分别为采样点年内总磷浓度在不同位点时间变异系数和年内逐月不同采样点的空间变异系数,分别反映采样点总磷浓度的时间和空间差异程度,8个监测点年内总磷变异系数变化范围在52.1%~90.1%之间,时间上表现出异质性,对应图 4a中年内总磷浓度呈现较大幅度变化;1-12月不同采样点之间的总磷浓度空间变异系数变化范围为0.3%~16.2%,空间上则差异不显著,对应图 4a中年内各个采样点总磷变化趋势基本一致。

图 4 石塘湖和外源输入性总磷浓度随时间变化图 (b、c图为变异系数在时间和空间上的变化。变异系数=标准差/平均值) Fig.4 Changes of total phosphorus concentration and exogenous input total phosphorus concentration over time in Lake Shitang (b and c the variation coefficient of dimension in time and space. Coefficient of variation=standard deviation/mean)
表 2 石塘湖及入湖河流总磷浓度相关性系数 Tab. 2 Correlation coefficient of total phosphorus concentration in Lake Shitang and rivers entering the lake

根据相关研究表明,土地利用面积大小与河流水质存在显著相关性[41],土地利用方式的变化改变了流域自然水文过程,是影响受纳水体水质的主要因素之一[42]。从石塘湖与入湖河流总磷相关性及时空变异性结果可以合理推断,石塘湖总磷主要来自入湖河流输入,入湖河流总磷的季节性变化是影响湖泊水质变化的主要成因,而入湖河流水质变化又受流域内高强度的人类活动影响,高强度的人类活动使流域中不同汇水单元间土地利用呈现异质性,这不仅会导致土地利用类型发生变化,而且其空间分布也会强烈改变,这种变化常常造成不同汇水单元内的土地利用类型对水质的影响效果产生较大差异[43]

农田和建设用地是湖泊面源污染的主要来源,充当着“源”;坑塘水面和绿地由于能有效滞留、吸收污染物,一般充当着污染物进入湖泊的预前“汇”,其中农田单位面积的污染程度甚至远高于其他用地类型[44]。根据Yong等[45]研究俄亥俄州的小迈阿密河流域地表径流时发现,耕地为流域水质磷污染的主要来源,建设用地次之。相关数据表明本研究区域石塘湖农业用地主要集中在沿湖岸缓冲带,对入湖径流水质影响最大,城镇建设用地单位面积产生的单位负荷量远小于农业用地,城市管网的不健全和雨季城市地表径流冲刷下大量磷汇集于城郊湖泊;而林地和草地则远离石塘湖湖岸线,总体对地表径流磷起削减作用。根据金春玲等[46]发现洱海西部雨季地表径流磷污染特征与2 km外缓冲区尺度下林地面积呈显著负相关,林地面积占比增加有利于径流水质的改善,林木、林下植物和土壤可通过削减暴雨径流、减少水土流失、吸收截留污染物,有效减少磷的流失量[47],从流域总磷浓度变化曲线来看,受纳于农业面源(占流域总面积15.9%)的入湖河流界河、余敦排洪河排口与受纳于城市建设用地(占流域总面积的30.8%)的入湖河流跃进大沟、606大沟总磷年内变化曲线相差不大,由此可知,石塘湖流域在单位面积下农业用地产生单位污染负荷量远高于城镇建设用地总磷负荷量。因此,对于城郊中小型湖泊治理策略上主要集中在流域控源、水文调控、生态修复和食物链调控等技术,包括实施更为彻底的控源截污、流域自净体系构建策略等[48],减少源头输入和入湖河口湿地建设是控制外源负荷输入最有效措施。

4 结论

1) 石塘湖水质理化指标丰、平、枯水期的聚类效果显著,TP浓度主导着石塘湖水质时空变化,其他指标对石塘湖水质贡献较小。丰水期石塘湖总磷浓度达到峰值,远超地表水Ⅲ类水标准;枯水期TP浓度处于低值,水质较好,主要呈现时间上的差异,空间差异性不大。

2) 闸控湖泊常常兼具湖泊和河流2种形态特征,对湖泊总磷时空分布影响存在较大差异,在平水期和枯水期主要表现出阻隔湖泊特征,受纳和净化入湖污染物,外源输入磷主要以沉积作用为主;而丰水期开闸后更接近于河流特征,入湖含磷微团在丰水期水动力驱动下向全湖迁移输运,石塘湖总磷浓度变化与入湖河流总磷输入同步且差异较小,湖区空间差异不明显;高强度的降雨和闸站水位调度导致石塘湖水动力变化是造成石塘湖总磷浓度在丰水期达到峰值的主要驱动因子。

3) 石塘湖总磷主要来源于高强度人类活动的城乡区域,相关性显著。不同土地利用类型单位面积产生单位总磷负荷量大小来看:农业用地>城市建设用地,而林地和草地能有效滞留吸收污染物,充当着污染物进入湖泊的预前“汇”,丰水期入湖河流输入是导致湖泊总磷过高的决定性因素。

4) 关于沿江城郊闸控过水性湖泊治理途径来看,严格控制外源输入的同时降低丰水期水动力强度,以减轻外源入湖河流输入总磷负荷的影响和提升湖泊营养盐的有效截留能力,可以较好地减少总磷等负荷对水质的影响,为此类湖泊污染治理提供理论支撑。

5 附录

附图 Ⅰ见电子版(DOI: 10.18307/2023.0427)。

附图Ⅰ 石塘湖周边土地利用分布 AttachedFig1 Distribution of land use around Lake Shitang
附图Ⅱ 石塘湖流域各土地利用类型占比 AttachedFig2 The proportion of land use types in Lake Shitang Basin
6 参考文献

[1]
Ma RH, Yang GS, Duan HT et al. China's lakes at present: Number, area and spatial distribution. Science China Earth Sciences, 2011, 54(2): 283-289. DOI:10.1007/s11430-010-4052-6
[2]
Yang HF, Yang SL, Xu KH et al. Erosion potential of the Yangtze Delta under sediment starvation and climate change. Scientific Reports, 2017, 7: 10535. DOI:10.1038/s41598-017-10958-y
[3]
Cui YH, Wang J. Influence of the sluice on water level and area of Yangtze River-connected lakes: A case study in Shengjin Lake. Journal of Water Resources and Water Engineering, 2018, 29(4): 47-52. [崔玉环, 王杰. 控水闸对通江湖泊水位及面积变化的影响分析——以升金湖为例. 水资源与水工程学报, 2018, 29(4): 47-52.]
[4]
吴涛哲. 城市浅水型湖泊富营养化特征及界面交换研究[学位论文]. 上海: 华东交通大学, 2021.
[5]
Correll DL. The role of phosphorus in the eutrophication of receiving waters: A review. Journal of Environmental Quality, 1998, 27(2): 261-266. DOI:10.2134/jeq1998.00472425002700020004x
[6]
Daniel TC, Sharpley AN, Lemunyon JL. Agricultural phosphorus and eutrophication: A symposium overview. Journal of Environmental Quality, 1998, 27(2): 251-257. DOI:10.2134/jeq1998.00472425002700020002x
[7]
Camargo JA, Alonso A. Ecological and toxicological effects of inorganic nitrogen pollution in aquatic ecosystems: A global assessment. Environment International, 2006, 32(6): 831-849. DOI:10.1016/j.envint.2006.05.002
[8]
Reynolds CS. Phosphorus and the eutrophication of lakes—a personal view. Ciba Foundation Symposium, 1977(57): 201-228. DOI:10.1002/9780470720387.ch11
[9]
彭近新, 陈慧君. 水质富营养化与防治. 北京: 中国环境出版社, 1988.
[10]
Karl D, Letelier R, Tupas L et al. The role of nitrogen fixation in biogeochemical cycling in the subtropical North Pacific Ocean. Nature, 1997, 388(6642): 533-538. DOI:10.1038/41474
[11]
Zhang YL, Qin BQ. Evolution and research progress of eutrophication in Taihu Lake. Shanghai Environmental Science, 2001, 20(6): 3. [张运林, 秦伯强. 太湖水体富营养化的演变及研究进展. 上海环境科学, 2001, 20(6): 3.]
[12]
卢雅婷. 三峡水库蓄水以来氮、磷营养盐时空变化特征研究[学位论文]. 武汉: 武汉大学, 2019.
[13]
Xu Q, Yang TX, Liu XD et al. Analysis and prediction of eutrophication for Miyun reservoir. Journal of Changchun University of Science and Technology, 2003, 33(3): 315-318. [徐清, 杨天行, 刘晓端等. 密云水库总磷的富营养化分析与预测. 吉林大学学报: 地球科学版, 2003, 33(3): 315-318.]
[14]
Xia J. Research on the current situation of nitrogen and phosphorus pollution in pond culture water in Henan Province. Henan Agriculture, 2021(22): 61-63. [夏杰. 河南省池塘养殖水体氮磷元素污染现状研究. 河南农业, 2021(22): 61-63.]
[15]
Wang KM. Study on ecological restoration of Pogang Lake Wetland in Anqing City. Anhui Agricultural Science Bulletin, 2015, 21(13): 79-80. [王康明. 安庆市破罡湖湿地生态恢复研究. 安徽农学通报, 2015, 21(13): 79-80. DOI:10.3969/j.issn.1007-7731.2015.13.041]
[16]
Zhang QH, Dong XH, Chen SY et al. Quantitative reconstruction of the history of total phosphorus concentrations in Lake Shitang, Anhui Province over the past 100 years. Journal of Ecology and Rural Environment, 2013, 29(4): 494-499. [张清慧, 董旭辉, 陈诗越等. 安徽省石塘湖过去百余年湖水总磷浓度的定量重建. 生态与农村环境学报, 2013, 29(4): 494-499. DOI:10.3969/j.issn.1673-4831.2013.04.016]
[17]
Chu AX, Ding AS. Pollution source analysis and ecological restoration for Anqing Shitanghu alternate water. Shanxi Architecture, 2018, 44(21): 172-173. [储安西, 丁爱生. 安庆市石塘湖备用水源地污染源分析与生态修复. 山西建筑, 2018, 44(21): 172-173. DOI:10.3969/j.issn.1009-6825.2018.21.093]
[18]
Bouguerne A, Boudoukha A, Benkhaled A et al. Assessment of surface water quality of Ain Zada Dam (Algeria) using multivariate statistical techniques. International Journal of River Basin Management, 2017, 15(2): 133-143. DOI:10.1080/15715124.2016.1215325
[19]
Qi JY, Ruan XH. Dam construction-induced environmental impact on riverine ecosystem. Journal of Hehai University: Natural Sciences, 2005, 33(1): 37-40. [祁继英, 阮晓红. 大坝对河流生态系统的环境影响分析. 河海大学学报: 自然科学版, 2005, 33(1): 37-40.]
[20]
Xu L, Li H, Chen YX et al. Research on seasonal variation of self-purification ability for small shallow lakes in south lake Taihu. Environmental Science, 2010, 31(4): 924-930. [许磊, 李华, 陈英旭等. 南太湖地区小型浅水湖泊自净能力季节变化研究. 环境科学, 2010, 31(4): 924-930. DOI:10.13227/j.hjkx.2010.04.029]
[21]
Zhang XY, Ma PM, Gao QH et al. Temporal and spatial variation characteristics of water quality in upstream and downstream of the three gorges dam. J Lake Sci, 2019, 31(3): 633-645. [张馨月, 马沛明, 高千红等. 三峡大坝上下游水质时空变化特征. 湖泊科学, 2019, 31(3): 633-645. DOI:10.18307/2019.0303]
[22]
Yao HM, Huang RT, Gan FX et al. Principal component analysis of the water quality evaluation in east lake. Geomatics and Information Science of Wuhan University, 2005, 30(8): 732-735. [姚焕玫, 黄仁涛, 甘复兴等. 用改进的主成分分析法对东湖的水质污染进行评价. 武汉大学学报·信息科学版, 2005, 30(8): 732-735.]
[23]
Yin W, Xin XK, Liang JK et al. Evaluation of water quality for tributaries of Danjiangkou reservoir with principal component analysis. Water Resources and Power, 2015, 33(1): 34-38. [尹炜, 辛小康, 梁建奎等. 基于主成分分析的丹江口水库支流水质评价. 水电能源科学, 2015, 33(1): 34-38.]
[24]
龚珞军. 五个湖泊浮游生物和主要理化指标的比较研究[学位论文]. 武汉: 华中农业大学, 2011.
[25]
Huang YL, Ji DB, Chen MX et al. Effect of pH value of water body on the generation and elimination of cyanobacterial blooms. Yangtze River, 2008, 39(2): 63-65. [黄钰铃, 纪道斌, 陈明曦等. 水体pH值对蓝藻水华生消的影响. 人民长江, 2008, 39(2): 63-65. DOI:10.16232/j.cnki.1001-4179.2008.02.008]
[26]
Ji PF, Xu H, Zhan X et al. Spatial-temporal variations and driving of nitrogen and phosphorus ratios in lakes in the middle and lower reaches of Yangtze River. Environmental Science, 2020, 41(9): 4030-4041. [季鹏飞, 许海, 詹旭等. 长江中下游湖泊水体氮磷比时空变化特征及其影响因素. 环境科学, 2020, 41(9): 4030-4041. DOI:10.13227/j.hjkx.202002017]
[27]
杜广龙. 太湖沉积物再悬浮及营养盐释放规律机制研究[学位论文]. 南京: 南京林业大学, 2013.
[28]
Zhu WJ, Pang WT, You QM et al. Phytoplankton community structure and the evaluation of water quality in spring, Huaihe River Basin. J Lake Sci, 2017, 29(3): 637-645. [朱为菊, 庞婉婷, 尤庆敏等. 淮河流域春季浮游植物群落结构特征及其水质评价. 湖泊科学, 2017, 29(3): 637-645. DOI:10.18307/2017.0312]
[29]
Zhang OY, Xiong W, Ding HL. Drainage connectivity characteristics and influential factors of Yangtze River Basin. Yangtze River, 2010, 41(1): 1-5, 78. [张欧阳, 熊文, 丁洪亮. 长江流域水系连通特征及其影响因素分析. 人民长江, 2010, 41(1): 1-5, 78.]
[30]
Zhu W, Hu SY, Feng GY et al. Effects of great floods on phosphorus in shallow lakes: A case study of Lake Taihu in 2016. J Lake Sci, 2020, 32(2): 325-336. [朱伟, 胡思远, 冯甘雨等. 特大洪水对浅水湖泊磷的影响: 以2016年太湖为例. 湖泊科学, 2020, 32(2): 325-336. DOI:10.18307/2020.0201]
[31]
Zhu GW. Eutrophic status and causing factors for a large, shallow and subtropical Lake Taihu, China. J Lake Sci, 2008, 20(1): 21-26. [朱广伟. 太湖富营养化现状及原因分析. 湖泊科学, 2008, 20(1): 21-26. DOI:10.18307/2008.0103]
[32]
Lv XY, Wu SQ, Zhang Y et al. Analysis on variation of main indicators of eutrophication and mutrition level in Taihu Lake. Journal of Water Resources and Water Engineering, 2014, 25(4): 1-6. [吕学研, 吴时强, 张咏等. 太湖富营养化主要指标及营养水平变化分析. 水资源与水工程学报, 2014, 25(4): 1-6. DOI:10.11705/j.issn.1672-643X.2014.04.001]
[33]
Zeng K, Wang JS, Zhang YC et al. Influence of water level fluctuation on water quality of Huayang River and lake group. Journal of Yangtze River Scientific Research Institute, 2020, 37(8): 49-53, 68. [曾凯, 王家生, 章运超等. 华阳河湖群水位变化对水质的影响. 长江科学院院报, 2020, 37(8): 49-53, 68. DOI:10.11988/ckyyb.20190419]
[34]
Lesheng AN, Kaihua L, Lei Z et al. Influence of river-lake isolation on the water level variations of Caizi Lake, lower reach of the Yangtze River. Journal of Geographical Sciences, 2021, 31(4): 14. DOI:10.1007/S11442-021-1858-4
[35]
Chen RT, Xie ZY. The evaluation of eutropher on cihu lake. Journal of Huangshi Polytechnic College, 2004, 20(4): 61-63. [陈瑞庭, 谢志勇. 磁湖富营养化评价. 黄石高等专科学校学报, 2004, 20(4): 61-63.]
[36]
Lu Z, Li Y, Ma X et al. Surface water quality assessment in period 2000-2009 and forecasting on changing tendency of pollution in Daye Lake. Environmental Science & Technology, 2012, 35(5): 174-178. [鲁珍, 李晔, 马啸等. 大冶湖2000-2009年地表水质评价及污染趋势预测. 环境科学与技术, 2012, 35(5): 174-178.]
[37]
Su JQ, Wang X, Yang ZF. Lake eutrophication modeling in considering climatic factors change: A review. Chinese Journal of Applied Ecology, 2012, 23(11): 3197-3206. [苏洁琼, 王烜, 杨志峰. 考虑气候因子变化的湖泊富营养化模型研究进展. 应用生态学报, 2012, 23(11): 3197-3206. DOI:10.13287/j.1001-9332.2012.0477]
[38]
Schindler DW. Recent advances in the understanding and management of eutrophication. Limnology and Oceanography, 2006, 51(1part2): 356-363. DOI:10.4319/lo.2006.51.1_part_2.0356
[39]
Zuo SP, Ye LT, Mei H. Driving effects of regional climate on dynamics of water environment in eastern half of Chaohu Lake. Bulletin of Soil and Water Conservation, 2011, 31(4): 28-31. [左胜鹏, 叶良涛, 梅慧. 气候对巢湖东半湖水体环境变化的驱动效应. 水土保持通报, 2011, 31(4): 28-31. DOI:10.1016/S1671-2927(11)60313-1]
[40]
Peng JP, Pang Y, Li YP et al. Change of phosphorus in lake after water dynamical condition and its contribution to eutrophication. Ecology and Environmnet, 2004, 13(4): 503-505. [彭进平, 逄勇, 李一平等. 水动力过程后湖泊水体磷素变化及其对富营养化的贡献. 生态环境, 2004, 13(4): 503-505.]
[41]
Ren WW, Zhong Y, Meligrana J et al. Urbanization, land use, and water quality in Shanghai: 1947-1996. Environment International, 2003, 29(5): 649-659. DOI:10.1016/S0160-4120(03)00051-5
[42]
Changnon SA, Demissie M. Detection of changes in streamflow and floods resulting from climate fluctuations and land use-drainage changes. Climatic Change, 1996, 32(4): 411-421. DOI:10.1007/BF00140354
[43]
Yang L, Ma KM, Guo QH et al. Zoning planning in non-point source pollution control in Hanyang district. Environmental Science, 2006, 27(1): 31-36.
[44]
Jiang QH, Liu YF, Huang PJ et al. The source-sink effect between urban diffuse pollution and Lake watershed: A case study in East Lake of Wuhan. Ecology and Environmental Sciences, 2013, 22(3): 469-474. [姜庆虎, 刘艳芳, 黄浦江等. 城市湖泊流域面源污染的源-汇效应研究——以武汉市东湖为例. 生态环境学报, 2013, 22(3): 469-474.]
[45]
Tong STY, Chen W. Modeling the relationship between land use and surface water quality. Journal of Environmental Management, 2002, 66(4): 377-393. DOI:10.1006/jema.2002.0593
[46]
Jin CL, Gao SJ, Ye BB et al. Nitrogen and phosphorus pollution characteristics of surface runoff and the impacts of land use on runoff water quality in rainy season in the western Erhai Lake Basin. Environmental Science Research, 2018, 31(11): 1891-1899. [金春玲, 高思佳, 叶碧碧等. 洱海西部雨季地表径流氮磷污染特征及受土地利用类型的影响. 环境科学研究, 2018, 31(11): 1891-1899. DOI:10.13198/j.issn.1001-6929.2018.03.25]
[47]
Floyd WC, Schoenholtz SH, Griffith SM et al. Nitrate-nitrogen, land use/land cover, and soil drainage associations at multiple spatial scales. Journal of Environmental Quality, 2009, 38(4): 1473-1482. DOI:10.2134/jeq2008.0099
[48]
Qin BQ. Lake eutrophication: Control countermeasures and recycling exploitation. Ecological Engineering, 2009, 35(11): 1569-1573. DOI:10.1016/j.ecoleng.2009.04.003