鄉(xiāng)村振興背景下廣東省農(nóng)業(yè)面源污染時空分析與防治對策

張明浩

（廣東國地規(guī)劃科技股份有限公司，廣州 510640）

2018年《中共中央國務(wù)院關(guān)于實施鄉(xiāng)村振興戰(zhàn)略的意見》指出，實施鄉(xiāng)村振興戰(zhàn)略，讓農(nóng)業(yè)成為有奔頭的產(chǎn)業(yè)，讓農(nóng)民成為有吸引力的職業(yè)，讓農(nóng)村成為安居樂業(yè)的美麗家園。中國農(nóng)村地域廣闊，農(nóng)村生態(tài)環(huán)境質(zhì)量直接關(guān)系到人民群眾的身體健康，關(guān)系到鄉(xiāng)村振興戰(zhàn)略的具體落實，關(guān)系到“美麗中國”的建設(shè)。然而，農(nóng)業(yè)面源污染已成為中國水體污染、土壤環(huán)境惡化、生態(tài)系統(tǒng)破壞等的最主要來源之一［1］。全國各污染物排放總量巨大，來自農(nóng)業(yè)污染源排放的化學(xué)需氧量（COD）、總磷（TP）、總氮（TN）的貢獻(xiàn)率分別達(dá)43.71%、67.27%和57.19%［2］，可見農(nóng)業(yè)面源污染已成為影響農(nóng)村增綠、農(nóng)業(yè)增效和農(nóng)民增收、制約鄉(xiāng)村振興戰(zhàn)略實施的關(guān)鍵因素［3］。

當(dāng)前，學(xué)術(shù)界從不同角度對農(nóng)業(yè)面源污染進(jìn)行了大量的研究。其內(nèi)容涵蓋農(nóng)業(yè)面源污染的形成機制［4］、驅(qū)動力［5］、防治政策［6，7］、與多主體行為關(guān)系研究［8-10］、污染現(xiàn)狀分布［11］與時空分析［12，13］等多方面。從已有研究看，采用不同的方法對農(nóng)業(yè)面源污染現(xiàn)狀進(jìn)行測度與分析是該領(lǐng)域研究的熱點問題［14］，但不難發(fā)現(xiàn)，當(dāng)前對農(nóng)業(yè)面源污染時空格局的研究多集中于中國中部、西部、北部地區(qū)或者全國整體層面，較少對廣東省的農(nóng)業(yè)面源污染空間特征以及形成原因等的深度分析。廣東省作為全國農(nóng)業(yè)大省之一，也是全國首個利用世行貸款實施農(nóng)業(yè)面源污染治理的省。近年來，廣東省農(nóng)業(yè)產(chǎn)業(yè)發(fā)展迅速，2010—2017年農(nóng)業(yè)總產(chǎn)值連續(xù)8年保持高速增長，農(nóng)業(yè)增效與農(nóng)民增收取得明顯效果。然而，隨著種植業(yè)、畜牧業(yè)、水產(chǎn)養(yǎng)殖業(yè)等的不斷發(fā)展，由此帶來的農(nóng)業(yè)面源污染日益增大，進(jìn)一步導(dǎo)致水污染、土壤污染等環(huán)境壓力日益增加［15］。據(jù)研究，2012年廣東省生豬與家禽等養(yǎng)殖污染排放量居全國第七位，而且全省47.4%的中小規(guī)模畜禽養(yǎng)殖污染得不到有效控制［16］。鄉(xiāng)村振興戰(zhàn)略明確提出加強農(nóng)業(yè)面源污染防治，開展農(nóng)業(yè)綠色發(fā)展行動，推進(jìn)農(nóng)村生態(tài)文明建設(shè)。因此，在鄉(xiāng)村振興背景下如何在保持農(nóng)業(yè)增效、農(nóng)民增收的基礎(chǔ)上實現(xiàn)農(nóng)村增綠，是廣東省農(nóng)村生態(tài)環(huán)境治理面臨的重要問題。

鑒于此，本研究估算廣東省2008—2017年農(nóng)業(yè)面源污染的COD、TN、TP污染負(fù)荷，并結(jié)合ArcGIS10.2可視化分析探討其時空分布，同時通過各污染物的地均排放量及其綜合指數(shù)對廣東省21個地市農(nóng)業(yè)面源污染強度進(jìn)行風(fēng)險分析并提出針對性措施。

1 材料與方法

1.1 研究區(qū)概況

廣東省位于北緯20°09′—25°31′、東經(jīng)109°45′—117°20′，地處中國東南沿海，交通條件便利，地理區(qū)位優(yōu)勢明顯。2017年廣東省實現(xiàn)地區(qū)生產(chǎn)總值89 879.23億元，比上年增長7.5%［17］。其中，第一產(chǎn)業(yè)增加值為3 792.40億元，增長3.5%。雖然2017年糧食作物播種面積為250萬hm2，比上年下降0.4%，但全年糧食總產(chǎn)量為1 365.13萬t，增長0.4%，這與廣東省農(nóng)業(yè)面源污染治理項目實施有一定關(guān)系。2017年全省農(nóng)村常住人口3 367.45萬人，比上年減少了20.24萬人，占年末常住人口的43.16%［18］。2017年畜禽養(yǎng)殖總產(chǎn)量總體呈下降趨勢。其中，豬肉產(chǎn)量262.19萬t，下降0.8%；禽肉產(chǎn)量134.27萬t，下降0.6%。全年水產(chǎn)品產(chǎn)量增長，其中海水產(chǎn)品增長1.6%，淡水產(chǎn)品增長1.3%。2017年化肥施用量（折純）為237.94萬t，農(nóng)藥施用量為9.46萬t［19］，對當(dāng)?shù)厣鷳B(tài)環(huán)境產(chǎn)生一定影響。總而言之，農(nóng)業(yè)面源污染與種植業(yè)、畜禽養(yǎng)殖業(yè)、水產(chǎn)養(yǎng)殖業(yè)的發(fā)展密不可分。

1.2 數(shù)據(jù)來源與研究方法

1.2.1 數(shù)據(jù)來源 本研究采用的關(guān)于廣東省21個地級市農(nóng)業(yè)面源污染中各種污染源（如化肥施用量、畜禽養(yǎng)殖、水產(chǎn)養(yǎng)殖等）的農(nóng)業(yè)面源污染物指標(biāo)等所有原始數(shù)據(jù)均來自2008—2017年的廣東省統(tǒng)計年鑒、廣東省統(tǒng)計公報、各地市的統(tǒng)計年鑒與統(tǒng)計公報等。廣東省市縣級行政區(qū)劃矢量數(shù)據(jù)由國家基礎(chǔ)地理信息中心提供。

1.2.2 研究方法

1）污染負(fù)荷估算。獲取全省種植業(yè)、畜禽養(yǎng)殖和水產(chǎn)養(yǎng)殖污染物排放總量，采用排污系數(shù)法［20，21］計算21個地市農(nóng)村畜禽養(yǎng)殖、種植業(yè)、水產(chǎn)養(yǎng)殖污染負(fù)荷。

式中，Pij為i污染源的污染物j在區(qū)域的總負(fù)荷量；j為區(qū)域中的污染物類型，包括COD、TN、TP共3種；mij為污染物排放系數(shù)；i為區(qū)域中的污染源類型，包括農(nóng)村種植業(yè)、畜禽養(yǎng)殖業(yè)、水產(chǎn)養(yǎng)殖業(yè)；N為區(qū)域中種植業(yè)化肥施用折純量、畜禽養(yǎng)殖總量、水產(chǎn)養(yǎng)殖總量［22］。從農(nóng)業(yè)面源污染物排放量的核算來看，主要核算污染物包括TN、TP、COD。根據(jù)現(xiàn)有文獻(xiàn)，這3類污染物中農(nóng)藥和農(nóng)膜的貢獻(xiàn)率很低，因此本研究僅選取化肥施用量。

從廣東省第一次全國污染源普查公報中農(nóng)業(yè)污染源排放情況的種植業(yè)、畜禽養(yǎng)殖業(yè)、水產(chǎn)養(yǎng)殖業(yè)的污染物排放數(shù)據(jù)獲取COD、TN和TP年總排放量，從廣東省2008年統(tǒng)計年鑒獲取全省畜禽養(yǎng)殖量、種植業(yè)化肥施用折純量及水產(chǎn)養(yǎng)殖量，將對應(yīng)污染物排放量/對應(yīng)施用量或養(yǎng)殖量作為排放系數(shù)，其中農(nóng)村種植業(yè)沒有對應(yīng)的COD年排放量統(tǒng)計數(shù)據(jù)（表1）。

表1 廣東省農(nóng)業(yè)面源污染物排放系數(shù)

2）污染風(fēng)險評估。由于不同區(qū)域污染物的排放不同，農(nóng)村土地面積也不一樣，無法在同一尺度上來比較21個地級市COD、TN、TP污染物的危害量比關(guān)系，因此在分析廣東省農(nóng)業(yè)面源污染強度時引入地均污染強度。用各地市的污染物實際排放量分別除以各自的農(nóng)村面積可以得出各地市的地均污染物排放量，以此來反映廣東省農(nóng)業(yè)面源污染環(huán)境狀態(tài)。其中，農(nóng)村面積分別為21個地市的總面積扣除城市面積，根據(jù)廣東省自然資源廳公布的2008—2017年的土地利用現(xiàn)狀數(shù)據(jù)作為各地級市的總面積，參考《土地利用現(xiàn)狀分類》（GB/T 21010—2007），其中城市、建制鎮(zhèn)、機場用地、港口碼頭用地面積歸類為城市面積，從而估算出農(nóng)村面積。

式中，Pij地為第i個市縣的第j個污染物地均排放量，kg/km2；Pij為第i個市縣的第j個污染物排放量，t/年；Si為第i個市縣的農(nóng)村土地面積，km2。

求得COD、TN、TP地均排放量后進(jìn)行離差標(biāo)準(zhǔn)化，使其結(jié)果落于［0，1］，再對這3種污染強度標(biāo)準(zhǔn)化值按相同權(quán)重加和得到污染綜合指數(shù)，從而確定農(nóng)業(yè)面源污染風(fēng)險排序。

2 結(jié)果與分析

2.1 廣東省2008—2017年農(nóng)業(yè)面源污染時空特征分析

2.1.1 2008—2017年COD污染物排放時空特征農(nóng)業(yè)面源污染主要包括畜禽養(yǎng)殖污染物、水產(chǎn)養(yǎng)殖污染物和種植業(yè)污染物。畜禽養(yǎng)殖業(yè)、水產(chǎn)養(yǎng)殖業(yè)是COD污染物的主要污染源。近10年廣東省畜禽養(yǎng)殖量呈遞減趨勢，從2008年的2 647.39萬只，減少至2017年的2 346.80萬只。水產(chǎn)養(yǎng)殖量逐年增長，2008年為680.41萬t，2017年為833.54萬t，其增長率為22.51%，由此導(dǎo)致水產(chǎn)養(yǎng)殖排放污染物逐年增多，對COD污染物排放的貢獻(xiàn)率逐年擴(kuò)大，成為COD污染物的重要污染源。

由公式（1）可得廣東省2008—2017年COD污染物排放量，詳見表2。2008—2017年，廣東省COD污染物排放量總體呈遞減趨勢，至2017年為104 467.16萬kg，減少7 456.75萬kg。不同經(jīng)濟(jì)地區(qū)COD污染物排放差異明顯，具體表現(xiàn)：①珠三角地區(qū)近10年來COD污染物排放量減少9 125.80萬kg，深圳、珠海、中山、東莞4市一直保持COD污染物低排放趨勢（圖1a）。肇慶、江門2市畜禽養(yǎng)殖量、水產(chǎn)養(yǎng)殖量均處于珠三角9市前列，其中2017年肇慶市COD污染物排放量為10 941.12萬kg，占珠三角地區(qū)COD污染物總排放量的35.01%，是該經(jīng)濟(jì)地區(qū)COD污染物的主要貢獻(xiàn)市。②東翼地區(qū)COD污染物排放呈先緩慢遞增后驟降的特征，由2016年的12 143.98萬kg驟然下降至2017年的10 244.95萬kg，同比下降15.64%。一方面由于汕尾市水產(chǎn)養(yǎng)殖量從2016年的65.223 5萬t減少至2017年的55.773 6萬t，降幅14.49%；另一方面東翼4市2016—2017年畜禽養(yǎng)殖量均大量減少，2016年東翼4市的畜禽養(yǎng)殖總量為241.66萬只，2017年為196.67萬只，同比減少44.99萬只，由此導(dǎo)致東翼4市COD污染物排放在2017年明顯下降（圖1b）。可見，畜禽養(yǎng)殖是該經(jīng)濟(jì)區(qū)域COD污染物排放的主要來源，控制畜禽養(yǎng)殖量是改善農(nóng)業(yè)面源污染的重要措施。③西翼地區(qū)3市COD污染物排放量處于7 000萬～17 000萬kg（圖1c），2017年茂名COD污染物排放量是陽江的近2倍，茂名市的COD污染物排放量位列廣東省21個地市首位；山區(qū)5市COD污染物排放量處于3 000萬～9 000萬kg（圖1d）。綜上所述，2017年COD污染區(qū)域前10名依次為茂名、湛江、肇慶、陽江、清遠(yuǎn)、江門、梅州、韶關(guān)、佛山、惠州，可見COD污染嚴(yán)重地區(qū)主要分布在粵西、粵北地區(qū)。

圖1 廣東省2008—2017年不同經(jīng)濟(jì)區(qū)域COD污染物排放量

表2 廣東省2008—2017年不同經(jīng)濟(jì)區(qū)域COD污染物排放量

2.1.2 2008—2017年TN、TP污染物排放時空特征 畜禽養(yǎng)殖污染物、水產(chǎn)養(yǎng)殖污染物和化肥施用量是TN、TP污染物的主要來源。總體上，2008—2017年廣東省化肥施用量（折純）呈緩慢遞增趨勢，但2016年廣東省化肥施用量（折純）為261.02萬t，驟降至2017年的237.94萬t，降幅8.84%，主要與湛江、云浮、陽江、廣州等市減少施用化肥有關(guān)（表3）。2017年化肥施用量（折純）超過10萬t的地市有湛江、茂名、肇慶、清遠(yuǎn)、梅州、河源、江門、揭陽。其中，湛江市為46.32萬t，茂名為35.15萬t，遠(yuǎn)遠(yuǎn)超過其他地區(qū)。河源市2017年的化肥施用量為12.72萬t，其對TN污染物排放量的貢獻(xiàn)達(dá)59%以上，對TP污染物排放量的貢獻(xiàn)率為50.18%。過度施用化肥是造成河源市農(nóng)業(yè)面源污染的主要原因。

表3 廣東省2008—2017年化肥施用量（折純）

由公式（1）可得廣東省2008—2017年TN、TP污染物排放量，詳見圖2和圖3。珠三角地區(qū)TN、TP污染物排放趨勢特征相似，深圳、珠海、中山、東莞4市TN、TP污染物排放量較少，TN、TP污染物排放量均分別處于300萬kg、50萬kg以下（圖2a、圖3a）。肇慶TN、TP污染物排放量遠(yuǎn)遠(yuǎn)超過其他地區(qū)，江門TN、TP污染物排放量從2008年遞增至2016年，而后大幅下降。廣州、惠州、佛山3市的TN、TP污染物排放均呈下降趨勢，農(nóng)業(yè)面源污染趨勢有所減緩。東翼地區(qū)4市TN污染物排放量處于300萬～800萬kg，TP污染物排放量處于40萬～110萬kg，其中揭陽市整體呈下降趨勢（圖2b、圖3b）。西翼地區(qū)3市TN、TP污染物排放量均分別超過1 000萬kg、150萬kg，其中湛江、茂名TN、TP污染物排放量處于廣東省前列（圖2c、圖3c）。山區(qū)5市TN、TP污染物排放量趨勢特征相似，清遠(yuǎn)、河源、梅州、云浮的TN、TP污染物排放量均在2016—2017年出現(xiàn)劇烈變化（圖2d、圖3d）。其中，清遠(yuǎn)、河源的TN污染物排放量分別突增163.41萬、136.62萬kg，TP污染物排放量分別突增27.15萬、13.94萬kg，與清遠(yuǎn)畜禽養(yǎng)殖量增加、河源化肥施用量增加密切相關(guān)。綜上所述，2017年TN污染區(qū)域排在前10名的依次為茂名、湛江、肇慶、清遠(yuǎn)、陽江、江門、梅州、韶關(guān)、河源、揭陽；TP污染區(qū)域排在前10名的依次為為茂名、湛江、肇慶、陽江、清遠(yuǎn)、江門、梅州、韶關(guān)、河源、揭陽，TN、TP污染區(qū)域基本一致。

圖2 廣東省2008—2017年不同經(jīng)濟(jì)區(qū)域TN污染物排放量

圖3 廣東省2008—2017年不同經(jīng)濟(jì)區(qū)域TP污染物排放量

2.2 廣東省2008—2017年農(nóng)業(yè)面源污染風(fēng)險分析

2.2.1 2008—2017年農(nóng)業(yè)面源污染強度分析 采用ArcGIS 10.2空間分析技術(shù)，繪制2008年和2017年廣東省各污染源的地均排放量空間分布圖。由圖4可知，2008—2017年廣東省農(nóng)業(yè)面源污染地均排放量總體呈下降趨勢，空間差異明顯。具體分析可知，2008—2017年，種植規(guī)模、畜禽養(yǎng)殖規(guī)模大是造成地區(qū)COD、TN、TP排放量居高不下的主要原因，樂觀的是，廣東省地均排放量呈下降趨勢，農(nóng)業(yè)面源污染有所改善。COD地均排放量由中部地區(qū)向粵東、粵北遞減，與2008年相比，佛山、廣州、江門等市2017年COD地均排放量均下降一個污染等級。2017年茂名、佛山、珠海、汕頭4市COD地均排放量處于12 000～16 000 kg/km2，其他大部分地區(qū)低于8 000 kg/km2。TN、TP地均排放量總體波動穩(wěn)定并逐漸下降，湛江、廣州、佛山、揭陽均呈改善趨勢；韶關(guān)、河源、東莞3市TN、TP地均排放量均處于最低污染等級。

圖4 廣東省2008—2017年COD、TN、TP污染物地均排放量

2.2.2 2008—2017年農(nóng)業(yè)面源污染風(fēng)險分析 經(jīng)計算可得廣東省農(nóng)業(yè)面源污染綜合指數(shù)，2017年21個地市COD、TN、TP的地均排放量標(biāo)準(zhǔn)化值及綜合指數(shù)結(jié)果見表4。綜合指數(shù)從大到小順序為茂名、汕頭、湛江、珠海、佛山、陽江、中山、江門、揭陽、肇慶、潮州、汕尾、廣州、清遠(yuǎn)、云浮、深圳、梅州、惠州、韶關(guān)、河源、東莞。

表4 廣東省2017年農(nóng)業(yè)面源污染綜合指數(shù)及順序

分析可知，農(nóng)業(yè)面源污染風(fēng)險與單個污染強度指標(biāo)分布規(guī)律較為一致。其中，茂名和汕頭所得綜合分值為農(nóng)業(yè)面源污染風(fēng)險最高區(qū)域；湛江、珠海、佛山、陽江各污染物地均排放量的標(biāo)準(zhǔn)化值較大，處于高度污染風(fēng)險區(qū)域；中山、江門、揭陽、肇慶、潮州、汕尾為中度污染風(fēng)險區(qū)域；廣州、深圳、惠州等市為低度污染風(fēng)險區(qū)域，其中東莞市COD、TN、TP地均排放量的標(biāo)準(zhǔn)化值均為0.00，東莞市所得綜合評價為農(nóng)業(yè)面源污染風(fēng)險最低。總體而言，廣東省農(nóng)業(yè)面源污染呈中西部污染風(fēng)險高、東中北污染風(fēng)險低的規(guī)律（圖5）。低風(fēng)險區(qū)域主要為城市化率較高、農(nóng)業(yè)比重小的深圳市、廣州市等地區(qū)；高風(fēng)險區(qū)域表現(xiàn)為農(nóng)地面積少、規(guī)模大、集約化程度低的茂名市、湛江市等地區(qū)。

圖5 廣東省2017年農(nóng)業(yè)面源污染風(fēng)險分布

3 討論

農(nóng)村農(nóng)業(yè)面源污染受到農(nóng)業(yè)經(jīng)濟(jì)活動的影響，土地的數(shù)量、質(zhì)量和結(jié)構(gòu)決定了種植業(yè)產(chǎn)量。隨著廣東省高產(chǎn)、高效農(nóng)業(yè)的發(fā)展，化肥施用對糧食增產(chǎn)、農(nóng)民增收起到了積極作用［23］，但化肥施用量越來越大，所貢獻(xiàn)的TN、TP污染越來越明顯，對農(nóng)村增綠產(chǎn)生了消極影響。廣東省自開展了世行貸款廣東農(nóng)業(yè)面源污染治理項目與制定了化肥零增長的行動方案以來，化肥施用量呈下降趨勢，特別是2016—2017年大部分地市的化肥施用量下降明顯。推進(jìn)精準(zhǔn)施肥、調(diào)整化肥施用結(jié)構(gòu)、改進(jìn)施肥方式、有機肥替代化肥等是從技術(shù)角度減少化肥施用量、保證化肥零增長、提高農(nóng)業(yè)生產(chǎn)效率的重要措施。

養(yǎng)殖業(yè)是COD污染的主要來源，也是廣東省農(nóng)業(yè)面源污染的主要因素。近10年廣東省加強畜禽養(yǎng)殖污染防治監(jiān)管，特別是嚴(yán)格執(zhí)行禁養(yǎng)政策，畜禽養(yǎng)殖規(guī)模有所減少。截至2017年底，廣東省畜禽禁養(yǎng)區(qū)內(nèi)的2萬多個養(yǎng)殖場戶已全部搬遷或關(guān)閉，畜禽養(yǎng)殖污染防治效果顯著。然而，隨著畜禽養(yǎng)殖污染防治力度加大，廣東省水產(chǎn)養(yǎng)殖規(guī)模逐年擴(kuò)大，進(jìn)而造成的污染也越來越明顯。西翼地區(qū)3市水產(chǎn)養(yǎng)殖規(guī)模均位于廣東省前列，其中2017年湛江市水產(chǎn)養(yǎng)殖量達(dá)120多萬噸，位居廣東省第一，由此導(dǎo)致的農(nóng)業(yè)面源污染風(fēng)險顯著提升。2019年《廣東省打贏農(nóng)業(yè)農(nóng)村污染治理攻堅戰(zhàn)實施方案》提出加強水產(chǎn)養(yǎng)殖污染防治和水生生態(tài)保護(hù)。因此，亟需在污染區(qū)域大力推進(jìn)水產(chǎn)養(yǎng)殖污染防治，推動水產(chǎn)養(yǎng)殖尾水達(dá)標(biāo)排放，加快推廣水產(chǎn)生態(tài)養(yǎng)殖技術(shù)和模式，從而達(dá)到水生生態(tài)保護(hù)的目標(biāo)。

廣東省已出臺關(guān)于農(nóng)業(yè)農(nóng)村污染治理多項政策，但政策執(zhí)行與具體落實有待提高，農(nóng)民真正受益與獲得感有待提升。雖然世行貸款廣東農(nóng)業(yè)面源污染治理項目效果明顯，但參加世行項目平均每個農(nóng)戶每年獲補貼458元對于純農(nóng)戶而言無疑為杯水車薪。目前，世行貸款廣東省農(nóng)業(yè)面源污染治理項目仍處于定點實施，所獲經(jīng)驗還未全省推廣，僅部分農(nóng)戶受益。因此對于未參加該項目的農(nóng)戶而言，在沒有政策支持下仍購買中高毒農(nóng)藥，如敵百蟲、殺蟲雙等，不會主動回收農(nóng)藥、肥料包裝廢棄物，環(huán)保意識淺薄。因此，加大對農(nóng)業(yè)農(nóng)村污染治理政策支持、技術(shù)指導(dǎo)，加大對農(nóng)戶農(nóng)業(yè)生產(chǎn)的經(jīng)濟(jì)補貼，促使農(nóng)戶積極主動參與農(nóng)業(yè)農(nóng)村污染治理，以農(nóng)民獲得實際收益的方式，引導(dǎo)農(nóng)戶朝著親環(huán)境型的生產(chǎn)方式轉(zhuǎn)變［24，25］，以期控制與治理農(nóng)業(yè)面源污染、實現(xiàn)生態(tài)文明建設(shè)的目標(biāo)。

本研究中，農(nóng)業(yè)面源污染源重點研究了種植業(yè)與養(yǎng)殖業(yè)所排放的COD、TN、TP具體化學(xué)指標(biāo)并選取地均污染強度及其綜合指數(shù)作為污染風(fēng)險分析指標(biāo)。但因COD、TN、TP對土壤的污染閾值數(shù)據(jù)獲取等原因，并未對上述污染類型閾值風(fēng)險進(jìn)行研究，今后應(yīng)加強對此方面的研究，使研究結(jié)果更為全面。同時，農(nóng)業(yè)面源污染具有隱蔽性、難測性、地域性等特點［26］，從省級層面宏觀分析難免忽視縣級局部污染現(xiàn)狀，今后亟需從縣級層次進(jìn)一步完善廣東省農(nóng)業(yè)面源污染時空分析，使研究結(jié)果更能反映實際情況。

4 小結(jié)

農(nóng)業(yè)面源污染時空分布差異造成不同地域的污染風(fēng)險不同，亟需因地制宜提出防治措施。采用排污系數(shù)法估算廣東省2008—2017年農(nóng)業(yè)面源污染的COD、TN、TP污染負(fù)荷，并結(jié)合ArcGIS 10.2可視化探討其時空分布；通過各污染物的地均排放量及其綜合指數(shù)對廣東省21個地市農(nóng)業(yè)面源污染強度進(jìn)行風(fēng)險分析。結(jié)果表明：①2017年廣東省農(nóng)村地區(qū)COD、TN及TP污染物排放量分別為104.47萬、16.87萬、2.45萬t；2008—2017年COD、TN及TP地均排放量整體呈下降趨勢，但2016—2017年部分地區(qū)變化幅度明顯，主要是由于部分市化肥施用量、畜禽養(yǎng)殖量大量增加；②21個地市農(nóng)業(yè)面源污染風(fēng)險大小順序為茂名、汕頭、湛江、珠海、佛山、陽江、中山、江門、揭陽、肇慶、潮州、汕尾、廣州、清遠(yuǎn)、云浮、深圳、梅州、惠州、韶關(guān)、河源、東莞。鄉(xiāng)村振興背景下農(nóng)業(yè)面源污染的時空分析與防治對策的提出可為廣東省鄉(xiāng)村振興、農(nóng)村生態(tài)環(huán)境治理提供借鑒。