葉宏萌，袁旭音，林 皓

（武夷學(xué)院 生態(tài)與資源工程學(xué)院，福建省生態(tài)產(chǎn)業(yè)綠色技術(shù)重點(diǎn)實(shí)驗(yàn)室，福建 武夷山 354300）

隨著經(jīng)濟(jì)發(fā)展和人口增加，全球范圍內(nèi)的水污染問題日益加劇。其中，水體氮、磷的超標(biāo)導(dǎo)致的富營養(yǎng)化問題最為嚴(yán)重和普遍[1]，已成為制約社會(huì)經(jīng)濟(jì)發(fā)展的重大環(huán)境問題。例如，歐洲、北美洲、南美洲和亞太地區(qū)發(fā)生富營養(yǎng)化水域的比例分別是53%、48%、41%和54%；而我國的700多條河流（包括長江、黃河和珠江等）中，就有70.6%的河流面臨富營養(yǎng)化的風(fēng)險(xiǎn)[2]。通常而言，氮、磷是控制湖泊或河流中藻類生長的關(guān)鍵因子[3]。隨著點(diǎn)源污染的控制，非點(diǎn)源污染成為了全球水體氮、磷超標(biāo)的主要來源[4-7]，被國際社會(huì)公認(rèn)為人類飲用水源的最大威脅[8]，近幾十年來受到越來越多的關(guān)注[9]。

據(jù)2010年《第一次全國污染源普查公報(bào)》顯示，我國農(nóng)業(yè)非點(diǎn)源污染排放的氮、磷的負(fù)荷量依次為270.46萬t和28.47萬t，分別占氮、磷排放總負(fù)荷量的53.1%和67.3%[7]。因此可以認(rèn)為，農(nóng)業(yè)非點(diǎn)源是我國水體富營養(yǎng)化的最主要來源。同時(shí)，氮磷流失過程是一個(gè)綜合氣候、地形、土壤、降雨、土地利用和管理方式等多種因素的復(fù)雜過程，具有顯著的時(shí)空差異性[10-11]。但是，當(dāng)前對(duì)農(nóng)業(yè)非點(diǎn)源氮磷及綜合流失的研究多集中于空間差異的研究，對(duì)時(shí)間分布差異的關(guān)注尚薄弱。因此，文中較為系統(tǒng)的介紹了農(nóng)業(yè)非點(diǎn)源污染的定義、產(chǎn)生原因和特點(diǎn)，并針對(duì)農(nóng)業(yè)非點(diǎn)源導(dǎo)致的氮磷流失風(fēng)險(xiǎn)進(jìn)行介紹，包括風(fēng)險(xiǎn)評(píng)估方法（辨識(shí)空間差異）和時(shí)間差異性的研究，指出已有研究的不足，以期更好的了解農(nóng)業(yè)非點(diǎn)源污染，對(duì)農(nóng)業(yè)非點(diǎn)源氮磷流失風(fēng)險(xiǎn)的時(shí)空差異評(píng)估和針對(duì)性治理研究提供參考依據(jù)。

1 農(nóng)業(yè)非點(diǎn)源污染概述

1.1 定義

非點(diǎn)源污染通常是在降水和徑流沖刷的作用下，溶解態(tài)或顆粒態(tài)的污染物從非特定范圍的地域通過徑流過程匯入受納水體而引起的污染[12]。其污染源范圍的界定就直接影響污染排放負(fù)荷量的計(jì)算，但由于污染管理水平的差異，不同國家或者地區(qū)，對(duì)其界定范圍往往存在差異[5]。

其中，農(nóng)業(yè)非點(diǎn)源污染是最為普遍和嚴(yán)重的非點(diǎn)源污染，長期以來受人們所關(guān)注。農(nóng)業(yè)非點(diǎn)源污染（又稱農(nóng)業(yè)面源污染）在學(xué)術(shù)界主要有兩種理解[13]，狹義的含義專指種植業(yè)農(nóng)田污染，而廣義上則涵蓋了農(nóng)村生活與農(nóng)業(yè)生產(chǎn)（包括種植業(yè)、畜禽和水產(chǎn)養(yǎng)殖業(yè)）兩個(gè)方面。

1.2 產(chǎn)生原因

導(dǎo)致農(nóng)業(yè)非點(diǎn)源污染的主要原因包括以下幾個(gè)方面：（1）農(nóng)藥、化肥的施用；（2）水土流失污染。其中，土壤侵蝕是規(guī)模最大、危害程度最嚴(yán)重的一種農(nóng)業(yè)非點(diǎn)源污染；（3）農(nóng)村生活污水的未處理或未達(dá)標(biāo)排放；（4）畜禽或水產(chǎn)養(yǎng)殖的糞便排放及剩余餌料等；（5）農(nóng)村生活垃圾的隨意堆放或處理；（6）灌溉污水和（7）大氣的干濕沉降等[5,12]。

1.3 污染特點(diǎn)

由于非點(diǎn)源污染主要由降雨徑流、土壤侵蝕、地表溶質(zhì)溶出和土壤溶質(zhì)滲漏4個(gè)過程形成，不同過程之間相互聯(lián)系、相互作用，導(dǎo)致了農(nóng)業(yè)非點(diǎn)源污染具有來源模糊性、發(fā)生隨機(jī)性、分布廣泛性、機(jī)理復(fù)雜性、時(shí)間潛伏性、污染隱蔽性和難監(jiān)測(cè)、難治理等特點(diǎn)[12-13]。

2 農(nóng)業(yè)非點(diǎn)源氮磷流失風(fēng)險(xiǎn)評(píng)估的研究綜述

2.1 主要氮磷流失風(fēng)險(xiǎn)評(píng)估方法

截止目前，非點(diǎn)源氮磷流失風(fēng)險(xiǎn)計(jì)算的主要方法有野外實(shí)測(cè)法、輸出系數(shù)法、物理模型法和污染指數(shù)法[14]。其中，野外實(shí)測(cè)法需要結(jié)合氣象條件，長期連續(xù)的同步監(jiān)測(cè)降雨量、徑流量及水質(zhì)的過程變化，準(zhǔn)確度較高，但兼具監(jiān)測(cè)工作量大、持續(xù)時(shí)間長且費(fèi)用消耗高等缺點(diǎn)，不便于該法的普及運(yùn)用。輸出系數(shù)法屬于經(jīng)驗(yàn)?zāi)Ｐ?，以?shí)地觀察或?qū)嶒?yàn)數(shù)據(jù)為基礎(chǔ)，具有結(jié)構(gòu)簡單和數(shù)據(jù)獲取容易等特點(diǎn)[15]，但忽略了非點(diǎn)源污染復(fù)雜的遷移轉(zhuǎn)化過程[16]，測(cè)算精度通常較低，更適用于大范圍的污染普查，尤其適用于農(nóng)業(yè)非點(diǎn)源研究初期以及資料缺乏地區(qū)，具有較強(qiáng)的實(shí)用性。物理模型法考慮了污染物的產(chǎn)生和遷移過程，并定量模擬不同污染物的組分含量，計(jì)算結(jié)果準(zhǔn)確可靠，但模型所需數(shù)據(jù)資料繁多、機(jī)理相對(duì)復(fù)雜、運(yùn)行周期長，計(jì)算效率低、對(duì)尺度較為敏感，更適合于中小流域[14,17]。污染指數(shù)法又稱概念模型，是一種半定量化評(píng)估指標(biāo)體系，該法包括源因子和遷移因子兩個(gè)子體系進(jìn)行評(píng)估氮、磷元素從各污染源流失至水體的風(fēng)險(xiǎn)程度。與傳統(tǒng)物理模型法的定量化評(píng)估不同，氮磷污染指數(shù)法的數(shù)據(jù)需求量少且容易獲得、機(jī)理簡單、操作簡便、運(yùn)行周期短、推廣實(shí)用性強(qiáng)，并可準(zhǔn)確評(píng)價(jià)區(qū)域內(nèi)氮磷流失的風(fēng)險(xiǎn)高低，便于篩選關(guān)鍵源區(qū)和對(duì)應(yīng)的采取針對(duì)性治理措施，在大尺度范圍的風(fēng)險(xiǎn)評(píng)估與關(guān)鍵源區(qū)識(shí)別方面比其他方法有著顯著的優(yōu)勢(shì)，近年來在我國得到了廣泛的應(yīng)用[18-19]。以下針對(duì)氮磷污染指數(shù)法進(jìn)行闡述。

氮磷污染指數(shù)法（氮磷指標(biāo)體系法）從評(píng)價(jià)元素的角度可以分為磷指標(biāo)體系法、氮指標(biāo)體系法及氮磷綜合評(píng)估體系[13]。

2.1.1 磷指標(biāo)體系法

最早的磷指標(biāo)體系（PI，Phosphorus Index）法來源于1993年，Lemunyon和Gilbert兩位學(xué)者以土壤本底養(yǎng)分含量、化肥施用量、土壤侵蝕及地表徑流為基本指標(biāo)構(gòu)建了土壤磷流失的敏感性評(píng)價(jià)體系，但只對(duì)基本指標(biāo)的影響因子進(jìn)行統(tǒng)一計(jì)算[20]。20世紀(jì)90年代末，人們開始將影響因子區(qū)分為源因子與遷移因子兩類，并以兩類因子指標(biāo)的乘法運(yùn)算得到最終的磷流失風(fēng)險(xiǎn)[21]，同時(shí)發(fā)現(xiàn)了磷流失嚴(yán)重的區(qū)域主要分布在兩類因子的高風(fēng)險(xiǎn)重疊區(qū)域，從此，開始了PI法對(duì)磷流失的關(guān)鍵源區(qū)識(shí)別的運(yùn)用[22]。2000年，Gburek和Sharpley對(duì)PI法進(jìn)行了第二次修正，引入了距河道有效距離因子（農(nóng)業(yè)區(qū)距離河道距離）[22]，進(jìn)一步提高了該法的準(zhǔn)確度和可操作性，利于普及運(yùn)用[23-26]。

國內(nèi)的學(xué)者也對(duì)PI法進(jìn)行了較為廣泛的改進(jìn)與應(yīng)用。最先引入的是，張淑榮等（2003）應(yīng)用該法估算了于橋水庫流域農(nóng)業(yè)非點(diǎn)源的磷流失風(fēng)險(xiǎn)[27]。后來，王麗華等（2006）根據(jù)我國的實(shí)際情況，增設(shè)了最佳管理措施因子和冰雪凍融因子，修正了磷肥施用量的計(jì)算，加強(qiáng)了PI法在我國的適用性[28]；李琪等（2008）基于PI分級(jí)方案[29]，也提出了修正的磷分級(jí)方案[30]；周惠平等（2008）以巢湖流域?yàn)槔?，引入了土壤磷吸持指?shù)和磷飽和指標(biāo)，提高了風(fēng)險(xiǎn)評(píng)估的準(zhǔn)確性[31]。李文超等（2014）以洱海源頭典型小流域鳳羽河流域?yàn)槔鶕?jù)中國高原特征對(duì)PI法進(jìn)行了簡化和修正，建立了中國高原農(nóng)業(yè)流域磷指數(shù)評(píng)價(jià)體系，并應(yīng)用于研究區(qū)農(nóng)業(yè)非點(diǎn)源磷（溶解態(tài)磷和顆粒態(tài)磷）流失風(fēng)險(xiǎn)的評(píng)價(jià)及關(guān)鍵源區(qū)的識(shí)別[32]。

2.1.2 氮指標(biāo)體系法

隨著PI法的發(fā)展，氮指標(biāo)體系（NLIs，Nitrogen Loss Index System）法也于20世紀(jì)九十年代中期得到了研究和發(fā)展[33]。根據(jù)指標(biāo)體系構(gòu)建重點(diǎn)的差異，NLIs法可以分為三類：源因子型[34-37]、遷移因子型[38-39]與綜合型指標(biāo)[33]。其中，源因子型的NLIs法的研究集中于對(duì)氮源（以氮素平衡指標(biāo)體系最為典型，即研究氮的輸入輸出平衡來得到氮的流失量）的分析構(gòu)建而成，應(yīng)用最為廣泛[35-37]。

就NLIs法在國內(nèi)的研究發(fā)展而言，單獨(dú)采用該法進(jìn)行氮流失風(fēng)險(xiǎn)評(píng)估的研究很少，主要集中于農(nóng)田、林地氮流失負(fù)荷的估算和氮形態(tài)組分流失模型的構(gòu)建[13]。例如，侯彥林等（2008）等監(jiān)測(cè)了農(nóng)田的氮素流失，并提出了農(nóng)田氮淋失過程的“物質(zhì)、動(dòng)力和因素”三環(huán)節(jié)的概念模型，制定了微觀或小區(qū)域非點(diǎn)源氮污染評(píng)價(jià)的若干標(biāo)準(zhǔn)與方法[40]。

2.1.3 氮磷綜合評(píng)估體系

隨著氮、磷污染指數(shù)體系法的發(fā)展，學(xué)者們逐步認(rèn)識(shí)到氮磷流失的差異性及分別單獨(dú)治理的矛盾性。例如，某些減少氮流失的管理措施可能會(huì)增加土壤中磷的含量，進(jìn)而加劇了土壤磷流失；而免耕措施雖然減少了土壤磷流失，卻可能促進(jìn)氮的淋溶，提高氮的流失風(fēng)險(xiǎn)[41]。因此，結(jié)合氮、磷兩元素進(jìn)行的氮磷綜合評(píng)估體系近年來受到了學(xué)者們的關(guān)注和研究[42]。

最早于2000年，Heathwaite基于對(duì)氮、磷單獨(dú)指標(biāo)體系的綜合分析，開始了構(gòu)建氮磷綜合評(píng)估體系的初探，首次提出了綜合源因子與遷移因子的農(nóng)田氮磷流失指標(biāo)體系的概念模型[43]。隨后，該模型體系得到了很快的發(fā)展，較具有影響力的代表學(xué)者有McDowell[44]、Bechmaim等[45]和Drewry等[18]研究者都相繼構(gòu)建了各種類型的綜合評(píng)估體系。

隨著氮磷綜合評(píng)估體系的建立和發(fā)展，該法近幾年在國內(nèi)也得到了較廣泛的應(yīng)用。張平等（2011）以密云水庫沿湖集約化農(nóng)區(qū)東莊小流域?yàn)槔?，發(fā)現(xiàn)該流域氮、磷不同元素流失風(fēng)險(xiǎn)的空間分布差異較大，進(jìn)一步說明了氮磷綜合指標(biāo)的重要性[46]。張汪壽等（2013）以土地利用、徑流量和距離因子等為參評(píng)指標(biāo)，利用改進(jìn)的理想解法確定指標(biāo)權(quán)重，最終建立了多準(zhǔn)則分析評(píng)價(jià)體系，應(yīng)用于北京懷柔區(qū)北宅小流域的農(nóng)業(yè)非點(diǎn)源氮磷流失風(fēng)險(xiǎn)評(píng)估與關(guān)鍵源區(qū)識(shí)別[47]。孫作雷（2015）以土壤侵蝕、地表徑流、土壤質(zhì)地、土壤滲透性及距離河道有效距離為參評(píng)指標(biāo)對(duì)苕溪流域農(nóng)業(yè)非點(diǎn)源氮磷污染風(fēng)險(xiǎn)進(jìn)行了氮、磷及氮磷綜合指標(biāo)的流失風(fēng)險(xiǎn)評(píng)價(jià)[13]。樊慧慧（2016）對(duì)合肥城郊二十辟河流域的某農(nóng)業(yè)匯水區(qū)進(jìn)行氮磷流失風(fēng)險(xiǎn)指數(shù)模型計(jì)算，表明氮磷流失的高風(fēng)險(xiǎn)區(qū)主要分布在距離河流較近的區(qū)域，表明遷移因子中距河道距離的重要性[48]。耿若楠（2017）則對(duì)淮北地區(qū)蚌埠市懷遠(yuǎn)縣某一人工林地土壤的氮磷流失風(fēng)險(xiǎn)進(jìn)行指數(shù)模型計(jì)算，結(jié)果表明該區(qū)人工林氮流失以中等風(fēng)險(xiǎn)為主，磷流失以低風(fēng)險(xiǎn)為主[49]。Ye et al.（2017）對(duì)苕溪流域進(jìn)行了氮、磷及綜合指標(biāo)評(píng)價(jià)，并篩選關(guān)鍵源區(qū)，提出對(duì)應(yīng)非點(diǎn)源污染防治對(duì)策[10]。

綜上所述，國內(nèi)學(xué)者對(duì)污染指數(shù)法的應(yīng)用越來越注重氮磷綜合風(fēng)險(xiǎn)的評(píng)估。但其應(yīng)用多數(shù)局限于農(nóng)田、林地等局部區(qū)域或者小流域農(nóng)業(yè)非點(diǎn)源氮磷綜合流失風(fēng)險(xiǎn)的研究，結(jié)合不同水文期的遷移因子進(jìn)行農(nóng)業(yè)非點(diǎn)源氮磷及綜合流失風(fēng)險(xiǎn)的時(shí)空分布特征的有關(guān)報(bào)道仍屬鮮見[10,18]。

2.2 氮磷流失的時(shí)間差異性

氮、磷流失過程是一個(gè)涉及土壤類型、水文條件、地形地貌、土地利用類型等多種因素的復(fù)雜非點(diǎn)源污染過程。與空間上的關(guān)鍵源區(qū)相似，農(nóng)業(yè)非點(diǎn)源污染負(fù)荷在時(shí)間上也存在一定的變化規(guī)律，與降雨徑流是農(nóng)業(yè)非點(diǎn)源污染輸出、遷移的最主要?jiǎng)恿o密相關(guān)。研究表明，年際或年內(nèi)水文過程的差異對(duì)氮磷流失的影響顯著[10]，其污染負(fù)荷輸出可能會(huì)集中于某一水文年或年內(nèi)某一水文時(shí)期或某一季節(jié)表現(xiàn)得更為嚴(yán)重[50]。

高超等（2005）表明旱地和水稻田用地在極端降水事件（高強(qiáng)度降水）期間隨地表徑流發(fā)生流失的養(yǎng)分含量占全年流失總量的50%以上[51]。張繼宗等（2009）對(duì)太湖平原農(nóng)田區(qū)域徑流氮磷污染監(jiān)測(cè)結(jié)果表明，降雨量越大，歷時(shí)越長，農(nóng)田氮磷流失濃度和流失量越大，因此，汛期的徑流氮磷流失幾率和流失強(qiáng)度均大于非汛期[52]。張皓天（2010）[53]、劉博等（2011）[14]、Francesconi,et al.（2016）[54]、Lee,et al.（2018）[55]利用SWAT（Soil and Water Assessment Tool）模型進(jìn)行不同流域非點(diǎn)源氮磷污染評(píng)估，結(jié)果都表明了，氮磷污染負(fù)荷量與降雨量呈正相關(guān)關(guān)系，汛期是氮磷負(fù)荷輸出的關(guān)鍵期。施練東等（2011）通過應(yīng)用SWAT模型，模擬出湯浦水庫流域農(nóng)業(yè)非點(diǎn)源氮磷污染負(fù)荷的季節(jié)變化規(guī)律呈不規(guī)則的“W”形態(tài)變化，雨季前期（5到6月）和作物收獲后（10月），流域內(nèi)氮磷的流失量會(huì)達(dá)到全年的最大值[56]。張桂軻（2016）以SWAT模型對(duì)長江流域非點(diǎn)源氮磷流失風(fēng)險(xiǎn)進(jìn)行模擬的結(jié)果表明，氮磷輸出負(fù)荷在年內(nèi)分布差異顯著，豐水期總氮和總磷的負(fù)荷強(qiáng)度遠(yuǎn)遠(yuǎn)高于枯水期，5月至9月總氮、總磷負(fù)荷量占全年總負(fù)荷量的比例分別為74%和85%。其中，氮素的輸出受到降雨、徑流動(dòng)力傳輸過程和氮素供給量的共同控制：農(nóng)業(yè)流域氮素供給充足，降雨、徑流過程的控制作用更加凸顯，因此，汛期的氮素濃度高于非汛期；而在非農(nóng)業(yè)流域，受氮素供給量的控制，汛期氮素濃度可能還低于非汛期[57]。

3 小結(jié)

綜上所述，在不同區(qū)域內(nèi)，各因素隨時(shí)空尺度及局部環(huán)境的變化對(duì)非點(diǎn)源氮磷流失產(chǎn)生的影響力亦隨之變化，因此，如何有效地模擬和估算全球各流域氮磷流失風(fēng)險(xiǎn)狀況仍是困擾研究者的一個(gè)現(xiàn)實(shí)難題[5]。同時(shí)，水文過程的差異對(duì)氮磷流失的影響顯著[10,18]，通過合理評(píng)估研究區(qū)域的農(nóng)業(yè)非點(diǎn)源氮磷污染負(fù)荷的時(shí)空分布差異，識(shí)別并篩選氮磷流失發(fā)生嚴(yán)重的關(guān)鍵源區(qū)和關(guān)鍵期，對(duì)于采取時(shí)空針對(duì)性的非點(diǎn)源污染的防治是極為必要的。但有關(guān)非點(diǎn)源氮磷流失關(guān)鍵期篩選的報(bào)道多運(yùn)用物理模型（主要是SWAT模型）進(jìn)行，尚缺乏氮磷綜合流失風(fēng)險(xiǎn)模型的運(yùn)用。而運(yùn)用氮磷及綜合污染指數(shù)法尚局限于關(guān)鍵源區(qū)的篩選，并多集中于農(nóng)業(yè)小流域或者局部的林業(yè)區(qū)，缺乏關(guān)鍵期的篩選[31]。因此，亟待結(jié)合不同水文期的差異加強(qiáng)氮磷及綜合流失風(fēng)險(xiǎn)的時(shí)空差異研究。