李發(fā)榮，邱學(xué)禮，周 璟，劉 永，劉淑春，郭懷成

1.昆明市環(huán)境監(jiān)測(cè)中心，云南 昆明 650028

2.云南省農(nóng)業(yè)科學(xué)研究院，云南 昆明 650205

3.北京大學(xué)環(huán)境科學(xué)與工程學(xué)院，教育部水沙科學(xué)重點(diǎn)實(shí)驗(yàn)室，北京 100871

4.云南省高原湖泊流域污染過(guò)程與管理重點(diǎn)實(shí)驗(yàn)室，云南 昆明 650034

5.北京林業(yè)大學(xué)，北京 100083

目前，水體富營(yíng)養(yǎng)化是許多國(guó)家密切關(guān)注的環(huán)境問(wèn)題［1-2］。非點(diǎn)源污染是造成水體富營(yíng)養(yǎng)化的重要原因，隨著點(diǎn)源污染治理程度的逐步加深，非點(diǎn)源污染成為關(guān)注的焦點(diǎn)［3-4］。土壤以及人為活動(dòng)產(chǎn)生的氮、磷等營(yíng)養(yǎng)元素受到降雨的沖刷，隨著地表徑流及泥沙遷移匯入受納水體，造成水質(zhì)退化并產(chǎn)生富營(yíng)養(yǎng)化，由此帶來(lái)的環(huán)境、經(jīng)濟(jì)及社會(huì)問(wèn)題已引起國(guó)內(nèi)外普遍關(guān)注［5-9］。研究者主要通過(guò)實(shí)地監(jiān)測(cè)［9-10］和模型模擬［4，11-12］等研究非點(diǎn)源營(yíng)養(yǎng)負(fù)荷的排放規(guī)律，但由于受到氣候、降雨、植被、土壤、地形、耕作方式等多個(gè)因素的影響，非點(diǎn)源污染具有隨機(jī)性強(qiáng)、涉及范圍廣、時(shí)空差異大等特點(diǎn)，導(dǎo)致其監(jiān)測(cè)和模擬十分困難。流域尺度的模型預(yù)測(cè)是今后研究和評(píng)價(jià)非點(diǎn)源污染的主要技術(shù)手段和發(fā)展趨勢(shì)［15］，但該方面的研究需要建立在詳盡的實(shí)地監(jiān)測(cè)資料的基礎(chǔ)之上。在流域內(nèi)開(kāi)展小尺度的源區(qū)監(jiān)測(cè)研究，有助于識(shí)別地表徑流產(chǎn)生的關(guān)鍵源區(qū)及其水質(zhì)污染特征［16］，探討湖泊等水體的外源輸入，進(jìn)而有針對(duì)性地采取管理措施。實(shí)地監(jiān)測(cè)結(jié)果能更準(zhǔn)確地反映匯入受納水體的地表徑流來(lái)源及其污染特征，為流域尺度上的模擬預(yù)測(cè)及實(shí)現(xiàn)流域生態(tài)系統(tǒng)管理提供基礎(chǔ)。

滇池流域是云南省人口最密集、經(jīng)濟(jì)最發(fā)達(dá)的地區(qū)。由于毗鄰昆明主城區(qū)，滇池水質(zhì)受到人為活動(dòng)的嚴(yán)重影響，隨著污染物產(chǎn)生量迅速增加，滇池富營(yíng)養(yǎng)化已成為主要的環(huán)境問(wèn)題。滇池位于流域下游，寶象河、馬料河等20多條河流匯入湖體。已有研究表明，流域地表徑流使得大量的氮、磷等進(jìn)入滇池［17］，加重水體富營(yíng)養(yǎng)化并導(dǎo)致藍(lán)藻水華頻繁發(fā)生［18］。進(jìn)入滇池的營(yíng)養(yǎng)負(fù)荷主要來(lái)自工業(yè)廢水、城市污水等點(diǎn)源及農(nóng)田地表徑流流失和磷礦開(kāi)采［19］。隨著滇池流域點(diǎn)源污染控制工程的實(shí)施，農(nóng)田地表徑流流失、磷礦開(kāi)采成為滇池水體富營(yíng)養(yǎng)化控制的重要內(nèi)容［20-21］，探索并研究其相關(guān)規(guī)律已成為滇池治理的重點(diǎn)。因此，開(kāi)展對(duì)滇池東南岸農(nóng)業(yè)區(qū)和富磷采礦區(qū)的入湖河流野外布點(diǎn)監(jiān)測(cè)，通過(guò)實(shí)地采樣、降雨徑流監(jiān)測(cè)，結(jié)合室內(nèi)分析，識(shí)別各入湖河流的地表徑流及其水質(zhì)污染特征，探索其污染物的時(shí)空分布特征及其變化規(guī)律，以期為滇池污染物入湖時(shí)空分布規(guī)律研究、流域非點(diǎn)源污染的控制、滇池環(huán)境保護(hù)管理決策提供科學(xué)依據(jù)。

1 實(shí)驗(yàn)部分

1.1 研究區(qū)概況

滇池流域位于云貴高原的中部，流域面積2 920 km2，其海拔為1 880.6～2 837.6 m。滇池湖體面積約為 309 km2，占整個(gè)流域面積的10.6%。流域?qū)賮啛釒駶?rùn)氣候，年平均氣溫為14.7℃，年平均降雨量為1 007 mm，年平均相對(duì)濕度為73%。一年內(nèi)干濕季節(jié)分明，80%以上的年降雨集中在5—10月，汛期主要集中在7、8月。11月至次年4月為干季，日照充足，雨量稀少。流域多年平均陸地徑流量約為9.7 m3/s［19］，常年匯入滇池的有盤(pán)龍江、寶象河、馬料河等20多條入湖河流。農(nóng)業(yè)用地約占整個(gè)流域面積的20%;流域內(nèi)最具代表性的富磷區(qū)域位于滇池南岸的晉寧縣上蒜鄉(xiāng)，區(qū)域常年集中開(kāi)采磷礦［21］。大規(guī)模的磷礦開(kāi)采在創(chuàng)造良好經(jīng)濟(jì)效益的同時(shí)，卻造成了磷礦開(kāi)發(fā)區(qū)生態(tài)系統(tǒng)的嚴(yán)重破壞和流域水環(huán)境污染。

1.2 布點(diǎn)及采樣方案

滇池東南岸是整個(gè)流域農(nóng)業(yè)用地的主要分布區(qū)，農(nóng)業(yè)非點(diǎn)源污染也主要集中在這一區(qū)域。如圖1所示，選擇滇池東南岸10條主要入湖河流，在其下游設(shè)置監(jiān)測(cè)斷面，斷面靠近入湖口并不受滇池回水影響。由表1可見(jiàn)，這10條河流的流域面積幾乎占整個(gè)流域面積的一半，且具有相同的特征，即徑流量較大、匯水面積大、水土流失嚴(yán)重、非點(diǎn)源污染突出。其中，上蒜礦區(qū)所帶來(lái)的污染主要通過(guò)次巷河(8#)匯入滇池，因此對(duì)次巷河的監(jiān)測(cè)綜合反映了農(nóng)業(yè)面源和磷礦開(kāi)采帶來(lái)的污染特征。

由于7—9月為滇池流域集中降雨時(shí)間，從2008年8月16日至9月25日對(duì)上述10條河流進(jìn)行監(jiān)測(cè)。在無(wú)降雨?duì)顟B(tài)下，每隔5或7 d進(jìn)行一次常規(guī)監(jiān)測(cè)，8月24—26日期間出現(xiàn)了連續(xù)降雨，因而在降雨前(8月24日)和降雨后(8月27日)分別對(duì)各河流加測(cè)一次，共進(jìn)行了8次常規(guī)監(jiān)測(cè)。由于降雨存在空間差異性，為了解污染物在降雨過(guò)程中的流失特征，選取了降雨過(guò)程相同的3條河流，即大河(7#)、次巷河(8#)、甸心河(10#)開(kāi)展降雨徑流分析。共監(jiān)測(cè)了2次降雨徑流全過(guò)程，分別發(fā)生在8月25日和8月26日，降雨量分別為35、24 mm。在降雨過(guò)程中，從降雨產(chǎn)生徑流開(kāi)始，每隔10 min或30 min取1次樣，取樣頻率取決于河流水位上漲或回落的快慢程度，直至徑流水質(zhì)趨于穩(wěn)定。樣品采集時(shí)，將采樣器放至河道中，采集二分之一水深處的水樣，然后將水樣放到預(yù)先準(zhǔn)備好的水樣瓶中，并根據(jù)分析項(xiàng)目，添加固定劑、貼上標(biāo)簽，帶回實(shí)驗(yàn)室分析。

圖1 入湖河流監(jiān)測(cè)斷面分布圖

表1 河流基本概況

1.3 樣品分析及指標(biāo)測(cè)定

同時(shí)進(jìn)行水量、水質(zhì)指標(biāo)的分析。水量指標(biāo)包括降雨量、河流斷面面積、水位、流速、流量。水位現(xiàn)場(chǎng)采用水文標(biāo)尺測(cè)量，同時(shí)用流速儀按水文監(jiān)測(cè)規(guī)范進(jìn)行記錄，流速測(cè)定采用信控流速儀(XHW-1型)，流量則根據(jù)流速與斷面面積計(jì)算。常規(guī)監(jiān)測(cè)的水質(zhì)分析指標(biāo)包括 pH、總氮(TN)、總 磷 (TP)、溶解性磷(DP)、氨氮硝酸鹽氮化學(xué)需氧量(COD)、懸浮物(SS);在進(jìn)行降雨徑流監(jiān)測(cè)時(shí)，水質(zhì)分析指標(biāo)包括 pH、TN、TP、COD、SS、氟化物(F－)、硫化物(S－)、砷(As)、鉛(Pb)。分析方法按照《水和廢水監(jiān)測(cè)分析方法(第四版)》進(jìn)行，河流采樣按照《水質(zhì)采樣技術(shù)指導(dǎo)》(HJ 494—2009)進(jìn)行。

1.4 數(shù)據(jù)處理與研究方法

分別分析常規(guī)監(jiān)測(cè)和降雨徑流監(jiān)測(cè)的數(shù)據(jù)，8次常規(guī)監(jiān)測(cè)共有80個(gè)樣品，一次降雨過(guò)程中每條監(jiān)測(cè)河流的取樣數(shù)為7～10個(gè)，2次降雨徑流監(jiān)測(cè)共采集樣品52個(gè)。對(duì)于常規(guī)監(jiān)測(cè)，首先通過(guò)系統(tǒng)聚類的方法探討水質(zhì)污染的空間分布特征，進(jìn)而通過(guò)因子分析方法分析水質(zhì)影響因子，識(shí)別污染源。數(shù)據(jù)預(yù)處理與統(tǒng)計(jì)分析過(guò)程采用 SPSS.19.0。

1.4.1 系統(tǒng)聚類分析

系統(tǒng)聚類分析是一種應(yīng)用廣泛的聚類分析方法，其實(shí)質(zhì)是先將每一個(gè)個(gè)體看做一類，然后將相近程度最高的2類進(jìn)行合并組成一個(gè)新類，再將該新類與相似度最高的類進(jìn)行合并，直到所有的個(gè)體都?xì)w為一類。此方法主要用于分析流域的水質(zhì)時(shí)空變化特征［22］。分析采用歐氏距離平方，計(jì)算方法采用離差平方和法(Ward法)，即以平方歐式距離作為2類之間的距離，先將集合中每個(gè)樣本自成一類，在進(jìn)行類別合并時(shí)，計(jì)算類重心間方差，將離差平方和增加的幅度最小的2類首先合并，再依次將所有類別逐級(jí)合并。

1.4.2 因子分析

因子分析法通過(guò)研究眾多變量之間的內(nèi)部依賴程度，探求觀測(cè)數(shù)據(jù)之間的基本結(jié)構(gòu)，從多個(gè)實(shí)測(cè)的原變量中提取出較少的、互不相關(guān)的抽象綜合指標(biāo)，即因子，用于反映原來(lái)眾多變量的主要信息，每個(gè)原變量可以用這些提取出的因子的線性組合表示。此方法主要用于提取水質(zhì)污染因子和識(shí)別污染源［23］。

2 結(jié)果與討論

2.1 常規(guī)監(jiān)測(cè)

2.1.1 入湖河流流量及其水質(zhì)污染特征

在8次常規(guī)監(jiān)測(cè)中，各河流的流量情況見(jiàn)圖2。由圖2可知，新寶象河的流量最大，其平均流量達(dá)到2.6 m3/s，占總?cè)牒髁康?6.5%;其次為洛龍河，平均流量為1.47 m3/s，占總?cè)牒髁康?5.0%。入湖流量最低的是甸心河，平均流量?jī)H為0.16 m3/s。

圖2 河流入湖流量監(jiān)測(cè)結(jié)果

表2 河流水質(zhì)指標(biāo)統(tǒng)計(jì)描述

2.1.2 聚類分析

10條入湖河流水質(zhì)在空間上大體可分為3組(圖3):海河、新寶象河和馬料河、其余河流。在其余河流中，洛龍河、東大河、甸心河可歸為一組，大河、次巷河、南沖河、撈魚(yú)河又可歸為一組。結(jié)合水質(zhì)污染程度可看出，10條河流水質(zhì)狀況具有明顯的空間分布特征。海河最靠近昆明主城區(qū)，人類活動(dòng)與工業(yè)活動(dòng)最為頻繁，水資源開(kāi)發(fā)強(qiáng)度大，水質(zhì)污染程度最為嚴(yán)重。新寶象河與馬料河區(qū)域，仍然有一定的人為活動(dòng)，但工業(yè)活動(dòng)已略有減少，因而這2條河的水質(zhì)要好于海河。其余河流受人類活動(dòng)和工業(yè)活動(dòng)的影響相對(duì)較少，主要以農(nóng)業(yè)面源污染為主，因而水質(zhì)狀況在10條河流中相對(duì)較好。

圖3 入湖河流水質(zhì)聚類樹(shù)狀圖

2.1.3 因子分析

按照特征值大于1的原則，提取3個(gè)公因子用以反映入湖河流的水質(zhì)特征，即 F1、F2、F3，累計(jì)貢獻(xiàn)率為82.573%(表3)，能較好地反映原始數(shù)據(jù)的基本信息。F1貢獻(xiàn)率為53.636%，其中 TN、TP、DP、NO－3-N所占的因子載荷較大，且均與F1呈較強(qiáng)的正相關(guān)，主要代表了水體中的N、P含量;F2貢獻(xiàn)率為16.601%，其中 NH+4-N、COD、SS所占的因子載荷較大，且均與 F2呈正相關(guān)關(guān)系，主要為水體中SS及COD含量水平;F3貢獻(xiàn)率為12.336%，只有 pH在其中占有較大載荷且與F3呈較強(qiáng)的正相關(guān)關(guān)系，主要代表水體的酸堿度。

因子得分反映了各監(jiān)測(cè)斷面的污染狀況。得分越高，說(shuō)明水質(zhì)越差。由表4可見(jiàn)，海河綜合得分最高，污染程度也最為嚴(yán)重。位于南岸的東大河，得分最低，其水質(zhì)狀況最好。這一結(jié)果也反映了水質(zhì)的空間分異性，上游的海河、新寶象河、馬料河主要受F1影響，含有較高的因子得分，頻繁的人類活動(dòng)會(huì)使大量的 N、P排入河道，加重污染。東岸的大河、次巷河、南沖河、撈魚(yú)河等主要受F2影響，這些地區(qū)農(nóng)耕地較為集中，以農(nóng)業(yè)面源污染為主。下游的甸心河，在常規(guī)監(jiān)測(cè)中其平均入湖流量為最低，但由于臨近滇池出口，受下游人為活動(dòng)的影響，其水質(zhì)污染程度在10條河流中排在第4位，可見(jiàn)人為活動(dòng)與生活污水對(duì)流域河流水質(zhì)的影響較大。

表3 空間分析的旋轉(zhuǎn)因子載荷矩陣

表4 監(jiān)測(cè)斷面因子得分

2.2 降雨徑流監(jiān)測(cè)

通過(guò)對(duì)大河(7#)、次巷河(8#)、甸心河(10#)3條河流2次降雨徑流過(guò)程監(jiān)測(cè)，探討污染物的流失特征，監(jiān)測(cè)結(jié)果見(jiàn)圖4。由圖4可見(jiàn)，在降雨徑流的初期，污染物濃度開(kāi)始升高，中期增長(zhǎng)迅速，當(dāng)流量到達(dá)峰值或接近峰值時(shí)，污染物也開(kāi)始接近其濃度最高值，此后濃度開(kāi)始下降并最終趨于穩(wěn)定。COD、SS的濃度變化比 TN、TP劇烈得多，4個(gè)指標(biāo)的濃度變化范圍均為甸心河＞次巷河＞大河，這也反映了3條河流不同的狀況。在常規(guī)監(jiān)測(cè)中次巷河的水質(zhì)狀況要好于大河，但由于受磷礦開(kāi)采區(qū)的影響，在降雨條件下，其污染物濃度變化比大河高許多，降雨量較大時(shí)，這種現(xiàn)象更加明顯。甸心河由于受到下游人為活動(dòng)的影響，其污染物濃度變化最大。值得注意的是，在25日降雨過(guò)程中，3條河流SS濃度在流量到達(dá)峰值前其濃度最高值已經(jīng)出現(xiàn)。8月25日的降雨量(35 mm)大于8月26日(24 mm)。可見(jiàn)，當(dāng)降雨量較大時(shí)，由于降雨徑流初期的溶解沖刷作用，地表徑流能攜帶大量的SS進(jìn)入河流并迅速達(dá)到最大值，徑流中后期稀釋作用逐漸占主導(dǎo)地位，使得SS濃度開(kāi)始下降。表5列出了甸心河在8月25日降雨過(guò)程中的污染物濃度變化，由于 N、P等營(yíng)養(yǎng)元素大多吸附于SS上，因而使得污染物濃度也不斷增加。

圖4 甸心河、次巷河、大河在2次降雨過(guò)程中污染物濃度隨地表徑流的流失過(guò)程

表5 8月25日降雨過(guò)程中甸心河污染物濃度的變化

3條河流降雨徑流中，流量、SS與其余指標(biāo)的相關(guān)性見(jiàn)表6。由表6可知，在次巷河，流量與SS呈極顯著相關(guān)，相關(guān)系數(shù)達(dá)到0.93。除COD、As、Pb外，次巷河的流量、SS與其余水質(zhì)指標(biāo)均呈正相關(guān)或顯著正相關(guān)，這也說(shuō)明在次巷河，過(guò)大的地表徑流量將會(huì)帶來(lái)更為嚴(yán)重的水質(zhì)污染。除農(nóng)業(yè)面源外，磷礦開(kāi)采過(guò)程中會(huì)產(chǎn)生大量的揚(yáng)塵，降雨時(shí)會(huì)隨地表徑流產(chǎn)生大量水土流失，并攜帶營(yíng)養(yǎng)物質(zhì)進(jìn)入到滇池湖體。而在以農(nóng)業(yè)面源為主的大河，流量與SS呈顯著正相關(guān)(相關(guān)系數(shù)為0.66)，但流量、SS與其余指標(biāo)的相關(guān)性均不強(qiáng)。其產(chǎn)生原因在于:①在降雨過(guò)程中可能沒(méi)有采集足夠的樣本以正確反映污染物流失過(guò)程;②2次降雨持續(xù)時(shí)間都較長(zhǎng)，因而使得降雨強(qiáng)度不大，沒(méi)有帶來(lái)更多的污染物;③還與該河流控制面積內(nèi)的土地利用(如林地)和土壤性質(zhì)有關(guān)。在受下游人為活動(dòng)影響的甸心河，除 S－、Pb外，SS與其余指標(biāo)均為顯著相關(guān)，說(shuō)明地表徑流過(guò)大仍然會(huì)產(chǎn)生較嚴(yán)重的水質(zhì)污染。

表6 降雨徑流監(jiān)測(cè)中河流流量、SS與其余指標(biāo)的相關(guān)性

3 結(jié)論與建議

3.1 結(jié)論

1)在8次常規(guī)監(jiān)測(cè)中，新寶象河的平均流量為2.6 m3/s，占總?cè)牒髁康?6.5%;南岸的甸心河入湖流量最低，平均流量?jī)H為0.16 m3/s。最為靠近昆明主城區(qū)的海河，除NH4+-N外，其余指標(biāo)濃度均比其他河流高。TN、TP、COD、SS是滇池的主要污染指標(biāo)，許多河流均已嚴(yán)重超標(biāo)。

2)河流水質(zhì)狀況在空間上可分為3類，反映了不同河流的水質(zhì)污染程度;3類結(jié)果包括:①海河，靠近昆明主城區(qū)，其水質(zhì)污染最為嚴(yán)重;②海河的新寶象河、馬料河;③其余河流，水質(zhì)狀況相對(duì)較好。

3)N、P含量是河流水質(zhì)污染的主要貢獻(xiàn)因子，TN、TP、DP 和 NO3－-N對(duì)水質(zhì)污染的貢獻(xiàn)率達(dá)到了53.636%;因子得分結(jié)果表明，受人為活動(dòng)和工業(yè)活動(dòng)影響較大的海河、新寶象河、馬料河，是研究區(qū)域中污染最嚴(yán)重的3條河流，東大河水質(zhì)最好。

4)降雨條件下COD、SS濃度增長(zhǎng)尤為迅速，磷礦開(kāi)采對(duì)河流水質(zhì)的影響在降雨條件下更加明顯，其SS濃度在降雨條件下比只受農(nóng)業(yè)面源影響的河流最高高出1.9倍。流量、SS與大多數(shù)水質(zhì)指標(biāo)均有相關(guān)性，當(dāng)降雨量較大時(shí)，由于降雨徑流初期的溶解沖刷作用，SS能迅速達(dá)到最大值。

3.2 建議

研究發(fā)現(xiàn)，TN、TP、COD、SS是滇池的主要污染負(fù)荷，降雨情況下COD、SS增長(zhǎng)尤為迅速，將帶來(lái)更多的水質(zhì)污染。結(jié)合滇池以山地為主的地形特征，農(nóng)業(yè)面源污染防治的主要措施應(yīng)在進(jìn)入滇池河道的徑流區(qū)建立生態(tài)型溝渠，同時(shí)在周邊建立生態(tài)隔離帶，以阻截地表徑流攜帶的泥沙和負(fù)荷，并在隔離帶內(nèi)種植植物以吸收徑流中的養(yǎng)分;其次是在進(jìn)入湖泊的下游建設(shè)緩沖池，以沉降N、P等營(yíng)養(yǎng)物質(zhì)［25-26］，達(dá)到控制地表徑流，減少養(yǎng)分排放的目的。雖然經(jīng)過(guò)了生態(tài)溝渠、生態(tài)攔截系統(tǒng)以及緩沖滯留池的阻截，但仍會(huì)有一部分養(yǎng)分進(jìn)入河道，可再利用滇池的陸生湖濱帶地區(qū)建設(shè)生態(tài)型人工濕地以充分阻攔泥沙、吸收養(yǎng)分，從而有效消減N、P入湖負(fù)荷。監(jiān)測(cè)發(fā)現(xiàn)，河流水質(zhì)狀況呈現(xiàn)明顯的空間差異性，因而今后在進(jìn)行水環(huán)境模型建立、污染物目標(biāo)控制等研究時(shí)，必須考慮水質(zhì)參數(shù)的時(shí)空變化。不同區(qū)域的河流，如靠近主城區(qū)、以農(nóng)業(yè)面源為主和綜合農(nóng)業(yè)面源、磷礦開(kāi)采為主的河流，其水質(zhì)特征差異顯著，需根據(jù)具體的特征考慮不同的水質(zhì)參數(shù)。

［1］Zhang H，Huang G H.Assessment of non-point source pollution using a spatial multicriteria analysis approach［J］.Ecological Modelling，2011，222:313-321.

［2］Schaffner M，Bader H P，Scheidegger R.Modeling the contribution of point sources and non-point sources to Thachin River water pollution［J］.Science of the Total Environment，2009，407:4 902-4 915.

［3］李苗，王曉，龐宗強(qiáng)，等.沛沿河流域農(nóng)田地表徑流氮流失的模擬研究［J］.環(huán)境污染與防治，2010，32(12):18-21.

［4］張展羽，司涵，孔莉莉.基于SWAT模型的小流域非點(diǎn)源氮磷遷移規(guī)律研究［J］.農(nóng)業(yè)工程學(xué)報(bào)，2013，29(2):93-100.

［5］張玉斌，鄭粉莉，武敏.土壤侵蝕引起的農(nóng)業(yè)非點(diǎn)源污染研究進(jìn)展［J］.水科學(xué)進(jìn)展，2007，18(1):123-132.

［6］李瑞玲，張永春，曾遠(yuǎn)，等.太湖流域丘陵地區(qū)暴雨條件下農(nóng)田氮素隨地表徑流遷移特征［J］.農(nóng)業(yè)環(huán)境科學(xué)學(xué)報(bào)，2009，28(6):1 185-1 190.

［7］焦平金，王少麗，許迪，等.次暴雨下作物植被類型對(duì)農(nóng)田氮磷徑流流失的影響［J］.水利學(xué)報(bào)，2009，40(3):296-302.

［8］Ongley E D，Zhang X L，Yu T.Current status of agricultural and rural non-point source pollution assessment in China［J］.Environmental Pollution，2010，158:1 159-1 168.

［9］王菊翠，仵彥卿，黨碧玲，等.基于統(tǒng)計(jì)分析的陜西段涇河水質(zhì)時(shí)空分布特征［J］.自然資源學(xué)報(bào)，2012，27(4):674-685.

［10］蔣銳，朱波，唐家良，等.紫色丘陵區(qū)典型小流域暴雨徑流氮磷遷移過(guò)程與通量［J］.水利學(xué)報(bào)，2009，40(6):659-666.

［11］程紅光，岳勇，楊勝天，等.黃河流域非點(diǎn)源污染負(fù)荷估算與分析［J］.環(huán)境科學(xué)學(xué)報(bào)，2006，26(3):384-391.

［12］Maillard P，Antonia N，Santos P.A spatial-statistical approach for modeling the effect of non-point source pollution on different water quality parameters in the Velhas river watershed—Brazil［J］.Journal of Environmental Management，2008，86:158-170.

［13］張忠明，周立軍，宋明順，等.太湖苕溪流域農(nóng)業(yè)面源污染評(píng)價(jià)及對(duì)策［J］.環(huán)境污染與防治，2012，34(3):105-109.

［14］郭鴻鵬，朱靜雅，楊印生.農(nóng)業(yè)非點(diǎn)源污染防治技術(shù)的研究現(xiàn)狀及進(jìn)展［J］.農(nóng)業(yè)工程學(xué)報(bào)，2008，24(4):290-295.

［15］賴格英，于革.流域尺度的營(yíng)養(yǎng)物質(zhì)輸移模型研究綜述［J］.長(zhǎng)江流域資源與環(huán)境，2005，14(5):574-578.

［16］李立青，朱仁肖，郭樹(shù)剛，等.基于源區(qū)監(jiān)測(cè)的城市地表徑流污染空間分異性研究［J］.環(huán)境科學(xué)，2010，31(12):2 896-2 904.

［17］Gao L，Zhou J M，Yang H，et al.Phosphorus fractions in sediment profiles and their potential contributions to eutrophication in Dianchi Lake［J］.Environmental Geology，2005，48:835-844.

［18］Wu D L，Qian B，He L H.Contributing factor analysis of eutrophication of Dianchi Lake［J］.Research of Environmental Sciences，1992，5:26-28.

［19］余曉飛，和樹(shù)莊，胡斌，等.滇池柴河小流域暴雨徑流中COD的輸移特征研究［J］.環(huán)境污染與防治，2012，34(10):61-70.

［20］湛方棟，傅志興，楊靜，等.滇池流域套作玉米對(duì)蔬菜農(nóng)田地表徑流污染流失特征的影響［J］.環(huán)境科學(xué)學(xué)報(bào)，2012，32(4):847-855.

［21］孔燕，和樹(shù)莊，胡斌，等.滇池流域富磷地區(qū)暴雨徑流中磷素的沉降及輸移規(guī)律［J］.環(huán)境科學(xué)學(xué)報(bào)，2012，32(9):2 160-2 166.

［22］徐華山，徐宗學(xué)，唐芳芳，等.漳衛(wèi)南運(yùn)河流域水質(zhì)時(shí)空變化特征及其污染源識(shí)別［J］.環(huán)境科學(xué)，2012，33(2):359-369.

［23］卜紅梅，劉文治，張全發(fā).多元統(tǒng)計(jì)方法在金水河水質(zhì)時(shí)空變化分析中的應(yīng)用［J］.資源科學(xué)，2009，31(3):42-434.

［24］GB 3838—2002 中華人民共和國(guó)地表水環(huán)境質(zhì)量標(biāo)準(zhǔn)［S］.

［25］劉永，陽(yáng)平堅(jiān)，盛虎，等.滇池流域水污染防治規(guī)劃與富營(yíng)養(yǎng)化控制戰(zhàn)略研究［J］.環(huán)境科學(xué)學(xué)報(bào)，2012，32(8):1 962-1 972.

［26］劉瑞志，朱麗娜，雷坤，等.滇池入湖河流“十一五”綜合整治效果分析［J］.環(huán)境污染與防治，2012，34(3):95-100.