太湖流域上游平原河網(wǎng)區(qū)水質(zhì)空間差異與季節(jié)變化特征*

張 濤,陳求穩(wěn),易齊濤,王 敏,黃 蔚,馮然然

(1:南京水利科學(xué)研究院生態(tài)環(huán)境研究中心,南京 210029)(2:太湖流域管理局水利發(fā)展研究中心，上海 200434)

太湖流域上游平原河網(wǎng)區(qū)水質(zhì)空間差異與季節(jié)變化特征*

張 濤1,2,陳求穩(wěn)1**,易齊濤1,王 敏1,黃 蔚1,馮然然1

(1:南京水利科學(xué)研究院生態(tài)環(huán)境研究中心,南京 210029)(2:太湖流域管理局水利發(fā)展研究中心，上海 200434)

在太湖流域上游的宜溧-洮滆水系主要河道設(shè)置67個(gè)監(jiān)測(cè)點(diǎn),分別于2014年1月(冬季)、4月(春季)、8月(夏季)、11月(秋季)進(jìn)行水質(zhì)監(jiān)測(cè),采用多元統(tǒng)計(jì)方法分析了水質(zhì)的空間差異性和季節(jié)性變化,并利用水質(zhì)標(biāo)識(shí)指數(shù)法對(duì)水環(huán)境質(zhì)量進(jìn)行評(píng)價(jià). 結(jié)果表明,宜溧-洮滆水系污染程度較嚴(yán)重,總氮(TN)、總磷(TP)和高錳酸鹽指數(shù)(CODMn)濃度年均值分別為4.93、0.26和7.63 mg/L；單因素多元方差分析和聚類分析顯示污染物濃度具有顯著時(shí)空差異性,時(shí)間上冬、春季污染程度較高而夏、秋季較低,空間上無(wú)錫和常州氮、磷污染較為嚴(yán)重,宜興和溧陽(yáng)市有機(jī)污染程度較高；水質(zhì)標(biāo)識(shí)評(píng)價(jià)結(jié)果顯示流域內(nèi)水質(zhì)基本為IV類或V類,其中TN、TP及CODMn是關(guān)鍵污染指標(biāo).

太湖流域；宜溧-洮滆水系；多元統(tǒng)計(jì)分析；水質(zhì)評(píng)價(jià)；水質(zhì)標(biāo)識(shí)指數(shù)法

近幾十年來(lái), 隨著我國(guó)工農(nóng)業(yè)和城市化的快速發(fā)展, 強(qiáng)烈的人類活動(dòng)導(dǎo)致湖泊內(nèi)營(yíng)養(yǎng)鹽不斷富集[1],富營(yíng)養(yǎng)化程度加劇[2],有研究表明我國(guó)85%以上的湖泊處于富營(yíng)養(yǎng)化的水平[3]. 湖泊富營(yíng)養(yǎng)化及由此引起的藍(lán)藻水華已嚴(yán)重影響湖泊的生態(tài)功能與飲水安全[4],成為制約經(jīng)濟(jì)和社會(huì)持續(xù)發(fā)展的重要因素. 我國(guó)過(guò)去二十多年中針對(duì)富營(yíng)養(yǎng)化嚴(yán)重的大型湖泊開(kāi)展了一系列的整治措施[2],部分重點(diǎn)湖泊水質(zhì)有所改善,但水環(huán)境狀況總體惡化的趨勢(shì)尚未得到根本遏制[5]. 湖泊富營(yíng)養(yǎng)化治理主要在于外源控制和內(nèi)源削減,其中外源控制是根本. 因此解析外源營(yíng)養(yǎng)物質(zhì)來(lái)源,有效降低外部營(yíng)養(yǎng)物質(zhì)的輸入,實(shí)現(xiàn)流域水系水質(zhì)目標(biāo)管理是當(dāng)前湖泊富營(yíng)養(yǎng)化治理的首要任務(wù)[6-7]. 入湖河流是營(yíng)養(yǎng)物質(zhì)輸運(yùn)入湖的關(guān)鍵路徑,其水質(zhì)對(duì)湖泊水環(huán)境狀況有著重要而直接的影響,因此研究入湖河流水質(zhì)的變化規(guī)律,控制入湖河流污染物濃度,對(duì)改善湖泊水體富營(yíng)養(yǎng)化狀況具有重要意義[8].

目前,針對(duì)脈絡(luò)較為清晰且結(jié)構(gòu)相對(duì)簡(jiǎn)單的樹狀河流的水質(zhì)時(shí)空變化特征研究比較多[9-13]；太湖流域是我國(guó)典型的平原河網(wǎng)區(qū),流域城鎮(zhèn)化、工業(yè)化程度高,污染負(fù)荷嚴(yán)重,河網(wǎng)水文水動(dòng)力條件及負(fù)荷來(lái)源復(fù)雜,河流水質(zhì)影響因素眾多[14-15]. 此前入湖河流水質(zhì)監(jiān)測(cè)與分析主要集中于環(huán)湖河口及周邊區(qū)域[16-17],缺乏入湖河網(wǎng)多斷面、大尺度的水質(zhì)時(shí)空分布研究. 本研究選取太湖上游流域中的宜溧-洮滆水系作為對(duì)象，優(yōu)化設(shè)置高密度的水質(zhì)監(jiān)測(cè)網(wǎng)點(diǎn)，對(duì)河流水質(zhì)指標(biāo)進(jìn)行測(cè)定，并采用水質(zhì)標(biāo)識(shí)指數(shù)法評(píng)價(jià)河流水質(zhì)狀況，旨在揭示典型平原河網(wǎng)水質(zhì)空間差異和季節(jié)變化特征，以期為太湖流域入湖河流水質(zhì)目標(biāo)管理及污染負(fù)荷控制提供相關(guān)依據(jù).

1 研究區(qū)域與方法

1.1 研究區(qū)概況與樣點(diǎn)布設(shè)

宜溧-洮滆水系主要由宜溧河水系和洮滆水系兩部分組成. 宜溧河水系也稱南河水系，發(fā)源于茅山山區(qū)和蘇浙皖三省交界處的丘陵山地[18]，途徑無(wú)錫宜興市、溧陽(yáng)市，土地利用類型以耕地為主，干流長(zhǎng)50 km，下游北與洮滆水系相連，宜溧河水系入湖水量約占太湖上游來(lái)水總量的25%. 洮滆水系是由山區(qū)河道和平原河道組成的河網(wǎng)，集鎮(zhèn)江、丹陽(yáng)、金壇一帶的丘陵崗坡徑流，經(jīng)洮湖、滆湖調(diào)蓄后由常州太滆運(yùn)河、漕橋河、殷村港等河流匯入太湖[19]，同時(shí)又以丹金溧漕河、扁擔(dān)河、武宜運(yùn)河等多條南北向河道與沿江水系相通形成東西逢源、南北交匯的網(wǎng)絡(luò)狀水系. 洮滆水系入湖水量約占太湖上游來(lái)水總量的20%左右[20].

在宜溧-洮滆水系共布設(shè)67個(gè)水質(zhì)監(jiān)測(cè)站點(diǎn)(圖1). 監(jiān)測(cè)點(diǎn)布設(shè)覆蓋主要4級(jí)以上河道，并充分考慮水系交匯結(jié)構(gòu). 其中宜溧水系共布設(shè)25個(gè)采樣點(diǎn)，分布在大浦港(Y1、Y2)、東氿(Y3)、西氿(Y4)、北溪河(Y5、Y19、Y20、Y22)、溧宜河(Y6～Y11)、埝徑河(Y13、Y16)、桃溪(Y12、Y14、Y15)、戴溧河(Y17)、趙村河(Y18)、丹金溧漕河(Y21)、淳溧河(Y23)、中河(Y24)和南河(Y25). 洮滆水系共布設(shè)42個(gè)采樣點(diǎn)，分布在丹金溧漕河(DJ1～DJ4)、夏溪河(XX1～XX3)、扁擔(dān)河(BD1、BD2)、孟津河(MJ1～MJ3)、湟里河(HL1、HL2)、北干河(BG1、BG2)、中干河(ZG1、ZG2)、北河(BH1、BH2)、武宜漕河(WY1～WY4)、漕橋河(CQ1、CQ2)、太滆運(yùn)河(TG1～TG4)、殷村港(YC1、YC2)、燒香河(SX1、SX2)、橫塘河(HT1)、錫溧漕河(XL1～XL4)和直湖港(ZH1～ZH3).

圖1 宜溧-洮滆水系水質(zhì)監(jiān)測(cè)點(diǎn)分布Fig.1 The distribution of water sampling sites in the Yili-Taoge river networks

1.2 樣品采集與分析

1.3 水質(zhì)時(shí)空變異分析與評(píng)價(jià)

1.3.1 時(shí)空變異分析 在水質(zhì)評(píng)價(jià)過(guò)程中，考慮到水質(zhì)在時(shí)間、空間上的差異性及相似性，采用方差分析(ANOVA)的方法對(duì)各評(píng)價(jià)指標(biāo)分別進(jìn)行季節(jié)間和空間上的顯著差異性檢驗(yàn)，以判斷各檢測(cè)指標(biāo)在不同斷面以及不同季節(jié)間是否存在顯著差異[22]. 方差分析由Fisher等提出，用于兩個(gè)及以上樣本均數(shù)差別的顯著性檢驗(yàn)，其基本原理是將全部觀測(cè)數(shù)據(jù)的總變異分解為組間變異和組內(nèi)變異，如果組間變異遠(yuǎn)遠(yuǎn)大于組內(nèi)變異，則認(rèn)為控制因素對(duì)因變量產(chǎn)生了影響[23]. 本研究的控制變量分別為時(shí)間和空間，屬于單因素多元ANOVA，顯著性檢驗(yàn)公式為：

(1)

式中，n為樣本數(shù)；B為組間離差平方和矩陣；W為組內(nèi)離差平方和矩陣；p為向量維數(shù)；m為水平數(shù).Λ統(tǒng)計(jì)量服從自由度為(p，n-m，m-1)的Λ分布，當(dāng)λ2/n值大于顯著水平α?xí)r，則表明在該控制變量下不同水平各總體均值不存在顯著性差異，反之，則存在顯著性差異.

1.3.2 空間層次聚類分析 聚類分析是根據(jù)某些數(shù)量特征將觀察對(duì)象進(jìn)行分類的一種統(tǒng)計(jì)方法，其實(shí)質(zhì)是根據(jù)樣品或變量之間的親疏程度，通過(guò)逐次聚合，將性質(zhì)最接近的對(duì)象結(jié)合在一起聚成一類[11]. 在水質(zhì)評(píng)價(jià)中，考慮到水質(zhì)在時(shí)間、空間上的差異性及相似性，本研究采取應(yīng)用較為廣泛的層次聚類分析，親疏程度的計(jì)算包括兩類：樣本間距離和組間距離，其中計(jì)算樣本間距離采用歐氏距離平方和法，即兩點(diǎn)間實(shí)際距離的平方和；計(jì)算組間距離采用離差平方和法，又稱Ward法，其基本思想是以平方歐氏距離為標(biāo)準(zhǔn)，先將集合中每個(gè)樣本自成一類，計(jì)算類重心間方差，將離差平方和增加幅度最小的兩類進(jìn)行合并，每次通過(guò)合并減少一類，此時(shí)離差平方和出現(xiàn)并逐步增大，選擇使離差平方和增加最小的兩類合，隨后依次將所有類別逐級(jí)合并[24]. 聚類分析選取的水質(zhì)指標(biāo)包括耗氧污染物類(DO、CODMn及NH3-N)和營(yíng)養(yǎng)鹽類(TN及TP). 時(shí)空變異分析和空間層次聚類分析均使用SPSS 19.0軟件完成.

1.3.3 綜合水質(zhì)評(píng)價(jià) 宜溧-洮滆水系多為漁業(yè)、工業(yè)及農(nóng)業(yè)用水，其水功能區(qū)目標(biāo)以Ⅲ類或Ⅳ類為主，水質(zhì)評(píng)價(jià)采用綜合水質(zhì)標(biāo)識(shí)指數(shù)法(Iwq)[25]，標(biāo)識(shí)指數(shù)由整數(shù)位和3位或4位小數(shù)位組成，其結(jié)構(gòu)為：

Iwq=X1.X2X3X4

(2)

式中，X1.X2由計(jì)算獲得，X3和X4根據(jù)比較結(jié)果得到. 其中，X1為河流總體的綜合水質(zhì)類別；X2為綜合水質(zhì)在X1類水質(zhì)變化區(qū)間內(nèi)所處位置，從而實(shí)現(xiàn)在同類水中進(jìn)行水質(zhì)優(yōu)劣比較；X3為參與綜合水質(zhì)評(píng)價(jià)的水質(zhì)指標(biāo)中，劣于水環(huán)境功能區(qū)目標(biāo)的單項(xiàng)指標(biāo)個(gè)數(shù)；X4為綜合水質(zhì)類別與水體功能區(qū)類別的比較結(jié)果，視綜合水質(zhì)的污染程度，一位或兩位有效數(shù)字.

(3)

式中,m為參加綜合水質(zhì)評(píng)價(jià)的水質(zhì)單項(xiàng)指標(biāo)的數(shù)目；P1、P2、Pm分別為第1、2、m個(gè)水質(zhì)因子的單因子水質(zhì)指數(shù)，為對(duì)應(yīng)單因子水質(zhì)標(biāo)識(shí)指數(shù)中的整數(shù)位和小數(shù)點(diǎn)后第1位(單因子水質(zhì)標(biāo)識(shí)指數(shù)中的X1.X2).

通過(guò)綜合水質(zhì)標(biāo)識(shí)指數(shù)X1.X2，結(jié)合《國(guó)家地表水環(huán)境質(zhì)量標(biāo)準(zhǔn)》的分類，可以判斷水體的水質(zhì)類別和污染程度：當(dāng)1.0≤X1.X2≤2.0時(shí)，水質(zhì)狀況為Ⅰ類；當(dāng)2.07.0時(shí)，水質(zhì)狀況為劣Ⅴ類且黑臭.

2 結(jié)果與分析

2.1 流域水質(zhì)時(shí)空分布特征

宜溧-洮滆水系水質(zhì)總體呈弱堿性，DO濃度季節(jié)均值范圍為3.75～6.51 mg/L，部分監(jiān)測(cè)點(diǎn)濃度低于Ⅳ類水質(zhì)標(biāo)準(zhǔn)(3.0 mg/L). CODMn和TP濃度季節(jié)均值范圍分別為5.98～8.78和0.22～0.31 mg/L(表1)，均屬于Ⅳ類水，TN濃度季節(jié)均值范圍為3.23～6.89 mg/L，遠(yuǎn)高于Ⅴ類水質(zhì)標(biāo)準(zhǔn).

表1 宜溧-洮滆水系監(jiān)測(cè)站點(diǎn)水質(zhì)指標(biāo)季節(jié)均值統(tǒng)計(jì)

2.1.1 CODMn的時(shí)空分布特征 各監(jiān)測(cè)點(diǎn)CODMn濃度總體呈現(xiàn)出夏、冬季較高，春秋季較低的特征(圖2). 春季變化范圍為5.60～9.30 mg/L，平均為7.01 mg/L. 秋季CODMn濃度為全年最低，其變化區(qū)間在3.30～12.20 mg/L之間，均值為5.98 mg/L. 夏季CODMn變化范圍為4.10～21.60 mg/L之間，平均為8.71 mg/L. 冬季CODMn變化范圍為5.00～15.40 mg/L，平均為8.78 mg/L，與夏季相差甚微. 流域空間上也存在一定的變異性. 整體來(lái)看城鎮(zhèn)區(qū)域周圍CODMn濃度明顯高于流域上游和以農(nóng)業(yè)活動(dòng)為主的區(qū)域，在城鎮(zhèn)集中、工業(yè)發(fā)達(dá)的錫溧漕河—武宜運(yùn)河、扁擔(dān)—孟津河、溧宜河及環(huán)太湖河口一線濃度較高，而農(nóng)業(yè)為主的區(qū)域洮湖及滆湖周圍濃度相對(duì)較低.

圖2 宜溧-洮滆水系污染物濃度的時(shí)空分布Fig.2 The spatio-temporal distributions of pollutant concentration in the Yili-Taoge river networks

2.1.2 TN和DIN的時(shí)空分布特征 宜溧-洮滆水系TN濃度季節(jié)變化規(guī)律明顯，其大小順序?yàn)槎?春季>夏秋季. 各監(jiān)測(cè)站點(diǎn)TN濃度春季變化范圍為1.60～9.40 mg/L，平均為5.20 mg/L(表1). 夏、秋季流域內(nèi)大部分監(jiān)測(cè)點(diǎn)濃度均較低，其中夏季TN濃度變化范圍為1.40～6.20 mg/L，均值為3.23 mg/L；秋季TN濃度略有升高，平均為4.42 mg/L. 冬季TN濃度顯著增加，最大值高達(dá)11.30 mg/L，平均為6.89 mg/L. TN濃度的空間變異性較大且呈現(xiàn)明顯的分布規(guī)律(圖2). 從流域上游向下游流動(dòng)的過(guò)程中TN濃度逐漸升高，同時(shí)在城鎮(zhèn)、工業(yè)化集中的區(qū)域如京杭運(yùn)河、錫溧漕河—武宜運(yùn)河、扁擔(dān)—孟津河、夏溪河等濃度較高，濃度最高的河流集中于常州的武進(jìn)區(qū)周圍的水系.

圖3 宜溧-洮滆水系各形態(tài)氮占無(wú)機(jī)氮比例的時(shí)空分布Fig.3 The spatio-temporal distributions of percentages of different forms nitrogen to inorganic nitrogen in the Yili-Taoge river networks

2.1.3 TP的時(shí)空分布特征 宜溧-洮滆水系TP濃度季節(jié)變化主要表現(xiàn)為冬季最大，秋、春季次之，夏季最小(圖2). 其中，夏季TP濃度范圍為0.10～0.61 mg/L，均值0.22 mg/L；冬季濃度范圍為0.07～0.70 mg/L，均值為0.31 mg/L(表1). 空間上主要是運(yùn)河明顯高于其他區(qū)域，與TN一樣，污染最為嚴(yán)重區(qū)域集中于武進(jìn)區(qū)周邊水系，如錫溧漕河—武宜運(yùn)河、扁擔(dān)—孟津河和夏溪河等河流，流域下游區(qū)域河流要高于上游區(qū)域河流.

2.2 數(shù)據(jù)時(shí)空尺度多元方差分析

監(jiān)測(cè)數(shù)據(jù)時(shí)空多元方差分析結(jié)果如表2、表3所示. 4種不同統(tǒng)計(jì)量中Pillai’s Trace、Hotelling-Lawley及Roy’s Greatest Root統(tǒng)計(jì)值越大，表明變量對(duì)結(jié)果的影響越大，而Wilks’ Lambda統(tǒng)計(jì)量則相反. 以時(shí)間多元方差分析中Wilks’ Lambda為例，P(Sig.=0.000)<0.05，說(shuō)明不同季節(jié)內(nèi)監(jiān)測(cè)數(shù)據(jù)差異顯著，偏Eta方值為0.430，說(shuō)明可以解釋變異的43.0%. 結(jié)果顯示4種不同統(tǒng)計(jì)量在時(shí)間尺度及空間尺度上P(Sig.=0.000)值均小于 0.05，因此在變量DO、CODMn、NH3-N、TN和TP的聯(lián)合分布下，不同監(jiān)測(cè)斷面在不同季節(jié)內(nèi)均存在顯著差異.

表2 采樣點(diǎn)時(shí)間多元方差分析*

*P<0.05，下同.

表3 采樣點(diǎn)空間多元方差分析

2.3 聚類分析結(jié)果

針對(duì)耗氧污染物類指標(biāo)(DO、CODMn、NH3-N)和營(yíng)養(yǎng)鹽類指標(biāo)(TN、TP)，利用Ward法對(duì)各監(jiān)測(cè)點(diǎn)進(jìn)行聚類分析，結(jié)果如表4和圖4所示. 耗氧污染物類指標(biāo)分為3組,其中,A1樣本量為12個(gè),主要聚集了CODMn較高樣本點(diǎn),主要分布在溧陽(yáng)市的部分宜溧水系、埝徑河及部分入湖河口附近；A2聚集了較高DO、同時(shí)NH3-N濃度較低的樣本點(diǎn),主要分布在流域上游的金壇市及丹陽(yáng)市內(nèi),即流域污染較小的區(qū)域；A3聚集了NH3-N濃度較高的樣本點(diǎn),其個(gè)數(shù)占所有樣本點(diǎn)的60%左右,表明流域內(nèi)NH3-N污染較為嚴(yán)重,空間分布主要在城市、工業(yè)化較為集中的區(qū)域.

營(yíng)養(yǎng)鹽類指標(biāo)(TN、TP)分為3組,污染程度順序?yàn)锽1

表4 不同聚類組水質(zhì)指標(biāo)均值

圖4 宜溧-洮滆水系不同水質(zhì)指標(biāo)聚類分析結(jié)果Fig.4 The results of cluster analysis of different water-quality indexs in the Yili-Taoge river networks

2.4 水質(zhì)評(píng)價(jià)與特征污染物分析

利用綜合水質(zhì)標(biāo)識(shí)指數(shù)法對(duì)宜溧-洮滆水系監(jiān)測(cè)期間內(nèi)水質(zhì)狀況進(jìn)行評(píng)價(jià),結(jié)果顯示春季除極個(gè)別站點(diǎn)達(dá)到Ⅱ類或Ⅲ類水標(biāo)準(zhǔn),其余各站點(diǎn)均為Ⅳ類或Ⅴ類,其中Ⅴ類水主要分布在流域下游的太滆運(yùn)河、漕橋河、殷村港、錫溧漕河及直湖港；夏、秋季水質(zhì)狀況有所改善,僅少數(shù)站點(diǎn)為Ⅴ類水,其余大部分監(jiān)測(cè)站點(diǎn)水質(zhì)類別為Ⅳ類,部分達(dá)到Ⅲ類,空間差異性較小；冬季水質(zhì)狀況較其他3季明顯惡化,僅少量站點(diǎn)為Ⅳ類,其余站點(diǎn)均為Ⅴ類(圖5). 根據(jù)不同季節(jié)綜合水質(zhì)標(biāo)識(shí)指數(shù)均值可以看出,流域內(nèi)水質(zhì)以Ⅳ類及Ⅴ類為主.

不同季節(jié)DO、CODMn、NH3-N、TN及TP各指標(biāo)單因子水質(zhì)標(biāo)識(shí)指數(shù)均值分布情況(圖6),DO及NH3-N指數(shù)值相對(duì)較低,以Ⅱ類、Ⅲ類水質(zhì)為主；而CODMn、TN及TP的指數(shù)值明顯偏高,其中TN指數(shù)均值在不同季節(jié)均超過(guò)5.0,總體處于Ⅴ類水質(zhì)標(biāo)準(zhǔn)以上. 耗氧污染類指標(biāo)中,CODMn平均超過(guò)Ⅳ類標(biāo)準(zhǔn),為主控指標(biāo)；營(yíng)養(yǎng)鹽指標(biāo)則需要對(duì)TP和TN進(jìn)行控制,其中TN控制是改善平原河網(wǎng)區(qū)整體水質(zhì)的關(guān)鍵所在.

3 討論

圖5 宜溧-洮滆水系綜合水質(zhì)標(biāo)識(shí)指數(shù)評(píng)價(jià)結(jié)果時(shí)空分布Fig.5 The spatio-temporal distributions of assessment result class using comprehensive water quality identification index in the Yili-Taoge river networks

圖6 宜溧-洮滆水系單因子水質(zhì)標(biāo)識(shí)指數(shù)季節(jié)分布Fig.6 Seasonal distribution of single-factor water identification index in the Yili-Taoge river networks

從季節(jié)變化來(lái)看,在2015年宜溧—洮滆水系污染程度整體呈現(xiàn)冬、春季較高,夏、秋季較低的特征. 不同季節(jié)水質(zhì)變化主要受其溫度、降水、水文條件等影響[35-36],在冬、春兩季宜溧—洮滆水系污染較為嚴(yán)重,其中在冬季水質(zhì)指標(biāo)CODMn、TN、TP濃度均值分別達(dá)到8.78、6.89、0.31 mg/L,基本處于Ⅴ類或劣Ⅴ類水質(zhì)標(biāo)準(zhǔn),春季污染物濃度較冬季有所降低,但仍然處于高值范圍,一方面冬、春兩季降雨量小導(dǎo)致河流水位下降,稀釋凈化能力相對(duì)降低,另一方面冬季微生物活動(dòng)的減弱,使得流域內(nèi)污水處理廠污染物處理效率降低,同時(shí)河流自凈能力減弱也是污染物濃度高的重要原因. 與此相反,夏、秋兩季豐水期對(duì)河流氮、磷污染物濃度產(chǎn)生了較強(qiáng)的稀釋作用,同時(shí)水溫增加，流域及河道內(nèi)微生物作用的增強(qiáng)，能夠降低水體營(yíng)養(yǎng)鹽濃度.

宜溧—洮滆水系水質(zhì)污染程度在空間上大體呈現(xiàn)北高南低、東高西低的分布特征，同時(shí)利用層次聚類的方法能夠很好表征水質(zhì)指標(biāo)的空間分布規(guī)律. 流域內(nèi)兩類指標(biāo)均聚成3類，具有顯著的空間差異性，同時(shí)營(yíng)養(yǎng)鹽類指標(biāo)與耗氧污染物類指標(biāo)空間分布具有一定相似性. 位于流域下游的錫溧漕河—武宜運(yùn)河、太滆運(yùn)河以及殷村港、漕橋河等入湖河口附近樣本點(diǎn)營(yíng)養(yǎng)鹽類指標(biāo)及NH3-N污染均最為嚴(yán)重，這些河流主要位于太湖流域城鎮(zhèn)化和工業(yè)化高度集中的區(qū)域(即無(wú)錫、常州以及宜興市的交匯區(qū)域)，體現(xiàn)了強(qiáng)烈的人為活動(dòng)影響. 土地利用是影響污染物排放量的重要因素[17]，其一方面通過(guò)不同人類活動(dòng)方式和強(qiáng)度影響營(yíng)養(yǎng)鹽的輸入量，另一方面則通過(guò)改變地表糙率影響地表徑流過(guò)程. 已有研究表明,林地面積的增加對(duì)水質(zhì)改善發(fā)揮積極作用，而農(nóng)田和城鎮(zhèn)用地面積的增加則使水質(zhì)趨于惡化[32]. 近些年，太湖流域土地利用情況發(fā)生明顯改變，與2000年相比，2008年農(nóng)田占地率從60%下降至51%，城鎮(zhèn)用地則由13%上升至22%[15]. 其中無(wú)錫、常州市的土地利用類型均以城鎮(zhèn)用地及耕地為主，林地、水域所占比例均較低，區(qū)域內(nèi)人口居住密度大，因此認(rèn)為大量的工業(yè)點(diǎn)源污染排放、生活污水排放及不完善的污水處理系統(tǒng)是造成該區(qū)域附近水域氮、磷污染嚴(yán)重的主要原因[37]. 同時(shí)由于受河流流向影響，污染物在入湖河口附近呈逐漸累積現(xiàn)象，從而導(dǎo)致流域下游區(qū)域污染物濃度偏高. 此外，距離監(jiān)測(cè)點(diǎn)最近的土地利用類型對(duì)水質(zhì)的影響較區(qū)域綜合土地利用情況對(duì)水質(zhì)的影響更大[38]，由于污染嚴(yán)重區(qū)域主要集中于入湖河口附近，流量大且污染物衰減距離短，也可能是限制太湖水質(zhì)提高的重要因素.

綜合水質(zhì)指數(shù)標(biāo)識(shí)法能夠較為客觀地反映河流的綜合水質(zhì)類別、與水功能區(qū)類別的比較結(jié)果等信息[39]. 評(píng)價(jià)結(jié)果表明宜溧—洮滆水系主要以Ⅲ～Ⅴ類水質(zhì)為主，總體能夠滿足區(qū)域水體功能目標(biāo). 單因子水質(zhì)標(biāo)識(shí)指數(shù)結(jié)果表明,TN、TP、CODMn為目前宜溧—洮滆水系特征污染物，以氮為主要控制指標(biāo)，與余輝等[28]研究結(jié)果一致，因此控氮是提高研究區(qū)域水環(huán)境質(zhì)量的關(guān)鍵所在. 同時(shí)，從太湖湖泊生態(tài)系統(tǒng)來(lái)看，氮的控制也具有重要意義. 傳統(tǒng)觀點(diǎn)認(rèn)為磷是湖泊水體生產(chǎn)力的主要限制因子，我國(guó)自1998年頒布對(duì)太湖禁磷政策以來(lái)，磷污染程度有所下降，對(duì)太湖富營(yíng)養(yǎng)化進(jìn)一步惡化起到延緩的作用，但近些年的研究表明太湖藻類生長(zhǎng)受到氮、磷的共同限制[40-43]，夏、秋季以氮限制為主，冬、春季以磷限制為主. 從流域輸入上可以看出，夏、秋季營(yíng)養(yǎng)鹽濃度較低而冬、春季較高，而湖泊水體中夏、秋季內(nèi)源磷負(fù)荷較大，加上氮的反硝化作用，使得水體受到氮限制，因此氮外源輸入控制對(duì)減輕太湖藍(lán)藻水華具有重要的意義.

4 結(jié)論

2)基于營(yíng)養(yǎng)鹽和耗氧污染物兩類指標(biāo)，利用單因素多元方差分析和聚類分析結(jié)果表明，宜溧—洮滆水系不同水質(zhì)指標(biāo)濃度均具有顯著的時(shí)空差異性. 季節(jié)上主要受溫度、降雨量等因素影響，水系污染程度總體呈現(xiàn)冬、春季較高，夏、秋季較低特征；空間上污染程度大體呈現(xiàn)北高南低、東高西低的分布特征，且空間聚集規(guī)律明顯；受土地利用類型和地理位置影響，無(wú)錫、常州、宜興交匯區(qū)域?yàn)榈?、磷污染最為?yán)重區(qū)域，也是污染控制的關(guān)鍵區(qū).

[1] Fan Ronggui, Zhu Dongnan, Deng Lan. Causes of eutrophication and progress of integrated management technology in lake.JournalofWaterResources&WaterEngineering, 2010,21(6): 48-52. [范榮桂, 朱東南, 鄧嵐. 湖泊富營(yíng)養(yǎng)化成因及其綜合治理技術(shù)進(jìn)展. 水資源與水工程學(xué)報(bào), 2010,21(6): 48-52.]

[2] Xu Qigong, Cao Jinling, Gao Rutaietal. Trend of water quality deterioration and eutrophication control phases partition in China.EnvironmentalScience&Technology, 2011,34(11): 147-151. [許其功, 曹金玲, 高如泰等. 我國(guó)湖泊水質(zhì)惡化趨勢(shì)及富營(yíng)養(yǎng)化控制階段劃分. 環(huán)境科學(xué)與技術(shù), 2011,34(11): 147-151.]

[3] Wu Yali, Xu Hai, Yang Guijunetal. Progress in nitrogen pollution research in Lake Taihu.JLakeSci, 2014,26(1): 19-28. DOI:10.18307/2014.0103. [吳雅麗, 許海, 楊桂軍等. 太湖水體氮素污染狀況研究進(jìn)展. 湖泊科學(xué), 2014,26(1): 19-28.]

[4] Xu Hai, Qin Boqiang, Zhu Guangwei. Nutrient limitation of cyanobacterial growth in different regions of Lake Taihu in summer.ChinaEnvironmentalScience, 2012,32(12): 2230-2236. [許海, 秦伯強(qiáng), 朱廣偉. 太湖不同湖區(qū)夏季藍(lán)藻生長(zhǎng)的營(yíng)養(yǎng)鹽限制研究. 中國(guó)環(huán)境科學(xué), 2012,32(12): 2230-2236.]

[5] Wang Wei. Study and application digital simulation on river water quality management by computational intelligence [Dissertation]. Harbin: Harbin Institute of Technology, 2013. [王祎. 基于計(jì)算智能方法的河流水質(zhì)管理數(shù)字模擬研究與應(yīng)用[學(xué)位論文]. 哈爾濱: 哈爾濱工業(yè)大學(xué), 2013.]

[6] Zhao Tingting, Li Qiuyan, Lu Liqiaoetal. Analysis of water environmental pollution in plain river network region in smallwatershed of Caoqiao River.TransactionsoftheCSAE, 2011,27(Supp.2): 170-175. [趙婷婷, 李秋艷, 陸麗巧等. 漕橋河小流域平原河網(wǎng)地區(qū)水環(huán)境污染分析. 農(nóng)業(yè)工程學(xué)報(bào), 2011,27(S2): 170-175.]

[7] Pan Hongbo. Lake eutrophication and its prevention measures.EnvironmentalScienceandTechnology, 2011,24(S1): 123-126. [潘紅波. 湖泊富營(yíng)養(yǎng)化問(wèn)題及其防治淺議. 環(huán)境科技, 2011,24(S1): 123-126.]

[8] Xu Mei, Ren Ruili, Liu Maosong. Annual changes of water quality in an upstream river of Taihu Lake.JournalofNanjingForestryUniversity:NaturalSciencesEdition, 2007,31(6): 121-124. [許梅,任瑞麗,劉茂松. 太湖入湖河流水質(zhì)指標(biāo)的年變化規(guī)律. 南京林業(yè)大學(xué)學(xué)報(bào): 自然科學(xué)版, 2007,31(6): 121-124.]

[9] Jing Hongwei, Zhang Zhigang, Guo Jing. Water pollution characteristics and pollution sources of Bei Canal river system in Beijing. China.EnvironmentalScience, 2013,33(2): 319-327. [荊紅衛(wèi), 張志剛, 郭婧. 北京北運(yùn)河水系水質(zhì)污染特征及污染來(lái)源分析. 中國(guó)環(huán)境科學(xué), 2013,33(2): 319-327.]

[10] Wang XZ, Cai QH, Ye Letal. Evaluation of spatial and temporal variation in stream water quality by multivariate statistical techniques: A case study of the Xiangxi River basin, China.QuaternaryInternational, 2012,282(2): 137-144.

[11] Wang Gang, Li Zhaofu, Wan Rongrongetal. Analysis of temporal and spatial variations in water quality of Xitiaoxi Watershed using multivariate statistical techniques.JournalofAgro-EnvironmentScience, 2015,34(9): 1797-1803. [王剛, 李兆富, 萬(wàn)榮榮等. 基于多元統(tǒng)計(jì)分析方法的西苕溪流域水質(zhì)時(shí)空變化研究. 農(nóng)業(yè)環(huán)境科學(xué)學(xué)報(bào), 2015,34(9): 1797-1803.]

[12] Xu Huashan, Xu Zongxue, Tang Fangfangetal. Spatio-temporal variation analysis and identification of water pollution sources in the Zhangweinan River Basin.EnvironmentalScience, 2012,33(2): 359-369. [徐華山, 徐宗學(xué), 唐芳芳等. 漳衛(wèi)南運(yùn)河流域水質(zhì)時(shí)空變化特征及其污染源識(shí)別. 環(huán)境科學(xué), 2012,33(2): 359-369.]

[13] Xie Ying, Li Xuyong, Wang Huiliangetal. Stream pollution analysis and hydro-chemical assessment of the upper Luanhe River Basin.ActaScientiaeCircumstantiae, 2012,32(3): 645-653. [解瑩, 李敘勇, 王慧亮等. 灤河流域上游地區(qū)主要河流水污染特征及評(píng)價(jià). 環(huán)境科學(xué)學(xué)報(bào), 2012,32(3): 645-653.]

[14] Li YP, Acharya K, Stone MCetal. Spatio-temporal patterns in nutrient loads, nutrient concentrations, and algal biomass in Lake Taihu, China.LakeandReservoirManagement, 2011,27(4): 298-309.

[15] Xu H, Yang LZ, Zhao GMetal. Anthropogenic impact on surface water quality in Taihu lake region, China.Pedosphere, 2009,19(6): 765-778.

[16] Lu Shaoyong, Jiao Wei, Wang Qiangetal. Spatial-temporal distribution characteristics of water quality of rivers around Taihu Lake.ResearchofEnvironmentalSciences, 2011,24(11): 1220-1225. [盧少勇, 焦偉, 王強(qiáng)等. 環(huán)太湖河流水質(zhì)時(shí)空分布特征. 環(huán)境科學(xué)研究, 2011,24(11): 1220-1225.]

[17] Niu Yong, Yu Hui, Zhang Minetal. Impact of land-using pattern on water quality of rivers in Lake Taihu watershed.ResourcesandEnvironmentintheYangtzeBasin, 2013,22(2): 205-211. [牛勇, 余輝, 張敏等. 太湖流域典型土地利用方式下入湖河流水質(zhì)污染特征研究. 長(zhǎng)江流域資源與環(huán)境, 2013,22(2): 205-211.]

[18] Shi Yunxiang. Water environment evolution tends and the pollutant transport process research of Yili River [Dissertation]. Nanjing: Nanjing University of Science and Technology, 2002. [史云祥. 宜溧河流域水環(huán)境演變趨勢(shì)以及污染物輸移過(guò)程研究[學(xué)位論文]. 南京: 南京理工大學(xué), 2002.]

[19] Gu Xiejun, Xu Dongtong, Cai Huanxingetal. Feature analysis and comparison on ecological risk of eutrophication in lakes and reservoirin Taoge River.EnvironmentalMonitoringandForewarning, 2011,3(3): 13-17. [顧謝軍, 徐東炯, 蔡煥興等. 洮滆水系湖庫(kù)富營(yíng)養(yǎng)化生態(tài)風(fēng)險(xiǎn)的特點(diǎn)與比較. 環(huán)境監(jiān)控與預(yù)警, 2011,3(3): 13-17.]

[20] Zhang Yiming, Zhang Yongchun, Gao Yuexiangetal. Water pollution control technology and strategy for pollution control of water body of River-Lake System.JournalofEcologyandRuralEnvironment, 2010,26(S1): 9-13. [張毅敏, 張永春, 高月香等. 河湖相連水系水體污染控制技術(shù)與策略. 生態(tài)與農(nóng)村環(huán)境學(xué)報(bào), 2010,26(S1): 9-13.]

[21] Ministry of Environmental Protection of the People’s Republic of China, Editorial Board of Water and Wastewater Monitoring and Analysis Methods eds. Water and Wastewater Monitoring and Analysis Methods: 4th edition. Beijing: China Environmental Science Press, 2002. [國(guó)家環(huán)境保護(hù)總局《水和廢水監(jiān)測(cè)分析方法》編委會(huì). 水和廢水監(jiān)測(cè)分析方法:第4版. 北京: 中國(guó)環(huán)境科學(xué)出版社, 2002.]

[22] Fu Tianyi, Zou Zhihong, Wang Xiaojing. Water quality assessment for Taizi River watershed in Liaoyang section based on multivariate statistical analysis and water quality identification index.ActaScientiaeCircumstantiae, 2014,34(2): 473-480. [富天乙, 鄒志紅, 王曉靜. 基于多元統(tǒng)計(jì)和水質(zhì)標(biāo)識(shí)指數(shù)的遼陽(yáng)太子河水質(zhì)評(píng)價(jià)研究. 環(huán)境科學(xué)學(xué)報(bào), 2014,34(2): 473-480.]

[23] Harrington PB, Vieira NF, Espinoza Jetal. Analysis of variance-principal component analysis: A soft tool for proteomic discovery.AnalyticaChimicaActa, 2005,544(1): 118-127.

[24] Tang Linan, Liu Yu, Pan Yuchunetal. Evaluation and zoning of rural regional multi-function based on BP model and Ward method: A case in the Pinggu district of Beijing city.ScientiaGeographica, 2016,36(10): 1514-1521. [唐林楠, 劉玉, 潘瑜春等. 基于BP模型和Ward 法的北京市平谷區(qū)鄉(xiāng)村地域功能評(píng)價(jià)與分區(qū). 地理科學(xué), 2016,36(10): 1514-1521.]

[25] Xu Zuxin. Comprehensive water quality identification index for environmental quality assessment of surface water.JournalofTongjiUniversity:NaturalScience, 2005,33(4): 482-488. [徐祖信. 我國(guó)河流綜合水質(zhì)標(biāo)識(shí)指數(shù)評(píng)價(jià)方法研究. 同濟(jì)大學(xué)學(xué)報(bào): 自然科學(xué)版, 2005,33(4): 482-488.]

[26] Cheng Bo, Zhang Ze, Chen Lingetal. Eutrophication of Taihu Lake and pollution from agricultural non-point sources in Lake Taihu Basin.JournalofAgro-EnvironmentScience, 2005,24(S1): 118-124. [程波, 張澤, 陳凌等. 太湖水體富營(yíng)養(yǎng)化與流域農(nóng)業(yè)面源污染的控制. 農(nóng)業(yè)環(huán)境科學(xué)學(xué)報(bào), 2005,24(S1): 118-124.]

[27] Wang Xue, Xu Hui, Yan Shuwenetal. Spatial distribution characteristics of up-stream river pollution in lakeTaihu basin.ResourcesandEnvironmentintheYangtzeBasin, 2012,21(3): 341-348. [王雪, 余輝, 燕姝雯等. 太湖流域上游河流污染空間分布特征研究. 長(zhǎng)江流域資源與環(huán)境, 2012,21(3): 341-348.]

[28] Yu Hui, Yan Shuwen, Xu Jun. Multivariate statistical analysis of water quality in the inflow and outflow rivers of Lake Taihu.ResourcesandEnvironmentintheYangtzeBasin, 2010,19(6): 696-702. [余輝, 燕姝雯, 徐軍. 太湖出入湖河流水質(zhì)多元統(tǒng)計(jì)分析. 長(zhǎng)江流域資源與環(huán)境, 2010,19(6): 696-702.]

[29] Jarvie HP, Whitton BA, Neal C. Nitrogen and phosphorus in east-coast British rivers: speciation, sources and biological significance.ScienceoftheTotalEnvironment, 1998, S210-211(1-6): 79-109.

[30] Ravichandran S. Hydrological influences on the water quality trends in Tamiraparani Basin, south India.EnvironmentalMonitoringandAssessment, 2003,87(3): 293-309.

[31] Xian G, Crane M, Junshan S. Analysis of urban development and its environmental impact on the Tampa Bay watershed.JournalofEnvironmentalManagement, 2007,85(4): 965-976.

[32] Hu Jian, Liu Maosong, Zhou Wenetal. Correlations between water quality and land use pattern in Taihu Lake basin.ChineseJournalofEcological, 2011,30(6): 1190-1197. [胡建, 劉茂松, 周文等. 太湖流域水質(zhì)狀況與土地利用格局的相關(guān)性. 生態(tài)學(xué)雜志, 2011,30(6): 1190-1197.]

[33] Xu Hai, Liu Zhaopu, Jiao Jiaguoetal. Nitrogen pollution status of various types of passing-by water bodies in upper reaches of Taihu Lake.ChineseJournalofEcological, 2008,27(1): 43-49. [許海, 劉兆普, 焦佳國(guó)等. 太湖上游不同類型過(guò)境水氮素污染狀況. 生態(tài)學(xué)雜志, 2008,27(1): 43-49.]

[34] Liu Yun, Yi Qitao, Chen Qiuwenetal. The spatial distribution of nitrogen species in surface sediment in the west river networks of the Taihu Lake Basin, China.ActaScientiaeCircumstantiae, 2015,35(12): 3890-3897. [劉蕓, 易齊濤, 陳求穩(wěn)等. 太湖西部河網(wǎng)中沉積物氮的空間分布特征. 環(huán)境科學(xué)學(xué)報(bào), 2015,35(12): 3890-3897.]

[35] Vega M, Pardo R, Barrado Eetal. Assessment of seasonal and polluting effects on the quality of river water by exploratory data analysis.WaterResearch, 1998,32(12): 3581-3592.

[36] Ouyang Y, Nkedi-kizza P, Wu QTetal. Assessment of seasonal variations in surface water quality.WaterResearch, 2006,40(20): 3800-3810.

[37] Xu Pengzhu, Qin Boqiang, Huang Wenyuetal. Assessment of water quality and nutrition to water bodies the Yilihe Watershed, Taihu Basin.JLakeSci, 2001,13(4): 315-321. DOI:10.18307/2001.0405. [許朋柱, 秦伯強(qiáng), 黃文鈺等. 太湖流域宜溧河地區(qū)水體水質(zhì)狀況及營(yíng)養(yǎng)狀態(tài)評(píng)價(jià). 湖泊科學(xué), 2001,13(4): 315-321.]

[38] Cao Fangfang, Li Xue, Wang Dong. Effects of land use structure on water quality in Xin’anjiang River.EnvironmentalScience, 2013,34(7): 2582-2587. [曹芳芳, 李雪, 王東. 新安江流域土地利用結(jié)構(gòu)對(duì)水質(zhì)的影響. 環(huán)境科學(xué), 2013,34(7): 2582-2587.]

[39] Xie Ying, Li Xuyong, Wang Huiliangetal. Stream pollution analysis and hydro-chemical assessment of the upper Luanhe River Basin.ActaScientiaeCircumstantiae, 2012,32(3): 645-653. [解瑩, 李敘勇, 王慧亮等. 灤河流域上游地區(qū)主要河流水污染特征及評(píng)價(jià). 環(huán)境科學(xué)學(xué)報(bào), 2012,32(3): 645-653.]

[40] Paerl HW, Hall NS, Calandrino ES. Controlling harmful cyanobacterial blooms in a world experiencing anthropogenic and climatic-induced change.ScienceoftheTotalEnvironment, 2011,409(10): 1739-1745.

[41] Paerl HW, Xu H, Hall NS. Nutrient limitation dynamics examined on a multiannual scale in Lake Taihu, China: implications for controlling eutrophication and harmful algal blooms.FreshwaterEcology, 2015,30(1): 5-24.

[42] Paerl HW, Xu H, McCarthy MJetal. Controlling harmful cyanobacterial blooms in a hyper-eutrophic lake (Lake Taihu, China): The need for a dual nutrient (N & P) management strategy.WaterResource, 2011,45(5): 1973-1983.

[43] Ma JR, Qin BQ, Wu Petal. Controlling cyanobacterial blooms by managing nutrient ratio and limitation in a large hypereutrophic lake: Lake Taihu, China.JournalofEnvironmentScience, 2015,27(1): 80-86.

SpatialandseasonalvariationsofwaterqualityintheupstreamplainrivernetworksoftheTaihuBasin

ZHANG Tao1,2, CHEN Qiuwen1**, YI Qitao1, WANG Min1, HUANG Wei1& FENG Ranran1

(1:CenterforEco-EnvironmentalResearch,NanjingHydraulicResearchInstitute,Nanjing210029,P.R.China)(2:WaterConservancyDevelopmentResearchCenter,TaihuBasinAuthority,Shanghai200434,P.R.China)

Sixty-seven water sampling sites in the Yili-Taoge plain river networks of the upstream area of the Taihu Basin were set to characterize its water quality with spatiotemporal distribution. Water sampling trips were conducted in January (winter), April (spring), August (summer), November (autumn) of 2014, and main water quality parameters were concentrated on nutrients and organic matters. Spatial and seasonal variations of water quality were analyzed by using the multivariate statistical methods and water quality were assessed by employing the comprehensive water quality identification index. The results indicated that pollution degree in the research area was severe. The averaged concentrations of total nitrogen (TN), total phosphorus (TP) and permanganate index (CODMn) were 4.93 mg/L, 0.26 mg/L and 7.63 mg/L, respectively. One-way multivariate analysis and cluster analysis showed significantly spatiotemporal patterns of pollutants. Pollution levels in winter and spring were higher than that in summer and autumn. Pollution of nitrogen and phosphorus were heavier in the Wuxi and Changzhou areas, whereas organic matter pollution was more serious in the Yixing and Liyang areas. Water quality in the upstream plain river networks can meet the criteria of Grade IV or V of National Standards for Surface Water Quality, characterized by pollution of TN, TP and CODMn.

Taihu Basin; Yili-Taoge water system; multivariate statistical analysis; water quality assessment; water quality identification index

*國(guó)家自然科學(xué)基金項(xiàng)目(51579149,51609142,51609144)和江蘇省水利廳水利科技重大項(xiàng)目(2015005,2016021,2016030)聯(lián)合資助. 2016-11-11收稿；2017-01-24收修改稿. 張濤(1993～)，女，碩士研究生; E-mail： zt7914375@163.com.

**通信作者; E-mail：qwchen@nhri.cn.