江西省農(nóng)田灰水足跡時(shí)空分析

傅春，陳毓迪，劉業(yè)忠，李帆，黃金燕

（南昌大學(xué)工程建設(shè)學(xué)院，南昌 330036）

隨著我國城鎮(zhèn)化發(fā)展進(jìn)程不斷加快，人均耕地正在逐年減少，為確?，F(xiàn)有耕地糧食生產(chǎn)能力不斷提升，大量施用農(nóng)藥和化肥成為了農(nóng)業(yè)生產(chǎn)過程中重要的保障措施。根據(jù)聯(lián)合國糧食及農(nóng)業(yè)組織（FAO）和《中國農(nóng)村統(tǒng)計(jì)年鑒》提供的相關(guān)數(shù)據(jù)，2017年全球化肥施用量為19 222.4萬t，同年我國農(nóng)業(yè)化肥施用量達(dá)到5 859.4萬t，約占全球化肥施用總量的30.5%，我國單位耕地面積化肥施用量也超過世界平均水平1倍多。由于農(nóng)田化肥的過量施用，我國面臨著巨大的環(huán)境污染風(fēng)險(xiǎn)，農(nóng)田面源污染的評估和治理迫在眉睫。

2008年，HOEKSTRA等首次提出灰水足跡概念，并將其定義為以水環(huán)境質(zhì)量標(biāo)準(zhǔn)為基準(zhǔn)，將水污染物稀釋至無害所需的淡水體積。相較于傳統(tǒng)的單因子指數(shù)評價(jià)法、模糊綜合評價(jià)法、水污染指數(shù)法等評價(jià)方法，灰水足跡概念的提出，有效地將水量和水污染聯(lián)系起來，通過水量直觀地表現(xiàn)出水污染嚴(yán)重程度，為水污染評價(jià)提供了一種新方法。在農(nóng)田面源污染方面，國內(nèi)外學(xué)者利用灰水足跡方法開展了相應(yīng)的研究：曹連海等計(jì)算內(nèi)蒙古河套灌區(qū)糧食生產(chǎn)灰水足跡后發(fā)現(xiàn)，節(jié)水灌溉技術(shù)的推廣可以有效降低糧食生產(chǎn)灰水足跡；班榮舶等對重慶市種植業(yè)灰水足跡核算后得出重慶市水環(huán)境保護(hù)要重視種植業(yè)產(chǎn)生的污染；付永虎等借助灰水足跡概念對洞庭湖糧食生產(chǎn)環(huán)境效應(yīng)進(jìn)行評價(jià)，得出區(qū)域水環(huán)境受糧食生產(chǎn)負(fù)面影響會持續(xù)增強(qiáng)的結(jié)論；何立新等通過研究河北省主要農(nóng)作物相應(yīng)的灰水足跡，提出了減少污染的相關(guān)措施；張鑫等針對汾河流域農(nóng)業(yè)灰水足跡進(jìn)行了時(shí)空和效率分析，并就汾河流域種植業(yè)污染依舊呈緩慢增長的變化趨勢提出了相應(yīng)的解決措施。

但目前灰水足跡研究中存在著淋溶系數(shù)均一化問題，即當(dāng)研究區(qū)域面積較大時(shí)，大多數(shù)學(xué)者通常采用統(tǒng)一的污染物淋溶系數(shù)來代指整個(gè)研究區(qū)域。污染物淋溶系數(shù)通常受到土地利用方式、地形坡度等因素的共同影響，采用統(tǒng)一的淋溶系數(shù)容易忽視不同地形地貌對污染物流失系數(shù)的影響，可能導(dǎo)致在污染物空間分布研究中存在一定誤差。

目前雖然有部分學(xué)者就江西省農(nóng)田污染進(jìn)行過相應(yīng)的研究，但大多數(shù)是以市級行政區(qū)為研究對象，而以江西省縣（區(qū)）級為研究區(qū)域，利用地理信息系統(tǒng)軟件ArcGIS對農(nóng)田面源污染進(jìn)行時(shí)空分析的研究鮮有報(bào)道。本研究運(yùn)用ArcGIS軟件相關(guān)功能，結(jié)合《全國農(nóng)田面源污染排放系數(shù)手冊》等相關(guān)資料，對鄱陽湖流域污染系數(shù)進(jìn)行空間化處理，不同地區(qū)采用不同污染系數(shù)，實(shí)現(xiàn)了將污染系數(shù)和坡度與土地利用方式結(jié)合，使不同地區(qū)種植業(yè)灰水足跡測算更加準(zhǔn)確科學(xué)，以期為江西省農(nóng)業(yè)綠色化發(fā)展和水環(huán)境保護(hù)提供一定的參考依據(jù)。

1 材料與方法

1.1 研究區(qū)域概況

江西省總面積為16.69萬km，區(qū)域內(nèi)地形地貌復(fù)雜，流域內(nèi)南部、西部、東部均被山地環(huán)繞，由外向里地貌逐漸變?yōu)榍鹆甑匦?，流域北部地貌則主要為平原地形。境內(nèi)水系發(fā)達(dá)、水資源豐富，由“五江一湖”構(gòu)成的水網(wǎng)連通全省各處，流域內(nèi)氣候主要為亞熱帶季風(fēng)氣候，年平均降水量為1 341～1 943 mm，降水多集中于每年的4—7月，東南豐富，西北較少，水資源在空間分布上呈現(xiàn)北多南少的格局。江西省作為我國南方地區(qū)重要的糧食生產(chǎn)基地之一，其糧食作物種植面積總計(jì)達(dá)5.6×10hm（2020年），年平均化肥施用量達(dá)192.75 kg·hm，是歐美國家的數(shù)倍以上，而我國平均氮肥實(shí)際利用率卻只有30%左右，磷肥利用率不足20%，大量未被植物吸收的氮、磷流入周邊水體，引發(fā)水體富營養(yǎng)化，造成嚴(yán)重水污染問題。

1.2 研究方法

1.2.1 污染系數(shù)空間化

當(dāng)研究區(qū)域范圍較大時(shí)，由于化肥淋溶率大小受種植方式、地形地貌、施肥量等因素的影響，不同條件下農(nóng)田內(nèi)的氮、磷流失量存在一定差異。江西省內(nèi)地形地貌復(fù)雜，多地區(qū)開展實(shí)驗(yàn)測量費(fèi)時(shí)費(fèi)力，而全區(qū)域采用統(tǒng)一的流失系數(shù)又容易忽視地區(qū)差異性。故文章運(yùn)用ArcGIS將鄱陽湖流域內(nèi)農(nóng)田區(qū)域按“農(nóng)區(qū)-坡度-種植”的分類原則進(jìn)行相應(yīng)劃分（圖1），并結(jié)合《全國農(nóng)田面源污染排放系數(shù)手冊》對不同種植模式下的農(nóng)區(qū)進(jìn)行污染系數(shù)賦值，以實(shí)現(xiàn)研究區(qū)域污染系數(shù)空間化。圖1中各數(shù)值代表的含義如表1所示。

表1 種植模式類型編號與名稱Table 1 Numbers and names of planting pattern types

圖1 鄱陽湖流域種植模式空間分布Figure 1 Spatial distribution of planting patterns in Poyang Lake watershed

根據(jù)《全國農(nóng)田面源污染排放系數(shù)手冊》中介紹的相關(guān)使用方法，依據(jù)中科院資源環(huán)境科學(xué)與數(shù)據(jù)中心數(shù)據(jù)平臺提供的2000—2020年江西土地利用相關(guān)數(shù)據(jù)，參考孫鋮等構(gòu)建的方法，用公式（1）計(jì)算出江西省多年平均氮、磷流失系數(shù)，用于灰水足跡時(shí)間層次分析：

式中：β為某一田塊下某農(nóng)作物種植面積占市級行政區(qū)農(nóng)作物總播種面積的比例；μ為該田塊的氮、磷流失系數(shù)。

同時(shí)選取2020年相關(guān)數(shù)據(jù)，依照手冊對不同種植模式進(jìn)行污染排放系數(shù)賦值。用相同的方法計(jì)算得到江西省各地區(qū)的氮、磷流失系數(shù)，用于灰水足跡空間分布分析。

1.2.2 灰水足跡

（1）種植業(yè)灰水足跡

種植業(yè)灰水足跡主要來源于農(nóng)業(yè)生產(chǎn)過程中農(nóng)藥與化肥的過度使用，未被植物吸收利用的化肥和農(nóng)藥在降水和地表徑流的作用下，通過淋溶的方式進(jìn)入地表水體。由于鉀離子能被土壤膠體離子吸附，因此鉀肥不容易被淋濾。因?yàn)槭?nèi)農(nóng)藥施用種類繁多，各類農(nóng)藥施用量統(tǒng)計(jì)數(shù)據(jù)難以獲取，各農(nóng)藥允許濃度不同，并且稀釋氮磷污染的水體和稀釋農(nóng)藥污染的水體會有重疊部分，故選用氮、磷為種植業(yè)污染物進(jìn)行計(jì)算。參考HOEKSTRA等的核算方法，核算公式為：

式中：GWF為種植業(yè)灰水足跡，m·a；L為第類污染物的排放負(fù)荷，kg；α為第種化學(xué)物質(zhì)進(jìn)入水體的比例，即淋溶率、徑流損失率，%；Appl為第種化學(xué)物質(zhì)在生產(chǎn)過程中的使用量，kg；為第種非點(diǎn)源污染環(huán)境允許最大濃度，kg·m；為第種非點(diǎn)源污染在環(huán)境中的本底濃度，kg·m。

（2）灰水足跡強(qiáng)度

農(nóng)業(yè)灰水足跡強(qiáng)度反映單位耕地面積承載的灰水足跡大小，體現(xiàn)了區(qū)域農(nóng)業(yè)水污染程度，其公式為：

式中：為農(nóng)業(yè)灰水足跡強(qiáng)度，m·hm；為耕地面積，hm。

1.2.3 相關(guān)數(shù)據(jù)及來源

本研究的基礎(chǔ)數(shù)據(jù)資料主要包括江西省100個(gè)縣市氮磷化肥折純量、各區(qū)域相關(guān)農(nóng)作物播種面積、氮磷流失系數(shù)及計(jì)算流失系數(shù)所需的相關(guān)數(shù)據(jù)。其中用于地形坡度劃分的DEM高程數(shù)據(jù)來源于地理空間數(shù)據(jù)云，30 m精度土地利用數(shù)據(jù)來源于中科院資源環(huán)境科學(xué)與數(shù)據(jù)中心數(shù)據(jù)平臺（http：//www.resdc.cn），水質(zhì)標(biāo)準(zhǔn)依據(jù)《地表水環(huán)境質(zhì)量標(biāo)準(zhǔn)》（GB 3838—2002），將Ⅲ類水的污染物濃度設(shè)為環(huán)境最大容許濃度，總磷（TP）為0.2 mg·L，總氮（TN）為1.0 mg·L，自然本底濃度以0計(jì)。各區(qū)域中相關(guān)農(nóng)作物種植面積和農(nóng)田化肥折純量來源于《江西統(tǒng)計(jì)年鑒》和各市級行政區(qū)統(tǒng)計(jì)年鑒，農(nóng)田中氮、磷肥總施用量由單一化肥和復(fù)合肥相加得到，復(fù)合肥中氮、磷、鉀比例參照江西省農(nóng)業(yè)農(nóng)村廳相關(guān)報(bào)告，取1∶1∶1。

2 結(jié)果與分析

2.1 農(nóng)田氮磷灰水足跡時(shí)空分析

2.1.1 時(shí)序分析

江西省化肥施用總量和氮磷灰水足跡變化趨勢如圖2所示。在2000—2020年期間，江西省農(nóng)田化肥施用量呈先增長后下降的變化趨勢，其中2000—2006年期間增長快速，年平均增長率為3.2%，2007—2015年期間波動變化，總體呈上升趨勢，2015年達(dá)到最高點(diǎn)，為142.71萬t，約占當(dāng)年全國化肥施用總量的2.4%；隨后，2016—2020年，江西省化肥施用總量快速下降，在2020年已下降到108.81萬t，下降率為23.75%，體現(xiàn)出近幾年來江西省綠色生態(tài)農(nóng)業(yè)建設(shè)取得較大進(jìn)展。對于各類化肥來說，在過去的21 a間，氮肥、磷肥和鉀肥占總化肥施用量的比例逐年減少，復(fù)合肥占比逐年增加，已成為了農(nóng)業(yè)生產(chǎn)化肥中最主要的部分，需要額外加以關(guān)注。

圖2 2000—2020年鄱陽湖流域化肥施用量及灰水足跡變化Figure 2 Changes in fertilizer application and grey water footprint in the Poyang Lake basin from 2000 to 2020

就氮肥灰水足跡來說，在2000—2020年期間，江西省農(nóng)田氮肥灰水足跡總體呈現(xiàn)先增后降的變化趨勢，其峰值出現(xiàn)在2011年，為7.27×10m，較2000年增加12.1%，隨后一直保持下降趨勢，直到2020年氮肥灰水足跡已降至5.35×10m。

就磷肥灰水足跡來說，峰值出現(xiàn)在2015年，為1.47×10m，相較2000年增幅較大，增加了33.73%。此后，在2015—2020年期間，磷肥灰水足跡逐漸下降，在2020年降至1.12×10m，下降幅度為23.56%。

2.1.2 空間分析

鄱陽湖流域是我國重要糧食生產(chǎn)基地，為保障農(nóng)業(yè)生產(chǎn)每年會使用大量化肥。區(qū)域內(nèi)各個(gè)市的行政面積和地形地貌存在明顯差異，為更加準(zhǔn)確分析鄱陽湖流域內(nèi)不同地區(qū)的種植業(yè)灰水足跡，本研究利用ArcGIS對區(qū)域內(nèi)100個(gè)縣區(qū)進(jìn)行空間分布特征分析，以更加準(zhǔn)確直觀地顯示區(qū)域氮磷灰水足跡狀況，結(jié)果見圖3。

圖3 2020年江西省農(nóng)田氮、磷灰水足跡空間分布Figure 3 Spatial distribution of nitrogen and phosphorus grey water footprint of farmland in Jiangxi Province in 2020

對于氮肥灰水足跡來說，流域內(nèi)總體呈現(xiàn)由外向內(nèi)逐漸增多的變化趨勢，高氮肥灰水足跡區(qū)域主要集中地區(qū)為鄱陽湖平原及周邊地區(qū)，例如上饒市的鄱陽縣和余干縣、宜春市的豐城市和樟樹市，以及南昌市的南昌縣和進(jìn)賢縣，其中南昌縣、豐城市和余干縣氮肥灰水足跡最高，分別為2.26×10、1.99×10m和1.84×10m；除以上縣區(qū)外，寧都縣、泰和縣、南豐縣、吉水縣、渝水區(qū)、臨川區(qū)、高安市、永修縣和新建區(qū)境內(nèi)氮肥灰水足跡也相對較高，值得加以關(guān)注。

磷肥灰水足跡空間分布與氮肥灰水足跡基本相同，江西中東部地區(qū)為高磷肥灰水足跡區(qū)域，南昌縣、余干縣、豐城縣、鄱陽縣為全省磷肥灰水足跡最高區(qū)域，分別為5.56×10、3.98×10、3.92×10m和3.82×10m。萍鄉(xiāng)市、景德鎮(zhèn)市、宜春市西部、贛州市南部、上饒市東部和吉安市磷肥灰水足跡相對其他地區(qū)較低，而撫州市北部、宜春市東部、南昌市和九江市南部等地區(qū)存在磷肥灰水足跡較高現(xiàn)象，需要加強(qiáng)綠色生態(tài)農(nóng)業(yè)建設(shè)，減少化肥使用。

2.2 農(nóng)田氮磷灰水足跡強(qiáng)度空間分析

灰水足跡僅反映了區(qū)域內(nèi)種植業(yè)污染量的多少，由于不同地區(qū)種植面積差異較大，單純從灰水足跡角度無法分辨地區(qū)農(nóng)田化肥污染程度，故采用灰水足跡強(qiáng)度來反映不同地區(qū)單位耕地在農(nóng)業(yè)生產(chǎn)過程中承載的面源污染程度，計(jì)算后結(jié)果等分為5個(gè)層次，分別為高污染強(qiáng)度地區(qū)、較高污染強(qiáng)度地區(qū)、中污染強(qiáng)度地區(qū)、較低污染強(qiáng)度地區(qū)、輕污染強(qiáng)度地區(qū)，結(jié)果如圖4所示。

圖4 2020年江西省農(nóng)田氮、磷灰水足跡強(qiáng)度空間分布Figure 4 Spatial distribution of nitrogen and phosphorus grey water footprint intensity of farmland in Jiangxi Province in 2020

從氮肥灰水足跡強(qiáng)度來看，鷹潭市和吉安市整體農(nóng)田氮肥污染程度較輕；中污染強(qiáng)度地區(qū)主要位于宜春市東部、九江市北部、撫州市南部和贛州市中部與西部地區(qū)。其余地區(qū)，如上饒市弋陽縣、萍鄉(xiāng)市上栗縣和新余市渝水縣均為中污染強(qiáng)度地區(qū)；鉛山縣、永修縣、樟樹市、奉新縣、珠江區(qū)和紅谷灘區(qū)為較高污染強(qiáng)度地區(qū)；全流域氮肥灰水足跡強(qiáng)度最大區(qū)域是撫州市南豐縣，為2 873.38 m·hm，為高污染強(qiáng)度地區(qū)。

從磷肥灰水足跡強(qiáng)度來看，吉安市、宜春市西部、萍鄉(xiāng)市北部、九江市西部和上饒市大部分地區(qū)為低和較低污染強(qiáng)度地區(qū)，其中吉安市低污染強(qiáng)度地區(qū)占全市面積的54%，說明當(dāng)?shù)剞r(nóng)業(yè)綠色生態(tài)化發(fā)展良好。中污染強(qiáng)度地區(qū)為上饒市鉛山縣、撫州市金溪縣、萍鄉(xiāng)市蓮花縣、宜春市奉新縣和樟樹市、贛州市會昌縣和崇義縣、南昌市安義縣和南昌縣；景德鎮(zhèn)市珠江區(qū)、南昌市紅谷灘區(qū)和九江市永修縣為較高污染強(qiáng)度地區(qū)，其磷肥灰水足跡強(qiáng)度依次為5 578.60、5 056.87 m·hm和4 651.24 m·hm；撫州市南豐縣為全區(qū)域磷肥灰水足跡強(qiáng)度最高地區(qū)，為6 885.59 m·hm。

3 討論

相較于其他研究方法，運(yùn)用灰水足跡理論來評估地區(qū)農(nóng)業(yè)化肥面源污染可有效地展示化肥污染對地區(qū)水資源的影響，是目前農(nóng)田面源污染研究的主要方法之一。

3.1 區(qū)域主要污染物及來源分析

本研究發(fā)現(xiàn)，在過去21 a江西省農(nóng)田化肥施用過程中，復(fù)合肥占總化肥施用量比例逐年增加，是目前農(nóng)田面源污染的主要來源。比較氮肥和磷肥灰水足跡可以發(fā)現(xiàn)，盡管氮肥施用量大于磷肥施用量，但由于自然水體對磷元素的最大容許濃度遠(yuǎn)小于氮元素的容許濃度，使得磷肥灰水足跡大于氮肥灰水足跡，根據(jù)灰水足跡中“木桶理論”，相較于施用量更大的氮肥，灰水足跡更高的磷肥才是農(nóng)田生產(chǎn)的主要污染來源，這與《江西省環(huán)境狀況公報(bào)》中的結(jié)論一致。因而在判斷主要污染物時(shí)，不僅需要比較污染物排放量的多少，同時(shí)也需要考慮水體對不同污染物耐受程度的不同。

3.2 區(qū)域重污染程度地區(qū)分析

農(nóng)田氮磷污染較嚴(yán)重區(qū)域大多集中在江西省中部平原地區(qū)，其中南昌市灰水足跡強(qiáng)度達(dá)到了3 686.23 m·hm（2010—2020年平均值），為江西省地級市中農(nóng)田污染最為嚴(yán)重的區(qū)域?？h區(qū)級方面，撫州市南豐縣為農(nóng)田化肥污染最為嚴(yán)重區(qū)域，比較南豐縣肥料施用情況發(fā)現(xiàn)，其復(fù)合肥施用量占全縣化肥施用量的44%，占撫州市復(fù)合肥施用總量的20%，遠(yuǎn)超過撫州市其他地區(qū)。景德鎮(zhèn)市珠山區(qū)、九江市永修縣和南昌市紅谷灘區(qū)為較高風(fēng)險(xiǎn)地區(qū)，其中景德鎮(zhèn)市珠山區(qū)用于蔬菜種植的耕地面積占其整體耕地面積的98%，在單一化種植模式下，當(dāng)?shù)剞r(nóng)民為保障作物產(chǎn)量而過量施用化肥，導(dǎo)致單位耕地面積下化肥使用量過高。九江市永修縣復(fù)合肥和磷肥施用量是周邊縣市的2～3倍，南昌市紅谷灘區(qū)與東湖區(qū)種植面積相近，但其各類化肥施用量是東湖區(qū)的數(shù)倍以上，兩地均存在化肥施用強(qiáng)度過高的情況，值得地區(qū)相關(guān)部門重視。

3.3 種植業(yè)面源污染防治對策分析

借鑒我國農(nóng)業(yè)面源污染防治的相關(guān)政策，結(jié)合研究區(qū)內(nèi)實(shí)際情況，本文從技術(shù)和政府管理兩方面提出相關(guān)的意見建議。

（1）技術(shù)方面

研究結(jié)果表明，占比不斷增加的復(fù)合肥是種植業(yè)面源污染的主要污染物來源。因此，可以在確保農(nóng)業(yè)生產(chǎn)不受影響的前提下，選用磷含量占比較低的復(fù)合肥品種，同時(shí)降低化肥的整體施用量，并嘗試搭配有機(jī)肥來改善耕地肥力，提高種植業(yè)生產(chǎn)水平，進(jìn)而從源頭上減少農(nóng)業(yè)化肥對水環(huán)境的影響。在坡度較高的地區(qū)，可以通過修建梯田來降低地形坡度，減少化肥水土流失風(fēng)險(xiǎn)；在平原或丘陵地區(qū)，可以在農(nóng)田水渠和周邊水體周圍種植具有過濾、吸附能力的植被，從傳輸路徑上減少化肥污染物的入河量。

（2）政府管理方面

相關(guān)環(huán)保部門要從宣傳和監(jiān)督兩方面開展工作，為農(nóng)戶開展科學(xué)種植知識的普及與教育，減少農(nóng)戶施肥的盲目性，提高農(nóng)民對面源污染的認(rèn)識程度。同時(shí)監(jiān)督各區(qū)域化肥施用情況，對農(nóng)田高污染的地區(qū)開展調(diào)研，找出污染的主要問題，推廣綠色生態(tài)農(nóng)業(yè)建設(shè)工作，緩解農(nóng)業(yè)生產(chǎn)和環(huán)境問題之間的矛盾，以達(dá)到農(nóng)業(yè)面源污染控制和管理目標(biāo)。

目前，農(nóng)田種植業(yè)面源污染依舊是地區(qū)面源污染的主要組成部分，是導(dǎo)致地表水富營養(yǎng)化的主要原因，對農(nóng)田周邊水環(huán)境構(gòu)成嚴(yán)重威脅。農(nóng)田面源污染主要是由當(dāng)?shù)剞r(nóng)業(yè)生產(chǎn)過程中化肥的過量施用和施肥不當(dāng)導(dǎo)致的，杜鵑等研究發(fā)現(xiàn)耕地總氮排放貢獻(xiàn)率大于人口和畜禽養(yǎng)殖，項(xiàng)頌等對洱海流域農(nóng)業(yè)面源污染聚類分析發(fā)現(xiàn)，總氮和總磷的排放量均主要來自于種植業(yè)。因此，盡管近幾年來農(nóng)業(yè)化肥施用量不斷降低，但依舊需要對農(nóng)田化肥污染保持較高的警惕。

4 結(jié)論

（1）鄱陽湖流域農(nóng)田灰水足跡呈現(xiàn)先增后降趨勢，且近幾年降低趨勢明顯，說明整體生態(tài)農(nóng)業(yè)建設(shè)取得較大進(jìn)展。

（2）境內(nèi)氮磷灰水足跡及強(qiáng)度分布特征相近，呈現(xiàn)中部地區(qū)污染嚴(yán)重、四周區(qū)域污染較輕的分布規(guī)律，其中污染較為嚴(yán)重的區(qū)域?yàn)槟喜?、宜春市東部地區(qū)、南豐縣、鉛山縣和永修縣。

（3）磷是鄱陽湖流域農(nóng)田生產(chǎn)主要污染物，占比逐年增加的復(fù)合肥是農(nóng)田污染的主要來源。建議選用配比更加合理的復(fù)合肥，采取科學(xué)施肥方式，提高化肥利用率，實(shí)現(xiàn)從源頭減少種植業(yè)污染。