農(nóng)業(yè)經(jīng)濟增長對農(nóng)業(yè)面源污染影響的實證研究

黃佳琦 劉舜佳

摘要根據(jù)1991—2014年我國31省（市、自治區(qū)）3項環(huán)境污染指標的面板數(shù)據(jù)，就農(nóng)業(yè)經(jīng)濟增長與農(nóng)業(yè)環(huán)境污染之間的關系，以環(huán)境3效應為基礎進行EKC檢驗。結(jié)果表明：農(nóng)業(yè)經(jīng)濟增長與農(nóng)業(yè)環(huán)境污染之間的關系正好與傳統(tǒng)環(huán)境庫茲涅茨假說的規(guī)律吻合，而我國目前仍處于EKC曲線的前期階段。在此基礎上對我國農(nóng)業(yè)經(jīng)濟水平與污染的現(xiàn)狀進行分析總結(jié)，提出相應建議。

關鍵詞農(nóng)業(yè)經(jīng)濟增長；農(nóng)業(yè)面源污染；面板數(shù)據(jù)

中圖分類號S-9文獻標識碼A文章編號0517-6611（2017）14-0216-02

AbstractBased on the panel data of three environmental pollution indicators in 31 provinces （municipalities and autonomous regions） in China from 1991 to 2014，the EKC test was carried out on the basis of environmental effects on the relationship between agricultural economic growth and agricultural environmental pollution.The results show that the relationship between agricultural economic growth and agricultural environmental pollution coincides with the traditional Kuznets hypothesis，and China is still in the early stage of EKC curve.On the basis of this，the present situation of agricultural economy and pollution in China was analyzed and summarized，and the corresponding suggestions were put forward.

Key wordsAgricultural economic growth；Agricultural nonpoint source pollution；Panel data

農(nóng)業(yè)經(jīng)濟得到一定發(fā)展，而農(nóng)業(yè)經(jīng)濟的發(fā)展與農(nóng)業(yè)生態(tài)環(huán)境的矛盾日漸嚴重。我國部分地區(qū)化肥、農(nóng)藥等不合理使用現(xiàn)象頻繁發(fā)生，導致水體污染、降低土壤肥力，成為水體富營養(yǎng)化的直接主要來源，嚴重破壞農(nóng)業(yè)生態(tài)環(huán)境。農(nóng)業(yè)環(huán)境問題對農(nóng)產(chǎn)品質(zhì)量安全問題具有關鍵性的影響，直接影響居民生產(chǎn)、生活水平。2016年中央一號文件提出加強資源保護與生態(tài)修復，推動綠色農(nóng)業(yè)發(fā)展。農(nóng)業(yè)生產(chǎn)對生態(tài)環(huán)境的負面效應引發(fā)了國內(nèi)外學者的思考。不少學者基于倒“U”型的環(huán)境庫茨涅茲曲線（EKC）分析驗證經(jīng)濟增長與環(huán)境質(zhì)量是否能夠同向變動以及兩者間的協(xié)整關系[1-3]。胡宗義等[4]研究了低碳經(jīng)濟背景下碳排放的環(huán)境庫茨涅茲曲線。李海鵬等[5]利用環(huán)境庫茲涅茨模型對我國經(jīng)濟增長與農(nóng)業(yè)面源污染進行實證分析。但ECK假定受到Arrow等批評，認為由于技術水平與農(nóng)業(yè)結(jié)構合理性的限制在低水平地區(qū)與發(fā)達地區(qū)的環(huán)境惡化狀況并不一致。因此，部分學者不再基于EKC曲線分析，而將視角轉(zhuǎn)向為Grossman等[6]提出的環(huán)境三效應分析框架：規(guī)模效應、結(jié)構效應、技術效應。葛繼紅等[7]、梁流濤等[8]從三效應出發(fā)分解污染演變的影響因素，構建其污染演變過程的內(nèi)在驅(qū)動機理模型以分析經(jīng)濟增長對農(nóng)業(yè)面源污染的影響。戴紅軍等[9]建立OLS多元線性模型，從三效應角度對農(nóng)村生態(tài)環(huán)境污染的影響因素進行動態(tài)分析。張智奎等[10]在結(jié)構效應視角下對經(jīng)濟發(fā)展與農(nóng)業(yè)面源污染關系進行協(xié)整檢驗，認為經(jīng)濟結(jié)構嚴重影響農(nóng)村生活污染源的排放。楊珂玲等[11]基于產(chǎn)業(yè)結(jié)構調(diào)整視角，強調(diào)從結(jié)構效應上制定農(nóng)業(yè)面源污染在種植業(yè)和畜牧業(yè)兩方面的控制政策。

鑒于我國農(nóng)業(yè)面源污染與經(jīng)濟增長之間關系的研究較少，筆者以1991—2014年我國31?。ㄊ?、自治區(qū)）面板數(shù)據(jù)為樣本，就兩者之間的關系構造EKC曲線，并借鑒Grossman等[6]環(huán)境三效應分析框架對所構造的曲線進行理論釋義。

1指標選取、數(shù)據(jù)來源與研究方法

1.1指標選取與數(shù)據(jù)來源筆者主要研究農(nóng)業(yè)經(jīng)濟增長對農(nóng)業(yè)環(huán)境污染的影響，衡量經(jīng)濟發(fā)展水平指標變量選取31省（市、自治區(qū)）農(nóng)林漁牧業(yè)生產(chǎn)總值。根據(jù)傳統(tǒng)EKC曲線理論，隨著農(nóng)林漁牧業(yè)生產(chǎn)總值不斷上升，初期階段農(nóng)業(yè)環(huán)境污染程度逐漸加深。而當農(nóng)業(yè)經(jīng)濟增長到某一程度后，污染程度開始下降，并隨經(jīng)濟增長繼續(xù)得到改善。根據(jù)理論假設，將農(nóng)、林、漁牧業(yè)生產(chǎn)總值作為實證模型的解釋變量，擬合二次線性函數(shù)驗證EKC曲線與我國農(nóng)業(yè)經(jīng)濟增長對農(nóng)業(yè)環(huán)境污染影響的適應情況。

衡量經(jīng)濟發(fā)展水平指標與農(nóng)業(yè)面源污染物指標數(shù)據(jù)主要來源于中國統(tǒng)計局網(wǎng)站、《中國農(nóng)村統(tǒng)計年鑒》，部分數(shù)據(jù)來自于中國農(nóng)業(yè)部、中國環(huán)境統(tǒng)計年報。為了避免消費物價指數(shù)CPI等宏觀經(jīng)濟因素隨市場價格波動的影響，保證模型設定的準確、可靠性，各地區(qū)各年份農(nóng)林漁牧業(yè)生產(chǎn)總值不按照名義產(chǎn)值，而按照價格定基指數(shù)以1991年為基年折算為GDP實際數(shù)據(jù)。由于重慶市1996年以前原始數(shù)據(jù)的部分缺失性，統(tǒng)計研究數(shù)據(jù)為非平衡的面板數(shù)據(jù)。

衡量農(nóng)業(yè)環(huán)境污染程度的污染性指標包括3項：農(nóng)業(yè)化肥（Agriferti）、氮肥施用量（N-ferti）、磷肥施用量（P-ferti）（折純量）。農(nóng)業(yè)化肥的大量施用造成土壤水溶性養(yǎng)分等隨雨水溶入地下水、河流、湖泊，水體富營養(yǎng)化；化肥中的氮、磷、等化學物質(zhì)在土壤中固結(jié)成鹽，土壤養(yǎng)分結(jié)構失衡，是農(nóng)業(yè)環(huán)境加劇的關鍵性因素。鑒于農(nóng)業(yè)施肥品種結(jié)構中氮肥、磷肥所占權重最大， P、N是導致水體污染，惡化土壤理化質(zhì)量的主要元素，氮肥、磷肥選作污染統(tǒng)計指標更具代表性。

1.2實證模型設立根據(jù)以上模型變量選取及農(nóng)業(yè)經(jīng)濟增長與農(nóng)業(yè)環(huán)境污染間關系的預期假設，建立實證回歸模型，得出各自變量指標對因變量的貢獻度。農(nóng)業(yè)經(jīng)濟增長對面源污染影響的模型如下：

式中，Pollutionit代表i地區(qū)第t年的污染指標；Agrigdpit表示i地區(qū)第t年的國內(nèi)生產(chǎn)總值；α、β1、β2分別表示模型截距項與解釋變量系數(shù)；ζit為隨機干擾項，表示無法觀測到的各項其他影響因素。在不考慮其他各項非關鍵性變量影響的條件下，對農(nóng)業(yè)經(jīng)濟增長對環(huán)境污染的影響進行重點分析。若有β1>0且β2<0，則有農(nóng)業(yè)經(jīng)濟增長與環(huán)境污染之間滿足預期假設的倒“U”型二次線性關系。經(jīng)過計算可得EKC曲線拐點處的地區(qū)GDP值滿足：

（2）

2結(jié)果與分析

2.1實證結(jié)果基于1991—2014年省際面板數(shù)據(jù)分別對Agriferti、N-ferti、P-ferti這3項污染指標進行變截距模型估計與檢驗，結(jié)果見表1。

表1給出了三項污染指標為基礎的模型估計結(jié)果，不難發(fā)現(xiàn)化肥施用量為污染指標的固定效應模型、隨機效應模型一次變量系數(shù)為均正數(shù)，二次項系數(shù)均為負數(shù)，且在0.01的水平上統(tǒng)計顯著?？傮w上看，可決系數(shù)R2分別為0.776 6、0.791 2，F(xiàn)統(tǒng)計量均在0.01的統(tǒng)計水平下顯著，可判斷以化肥施用總量為污染指標實證得到的模型整體擬合效果較好，證明了農(nóng)業(yè)經(jīng)濟對農(nóng)業(yè)環(huán)境污染的影響表現(xiàn)為倒“U”型曲線。除了化肥施用模型之外，其他2個污染指標構建的模型變量系數(shù)同樣滿足β1>0且β2<0，各項指標很好的通過檢驗，進一步證實農(nóng)業(yè)經(jīng)濟與環(huán)境污染之間的關系。綜上，關于農(nóng)業(yè)經(jīng)濟增長對農(nóng)業(yè)環(huán)境污染影響估計得到的模型與預期假設完全一致，表明我國農(nóng)業(yè)經(jīng)濟增長對農(nóng)業(yè)環(huán)境污染表現(xiàn)出的影響與Grossman等[6]提出的環(huán)境庫茲涅茨理論正好吻合。

表1還顯示，在3個不同指標設定的模型下計算出的農(nóng)業(yè)經(jīng)濟增長“拐點”，各地區(qū)農(nóng)業(yè)生產(chǎn)總值大概在到達9 000～13 000億元區(qū)間之后，農(nóng)業(yè)環(huán)境污染便開始呈明顯下降趨勢。目前來看，僅有山東、上海、江西等少數(shù)?。ㄊ校┺r(nóng)業(yè)經(jīng)濟水平將達到或剛達到EKC“拐點”，新疆、海南、四川、云南、寧夏等地區(qū)離經(jīng)濟“拐點”差距較大，整體上各地區(qū)均在努力提高農(nóng)業(yè)發(fā)展水平以達到“拐點”處農(nóng)業(yè)經(jīng)濟水平。

2.2EKC曲線機理分析實證檢驗結(jié)果顯示，我國農(nóng)業(yè)面源與農(nóng)業(yè)經(jīng)濟增長之間呈現(xiàn)經(jīng)典的倒“U”型曲線，這可以用Grossman等[6]的環(huán)境三效應理論進行釋義：家庭聯(lián)產(chǎn)承包責任制在制度上解決了對農(nóng)業(yè)生產(chǎn)力的束縛，推動農(nóng)業(yè)生產(chǎn)規(guī)模擴大；對外開放為我國農(nóng)業(yè)經(jīng)濟增長注入了新鮮的生產(chǎn)要素——化肥和農(nóng)藥的施用，生產(chǎn)規(guī)模效應開始顯現(xiàn)，加大了對自然資源開采，破壞生態(tài)環(huán)境，尤其是化肥、農(nóng)藥的過量施用使得前期農(nóng)業(yè)環(huán)境污染程度逐漸嚴重（規(guī)模效應）。農(nóng)業(yè)經(jīng)濟增長使得人民收入水平提升，以至于居民偏好優(yōu)質(zhì)生態(tài)環(huán)境，開始顯現(xiàn)對污染治理的努力。普及農(nóng)業(yè)環(huán)境友好型生產(chǎn)技術，旨在減少面源污染的項目，發(fā)展測土配方施肥技術（技術效應）。隨著技術服務的普及，農(nóng)業(yè)生產(chǎn)性服務業(yè)開始出現(xiàn)，這些新興行業(yè)主要為農(nóng)業(yè)綠色轉(zhuǎn)型提供技術支持，它們的存在有力地促進了農(nóng)業(yè)產(chǎn)業(yè)結(jié)構的進一步優(yōu)化（結(jié)構效應），從而農(nóng)業(yè)環(huán)境狀況開始隨著農(nóng)業(yè)經(jīng)濟的發(fā)展得到不斷的改善。

3結(jié)論與建議

該研究以1991—2014年我國31省（市、自治區(qū)）的面板數(shù)據(jù)為研究對象，對農(nóng)業(yè)環(huán)境污染演變的內(nèi)在機理進行實證分析。結(jié)果表明，農(nóng)業(yè)經(jīng)濟增長是影響農(nóng)業(yè)環(huán)境污染的重要原因，我國農(nóng)業(yè)經(jīng)濟增長與環(huán)境污染之間的關系符合環(huán)境EKC曲線的走勢，一定時期內(nèi)農(nóng)業(yè)經(jīng)濟增長與農(nóng)業(yè)環(huán)境污染表現(xiàn)為同方向的變動，隨著農(nóng)業(yè)經(jīng)濟增長，農(nóng)業(yè)環(huán)境污染在經(jīng)濟達到某一程度后得到改善，并將呈現(xiàn)出污染程度逐步下降的趨勢。

針對我國目前農(nóng)業(yè)經(jīng)濟與環(huán)境的總體情況，提出下列建議：改善產(chǎn)業(yè)結(jié)構，實現(xiàn)農(nóng)業(yè)經(jīng)濟低污染、少消耗、高產(chǎn)出。主要從農(nóng)業(yè)結(jié)構出發(fā)，調(diào)整農(nóng)、林、漁、牧業(yè)之間的比例，把握種植業(yè)與養(yǎng)殖業(yè)的均衡協(xié)調(diào)發(fā)展，對于對環(huán)境質(zhì)量影響較大的產(chǎn)業(yè)進行產(chǎn)業(yè)項目升級或轉(zhuǎn)移。將科學技術用于農(nóng)業(yè)生產(chǎn)，鼓勵科技創(chuàng)新。改變傳統(tǒng)農(nóng)業(yè)的經(jīng)營管理模式，運用生物技術、信息技術、新材料技術改造農(nóng)業(yè)，通過高新技術研制可降解度較高的農(nóng)藥、化肥、農(nóng)用薄膜等，在推動農(nóng)業(yè)經(jīng)濟發(fā)展的同時達到環(huán)境影響最小化。

參考文獻

[1] 潘培，楊順順，欒勝基.我國農(nóng)村居民消費結(jié)構變化及其環(huán)境影響分析[J].安徽農(nóng)業(yè)科學，2009，37（26）：12732-12735，12772.

[2] 胡浩，張暉，黃士新.規(guī)模養(yǎng)殖戶健康養(yǎng)殖行為研究：以上海市為例[J].農(nóng)業(yè)經(jīng)濟問題，2009（8）：25-32.

[3] 覃巍，丁慧.基于環(huán)境庫茲涅茨理論的廣西經(jīng)濟增長與環(huán)境質(zhì)量協(xié)調(diào)性研究[J].學術論壇，2011（10）：142-148.

[4] 胡宗義，劉亦文，唐李偉.低碳經(jīng)濟背景下碳排放的庫茲涅茨曲線研究[J].統(tǒng)計研究，2013，30（2）：73-79.

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

[6] GROSSMAN G M，KRUEGER A B.Economic growth and the environment[J].Quarterly joural of economics，1995，110（2）：353-377.

[7] 葛繼紅，周曙東.農(nóng)業(yè)面源污染的經(jīng)濟影響因素分析：基于1978～2009年的江蘇省數(shù)據(jù)[J].中國農(nóng)村經(jīng)濟，2011（5）：72-81.

[8] 梁流濤，曲福田，馮淑怡.經(jīng)濟發(fā)展與農(nóng)業(yè)面源污染：分解模型與實證研究[J].長江流域資源與環(huán)境，2013，22（10）：1369-1374.

[9] 戴紅軍，孫濤.城鄉(xiāng)一體化進程中農(nóng)村生態(tài)環(huán)境污染影響因素分析：基于江蘇省1990-2009年的數(shù)據(jù)[J].生態(tài)經(jīng)濟，2013（11）：177-181.

[10] 張智奎，肖新成.經(jīng)濟發(fā)展與農(nóng)業(yè)面源污染關系的協(xié)整檢驗：基于三峽庫區(qū)重慶段1992-2009年數(shù)據(jù)的分析[J].中國人口·資源與環(huán)境，2012，22（1）：57-61.

[11] 楊珂玲，張宏志.基于產(chǎn)業(yè)結(jié)構調(diào)整視角的農(nóng)業(yè)面源污染控制政策研究[J].生態(tài)經(jīng)濟，2015，31（3）：89-92.

[12] 包群，彭水軍.經(jīng)濟增長與環(huán)境污染：基于面板數(shù)據(jù)的聯(lián)立方程估計[J].世界經(jīng)濟，2006（11）：48-58.

[13] 鐘茂初，張學剛.環(huán)境庫茲涅茨曲線理論及研究的批評綜述[J].中國人口·資源與環(huán)境，2010，20（2）：62-67.