三峽庫區(qū)面源污染的時空特征及EKC分析

馮 琳,張婉婷,張鈞珂,毛馨卉,于 鑫,葉成松,常天陽

三峽庫區(qū)面源污染的時空特征及EKC分析

馮 琳1*,張婉婷1,張鈞珂1,毛馨卉2,于 鑫3,葉成松3,常天陽1

(1.中國人民大學環(huán)境學院,北京 100872；2.中國自然資源經(jīng)濟研究院,北京 101149；3.廈門大學環(huán)境與生態(tài)學院,福建 廈門 361102)

在全面核算2008～2018年三峽庫區(qū)19個區(qū)縣農(nóng)業(yè)面源污染TN、TP排放量的基礎(chǔ)上,分析了其時空變化特征.基于環(huán)境庫茲涅茨曲線(EKC)理論,構(gòu)建了基于面板數(shù)據(jù)的回歸模型,探究了庫區(qū)農(nóng)業(yè)經(jīng)濟發(fā)展與面源污染排放強度的演替關(guān)系.結(jié)果顯示,庫區(qū)農(nóng)業(yè)TN排放波動減少,TP波動增加.各區(qū)縣的TN和TP年均排放量分別在374～6046t和105～1267t之間.其中,農(nóng)田化肥與畜禽養(yǎng)殖單元的總產(chǎn)污貢獻率達80%以上.庫區(qū)TN排放強度、畜禽養(yǎng)殖與農(nóng)村生活單元的TN、TP排放強度均存在顯著的“倒U型”EKC關(guān)系,目前已跨越拐點.農(nóng)田化肥TP排放強度、水產(chǎn)養(yǎng)殖與農(nóng)田固廢單元的TN、TP排放強度呈現(xiàn)顯著的“直線型”EKC關(guān)系,處于與經(jīng)濟同步增長的階段.建議重點升級農(nóng)田化肥單元的污染防控能力,以配套推進農(nóng)村人居環(huán)境的改善,促進區(qū)域氮磷減排.

農(nóng)業(yè)面源污染；單元分析法；環(huán)境庫茲涅茨曲線(EKC)；三峽庫區(qū)

全國污染源普查數(shù)據(jù)顯示,農(nóng)業(yè)面源已超過工業(yè)源和生活源,成為我國水體污染的第一大來源[1].與點源相比,農(nóng)業(yè)面源分散、面廣、來源和途徑復雜,防治難度更大,已成為許多地區(qū)水環(huán)境問題的重要原因[2-4].

三峽庫區(qū)是長江經(jīng)濟帶的重要組成部分,戰(zhàn)略地位重要.它集農(nóng)村、山區(qū)、移民區(qū)為一體,生態(tài)環(huán)境保護與移民安穩(wěn)致富等問題交叉[5-6],面源污染防治任重而道遠.根據(jù)全國生態(tài)環(huán)境狀況公報,2019年庫區(qū)主要支流的77個監(jiān)測斷面中,富營養(yǎng)、中營養(yǎng)狀態(tài)分別占了20.8%、77.9%.在2021年中央一號文件、《全國農(nóng)業(yè)可持續(xù)發(fā)展規(guī)劃(2015-2030年)》中,三峽庫區(qū)均被列為開展面源污染綜合防治示范建設(shè)的重點地區(qū)之一.

關(guān)于農(nóng)業(yè)面源污染負荷的評估,主要有經(jīng)驗過程模型和機理模型兩類方法.前者包括統(tǒng)計法、水文分離法、輸出系數(shù)法、單元分析法[7-8]等;后者包括AnnAGNPS、SWAT方法等[9].目前,針對三峽庫區(qū)蘭陵溪、香溪河、筍溪河等小流域[10-12],以及重慶段、湖北段[13-15]分段負荷評估的研究較豐富,為揭示庫區(qū)農(nóng)業(yè)面源污染的部分排放特征作出了重要貢獻,但已有研究空間上多局限于重慶段或湖北段,時間上多關(guān)注前三峽時期(2008年以前)或者后三峽時期的某個具體年份,尚缺乏將重慶、湖北庫區(qū)作為整體,并貫穿后三峽時期的面源負荷評估.

根據(jù)污染負荷的評估值,探究經(jīng)濟增長與面源污染水平的內(nèi)在聯(lián)系,有助于識別各污染單元潛在的問題,為配套設(shè)計面源污染防治政策提供科學依據(jù)[16-18].“環(huán)境庫茲涅茨曲線”(EKC)為環(huán)境-經(jīng)濟關(guān)系的實證研究提供了重要理論[19-21].雖然傳統(tǒng)的EKC曲線形狀一般為“倒U型”,但從全國和地區(qū)的實證研究結(jié)果看,農(nóng)業(yè)面源污染的EKC曲線有“倒U型”[22-23]、“U型”[24-25]、甚至是“N型”、“倒N型”[26-27]、線性增長[28]等多種可能[29].目前關(guān)于三峽庫區(qū)面源EKC的研究,分析維度單一,多用時間序列數(shù)據(jù)[30-32].如果采用兼具時空信息的面板數(shù)據(jù),將可增加數(shù)據(jù)的自由度并降低解釋變量間的共線性,彌補時序和截面數(shù)據(jù)的不足[33],提高EKC模型估計的有效性.

基于以上背景,本文拓展評估的時空尺度,全面估算庫區(qū)(重慶段及湖北段)在后三峽時期(2008～ 2018年)的農(nóng)業(yè)面源污染負荷,識別重點防控地域和污染單元.基于評估獲得的面板數(shù)據(jù)求得污染排放強度,借助EKC理論,探究庫區(qū)農(nóng)業(yè)面源污染與農(nóng)業(yè)經(jīng)濟增長的關(guān)系,闡釋EKC的現(xiàn)實意義,以期為庫區(qū)和長江流域農(nóng)業(yè)面源污染防治提供有益的參考.

1 研究方法

1.1 面源污染排放量核算

由于三峽庫區(qū)涉及的行政區(qū)縣多,地理尺度較大,從數(shù)據(jù)的可獲得性角度,本文擬采用官方統(tǒng)計數(shù)據(jù),運用單元分析法[12]對庫區(qū)農(nóng)業(yè)面源TN、TP污染的排放量進行估算.計算公式如下:

式中:ANP為進入水系的農(nóng)業(yè)面源污染物的排放量;N為單元指標統(tǒng)計數(shù);為單元污染物的產(chǎn)污系數(shù);C(N,,)表示單元污染物的入河系數(shù),由單元特性(N)、資源利用率()和區(qū)域環(huán)境地理特征()共同決定.

各產(chǎn)污單元TN、TP排放量的具體公式和系數(shù)說明如表1所示.

表1 各產(chǎn)污單元TN、TP污染排放量計算公式與系數(shù)

注:*《第一次全國污染源普查水產(chǎn)養(yǎng)殖業(yè)污染物產(chǎn)排污系數(shù)手冊》.

1.2 EKC模型

利用面板數(shù)據(jù)進行EKC模型回歸時,由于各區(qū)域在地理環(huán)境等方面存在異質(zhì)性,所以被解釋變量宜采用強度變量.考慮到EKC曲線的形狀有多種可能性,先將其設(shè)為三次曲線,若不顯著,再設(shè)為二次曲線,若還不顯著,則擬合成線性模型[37].為減少異方差,本文對各指標進行對數(shù)處理,具體模型構(gòu)建如下:

參與回歸的具體指標設(shè)置如表2.為消除通脹影響,借助農(nóng)林牧漁總產(chǎn)值指數(shù),將農(nóng)業(yè)產(chǎn)值調(diào)整為以2008年為基期的實際值.面板數(shù)據(jù)的估計模型主要包括混合模型、隨機效應模型及固定效應模型[28].依次利用Eviews軟件中的F-test、Hausman-test、LM-test功能對3類模型進行選擇.F-test用于檢驗固定效應與混合效應的適用性,Hausman-test用于檢驗固定效應與隨機效應的適用性,LM-test用于檢驗隨機效應與混合效應的適用性.F-test如果值小于0.05,則拒絕混合效應,固定效應模型更優(yōu). Hausman-test如果值小于0.05,則拒絕隨機效應,固定效應模型更優(yōu).LM-test如果值小于0.05,則拒絕混合效應,隨機效應模型更優(yōu).

表2 農(nóng)業(yè)面源污染指標與經(jīng)濟指標設(shè)置情況

1.3 數(shù)據(jù)來源

本文的研究對象涉及三峽庫區(qū)重慶段的15個區(qū)縣—巫溪、開州、巫山、云陽、奉節(jié)、萬州、忠縣、石柱、豐都、長壽、武隆、涪陵、巴南、渝北、江津,以及湖北段的4個區(qū)縣—夷陵、秭歸、興山、巴東.各區(qū)縣的化肥折純施用量、畜禽飼養(yǎng)量、各作物產(chǎn)量、水產(chǎn)品產(chǎn)量、農(nóng)村常住人口、農(nóng)林牧漁總產(chǎn)值和農(nóng)村居民人均可支配收入數(shù)據(jù),來源于2009～2019年重慶各區(qū)縣的統(tǒng)計年鑒和《湖北農(nóng)村統(tǒng)計年鑒》.各區(qū)縣的農(nóng)林牧漁總產(chǎn)值指數(shù),來源于歷年的《重慶市統(tǒng)計年鑒》和《湖北省統(tǒng)計年鑒》;地表水資源量由歷年的重慶市水資源公報和湖北省水資源公報獲得.本文其余數(shù)據(jù),如各污染單元的產(chǎn)污系數(shù)、入河系數(shù)等則主要參考相關(guān)文獻[12,18,39-41].

2 結(jié)果與討論

2.1 三峽庫區(qū)農(nóng)業(yè)面源污染排放時空特征

2.1.1 時間尺度 利用1.1部分的模型、參數(shù)和數(shù)據(jù),得到2008～2018年三峽庫區(qū)農(nóng)業(yè)面源的TN和TP排放量(圖1).TN從2008年的41248.62t減少到2018年的38527.86t,年均排放量39770.55t.11a間TN排放整體呈顯著下降趨勢,擬合曲線的斜率為-166.9t/ a (sig=0.000).TP從2008年的7347.19t增加到2018年的8210.21t,年均排放量為8795.23t.11a間TP排放先升后降,但總體呈上升趨勢,擬合曲線的斜率為208.03t/a (sig=0.000).其中,2014年TP排放的年增幅最大,為31.13%.開州的磷肥施用量在當年驟然增加,是造成這一結(jié)果的主要原因.

如圖2所示,5個污染單元對庫區(qū)農(nóng)業(yè)面源TN、TP的排放貢獻由高到低依次為:農(nóng)田化肥>畜禽養(yǎng)殖>農(nóng)村生活>水產(chǎn)養(yǎng)殖>農(nóng)田固廢.該結(jié)果與江西省面源污染貢獻率“農(nóng)村生活>畜禽養(yǎng)殖>種植業(yè)>水產(chǎn)養(yǎng)殖”[37]有所不同,也與漢江流域荊門段TN貢獻率“農(nóng)村生活>農(nóng)業(yè)用地>畜禽養(yǎng)殖”,TP貢獻率“畜禽養(yǎng)殖>農(nóng)村生活>農(nóng)業(yè)用地”[38]有著較大的差異.三峽庫區(qū)農(nóng)田化肥和畜禽養(yǎng)殖兩個單元對TN與TP的貢獻之和大于80%,該結(jié)果與張廣納等[14]、蔡金洲等[15]學者的發(fā)現(xiàn)一致.

各單元污染排放量隨時間變化的特點為:

農(nóng)田化肥的TN排放整體上波動減少,從2008年30407.67t下降到2018年28006.89t,降幅約8%;TP排放雖然總體增幅達13.3%,但在2015～2018年期間持續(xù)下降.以上結(jié)果可能與國家2015年頒布《到2020年化肥使用量零增長行動方案》,庫區(qū)隨后積極開展橘渣飼料化、肥料化和沼液利用等技術(shù)研究,減施化肥,增施有機肥有關(guān).

畜禽養(yǎng)殖的TN排放略有減少,2018年減至5163.34t,降幅約0.58%;TP排放波動上升,2018年增至1908.91t,總體增幅約12.4%.這與十二五期間豐都曾大力發(fā)展肉牛養(yǎng)殖業(yè),江津、開州、云陽和萬州新開了一些大型生豬養(yǎng)殖廠有關(guān).

農(nóng)村生活的TN排放變化不大,TP排放略有下降.水產(chǎn)養(yǎng)殖、農(nóng)田固廢的TN與TP排放略有上升.

2.1.2 空間尺度 2008～2018年間,三峽庫區(qū)的19個區(qū)縣農(nóng)業(yè)面源TN、TP的年均排放量空間差異明顯(圖3).

各區(qū)縣TN年均排放量在374～6046t之間.其中,開州最高(6046t),其次為涪陵、江津、萬州(3000～4000t),渝北、巴南和興山最少(小于1000t).渝北、巴南在重慶主城區(qū)外圍,城鎮(zhèn)化建設(shè)較為發(fā)達,農(nóng)業(yè)生產(chǎn)活動較少;而湖北興山地形復雜,“八分半山一分田”,種植業(yè)較少,所以這3個區(qū)縣的TN排放量比較低.

各區(qū)縣TP年均排放量在105～1267t.其中,開州最高(1267t),其次為夷陵、萬州、江津、云陽(600～ 800t),巴南和興山最少(小于200t).相對于TN排放量的地區(qū)排序,夷陵、渝北的TP次序前移,其農(nóng)用磷肥施用量較多可能是主要原因.

圖3 三峽庫區(qū)2008～2018年農(nóng)業(yè)面源污染TN及TP年均排放量

5個污染單元對各個區(qū)縣TN和TP的年均排放貢獻,如圖4所示.除了巴南的農(nóng)村生活單元對TN的貢獻率大于畜禽養(yǎng)殖單元,奉節(jié)、巴南的畜禽養(yǎng)殖單元對TP的貢獻率大于農(nóng)田化肥以外,其他區(qū)縣5個污染單元對農(nóng)業(yè)面源TN、TP排放量的貢獻率由大到小的順序均為:農(nóng)田化肥>畜禽養(yǎng)殖>農(nóng)村生活>水產(chǎn)養(yǎng)殖>農(nóng)田固廢.

庫區(qū)絕大多數(shù)區(qū)縣農(nóng)田化肥單元對TN、TP的年均貢獻率超過了50%.其中,化肥單元對TN貢獻率最高的4個區(qū)縣為:涪陵(81.88%)、開州(81.61%)、秭歸(81.50%)和夷陵(80.07%);對TP貢獻率最高的4個區(qū)縣為:秭歸(84.60%),夷陵(81.11%)、開州(77.84%)、巫溪(76.65%).畜禽養(yǎng)殖單元對TN年均貢獻率最高的4個區(qū)縣為:豐都(32.95%)、興山(30.70%)、忠縣(26.59%)、巴南(25.35%);對TP貢獻率最高的4個區(qū)縣為:奉節(jié)(44.05%),涪陵(38.63%)、豐都(38.34%)、興山(36.12%).

各區(qū)縣在農(nóng)田化肥與畜禽養(yǎng)殖這兩個主要產(chǎn)污單元TN、TP的年均排放量,如圖5所示.可以看出,在19個區(qū)縣中,開州的化肥TN、TP年均排放量最多,豐都的畜禽養(yǎng)殖TN、TP年均排放量最多.庫區(qū)的農(nóng)業(yè)面源污染防治工作,對這2個區(qū)縣需要給予較多的關(guān)注.

圖4 2008～2018年污染單元對三峽各區(qū)縣TN、TP排放的貢獻

審圖號:GS(2021)6024號

2.2 三峽庫區(qū)農(nóng)業(yè)面源污染的EKC曲線分析

2.2.1 面板單位根檢驗 面板數(shù)據(jù)兼具截面數(shù)據(jù)和時間序列數(shù)據(jù)的雙重特性,為了識別面板數(shù)據(jù)中各組變量的平穩(wěn)性、避免偽回歸現(xiàn)象的存在,常采用LLC、IPS、Fisher等方法進行單位根檢驗.本文對庫區(qū)及各污染單元的TNS、TPS和AGDP數(shù)據(jù)進行了上述檢驗,得到的平穩(wěn)性結(jié)論一致,限于篇幅文中僅列出LLC檢驗結(jié)果,如表3所示.所有變量的原始序列均拒絕了原假設(shè),各變量均沒有單位根,說明數(shù)據(jù)是長期平穩(wěn)的,無需進行協(xié)整檢驗,可以構(gòu)建面板模型.

2.2.2 庫區(qū)及各產(chǎn)污單元的EKC擬合結(jié)果 由于各組數(shù)據(jù)的三次曲線模型均未通過10%的顯著水平檢驗,故先設(shè)為二次曲線.如不顯著,再設(shè)為一次曲線,回歸結(jié)果如表4所示.

TNS、CQTNS、CQTPS、SHTNS、SHTPS與經(jīng)濟發(fā)展呈顯著的“倒U型”關(guān)系,總體上模型擬合效果較好,說明伴隨著農(nóng)業(yè)發(fā)展,庫區(qū)農(nóng)業(yè)面源的TN、畜禽養(yǎng)殖、農(nóng)村生活的污染強度表現(xiàn)為“升-降”趨勢.TNS擬合的2值相對較小,可能是由于面源污染產(chǎn)生的機制較為復雜,僅將AGDP作為解釋變量,對該模型回歸結(jié)果的解釋能力不足,種植結(jié)構(gòu)、農(nóng)業(yè)機械投入強度等因素或許應增設(shè)為該模型重要的控制變量[33].

表3 三峽庫區(qū)農(nóng)業(yè)面源及各污染單元TN、TP排放強度的面板單位根(LLC)檢驗結(jié)果

注:*、**和***分別表示在10%、5%和1%水平顯著;括號內(nèi)數(shù)據(jù)為值.

表4 三峽庫區(qū)農(nóng)業(yè)面源及各產(chǎn)污單元TN、TP排放強度的EKC回歸結(jié)果

注:*、**和***分別表示在10%、5%和1%水平顯著.

經(jīng)過一階求導和自然底數(shù)的指數(shù)函數(shù)求值,得到TNS曲線拐點對應的AGDP為6933.96元, CQTNS、CQTPS的拐點分別為7464.05、7670.93元,SHTNS、SHTPS的拐點為5268.38元.對比三峽各年份的AGDP,庫區(qū)于2011年跨越了TNS的拐點,2012年跨越了CQTNS、CQTPS的拐點,2008年跨越了SHTNS、SHTPS的拐點.說明庫區(qū)在畜禽養(yǎng)殖與農(nóng)村生活單元的環(huán)境污染與經(jīng)濟發(fā)展已基本脫鉤.究其原因,一方面是由于當?shù)剞r(nóng)民對畜禽養(yǎng)殖收入的依賴程度降低,以及外出打工使得農(nóng)村常住人口減少.另一方面,庫區(qū)兩級政府自2008年實施了農(nóng)村居民點污水處理示范項目,對這兩個單元的拐點跨越作出了積極貢獻.

HFTPS、SCTNS、SCTPS、GFTNS、GFTPS與經(jīng)濟發(fā)展呈顯著的“直線型”關(guān)系,說明庫區(qū)農(nóng)田化肥的TP、水產(chǎn)養(yǎng)殖、農(nóng)田固廢的污染強度處于與經(jīng)濟同步增長的耦合階段,需加以重視并科學防控.農(nóng)用化肥單元TPS的一次曲線雖然模型顯著,但擬合的2值較小.這可能是由于各區(qū)縣在種植結(jié)構(gòu)、磷肥施用種類、農(nóng)業(yè)政策等方面存在差異,產(chǎn)生了門檻效應[29],從而使得ADGP在跨越門檻值前后對農(nóng)田化肥TPS的影響有所不同,進而呈現(xiàn)出不完全的線性關(guān)系. TPS、HFTNS與經(jīng)濟發(fā)展之間的關(guān)系不顯著.

對于模型的穩(wěn)健性檢驗,本文選用變量替換法[39].鑒于人均收入是衡量經(jīng)濟發(fā)展的常用指標之一,且曾用于世行EKC公式的推導以及北美自由貿(mào)易協(xié)定的環(huán)境效應分析[19],本文將解釋變量替換為農(nóng)村居民人均可支配收入,重新進行回歸分析.結(jié)果顯示新模型的回歸系數(shù),在顯著性、方向和大小關(guān)系上與原回歸模型基本一致,說明本研究的模型估計結(jié)果具有一定的穩(wěn)健性.

3 結(jié)論及建議

3.1 結(jié)論

3.1.1 面源排放 2008～2018年,庫區(qū)TN排放波動減少,年均排放量39770.55t;TP排放波動增加,年均排放量為8795.23t.各區(qū)縣的TN、TP年均排放量,分別在374～6046t和105～1267t之間.農(nóng)田化肥和畜禽養(yǎng)殖的總產(chǎn)污貢獻率達80%以上.在庫區(qū)19個區(qū)縣中,TN年均排放量最多的前4名區(qū)縣為開州、涪陵、江津和萬州,TP年均排放量最多的前4名為開州、夷陵、萬州和江津.開州與豐都分別是農(nóng)田化肥污染、畜禽養(yǎng)殖污染最嚴重的區(qū)縣.

3.1.2 面源排放強度的EKC特征 基于2008～ 2018年的面板數(shù)據(jù),三峽庫區(qū)的TN排放強度、畜禽養(yǎng)殖與農(nóng)村生活單元的TN、TP排放強度,與農(nóng)業(yè)經(jīng)濟發(fā)展呈顯著的“倒U型”關(guān)系,跨越拐點的時間分別為2011年、2012年、2008年.畜禽養(yǎng)殖與農(nóng)村生活單元的環(huán)境污染與經(jīng)濟發(fā)展已基本脫鉤.農(nóng)田化肥的TP排放強度、水產(chǎn)養(yǎng)殖與農(nóng)田固廢單元的TN、TP排放強度與經(jīng)濟發(fā)展呈顯著的“直線型”關(guān)系,它們處于與經(jīng)濟同步增長的階段,需加以重視并科學防控.

3.2 建議

農(nóng)業(yè)面源污染的控制,不僅是我國后續(xù)水污染防治工作的關(guān)鍵,還是實現(xiàn)鄉(xiāng)村振興戰(zhàn)略中“生態(tài)宜居”要求的重要組成部分.綜合以上分析和背景,建議重點升級三峽庫區(qū)農(nóng)田化肥的污染防控能力,以配套推進農(nóng)村人居環(huán)境的改善.可采取的措施包括:適當調(diào)整種植結(jié)構(gòu),因地制宜地推廣精準施肥、高效緩釋肥、測土配方施肥、有機肥,研發(fā)推廣適用施肥設(shè)備,機械深施,提高化肥利用率.

[1] 生態(tài)環(huán)境部,國家統(tǒng)計局,農(nóng)業(yè)農(nóng)村部.第二次全國污染源普查公報[R]. 2020.6.

Ministry of Ecology and Environment, National Bureau of Statistics, Ministry of Agricutural and Rural Affairs. Second National Pollution Source Census Bulletin [R]. 2020.6.

[2] Han D M, Currell M J, Cao G L. Deep challenges for China's war on water pollution [J]. Environmental Pollution, 2016,218:1222-1233.

[3] 歐陽威,鞠欣妍,高 翔,等.考慮面源污染的農(nóng)業(yè)開發(fā)流域生態(tài)安全評價研究 [J]. 中國環(huán)境科學, 2018,38(3):1194-1200.

Ou Y W, Ju X Y, Gao X, et al. Ecological security assessment of agricultural development watershed considering non-point source pollution [J]. China Environmental Science, 2018,38(3):1194-1200.

[4] 楊金鳳,馮愛萍,王雪蕾,等.海河流域農(nóng)業(yè)面源污染潛在風險識別方法[J]. 中國環(huán)境科學, 2021,41(10):4782-4791.

Yang J F, Feng A P, Wang X L, et al. An identification method of potential risk for agricultural non-point source pollution in the Haihe River Basin [J]. China Environmental Science, 2021,41(10):4782- 4791.

[5] Feng L, XuJ Y. Farmers’ willingness to participate in the next-stage grain-for-green project in the Three Gorges Reservoir Area, China [J]. Environmental Management, 2015,56:505-518.

[6] 馮 琳,孫 平,李 丁,等.三峽庫區(qū)肉牛養(yǎng)殖環(huán)境承載力研究[J]. 水生態(tài)學雜志, 2016,37(3):26-33.

Feng L, Sun P, Li D, et al. Environmental carrying capacity of beef cattle in the Three Gorges Reservoir Area [J]. Journal of Hydroecology, 2016,37(3):26-33.

[7] 徐 浩,左新宇,蘭 峰.基于環(huán)境庫茲涅茨曲線方法的三峽庫區(qū)重慶段水質(zhì)演變預測[J]. 水利水電快報, 2019,40(8):40-44,56.

Xu H, Zuo X Y, Lan F. Prediction of water quality evolution in Chongqing section of the Three Gorges Reservoir Area based on environmental Kuznets curve method [J]. Express Water Resources & Hydropower Information, 2019,40(8):40-44,56.

[8] 賴斯蕓,杜鵬飛,陳吉寧.基于單元分析的非點源污染調(diào)查評估方法[J]. 清華大學學報(自然科學版), 2004,(9):1184-1187.

Lai S Y, Du P F, Chen J N. Evaluation of non-point source pollution based on unit analysis [J]. Journal of Tsinghua University (Science and Technology), 2004,(9):1184-1187.

[9] 李玉華.基于SWAT模型的三峽庫區(qū)徑流模擬研究[D]. 重慶:西南大學, 2010.

Li Y H. Runoff simulation using SWAT Model in the Three Gorges Reservoir [D]. Chongqing: Master's Dissertation, Southwestern University, 2010.

[10] 胡文杰,王曉榮,付 甜,等.三峽庫區(qū)蘭陵溪小流域農(nóng)業(yè)面源污染排放特征解析[J]. 環(huán)境污染與防治, 2021,43(5):568-573.

Hu W J, Wang X R, Fu T, et al.Analysis of agricultural non-point source pollution emission characteristics of Lanlingxi small watershed in Three Gorges Reservoir Area [J]. Environmental Pollution & Control, 2021,43(5):568-573.

[11] 鄭永林,王海燕,秦倩倩,等.三峽庫區(qū)筍溪河流域水質(zhì)評價及其空間差異[J]. 水生態(tài)學雜志, 2021,42(4):1-9.

Zheng Y L, Wang H Y, Qin Q Q, et al. Assessment and spatial variation of water quality in the Sunxi River Basin in the Three Gorges Reservoir Area [J]. Journal of Hydroecology, 2021,42(4):1-9.

[12] 蘇青青,劉德富,劉綠波,等.三峽水庫蓄水期支流水體營養(yǎng)鹽來源估算[J]. 中國環(huán)境科學, 2018,38(10):3925-3932.

Su Q Q, Liu D F, Liu L B, et al. Analysis of the nutrient supply of tributaries in the Three Gorges Reservoir during impounding period [J]. China Environmental Science, 2018,38(10):3925-3932.

[13] 胡靜鋒.重慶市農(nóng)業(yè)面源污染測算與空間特征解析[J]. 中國農(nóng)業(yè)資源與區(qū)劃, 2017,38(1):135-144.

Hu J F. Measurement and spatial feature analysis of agricultural non-point source pollution in Chongqing [J]. Chinese Journal of Agricultural Resources and Regional Planning, 2017, 38(1):135-144.

[14] 張廣納,邵景安,王金亮,等.三峽庫區(qū)重慶段農(nóng)村面源污染時空格局演變特征[J]. 自然資源學報, 2015,30(7):1197-1209.

Zhang G N, Shao J A, Wang J L, et al. Spatial and temporal variations of agricultural non-point source pollution in the Three Gorges Reservoir Area of Chongqing [J]. Journal of Natural Resources, 2015,30(7):1197-1209.

[15] 蔡金洲,范先鵬,黃 敏,等.湖北省三峽庫區(qū)農(nóng)業(yè)面源污染解析[J]. 農(nóng)業(yè)環(huán)境科學學報, 2012,31(7):1421-1430.

Cai J Z, Fan X P, Huang M, et al. Sources analysis of agricultural non-point source pollution in the Three Gorges Reservoir Area of Hubei Province, China [J]. Journal of Agro-Environment Science, 2012,31(7):1421-1430.

[16] 于 驥,蒲 實,周 靈.四川省農(nóng)業(yè)面源污染與農(nóng)業(yè)增長的實證分析[J]. 農(nóng)村經(jīng)濟, 2016,(9):56-60.

Yu J, Pu S, Zhou L. Empirical analysis of agricultural surface source pollution and agricultural growth in Sichuan Province [J]. Rural Economy, 2016,(9):56-60.

[17] 閆桂權(quán),何玉成,張曉恒.綠色技術(shù)進步、農(nóng)業(yè)經(jīng)濟增長與污染空間溢出——來自中國農(nóng)業(yè)水資源利用的證據(jù)[J]. 長江流域資源與環(huán)境, 2019,28(12):2921-2935.

Yan G Q, He Y C, Zhang X H. Green technology progress, agricultural economic growth and pollutionspace spillover effect: Evidence of agricultural water utilization process in China [J]. Resources and Environment in the Yangtze Basin, 2019,28(12): 2921-2935.

[18] 彭甲超,肖建忠,李 綱,等.長江經(jīng)濟帶農(nóng)業(yè)廢水面源污染與農(nóng)業(yè)經(jīng)濟增長的脫鉤關(guān)系[J]. 中國環(huán)境科學, 2020,40(6):2770-2784.

Peng J C, Xiao J Z, Li G, et al. Decoupling relationship between agricultural wastewater non-point source pollution and agricultural economic growth in the Yangtze River Economic Belt [J]. China Environmental Science, 2020,40(6):2770-2784.

[19] Grossman G M, Krueger A B. Environmental impacts of a North American Free Trade Agreement [R]. Cambridge, MA: National Bureau of Economic Research, 1991.

[20] 李 競,侯麗朋,唐立娜.基于環(huán)境庫茲涅茨曲線的我國大氣污染防治重點區(qū)域環(huán)境空氣質(zhì)量與經(jīng)濟增長關(guān)系研究 [J]. 生態(tài)學報, 2021,41(22):8845-8859.

Li J, Hou L P, Tang L N. Relationship between air quality and economic growth in key areas ofair pollution control in China based on the Environmental Kuznets Curve [J]. Acta Ecologica Sinica, 2021,41(22):8845-8859.

[21] 王 遠,王春春,趙 靜,等.江蘇省SO2排放環(huán)境庫茨涅茨關(guān)系研究——基于半?yún)?shù)面板數(shù)據(jù)模型分析[J]. 中國環(huán)境科學, 2016, 36(10):3143-3149.

Wang Y, Wang C C, Zhao J, et al. Investigating the environmental Kuznets curve hypothesis for SO2emissions in Jiangsu Province: A semi-parametric panel data analysis [J]. China Environmental Science, 2016,36(10):3143-3149.

[22] 孫大元,楊祁云,張景欣,等.廣東省農(nóng)業(yè)面源污染與農(nóng)業(yè)經(jīng)濟發(fā)展的關(guān)系[J]. 中國人口·資源與環(huán)境, 2016,26(S1):102-105.

Sun D Y, Yang Q Y, Zhang J X, et al.Positive analysis on the relationship between agriculture non-point source pollution and economic growth of Guangdong Province [J]. China Population, Resources and Environment, 2016,26(S1):102–105.

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

Jie C L, Wang J L, Pang Y N. Agricultural growth and fertilizer surface-source pollution in China: Does the environmental Kuznets curve exist [J].Rural Economy, 2018,(11):110–117.

[24] 沈 能,張 斌.農(nóng)業(yè)增長能改善環(huán)境生產(chǎn)率嗎?——有條件“環(huán)境庫茲涅茨曲線”的實證檢驗[J]. 中國農(nóng)村經(jīng)濟, 2015,(7):17-30.

Shen N, Zhang B. Can agricultural growth improve environmental productivity? --An empirical test of the conditional "environmental Kuznets curve"[J].Chinese Rural Economy, 2015,(7):17–30.

[25] Martin Wagner. The carbon Kuznets curve: A cloudy picture emitted by bad econometrics? [J]. Resource and Energy Economics, 2008, 30(3):388–408.

[26] 劉 揚,陳劭鋒,張云芳.中國農(nóng)業(yè)EKC研究:以化肥為例[J]. 中國農(nóng)學通報, 2009,25(16):263-267.

Liu Y, Chen S F, Zhang Y F.Study on Chinese agricultural EKC: Evidence from fertilizer [J].Chinese Agricultural Science Bulletin, 2009,25(16):263–267.

[27] 尚 杰,李 新,鄧雁云.基于EKC的農(nóng)業(yè)經(jīng)濟增長與農(nóng)業(yè)面源污染的關(guān)系分析——以黑龍江省為例[J]. 生態(tài)經(jīng)濟, 2017,33(6):157- 160,166.

Shang J, Li X, Deng Y Y. Analysis of the relationship between agricultural economic growth and agricultural non-point source pollution based on EKC: Taking HeilongjiangProvince as an example [J]. Ecological Economy, 2017,33(6):157-160,166.

[28] 陳 勇,馮永忠,楊改河.農(nóng)業(yè)非點源污染源的環(huán)境庫茲涅茨曲線實證研究——基于陜西省農(nóng)業(yè)投入和廢棄物排放的研究[J]. 干旱地區(qū)農(nóng)業(yè)研究, 2010,28(3):191-198.

ChenY, Feng Y Z, Yang G H.Empirical study on Environmental Kuznets Curve of non-point sourcepollution——A case study of agricultural input and agricultural waste'soutput in Shaanxi Province [J]. Agricultural Research in the Arid Areas, 2010,28(3):191–198.

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

Hou M Y, Yao S B. EKC Retest of fertilizer non-point source pollution emissionunderheterogeneous conditions——grouping based on Panel Threshold Model [J]. Journal of Agrotechnical Economics, 2019,(4): 104-118.

[30] 劉志欣,邵景安,李陽兵.重慶市農(nóng)業(yè)面源污染源的EKC實證分析[J]. 西南師范大學學報(自然科學版), 2015,40(11):94-101.

Liu Z X, Shao J A, Li Y B. The empirical analysis of EKC on agricultural non-point source pollution in Chongqing [J].Journal of Southwest China Normal University (Natural Science Edition), 2015, 40(11):94-101.

[31] 牟新利,祁俊生,黃 宇,等.重慶農(nóng)業(yè)面源污染的環(huán)境庫茲涅茨曲線評價[J]. 貴州農(nóng)業(yè)科學, 2011,39(1):228-230.

Mu X Li, Qi J S, Huang Y, et al.Pollution of agricultural non-point source based on Environmental Kuznets Curve Model in Chongqing [J]. Guizhou Agricultural Sciences, 2011,39(1):228-230.

[32] Zhang T, Ni J P, Xie D T. Assessment of the relationship between rural non-point source pollution and economic development in the Three Gorges Reservoir Area. Environmental Science and Pollution Research [J]. 2016,23(8):8125-8132.

[33] 雷俊華,蘇時鵬,余文夢,等.中國省域化肥面源污染時空格局演變與分組預測[J]. 中國生態(tài)農(nóng)業(yè)學報(中英文), 2020,28(7):1079-1092.

Lei J H, Su S P, Yu W M, et al. Temporal and spatial pattern evolution and grouping prediction of non-point source pollution of chemical fertilizers in China [J]. Chinese Journal of Eco-Agriculture, 2020, 28(7):1079-1092.

[34] 陳玉成,楊志敏,陳慶華,等.基于“壓力-響應”態(tài)勢的重慶市農(nóng)業(yè)面源污染的源解析[J]. 中國農(nóng)業(yè)科學, 2008,41(8):2362-2369.

Chen Y C, Yang Z M, Chen Q H, et al. Source apportionment of agricultural non-point source pollution in Chongqing based on pressure-response system [J]. Scientia Agricultura Sinica, 2008,41(8): 2362-2369.

[35] 國家環(huán)境保護總局.全國規(guī)?；笄蒺B(yǎng)殖業(yè)污染情況調(diào)查及防治對策[M]. 北京:中國環(huán)境科學出版社, 2002,25:77-78.

State Environmental Protection Administration. Survey of pollution in the national large-scale livestock and poultry farming industry and countermeasures for prevention and control [M]. Beijing: China Environmental Science Press, 2002,25:77-78.

[36] 葛繼紅.江蘇省農(nóng)業(yè)面源污染及治理的經(jīng)濟學研究[D]. 南京:南京農(nóng)業(yè)大學, 2011.

Ge J H. Economics of agricultural surface source pollution and management in Jiangsu Province [D]. Nanjing: Nanjing Agricultural University, 2011.

[37] 熊昭昭,王書月,童 雨,等.江西省農(nóng)業(yè)面源污染時空特征及污染風險分析[J]. 農(nóng)業(yè)環(huán)境科學學報, 2018,37(12):2821-2828.

Xiong Z Z, Wang S Y, TONG Y, et al. Analysis on spatio-temporal characteristics of agricultural non-point sources pollution and pollution risk in Jiangxi Province, China [J]. Journal of Agro- Environment Science, 2018,37(12):2821-2828.

[38] 許 策,李 曄,束繼年等.漢江流域荊門段面源污染負荷時空分布與污染現(xiàn)狀評價[J]. 水土保持通報, 2017,37(4):63-68.

XU C, Li Y, SHU J N, et al. Spatial distribution and evaluation of non-point pollution source in Jingmen Region of Hanjiang Watershed [J]. Bulletin of Soil and Water Conservation, 2017,37(4):63-68.

[39] 孫傳旺,羅 源,姚 昕.交通基礎(chǔ)設(shè)施與城市空氣污染——來自中國的經(jīng)驗證據(jù) [J]. 經(jīng)濟研究, 2019,54(8):136-151.

Sun C W, Luo Y, Yao X.Transportation Infrastructure and urban air pollution - Empirical evidence from China [J]. Economic Research Journal, 2019,54(8):136-151.

Analysis on spatial-temporal characteristics and Environmental Kuznets Curve of non-point source pollution in the Three Gorges Reservoir Area.

FENG Lin1*, ZHANG Wan-ting1, ZHANG Jun-ke1, MAO Xin-hui2, YU Xin3, YE Cheng-song3, CHANG Tian-yang1

(1.School of Environment & Natural Resources, Renmin University of China, Beijing 100872, China；2.Chinese Academy of Natural Resources Economics,Beijing 101149, China；3.College of the Environment & Ecology, Xiamen University, Xiamen361102, China)., 2022,42(7)：3325～3333

Total Nitrogen(TN) and total phosphorous(TP) emissions of agricultural non-pointsource pollution in 19 regions of the Three Gorges Reservoir Area(TGRA) were estimated from 2008 to 2018. Their spatial-temporal variation characteristics were further analyzed. According to Environmental Kuznets Curve (EKC) theory, a regression model based on panel data was constructed to explore the relationship between agricultural economy development and non-point source pollution emission intensity in TGRA. Our study found that TN emissions of agriculture decreased with fluctuation, and TP increased with fluctuation in TGRA. The average annual emissions of TN and TP in each region ranged from 374～6046 tons and 105～1267 tonsrespectively. Agricultural fertilizer and livestock breeding units contributed to more than 80% of total TN and TP. There were significant "inverted U-shape" EKC relationships to TN emission intensity, TN/TP emission intensity in livestock & poultry breeding unit and rural living unit, indicating that the turning points have been showed in these curves. There were significant "linear" EKC relationships to TP emission intensity in agricultural fertilizer unit, TN/TP emission intensity in aquaculture unit and farmland solid-waste unit, indicating a simultaneous growth with economy. The research finally proposed that agricultural fertilizer unit should upgrade the pollution prevention capacity, so as to improve rural living environment and reduce regional TN and TP emission.

agricultural non-point source pollution；unit analysis method；Environmental Kuznets Curve (EKC)；Three Gorges reservoir area(TGRA)

X501

A

1000-6923(2022)07-3325-09

馮 琳(1979-),女,江蘇泰州人,副教授,博士,主要研究方向為環(huán)境規(guī)劃與管理.發(fā)表論文30余篇.

2021-12-17

中國人民大學2022年度“中央高校建設(shè)世界一流大學(學科)和特色發(fā)展引導專項資金”

* 責任作者, 副教授, feng.lin@163.com