江蘇省級(jí)區(qū)域空氣質(zhì)量數(shù)值預(yù)報(bào)模式效果評(píng)估

朱莉莉，晏平仲，王自發(fā)，李 杰，張祥志，湯莉莉，李健軍，劉 冰

1．中國(guó)科學(xué)院大氣物理研究所，大氣邊界層物理和大氣化學(xué)國(guó)家重點(diǎn)實(shí)驗(yàn)室，北京 100029

2．中國(guó)科學(xué)院大學(xué)，北京 100049

3．江蘇省環(huán)境監(jiān)測(cè)中心，江蘇 南京 210036

4．中國(guó)環(huán)境監(jiān)測(cè)總站，國(guó)家環(huán)境保護(hù)環(huán)境監(jiān)測(cè)質(zhì)量控制重點(diǎn)實(shí)驗(yàn)室，北京 100012

隨著中國(guó)城市化進(jìn)程的加快，區(qū)域灰霾污染現(xiàn)象頻繁發(fā)生，引起社會(huì)的廣泛關(guān)注。在巨大的能源消耗和污染物排放背景下，相當(dāng)長(zhǎng)的一段時(shí)間內(nèi)，中國(guó)大部分地區(qū)的空氣質(zhì)量仍將處于重污染頻發(fā)期，區(qū)域大氣污染聯(lián)防聯(lián)控成為一項(xiàng)長(zhǎng)期而艱巨的任務(wù)，需要建立完善的空氣質(zhì)量預(yù)報(bào)預(yù)警系統(tǒng)，及時(shí)有效地支持重污染天氣應(yīng)急減排，緩解空氣污染對(duì)公眾健康、交通出行等多方面的影響。

2010—2012年，江蘇省 PM2.5年均值均超過(guò)國(guó)家二級(jí)標(biāo)準(zhǔn)［1］50%以上，大氣污染問(wèn)題已非常嚴(yán)峻。2013年，國(guó)務(wù)院印發(fā)了《大氣污染防治行動(dòng)計(jì)劃》，明確要求“建立監(jiān)測(cè)預(yù)警應(yīng)急體系，妥善應(yīng)對(duì)重污染天氣”［2］。2014年夏季，南京舉辦第二屆夏季青年奧林匹克運(yùn)動(dòng)會(huì)，對(duì)江蘇省空氣質(zhì)量保障工作提出了更高的要求。

空氣質(zhì)量數(shù)值預(yù)報(bào)以大氣動(dòng)力學(xué)理論為基礎(chǔ)，在給定的氣象場(chǎng)、源排放以及初始和邊界條件下，通過(guò)一套復(fù)雜的偏微分方程組描述大氣污染物在空氣中的各種物理化學(xué)過(guò)程(輸送、擴(kuò)散、轉(zhuǎn)化、清除等)，并利用高性能計(jì)算集群進(jìn)行數(shù)值求解計(jì)算，預(yù)報(bào)污染物濃度動(dòng)態(tài)分布和變化趨勢(shì)。研究采用中國(guó)科學(xué)院大氣物理研究所嵌套網(wǎng)格空氣質(zhì)量模式系統(tǒng)(NAQPMS)空氣質(zhì)量模式，建立江蘇省級(jí)區(qū)域空氣質(zhì)量數(shù)值預(yù)報(bào)系統(tǒng)，實(shí)現(xiàn)對(duì)PM2.5質(zhì)量濃度的未來(lái)24 h預(yù)報(bào)和未來(lái)7 d趨勢(shì)預(yù)測(cè)。并利用2013年6—8月PM2.5質(zhì)量濃度觀測(cè)資料評(píng)估了該系統(tǒng)對(duì)江蘇省PM2.5質(zhì)量濃度未來(lái)24 h預(yù)報(bào)和7 d潛勢(shì)預(yù)測(cè)效果。

1 模式的業(yè)務(wù)化及評(píng)估指標(biāo)的確立

1.1數(shù)值模式及設(shè)置

NAQPMS［3］是中國(guó)科學(xué)院大氣物理研究所基于“一個(gè)大氣”的理念自主開發(fā)的三維歐拉化學(xué)傳輸模式，廣泛應(yīng)用于區(qū)域大氣污染數(shù)值模擬研究。NAQPMS模式考慮氣溶膠及其前體物排放后在大氣中的平流、擴(kuò)散、干沉降和濕沉降等物理過(guò)程以及氣相化學(xué)、液相化學(xué)和氣溶膠非均相化學(xué)等化學(xué)過(guò)程，模式方案詳細(xì)介紹見文獻(xiàn)［4］。目前，北京、上海、廣州等多個(gè)城市已采用NAQPMS模式建成數(shù)值預(yù)報(bào)業(yè)務(wù)系統(tǒng)，并在北京奧運(yùn)會(huì)［5-6］、上海世博會(huì)［7］和廣州亞運(yùn)會(huì)［8］成功應(yīng)用。

模式區(qū)域覆蓋中國(guó)中東部地區(qū)，中心經(jīng)緯度為115°E、25°N，水平網(wǎng)格格距為 45 km，網(wǎng)格數(shù)為120(東西向)×162(南北向)。采用Sigma-Z地形追隨坐標(biāo)系，垂直方向不均勻地劃分為20層，其中2 km高度以下劃分8層，近地層中心高度約47 m，模式層頂海拔高度為20 km。

采用中尺度氣象模式WRFv 3.5預(yù)報(bào)未來(lái)天氣過(guò)程，以驅(qū)動(dòng)空氣質(zhì)量模式預(yù)測(cè)未來(lái)區(qū)域污染形勢(shì)和城市空氣質(zhì)量。氣象預(yù)報(bào)的初始、邊界條件取自美國(guó)國(guó)家環(huán)境預(yù)報(bào)中心(NCEP)的GFS全球氣象預(yù)報(bào)資料，時(shí)間間隔為24 h，空間分辨率為0.5°×0.5°。以前日北京時(shí)20:00為模式起始時(shí)刻，每日計(jì)算未來(lái)192 h。WRF模式輸出逐時(shí)氣象要素，包括氣壓場(chǎng)、風(fēng)場(chǎng)、溫度場(chǎng)、濕度場(chǎng)、云量和降水等。

人為排放清單采用日本國(guó)立環(huán)境研究所編制的亞洲區(qū)域大氣污染物排放清單REAS 2.0［9］，空間分辨率為0.25°×0.25°，基準(zhǔn)年為2007年，排放物種類包括二氧化硫(SO2)、氮氧化物(NOx)、非甲烷揮發(fā)性有機(jī)物(NMVOC)、黑碳(BC)、一氧化碳(CO)和一次排放的細(xì)顆粒物(PMF)等。碳?xì)浠衔锏淖匀辉磁欧湃∽匀蛏l(fā)物質(zhì)名錄行動(dòng)［10］。圖1給出了中國(guó)中東部地區(qū)SO2排放速率的水平分布。

圖1 江蘇省SO2排放速率與空氣質(zhì)量監(jiān)測(cè)站點(diǎn)分布圖(圖中圓點(diǎn)表示監(jiān)測(cè)站點(diǎn)位置)

由圖1可以看出，華北平原和長(zhǎng)江中下游平原為SO2排放速率的高值區(qū)。江蘇省北部與華北平原高排放區(qū)相連，南部處于長(zhǎng)三角排放高值區(qū)。江蘇南部SO2排放速率可達(dá)1 μg/(m2·s);江蘇北部的徐州、宿遷等部分地區(qū)SO2排放速率也在0.4 μg/(m2·s)以上。江蘇相鄰省份中，南部的浙江省和西部的安徽省SO2排放速率相對(duì)較低，北部的山東省SO2排放速率較高，國(guó)內(nèi)外廣泛應(yīng)用的其他幾個(gè)排放源清單均顯示類似的分布特征［11］。

選取2013年6—8月江蘇省13個(gè)城市未來(lái)7 d PM2.5的預(yù)報(bào)預(yù)測(cè)結(jié)果進(jìn)行效果評(píng)估，每日預(yù)報(bào)以前日北京時(shí)20:00為起始時(shí)刻，預(yù)測(cè)未來(lái)192 h的污染物濃度的時(shí)空分布，未來(lái)7 d預(yù)報(bào)預(yù)測(cè)值分別用 F1、F2、F3、F4、F5、F6 和 F7 表示。預(yù)報(bào)系統(tǒng)從5月1日開始模式初始化(spin-up)，每日預(yù)報(bào)初始條件取自前日模式結(jié)果，區(qū)域邊界條件設(shè)置為模式缺省值。模式的PM2.5組分包含硫酸鹽(SO42－)、硝酸鹽(NO3－)、銨鹽(NH4+)、BC、有機(jī)物(OM)、PMF、沙塵和海鹽等。

1.2 觀測(cè)資料

PM2.5觀測(cè)資料來(lái)源于全國(guó)城市空氣質(zhì)量實(shí)時(shí)發(fā)布平臺(tái)［12］的江蘇省72個(gè)國(guó)控站小時(shí)濃度數(shù)據(jù)。觀測(cè)站點(diǎn)位置信息見圖1?？紤]到模式水平網(wǎng)格格距為45 km，多個(gè)觀測(cè)站點(diǎn)可能位于同一模式網(wǎng)格，采用網(wǎng)格內(nèi)觀測(cè)數(shù)據(jù)的站點(diǎn)平均值與網(wǎng)格預(yù)報(bào)值進(jìn)行驗(yàn)證。參考《環(huán)境空氣質(zhì)量標(biāo)準(zhǔn)》(GB 3095—2012)［1］，采用每日 00:00 開始的24 h濃度平均值進(jìn)行驗(yàn)證。在觀測(cè)資料預(yù)處理時(shí)，首先選取位于相同模式網(wǎng)格的所有站點(diǎn)的每小時(shí)濃度的平均值作為該網(wǎng)格的小時(shí)濃度，然后計(jì)算該網(wǎng)格每日24 h濃度平均值。每日參與計(jì)算的小時(shí)濃度數(shù)據(jù)少于20 h時(shí)，當(dāng)日濃度平均值記為缺測(cè)。

1.3 評(píng)估指標(biāo)的確定

在氣象、環(huán)境的模式研究領(lǐng)域通用的統(tǒng)計(jì)指標(biāo)中，選取相關(guān)系數(shù)(r)、均方根誤差(RMSE)、平均偏差(MB)、平均誤差(ME)、標(biāo)準(zhǔn)化分?jǐn)?shù)偏差(MFB)、標(biāo)準(zhǔn)化分?jǐn)?shù)誤差(MFE)、標(biāo)準(zhǔn)化平均偏差(NMB)和標(biāo)準(zhǔn)化平均誤差(NME)共8項(xiàng)指標(biāo)對(duì)PM2.5數(shù)值預(yù)報(bào)效果進(jìn)行評(píng)估，具體指標(biāo)定義和計(jì)算公式參見文獻(xiàn)［13］。其中，r反映預(yù)報(bào)值與觀測(cè)值隨時(shí)間變化趨勢(shì)的相似程度，是一個(gè)無(wú)量綱數(shù);r＞0表示預(yù)報(bào)值與觀測(cè)值正相關(guān)，r＜0表示預(yù)報(bào)值與觀測(cè)值負(fù)相關(guān)，通過(guò)統(tǒng)計(jì)顯著性檢驗(yàn)判定預(yù)報(bào)值和觀測(cè)值的相關(guān)性。此外，RMSE、MB和ME 3項(xiàng)參數(shù)的單位是μg/m3，使用時(shí)需要結(jié)合觀測(cè)濃度大小進(jìn)行評(píng)估，這3項(xiàng)指標(biāo)只適用于觀測(cè)濃度值接近的情況。而NMB和NME 2項(xiàng)指標(biāo)克服了上述3項(xiàng)指標(biāo)在觀測(cè)濃度差異較大時(shí)難以直接比較的缺點(diǎn)。MFB和MFE數(shù)值的變化范圍分別為(－100%，∞)和(0，∞)。但是NMB在預(yù)報(bào)值低估與高估時(shí)的數(shù)值大小并不對(duì)稱(如預(yù)報(bào)值僅為觀測(cè)值的二分之一時(shí)NMB等于－50%，但預(yù)報(bào)值為觀測(cè)值的2倍時(shí)NMB等于200%)。NME也存在類似的問(wèn)題。

Boylan等［13］認(rèn)為 MFB和 MFE指標(biāo)更適用于評(píng)估顆粒物的模式預(yù)報(bào)效果。這2項(xiàng)參數(shù)克服了NMB和MFE在預(yù)報(bào)值低估與高估時(shí)的數(shù)值大小不對(duì)稱的問(wèn)題，并且考慮到觀測(cè)值的不確定性，采用觀測(cè)值和預(yù)報(bào)值的平均值作為效果評(píng)估的參考值。同時(shí)，給出了模式預(yù)報(bào)的“合理”范圍為－60%≤MFB≤60%、MFE≤75%，模式預(yù)報(bào)效果達(dá)到“理想”水平的范圍為 －30%≤MFB≤30%、MFE≤50%。

2 預(yù)報(bào)預(yù)測(cè)效果評(píng)估

2.1 未來(lái)24 h濃度預(yù)報(bào)效果評(píng)估

圖2(a)展示了中國(guó)中東部2013年夏季預(yù)報(bào)和觀測(cè)的PM2.5質(zhì)量濃度日均值的空間分布，為清晰展示研究區(qū)域濃度分布，圖2(b)將江蘇省區(qū)域放大，其中預(yù)報(bào)值為未來(lái)24 h(即模式預(yù)報(bào)的第5～28 h)的預(yù)報(bào)結(jié)果。

圖2 2013年夏季PM2.5質(zhì)量濃度預(yù)報(bào)及觀測(cè)結(jié)果空間分布

圖2(a)顯示，華北平原出現(xiàn)高污染區(qū)，PM2.5質(zhì)量濃度超過(guò) 75 μg/m3，達(dá)到輕度污染［1］，江蘇北部位于該污染區(qū)南部;江蘇南部位于長(zhǎng)三角次高污染區(qū)，PM2.5質(zhì)量濃度超過(guò) 60 μg/m3，超過(guò)WHO 日平均指導(dǎo)值(25 μg/m3)［14］。

圖2(b)顯示，江蘇省內(nèi)，PM2.5質(zhì)量濃度預(yù)報(bào)結(jié)果呈現(xiàn)南北高、中部低的分布特征，與排放源的分布相對(duì)應(yīng)。江蘇省多數(shù)城市的PM2.5質(zhì)量濃度預(yù)報(bào)值與觀測(cè)值較為吻合，南部部分城市預(yù)報(bào)結(jié)果有所高估。

2013年夏季，南京市PM2.5質(zhì)量濃度平均值為 45 μg/m3(標(biāo)準(zhǔn)差 18 μg/m3)，定義 PM2.5質(zhì)量濃度超過(guò)國(guó)家二級(jí)標(biāo)準(zhǔn)為超標(biāo)，則超標(biāo)率為7%。選取無(wú)錫市和徐州市分別作為江蘇省南部和北部的代表城市進(jìn)行預(yù)報(bào)效果評(píng)估。無(wú)錫市PM2.5質(zhì)量濃度平均值為 43 μg/m3(標(biāo)準(zhǔn)差 22 μg/m3)，超標(biāo)率為8%;徐州市PM2.5質(zhì)量濃度平均值為52 μg/m3(標(biāo)準(zhǔn)差 46 μg/m3)，超標(biāo)率為 17%。南京與無(wú)錫的情況較為接近，徐州市的PM2.5質(zhì)量濃度和污染超標(biāo)率較高，并且徐州市的樣本離散度也較高，濃度值波動(dòng)更為顯著。

圖3為2013年夏季南京、無(wú)錫和徐州PM2.5質(zhì)量濃度預(yù)報(bào)值與觀測(cè)值時(shí)間序列與散點(diǎn)圖，其中預(yù)報(bào)值為未來(lái)24 h預(yù)報(bào)結(jié)果。

圖3 2013年夏季南京、無(wú)錫和徐州PM2.5質(zhì)量濃度預(yù)報(bào)結(jié)果與觀測(cè)時(shí)間序列、散點(diǎn)圖

由圖3可見，南京市PM2.5質(zhì)量濃度的預(yù)報(bào)值與觀測(cè)值的變化趨勢(shì)較為一致，r值達(dá)到0.69，預(yù)報(bào)平均值為63 μg/m3，與觀測(cè)值相比高估18 μg/m3，高估現(xiàn)象較多發(fā)生在 7、8月。無(wú)錫市PM2.5質(zhì)量濃度的預(yù)報(bào)值與觀測(cè)值的變化趨勢(shì)較為吻合，r值達(dá)到0.82，預(yù)報(bào)平均值為50 μg/m3，與觀測(cè)值相比略有偏高。徐州市7、8月的PM2.5質(zhì)量濃度預(yù)報(bào)值出現(xiàn)了明顯的高估現(xiàn)象，而6月的PM2.5質(zhì)量濃度預(yù)報(bào)值出現(xiàn)嚴(yán)重低估。

預(yù)報(bào)值處于觀測(cè)值的0.5～2倍范圍內(nèi)即為合理。散點(diǎn)圖中采用FAC表示預(yù)報(bào)值落于觀測(cè)值的0.5～2倍范圍內(nèi)的比例。南京市預(yù)報(bào)值的FAC為90%，無(wú)錫市 FAC高達(dá)97%，而徐州市FAC僅為55%。根據(jù)衛(wèi)星火點(diǎn)資料推測(cè)［15］，6月中旬江蘇省北部大部分地區(qū)出現(xiàn)大量火點(diǎn)，而預(yù)報(bào)模式的排放源尚未動(dòng)態(tài)考慮這些生物質(zhì)燃燒排放，可能是導(dǎo)致江蘇北部城市(如徐州)預(yù)報(bào)低估的重要原因。

表1是2013年夏季江蘇省13個(gè)城市PM2.5質(zhì)量濃度預(yù)報(bào)效果統(tǒng)計(jì)參數(shù)表，根據(jù)地理氣候特征將所有城市分為南部(沿江8市)和北部(蘇北5市)，其中MM和MO分別是模式模擬和觀測(cè)的PM2.5日均濃度的季度平均值。

表1 2013年夏季江蘇省各城市PM2.5質(zhì)量濃度預(yù)報(bào)統(tǒng)計(jì)參數(shù)

由表1可見，沿江8市夏季PM2.5質(zhì)量濃度預(yù)報(bào)與觀測(cè)的MB大部分都在±15 μg/m3范圍內(nèi)，ME 均在 30 μg/m3以下，RMSE 均在 35 μg/m3以下，r值大部分在0.5以上，均通過(guò)99%顯著性檢驗(yàn)。可見，對(duì)于沿江8市的預(yù)報(bào)系統(tǒng)性偏差較小，并且PM2.5質(zhì)量濃度的預(yù)測(cè)值與觀測(cè)值的變化趨勢(shì)大致吻合。蘇北5市的MB全部為負(fù)值，即對(duì)蘇北5市2013年夏季的PM2.5質(zhì)量濃度的預(yù)報(bào)整體偏低，并且大部分城市低估超過(guò)25 μg/m3;ME最大值可達(dá)到49 μg/m3;RMSE值大部分超過(guò)50 μg/m3;r值普遍低于0.35?？梢妼?duì)于蘇北5市的PM2.5濃度值和濃度變化趨勢(shì)的預(yù)報(bào)偏差較大。表1中MFB和MFE的計(jì)算結(jié)果顯示，沿江8市PM2.5預(yù)報(bào)結(jié)果均達(dá)到“理想”水平，而蘇北5市大部分城市的預(yù)報(bào)結(jié)果未達(dá)到“理想”水平，但所有城市的預(yù)報(bào)結(jié)果均在“合理”范圍內(nèi)?？傮w看來(lái)，沿江8市的預(yù)報(bào)效果普遍優(yōu)于蘇北5市。

目前，國(guó)內(nèi)關(guān)于長(zhǎng)三角地區(qū)的PM2.5質(zhì)量濃度的預(yù)報(bào)效果評(píng)估研究較少。Heather Simon等［16］對(duì)2006—2012年3月北美地區(qū)PM2.5質(zhì)量濃度的數(shù)值預(yù)報(bào)進(jìn)行了相關(guān)研究，其中r2為0～0.85，中位數(shù)為0.34左右;NMB范圍為±0.6;NME范圍為－0.35～0.75，中位數(shù)分別為0.01和0.42;MFB和MFE范圍為－0.55～0.25和0.25～0.7，中位數(shù)分別為0.15和0.5左右。江蘇省各個(gè)城市2013年夏季PM2.5質(zhì)量濃度日均值的預(yù)報(bào)結(jié)果均在上述范圍內(nèi)，其中南部城市的預(yù)報(bào)效果與北美的預(yù)報(bào)效果較為接近，北部城市的預(yù)報(bào)效果較差。

2.2 未來(lái)7 d污染趨勢(shì)預(yù)測(cè)效果評(píng)估

在環(huán)境空氣質(zhì)量數(shù)值預(yù)報(bào)中，污染源短時(shí)間并未發(fā)生顯著變化時(shí)，天氣形勢(shì)和氣象要素的預(yù)報(bào)性能直接影響著污染物濃度的預(yù)報(bào)效果。氣象模式的預(yù)報(bào)誤差會(huì)隨著時(shí)間推移而不斷增大，因此不同預(yù)報(bào)起始時(shí)間的PM2.5質(zhì)量濃度的預(yù)報(bào)效果也會(huì)存在差異。采用RMSE、NMB、MFB和MFE 4項(xiàng)指標(biāo)評(píng)價(jià)不同預(yù)報(bào)時(shí)效F1～F7對(duì)江蘇省13個(gè)城市PM2.5質(zhì)量濃度預(yù)測(cè)預(yù)報(bào)效果，如圖4所示。

由圖4可見，盡管大部分城市F1的RMSE和NMB明顯優(yōu)于 F2～F7，F(xiàn)2～F7的 MFB和 MFE參數(shù)值也都處于“合理”范圍，特別是F2～F5的MFE數(shù)值也同樣達(dá)到“理想”范圍。同時(shí)F2～F5的RMSE和NMB水平較為接近。其中，F(xiàn)2～F5的RMSE超過(guò)F1約1.5倍，F(xiàn)2～F5的NMB超過(guò)F1 50%以上，這表明未來(lái)可以構(gòu)建基于時(shí)間的集合預(yù)報(bào)產(chǎn)品。

利用F1～F7的7個(gè)預(yù)報(bào)值，可構(gòu)造出每日PM2.5質(zhì)量濃度預(yù)報(bào)值分布區(qū)間，提供更多預(yù)報(bào)信息。以無(wú)錫市(圖5)為例，當(dāng)F1～F7較接近、離散度較小時(shí)，預(yù)報(bào)值與觀測(cè)值較一致，預(yù)報(bào)效果較優(yōu);當(dāng)F1～F7差異較大、離散度較大時(shí)，預(yù)報(bào)效果有所降低。定義觀測(cè)值落于F1～F7預(yù)報(bào)值區(qū)間內(nèi)為預(yù)報(bào)命中。2013年夏季，江蘇省13地市的預(yù)報(bào)命中率的平均值為50%，其中蘇州市預(yù)報(bào)命中率最高，達(dá)到70%。無(wú)錫市預(yù)報(bào)命中率為江蘇省平均預(yù)報(bào)命中水平50%。比較F1～F7發(fā)現(xiàn)，92 d中有26 d的F1預(yù)報(bào)值最接近觀測(cè)值，而65 d為F2～F7中某個(gè)預(yù)報(bào)值更接近觀測(cè)值。

圖4 NAQPMS對(duì)江蘇省各城市PM2.5質(zhì)量濃度的7 d預(yù)測(cè)結(jié)果驗(yàn)證參數(shù)箱線圖

圖5 無(wú)錫市2013年夏季PM2.5質(zhì)量濃度預(yù)報(bào)和觀測(cè)結(jié)果的時(shí)間序列圖

3 結(jié)論

采用中國(guó)科學(xué)院大氣物理研究所NAQPMS空氣質(zhì)量模式，在建立的江蘇省級(jí)區(qū)域空氣質(zhì)量數(shù)值預(yù)報(bào)系統(tǒng)初步實(shí)現(xiàn)對(duì)PM2.5質(zhì)量濃度的未來(lái)24 h預(yù)報(bào)和未來(lái)7 d趨勢(shì)預(yù)測(cè)基礎(chǔ)上，利用2013年6—8月PM2.5質(zhì)量濃度觀測(cè)資料評(píng)估了系統(tǒng)對(duì)江蘇省PM2.5質(zhì)量濃度未來(lái)24 h預(yù)報(bào)和7 d潛勢(shì)預(yù)測(cè)效果。結(jié)果表明，該系統(tǒng)可合理預(yù)測(cè)出江蘇省2013年夏季PM2.5質(zhì)量濃度南北高、中部低的分布特征。對(duì)省屬13個(gè)地市的PM2.5質(zhì)量濃度未來(lái)24 h預(yù)報(bào)值達(dá)到“理想”水平。其中，南部沿江8市的PM2.5質(zhì)量濃度預(yù)報(bào)效果優(yōu)于蘇北5市。系統(tǒng)對(duì)PM2.5質(zhì)量濃度未來(lái)7 d趨勢(shì)預(yù)測(cè)效果達(dá)到“合理”水平，部分指標(biāo)接近“理想”水平，可較好地為污染預(yù)警應(yīng)急減排提供及時(shí)有效的技術(shù)支持。不同時(shí)段預(yù)測(cè)預(yù)報(bào)值區(qū)間對(duì)江蘇省13個(gè)城市的平均預(yù)報(bào)命中率可達(dá)50%，最高可達(dá)70%。今后可以通過(guò)提高模式分辨率、精細(xì)排放源清單(特別是本地源清單)、資料同化以及構(gòu)建集合預(yù)報(bào)等方法進(jìn)一步提高預(yù)報(bào)效果。

［1］GB 3095—2012 環(huán)境空氣質(zhì)量標(biāo)準(zhǔn)［S］．

［2］國(guó)務(wù)院辦公廳．大氣污染防治行動(dòng)計(jì)劃［R］．北京:國(guó)務(wù)院辦公廳，2013．

［3］王自發(fā)，謝付瑩，王喜全，等．嵌套網(wǎng)格空氣質(zhì)量預(yù)報(bào)模式系統(tǒng)的發(fā)展與應(yīng)用［J］．大氣科學(xué)，2006，30(5):778-790．

［4］Li J，Wang Z F，Zhuang G，et al．Mixing of Asian mineral dust with anthropogenic pollutants over East Asia:a model case study of a super-duststorm in March 2010［J］．Atmos Chem Phys，2012，12:7 591-7 607．

［5］王自發(fā)，吳其重，Gbaguidi A，等．北京空氣質(zhì)量多模式集成預(yù)報(bào)系統(tǒng)的建立及初步應(yīng)用［J］．南京信息工程大學(xué)學(xué)報(bào)(自然科學(xué)版)，2009，1(1):19-26．

［6］吳其重，王自發(fā)，徐文帥，等．多模式模擬評(píng)估奧運(yùn)賽事期間可吸入顆粒物減排效果［J］．環(huán)境科學(xué)學(xué)報(bào)，2010，30(9):1 739-1 748．

［7］王茜，伏晴艷，王自發(fā)，等．集合數(shù)值預(yù)報(bào)系統(tǒng)在上海市空氣質(zhì)量預(yù)測(cè)預(yù)報(bào)中的應(yīng)用研究［J］．環(huán)境監(jiān)控與預(yù)警，2010，2(4):1-6，11．

［8］陳煥盛，王自發(fā)，吳其重，等．空氣質(zhì)量多模式系統(tǒng)在廣州應(yīng)用及對(duì)PM10預(yù)報(bào)效果評(píng)估［J］．氣候與環(huán)境研究，2013，18(4):427-435．

［9］Kurokawa J，Ohara T，Morikawa T，et al．Emissions of air pollutants and greenhouse gases over Asian regions during 2000-2008:Regional Emission inventory in Asia(REAS)version 2［J］．Atmospheric Chemistry and Physics Discussions，2013，13(4):10 049-10 123．

［10］Guenther A，Hewitt N，Erickson D，et al．A global model of natural volatile organic compound emissions［J］．Journal of Geophysical Research，1995，100(D5):8 873-8 892．

［11］楊文夷，李杰，朱莉莉，等．我國(guó)空氣污染物人為源排放清單對(duì)比［J］．環(huán)境科學(xué)研究，2013，26(7):703-711．

［12］全國(guó)城市空氣質(zhì)量實(shí)時(shí)發(fā)布平臺(tái)．http://113.108．142.147:20035/emcpublish/．

［13］Boylan J W，Russell A G．PM and light extinction model performance metrics，goals，and criteria for threedimensional airquality models［J］．Atmospheric Environment，2006，40(23):4 946-4 959．

［14］E87950．World Health Organization Air Quality Guidelines Global Update［S］．

［15］MODIS Land ProductSubsets．http://lance-modis．eosdis．nasa．gov/imagery/subsets/?project=fas ＆subset=FAS_China4.2013164．terra．1km．

［16］Heather S， Kirk R， Sharon P．Compliation and interpretation ofphotochemicalmodelperformance statistics published between 2006 and 2012［J］．Atmospheric Environment，2012，60:124-139．