農(nóng)業(yè)增長與環(huán)境污染的動態(tài)分析*
——基于2006—2015年面板數(shù)據(jù)的驗(yàn)證

高如夢，杜 江，李曉濤

(武漢輕工大學(xué)經(jīng)濟(jì)與管理學(xué)院，武漢 430023)

0 引言

“三農(nóng)”問題關(guān)系到國計(jì)民生與社會的和諧穩(wěn)定發(fā)展，我國歷年來對農(nóng)業(yè)農(nóng)村發(fā)展問題都高度重視?！爸醒胍惶栁募边B續(xù)14年聚焦“三農(nóng)”問題， 2017年該文件指出提高糧食生產(chǎn)能力將繼續(xù)成為農(nóng)業(yè)經(jīng)濟(jì)發(fā)展的基本目標(biāo)?！笆晃濉薄ⅰ笆濉逼陂g，我國農(nóng)業(yè)農(nóng)村經(jīng)濟(jì)得到了快速的發(fā)展，生產(chǎn)實(shí)現(xiàn)了“十連增”，農(nóng)業(yè)總產(chǎn)值從2006年的2.152 228萬億元增加到2015年的3.133 746萬億元(以2006年為基期的不變價格計(jì)算)，年均增長幅度為4.56%，為國民經(jīng)濟(jì)持續(xù)穩(wěn)定發(fā)展做出了重要貢獻(xiàn)。在農(nóng)業(yè)總產(chǎn)值增加的同時，農(nóng)業(yè)面源污染也在不斷加劇嚴(yán)重影響了我國農(nóng)業(yè)的可持續(xù)發(fā)展。目前，我國已經(jīng)成為世界上最大的農(nóng)藥、化肥生產(chǎn)國和消費(fèi)國，塑料薄膜的使用量也居世界首位。其中，農(nóng)藥使用量由2006年的153.71萬t增加到2015年的178.12萬t，使用強(qiáng)度由2006年的11.82kg/hm2增加到2015 年的13.70kg/hm2，年均增長幅度為1.59%?；适褂昧坑?006年的4 927.7萬t增加到2015年的6 022.5萬t，使用強(qiáng)度由2006年的378.93kg/hm2增加到2015 年的463.12kg/hm2，年均增長幅度為2.22%。塑料薄膜使用量由2006年的184.55萬t增加到2015年的260.36萬t，使用強(qiáng)度由2006年的14.19kg/hm2增加到2015 年的20.02kg/hm2，年均增長幅度為4.11%。據(jù)有關(guān)資料顯示， 2015年我國水稻、玉米、小麥三大糧食作物化肥利用率為35.2%、農(nóng)藥利用率為36.6%、塑料薄膜回收率為30%，約有70%左右未得到利用。未利用的部分在很大程度上會造成大氣質(zhì)量下降、水體富營養(yǎng)化、土壤肥力下降，這是造成我國農(nóng)業(yè)面源污染的一個主要原因。而歐美等發(fā)達(dá)國家小麥、玉米等糧食作物的農(nóng)藥、化肥等利用率在50%～60%，比我國高15～25個百分點(diǎn)，可見我國在農(nóng)藥、化肥、塑料薄膜的使用上存在過量、不合理等問題，這是當(dāng)下我們面臨的重要問題?！笆濉逼陂g我國農(nóng)業(yè)生產(chǎn)重心由“增產(chǎn)”逐漸轉(zhuǎn)移到“增產(chǎn)”和“質(zhì)量”兩個重心上來，在提高農(nóng)業(yè)增長的基礎(chǔ)上開始更多的關(guān)注農(nóng)業(yè)可持續(xù)發(fā)展。

2008—2010， 2013—2015年的“中央一號文件”都明確提到農(nóng)業(yè)化肥污染治理問題，并且“十二五”規(guī)劃將防治農(nóng)業(yè)面源污染列為“十二五”重點(diǎn)工作。2015年“中央一號文件”對“加強(qiáng)農(nóng)業(yè)生態(tài)治理”作出專門部署，強(qiáng)調(diào)要加強(qiáng)農(nóng)業(yè)面源污染治理，并且當(dāng)年政府工作報告也提出了加強(qiáng)農(nóng)業(yè)面源污染治理的重大任務(wù)。目前我國農(nóng)業(yè)生產(chǎn)面臨的最大問題在于如何權(quán)衡農(nóng)業(yè)發(fā)展、資源利用與生態(tài)環(huán)境三者的關(guān)系，因此研究農(nóng)業(yè)增長與環(huán)境污染的關(guān)系，對于實(shí)現(xiàn)“十三五”時期產(chǎn)業(yè)結(jié)構(gòu)升級和“一控兩減三基本”的目標(biāo)具有重要的現(xiàn)實(shí)意義。

農(nóng)業(yè)增長的源泉主要來自于傳統(tǒng)要素的投入，其中農(nóng)藥、化肥投入對農(nóng)業(yè)經(jīng)濟(jì)增長具有重要的貢獻(xiàn)[1](lin， 1992)。農(nóng)藥的合理使用可以有效預(yù)防各種病蟲害，提高農(nóng)產(chǎn)品的產(chǎn)量； 化肥的合理使用可以有效保持土壤養(yǎng)分與肥力。反之不僅增加農(nóng)業(yè)生產(chǎn)成本，而且會致使農(nóng)產(chǎn)品產(chǎn)量下降，成為農(nóng)業(yè)面源污染的主要來源。楊晶[2](2016)在研究農(nóng)業(yè)面源污染與經(jīng)濟(jì)增長的關(guān)系時發(fā)現(xiàn)四川省農(nóng)田污染主要來源于農(nóng)藥化肥的不合理使用。目前我國農(nóng)業(yè)整體面臨著農(nóng)藥使用強(qiáng)度大、施肥方法不科學(xué)、利用率低等問題。因此，減少農(nóng)藥、化肥使用量有利于降低農(nóng)業(yè)面源污染，對于促進(jìn)農(nóng)業(yè)綠色可持續(xù)發(fā)展，具有顯著的經(jīng)濟(jì)效益、環(huán)境效益和社會效益[3](金書秦等， 2016)。

以往的研究常常借助環(huán)境庫茨涅茲曲線(EKC)來分析經(jīng)濟(jì)增長與環(huán)境污染的關(guān)系。國內(nèi)外學(xué)者運(yùn)用不同國家或部門的數(shù)據(jù)，采用不同的方法，分別驗(yàn)證農(nóng)藥、化肥等化學(xué)品投入與收入之間EKC關(guān)系的存在性、探討拐點(diǎn)的位置，并分析EKC 曲線形成的內(nèi)外在機(jī)理。Crossman 和Krueger(1991)[4]對42 個國家的環(huán)境與經(jīng)濟(jì)數(shù)據(jù)進(jìn)行分析后發(fā)現(xiàn)，經(jīng)濟(jì)增長與環(huán)境污染之間存在類似于倒U 型的曲線關(guān)系。Kuznets在1995年發(fā)現(xiàn)收入分配與經(jīng)濟(jì)增長之間存在類似于倒U形的曲線關(guān)系，并將其命名為庫茨涅茲曲線。Panayotou[5](1996)借用1995年庫茲涅茨界定的人均收入與收入不均等之間的倒 U型曲線，首次將這種環(huán)境質(zhì)量與人均收入間的關(guān)系稱為環(huán)境庫茲涅茨曲線(EKC)。EKC理論假說提出后，相關(guān)領(lǐng)域?qū)嵶C研究得到不斷深化，Shafik 和 Bandyopadhyay[6]、Selden和Song[7]、Crossman 和Krueger[8]的研究都在一定程度上支持了EKC的存在性，并且指出發(fā)達(dá)國家和發(fā)展中國家的庫茨涅茲曲線有所差異，這便引起了國內(nèi)相關(guān)學(xué)者的重視。

國內(nèi)關(guān)于農(nóng)業(yè)增長與環(huán)境污染的研究主要集中在兩個方面。一是定性研究，主要體現(xiàn)在對農(nóng)業(yè)面源污染現(xiàn)狀、成因的分析和相關(guān)對策的研究。二是定量研究，通過理論和實(shí)證分析來探討農(nóng)業(yè)EKC的存在性。王建兵等(2008)[9]、李秀芬等(2010)[10]分析了我國農(nóng)業(yè)面源污染的現(xiàn)狀。梁流濤等(2010)[11]、杜江等(2013)[12]對農(nóng)業(yè)面源污染的成因進(jìn)行了分析。饒靜等(2011)[13]對農(nóng)業(yè)面源污染的現(xiàn)狀和成因進(jìn)行了研究。朱兆良等(2005)[14]和崔鍵等(2006)[15]對農(nóng)業(yè)面源污染的防治對策進(jìn)行了研究，并從政策制定、環(huán)境立法、技術(shù)體系等方面提出了政策建議。以上學(xué)者是基于定性角度研究環(huán)境污染與經(jīng)濟(jì)增長的關(guān)系，定性研究由于無數(shù)據(jù)的支撐，所以相關(guān)結(jié)論也缺乏說服力。張暉等[16](2009)研究發(fā)現(xiàn)江蘇省經(jīng)濟(jì)發(fā)展水平與過剩氮污染之間存在顯著的倒U型曲線關(guān)系。劉揚(yáng)等[17](2009)和李太平等[18](2011)以化肥為例、李海鵬等[19](2009)以化肥、農(nóng)藥投入密度和畜禽糞尿排泄物密度為例、杜江等[20](2013)和于冀等[21](2016)以化肥、農(nóng)藥、農(nóng)用塑料薄膜為例、沈能等[22](2016)以農(nóng)藥投入為例，研究發(fā)現(xiàn)農(nóng)業(yè)增長與環(huán)境污染呈現(xiàn)倒U型曲線關(guān)系，但不同區(qū)域的EKC曲線轉(zhuǎn)折點(diǎn)和所處的階段不同。對于同一污染物研究區(qū)域、測算方法不同時，結(jié)果有所差異。杜江等[23](2009)研究得出人均農(nóng)業(yè)產(chǎn)值與農(nóng)藥投入間存在倒 U 型關(guān)系，與化肥投入間存在倒 N 型關(guān)系。沈能等[24](2015)從生產(chǎn)率的角度出發(fā)，研究發(fā)現(xiàn)農(nóng)業(yè)額增長與環(huán)境生產(chǎn)率之間存在穩(wěn)定的正 U 型曲線關(guān)系。馬乃毅等[25](2016)研究發(fā)現(xiàn)西部地區(qū)化肥的排放和經(jīng)濟(jì)增長之間存在“U”形關(guān)系，而農(nóng)藥和農(nóng)膜的排放與經(jīng)濟(jì)增長呈現(xiàn)線性遞增關(guān)系。以上學(xué)者主要是基于定量角度研究環(huán)境污染與經(jīng)濟(jì)增長的關(guān)系。目前，關(guān)于農(nóng)業(yè)面源污染與經(jīng)濟(jì)增長的研究尚存在一定的不足。一是研究指標(biāo)沒有統(tǒng)一的選取標(biāo)準(zhǔn)，指標(biāo)選取具有多樣性，得出的結(jié)論不具統(tǒng)一性。二是樣本的選取隨意性較大，得出的結(jié)論不具代表性。三是研究范圍主要集中在某個省份或地區(qū)，針對全國性的研究較少且一般選取的是時間序列數(shù)據(jù)，面板數(shù)據(jù)的研究相對匱乏。四是在以往的研究中普遍采用農(nóng)藥、化肥等投入總量，忽略了并不是所有要素的投入都會帶來環(huán)境污染，只是要素的流失部分才會導(dǎo)致農(nóng)業(yè)面源污染。該文采用各個要素的流失量作為污染指標(biāo)，試圖從化肥、農(nóng)藥、農(nóng)用塑料薄膜出發(fā)，利用2006—2015年省際面板數(shù)據(jù)進(jìn)行理論和實(shí)證分析來探討我國“十一五”、“十二五”時期農(nóng)業(yè)面源污染與經(jīng)濟(jì)增長的關(guān)系。該研究在一定程度上能夠?yàn)槲覈r(nóng)業(yè)可持續(xù)發(fā)展和環(huán)境治理提出相關(guān)政策建議，有利于我國“十三五”時期農(nóng)業(yè)發(fā)展目標(biāo)的實(shí)現(xiàn)，同時也是相關(guān)領(lǐng)域已有研究的一個補(bǔ)充。

該文試圖回答以下3個問題：一是農(nóng)業(yè)增長與面源污染是怎樣的關(guān)系。二是要什么樣的經(jīng)濟(jì)增長。三是如何實(shí)現(xiàn)在經(jīng)濟(jì)增長中控制或者消除農(nóng)業(yè)面源污染。在此基礎(chǔ)上，提出未來我國化肥、農(nóng)藥、農(nóng)用塑料薄膜科學(xué)使用的對策與建議。

1 理論框架與研究設(shè)想

在研究農(nóng)業(yè)EKC時，大多數(shù)學(xué)者往往采用簡約式模型來考查環(huán)境污染與經(jīng)濟(jì)增長之間的關(guān)系，理論模型如下：

yit=αit+β1xit+β2xit2+β3xit3+β4zit+εit

(1)

式(1)中，y和x分別表示環(huán)境質(zhì)量與收入變量，z為其他影響環(huán)境質(zhì)量的變量，即控制變量；i表示國家或地區(qū)，t表示時間；βk表示待估參數(shù)(k=1， 2， 3， 4)；εit表示隨機(jī)干擾項(xiàng)。通過實(shí)證分析可能得到以下7種x與y的關(guān)系，見表1。

表1 實(shí)證分析可能結(jié)果

可能情況系數(shù)滿足條件曲線形狀拐點(diǎn)1β1=β2=β3=0沒有必然聯(lián)系無2 β1>0,β2=β3=0單調(diào)遞增關(guān)系無3β1<0,β2=β3=0單調(diào)遞減關(guān)系無4β1>0,β2<0,β3=0倒U型曲線(EKC)-β1/(2β2)5β2>0,β1<0,β3=0U型關(guān)系β1/(2β2)6β1>0,β2<0,β3>0N型曲線-β2/(3β3)7β1<0,β2>0,β3<0倒N型曲線β2/(3β3)

2 模型構(gòu)建、指標(biāo)選取與數(shù)據(jù)說明

2.1 模型構(gòu)建

借鑒已有的研究，該文采用簡約式模型來探討農(nóng)業(yè)環(huán)境污染與經(jīng)濟(jì)增長的關(guān)系。建立如下方程：

yit=αit+β1agricultureit+β2(agricultureit)2+β3(agricultureit)3+∑Di+εit

(2)

yit分別包含fertilizerit、pesticideit、plasticit，分別表示t時期i省份化肥流失量、農(nóng)藥流失量、塑料薄膜殘留量。agricultureit表示t時期i省份農(nóng)業(yè)總產(chǎn)值，以下變量machineit、irrigateit、urbanit為控制變量Di，分別表示t時期i省份機(jī)械擁有量、有效灌溉率、城鎮(zhèn)化水平。

實(shí)證中的污染方程有3個，分別是以化肥、農(nóng)藥、塑料薄膜衡量。在不考慮控制變量的情況下估計(jì)方程，初步探討污染性要素投入與農(nóng)業(yè)增長之間的關(guān)系，在方程的估計(jì)過程中分別對包括和不包括agri3兩種情況進(jìn)行檢驗(yàn)，最終根據(jù)豪斯曼檢驗(yàn)結(jié)果選擇最為合適的模型在文章中展示。

2.2 指標(biāo)選取

該文主要選取2006—2015年全國31個省(省份、自治區(qū))的面板數(shù)據(jù)，探討“十一五”以來農(nóng)業(yè)發(fā)展是否已經(jīng)突破拐點(diǎn)處于庫茨涅茲曲線下降階段，檢驗(yàn)農(nóng)業(yè)各項(xiàng)減排技術(shù)是否得到推廣、各項(xiàng)政策是否得到落實(shí)。該文數(shù)據(jù)來源于《中國統(tǒng)計(jì)年鑒》《中國環(huán)境統(tǒng)計(jì)年鑒》《中國農(nóng)村統(tǒng)計(jì)年鑒》。采用的主要經(jīng)濟(jì)指標(biāo)為：

(1)選取各省農(nóng)業(yè)總產(chǎn)值作為經(jīng)濟(jì)增長的指標(biāo)。此處農(nóng)業(yè)為狹義農(nóng)業(yè)，即種植業(yè)。該文選取農(nóng)業(yè)總產(chǎn)值作為經(jīng)濟(jì)增長的指標(biāo)，單位為億元，并以2006年為基期，對后續(xù)年份農(nóng)業(yè)總產(chǎn)值進(jìn)行了處理，以消除通貨膨脹的影響。

(2)衡量污染性要素投入的指標(biāo)為化肥流失量、農(nóng)藥流失量、塑料薄膜殘留量，單位均為萬t。其中化肥流失量=化肥施用總量(折純量)×(1-化肥利用率)，農(nóng)藥流失量=農(nóng)藥施用總量(折純量)×(1-農(nóng)藥利用率)，塑料薄膜殘留量=塑料薄膜使用總量(折純量)×(1-塑料薄膜回收率)。根據(jù)已有文獻(xiàn)(邵振潤， 2014)[26]，農(nóng)藥利用率為30%，化肥利用率為35%，塑料薄膜回收率為30%，考慮到流失系數(shù)的影響因素較多并難以獲得，根據(jù)相關(guān)文獻(xiàn)的處理方法，所以該文也對以上比例進(jìn)行換算進(jìn)而得出各項(xiàng)指標(biāo)。

(3)控制變量主要有機(jī)械擁有量(萬kW/h，mac)，有效灌溉率(%，irr)，城鎮(zhèn)化水平(%，urban)。有效灌溉率=有效灌溉面積/當(dāng)年總耕地面積，城鎮(zhèn)化水平=城鎮(zhèn)人口/總?cè)丝?。機(jī)械擁有量和有效灌溉率是農(nóng)業(yè)生產(chǎn)條件的表征，是技術(shù)進(jìn)步的體現(xiàn)。技術(shù)進(jìn)步會帶來先進(jìn)的生產(chǎn)設(shè)備和生產(chǎn)技術(shù)，可以極大程度的改善農(nóng)業(yè)生產(chǎn)條件。城鎮(zhèn)化的發(fā)展是農(nóng)村人口向城市聚集和農(nóng)村地區(qū)逐漸轉(zhuǎn)變?yōu)槌鞘械貐^(qū)的復(fù)雜過程，表現(xiàn)為二、三產(chǎn)業(yè)向城市集中。城市通過增長極作用將經(jīng)濟(jì)增長的利益擴(kuò)散到農(nóng)村地區(qū)，為農(nóng)業(yè)的可持續(xù)發(fā)展提供大量物質(zhì)基礎(chǔ)與技術(shù)支撐，城市化的健康發(fā)展有利于降低農(nóng)業(yè)生產(chǎn)對環(huán)境造成的污染。由于機(jī)械擁有量、有效灌溉率、城鎮(zhèn)化水平在某種程度上都會影響農(nóng)業(yè)面源污染，所以選取為控制變量。但良種率、測土配方作為環(huán)保農(nóng)業(yè)技術(shù)中的關(guān)鍵因素，由于數(shù)據(jù)無法精確獲得，在此不再考慮。

2.3 數(shù)據(jù)說明

在進(jìn)行實(shí)證分析之前首先運(yùn)用stata13對上述相關(guān)數(shù)據(jù)做描述性統(tǒng)計(jì)。結(jié)果見表2。

表2 變量描述性統(tǒng)計(jì)分析結(jié)果

變量變量描述均值標(biāo)準(zhǔn)差最小值最大值樣本fertilizer化肥流失量116.82293.8942.86465.465N=310agriculture農(nóng)業(yè)總產(chǎn)值847.616663.79631.83 234.435N=310machine機(jī)械擁有量2 877.8372 769.16495.3133 530N=310urban城鎮(zhèn)化水平51.42914.61121.1389.6N=310pesticide農(nóng)藥流失量3.9023.1430.0612.142|N=310plastic塑料薄膜殘留量5.0374.5930.02824.047N=310irrigate有效灌溉率54.34727.25415.057322.689N=310

3 實(shí)證結(jié)果與分析

3.1 實(shí)證結(jié)果

(1)方程估計(jì)結(jié)果

表3 方程估計(jì)結(jié)果

pesticidefertilizerplasticFEREFEREFEREagriculture 0.00421??? 0.00498???0.120???0.142??? 0.0120??? 0.0106???(0.000485)(0.000454)(0.0112)(0.0102)(0.000957)(0.000846)agriculture2-0.000000964???-0.00000109???-0.0000184???-0.0000220???-0.00000236???-0.00000198???(0.000000127)(0.000000124)(0.00000293)(0.00000282)(0.000000251)(0.000000235)_cons 1.454???0.32836.57???27.40???-2.405???-2.030??(0.282)(0.378)(6.495)(8.042)(0.556)(0.642)調(diào)整R20.70680.7210.83960.83860.44290.4553F38.25—101.02—94.48—wald X2—130.9—344.45—214.15hausman檢驗(yàn)7.26(0.0265) 15.00(0.0006)6.51(0.0386)模型選擇FEFEFE曲線形式倒U型倒U型倒U型轉(zhuǎn)折點(diǎn)2 1833 2602 542 注:Standard errors in parentheses,? p<0.05,??p<0.01,???p<0.001。調(diào)整R2值用來衡量方程整體的擬合優(yōu)度, F值和waldX2用來檢驗(yàn)回歸方程的顯著性

在考慮和不考慮三次方項(xiàng)的情況下對方程分別做出實(shí)證分析，由于三次方項(xiàng)回歸系數(shù)均不顯著，所以初步判斷不存在N型或倒N型曲線關(guān)系。為了避免文章篇幅過長，在此只列出隨機(jī)效應(yīng)模型、固定效應(yīng)模型估計(jì)結(jié)果，在后續(xù)的穩(wěn)健性檢驗(yàn)中也是如此。方程隨機(jī)效應(yīng)模型、固定效應(yīng)模型估計(jì)結(jié)果，見表3。

根據(jù)表3可知Hausman檢驗(yàn)結(jié)果均小于0.1，所以拒絕原假設(shè)，均選用固定效應(yīng)模型?；?、農(nóng)藥、塑料薄膜與農(nóng)業(yè)總產(chǎn)值方程各回歸系數(shù)均顯著，兩者間存在顯著的倒U 型曲線關(guān)系。轉(zhuǎn)折點(diǎn)分別在農(nóng)業(yè)總產(chǎn)值為3 260億元、2 183億元、2 542億元達(dá)到。從樣本總產(chǎn)值整體來看，從2005年開始山東省已經(jīng)達(dá)到農(nóng)藥方程的拐點(diǎn)，從2010年開始達(dá)到塑料薄膜的拐點(diǎn)，河南省在2009年達(dá)到農(nóng)藥的拐點(diǎn)，在2012年達(dá)到塑料薄膜的拐點(diǎn)，但是均未達(dá)到化肥的拐點(diǎn)，這意味著“十三五”時期這兩個省份化肥的施用量將會繼續(xù)增加，農(nóng)藥、塑料薄膜的使用量將會降低，處于庫茨涅茲曲線下降階段。其他省份隨著經(jīng)濟(jì)的增長化肥、農(nóng)藥、塑料薄膜的使用量會繼續(xù)增加，環(huán)境質(zhì)量將未能得到有效改善。轉(zhuǎn)折點(diǎn)的出現(xiàn)是受多種因素綜合影響，山東省農(nóng)業(yè)增加值長期穩(wěn)居我國第一位，且山東、河南兩地區(qū)農(nóng)業(yè)總產(chǎn)值在生產(chǎn)總值中所占的比重、人均農(nóng)業(yè)總產(chǎn)值均位于前列。農(nóng)民生活水平更加富裕，農(nóng)民整體環(huán)保意識更強(qiáng)，政府出臺相關(guān)環(huán)保政策的實(shí)施效果，農(nóng)業(yè)技術(shù)水平和城鎮(zhèn)化水平也可能是導(dǎo)致山東、河南先于糧食主產(chǎn)區(qū)的其他省份出現(xiàn)轉(zhuǎn)折點(diǎn)的原因。作為經(jīng)濟(jì)發(fā)達(dá)地區(qū)的上海、浙江、廣東等省份，雖然人均總產(chǎn)值較高，但其人均農(nóng)業(yè)總產(chǎn)值水平較低。這些省份主要是以發(fā)展工業(yè)為主，對農(nóng)業(yè)的重視程度較低，使得農(nóng)民整體環(huán)保意識薄弱，各項(xiàng)減排政策落實(shí)不到位，農(nóng)業(yè)生產(chǎn)帶來了更多的環(huán)境污染，所以達(dá)到拐點(diǎn)還需一段時間。作為經(jīng)濟(jì)發(fā)展水平較落后的西部地區(qū)，農(nóng)民收入的增加主要來源于農(nóng)業(yè)。農(nóng)民為了增加其收入，在農(nóng)業(yè)生產(chǎn)中大量增加化肥、農(nóng)藥等物質(zhì)的投入，機(jī)械化水平的落后也會導(dǎo)致環(huán)境污染進(jìn)一步增加，所以其達(dá)到拐點(diǎn)還需很長一段時間。

(2)穩(wěn)健性檢驗(yàn)結(jié)果

表4 穩(wěn)健性檢驗(yàn)結(jié)果

pesticidefertilizerplasticFEREFEREFEREagriculture0.00406???0.00545???0.125???0.160???0.00911???0.00910???(0.00085)(0.000587)(0.0193)(0.0121)(0.00166)(0.00103)agriculture2-0.000000936???-0.00000114???-0.0000190???-0.0000240???-0.00000193???-0.00000179???(0.000000165)(0.000000133)(0.00000376)(0.00000287)(0.000000323)(0.000000247)machine0.00005820.0000678?0.00192?0.00222??0.00008330.000124(0.0000327)(0.000033)(0.000744)(0.000743)(0.0000639)(0.0000643)irrigate0.0004940.000916-0.0416-0.0398-0.00503 -0.00407(0.00147)(0.0015)(0.0335)(0.0338)(0.00288)(0.00292)Urban-0.00315-0.0246?-0.256-0.867???0.0578?0.0370?(0.0144)(0.00978)(0.328)(0.200)(0.0282)(0.0171)_cons1.510???1.011?42.94???54.49???-3.397???-3.024???(0.399)(0.461)(9.076)(9.612)(0.780)(0.825)調(diào)整R20.7260.72470.86430.86160.632206433F15.96—43.51—40.35—wald X2—143.95—422.45—228.69hausman檢驗(yàn)6.65(0.0557)4.96(0.0219) 14.24(0.0066)模型選擇FEFEFE曲線形式倒U型倒U型倒U型轉(zhuǎn)折點(diǎn)2 390.353 333.332 360 注:Standard errors in parentheses,? p<0.05,??p<0.01,???p<0.001

該文選用增加對被解釋變量有一定影響的控制變量來檢驗(yàn)?zāi)Ｐ偷姆€(wěn)健性。加入控制變量后，對方程進(jìn)行回歸，以便探討EKC模型的穩(wěn)健性，即控制變量的加入是否會對污染性要素投入與農(nóng)業(yè)增長的關(guān)系產(chǎn)生影響，穩(wěn)健性檢驗(yàn)結(jié)果見表4。

根據(jù)表4可知Hausman檢驗(yàn)結(jié)果均小于0.1，所以拒絕原假設(shè)，均選用固定效應(yīng)模型?；?、農(nóng)藥、塑料薄膜與農(nóng)業(yè)總產(chǎn)值方程各回歸系數(shù)均顯著，兩者間存在顯著的倒U 型曲線關(guān)系，轉(zhuǎn)折點(diǎn)分別在3 333.33億元、2 390.35億元、2 360億元時。由此可見，加入控制變量不影響模型的顯著性，只是改變了轉(zhuǎn)折點(diǎn)的出現(xiàn)。

3.2 結(jié)果分析

加入控制變量后，化肥、農(nóng)藥、塑料薄膜與農(nóng)業(yè)總產(chǎn)值方程系數(shù)均顯著，兩者間仍存在顯著的倒U 型曲線關(guān)系，符合EKC模型?？梢娀?、農(nóng)藥、塑料薄膜與農(nóng)業(yè)總產(chǎn)值EKC模型穩(wěn)健性較好，不易受其他因素的影響。此時轉(zhuǎn)折點(diǎn)分別在3 333.33億元、2 390.35億元、2 360億元達(dá)到。相比于未加控制變量時，化肥、農(nóng)藥均推遲了轉(zhuǎn)折點(diǎn)的出現(xiàn)，塑料薄膜則提前了出現(xiàn)了轉(zhuǎn)折點(diǎn)。按照我國2006—2015年農(nóng)業(yè)產(chǎn)值空間分布來看，除了山東、河南部分年份到達(dá)了拐點(diǎn)之外，其余所有省份農(nóng)業(yè)總產(chǎn)值均低于拐點(diǎn)值。表明我國整體農(nóng)業(yè)發(fā)展仍然處于農(nóng)業(yè)環(huán)境庫涅茲曲線的上升階段，未來的農(nóng)業(yè)增長勢必會帶來大量的農(nóng)業(yè)污染問題。目前農(nóng)民增收、糧食增長依然是奮斗目標(biāo)，加上我國農(nóng)業(yè)生產(chǎn)正處于轉(zhuǎn)型時期，各項(xiàng)技術(shù)措施并不成熟、各項(xiàng)政策尚未完善。未來“十三五”時期的農(nóng)業(yè)發(fā)展將會繼續(xù)依賴傳統(tǒng)要素的投入，農(nóng)業(yè)面源污染將會繼續(xù)。

下面探討一下控制變量的加入對農(nóng)業(yè)污染排放的影響。有效灌溉率這一變量的系數(shù)為正但未通過顯著性檢驗(yàn)，不能說明有效灌溉是否可以減少農(nóng)業(yè)污染。在對農(nóng)作物施肥后，及時有效灌溉一方面可以促進(jìn)化肥的吸收，另一方面也會導(dǎo)致化肥的流失，造成農(nóng)業(yè)環(huán)境污染。農(nóng)業(yè)機(jī)械化水平的提高會顯著增加化肥、農(nóng)藥的使用量，加劇面源污染程度。這一結(jié)果與農(nóng)業(yè)技術(shù)進(jìn)步有利于降低農(nóng)業(yè)面源污染理論相矛盾，可能的原因是中國農(nóng)戶的經(jīng)營規(guī)模小、各項(xiàng)技術(shù)推廣并未得到很好的應(yīng)用，使得農(nóng)業(yè)技術(shù)進(jìn)步環(huán)境效應(yīng)并未得到充分發(fā)揮。城鎮(zhèn)化水平可以減少農(nóng)藥化肥的使用，從而降低農(nóng)業(yè)污染，這與杜江等(2009)[27]、曹大宇(2011)[28]研究結(jié)論相一致。該文發(fā)現(xiàn)城鎮(zhèn)化水平的提高會增加塑料薄膜的使用，我國正處于城鎮(zhèn)化的快速發(fā)展階段，雖然在某種程度上減少了農(nóng)藥、化肥的使用，但是卻增加了固體廢棄物的排放。

4 政策建議

理論上農(nóng)藥、化肥、塑料薄膜的使用均能顯著的促進(jìn)農(nóng)業(yè)經(jīng)濟(jì)增長，但我們不能繼續(xù)走“先污染，后治理”的農(nóng)業(yè)發(fā)展老路子，要走綠色、循環(huán)、可持續(xù)發(fā)展的農(nóng)業(yè)現(xiàn)代化新道路。積極采取相關(guān)措施去改善農(nóng)業(yè)生產(chǎn)條件，降低農(nóng)業(yè)發(fā)展帶來的污染。

“十一五”、“十二五”時期我國雖然已經(jīng)高度重視農(nóng)業(yè)面源污染問題，但是目前污染現(xiàn)狀還是未能得到很大程度的緩解，依然處于環(huán)境庫茨涅茲曲線上升階段。要想促使U型曲線拐點(diǎn)提前出現(xiàn)，政府應(yīng)針對我國不同地區(qū)的實(shí)際情況，制定差異化的農(nóng)業(yè)產(chǎn)業(yè)和環(huán)保政策，通過多種方式的靈活應(yīng)用加快實(shí)現(xiàn)農(nóng)業(yè)由傳統(tǒng)生產(chǎn)方式向現(xiàn)代生產(chǎn)方式的轉(zhuǎn)變。該文給出以下幾點(diǎn)政策啟示。

(1)要借鑒中部地區(qū)山東省、河南省的成功發(fā)展經(jīng)驗(yàn)。其次是農(nóng)業(yè)研究要有新突破，不斷加大科技攻關(guān)，開發(fā)減少農(nóng)業(yè)廢棄物生成等農(nóng)業(yè)環(huán)保技術(shù)。積極支持農(nóng)田水利建設(shè)，研發(fā)和推廣節(jié)約型施肥、施藥、灌溉農(nóng)業(yè)。

(2)政策制定要具有針對性而不應(yīng)停留在維持和補(bǔ)救上。首先，根據(jù)生態(tài)系統(tǒng)內(nèi)物質(zhì)循環(huán)和能量轉(zhuǎn)化規(guī)律，在生產(chǎn)過程中建立和諧互補(bǔ)的耕作制度，把化肥和農(nóng)藥的投入降低到最低水平。其次，建立高效、環(huán)境友好型農(nóng)業(yè)技術(shù)與推廣體系，加大農(nóng)業(yè)先進(jìn)技術(shù)、測土配方施肥等環(huán)保農(nóng)業(yè)技術(shù)的推廣，使農(nóng)民真正意識到過量使用化肥、農(nóng)藥、塑料薄膜所帶來的不利影響。

(3)淘汰性能落后的機(jī)械，引進(jìn)先進(jìn)的農(nóng)業(yè)機(jī)械。目前，我國農(nóng)業(yè)基礎(chǔ)設(shè)施薄弱，農(nóng)用機(jī)械綜合生產(chǎn)能力弱，機(jī)械化水平有待進(jìn)一步提升?，F(xiàn)代化農(nóng)業(yè)需要更高水平的機(jī)械，農(nóng)機(jī)產(chǎn)品需要與現(xiàn)代先進(jìn)的信息技術(shù)、電子技術(shù)、計(jì)算機(jī)技術(shù)相結(jié)合，推動農(nóng)業(yè)機(jī)械設(shè)備向智能化一體化發(fā)展?！笆濉睍r期是全面建成小康社會的決勝階段，應(yīng)該將農(nóng)業(yè)重心轉(zhuǎn)移到依靠科技創(chuàng)新和提高勞動者素質(zhì)上來。