高齊圣 董志強(qiáng)

摘? ? 要：探究新型城鎮(zhèn)化如何影響農(nóng)業(yè)面源污染，對(duì)推進(jìn)農(nóng)業(yè)農(nóng)村綠色發(fā)展，全面推進(jìn)鄉(xiāng)村振興戰(zhàn)略具有重要意義。利用2007—2019年我國159個(gè)地級(jí)城市的面板數(shù)據(jù)，借助多期雙重差分方法，可評(píng)估新型城鎮(zhèn)化試點(diǎn)政策對(duì)農(nóng)業(yè)面源污染的影響效應(yīng)。研究發(fā)現(xiàn)，新型城鎮(zhèn)化試點(diǎn)政策能夠有效減少農(nóng)業(yè)面源污染排放。影響機(jī)制分析表明，新型城鎮(zhèn)化試點(diǎn)政策通過加強(qiáng)環(huán)境規(guī)制、促進(jìn)農(nóng)業(yè)綠色技術(shù)進(jìn)步來緩解農(nóng)業(yè)面源污染。沖擊分析表明，我國東、中、西部地區(qū)以及大、中、小規(guī)模城市對(duì)新型城鎮(zhèn)化試點(diǎn)的沖擊表現(xiàn)出較大差異。目前實(shí)施的新型城鎮(zhèn)化試點(diǎn)政策能有效減少東部城市、大規(guī)模城市的農(nóng)業(yè)面源污染排放，但對(duì)中西部城市、中小規(guī)模城市的農(nóng)業(yè)面源污染排放不構(gòu)成顯著影響。

關(guān)鍵詞：農(nóng)業(yè)面源污染；新型城鎮(zhèn)化；準(zhǔn)自然實(shí)驗(yàn)；多期DID；內(nèi)在機(jī)制

中圖分類號(hào)：F299.2；F323；X501 文章標(biāo)志碼：A 文章編號(hào)：1006-6152（2023）04-0005-09

DOI：10.16388/j.cnki.cn42-1843/c.2023.04.001

一、問題的提出

“綠水青山就是金山銀山”，良好的生態(tài)環(huán)境既是自然財(cái)富，也是經(jīng)濟(jì)財(cái)富。2022年中央一號(hào)文件指出，加強(qiáng)農(nóng)業(yè)面源污染綜合治理，對(duì)推進(jìn)鄉(xiāng)村生態(tài)振興、農(nóng)業(yè)農(nóng)村綠色發(fā)展具有重要意義?；释度胧寝r(nóng)業(yè)面源污染的重要來源之一，化肥的過量施用使得耕地質(zhì)量不斷下降、地下水污染日益嚴(yán)重［1］，給農(nóng)村生態(tài)環(huán)境帶來了巨大影響［2］。改革開放以來，我國城鎮(zhèn)化發(fā)展迅速，常住人口城鎮(zhèn)化率由1978年的17.9%增長(zhǎng)至2021年的63.9%①。然而，城鎮(zhèn)化推進(jìn)的盲目性，使得農(nóng)業(yè)面源污染現(xiàn)象日益嚴(yán)重［3］。在此背景下，更加強(qiáng)調(diào)“生態(tài)文明”的新型城鎮(zhèn)化概念被提出。新型城鎮(zhèn)化的本質(zhì)是“以人為本”［4］，核心在于不以犧牲農(nóng)業(yè)和糧食、生態(tài)和環(huán)境為代價(jià)，力圖實(shí)現(xiàn)由片面追求經(jīng)濟(jì)增長(zhǎng)的非生態(tài)城鎮(zhèn)化向人與自然和諧共生的生態(tài)城鎮(zhèn)化轉(zhuǎn)型。2014年，《國家新型城鎮(zhèn)化規(guī)劃（2014—2020年）》正式發(fā)布，文件對(duì)我國新型城鎮(zhèn)化發(fā)展進(jìn)行了總體部署。在新型城鎮(zhèn)化試點(diǎn)政策的實(shí)施中，國家共分三批次發(fā)布了試點(diǎn)地區(qū)名單，試點(diǎn)地區(qū)包含省、市、縣、鎮(zhèn)四個(gè)不同的實(shí)施層級(jí)，涵蓋了我國東、中、西、東北四個(gè)不同區(qū)域。2015年2月，國家發(fā)改委印發(fā)《國家新型城鎮(zhèn)化綜合試點(diǎn)方案》，將江蘇、安徽兩省和寧波等62個(gè)城市（鎮(zhèn)）列為國家新型城鎮(zhèn)化綜合試點(diǎn)地區(qū)；2015年11月，國家另外選取了73個(gè)城市（鎮(zhèn)）作為第二批試點(diǎn)地區(qū)；2016年12月，國家再另外選取了111個(gè)城市（鎮(zhèn)）作為第三批試點(diǎn)地區(qū)。在新的時(shí)代背景下，探究新型城鎮(zhèn)化試點(diǎn)對(duì)農(nóng)業(yè)面源污染的直接效應(yīng)和間接效應(yīng)，對(duì)實(shí)現(xiàn)農(nóng)村生態(tài)環(huán)境改善、推進(jìn)農(nóng)業(yè)農(nóng)村綠色發(fā)展具有重要的現(xiàn)實(shí)意義。

關(guān)于城鎮(zhèn)化與農(nóng)業(yè)面源污染之間的關(guān)系，現(xiàn)有研究主要分為三類。第一類研究認(rèn)為城鎮(zhèn)化會(huì)增加農(nóng)業(yè)面源污染排放。城鎮(zhèn)化的發(fā)展一方面會(huì)吸引大量?jī)?yōu)質(zhì)勞動(dòng)力非農(nóng)轉(zhuǎn)移，農(nóng)村老齡人口的增加導(dǎo)致了化肥的過量施用［5］；另一方面會(huì)占用大量農(nóng)地，而耕地的減少導(dǎo)致農(nóng)戶增加化肥投入以擴(kuò)大產(chǎn)出［6］。第二類研究認(rèn)為城鎮(zhèn)化的發(fā)展會(huì)緩解農(nóng)業(yè)面源污染問題。城鎮(zhèn)化的發(fā)展一方面促進(jìn)了綠色技術(shù)進(jìn)步［7］，另一方面增加了綠色農(nóng)產(chǎn)品的需求量［8］，進(jìn)而減少農(nóng)業(yè)面源污染。第三類研究認(rèn)為，城鎮(zhèn)化與農(nóng)業(yè)面源污染之間存在“倒U型”曲線關(guān)系或具有門檻特征。薛蕾等研究發(fā)現(xiàn)，農(nóng)業(yè)面源污染因城鎮(zhèn)化水平的提高會(huì)呈現(xiàn)出“倒Ｕ型”趨勢(shì)［9］；欒健等基于門檻效應(yīng)模型得出，城鎮(zhèn)化對(duì)化肥面源污染的影響會(huì)因人均收入水平的高低產(chǎn)生完全相反的作用效果［10］。

綜上可見，學(xué)者們對(duì)城鎮(zhèn)化和農(nóng)業(yè)面源污染進(jìn)行了大量有價(jià)值的研究，但鮮有學(xué)者從新型城鎮(zhèn)化角度探究它對(duì)農(nóng)業(yè)面源污染的影響。為此，本文將考察新型城鎮(zhèn)化對(duì)農(nóng)業(yè)面源污染的影響，并探索其作用機(jī)制。

二、理論分析與研究假說

（一）新型城鎮(zhèn)化對(duì)農(nóng)業(yè)面源污染的內(nèi)在影響機(jī)制分析

首先，新型城鎮(zhèn)化建設(shè)將促進(jìn)環(huán)境規(guī)制強(qiáng)度的提升，進(jìn)而減少農(nóng)業(yè)面源污染排放。許多研究認(rèn)為，環(huán)境規(guī)制可以有效降低環(huán)境污染［11-12］。一方面，地方政府通過對(duì)化肥、農(nóng)藥等污染源提岀嚴(yán)格治理措施，從源頭上治理農(nóng)業(yè)面源污染；另一方面，地方政府采用綠色農(nóng)業(yè)補(bǔ)貼的形式，比如“以獎(jiǎng)促治”“以獎(jiǎng)代補(bǔ)”等，鼓勵(lì)農(nóng)民向綠色友好型生產(chǎn)方式轉(zhuǎn)變，控制農(nóng)業(yè)面源污染的產(chǎn)生。新型城鎮(zhèn)化試點(diǎn)的環(huán)境規(guī)制效應(yīng)具體可能表現(xiàn)在政府通過明確污染物排放的具體目標(biāo)、提高污染物排放標(biāo)準(zhǔn)以減少污染物的排放［13］。

其次，新型城鎮(zhèn)化建設(shè)將促進(jìn)農(nóng)業(yè)綠色技術(shù)進(jìn)步，進(jìn)而減少農(nóng)業(yè)面源污染排放。許多研究認(rèn)為，技術(shù)進(jìn)步在農(nóng)業(yè)面源污染治理中能夠發(fā)揮重要作用［7］［14］。新型城鎮(zhèn)化沖擊的技術(shù)效應(yīng)可能表現(xiàn)在新的技術(shù)進(jìn)步促進(jìn)了節(jié)能減排。目前我國面源污染控制技術(shù)在效率、規(guī)模方面存在嚴(yán)重不足［15］。但是，新型城鎮(zhèn)化建設(shè)一方面會(huì)促使企業(yè)開發(fā)出更加清潔、利用效率更高的生產(chǎn)技術(shù)［16］；另一方面，技術(shù)進(jìn)步有利于污染防護(hù)和治理水平的提高［13］，可以更好地實(shí)現(xiàn)農(nóng)業(yè)面源污染減排的目標(biāo)。

（二）地區(qū)農(nóng)業(yè)面源污染對(duì)新型城鎮(zhèn)化沖擊反應(yīng)分析

首先，對(duì)待共同的新型城鎮(zhèn)化試點(diǎn)沖擊，各個(gè)地區(qū)或城市可能會(huì)因地理情況呈現(xiàn)不同反應(yīng)而凸顯差異。不同區(qū)域城市在基礎(chǔ)設(shè)施水平、經(jīng)濟(jì)環(huán)境、人才供給狀況等方面往往存在較大差異。一般認(rèn)為，東部城市在技術(shù)進(jìn)步、人才集聚等方面相比于中西部城市具有明顯的優(yōu)勢(shì)［17］，能綜合利用不同類型的環(huán)境規(guī)制推動(dòng)清潔生產(chǎn)技術(shù)的革新，試點(diǎn)政策可能更有助于減少東部城市的農(nóng)業(yè)面源污染排放。其次，新型城鎮(zhèn)化試點(diǎn)的效果可能會(huì)因城市規(guī)模的不同而有所差異。由于大規(guī)模城市的市場(chǎng)需求多樣、科研基礎(chǔ)較好，對(duì)創(chuàng)新人才、高新企業(yè)往往有更強(qiáng)的吸引力，因而相較于中小規(guī)模城市能更高效地配置大量資源，促進(jìn)企業(yè)進(jìn)行綠色技術(shù)創(chuàng)新［18］，因此，新型城鎮(zhèn)化試點(diǎn)政策可能更有助于緩解大規(guī)模城市的農(nóng)業(yè)面源污染。

以上述理論分析為基礎(chǔ)，本文提出以下3個(gè)研究假設(shè)：

假設(shè)1：新型城鎮(zhèn)化試點(diǎn)政策能夠通過提升環(huán)境規(guī)制強(qiáng)度減少農(nóng)業(yè)面源污染排放。

假設(shè)2：新型城鎮(zhèn)化試點(diǎn)政策能夠促進(jìn)農(nóng)業(yè)綠色技術(shù)進(jìn)步，進(jìn)而緩解農(nóng)業(yè)面源污染問題。

假設(shè)3：新型城鎮(zhèn)化試點(diǎn)的農(nóng)業(yè)面源污染減排效應(yīng)會(huì)因城市所在地區(qū)以及城市規(guī)模的不同存在差異。

三、研究設(shè)計(jì)

（一）模型設(shè)計(jì)

本文研究地級(jí)市層面的新型城鎮(zhèn)化試點(diǎn)政策對(duì)農(nóng)業(yè)面源污染的影響效應(yīng)，根據(jù)國家發(fā)展和改革委員會(huì)公布的三批新型城鎮(zhèn)化試點(diǎn)城市名單，將試點(diǎn)城市定義為實(shí)驗(yàn)組，非試點(diǎn)城市定義為控制組。其中，對(duì)實(shí)驗(yàn)組和控制組的選取做了如下處理：剔除了只將地級(jí)市內(nèi)的某個(gè)縣級(jí)市或區(qū)作為試點(diǎn)的地級(jí)市；將數(shù)據(jù)缺失嚴(yán)重的地級(jí)市在樣本中剔除。最終選擇了75個(gè)地級(jí)城市進(jìn)入實(shí)驗(yàn)組，84個(gè)地級(jí)城市進(jìn)入控制組。考慮到政策實(shí)施需要時(shí)間，本文將2015年、2016年和2017年作為試點(diǎn)政策實(shí)施的時(shí)間節(jié)點(diǎn)。

傳統(tǒng)DID方法通過構(gòu)建政策分組虛擬變量和政策實(shí)施時(shí)間虛擬變量的交互項(xiàng)，來估計(jì)政策實(shí)施的凈效應(yīng)，這種方法只適用于政策是一次實(shí)施的，而新型城鎮(zhèn)化政策的實(shí)施是先試點(diǎn)后推廣，因此本文將現(xiàn)行的三批新型城鎮(zhèn)化試點(diǎn)政策作為一項(xiàng)準(zhǔn)自然實(shí)驗(yàn)，采用多期DID方法［19-20］，通過比較實(shí)驗(yàn)組和控制組在政策實(shí)施前后的差異，評(píng)估新型城鎮(zhèn)化試點(diǎn)政策對(duì)農(nóng)業(yè)面源污染的影響效應(yīng)。具體模型設(shè)定如下：

[TEit=β0+β1Dit+β2Xit+μi+γt+εit] （1）

這個(gè)模型中，[TEit]為被解釋變量；[X]表示控制變量；[μi]和[γt]分別代表控制城市和時(shí)間的固定效應(yīng)；[εit]為隨機(jī)誤差項(xiàng)，[i]表示城市，[t]為時(shí)間。[Dit]表示個(gè)體[i]在[t]期的處理狀態(tài)，接受處理時(shí)值為1，未接受處理時(shí)值為0。實(shí)驗(yàn)組在政策實(shí)施前的樣本均值為[β0+β2?E（X）]，政策實(shí)施后的樣本均值為[β0+β1+β2?E（X）]，所以實(shí)驗(yàn)組在政策實(shí)施前后的差異為[β1]，控制組在政策實(shí)施前后的樣本均值均為[β0+β2?E（X）]，故控制組政策實(shí)施前后的差異為0，所以平均處理效應(yīng)即為[β1-0=β1]，它代表了新型城鎮(zhèn)化試點(diǎn)政策對(duì)農(nóng)業(yè)面源污染的凈效應(yīng)。若[β1]小于零，則說明新型城鎮(zhèn)化政策的實(shí)施有利于減少農(nóng)業(yè)面源污染。

（二）變量設(shè)定

1. 被解釋變量：農(nóng)業(yè)面源污染（TE）

本文運(yùn)用應(yīng)用較為廣泛的清單分析法［21-22］測(cè)算我國159個(gè)地級(jí)城市的農(nóng)業(yè)面源污染排放情況。計(jì)算的主要污染物為總氮（TN）、總磷（TP），具體計(jì)算公式如下：

[TE=Ti×λij×ωi] （2）

式（2）中，TE為農(nóng)業(yè)面源污染排放總量（萬噸）；[Ti]為產(chǎn)污單元i的指標(biāo)統(tǒng)計(jì)量，本文的產(chǎn)污單元為氮肥、磷肥、復(fù)合肥；[λij]為單元[i]中污染物[j]的產(chǎn)污系數(shù)，產(chǎn)污系數(shù)根據(jù)化肥折純量的化學(xué)成分計(jì)算，氮肥、磷肥、復(fù)合肥的總氮產(chǎn)污系數(shù)分別為1.00、0、0.33，總磷產(chǎn)污系數(shù)分別為0、0.44、0.15 ［22］；[ωi]為[i]類化肥的流失率，氮肥、磷肥流失率參照史常亮等［5］整理出的中國各省區(qū)化肥流失情況。

2. 核心解釋變量：新型城鎮(zhèn)化試點(diǎn)（DID）

本文的核心解釋變量為新型城鎮(zhèn)化試點(diǎn)，我們將其定義為：某城市實(shí)施新型城鎮(zhèn)化試點(diǎn)的當(dāng)年及之后的年份取值為1，否則為0。

3. 中介變量

環(huán)境規(guī)制（ER）：基于李虹等［23］的做法，采用熵值法計(jì)算環(huán)境規(guī)制綜合指標(biāo)。為此，選取工業(yè)廢水排放量、工業(yè)SO2排放量、工業(yè)煙（粉）塵排放量三個(gè)單項(xiàng)指標(biāo)。首先對(duì)各個(gè)指標(biāo)進(jìn)行標(biāo)準(zhǔn)化處理，然后確定指標(biāo)權(quán)重，最后根據(jù)權(quán)重和標(biāo)準(zhǔn)化數(shù)值計(jì)算環(huán)境規(guī)制強(qiáng)度。

農(nóng)業(yè)綠色技術(shù)進(jìn)步（AGTP）：現(xiàn)有文獻(xiàn)主要通過構(gòu)建環(huán)境效率指標(biāo)［24］、選取綠色全要素生產(chǎn)率［25］衡量綠色技術(shù)進(jìn)步。本文采用第二種方法。針對(duì)農(nóng)業(yè)綠色全要素生產(chǎn)率的度量，學(xué)術(shù)界廣泛采用基于松弛的效率模型（Slacks-Based Measure，SBM）［14］。然而由于ML指數(shù)不具備傳遞性，借鑒吳傳清等［26］的研究采用SBM模型、GML指數(shù)測(cè)度農(nóng)業(yè)綠色全要素生產(chǎn)率。使用MaxDEA軟件進(jìn)行測(cè)度，所使用的數(shù)據(jù)說明如表1。

4. 控制變量

為控制影響農(nóng)業(yè)面源污染的因素，參考相關(guān)文獻(xiàn)，引入如下控制變量：產(chǎn)業(yè)結(jié)構(gòu)（STR）［27］（第一產(chǎn)業(yè)增加值/（第二產(chǎn)業(yè)增加值+第三產(chǎn)業(yè)增加值））；農(nóng)村勞動(dòng)力轉(zhuǎn)移（LT）［7］（第一產(chǎn)業(yè)從業(yè)人員數(shù)量）；化肥施用技術(shù)水平（FTL）［10］（化肥施用折純量/農(nóng)業(yè)總產(chǎn)值）；種植結(jié)構(gòu)（CPS）［28］（糧食作物播種面積/（農(nóng)作物總播種面積-糧食作物播種面積））；富裕程度（INC）［5］（農(nóng)業(yè)總產(chǎn)值/農(nóng)作物總播種面積）。

（三）數(shù)據(jù)來源

以上所運(yùn)用的樣本數(shù)據(jù)來自各省統(tǒng)計(jì)年鑒、各市統(tǒng)計(jì)年鑒、各市國民經(jīng)濟(jì)和社會(huì)發(fā)展統(tǒng)計(jì)公報(bào)、《中國城市統(tǒng)計(jì)年鑒》和《中國區(qū)域經(jīng)濟(jì)統(tǒng)計(jì)年鑒》。本文所運(yùn)用到的農(nóng)業(yè)總產(chǎn)值數(shù)據(jù)利用農(nóng)業(yè)總產(chǎn)值指數(shù)（上年=100）將其折算成以2007年為基準(zhǔn)的不變價(jià)?；鶞?zhǔn)回歸模型中勞動(dòng)力轉(zhuǎn)移（LT）為第一產(chǎn)業(yè)從業(yè)人數(shù)做對(duì)數(shù)處理后的數(shù)值。本文部分缺失數(shù)據(jù)采用線性插值法將其補(bǔ)齊。變量的描述性統(tǒng)計(jì)見表2。

四、實(shí)證分析

（一）基準(zhǔn)回歸分析

根據(jù)基準(zhǔn)模型，本節(jié)量化分析新型城鎮(zhèn)化試點(diǎn)政策的出臺(tái)對(duì)農(nóng)業(yè)面源污染的政策效應(yīng)。由表3可知，在未納入控制變量的情況下（模型1），DID估計(jì)系數(shù)顯示新型城鎮(zhèn)化試點(diǎn)對(duì)農(nóng)業(yè)面源污染的作用在10%的水平下負(fù)向顯著。依次加入控制變量后（模型2—模型6），DID估計(jì)系數(shù)由-0.104變?yōu)?0.115，顯著性從10%提升至5%，說明新型城鎮(zhèn)化試點(diǎn)政策有效緩解了農(nóng)業(yè)面源污染的增加。在控制變量方面，種植結(jié)構(gòu)（CPS）、化肥施用技術(shù)水平（FTL）、農(nóng)村勞動(dòng)力轉(zhuǎn)移（LT）對(duì)農(nóng)業(yè)面源污染的影響顯著為正，說明優(yōu)化作物種植結(jié)構(gòu)、減少化肥施用量、促進(jìn)農(nóng)村地區(qū)勞動(dòng)力轉(zhuǎn)移會(huì)抑制農(nóng)業(yè)面源污染的提升。

（二）穩(wěn)健性檢驗(yàn)

1. 平行趨勢(shì)檢驗(yàn)

多期雙重差分法的前提是在政策實(shí)施之前，實(shí)驗(yàn)組和控制組之間滿足平行趨勢(shì)假設(shè)，即在未試點(diǎn)新型城鎮(zhèn)化政策前，實(shí)驗(yàn)組與控制組城市的農(nóng)業(yè)面源污染變化趨勢(shì)應(yīng)該是平行的或不存在顯著的差異。為此，運(yùn)用事件分析法檢驗(yàn)新型城鎮(zhèn)化試點(diǎn)對(duì)農(nóng)業(yè)面源污染的影響是否滿足平行趨勢(shì)假設(shè)，參照Beck等［19］具體模型設(shè)定如下：

[TEit=δ0+δ1did-9it+δ2did-8it+…+δ14did+4it+φXit+μi+γt+εit]? （3）

式（3）中，[didKit]（K=-9，-8，···，+4）表示新型城鎮(zhèn)化試點(diǎn)批準(zhǔn)這一事件的虛擬變量；K的符號(hào)為負(fù)代表城市[i]試點(diǎn)新型城鎮(zhèn)化之前的第K年，K的符號(hào)為正代表批復(fù)新型城鎮(zhèn)化試點(diǎn)后的第K年，本文以政策實(shí)施前的第10年為基準(zhǔn)組，選取新型城鎮(zhèn)化試點(diǎn)前9年與試點(diǎn)當(dāng)年及之后的4年進(jìn)行估計(jì)。其他變量與式（1）含義相同。檢驗(yàn)結(jié)果如圖1所示。

從圖1可知，在新型城鎮(zhèn)化試點(diǎn)之前，系數(shù)[δk]不具有統(tǒng)計(jì)顯著性，說明在試點(diǎn)開始前，實(shí)驗(yàn)組與控制組城市的農(nóng)業(yè)面源污染變化趨勢(shì)不存在顯著的差異，滿足平行趨勢(shì)假設(shè)；試點(diǎn)開始后，除第2年以外，[δk]均負(fù)向顯著，表明新型城鎮(zhèn)化試點(diǎn)政策的實(shí)施有助于抑制農(nóng)業(yè)面源污染水平的提高。

[-9][-8][-7][-6][-5][-4][-3][-2][-1][0][1][2][3][4][政策時(shí)點(diǎn)][政策動(dòng)態(tài)經(jīng)濟(jì)效應(yīng)][-0.6][-0.5][-0.4][-0.3][-0.2][-0.1][0][0.1][0.2]圖1? ? 平行趨勢(shì)檢驗(yàn)

2. PSM-DID檢驗(yàn)

為避免實(shí)驗(yàn)組與控制組的變動(dòng)趨勢(shì)存在“選擇性偏差”，本文使用PSM—DID方法進(jìn)行穩(wěn)健性檢驗(yàn)。以樣本期內(nèi)75個(gè)試點(diǎn)城市作為實(shí)驗(yàn)組，用控制變量代表協(xié)變量對(duì)實(shí)驗(yàn)組進(jìn)行重新匹配，匹配后的平衡性檢驗(yàn)結(jié)果顯示大多數(shù)協(xié)變量的標(biāo)準(zhǔn)化偏差在5%以下，將無法匹配成功的控制組城市樣本剔除，再按照式（1）進(jìn)行估計(jì)。匹配之后的新樣本回歸結(jié)果顯示，解釋變量DID的估計(jì)系數(shù)并無明顯變化，變量DID估計(jì)系數(shù)仍在5%水平下顯著為負(fù)（表4），這與基準(zhǔn)回歸的結(jié)論基本一致，驗(yàn)證了基準(zhǔn)回歸結(jié)果的穩(wěn)健性。

3. 安慰劑檢驗(yàn)

為進(jìn)一步檢驗(yàn)新型城鎮(zhèn)化試點(diǎn)對(duì)農(nóng)業(yè)面源污染的抑制效應(yīng)是否是由于其他不可觀測(cè)的因素導(dǎo)致的，本文隨機(jī)選擇虛擬實(shí)驗(yàn)組與控制組進(jìn)行安慰劑檢驗(yàn)。具體方法為，從所有樣本城市中隨機(jī)選擇75個(gè)城市作為試點(diǎn)新型城鎮(zhèn)化政策的實(shí)驗(yàn)組，并隨機(jī)為被選中的城市設(shè)定政策年份，設(shè)置偽政策虛擬變量，在此基礎(chǔ)上利用新的樣本組對(duì)（1）式進(jìn)行回歸，將上述步驟重復(fù)500次，得到500個(gè)DID的估計(jì)系數(shù)，將這些系數(shù)與實(shí)際DID系數(shù)進(jìn)行對(duì)比。結(jié)果分布如圖2所示。從圖2可知，偽政策虛擬變量的估計(jì)值集中分布在0附近，并且大部分在10%的水平下不顯著，說明本文的實(shí)證結(jié)果較為穩(wěn)健，其他不可觀測(cè)因素不會(huì)干擾新型城鎮(zhèn)化試點(diǎn)政策對(duì)農(nóng)業(yè)面源污染的影響。

4. 替換被解釋變量

為避免計(jì)量結(jié)果存在偏誤，本文進(jìn)一步采用替換被解釋變量的方法對(duì)本文得出的結(jié)論進(jìn)行穩(wěn)健性檢驗(yàn)?？紤]到新型城鎮(zhèn)化試點(diǎn)可能會(huì)減少農(nóng)業(yè)用地，進(jìn)而減少農(nóng)業(yè)面源污染總量，選取相對(duì)指標(biāo)化肥投入密度（FER）、農(nóng)藥投入密度（PES），前者為單位農(nóng)作物播種面積的化肥投入量，后者為單位農(nóng)作物播種面積的農(nóng)藥投入量，作為農(nóng)業(yè)面源污染排放量的代理變量進(jìn)行穩(wěn)健性分析，回歸結(jié)果見表5。核心解釋變量的系數(shù)符號(hào)和顯著性與前文保持一致，說明新型城鎮(zhèn)化試點(diǎn)顯著減少了農(nóng)業(yè)面源污染。新型城鎮(zhèn)化是以人為核心的城鎮(zhèn)化，它一方面鼓勵(lì)農(nóng)民到城市生存，另一方面又保留了農(nóng)村的耕地和宅基地，因此農(nóng)用土地實(shí)質(zhì)上并沒有發(fā)生顯著變化，因此新型城鎮(zhèn)化通過減少農(nóng)業(yè)用地進(jìn)而減少農(nóng)業(yè)面源污染，這一影響基本不存在，這在下面的實(shí)證結(jié)果中也得到了驗(yàn)證，進(jìn)一步證實(shí)了本文結(jié)論的穩(wěn)健性。

（三）影響機(jī)制的中介效應(yīng)檢驗(yàn)

通過前文理論分析，新型城鎮(zhèn)化可通過加強(qiáng)環(huán)境規(guī)制、促進(jìn)農(nóng)業(yè)綠色技術(shù)進(jìn)步來減少農(nóng)業(yè)面源污染排放，為進(jìn)一步驗(yàn)證上述影響機(jī)制是否存在，本文借鑒溫忠麟等［29］的研究成果構(gòu)建如下中介效應(yīng)檢驗(yàn)?zāi)Ｐ停?/p>

[TEit=α0+α1DIDit+α2Controlit+μi+δt+εit] （4）

[Mit=β0+β1DIDit+β2Controlit+μi+δt+εit] （5）

[TEit=γ0+γ1DIDit+γ2Mit+γ3Controlit+μi+δt+εit] （6）

上述模型中，[Mit]代表可能的中介變量；[μi]和[δt]分別代表控制城市和年份的固定效應(yīng)；[εit]為隨機(jī)擾動(dòng)項(xiàng)。其他變量與式（1）中的含義相同。

中介效應(yīng)檢驗(yàn)結(jié)果見表6，模型7為基準(zhǔn)回歸結(jié)果，本部分不再贅述；模型8表明新型城鎮(zhèn)化試點(diǎn)政策對(duì)環(huán)境規(guī)制具有顯著正向影響，但由于模型9中變量ER的估計(jì)系數(shù)不顯著，根據(jù)溫忠麟等［29］提出的中介效應(yīng)檢驗(yàn)流程，運(yùn)用bootstrap檢驗(yàn)法進(jìn)行系數(shù)乘積的檢驗(yàn)，即檢驗(yàn)原假設(shè)[H0： β1γ2=0]，重復(fù)估計(jì)500次的檢驗(yàn)結(jié)果顯示，系數(shù)乘積在10%的顯著性水平下的置信區(qū)間為［-0.0125，-0.0004］，說明環(huán)境規(guī)制在新型城鎮(zhèn)化抑制農(nóng)業(yè)面源污染中存在中介效應(yīng)，假說1得到驗(yàn)證。由于[γ1]顯著，故存在的是“部分中介效應(yīng)”。模型10顯示新型城鎮(zhèn)化試點(diǎn)顯著促進(jìn)了農(nóng)業(yè)綠色技術(shù)進(jìn)步，模型11表明農(nóng)業(yè)綠色技術(shù)進(jìn)步的提高會(huì)減少農(nóng)業(yè)面源污染排放，且由于[γ1]顯著，表明綠色技術(shù)進(jìn)步是新型城鎮(zhèn)化試點(diǎn)減少農(nóng)業(yè)面源污染排放的部分中介，假說2得到驗(yàn)證。以上結(jié)果表明，新型城鎮(zhèn)化試點(diǎn)政策可通過加強(qiáng)環(huán)境規(guī)制、促進(jìn)農(nóng)業(yè)綠色技術(shù)進(jìn)步減少農(nóng)業(yè)面源污染。

（四）沖擊反應(yīng)的異質(zhì)性分析

新型城鎮(zhèn)化試點(diǎn)作為一項(xiàng)沖擊性政策，對(duì)農(nóng)業(yè)面源污染減排效應(yīng)是否會(huì)在不同區(qū)域或者不同規(guī)模城市表現(xiàn)出差異性？在接下來的異質(zhì)性分析中，我們檢驗(yàn)在新型城鎮(zhèn)化試點(diǎn)引入后，不同地理區(qū)域和不同規(guī)模城市的農(nóng)業(yè)面源污染減排效應(yīng)的差異性。

1. 區(qū)域異質(zhì)性

中國地區(qū)間的經(jīng)濟(jì)狀況、資源稟賦等存在較大差距。本文首先根據(jù)國家政策上的劃分，將樣本城市分為東部城市、中部城市和西部城市，進(jìn)行區(qū)域異質(zhì)性分析，具體結(jié)果如表7所示。新型城鎮(zhèn)化試點(diǎn)政策的實(shí)施有效減少了東部城市的農(nóng)業(yè)面源污染排放，但對(duì)中、西部城市的農(nóng)業(yè)面源污染排放不存在顯著影響。原因可能在于兩方面，一方面東部地區(qū)相較于中西部地區(qū)在經(jīng)濟(jì)實(shí)力、人力資本、技術(shù)基礎(chǔ)等方面具有顯著優(yōu)勢(shì)，通過技術(shù)效應(yīng)來實(shí)現(xiàn)污染減排的激勵(lì)更強(qiáng)［30］。另一方面可能是由于東部城市較中西部城市人均收入水平較高，隨著人均收入水平的提升，城鎮(zhèn)化的推進(jìn)對(duì)農(nóng)業(yè)面源污染的影響會(huì)由擴(kuò)張效應(yīng)轉(zhuǎn)為質(zhì)量效應(yīng)［10］。

2. 城市規(guī)模異質(zhì)性

本文參考《國務(wù)院關(guān)于調(diào)整城市規(guī)模劃分標(biāo)準(zhǔn)的通知》以及2020年《中國城市建設(shè)統(tǒng)計(jì)年鑒》統(tǒng)計(jì)的城區(qū)人口數(shù)，將城市規(guī)模等級(jí)劃分為大規(guī)模城市（城區(qū)人口在100萬以上）、中等規(guī)模城市（城區(qū)人口在50—100萬之間）和小規(guī)模城市（城區(qū)人口在50萬以下），從規(guī)模異質(zhì)性視角，探討新型城鎮(zhèn)化試點(diǎn)政策對(duì)農(nóng)業(yè)面源污染的影響，具體結(jié)果如表7所示。估計(jì)結(jié)果顯示，新型城鎮(zhèn)化試點(diǎn)政策的實(shí)施能夠顯著減少大規(guī)模城市的農(nóng)業(yè)面源污染排放，但對(duì)中、小規(guī)模城市存在不顯著的抑制作用。這可能是因?yàn)榇笠?guī)模城市的產(chǎn)業(yè)集聚水平較高，能夠更好地對(duì)污染物進(jìn)行規(guī)?；幚?，從而提高環(huán)境治理效率和技術(shù)水平，因此，大規(guī)模城市的環(huán)境改善效應(yīng)更大［31］。至此，假說3得到了驗(yàn)證。

五、結(jié)論與建議

“新型城鎮(zhèn)化”是以人為核心的城鎮(zhèn)化，是以生態(tài)宜居、和諧發(fā)展等為基本特征的城鎮(zhèn)化。本文利用2007—2019年我國159個(gè)地級(jí)城市的面板數(shù)據(jù)為樣本，基于清單分析法計(jì)算出我國159個(gè)地級(jí)城市的農(nóng)業(yè)面源污染排放指標(biāo)。將新型城鎮(zhèn)化試點(diǎn)作為一項(xiàng)準(zhǔn)自然實(shí)驗(yàn)，采用多期DID方法研究了新型城鎮(zhèn)化試點(diǎn)對(duì)農(nóng)業(yè)面源污染的沖擊效應(yīng)、影響機(jī)制和異質(zhì)性。研究發(fā)現(xiàn)：第一，我國新型城鎮(zhèn)化試點(diǎn)政策的實(shí)施能有效減少農(nóng)業(yè)面源污染排放，這對(duì)當(dāng)前鄉(xiāng)村振興戰(zhàn)略實(shí)施具有現(xiàn)實(shí)意義；第二，影響機(jī)制檢驗(yàn)發(fā)現(xiàn)，新型城鎮(zhèn)化試點(diǎn)主要依靠促進(jìn)環(huán)境規(guī)制強(qiáng)度的提升、促進(jìn)農(nóng)業(yè)綠色技術(shù)進(jìn)步來減少農(nóng)業(yè)面源污染排放；第三，異質(zhì)性分析發(fā)現(xiàn)，不同地區(qū)、不同規(guī)模城市在新型城鎮(zhèn)化試點(diǎn)對(duì)農(nóng)業(yè)面源污染的沖擊上表現(xiàn)出不同的差異，目前東部地區(qū)和大規(guī)模城市的政策效應(yīng)明顯，但對(duì)中西部地區(qū)、中小規(guī)模城市的政策效應(yīng)尚不明顯。

基于上述結(jié)論，本文提出深化農(nóng)業(yè)面源污染治理的政策建議：第一，進(jìn)一步擴(kuò)大新型城鎮(zhèn)化政策試點(diǎn)范圍。研究發(fā)現(xiàn)，新型城鎮(zhèn)化試點(diǎn)對(duì)農(nóng)業(yè)面源污染排放具有抑制作用，因此，可以在現(xiàn)有試點(diǎn)城市基礎(chǔ)上納入新的指標(biāo)體系或新的發(fā)展理念等，擴(kuò)大試點(diǎn)范圍。第二，在制定新型城鎮(zhèn)化政策時(shí)，應(yīng)當(dāng)將促進(jìn)綠色技術(shù)進(jìn)步作為重要目標(biāo)。理論分析和實(shí)證結(jié)果表明，新型城鎮(zhèn)化試點(diǎn)政策能夠通過促進(jìn)綠色技術(shù)進(jìn)步減少農(nóng)業(yè)面源污染排放，因此，促進(jìn)農(nóng)業(yè)綠色技術(shù)進(jìn)步是減少農(nóng)業(yè)面源污染排放的有效手段。第三，農(nóng)業(yè)面源污染治理是一個(gè)長(zhǎng)期過程。實(shí)證結(jié)果表明，新型城鎮(zhèn)化試點(diǎn)政策通過加強(qiáng)環(huán)境規(guī)制來減少農(nóng)業(yè)面源污染排放。因此，要因地制宜地制定和完善環(huán)境規(guī)制政策，地方政府要基于本地區(qū)的實(shí)際情況，制定和完善農(nóng)業(yè)環(huán)境標(biāo)準(zhǔn)，不斷強(qiáng)化環(huán)境規(guī)制約束，推行標(biāo)準(zhǔn)化生產(chǎn)，引導(dǎo)農(nóng)戶按照地區(qū)標(biāo)準(zhǔn)使用化肥等化學(xué)投入品，控制化肥、農(nóng)藥等化學(xué)投入品的施用，驅(qū)動(dòng)農(nóng)業(yè)生產(chǎn)方式朝著綠色友好型方向轉(zhuǎn)變。第四，新型城鎮(zhèn)化建設(shè)要因地制宜。實(shí)證結(jié)果表明，新型城鎮(zhèn)化試點(diǎn)政策對(duì)東部城市、大規(guī)模城市的農(nóng)業(yè)面源污染治理效果更強(qiáng)。因此，在制定新型城鎮(zhèn)化政策時(shí)，應(yīng)充分認(rèn)識(shí)試點(diǎn)城市在地理位置、人口規(guī)模等方面的差異，采取不同的差異化預(yù)案。

注釋：

①? ? 常住人口城鎮(zhèn)化率數(shù)據(jù)來源于國家統(tǒng)計(jì)局網(wǎng)站。

參考文獻(xiàn)：

［1］? ? Fleming P， Lichtenberg E， Newburn D A.Evaluating Impacts of Agricultural Cost Sharing on Water Quality：Additionality， Crowding in， and Slippage［J］.Journal of Environmental Economics and Management， 2018（92）：1-19.

［2］? ? 張艦，亞伯拉罕·艾賓斯坦，瑪格麗特·麥克米倫，等.農(nóng)村勞動(dòng)力轉(zhuǎn)移、化肥過度使用與環(huán)境污染［J］.經(jīng)濟(jì)社會(huì)體制比較，2017（3）：149-160.

［3］? ? 馬道明.城鎮(zhèn)化背景下農(nóng)村環(huán)境持續(xù)惡化的內(nèi)生性探析：基于蘇南Q鎮(zhèn)的調(diào)查［J］.學(xué)海，2018（4）：53-58.

［4］? ? 史桂芬，沈淘淘.新型城鎮(zhèn)化背景下農(nóng)業(yè)轉(zhuǎn)移人口社會(huì)融合路徑［J］.東北師大學(xué)報(bào)（哲學(xué)社會(huì)科學(xué)版），2021（3）：103-114.

［5］? ? 史常亮，李赟，朱俊峰.勞動(dòng)力轉(zhuǎn)移、化肥過度使用與面源污染［J］.中國農(nóng)業(yè)大學(xué)學(xué)報(bào)，2016（5）：169-180.

［6］? ? 廖煒，李璐，楊偉，等.城鎮(zhèn)化過程中的流域面源污染時(shí)空變化［J］.長(zhǎng)江流域資源與環(huán)境，2018（8）：1776-1783.

［7］? ? 羅海平，何志文，胡學(xué)英.城鎮(zhèn)化對(duì)種植業(yè)面源污染影響的中介效應(yīng)分析［J］.中國生態(tài)農(nóng)業(yè)學(xué)報(bào)， 2021（9）：1625-1635.

［8］? ? 徐承紅，薛蕾.農(nóng)業(yè)產(chǎn)業(yè)集聚與農(nóng)業(yè)面源污染：基于空間異質(zhì)性的視角［J］.財(cái)經(jīng)科學(xué)，2019（8）：82-96.

［9］? ? 薛蕾，廖祖君，王理.城鎮(zhèn)化與農(nóng)業(yè)面源污染改善：基于農(nóng)民收入結(jié)構(gòu)調(diào)節(jié)作用的空間異質(zhì)性分析［J］.農(nóng)村經(jīng)濟(jì)，2019（7）：55-63.

［10］? 欒健，韓一軍.城鎮(zhèn)化會(huì)加劇化肥面源污染嗎：基于門檻效應(yīng)與空間溢出的雙重視角［J］.中國農(nóng)業(yè)大學(xué)學(xué)報(bào)，2020（5）：174-186.

［11］? Chen X， Chang C P.Fiscal Decentralization， Environmental Regulation， and Pollution：A Spatial Investigation［J］.Environmental Science and Pollution Research， 2020（25）：31946-31968.

［12］? 秦天，彭玨，鄧宗兵，等.環(huán)境分權(quán)、環(huán)境規(guī)制對(duì)農(nóng)業(yè)面源污染的影響［J］.中國人口·資源與環(huán)境， 2021（2）：61-70.

［13］? 王華星，石大千.新型城鎮(zhèn)化有助于緩解霧霾污染嗎：來自低碳城市建設(shè)的經(jīng)驗(yàn)證據(jù)［J］.山西財(cái)經(jīng)大學(xué)學(xué)報(bào)，2019（10）：15-27.

［14］? 閆桂權(quán)，何玉成，張曉恒.綠色技術(shù)進(jìn)步，農(nóng)業(yè)經(jīng)濟(jì)增長(zhǎng)與污染空間溢出：來自中國農(nóng)業(yè)水資源利用的證據(jù)［J］.長(zhǎng)江流域資源與環(huán)境，2019（12）：2921-2935.

［15］? 王一格，王海燕，鄭永林，等.農(nóng)業(yè)面源污染研究方法與控制技術(shù)研究進(jìn)展［J］.中國農(nóng)業(yè)資源與區(qū)劃，2021（1）：25-33.

［16］? 陳海波，姜娜娜，劉潔.新型城鎮(zhèn)化試點(diǎn)政策對(duì)區(qū)域生態(tài)環(huán)境的影響：基于PSM-DID的實(shí)證檢驗(yàn)［J］.城市問題，2020（8）：33-41.

［17］? 王巧，佘碩.城市異質(zhì)性視角下中國低碳試點(diǎn)政策的綠色增長(zhǎng)效應(yīng)評(píng)估［J］.軟科學(xué)，2020（9）：1-8.

［18］? 胡求光，馬勁韜.低碳城市試點(diǎn)政策對(duì)綠色技術(shù)創(chuàng)新效率的影響研究：基于創(chuàng)新價(jià)值鏈視角的實(shí)證檢驗(yàn)［J］.社會(huì)科學(xué)，2022（1）：62-72.

［19］? Beck T， Levine R， Levkov A.Big Bad Banks：The Winners and Losers from Bank Deregulation in The United States［J］.The Journal of Finance， 2010（5）：1637-1667.

［20］? Cheng J， Yi J， Dai S， et al.Can Low-carbon City Construction Facilitate Green Growth：Evidence from Chinas Pilot Low-carbon City Initiative［J］.Journal of Cleaner Production， 2019（5）：1158-1170.

［21］? 陳敏鵬，陳吉寧，賴斯蕓.中國農(nóng)業(yè)和農(nóng)村污染的清單分析與空間特征識(shí)別［J］.中國環(huán)境科學(xué)， 2006（6）：751-755.

［22］? 賴斯蕓，杜鵬飛，陳吉寧.基于單元分析的非點(diǎn)源污染調(diào)查評(píng)估方法［J］.清華大學(xué)學(xué)報(bào)（自然科學(xué)版），2004（9）：1184-1187.

［23］? 李虹，鄒慶.環(huán)境規(guī)制，資源稟賦與城市產(chǎn)業(yè)轉(zhuǎn)型研究：基于資源型城市與非資源型城市的對(duì)比分析［J］.經(jīng)濟(jì)研究，2018（11）：182-198.

［24］? Levidow L， Lindgaard-jorgensen P， Nilsson A， et al.Process Eco-innovation：Assessing Meso-level Eco-efficiency in Industrial Water-service Systems［J］.Journal of Cleaner Production， 2015（1）：54-65.

［25］? Beltran-esteve M， Picazo-tadeo A J.Assessing Environmental Performance Trends in the Transport Industry：Eco-innovation or Catching-up［J］.Energy Economics， 2015（8）：570-580.

［26］? 吳傳清，宋子逸.長(zhǎng)江經(jīng)濟(jì)帶農(nóng)業(yè)綠色全要素生產(chǎn)率測(cè)度及影響因素研究［J］.科技進(jìn)步與對(duì)策， 2018（17）：35-41.

［27］? 姜松，周潔，邱爽.適度規(guī)模經(jīng)營是否能抑制農(nóng)業(yè)面源污染：基于動(dòng)態(tài)門檻面板模型的實(shí)證［J］.農(nóng)業(yè)技術(shù)經(jīng)濟(jì)，2021（7）：33-48.

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

［29］? 溫忠麟，葉寶娟.中介效應(yīng)分析：方法和模型發(fā)展［J］.心理科學(xué)進(jìn)展，2014（5）：731-745.

［30］? 臧傳琴，孫鵬.低碳城市建設(shè)促進(jìn)了地方綠色發(fā)展嗎：來自準(zhǔn)自然實(shí)驗(yàn)的經(jīng)驗(yàn)證據(jù)［J］.財(cái)貿(mào)研究， 2021（10）：27-40.

［31］? 劉習(xí)平，宋德勇.城市產(chǎn)業(yè)集聚對(duì)城市環(huán)境的影響［J］.城市問題，2013（3）：9-15.

責(zé)任編輯：倪貝貝

（E-mail：shellni@ 163. com）

Impact of New-type Urbanization Pilot on Agricultural Non-point Source

Pollution： Verification Using DID with Multiple Time Periods

GAO Qisheng， DONG Zhiqiang

（The School of Economics， Qingdao University， Qingdao Shandong 266061）

Abstract： To promote the green development in agriculture and rural area， as well as the rural revitalization strategy， it is of great significance to explore how new type of urbanization affects agricultural non-point source pollution（ANSP）. Based on the panel data of 159 prefecture cities of China between 2007 and 2019， the influential effect of new-type urbanization pilot policies on ANSP is evaluated using Difference-in-Differences （DID） with Multiple Time Periods. The result shows that the new-type urbanization pilot policies can effectively reduce ANSP discharge. An analysis of the influencing mechanism reveals that the policies relieve ANSP discharge by enhancing environmental regulation and agricultural green technology. An impact analysis reveals that the impact on new-type urbanization pilot differs significantly depending on the citys size and location in the eastern， middle， or western part of the country. The current new-type urbanization pilot policies can effectively reduce the ANSP discharge of eastern cities and large cities， while the influence is insignificant on middle and western cities and cities of medium and small sizes.

Key words： agricultural non-point source pollution（ANSP）； new type of urbanization； quasi-natural experiment； DID with Multiple Time Periods； internal mechanism