江求川，張錕，湯凱

（鄭州大學(xué) 商學(xué)院，河南 鄭州 450001）

1 問題的提出

隨著我國探索生態(tài)優(yōu)先與綠色發(fā)展協(xié)同推進(jìn)路子的步伐不斷加快，我國北部地區(qū)東中西板塊區(qū)域發(fā)展分散、區(qū)域治理體系無法有效聯(lián)動、區(qū)域間發(fā)展差異較大等問題突顯，整體區(qū)域亟須一個生態(tài)經(jīng)濟(jì)帶國家戰(zhàn)略，黃河經(jīng)濟(jì)帶開始被廣大學(xué)者關(guān)注[1-2]。2019 年3 月全國兩會期間，河南、山東等九省份政協(xié)主席聯(lián)名向全國政協(xié)提交提案，建議將黃河生態(tài)經(jīng)濟(jì)帶上升為國家戰(zhàn)略。2019年9 月，習(xí)近平總書記在河南鄭州主持黃河流域生態(tài)保護(hù)和高質(zhì)量發(fā)展座談會并發(fā)表重要講話時指出，要共同抓好大保護(hù)，推動整個黃河流域高質(zhì)量發(fā)展。2020 年5月，李克強(qiáng)總理在政府工作報告中提出，計劃編制黃河流域生態(tài)保護(hù)和高質(zhì)量發(fā)展規(guī)劃綱要。因此，針對黃河經(jīng)濟(jì)帶生態(tài)保護(hù)和綠色發(fā)展進(jìn)行深入研究，對于加快推進(jìn)黃河經(jīng)濟(jì)帶高質(zhì)量發(fā)展具有重要現(xiàn)實意義。

黃河經(jīng)濟(jì)帶與長江經(jīng)濟(jì)帶的最大區(qū)別在于，長江經(jīng)濟(jì)帶覆蓋區(qū)域可以憑借長江這個“黃金水道”便捷聯(lián)通世界各地，極大提高了整體區(qū)域開放程度。然而到21 世紀(jì)初，黃河由于北方降水量較少、水土流失嚴(yán)重、河道落差大等原因，造成河道淤積斷流，很多地區(qū)已不具備通航能力[2]。王海江等[3]提出運(yùn)用現(xiàn)代交通工具如航空、高鐵等改善和替代流域內(nèi)的通航能力。但是，還未將航空運(yùn)輸在經(jīng)濟(jì)帶的地位提升到一定高度。本文提出：以臨空經(jīng)濟(jì)區(qū)和主要港口為核心點和綜合交通樞紐，以高速鐵路和高速公路為支撐和連接線，構(gòu)建區(qū)域立體式現(xiàn)代交通運(yùn)輸體系，以便助推產(chǎn)業(yè)優(yōu)化升級和綠色產(chǎn)業(yè)集群形成，最終協(xié)同生態(tài)保護(hù)推動全區(qū)域經(jīng)濟(jì)高質(zhì)量發(fā)展。整個體系構(gòu)建的問題核心點集中在利用空港資源發(fā)展臨空經(jīng)濟(jì)，強(qiáng)化黃河經(jīng)濟(jì)帶集疏運(yùn)服務(wù)功能，這是否能夠促進(jìn)整個經(jīng)濟(jì)帶的綠色發(fā)展？

臨空經(jīng)濟(jì)對黃河經(jīng)濟(jì)帶本地城市和周邊城市綠色發(fā)展效率是否有影響？其作用效應(yīng)是正還是負(fù)？本文針對這些問題研究黃河經(jīng)濟(jì)帶臨空經(jīng)濟(jì)與綠色發(fā)展效率之間的因果關(guān)系，并考察兩者的空間溢出效應(yīng)，以期為政府科學(xué)規(guī)劃黃河經(jīng)濟(jì)帶生態(tài)保護(hù)和高質(zhì)量發(fā)展提供一定理論參考。

在現(xiàn)存文獻(xiàn)對黃河經(jīng)濟(jì)帶生態(tài)保護(hù)及其高質(zhì)量發(fā)展研究基礎(chǔ)上，本文在以下幾個方面做出創(chuàng)新。首先，鑒于黃河經(jīng)濟(jì)帶不同于長江經(jīng)濟(jì)帶擁有“黃金水道”，無法支撐黃河經(jīng)濟(jì)帶可持續(xù)地高質(zhì)量發(fā)展，因此，本文提出構(gòu)建以臨空經(jīng)濟(jì)和主要港口為綜合樞紐，以高速鐵路和高速公路為支撐和連接線，以倉儲業(yè)集裝箱物流強(qiáng)化全經(jīng)濟(jì)帶疏運(yùn)能力的區(qū)域立體式現(xiàn)代化交通運(yùn)輸體系。其次，由于以往衡量臨空經(jīng)濟(jì)發(fā)展指標(biāo)的航空客運(yùn)量和貨運(yùn)量會受到當(dāng)?shù)亟?jīng)濟(jì)和綠色發(fā)展水平的影響，這種反向因果形成的內(nèi)生性問題會造成后續(xù)模型估計結(jié)果有偏。因此，本文根據(jù)Nunn &Qian[4]和Qian[5]的雙重差分設(shè)計方法，運(yùn)用2008 年11 月我國政府為應(yīng)對金融危機(jī)、擴(kuò)大內(nèi)需推出的一系列包括加快機(jī)場等公共基礎(chǔ)設(shè)施建設(shè)在內(nèi)的政策效應(yīng)沖擊解決航空客運(yùn)量和運(yùn)輸量與區(qū)域綠色發(fā)展效率之間的內(nèi)生性問題，當(dāng)然文中也控制了2008 年金融危機(jī)的沖擊。再次，以往文獻(xiàn)衡量綠色發(fā)展效率主流做法是使用非期望產(chǎn)出的超效率SBM（Slack Based Model）[6-7]。但是，其期望產(chǎn)出僅包括地區(qū)生產(chǎn)總值，對于人民生活水平、政府財政收入、區(qū)域產(chǎn)業(yè)結(jié)構(gòu)優(yōu)化等產(chǎn)出未加考慮。根據(jù)綠色發(fā)展內(nèi)涵[8]，將以上因素綜合考慮加入SBM 模型中測度黃河經(jīng)濟(jì)帶綠色發(fā)展效率。最后，本文構(gòu)建了空間面板Durbin 模型和將DID 的設(shè)計納入其中形成空間DID 模型，討論了黃河經(jīng)濟(jì)帶臨空經(jīng)濟(jì)區(qū)對區(qū)域綠色發(fā)展的影響效應(yīng)及其空間溢出效應(yīng)。鑒于可能存在的空間異質(zhì)性問題，本文采用地理加權(quán)回歸（GWR）方法緩解這一問題。

2 理論分析和假說提出

臨空經(jīng)濟(jì)是指以航空運(yùn)輸為核心，吸引生產(chǎn)、資本、技術(shù)、人才等要素不斷集聚，具有自我造血機(jī)制的多功能新型經(jīng)濟(jì)形態(tài)[9]。Conway[10]最先提出臨空經(jīng)濟(jì)的概念和內(nèi)涵。在此基礎(chǔ)上，一些研究發(fā)現(xiàn)大型運(yùn)輸機(jī)場通過其集散功能和乘數(shù)效應(yīng)吸引著與航空運(yùn)輸相關(guān)聯(lián)的產(chǎn)業(yè)向機(jī)場周邊區(qū)域集中。對樞紐機(jī)場的研究認(rèn)為航空運(yùn)輸能夠顯著推動周邊產(chǎn)業(yè)的發(fā)展[11]。Lindsay &Kasarda[12]更是提出了航空大都市的理念，指出臨空經(jīng)濟(jì)發(fā)展必須注重與周邊地區(qū)的協(xié)同關(guān)系。在整個城市化的進(jìn)程中，臨空經(jīng)濟(jì)將會成為新的增長極和區(qū)域經(jīng)濟(jì)發(fā)展的催化劑[13-14]，并且推動所在城市向多功能航空都市轉(zhuǎn)變[15]。國內(nèi)對于臨空經(jīng)濟(jì)的研究較晚，如魏曉芳等[16]研究了臨空經(jīng)濟(jì)區(qū)規(guī)劃的產(chǎn)業(yè)選擇和產(chǎn)業(yè)空間布局問題。另外，不同類型臨空經(jīng)濟(jì)對其所在區(qū)域經(jīng)濟(jì)增長的影響效應(yīng)也不同[17]。基礎(chǔ)設(shè)施建設(shè)水平、人力資本以及技術(shù)進(jìn)步對臨空經(jīng)濟(jì)的發(fā)展有著顯著促進(jìn)作用[18]。班奕[19]發(fā)現(xiàn)臨空經(jīng)濟(jì)與自貿(mào)區(qū)之間存在明顯的耦合關(guān)系。

以上的文獻(xiàn)基本都是研究臨空經(jīng)濟(jì)對所在區(qū)域的產(chǎn)業(yè)和經(jīng)濟(jì)增長的影響。然而，研究臨空經(jīng)濟(jì)區(qū)是否影響區(qū)域綠色發(fā)展效率的文獻(xiàn)還很少，更鮮有文獻(xiàn)研究臨空經(jīng)濟(jì)區(qū)對周邊城市綠色發(fā)展的溢出效應(yīng)。不僅如此，在現(xiàn)有的文獻(xiàn)中大多采用臨空經(jīng)濟(jì)區(qū)內(nèi)航空客運(yùn)量和貨運(yùn)量來代替臨空經(jīng)濟(jì)區(qū)的發(fā)展水平。然而，綠色發(fā)展效率較高的地區(qū)其臨空經(jīng)濟(jì)發(fā)展水平也相應(yīng)提升，造成嚴(yán)重的內(nèi)生性問題。本文正是基于這些問題以及研究黃河經(jīng)濟(jì)帶生態(tài)保護(hù)和高質(zhì)量發(fā)展的現(xiàn)實意義下決定深入探討該經(jīng)濟(jì)帶的綠色發(fā)展變動問題和通過雙重差分的設(shè)計解決可能存在的內(nèi)生性問題。

臨空經(jīng)濟(jì)作為城市內(nèi)新型經(jīng)濟(jì)形態(tài)，在三個方面影響本地城市經(jīng)濟(jì)的綠色發(fā)展。其一，臨空經(jīng)濟(jì)通過重構(gòu)區(qū)域產(chǎn)業(yè)結(jié)構(gòu)，淘汰高污染高耗能產(chǎn)業(yè)，推動區(qū)域經(jīng)濟(jì)綠色發(fā)展。航空運(yùn)輸貨物特點主要為價值高、重量輕、時效需要高，就要求周邊產(chǎn)業(yè)逐步優(yōu)化轉(zhuǎn)型為高新技術(shù)產(chǎn)業(yè)和時間偏好型產(chǎn)業(yè)[9]。這兩類產(chǎn)業(yè)能較好協(xié)同生態(tài)保護(hù)，提高區(qū)域綠色發(fā)展效率。其二，臨空經(jīng)濟(jì)具有較強(qiáng)溢出效應(yīng)和極化效應(yīng)，引導(dǎo)生產(chǎn)、技術(shù)、人才、資本集聚，提高區(qū)域生產(chǎn)效率，推動區(qū)域高質(zhì)量發(fā)展，即集聚促進(jìn)效應(yīng)。其三，這種集聚效應(yīng)也可能會造成擁堵效應(yīng)，阻礙區(qū)域的綠色發(fā)展[20]。人口和生產(chǎn)的空間集聚造成城市污染排放和碳排放在一定程度上加劇，加速城市集聚區(qū)資源消耗，造成當(dāng)?shù)厣鷳B(tài)環(huán)境惡化[21-22]。

同時，臨空經(jīng)濟(jì)區(qū)從兩個路徑影響周邊城市綠色發(fā)展。一方面，臨空經(jīng)濟(jì)城市通過高速鐵路和高速公路推動周邊城市產(chǎn)業(yè)轉(zhuǎn)型升級，促進(jìn)周邊城市綠色發(fā)展[23]。但是，另一方面，如果本地城市與鄰接城市產(chǎn)業(yè)間構(gòu)成競爭關(guān)系，臨空經(jīng)濟(jì)城市憑借靠近臨空經(jīng)濟(jì)區(qū)的優(yōu)勢打壓鄰接城市競爭性產(chǎn)業(yè)[24]，造成鄰接城市綠色產(chǎn)業(yè)和優(yōu)質(zhì)資源向臨空經(jīng)濟(jì)城市集聚，抑制周邊城市生產(chǎn)效率的提高[25]。綜上所述，臨空經(jīng)濟(jì)區(qū)對本地城市和周邊城市綠色發(fā)展的影響效應(yīng)方向和顯著性并不確定，沒有一個明確的結(jié)論?；谝陨戏治?，因此本文提出：

假說一：臨空經(jīng)濟(jì)對本地區(qū)域綠色發(fā)展效率有正向影響。

假說二：臨空經(jīng)濟(jì)遏制鄰接城市綠色發(fā)展。

3 研究策略

3.1 綠色發(fā)展效率的測度方法

Tone[26-27]提出將非零松弛考慮到效率值的計算中構(gòu)建非徑向和非角度的SBM 模型可以彌補(bǔ)傳統(tǒng)前沿分析方法的缺陷；并且非期望產(chǎn)出問題也被考慮到SBM模型中。本文參照以上文獻(xiàn)方法，建立Super-SBMUndesirable 測度黃河經(jīng)濟(jì)帶綠色發(fā)展效率。

基于以上等式（2）條件，可以定義SBM 為：

式中：目標(biāo)函數(shù)值ρ*是關(guān)于s-、sg、sb嚴(yán)格遞減的，且ρ*∈[0,1]。當(dāng)且僅當(dāng)ρ*=1 時認(rèn)為生產(chǎn)決策單元是有效率的，ρ*＜1 時，生產(chǎn)決策單元未達(dá)到最優(yōu)，投入產(chǎn)出結(jié)構(gòu)有必要再改進(jìn)。

等式（3）為非線性規(guī)劃模型，可以通過Charnes-Cooper 方法轉(zhuǎn)變?yōu)榫€性規(guī)劃問題后求解[28]，因此定義目標(biāo)函數(shù)τ的最優(yōu)解作為的超效率值，形式如下（t為未知參數(shù)，為超效率得分）：

3.2 模型構(gòu)建

本文運(yùn)用DID 的設(shè)計將2009 年政府對于機(jī)場、公路等基礎(chǔ)設(shè)施建設(shè)的大規(guī)模投資作為臨空經(jīng)濟(jì)發(fā)展水平的外生沖擊來解決潛在的內(nèi)生性問題，構(gòu)建基礎(chǔ)模型為：

式中：gdet為第t期綠色發(fā)展效率（green development efficiency）的n×1 矩陣，n為空間單元數(shù)；IT×1為T×1階單位矩陣，這里T為總時期數(shù)；ιn×1為n×1 階單位矩陣；timet設(shè)定為一個虛擬變量，當(dāng)?shù)趖期樣本處于2009 年之后（包括2009 年）為1，否則為0；treatmentt設(shè)定為一個虛擬變量，當(dāng)?shù)趖期樣本點為臨空經(jīng)濟(jì)區(qū)取值為1，否則為0；Xt為第t期n×k階控制變量矩陣，k為控制變量數(shù)；τT×1為T×1 階時間固定效應(yīng)矩陣；μn×1為n×1 階空間固定效應(yīng)矩陣；為隨機(jī)擾動項矩陣，滿足零條件均值假設(shè)；α0、λ0為待估參數(shù)，β0為k×1 階待估參數(shù)矩陣。這里，“”為克羅內(nèi)克積法則，“ ”為哈達(dá)瑪積法則。

在空間Durbin 模型的基礎(chǔ)上[29]，本文構(gòu)建如下空間DID 設(shè)計模型[30-31]，見模型（6）：

式中：W為n×n階空間權(quán)重矩陣；IT為T×T階單位矩陣；為隨機(jī)擾動項矩陣，滿足零條件均值假設(shè)；δ1、α1、λ1、φ1為待估參數(shù)，β1和θ1為k×1 階待估參數(shù)矩陣。其他變量和參數(shù)同模型（5）所注。

由于實際估計中各變量特征在空間分布上呈現(xiàn)異質(zhì)性，造成估計結(jié)果非穩(wěn)健性和有偏的問題，可以通過地理加權(quán)回歸的方法緩解這一問題[32]。為了緩解空間異質(zhì)性的問題，本文在模型（5）的基礎(chǔ)上，建立如下模型形式：

式中：(ν,u)為根據(jù)坐標(biāo)調(diào)整系數(shù)權(quán)重的n×1 階矩陣；為隨機(jī)擾動項，滿足零條件均值假設(shè)。α2、λ2為待估參數(shù)，β2為k×1 階待估參數(shù)矩陣；其他變量和參數(shù)同模型（5）和模型（6）所注。

GWR 回歸較特殊的一點是其必須確定一個帶寬，通過帶寬窗口來掃描全局空間單元，給出參數(shù)估計值中各分位點的值；根據(jù)給定帶寬和衰減函數(shù)確定某一空間單元回歸的空間權(quán)重，以此調(diào)整估計參數(shù)的權(quán)重。一般有高斯（Gauss）衰減函數(shù)[見式（8）]，權(quán)重函數(shù)形式如下：

式中：dij為空間單元i和j之間的距離；bandwidth為帶寬長度，一般通過交叉驗證CV 的最小值確定最優(yōu)帶寬長度。

4 數(shù)據(jù)來源和變量處理

為研究黃河經(jīng)濟(jì)帶臨空經(jīng)濟(jì)與綠色發(fā)展效率的因果關(guān)系，選取研究樣本為2003—2017 年黃河經(jīng)濟(jì)帶82 個地級市，共計1 230 個樣本。由于海東市2013 年撤地建市，其2013 年及以前數(shù)據(jù)缺失嚴(yán)重，故將海東市刪除。本文收集數(shù)據(jù)主要來源于2002—2018 年的《中國城市統(tǒng)計年鑒》。部分指標(biāo)數(shù)據(jù)存在缺失，比如2017 年全社會固定資產(chǎn)投資額、名義GDP、按三次產(chǎn)業(yè)分類的產(chǎn)值比重等一些指標(biāo)，這些指標(biāo)數(shù)據(jù)從2018 年各省份統(tǒng)計年鑒中獲取。除此之外，還有部分缺失值，本文采取插值法補(bǔ)充。另外，后續(xù)研究分析還涉及臨空經(jīng)濟(jì)城市的航空運(yùn)輸量和貨運(yùn)量指標(biāo)，數(shù)據(jù)從中國民航局2005—2017 年民航機(jī)場生產(chǎn)統(tǒng)計年報中獲得。

針對研究的被解釋變量——綠色發(fā)展效率，本文運(yùn)用Super-SBM-Undesirable 進(jìn)行測度，這一測度方法需要確定投入要素矩陣、期望產(chǎn)出矩陣以及非期望產(chǎn)出矩陣，具體的指標(biāo)選取和變量處理說明如下：

（1）投入要素矩陣。資本存量作為其中一個投入要素，采用永續(xù)盤存法計算，其計算公式以及基期資本存量估計公式如（9）[33]：

式中：kit為第t期i區(qū)域的資本存量；ki0為基期i區(qū)域的資本存量；ψit為第t期i區(qū)域的固定資產(chǎn)折舊率；ψi0為基期i區(qū)域的固定資產(chǎn)折舊率；Iit為第t期i區(qū)域的全社會固定資產(chǎn)投資額；Ii0為基期i區(qū)域的全社會固定資產(chǎn)投資額；Pit為第t期i區(qū)域的固定資產(chǎn)投資價格指數(shù)（以2000 年=100 進(jìn)行平減）；為i區(qū)域t1到t2時期固定資產(chǎn)投資的平均幾何增長率。最后將計算得到的資本存量單位換算為百億元。

其他投入要素為：勞動投入使用各地市年末城鎮(zhèn)單位就業(yè)人數(shù)衡量，單位為十萬人；固定資產(chǎn)折舊額是根據(jù)分省固定資產(chǎn)折舊率計算各地級市的固定資產(chǎn)折舊額，單位為百億元；科技投入為地區(qū)財政預(yù)算內(nèi)支出項目中科技支出項目獲得，單位為千萬元。以上數(shù)據(jù)指標(biāo)從《中國城市統(tǒng)計年鑒》獲得。

（2）期望產(chǎn)出矩陣。本文綜合考慮地區(qū)綠色發(fā)展指標(biāo)，使用實際地區(qū)生產(chǎn)總值衡量地市經(jīng)濟(jì)發(fā)展總體實力，單位為百億元；使用實際人均地區(qū)生產(chǎn)總值衡量人民生活水平，單位為萬元/人；使用一般公共預(yù)算內(nèi)收入反映政府財政能力，單位為十億元；以及用第三產(chǎn)業(yè)與第二產(chǎn)業(yè)產(chǎn)值比重反映地區(qū)產(chǎn)業(yè)結(jié)構(gòu)。最終，將以上指標(biāo)作為期望產(chǎn)出。

（3）非期望產(chǎn)出矩陣。大部分文獻(xiàn)采用了工業(yè)廢水、工業(yè)二氧化硫、工業(yè)廢氣排放量[7]作為非期望產(chǎn)出，本文考慮到三個指標(biāo)的一些問題，比如地區(qū)環(huán)境治理因素沒有考慮到指標(biāo)體系中和一些行政或者企業(yè)的干擾等問題，因此，本文選用PM2.5（單位ug/m3）指標(biāo)衡量地市環(huán)境污染水平，作為非期望產(chǎn)出，由于PM2.5數(shù)據(jù)來自衛(wèi)星采集現(xiàn)實觀測數(shù)據(jù)，排除了人為干擾問題和考慮到環(huán)境治理因素，該數(shù)據(jù)來自達(dá)爾豪斯大學(xué)大氣成分分析組，數(shù)據(jù)采用NASA 衛(wèi)星的氣溶膠數(shù)據(jù)，為一年期均值[34]。

最終將投入要素矩陣、期望產(chǎn)出矩陣、非期望產(chǎn)出矩陣代入綠色發(fā)展效率測度模型后獲得2003—2017 年黃河經(jīng)濟(jì)帶82 個地級市的綠色發(fā)展效率。

本文核心變量選取2009 年我國政府為應(yīng)對金融危機(jī)出臺一攬子包括大規(guī)模投資機(jī)場等基礎(chǔ)設(shè)施建設(shè)的政策沖擊作為臨空經(jīng)濟(jì)發(fā)展水平的代理變量。一攬子經(jīng)濟(jì)刺激計劃包括安排中西線干線機(jī)場和支線機(jī)場建設(shè)，大部分干線城市獲批以及購置更高規(guī)格的飛機(jī)，推動中西部城市臨空經(jīng)濟(jì)的快速發(fā)展。最終，通過考察黃河經(jīng)濟(jì)帶82 個地級市處理組是否設(shè)置或者規(guī)劃臨空經(jīng)濟(jì)區(qū)以及臨空經(jīng)濟(jì)區(qū)發(fā)展程度，選取12 個臨空經(jīng)濟(jì)城市，包括西安、鄭州、青島、濟(jì)南、蘭州、太原、呼和浩特、銀川、煙臺、西寧、鄂爾多斯、臨沂，控制組為其他70個地級城市；設(shè)置虛擬變量treatmentt，如果為臨空經(jīng)濟(jì)城市，取值為1，否則為0。同時設(shè)置虛擬變量timet，如果樣本在2009 年發(fā)生以后，取值為1，否則取值為0。

本文選取控制變量為：產(chǎn)業(yè)集聚度（即diat），一般的思路是用市區(qū)的規(guī)模以上工業(yè)企業(yè)數(shù)與全市規(guī)模以上工業(yè)企業(yè)數(shù)之比衡量，但考慮到地級市工業(yè)企業(yè)數(shù)規(guī)模差異造成的內(nèi)生性沖擊，本文使用區(qū)位熵來緩解這一問題。其形式如下：

式中：diait為第t期i區(qū)域的產(chǎn)業(yè)集聚度；uinit為第t期i區(qū)域的市區(qū)工業(yè)企業(yè)數(shù)（urban industrial number）；tinit為第t期i區(qū)域的全地區(qū)工業(yè)企業(yè)數(shù)（total industrial number）。

本文還控制了人均實際GDP 及其平方項，人均實際GDP 為人均名義GDP 按GDP 平減指數(shù)（以2000 年=100）換算，單位為萬元/人；產(chǎn)業(yè)高級化指標(biāo)為第三產(chǎn)業(yè)與第二產(chǎn)業(yè)產(chǎn)值比重；產(chǎn)業(yè)合理化指數(shù)可以用下式泰爾指數(shù)求倒數(shù)得到[35]。

式中：RISit為第t期i區(qū)域的產(chǎn)業(yè)合理化指數(shù)；TLit為第t期i區(qū)域的泰爾指數(shù)；Vitq表示第t期i區(qū)域的第q產(chǎn)業(yè)產(chǎn)值；GDPit表示第t期i區(qū)域的地區(qū)生產(chǎn)總值；Litq為第t期i區(qū)域的第q產(chǎn)業(yè)城鎮(zhèn)單位就業(yè)職工數(shù)；Lit為第t期i區(qū)域城鎮(zhèn)單位就業(yè)職工總數(shù)。

另外，其余控制變量中創(chuàng)新水平用科技支出占地區(qū)GDP 的比重衡量；考慮到政府干預(yù)對當(dāng)?shù)鼐G色發(fā)展效率有較大影響[36]，因此還控制了政府干預(yù)度，用地方財政預(yù)算內(nèi)支出占GDP 的比重衡量[20]；由于2009 年政策效應(yīng)一定程度上受到2008 年金融危機(jī)的影響，本文通過添加經(jīng)濟(jì)沖擊指標(biāo)來控制，使用地市本期相對于上一期實際GDP 增長率變動量來衡量，人口增長率選用年末戶籍人數(shù)的增長率衡量。

本文設(shè)定空間權(quán)重矩陣為逆經(jīng)濟(jì)距離和逆地理距離混合權(quán)重矩陣①本文還對逆經(jīng)濟(jì)距離權(quán)重矩陣、鄰接權(quán)重矩陣以及逆地理距離權(quán)重矩陣進(jìn)行了檢驗，發(fā)現(xiàn)與本文結(jié)果一致。，其矩陣形式如下：

5 研究結(jié)果

5.1 主要結(jié)果

為緩解綠色發(fā)展效率與臨空經(jīng)濟(jì)發(fā)展水平的反向因果而造成的內(nèi)生性問題，本文引入了2009 年大規(guī)模機(jī)場等基建投資政策作為臨空經(jīng)濟(jì)發(fā)展外生沖擊的代理變量，以解決內(nèi)生性造成的估計偏誤問題。圖1 展示了2003—2017 年黃河經(jīng)濟(jì)帶按照是否規(guī)劃臨空經(jīng)濟(jì)區(qū)分組后綠色發(fā)展效率的走勢圖，可以得出臨空經(jīng)濟(jì)城市與非臨空經(jīng)濟(jì)城市的綠色發(fā)展效率均值變動趨勢基本一致，整體趨勢基本滿足共同趨勢假定。后續(xù)穩(wěn)健性檢驗中安慰劑檢驗結(jié)果表明，2009 年以前的年份虛擬變量與處理組的交互項系數(shù)并不顯著，說明2009 年政策發(fā)生之前不同組別間綠色發(fā)展效率走勢差異并無明顯變動，同樣證明共同趨勢假定成立。在隨后估計模型中加入固定效應(yīng)和考慮經(jīng)濟(jì)沖擊后，保證整個估計結(jié)果更加精確。

圖1 黃河經(jīng)濟(jì)帶按臨空經(jīng)濟(jì)分綠色發(fā)展效率走勢圖

圖2 為2005—2017 年黃河經(jīng)濟(jì)帶臨空經(jīng)濟(jì)城市航空客運(yùn)量和貨運(yùn)量的走勢圖。臨空經(jīng)濟(jì)城市的航空貨運(yùn)量和航空客運(yùn)量曲線整體呈凹形增長，在2007—2008年之間增長較為平緩，到2009 年之后整個臨空經(jīng)濟(jì)區(qū)的航空客運(yùn)量和貨運(yùn)量出現(xiàn)較大幅度增長，尤其是在航空貨運(yùn)量上表現(xiàn)得尤為顯著。基本可以了解到2009 年大規(guī)模的基建投資以后，臨空經(jīng)濟(jì)城市的航空客運(yùn)量和貨運(yùn)量不斷攀升，而且增速不斷加快。因此，可以認(rèn)為2009 年的基建投資政策對臨空經(jīng)濟(jì)航空客運(yùn)量和貨運(yùn)量有著一個外生沖擊。

圖2 黃河經(jīng)濟(jì)帶臨空經(jīng)濟(jì)城市航空客運(yùn)量和貨運(yùn)量走勢圖

在對基礎(chǔ)模型（5）估計后，表1 報告了OLS 回歸結(jié)果。第二列展示未控制時間和空間固定效應(yīng)的模型估計結(jié)果。結(jié)果表明，臨空經(jīng)濟(jì)區(qū)能夠促進(jìn)本地城市加快綠色發(fā)展，驗證假說一的成立。考慮到城市間發(fā)展水平、地方政策等方面差異和城市發(fā)展在不同時間段也存在差異。因此，本文控制了空間固定效應(yīng)和時間固定效應(yīng)并再次估計模型，結(jié)果顯示，臨空經(jīng)濟(jì)區(qū)仍然對本地區(qū)域綠色發(fā)展效率有正向影響，支持假說一。但是，以上分析僅能驗證假說一，為了驗證假說二成立的可能性，本文繼續(xù)構(gòu)造并估計基于DID 設(shè)計的空間Durbin 模型。

表1 模型（5）的OLS估計結(jié)果

表2 中展示了模型（6）空間Durbin 模型的估計結(jié)果。加入空間固定效應(yīng)和時間固定效應(yīng)后，并且考慮到固定效應(yīng)會影響參數(shù)的一致性問題后進(jìn)行有偏糾正[37]，可以發(fā)現(xiàn)：臨空經(jīng)濟(jì)對本地區(qū)域綠色發(fā)展依然有顯著促進(jìn)作用，同樣證實了假說一的成立。同時，城市臨空經(jīng)濟(jì)區(qū)發(fā)展水平的提高會抑制周邊城市的綠色發(fā)展，因此證實假說二成立。主要由于臨空經(jīng)濟(jì)城市與鄰接城市產(chǎn)業(yè)間構(gòu)成競爭關(guān)系，臨空經(jīng)濟(jì)城市憑借靠近臨空經(jīng)濟(jì)區(qū)的優(yōu)勢打壓鄰接城市競爭性產(chǎn)業(yè)。具體表現(xiàn)為，臨空經(jīng)濟(jì)區(qū)所在城市吸引城市或者區(qū)域外的高新技術(shù)等綠色產(chǎn)業(yè)加快流入。還有一些發(fā)現(xiàn)，區(qū)域或者城市內(nèi)的產(chǎn)業(yè)集聚會抑制本地綠色發(fā)展，但是不會影響周圍城市的綠色發(fā)展?？紤]到城市的產(chǎn)業(yè)集聚度提高和競爭不斷加大，區(qū)域間企業(yè)壁壘和地方保護(hù)主義增強(qiáng)，削弱了本地產(chǎn)業(yè)集聚的外溢效應(yīng)。綠色發(fā)展效率與人均GDP 成正比，與其平方項成反比，意味著庫茲涅茲曲線假說成立[38]。產(chǎn)業(yè)結(jié)構(gòu)升級中產(chǎn)業(yè)結(jié)構(gòu)高級化會推動地區(qū)綠色發(fā)展，而產(chǎn)業(yè)結(jié)構(gòu)合理化對城市的綠色發(fā)展效率并無太大影響。政府干預(yù)會極大促進(jìn)區(qū)域的綠色發(fā)展，推動地區(qū)環(huán)境向好轉(zhuǎn)變。同時政府干預(yù)的空間溢出效應(yīng)不是太顯著，主要由于我國行政管理上的隔離，某一地級市的政策手段無法干預(yù)其他地級市的生產(chǎn)活動。

表2 模型（6）的空間Durbin模型估計結(jié)果

5.2 穩(wěn)健性檢驗結(jié)果

為了進(jìn)一步研究各個時期可能影響到臨空經(jīng)濟(jì)發(fā)展水平的政策和事件的作用效果，本文將各年份單獨(dú)設(shè)置為虛擬變量放入模型，借此識別各年份政策或事件對臨空經(jīng)濟(jì)的影響效應(yīng)，設(shè)定模型（13）和（14）如下：

表3 中報告了分年份政策或者事件影響效應(yīng)的估計結(jié)果。總的來說，2009 年之前，一些可能性政策或事件對區(qū)域綠色發(fā)展并無顯著影響，而2009 年作為一個時期標(biāo)志，其后年份基本表現(xiàn)出不同程度的正向影響效應(yīng)，這直接說明了本文DID 設(shè)計的穩(wěn)健性。當(dāng)然也說明了2009 年后一些年份（包括2009 年）臨空經(jīng)濟(jì)發(fā)展水平的政策或事件沖擊對本地綠色發(fā)展有著顯著的正向影響，而對鄰居地區(qū)的綠色發(fā)展有明顯的抑制作用，驗證了主要結(jié)果的穩(wěn)健性，可以認(rèn)為假說一和假說二仍然成立。此外，2009 年以前的政策或事件對區(qū)域綠色發(fā)展并無明顯影響。到2009 年及其以后的結(jié)果中，可以發(fā)現(xiàn)政策沖擊對本地綠色發(fā)展效率存在顯著的正向影響?？傮w上參數(shù)估計值依然支持假說一成立。因此，相對于2009 年以前年份事件效應(yīng)的不顯著，2009 年及以后年份顯著的影響效應(yīng)證實空間DID 設(shè)計的穩(wěn)健性。

表3 模型（13）和（14）估計結(jié)果

為了緩解空間異質(zhì)性問題，利用地理加權(quán)回歸（GWR）的方法構(gòu)造模型（10），以驗證臨空經(jīng)濟(jì)區(qū)是否仍對區(qū)域綠色發(fā)展有促進(jìn)作用。由于地理加權(quán)回歸是以帶寬窗口掃描全局空間單元，因此我們僅控制時間固定效應(yīng)。表4 報告了不同帶寬選擇下模型（7）的GWR 估計結(jié)果。在25%、50%、75%分位處，臨空經(jīng)濟(jì)區(qū)對區(qū)域綠色發(fā)展仍保持正向影響，這與主要結(jié)果中的結(jié)論基本一致。

表4 不同帶寬選擇下模型（7）的GWR估計結(jié)果

6 結(jié)論與政策建議

本文通過在空間Durbin 模型基礎(chǔ)上構(gòu)建DID 設(shè)計緩解臨空經(jīng)濟(jì)發(fā)展水平和綠色發(fā)展效率之間的內(nèi)生性問題，使用區(qū)位熵消除產(chǎn)業(yè)集聚指標(biāo)中不同地區(qū)規(guī)模差異的影響，最終研究發(fā)現(xiàn)：首先，臨空經(jīng)濟(jì)會促進(jìn)自身地區(qū)的綠色發(fā)展，并且臨空經(jīng)濟(jì)長期內(nèi)會通過競爭機(jī)制遏制周邊地區(qū)的綠色發(fā)展，主要原因在于未合理配置臨空經(jīng)濟(jì)區(qū)所在城市和周邊區(qū)域的互補(bǔ)性產(chǎn)業(yè)資源，造成區(qū)域間不必要的競爭損失。其次，本文研究還發(fā)現(xiàn)產(chǎn)業(yè)集聚會抑制本地綠色發(fā)展，但是不會影響周圍城市的綠色發(fā)展。產(chǎn)業(yè)結(jié)構(gòu)升級中產(chǎn)業(yè)結(jié)構(gòu)高級化會顯著促進(jìn)自身加快綠色發(fā)展，并且有正的溢出效應(yīng)；而產(chǎn)業(yè)結(jié)構(gòu)合理化對區(qū)域綠色發(fā)展并無太大影響，主要是由于當(dāng)前絕大部分黃河經(jīng)濟(jì)帶地級市處于經(jīng)濟(jì)高速發(fā)展過程中，政府更多強(qiáng)調(diào)第三產(chǎn)業(yè)占比，因此造成當(dāng)?shù)鼐G色發(fā)展效率對產(chǎn)業(yè)結(jié)構(gòu)高級化的變動表現(xiàn)更加敏感，而這些地市政府并未太多關(guān)注產(chǎn)業(yè)結(jié)構(gòu)合理化配置，所以其影響效應(yīng)并不顯著。另外，政府干預(yù)可以有效抑制環(huán)境污染和改善整體生態(tài)水平，推動區(qū)域加速綠色發(fā)展，但由于行政等政治因素，其并沒表現(xiàn)出明顯的空間溢出效應(yīng)。最后，本文對主回歸進(jìn)行穩(wěn)健性檢驗，都支持假說一和假說二。對各年份政策或者事件的影響效應(yīng)單獨(dú)識別后發(fā)現(xiàn)，本文的空間DID 設(shè)計的政策處理時點選取具有穩(wěn)健性，同時臨空經(jīng)濟(jì)區(qū)對本地和周邊區(qū)域綠色發(fā)展水平的影響方向同主要回歸結(jié)果基本一致。另外，在通過GWR 方法緩解空間異質(zhì)性問題后，仍然發(fā)現(xiàn)臨空經(jīng)濟(jì)區(qū)對區(qū)域綠色發(fā)展效率的作用方向與主回歸結(jié)果保持一致。

從政策角度來看，綜合以上研究結(jié)論，本文提出以下政策建議：其一，在建設(shè)臨空經(jīng)濟(jì)以推動區(qū)域綠色發(fā)展過程中，要合理規(guī)劃臨空經(jīng)濟(jì)城市和選擇空間區(qū)位，注意處理同周邊區(qū)域產(chǎn)業(yè)資源配置問題，引導(dǎo)臨空經(jīng)濟(jì)城市向下游高新綠色產(chǎn)業(yè)發(fā)展而周邊地區(qū)向產(chǎn)業(yè)鏈上游供應(yīng)端發(fā)展，中心區(qū)和外圍區(qū)產(chǎn)業(yè)資源相互互補(bǔ)，最終帶動臨空經(jīng)濟(jì)城市輻射區(qū)域整體經(jīng)濟(jì)快速騰飛。其二，在城市化的進(jìn)程中要特別注重產(chǎn)業(yè)集聚帶來的負(fù)面影響，合理布局產(chǎn)業(yè)在城市區(qū)域內(nèi)的空間分布，加快產(chǎn)業(yè)結(jié)構(gòu)高級化進(jìn)度，推動區(qū)域產(chǎn)業(yè)轉(zhuǎn)型升級，避免當(dāng)?shù)刭Y源狀況惡化；另外，要注重破除區(qū)域間企業(yè)壁壘和地方保護(hù)主義的阻礙，整體推進(jìn)經(jīng)濟(jì)帶地市間相互促進(jìn)和協(xié)調(diào)發(fā)展。其三，應(yīng)適當(dāng)加大政府干預(yù)和管控，制定環(huán)境規(guī)制政策，監(jiān)督企業(yè)生產(chǎn)過程和督促污染企業(yè)盡快整改，快速推動整個黃河經(jīng)濟(jì)帶生態(tài)保護(hù)和高質(zhì)量發(fā)展。