新鄉(xiāng)市空氣質(zhì)量數(shù)值預(yù)報(bào)系統(tǒng)及重污染過程模擬分析

馬 琳，張穩(wěn)定，魏 巍，晏平仲，3

（1.新鄉(xiāng)市環(huán)境保護(hù)監(jiān)測(cè)站，河南新鄉(xiāng)453003；2.中科三清科技有限公司，北京100029；3.中國(guó)科學(xué)院大氣物理研究所，北京100029）

自2000年開始，國(guó)家組織環(huán)境保護(hù)重點(diǎn)城市開展了城市環(huán)境空氣質(zhì)量日?qǐng)?bào)和預(yù)報(bào)工作。近幾年來，隨著大氣環(huán)境污染問題日漸突出，環(huán)境空氣質(zhì)量預(yù)報(bào)逐漸從幕后走到了臺(tái)前，開始為國(guó)家制定經(jīng)濟(jì)發(fā)展規(guī)劃以及人民群眾日常工作和生產(chǎn)生活提供指導(dǎo)。2013年12月，我國(guó)中東部出現(xiàn)長(zhǎng)時(shí)間大范圍高強(qiáng)度的細(xì)顆粒物污染事件，全國(guó)20多個(gè)省份、100多個(gè)城市顆粒物污染達(dá)重度至嚴(yán)重污染水平，引起了公眾的高度關(guān)注。近年來，我國(guó)中東部地區(qū)類似的重污染事件日益增多，呈現(xiàn)污染程度重、持續(xù)時(shí)間長(zhǎng)、污染面積大等特征。2014年，中國(guó)環(huán)境監(jiān)測(cè)總站印發(fā)《環(huán)境空氣質(zhì)量預(yù)報(bào)預(yù)警業(yè)務(wù)工作指南（暫行）》[1]，用于指導(dǎo)各省、自治區(qū)、直轄市及相關(guān)城市環(huán)境監(jiān)測(cè)中心（站）的空氣質(zhì)量預(yù)報(bào)預(yù)警業(yè)務(wù)工作。對(duì)于預(yù)報(bào)方法，可選擇數(shù)值預(yù)報(bào)、統(tǒng)計(jì)預(yù)報(bào)方法或2者的集合預(yù)報(bào)。以第3代空氣質(zhì)量模型嵌套網(wǎng)格空氣質(zhì)量數(shù)值模式系統(tǒng)（NAQPMS）、多尺度空氣質(zhì)量模式（CMAQ）、復(fù)合型空氣質(zhì)量模型（CAMx）以及氣象研究和預(yù)報(bào)化學(xué)耦合模型（WRFCHEM）等數(shù)值預(yù)報(bào)模型為主要工具，在科學(xué)研究[2-5]和空氣質(zhì)量預(yù)報(bào)領(lǐng)域[6-9]發(fā)揮出重要的作用。

2015年，河南省環(huán)境保護(hù)廳下發(fā)《河南省環(huán)境保護(hù)廳關(guān)于印發(fā)河南省環(huán)境空氣質(zhì)量預(yù)報(bào)預(yù)警試行方案的通知》[10]，河南省各市陸續(xù)建立環(huán)境空氣質(zhì)量預(yù)報(bào)預(yù)警業(yè)務(wù)平臺(tái)，為響應(yīng)國(guó)家大氣污染防治政策，積極應(yīng)對(duì)空氣污染問題。2016年，新鄉(xiāng)市環(huán)保局啟動(dòng)了空氣質(zhì)量預(yù)報(bào)預(yù)警平臺(tái)的建設(shè)工作，選用NAQPMS和CMAQ模型作為空氣質(zhì)量數(shù)值預(yù)報(bào)的工具。項(xiàng)目建設(shè)單位基于多個(gè)預(yù)報(bào)預(yù)警系統(tǒng)項(xiàng)目的建設(shè)經(jīng)驗(yàn)，結(jié)合新鄉(xiāng)實(shí)際情況，設(shè)計(jì)開發(fā)出適合新鄉(xiāng)本地的空氣質(zhì)量數(shù)值預(yù)報(bào)系統(tǒng)，輔助預(yù)報(bào)人員做好新鄉(xiāng)市空氣質(zhì)量預(yù)報(bào)預(yù)警工作。

1 系統(tǒng)的設(shè)計(jì)

空氣質(zhì)量數(shù)值預(yù)報(bào)系統(tǒng)是新鄉(xiāng)市空氣質(zhì)量調(diào)控綜合決策支撐平臺(tái)的重要組成部分，是開展空氣質(zhì)量預(yù)報(bào)預(yù)警業(yè)務(wù)的核心。該系統(tǒng)以大氣污染物排放清單和氣象預(yù)報(bào)數(shù)據(jù)為基礎(chǔ)，通過三維歐拉空氣質(zhì)量模型計(jì)算出未來一段時(shí)間內(nèi)大氣中多種污染物的濃度，結(jié)合后處理和業(yè)務(wù)處理流程，為空氣質(zhì)量預(yù)報(bào)預(yù)警工作平臺(tái)提供各類預(yù)報(bào)產(chǎn)品。

本研究選擇國(guó)內(nèi)先進(jìn)的中國(guó)科學(xué)院大氣物理研究所自主研發(fā)的NAQPMS以及國(guó)際主流的美國(guó)環(huán)保署Model-3-CMAQ模式作為空氣質(zhì)量數(shù)值預(yù)報(bào)核心模式，采用廣泛應(yīng)用的中尺度氣象模式WRF（weather research and forecast）提供輸入氣象場(chǎng)，耦合清華大學(xué)研究建立的用于空氣質(zhì)量模式的網(wǎng)格化全國(guó)污染源清單，基于模式不確定性分析評(píng)估各輸入數(shù)據(jù)不確定性對(duì)模式預(yù)報(bào)誤差的影響，采用新鄉(xiāng)市及周邊區(qū)域空氣質(zhì)量集合預(yù)報(bào)優(yōu)化方法，利用包括多元回歸、偏差訂正、權(quán)重因子和神經(jīng)網(wǎng)絡(luò)等方法對(duì)空氣質(zhì)量集合預(yù)報(bào)結(jié)果進(jìn)行統(tǒng)計(jì)集成，滿足新鄉(xiāng)市環(huán)境空氣質(zhì)量預(yù)報(bào)預(yù)警業(yè)務(wù)準(zhǔn)確性需求。

1.1 系統(tǒng)框架

新鄉(xiāng)市空氣質(zhì)量數(shù)值模式系統(tǒng)框架圖如圖1。

圖1 新鄉(xiāng)市環(huán)境空氣質(zhì)量多模式預(yù)報(bào)預(yù)警平臺(tái)系統(tǒng)流程圖Fig.1 System flowchart of environmental air quality multi-model forecast and warning platform in Xinxiang city

該系統(tǒng)采用多模式集合預(yù)報(bào)系統(tǒng)，通過統(tǒng)一的區(qū)域設(shè)置、排放清單、氣象場(chǎng)驅(qū)動(dòng)，空氣質(zhì)量數(shù)值模式采用國(guó)內(nèi)應(yīng)用成熟、業(yè)務(wù)化程度高的NAQPMS、CMAQ等，并提供污染資料及氣象資料同化功能。系統(tǒng)可利用多元回歸、偏差訂正、神經(jīng)網(wǎng)絡(luò)等方法對(duì)空氣質(zhì)量集合預(yù)報(bào)結(jié)果進(jìn)行統(tǒng)計(jì)集成，實(shí)現(xiàn)多模式的集合預(yù)報(bào)。

相對(duì)于單模式預(yù)報(bào)，多模式集合預(yù)報(bào)可以集成不同模式的優(yōu)點(diǎn)，使得通過合適的集成方法來提高確定性預(yù)報(bào)的準(zhǔn)確率成為可能。

本系統(tǒng)設(shè)計(jì)符合《環(huán)境空氣質(zhì)量預(yù)報(bào)預(yù)警方法技術(shù)指南》[11]、《全國(guó)環(huán)境空氣質(zhì)量預(yù)報(bào)信息交換技術(shù)指南（試行）》[12]、《環(huán)境空氣質(zhì)量可視化預(yù)報(bào)會(huì)商技術(shù)指南（試行）》[13]、《環(huán)境空氣質(zhì)量數(shù)值模式源清單技術(shù)指南（試行）》[14]、《大氣細(xì)顆粒物（PM2.5）源排放清單編制技術(shù)指南（試行）》（征求意見稿）[15]等方面要求。

1.2 數(shù)值模式簡(jiǎn)介

1.2.1 NAQPMS模式

NAQPMS的設(shè)計(jì)是以我國(guó)當(dāng)前計(jì)算硬件條件和業(yè)務(wù)水平為出發(fā)點(diǎn)，結(jié)合我國(guó)城市群大氣復(fù)合污染的排放、輸送、演變特點(diǎn)，綜合評(píng)估多個(gè)有代表性的數(shù)值模式。NAQPMS主要由氣象處理、排放源處理、空氣質(zhì)量模式及模式輸出等4個(gè)主要部分構(gòu)成。NAQPMS采用開放式氣象驅(qū)動(dòng)場(chǎng)，可利用MM5、WRF等中尺度氣象模式輸出的氣象要素場(chǎng)作為模式的動(dòng)力驅(qū)動(dòng)，而中尺度氣象模式預(yù)報(bào)所需的初始條件和邊界條件可由美國(guó)國(guó)家環(huán)境預(yù)報(bào)中心（NCEP）-美國(guó)國(guó)家大氣研究中心（NCAR）再分析數(shù)據(jù)、歐洲數(shù)值預(yù)報(bào)中心的ECMWF數(shù)據(jù)等全球數(shù)據(jù)提供。結(jié)合SMOKE模型實(shí)時(shí)輸出的排放源，NAQPMS可以對(duì)大氣中主要化學(xué)成分的分布狀況、輸送態(tài)勢(shì)、沉降特征進(jìn)行數(shù)值模擬，從而使得模式系統(tǒng)能夠合理反映大氣化學(xué)成分在輸送過程中的物理化學(xué)特性變化。

1.2.2 CMAQ模式

Model-3-CMAQ系統(tǒng)是美國(guó)環(huán)境保護(hù)局開發(fā)的第3代空氣質(zhì)量模式，以“一個(gè)大氣”為設(shè)計(jì)理念，擬將所有的大氣問題均考慮進(jìn)模式之中。Model-3模式系統(tǒng)由中尺度氣象模式、排放模式模式以及通用多尺度模式（CMAQ）等3大模塊組成。中尺度氣象模式采用MM5和WRF，SMOKE排放模式提供污染源的排放清單，估算污染源位置和產(chǎn)生量隨時(shí)間的變化；CMAQ為空氣質(zhì)量系統(tǒng)的核心部分，設(shè)計(jì)為多重網(wǎng)格嵌套模式，可以模擬多種污染物的輸送和轉(zhuǎn)化過程。

1.3 數(shù)值預(yù)報(bào)方案

1.3.1 模式區(qū)域和預(yù)報(bào)時(shí)間

對(duì)多模式預(yù)報(bào)系統(tǒng)的區(qū)域網(wǎng)格進(jìn)行了統(tǒng)一的設(shè)置，本系統(tǒng)中模式計(jì)算采用三重嵌套區(qū)域設(shè)置。第1區(qū)域?yàn)橹袊?guó)中東部區(qū)域，水平分辨率為27 km；第2區(qū)域?yàn)楹幽鲜〖捌渲苓厖^(qū)域，水平分辨率為9 km；第3區(qū)域?yàn)樾锣l(xiāng)市全境及周邊，水平分辨率為新鄉(xiāng)市區(qū)為1 km；其他區(qū)域?yàn)? km。模式計(jì)算垂直范圍從地面到20 km高度，垂直分層不少于20層，邊界層內(nèi)分為10層。

模式以前日北京時(shí)間20時(shí)為起點(diǎn)，計(jì)算未來168 h的氣象場(chǎng)和污染物濃度場(chǎng)，以實(shí)現(xiàn)未來3 d可用的區(qū)域空氣質(zhì)量預(yù)報(bào)和未來7 d可供參考的區(qū)域污染趨勢(shì)預(yù)測(cè)，預(yù)報(bào)輸出結(jié)果的時(shí)間分辨率為1 h。

1.3.2 氣象場(chǎng)預(yù)報(bào)及設(shè)置

空氣質(zhì)量模式的氣象驅(qū)動(dòng)場(chǎng)由NCEP、NCAR等科研機(jī)構(gòu)和大學(xué)聯(lián)合開發(fā)的新一代中尺度氣象模式WRF提供，采用NCEP的數(shù)值天氣預(yù)報(bào)中心GFS數(shù)據(jù)集全球預(yù)報(bào)分析資料作為WRF模式運(yùn)行的初始及邊界條件，輸出頻率為每小時(shí)1次。WRF模式垂直方向上采用地形跟隨質(zhì)量坐標(biāo)系，垂直分層不少20層。

1.3.3 排放清單及其處理

研究的污染物主要涉及無機(jī)污染物和有機(jī)污染物兩大類。無機(jī)污染物包括顆粒物PM10、PM2.5，氣態(tài)污染物SO2、NOx、CO和NH3等，其中PM10和PM2.5特指由污染源直接排放的一次顆粒物；有機(jī)污染物主要指VOC。

大氣污染物排放清單以實(shí)測(cè)為基礎(chǔ)，結(jié)合重點(diǎn)污染源核查等專項(xiàng)工作，建立了與行業(yè)或通用設(shè)備相關(guān)聯(lián)的本地污染物排放系數(shù)，同時(shí)，還借鑒了國(guó)內(nèi)《工業(yè)污染物產(chǎn)生和排放系數(shù)于冊(cè)》中部分排放系數(shù)、歐美等發(fā)達(dá)國(guó)家清單統(tǒng)計(jì)的規(guī)范和方法（如美國(guó)EPA的Air Chief、NEI和NIF、歐盟的CORINAIR、英國(guó)的NAEI 96年排放清單等）等。

1.3.4 模式參數(shù)化方案設(shè)置

本研究首先根據(jù)前期研究經(jīng)驗(yàn)初步預(yù)設(shè)多模式集合預(yù)報(bào)系統(tǒng)4個(gè)核心模式的物理化學(xué)方案，通過模擬評(píng)估進(jìn)一步優(yōu)化方案參數(shù)，具體方案見表1。

表1 模式應(yīng)用方案Tab.1 Application schemes of models

1.3.5 模式自動(dòng)化運(yùn)行技術(shù)

模式自動(dòng)化運(yùn)行是系統(tǒng)的技術(shù)核心，本研究中基于Linux SHELL腳本和Fortran程序?qū)崿F(xiàn)多模式系統(tǒng)的自動(dòng)化運(yùn)行，自動(dòng)生成預(yù)報(bào)所需圖片和文本文件，每日定期備份關(guān)鍵數(shù)據(jù)，并且清理模式計(jì)算產(chǎn)生的臨時(shí)文件，滿足自動(dòng)、穩(wěn)定、無人值守是空氣質(zhì)量預(yù)報(bào)系統(tǒng)的基本要求。

以NAQPMS模式系統(tǒng)的自動(dòng)化為例，該系統(tǒng)由一個(gè)主SHELL腳本整體控制，負(fù)責(zé)調(diào)用控制NAQPMS各個(gè)部分的子SHELL腳本，每個(gè)子SHELL腳本控制相應(yīng)模塊的自動(dòng)運(yùn)行。自動(dòng)化運(yùn)行的時(shí)間由主SHELL腳本產(chǎn)生，每天0時(shí)，系統(tǒng)第1個(gè)子部分在相應(yīng)SHELL腳本的控制下自動(dòng)下載NCEP的全球預(yù)報(bào)資料，然后依次運(yùn)行WRF模式的WPS模塊和WRF模塊SHELL腳本，接著第2個(gè)子部分運(yùn)行wrf2apops氣象-化學(xué)接口SHELL腳本，將結(jié)果輸入空氣質(zhì)量模式，運(yùn)行real-NAQPMS的SHELL腳本進(jìn)行污染場(chǎng)的計(jì)算，完成后用后處理SHELL腳本將氣象場(chǎng)和污染場(chǎng)的預(yù)報(bào)結(jié)果繪圖，生成多種預(yù)報(bào)產(chǎn)品，并且調(diào)用Fortran程序計(jì)算AQI等表格，生成網(wǎng)頁發(fā)布。此系統(tǒng)充分利用了0—7時(shí)這段相對(duì)空余的時(shí)間運(yùn)行，呈現(xiàn)未來幾天的詳細(xì)天氣走勢(shì)以及不同高度污染物的分布，為空氣污染的預(yù)報(bào)預(yù)警和分析預(yù)報(bào)提供有力的支持。

系統(tǒng)每日0時(shí)自動(dòng)下載全球氣象預(yù)報(bào)場(chǎng)數(shù)據(jù)，整個(gè)計(jì)算時(shí)間約為4 h，保證預(yù)報(bào)結(jié)果在每日7時(shí)之前輸入預(yù)報(bào)平臺(tái)，滿足業(yè)務(wù)預(yù)報(bào)時(shí)效性需求。

2 重污染過程模擬分析

2.1 重污染過程回顧

2016年12月16—20日，受采暖季污染物排放增強(qiáng)和不利氣象條件的綜合影響，我國(guó)華北地區(qū)發(fā)生了大范圍的細(xì)顆粒物持續(xù)重污染過程，20多個(gè)城市同步啟動(dòng)了重污染紅色預(yù)警。基于新鄉(xiāng)市空氣質(zhì)量數(shù)值模式業(yè)務(wù)預(yù)報(bào)系統(tǒng)和多種觀測(cè)資料，本文中初步分析了新鄉(xiāng)市本次污染過程的成因，評(píng)估了新鄉(xiāng)市數(shù)值模式對(duì)重污染過程的預(yù)報(bào)效果。

圖2展示了新鄉(xiāng)市觀測(cè)（柱狀）的AQI和細(xì)顆粒物PM2.5的時(shí)間序列圖。細(xì)顆粒物PM2.5濃度從15日開始累積，17日超過了輕度污染等級(jí)，18日濃度飆升到300 μg/m3，19日達(dá)到最高峰（PM2.5濃度超過了450 μg/m3），20日擴(kuò)散條件改善，濃度開始下降，22日空氣污染等級(jí)降到良。對(duì)比NAQPMS和CMAQ的24 h預(yù)報(bào)結(jié)果，2個(gè)模式很好的預(yù)報(bào)出了此次過程PM2.5的變化趨勢(shì)，但在污染較重的幾天（18—21日）同觀測(cè)值相比，有一定偏差。NAQPMS和CMAQ24 h預(yù)報(bào)的PM2.5值的相關(guān)系數(shù)r分別為0.81和0.63，2個(gè)模式24 h預(yù)報(bào)時(shí)效同觀測(cè)的相關(guān)性較強(qiáng)。

圖2 AQI和PM2.5濃度模擬和觀測(cè)對(duì)比圖Fig.2 Comparison of predicted and observed for AQI and PM2.5

2016年12月17—22日全國(guó)空氣質(zhì)量監(jiān)測(cè)站點(diǎn)PM2.5空間分布顯示，17日PM2.5主要集聚在長(zhǎng)三角及京津冀地區(qū)，河南全省污染等級(jí)不高；18日華北地區(qū)污染急劇加重，河南全省市達(dá)到重污染等級(jí)；19—20日重污染等級(jí)地區(qū)擴(kuò)大至華北至華中地區(qū)；21—22日在偏北西風(fēng)的影響下，顆粒物往東部沿海地區(qū)輸送，河南省的污染減輕。

2.2 重污染過程的氣象條件分析

2016年12月17—20日，華北地區(qū)高空主要受緯向型環(huán)流控制，冷空氣活動(dòng)不活躍，天氣情況較為穩(wěn)定。17日，新鄉(xiāng)上游地區(qū)有弱短波槽，并不斷東移；19—20日，新鄉(xiāng)市受弱槽前偏南氣流影響，空中增溫增濕，在此期間，起源于內(nèi)蒙中北部的槽的不斷發(fā)展加深；至22日，新鄉(xiāng)市受槽后強(qiáng)大西北氣流控制，有利于污染物的擴(kuò)散。

地面上，12月17日8時(shí)，新鄉(xiāng)受華北小低壓控制，有利于污染物的輻合；19日8時(shí)，新鄉(xiāng)市逐漸由前一日的高后倒槽控制轉(zhuǎn)變?yōu)槿醺邏旱撞靠刂?，且處于北部弱高前輻合線的南側(cè)，有弱輻合，地面風(fēng)速較小，污染持續(xù)加重；21日，新鄉(xiāng)市再次為倒槽控制，22日，受蒙古冷高壓迅速南下影響，新鄉(xiāng)市空氣質(zhì)量迅速轉(zhuǎn)好。

圖3為此次重污染過程每日8時(shí)鄭州的探空曲線圖。由圖可見，此次污染地區(qū)的深厚逆溫層，亦不利于污染物擴(kuò)散。采用和新鄉(xiāng)市距離較近的鄭州市的探空數(shù)據(jù)可看到，從17日8時(shí)起，鄭州地區(qū)自地面至高空850 hPa存在逆溫層，并且一直持續(xù)到21日20時(shí)，直到22日，受近地面層強(qiáng)冷平流的影響，逆溫層才逐漸消失。鄭州地區(qū)水汽主要聚積在近地面層，垂直方向上呈上干下濕分布，有利于污染物累積。

綜上所述，此次污染過程中，主要受緯向型環(huán)流控制，冷空氣不活躍，天氣較為穩(wěn)定；地面多受弱氣壓場(chǎng)影響，水平風(fēng)速較??；且垂直方向上水汽呈上干下濕分布，并有深厚逆溫層存在，導(dǎo)致水平和垂直擴(kuò)散條件均較差，從18日開始，新鄉(xiāng)地區(qū)受此天氣形勢(shì)影響，污染物的濃度逐漸累積，直至22日冷空氣大舉南下，此次污染過程結(jié)束。

圖3 2016年12月17—22日8時(shí)鄭州（57083）探空曲線Fig.3 Sounding curve in Zhengzhou airport at 8∶00 from December 17 to 22，2016

2.3 新鄉(xiāng)市預(yù)報(bào)系統(tǒng)的模擬評(píng)估及偏差原因分析

預(yù)報(bào)系統(tǒng)不但預(yù)報(bào)出了新鄉(xiāng)本地的污染過程，對(duì)整個(gè)區(qū)域的顆粒物污染也有較好的反應(yīng)。2016年12月17—22日8時(shí)NAQPMS 24 h時(shí)效預(yù)報(bào)PM2.5空間分布顯示，17日PM2.5預(yù)報(bào)高值主要集聚在長(zhǎng)三角及京津冀地區(qū)，河南省和新鄉(xiāng)市污染等級(jí)不高；18日華北地區(qū)污染開始加重，河南省和新鄉(xiāng)市達(dá)到重污染等級(jí)；19—20日重污染等級(jí)地區(qū)逐漸擴(kuò)大；21—22日顆粒物緩慢往北輸送，河南省的污染有所減輕。整個(gè)過程N(yùn)AQPMS預(yù)報(bào)結(jié)果和監(jiān)測(cè)結(jié)果基本一致。

3 結(jié)論

新鄉(xiāng)市空氣質(zhì)量數(shù)值業(yè)務(wù)系統(tǒng)分析模式對(duì)重污染過程的預(yù)報(bào)效果說明，預(yù)報(bào)系統(tǒng)較好地預(yù)報(bào)出了污染過程，尤其是NAQPMS模式，較為準(zhǔn)確地預(yù)報(bào)出了污染時(shí)段和污染物的濃度水平，CMAQ也準(zhǔn)確地預(yù)報(bào)出了污染趨勢(shì)，但是對(duì)于細(xì)顆粒物的峰值濃度預(yù)報(bào)不足。這可能是由于NAQPMS的物理化學(xué)機(jī)制更適合我國(guó)實(shí)際污染過程，而CMAQ模式在本地更好的應(yīng)用可能還需要更多的本地化工作。

致謝：在論文撰寫期間以及后期修改過程中得到了中國(guó)科學(xué)院大氣物理研究所吳劍斌老師的熱心指導(dǎo)和幫助，在此表示衷心感謝。

（

）：

[1]中國(guó)環(huán)境監(jiān)測(cè)總站.關(guān)于印發(fā)《環(huán)境空氣質(zhì)量預(yù)報(bào)預(yù)警業(yè)務(wù)工作指南（暫行）》的通知[EB/OL].（2014-03-11）[2017-08-01].http://www.cnemc.cn/publish/117/news/news_40461.htm l.

[2]WU J B，WANG Z F，WANG Q，et al.Development of an on-line source-tagged model for sulfate，nitrate and ammonium:a modeling study for highly polluted periods in Shanghai，China[J].Environmental Pollution，2017，221:168-179.

[3]WANG L T，WEI Z，WEI W，et al.Source apportionment of PM2.5in top polluted cities in Hebei，China using the CMAQ model[J].Atmospheric Environment，2015，122:723-736.

[4]BONYOUNG K，ELADIO K，GREG Y.1.5-Dimensional volatility basis set approach for modeling organic aerosol in CAMx and CMAQ[J].Atmospheric Environment，2014，95:158-164.

[5]LI N，HE Q Y，TIE X X，et al.Quantifying sources of elemental carbon over the Guanzhong Basin of China:a consistent network of measurements and WRF-Chem modeling[J].Environmental Pollution，2016，214:86-93.

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

[7]徐文帥，李云婷，吳其重，等.AQI標(biāo)準(zhǔn)下北京市空氣質(zhì)量數(shù)值預(yù)報(bào)系統(tǒng)及其在重大活動(dòng)保障中的應(yīng)用[J].中國(guó)環(huán)境監(jiān)測(cè)，2016，32（2）：11-19.

[8]程念亮，李紅霞，孟凡，等.山東省空氣質(zhì)量預(yù)報(bào)平臺(tái)設(shè)計(jì)及其預(yù)報(bào)效果評(píng)估[J].環(huán)境污染與防治，2015，37（9）：92-99.

[9]蘇樅樅.應(yīng)用嵌套網(wǎng)格模式系統(tǒng)對(duì)沈陽市環(huán)境空氣質(zhì)量的數(shù)值模擬研究[J].北方環(huán)境，2013，29（4）：58-61.

[10]河南省環(huán)境保護(hù)廳.河南省環(huán)境保護(hù)廳關(guān)于印發(fā)河南省環(huán)境空氣質(zhì)量預(yù)報(bào)預(yù)警試行方案的通知[EB/OL].（2015-04-16）[2017-08-01].http://www.hnep.gov.cn/tabid/435/InfoID/15739/Default.aspx.

[11]中國(guó)環(huán)境監(jiān)測(cè)總站，環(huán)境空氣質(zhì)量預(yù)報(bào)預(yù)警方法技術(shù)指南[M].北京：中國(guó)環(huán)境出版社，2014.

[12]中國(guó)環(huán)境監(jiān)測(cè)總站.環(huán)境空氣質(zhì)量預(yù)報(bào)信息交換技術(shù)指南（試行）[EB/OL].（2015-03-09）[2017-08-01].http://www.cnemc.cn/publish/117/news/news_44295.htm l.

[13]中國(guó)環(huán)境監(jiān)測(cè)總站.環(huán)境空氣質(zhì)量可視化預(yù)報(bào)會(huì)商技術(shù)指南（試行）[EB/OL].（2015-03-09）[2017-08-01].http://www.cnemc.cn/publish/117/news/news_44299.htm l.

[14]中國(guó)環(huán)境監(jiān)測(cè)總站.環(huán)境空氣質(zhì)量數(shù)值預(yù)報(bào)模式源清單技術(shù)指南（試行）[EB/OL].（2015-03-09）[2017-08-01].http://www.cnemc.cn/publish/totalWebSite/news/news_44325.htm l.

[15]環(huán)境保護(hù)部辦公廳.關(guān)于征求《大氣細(xì)顆粒物（PM2.5）源排放清單編制技術(shù)指南（試行）》（征求意見稿）等8項(xiàng)環(huán)境保護(hù)部文件意見的函[EB/OL].（2014-01-17）[2017-08-01].http://www.zhb.gov.cn/gkm l/hbb/bgth/201401/t20140124_266910.htm l.