王思如,楊大文,孫金華,唐莉華,王朋杰，盧少勇

(1.清華大學(xué)土木水利學(xué)院，北京 100084； 2.清華大學(xué)水沙科學(xué)與水利水電工程國(guó)家重點(diǎn)實(shí)驗(yàn)室,北京 100084；3.南京水利科學(xué)研究院水文水資源與水利工程科學(xué)國(guó)家重點(diǎn)實(shí)驗(yàn)室,江蘇 南京 210029；4.中國(guó)環(huán)境科學(xué)研究院環(huán)境基準(zhǔn)與風(fēng)險(xiǎn)評(píng)估國(guó)家重點(diǎn)實(shí)驗(yàn)室, 北京 100012；5.中國(guó)環(huán)境科學(xué)研究院國(guó)家環(huán)境保護(hù)湖泊污染控制重點(diǎn)實(shí)驗(yàn)室，北京 100012；6.中國(guó)環(huán)境科學(xué)研究院國(guó)家環(huán)境保護(hù)洞庭湖科學(xué)觀測(cè)研究站,北京 100012)

水污染嚴(yán)重影響飲用水、糧食、生態(tài)和人居環(huán)境的安全[1-2]，成為我國(guó)社會(huì)經(jīng)濟(jì)可持續(xù)發(fā)展的瓶頸之一[3]。隨著工業(yè)和城市生活點(diǎn)源污染控制措施逐步取得成效，范圍更廣、治理難度更大的農(nóng)業(yè)面源污染(agricultural non-point source pollution， ANP)成為我國(guó)水污染的另一重要來(lái)源[4-7]。ANP是指農(nóng)業(yè)生產(chǎn)過(guò)程中，農(nóng)田化肥流失、畜禽及水產(chǎn)養(yǎng)殖廢水排放等導(dǎo)致污染物進(jìn)入水體造成的污染[8-9]。在太湖流域污染負(fù)荷中，83%的總氮和84%的總磷來(lái)自ANP[10-11]；ANP對(duì)洞庭湖總氮和總磷污染貢獻(xiàn)率分別達(dá)到61%和80%[12]。因此，掌握我國(guó)ANP現(xiàn)狀與特征對(duì)指導(dǎo)水資源保護(hù)與水環(huán)境治理十分必要。

相關(guān)研究通常采用Johnes[13]提出的輸出系數(shù)法估算ANP排放量。該方法假定單位數(shù)量(單位面積)同一污染源的輸出系數(shù)相同，相比于考慮流域內(nèi)污染物形成、運(yùn)移和轉(zhuǎn)化過(guò)程的機(jī)理模型[14]，輸出系數(shù)法簡(jiǎn)單方便，不依賴大量觀測(cè)數(shù)據(jù)來(lái)率定模型參數(shù)。當(dāng)觀測(cè)數(shù)據(jù)有限時(shí)，輸出系數(shù)法是可行的方法，因此得到了廣泛使用[15-16]。我國(guó)針對(duì)ANP排放現(xiàn)狀的研究側(cè)重于對(duì)局部區(qū)域的研究[17-20]，在已有少數(shù)針對(duì)全國(guó)的研究中[21-22]，沒(méi)有考慮污染物輸出系數(shù)在不同區(qū)域的差異。此外，國(guó)內(nèi)學(xué)者在運(yùn)用該方法時(shí)極少考慮水產(chǎn)養(yǎng)殖污染排放，我國(guó)有71.20%的水產(chǎn)品來(lái)自于池塘養(yǎng)殖[9]，有必要在計(jì)算ANP排放量時(shí)考慮水產(chǎn)養(yǎng)殖污染。

本文基于2007年全國(guó)第一次污染源普查資料，考慮輸出系數(shù)在不同區(qū)域的差異性，根據(jù)2016年農(nóng)田種植、畜禽養(yǎng)殖和水產(chǎn)養(yǎng)殖統(tǒng)計(jì)數(shù)據(jù)，計(jì)算了2016年我國(guó)31個(gè)省、自治區(qū)、直轄市(以下簡(jiǎn)稱省市，港澳臺(tái)地區(qū)數(shù)據(jù)缺失)ANP中的總氮、總磷排放量；采用單位面積污染負(fù)荷系數(shù)和等標(biāo)污染負(fù)荷指數(shù)評(píng)價(jià)了各省、市環(huán)境影響程度與污染嚴(yán)重程度，分析了ANP的空間特征及不同類型ANP比例，指出ANP的關(guān)鍵污染源，并利用SPSS聚類分析識(shí)別各省市污染源特征。

1 研究方法與數(shù)據(jù)來(lái)源

1.1 研究方法

1.1.1污染排放量計(jì)算

計(jì)算ANP排放量采用改進(jìn)輸出系數(shù)法[23]，即在原輸出系數(shù)法基礎(chǔ)上考慮不同類型污染源輸出系數(shù)的差異。以省市為計(jì)算單元，計(jì)算農(nóng)田種植、畜禽養(yǎng)殖和水產(chǎn)養(yǎng)殖的氮磷排放量。計(jì)算公式為

TAL=FL+LL+AL

(1)

其中

FL=10-3Sαf+B

式中：TAL為ANP年總排放量，t；FL為農(nóng)田種植氮磷污染年排放量，t；LL為畜禽養(yǎng)殖氮磷污染年排放量，t；AL為水產(chǎn)養(yǎng)殖氮磷污染年排放量，t；S為各省市農(nóng)作物氮磷肥施用折純量，t；αf為農(nóng)田種植排放/流失系數(shù)，g/kg；B為農(nóng)田種植氮磷基礎(chǔ)流失量，t；i為某種養(yǎng)殖類型序號(hào)，如i=1，2，3，4分別為豬、牛、羊、禽；n為畜禽或水產(chǎn)養(yǎng)殖類型數(shù)目；Hi為不同畜禽養(yǎng)殖類型牲畜出欄量，萬(wàn)頭；αli為i類型畜禽養(yǎng)殖氮磷排污系數(shù)，kg/(頭·a)；Yi為不同水產(chǎn)養(yǎng)殖類型養(yǎng)殖產(chǎn)量，t；αai為i類型水產(chǎn)養(yǎng)殖氮磷排污系數(shù)，g/kg。

假設(shè)農(nóng)田種植基礎(chǔ)流失量隨時(shí)間不變，肥料輸出系數(shù)等于《第一次污染源普查數(shù)據(jù)集》(以下簡(jiǎn)稱“一普”)肥料輸出量除以相應(yīng)年份施肥量；不同類型畜禽養(yǎng)殖的污染輸出系數(shù)依據(jù)文獻(xiàn)[15-16,19,24-44]并結(jié)合一普的畜禽養(yǎng)殖污染排放量與相應(yīng)年份各類型畜禽養(yǎng)殖規(guī)模進(jìn)行修正后確定；水產(chǎn)養(yǎng)殖輸出系數(shù)依據(jù)一普附件《水產(chǎn)養(yǎng)殖業(yè)污染源產(chǎn)排污系數(shù)手冊(cè)》中現(xiàn)場(chǎng)試驗(yàn)數(shù)據(jù)，并根據(jù)各省市不同養(yǎng)殖模式比例，綜合考慮該類養(yǎng)殖品種養(yǎng)殖模式可能性，計(jì)算其采取某種養(yǎng)殖模式的權(quán)重，得到該種類型水產(chǎn)養(yǎng)殖的綜合輸出系數(shù)，并根據(jù)一普的水產(chǎn)養(yǎng)殖污染排放量與相應(yīng)年份各類型水產(chǎn)養(yǎng)殖規(guī)模對(duì)輸出系數(shù)進(jìn)行修正。本文暫不考慮年降雨差異和流域損失對(duì)污染輸出系數(shù)的影響。計(jì)算過(guò)程中并未考慮輸出系數(shù)隨時(shí)間發(fā)生變化，可能會(huì)對(duì)污染排放量存在一定程度的高估[45]。

1.1.2污染程度評(píng)價(jià)

a.單位面積污染負(fù)荷系數(shù)法。采用單位面積污染負(fù)荷系數(shù)評(píng)價(jià)污染物對(duì)某省市的環(huán)境影響程度[46]。計(jì)算公式為

(2)

其中

式中：K為單位面積污染負(fù)荷系數(shù)；qi為某省市第i種類型污染的單位面積年排放量，t/km2；Qi為某省市第i種類型污染的年排放量，t；Si為某省市行政面積，km2。

根據(jù)K值評(píng)價(jià)總氮、總磷對(duì)各省市環(huán)境影響程度[47]：K<60，對(duì)環(huán)境不構(gòu)成威脅(Ⅰ級(jí))；60≤K<80，對(duì)環(huán)境稍有威脅(Ⅱ級(jí))；80≤K<100，對(duì)環(huán)境構(gòu)成威脅(Ⅲ級(jí))；K≥100，對(duì)環(huán)境構(gòu)成嚴(yán)重威脅(Ⅳ級(jí))。

b.等標(biāo)污染負(fù)荷指數(shù)法。采用等標(biāo)污染負(fù)荷指數(shù)來(lái)反映各省市相對(duì)污染嚴(yán)重程度，即某省市等標(biāo)污染負(fù)荷總量除以其水資源總量[48]，等標(biāo)污染負(fù)荷總量為把污染物全部稀釋到評(píng)價(jià)標(biāo)準(zhǔn)所需的介質(zhì)量[49]，表征污染物排放濃度超過(guò)某一基準(zhǔn)排放標(biāo)準(zhǔn)的程度。計(jì)算公式為

(3)

其中

式中:λij為第j種污染源中第i種污染物的等標(biāo)污染指數(shù)；Pij為第j種污染源中第i種污染物的等標(biāo)污染負(fù)荷總量，m3；Wi為該污染物所處省市的水資源總量，m3；Cij為第j種污染源中第i種污染物的年排放量，t；ρi0為GB 3838—2002《地表水環(huán)境質(zhì)量標(biāo)準(zhǔn)》中第i種污染物的水質(zhì)標(biāo)準(zhǔn)，mg/L。本文按照水質(zhì)標(biāo)準(zhǔn)為Ⅲ類進(jìn)行計(jì)算，總氮為1 mg/L，總磷為 0.2 mg/L。

1.1.3污染特征識(shí)別

采用SPSS軟件將農(nóng)田種植、畜禽養(yǎng)殖與水產(chǎn)養(yǎng)殖的氮磷等標(biāo)污染負(fù)荷總量作為輸入，利用分層聚類方法對(duì)污染排放特征進(jìn)行聚類分析。ANP變化的主導(dǎo)因素為對(duì)ANP氮磷排放變化貢獻(xiàn)率最大的行業(yè)，貢獻(xiàn)率按各行業(yè)氮磷排放量值占ANP氮磷總排放量值的百分比進(jìn)行計(jì)算。

1.2 數(shù)據(jù)來(lái)源

我國(guó)各省市農(nóng)林牧漁業(yè)總產(chǎn)值數(shù)據(jù)來(lái)源于《中國(guó)農(nóng)業(yè)統(tǒng)計(jì)資料》；農(nóng)作物氮肥、磷肥施用折純量和畜禽養(yǎng)殖出欄量均源自《中國(guó)農(nóng)業(yè)年鑒》；水產(chǎn)養(yǎng)殖產(chǎn)量數(shù)據(jù)來(lái)自《中國(guó)漁業(yè)年鑒》；污染輸出系數(shù)計(jì)算數(shù)據(jù)來(lái)自《第一次全國(guó)污染源普查數(shù)據(jù)集》及其附件《水產(chǎn)養(yǎng)殖業(yè)污染源產(chǎn)排污系數(shù)手冊(cè)》；水資源總量數(shù)據(jù)來(lái)源于各年《中國(guó)水資源公報(bào)》。

2 ANP排放量及特征分析

2.1 污染排放量

由表1可見(jiàn)，2016年，我國(guó)ANP總氮排放量為294.3萬(wàn)t，其中，農(nóng)田種植、畜禽養(yǎng)殖、水產(chǎn)養(yǎng)殖排放分別占53.8%、41.9%和4.3%。全國(guó)總氮排放量各省市平均值為9.5萬(wàn)t，高于該值的省市有13個(gè)，均位列該年農(nóng)林牧漁業(yè)總產(chǎn)值的前13位，說(shuō)明地區(qū)的ANP總氮排放量與當(dāng)?shù)剞r(nóng)業(yè)發(fā)展水平直接相關(guān)。山東、河南、天津、江蘇4省市的總氮排放強(qiáng)度均高于全國(guó)平均值的2倍，表明上述省市在產(chǎn)業(yè)發(fā)展過(guò)程中需更加注重環(huán)境承載能力與可持續(xù)發(fā)展。由表1中各行業(yè)總氮排放貢獻(xiàn)率可知，不同省市的ANP總氮排放源差異較大。13個(gè)總氮排放量較大的省市中，山東、河北、湖北、江蘇、安徽、廣西、四川7省市由農(nóng)田種植主導(dǎo)，河南、廣東、湖南、福建、黑龍江、遼寧6省市由畜禽養(yǎng)殖主導(dǎo)，這主要與各省市的農(nóng)業(yè)生產(chǎn)結(jié)構(gòu)有重要關(guān)系。

由表1可見(jiàn)，2016年，我國(guó)ANP總磷排放量為33.1萬(wàn)t。其中，農(nóng)田種植、畜禽養(yǎng)殖、水產(chǎn)養(yǎng)殖排放分別占32.7%、60.1和7.2%。全國(guó)ANP總磷排放量各省市平均值為1.1萬(wàn)t，高于該值的省市有12個(gè)，分別是山東、河南、廣東、福建、江蘇、湖北、湖南、河北、遼寧、廣西、四川和安徽。上述省市總氮排放量也高于全國(guó)平均值，說(shuō)明不同省市總磷與總氮排放量排序較為一致(相關(guān)系數(shù)為0.97)，表明農(nóng)業(yè)發(fā)展對(duì)氮、磷排放的貢獻(xiàn)同步。此外，除福建外11個(gè)省市2016年農(nóng)林牧漁業(yè)總產(chǎn)值均位列全國(guó)前12位，說(shuō)明地區(qū)的ANP總磷排放量也與當(dāng)?shù)氐霓r(nóng)業(yè)發(fā)展水平直接相關(guān)。山東、河南、天津、江蘇、福建5省市的總磷排放強(qiáng)度均高于全國(guó)平均值的2倍，此結(jié)論與總氮類似，說(shuō)明了我國(guó)ANP氮磷排放問(wèn)題的一致性。由表1中的各行業(yè)總磷排放貢獻(xiàn)率可知，12個(gè)總磷排放量較大的省市中，山東、河南、廣東、福建、江蘇、遼寧等10省市主要來(lái)自畜禽養(yǎng)殖，尤其是遼寧省，其畜禽養(yǎng)殖貢獻(xiàn)率占到96%以上，廣西和四川2省主要來(lái)自農(nóng)田種植。

表1 2016年全國(guó)各省市ANP排放及各行業(yè)貢獻(xiàn)率

2.2 污染程度分析

2.2.1單位面積污染負(fù)荷系數(shù)分析

根據(jù)2016年全國(guó)各省市單位面積負(fù)荷系數(shù)分析結(jié)果，除了重慶和北京以外，其余各省市總氮和總磷的分級(jí)空間分布基本一致(表2)。東部、中部較西部地區(qū)的環(huán)境威脅更嚴(yán)重，而前者農(nóng)業(yè)開(kāi)發(fā)強(qiáng)度較后者高，導(dǎo)致前者單位面積污染負(fù)荷系數(shù)(污染程度)更高。氮磷排放量在全國(guó)平均值以上的省市中，只有黑龍江和四川對(duì)環(huán)境不構(gòu)成威脅，廣西對(duì)環(huán)境稍有威脅，其余省市均對(duì)環(huán)境構(gòu)成嚴(yán)重威脅。

表2 2016年全國(guó)各省市ANP單位面積污染負(fù)荷系數(shù)與等標(biāo)污染負(fù)荷指數(shù)

2.2.2等標(biāo)污染負(fù)荷指數(shù)分析

2016年，全國(guó)各省市氮磷等標(biāo)污染負(fù)荷指數(shù)與污染排放量的排序結(jié)果截然不同，氮磷等標(biāo)污染負(fù)荷指數(shù)在各省市間差異很大，且各省市總磷與總氮的等標(biāo)污染負(fù)荷指數(shù)排序具有一致性(相關(guān)系數(shù)為0.97)。對(duì)于總氮和總磷等標(biāo)污染負(fù)荷指數(shù)，超過(guò)全國(guó)平均值的有寧夏、山東、天津、河南、河北、山西和遼寧，上述省市均值為其他24個(gè)省市均值的10倍以上(表2)。寧夏的氮磷排放量和排放強(qiáng)度均顯著低于全國(guó)平均值，其等標(biāo)污染負(fù)荷指數(shù)卻均位于全國(guó)前五，與其水資源總量小、水環(huán)境容量低直接相關(guān)。廣東、湖北、湖南、江蘇、福建的氮磷排放量與排放強(qiáng)度均高于全國(guó)平均值，其等標(biāo)污染指數(shù)均顯著低于全國(guó)平均值，可見(jiàn)豐沛的水資源量對(duì)于減輕污染程度具有至關(guān)重要的作用。

2.3 污染特征識(shí)別

根據(jù)污染排放特征聚類分析結(jié)果(圖1)，將全國(guó)各省市劃分為4類：①農(nóng)田種植對(duì)ANP貢獻(xiàn)占絕對(duì)主導(dǎo)的省市，包括貴州、云南和西藏，其農(nóng)田種植等標(biāo)污染負(fù)荷占比分別為77.0%、84.0%和98.8%；②畜禽養(yǎng)殖對(duì)ANP貢獻(xiàn)占絕對(duì)主導(dǎo)的省市，包括福建和遼寧，其畜禽養(yǎng)殖等標(biāo)污染負(fù)荷占比分別為69.6%和82.8%；③農(nóng)田種植和畜禽養(yǎng)殖對(duì)ANP貢獻(xiàn)為主要因素，且農(nóng)田種植占比更大的省市，包括山西、甘肅、新疆、廣西、北京、上海、浙江、河北、重慶、四川、安徽、陜西和海南；④農(nóng)田種植和畜禽養(yǎng)殖對(duì)ANP貢獻(xiàn)為主要因素，且畜禽養(yǎng)殖占比更大的省市，包括山東、河南、吉林、青海、黑龍江、寧夏、湖南、江蘇、江西、內(nèi)蒙古、廣東、天津、湖北。我國(guó)ANP污染源主要是畜禽養(yǎng)殖和農(nóng)田種植，雖然在各省市污染源占比上不盡相同，但是兩者之和占比大多數(shù)超過(guò)90%，可見(jiàn)，加強(qiáng)畜禽養(yǎng)殖行業(yè)管理與污染治理、生態(tài)種植與水土保持等措施對(duì)解決我國(guó)農(nóng)業(yè)面源污染問(wèn)題至關(guān)重要。

圖1 等標(biāo)污染負(fù)荷量聚類樹(shù)形圖

3 ANP排放量變化及特征分析

3.1 污染排放變化分析

2007—2016年，全國(guó)ANP總氮排放量增長(zhǎng)了23.8萬(wàn)t，占2007年排放總量的8.8%(表3)。總氮排放增長(zhǎng)的省市有26個(gè)，其中絕對(duì)量增長(zhǎng)最大的5個(gè)省市依次為福建、黑龍江、湖北、河南和遼寧，增長(zhǎng)幅度最大的5個(gè)省市依次為福建、寧夏、新疆、黑龍江和內(nèi)蒙古，排放減幅前5位依次為北京、上海、浙江、河北、山東?？偟欧旁鲩L(zhǎng)絕對(duì)量或幅度位于全國(guó)前5位的省市，其農(nóng)林牧漁業(yè)總產(chǎn)值的增幅均超過(guò)了100%，說(shuō)明農(nóng)業(yè)發(fā)展水平提升會(huì)造成更多總氮排放?？偟欧艤p少的省市中，上海、北京和山東的農(nóng)林牧漁業(yè)總產(chǎn)值增幅不足100%，但浙江和河北的總產(chǎn)值增幅在全國(guó)位列前2位，這兩個(gè)以農(nóng)田種植為ANP主要來(lái)源的省市在大幅提升農(nóng)業(yè)發(fā)展水平的同時(shí)能降低總氮排放量，可能是其農(nóng)業(yè)生產(chǎn)技術(shù)得以革新，提高了農(nóng)田氮素利用效率。

表3 2007年至2016年全國(guó)各省市ANP排放量變化及其占比情況

2007—2016年，全國(guó)有28個(gè)省市的ANP總磷排放量出現(xiàn)了增長(zhǎng)，有3個(gè)省市出現(xiàn)了下降。總體上，總磷排放量增長(zhǎng)了4.7萬(wàn)t，相比于2007年的排放總量上升了16.4%(表3)。其中，總磷排放絕對(duì)量增長(zhǎng)最大的5個(gè)省市依次為山東、福建、河南、湖北和遼寧，總磷增幅最大的5個(gè)省市依次為福建、新疆、黑龍江、內(nèi)蒙古和湖北，總磷排放減幅前3位依次為北京、上海、浙江。福建、河南、湖北和遼寧4省市總氮和總磷排放的絕對(duì)量增長(zhǎng)均位列全國(guó)前5位，而福建、新疆、黑龍江和內(nèi)蒙古4省市的總氮和總磷排放量增幅也位列全國(guó)前5位，說(shuō)明大部分省市農(nóng)業(yè)生產(chǎn)中的總氮和總磷排放量的變化同步。然而，相比于同期總氮排放量的變化，總磷排放增長(zhǎng)的省市更多，其中山東和河北雖然總氮排放量下降，但其總磷排放量上升，體現(xiàn)了這兩省農(nóng)業(yè)發(fā)展中不同行業(yè)發(fā)展變化方向的差異性。

3.2 污染排放變化主導(dǎo)因素

根據(jù)2007—2016年各行業(yè)氮磷排放量變化情況，計(jì)算得到導(dǎo)致總氮、總磷排放量變化的主導(dǎo)因素(表4)。農(nóng)田種植、畜禽養(yǎng)殖和水產(chǎn)養(yǎng)殖變化對(duì)總氮排放增長(zhǎng)的貢獻(xiàn)率分別為42.5%、44.6%和12.9%。83.8%的省市總氮排放量增長(zhǎng)，其中，69.2%的省市由畜禽養(yǎng)殖排污增長(zhǎng)主導(dǎo)，23.1%的省市(吉林、陜西、海南、云南、新疆、西藏)由農(nóng)田種植排污增長(zhǎng)主導(dǎo)，其余7.7%的省市(貴州、廣東)由水產(chǎn)養(yǎng)殖排污增長(zhǎng)主導(dǎo)。16.2%的省市總氮排放量降低，其中40.0%的省市(山東、河北)由農(nóng)田種植排污減少主導(dǎo)，60.0%的省市(浙江、北京、上海)由畜禽養(yǎng)殖排污減少主導(dǎo)。農(nóng)田種植、畜禽養(yǎng)殖和水產(chǎn)養(yǎng)殖對(duì)總磷排放增長(zhǎng)的貢獻(xiàn)率分別為9.4%、71.7%和18.9%。90.3%的省市總磷排放量增長(zhǎng)，其中，67.8%由畜禽養(yǎng)殖排污增長(zhǎng)主導(dǎo)，17.9%的省市(陜西、海南、云南、新疆、西藏)由農(nóng)田種植排污增長(zhǎng)主導(dǎo)，其余14.3%(廣東、湖北、貴州、寧夏)由水產(chǎn)養(yǎng)殖排污增長(zhǎng)主導(dǎo)。9.7%的省市總磷排放量降低，其中，33.3%(上海)由農(nóng)田種植排污減少主導(dǎo)，66.7%(浙江、北京)由畜禽養(yǎng)殖排污減少主導(dǎo)。綜上所述，盡管農(nóng)田種植和水產(chǎn)養(yǎng)殖存在一定影響，2007—2016年影響我國(guó)ANP氮磷排放的核心污染源是畜禽養(yǎng)殖，因此，抑制ANP增長(zhǎng)應(yīng)著重考慮削減畜禽養(yǎng)殖污染。

表4 2007年至2016年全國(guó)各省市ANP變化主導(dǎo)因素

4 結(jié) 語(yǔ)

我國(guó)2016年ANP總氮、總磷排放量分別為294.3萬(wàn)t、33.1萬(wàn)t，農(nóng)田種植和畜禽養(yǎng)殖污染是ANP排放的主要來(lái)源。ANP排放對(duì)我國(guó)東部、中部的環(huán)境影響較西部地區(qū)更嚴(yán)重。受各省市的水資源自然稟賦影響，對(duì)水資源匱乏的省市而言，雖然ANP排放總量與強(qiáng)度不高，但是水環(huán)境污染程度仍然較高。相較于2007年，2016年全國(guó)ANP排放的總氮、總磷分別增長(zhǎng)了8.8%和16.4%，增長(zhǎng)量主要來(lái)自畜禽養(yǎng)殖的快速發(fā)展。抑制我國(guó)ANP增長(zhǎng)，應(yīng)加強(qiáng)對(duì)畜禽養(yǎng)殖的管理，尤其是東部與中部地區(qū)的畜禽養(yǎng)殖管理。