土地經(jīng)營規(guī)模與農(nóng)業(yè)面源污染的關(guān)系研究

摘 要：面源污染這一環(huán)境問題日益嚴(yán)重，現(xiàn)已成為國內(nèi)外專家學(xué)者研究的活躍領(lǐng)域。本文以農(nóng)業(yè)生產(chǎn)中污染性要素（化肥、農(nóng)藥和農(nóng)用塑料薄膜）的投入為例，利用2000—2012年省際面板數(shù)據(jù)，運(yùn)用簡約式與結(jié)構(gòu)式模型，實(shí)證分析土地經(jīng)營規(guī)模對(duì)農(nóng)業(yè)面源污染的影響。研究發(fā)現(xiàn)，我國化肥、農(nóng)藥和農(nóng)膜的投入均與人均收入之間存在倒U型曲線關(guān)系，倒U型曲線轉(zhuǎn)折點(diǎn)對(duì)應(yīng)的人均收入水平大于全國大部分地區(qū)及全國平均水平，且規(guī)模經(jīng)營的控制變量使各曲線轉(zhuǎn)折點(diǎn)對(duì)應(yīng)的人均收入，意味著面源污染負(fù)荷強(qiáng)度會(huì)隨著人均收入的增加而增加化肥的投入與人均耕地面積、非農(nóng)就業(yè)水平、農(nóng)業(yè)機(jī)械化水平顯著正相關(guān)，與勞均農(nóng)作物播種面積顯著負(fù)相關(guān)；農(nóng)藥的投入與人均耕地面積和農(nóng)業(yè)機(jī)械化水平顯著正相關(guān)，與勞均農(nóng)作物播種面積顯著負(fù)相關(guān)。

關(guān)鍵詞：規(guī)模經(jīng)營 面源污染

一、引言

國外研究資料表明，地球約有30-50%的受到面源污染。我國面源污染由于其分散廣、很難監(jiān)測和量化等特點(diǎn)，長期受到忽視，面源污染帶來的水域富營養(yǎng)化是污染重要表現(xiàn)，其程度和廣度甚至超過許多發(fā)達(dá)國家，其潛在的威脅更是其他國家無法比擬的?，F(xiàn)在，我國農(nóng)業(yè)已成為最大的面源污染產(chǎn)業(yè)，而非工業(yè)，農(nóng)業(yè)面源污染形勢嚴(yán)峻。但是，農(nóng)業(yè)在我國產(chǎn)業(yè)結(jié)構(gòu)中扮演著重要的角色，從發(fā)達(dá)國家農(nóng)業(yè)改革的實(shí)踐經(jīng)驗(yàn)可以看出，農(nóng)業(yè)生產(chǎn)效率的提高和農(nóng)業(yè)現(xiàn)代化的一個(gè)重要方式就是實(shí)現(xiàn)農(nóng)地的適度規(guī)模經(jīng)營。2014 年中央一號(hào)文件指出構(gòu)建新型農(nóng)業(yè)經(jīng)營體系要發(fā)展多種形式規(guī)模經(jīng)營，鼓勵(lì)有條件的農(nóng)戶流轉(zhuǎn)承包土地的經(jīng)營權(quán)，推動(dòng)土地經(jīng)營權(quán)流轉(zhuǎn)市場的不斷健全?？梢娢覈鴮?duì)農(nóng)地適度規(guī)模經(jīng)營問題非常重視，因地制宜發(fā)展土地經(jīng)營規(guī)模是我國農(nóng)業(yè)發(fā)展的方向。

相對(duì)于現(xiàn)實(shí)的緊迫性，有關(guān)土地經(jīng)營規(guī)模與農(nóng)業(yè)面源污染兩者關(guān)系的學(xué)術(shù)研究卻嚴(yán)重滯后，使得更為準(zhǔn)確地認(rèn)識(shí)兩者的關(guān)系缺乏經(jīng)驗(yàn)支撐和理論基礎(chǔ)。因此，利用理論分析和實(shí)證模型分析兩者的關(guān)系刻不容緩。本文以 EKC為研究基礎(chǔ)，通過建立EKC模型，利用我國 2000年-2012年31個(gè)?。ㄗ灾螀^(qū)、直轄市）的面板數(shù)據(jù)，以農(nóng)業(yè)生產(chǎn)中化肥、農(nóng)藥、農(nóng)用塑料薄膜投入為例，分析土地經(jīng)營規(guī)模各要素與污染性要素投入之間的關(guān)系。

二、國內(nèi)外研究回顧

目前，環(huán)境庫茲涅茨曲線（EKC）主要用于分析探討經(jīng)濟(jì)發(fā)展中的環(huán)境污染的演變趨勢，研究對(duì)象主要是工業(yè)污染或者城市污染。隨著農(nóng)業(yè)污染問題的日益嚴(yán)重，國外學(xué)者開始將EKC曲線運(yùn)用到農(nóng)業(yè)面源污染與經(jīng)濟(jì)增長關(guān)系的問題，國內(nèi)學(xué)者主要以農(nóng)資如化肥、農(nóng)藥和農(nóng)膜等為變量，或者以化學(xué)產(chǎn)物如氮、氨和農(nóng)業(yè)化學(xué)需氧量為變量研究其與經(jīng)濟(jì)增長之間的關(guān)系，驗(yàn)證是否符合EKC簡約式或結(jié)構(gòu)式模型等。

綜上所述，利用EKC模型分析農(nóng)業(yè)問題的研究仍不夠豐富，尤其是針對(duì)土地經(jīng)營規(guī)模所帶來的規(guī)模、結(jié)構(gòu)和技術(shù)效應(yīng)與農(nóng)業(yè)面源污染的關(guān)系的研究仍然處于空白。鑒于此，本文將簡約式與結(jié)構(gòu)式分析相結(jié)合，對(duì)土地經(jīng)營規(guī)模與以污染性要素投入為表征的農(nóng)業(yè)面源環(huán)境問題進(jìn)行實(shí)證研究，并分析規(guī)模經(jīng)營等控制變量的影響。

三、研究數(shù)據(jù)與研究方法

本文數(shù)據(jù)來源于《新中國60年統(tǒng)計(jì)資料匯編：1949-2008》，《中國統(tǒng)計(jì)年鑒》，《中國農(nóng)村統(tǒng)計(jì)年鑒》以及各省的統(tǒng)計(jì)年鑒。人均收入（I）。以各地區(qū)2000年農(nóng)林牧漁從業(yè)人員數(shù)為基期的人均農(nóng)林牧漁總產(chǎn)值來表示。

1.土地經(jīng)營規(guī)模指標(biāo)。本文采用人均耕地面積、勞均農(nóng)作物播種面積、單位面積農(nóng)業(yè)機(jī)械動(dòng)力、有效灌溉比例來衡量土地經(jīng)營規(guī)模狀況。人均耕地面積（Pl）。以各省農(nóng)村居民家庭經(jīng)營耕地面積表示。人均耕地面積的大小將直接影響耕地市場的價(jià)格和流轉(zhuǎn)租金，也是耕地資源的稀缺程度的表現(xiàn)。而較高的農(nóng)地流轉(zhuǎn)價(jià)格可能使得供需難以達(dá)到平衡，阻礙土地經(jīng)營規(guī)模的發(fā)展；當(dāng)然，人均耕地面積越小，也可能意味著耕地細(xì)碎且分散，這樣的結(jié)果就是耕地流轉(zhuǎn)所附帶的交易成本增加，耕地的細(xì)碎化不利于規(guī)模經(jīng)營主體進(jìn)行現(xiàn)代化農(nóng)業(yè)生產(chǎn)，這樣一來，其經(jīng)營規(guī)模的意愿就隨之降低，從而影響土地經(jīng)營規(guī)模狀況。

勞均農(nóng)作物播種面積（Lp）。以各省鄉(xiāng)村農(nóng)林牧漁從業(yè)人員為基期的人均農(nóng)作物播種總面積。即：勞均農(nóng)作物播種面積=農(nóng)作物播種面積/農(nóng)林牧漁從人員數(shù)。勞均農(nóng)作物播種面積越大，則表明勞動(dòng)力的生產(chǎn)率越高，意味著規(guī)模經(jīng)營水平越高。

非農(nóng)就業(yè)水平（Na）。鄉(xiāng)村勞動(dòng)力中非農(nóng)就業(yè)人數(shù)所占的比例。非農(nóng)就業(yè)人數(shù)是除農(nóng)林牧漁從業(yè)人員外從事其他產(chǎn)業(yè)的勞動(dòng)力。非農(nóng)就業(yè)水平=非農(nóng)就業(yè)人數(shù)/勞動(dòng)力人數(shù)。從事第一產(chǎn)業(yè)的農(nóng)村勞動(dòng)力向第二三產(chǎn)業(yè)轉(zhuǎn)移，使大量土地?zé)o勞動(dòng)力耕種，閑置的土地有助于土地集中和發(fā)展規(guī)模經(jīng)營。

農(nóng)業(yè)機(jī)械化水平（Am）。本文采用單位農(nóng)作物播種面積的農(nóng)業(yè)機(jī)械投入量來表示。即農(nóng)業(yè)機(jī)械化水平=農(nóng)業(yè)機(jī)械總動(dòng)力/農(nóng)作物播種面積。農(nóng)業(yè)機(jī)械的投入是保證規(guī)模經(jīng)營的基礎(chǔ)，農(nóng)業(yè)機(jī)械化水平的變化，可能帶來土地經(jīng)營規(guī)模狀況的改變。

2.農(nóng)業(yè)面源污染指標(biāo)。農(nóng)業(yè)面源污染程度主要采用化肥污染、農(nóng)藥污染、農(nóng)膜污染污染指標(biāo)。對(duì)于化肥、農(nóng)藥和農(nóng)膜等農(nóng)資污染物排放量指標(biāo)，用投入密度來表示，即用單位播種面積投入量作為農(nóng)業(yè)面源污染變量，F(xiàn)er、Pes、Pla分別表示單位農(nóng)作物播種面積化肥、農(nóng)藥和農(nóng)膜使用量。

本文采用農(nóng)業(yè)面源污染產(chǎn)生總量來衡量面源污染程度。農(nóng)村面源污染產(chǎn)生總量核算公式如下：

式（2.1）和式（2.2）為農(nóng)業(yè)面源污染產(chǎn)生總量核算公式。其中，表示i地區(qū)第j年農(nóng)業(yè)是面源污染第K種污染物排放總量；表示i地區(qū)第j年單位農(nóng)作物播種面積上污染負(fù)荷量（簡稱污染負(fù)荷強(qiáng)度）；表示i地區(qū)第j年農(nóng)作物播種面積。

3.模型建立。本文在EKC模型的基礎(chǔ)上，選擇農(nóng)業(yè)面源污染排放強(qiáng)度作為因變量，實(shí)際人均收入和土地經(jīng)營規(guī)模水平的相關(guān)指標(biāo)作為自變量來擬合方程，構(gòu)建土地經(jīng)營規(guī)模對(duì)農(nóng)業(yè)面源污染影響的理論分析框架，即：

其中，AM為污染負(fù)荷強(qiáng)度，I為實(shí)際人均收入，Pl為人均耕地面積，Lp為勞均農(nóng)作物播種面積，Na為非農(nóng)就業(yè)水平，Am為農(nóng)業(yè)機(jī)械化水平。本文將利用2000年至2012年我國各省的面板數(shù)據(jù)先估計(jì)（2.3）式并計(jì)算轉(zhuǎn)折點(diǎn)所對(duì)應(yīng)的收入水平，再對(duì)（2.4）式進(jìn)行估計(jì)并與（2.3）式進(jìn)行對(duì)比，分析EKC的穩(wěn)健性和控制變量的影響。

上述各變量統(tǒng)計(jì)性描述見表1。

四、研究結(jié)果與分析

1.對(duì)（2.3）式的估計(jì)。本文利用stata12.0對(duì)面板數(shù)據(jù)進(jìn)行回歸分析，采用固定效益模型和隨機(jī)效益模型進(jìn)行估計(jì)，并通過Hausman檢驗(yàn)確定為計(jì)算EKC轉(zhuǎn)折點(diǎn)所對(duì)應(yīng)的收入水平，對(duì)（2.3）式進(jìn)行估計(jì)結(jié)果見表2。

根據(jù)表2，化肥-收入、農(nóng)藥-收入選用隨機(jī)效應(yīng)模型，農(nóng)膜-收入方程選用固定效應(yīng)模型。各方程系數(shù)均顯著，兩者都存在顯著的倒U型曲線關(guān)系，表3說明化肥-收入、農(nóng)藥-收入和農(nóng)膜-收入方程轉(zhuǎn)折點(diǎn)分別在人均收入為17289.72元、17405.46元和21799.63元達(dá)到。

表4為2012年人均農(nóng)業(yè)總產(chǎn)值前9位的地區(qū)以及全國的均值，除上海市、江蘇省人均農(nóng)業(yè)產(chǎn)值高于17300元外，其余地區(qū)均低于17300元，且全國均值僅9474.03元，遠(yuǎn)未達(dá)到化肥、農(nóng)藥和農(nóng)膜的倒U曲線轉(zhuǎn)折點(diǎn)。全國仍有29個(gè)個(gè)地區(qū)未達(dá)到化肥的倒U曲線轉(zhuǎn)折點(diǎn)，30個(gè)地區(qū)未達(dá)到農(nóng)藥的倒U曲線轉(zhuǎn)折點(diǎn)，31個(gè)地區(qū)未達(dá)到農(nóng)膜倒U曲線的轉(zhuǎn)折點(diǎn)。因此，從各地區(qū)的人均農(nóng)業(yè)總產(chǎn)值來看，中國農(nóng)業(yè)經(jīng)濟(jì)仍然處于 EKC 曲線的左半段，化肥、農(nóng)藥和農(nóng)膜的投入量隨著農(nóng)業(yè)的經(jīng)濟(jì)增長會(huì)呈現(xiàn)出上升趨勢，如不加以治理，面源污染問題仍將持續(xù)惡化。

2.對(duì)（2.4）式進(jìn)行估計(jì)與對(duì)比。

根據(jù)表5，加入控制變量后，化肥和農(nóng)藥投入的曲線特征沒有改變，農(nóng)膜曲線關(guān)系改變，未呈現(xiàn)倒U型曲線關(guān)系。（2.3）式和（2.4）式估計(jì)結(jié)果的對(duì)比見表6，控制變量的加入使Fer和Pes轉(zhuǎn)折點(diǎn)對(duì)應(yīng)的人均收入均上升。Fer的轉(zhuǎn)折點(diǎn)由17289.72元上升到22584.27元；Pes的轉(zhuǎn)折點(diǎn)由17405.46元上升到29850.75元。對(duì)照表4，實(shí)際人均收入水平與轉(zhuǎn)折點(diǎn)還有較大差距，全部地區(qū)都處于各方程曲線的左半段。

2.1人均耕地面積（Pl）。本文用人均耕地面積表征耕地資源的稀缺程度，結(jié)果Pl的系數(shù)在各方程中為正，表明耕地資源的稀缺程度與化肥、農(nóng)藥和農(nóng)膜使用顯著正相關(guān)。1995—2012年我國人均耕地面積從0.142公頃/人減少到0.139公頃/人，總體上處于減少趨勢。人均耕地的增加主要對(duì)土地經(jīng)營規(guī)模產(chǎn)生兩種影響：（1）耕地資源變豐富，其市場價(jià)格和流轉(zhuǎn)租金下降，較低的農(nóng)地流轉(zhuǎn)價(jià)格可能使得供需更易達(dá)到平衡，促進(jìn)土地經(jīng)營規(guī)模的發(fā)展。（2）耕地資源集約化化程度提高，減少了土地經(jīng)營規(guī)模所需的交易成本，更利于規(guī)模經(jīng)營主體的現(xiàn)代化生產(chǎn)，因此其進(jìn)一步流入土地?cái)U(kuò)大經(jīng)營規(guī)模的意愿越來越高。根據(jù)本文實(shí)證結(jié)果可知，隨著我國農(nóng)村居民家庭經(jīng)營人均耕地面積的增加，規(guī)模經(jīng)營成本下降，單位面積化肥施用量卻增加，化肥污染程度加深。

2.2勞均農(nóng)作物播種面積（Lp）。各投入方程的Lp系數(shù)均顯著，且農(nóng)均農(nóng)作物播種面積與化肥、農(nóng)藥和農(nóng)膜投入均為顯著負(fù)相關(guān)。2000年到2012年我國農(nóng)均農(nóng)作物播種面積呈現(xiàn)逐步增加趨勢，從2000年的0.48公頃/人增加到2012年的0.61公頃/人，增加了0.13公頃/人，增幅達(dá)27.08%?？梢?，隨著我國農(nóng)均播種面積的不斷增加，化肥、農(nóng)藥和農(nóng)膜的投入將下降。

2.3非農(nóng)就業(yè)水平（Na）。Na系數(shù)在化肥和農(nóng)膜投入方程中顯著（在農(nóng)藥投入方程中不顯著），且系數(shù)為正，表明非農(nóng)就業(yè)水平與化肥和農(nóng)膜投入顯著正相關(guān)。非農(nóng)就業(yè)水平越高，說明農(nóng)村勞動(dòng)力更普遍地向非農(nóng)產(chǎn)業(yè)轉(zhuǎn)移，從事農(nóng)業(yè)生產(chǎn)的家庭勞動(dòng)力減少，農(nóng)業(yè)產(chǎn)能不足。農(nóng)業(yè)生產(chǎn)的效益明顯低于第二、三產(chǎn)業(yè)，農(nóng)業(yè)勞動(dòng)力的機(jī)會(huì)成本越來越高。對(duì)于從事農(nóng)業(yè)生產(chǎn)的勞動(dòng)力來說，提高單位面積生產(chǎn)要素的投入（化肥、農(nóng)藥和農(nóng)膜等）是簡單可行的辦法，這不可避免的加劇農(nóng)業(yè)面源污染程度。我國非農(nóng)就業(yè)水平從2000年到2012年呈現(xiàn)上升趨勢，表明我國化肥和農(nóng)膜的投入將隨著非農(nóng)就業(yè)水平的上升而增加。

2.4農(nóng)業(yè)機(jī)械化水平（Am）。Am系數(shù)在化肥和農(nóng)藥投入方程中顯著（在農(nóng)膜投入方程中不顯著），且系數(shù)均為正，表明農(nóng)業(yè)機(jī)械化水平與化肥和農(nóng)藥投入顯著正相關(guān)。由于現(xiàn)階段我國農(nóng)業(yè)機(jī)械制造水平不高，農(nóng)林牧漁的各種動(dòng)力機(jī)械在使用過程中容易造成生產(chǎn)要素的流失和浪費(fèi)，導(dǎo)致農(nóng)民不得不追加施用化肥和農(nóng)藥。我國農(nóng)業(yè)機(jī)械化水平從2000年到2013年呈顯著上升趨勢。實(shí)證結(jié)果也說明，農(nóng)業(yè)機(jī)械化水平越高，化肥和農(nóng)藥的污染將顯著增加。

五、結(jié)語

本文的實(shí)證分析發(fā)現(xiàn)：

1.化肥、農(nóng)藥和農(nóng)膜與收入之間表現(xiàn)出穩(wěn)定的倒U型曲線，控制變量的加入對(duì)化肥和農(nóng)藥與收入曲線形狀沒有影響， 曲線轉(zhuǎn)折點(diǎn)對(duì)應(yīng)的收入水平提高，使農(nóng)膜與收入的曲線不再呈現(xiàn)倒U型；

2.各曲線轉(zhuǎn)折點(diǎn)所對(duì)應(yīng)的收入水平說明我國仍處于化肥、農(nóng)藥和農(nóng)膜EKC曲線的左半端，面源污染負(fù)荷強(qiáng)度會(huì)隨著收入的增長而增加；

3.各土地經(jīng)營規(guī)模指標(biāo)對(duì)不同面源污染負(fù)荷強(qiáng)度指標(biāo)產(chǎn)生不同的影響。化肥的投入與人均耕地面積、非農(nóng)就業(yè)水平、農(nóng)業(yè)機(jī)械化水平顯著正相關(guān)，與勞均農(nóng)作物播種面積顯著負(fù)相關(guān)；農(nóng)藥的投入與人均耕地面積和農(nóng)業(yè)機(jī)械化水平顯著正相關(guān)，與勞均農(nóng)作物播種面積顯著負(fù)相關(guān)；

參考文獻(xiàn)：

[1]DENNISLC，PETER JV，KEITHL. Modeling non-point source pollutants in vadose zone with GIS[J].Environmental Science and Technology，1997，31（8）：2157-217.

[2]張維理，武淑霞，冀宏杰，等.中國農(nóng)業(yè)面源污染形勢估計(jì)及控制對(duì)策：21世紀(jì)初期中國農(nóng)業(yè)面源污染的形勢估計(jì)[J]中國農(nóng)業(yè)科學(xué)，2004，37（7）：1008-1017.

[3]張?zhí)伊?農(nóng)業(yè)部副部長：農(nóng)業(yè)已超工業(yè)成最大面源污染產(chǎn)業(yè)[OL]http：//env.people.com.cn/n/2015/0415/c1010-26846588.html，20150414國新辦發(fā)布會(huì)

[4]董雪嬌，湯惠君.國內(nèi)外土地經(jīng)營規(guī)模述評(píng)[J].中國農(nóng)業(yè)資源與區(qū)劃，2015，（03）： 62～71.

[5] DINDA S.Environmental Kuznets curve hypothesis：A survey[J].Ecological Economics，2004，49：431-455.

[6]MANAGI S.Are there increasing returns to pollution abatement？ Empirical analytics of the environmental Kuznets curve in pesticides[I].Ecological Economics，2006，58（3）：617-636.

[7]TSUZUKI Y.An index directly indicates land-based pollutant load contribution of domestic wastewater to the water pollution and its application[J].Science of the Total Environment，2006，370：425-440.

[8]李海鵬，張俊飚.中國農(nóng)業(yè)面源污染與經(jīng)濟(jì)發(fā)展關(guān)系的實(shí)證研究[J].長江流域資源與環(huán)境， 2009，18（6）：585-590.

[9]張暉，胡浩.農(nóng)業(yè)面源污染的環(huán)境庫茲涅茨曲線驗(yàn)證——基于江蘇省時(shí)序數(shù)據(jù)的分析[J].中國農(nóng)村經(jīng)濟(jì)， 2009（04）：48-53.

[10]杜江，羅珺.農(nóng)業(yè)經(jīng)濟(jì)增長與污染性要素投入——基于簡約式及結(jié)構(gòu)式模型的實(shí)證分析[J].經(jīng)濟(jì)評(píng)論，2013（03）：56-65.

作者簡介：劉韜（1990-）男，湖南邵東人，中南民族大學(xué)土地資源管理專業(yè)研究生。