【摘要】本文基于1994～2014年間的省級面板數(shù)據(jù)，通過構(gòu)建土壤環(huán)境與農(nóng)業(yè)經(jīng)濟增長的VAR模型，在區(qū)域?qū)用鏅z驗與分析了我國土壤環(huán)境與農(nóng)業(yè)經(jīng)濟增長之間的相互關(guān)系。研究結(jié)果表明：農(nóng)業(yè)經(jīng)濟增長對土壤環(huán)境的影響以及土壤環(huán)境對農(nóng)業(yè)經(jīng)濟增長的影響都存在明顯的區(qū)域差異，對于前者而言，西北、東北地區(qū)所受影響最大，西南地區(qū)所受影響最小，驗證了西北、東北地區(qū)的農(nóng)業(yè)經(jīng)濟增長對土壤環(huán)境的破壞力度相對較高，西南地區(qū)破壞力度相對較小。對于后者而言，西南地區(qū)所受影響最大，華中、華南地區(qū)所受影響較小，這說明西南地區(qū)目前土壤環(huán)境比較優(yōu)良，對農(nóng)業(yè)經(jīng)濟可持續(xù)發(fā)展能帶來更大的增長空間，而華中、華南地區(qū)土壤環(huán)境已經(jīng)欠佳，能為地區(qū)農(nóng)業(yè)經(jīng)濟增長提供的動力較小。

【關(guān)鍵詞】土壤環(huán)境 農(nóng)業(yè)經(jīng)濟增長 VAR模型

一、引言

隨著農(nóng)村環(huán)境問題越來越嚴峻，關(guān)于農(nóng)業(yè)生產(chǎn)環(huán)境與農(nóng)業(yè)經(jīng)濟增長之間的關(guān)系，學術(shù)界已經(jīng)展開了較為深入的研究，目前主要有以下幾類，一類是沿用Grossman和Krueger（1993）等研究經(jīng)濟增長與工業(yè)污染之間關(guān)系所運用的環(huán)境庫茲涅茲曲線（EKC）驗證方法，來探究農(nóng)業(yè)面源污染與農(nóng)業(yè)經(jīng)濟增長之間的關(guān)系，如張暉、胡浩（2009）通過建立江蘇省經(jīng)濟發(fā)展水平與過剩氮污染的經(jīng)濟計量模型，對二者的關(guān)系進行了環(huán)境庫茲涅茲曲線驗證，研究結(jié)果表明，人均過剩氮排放量和經(jīng)濟增長之間存在顯著的“倒U型”曲線關(guān)系。[1]還有一類是借鑒Islam的研究思路，通過構(gòu)建農(nóng)業(yè)面源污染影響因素分解的理論模型，進行實證研究。如梁流濤、曲福天等（2013）利用計量經(jīng)濟模型在農(nóng)業(yè)面源污染影響因素分解理論模型的基礎(chǔ)上探尋了農(nóng)業(yè)面源污染演變的內(nèi)在驅(qū)動機制，得出經(jīng)濟作物比重提高和畜禽養(yǎng)殖業(yè)產(chǎn)值比重的提高都會導致農(nóng)業(yè)面源污染加劇等結(jié)論。[2]第三類是以沈霞林（2011）為代表，選取多個因素作為農(nóng)業(yè)經(jīng)濟發(fā)展指標，采用多元回歸分析法，測算區(qū)域農(nóng)業(yè)面源污染水平與區(qū)域農(nóng)業(yè)經(jīng)濟發(fā)展水平之間的關(guān)系，以給出區(qū)域農(nóng)業(yè)面源污染差異的經(jīng)濟解釋。[3]本文正是在包含以上這些文獻的基礎(chǔ)上，展開了對農(nóng)業(yè)生產(chǎn)環(huán)境與農(nóng)業(yè)經(jīng)濟增長之間關(guān)系的實證研究。

二、研究方法與數(shù)據(jù)說明

研究方法上，本文運用的是VAR模型，數(shù)據(jù)方面，以1994～2014年這20年間的中國省級面板數(shù)據(jù)為依托，以農(nóng)業(yè)一項的產(chǎn)值作為農(nóng)業(yè)經(jīng)濟增長的指標，以化肥、農(nóng)藥及地膜的投入密度為基礎(chǔ)計算出的土壤環(huán)境污染綜合值[4]作為土壤污染測算指標，基礎(chǔ)數(shù)據(jù)來源于《中國統(tǒng)計年鑒》與《中國農(nóng)村統(tǒng)計年鑒》。

另外，考慮到我國國土面積廣大，不同區(qū)域經(jīng)濟發(fā)展狀況與農(nóng)業(yè)發(fā)展水平都有所不同，所以本文對所研究的區(qū)域進行了劃分。依據(jù)經(jīng)濟社會發(fā)展程度以及地理方位的不同，將研究的31個?。ㄊ小⒆灾螀^(qū)）劃分為東北、西北、華北、華中、華東、華南、西南一共七個區(qū)域進行對比研究[5]。在基礎(chǔ)數(shù)據(jù)處理上，為避免物價波動的干擾，本文將每年農(nóng)業(yè)總產(chǎn)值都折算為1994年可比價，同時，為消除數(shù)據(jù)中存在的異方差，將代表各地區(qū)農(nóng)業(yè)經(jīng)濟增長的農(nóng)業(yè)總產(chǎn)值以及代表土壤環(huán)境的土壤環(huán)境污染綜合值再分別進行對數(shù)化處理，其對應形式分別為lngdp、lnsoil，而下文中的Dlngdp、Dlnsoil與D（lngdp，2）、D（lnsoil，2）分別表示一階和二階差分序列。

三、實證分析

（一）單位根檢驗

對時間序列進行分析的前提是保證序列的平穩(wěn)性，因此一般情況下進行分析之前，需要對原變量序列進行單位根檢驗，判斷序列的平穩(wěn)性。[6]本文選取較為常用的ADF單位根檢驗方法對處理過的數(shù)據(jù)進行檢驗。通過檢驗結(jié)果可以看出七個地區(qū)中，除了華東地區(qū)的lngdp與lnsoil是在二階差分形式下同時平穩(wěn)外，其他地區(qū)的lngdp與lnsoil都是在一階差分形式下同時平穩(wěn)，所以，華東地區(qū)是二階單整序列，其他地區(qū)都是一階單整序列，滿足協(xié)整檢驗的前提。

（二）協(xié)整檢驗

協(xié)整檢驗是為了進一步判斷土壤環(huán)境與農(nóng)業(yè)經(jīng)濟增長之間是否存在長期的均衡關(guān)系，為此本文采取常用的協(xié)整檢驗方法Johansen協(xié)整檢驗，做進一步探究。協(xié)整檢驗結(jié)果表明，在最佳滯后期，華東、西南、華北、東北地區(qū)的lngdp與lnsoil之間有協(xié)整關(guān)系，而華中、華南以及西北地區(qū)lngdp與lnsoil之間沒有協(xié)整關(guān)系。這說明，前者的土壤環(huán)境與農(nóng)業(yè)經(jīng)濟增長具有長期協(xié)整關(guān)系，而后者兩者之間的長期協(xié)整關(guān)系在最佳滯后期并不顯著。不顯著的原因可能是因為這三個地區(qū)的土壤環(huán)境或農(nóng)業(yè)經(jīng)濟增長已經(jīng)基本達到了飽和，因此造成土壤環(huán)境的變化對農(nóng)業(yè)經(jīng)濟增長或農(nóng)業(yè)經(jīng)濟增長對土壤環(huán)境的變化沒有顯著的反應。

（三）Granger因果關(guān)系檢驗

協(xié)整檢驗只反映了變量之間存在長期因果關(guān)系，并不能明確兩者之間因果關(guān)系的具體方向。為了進一步明確二者之間的關(guān)系，需要在此基礎(chǔ)上進行Granger因果關(guān)系檢驗。因為之前的協(xié)整檢驗結(jié)果表明，華北、東北、西南地區(qū)的土壤環(huán)境與農(nóng)業(yè)經(jīng)濟增長之間都具有長期的協(xié)整關(guān)系，所以在此使用原序列進行長短期的Granger因果檢驗，而華中、華南與西北地區(qū)沒有長期協(xié)整關(guān)系，所以在此用一階差分形式來做Granger因果關(guān)系檢驗，此外，華東地區(qū)因為是二階單整，所以也用一階差分形式來做Granger因果關(guān)系檢驗。檢驗結(jié)果顯示，在最佳滯后期，華北、東北地區(qū)存在lnsoil到lngdp的單項因果關(guān)系，西北也存在Dlnsoil到Dlngdp的單項因果關(guān)系，而西南地區(qū)存在lngdp到lnsoil的單項因果關(guān)系，華東也存在Dlngdp到Dlnsoil的單項因果關(guān)系，但華中、華南地區(qū)的Dlngdp和Dlnsoil并不存在Granger因果關(guān)系。這說明，華北、東北、西北地區(qū)前期土壤環(huán)境的變化對隨后的農(nóng)業(yè)經(jīng)濟增長的程度有較大影響，而西南地區(qū)和華東地區(qū)則是前期的農(nóng)業(yè)經(jīng)濟增長的程度對隨后的土壤環(huán)境變化有較大影響，華中、華南地區(qū)土壤環(huán)境與農(nóng)業(yè)經(jīng)濟增長之間并不存在顯著的Granger因果關(guān)系。

（四）方差分解

為了更清楚地刻畫和度量土壤環(huán)境與農(nóng)業(yè)經(jīng)濟增長的相互影響程度，本文進一步采用方差分解的方法，獲得不同方程的沖擊反應對各個變量波動的方差貢獻率構(gòu)成。下表1-1、表1-2給出了各個地區(qū)預測期的方差分解結(jié)果。

表1-1、表1-2的方差結(jié)果表明，（1）各地區(qū)經(jīng)過5個預測期以后，系統(tǒng)就已經(jīng)基本穩(wěn)定。（2）土壤環(huán)境與農(nóng)業(yè)經(jīng)濟增長的波動均主要來自于自身，農(nóng)業(yè)經(jīng)濟增長對自身波動的貢獻比率都在80%以上，土壤環(huán)境對自身波動的貢獻比率除了西南地區(qū)外也都在60%以上。（3）農(nóng)業(yè)經(jīng)濟增長對土壤環(huán)境的影響在2%-15%之間，其中西北、東北地區(qū)所受影響最大占到了15%左右，其次為華北、華南、華東、華中，影響最小的是西南地區(qū)，僅為2%。（4）土壤環(huán)境對農(nóng)業(yè)經(jīng)濟增長的影響方面，西南地區(qū)影響最大，達到了60%，其次是華東、東北、華北、西北地區(qū)，影響程度在30%左右，而華中、華南地區(qū)土壤環(huán)境對農(nóng)業(yè)經(jīng)濟增長的影響較小，分別為13%何9%左右。

四、結(jié)論

華北、東北、華東以及西南地區(qū)的土壤環(huán)境與農(nóng)業(yè)經(jīng)濟增長之間在各自的最佳滯后期內(nèi)具有長期協(xié)整關(guān)系，而華中、華南以及西北三個地區(qū)長期協(xié)整關(guān)系在最佳滯后期內(nèi)并不顯著。

華北、東北、西北地區(qū)土壤環(huán)境前期的變化對隨后的農(nóng)業(yè)經(jīng)濟增長程度有較大影響，而西南地區(qū)和華東地區(qū)則是前期的農(nóng)業(yè)經(jīng)濟增長的程度對隨后的土壤環(huán)境變化有較大影響，華中、華南地區(qū)土壤環(huán)境與農(nóng)業(yè)經(jīng)濟增長之間并不存在明顯的Granger因果關(guān)系。

方差分解的檢驗結(jié)果更進一步表明，基本上經(jīng)過5個預測期后，各地區(qū)的系統(tǒng)基本趨于穩(wěn)定，土壤環(huán)境與農(nóng)業(yè)經(jīng)濟增長的波動均主要來自于自身，農(nóng)業(yè)經(jīng)濟增長對土壤環(huán)境的影響僅在2%-15%之間，其中西北、東北地區(qū)所受影響最大，其次為華北、華南、華東、華中，影響最小的是西南地區(qū)。這說明西北、東北地區(qū)的農(nóng)業(yè)經(jīng)濟增長對土壤環(huán)境的破壞力度相對較高，西南地區(qū)破壞力度相對較小。土壤環(huán)境對農(nóng)業(yè)經(jīng)濟增長的影響方面，西南地區(qū)最大，其次是華東、東北、華北、西北地區(qū)，而華中、華南地區(qū)土壤環(huán)境對農(nóng)業(yè)經(jīng)濟增長的影響較小。這也側(cè)面反映了西南地區(qū)目前土壤環(huán)境比較優(yōu)良，對農(nóng)業(yè)經(jīng)濟可持續(xù)發(fā)展而言能帶來更大的增長空間，而華中、華南地區(qū)土壤環(huán)境已經(jīng)欠佳，能為地區(qū)農(nóng)業(yè)經(jīng)濟增長提供的動力較小。

注釋

{1}土壤環(huán)境污染綜合值是以94年為基期，算出1994-2014年各地區(qū)化肥、農(nóng)藥以及地膜投入密度與當年全國平均值的比例，將這三項比例乘以相同權(quán)重再相加在一起得出。

{2}東北地區(qū)包括：黑龍江、吉林、遼寧3個省（市、自治區(qū)）；西北地區(qū)包括：寧夏、青海、新疆、陜西、甘肅5個?。ㄊ小⒆灾螀^(qū)）；華北地區(qū)包括：北京、天津、河北、山西、內(nèi)蒙古5個省（市、自治區(qū)）；華中地區(qū)包括：湖北、湖南、河南3個?。ㄊ?、自治區(qū)）；華東地區(qū)包括：山東、江蘇、安徽、江西、浙江、福建、上海7個?。ㄊ?、自治區(qū)）；華南地區(qū)包括：廣東、廣西、海南3個?。ㄊ小⒆灾螀^(qū)）；西南地區(qū)包括：四川、重慶、貴州、云南、西藏5個?。ㄊ?、自治區(qū)）。

參考文獻

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

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

[3]沈霞林.中國區(qū)域農(nóng)業(yè)面源污染與經(jīng)濟發(fā)展關(guān)系的實證分析[D].南京林業(yè)大學，2011.

[4]李嫣怡、劉榮、丁維岱等.EViews統(tǒng)計分析與應用[M].電子工業(yè)出版社.2013：312-319.

作者簡介：馬鑫然（1991-），女，漢，河南新鄉(xiāng)人，就讀于貴州大學經(jīng)濟學院，研究方向：國民經(jīng)濟。