滕 麗，盧 君

(廣州大學(xué) 地理科學(xué)學(xué)院，廣東 廣州 510006)

?

珠江三角洲城市空氣質(zhì)量的區(qū)域影響分析

滕 麗，盧 君

(廣州大學(xué) 地理科學(xué)學(xué)院，廣東 廣州 510006)

基于廣東省環(huán)境信息綜合發(fā)布平臺(tái)發(fā)布的珠三角9城市空氣質(zhì)量日?qǐng)?bào)數(shù)據(jù)構(gòu)建VAR模型，通過格蘭杰檢驗(yàn)和脈沖響應(yīng)分析，探討2002～2012年珠三角9城市空氣質(zhì)量的區(qū)域影響特征。結(jié)果顯示：珠三角9城市間空氣質(zhì)量的區(qū)域影響包括自我影響和跨域影響。且城市個(gè)體的自我影響是負(fù)效應(yīng)，其強(qiáng)度遠(yuǎn)遠(yuǎn)超過跨域影響。在跨域影響中單向性影響占比84%，雙向性影響占比16%，鄰域影響占42%，非鄰域影響占比58%。南部沿海城市影響范圍普遍大于北部城市的影響范圍。具有因果關(guān)系的城市間的跨域影響中，65%表現(xiàn)為正效應(yīng)，22%為負(fù)效應(yīng)，有13%顯示為正負(fù)效應(yīng)轉(zhuǎn)換。而正負(fù)效應(yīng)轉(zhuǎn)換的情況說明存在兩種不同的影響模式。

空氣質(zhì)量；API指數(shù)；跨域影響；向量自回歸模型

0 引言

隨著區(qū)域經(jīng)濟(jì)發(fā)展，資源環(huán)境問題對(duì)區(qū)域經(jīng)濟(jì)增長的限制性影響加劇。珠江三角洲是中國自然地理環(huán)境獨(dú)特，人口經(jīng)濟(jì)高度密集的地區(qū)之一。區(qū)域資源開發(fā)強(qiáng)度大，發(fā)展速度快。人地之間相互作用劇烈，資源環(huán)境問題對(duì)該地區(qū)經(jīng)濟(jì)發(fā)展的限制性凸現(xiàn)。其中大氣污染問題是珠三角面臨的主要環(huán)境問題之一。而城市空氣質(zhì)量的區(qū)域影響是建立跨域生態(tài)補(bǔ)償機(jī)制的重要基礎(chǔ)，對(duì)區(qū)域環(huán)境治理政策制定具有重要意義。

城市空氣質(zhì)量的區(qū)域影響是大氣污染物空間擴(kuò)散的結(jié)果。在環(huán)境科學(xué)領(lǐng)域，對(duì)大氣污染擴(kuò)散的研究主要基于大氣動(dòng)力學(xué)理論，通常采用高斯模式及其變形模式，統(tǒng)計(jì)模式、大氣壓擴(kuò)散相似模式或K模式等精確模擬污染物濃度的時(shí)空分布及其變化[1-3]。對(duì)珠三角地區(qū)的實(shí)證研究主要有余緯等采用數(shù)值模擬分析了珠三角一次重空氣污染過程的特征[4]。胡曉宇等利用 Models-3/CMAQ 模擬系統(tǒng)，對(duì)珠江三角洲地區(qū)2006 年10 月的大氣PM10 污染進(jìn)行模擬研究。通過敏感性分析獲得各城市PM10 濃度隨不同污染源削減的變化情況，量化出城市間空氣污染的相互影響，結(jié)果顯示城市間輸送已成為造成珠三角地區(qū)PM10 污染的重要因素[5]。王淑蘭的數(shù)值模擬結(jié)果也顯示珠三角城市之間大氣污染物存在相互輸送的關(guān)系，其中廣州市大氣污染最容易與周邊城市發(fā)生相互作用[6]。這些研究為本文的研究結(jié)果提供了較好的對(duì)照和驗(yàn)證。作為珠三角空氣質(zhì)量相對(duì)較差的城市，廣州空氣質(zhì)量問題也得到了較多討論。如湯惠君分析了廣州市1981～1997年的大氣監(jiān)測(cè)資料，指出廣州大氣污染物的時(shí)空分布特征與其氣象條件、污染源的分布、排放強(qiáng)度的相關(guān)性較大[7]。沈佳芬等利用2001～2004年廣州市空氣自動(dòng)監(jiān)測(cè)點(diǎn)監(jiān)測(cè)的PM10日平均濃度數(shù)據(jù)，分析了廣州市PM10的時(shí)空變化特征。揭示了單一城市多種污染物的時(shí)空分布特征[8]。馬民濤采用傳統(tǒng)的空氣污染負(fù)荷方法揭示廣州市大氣環(huán)境質(zhì)量演變特征與經(jīng)濟(jì)發(fā)展的相關(guān)性[9]。這些研究揭示了在特定時(shí)間段廣州市的大氣污染分布狀況，并對(duì)其形成原因提供了多角度的理解，但是單個(gè)城市的研究往往對(duì)外源影響因素考慮不夠，所以有必要將城市置于區(qū)域的背景中去開展研究。

總之，數(shù)值模擬適合對(duì)短期內(nèi)污染物空間擴(kuò)散的微觀機(jī)制做出解釋，但是數(shù)值模擬一般將空間劃分為若干個(gè)環(huán)境單元，環(huán)境單元與行政單元并非完全耦合，在應(yīng)用到以行政單元界定的跨域環(huán)境管理決策時(shí)往往有一定的局限性[10，11]。在屬地管理的權(quán)限約束下地方政府對(duì)跨域空氣污染治理失效，有必要開展基于行政單元的空氣質(zhì)量的區(qū)域影響研究，為建立城市間的生態(tài)補(bǔ)償機(jī)制和環(huán)境政策決策提供參考。本文基于VAR模型以珠三角9城市為空間單元，量化分析珠三角9城市2002～2012年間空氣質(zhì)量的區(qū)域影響特征。

1 模型原理

向量自回歸(簡(jiǎn)稱VAR)模型是一種不以經(jīng)濟(jì)理論為基礎(chǔ)，而是把系統(tǒng)中的每一個(gè)內(nèi)生變量作為系統(tǒng)中所有內(nèi)生變量的滯后值的函數(shù)來構(gòu)造模型，從而估計(jì)出全部?jī)?nèi)生變量的動(dòng)態(tài)關(guān)系[12]。VAR模型是一種基于數(shù)據(jù)構(gòu)建的統(tǒng)計(jì)學(xué)模型，將每個(gè)進(jìn)入模型的變量都視為平等的。自西姆斯1980年將該模型引入經(jīng)濟(jì)學(xué)之后，模型常常用于預(yù)測(cè)相互聯(lián)系的時(shí)間序列系統(tǒng)和分析隨機(jī)擾動(dòng)對(duì)變量系統(tǒng)的影響[13][14]。VAR(p)模型的基本形式是：

yt=φ1yt-1+…+φpyt-p+εtt=1，2…T

(1)

式中：yt是k維內(nèi)生變量向量，p是滯后階數(shù)，T是樣本個(gè)數(shù)，εt是k維擾動(dòng)列向量。

VAR模型除了能考察變量來自于自身的影響，還能考察來自于其他變量的影響。在理論上，一個(gè)城市的空氣質(zhì)量不僅與其前期空氣質(zhì)量有關(guān)，也可能受到其他城市空氣質(zhì)量的影響。所以將VAR模型運(yùn)用于珠三角9城市的空氣質(zhì)量的區(qū)域影響研究是適合的。

運(yùn)用VAR模型的方法策略(圖1)首先針對(duì)珠三角城市空氣質(zhì)量，建立VAR模型，每個(gè)城市的空氣質(zhì)量都是一個(gè)序列，9個(gè)城市共有9個(gè)序列進(jìn)入模型，最后得到9個(gè)方程。考慮到大氣污染物的滯留時(shí)間有限，模型的滯后階段設(shè)定為2。用最小二乘法運(yùn)行VAR(2)模型并采用AR根圖檢驗(yàn)其穩(wěn)定性，若穩(wěn)定則代表珠三角城市間空氣質(zhì)量存在長期均衡關(guān)系。珠三角城市間空氣質(zhì)量的相互影響不是一種隨機(jī)現(xiàn)象，而是具有統(tǒng)計(jì)規(guī)律的。其次基于VAR模型開展珠三角城市對(duì)的格蘭杰因果檢驗(yàn)。由于格蘭杰因果關(guān)系判斷具有序列平穩(wěn)且具有協(xié)整性的要求，只有通過單位根檢驗(yàn)和協(xié)整檢驗(yàn)的城市對(duì)才進(jìn)一步進(jìn)行格蘭杰因果關(guān)系檢驗(yàn)。最后為了一步判斷城市間空氣質(zhì)量相互影響強(qiáng)度。該方法策略說明了分析的邏輯過程是客觀的。即只對(duì)滿足一定條件的城市對(duì)進(jìn)行空氣質(zhì)量相互影響判斷。因?yàn)榭諝赓|(zhì)量的區(qū)域影響是一種非常復(fù)雜的自然現(xiàn)象，在缺少大氣動(dòng)力機(jī)制基礎(chǔ)的情況下，這樣的方法策略可以防止對(duì)任意兩個(gè)變量進(jìn)行主觀判斷。

圖1 方法策略

2 數(shù)據(jù)來源與處理

數(shù)據(jù)來源于2002～2012年廣東省環(huán)境信息綜合發(fā)布平臺(tái)(www-app.gdepb.gov.cn)發(fā)布的珠江三角9大城市的空氣污染指數(shù)日?qǐng)?bào)數(shù)據(jù)?？諝馕廴局笖?shù)(Air Pollution Index，簡(jiǎn)稱API)是將監(jiān)測(cè)的幾種空氣污染物濃度簡(jiǎn)化成為單一的概念性指數(shù)值形式，并分級(jí)表征空氣污染程度和空氣質(zhì)量狀況的指標(biāo)。直接使用空氣污染指數(shù)來研究城市間空氣質(zhì)量的相互溢出關(guān)系具有不合理性，因?yàn)楦魑廴疚锏臄U(kuò)散條件、介質(zhì)和特征各不相同，基于空氣污染指數(shù)的城市間相互溢出關(guān)系的分析難以從機(jī)制上深入解釋。所以有必要將空氣污染指數(shù)換算成首要污染物濃度。一方面將合成指標(biāo)轉(zhuǎn)換為單一指標(biāo)，消除了機(jī)制解釋的障礙，另一方面將單一指標(biāo)替代合成指標(biāo)所面臨的可代表性問題降到最低。通過對(duì)珠三角個(gè)城市2002～2012年首要污染物，可吸入顆粒物(即PM10)的頻數(shù)統(tǒng)計(jì)，在2002～2012年，珠三角9城市的PM10為首要污染物出現(xiàn)的頻數(shù)均在80%以上，所以用PM10濃度來表征珠三角城市的空氣質(zhì)量較為合適(表1)。

表1 珠江三角洲各城市首要污染物為可吸入顆粒物的頻數(shù)(2002～2012年)

因此通過下式將各城市2002～2012年的空氣污染指數(shù)轉(zhuǎn)換為PM10濃度。為了降低短時(shí)、隨機(jī)影響因素的干擾，同時(shí)考慮PM10在空氣中滯留或擴(kuò)散的時(shí)間較長，所以用于VAR建模的PM10濃度指標(biāo)采用了周平均濃度。

(1)

式中：I為污染物的污染指數(shù)，C為污染物的濃度，C大和C小分別為C所處的濃度段的上下值范圍，I大和I小則為該濃度段對(duì)應(yīng)的API值的上下限。

3 計(jì)算結(jié)果與分析

3.1 VAR模型回歸結(jié)果

從VAR系統(tǒng)整體來看，依據(jù)AIC和SC信息量取值最小原則，模型階數(shù)為2時(shí)，AIC值(-70.684)小于SC值(-66.714)，滿足統(tǒng)計(jì)合理性。9個(gè)子方程的可決系數(shù)R2除了廣州稍低(0.58)以外，其他均在0.84以上，說明各變量所形成的結(jié)構(gòu)性相關(guān)較強(qiáng)(表2)。AR根圖的結(jié)果顯示所有根的模的倒數(shù)均位于單位圓內(nèi)，模型是穩(wěn)定的。珠江三角洲9城市的空間質(zhì)量相互影響關(guān)系是長期穩(wěn)定的(圖2)。

圖2 珠三角9城市VAR模型的AR根圖

結(jié)果城市廣州東莞佛山惠州江門深圳珠海肇慶中山R20.6480.890.8890.870.9190.9090.8650.8740.921可調(diào)R20.5870.8720.8690.8470.9060.8930.8420.8520.907F10.64547.05746.18838.64266.11357.50237.07840.05967.562最大似然354.735476.061475.661468.087472.863484.502485.762462.955485.097AIC-5.15-7.016-7.01-6.894-6.967-7.146-7.166-6.815-7.155SC-4.709-6.575-6.569-6.452-6.526-6.705-6.724-6.374-6.714行列式殘差協(xié)方差(dofadj.)4.57E-43行列式殘差協(xié)方差1.02E-43最大似然4774.464AIC-70.684SC-66.714

3.2 格蘭杰檢驗(yàn)結(jié)果

格蘭杰因果檢驗(yàn)結(jié)果表明在珠三角各個(gè)城市間的空氣質(zhì)量的區(qū)域影響差異比較明顯。空氣質(zhì)量具有雙向因果關(guān)系的3對(duì)城市分別是廣州-惠州，廣州-深圳，深圳-惠州，占比16%。具有單向因果關(guān)系的17對(duì)城市占比84%。所以珠三角地區(qū)城市間空氣質(zhì)量以單向影響為主(表3)。

由于環(huán)境影響是市場(chǎng)失靈的，所以地方政府具有將一些有污染的產(chǎn)業(yè)布局到行政區(qū)邊緣地帶的傾向。接近行政邊界附近的產(chǎn)業(yè)布局可能給鄰近區(qū)域帶來環(huán)境影響。所以從具有因果關(guān)系的城市間是否擁有公共地理邊界的角度觀察，珠三角各城市空氣質(zhì)量鄰域影響的占42%，非鄰域影響的占58%。

從影響范圍的大小看，深圳市的空氣質(zhì)量影響范圍最廣，達(dá)到5個(gè)城市。其次是中山、珠海，其影響的城市分別為4和3個(gè)，其余城市的影響范圍均在1～2個(gè)城市。所以珠三角地區(qū)PM10濃度值較高的廣州、佛山和東莞空氣質(zhì)量的影響有限，而位于沿海的深圳、中山、珠海等城市的空氣質(zhì)量影響范圍較廣。之所以形成這樣的格局，是因?yàn)橹槿谴笾率且粋€(gè)北高南低，開口在南的馬蹄形地貌，南部城市的大氣污染物會(huì)隨著東南風(fēng)或西南風(fēng)向北擴(kuò)散，而廣佛莞北部受山地影響，污染物擴(kuò)散條件較差，影響范圍受限。

表3 空氣質(zhì)量具有格蘭杰因果關(guān)系的城市

注：*城市A是城市B的格蘭杰原因；[]中的數(shù)字為格蘭杰因果檢驗(yàn)的置信度P值；(Y)表示與城市A有公共邊界；(N)表示與城市A沒有公共邊界.

3.3 脈沖響應(yīng)結(jié)果

脈沖響應(yīng)分析就是用來衡量來自隨機(jī)擾動(dòng)項(xiàng)的一個(gè)新息沖擊對(duì)內(nèi)生變量當(dāng)前和未來值之間的動(dòng)態(tài)交互作用及其效應(yīng)。運(yùn)用廣義脈沖響應(yīng)函數(shù)分析空氣質(zhì)量具有因果關(guān)系的城市之間的沖擊響應(yīng)，沖擊響應(yīng)期設(shè)定為20期。脈沖響應(yīng)曲線圖的縱坐標(biāo)分別表示一個(gè)單位沖擊引起的波動(dòng)響應(yīng)值，橫坐標(biāo)表示波動(dòng)持續(xù)的響應(yīng)期。因PM10濃度值是表征空氣質(zhì)量的反向指標(biāo)，值越大空氣質(zhì)量越差，所以當(dāng)波動(dòng)響應(yīng)值為正，表示受影響城市PM10濃度值提高，是一種負(fù)效應(yīng)。相反，波動(dòng)響應(yīng)值為負(fù)，表示受影響城市PM10濃度值降低，是一種正效應(yīng)。

圖3顯示所有城市PM10濃度的變化隨自身沖擊的響應(yīng)均呈非線性函數(shù)的變化趨勢(shì)，大體上經(jīng)過10期之后收斂到0。除了惠州收斂速度較為均勻之外，其他城市在前4期收斂較快，第4期之后收斂較慢。這說明初始階段所有城市的PM10無一例外都受到城市自身的影響，這種影響是一種負(fù)效應(yīng)，也就是一個(gè)新息沖擊會(huì)導(dǎo)致城市空氣質(zhì)量的下降，隨后這種內(nèi)源作用下降，外圍城市的影響逐漸顯現(xiàn)。對(duì)比各城市內(nèi)源作用的程度，內(nèi)源作用最強(qiáng)的是廣州，響應(yīng)系數(shù)初始值為0.17，其他8個(gè)城市的自我響應(yīng)系數(shù)初始值均在0.06左右。

跨域影響有正負(fù)以及正負(fù)轉(zhuǎn)換3種效應(yīng)，在珠三角9城市空氣質(zhì)量具有格蘭杰因果關(guān)系的城市對(duì)中，對(duì)脈沖新息的響應(yīng)出現(xiàn)了15個(gè)正效應(yīng)(65%)，5個(gè)負(fù)效應(yīng)(22%)，3個(gè)正負(fù)效應(yīng)轉(zhuǎn)換的情況(13%)(表4)。據(jù)李婕等研究表明珠三角9城市在2002～2012年單個(gè)城市的空氣質(zhì)量上升[15]?？缬虺鞘虚g普遍具有正向的影響效應(yīng)，說明城市個(gè)體空氣質(zhì)量的改善對(duì)區(qū)域其他城市空氣質(zhì)量的改善具有影響。但從橫向?qū)Ρ瓤矗L期以來，廣州，東莞和佛山是珠三角區(qū)域空氣質(zhì)量相對(duì)較差的地區(qū)，該地區(qū)跨域影響的情況如何？脈沖響應(yīng)分析的結(jié)果顯示該區(qū)域?qū)ν猱a(chǎn)生的是一種正效應(yīng)，說明這3個(gè)城市空氣質(zhì)量的一個(gè)殘差標(biāo)準(zhǔn)差的變化，會(huì)降低相應(yīng)城市的PM10濃度。廣佛莞3個(gè)城市的地形和氣象等不利于大氣污染物的擴(kuò)散，相反該區(qū)域還有吸納外源污染物的可能，從而呈現(xiàn)一種正的影響效果。

圖3 珠三角9城市PM10濃度的自我響應(yīng)

影響城市響應(yīng)城市響應(yīng)系數(shù)初始值響應(yīng)系數(shù)極值持續(xù)期數(shù)實(shí)際效應(yīng)廣州惠州深圳-0.02-0.038-0.0260.0052014正先正后負(fù)佛山東莞0.0120.0128負(fù)東莞江門肇慶-0.012-0.008-0.012-0.012613正正肇慶廣州佛山江門-0.006-0.0005-0.0014-0.0060.0001-0.001451613正先負(fù)后正正中山深圳惠州肇慶廣州0.002-0.00250.0002-0.0060.002-0.00250.0014-0.006916105負(fù)正負(fù)正影響城市響應(yīng)城市響應(yīng)系數(shù)初始值響應(yīng)系數(shù)極值持續(xù)期數(shù)實(shí)際效應(yīng)江門佛山-0.012-0.01212正惠州廣州深圳-0.0050.0002-0.0050.00022013正正深圳珠?；葜輳V州東莞肇慶0.0250.002-0.1-0.00140.00010.025-0.0018-0.1-0.00150.0122071498負(fù)先負(fù)后正正正負(fù)珠海廣州東莞肇慶-0.008-0.0012-0.001-0.008-0.0012-0.001202010正正正

值得特別關(guān)注的是，局部城市對(duì)之間的負(fù)影響效應(yīng)(表4)。具體是東莞對(duì)佛山，深圳對(duì)中山，肇慶對(duì)中山，珠海對(duì)深圳，肇慶對(duì)深圳的脈沖響應(yīng)為負(fù)。這種負(fù)的單向影響往往發(fā)生在有一定距離的兩個(gè)城市而不是相鄰城市。所以珠三角城市空氣質(zhì)量的負(fù)影響并不是嚴(yán)格按照傳染型擴(kuò)散模式而發(fā)生。這說明PM10的遠(yuǎn)程輸送也是存在的。

另外看3對(duì)正負(fù)效應(yīng)轉(zhuǎn)化的城市，其中深圳對(duì)廣州的脈沖響應(yīng)是先正后負(fù)，佛山對(duì)肇慶，惠州對(duì)深圳的響應(yīng)則是先負(fù)后正。這種效應(yīng)轉(zhuǎn)化的情況說明珠三角城市空氣質(zhì)量的影響出現(xiàn)了兩種模式。其中深圳對(duì)廣州的響應(yīng)關(guān)系代表了PM10高濃度值城市對(duì)低濃度值城市的影響(定義為模式I)。相應(yīng)的事實(shí)理解為廣州先吸納外部城市的污染物，然后再對(duì)外擴(kuò)散加重對(duì)應(yīng)城市污染物濃度。所以對(duì)外影響表現(xiàn)為先正后負(fù)。佛山對(duì)肇慶的響應(yīng)關(guān)系代表了PM10低濃度值城市對(duì)高濃度值城市的影響(定義為模式II)，相應(yīng)的事實(shí)理解為肇慶污染物擴(kuò)散在先，造成下風(fēng)向城市污染物濃度上升，隨后以更大的環(huán)境容量稀釋佛山擴(kuò)散出的污染物，表現(xiàn)為先負(fù)后正的效果。

4 結(jié)論

通過珠三角9城市空氣質(zhì)量的VAR模型構(gòu)建和分析，得到的主要結(jié)論有：

(1)2002～2012年可吸入顆粒物(PM10)是珠三角地區(qū)城市空氣污染首要污染物，用PM10濃度表征城市空氣質(zhì)量具有典型意義。采用PM10濃度周平均值作為研究城市空氣質(zhì)量的區(qū)域影響是適宜的。

(2)珠三角9城市在2002～2012年空氣質(zhì)量的動(dòng)態(tài)變化特征是城市個(gè)體的自我影響強(qiáng)度遠(yuǎn)遠(yuǎn)超過跨域影響強(qiáng)度。9城市的自我響應(yīng)是負(fù)效應(yīng)，但城市間具有正，負(fù)或正負(fù)轉(zhuǎn)換三種效應(yīng)。

(3)總體而言，珠三角9城市間空氣質(zhì)量的影響特征是以單向性，非鄰近性影響為主。南部沿海城市影響范圍普遍大于北部城市的影響范圍。

(4)部分PM10濃度高值或低值城市的對(duì)外影響效應(yīng)出現(xiàn)了正負(fù)轉(zhuǎn)化的模式。PM10濃度高值城市向低值城市的影響是先正后負(fù)，PM10濃度低值城市向高值城市的影響是先負(fù)后正。前者揭示了PM10濃度高值城市先吸納外源污染物，然后再擴(kuò)散，而后者則相反，PM10濃度低值城市先擴(kuò)散污染物，然后再吸納外源污染物。結(jié)合兩種模式可以說明，在一定的區(qū)域范圍內(nèi)，城市個(gè)體的空氣質(zhì)量難以獨(dú)善其身。

(5)對(duì)照珠三角城市空氣質(zhì)量相互作用的原有研究，珠三角城市的跨域影響已經(jīng)成為共識(shí)，污染物濃度處于高位的城市傾向于對(duì)外產(chǎn)生不良影響，例如廣州。但是本文的研究卻發(fā)現(xiàn)并不是所有城市間都會(huì)發(fā)生空氣質(zhì)量相互作用，并且廣州的自我影響強(qiáng)度居珠三角各城市之首，對(duì)外的空氣質(zhì)量影響范圍并不大，甚至對(duì)外產(chǎn)生了正的影響效果。之所以出現(xiàn)這樣的出入，是因?yàn)楸疚难芯克\(yùn)用的VAR模型揭示的變量間相互關(guān)系是一種宏觀層面的長期均衡關(guān)系，而基于數(shù)值模擬的變量間相互關(guān)系往往是一種微觀層面的短期關(guān)系。總之，珠三角城市的相對(duì)位置和局部地理?xiàng)l件存在差異，對(duì)珠三角城市空氣質(zhì)量的區(qū)域性影響不能一概而論。

[1]ByunBW，ChingTKS.SciencealgorithmsoftheEPAmodels3communitymultiscaleairquality(CMAQ)modelingsystem[M].UnitedStatesEnvironmentalProtectionAency，1999，EPA.

[2]徐祥德，周秀驥，施曉暉.城市群落大氣污染源影響的空間結(jié)構(gòu)及尺度特征[J].中國科學(xué)D輯：地球科學(xué)，2005，35(增刊Ⅰ)：1-19.

[3]YeBM，JiXL，YangHZ，etal.ConcentrationandchemicalcompositionofPM2.5inShanghaifora1-yearperiod[J].AtmosphericEnvironment，2003，37(4)：499-510.

[4]余緯，羅栩羽，范紹佳，等.珠三角一次重空氣污染過程特征分析及數(shù)值模擬[J].環(huán)境科學(xué)研究，2011，24(6)：645-653.

[5]胡曉宇，李云鵬，李金鳳，等.珠江三角洲城市群PM10的相互影響研究[J].北京大學(xué)學(xué)報(bào)：自然科學(xué)版，2011，47(3)：519-524.

[6]王淑蘭，張遠(yuǎn)航，鐘流舉，等.珠江三角洲城市間空氣污染的相互影響[J].中國環(huán)境科學(xué)，2005，25(2)：133-137.

[7]湯惠君.廣州市大氣污染分布規(guī)律[J].地理研究，2004，23(4)：495-503.

[8]沈家芬，馮建軍，謝春玲，等.廣州市PM10的時(shí)空變化特征分析[J].生態(tài)環(huán)境，2008，17(2)：553-559.

[9]馬民濤，劉潔，廉婕，等.廣州市大氣環(huán)境演變過程及其控制要素的因子分析[J].環(huán)境工程，2013，31(6)：62-66。

[10]FacundoA，MatthewAC，RobertJRE，etal.InSearchofEnvironmentalSpillovers[J].(2009)TheWorldEconomy.doi：10.1111/j.1467-9701，2008，01160.x.，136-163.

[11]HaoJiming，HeKebin，DuanLei，etal.AirpollutionanditscontrolinChina[J].Front.Environ.Sci.Engin.China，2007，1(2)：129-142.

[12]張曉峒.計(jì)量經(jīng)濟(jì)學(xué)[M].天津：南開大學(xué)出版社，2014.

[13]趙進(jìn)文，范繼濤.經(jīng)濟(jì)增長與能源消費(fèi)內(nèi)在依從關(guān)系的實(shí)證研究[J].經(jīng)濟(jì)研究，2007(8)：31-42.

[14]劉海兵，劉麗.工業(yè)化、城鎮(zhèn)化與農(nóng)業(yè)現(xiàn)代化互動(dòng)關(guān)系實(shí)證研究[J].統(tǒng)計(jì)與決策，2014(22)：98-101.

[15]李婕，滕麗.珠三角城市空氣質(zhì)量的時(shí)空變化特征及影響因素[J].城市觀察，2014(5)：85-95.

ANALYSIS ON REGIONAL EFFECTS OF AIR QUALITY INTERACTION ACROSS CITIES IN PEARL RIVER DELTA

TENG Li，LU Jun

(SchoolofGeographicalScience，GuangzhouUniversity，Guangzhou510006，Guangdong，China)

The Interaction of air quality across the cities exists objectively as a consequence of air pollution spread around.It is meaningful to do some research on regional effects in which both cross-border and intra-city effect involved at the level of administrative units for government making the environmental governance policy and building the mechanism of ecological compensation between the different administrative districts.The article builds the VAR model composed of nine cities in the Pearl River Delta，based on the daily data of API of the cities from 2002 to 2012 year，which is from the Guangdong Province Environmental Information Integrated publishing platform.The results of the Granger Causality test and the impulse response analysis show there exist interaction of air quality intra-city and inter-city.The interaction degree of intra-city is generally stronger than that of inter-city.It is seen that unidirectional influence is high up to 84% and the bi-directional effect only accounts for 16%.The cross-border effects of air quality are happened between the two cities no matter whether they share the common geographical boundary or not.However，such effects occur much more in non-adjacent domain than in adjacent domain.The proportion of the former is 58% while of the later is 42%.Much more cities are influenced by the cities on the south of Pearl River Delta，rather than by those on the north of the region.The regional effects of those cities which satisfied the test of Granger Causality are identified to three kinds which positive，negative or transforming between positive and negative effect.and their proportions are 65%，22% 13% in proper sequence.Specially，the transforming between positive and negative effect reveals the two different interaction processes.

air quality；air pollution index；cross-border effect；VAR model

2015-10-11；

2015-11-18．

教育部人文社會(huì)科學(xué)項(xiàng)目：“資源溢出下的區(qū)域經(jīng)濟(jì)增長機(jī)理和對(duì)策”(12YJC790176)資助.

滕麗(1972-)女，云南省大理人，副教授，博士，主要研究方向?yàn)榻?jīng)濟(jì)地理與區(qū)域發(fā)展.

X513

A

1001-7852(2015)06-0001-07