环境科学研究  2020, Vol. 33 Issue (5): 1170-1178  DOI: 10.13198/j.issn.1001-6929.2020.04.22

引用本文  

范中亚, 王文才, 蒋锦刚, 等. 华阳河湖群沉积物内源磷释放风险及控制策略[J]. 环境科学研究, 2020, 33(5): 1170-1178.
FAN Zhongya, WANG Wencai, JIANG Jingang, et al. Risk and Control Strategy of Internal Phosphorus Release from Sediments in Huayang Lakes[J]. Research of Environmental Sciences, 2020, 33(5): 1170-1178.

基金项目

长江生态环境保护修复联合研究项目(第一期)(No.2019-LHYJ-01-0212,2019-LHYJ-01-0212-17);长江生态保护与修复安全驻点研究项目(No.2018CJA030301-014)
Supported by Joint Research Project for the Yangtze River Conservtion (Phase Ⅰ), China (No.2019-LHYJ-01-0212, 2019-LHYJ-01-0212-17); National Yangtze River Conservation and Restoration Joint Research Project China (No.2018CJA030301-014)

责任作者

蒋锦刚(1981-), 男, 安徽桐城人, 助理研究员, 博士, 主要从事湖泊环境保护研究, gangzg@ipp.ac.cn; 曾凡棠(1963-), 男, 广东河源人, 研究员, 硕士, 主要从事流域水环境治理研究, zengfantang@scies.org.

作者简介

范中亚(1985-), 男, 安徽舒城人, 高级工程师, 硕士, 主要从事流域水环境治理研究, fanzhongya@scies.org

文章历史

收稿日期:2020-02-15
修订日期:2020-04-13
华阳河湖群沉积物内源磷释放风险及控制策略
范中亚1, 王文才1, 蒋锦刚2, 马金玉1, 王钟1, 曾凡棠1    
1. 生态环境部华南环境科学研究所, 国家环境保护水环境模拟与污染控制重点实验室, 广东 广州 510655;
2. 中国科学院合肥物质科学研究院技术生物与农业工程研究所, 安徽 合肥 230031
摘要:为研究沉积物内源磷释放对浅水湖泊水体总磷的影响,以华阳河湖群为对象,根据不同季节全湖水质监测数据,分析总磷的空间和季节性分布特征;结合沉积物间隙水磷浓度数据,采用浓度梯度扩散模型估算湖泊沉积物内源磷静态释放量;基于2019年在线水质监测站连续数据,分析典型沉积物再悬浮过程内源磷释放对水体总磷贡献及估算其年累积释放量.结果表明:①华阳河湖群现场监测总磷浓度变化范围为0.014~0.297 mg/L,龙感湖与黄大湖总磷浓度空间变幅大于泊湖;②典型再悬浮过程中内源磷动态释放对水体总磷的贡献可达75%,华阳河湖群内源磷的年释放量在1.129×104~1.684×104 t之间;③静态的内源磷年释放量为12.92 t,仅占动态释放量0.1%.研究显示:华阳河湖群总磷浓度时空差异较大,内源磷释放对总磷短期波动具有重要作用,静态释放量对全年内源磷释放量的贡献几乎可以忽略;严格控制入湖营养盐通量的同时,恢复草型湖泊生态系统以控制内源释放已成为解决华阳河湖群水环境问题的重要管理决策方向.
关键词华阳河湖群    内源释放    沉积物再悬浮    湖泊水质保护策略    
Risk and Control Strategy of Internal Phosphorus Release from Sediments in Huayang Lakes
FAN Zhongya1, WANG Wencai1, JIANG Jingang2, MA Jinyu1, WANG Zhong1, ZENG Fantang1    
1. National Key Laboratory of Water Environmental Simulation and Pollution Control, South China Institute of Environmental Sciences, Ministry of Ecology and Environment, Guangzhou 510655, China;
2. Institute of Technical Biology & Agriculture Engineering, Hefei Institutes of Physical Science, Chinese Academy of Sciences, Hefei 230031, China
Abstract: In order to investigate the effect of internal phosphorus release on total phosphorus concentration in shallow lakes, this article takes Huayang Lakes as the research object, the spatial and seasonal distribution characteristics of total phosphorus were analyzed according to the monitoring data of the whole lakes in different seasons. The amount of static internal phosphorus release in lake sediment was estimated by the concentration gradient diffusion model with phosphorous concentration in pore water and in sediment. The contribution of internal phosphorous release to total phosphorous and annual cumulative release during sediment resuspension process was analyzed based on the time series data on-line water quality monitoring station in 2019. The results showed: (1) The spatial temporal monitored total phosphorus concentration of the Huayang Lakes carried in a range of 0.014-0.297 mg/L, and the spatial variation of total phosphorus concentration in Longgan Lake and Huangda Lake was higher than that in Po Lake. (2) The contribution of dynamic internal release of total phosphorus could reach 75% in a typical resuspension process, and the annual internal phosphorus release in Huayang Lakes was about 1.129×104-1.684×104 t. (3) The static release was about 12.92 t/a, accounting for only 0.1% of the dynamic internal release. This study revealed that the spatial temporal total phosphorus concentration in Huayang Lakes changed remarkably, internal phosphorus release was significance for short-term sequence variation of total phosphorus, and the contribution of static internal phosphorus was limited. Strictly controlling the nutrient flux from rivers into lakes and restoring macrophytic lake ecosystem have become the important measures for solving the water environmental problems of Huayang Lakes.
Keywords: Huayang Lakes    internal release    sediment resuspension    protection strategy of lake water quality    

浅水湖泊受风浪扰动下易发生沉积物再悬浮及内源磷释放现象;而在动力扰动较弱时,依然存在浓度梯度作用下的静态释放风险.沉积物再悬浮引起的水下光场改变与内源磷释放引起的水体磷浓度改变,会对湖泊初级生产力(水生植被、浮游藻类)的生长产生重要影响[1].因此,揭示内源磷释放规律、准确估算内源磷释放规模及对其上覆水的影响,对于制定科学、有效的湖泊管理策略具有重要意义,也是湖泊富营养化相关研究的热点.

沉积物内源释放主要有静态释放和动态释放两种模式.静态释放依靠沉积物间隙水与上覆水营养盐浓度差的分子扩散作用,作为风浪扰动较小条件下营养盐释放的主要方式,释放通量会受到氧化还原环境、pH、温度以及铁、锰含量等影响[2];动态释放指受动力扰动时在沉积物-水界面切应力作用下发生的沉积物再悬浮及营养盐释放,释放规模会受到风浪强度、水深、底质特性影响.研究表明,太湖沉积物内源磷动态释放量达2×104 t/a,远高于外源输入和静态释放量[1];沉积物再悬浮及营养盐释放显著影响了太湖富营养化湖泊水环境特征及水华生消过程[3].在湖泊沉积物再悬浮状态下的内源磷动态释放以颗粒态为主,溶解态磷占比较小[4-6].已有研究多以室内或短期高频野外原位试验为基础,难以解耦风速风向不稳定下的长期内源磷动态释放规律和规模.大量在线水质观测站的建设为研究人员利用长期水质监测数据分析内源磷动态释放对水体总磷的贡献奠定了基础.

长江中下游地区的湖泊约占我国淡水湖泊总面积的60%,且绝大多数为浅水湖泊[7],基本没有明显的热力学分层现象,湖水会经常在垂向上发生交换作用[8].近年来,环保基础设施的滞后以及湖泊资源过度开发造成大量营养盐进入湖泊水体和蓄积在湖底沉积物之中,导致湖泊富营养化成为长江中下游湖泊面临的主要问题,且众多湖泊表现出总磷浓度偏高的现象.因此,研究这些浅水湖泊的内源磷释放规律、规模和影响具有十分重要的意义.华阳河湖群是长江蓄滞洪区之一,季节性水位波动大,湖泊总磷本底含量高[9].近年来,华阳河湖群富营养呈增长趋势,且湖泊总磷季节性波动幅度较大[10-11],总磷浓度是控制湖泊富营养化趋势的主要因子之一.开展华阳河湖群沉积物内源磷释放风险研究,探讨内源释放对湖泊水体磷的影响,以期为华阳河湖群水环境管理提供决策依据.

1 材料与方法 1.1 流域概况

华阳河湖群流域面积为5 511.4 km2,属于北亚热带湿润气候,年均温度为16.6 ℃,年降水量为1 410 mm.华阳河湖群包括龙感湖、黄大湖、泊湖,其中,龙感湖、黄大湖、泊湖面积分别为420、299、180 km2,湖底高程为9.5~10.5 m,龙感湖、黄大湖、泊湖湖底高程呈逐渐下降趋势,湖泊多年平均水位为12.93 m.湖泊北部分布着三条主要入湖河流,其中两条(新县河、二郎河)汇入龙感湖,一条(凉亭河)汇入泊湖.随着湖区农业及水产业等发展,湖泊湿地面积减少了37%,水生植被覆盖度下降明显.湖区大面积的围栏养殖造成沉积物内源负荷不断增加. 2016年龙感湖、黄大湖总磷年均值分别为0.06、0.04 mg/L,季节性波动幅度较大. 2018年龙感湖已接近中度富营养化[12]. 2018年华阳河湖群的网箱全部被拆除,养殖投饵产生的污染负荷不再增加,但蓄积在湖泊沉积物中的营养盐还没有消除,一旦发生沉积物再悬浮,水体中的营养负荷将大幅增加.

1.2 数据来源

该研究所用数据分为四部分:①2014年9月实测的华阳河湖群空间27个点位沉积物间隙水与上覆水营养盐浓度,用于计算磷的静态释放通量;②2018年8月(夏季)、2018年10月(秋季)、2019年1月(冬季)、2019年4月(春季)现场监测点如图 1所示,其监测的湖泊水质数据用于分析水质指标的季节性和空间变化特征,监测指标包括悬浮物(SS)、总磷(TP)等;③2019年水质在线监测站(L4监测点)水质数据,水质数据来源于龙感湖水质监测自动站,包括浊度、总磷等指标;④气象数据主要包括逐日的风速和风向,来源于国家气象数据共享平台(http://data.cma.cn/site/index.html).

注:L1~L7为龙感湖监测点,H1~H11为黄大湖监测点,P1~P4为泊湖监测点. 图 1 华阳河湖群流域监测点位置 Fig.1 Location map of monitoring points in Huayang Lakes
1.3 样品采集与处理

沉积物柱状样采用柱状采样器采集,采集过程中尽量不扰动沉积物.柱状沉积物现场按5 cm进行分层并分装于聚乙烯瓶中,同步采集上覆水储存于用5%的HNO3处理过的聚乙烯瓶中.样品运回实验室做进一步处理,取部分沉积物鲜样用烘干法测定容重和孔隙度,剩余沉积物用离心法获取间隙水.上覆水和间隙水样品经孔径为0.45 μm的玻璃纤维滤膜过滤,并在24 h之内利用钼酸铵分光光度法完成溶解性正磷酸盐(PO43--P)、可溶性总磷酸盐(DTP)和可溶性无机磷(DIP)的测定.

华阳河湖群水质现场调查于2018年8月、2018年10月、2019年1月、2019年4月在华阳河湖群22个点位(见图 1)采集表层水样,并在现场利用多参数水质分析仪监测水温、pH、浊度和氧化还原电位.

1.4 分析方法 1.4.1 静态释放

采用修正的Fick′s第一扩散定律[13]对华阳河湖群沉积物-水界面磷扩散通量进行计算,计算公式:

$ F = - \varphi \times {D_{\rm{S}}}\left( {{\delta _{\rm{C}}}/{\delta _{\rm{S}}}} \right) $ (1)

式中:F为沉积物-水界面营养盐扩散速率,F为正值代表营养盐通量是从沉积物向上覆水扩散,其为负值代表营养盐通量是从上覆水向沉积物扩散,mg/(d·m2);φ为表层沉积物的孔隙度,%;δC/δS为浓度梯度,mg/(L·cm);DS为沉积物中溶质的分子扩散系数,cm2/s, 由于DS是考虑了沉积物的弯曲效应,通常通过其与孔隙度之间的关系来进行推导[14]

$ {D_{\rm{S}}} = \varphi \times {D_0}\quad (\varphi < 0.7) $ (2)
$ {D_{\rm{s}}} = {\varphi ^2} \times {D_0}\quad (\varphi \ge 0.7) $ (3)

华阳河湖群表层沉积物φ在0.37~0.53之间,取DS=φ×D0D0为理想溶液的扩散系数,该研究取值为5×10-6 cm2/s.

1.4.2 风浪引起的动态释放过程分析 1.4.2.1 沉积物再悬浮过程内源占比分析

研究[1-2]表明,风浪扰动下沉积物再悬浮会引起内源营养盐的动态释放,显著增加水体中的颗粒态磷浓度.一次完整的沉积物再悬浮及内源释放过程会随着风浪扰动的增强而增强,呈现悬浮颗粒物浓度先升后降的情况.因此,可以将风浪扰动前相对平静状态下的总磷浓度作为背景浓度,悬浮颗粒物和总磷浓度明显升高的时刻作为沉积物再悬浮及动态释放的起始时刻;风浪扰动后相对平静,悬浮颗粒物、总磷浓度相对平稳或恢复到扰动前水平的时刻作为沉积物再悬浮及动态释放的结束时刻.据此计算在一个完整的沉积物再悬浮过程中总磷动态释放对水体中总磷浓度的贡献(内源占比).

$ \omega (t) = \frac{{{c_{{\rm{TP}}}}(t) - {c_{{\rm{TP}}}}^\prime }}{{{c_{{\rm{TP}}}}(t)}} \times 100\% $ (4)
$ {c_{{\rm{TP}}}}^\prime = \frac{1}{n}\sum\limits_{i = 1}^n {{c_i}} $ (5)

式中:ω(t)为t时刻内源动态磷释放量占当前时刻水体中磷浓度的权重百分比,%;cTP′为风浪扰动引起沉积物再悬浮前的总磷背景平均浓度,mg/L;cTP(t)为t时刻水体中总磷浓度,mg/L;ci为风浪扰动前用于计算平均背景浓度时刻的总磷浓度,mg/L; i为风浪扰动前被选取用于计算平均背景浓度的数据序号.

1.4.2.2 年再悬浮动态释放量计算

该研究将总磷浓度、浊度值在小时尺度同步增高,且总磷浓度增高值大于某一阈值的时刻判定为再悬浮释放状态,即同时满足式(6)~(9).

$ {\Delta {c_{{\rm{TP}}}}(\Delta t) = {c_{{\rm{TP}}}}(t) - {c_{{\rm{TP}}}}(t - \Delta t) > 0} $ (6)
$ {\Delta {c_{{\rm{TP}}}}(2\Delta t) = {c_{{\rm{TP}}}}(t) - {c_{{\rm{TP}}}}(t - 2\Delta t) > {c_{{\rm{thr}}}}} $ (7)
$ {\Delta {\mathop{\rm Turb}\nolimits} (\Delta t) = {\mathop{\rm Turb}\nolimits} (t) - {\mathop{\rm Turb}\nolimits} (t - \Delta t) > 0} $ (8)
$ {\Delta {\mathop{\rm Turb}\nolimits} (2\Delta t) = {\mathop{\rm Turb}\nolimits} (t) - {\mathop{\rm Turb}\nolimits} (t - 2\Delta t) > 0} $ (9)

式中:ΔcTPt)为总磷浓度增高值,mg/L;cTP(tt)为前一时刻的总磷浓度,mg/L;ΔcTP(2Δt)为间隔两个时刻的总磷浓度增高值,mg/L;cTP(t-2Δt)为前两个时刻的总磷浓度,mg/L; Turb(t)为t时刻的浊度,NTU;ΔTurb(Δt)为间隔一个时刻的水体浊度增高值,NTU;Turb(tt)为前一时刻的水体浊度值;ΔTurb(2Δt)为间隔两个时刻的水体浊度增高值,NTU;Turb(t-2Δt)为前两个时刻的水体浊度值, NTU; cthr为总磷变化阈值,mg/L,将其设定为0~0.05区间,按照0.001间隔分别计算.

在线监测点年再悬浮量计算可根据再悬浮状态下总磷浓度相比前一采测时刻的增量与单位面积上水体体积的乘积累加,计算公式:

$ {\rm{PRRA}} = \sum\limits_{j = 1}^m {\left[ {\Delta {c_{{\rm{TP}}}}\left( {\Delta {t_j}} \right) \times d\left( {{t_j}} \right)} \right]} $ (10)

式中:PRRA为每m2沉积物每年总磷释放量,g/(m2·a);m为一年中所有处于再悬浮释放状态的事件数量;ΔcTPtj)为第j个再悬浮状态下相比前一时刻的总磷浓度增量,mg/L;d(tj)为第j个再悬浮状态对应时刻的水深,m.

2 结果与分析 2.1 湖泊水体磷分布

图 2可见:华阳河湖群总磷浓度季节性变化较大,表现为秋季 < 夏季 < 春季 < 冬季,变化范围为0.014~0.297 mg/L,龙感湖与黄大湖总磷浓度空间变幅大于泊湖.夏季与秋季总磷浓度较低,夏季总磷浓度在0.023~0.110 mg/L之间,平均值为0.058 mg/L,龙感湖、黄大湖和泊湖秋季总磷浓度平均值分别为0.048、0.022和0.022 mg/L;冬季与春季总磷浓度较高,龙感湖、黄大湖和泊湖冬季总磷浓度平均值分别为0.183、0.164和0.087 mg/L, 春季龙感湖、黄大湖和泊湖平均值分别为0.081、0.089和0.078 mg/L.

图 2 华阳河湖群总磷浓度季节性分布特征 Fig.2 Seasonal distribution of total phosphorus in Huayang Lakes
2.2 静态释放量估算

基于4次现场监测结果发现:华阳河湖群沉积物间隙水中可溶性总磷(DTP)浓度在0.04~0.19 mg/L之间,平均值为0.10 mg/L;间隙水无机磷(DIP)浓度在0.008~0.059 mg/L之间,平均值为0.02 mg/L.华阳河湖群沉积物中的可溶性磷浓度高于上覆水的1~3倍,这也为沉积物中的磷释放到上覆水提供了必要条件.

华阳河湖群表层沉积物磷扩散速率在0.025 3~0.064 8 mg/(d·m2)之间,日释放通量为(0.035 4±0.025 9)t(见表 1).总体上,表层沉积物中的磷有从沉积物向上覆水释放的可能性,即表层沉积物表现为磷的“源”.各湖区磷扩散通量平均值大小顺序为黄大湖>泊湖>龙感湖,其中黄大湖的日释放通量约为龙感湖和泊湖的2.4倍.磷扩散速率的最大值和最小值分别出现在黄大湖和龙感湖,磷扩散速率相对较高的地区主要集中在黄大湖西北部.经估算,华阳河湖群磷的年静态释放通量约为12.92 t.

表 1 黄大湖、龙感湖、泊湖沉积物-水界面磷扩散速率及日释放通量 Table 1 Phosphorus diffusion rate and daily release flux at the sediment-water interface of Huangda Lake, Longgan Lake and Po Lake
2.3 动态释放量计算 2.3.1 风速、风向、浊度与总磷的关系

图 3为2019年9月15日—12月30日风速、风向、浊度和总磷浓度变化的时间序列图.由图 3可见:浊度与总磷浓度存在明显的同步性,总磷浓度随浊度同步变化;观测期浊度有几个明显的先升后降的过程,风速同样呈先升后降的特征,且风向相对稳定为(东北风).在一个相对稳定的风向下,浊度和总磷浓度随风速的增加而略有滞后的上升;风速减小后,浊度与总磷浓度会随之下降,表现出明显的再悬浮及磷的释放过程.

注:Ⅰ、Ⅱ、Ⅲ、Ⅳ表示4个典型再悬浮过程. 图 3 华阳河湖群风速、风向、浊度与总磷浓度随时间变化 Fig.3 Temporal variety of wind speed, wind direction, turbidity and total phosphorus in Huayang Lakes

图 4为华阳河湖群4个典型再悬浮过程中风速、风向频率及风速与总磷浓度的关系.由图 4可见:4个再悬浮过程中,风向主要集中在45°方向,风速与总磷浓度变化存在时间滞后性,相关系数峰值滞后时间为10~40 h,滞后时长主要受风向、湖泊水位等影响;另外,4个典型再悬浮过程时段内浊度与总磷浓度的相关系数为0.939(P < 0.01).

图 4 华阳河湖群风速、风向及风速与总磷浓度的关系 Fig.4 Wind speed, wind direction and correlation between wind speed and total phophorus in Huayang Lakes
2.3.2 再悬浮过程中磷的释放动态过程

图 5为4个典型再悬浮过程中水体总磷浓度及内源释放贡献率的变化过程.由图 5可见:再悬浮期间内源释放的总磷浓度最高可占水体总磷浓度的75%左右,后期逐渐沉降后出现贡献率为负的现象.典型再悬浮过程Ⅰ和Ⅱ是两个相对独立的再悬浮及磷释放过程,呈总磷浓度先升后降的单峰特征,在两个相对独立的再悬浮过程后期,总磷贡献率为负的情况不明显;典型再悬浮过程Ⅲ和Ⅳ是两个相对连续的再悬浮及磷释放过程,总磷浓度变化规律复杂,在第2个再悬浮过程后期,总磷浓度均出现下降趋势,且总磷贡献率为负的情况较明显.

图 5 典型再悬浮过程总磷浓度及再悬浮磷对总磷贡献随时间变化 Fig.5 Temporal variety of total phosphorus and contribution of internal phosphorous release during typical resuspension process
2.3.3 全年再悬浮磷释放量计算

图 6(a)可见,单位面积沉积物全年内源磷动态释放量随总磷浓度变化阈值的增大而减小,当总磷浓度变化阈值介于0~0.05 mg/L时,内源磷动态释放量在12.558 4~18.735 5 g/(m2·a)之间.研究[15]发现,华阳河湖群各湖区沉积物均以细颗粒泥沙为主,各湖区表层沉积物全磷含量总体差别较小,采用龙感湖在线站再悬浮结果对全湖进行估算有一定的参考,按899 km2的湖泊面积估算,华阳河湖群磷的年再悬浮动态释放量约为1.129×104~1.684×104 t.由图 6(b)可见,总磷浓度变化阈值为0.02 mg/L时,春季和冬季(枯水期)内源磷动态释放量占比明显高于夏季和秋季(丰水期),春季与冬季内源磷动态释放占全年内源磷动态释放量的84.3%,与全湖现场监测的总磷浓度季节性变化基本一致.

图 6 单位面积沉积物内源磷动态释放量及占比 Fig.6 Annual internal phosphorus release per unit area from sediment
3 讨论 3.1 沉积物内源释放对湖泊水质的影响

华阳河湖群表层沉积物磷扩散速率在0.025 3~0.064 8 mg/(d·m2)之间,日释放通量为(0.035 4±0.025 9)t,年释放通量约为12.92 t,与长江流域其他浅水湖泊单位面积磷静态释放通量相比[2, 16],磷静态释放量呈太湖>华阳河湖群>巢湖的特征,太湖和巢湖都是富营养化浅水湖泊,而华阳河湖群也面临富营养化趋势.由此可见,华阳河湖群内源静态释放量不容忽视,是浅水湖泊富营养化控制要考虑的因素之一.

图 3中显示在东北风条件下,华阳河湖群悬浮物浓度、总磷浓度的变化滞后于风速的变化,悬浮物与总磷浓度存在明显的正相关.风强迫作用下的波切应力是浅水湖泊沉积物再悬浮的主要驱动力[17-18];野外观测研究[4, 6, 19]发现,水体中悬浮物及总磷浓度受风浪引起的沉积物再悬浮影响波动明显,再悬浮过程中主要是颗粒态磷浓度的变化;室内波浪水槽模拟结果[20-21]显示,沉积物再悬浮过程水体中溶解态磷变化不显著.浊度与悬浮颗粒物浓度通常具有较好的线性回归关系[22],龙感湖风速、风向与浊度和总磷浓度的关系也反映了风浪对沉积物再悬浮及磷释放的影响. 4个典型再悬浮过程中内源磷释放量对水体总磷浓度贡献率较高,最高达75%;而再悬浮后期会出现贡献率为负的现象,水体中总磷浓度比背景磷浓度低,并逐渐下降.该研究计算的华阳河湖群总磷释放量为1.129×104~1.684×104 t,与秦伯强等[2]对太湖全年总磷释放量估算结果(2.1×104 t)在一个数量级,也表明华阳河湖群沉积物再悬浮引起的总磷释放是引起湖泊水体总磷浓度变化的主导因素.由于未考虑各湖区动力条件不同引起的沉积物再悬浮及总磷动态释放时空变异性,全湖总磷释放量的计算结果存在一定的偏差,但基本能反映沉积物再悬浮下总磷动态释放量的量级,且远高于其静态释放量.

季节性沉积物内源磷动态释放对比显示,春、冬两季(枯水期)沉积物内源磷释放量远高于夏、秋两季(丰水期),枯水期较低的水深会加重内源磷释放对水体总磷的影响,全湖不同季节现场监测结果也证实了春、冬两季水体总磷浓度明显高于夏、秋两季. 4个典型沉积物再悬浮及内源磷释放过程分析结果显示,沉积物内源磷动态释放对上覆水体风浪期的短期总磷浓度具有决定性的影响,因此风浪刚结束的冬季全湖采样监测的总磷浓度明显高于其他季节.研究[6, 20, 23]表明,风速、风向(风区长度)、水深、沉积物特性都会对沉积物再悬浮强度造成影响;另外,由于不同点位受到的扰动强度不同,也造成空间湖泊总磷浓度在冬季和春季的空间差异性更明显,在风浪扰动较小的夏季和秋季,湖泊总磷浓度空间差异相对较小,这时总磷浓度变化主要受控于入湖河流及浮游植物等初级生产力的影响.

3.2 华阳河湖群水质保护策略分析

华阳河湖群主要的水质问题是季节性总磷浓度偏高.由于华阳河湖群属于长江古河道,湖泊沉积物本底含量高[9];另外,冬季与春季湖泊水体总磷浓度普遍高于夏季,与吴奇丽等[10-11]研究结论一致.因此,华阳河湖群水质保护核心之一是控制湖泊水体总磷浓度波动,特别是解决春季与冬季水体总磷浓度偏高的问题.

湖泊总磷浓度控制分外源与内源污染源控制.外源主要来自入湖河流及沿岸面源,雨季的外源输入磷通量高于旱季,然而湖泊水体却表现为旱季总磷浓度高于雨季,这表明华阳河湖群低水位季节总磷浓度波动主要受内源释放影响.湖泊内源污染控制要结合湖泊面积、水文、水体营养盐水平、沉积物内源释放特征等,科学制定综合治理措施.华阳河湖群内源控制措施又分两种,分别为拆除湖泊围栏网和控制湖泊生态水位.自2018年以来,宿松县拆除湖泊围栏网并禁止投饵,切断了每年大量外源性饵料的投放以及鱼类排泄物和残留物;同时,调整实施华阳河湖群生态水位管理办法,将10月—翌年2月水位普遍提升0.5 m,旨在抑制沉积物再悬浮,减少内源引起的磷浓度波动.经两项措施的实施,短期内缓解了枯水期华阳河湖群总磷浓度偏高的问题.然而,并不能夸大水位控制在内源污染控制中的作用.研究[24-25]表明,打破通江湖泊原有的水文规律维持高水位将破坏水生植物的生境,特别是春季大部分水生植被萌芽期需要低水位条件.受1998年洪水影响,长江中下游大部分浅水湖泊水生植物消亡,造成了湖泊营养盐含量升高[26-27].

对于华阳河湖群这样的浅水湖泊,实施湖泊生态恢复工程,恢复20年前的草型湖泊生态系统非常重要.事实证明,在沉水植物为主的草型湖泊生态系统,水体营养盐普遍较低,也不易发生水华,东太湖、梁子湖等以沉水植物为优势的草型生态系统富营养化指数明显较低[28-29].在浅水湖泊,水生植物在吸收水体与沉积物中营养盐的同时,还有削减风浪、固结沉积物、抑制沉积物再悬浮过程及释放营养盐的重要作用.大面积湖区水生植物种植的成功与否受水文规律、营养盐、风浪扰动、透明度、沉积物底质等影响.近年来,在太湖、梁子湖等湖泊实施了水生植被恢复工程,但效果不一,究其原因是水生植被生长、发育过程与相关限制性因素的关系不够明确[30].华阳河湖群有显明的通江湖泊水位变化规律,其水生植被恢复模式、试验区域等需要进一步评估.

外源污染控制始终是华阳河湖群区域生态环境管理的主要策略.第二次全国污染源普查数据显示,研究区域内农业源磷贡献率最大.在严格控制流域内农药、化肥等施用量同时,也需掌握流域尺度营养盐输移和生物地球化学循环过程,以及研究实施入湖河流及面源污染控制措施,特别是磷的颗粒态的特性.无论是自然湖滨缓冲带,还是人工湿地等,都可有效拦截颗粒态磷进入湖泊.华阳河湖群有广袤的湖滨带,充分利用湖滨带拦截雨季悬浮物入湖将会有较好的效果.

4 结论

a) 华阳河湖群沉积物内源磷静态年释放量为12.92 t,动态磷年释放量为1.129×104~1.684×104 t,静态释放量约占动态释放量的0.1%;再悬浮期间内源磷释放占水体总磷比例可达75%;沉积物内源磷静态释放对短期水体总磷影响相对较小.

b) 春季和冬季(枯水期)沉积物内源磷动态释放量占全年释放量的比例在80%左右,是造成湖泊总磷浓度季节性差异的主要原因,总磷浓度的空间差异性可能与入湖河流分布及湖底地形变化有关.

c) 严格控制入湖营养盐通量的同时,恢复草型湖泊生态系统以控制内源磷释放是华阳河湖群水环境质量改善的重要措施之一.

d) 该研究主要是对实测数据的解耦分析,但是未对区域风速、风向、地形、波流作用、沉积物再悬浮与磷释放机制做出定量分析,因此还有待开展进一步的深入研究.

参考文献
[1]
秦伯强, 朱广伟, 张路, 等. 大型浅水湖泊沉积物内源营养盐释放模式及其估算方法:以太湖为例[J]. 中国科学:地球科学, 2005, 35(S2): 33-44. (0)
[2]
秦伯强, 朱广伟. 长江中下游地区湖泊水和沉积物中营养盐的赋存, 循环及其交换特征[J]. 中国科学:地球科学, 2005, 35(S2): 1-10. (0)
[3]
秦伯强, 杨桂军, 马健荣, 等. 太湖蓝藻水华"暴发"的动态特征及其机制[J]. 科学通报, 2016, 61(7): 759-770.
QIN Boqiang, YANG Guijun, MA Jianrong, et al. Dynamics of variability and mechanism of harmful cyanobacteria bloom in Lake Taihu, China[J]. Chinese Science Bulletin, 2016, 61(7): 759-770. (0)
[4]
朱广伟, 秦伯强, 高光. 风浪扰动引起大型浅水湖泊内源磷暴发性释放的直接证据[J]. 科学通报, 2005, 50(1): 66. DOI:10.3321/j.issn:0023-074X.2005.01.013 (0)
[5]
WU Dan, HUA Zulin. The effect of vegetation on sediment resuspension and phosphorus release under hydrodynamic disturbance in shallow lakes[J]. Ecological Engineering, 2014, 69(1): 55-62. (0)
[6]
TANG Chunyan, LI Yiping, HE Cao, et al. Dynamic behavior of sediment resuspension and nutrients release in the shallow and wind-exposed meiliang bay of lake Taihu[J]. Science of the Total Environment, 2020, 708: 131-135. (0)
[7]
王苏民, 窦鸿身. 中国湖泊志[M]. 北京: 科学出版社, 1998. (0)
[8]
NIXDORF B, DENEKE R. Why very shallow lakes are more successful opposing reduced nutrient loads[J]. Hydrobiologia, 1997, 342(1): 269-284. (0)
[9]
羊向东, 沈吉, 董旭辉, 等. 长江中下游浅水湖泊历史时期营养态演化及湖泊生态响应[J]. 中国科学:地球科学, 2005, 35(S2): 45-54. (0)
[10]
吴奇丽.水位波动对华阳河湖群浮游植物群落结构的影响[D].合肥: 安徽大学, 2018: 31-36. http://cdmd.cnki.com.cn/Article/CDMD-10357-1018087775.htm (0)
[11]
曾凯.基于不同目标的典型闸控河湖生态需水计算研究[D].武汉: 长江科学院, 2019: 83-94. http://cdmd.cnki.com.cn/Article/CDMD-82305-1019849938.htm (0)
[12]
安徽省生态环境厅.安徽省生态环境状况公报: 2018[R].合肥: 安徽省生态环境厅, 2019. (0)
[13]
FICK A. Ueber diffusion[J]. Annalen der Physik, 1855, 170(1): 59-86. (0)
[14]
ULLMAN W J, SANDSTROM M W. Dissolved nutrient fluxes from the nearshore sediments of Bowling Green Bay, Central Great Barrier Reef Lagoon (Australia)[J]. Estuarine Coastal & Shelf Science, 1987, 24(3): 289-303. (0)
[15]
杜宏伟, 张恒军, 范中亚, 等. 华阳河湖群底泥沉积物特性研究[J]. 环境科学与技术, 2015, 38(4): 128-132.
DU Hongwei, ZHANG Hengjun, FAN Zhongya, et al. Study on characteristics of sediments in Huayang River-Lake Group[J]. Environmental Science & Technology (China), 2015, 38(4): 128-132. (0)
[16]
夏守先, 杨丽标, 张广萍, 等. 巢湖沉积物-水界面磷酸盐释放通量研究[J]. 农业环境科学学报, 2011, 30(2): 322-327.
XIA Shouxian, YANG Libiao, ZHANG Guangping, et al. Benthic flux of phosphate across sediment-water interface in Lake Chaohu[J]. Journal of Agro-Environment Science, 2011, 30(2): 322-327. (0)
[17]
罗潋葱, 秦伯强. 太湖波浪与湖流对沉积物再悬浮不同影响的研究[J]. 水文, 2003, 23(3): 1-4.
LUO Liancong, QIN Boqiang. Comparison between wave effects and current effects on sediment resuspension in Lake Taihu[J]. Hydrology, 2003, 23(3): 1-4. (0)
[18]
罗潋葱, 秦伯强, 胡维平, 等. 不同水动力扰动下太湖沉积物的悬浮特征[J]. 湖泊科学, 2004, 16(3): 273-276.
LUO Liancong, QIN Boqiang, HU Weiping, et al. Sediment re-suspension under different hydrodynamic disturbances in Lake Taihu[J]. Journal of Lake Sciences, 2004, 16(3): 273-276. DOI:10.3321/j.issn:1003-5427.2004.03.013 (0)
[19]
LI Yiping, TANG Chunyan, WANG Jianwei, et al. Effect of wave-current interactions on sediment resuspension in large shallow lake Taihu, China[J]. Environmental Science & Pollution Research, 2016, 24(4): 1-11. (0)
[20]
许遐祯, 陶蓉茵, 赵巧华, 等. 大型浅水湖泊太湖波浪特征及其对风场的敏感性分析[J]. 湖泊科学, 2013, 25(1): 55-64.
XU Xiazhen, TAO Rongyin, ZHAO Qiaohua, et al. Wave characteristics and sensitivity analysis of the wind field in a large shallow lake:Lake Taihu[J]. Journal of Lake Sciences, 2013, 25(1): 55-64. DOI:10.3969/j.issn.1003-5427.2013.01.008 (0)
[21]
胡维平, 胡春华, 张发兵, 等. 太湖北部风浪波高计算模式观测分析[J]. 湖泊科学, 2005, 17(1): 41-46.
HU Weiping, HU Chunhua, ZHANG Fabing, et al. On the emperical wind-wave height model in northern Lake Taihu, China[J]. Journal of Lake Sciences, 2005, 17(1): 41-46. DOI:10.3321/j.issn:1003-5427.2005.01.007 (0)
[22]
高小孟, 李一平, 杜薇, 等. 太湖梅梁湾湖区悬浮物动态沉降特征的野外观测[J]. 湖泊科学, 2017, 29(1): 52-58.
GAO Xiaomeng, LI Yiping, DU Wei, et al. Suspended solid settlement characteristics in Meiliang Bay of Lake Taihu:in situ observation[J]. Journal of Lake Sciences, 2017, 29(1): 52-58. (0)
[23]
CHUNG E G, BOMBARDELLI F A, SCHLADOW S G. Sediment resuspension in a shallow lake[J]. Water Resources Research, 2009. DOI:10.1029/2007WR006585 (0)
[24]
徐治国, 何岩, 闫百兴, 等. 营养物及水位变化对湿地植物的影响[J]. 生态学杂志, 2006, 25(1): 87-92.
XU Zhiguo, HE Yan, YAN Baixing, et al. Effects of nutrients and water level fluctuation on wetland plants[J]. Chinese Journal of Ecology, 2006, 25(1): 87-92. DOI:10.3321/j.issn:1000-4890.2006.01.017 (0)
[25]
刘永, 郭怀成, 周丰, 等. 湖泊水位变动对水生植被的影响机理及其调控方法[J]. 生态学报, 2006, 26(9): 3117-3126.
LIU Yong, GUO Huaicheng, ZHOU Feng, et al. Role of water level fluctuation on aquatic vegetation in lakes[J]. Acta Ecologica Sinica, 2006, 26(9): 3117-3126. DOI:10.3321/j.issn:1000-0933.2006.09.042 (0)
[26]
陈中义, 雷泽湘, 周进, 等. 梁子湖六种沉水植物种群数量和生物量周年动态[J]. 水生生物学报, 2000, 24(6): 582-588.
CHEN Zhongyi, LEI Zexiang, ZHOU Jin, et al. Monthly quantitative and biomass dynamics of six submerged macrophytes populations in Liangzi Lake[J]. Acta Hydrobiologica Sinica, 2000, 24(6): 582-588. DOI:10.3321/j.issn:1000-3207.2000.06.002 (0)
[27]
李伟, 刘贵华, 熊秉红, 等. 1998年特大洪水后鄱阳湖自然保护区主要湖泊水生植被的恢复[J]. 武汉植物学研究, 2004, 22(4): 301-306.
LI Wei, LIU Guihua, XIONG Binghong, et al. The restoration of aquatic vegetation in lakes of Poyang Lake nature reserve after catastrophic flooding in 1998[J]. Journal of Wuhan Botanical Research, 2004, 22(4): 301-306. DOI:10.3969/j.issn.2095-0837.2004.04.005 (0)
[28]
查慧铭, 朱梦圆, 朱广伟, 等. 太湖出入湖河道与湖体水质季节差异分析[J]. 环境科学, 2018, 39(3): 1102-1112.
ZHA Huiming, ZHU Mengyuan, ZHU Guangwei, et al. Seasonal difference in water quality between lake and inflow/outflow rivers of Lake Taihu, China[J]. Environmental Science, 2018, 39(3): 1102-1112. (0)
[29]
秦云, 李艳蔷, 吴丽秀, 等. 梁子湖水质时空格局分析[J]. 湖泊科学, 2016, 28(5): 994-1003.
QIN Yun, LI Yanqiang, WU Lixiu, et al. Analysis on the spatial-temporal patterns of water quality in Lake Liangzi[J]. Journal of Lake Sciences, 2016, 28(5): 994-1003. (0)
[30]
王志强, 崔爱花, 缪建群, 等. 淡水湖泊生态系统退化驱动因子及修复技术研究进展[J]. 生态学报, 2017, 37(18): 6253-6264.
WANG Zhiqiang, CUI Aihua, MIU Jianqun, et al. Research progress on the driving factors of freshwater lake ecosystem degradation and associated restoration techniques[J]. Acta Ecologica Sinica, 2017, 37(18): 6253-6264. (0)