農(nóng)業(yè)非點(diǎn)源污染與經(jīng)濟(jì)增長關(guān)系的實(shí)證研究
——基于環(huán)境庫茲涅茨曲線檢驗(yàn)

洪業(yè)應(yīng)

(涪陵行政學(xué)院科研管理與區(qū)情研究室，重慶涪陵 408000)

1 研究背景

非點(diǎn)源污染即面源污染，是指溶解性或固體污染物在大面積降水和徑流沖刷作用下匯入地表水體而引起的污染。農(nóng)業(yè)非點(diǎn)源污染則是指在農(nóng)業(yè)生產(chǎn)活動中，氮、磷等營養(yǎng)物質(zhì)，農(nóng)藥、重金屬等有機(jī)和無機(jī)污染物，土壤顆粒等沉積物，通過地表徑流和地下滲透，造成環(huán)境尤其是水域環(huán)境的污染［1］。國內(nèi)學(xué)者對農(nóng)業(yè)非點(diǎn)源污染和經(jīng)濟(jì)增長的關(guān)系進(jìn)行了分析、檢驗(yàn)，其主要集中在環(huán)境庫茲涅茨曲線(Environmental Kuznets Curve，EKC)是否具有普遍性及轉(zhuǎn)折點(diǎn)的確定，并取得了豐富的研究成果［2－3］，目前對農(nóng)業(yè)非點(diǎn)源污染的EKC假說的實(shí)證研究多數(shù)采用面板數(shù)據(jù)分析，該方法的假設(shè)是所有國家和區(qū)域都遵循相同的EKC發(fā)展路徑［4］;但是根據(jù)相關(guān)研究得出，時間序列數(shù)據(jù)能夠較好地擬合EKC模型，而面板數(shù)據(jù)卻不能夠得出理想的模擬效果［5］;在選取經(jīng)濟(jì)發(fā)展指標(biāo)時，諸多學(xué)者主要選用人均GDP、農(nóng)民人均純收入等指標(biāo)來表示［6－8］，然而化肥、農(nóng)膜和農(nóng)藥作為重要的農(nóng)業(yè)生產(chǎn)資料，主要用于農(nóng)業(yè)活動，尤其是農(nóng)業(yè)種植業(yè)，因此選取人均GDP、農(nóng)民人均純收入未免存在計算結(jié)果的誤差;同時國內(nèi)研究成果更多的是從宏觀層面和定性角度加以研究，研究的區(qū)域主要集中于東、中部區(qū)域，對于西部尤其是山區(qū)、巖溶地貌區(qū)域的農(nóng)業(yè)非點(diǎn)源污染是否符合“EKC”假說涉及較少，對于巖溶山區(qū)農(nóng)業(yè)非點(diǎn)源污染與經(jīng)濟(jì)增長是否具有類似“倒U型”關(guān)系，是否與前述研究一致，還有待于進(jìn)一步研究。

重慶市地處我國西南巖溶山區(qū)，其巖溶面積3.28 ×105km2，占全市國土面積的39.71%;山地、丘陵地貌占全市國土面積的94.00%;2011年年底全市常住人口為2 919萬人，人口密度為354人/km2，是全國平均水平(全國143人/km2)的2倍多。重慶人多地少、山高坡陡，“巴掌田”、“雞窩地”等居多，農(nóng)業(yè)生產(chǎn)條件差，生態(tài)環(huán)境脆弱等問題仍然十分突出;因地處長江上游三峽庫區(qū)，其生態(tài)地位十分重要。在重慶直轄以來，農(nóng)業(yè)取得了長足發(fā)展，但也由此帶來了農(nóng)業(yè)資源與環(huán)境問題的凸顯，其中農(nóng)業(yè)非點(diǎn)源污染在相當(dāng)時期內(nèi)已成為制約農(nóng)業(yè)循環(huán)經(jīng)濟(jì)發(fā)展的重要因素之一［9］。發(fā)達(dá)國家的經(jīng)驗(yàn)表明，農(nóng)村生產(chǎn)與生活引起的非點(diǎn)源污染已經(jīng)開始取代點(diǎn)源污染，逐漸成為了水體污染的重要因素。為此，本文基于重慶市1996—2011年時間序列數(shù)據(jù)，采用Eviews 6.0軟件，運(yùn)用計量分析方法研究農(nóng)用化肥、農(nóng)膜和農(nóng)藥單位播種面積施用量與農(nóng)業(yè)經(jīng)濟(jì)增長水平變量關(guān)系，探討其關(guān)系是否符合EKC模型，對山區(qū)農(nóng)業(yè)經(jīng)濟(jì)增長與農(nóng)業(yè)非點(diǎn)源污染之間的演變關(guān)系進(jìn)行系統(tǒng)分析。該研究結(jié)論也可為相關(guān)部門制定農(nóng)業(yè)環(huán)境政策和農(nóng)業(yè)的可持續(xù)發(fā)展提供一定的參考，提高人們對區(qū)域農(nóng)業(yè)非點(diǎn)源污染的重視。

2 變量選取與數(shù)據(jù)來源

2.1 變量選取

根據(jù)重慶市農(nóng)業(yè)面源污染的相關(guān)研究［9］，本文選取投入農(nóng)業(yè)生產(chǎn)的主要農(nóng)用物資化肥、農(nóng)膜和農(nóng)藥作為農(nóng)業(yè)非點(diǎn)源污染源變量，為更好反映農(nóng)戶施用水平，因此選以單位播種面積化肥、農(nóng)膜和農(nóng)藥投入施用量為指標(biāo)來表示農(nóng)業(yè)非點(diǎn)源污染強(qiáng)度［10］，用yi表示(i=1，2，3)，即 y1為化肥投入強(qiáng)度(kg/hm2)、y2為農(nóng)膜投入強(qiáng)度(kg/hm2)、y3為農(nóng)藥投入強(qiáng)度(kg/hm2)。反映農(nóng)村經(jīng)濟(jì)發(fā)展的指標(biāo)有很多，諸如農(nóng)業(yè)人均GDP、農(nóng)民人均純收入等指標(biāo)，然而化肥、農(nóng)膜和農(nóng)藥作為主要的農(nóng)業(yè)生產(chǎn)資料，其主要用于農(nóng)業(yè)活動，尤其是農(nóng)業(yè)種植業(yè)，對于我國目前仍實(shí)行的家庭聯(lián)產(chǎn)承包經(jīng)營的小戶來說，農(nóng)業(yè)種植業(yè)很重要的功能就是家庭食用品的供給保障，所以，農(nóng)戶施肥量、防蟲等決策主要受農(nóng)業(yè)種植業(yè)產(chǎn)出水平的影響。因此，選取x1為人均大農(nóng)業(yè)產(chǎn)值(元)和x2為人均種植業(yè)產(chǎn)值(元)來表示農(nóng)業(yè)經(jīng)濟(jì)增長水平更具有重要的現(xiàn)實(shí)意義。

2.2 數(shù)據(jù)來源與說明

為探究重慶直轄后的農(nóng)業(yè)非點(diǎn)源污染狀況，選取1996—2011年為研究時段，分別研究農(nóng)用化肥、農(nóng)膜和農(nóng)藥單位播種面積投入量與人均大農(nóng)業(yè)產(chǎn)值(元)、人均種植業(yè)產(chǎn)值(元)的關(guān)系。為消除價格因素對經(jīng)濟(jì)變量的影響，經(jīng)濟(jì)變量分別以1995年為基期計算成可比價。分析研究的指標(biāo)數(shù)據(jù)，依據(jù)《重慶統(tǒng)計年鑒1996—2012》、《2011、2012 年重慶市國民經(jīng)濟(jì)和社會發(fā)展統(tǒng)計公報》等公開發(fā)布的數(shù)據(jù)為基礎(chǔ)，見表1。本文中的人均大農(nóng)業(yè)產(chǎn)值以農(nóng)林牧漁業(yè)總產(chǎn)值為基礎(chǔ)計算，人均種植業(yè)產(chǎn)值是以狹義的農(nóng)業(yè)產(chǎn)值為基礎(chǔ)計算(以下稱人均小農(nóng)業(yè)產(chǎn)值)。

3 農(nóng)業(yè)非點(diǎn)源污染源與農(nóng)業(yè)經(jīng)濟(jì)增長關(guān)系的實(shí)證研究

3.1 模型構(gòu)建

根據(jù)環(huán)境庫茲涅茨曲線(EKC)相關(guān)理論，構(gòu)建農(nóng)業(yè)非點(diǎn)源污染的庫茲涅茨曲線模型。為檢驗(yàn)農(nóng)業(yè)非點(diǎn)源污染與農(nóng)業(yè)經(jīng)濟(jì)增長關(guān)系是否符合環(huán)境庫茲涅茨假說，本文以單位播種面積化肥、農(nóng)膜和農(nóng)藥投入施用量為因變量，以人均大農(nóng)業(yè)產(chǎn)值(元)和人均小農(nóng)業(yè)產(chǎn)值(元)作為自變量，分別研究每個因變量與各個自變量的關(guān)系。根據(jù)每個因變量與各個自變量的散點(diǎn)圖，分別用非“倒U型”特征一次函數(shù)和具有“倒U型”特征的二次或三次函數(shù)型［7］進(jìn)行模擬，檢驗(yàn)農(nóng)業(yè)非點(diǎn)源污染指標(biāo)與農(nóng)業(yè)經(jīng)濟(jì)增長變量是否符合環(huán)境庫茲涅茨曲線(EKC)假說。用于模擬每個因變量與各個自變量的模型如下:

表1 重慶市1996—2011年農(nóng)業(yè)非點(diǎn)源污染源與農(nóng)業(yè)經(jīng)濟(jì)增長變量相關(guān)匯總表Table 1 Data of agricultural NPS pollution and agricultural economic growth indicators of Chongqing from 1996 to 2011

式中:yi表示各種農(nóng)業(yè)非點(diǎn)源污染強(qiáng)度指標(biāo);β0表示常數(shù);βi表示回歸方程參數(shù);xi表示經(jīng)濟(jì)增長水平;μ為隨機(jī)影響因素。

3.2 實(shí)證檢驗(yàn)與分析

使用Eviews 6.0軟件用OLS估計方法，以表1數(shù)據(jù)用式(1)研究每個因變量與各個自變量關(guān)系。依據(jù)F檢驗(yàn)值來確定回歸模型的總體擬合情況，依據(jù)t統(tǒng)計量來判定回歸系數(shù)是否顯著，依據(jù)校正R2和估計標(biāo)準(zhǔn)誤差大小選擇最優(yōu)估計模型。

3.2.1 化肥投入強(qiáng)度與農(nóng)業(yè)經(jīng)濟(jì)增長的“EKC”檢驗(yàn)

采用OLS估計方法對單位播種面積化肥投入施用量與人均大農(nóng)業(yè)產(chǎn)值和人均小農(nóng)業(yè)產(chǎn)值關(guān)系，用式(1)進(jìn)行模擬，其最優(yōu)回歸估計模型見表2。由表2看到，在R2和F統(tǒng)計量均通過1%顯著水平情況下，單位播種面積化肥施用量與人均大農(nóng)業(yè)產(chǎn)值和人均小農(nóng)業(yè)產(chǎn)值關(guān)系均呈二次函數(shù)型，即呈現(xiàn)“倒U型”曲線形狀，具有“EKC”假設(shè)圖形的趨勢特征，見圖1。

圖1 單位播種面積化肥施用量與農(nóng)業(yè)經(jīng)濟(jì)增長變量擬合曲線Fig.1 Fitted curves of the relationship between the amount of chemical fertilizer per unit sown area and the agricultural economic growth indicators

與表1時間序列數(shù)據(jù)比較，“倒U型”曲線的轉(zhuǎn)折點(diǎn)在2009年與2011年之間，但更靠近2009年數(shù)據(jù)。由表1看到，重慶市單位播種面積化肥施用量自1996—2005年期間增速較快，年均增長率為2.58%，到2006—2011年增長減慢，年均增長率為0.8%并呈現(xiàn)向下的拐點(diǎn)(2010年左右);從全市單位播種面積化肥施用量分析結(jié)果看，以2009年的275.58 kg/hm2為最大值，2011年比2009年下降2.51 kg/hm2。綜上所述，重慶市單位播種面積化肥施用量與農(nóng)村經(jīng)濟(jì)發(fā)展變量的關(guān)系，具有“EKC”假設(shè)圖形的趨勢特征，并以2009年為拐點(diǎn)。單位播種面積化肥投入施用量與農(nóng)業(yè)經(jīng)濟(jì)增長變量呈現(xiàn)“倒U型”曲線，說明農(nóng)業(yè)經(jīng)濟(jì)增長對化肥投入施用量需求的規(guī)模效應(yīng)在2009年后開始脫鉤，農(nóng)業(yè)經(jīng)濟(jì)增長方式已從2009年前對化肥投入施用量規(guī)模效應(yīng)向投入要素的結(jié)構(gòu)效應(yīng)轉(zhuǎn)變。農(nóng)業(yè)經(jīng)濟(jì)增長方式正在由粗放型向集約型過渡，農(nóng)業(yè)產(chǎn)業(yè)結(jié)構(gòu)調(diào)整穩(wěn)步推進(jìn)，施肥結(jié)構(gòu)逐漸合理化，引進(jìn)化肥新品種，逐步擴(kuò)大實(shí)行測土配方施肥。

3.2.2 農(nóng)膜投入強(qiáng)度與農(nóng)業(yè)經(jīng)濟(jì)增長的“EKC”檢驗(yàn)

采用OLS估計方法對單位播種面積農(nóng)膜用量與人均大農(nóng)業(yè)增加值、人均小農(nóng)業(yè)增加值、人均大農(nóng)業(yè)產(chǎn)值和人均小農(nóng)業(yè)產(chǎn)值關(guān)系，用式(1)進(jìn)行模擬，其最優(yōu)回歸估計模型見表3。由表3看到，在R2和F統(tǒng)計量均通過1%顯著水平情況下，單位播種面積農(nóng)膜施用量與人均大農(nóng)業(yè)產(chǎn)值和人均小農(nóng)業(yè)產(chǎn)值關(guān)系均呈二次函數(shù)型，即呈現(xiàn)“倒U型”曲線形狀，具有“EKC”假設(shè)圖形的趨勢特征，見圖2。

圖2 單位播種面積農(nóng)膜施用量與農(nóng)業(yè)經(jīng)濟(jì)增長變量擬合曲線Fig.2 Fitted curves of the relationship between the amount of farming film per unit sown area and the agricultural economic growth indicators

表2 單位播種面積化肥施用量與農(nóng)業(yè)經(jīng)濟(jì)增長變量模型估計的最優(yōu)結(jié)果Table 2 Model-estimated optimal results of the relationship between the amount of fertilizer per unit sown area and the agricultural economic growth indicators

表3 單位播種面積農(nóng)膜施用量與農(nóng)業(yè)經(jīng)濟(jì)增長變量模型估計的最優(yōu)結(jié)果Table 3 Model-estimated optimal results of the relationship between the amount of farming film per unit sown area and the agricultural economic growth indicators

與表1時間序列數(shù)據(jù)比較，“倒U型”曲線的轉(zhuǎn)折點(diǎn)在2010年左右，靠近2009年數(shù)據(jù)，因此以2009年作為“倒U型”的拐點(diǎn)。因此，認(rèn)為重慶市單位播種面積農(nóng)膜用量與人均小農(nóng)業(yè)產(chǎn)值的關(guān)系，具有“EKC”假設(shè)圖形的趨勢特征，并以2009年為拐點(diǎn)。單位播種面積農(nóng)膜用量與農(nóng)業(yè)經(jīng)濟(jì)增長變量具有“EKC”假設(shè)的趨勢特征，說明農(nóng)業(yè)經(jīng)濟(jì)增長對農(nóng)膜用量需求的規(guī)模效應(yīng)在2009年后開始脫鉤。其主要原因大致歸結(jié)為農(nóng)膜利用效率提高和生態(tài)農(nóng)業(yè)的發(fā)展。提高農(nóng)膜利用效率途徑，一是通過土地整理，使土地小塊經(jīng)營向規(guī)模經(jīng)營轉(zhuǎn)變，加上農(nóng)膜正確安裝和使用;二是正確選用和科學(xué)使用各種專用膜產(chǎn)品。單位播種面積農(nóng)膜用量與大農(nóng)業(yè)經(jīng)濟(jì)變量的擬合度要優(yōu)于與小農(nóng)業(yè)的經(jīng)濟(jì)變量的擬合度。這主要是農(nóng)膜不僅使用于小農(nóng)業(yè)，在林業(yè)、畜牧業(yè)上也正在廣泛使用，如用于苗木的繁育，大棚養(yǎng)殖雞、鴨、鵝等。

3.2.3 農(nóng)藥投入強(qiáng)度與農(nóng)業(yè)經(jīng)濟(jì)增長的“EKC”檢驗(yàn)

采用OLS估計方法對單位播種面積農(nóng)藥施用量與人均大農(nóng)業(yè)產(chǎn)值和人均小農(nóng)業(yè)產(chǎn)值關(guān)系，用式(1)進(jìn)行模擬，其最優(yōu)回歸估計模型見表4。由表4看到，在R2和F統(tǒng)計量均通過1%顯著水平情況下，均呈二次函數(shù)型，即呈現(xiàn)“倒U型”曲線形狀，具有“EKC”假設(shè)圖形的趨勢特征，見圖3。

圖3 單位播種面積農(nóng)藥施用量與農(nóng)業(yè)經(jīng)濟(jì)增長變量擬合曲線Fig.3 Fitted curves of the relationship between the amount of pesticide per unit sown area and the agricultural economic growth indicators

與表1時間序列數(shù)據(jù)比較，“倒U型”曲線的轉(zhuǎn)折點(diǎn)均在2009年左右，靠近2009年數(shù)據(jù)，因此以2009年作為“倒U型”的拐點(diǎn)。因此，認(rèn)為重慶市單位播種面積農(nóng)藥施用量與農(nóng)業(yè)經(jīng)濟(jì)增長變量呈現(xiàn)“倒U型”曲線，具有“EKC”假設(shè)圖形的趨勢特征，并以2009年為拐點(diǎn)。說明農(nóng)業(yè)經(jīng)濟(jì)增長對農(nóng)藥用量需求的規(guī)模效應(yīng)在2009年后開始脫鉤。同時，單位播種面積農(nóng)藥施用量與大農(nóng)業(yè)經(jīng)濟(jì)變量的擬合度要優(yōu)于與小農(nóng)業(yè)的經(jīng)濟(jì)變量的的擬合度，這與農(nóng)藥不僅用于小農(nóng)業(yè)，而且還用于林業(yè)、草業(yè)病蟲害防治有關(guān)。

4 結(jié)論與討論

通過對重慶市1996—2011年農(nóng)業(yè)經(jīng)濟(jì)增長變量分別與農(nóng)業(yè)非點(diǎn)源污染源單位播種面積投入強(qiáng)度的統(tǒng)計資料，采用環(huán)境庫茲涅茨曲線模型，探討重慶市農(nóng)業(yè)經(jīng)濟(jì)增長與農(nóng)業(yè)非點(diǎn)源污染源單位播種面積投入強(qiáng)度的相關(guān)關(guān)系。

(1)單位播種面積化肥、農(nóng)膜、農(nóng)藥施用量與農(nóng)業(yè)經(jīng)濟(jì)增長變量關(guān)系均呈二次函數(shù)型，即呈現(xiàn)“倒U型”曲線形狀，具有“EKC”假設(shè)圖形的趨勢特征，并以2009年為拐點(diǎn)，從而證明重慶市農(nóng)業(yè)非點(diǎn)源污染與農(nóng)業(yè)經(jīng)濟(jì)增長存在“強(qiáng)相關(guān)—減排—脫鉤”的關(guān)系，符合環(huán)境庫茲涅茨(EKC)理論。重慶市農(nóng)村經(jīng)濟(jì)發(fā)展水平低，單位播種面積化肥、農(nóng)藥施用量與農(nóng)業(yè)經(jīng)濟(jì)增長變量開始脫鉤。

(2)從污染結(jié)構(gòu)來看，農(nóng)資污染尤其是農(nóng)膜污染在較長時間內(nèi)呈現(xiàn)持續(xù)上升趨勢，應(yīng)成為重慶市農(nóng)業(yè)非點(diǎn)源的管理重心。

(3)單位播種面積化肥、農(nóng)膜、農(nóng)藥施用量與大農(nóng)業(yè)經(jīng)濟(jì)變量的擬合度要明顯優(yōu)于與小農(nóng)業(yè)經(jīng)濟(jì)變量的的擬合度，說明大農(nóng)業(yè)的產(chǎn)業(yè)組合效應(yīng)或交互效應(yīng)促進(jìn)了單位播種面積化肥、農(nóng)膜、農(nóng)藥施用量與大農(nóng)業(yè)經(jīng)濟(jì)變量關(guān)系更緊密。

(4)在“大農(nóng)村、大庫區(qū)”市情背景下，地處低山丘陵巖溶區(qū)，單位播種面積化肥、農(nóng)膜、農(nóng)藥施用量與農(nóng)業(yè)經(jīng)濟(jì)增長關(guān)系呈現(xiàn)“倒U型”曲線形狀，說明走科學(xué)發(fā)展之路，轉(zhuǎn)變經(jīng)濟(jì)發(fā)展方式，調(diào)整農(nóng)業(yè)產(chǎn)業(yè)結(jié)構(gòu)，實(shí)施“沃土工程”建設(shè)，大力發(fā)展生態(tài)循環(huán)農(nóng)業(yè)，農(nóng)業(yè)生態(tài)園區(qū)等是促進(jìn)農(nóng)用化肥、農(nóng)膜、農(nóng)藥“減量化”的有效途徑。如何鞏固這一來之不易的成果，進(jìn)一步促進(jìn)農(nóng)用化肥、農(nóng)膜和農(nóng)藥“減量化”成為重慶市農(nóng)業(yè)生產(chǎn)方式轉(zhuǎn)變的重要內(nèi)容。

表4 單位播種面積農(nóng)藥施用量與農(nóng)業(yè)經(jīng)濟(jì)增長變量模型估計的最優(yōu)結(jié)果Table 4 Model-estimated optimal results of the relationship between the amount of pesticide per unit sown area and the agricultural economic growth indicators

(5)在治理農(nóng)業(yè)非點(diǎn)源污染時，首先應(yīng)進(jìn)一步加大對農(nóng)業(yè)科技創(chuàng)新投入，大力推廣資源節(jié)約型、環(huán)境友好型農(nóng)業(yè)新技術(shù)，進(jìn)一步加大力度推廣測土配方施肥技術(shù)擴(kuò)大應(yīng)用面積，科學(xué)選用和正確使用農(nóng)膜，科學(xué)施用農(nóng)藥，減輕副作用發(fā)生，從多個方面采取多種措施，從而為農(nóng)民合理施肥，減少對環(huán)境的危害;其次，重慶市處于三峽庫區(qū)上游，其農(nóng)業(yè)非點(diǎn)源污染將會對環(huán)境造成巨大危害，因此應(yīng)構(gòu)建農(nóng)業(yè)生態(tài)補(bǔ)償機(jī)制，引導(dǎo)農(nóng)業(yè)產(chǎn)業(yè)的升級，轉(zhuǎn)向生態(tài)農(nóng)業(yè)、生態(tài)園區(qū)等方向發(fā)展;第三，加大農(nóng)村環(huán)境保護(hù)宣傳力度，提高農(nóng)民的現(xiàn)代文明素質(zhì)。通過舉辦知識講座、科技人員現(xiàn)場指導(dǎo)、輿論媒體的廣泛宣傳，使農(nóng)戶相信科學(xué)，充分認(rèn)識到過量施用化肥、不合理使用農(nóng)膜、盲目濫用農(nóng)藥的危害性和安全用藥的必要性、迫切性。以實(shí)施“科技入戶”工程為載體，創(chuàng)新農(nóng)技服務(wù)方式，整合培訓(xùn)資源，加強(qiáng)農(nóng)民科技培訓(xùn)，全面提高農(nóng)民的現(xiàn)代文明素質(zhì)。

［1］張維理，徐愛國，冀宏杰，等.中國農(nóng)業(yè)面源污染形勢估計及控制對策Ⅲ:中國農(nóng)業(yè)面源污染控制中存在問題分析［J］.中國農(nóng)業(yè)科學(xué)，2004，37(7):1026 －1033.(ZHANG Wei-li，XU Ai-guo，JI Hong-jie，et al.Estimation of Agricultural Non-point Source Pollution in China and the Alleviating StrategiesⅢ.Status of Agricultural Non-point Source Pollution and the Alleviating Strategies in China［J］.Scientia Agricultura Sinica，2004，37(7):1026 －1033.(in Chinese))

［2］胡 聘，許開鵬，楊建新，等.經(jīng)濟(jì)發(fā)展對環(huán)境質(zhì)量的影響——環(huán)境庫茲涅茨曲線國內(nèi)外研究進(jìn)展［J］.生態(tài)學(xué)報，2004，24(6):1260 － 1265.(HU Pin，XU Kaipeng，YANG Jian-xin，et al.How Does Economic Development Influence on the Environment Quality:The Environmental Kuznets Curve’s Research Progress Home and Abroad［J］.Acta Ecologica Sinica，2004，24(6):1260 －1265.(in Chinese))

［3］宋 濤，鄭挺國，佟連軍.基于面板協(xié)整的環(huán)境庫茲茨曲線的檢驗(yàn)與分析［J］.中國環(huán)境科學(xué)，2007，27(4):572 －575.(SONG Tao，ZHENG Ting-guo，TONG Lianjun.Test and Analysis on the Environment Kuznets Curve Based on Panel Co-integration［J］.Journal of China Environmental Science，2007，27(4):572 － 575.(in Chinese))

［4］劉 楊，陳邵峰，張云芳.中國農(nóng)業(yè)EKC研究:以化肥為例［J］.中國農(nóng)學(xué)通報，2009，25(16):260－265.(LIU Yang，CHEN Shao-feng，ZHANG Yun-fang.Research on China Agricultural EKC with Chemical Fertilizer for example［J］.Chinese Agricultural Science Bulletin，2009，25(16):260 －265.(in Chinese))

［5］CHIMELI A B.Growth and the Environment:Are We Looking at the Right Data［J］.Economics Letters，2007，96(1):90－95.

［6］彭水軍，包 群.經(jīng)濟(jì)增長與環(huán)境污染——環(huán)境庫茲涅茨曲線假說的中國檢驗(yàn)［J］.財經(jīng)問題研究，2006，(8):3 － 15.(PENG Shui-jun，BAO Qun.Economic Growth and Environmental Pollution:An Empirical Test for the Environmental Kuznets Curve Hypothesis in China［J］.Research on Financial and Economic Issues，2006，(8):3 －15.(in Chinese))

［7］李海鵬，張俊飚.中國農(nóng)業(yè)面源污染與經(jīng)濟(jì)發(fā)展關(guān)系的實(shí)證研究［J］.長江流域資源與環(huán)境，2009，18(6):586 －590.(LI Hai-peng，ZHANG Jun-biao.An Empirical Test on the“EKC”Relationship between the Agricultural Non-point Source Pollution and Economic Development in China［J］.Resources and Environment in the Yangtze Basin，2009，18(6):586 －590.(in Chinese))

［8］牟新利，祁俊生，黃 宇，等.重慶農(nóng)業(yè)面源污染的環(huán)境庫茲涅茨曲線評價［J］.貴州農(nóng)業(yè)科學(xué)，2011，39(1):228 － 230.(MOU Xin-li，QI Jun-sheng，HUANG Yu， etal.Pollution ofAgriculturalEnvironmental Kuznets Non-point Source Based on Curve Model in Chongqing［J］.Guizhou Agricultural Sciences，2011，39(1):228 －230.(in Chinese))

［9］楊志敏，陳玉成，魏世強(qiáng)，等.重慶市農(nóng)業(yè)面源污染影響因子的系統(tǒng)分析［J］.農(nóng)業(yè)環(huán)境科學(xué)學(xué)報，2009，28(5):1001 － 1004.(YANG Zhi-min，CHEN Yu-cheng，WEI Shi-qiang，et al.System Analysis of Chongqing Agricultural Non-point Source Pollution Influence Factor［J］.Journal of Agro-environment Science，2009，28(5):1001 －1004.(in Chinese))

［10］鞏前文，田志宏.農(nóng)村經(jīng)濟(jì)發(fā)展與化肥使用量的庫茲涅茨曲線假說及驗(yàn)證［J］.中國農(nóng)業(yè)大學(xué)學(xué)報(社會科學(xué)版)，2011，27(4):157 －160.(GONG Qian-wen，Tian Zhi-hong.The Kuznets Curve’s Hypothesis and Validation between Rural Economic Development and Fertilizer Usage［J］.Journal of China Agricultural University(Social Science)，2011，27(4):157 －160.(in Chinese ))