翟淑華，周 婭**，程媛華，蔡 杰，胡燕雯

(1：太湖流域水資源保護局，上海 200434) (2：香港理工大學(xué)，香港 999077)

2007年太湖藍藻暴發(fā)引發(fā)無錫供水危機后，太湖流域水環(huán)境綜合治理的力度進一步加大，連續(xù)十一年實現(xiàn)太湖“確保飲用水安全、確保太湖水體不發(fā)生大面積水質(zhì)黑臭”的目標. 太湖氨氮濃度和總氮(TN)濃度明顯下降，但總磷(TP)濃度于2015-2016年重回升勢，相應(yīng)2017年太湖出現(xiàn)歷史上最大面積的藍藻水華，太湖藍藻暴發(fā)勢頭未得到有效控制，富營養(yǎng)化形勢依然嚴峻[1-5]. 太湖TP濃度波動上升和藍藻大面積暴發(fā)的原因，再次成為政府各級部門和社會關(guān)注的焦點，但尚未達成統(tǒng)一認識. Qin等[2]的研究認為近年來TP濃度波動上升的原因是多方面的，是污水排放標準較低、面源控制力度不足、內(nèi)源措施不夠經(jīng)濟有效、氣候變化等多因素共同作用的結(jié)果. 張民等[3]認為，對藍藻水華生物量的相關(guān)驅(qū)動因素中，營養(yǎng)鹽(TP或氮磷比)可單獨解釋藍藻生物量變化的29%. 朱廣偉等[4]重點對太湖北部湖區(qū)開展了研究，認為藻類生物體含有的營養(yǎng)鹽是水體營養(yǎng)鹽的重要組成部分，水華堆積之后的耗氧效應(yīng)大大加快了底泥磷的釋放，造成夏季水體磷濃度升高，同時近年來北部湖區(qū)葉綠素a濃度的波動很大程度上受水文氣象因子的影響. 王華等[5]基于2010-2017年的監(jiān)測數(shù)據(jù)及其相關(guān)性分析，認為太湖2014-2017年TP濃度總體呈上升趨勢，尤其是2016年為8年間最高值，并從長系列趨勢分析的角度分析了入湖河流污染負荷量、水生植被和藍藻變化與TP濃度的關(guān)系，但并未對2015-2016年大洪水條件對TP的影響進行深入分析.

太湖流域2015-2016年發(fā)生流域性大洪水，特別是2016年太湖最高水位達4.87 m，僅低于1999年的4.97 m的歷史最高水位0.10 m，太湖TP濃度在高水位條件下不降反升. 朱偉等[6]認為，太湖梅梁湖等北部湖灣2015-2017年TP濃度相較于2012-2014年提高了20%～25%左右，為藍藻水華形成提供了豐富的物質(zhì)基礎(chǔ)，并認為TP濃度波動的主要原因為：大洪水攜帶大量營養(yǎng)鹽進入太湖，由于磷大多數(shù)以顆粒態(tài)存在，洪水消退后仍沉積在湖泊中，并在適宜條件下釋放進入水體. 但該研究僅通過2015-2016年入湖TP負荷總量偏大進行推斷，詳細論證主要針對北部湖灣. 因此，有必要進一步對2015-2016年進出全太湖的TP來源開展更為全面、深入的定量計算和詳細論證. 研究表明，環(huán)太湖河流是影響太湖水環(huán)境的重要因素，其進出太湖的水量、水質(zhì)和污染負荷量一直以來得到廣泛關(guān)注，有對長系列數(shù)據(jù)的變化趨勢分析[7]，也有對個別年份的深入探討[8-10]，但針對流域性大洪水條件下環(huán)太湖河道進出太湖的污染負荷量計算分析及其對太湖TP變化影響，以太湖水量和質(zhì)量平衡為基礎(chǔ)開展太湖TP來源與影響分析的研究還未見報道.

經(jīng)過太湖流域水環(huán)境綜合治理，入湖河道TP濃度呈改善趨勢，但太湖TP濃度卻于2015-2016年出現(xiàn)較為明顯的上升. 因此，鑒于2015-2016年水文條件的特殊性，不僅要考慮水質(zhì)濃度的變化，還需要綜合考慮TP入湖負荷量的變化. 本文從2015-2016年環(huán)太湖河道的進出太湖水量、TP負荷量計算入手，結(jié)合雨情、水情、太湖調(diào)蓄以及人為影響等各方面因素，分別開展出入太湖的水量和TP負荷質(zhì)量的平衡分析. 在此基礎(chǔ)上，開展太湖TP的外源、內(nèi)源變化趨勢及來源分析，探討2015-2016年太湖TP濃度升高的原因，提出有關(guān)今后太湖控磷措施的重點方向.

1 方法與計算

自1998年起，太湖流域水資源保護局連續(xù)每年組織開展環(huán)太湖河道進出太湖的污染負荷量計算. 為保持數(shù)據(jù)的連續(xù)性，本文采用的水量資料為整編后的環(huán)太湖水文巡測資料，水質(zhì)資料為流域機構(gòu)監(jiān)測的環(huán)湖河流及水功能區(qū)水質(zhì)監(jiān)測資料. 計算望亭(立交)站入湖污染負荷時，采用“引江濟太”每日監(jiān)測的水質(zhì)數(shù)據(jù)；計算其他河流污染負荷時，采用當月環(huán)湖河流水質(zhì)數(shù)據(jù)；計算出湖污染負荷時，采用太湖沿岸區(qū)相應(yīng)監(jiān)測點的水質(zhì)數(shù)據(jù). 計算采用的環(huán)太湖水文巡測段、水質(zhì)監(jiān)測站點分布及對應(yīng)關(guān)系詳見圖1.

圖1 環(huán)太湖水文巡測段、水質(zhì)監(jiān)測站點分布及對應(yīng)關(guān)系示意圖Fig.1 Sketch of the hydrological survey sections and water quality monitoring stations and their corresponding relationship

1.1 太湖水量平衡要素與計算方法

1.1.1 水量平衡計算方法 采用翟淑華等[11]的方法，基于2015-2016年環(huán)太湖河道進出湖水量計算，結(jié)合降雨量、蒸發(fā)量、太湖調(diào)蓄水量、取用水等要素，開展進出太湖的水量平衡計算分析，公式為：

WI-WO+P-E-WC+D-ΔV=ΔE1

(1)

式中，WI、WO為環(huán)太湖河流進、出湖水量，根據(jù)環(huán)太湖水文巡測整編資料成果統(tǒng)計；P為湖面降水量，即太湖年降雨量與太湖平均水位下對應(yīng)的水面面積的乘積；E為湖面蒸發(fā)量，即太湖年蒸發(fā)量與太湖平均水位下對應(yīng)的水面面積的乘積；D為陸地產(chǎn)水入湖量，即環(huán)太湖水文巡測線范圍內(nèi)的陸地區(qū)域的產(chǎn)水入湖量；WC為取水戶直接取水量；ΔV為太湖蓄變量，太湖年末與年初蓄水量的差值；ΔE1為水量平衡計算的絕對誤差.

水量平衡計算的相對誤差ΔE2的計算公式為：

ΔE2=2ΔE1/(WI+WO+P+E+WC+D+|ΔV|)

(2)

1.1.2 水量平衡要素計算 1)環(huán)太湖河道進出湖水量計算. 根據(jù)實測水量整編數(shù)據(jù)計算得到，2015、2016年環(huán)太湖河道入湖水量分別為119.08億和159.85億m3，出湖水量分別為118.47億和167.27億m3.

2)降雨量、蒸發(fā)量計算. 2015年太湖全年平均水位3.42 m，對應(yīng)的太湖水面面積為2338.6 km2，太湖湖面降水量1523.3 mm，折算為入湖量35.62億m3，湖面蒸發(fā)量為914.1 mm，折算為出湖量20.74億m3；2016年太湖全年平均水位3.58 m，對應(yīng)的太湖水面面積為2351 km2，太湖湖面降水量1861.1 mm，折算為入湖量43.75億m3，湖面蒸發(fā)量為1061.8 mm，折算為出湖量24.96億m3.

3)太湖調(diào)蓄水量計算. 2015年太湖年初水位為3.18 m，年末水位為3.50 m，根據(jù)太湖水位-庫容關(guān)系得到太湖蓄水量年初為49.00億m3，年末為56.52億m3，蓄變量為7.52億m3. 2016年太湖年初水位為3.41 m，年末水位為3.28 m，根據(jù)太湖水位-庫容關(guān)系得到太湖蓄水量年初為54.37億m3，年末為51.28億m3，蓄變量為-3.09億m3.

4)其他水量要素計算. 根據(jù)太湖流域管理局發(fā)布的水資源公報，2015年直接從太湖取水水量約為10.02億m3，2016年直接從太湖取水水量約為10.86億m3. 經(jīng)估算，環(huán)太湖巡測線所圍面積為3128.0 km2，扣除農(nóng)業(yè)灌溉耗水量后，2015年陸地產(chǎn)水入湖量大致為1.82億m3，2016年陸地產(chǎn)水入湖量大致為1.89億m3.

1.1.3 水量平衡計算結(jié)果 計算結(jié)果見表1，可見水量平衡誤差在規(guī)范要求的誤差范圍內(nèi).

表1 2015-2016年太湖水量平衡(單位：億m3)

1.2 太湖總磷質(zhì)量平衡計算

1.2.1 質(zhì)量平衡計算方法 采用翟淑華等[11]的方法，基于2015-2016年環(huán)太湖河道進出湖TP負荷量計算，結(jié)合大氣干濕沉降、取用水、水生生物收獲、太湖年初和年末存量等要素，開展進出太湖的TP質(zhì)量平衡計算分析，公式為：

W河入+W降雨+W降塵+W污染源+W年初=W河出+W取水+W收獲+W年末+W滯留量

(3)

式中，W河入為通過河流直接進入太湖的量，W降雨、W降塵為通過降雨、降塵直接進入太湖的量，W污染源為巡測線內(nèi)各污染源進入太湖的量；W年初為年初太湖水體中存量；W河出為通過河流直接出太湖的量；W取水為通過水廠取水帶出太湖的量；W收獲為通過水生生物收獲帶出太湖的量，包括藍藻打撈、水草收割、魚類捕獲等；W年末為年末太湖水體中存量；W滯留量為太湖扣除水體存量變化以后滯留在湖體中的量.

太湖地區(qū)污染源排放的污染負荷量主要通過環(huán)太湖河道進入太湖，暫未考慮巡測線內(nèi)污染源直接排放入湖量. 太湖物料平衡計算公式簡化后得到滯留的污染負荷量為：

W滯留量=(W河入+W降雨+W降塵)-(W河出+W取水+W收獲)-(W年末-W年初)

(4)

質(zhì)量平衡計算得到的平衡差主要為滯留量及資料統(tǒng)計帶來的誤差，當假設(shè)統(tǒng)計誤差為零時，則平衡差即認為是當年度太湖滯留量.

1.2.2 質(zhì)量平衡要素計算 1)環(huán)太湖河道出入湖TP負荷量計算. 環(huán)太湖水文巡測段水量資料有每天一次的連續(xù)監(jiān)測數(shù)據(jù)，而環(huán)太湖水質(zhì)監(jiān)測資料多為每月一次. 針對水量水質(zhì)監(jiān)測頻次不匹配的問題，時段通量的估算主要有兩種方法：第1種是將每月一次水質(zhì)監(jiān)測的當天負荷量作為當月的平均負荷量，適用于污染物量受流量影響較小的情況；第2種是將每月一次的水質(zhì)數(shù)據(jù)作為當月的平均水質(zhì)濃度，適用于污染物量受流量影響較大的情況[12]. 根據(jù)環(huán)太湖河道進出太湖的逐日流量統(tǒng)計，河流的日流量變化較大，水量影響明顯大于水質(zhì)的影響，故選用第2種方法進行環(huán)太湖出入湖污染負荷量的計算. 采用與水量巡測段相匹配的河流水質(zhì)資料，計算得到2015和2016年通過環(huán)太湖河流帶入TP的負荷量分別為2209和2594 t，帶出分別為643和906 t，通過河流凈入湖的TP負荷量為1565和1688 t.

2)大氣干濕沉降入湖量估算. 研究表明，氮、磷大氣干濕沉降是太湖營養(yǎng)鹽輸入的重要來源之一[13-17]. 本文綜合考慮監(jiān)測時間、點位、頻次和方法，在各研究結(jié)果中選取適宜的沉降率參數(shù)，其中，濕沉降率采用太湖流域管理局《太湖大氣干濕沉降監(jiān)測試驗分析單元工程成果報告》中的數(shù)據(jù)，干沉降率采用劉濤等[16]的研究結(jié)果. 結(jié)合流域降水情況，計算得到2015年太湖TP濕沉降量和干沉降量分別為182 和562 t，2016年太湖TP濕沉降量和干沉降量分別為222 和390 t.

3)其他途徑出湖量估算. 水廠直接取水帶出的TP量根據(jù)當年的取水量和原水水質(zhì)濃度估算，2015和2016年分別為65和65 t.W收獲分為W藍藻、W水草和W魚類，分別為藍藻打撈、水草收割及魚類捕撈帶出量. 藍藻打撈帶出的TP量根據(jù)江蘇省水資源公報中的藍藻打撈量和太湖健康狀況報告中藍藻的氮磷含量估算，2015和2016年分別為164和160 t. 水草收割帶出的TP量根據(jù)江蘇省水資源公報中的水草打撈量和王強等[18]所測水草中氮磷含量估算，2015和2016年分別為933和714 t. 漁業(yè)統(tǒng)計資料顯示，太湖魚類捕撈產(chǎn)量總體呈不斷增長趨勢[19]. 魚類捕撈帶出TP量根據(jù)江蘇省太湖漁業(yè)管理委員會統(tǒng)計的魚類捕撈量和太湖污染物自凈能力與適宜換水周期研究項目所測魚體中氮磷含量估算，2015和2016年分別為183和226 t.

4)太湖水體年初和年末存量計算. 根據(jù)太湖年初、年末蓄水量和年初、年末水質(zhì)進行計算，其中蓄水量按照年初、年末水位條件下的水位-庫容曲線插值計算. 計算得到2015年年初存量為353 t，年末存量為607 t；2016年年初存量為288 t，年末存量為462 t.

1.2.3 質(zhì)量平衡計算結(jié)果 依據(jù)上述方法的數(shù)據(jù)，開展2015-2016年太湖TP質(zhì)量平衡計算，計算結(jié)果見表2，2015和2016年太湖TP的滯留量分別為711和961 t.

表2 2015-2016年進出太湖的總磷負荷量

2 分析與討論

2.1 2015-2016年太湖總磷來源分析

一般情況下，輸入湖泊系統(tǒng)的磷大多來自河流，這部分磷主要是人類活動產(chǎn)生的污水帶來的，其他輸入途徑還有來自大氣的干濕沉降. 進入水體的TP在經(jīng)過河流帶出、取水帶出、水生生物收獲帶出后，剩余的TP滯留于湖體，一部分造成水體TP存量的變化，另一部分滯留于底泥中或水生生物體內(nèi). 2015和2016年太湖各途徑進出湖TP負荷量分布見圖2.

2015年，太湖水體年初存量為353 t，外源輸入量為2953 t，合計總負荷量為3306 t；經(jīng)各途徑帶出1988 t后，剩余TP量為1318 t，其中607 t為水體中年末存量，其余711 t即為滯留量. 在各帶入途徑中，通過河道帶入的TP量為2209 t，占總負荷量的66.8%；通過大氣干、濕沉降帶入的TP量分別為562和182 t，占比分別為17.0%和5.5%. 在各帶出途徑中，通過河道帶出的TP量為643 t，占總負荷量的19.4%；通過打撈藍藻、水草和魚類等水生生物帶出的TP量為1280 t，占比38.7%；通過取水帶出的TP量為65 t，占比2.0%. 相應(yīng)計算得出，太湖水體中增加的TP為254 t，而滯留于底泥和被水生生物吸收轉(zhuǎn)化的量則達到711 t，占總負荷量的21.5%(圖2a).

2016年，太湖水體年初存量為288 t，外源輸入量為3206 t，合計總負荷量為3494 t；經(jīng)各途徑帶出2071 t后，剩余TP量為1423 t，其中462 t為水體中年末存量，其余961 t即為滯留量. 在各帶入途徑中，通過河道帶入的TP量為2594 t，占總負荷量的74.2%；通過大氣干、濕沉降帶入的TP量分別為390和222 t，占比分別為11.2%和6.4%. 在各帶出途徑中，通過河道帶出的TP量為906 t，占總負荷量的25.9%；通過打撈藍藻、水草和魚類等水生生物帶出的TP量為1100 t，占比31.5%；通過取水帶出的TP量為65 t，占比1.9%. 相應(yīng)計算得出，太湖水體中增加的TP為174 t，而滯留于底泥和被水生生物吸收轉(zhuǎn)化的量則達到961 t，占總負荷量的27.5%(圖2b).

圖2 2015年(a)和2016年(b)各途徑進出太湖的總磷負荷量分布Fig.2 Distribution of total phosphorus flux in and out of Lake Taihu in different ways in 2015 (a) and 2016 (b)

綜合上述分析認為，影響太湖TP濃度變化的主要因素包括：1)環(huán)太湖河道是TP進入太湖最主要的途徑，尤其是2016年的特大洪水，大量營養(yǎng)鹽通過河道帶入太湖，占比74.2%；2)水生生物收獲對太湖水生態(tài)環(huán)境具有重要影響，2015年受太湖水草大面積收割及太湖持續(xù)高水位的影響，水生生物帶出TP量占比為38.7%，可能導(dǎo)致水生生物對磷的吸收利用減少. 3)逐年累積的內(nèi)源污染是太湖TP濃度升高的又一重要原因，底泥中的滯留量和生物中吸收轉(zhuǎn)化的TP量占比為21.5%～27.5%，遠超水體中TP存量的變化，而這部分磷在適宜條件下很容易再次釋放到水體中. 與王華等[5]從定性角度推斷上述因素為太湖TP波動的初步原因不同，TP質(zhì)量平衡分析對大洪水條件下各因素對太湖TP存量的貢獻率進行了定量. 下文將重點研究2015-2016年環(huán)太湖河道、水生植物及底泥變化情況，尤其是通過與多年平均的比較開展深入細化分析.

圖3 2015和2016年環(huán)太湖河道凈入湖總磷負荷量與多年平均值(2010-2017年)的比較Fig.3 Comparison of the net input of total phosphorus to Lake Taihu from 2015 to 2016 and the annual average from 2010 to 2017

2.2 環(huán)太湖河道進出太湖總磷負荷量變化的影響

將2015和2016年TP凈入湖量與多年平均值(2010-2017年)相比較，結(jié)果見圖3. 與多年平均值相比，2015年凈入湖量差別不大. 從環(huán)太湖河道所在行政區(qū)分析，常州略偏少，宜興和長興略偏多，蘇州出湖偏多；2016年除常州、無錫凈入湖量與多年平均值基本持平以外，其余地區(qū)均大幅偏多，其中宜興比多年平均偏多233.2 t，湖州為多年平均值的7倍，長興為多年平均值的2.7倍. 從空間分布上，2015和2016年TP的入湖負荷量以宜興地區(qū)入湖河道帶入為最多，所占比重分別約為65%和60%.

從具體環(huán)太湖河道分析，環(huán)太湖水文巡測段的入湖、出湖及凈入湖TP負荷見表3. 結(jié)合圖1可以看出，太湖西北部陳東港橋段和浯溪橋段的入湖TP負荷量最大，2015和2016年兩者負荷量之和分別占所有河道入湖量的77%和71%；而望虞河由于受到太湖高水位的影響，則由引水入湖為主改為排泄太湖洪水為主，兩年均為凈出湖.

環(huán)太湖河道TP負荷量年內(nèi)變化則與流域水情密切相關(guān). 2015和2016年太湖流域年降水量分別位列1951年以來第3位和第1位. 時間分布上，2015年在6月出現(xiàn)明顯的峰值，遠高于其他月份；2016年最高值同樣出現(xiàn)在6月，但9-10月降雨也異常偏多，太湖流域發(fā)生秋汛. 相應(yīng)于雨情，環(huán)太湖河道出入太湖的水量也呈現(xiàn)大進大出的格局，2015年入湖119.1億m3，出湖118.5億m3；2016年入湖159.9億m3，出湖167.3億m3. 時間分布上，2015年入湖水量主要集中在6月和7月，兩月的水量約占全年入湖總水量的35%；2016年入湖水量最大值同樣出現(xiàn)在6月和7月，兩月的水量約占全年入湖總水量的33%. 相應(yīng)地，環(huán)太湖河道入湖TP量也呈現(xiàn)類似的規(guī)律，與流域降雨量的變化呈現(xiàn)明顯的正相關(guān)關(guān)系(圖4). 2015年入湖TP量主要集中分布在6月和7月，兩月的入湖量約占全年入湖總量的37%；2016年入湖TP量最大值同樣出現(xiàn)在6月和7月，兩月的入湖量約占全年入湖總量的32%.

表3 2015和2016年環(huán)太湖河道出入湖總磷負荷量分布

圖4 2015和2016年環(huán)太湖河道入湖總磷量及流域降雨量的年內(nèi)變化Fig.4 Variation of total phosphorus flux from rivers into Lake Taihu and rainfall of the basin in 2015 and 2016

圖5 太湖水生植物分布面積(5月)與太湖總磷濃度的關(guān)系 Fig.5 Correlation of the total phosphorus concentration and the area (in May) with macrophytes in Lake Taihu

2.3 太湖沉水植物覆蓋度變化的影響

針對太湖沉水植物覆蓋度變化，分春、夏、秋3季連續(xù)6年，采用遙感方法結(jié)合現(xiàn)場調(diào)查進行反演計算[20]，開展了太湖沉水植物覆蓋度調(diào)查研究，鑒于資料的連續(xù)性和太湖水體黑臭出現(xiàn)的時間最早為5月，故采用5月數(shù)據(jù)進行趨勢分析. 2012-2016年太湖沉水植物5月平均分布面積與太湖TP濃度的關(guān)系見圖5.

2014年之前太湖東部湖區(qū)水草(沉水植物)分布廣泛，生長情況良好，沉水植物的根莖葉能吸收底泥和水體中的氮、磷營養(yǎng)物質(zhì)，促進水體懸浮物質(zhì)沉淀，抑制沉積物再懸浮，有利于水質(zhì)的改善；多數(shù)沉水植物葉片長至水面，能夠有效抑制并阻攔藍藻的生長和漂移. 2012-2014年太湖水生植物面積維持在250～350 km2，沉水植物面積維持在200～300 km2左右. 但因大規(guī)模打撈水草的過度收割作業(yè)以及連續(xù)高水位的不利影響，2015年5月太湖沉水植物面積驟減為27.65 km2，較2014年5月的244.31 km2大幅減少，降幅達到88.7%. 2016年5月面積略有上升，為47.33 km2，但仍比2014年同期少80.6%. 相對于分布面積的變化，太湖沉水植物年均單位面積生物量變化幅度相對較小，2012-2016年年際變化幅度為7%～31%.

通常，湖泊水體中的生物有效磷很低，但沉積物間隙水的磷酸鹽濃度是上覆水的9～600倍，有根大型水生植物體內(nèi)85%的磷是從沉積物間隙水中獲得的[21]. 張來甲等[22]研究發(fā)現(xiàn)，苦草腐解使3月水體 TP 總量較之前增長了144.40%, 而底泥中增長了19.99%, 可見植物腐解對水體的影響遠大于底泥. 當沉水植物腐爛后，植物體內(nèi)的磷主要釋放到水體中，而不是底泥中，因此，適時開展對水生生物的收割是十分必要的，但同時保持湖泊中適量的沉水植物對于減少磷含量起著重要作用. 而2015年大量收割太湖湖體中生長的沉水植物，湖體中的沉水植物分布面積和總生物量急劇減少，雖然有效帶出了TP，但湖中沉水植物的減少相應(yīng)減少了水生生物對磷的利用轉(zhuǎn)化，也增加了底泥向上覆水釋放磷的幾率.

2.4 底泥對太湖水體磷的影響

除水體中TP存量變化以外，滯留的TP一部分被水生生物吸收轉(zhuǎn)化，另一部分則沉降至底泥中. 底泥對水體磷的貢獻取決于3方面，即底泥中磷的含量、底泥受擾動的程度以及藍藻的泵吸作用.

2015年太湖沉水植物大幅下降，使得可被水生生物吸收轉(zhuǎn)化的TP量相應(yīng)減少，從而間接地增加了底泥中的TP存量. 研究表明，水生植物生長密集區(qū)底泥中TN、TP、總有機碳含量均顯著低于水生植物零星生長區(qū)[23]. 2015-2017年太湖底泥TP質(zhì)量比約為650 mg/kg，較1980s增加了約2.5倍，西部湖區(qū)的底泥TP質(zhì)量比更是達到1020 mg/kg，是1980s的5倍[24]. 不斷累積在底泥中的磷是太湖水體TP的重要內(nèi)源.

2015和2016年的大洪水和驟減的水生植物共同增強了底泥的受擾動程度，加劇了底泥中的磷向水體的遷移釋放. 在對底泥開展的靜態(tài)釋放和動態(tài)釋放試驗中發(fā)現(xiàn)，擾動可促使底泥發(fā)生大量釋放，而風(fēng)浪擾動是淺水湖泊沉積物懸浮和營養(yǎng)鹽釋放的主要驅(qū)動力之一，半天的強風(fēng)浪擾動可導(dǎo)致水體營養(yǎng)鹽濃度增加近1倍. 理論上，在大洪水條件下，底泥中高濃度的TP受到強風(fēng)浪擾動，更容易釋放到水體中.

依據(jù)遙感方法開展的太湖藍藻覆蓋面積數(shù)據(jù)，2015和2016年太湖藍藻暴發(fā)面積最大分別為1091.4和936.4 km2. 已有研究成果表明[25]，藍藻的快速生長導(dǎo)致其從底泥中泵取大量的磷，使得水體中顆粒磷的含量增加，水體與沉積物之間的磷平衡被打破，大量磷由沉積物進入水體，太湖中的磷含量與藍藻生長形成不斷放大的正向反饋. 2015和2016年太湖藍藻暴發(fā)與TP濃度的升高也進一步說明其正向反饋作用.

3 結(jié)論和建議

1)外源輸入TP負荷仍然是太湖水體中TP的主要來源. 2015-2016年流域為豐水年，尤其在2016年發(fā)生特大洪水的情況下，雖然環(huán)湖河道入湖水質(zhì)是近年來最好的一年，但由于總水量遠超常年，環(huán)湖河道進入太湖的TP負荷量仍為近年來的第2高值，達2594 t，占太湖TP負荷總量的74.2%. 而通過環(huán)湖河道帶入太湖的磷，絕大部分難以再通過河流帶出，只能通過人工收獲帶出、沉降、生物吸收轉(zhuǎn)化等途徑得以消減. 因此，對于累積效應(yīng)明顯的磷來說，河道輸入TP負荷量是太湖TP濃度升高的主要外源.

2)太湖沉水植物驟減導(dǎo)致湖體對磷的吸收轉(zhuǎn)化能力下降. 2015和2016年水生生物帶出TP負荷量為1280和1100 t，占比分別為38.7%和31.5%，雖然帶出了大量的磷負荷，但水生植物的減少一方面使水生植物吸收轉(zhuǎn)化磷的能力急劇下降，另一方面也增加了底泥向上覆水釋放的幾率. 不合理的過度收割水草，以及太湖持續(xù)高水位對水草生長的不利影響，成為太湖TP濃度升高的重要影響因素.

3)滯留在底泥中的TP在適宜條件下可進入水體，成為太湖水體中TP潛在的內(nèi)部來源. 2015和2016年，大量外源負荷進入太湖造成的TP滯留量達711和961 t，分別占總負荷的21.5%和27.5%，大洪水條件下加劇了底泥的受擾動程度，同時藍藻的快速生長與底泥中滯留的磷形成不斷放大的正向反饋，成為太湖TP濃度升高的又一重要影響因素.

4)2015-2016年太湖TP濃度升高，雖然體現(xiàn)了大洪水對湖泊生態(tài)系統(tǒng)的巨大沖擊，增加了TP外源的入湖以及內(nèi)源的擾動釋放，但氣象條件無法人為控制，對太湖TP濃度的控制主要還是要控源截污，從源頭減少TP輸入. 在今后的太湖水環(huán)境綜合治理中，應(yīng)進一步針對太湖TP的來源采取內(nèi)外兼控的控磷措施，強化對環(huán)太湖河流TP來源的排摸及控制、深化對底泥的影響及創(chuàng)新處置措施研究、科學(xué)合理開展水生生物收割，修復(fù)以沉水植被為主的太湖自然生態(tài)系統(tǒng).