長江經(jīng)濟(jì)帶城鎮(zhèn)化發(fā)展對霧霾污染影響的門檻效應(yīng)研究——基于居民消費(fèi)水平的視角

■彭迪云，劉 暢，周依仿

2013年1月，接連4次嚴(yán)重等級的霧霾侵蝕了中國大部分疆域，在眾多經(jīng)濟(jì)發(fā)達(dá)的一、二線城市中，只有不到1%的城市達(dá)到世界衛(wèi)生組織推薦的最低標(biāo)準(zhǔn)。霧霾污染給人們的生活帶來嚴(yán)重的影響，并危害著公眾的健康。2015年2月28日，柴靜在互聯(lián)網(wǎng)上公開了一部自費(fèi)拍攝的紀(jì)錄片 《穹頂之下》，引起了全民關(guān)于霧霾的廣泛討論。在互聯(lián)網(wǎng)時(shí)代強(qiáng)大的信息傳播作用帶動(dòng)下，無數(shù)人重新審視了霧霾污染的嚴(yán)重性和危害性。因此，盡管現(xiàn)階段我國經(jīng)濟(jì)增長已取得舉世矚目的成績，但這種高速增長背后各種環(huán)境問題日益突出，霧霾污染問題不容忽視。城鎮(zhèn)化是現(xiàn)代化建設(shè)中的重要一環(huán)。對城鎮(zhèn)化和霧霾污染之間的關(guān)系進(jìn)行深入探究，對于正確認(rèn)識城鎮(zhèn)化發(fā)展進(jìn)程中霧霾污染產(chǎn)生的機(jī)理和應(yīng)對策略具有重要的理論價(jià)值和現(xiàn)實(shí)指導(dǎo)意義。

一、相關(guān)研究回顧及評價(jià)

關(guān)于霧霾的研究大多停留在自然科學(xué)階段。朱成章（2013）研究了我國霧霾的形成原因，認(rèn)為煤炭和石油是導(dǎo)致我國霧霾嚴(yán)重的重要原因，提出防治霧霾應(yīng)先解決好化石能源利用中的環(huán)境問題。顧為東（2014）研究認(rèn)為中國工業(yè)等污染和廣大農(nóng)村的土壤、水源嚴(yán)重污染的疊加效應(yīng)在工業(yè)化進(jìn)程中成為中國嚴(yán)重霧霾形成的特殊機(jī)理。呂效譜等（2013）針對我國2013年1月份大范圍霧霾的成因及特點(diǎn)進(jìn)行分析，認(rèn)為霧霾是以特殊氣象條件為主導(dǎo)的機(jī)動(dòng)車尾氣及煤煙型復(fù)合污染引起的大范圍污染現(xiàn)象。然而關(guān)于霧霾的經(jīng)濟(jì)思路研究比較缺乏，對其的定量分析更是罕見。通過閱讀歷史文獻(xiàn)發(fā)現(xiàn)，發(fā)達(dá)國家在發(fā)展過程中也曾飽受霧霾污染的困擾，如倫敦霧霾事件等?？梢哉f，霧霾是后工業(yè)時(shí)代經(jīng)濟(jì)發(fā)展過程中的必經(jīng)階段，霧霾的產(chǎn)生過程有一定的自然規(guī)律。

在城鎮(zhèn)化與霧霾污染的定量關(guān)系方面，國內(nèi)外幾乎罕有研究，相關(guān)文獻(xiàn)主要圍繞城鎮(zhèn)化與碳排放展開。霧霾不等于碳排放，但兩者關(guān)系十分緊密，霧霾治理中的重要一環(huán)就是發(fā)展低碳經(jīng)濟(jì)，節(jié)能減排。因此，碳排放的研究能為本文霧霾的定量研究提供重要借鑒價(jià)值。20世紀(jì)70年代美國生態(tài)學(xué)家Ehrlich等提出IPAT模型以估計(jì)環(huán)境壓力（I），該模型表明人口（P）、人均財(cái)富量（A）以及技術(shù)（T）三者間的相互作用對環(huán)境具有重要影響。PARIKH等（1995）采用83個(gè)發(fā)展中國家的截面數(shù)據(jù)樣本進(jìn)行研究，發(fā)現(xiàn)城鎮(zhèn)化對碳排放具有顯著正影響。YORK等（2003）在截面計(jì)量模型中引入城鎮(zhèn)化二次項(xiàng)，利用1995年137個(gè)國家的數(shù)據(jù)研究了碳排放與城鎮(zhèn)化之間的非線性關(guān)系。隨后的眾多研究表明，截面模型無法對個(gè)體異質(zhì)性進(jìn)行控制，因此學(xué)術(shù)界逐漸轉(zhuǎn)變?yōu)槔妹姘鍞?shù)據(jù)對經(jīng)濟(jì)、環(huán)境的關(guān)系進(jìn)行研究。FAN等（2006）認(rèn)為，收入水平會(huì)對碳排放量與城鎮(zhèn)化之間的關(guān)系造成影響，按照世界銀行的收入分組標(biāo)準(zhǔn)，他將198個(gè)國家的研究樣本分類，發(fā)現(xiàn)在低收入組中，城鎮(zhèn)化對碳排放影響是正向的，其他收入組對城鎮(zhèn)化影響均為負(fù)向?；谶@一思路，文章運(yùn)用面板門檻模型對該議題進(jìn)行研究。關(guān)于“門檻模型”的應(yīng)用，國內(nèi)外學(xué)者近年來進(jìn)行了大量嘗試。如李子豪、劉輝煌（2012）利用門檻面板回歸，以收入和人力資本為門檻變量，檢驗(yàn)了FDI對環(huán)境影響的門檻效應(yīng)；Roberto Duncan（2015）則以時(shí)間為門檻變量，運(yùn)用門檻模型研究了美國經(jīng)常賬戶失衡的現(xiàn)象。

本文基于生態(tài)經(jīng)濟(jì)學(xué)視角，認(rèn)為可能存在某個(gè)經(jīng)濟(jì)拐點(diǎn)，跨越拐點(diǎn)后霧霾污染將會(huì)大幅度減少。為探索拐點(diǎn)的存在，本文選用“門檻模型”進(jìn)行分析。文章以長江經(jīng)濟(jì)帶11個(gè)省市為研究對象，通過借鑒FAN等的研究思路，結(jié)合國內(nèi)現(xiàn)有 “門檻模型”實(shí)證，創(chuàng)新性地研究居民消費(fèi)水平視角下城鎮(zhèn)化發(fā)展對霧霾污染影響的門檻效應(yīng)。在這里，門檻效應(yīng)是指當(dāng)一個(gè)經(jīng)濟(jì)參數(shù)達(dá)到特定的數(shù)值后，引起另外一個(gè)經(jīng)濟(jì)參數(shù)發(fā)生突然轉(zhuǎn)向其他發(fā)展形式的現(xiàn)象，而作為原因現(xiàn)象的臨界值稱為門檻值。

二、模型設(shè)定與估計(jì)方法

（一）模型設(shè)定

為研究城鎮(zhèn)化發(fā)展?fàn)顩r與霧霾污染間的相互關(guān)系，本文以居民消費(fèi)水平作為城鎮(zhèn)化發(fā)展對霧霾污染影響的門檻變量，即考察長江經(jīng)濟(jì)帶各地區(qū)居民消費(fèi)水平是否達(dá)到城鎮(zhèn)化發(fā)展影響霧霾污染的門檻水平，在此基礎(chǔ)上可以通過進(jìn)一步檢驗(yàn)，判斷是否存在多個(gè)門檻變量。傳統(tǒng)門檻效應(yīng)檢驗(yàn)方法分為分組檢驗(yàn)法和交叉項(xiàng)檢驗(yàn)法，但在實(shí)際應(yīng)用過程中，前者具有較強(qiáng)主觀性；后者雖然能夠具體評估門檻值，但苛刻的線性要求無法保證門檻值的正確性。目前學(xué)術(shù)界研究門檻效應(yīng)常用的是Hansen（1999）提出的非動(dòng)態(tài)面板門檻模型，它彌補(bǔ)了上述兩種方法的不足，不僅可以估計(jì)門檻值，而且對門檻值的正確性和“內(nèi)生門檻效應(yīng)”進(jìn)行檢驗(yàn)。基本思想是將門檻值作為位置變量納入實(shí)證模型中，構(gòu)建關(guān)鍵變量的分段函數(shù)，并對門檻值及“門檻效應(yīng)”進(jìn)行一系列估計(jì)和檢驗(yàn)，最終得出存在的“門檻效應(yīng)”。

本研究以生態(tài)經(jīng)濟(jì)領(lǐng)域廣泛使用的STIRPAT模型為基礎(chǔ)，構(gòu)建檢驗(yàn)城鎮(zhèn)化發(fā)展與霧霾污染之間的實(shí)證模型。STIRPAT模型是經(jīng)典IPAT等式改造而成隨機(jī)形式，其基本形式為：I=aPβ1Aβ2Tβ3ε。 式中，I、P、A和T分別為環(huán)境要素、人口要素、財(cái)富要素和技術(shù)要素，a為系數(shù)值，ε為隨機(jī)誤差項(xiàng)。研究霧霾的產(chǎn)生，產(chǎn)業(yè)結(jié)構(gòu)關(guān)系密切，結(jié)合具體應(yīng)用，本文對STIRPAT模型進(jìn)行了拓展，選取第二、三產(chǎn)業(yè)占比和城鎮(zhèn)化水平兩個(gè)因素作為技術(shù)因素的表征。故得到霧霾污染影響因素的基礎(chǔ)數(shù)學(xué)模型：

式中，Iit為霧霾嚴(yán)重性程度；a為常數(shù)；Pit為人口要素的表征值，用“常住人口數(shù)”表示；Ait為財(cái)富要素的表征值，用“人均GDP”表示；Sit為產(chǎn)業(yè)結(jié)構(gòu)狀況，即二、三產(chǎn)業(yè)占比；Uit為城鎮(zhèn)化水平，用城鎮(zhèn)人口占常住人口比值衡量；ε為隨機(jī)擾動(dòng)項(xiàng)。

為通過回歸分析得到相應(yīng)參數(shù)，對（1）式兩邊取對數(shù)，得

模型（2）沒有考慮居民消費(fèi)水平的“門檻效應(yīng)”。借鑒Hansen非動(dòng)態(tài)面板門檻回歸模型思路，假設(shè)存在“單門檻效應(yīng)”，在模型（2）的基礎(chǔ)上構(gòu)建單門檻模型：

其中，居民消費(fèi)水平C為門檻變量，η為帶估算門檻值，L（·）為一指標(biāo)函數(shù)。以上為假設(shè)僅存在“單門檻效應(yīng)”模型構(gòu)建，但實(shí)際應(yīng)用中，可能存在多個(gè)門檻值?，F(xiàn)列舉“雙門檻效應(yīng)”模型進(jìn)行說明，“多門檻效應(yīng)”模型構(gòu)建時(shí)可以此類推：

（二）估計(jì)方法與檢驗(yàn)

為敘述的簡練性，文章將以單門檻模型為例進(jìn)行說明，雙（多）門檻模型估計(jì)方法與檢驗(yàn)方法與此類似。門檻分析的關(guān)鍵是為了解決：一是聯(lián)合估計(jì)門檻值η和變量系數(shù)θ；二是進(jìn)行門檻效應(yīng)的相關(guān)檢驗(yàn)。下面針對這兩個(gè)目標(biāo)進(jìn)行具體說明。

對聯(lián)合門檻值門檻估計(jì)的過程中，選定一個(gè)門檻初始值β0，利用最小二乘法估計(jì)出解釋變量系數(shù)，并求出殘差平方和，最小殘差平方和所對應(yīng)的門檻值η即為所求。為了對模型參數(shù)進(jìn)行估計(jì)，一般常用格點(diǎn)搜索法（Grid Search），求出殘差平方和最小的門檻值η后，利用常規(guī)估計(jì)方法對各變量進(jìn)行參數(shù)估計(jì)。

門檻效應(yīng)的相關(guān)檢驗(yàn)主要包括：（1）門檻效應(yīng)的顯著性水平；（2）門檻估計(jì)值的真實(shí)性檢驗(yàn)。前者原假設(shè)為 H0：θ1=θ2，檢驗(yàn)統(tǒng)計(jì)量為：

其中，S0為原假設(shè)下進(jìn)行參數(shù)估計(jì)后得到的殘差平方和，δ2（η）為備擇假設(shè)下進(jìn)行參數(shù)估計(jì)后得到的殘差平方。原假設(shè)下門檻值η是不確定的，因此統(tǒng)計(jì)量的分布為非標(biāo)準(zhǔn)分布，但通過采用“自抽樣”方法模擬其漸進(jìn)分布，構(gòu)建對應(yīng)P值。第二個(gè)檢驗(yàn)原假設(shè)為H0：η*=η0，相應(yīng)的似然比檢驗(yàn)統(tǒng)計(jì)量為：LR1（η）=（S1（η）-S1（η*））/δ2（η*）。 統(tǒng)計(jì)量 LR1的分布是非標(biāo)準(zhǔn)的，因此Hansen提供了一個(gè)簡單公式計(jì)算拒絕域，即當(dāng)）時(shí)，拒絕原假設(shè)，其中a為顯著性水平。

（三）樣本變量的選取及描述

文章采用2001年至2013年長江經(jīng)濟(jì)帶11個(gè)省市的面板數(shù)據(jù)進(jìn)行實(shí)證研究。數(shù)據(jù)來自歷年《中國統(tǒng)計(jì)年鑒》與各省統(tǒng)計(jì)年鑒，以及《新中國60年統(tǒng)計(jì)資料匯編》、《中國環(huán)境統(tǒng)計(jì)年鑒》。個(gè)別缺失數(shù)據(jù)采用灰色預(yù)測法獲得。

1.被解釋變量的確定

根據(jù)文章討論話題，被解釋變量為霧霾水平。由于缺乏直接數(shù)據(jù)反映霧霾程度，因此，文章運(yùn)用因子分析方法為各地區(qū)霧霾狀況進(jìn)行打分，將得分值作為霧霾水平指標(biāo)。

通過相關(guān)文獻(xiàn)閱讀可以歸納出，造成霧霾的原因可以歸類為工廠廢氣排放、工業(yè)產(chǎn)生的二氧化硫及煙塵、煤炭使用產(chǎn)生的可吸入灰塵、汽車尾氣排放等?；跀?shù)據(jù)可獲得性，將上述指標(biāo)指代如下：x1表示工業(yè)廢氣排放量，x2表示煙塵排放量，x3表示二氧化硫排放量，x4表示煤炭消費(fèi)量，x5表示居民汽車擁有量。所有數(shù)據(jù)經(jīng)過無量綱化處理。

（1）因子分析適合度檢驗(yàn)

文章應(yīng)用SPSS18.0進(jìn)行因子分析。數(shù)據(jù)經(jīng)KMO與巴特利特球形檢驗(yàn)后得到KMO值為0.691，Bartlett球形檢驗(yàn)統(tǒng)計(jì)值為505.453，自由度為10，P值為0.000，因此數(shù)據(jù)適合做因子分析。

（2）公共因子的提取

運(yùn)用SPSS18.0進(jìn)行因子分析，依據(jù)特征值大于1，累計(jì)方差貢獻(xiàn)率大于85%的標(biāo)準(zhǔn)提取公共因子，解釋方差數(shù)據(jù)如表1所示。

表1 解釋的總方差

由表1可以看出前兩個(gè)公共因子的特征值分別為2.756和1.522，均大于1，并且它們的累計(jì)方差貢獻(xiàn)率為85.548%，大于85%。根據(jù)以上分析，我們可以認(rèn)為前兩個(gè)公共因子基本描述了5個(gè)指標(biāo)變量的變化，因此提取的公共因子為2個(gè)。

（3）因子分析結(jié)果

在因子分析時(shí)，以主成分分析法進(jìn)行提取，用方差最大正交旋轉(zhuǎn)法進(jìn)行旋轉(zhuǎn)，所得到的因子載荷矩陣如下表所示：

表2 因子載荷矩陣

通過上表可以看出，工業(yè)廢氣排放量、煤炭消費(fèi)量和汽車擁有量的載荷值分別為0.97、0.930和0.933，計(jì)入公共因子1；煙塵排放量和二氧化硫排放量的載荷值分別為0.866和0.851，計(jì)入公共因子2。

（4）分子得分計(jì)算

運(yùn)用SPSS17.0對提取的兩個(gè)因子進(jìn)行評分及排名，然后將兩個(gè)公共因子的特征值代入公式Wi=，得到W1和W2分別為0.6442和0.3557。根據(jù)公式WiGi(其中Gi為因子得分)計(jì)算因子綜合得分。

2.解釋變量的描述

單門檻模型的形式如（3）所示，由于數(shù)據(jù)進(jìn)行對數(shù)處理，則無需再進(jìn)行標(biāo)準(zhǔn)化處理。需要特別指出的是，基于因子分析綜合得分的計(jì)算方式，所得得分區(qū)間為（-2，2），為保證對數(shù)處理的可行性，霧霾得分?jǐn)?shù)據(jù)均平移2個(gè)單位，但這并不影響霧霾得分的趨勢，因此對文章計(jì)算結(jié)果不會(huì)產(chǎn)生影響。

居民消費(fèi)水平指數(shù)反映不同時(shí)期居民消費(fèi)水平變動(dòng)程度，是一個(gè)國家或地區(qū)居民的生活狀況的表現(xiàn)，可以體現(xiàn)一個(gè)國家或地區(qū)經(jīng)濟(jì)建設(shè)的成果。因此，模型中的門檻變量設(shè)定為居民消費(fèi)水平C。樣本的描述性統(tǒng)計(jì)量如下（表3）。

表3 描述性統(tǒng)計(jì)量

三、實(shí)證結(jié)果及分析

（一）門檻效應(yīng)的顯著性檢驗(yàn)及置信區(qū)間估計(jì)

本文利用stata12.0軟件對門檻模型進(jìn)行估計(jì)和檢驗(yàn)。以居民消費(fèi)水平作為門檻變量建立門檻回歸模型，為確定門檻的個(gè)數(shù)，分別在單門檻、雙門檻和三門檻假設(shè)下分析城鎮(zhèn)化水平的門檻效應(yīng)。表4為居民消費(fèi)水平門檻的顯著性檢驗(yàn)、門檻估計(jì)值及其置信區(qū)間，從表中可以看出，在1%的顯著性水平上城鎮(zhèn)化水平對霧霾存在單門檻效應(yīng)，在5%的顯著性水平下存在雙重門檻效應(yīng)，而三門檻檢驗(yàn)只在10%下顯著，本文不予考慮。因此，我們將討論城鎮(zhèn)化水平與霧霾之間的雙門檻效應(yīng)，兩個(gè)門檻值分別為8190和15895。同時(shí)，門檻1估計(jì)值處于[7683,8700]區(qū)間內(nèi)且門檻2估計(jì)值處于 [15500,17000]區(qū)間內(nèi)時(shí)，似然比值小于5%顯著性水平下的臨界值，在原假設(shè)接受域內(nèi)，即兩個(gè)門檻值都與實(shí)際門檻值相等。

表4 門檻變量的顯著性檢驗(yàn)和置信區(qū)間

（二）門檻模型估計(jì)結(jié)果及分析

1.回歸結(jié)果與分析

確定了門檻值后，對門檻回歸模型進(jìn)行參數(shù)估計(jì)結(jié)果見表5?；貧w表明，各解釋標(biāo)量對霧霾都會(huì)正的作用，其中人口和二三產(chǎn)業(yè)占比對霧霾產(chǎn)生比較大的作用，這與實(shí)際相符。因?yàn)槿丝谶^多會(huì)引發(fā)一系列環(huán)境問題，而二三產(chǎn)業(yè)的發(fā)展，尤其是第二產(chǎn)業(yè)是造成環(huán)境破壞的主要原因。而通過前文分析，城鎮(zhèn)化水平對霧霾的影響存在著居民消費(fèi)水平門檻效應(yīng)，具體分析如下：

城鎮(zhèn)化發(fā)展對長江經(jīng)濟(jì)帶霧霾程度的雙門檻模型中，門檻值分別為8190和15895。從表5的參數(shù)估計(jì)結(jié)果看，當(dāng)長江經(jīng)濟(jì)帶的居民消費(fèi)水平達(dá)到8190元時(shí)，城鎮(zhèn)化發(fā)展對霧霾程度的影響有正向作用，其影響系數(shù)為0.1054,；當(dāng)跨越第一門檻值后，即居民消費(fèi)水平超過8190元時(shí)，城鎮(zhèn)化發(fā)展對霧霾的影響依然為正向作用，但影響系數(shù)降到0.06；當(dāng)跨越第二門檻值時(shí)，即居民消費(fèi)水平超過15895元時(shí)，城鎮(zhèn)化發(fā)展對霧霾的影響轉(zhuǎn)變?yōu)樨?fù)向作用，影響系數(shù)為-0.2407，即此時(shí)的城鎮(zhèn)化發(fā)展有利于減輕霧霾污染。這種現(xiàn)象的原因可能是，經(jīng)濟(jì)發(fā)展初期，居民消費(fèi)水平較低，城鎮(zhèn)化建設(shè)需要大量的基礎(chǔ)設(shè)施建設(shè)，以及外來人口涌入等人口作用，由此產(chǎn)生了巨大的環(huán)境破壞；當(dāng)消費(fèi)水平進(jìn)入一個(gè)階段，即邁過第一門檻值時(shí)，基礎(chǔ)設(shè)施等逐步完善，城市開始由經(jīng)濟(jì)發(fā)展重心向綠色發(fā)展轉(zhuǎn)移，城市建設(shè)開始注重環(huán)境保護(hù)，由此雖然依然有外來人口遷入等一系列問題，但環(huán)境的破壞程度開始減弱；當(dāng)消費(fèi)水平邁過第二門檻值時(shí)，經(jīng)濟(jì)已高度發(fā)達(dá)，城市建設(shè)更加注重人與自然和諧統(tǒng)一，此時(shí)城鎮(zhèn)化建設(shè)中將對環(huán)境造成破壞的項(xiàng)目進(jìn)行拆除，對環(huán)境保護(hù)、城市綠化等更加注重，由此霧霾程度也逐步下降。

表5 門檻模型參數(shù)的估計(jì)結(jié)果

2.拓展分析

長江經(jīng)濟(jì)帶城鎮(zhèn)化發(fā)展霧霾程度影響的居民消費(fèi)水平門檻分析表明，居民消費(fèi)水平影響城鎮(zhèn)化發(fā)展對霧霾程度的作用有正、負(fù)兩種，正向又分為高水平和低水平兩類。居民消費(fèi)水平較低時(shí)，城鎮(zhèn)化發(fā)展引起較高水平的霧霾，隨著居民消費(fèi)水平提升，城鎮(zhèn)化發(fā)展對霧霾的作用逐漸減少，當(dāng)居民消費(fèi)水平足夠高時(shí)，城鎮(zhèn)化發(fā)展與霧霾呈負(fù)向關(guān)系。將研究期間各省市消費(fèi)水平與門檻值進(jìn)行比較，將11個(gè)城市劃分為三類。即居民消費(fèi)低水平區(qū)域（C≤8190），居民消費(fèi)中水平區(qū)域（8190〈C≤15895）和居民消費(fèi)高水平區(qū)域（C〉15895）。根據(jù)所求門檻值及各地區(qū)歷年居民消費(fèi)水平狀況，畫出2001年至2013年長江經(jīng)濟(jì)帶11個(gè)省市居民消費(fèi)水平走勢與門檻值比較圖（圖1），同時(shí)得到各年份不同區(qū)域省份個(gè)數(shù)分布與狀況（表5）及2013年各省市分布情況（表6）。

圖1 長江經(jīng)濟(jì)帶各省市歷年消費(fèi)水平與門檻值比較圖

表6 2001～2013年長江經(jīng)濟(jì)帶居民消費(fèi)水平區(qū)域省份個(gè)數(shù)分布表

表7 2013年居民消費(fèi)情況表

通過上文分析可以看出，進(jìn)入2013年以來，所有省市均跨過了第一門檻值，其中上海、江蘇、浙江已跨過第二門檻值。對比2013年霧霾情況來看，上海、江蘇、浙江霧霾最嚴(yán)重，而中部地區(qū)較差，云貴地區(qū)稍好?？缭降谝婚T檻值后，霧霾的產(chǎn)生速度會(huì)減緩，但城鎮(zhèn)化發(fā)展依舊與霧霾產(chǎn)生成正比，因此東部地區(qū)在跨越第二門檻值時(shí)，處于霧霾程度峰值，因此霧霾情況最為嚴(yán)重，但預(yù)計(jì)未來情況會(huì)有所好轉(zhuǎn)；而中西部地區(qū)僅跨越第一門檻值，霧霾產(chǎn)生速度會(huì)有所減緩，但是隨著時(shí)間推移，預(yù)計(jì)未來霧霾狀況還會(huì)有所嚴(yán)重。同時(shí)，霧霾屬于大氣污染，而大氣處于一個(gè)循環(huán)流動(dòng)過程中，東部地區(qū)自身產(chǎn)生霧霾在未來會(huì)逐步減少，但仍不能排除中西部地區(qū)產(chǎn)生的污染物進(jìn)入隨季風(fēng)等原因進(jìn)入東部地區(qū)。根據(jù)文章結(jié)論，初步繪制出霧霾與城鎮(zhèn)化發(fā)展的消費(fèi)水平門檻效應(yīng)圖（圖 2）。

圖2 霧霾與城鎮(zhèn)化發(fā)展的居民消費(fèi)水平門檻效應(yīng)圖

四、結(jié)論與政策建議

文章基于“門檻回歸”技術(shù)為代表的非線性計(jì)量經(jīng)濟(jì)學(xué)理論，采用長江經(jīng)濟(jì)帶2001～2013年的省級面板數(shù)據(jù)構(gòu)建門檻回歸模型，并以居民消費(fèi)水平作為門檻變量，實(shí)證研究了城鎮(zhèn)化發(fā)展與霧霾污染之間的關(guān)系。對城鎮(zhèn)化門檻效應(yīng)的檢驗(yàn)表明，在居民消費(fèi)水平較低時(shí)，城鎮(zhèn)化發(fā)展對霧霾污染有較大的正相關(guān)，當(dāng)居民消費(fèi)水平跨越第一道門檻時(shí)，城鎮(zhèn)化發(fā)展對霧霾污染的作用開始減弱，但仍呈現(xiàn)正相關(guān)關(guān)系，當(dāng)居民消費(fèi)水平跨越第二道門檻時(shí)，城鎮(zhèn)化發(fā)展與霧霾污染開始呈現(xiàn)負(fù)相關(guān)關(guān)系，城鎮(zhèn)化發(fā)展有利于霧霾的減少。由于第二道門檻之前，城鎮(zhèn)化發(fā)展與霧霾污染均呈現(xiàn)正相關(guān)關(guān)系，因此在跨越第二道門檻時(shí)，霧霾污染程度達(dá)到峰值，這與近些年來東部地區(qū)污染普遍較中西部地區(qū)嚴(yán)重的現(xiàn)象不謀而合。預(yù)測在未來，隨著城鎮(zhèn)化發(fā)展，居民消費(fèi)水平的不斷提高，霧霾污染將會(huì)逐步減少。因此，可以認(rèn)為霧霾是城鎮(zhèn)化發(fā)展過程中必經(jīng)的現(xiàn)象。對比發(fā)達(dá)國家發(fā)展經(jīng)驗(yàn)，這一結(jié)論也與發(fā)達(dá)國家發(fā)展歷史較為一致。然而，由于霧霾屬于大氣污染，流動(dòng)性強(qiáng)，影響因素涉及自然、經(jīng)濟(jì)、社會(huì)等方方面面，因此僅從經(jīng)濟(jì)學(xué)視角進(jìn)行探討具有一定局限性，研究結(jié)論還需要綜合多方學(xué)者共同探討。

打造 “長江經(jīng)濟(jì)帶”是我國未來發(fā)展的重要戰(zhàn)略，研究長江經(jīng)濟(jì)帶城鎮(zhèn)化發(fā)展對霧霾污染的影響極具現(xiàn)實(shí)意義，有助于幫助制定相關(guān)政策，實(shí)現(xiàn)長江經(jīng)濟(jì)帶乃至全國范圍內(nèi)的綠色發(fā)展?；诖?，文章提出以下建議。

首先，加強(qiáng)區(qū)域間合作交流，建立區(qū)域信息互通機(jī)制，提高居民消費(fèi)水平。區(qū)域間合作有利于區(qū)域經(jīng)濟(jì)一體化的發(fā)展，幫助區(qū)域協(xié)調(diào)發(fā)展。長江經(jīng)濟(jì)帶內(nèi)各地區(qū)加強(qiáng)合作，減少重復(fù)建設(shè)，通過信息互通機(jī)制，為區(qū)域發(fā)展?fàn)I造良好共建的氛圍，幫助區(qū)域經(jīng)濟(jì)健康平穩(wěn)發(fā)展。拉動(dòng)內(nèi)需，增加消費(fèi)，提高居民消費(fèi)水平，助推區(qū)域快速邁過第二門檻值，實(shí)現(xiàn)霧霾污染速度下降。

其次，轉(zhuǎn)變東部地區(qū)經(jīng)濟(jì)發(fā)展方式，提高中西部地區(qū)的招商引資能力。東部經(jīng)濟(jì)增長在過去主要靠投資和出口拉動(dòng)，為了保障經(jīng)濟(jì)持續(xù)穩(wěn)健增長，必須重視消費(fèi)主導(dǎo)的經(jīng)濟(jì)模式，充分利用優(yōu)惠政策及配套管理措施，以加快地區(qū)消費(fèi)的最優(yōu)化進(jìn)程,使經(jīng)濟(jì)增長模式逐步向消費(fèi)為主轉(zhuǎn)變，以使消費(fèi)達(dá)到較高水平，減弱城鎮(zhèn)化對霧霾污染的助推作用。而對于中部地區(qū)，應(yīng)進(jìn)一步提高招商引資能力，并改善投資環(huán)境，加強(qiáng)產(chǎn)業(yè)建設(shè)，并提高購居民消費(fèi)買力，使消費(fèi)拉動(dòng)經(jīng)濟(jì)增長，并達(dá)到較高水平以減弱城鎮(zhèn)化對霧霾的助推作用。

第三，促進(jìn)城鎮(zhèn)化綠色發(fā)展，實(shí)現(xiàn)新型城鎮(zhèn)化。近年來，城鎮(zhèn)化發(fā)展越來越注重生態(tài)優(yōu)先、綠色發(fā)展，在不破壞生態(tài)的前提下實(shí)現(xiàn)城鎮(zhèn)化。完善企業(yè)監(jiān)督管理機(jī)制，加強(qiáng)對產(chǎn)業(yè)項(xiàng)目的事前問責(zé)、事后監(jiān)督，對于影響生態(tài)的企業(yè)環(huán)保部門應(yīng)被賦予更大的權(quán)力進(jìn)行監(jiān)管。在城鎮(zhèn)化過程中，非農(nóng)人員的流入要做好安置，基礎(chǔ)設(shè)施配套完善。另一方面，加大產(chǎn)業(yè)集群，促進(jìn)產(chǎn)業(yè)結(jié)構(gòu)的轉(zhuǎn)型升級，轉(zhuǎn)變經(jīng)濟(jì)發(fā)展方式，使經(jīng)濟(jì)體綠色化、效益化、生態(tài)化。

第四，針對霧霾產(chǎn)生的原因進(jìn)行系統(tǒng)性防范，綜合治理環(huán)境污染。優(yōu)化能源結(jié)構(gòu)，減少煤炭等能源的消費(fèi)，提高天然氣普及程度，有管制的放開能源市場，實(shí)現(xiàn)能源市場化。對工業(yè)企業(yè)進(jìn)行24小時(shí)排污監(jiān)控，通過設(shè)立排污費(fèi)等方式減少工業(yè)廢氣等有害氣體的排放。另外，鑒于各地區(qū)霧霾污染程度與城鎮(zhèn)化發(fā)展的差異，可以實(shí)行區(qū)域生活廢氣排放平衡與監(jiān)督機(jī)制，將廢氣減放目標(biāo)分配給各個(gè)區(qū)域，鼓勵(lì)各區(qū)域自行探索自身城鎮(zhèn)化建設(shè)中的霧霾形成原因。

最后，加快實(shí)施國家主體功能區(qū)戰(zhàn)略，大力推進(jìn)生態(tài)文明建設(shè)，形成有序開發(fā)格局。為此，必須構(gòu)建完善的利益補(bǔ)償機(jī)制、科學(xué)高效的考核評價(jià)機(jī)制和統(tǒng)一的監(jiān)管機(jī)構(gòu)，以及長效的地方政府聯(lián)席和對話機(jī)制，以有效兼顧發(fā)展與民生，實(shí)現(xiàn)區(qū)域協(xié)調(diào)發(fā)展和共同富裕。

總之，我國正處于工業(yè)化和新型城鎮(zhèn)化建設(shè)的重要階段，發(fā)展的同時(shí)環(huán)境的保護(hù)也越發(fā)不容忽視。目前，我國大部分省市還沒有跨越第二道門檻值，在當(dāng)前及以后很長一段時(shí)間內(nèi)，我國城鎮(zhèn)化的發(fā)展還將伴隨著霧霾污染的困擾，但隨著城鎮(zhèn)化發(fā)展質(zhì)量提升階段的到來，隨著居民生活及企業(yè)生產(chǎn)資源更加合理、高效地利用，所產(chǎn)生的霧霾污染物會(huì)越來越少，最終將實(shí)現(xiàn)經(jīng)濟(jì)社會(huì)的可持續(xù)發(fā)展。

[1]朱成章.我國防止霧霾污染的對策與建議[J].中外能源,2013，(6):1～4.

[2]顧為東.中國霧霾特殊形成機(jī)理研究[J].宏觀經(jīng)濟(jì)研究,2014,(6):3～7.

[3]呂效譜,成海容,王祖武,張帆.中國大范圍霧霾期間大氣污染特征分析[J].湖南科技大學(xué)學(xué)報(bào)(自然科學(xué)版),2013,(03):104～110.

[4]Ehrlich P R,Holden J P.Impact of population growth[J].Science,1971,(171):1212～1217.

[5]Ehrlich P R,Holden J P.One dimensional economy[J].Bulletin of Atomic Scientists,1972,(16):18～27.

[6]PARIKH J,SHUKLA V.Urbanization,Energy Use and Greenhouse Effects in Economic Development:Result from a Cross-National Study of Developing Countries[J].Global Environmental Change,1995,5(2):87～103.

[7]YORK R,ROSA E A,DIETZ T.STIRPAT,IPAT and Impacts:Analytic Tools for Unpacking the Driving Forces of Environmental Impacts[J].Ecological Economics,2003,46(3):351-365.

[8]FAN Y,LIU L C,WU G,et al.Analyzing Impact Factors of CO2 Emissions Using the STIRPAT Model[J].Environmental Impact Assessment Review,2006,26(4)：377～395.

[9]李子豪,劉輝煌.FDI對環(huán)境的影響存在門檻效應(yīng)嗎——基于中國220個(gè)城市的檢驗(yàn)[J].財(cái)貿(mào)經(jīng)濟(jì),2012,09:101～108.

[10]Duncan R.A threshold model of the US current account[J].Economic Modeling,2014.

[11]李梅,柳士昌.對外直接投資逆向技術(shù)溢出的地區(qū)差異和門檻效應(yīng)——基于中國省際面板數(shù)據(jù)的門檻回歸分析[J].管理世界,2012(1):21～32.

[12]Hansen B.E.,1999,”Threshold Effects in Non-dynamic Panels:Estimation,Testing and Inference”,Journal of Econometrics,Vol.93,pp.345～368.

[13]DIETZ T,ROSA E A.Rethinking the Environmental Impacts of Population,Affluence and Technology[J].Human Ecology Review,1994,1(2):277～300.

[14]劉敏,曹衷陽.外商直接投資對經(jīng)濟(jì)發(fā)展影響的門檻效應(yīng)研究——基于居民相對消費(fèi)水平視角[J].工業(yè)技術(shù)經(jīng)濟(jì),2011,12:134～142.