農(nóng)業(yè)面源污染EKC再檢驗

馬子龍 閔晨晨

摘? 要：污染防治攻堅戰(zhàn)的提出督促著學術(shù)界對促進經(jīng)濟與自然的和諧統(tǒng)一展開進一步的研究?；?003-2017年中國省級面板數(shù)據(jù)，采用動態(tài)面板門檻模型對農(nóng)業(yè)不同發(fā)展水平下的化肥污染EKC進行檢驗后發(fā)現(xiàn)：化肥污染與人均農(nóng)業(yè)實際產(chǎn)值間存在以農(nóng)業(yè)發(fā)展水平為門檻變量的EKC關系。當農(nóng)業(yè)實際產(chǎn)值低于611.55億元時，兩者呈“N”型關系;當農(nóng)業(yè)實際產(chǎn)值高于611.55億元時，兩者呈倒“U”型關系。

關鍵詞：農(nóng)業(yè)污染;EKC檢驗;動態(tài)面板門檻回歸

一、引言

黨的第十九屆全國人民代表大會提出：打贏防范化解重大風險、精準脫貧和污染防治的“三大攻堅戰(zhàn)”，是新時代決勝建成全面小康社會的關鍵。自2020年底中國實現(xiàn)全面脫貧后，污染防治問題的重要性越發(fā)凸顯。在此背景下，加強生態(tài)文明建設，踐行“綠水青山就是金山銀山”理念成為農(nóng)業(yè)污染治理工作的重要理論指導。

在現(xiàn)有文獻中，環(huán)境庫茲涅茨曲線（EKC）是用來描述環(huán)境和經(jīng)濟關系最常見的理論。大量的研究驗證了兩者之間的密切相關性。那么，EKC假說是否適用于農(nóng)業(yè)污染領域？農(nóng)業(yè)污染與農(nóng)業(yè)增長的關系曲線為何種形狀？厘清上述問題不僅有利于探究農(nóng)業(yè)污染性質(zhì)，科學把握污染與經(jīng)濟的內(nèi)在聯(lián)系;同時對實現(xiàn)農(nóng)業(yè)面源污染治理和促進農(nóng)業(yè)現(xiàn)代化也有著重要的現(xiàn)實意義。

二、研究設計與實證分析

（一）模型構(gòu)建

根據(jù)EKC假說，經(jīng)濟增長和環(huán)境污染之間一般存在著倒“U”型關系。本文采取雙對數(shù)形式來構(gòu)建基本模型：

其中， 表示農(nóng)業(yè)化肥污染; 表示人均農(nóng)業(yè)GDP， 與 為其高次項;個體固定效應μi代表著個省份不隨時間變化的特征，年份固定效應γt則表示所有省份共有的時間因素。εi，t為隨機擾動項，i表示個體，t表示時間。

進一步的，由于農(nóng)業(yè)污染水平一般存在著顯著的“路徑依賴”，即當前的農(nóng)業(yè)污染水平往往會受到往期農(nóng)業(yè)污染的影響。為了考察農(nóng)業(yè)污染的長期效應，本文在模型（1）中加入農(nóng)業(yè)污染的滯后項：

另一方面，一些學者的結(jié)論表明農(nóng)業(yè)污染與農(nóng)業(yè)增長的關系可能會受到其他變量的門檻作用，為了檢驗不同經(jīng)濟水平下，中國農(nóng)業(yè)污染與農(nóng)業(yè)增長的EKC差異，本文在模型（3）的基礎上加入農(nóng)業(yè)實際GDP作為門檻變量，構(gòu)建動態(tài)面板門檻回歸模型：

模型（3）中Q表示實際農(nóng)業(yè)產(chǎn)值的門檻值，通過觀測實際農(nóng)業(yè)產(chǎn)值在不同區(qū)間時農(nóng)業(yè)污染及其高次項系數(shù)，便可檢驗不同農(nóng)業(yè)發(fā)展水平下污染與增長的EKC形狀。

（二）數(shù)據(jù)及變量說明

本文以2003-2017年中國30個省市面板數(shù)據(jù)為樣本，對當前農(nóng)業(yè)面源污染的EKC形狀進行檢驗。數(shù)據(jù)來源于《中國統(tǒng)計年鑒》、《中國農(nóng)村統(tǒng)計年鑒》以及國家統(tǒng)計局數(shù)據(jù)網(wǎng)站。部分缺失數(shù)據(jù)用插值法進行插補。各變量所表示的具體內(nèi)容如下所示：

1.本文解釋變量為農(nóng)業(yè)人均GDP，為了剔除價格變動帶來的影響，根據(jù)已有研究，采取1990年不變價格計算出實際農(nóng)業(yè)產(chǎn)值后，再結(jié)合農(nóng)業(yè)勞動力人數(shù)計算得出。被解釋變量為農(nóng)業(yè)污染，選取中國最主要的農(nóng)業(yè)污染源化肥污染作為研究對象，在現(xiàn)有文獻的基礎上，采用化肥施用密度作為衡量。

2.一此外，些學者認為：EKC形狀在不同的收入水平存在著差異，因此，本文選取了農(nóng)業(yè)發(fā)展水平作為門檻變量以考察在異質(zhì)性條件下農(nóng)業(yè)污染與增長的EKC關系。

3.城鎮(zhèn)化率和產(chǎn)業(yè)結(jié)構(gòu)被認為是生態(tài)環(huán)境的重要影響因素，故本文將其作為控制變量納入模型。

（三）結(jié)果展示

表1的第一列結(jié)果顯示了農(nóng)業(yè)產(chǎn)值處于低閾值區(qū)間時的回歸結(jié)果。其中：農(nóng)業(yè)人均產(chǎn)值三次項的系數(shù)為正，在5%的水平上顯著，表明當農(nóng)業(yè)實際GDP小于611.55億元時，農(nóng)業(yè)增長與化肥污染呈“N”型關系，即農(nóng)業(yè)人均產(chǎn)值的增加最終會加劇農(nóng)業(yè)污染的嚴重程度。第二列結(jié)果則給出了農(nóng)業(yè)產(chǎn)值處于高閾值區(qū)間時的回歸結(jié)果。結(jié)果顯示：農(nóng)業(yè)人均產(chǎn)值二次項的系數(shù)在1%的水平上顯著為負，三次項系數(shù)不顯著，表明當農(nóng)業(yè)實際產(chǎn)值大于611.55億元時，農(nóng)業(yè)增長與化肥污染呈倒“U”型關系。

三、結(jié)論與建議

科學把握環(huán)境與經(jīng)濟的關系，對我國推動農(nóng)業(yè)綠色轉(zhuǎn)型，打造現(xiàn)代化農(nóng)業(yè)，實現(xiàn)科學發(fā)展有著重要的理論和現(xiàn)實意義。鑒于當前研究的主要缺陷，本文基于2003-2017年中國省級面板數(shù)據(jù)，采用動態(tài)面板門檻回歸模型對中國農(nóng)業(yè)化肥污染的EKC關系及形成機制進行了檢驗，研究主要結(jié)論如下：化肥污染與人均農(nóng)業(yè)實際產(chǎn)值間存在以農(nóng)業(yè)發(fā)展水平為門檻變量的EKC關系。當農(nóng)業(yè)實際產(chǎn)值低于611.55億元時，兩者呈“N”型;當農(nóng)業(yè)實際產(chǎn)值高于611.55億元時，兩者呈倒“U”型。

根據(jù)研究結(jié)論，政策建議如下：

大力發(fā)展農(nóng)業(yè)現(xiàn)代化，促進農(nóng)業(yè)生產(chǎn)專業(yè)化和高度化。應增大資本投入，增強農(nóng)業(yè)產(chǎn)出能力，加快農(nóng)業(yè)產(chǎn)值由低到高過渡;提高居民環(huán)保意識，促進清潔生產(chǎn)技術(shù)的改進和創(chuàng)新;積極調(diào)整生產(chǎn)結(jié)構(gòu)，引領種植結(jié)構(gòu)和要素投入合理化轉(zhuǎn)變。

參考文獻

[1]揭昌亮，王金龍，龐一楠.中國農(nóng)業(yè)增長與化肥面源污染：環(huán)境庫茲涅茨曲線存在嗎？[J].農(nóng)村經(jīng)濟，2018（11）：110-117.

[2]Jiang S， Qiu S， Zhou J. Re-Examination of the Relationship between Agricultural Economic Growth and Non-Point Source Pollution in China： Evidence from the Threshold Model of Financial Development[J]. Water， 2020， 12（9）： 2609.

[3]侯孟陽，姚順波.異質(zhì)性條件下化肥面源污染排放的EKC再檢驗——基于面板門檻模型的分組[J].農(nóng)業(yè)技術(shù)經(jīng)濟，2019（04）：104-118.

[4] Song J， Guo Y， Wu P， et al. The agricultural water rebound effect in China[J]. Ecological Economics， 2018， 146： 497-506.