徐麗賢, 梁新強(qiáng),, 周柯錦, 林麗敏,王知博, 華桂芬, 王欣煜, 王志榮, 黃 武

(1.浙江大學(xué) 環(huán)境與資源學(xué)院, 浙江 杭州 310058; 2.浙江省水體污染控制與環(huán)境安全技術(shù)重點實驗室,浙江 杭州 310058; 3.浙江省環(huán)境檢測中心, 浙江 杭州 310015; 4. 浙江省農(nóng)業(yè)生態(tài)與能源辦公室, 浙江 杭州 310012)

?

山地農(nóng)業(yè)小流域非點源氮磷輸出特征及來源

徐麗賢1, 梁新強(qiáng)1,2, 周柯錦3, 林麗敏1,王知博1, 華桂芬2, 王欣煜1, 王志榮4, 黃 武4

(1.浙江大學(xué) 環(huán)境與資源學(xué)院, 浙江 杭州 310058; 2.浙江省水體污染控制與環(huán)境安全技術(shù)重點實驗室,浙江 杭州 310058; 3.浙江省環(huán)境檢測中心, 浙江 杭州 310015; 4. 浙江省農(nóng)業(yè)生態(tài)與能源辦公室, 浙江 杭州 310012)

摘要：[目的] 分析非點源污染物輸出特征及來源,為相似農(nóng)業(yè)小流域水質(zhì)研究以及非點源污染控制提供參照基礎(chǔ)。[方法] 通過ArcGIS軟件對集水流域劃分,監(jiān)測2014年斷面及降雨水質(zhì),結(jié)合平均濃度法及輸出系數(shù)法,建立考慮降雨攜帶輸出的流域非點源負(fù)荷輸出系數(shù)法計算模型,并進(jìn)一步采用最優(yōu)化數(shù)學(xué)方法對流域內(nèi)污染物的來源進(jìn)行分析。[結(jié)果] 核算得到耕地、林地、城鎮(zhèn)村及工礦用地、農(nóng)村生活和畜禽養(yǎng)殖的總氮輸出系數(shù)分別為15.87,6.33,6.27 kg/(hm2·a),0.20 kg/(人·a),0.83 kg/(頭·a),總磷輸出系數(shù)分別為0.46,0.39,0.67 kg/(hm2·a),0.10 kg/(人·a),0.16 kg/(頭·a)。該研究區(qū)域徑流中-N是氮流失的主要形式,降雨高峰期徑流中顆粒態(tài)磷流失嚴(yán)重,-N及溶解態(tài)磷是該流域雨水中氮磷的主要存在形態(tài)。[結(jié)論] 該流域非點源污染輸出以降雨、林地、農(nóng)村生活輸出為主。

關(guān)鍵詞：非點源污染; 降雨輸出; 系數(shù)核算; 氮磷流失特征

文獻(xiàn)參數(shù)： 徐麗賢, 梁新強(qiáng), 周柯錦, 等.山地農(nóng)業(yè)小流域非點源氮磷輸出特征及來源[J].水土保持通報，2016,36(3)：30-37.DOI:10.13961/j.cnki.stbctb.2016.03.007

農(nóng)業(yè)非點源污染是影響水質(zhì)的關(guān)鍵因素，對其輸出特征及來源進(jìn)行分析是控制污染的前提。自2001年出臺畜禽飼養(yǎng)相關(guān)規(guī)定后，非點源污染受到人們的關(guān)注越來越大，其中，農(nóng)業(yè)非點源污染被認(rèn)為是水生態(tài)環(huán)境的重要影響因素[1]，氮、磷被普遍認(rèn)為是農(nóng)業(yè)非點源的特征污染物[2-3]。不同于城鎮(zhèn)生活污水與工業(yè)廢水污染，農(nóng)業(yè)非點源污染具有發(fā)生的隨機(jī)性和間歇性，復(fù)雜的機(jī)制與過程，不確定的排放渠道及量值，多變的空間及時間污染負(fù)荷，監(jiān)測、估算及控制的困難性等一系列特點[4]。對國內(nèi)外流域水環(huán)境模型，即非點源污染模型進(jìn)行分類，按照物理過程描述的不同分為經(jīng)驗型模型以及機(jī)理機(jī)制型模型；按照空間離散程度或分辨率大小的不同分為集總式和分布式模型[5]。其中，輸出系數(shù)模型即為經(jīng)驗型集總式模型，該模型最大的特點是利用相對容易得到的研究流域土地利用狀況、人類生產(chǎn)活動情況等資料直接建立各因素與收納水體非點源污染負(fù)荷的關(guān)系[6]。研究表明，空氣沉降等因素也可進(jìn)一步被考慮添加進(jìn)模型中，提高模型應(yīng)用過程中流域尺度非點源污染估算和預(yù)測的準(zhǔn)確性[7]。同時，此類模型通過結(jié)合GIS等工具也被廣泛應(yīng)用于非點源污染負(fù)荷核算中[8]。對于缺乏水文水質(zhì)資料的特定流域中，若參考其他流域輸出系數(shù)的方法，往往不符合研究區(qū)域的實際情況，容易給結(jié)果造成一定偏差。不同土地利用等因素也造成非點源污染輸出特征的差異性[9-10]，因此，清晰辨別研究區(qū)域污染源以及確定各污染源的輸出系數(shù)值是應(yīng)用該方法的關(guān)鍵。輸出系數(shù)法雖常被應(yīng)用于大尺度流域，鑒于本研究區(qū)域為集水流域可劃分、農(nóng)業(yè)污染源明確以及無工業(yè)點源污染的流域，且各個監(jiān)測斷面形狀規(guī)則利于監(jiān)測，本文通過定點監(jiān)測與輸出系數(shù)模型相結(jié)合的方法，在全城塢小流域進(jìn)行2014年1—12月各斷面水質(zhì)及大氣濕沉降(降雨)實地監(jiān)測的前提下，研究該地區(qū)年際非點源氮磷輸出的主要特征并核算驗證各污染源非點源氮磷污染輸出系數(shù)。以期為相似農(nóng)業(yè)小流域水質(zhì)研究以及非點源污染控制提供參照基礎(chǔ)。

1研究區(qū)概況

全城塢村小流域(119°41′50″—119°43′51″,30°27′37″—30°29′24″)位于浙江省杭州市余杭區(qū)鸕鳥鎮(zhèn)，東鄰黃湖，南連臨安市橫畈，北接百丈，西接安吉縣山川鄉(xiāng)，占地面積約659.92 hm2，人口為1 028人。調(diào)查研究表明，該小流域內(nèi)在監(jiān)測時段中無點源污染，因此將該流域的污染均視為非點源污染。同時經(jīng)實地考察，全城塢流域內(nèi)非點源污染來源主要包括：耕地、林地、城鎮(zhèn)村及工礦用地、農(nóng)村生活、畜禽養(yǎng)殖、降雨輸出。通過ArcGIS 10.2軟件生成研究區(qū)域的水系及集水流域，并進(jìn)一步依據(jù)輸出系數(shù)模型所需參數(shù)個數(shù)，采用嵌套式布點的方法在河流中布設(shè)水質(zhì)監(jiān)測斷面，分別位于該地區(qū)劃分得到集水流域的出口，其編號分別為D1，D2，D3，D4，D5。

本研究所在地區(qū)的耕地土壤屬于潴育水稻土(黃斑田)，林地土壤屬于黃紅壤(黃泥沙土)，土壤氮磷背景值如表1所示。通過對該研究區(qū)域的實地調(diào)查獲得農(nóng)村生活及畜禽養(yǎng)殖數(shù)據(jù)，解譯2014年的OLI-TIRS影像數(shù)據(jù)，并根據(jù)分類規(guī)范(GB/T21010-2007)獲得該全城塢流域的土地利用情況。進(jìn)一步用Arc-Map中的區(qū)域分析功能對各斷面所代表的集水流域進(jìn)行耕地面積、林地面積、城鎮(zhèn)村及工礦用地面積進(jìn)行統(tǒng)計，得到如表2所示信息。

表1　研究區(qū)土壤性質(zhì)

表2　研究區(qū)集水流域基本信息

注：畜禽養(yǎng)殖的頭數(shù)通過豬當(dāng)量換算得到，其中肉雞、蛋雞、羊、肉牛、奶牛的豬當(dāng)量折算系數(shù)分別為1/60，1/30，1/3，5，10。由于解譯精度受限，同時該研究區(qū)域農(nóng)用地(耕地+林地)的面積占93%，交通運(yùn)輸用地等土地利用臨近城鎮(zhèn)村及工礦用地且占地面積小，因此在本研究解譯及統(tǒng)計過程中將其歸于該流域的城鎮(zhèn)村及工礦用地中。D1，D2，D3，D4，D5表示5個斷面。下同。

2材料與方法

2.1樣品測試分析方法

本研究選取2014年1—12月中每月代表日8:00至翌日8:00進(jìn)行該流域水樣采集及流量的測定，其中水樣每30 min采集1次。采集完畢后以每個水樣采集時段的累計流量為體積權(quán)重制備當(dāng)次監(jiān)測混合水樣；流量采用STARFLOW超聲波多普勒流量計，在各斷面監(jiān)測現(xiàn)場進(jìn)行實時測定，每5 s記錄1次讀數(shù)，本研究將每月15日所監(jiān)測的數(shù)值用作于輸出系數(shù)模型的核算，每月5日所監(jiān)測的數(shù)值用作于該模型的驗證。

參照《大氣降水采樣與分析方法》GB13 580.2-92在該研究區(qū)域布設(shè)2個點，利用聚乙烯塑料小桶放置監(jiān)測斷面D2，D3附近相對高度為1.5 m的空曠地段進(jìn)行采樣，視降雨頻次而定(2014年每月監(jiān)測降雨2～7場)，監(jiān)測指標(biāo)以各場降雨量為權(quán)重計算得到每月降雨平均濃度，每次水樣取3個平行樣，在現(xiàn)場進(jìn)行預(yù)處理后，統(tǒng)一由高密度聚乙烯瓶置于移動冰箱冷藏保存后，帶回實驗室立即進(jìn)行各個指標(biāo)的測試。

2.2輸出負(fù)荷核算方法

通過實地調(diào)查，綜合考慮土地利用形式、農(nóng)村生活、畜禽養(yǎng)殖及降雨輸出等因素，利用ArcGIS 10.2軟件對研究區(qū)域進(jìn)行集水流域劃分，進(jìn)一步在集水流域出口設(shè)立監(jiān)測斷面，并通過斷面一整年的監(jiān)測數(shù)據(jù)，在Johns等[7]提出的輸出系數(shù)模型基礎(chǔ)上，以軟件劃分獲得的集水流域為基本計算單元建立該流域考慮降雨輸出的非點源負(fù)荷輸出系數(shù)法計算模型，并通過最優(yōu)化數(shù)學(xué)方法進(jìn)行參數(shù)率定，得到流域非點源各污染物的輸出系數(shù)，同時對各斷面污染物輸出特征及研究區(qū)域非點源污染物來源進(jìn)行分析。其中，以各斷面所控制的集水流域為基本計算單元(各集水流域出口設(shè)置監(jiān)測斷面)，根據(jù)每次各斷面的水量、水質(zhì)同步監(jiān)測數(shù)據(jù)，以徑流量為權(quán)重，求出加權(quán)平均濃度近似作為斷面徑流的平均濃度，乘以斷面徑流量得到非點源污染負(fù)荷量[11]，并通過核算時段時間轉(zhuǎn)換系數(shù)得到年非點源負(fù)荷量，其計算表達(dá)式為：

(1)

2.3考慮降雨輸出的模型建立

確定各污染源合理的輸出系數(shù)值是輸出系數(shù)模型在特定流域中應(yīng)用時的重點。影響流域非點源污染輸出系數(shù)的因素包括地形地貌、水文氣候等自然因素，也包括管理措施以及人類活動等人為因素[6]。本研究斷面水質(zhì)水量同步監(jiān)測的基礎(chǔ)下，同時對該研究區(qū)域的降雨情況進(jìn)行了監(jiān)測，將降雨監(jiān)測得到污染物濃度與斷面流量的乘積作為降雨污染物輸出的負(fù)荷(通常，降雨輸出為雨水濃度、降雨量及徑流系數(shù)三者的乘積，本研究中將斷面流量視為該斷面所覆蓋區(qū)域的降雨量與徑流系數(shù)的乘積)。確定污染物的類型和污染源的種類是決定模型的基本結(jié)構(gòu)和輸出系數(shù)個數(shù)的前提[6]，進(jìn)一步利用監(jiān)測及統(tǒng)計數(shù)據(jù)進(jìn)行參數(shù)率定得到全城塢流域各個非點源污染源的輸出系數(shù)，模型的一般表達(dá)式為[12]：

(2)式中：k——集水流域的類別；j——核算流域中污染源的種類，共n個；i——斷面及降雨監(jiān)測次數(shù)，共m次；Lk——污染物在流域k中的總負(fù)荷量(kg/a)；Ej——污染物在第j種污染源的輸出系數(shù)〔kg/(hm2·a)〕，其中農(nóng)村生活的單位為〔kg/(人·a)〕，畜禽養(yǎng)殖的單位為〔kg/(頭·a)〕；Aj——第j種污染源的數(shù)量；Ci——第i月降雨中污染物平均濃度(kg/m3)；Ka——降雨負(fù)荷核算時段時間轉(zhuǎn)換系數(shù)。

根據(jù)該流域的特點，識別主要污染物的類型有TN，TP，主要的污染源種類有耕地、林地、城鎮(zhèn)村及工礦用地、農(nóng)村生活、畜禽養(yǎng)殖及降雨輸出。當(dāng)降雨輸出負(fù)荷已知時，計算模型需要確定的參數(shù)的個數(shù)是5。因此將該流域劃分為5個集水流域，根據(jù)公式可以列出5個方程，組成線性方程如下：

(3)式中：L——集水流域的某種污染物的非點源污染負(fù)荷，利用平均濃度法[12]可以計算出其值；A——污染源的數(shù)量，根據(jù)調(diào)查資料為已知值；以月為單位，斷面及降雨同步監(jiān)測共計12次，公式中共有5個未知數(shù)，解此線性方程組可求得輸出系數(shù)，即E1，E2，E3，E4，E5。

3結(jié)果與分析

3.1流域降雨量及徑流量變化特征

降雨產(chǎn)生的徑流是往往是造成流域氮磷流失的重要因素[13]。自2014年1月1日至2014年12月31日，全城塢流域全年累積降雨量1 438.5 mm，其中5—8月的降雨量占全年降雨量的54%(全年最大的次降雨事件發(fā)生在7月27日，降雨量為87 mm，最大雨強(qiáng)為30 mm/h)。本研究區(qū)域降雨較多且雨量充沛。在該年的1，2，11月降雨較少，該3個月累計降雨量僅占全年降雨量的18%。

監(jiān)測結(jié)果顯示(圖1)，各斷面流量波動情況具有相似性。自6月份起，由于梅雨季節(jié)造成的持續(xù)性降雨，導(dǎo)致該區(qū)域在后3個月中各斷面流量增大，其中7月份各斷面流量達(dá)到峰值，因此降雨也是導(dǎo)致流量變化的關(guān)鍵因素，降雨量年內(nèi)分布不均造成流量呈現(xiàn)出季節(jié)性變化，流量高峰主要集中在夏季。各斷面間的流量差異性主要由于集雨面積以及農(nóng)村生活排水造成，D1—D5集雨面積不斷增大以及農(nóng)村生活排水量不斷增加，造成：D5流量>D4流量>D3流量>D2流量>D1流量的結(jié)果。

圖1　2014年全城塢流域日降雨量及各監(jiān)測斷面流量變化情況

3.2降雨及斷面監(jiān)測氮磷變化特征

從圖3可以得出，TP含量為0.03～0.14 mg/L，其平均值為0.08 mg/L，處于地表Ⅱ類或Ⅲ類水標(biāo)準(zhǔn)(以P計)。斷面監(jiān)測結(jié)果表明，在降雨高峰期斷面監(jiān)測TP濃度波動大，且主要以顆粒態(tài)磷(DTP)流失為主，其中斷面1及斷面4全年磷流失以顆粒態(tài)為主，其原因可能是該研究區(qū)域磷含量普遍較低，該監(jiān)測斷面又受區(qū)域農(nóng)村生活污水及河道淤泥擾動的影響大；降雨水質(zhì)監(jiān)測結(jié)果表明，該地區(qū)降雨中磷流失主要以溶解態(tài)形式為主。

將監(jiān)測得到2014年各月的污染物濃度(圖2及圖3)及流量(圖1)水文數(shù)據(jù)代入污染物年輸出負(fù)荷公式(1)，其監(jiān)測斷面的污染物年輸出負(fù)荷如表4所示。監(jiān)測斷面D1—D5分別代表5個集水流域的出口，其中斷面1及斷面2代表的集水流域相對獨立，其余監(jiān)測斷面所代表的集水流域相互嵌套。由于隨著流域面積的增加，各污染源的數(shù)量也隨之增加，監(jiān)測斷面D3—D5污染物總流失量呈逐漸遞增的趨勢。研究結(jié)果表明，降雨輸出是非點源污染物流失的重要組成部分，監(jiān)測研究區(qū)域降雨中的污染物水平，結(jié)合各斷面的流量數(shù)據(jù)得到各集水流域的降雨輸出負(fù)荷，降雨中的TN水平較TP水平相比普遍較高。核算得到D1—D5降雨輸出負(fù)荷分別為TN：2 540.20，4 751.67，6 208.18，13 284.94，14 410.61 kg；TP：15.04，28.15，36.78，78.69，85.00 kg。

3.3參數(shù)率定及驗證

將2014年各監(jiān)測斷面非點源污染物總流失負(fù)荷(表4)、降雨輸出負(fù)荷與各污染源統(tǒng)計信息(表2)代入輸出系數(shù)模型計算公式(3)中，通過最優(yōu)化的數(shù)學(xué)方法得到耕地、林地、城鎮(zhèn)村及工礦用地、農(nóng)村生活、畜禽養(yǎng)殖的TN，TP輸出系數(shù)核算結(jié)果(表3)，并在表3中統(tǒng)計了其他學(xué)者的輸出系數(shù)研究成果。耕地的TN及TP輸出系數(shù)相比其他流域的研究結(jié)果數(shù)值偏小，可能是由于該地區(qū)的耕地的種植主要以單季稻為主，在其他流域耕地的種植主要以多季稻等作物為主。農(nóng)村生活的TN，TP輸出系數(shù)明顯高于其他流域，通過實地調(diào)查，該研究區(qū)域農(nóng)村生活污水的處理相對簡單，主要以“三格式”化糞池或者直排的方式入河，這是造成農(nóng)村生活輸出系數(shù)偏高的主要原因。

圖2　全城塢小流域2014年非點源TN月流失量

通過數(shù)據(jù)核算得到污染物流失總負(fù)荷以及降雨輸出負(fù)荷，結(jié)合各斷面所代表的集水流域基本信息(表2)，根據(jù)各污染源的輸出系數(shù)(表3)及各污染源的統(tǒng)計信息(表2)相乘得到輸出負(fù)荷占比情況如表3。其中，降雨輸出負(fù)荷是通過實測降雨資料結(jié)合公式(2)中的降雨輸出部分核算獲得。結(jié)合表3可知，全城塢流域各污染源TN輸出負(fù)荷由大到小分別為：降雨輸出>林地>城鎮(zhèn)村及工礦用地>耕地>農(nóng)村生活>畜禽養(yǎng)殖，其中降雨中TN輸出負(fù)荷所占比例較大為75.82%，由此可見該地區(qū)降雨中TN的流失不容忽視，其次為林地輸出，占20.01%，實地考察得到該地區(qū)以林地為主且面積較大，部分林地施用有機(jī)肥是造成林地TN輸出較大的主要原因；各污染源的TP輸出負(fù)荷由大到小分別為：林地>農(nóng)村生活>降雨輸出>城鎮(zhèn)村及工礦用地>畜禽養(yǎng)殖>耕地，其中林地TP輸出占49.16%，其次為農(nóng)村生活TP輸出占總輸出的21.35%，因此，對林地、農(nóng)村生活中TP的輸出控制是治理該地區(qū)非點源污染物的關(guān)鍵。

圖3　全城塢小流域2014年非點源TP月流失量

污染源 輸出系數(shù)TNTP其他流域文獻(xiàn)參考系數(shù)TNTP負(fù)荷/kgTNTP百分比/%TNTP耕地/(kg·hm-2·a-1)15.870.4619.40～29.00[14-18]0.90～2.26[14-18]210.916.111.111.28林地/(kg·hm-2·a-1)6.330.390.43～11.90[14-18]0.12～1.53[14-18]3800.91234.1820.0149.16城鎮(zhèn)村及工礦用地/(kg·hm-2·a-1)6.270.677.17～17.14[14-15,17-18]0.24～7.33[14-15,17-18]289.4930.931.526.49農(nóng)村生活〔kg/(人·a)〕0.200.100.01～0.04[19-21]0.01～0.06[19-21]205.60102.801.0621.35畜禽養(yǎng)殖〔kg/(頭·a)〕0.830.160.27～4.42[22-24]0.21～2.53[22-24]91.3017.600.483.79降雨輸出————14410.6185.0075.8217.92

注：“—”表示未核算相應(yīng)結(jié)果。

將2014年1—12月每月5日的監(jiān)測結(jié)果代入負(fù)荷核算表達(dá)式(1)對該輸出系數(shù)計算模型進(jìn)行驗證，即為表4中的監(jiān)測值。

將輸出系數(shù)、降雨監(jiān)測數(shù)據(jù)代入公式(3)中獲得表4中計算值，由于輸出系數(shù)即為各斷面負(fù)荷值結(jié)合最優(yōu)化數(shù)學(xué)方法計算獲得，因此各斷面TN，TP計算值即為各斷面污染物總流失負(fù)荷。如表4所示，該輸出系數(shù)計算模型在各斷面平均模擬精度TN，TP分別達(dá)到80.52%以及82.90%。研究表明，該輸出系數(shù)計算模型具有一定的模擬精度，進(jìn)一步說明通過最優(yōu)化數(shù)學(xué)方法在該類型小流域核算得到的各污染源輸出系數(shù)具有一定的參考價值，在下一步試驗中增加采樣頻次以及監(jiān)測斷面有助于提高該方法在該流域中模擬精度。

表4　全城塢流域非點源污染輸出系數(shù)計算模型模擬精度

4結(jié) 論

(3) 基于實測資料以及輸出系數(shù)法計算模型，通過最優(yōu)化數(shù)學(xué)方法核算得到各污染源TN，TP輸出系數(shù)。該輸出模型具有一定的模擬精度，通過驗證得到TN，TP負(fù)荷的模擬精度分別為80.52%以及82.90%。說明此方法在該類型小流域的適用性，核算所得輸出系數(shù)為相似流域提供參考。

(4) 該流域TN輸出以降雨輸出為主占75.82%，其次為林地占20.01%；流域TP輸出中，林地輸出最大，為49.16%，農(nóng)村生活TP輸出占非點源污染輸出的21.35%。因此，為控制非點源污染流失，針對該小流域林地以種植毛竹為主的特點，通過竹腔施肥[25]等技術(shù)實現(xiàn)林地的合理施肥以及對農(nóng)村生活污水進(jìn)行納管收集等措施。

[參考文獻(xiàn)]

[1]Ongley E D, Zhang Xiaolan, Yu Tao. Current status of agricultural and rural non-point source Pollution assessment in China[J]. Environment Pullution, 2010,158(5):1159-1168.

[2]崔健,馬友華,趙艷萍,等.農(nóng)業(yè)面源污染的特性及防治對策[J].中國農(nóng)學(xué)報,2006,22(1):335-340.

[3]Weld J L, Parsons R L, Beelgle D B. Evaluation of phosphorus based nutrient management strategies in Pennsylvania[J]. Journal of Soil and Water Conservation, 2002,57(6):448-454.

[4]Shen Zhenyao, Liao Qian, Hong Qian, et al. An overview of research on agricultural non-point source pollution modelling in China[J]. Separation and Purification Technology, 2012,84(9):104-111.

[5]朱瑤,梁志偉,李偉,等.流域水環(huán)境污染模型及其應(yīng)用研究綜述[J].應(yīng)用生態(tài)學(xué)報,2013,24(10):3012-3018.

[6]黃永剛,付玲玲,胡筱敏.基于河流斷面監(jiān)測資料的非點源負(fù)荷估算輸出系數(shù)法的研究和應(yīng)用[J].水力發(fā)電學(xué)報,2012,31(5):159-162.

[7]Johns P J. Evaluation and management of the impact of land use on the Nitrogen and Phosphorus load delivered to surface water[J]. Journal of Hydrology, 1996,183(3):323-349.

[8]郝旭,張乃明,史靜.昆明市云龍水庫徑流區(qū)氮磷非點源污染負(fù)荷分析[J].水土保持通報,2013,33(6):274-278.

[9]李艷利,徐宗學(xué),李艷粉,等.遼河流域水質(zhì)狀況及其對土地利用/覆被變化的響應(yīng)[J].水土保持通報,2013,33(2):72-77.

[10]張福平,趙沙,周正朝,等.灃河流域土地利用格局與水質(zhì)變化的關(guān)系[J].水土保持通報,2014,34(4):308-312.

[11]杜娟,李懷恩,李家科.基于實測資料的輸出系數(shù)分析與陜西灃河流域非點源負(fù)荷來源探討[J].農(nóng)業(yè)環(huán)境科學(xué)學(xué)報,2013,32(4):827-837.

[12]李懷恩,李家科.流域非點源污染負(fù)荷定量化方法研究與應(yīng)用[M].北京:科學(xué)出版社,2013.

[13]田平.基于GIS杭嘉湖地區(qū)農(nóng)田氮磷徑流流失研究[D].杭州:浙江大學(xué),2006.

[14]廖義善,卓慕寧,李定強(qiáng),等.基于“徑流—地類”參數(shù)的非點源氮磷負(fù)荷估算方法[J].環(huán)境科學(xué)學(xué)報,2014,34(8):2126-2132.

[15]黃金良,洪華生,張珞平,等.基于GIS的九龍江流域農(nóng)業(yè)非點源氮磷負(fù)荷估算研究[J].農(nóng)業(yè)環(huán)境科學(xué)學(xué)報,2004,23(5):866-871.

[16]趙廣舉,田鵬,穆興民,等.基于PCRaster的流域非點源氮磷負(fù)荷估算[J].水科學(xué)進(jìn)展,2012,23(1):80-86.

[17]劉瑞民,楊志峰,丁曉雯,等.土地利用/覆蓋變化對長江上游非點源污染影響研究[J].環(huán)境科學(xué),2006,27(12):2407-2414.

[18]王歡歡,陳世儉.土地利用結(jié)構(gòu)變化對非點源污染的影響研究[J].環(huán)境科學(xué)與技術(shù),2011,34(S2):25-28.

[19]譚平,馬太玲,趙立欣,等.巢湖農(nóng)村生活污水產(chǎn)排污系數(shù)測算及處理模式分析[J].中國給水排水,2012,28(13):88-91.

[20]嚴(yán)婷婷,王紅華,孫治旭,等.滇池流域農(nóng)村生活污水產(chǎn)排污系數(shù)研究[J].環(huán)境科學(xué)導(dǎo)刊,2010,29(4):46-48.

[21]蔡金洲,范先鵬,黃敏,等.湖北省三峽庫區(qū)農(nóng)業(yè)面源污染解析[J].農(nóng)業(yè)環(huán)境科學(xué)學(xué)報,2012,31(7):1421-1430.

[22]蘇文幸.生豬養(yǎng)殖業(yè)主要污染源產(chǎn)排污量核算體系研究[D].長沙:湖南師范大學(xué),2012.

[23]何志平,曾凱,李正確,等.四川規(guī)模豬場產(chǎn)排污系數(shù)測定[J].中國沼氣,2010,28(4):10-14.

[24]汪開英,劉健,陳小霞,等.浙江省畜禽業(yè)產(chǎn)排污測算與土地承載力分析[J].應(yīng)用生態(tài)學(xué)報,2009,20(12):3043-3048.

[25]林自強(qiáng),江建煌.竹腔施肥技術(shù)[J].安徽林業(yè)科技,2011,37(5):79-80.

收稿日期：2015-05-14修回日期：2015-07-13

通訊作者：梁新強(qiáng)(1979—),男(漢族),浙江省杭州市人,博士,副教授,主要從事農(nóng)業(yè)面源污染控制與水質(zhì)評估研究。E-mail:liang410@zju.edu.cn。

文獻(xiàn)標(biāo)識碼：A

文章編號：1000-288X(2016)03-0030-08

中圖分類號：X131.2

Export Characteristics and Source Analysis of Non-point Source Nitrogen and Phosphorus in Mountain Agricultural Small Catchment

XU Lixian1, LIANG Xinqiang1,2, ZHOU Kejin3, LIN Limin1,WANG Zhibo1, HUA Guifen2, WANG Xinyu1, WANG Zhirong4, HUANG Wu4

(1.CollegeofEnvironmentandResourceSciences,ZhejiangUniversity,Hangzhou,Zhejiang310058,China; 2.KeyLaboratoryforWaterPollutionControlandEnvironmentSafetyofZhejiangProvince,Hangzhou,Zhejiang310058,China; 3.EnvironmentalMonitoringCenterofZhejiangProvince,Hangzhou,Zhejiang310015,China; 4.OfficeofAgriculturalEcologyandEnergyofZhejiangProvince,Hangzhou,Zhejiang310012,China)

Abstract：[Objective] To analyze the export characteristics and source of non-point source in order to provide reference basis for the study of the similar agricultural watershed water quality and non-point source pollution control. [Methods] This area was firstly divided into certain catchments with ArcGIS software. The water quality and precipitation was monitored across sections in 2014. Combined with the mean concentration method and export coefficient model, a calculation model was set up to estimate non-point source pollution export , and the sources for various non-point pollutants was analyzed with optimized mathematic method. [Results] The total nitrogen export coefficients of agricultural land, woodland, residential land, domestics and livestock was 15.87 kg/(hm2·a), 6.33 kg/(hm2·a), 6.27 kg/(hm2·a), 0.20 kg/(person·a), 0.83 kg/(anamal·a), respectively, while the total phosphorus export coefficients was 0.46 kg/(hm2·a), 0.39 kg/(hm2·a), 0.67 kg/(hm2·a), 0.10 kg/(person·a), 0.16 kg/(anamal·a), respectively. Nitrate nitrogen was the main form of nitrogen loss in the runoff and a large amount of particulate phosphorus run off in the periods of peak rainfall. In the precipitation, ammonia nitrogen and dissolved phosphorus was main loss form. [Conclusion] In this small catchment, precipitation, woodland and domestic waste are the main pollution sources.

Keywords:non-point source pollution; precipitation; coefficient estimation; nitrogen and phosphorus losses characteristics

資助項目：國家水體污染控制與治理科技重大專項“苕溪流域農(nóng)村污染治理技術(shù)集成與規(guī)模化工程示范”(2014ZX07101-012); 國家自然科學(xué)基金項目“有機(jī)肥施用對稻田土壤膠體磷賦存及運(yùn)移的影響機(jī)理”(41271314)

第一作者：徐麗賢(1991—),女(漢族),浙江省湖州市人,碩士研究生,研究方向為農(nóng)業(yè)面源污染控制與治理。E-mail:21314073@zju.edu.cn。