磷是湖泊富营养化的重要营养元素^[1-4]，是营养状态评价的控制性指标^[5]。2007年太湖蓝藻暴发引发无锡供水危机后，太湖流域综合治理各项措施逐步落地^[6]，太湖及主要入湖河道水质呈改善趋势^[7]，但太湖总磷仍处于高位波动^[8]，驱动太湖蓝藻生物量扩张^[2]，湖体富营养状况未得到根本扭转^[9]。外源输入的总磷负荷是太湖总磷的主要来源^[10]，入湖总磷负荷逐年持续累积，远超水体自净能力^{[6, 11]}，是太湖总磷居高不下的重要原因。太湖地处长江中下游，区域地理磷本底高^[12]，同时太湖作为大型浅水湖泊，磷内源供给通量大、速度快、效率高^[13]，也对太湖总磷浓度具有不可忽视的影响^[14]。

近年来，在当前全球气候普遍变暖，内源污染和外源污染仍未得到有效遏制等多重因素综合影响下，太湖蓝藻密度和水华强度总体呈上升趋势，防控形势较为严峻。氮磷等营养物质的输入不断累积是导致湖泊水体富营养化的主要原因，湖泊富营养化引起了藻类和水生生物的过度繁殖，继而可能引发水质严重恶化和藻类水华等水生态水环境灾害^[15]。科学探讨太湖总磷浓度处于高位波动的原因，对管理者科学施策、净化太湖水质、抑制蓝藻水华的发生、改善太湖水体生态环境和保障安全供水具有重要意义。本研究应用2007年以来太湖水质监测资料和调查数据，并基于太湖总磷污染负荷质量平衡估算和滞留量分析、占比动态变化等，对影响太湖总磷输入输出等各要素进行量化研究，期望从总磷污染负荷动态平衡角度分析太湖总磷主要来源与总磷浓度高位波动的原因，为太湖水环境的治理和保护工作提供参考。

1 材料与方法 1.1 数据来源

本文采用的太湖及环太湖河流相关水质数据为太湖流域水文水资源监测中心2007年以来每月一次的长时序监测资料，其中总磷测定参照《水质总磷的测定钼酸铵分光光度法》(GB 11893—1989)，太湖和22条主要入湖河道监测站点分布见图 1。环湖巡测段水量数据来自太湖流域管理局组织的江苏、浙江水文部门开展的环太湖水文巡测资料整编成果；蓝藻和水草打捞量、水厂取水量数据摘自历年《江苏省水资源公报》，鱼类捕捞量数据来自江苏省太湖渔业管理委员会。其它数据摘自文献并在文中逐一标注说明。

1.2 方法 1.2.1 质量平衡计算

影响太湖总磷污染负荷质量平衡计算的因子众多，转化机理复杂。近年来，关于太湖磷污染负荷的研究较多，朱广伟等^[8]基于2005—2018年长时序监测资料，重点关注水体中不同形态磷赋存量的季度变化，估算出太湖水体总磷赋存量；王华等^[6]和吴浩云等^[16]分别分析了2010—2017、1980—2020年入湖河道输入、输出磷污染负荷，但并未考虑干湿沉降输入和水生动植物捕捞等输出影响；翟淑华等^[10]依据太湖总磷质量平衡，考虑出入湖各要素，计算了2015—2016年大水年太湖总磷的滞留量约分别占当年磷总负荷的21.5%和27.5%，但分析数据系列有限，结论有一定的局限性。

本研究参考翟淑华等^[17]的方法，对太湖总磷污染负荷主要影响因子进行分析，基于太湖水量平衡，建立质量平衡计算公式对太湖总磷污染负荷质量平衡进行估算分析。入湖污染负荷估算中，因太湖现已全面取缔入湖排污口，入湖河道和大气干湿沉降是太湖主要污染负荷来源，同时受监测条件、资料所限，其它方式出入湖总磷负荷暂不纳入估算；出湖污染负荷中，综合考虑了出湖河道、水厂取水、生物收获带出量等因子。

水量平衡计算公式为：

$ W_{\mathrm{I}}+P+D \pm \Delta V=W_{\mathrm{o}}+E+W_{\mathrm{c}}+\Delta E $

(1)

式中，W_I、W_o分别为依据逐年环太湖水文巡测整编资料统计的入湖、出湖水量，P为湖面降水量，以环太湖周边五站(西山、大浦口、夹浦、望亭立交、小梅口)平均降水量计算太湖湖面降雨深，以年太湖平均水位计算太湖水面面积，相乘得出；D为陆地产水入湖量，ΔV为太湖蓄变量，W_c为取水户直接取水量；E为湖面蒸发量。

质量平衡计算公式为：

$ W_{\text {入湖 }}+W_{\text {干湿沉降 }}+W_{\text {年初 }}=W_{\text {滞留量 }}+W_{\text {出湖 }}+W_{\text {取水 }}+W_{\text {水生生物捕捞 }}+W_{\text {年末 }} $

(2)

因此，

$W_{\text {滞留量 }}=W_{\text {入湖 }}+W_{\text {干湿沉降 }}+W_{\text {年初 }}-\left(W_{\text {取水 }}+W_{\text {水生生物捕捞 }}+W_{\text {出湖 }}+W_{\text {年末 }}\right) $

(3)

式中，W_入湖和W_出湖分别为环太湖河道带入和带出太湖的磷负荷；W_干湿沉降为通过干沉降和湿沉降带入太湖的磷负荷；W_年初和W_年末分别为年初和年末太湖水体中磷存量；W_取水为通过水厂取水带出磷负荷；W_{水生生物捕捞}为通过捕捞带出磷负荷，包括蓝藻、水草打捞和鱼类捕捞等；W_滞留量为太湖总磷负荷滞留量，主要为滞留于底泥和被水生生物吸收转化的总磷负荷量。

1.2.2 总磷质量平衡各要素计算

1) 水质浓度: 水质浓度根据太湖33个监测站点代表的水面面积加权平均计算得出，监测站点代表面积按照泰森多边形法确定。入湖河流平均浓度为《太湖流域管理条例》中明确的22条主要入湖河道监测站点的数学平均浓度。

2) 出入湖污染负荷: 环太湖出入湖水量资料是以水文巡测段或单站为最小单元进行整编，为了计算匹配，出入湖污染负荷以水文巡测段或单站为单元进行计算统计。月度环太湖出入湖磷污染负荷计算采用水质监测断面的水质浓度与其巡测段或单站月平均流量乘积计算得出，年度出入湖磷负荷总量为所有巡测段(站)年度磷负荷相加的和。

3) 干湿沉降入湖污染负荷: 杨龙元等^[18]、刘涛等^[15]、翟淑华等^[10]、许志波等^[19]分别研究了2002年7月-2003年6月、2011年、2015—2016年和2018年的太湖总磷干湿沉降成果，依据这些研究成果中的太湖干、湿沉降率并结合太湖水雨情分别估算其相近年份的2007—2010、2012—2014、2017、2019—2020年总磷干湿沉降负荷，估算公式为：

$ Q=Q_{\mathrm{d}}+Q_{\mathrm{w}}=F_{\mathrm{d}} \times d \times S+F_{\mathrm{w}} \times S $

(4)

其中，

$F_{\mathrm{w}}=k_{\mathrm{w}} \times C \times R $

(5)

式中，Q为沉降量(kg)；Q_d、Q_W分别为干、湿沉降量(kg)；F_d、F_w为干、湿沉降率(kg/km²)；d为一年干沉降天数(d)，S为太湖水面积(km²)；k_w为换算系数，为1 L ·kg/(mg ·mm ·km² ·365 d), C为降雨总磷质量浓度(mg/L)，R为全年太湖湖面降水量(mm)。

4) 水生生物捕捞磷负荷: 水生生物捕捞输出磷负荷为捕捞量与总磷含量的乘积。其中，打捞的蓝藻藻水中含藻率约为0.5%，蓝藻体内总磷含量约为0.68%，蓝藻总磷含量为蓝藻藻水的0.0034%；根据相关研究，大型水生植物体中总磷含量约为0.044%(湿重)^[20-22]；鱼类体内总磷含量约为0.31%^[23]。

5) 取水输出磷负荷: 水厂直接取水带出的磷负荷量按照当年的取水量和原水水质浓度估算。水厂取水水质采用与水厂取水口位置邻近的沙墩港、渔业村、贡湖、胥口、漫山和庙港6个监测点浓度的算术平均值。

6) 年初原有量和年末残留量: 太湖水体中污染负荷年初原有量和年末残留量分别考虑太湖蓄水量、分区水面积和分区水质，其中蓄水量和分区水面积按照年初年末水位条件下水位-库容曲线和水位-面积曲线计算，计算公式为：

$ W=Q \times \sum\limits_{i=1}^9 \frac{S_i}{S} \times C_i $

(6)

式中，W为年初(年末)太湖水体中原有(残留)量；Q为年初(年末)太湖蓄水量，S_i为太湖分区水面积，S为太湖水面积；C_i为太湖第i分区年初(年末)总磷浓度。考虑到监测数据连续性影响，为保证系列成果衔接，年初原有量与上年度年末残留量均计为此两项的均值。

2 结果 2.1 太湖总磷污染负荷各项影响因子输入输出量估算 2.1.1 环湖河道出入湖总磷污染负荷

2007年以来入湖河道入湖总磷污染负荷量为1835~2799 t，未见明显下降趋势，其中以2010年最高，2016和2011年其次，近4年基本维持在2000 t左右。出湖总磷污染负荷量为516~906 t，其中以2016年最高，其次是2020年。净入湖总磷污染负荷量为1090~2115 t，以2010年最高，之后有所下降，2019年最低，近4年基本维持在1200 t左右(图 2)。

2007年以来，22条主要入湖河道总磷平均浓度由0.190 mg/L下降到2020年的0.128 mg/L(图 3), 但输入总磷污染负荷却未见明显下降趋势，总磷入湖污染负荷增加与入湖水量增加有关^[16]。从长时序水文资料看，2007年后太湖流域处于相对丰水期，季海萍等^[24]对太湖流域年降水量变化采用Mann-Kendall趋势检验法分析，发现2008—2017年较1986—2007年系列平均降水量多150 mm，增加了13%。2007—2020年环太湖年均入湖水量约为114.0亿m³，较1986—2006年均值增加了33.0亿m³，增幅近41%。

2.1.2 大气干湿沉降输入

以往分析表明，太湖外源输入总磷负荷除环太湖河道的径流输入外，大气干湿沉降也是太湖磷等营养盐输入的重要途径^[25]。工业生产中化工染料的大量使用、农业化学肥料的大量施用和畜牧业的快速发展均向大气中排放了大量的氮、磷污染物，这些污染物可通过干沉降(空气尘粒)和湿沉降(雨、雪、雾等)的方式进入太湖，增加太湖污染物总量^{[15, 19]}。根据估算结果，2007—2017年，太湖干湿沉降总磷污染负荷量呈先上升后下降的变化趋势，均值约为935 t，并以干沉降为主^[15-19]，其中干沉降总磷污染负荷量均值为660 t，湿沉降均值为275 t。2018年6月，国务院实施“打赢蓝天保卫战三年行动计划”，推动产业、能源等重点领域结构优化，大气污染物排放总量大幅减少，在一定程度上有效降低了入太湖干湿沉降负荷，2018年起太湖干湿沉降总磷污染负荷量大幅下降，2018—2020年干湿沉降总磷污染负荷量在300~350 t之间，其中干沉降均值为135 t，湿沉降均值为200 t(附表Ⅰ)。

2020年6—7月，太湖流域水文水资源监测中心曾对太湖湖面降水进行5次抽样监测，结果显示仅1次总磷浓度超过0.01 mg/L的检出限，为0.012 mg/L，浓度远低于同期太湖总磷浓度(0.068~0.077 mg/L)。若以0.012 mg/L作为当年降雨总磷浓度的参考，结合当年水雨情估算湿沉降总磷负荷输入量仅约为42.9 t；本文研究时因缺乏系统、完善的干沉降量资料和湿沉降总磷浓度资料，大气干湿沉降负荷量结果为估算值，下一步有待开展深入专题研究。

2.1.3 水生动植物的捕捞输出

水生动植物均能吸收转化水中的氮、磷污染物，水草打捞、藻类打捞和鱼类捕捞等能够将氮、磷等带出水体，是湖泊水体总磷污染物总量的主要影响因子之一。

2007年无锡市供水危机以来，为应对蓝藻暴发，无锡市、苏州市和常州市均开展了蓝藻打捞工作，蓝藻藻水打捞量持续上升。2012年以前，太湖蓝藻打捞量在100万t以下，蓝藻打捞带出总磷污染负荷量约6.8~33.0 t；2012年起，太湖蓝藻打捞量维持在120万t以上，2015年起，上升至150万t以上，近3年均超过170万t，每年通过蓝藻打捞从太湖带出总磷污染负荷量超过57.8 t^[26]。

大型水生植物可以通过根系吸收沉积物中的生物有机磷，同时能够抑制底泥再悬浮，有效减少沉积物间隙水中的磷向上覆水的扩散量，还可通过茎叶拦截、吸附水中的颗粒物质来吸附水中的可溶性磷^[27]。但其死亡腐烂会造成氮、磷等营养物质的再次释放，污染水质。2012年起，江苏省对贡湖、东太湖等部分湖区实施区域水草打捞。近年，太湖水草打捞量总体呈先上升后下降的趋势，以2014年最高，当年水草打捞量(湿重)为20万t，估算从太湖带出总磷污染负荷量为88 t；2016年后，为养护太湖水草，环湖各地市逐年减少了收割打捞量，2017—2020年均打捞量约为6万~7万t，带出总磷污染负荷量约为30 t。

2007—2014年太湖捕捞量参考谷孝鸿等的报道^[28]，之后采用江苏省太湖渔业管理委员会的统计数据。2007—2019年，太湖鱼类捕捞量在3.5万~7.2万t左右，总体呈上升趋势；鱼类捕捞带出太湖的总磷负荷量约为108~226 t，均值为169 t。2020年江苏省农业厅发布太湖“十年禁捕”通告，当年10月1日起全面停止捕捞作业，估算2020年捕捞量为5.4万t，带出总磷量约为167.4 t。

2007—2020年水生动植物捕捞量及其带出太湖的总磷污染负荷量详见附表Ⅱ。太湖水生动植物捕捞带出总磷量呈阶段式上升趋势，其中2007—2011年由100 t逐渐上升至200 t左右；水草打捞开始后，2012年起增加至200 t以上；2015年起均超过了300 t，并以2016年为最高，之后水草打捞减少、鱼类禁捕后有所降低。

2.1.4 取水输出

2007年以来，太湖各水厂取水量有所上升，取水总磷平均浓度总体呈下降趋势。水厂取水输出总磷污染负荷均值为60 t左右，年际间呈现不同的变化，近5年均值在65 t左右。2007年以来，太湖水厂输出总磷污染负荷量详见附表Ⅲ。

2.1.5 太湖水体总磷负荷年初原有量和年末残留量估算

由2007—2020年太湖水体中总磷污染负荷量变化(附表Ⅳ)可见，近年太湖水体中总磷污染负荷约为300~500 t，年内水体中总磷污染负荷量变化有增有减，除2020年年末较年初减少168 t外，其余年份基本在150 t之内变化。太湖水体中总磷污染负荷量除与水体中总磷浓度有关外，也受太湖蓄水量直接影响。

2.2 太湖总磷污染负荷质量平衡计算结果

近年太湖总磷污染负荷输入量约为2155~3771 t，均值为2889 t；输出量基本在741~1309 t之间，均值为940 t；总磷输入量远大于输出量，多出的总磷污染负荷量一部分为太湖水体中总磷存量变化，其余为滞留于底泥中或水生生物体内的总磷负荷量^[10]，质量平衡计算结果详见表 1。其中，由于干湿沉降量缺少系列实测数据，部分年份采用文献数据结合水情相近年份进行估算，对质量平衡计算结果精度有一定影响；因未查询到相关数据，根据历年数据对2012年水草打捞量和2020年生物收获量进行了估算，估算值与真实值有一定误差，但水草打捞和鱼类捕捞均在出湖污染负荷总量中占比较小，对当年质量平衡计算结果的影响有限。

根据以上计算结果，2007—2020年太湖总磷年负荷滞留量为1080~3036 t，年均值为1959 t，表明目前太湖底泥仍是磷“汇”，对总磷的吸附和沉积作用要远远大于释放。总磷滞留量约占太湖总磷污染负荷量的41%~74%，总体呈先上升后下降的变化趋势，以2011年最高，2012年之后基本以2000 t以下为主，2018年开始年负荷滞留量明显减少，连年处于最低的水平。滞留于太湖中的磷主要吸附于底泥中或被水生生物吸收，大量总磷负荷逐年累积，成为影响太湖水体总磷浓度的重要内源，这也是近年来太湖入湖河流总磷浓度持续改善，但太湖总磷浓度持续高位波动的主要原因之一。

3 分析与讨论 3.1 太湖输入、输出总磷负荷分析

以环湖河道入湖总磷污染负荷量、干湿沉降入湖量和年初水体原有量合计估算太湖总磷污染负荷量。2007年以来，太湖总磷污染负荷量总体呈先上升后下降的变化趋势(图 4a)，以2011年最高(4119 t)，2017年以后基本维持在2600~3000 t。入太湖河道输入总磷污染负荷仍是太湖总磷的最主要来源，其占比达到55%~73%，2018年起基本维持在70%左右。干湿沉降输入是太湖总磷污染负荷输入的第二大来源，其输入总磷负荷占12%~38%，2018年起基本维持在12%左右。年初太湖水体中原有总磷负荷量占8%~15%，2015年之后基本维持在11%~15%之间。

太湖各途径输出总磷污染负荷量分布(图 4b)显示，水厂取水输出总磷污染负荷量仅占太湖负荷量的2%~3%左右；水生动植物捕捞约占4%~12%，近年基本维持在10%左右，2020年有所下降；环湖河道出湖带出总磷污染负荷是太湖总磷输出的主要途径，其占比总体呈上升趋势，2015年起基本维持在13%~30%之间，以2020年最高为29.8%；年末水体残留量基本以8%~15%为主，以2020年最低(8.6%)。滞留在底泥中和被吸收、利用分解的总磷污染负荷量约占太湖总磷负荷量的41%~74%，总体呈先升后下降的趋势，2012年之后逐步降低，2015年起维持在41%~55%左右。

3.2 滞留污染负荷对总磷浓度变化的影响

《太湖流域水环境综合治理总体方案》中明确太湖总磷的纳污能力为514 t/a，近年太湖总磷净输入污染负荷均值为1959 t，约为太湖总磷纳污能力的3.81倍，大量总磷污染负荷以太湖底泥的形式存在于湖体中^[10]。太湖流域管理局2013年和2018年的太湖底泥调查表明，太湖底泥中总磷含量明显上升，全湖平均由2003年的598.1 mg/kg上升为674.8 mg/kg^[29]，沉积的磷在一定条件下会逐步释放，成为影响水质的主导因子^[30]。底泥是磷元素累积和再生的重要场所，沉积在底泥中的磷易于泛起再悬浮，通过各种复杂的过程再次释放到上覆水中，造成水体中总磷浓度的升高^[31]。据逄勇等报道^[32]，每年太湖内源年均进入水体的净底泥量为19.03万t，总磷约275.5 t。

已有大量研究表明，底泥中磷的释放主要与水体中的溶解氧、水温、pH和水体扰动有关，溶解氧浓度的降低、水温的升高、pH值的上升以及湖流速度的加大都有利于底泥中磷的释放^{[29-30, 33-35]}。2007年后，太湖溶解氧浓度年均值基本在9.0~10.0 mg/L之间波动，pH年均值在8.11~8.52之间波动，均无明显的趋势性变化。吴凡等的研究表明，近年来太湖地区的气温有明显升高趋势，增幅超过了0.2℃/10 a^[36-37]，太湖流域水文水资源监测中心对太湖水温监测资料也表明，2007—2013年太湖平均水温较2014—2020年上升了0.9℃(图 5a)，而冬春季(12月-次年4月)太湖平均水温上升了2.0℃(图 5b)。一方面，太湖水温升高促进底泥中磷的交换，有利于底泥中磷的释放；另一方面，水温是影响太湖蓝藻水华发生的主要气象因子之一，尤其是冬春季水温偏高，有利于蓝藻越冬、复苏。蓝藻生长的“泵吸作用”以及水华堆积之后的“耗氧效应”，会加快底泥磷的释放，底泥中释放的磷除了满足藻类生产需求外，还能通过生物富集，使水体中总磷浓度升高^[38]。

太湖作为浅水型湖泊，不同于深水湖泊靠浓度梯度释放营养物质^[39]，风浪等动力扰动引起沉积物再悬浮对浅水湖泊内源释放影响较大。水流流速加大也会增加水体扰动，使沉积在底泥中的磷易于扩散到上覆水中，提高水体中的总磷浓度。2007年以来，太湖流域处于相对丰水期，加之“引江济太”工程的实施，促进了太湖水体的循环，环太湖出入湖水量均呈上升趋势^[40]，太湖换水周期明显缩短，太湖及其周边水体的流动速度加快，一定程度上促进了底泥间隙水中的营养物质的加速释放，同时水动力作用使脱离底泥的颗粒在数量上大过絮凝沉淀的颗粒，可导致总磷浓度升高。

4 结论

2007年以来太湖总磷污染负荷的年均输入量为2889 t，年均输出量为940 t，近年来年输入量已有明显下降，但远大于输出量的状况仍在持续；年均约有1959 t的总磷滞留在湖体的底泥中或水生生物体内，约为太湖总磷纳污能力的3.81倍，是影响太湖水体总磷浓度的重要内源。2007年以来，总磷输入污染负荷年际间总体为先升后降，2018年后基本维持在2200 t左右。入湖河道总磷污染负荷是太湖总磷最主要的输入来源，占太湖总磷负荷的55%~73%，仍未见明显下降趋势；湖面干湿沉降也是太湖总磷输入的主要来源之一。2007年以来，太湖总磷出湖污染负荷年均值约为940 t，总体呈上升趋势，仅占太湖总磷输入污染负荷的1/3左右。出湖河道带出总磷污染负荷是太湖总磷输出最主要的途径，占太湖总磷负荷量的13%~30%；年际间通过水草和蓝藻打捞、鱼类捕捞、水厂取水等带出的总磷负荷量较为稳定，也是太湖总磷污染负荷输出的重要途径。近年来，太湖地区气温升高、太湖水体流动速度加快一定程度上促进了底泥中总磷的释放，对太湖水体中总磷浓度升高有一定的影响。

5 附录

附表Ⅰ~Ⅳ见电子版(DOI: 10.18307/2023.0520)。