錢俊龍，牛文勝，查書瑤，王璐璐，陳瀟瀟，周彬

(1.無錫市氣象局，江蘇無錫214000；2.中國氣象局交通氣象重點開放實驗室，江蘇南京210000)

1 引 言

霾，是指大量極細微的干塵粒等均勻地浮游在空中，使水平能見度小于10 km，造成空氣普遍渾濁的天氣現象。統(tǒng)計表明，隨著社會經濟的發(fā)展，城市霧日在減少，霾日卻在增加[1-3]。霾天氣歸根結底是細顆粒物PM2.5(尤其是PM1.0)的污染事件[4-6]，PM2.5主要來源為燃煤、揚塵、機動車排放、建筑塵、生物質燃燒、二次硫酸鹽和硝酸鹽及有機物[7]。研究發(fā)現，PM2.5對老年人和兒童危害最大，長期暴露于燃燒產生的細顆粒，可導致心肺疾病和肺癌死亡率增加[8-10]。

在日益嚴重的大氣污染形勢下，國家頒布了《大氣十條》以期治理大氣污染，越來越多的專家學者也投入到大氣環(huán)境方向研究中。Chen等[11]分析了1960—2012年間華北地區(qū)霾的變化，發(fā)現北風分量的減弱、對流層底層逆溫增強、中層東亞槽的減弱等因素導致了該地區(qū)1984—1992年霾日的增加(與1993—2001年對比)。Qu等[12]研究發(fā)現華東地區(qū)能見度每十年降低0.68 km(1973—2012年)，這可能與每年寒潮發(fā)生次數減少有關。Tsai等[13]通過后向軌跡、天氣圖分析、區(qū)域沙塵模式等方法分析了ACE-Asia項目獲得的氣溶膠數據，研究東亞沙塵傳輸機制，發(fā)現沙塵傳輸距離與源地的天氣形勢密切相關，源地位于槽前更有利于沙塵傳輸。戴竹君等[14]對江蘇秋冬季重度霾的環(huán)流背景、邊界層特征、熱力條件、動力條件及氣流軌跡進行了分析，發(fā)現低層輻合下沉區(qū)域霾常加重，冷鋒前部型氣流輸送的氣溶膠主要來自西北-華北地區(qū)。陳一娜等[15]研究發(fā)現，氣溶膠散射作用占環(huán)境總消光作用的94%以上，在夏秋季相對濕度可以使氣溶膠的散射能力提升70%～80%。張蓬勃等[16]在研究中指出，江蘇冬季霾日數對氣溫增暖存在可能的顯著正響應。于興娜等[17]通過分析顆粒物濃度、氣溶膠細粒子譜分布等觀測數據，發(fā)現0.5～2.0 μm細粒子及相對濕度是導致南京北郊秋、冬季大氣能見度下降的主要因素。薛文博等[18]利用CAMx模式基于顆粒物來源追蹤技術(PSAT)，模擬了全國PM2.5的跨區(qū)域輸送規(guī)律，指出江蘇省PM2.5年均濃度受省外源貢獻超過45%。

數值預報正成為天氣預報、空氣質量預報中的主要參考資料，數值模式能夠定量、客觀地反映要素場的時空變化。李鋒等[19]運用WRF-CMAQ模擬了2013年12月長江三角洲地區(qū)一次重霾污染過程，探究重霾天氣下PM2.5的時空分布，并定量研究了外來輸送對長江三角洲地區(qū)PM2.5的貢獻。張小玲等[20]利用WRF-Chem對2013年1月華北平原一次霧霾過程進行模擬，模式對天氣系統(tǒng)演變和PM2.5質量濃度的空間分布及高濃度持續(xù)時間、消散減弱等過程做出了較好的預報。吳洛林等[21]利用WRF-Chem模式研究發(fā)現，不利氣象條件及本地排放源是造成肇慶冬季重污染過程發(fā)生的主要原因，肇慶境內的主要污染物輸送通道呈東南-西北走向，外地輸入的污染物也通過此通道影響肇慶的空氣質量。

每年冬季，江蘇地區(qū)都會迎來一輪又一輪污染過程，雖然各地政府都在加強環(huán)保監(jiān)督，大力減排，但是關于污染來自本地排放還是外來輸送的爭論卻是一直存在。2016年12月中下旬(14—24日)，江蘇省遭遇了一次長時間大范圍霾天氣過程，本文利用地面和高空氣象觀測資料、環(huán)境監(jiān)測資料結合數值模式(WRF-Chem)對此輪過程進行分析，定量研究外來輸送對江蘇省PM2.5的污染貢獻。

2 資料與方法

2.1 數據來源

本文采用的地面氣象觀測資料來自江蘇省13個市基本觀測站，PM2.5濃度資料來自環(huán)保部72個國控點監(jiān)測數據，氣象觀測資料和PM2.5資料時間分辨率均為1 h。NCEP再分析資料水平分辨率為1°×1°，垂直分層為26層，每天4個時次。數值模式所用的排放源資料來自清華大學MEIC項目[22]公開的2010年12月排放清單，該清單包含了農業(yè)、工業(yè)、火電廠、交通以及居民生活等人為排放源。

2.2 數值模擬方案設計

WRF-Chem模式是美國最新發(fā)展的區(qū)域大氣動力-化學耦合模式，最大優(yōu)點是氣象模式與化學傳輸模式在時間和空間分辨率上完全耦合[23]。本文中的數值模擬采用WRF-Chem模式V3.6版本，模擬采用兩層嵌套，最內層網格距為9 km，垂直分層32層，排放源在底部5層。模式模擬范圍最內層覆蓋了江蘇全境 (圖1)。模擬時間為2016年12月13日20:00—24日20:00(北京時間，下同)，時間分辨率為1 h。模式采用的主要參數化方案見表1。

表1 WRF-Chem參數化方案

為定量研究省內和省外排放源對江蘇PM2.5的污染貢獻，本文做了1組基礎實驗(CASE0)和6組對比實驗 (CASE1～6)，CACE1～6實驗中分別關閉了江蘇省、京津冀地區(qū)、山東省、山西省、河南省、安徽省的排放源。對比實驗與基礎實驗之差看作是該省排放源對江蘇的污染貢獻。

3 分析與討論

冬季，受北方采暖污染排放影響，來自內陸的西北風會將山東、河北一帶的污染物輸送至江蘇，造成PM2.5濃度增長[24-25]。

3.1 污染過程回顧

2016年12月中下旬江蘇省經歷了一輪長時間大范圍的重污染天氣過程，期間有兩次不同強度的冷空氣過程活動并伴隨PM2.5濃度顯著增長。第一次過程(15—17日)，前期上游污染輸送造成全省PM2.5濃度迅速增長，后期靜穩(wěn)的邊界層條件及風向的轉變(污染回流)共同造成此次重污染過程；第二次過程(22—23日)，平均風力較大，上游污染輸送是造成這次嚴重污染的主因[26]。本文主要討論15—17日和22—23日這兩次污染過程。

從高空形勢來看(圖2a、2b)，15日受高空槽后西北氣流控制，下沉氣流不利于污染物的垂直擴散；22日500 hPa高空槽前西南氣流，上升氣流有利于污染物遠距離傳輸，850 hPa為西北氣流控制，有利于高空污染向地面輸送。從地面形勢看(圖2c、2d)，17日早晨地面受高壓中心控制，風力較低，風向由前期的偏北風轉偏南風，造成污染物回流。22日開始地面處于高壓前部西北風場中，大量上游污染物被輸送至江蘇。

3.2 數值模擬

本文利用WRF-Chem模式對14—24日過程進行數值模擬，并對模擬結果進行模擬觀測檢驗。江蘇省13個市的PM2.5模擬值與觀測值的相關系數在0.51～0.75，相關性均通過0.01顯著性檢驗，標準化平均偏差(NMB)在-0.40～0.06，模擬值總體略偏低，均方根差(RMSE)在39～50(徐州90)。圖3分別給出了江蘇省13市和無錫 (代表蘇南)、泰州(代表蘇中)、宿遷(代表蘇北)的PM2.5模擬/觀測散點分布。總體來說，模式能夠較好地反映江蘇省PM2.5的時空分布。

圖4為15日和22日宿遷、揚州、無錫三站上空PM2.5濃度疊加風場的垂直剖面圖，PM2.5濃度峰值出現時刻順序為：宿遷-揚州-無錫(自西北向東南)，峰值均出現在西北風場中，等到風向轉為偏北風時，在來自山東半島以及海上的清潔氣流稀釋下，峰值逐漸減弱。22日的高值區(qū)范圍強度明顯強于15日，重度污染層厚度達到900～1 500 m，且持續(xù)時間較長。

3.3 各省排放源貢獻

本文通過數值模擬實驗，獲得各省排放源對江蘇PM2.5的污染貢獻(圖5和表2)?？煽闯?，第一次過程中江蘇省內排放源貢獻率在23%(連云港)～79%(蘇州)，自北向南逐漸增加，這與省內排放源分布相關，蘇北外來輸送貢獻率超過50%，蘇南仍以本地排放源污染為主。外來輸送山東省貢獻率最高，其次為京津冀地區(qū)，而其余周邊省份的貢獻率明顯低于兩者，安徽省對江蘇西部城市 (南京、徐州、宿遷)的貢獻率較高。第二次過程中江蘇省內貢獻率在5%(連云港)～32%(蘇州)，山東省和京津冀地區(qū)排放源是主要污染來源，二者貢獻之和在44%(南京)～70%(淮安)。從輸送量絕對值來看(圖6)，15日宿遷、揚州、無錫的PM2.5日平均外來 輸 送 量 分 別 為 40.2 μg/m3、20.9 μg/m3、11.1 μg/m3，外來輸送自北向南遞減；22日外來輸送明顯高于15日，PM2.5日平均外來輸送量分別為105.9 μg/m3、83.1 μg/m3、64.8 μg/m3，宿遷、揚州、無錫先后迎來外來輸送的峰值。圖7為京津冀地區(qū)和山東省的排放源對江蘇省PM2.5貢獻，橫向比較發(fā)現，第二次過程中京津冀地區(qū)和山東省排放源貢獻率均高于第一次過程，這是由于第二次過程西北風力更大，同時高空形勢更有利于污染物的遠距離傳輸。縱向比較發(fā)現，兩次過程中山東省貢獻率均高于京津冀地區(qū)，山東省貢獻率自北向南遞減，而京津冀地區(qū)貢獻率無此特征(不明顯)。

表2 兩次過程中各省排放源對江蘇各市PM2.5貢獻量 單位：μg/m3。

4 結 論

(1)2016年12月中下旬江蘇省經歷一次持續(xù)性霾天氣過程，期間有兩次冷空氣過程活動并伴隨PM2.5濃度顯著增長。PM2.5濃度自西北向東南(宿遷-揚州-無錫)先后出現峰值，峰值均出現在西北風場中，等到風向轉為偏北風時峰值逐漸減弱，第二次過程地面風力更大，高空形勢更有利于遠距離輸送，高值區(qū)范圍強度明顯強于第一次，重度污染層厚度達到900～1 500 m，且持續(xù)時間較長。

(2)江蘇省13個市的PM2.5模擬值與觀測值的相關系數在0.51～0.75，相關性均通過0.01顯著性檢驗，標準化平均偏差(NMB)在-0.40～0.06，模擬值總體略偏低，均方根差(RMSE)在39～50(徐州90)。模式能夠較好地反映江蘇省PM2.5的時空分布。

(3)第一次過程中宿遷、揚州、無錫的PM2.5外來輸送貢獻分別為 40.2 μg/m3、20.9 μg/m3、11.1 μg/m3，第二次過程中分別為 105.9 μg/m3、83.1 μg/m3、64.8 μg/m3。第二次過程外來輸送明顯高于第一次，外來輸送自北向南遞減，蘇北受污染輸送影響更大。

(4)第一次過程中江蘇省內排放源貢獻率在23%～79%，蘇北外來輸送貢獻率超過50%，蘇南仍以本地排放源污染為主。第二次過程中江蘇省內貢獻率在5%～32%，山東省和京津冀地區(qū)排放源是主要污染來源，二者貢獻之和在44%～70%。外來輸送山東省貢獻率高于京津冀地區(qū)，而其余周邊省份的貢獻率明顯低于兩者，安徽省對江蘇西部城市的貢獻率較高。