摘 要:本文基于2004—2023年大汶河流域內(nèi)各地市的面板數(shù)據(jù),以該流域生態(tài)補(bǔ)償實(shí)際情況為研究對(duì)象,運(yùn)用雙重差分模型探究流域生態(tài)補(bǔ)償機(jī)制的環(huán)境、經(jīng)濟(jì)效應(yīng)。研究發(fā)現(xiàn),大汶河生態(tài)補(bǔ)償對(duì)該流域的環(huán)境狀況改善和整體經(jīng)濟(jì)發(fā)展都存在正面影響,但其影響強(qiáng)度隨著政策實(shí)施周期推進(jìn)逐漸降低,應(yīng)當(dāng)延續(xù)并優(yōu)化此類流域生態(tài)補(bǔ)償政策,推動(dòng)大汶河流域可持續(xù)發(fā)展。
關(guān)鍵詞:流域生態(tài)補(bǔ)償;經(jīng)濟(jì)效應(yīng);環(huán)境效益;雙重差分模型
中圖分類號(hào):TP 79 文獻(xiàn)標(biāo)志碼:A
隨著經(jīng)濟(jì)發(fā)展與生態(tài)保護(hù)之間的矛盾凸顯,以促進(jìn)流域高質(zhì)量發(fā)展、上下游流域協(xié)同發(fā)展的生態(tài)補(bǔ)償制度逐漸清晰。黨中央和國(guó)務(wù)院對(duì)流域生態(tài)保護(hù)的關(guān)注度日益增加,并將健全我國(guó)生態(tài)保護(hù)補(bǔ)償機(jī)制作為重要任務(wù)。狹義上,流域生態(tài)補(bǔ)償就是由下游生態(tài)受益區(qū)居民對(duì)流域內(nèi)上游經(jīng)濟(jì)主體的生態(tài)保護(hù)行為進(jìn)行激勵(lì)和補(bǔ)貼,以達(dá)到生態(tài)保護(hù)的目的,我國(guó)的相關(guān)學(xué)術(shù)探索工作起源于20世紀(jì)80年代。從生態(tài)補(bǔ)償?shù)膬?nèi)涵來(lái)說(shuō),門寶輝等[1]認(rèn)為流域生態(tài)補(bǔ)償是指對(duì)因人類活動(dòng)而產(chǎn)生的流域生態(tài)環(huán)境負(fù)外部效應(yīng)進(jìn)行補(bǔ)償?shù)男袨?。牛坤玉等[2]提出生態(tài)補(bǔ)償也就是生態(tài)系統(tǒng)付費(fèi)服務(wù)。是人類從自然界獲得的福利。從生態(tài)補(bǔ)償標(biāo)準(zhǔn)來(lái)看,羅萬(wàn)云等[3]提出在生態(tài)公平視域下理清權(quán)責(zé)關(guān)系,確定流域補(bǔ)償標(biāo)準(zhǔn)。鄧夢(mèng)華等[4]基于學(xué)習(xí)模型XGBoost-SHAP,將流域居民補(bǔ)償意愿納入其中。綜上,學(xué)者們對(duì)流域生態(tài)補(bǔ)償進(jìn)行規(guī)范化研究,但該研究還有待進(jìn)一步深入。本文以大汶河流域的生態(tài)補(bǔ)償機(jī)制為研究對(duì)象,為流域協(xié)調(diào)發(fā)展提供依據(jù)。
1 模型構(gòu)建與變量選取
本文采用雙重差分法(DID)對(duì)大汶河流域生態(tài)補(bǔ)償政策進(jìn)行評(píng)估。雙重差分法(DID)雙重差分法最早由ASHENFELTER引入經(jīng)濟(jì)學(xué)研究中,是基于各地市面板數(shù)據(jù)進(jìn)行機(jī)制評(píng)估的方法,本文建立的模型如公式(1)所示。
Yit=β0+β1didit+β2x+δi+θt+εit (1)
式中:Yit為表征環(huán)境效應(yīng)、經(jīng)濟(jì)效應(yīng)的被解釋變量;didit為核心解釋變量;X為表征環(huán)境效應(yīng)、經(jīng)濟(jì)效應(yīng)的控制變量組;δi為個(gè)體效應(yīng);θt為時(shí)間效應(yīng);εit為隨機(jī)干擾項(xiàng);β1為所求政策效應(yīng)值,如果系數(shù)符號(hào)為正,就說(shuō)明生態(tài)補(bǔ)償機(jī)制將改善環(huán)境,促進(jìn)經(jīng)濟(jì)發(fā)展,如果系數(shù)符號(hào)為負(fù),就說(shuō)明生態(tài)補(bǔ)償機(jī)制將不利于環(huán)境改善和經(jīng)濟(jì)發(fā)展;β2為控制變量組系數(shù)。
1.1 被解釋變量
廢水排放量:減少流域內(nèi)污染物排放,促進(jìn)生態(tài)系統(tǒng)良性發(fā)展是流域生態(tài)補(bǔ)償?shù)闹匾康?。本文通過(guò)大汶河流域各地市的年工業(yè)廢水排放總量來(lái)衡量[4]。
地區(qū)生產(chǎn)總值:GDP是國(guó)民經(jīng)濟(jì)核算的核心指標(biāo),本文通過(guò)引用大汶河流域地市生產(chǎn)總值來(lái)反映流域內(nèi)生態(tài)補(bǔ)償政策實(shí)施后帶來(lái)的經(jīng)濟(jì)效果影響。
1.2 核心解釋變量
本文在模型中引入地區(qū)與時(shí)間虛擬變量的交互項(xiàng),來(lái)檢驗(yàn)試驗(yàn)組與對(duì)照組由大汶河生態(tài)補(bǔ)償機(jī)制而引起的流域內(nèi)經(jīng)濟(jì)和生態(tài)環(huán)境的變化情況。地區(qū)虛擬變量(treat)試驗(yàn)組取1,對(duì)照組取0;時(shí)間虛擬變量將大汶河流域機(jī)制實(shí)施的時(shí)間(t≥2008)取1,機(jī)制實(shí)施之前的年份(t<2008)取0。
1.3 控制變量
產(chǎn)業(yè)結(jié)構(gòu):區(qū)域產(chǎn)業(yè)結(jié)構(gòu)的調(diào)整是區(qū)域內(nèi)經(jīng)濟(jì)發(fā)展水平和生態(tài)環(huán)境變化的重要體現(xiàn)指標(biāo)之一。城市現(xiàn)代化的發(fā)展過(guò)程當(dāng)中,第二、三產(chǎn)業(yè)的發(fā)展情況對(duì)經(jīng)濟(jì)進(jìn)步和環(huán)境變化起到重要影響,因此本文取流域內(nèi)各地市的第三產(chǎn)業(yè)和第二產(chǎn)業(yè)增加值比重的對(duì)數(shù)來(lái)表示。
科技支出:創(chuàng)新會(huì)帶來(lái)科技水平的進(jìn)步,是發(fā)展的源泉和第一動(dòng)力,而科技進(jìn)步對(duì)轉(zhuǎn)換城市經(jīng)濟(jì)增長(zhǎng)動(dòng)力、提高各產(chǎn)業(yè)水資源利用率和提升環(huán)境污染治理能力等方面均能起到積極作用,促進(jìn)生態(tài)環(huán)境與經(jīng)濟(jì)建設(shè)可持續(xù)發(fā)展。因此,本文以各地區(qū)年政府科技支出的對(duì)數(shù)來(lái)表示。
城鎮(zhèn)化率:改革開放以來(lái),我國(guó)的城鎮(zhèn)化率大大提高,包括大汶河流域在內(nèi),為滿足日益增長(zhǎng)的人口對(duì)生態(tài)資源的需要,人類社會(huì)的發(fā)展在一定程度上對(duì)該流域生態(tài)環(huán)境施加了更大的壓力。因此本文以各地市的城鎮(zhèn)常住人口占地區(qū)常住人口的比例來(lái)表示。
教育支出:教育支出是影響該地區(qū)市民的受教育水平的重要指標(biāo),受教育程度越高,越能夠?qū)斫馍鷳B(tài)保護(hù)的重要性并積極配合政府相關(guān)的環(huán)境政策和環(huán)保號(hào)召。因此本文取教育和文化財(cái)政支出的對(duì)數(shù)來(lái)表示。
生態(tài)支出:生態(tài)支出是政府購(gòu)買的重要組成部分,主要用于政府生態(tài)修復(fù)工程的投資、相關(guān)生態(tài)產(chǎn)品與服務(wù)的購(gòu)買支出、有關(guān)生態(tài)金融衍生品和政策的補(bǔ)貼與補(bǔ)助等方面。因此本文取各地區(qū)歷年的農(nóng)業(yè)、林業(yè)和水資源這三類環(huán)境保護(hù)支出總和的對(duì)數(shù)來(lái)表示。
2 大汶河流域生態(tài)補(bǔ)償政策效應(yīng)評(píng)估
2.1 描述性統(tǒng)計(jì)
表1為大汶河流域生態(tài)補(bǔ)償機(jī)制的環(huán)境效應(yīng)和經(jīng)濟(jì)效應(yīng)的相關(guān)變量的描述性統(tǒng)計(jì)。樣本總數(shù)為270,其中2004—2023年該地區(qū)工業(yè)廢水排放量平均數(shù)為7543萬(wàn)t,標(biāo)準(zhǔn)差為5755萬(wàn)t,最小值1210萬(wàn)t,中位數(shù)為5772萬(wàn)t,最大值為21212萬(wàn)t,說(shuō)明該流域工業(yè)廢水排放量分布不均衡,也表明少部分地區(qū)工業(yè)較發(fā)達(dá),大部分地區(qū)工業(yè)化水平較落后。2004—2023年該地區(qū)生產(chǎn)總值平均數(shù)為2828億元,標(biāo)準(zhǔn)差為3259億元,最小為值223.9億元,中位數(shù)為2233億元,最大值為22901億元,說(shuō)明經(jīng)濟(jì)水平發(fā)展不均衡,生產(chǎn)總值最小地區(qū)只占生產(chǎn)總值最大地區(qū)的10%左右。2004—2023年該地區(qū)產(chǎn)業(yè)結(jié)構(gòu)中第三產(chǎn)業(yè)和第二產(chǎn)業(yè)增加值的比重平均數(shù)為0.810,標(biāo)準(zhǔn)差為0.330,最小值為0.380,最大值為1.680,說(shuō)明該流域內(nèi)不同地市在產(chǎn)業(yè)發(fā)展方向有所側(cè)重,并未全部發(fā)展工業(yè)或服務(wù)業(yè),產(chǎn)業(yè)發(fā)展結(jié)構(gòu)相對(duì)合理。
2.2 平行趨勢(shì)檢驗(yàn)
雙重差分法的應(yīng)用前提是試驗(yàn)組與對(duì)照組在機(jī)制的實(shí)施前具有相同的發(fā)展趨勢(shì)[5]。因此,在機(jī)制效應(yīng)評(píng)估前,本文構(gòu)建模型,如公式(2)所示。
(2)
式中:Yit為表征環(huán)境效應(yīng)、經(jīng)濟(jì)效應(yīng)的被解釋變量,分別代入廢水排放量和地區(qū)生產(chǎn)總值進(jìn)行檢驗(yàn);β0為生態(tài)補(bǔ)償系數(shù);δi為個(gè)體效應(yīng);θt為時(shí)間效應(yīng);εit為隨機(jī)干擾項(xiàng);縱坐標(biāo)為流域污水排放量和地區(qū)生產(chǎn)總值回歸系數(shù);prej和lastl分別為政策實(shí)施前、后的時(shí)間虛擬變量。
如果大汶河流域生態(tài)補(bǔ)償機(jī)制的實(shí)施年份為2008年,那么pre2代表政策實(shí)施的兩年前,即2006年;同理,post2為政策實(shí)施后第二年,即2010年。結(jié)果如圖1、圖2所示,在政策執(zhí)行前,對(duì)生態(tài)補(bǔ)償機(jī)制系數(shù)β的置信區(qū)間均為零,這說(shuō)明在政策執(zhí)行前,資料處理試驗(yàn)組與對(duì)照組的廢水排放量和地區(qū)生產(chǎn)總值并未出現(xiàn)明顯區(qū)別,因此可以將其視為平行。在2008年流域相關(guān)保護(hù)政策正式執(zhí)行后,政策系數(shù)β的置信區(qū)間在往后連續(xù)的四期內(nèi)都不為零,說(shuō)明政府制定并執(zhí)行的生態(tài)政策已經(jīng)對(duì)現(xiàn)實(shí)情況中的處理廢水排放量和地區(qū)生產(chǎn)總值等相關(guān)研究要素形成了一定限制與影響,該結(jié)果滿足平行趨勢(shì)假設(shè)成立所需要求,因此可以運(yùn)用雙重差分法展開進(jìn)一步的分析研究。
2.3 雙重差分檢驗(yàn)結(jié)果
本文借助Stata17軟件進(jìn)行基本回歸分析。雙重查分檢驗(yàn)結(jié)果見表2。其中(1)、(3)列為廢水排放量環(huán)境效應(yīng)評(píng)估模型,(2)、(4)列為地區(qū)生產(chǎn)總值為經(jīng)濟(jì)效應(yīng)評(píng)估模型。并且(1)列為對(duì)照組,(3)列為試驗(yàn)組,試驗(yàn)組相較對(duì)照組增加了上文所描述與統(tǒng)計(jì)的若干解釋變量。從該模型檢驗(yàn)的結(jié)果可知,生態(tài)補(bǔ)償機(jī)制的環(huán)境效應(yīng)系數(shù)為-4.273,在1%的水平上顯著為負(fù),表明其在減少該流域工業(yè)廢水排放量上起到積極作用,即減少了該流域工業(yè)廢水排放總量,顯著改善了流域水資源的生態(tài)環(huán)境。同時(shí)由表2可得,生態(tài)補(bǔ)償機(jī)制經(jīng)濟(jì)效應(yīng)系數(shù)為4.323,在1%的水平上顯著為正,表明生態(tài)補(bǔ)償政策對(duì)流域的經(jīng)濟(jì)發(fā)展起到了明顯作用,即制定實(shí)施相關(guān)政策顯推動(dòng)了GDP等經(jīng)濟(jì)指標(biāo)增長(zhǎng),提高了該區(qū)域的整體經(jīng)濟(jì)發(fā)展水平。
2.4 穩(wěn)健性檢驗(yàn)
為進(jìn)一步驗(yàn)證結(jié)果的準(zhǔn)確性,本文采用滯后一期實(shí)施有關(guān)流域生態(tài)政策的檢驗(yàn)方法對(duì)上述結(jié)果進(jìn)行穩(wěn)健性檢驗(yàn)?;貧w結(jié)果見表3,其中(1)列代表廢水排放量,(2)列代表生產(chǎn)總值。從表3中可得,就業(yè)率在1%的水平上對(duì)廢水排放量指標(biāo)顯著為負(fù),對(duì)地區(qū)生產(chǎn)總值指標(biāo)顯著為正,說(shuō)明該政策的實(shí)施通過(guò)在生態(tài)保護(hù)行業(yè)吸納就業(yè)人口,達(dá)到了明顯降低了污水排放量和提高經(jīng)濟(jì)社會(huì)發(fā)展的目的,同時(shí)該流域生態(tài)政策的環(huán)境機(jī)制效應(yīng)系數(shù)為-2.403,仍在1%的水平上顯著為負(fù),即延后實(shí)施一期比前一期工業(yè)污水排放量明顯下降,經(jīng)濟(jì)機(jī)制效應(yīng)系數(shù)為2.126,在5%的水平上顯著為正,即延后實(shí)施一期比前一期經(jīng)濟(jì)發(fā)展水平明顯上升,與本文雙重差分法的檢驗(yàn)結(jié)果保持一致。
3 結(jié)語(yǔ)
本文經(jīng)過(guò)相關(guān)研究結(jié)論如下:大汶河流域生態(tài)補(bǔ)償政策實(shí)施后起到了改善當(dāng)?shù)厣鷳B(tài)環(huán)境的作用,并且?guī)?dòng)了當(dāng)?shù)亟?jīng)濟(jì)的發(fā)展;但是大汶河流域生態(tài)補(bǔ)償機(jī)制的環(huán)境效應(yīng)隨時(shí)間發(fā)展其相關(guān)性逐漸減弱,最明顯的就是政策實(shí)施后的第三年效果最顯著。
大汶河流域的生態(tài)補(bǔ)償政策本身和其他方面仍存在不足。1)大汶河流域生態(tài)補(bǔ)償?shù)膮⑴c主體基本是政府,這在一定程度上帶來(lái)了轉(zhuǎn)移支付的財(cái)政壓力和政府部門為實(shí)現(xiàn)政策目標(biāo)而進(jìn)行人員配置等環(huán)節(jié)的執(zhí)行成本,在后疫情時(shí)代對(duì)財(cái)政相對(duì)緊張的各級(jí)政府來(lái)說(shuō)是很難繼續(xù)保證政策實(shí)施成效的。2)大汶河流域的經(jīng)濟(jì)效益提升更多來(lái)自社會(huì)生產(chǎn)力的發(fā)展和產(chǎn)業(yè)升級(jí),對(duì)生態(tài)環(huán)境本身經(jīng)濟(jì)價(jià)值的探索還停留在初級(jí)階段,沒(méi)有很好地發(fā)揮大汶河流域生態(tài)的經(jīng)濟(jì)效應(yīng)。3)大汶河生態(tài)補(bǔ)償機(jī)制的手段不夠豐富,基本依靠較直接的政府轉(zhuǎn)移支付和政府投資來(lái)推動(dòng)生態(tài)修復(fù)和拉動(dòng)經(jīng)濟(jì)發(fā)展,在更復(fù)雜的產(chǎn)權(quán)交易方面缺少探索和實(shí)踐,很難形成由市場(chǎng)為主導(dǎo)驅(qū)動(dòng)的生態(tài)產(chǎn)業(yè)鏈。
參考文獻(xiàn)
[1]門寶輝,劉菁蘋.基于最小數(shù)據(jù)法的黃河中游延安市生態(tài)補(bǔ)償標(biāo)準(zhǔn)研究[J].華北水利水電大學(xué)學(xué)報(bào)(自然科學(xué)版),2022,43(5):101-108.
[2]牛坤玉,劉靜,郭靜利,等.農(nóng)業(yè)生態(tài)補(bǔ)償:內(nèi)涵、要素特征與政策創(chuàng)設(shè)[J].中國(guó)農(nóng)業(yè)資源與區(qū)劃,2024(2):1-10.
[3]羅萬(wàn)云,周楊,王小娟.生態(tài)公平視域下額爾齊斯河流域生態(tài)補(bǔ)償標(biāo)準(zhǔn)及空間選擇研究[J].生態(tài)學(xué)報(bào),2024,(21):1-16.
[4]鄧夢(mèng)華,張?zhí)焓妫愜婏w.基于XGBoost-SHAP模型的太湖流域居民生態(tài)補(bǔ)償支付意愿影響因素研究[J].水利經(jīng)濟(jì),2024,42(2):44-50.
[5]于冰,史穎.新安江生態(tài)補(bǔ)償環(huán)境經(jīng)濟(jì)效應(yīng)及影響機(jī)制研究[J].生態(tài)經(jīng)濟(jì),2023,39(7):157-164.