江曲圖,何俊昱,王占山,葉觀瓊*,陳 倩,肖 璐

(1.浙江大學(xué)海島海岸研究所,浙江 舟山 316021；2.北京市環(huán)境保護(hù)監(jiān)測(cè)中心,北京 100048；3.江西師范大學(xué)地理與環(huán)境學(xué)院,江西 南昌 330022)

基于LUR/BME的海岸帶地區(qū)PM2.5時(shí)空特性研究

江曲圖1,何俊昱1,王占山2,葉觀瓊1*,陳 倩3,肖 璐1

(1.浙江大學(xué)海島海岸研究所,浙江 舟山 316021；2.北京市環(huán)境保護(hù)監(jiān)測(cè)中心,北京 100048；3.江西師范大學(xué)地理與環(huán)境學(xué)院,江西 南昌 330022)

以高程、距海距離、道路、歸一化植被指數(shù)構(gòu)建大區(qū)域土地利用回歸模型LUR,并結(jié)合貝葉斯最大熵BME對(duì)LUR模型的殘差進(jìn)行時(shí)空分析,得到2015年中國(guó)沿海部分省市PM2.5的時(shí)空分布.交叉驗(yàn)證結(jié)果表明在引入BME模型后R2由0.36提高至0.85,均方根誤差RMSE由23.53μg/m3降低至 11.08μg/m3;整體海岸帶地區(qū)以長(zhǎng)江三角洲為界PM2.5濃度呈現(xiàn)南低北高,且以京津冀及山東內(nèi)陸區(qū)域秋冬季污染最為嚴(yán)重,同時(shí)以山東省為例進(jìn)行各市室外人口空氣污染暴露分析,表明沿內(nèi)陸至近海,人均PM2.5暴露濃度逐步遞減,以濟(jì)南85.5μg/m3最高,沿海區(qū)域煙臺(tái),威海等地較低.

土地利用回歸；貝葉斯最大熵；PM2.5；人體暴露

隨著我國(guó)城市化進(jìn)程的加快,各地區(qū)工業(yè)發(fā)展迅速,工業(yè)排放及汽車保有量逐年增加,導(dǎo)致我國(guó)部分地區(qū)PM2.5空氣污染嚴(yán)重,灰霾天氣明顯增加[1-2].而我國(guó)自2012年底才開展PM2.5全面監(jiān)測(cè),其覆蓋點(diǎn)位稀疏且設(shè)備昂貴,如何獲得高時(shí)空分辨率的PM2.5分布對(duì)于環(huán)境管理,空間流行病研究等意義重大[3-4].

土地利用回歸模型(LUR)利用觀測(cè)站點(diǎn)周圍的多種自然環(huán)境因子進(jìn)行線性回歸建模,得到未觀測(cè)區(qū)域的空氣污染物濃度,其作為一種簡(jiǎn)單易操作的數(shù)學(xué)模型在污染物濃度空間制圖方面廣為使用[2,5-8].國(guó)外利用LUR模型對(duì)空氣污染進(jìn)行空間制圖已在多倫多(NO2)[9],加利福利亞(PM2.5)[10],德黑蘭(SO2)[11]等眾多城市得到實(shí)踐.國(guó)內(nèi)[12]首次在2009年采用了不同緩沖區(qū)道路長(zhǎng)度,土地利用,人口密度等構(gòu)建LUR模型對(duì)天津市PM2.5及NO2進(jìn)行了空間分布制圖.隨后相關(guān)學(xué)者分別在濟(jì)南(SO2,NO2,PM10)[13],北京(PM2.5)[6],長(zhǎng)沙(NO2,PM10)[14],上海(NO2,PM2.5)[15]等城市構(gòu)建LUR模型對(duì)空氣污染空間分布制圖.但LUR模型存在空間遷移能力較差,適用于小尺度研究區(qū)域及無(wú)法充分考慮時(shí)間變異等缺陷[16-19].

現(xiàn)代時(shí)空地統(tǒng)計(jì)學(xué)貝葉斯最大熵(BME)作為一種嚴(yán)謹(jǐn)科學(xué)的時(shí)空數(shù)據(jù)分析方法則能很好的綜合各種數(shù)據(jù)來(lái)源達(dá)到更精細(xì)化的時(shí)空制圖.傳統(tǒng)空間插值方法包括泰森多邊形法,反距離權(quán)重法(IDW),克里金(kriging)等.而其中應(yīng)用最廣效果最好的則是克里金插值方法,然而克里金一般限于對(duì)實(shí)測(cè)數(shù)據(jù)空間尺度上的插值,忽略了其時(shí)間序列上的相關(guān)性,且假設(shè)數(shù)據(jù)無(wú)不確定性,正態(tài)分布等,這些都限制了數(shù)據(jù)利用的多樣性及制圖的精確性.BME則能考慮到數(shù)據(jù)的不確定性及時(shí)間上的關(guān)聯(lián),與傳統(tǒng)地統(tǒng)計(jì)學(xué)相比其是非線性插值,不要求人為的正態(tài)分布且能與各種來(lái)源數(shù)據(jù)及方法相結(jié)合,通過(guò)對(duì)臺(tái)北地區(qū)PM2.5分布制圖LUR模型中引入BME,平均誤差從2.78降低至2.15[17];此外,通過(guò)在整個(gè)美國(guó)PM2.5時(shí)空特性研究中對(duì)LUR輔助變量是否加入遙感數(shù)據(jù)及其與BME結(jié)合方法的比較,結(jié)果表明在引入BME后,均方根誤差(RMSE)由15.91降低至4.63[20].

我國(guó)目前利用LUR模型對(duì)PM2.5分布制圖研究仍限于中小城市,且較少考慮時(shí)間變異,對(duì)于預(yù)測(cè)大范圍區(qū)域的PM2.5質(zhì)量濃度分布精度不高.本文以我國(guó)部分海岸帶地區(qū)為研究區(qū)域,采用貝葉斯最大熵BME時(shí)空數(shù)據(jù)分析方法,通過(guò)易獲取的數(shù)據(jù)源,結(jié)合LUR模型殘差進(jìn)行時(shí)空分析,對(duì)中國(guó)沿海部分省市PM2.5時(shí)空分布進(jìn)行研究,得到高精度的格點(diǎn)分布數(shù)據(jù),并以山東省為例,對(duì)各市室外PM2.5人體暴露進(jìn)行分析對(duì)比,進(jìn)一步驗(yàn)證該方法的應(yīng)用價(jià)值.

1 數(shù)據(jù)與方法

1.1 研究區(qū)域及數(shù)據(jù)

本研究選取我國(guó)部分沿海省市,包括北京、天津、河北、山東、江蘇、上海、浙江、福建、廣東由北至南沿海的3個(gè)直轄市和6個(gè)省份(圖1),縱跨渤海,黃海,東海及南海,包括珠江三角洲,長(zhǎng)江三角洲及京津冀,這幾個(gè)區(qū)域是我國(guó)人口密度極高,工商業(yè)最為發(fā)達(dá)的地區(qū),其環(huán)境健康問(wèn)題備受關(guān)注,研究該區(qū)域海岸帶城市PM2.5時(shí)空變化對(duì)于指導(dǎo)環(huán)境治理及開展我國(guó)環(huán)境健康與人體暴露等研究意義重大.

圖1 研究區(qū)域內(nèi)PM2.5監(jiān)測(cè)站點(diǎn)空間分布Fig.1 Spatial distribution of PM2.5monitoring sites in this study

研究數(shù)據(jù)包括:(1)2015年P(guān)M2.5月平均數(shù)據(jù),來(lái)源于中國(guó)環(huán)境監(jiān)測(cè)總站網(wǎng)上發(fā)布平臺(tái)(http://www.cnemc.cn),共448個(gè)國(guó)家站點(diǎn)PM2.5質(zhì)量濃度小時(shí)觀測(cè)值(部分站點(diǎn)數(shù)據(jù)缺失已移除),以日平均值算得448個(gè)站點(diǎn)5376個(gè)月平均PM2.5質(zhì)量濃度觀測(cè)值,由箱線圖(圖2)可知各站點(diǎn)觀測(cè)數(shù)值展現(xiàn)了較大的空間變異及時(shí)間變異特性;(2)DEM高程數(shù)據(jù),來(lái)源于SRTM全球90m分辨率高程數(shù)據(jù)(http://srtm.csi.cgiar.org);(3)歸一化植被指數(shù)NDVI數(shù)據(jù),從中國(guó)科學(xué)院計(jì)算機(jī)網(wǎng)絡(luò)信息中心國(guó)際科學(xué)數(shù)據(jù)鏡像網(wǎng)站(http://www.gscloud.cn)上下載全國(guó)2015年各月NDVI;(4)研究區(qū)域道路矢量圖來(lái)源于Wiki世界地圖數(shù)據(jù)庫(kù)(http://www.openstreetmap.org);(5)人口密度數(shù)據(jù)來(lái)源于哥倫比亞大學(xué)社會(huì)經(jīng)濟(jì)數(shù)據(jù)應(yīng)用中心網(wǎng)站(http://beta.sedac.ciesin.columbia. edu);(6)研究區(qū)域沿海各市煙粉塵排放量及汽車擁有量,數(shù)據(jù)來(lái)源于我國(guó)國(guó)家統(tǒng)計(jì)局發(fā)布的最新統(tǒng)計(jì)年鑒及沿海各省統(tǒng)計(jì)年鑒(http://www. stats.gov.cn).以上地圖數(shù)據(jù)處理均基于ArcGis10.2,采用A lbers等面積投影.

圖2 2015年1～12月各站點(diǎn)PM2.5觀測(cè)濃度Fig.2 Box plot of monthly average concentration PM2.5(fromJan. to Dec. 2015)

1.2 LUR模型

LUR模型的構(gòu)建主要包括在ArcGis10.2軟件中各種變量(variable,以下因子簡(jiǎn)稱v)數(shù)據(jù)的提取和在SPSS22.0軟件中采用逐步多元回歸對(duì)因子進(jìn)行篩選.變量的獲取過(guò)程如下,(1)監(jiān)測(cè)站點(diǎn)處高程值(v1)由ArcGis10.2軟件根據(jù)站點(diǎn)經(jīng)緯度提取高程影像值得到;(2)根據(jù)海岸線及道路矢量圖層計(jì)算得到監(jiān)測(cè)站點(diǎn)距海距離(v2)及距主要道路距離(v3)和距次要道路距離(v4);(3)以站點(diǎn)位置為中心設(shè)置500m、1km、3km、5km、7km、10km緩沖圓,計(jì)算個(gè)緩沖圓內(nèi)道路總長(zhǎng)度(v5～v10)及人口密度(v11～v16);(4)由于最終圖像分辨率為8km×8km,對(duì)于歸一化植被指數(shù)NDVI取監(jiān)測(cè)站點(diǎn)4km緩沖區(qū)內(nèi)NDVI平均值.模型構(gòu)建利用數(shù)據(jù)分析軟件SPSS22.0實(shí)現(xiàn),采用逐步多元回歸,最終剩余4類因子,多元回歸模型如下:

式中:x1代表監(jiān)測(cè)站點(diǎn)所在高程值,km;x2表示監(jiān)測(cè)站點(diǎn)離海岸線距離,km;x3表示5km緩沖區(qū)內(nèi)道路總長(zhǎng)度,km;x4表示4km緩沖半徑內(nèi)歸一化植被指數(shù)NDVI平均值.

1.3 BME時(shí)空分析

貝葉斯最大熵(BME)是由George Christakos在1990年提出的一種基于認(rèn)知學(xué)的時(shí)空數(shù)據(jù)分析方法[21-22],其最大的優(yōu)勢(shì)在于能夠綜合利用各種來(lái)源不同的不確定性數(shù)據(jù),顯著地提高未采樣點(diǎn)的估計(jì)精度[23].

BME方法的實(shí)施基于Matlab的Seks-gui時(shí)空分析軟件包(http://140.112.63.249/SEKSGUI/SEKSHome.html).用一個(gè)時(shí)空隨機(jī)場(chǎng)(space/time randomfield, S/TRF)來(lái)描述PM2.5的時(shí)空變異,用Z(p)=Z(s,t)表示時(shí)空隨機(jī)場(chǎng)中PM2.5的月均濃度,s表示觀測(cè)點(diǎn)所在經(jīng)緯度,t表示月份.那么最終PM2.5質(zhì)量濃度的估計(jì)值可由式(2～3)計(jì)算得出.

式中:Rs,t表示實(shí)際觀測(cè)值與LUR模型預(yù)測(cè)值的差值,即殘差;Z(s,t)表示最終的預(yù)測(cè)值, ZLUR(pre)表示LUR模型預(yù)測(cè)值,Os,t表示實(shí)際觀測(cè)值.

假設(shè)LUR模型的殘差Rs,t同樣能用S/TRF來(lái)描述其時(shí)空變異,利用BME對(duì)其進(jìn)行時(shí)空分析.本研究中,一般知識(shí)由時(shí)空平均趨勢(shì)函數(shù)mr(p)=E[R(p)]和協(xié)方差函數(shù)CR(p,p')=E{[R(p)-mr(p)][R(p')-mr(p')]}組成,特定知識(shí)由計(jì)算各個(gè)站點(diǎn)的殘差Rs,t組成.BME估計(jì)可由式(4)獲得:

其中:fG( r)是由一般知識(shí)獲得的關(guān)于r的先驗(yàn)概率密度函數(shù),fk( rk)是r在估計(jì)點(diǎn)pk的后驗(yàn)概率密度函數(shù),A是歸一化常數(shù).根據(jù)后驗(yàn)概率密度函數(shù),可以計(jì)算得出r在估計(jì)點(diǎn)pk處的估計(jì)均值:

式中:rk,mean是r在pk估計(jì)點(diǎn)處的估計(jì)均值,帶回式(3)即可計(jì)算得出PM2.5質(zhì)量濃度的估計(jì)值.

2 結(jié)果及討論

2.1 LUR/BME模型

本文LUR模型最終剩余的顯著性因子包括高程,距海距離,道路長(zhǎng)度,NDVI(表1).高程體現(xiàn)了由于研究區(qū)域廣闊,地形變換復(fù)雜所導(dǎo)致的PM2.5空間分布上的差異性,許多研究表明,一般來(lái)說(shuō)大部分地區(qū)高程與PM2.5質(zhì)量濃度呈負(fù)相關(guān)[6,24-25],在本文研究區(qū)域最終模型也是如此.同時(shí)本文研究區(qū)域主要是沿海省市,地處海陸交互作用地帶,相比于其它內(nèi)陸地區(qū)LUR模型的研究,距海距離在本文最終模型中貢獻(xiàn)很大,反映了距海越遠(yuǎn),PM2.5質(zhì)量濃度越高的宏觀趨勢(shì).而城市微觀方面,汽車尾氣及工業(yè)廢氣的排放是PM2.5主要來(lái)源,由于缺乏具體的排放清單及交通數(shù)據(jù),本文以道路長(zhǎng)度體現(xiàn)污染源的排放情況,由最終模型參數(shù)可知,道路長(zhǎng)度與PM2.5質(zhì)量濃度呈正相關(guān),符合科學(xué)認(rèn)知.但本文LUR模型最終擬合度R2(0.36)小于Reyes等[20]對(duì)于全美國(guó)PM2.5研究中的LUR模型擬合度R2(0.53),原因之一可能是由于Reyes構(gòu)建模型中引入了詳細(xì)的固定點(diǎn)污染源(NOx,SO2,NH3)排放數(shù)據(jù)及機(jī)動(dòng)車排放數(shù)據(jù),而Beckerman等[22]對(duì)于全美國(guó)PM2.5研究中的LUR模型構(gòu)建中引入遙感數(shù)據(jù),最終擬合度R2(0.21),本文LUR模型R2高于Beckerman的研究,雖然主要原因得益于本文相比于Beckerman的研究區(qū)域面積來(lái)說(shuō)較小,但說(shuō)明本文的模型一定程度上反應(yīng)了PM2.5的空間分布特征.另一方面許多研究表明植被覆蓋即地區(qū)的綠化程度對(duì)于空氣污染具有緩解作用[6,11],最終模型中NDVI與PM2.5質(zhì)量濃度呈顯著負(fù)相關(guān).

由于第1階段的LUR模型不足以充分解釋大范圍區(qū)域PM2.5的時(shí)空分布,本文在第2階段引入BME模型,BME作為一種有效的時(shí)空插值方法,能充分考慮空間及時(shí)間上的變異(圖3),對(duì)LUR模型得到的殘差進(jìn)行時(shí)空插值,得到全領(lǐng)域殘差的時(shí)空分布,最后將LUR估測(cè)值與BME對(duì)殘差的估測(cè)值相加,得到全領(lǐng)域的PM2.5時(shí)空分布.

表1 回歸分析參數(shù)Table 1 Parameters of correlation analysis

圖3 BME模型協(xié)方差擬合Fig.3 Covariance model fitting of BME

2.2 海岸帶PM2.5時(shí)空制圖精度及討論

本文研究區(qū)域較廣,僅以LUR模型并不能較好的展現(xiàn)全領(lǐng)域污染物分布的特性,且LUR模型在表現(xiàn)時(shí)間變異上也較為欠缺,而BME作為科學(xué)性的時(shí)空數(shù)據(jù)分析方法則能較好的與LUR相結(jié)合.本研究從5376個(gè)月平均觀測(cè)值中抽取1000個(gè)真實(shí)觀測(cè)值作為驗(yàn)證,用其他觀測(cè)值進(jìn)行模型預(yù)測(cè),結(jié)果也表明在引入BME后,預(yù)測(cè)精度得到了顯著提高,均方根誤差RMSE由23.53μg/ m3降低到11.08μg/m3.Akita等[18]對(duì)于西班牙加泰羅尼亞NO2分布估測(cè)的研究中也有結(jié)論表明LUR相比于傳統(tǒng)地統(tǒng)計(jì)學(xué)(克里金)更加精確.但Akita的研究也表明LUR與其它方法如大氣擴(kuò)散模型等的結(jié)合更能全面反映污染物分布的空間變異,而這一點(diǎn)也得到了許多研究的證實(shí)[17,20,26-27].同時(shí)相較于傳統(tǒng)地統(tǒng)計(jì)學(xué),該類方法的使用也避免了大范圍平滑效應(yīng)(圖3).鄒濱等[2]以美國(guó)休斯頓為研究區(qū)域比較了LUR與普通克里格對(duì)PM2.5分布作圖的精度,證實(shí)LUR模型在預(yù)測(cè)精度與圖像信息熵方面更優(yōu)于普通克里格.

表2 模型比較Table 2 Comparison between LUR and LUR-BME

圖4 2015年沿海部分省份PM2.5質(zhì)量濃度時(shí)空分布Fig.4 Spatiotemporal distribution of PM2.5concentration in Coastal Areas of China, 2015

整體上看(圖4),空間尺度上,海岸帶PM2.5質(zhì)量濃度由南至北,以長(zhǎng)江三角洲為分界線,宏觀上呈現(xiàn)南低北高,京津冀地區(qū)作為重點(diǎn)污染地區(qū)空氣狀況仍然形勢(shì)嚴(yán)峻,河北西南部,山東西部及北京城區(qū)屬嚴(yán)重污染區(qū)域,一方面與該部分地區(qū)經(jīng)濟(jì)工業(yè)發(fā)展及氣候有關(guān)[12],另一方面薛文博等[28]利用顆粒物來(lái)源追蹤技術(shù)PSAT定量模擬我國(guó)PM2.5城市來(lái)源研究表明,北京,天津及石家莊PM2.5受外省影響較大,分別達(dá)37%,42%,33%.而長(zhǎng)江三角洲地區(qū)2015年P(guān)M2.5質(zhì)量濃度也普遍較高,福建,廣東等地空氣質(zhì)量則相對(duì)較好.時(shí)間尺度上,整個(gè)研究區(qū)域秋冬2季PM2.5質(zhì)量濃度普遍偏高,夏季PM2.5質(zhì)量濃度較低.

本文最終得到的PM2.5質(zhì)量濃度分布分辨率為8km×8km均勻網(wǎng)格,為比較各市PM2.5污染狀況,我們計(jì)算各省市行政區(qū)劃內(nèi)所有網(wǎng)格平均值代表其四季PM2.5平均質(zhì)量濃度(圖5).各省市四季PM2.5質(zhì)量濃度差異較大,其中冬季空氣污染最為嚴(yán)重,PM2.5質(zhì)量濃度普遍較高,最高天津市達(dá)94μg/m3,最低福建省23.4μg/m3;而各省市夏季PM2.5質(zhì)量濃度普遍較低,最低福建省7.1μg/m3.

圖5 沿海部分省市四季PM2.5平均濃度Fig.5 Seasonal mean PM2.5concentration in Coastal A reas of China, 2015

2.3 山東省空氣污染人體暴露研究

山東作為我國(guó)北方人口基數(shù)較大的省份,其空氣污染問(wèn)題非常嚴(yán)峻[29],為進(jìn)一步驗(yàn)證本方法的應(yīng)用價(jià)值,本文以山東省為典例將LUR/ BME的空間數(shù)據(jù)分析方法應(yīng)用于室外PM2.5人體暴露分析,由于人口普查數(shù)據(jù)是以行政區(qū)劃來(lái)統(tǒng)計(jì),例如高山,水域等地人口較少,如此不利于計(jì)算操作且結(jié)果不精確,本文采用第4代全球人口密度1km分辨率柵格數(shù)據(jù),通過(guò)以下公式計(jì)算[6]:

式中:Cpe表示室外平均人口暴露PM2.5質(zhì)量濃度, μg/m3;pop(i,j)表示第i行第j列柵格代表的人口密度,人/km2;c(i,j)表示第i行第j列柵格代表的PM2.5質(zhì)量濃度,μg/m3.

山東省年均PM2.5分布呈現(xiàn)明顯的階梯型分布(圖6),內(nèi)陸城市污染嚴(yán)重,而沿海城市則較為改善,以人口密度柵格數(shù)據(jù)計(jì)算人均暴露濃度,減少了PM2.5分布不均而城市地區(qū)人口密度較高所產(chǎn)生的誤差.

圖6 2015年山東省PM2.5年平均濃度空間分布Fig.6 Annual mean PM2.5concentration distribution in Shandong Province, 2015

圖7 山東省2015年P(guān)M2.5人均暴露濃度Fig.7 Per capita outdoor population exposure to PM2.5in Shandong Province, 2015

由內(nèi)陸至近海,山東省PM2.5人均暴露質(zhì)量濃度逐步遞減(圖7),濟(jì)南市作為山東省會(huì)城市,工商業(yè)最為發(fā)達(dá)且人口集中,最高達(dá)85.5μg/m3,而近海城市日照,青島,煙臺(tái),威海則相對(duì)較低,但大部分人口較為集中的地區(qū)仍高于我國(guó)現(xiàn)行的PM2.5質(zhì)量濃度年平均限值35μg/m3.一方面可能與山東內(nèi)陸及周邊城市工廠排放量較高有關(guān),另一方面關(guān)于近海氣候因素對(duì)空氣污染的影仍需進(jìn)一步研究.

3 結(jié)論

3.1 本文通過(guò)站點(diǎn)實(shí)測(cè)數(shù)據(jù)及與空氣污染相關(guān)的地理因子,采用LUR與BME相結(jié)合的空間數(shù)據(jù)分析方法,得到了海岸帶部分地區(qū)高精度的PM2.5質(zhì)量濃度格點(diǎn)分布. 結(jié)果顯示各省市四季PM2.5質(zhì)量濃度差異較大,其中冬季空氣污染普遍較為嚴(yán)重,PM2.5質(zhì)量濃度最高為天津市達(dá)94μg/m3,最低福建省23.4μg/m3;而各省市夏季PM2.5質(zhì)量濃度普遍較低,最低福建省7.1μg/m3.

3.2 將該方法應(yīng)用到山東省PM2.5室外人體暴露分析的結(jié)果表明,由內(nèi)陸至近海人均PM2.5暴露質(zhì)量濃度逐步遞減,最高濟(jì)南市85.5μg/m3,煙臺(tái)、威海等沿海城市相對(duì)較低.

3.3 由于PM2.5在時(shí)間及空間上受許多因素影響,如氣象因素,污染源等,本文并未納入到研究中,導(dǎo)致結(jié)果可能受某一因子影響較大,局部出現(xiàn)較大誤差,且本文研究是基于月平均PM2.5質(zhì)量濃度,在后續(xù)的研究中可考慮提高時(shí)間分辨率[30],或結(jié)合遙感數(shù)據(jù)得到大范圍PM2.5高分辨率時(shí)空分布.

[1] Yang Y, Christakos G, Huang W, et al. Uncertainty assessment of PM2.5contamination mapping using spatiotemporal sequential indicator simulations and multi-temporal monitoring data [J]. Science Reports, 2016,6:24-35.

[2] Zou B, Luo Y Q, Wan N, et al. Performance comparison of LUR and OK in PM2.5concentration mapping: a multidimensional perspective [J]. Science Reports, 2015,5:86-98.

[3] Lv B, Hu Y, Chang H H, et al. Improving the Accuracy of Daily PM2.5Distributions Derived fromthe Fusion of Ground-Level Measurements with Aerosol Optical Depth Observations, a Case Study in North China [J]. Environmental Science & Technology, 2016,50(9):4752-4759.

[4] Zhang W, Capps S L, Hu Y, et al. Development of the high-order decoupled direct method in three dimensions for particulate matter: enabling advanced sensitivity analysis in air quality models [J]. Geoscientific Model Development Discussions, 2011,4(4):355-368.

[5] 焦利民,許 剛,趙素麗,等.基于LUR的武漢市PM2.5濃度空間分布模擬 [J]. 武漢大學(xué)學(xué)報(bào)(信息科學(xué)版), 2015,40(8):1088-1094.

[6] Wu J S, Li J C, Peng J, et al. Applying land use regression model to estimate spatial variation of PM2.5in Beijing, China [J]. Environmental Science Pollution Research, 2015,22(9):7045-7061.

[7] Briggs D J, Collins S, Elliott P, et al. Mapping urban air pollution using GIS: a regression-based approach [J]. International Journal of Geographical Information Science, 1997,11(7):699-718.

[8] Henderson S B, Beckerman B, Jerret M, et al. Application of land use regression to estimate long-termconcentrations of traffic-related nitrogen oxides and fine particulate matter [J]. Environmental Science & Technology, 2007,41(7):2422-2428.

[9] Miyawaki J, Kamei S, Sakayama K, et al. Modeling the Intraurban Variability of Ambient Traffic Pollution in Toronto, Canada [J]. Journal of Toxicology & Environmental Health Part A, 2007,70(3/4):200-212.

[10] Beckerman B S, Jerrett M, Martin R V, et al. Application of the deletion/substitution/addition algorithmto selecting land use regression models for interpolating air pollution measurements in California [J]. Atmospheric Environment, 2013,77(68):172-177.

[11] Hassan A, Seyed Mahmood T S, Henderson S B, et al. Land use regression models to estimate the annual and seasonal spatial variability of sulfur dioxide and particulate matter in Tehran, Iran. [J]. Science of the Total Environment, 2014,489(488-489C): 343-353.

[12] 陳 莉,白志鵬,蘇 笛,等.利用LUR模型模擬天津市大氣污染物濃度的空間分布 [J]. 中國(guó)環(huán)境科學(xué), 2009,29(7):685-691.

[13] Li C, Du S Y, Bai Z P, et al. Application of land use regression for estimating concentrations of major outdoor air pollutants in Jinan, China [J]. Journal of Zhejiang University-SCIENCE A, 2010, 11(11):857-867.

[14] Liu W, Li X, Chen Z, et al. Land use regression models coupled with meteorology to model spatial and temporal variability of NO2, and PM10, in Changsha, China [J]. A tmospheric Environment, 2015,116:272-280.

[15] Liu C, Henderson B H, Wang D, et al. A land use regression application into assessing spatial variation of intra-urban fine particulate matter PM2.5and nitrogen dioxide NO2concentrations in City of Shanghai, China. [J]. Science of the Total Environment, 2016,565:607-615.

[16] Akita Y, Baldasano J M, Beelen R, et al. Large scale air pollution estimation method combining land use regression and chemical transport modeling in a geostatistical framework. [J]. Environmental Science & Technology, 2014,48(8):4452-4459.

[17] Yu H L, Wang C H, Liu MC, et al. Estimation of Fine Particulate Matter in Taipei Using Landuse Regression and Bayesian MaximumEntropy Methods [J]. International Journal ofEnvironmental Research & Public Health, 2011,8(6):2153-2169.

[18] Reyes J, Serre ML. An LUR/BME Framework to Estimate PM2.5Explainedby on Road Mobile and StationarySources [J]. Environmental Science & Technology, 2014,48(3):1736-1744.

[19] Matheron G. The Theory of Regionalized Variables and Its Application [J]. 1971,5:218.

[20] Beckerman B S, Jerrett M, Serre M, et al. A Hybrid Approach to Estimating National Scale Spatiotemporal Variability of PM2.5in the Contiguous United States [J]. Environmental Science & Technology, 2013,47(13):7233-7241.

[21] Christakos G. A Bayesian maximum-entropy viewto the spatial estimation problem[J]. Mathematical Geosciences, 1990,22(7): 763-777.

[22] Christakos G, Serre ML. BME analysis of spatiotemporal particulate matter distributions in North Carolina [J]. Atmospheric Environment, 2000,34(20):3393-3406.

[23] 楊 勇,張若兮.貝葉斯最大熵地統(tǒng)計(jì)方法研究與應(yīng)用進(jìn)展 [J].土壤, 2014(3):402-406.

[24] Henderson S B, Beckerman B, Jerrett M, Brauer M. Application of Land Use Regression to Estimate Long-TermConcentrations of Traffic-Related Nitrogen Oxides and Fine Particulate Matter [J]. Environmental Science & Technology, 2007,41(7):2422-2428.

[25] Hart J E, Yanosky J D, Puett R C, et al. Spatial Modeling of PM10and NO2in the Continental United States, 1985～2000 [J]. Environmental Health Perspectives, 2009,117(11):1690-1696.

[26] Wang M, Sampson P D, Hu J, et al. Combining Land-Use Regression and Chemical Transport Modeling in a Spatiotemporal Geostatistical Model for Ozone and PM2.5. [J]. Environmental Science & Technology, 2016,50(10):5111-5118.

[27] Bechle MJ, Millet D B, Marshall J D. National Spatiotemporal Exposure Surface for NO2: Monthly Scaling of a SatelliteDerived Land-Use Regression, 2000～2010. [J]. Environmental Science & Technology, 2015,49(20):12297-12305.

[28] 薛文博,付 飛, 王金南,等.中國(guó)PM2.5跨區(qū)域傳輸特征數(shù)值模擬研究 [J]. 中國(guó)環(huán)境科學(xué), 2014,34(6):1361-1368.

[29] Yang Y, Christakos G. Spatiotemporal Characterization of Ambient PM2.5Concentrations in Shandong Province (China). [J]. Environmental Science & Technology, 2015,49(22):13431-13438.

[30] 郭 宇,鄒 濱,鄭 忠,等.高時(shí)空分辨率PM2.5濃度土地利用回歸模擬與制圖 [J]. 遙感信息, 2015(5):94-101.

Spatiotemporal analysis of PM2.5in large coastal domains by combining Land Use Regression and Bayesian MaximumEntropy.

JIANG Qu-tu1, HE Jun-yu1, WANG Zhan-shan2, YE Guan-qiong1*, CHEN Qian3, XIAO Lu1
(1.Institute of Island & Coastal Ecosystems, Zhejiang University, Zhoushan 316021, China；2.Beijing Municipal Environmental Monitoring Center, Beijing 100048, China；3.School of Geographic and Environment, Jiangxi Normal University, Nanchang 330022, China). China Environmental Science, 2017,37(2)：424～431

By combining Land Use Regression (LUR) and Bayesian MaximumEntropy (BME), this study constructed a LUR model based on the parameters of elevation, distance to sea, length of roads and Normalized Difference Vegetation Index (NDVI) to generate a global map of PM2.5distribution in a large costal area in 2015, china. The Bayesian MaximumEntropy was further introduced in the interpolation of LUR space-time residuals. Because of the introduction of BME, the cross-validation results showed that the R2increased from0.36 to 0.85, and the root-mean-square error (RMSE) decreased from23.53μg/m3to 11.08μg/m3. The average concentration of PM2.5in the northern coastal areas was higher than that of the southern areas, and the highest concentration of PM2.5appeared in the inland area of Beijing, Tianjin, Hebei and Shandong provinces during winter times. The annual spatial distribution of PM2.5was further integrated with population density in Shandong province for risk exposure analysis. The outcome showed that the outdoor population exposure of PM2.5decreased frominland to sea, and the highest Per capita outdoor exposure value occurred in the central city, Jinan (85.5μg/m3), while the lowest value occurred in coastal areas of Yantai and Weihai.

land use regression；bayesian maximumentropy；PM2.5；human exposure

X51

A

1000-6923(2017)02-0424-08

江曲圖(1992-),男,湖北咸寧人,浙江大學(xué)海洋學(xué)院碩士研究生,主要從事海岸帶環(huán)境健康研究.

2016-06-18

科技部國(guó)際合作項(xiàng)目(2015DFA01410)

* 責(zé)任作者, 講師, gqy@zju.edu.cn