一次大氣污染物的大氣環(huán)流輸送規(guī)律分析與模擬

傅云燕 王金強 蔣翠花 袁志勇 王文清

摘要：以２００５年１月６日遼寧中部城市群出現(xiàn)的區(qū)域污染事件為例，利用ＮＣＥＰ／ＮＣＡＲ全球再分析資料對２００５年冬季影響我國北部的大尺度環(huán)流系統(tǒng)進行了客觀分析，同時采用中尺度數(shù)值模式ＭＭ５對不同尺度大氣邊界層進行數(shù)值模擬研究。結(jié)果表明，此次污染過程中地形影響形成的氣流輸送匯是造成區(qū)域污染的關(guān)鍵原因，而較強的逆溫、逆濕等典型的穩(wěn)定邊界層的特征則加重了局地污染。

關(guān)鍵詞：大氣環(huán)流；冬季；污染物輸送；數(shù)值模擬

中圖分類號：Ｘ５１；Ｐ４３４＋．５；Ｐ４３７文獻標識碼：Ａ 文章編號：0439－８114（２０12）06－1119-04

Ｏｂｊｅｃｔｉｖｅ Ａｎａｌｙｓｉｓ ａｎｄ Ｓｉｍｕｌａｔｉｏｎ ｏｆ Ａ Ｃａｓｅ ｏｆ Ａｔｍｏｓｐｈｅｒｉｃ Ｐｏｌｌｕｔａｎｔｓ Ｔｒａｎｓｐｏｒｔａｔｉｏｎ

ＦＵ Ｙｕｎ－ｙａｎ１，２，ＷＡＮＧ Ｊｉｎ－ｑｉａｎｇ２，ＪＩＡＮＧ Ｃｕｉ－ｈｕａ２，ＹＵＡＮ Ｚｈｉ－ｙｏｎｇ２，ＷＡＮＧ Ｗｅｎ－ｑｉｎｇ２

（１． Ｄｅｐａｒｔｍｅｎｔ ｏｆ Ａｔｍｏｓｐｈｅｒｉｃ Ｓｃｉｅｎｃｅｓ，Ｎａｎｊｉｎｇ Ｕｎｉｖｅｒｓｉｔｙ，Ｎａｎｊｉｎｇ ２１００９３，Ｃｈｉｎａ；

２． Ｓｕｑｉａｎ Ｍｅｔｅｏｒｏｌｏｇｉｃａｌ Ｂｕｒｅａｕ，Ｓｕｑｉａｎ ２２３８００，Ｊｉａｎｇｓｕ，Ｃｈｉｎａ）

Ａｂｓｔｒａｃｔ： Ｔｏ ｅｘｐｌｏｒｅ ｔｈｅ ｍｅｃｈａｎｉｓｍ ｏｆ ａｔｍｏｓｐｈｅｒｉｃ ｐｏｌｌｕｔａｎｔｓ ｉｎ Ｌｉａｏｎｉｎｇ， Ｃｈｉｎａ， ｉｎ Ｊａｎｕａｒｙ ６， ２００５， ｔｈｅ ａｔｍｏｓｐｈｅｒｉｃ ｃｉｒｃｕｌａｔｉｏｎ ｉｎ ｔｈｅ ｌａｒｇｅ－ｓｃａｌｅ of North China in the winter of 2005 ｗａｓ ａｎａｌｙzed ｂｙ ｕｓｉｎｇ ＮＣＥＰ／ＮＣＡＲ Ｒｅａｎａｌｙｓｉｓ Ｄａｔａ． Ｍｅａｎｗｈｉｌｅ ｔｈｅ ｎｕｍｅｒｉｃａｌ ｓｉｍｕｌａｔｉｏｎ ｒｅｓｅａｒｃｈ on ｔｈｅ dｉｆｆｅｒｅｎｔ－sｃａｌｅ ｏｆ ａｔｍｏｓｐｈｅｒｉｃ ｂｏｕｎｄａｒｙ ｌａｙｅｒ ａｃｃｏｒｄｉｎｇ ｔｏ Ｍｅｓｏｓｃａｌｅ Ｍｏｄｅｌ ＭＭ５ ｗａｓ ｃａｒｒｉｅｄ ｏｕｔ． Ｔｈｅ ｒｅｓｕｌｔ ｓｈｏｗed ｔｈａｔ ｔｈｅ ｍａｉｎ ｒｅａｓｏｎ ｏｆ ｔｈｉｓ ｒｅｇｉｏｎａｌ ｐｏｌｌｕｔｉｏｎ ｗａｓ ｔｈｅ ａｉｒｓｔｒｅａｍ ｃｏｎｖｅｒｇｅｎｃｅ ａｆｆｅｃｔｅｄ ｂｙ ｔｅｒｒａｉｎ， ａｎｄ ｔｈｅ ｔｙｐｉｃａｌ ｃｈａｒａｃｔｅｒｉｓｔｉｃｓ ｏｆ ｔｈｅ ｓｔａｂｉｌｉｚｅｄ ｂｏｕｎｄａｒｙ ｌａｙｅｒ of ａ ｓｔｒｏｎｇ tｅｍｐｅｒａｔｕｒｅ iｎｖｅｒｓｉｏｎ ａｎｄ iｎｖｅｒｓｅ hｕｍｉｄｉｔｙ ａｇｇｒａｖａｔｅｄ ｔｈｅ ｌｏｃａｌ ｐｏｌｌｕｔｉｏｎ．

Ｋｅｙ ｗｏｒｄｓ： ａｔｍｏｓｐｈｅｒｉｃ ｃｉｒｃｕｌａｔｉｏｎ； ｗｉｎｔｅｒ； ｐｏｌｌｕｔａｎｔ ｔｒａｎｓｐｏｒｔａｔｉｏｎ； ｎｕｍｅｒｉｃａｌ ｓｉｍｕｌａｔｉｏｎ

大氣環(huán)流是不同尺度天氣系統(tǒng)發(fā)生、發(fā)展和演變的背景條件［１］。位于對流層最底部的大氣邊界層受大尺度、中小尺度天氣系統(tǒng)以及局地系統(tǒng)的直接影響，生產(chǎn)活動產(chǎn)生的污染物在大氣邊界層中排放和輸送，必然受到不同尺度大氣環(huán)流特征的影響。近年來，大氣污染物輸送過程受到國內(nèi)外科學工作者的廣泛關(guān)注。任陣海等［２］利用觀測數(shù)據(jù)研究了北京的一次重污染過程中不同尺度大氣系統(tǒng)對該過程邊界層的影響，分析了邊界層內(nèi)粒子濃度輸送的動態(tài)變化，提出北京地區(qū)重污染的形成原因。江燕如等［３］用６年時間收集到的氣象資料進行多層次的大氣污染物輸送場的研究，初步得到我國大氣污染物輸送場的時空分布特征。蔡旭暉等［４］采用拉格朗日隨機粒子擴散模式模擬分析了珠江三角洲東部地區(qū)低層大氣流動和氣載污染物的中尺度擴散情況。吳澗等［５］運用三維非靜力Ｅ－ε閉合模式模擬了山谷地形下的流場，發(fā)現(xiàn)在近地面源條件下，穩(wěn)定層間如果出現(xiàn)小風會造成一定范圍的嚴重污染。Ｔｅｒｍｏｎｉａ等［６］提出一個表征水平和垂直輸送距離的指數(shù)來預報不利于污染物擴散時的氣象條件，研究表明在小風和大氣輸送很弱的穩(wěn)定條件下，該指數(shù)可以準確預報由持續(xù)熱源所造成的極端污染事件。Ｅｇａｎ等［７］利用一個水平對流擴散的空氣污染物輸送模式，發(fā)現(xiàn)垂直風與水平風對污染物擴散的共同作用要比單純的垂直擴散大得多。

利用全球再分析資料對２００５年冬季影響我國北部的大尺度環(huán)流系統(tǒng)進行客觀分析，再利用中尺度數(shù)值模式ＭＭ５對發(fā)生在遼寧的一次污染過程進行大氣邊界層不同尺度的數(shù)值模擬研究，克服了再分析資料在時間、空間上精度較差的缺點，以期進一步認識遼寧地區(qū)污染物的大氣環(huán)流驅(qū)動機制。

１資料及計算方案

２００５年１月６日，遼寧出現(xiàn)了一次比較嚴重的的污染過程，全省各城市污染物濃度都比較高，其中沈陽的ＳＯ２、ＮＯ２和ＰＭ１０日均濃度分別為０．２４、０．１０和０．３７ ｍｇ／ｍ３，ＳＯ２和ＰＭ１０指數(shù)分別為１１３和２３４，空氣質(zhì)量達到輕度污染級別。當天遼寧地區(qū)受蒙古高壓控制，天氣晴朗，近地面氣壓梯度小，大氣層結(jié)穩(wěn)定。一股來自貝加爾湖的氣流在內(nèi)蒙古地區(qū)開始輻合并繼續(xù)向東南方向傳播，在遼寧邊界轉(zhuǎn)向東北方向，與同樣來自西北方向的氣流交匯形成一條較強的輻合帶，覆蓋遼寧大部分地區(qū)。為了解冬季污染的形成原因，首先利用每日１次、２．５°×２．５°的ＮＣＥＰ／ＮＣＡＲ全球再分析資料作為大尺度實況場，然后采用中尺度數(shù)值模式ＭＭ５反映遼寧區(qū)域氣象場分布，此次采用的是其最新的Ｖ ３．７；采用每日４次、１°×１°的ＮＣＥＰ／ＮＣＡＲ全球再分析資料作為模式初始場。區(qū)域中心點為４３°Ｎ、１２０°Ｅ，設(shè)置３重模擬區(qū)域，格距從外至內(nèi)分別為２７、９、３ ｋｍ（圖１）。各重模擬區(qū)域選取一致的參數(shù)化方案：Ｇｒｅｌｌ積云參數(shù)化方案、Ｍｉｘｅｄ－Ｐｈａｓｅ顯式降水方案、高分辨率的行星邊界層方案。以２００５年１月作為冬季試驗月進行大尺度環(huán)流的客觀分析，選擇發(fā)生在遼寧地區(qū)的典型污染個例進行數(shù)值模擬研究。

２我國北部冬季環(huán)流概況

由再分析資料得到２００５年１月影響我國北部的大氣環(huán)流形勢（圖２）。整個中國大陸冬季都在西風環(huán)流控制之下，我國上空基本氣流是西北風；地面上，蒙古的冷高壓強度達全年最高值，中心位于貝加爾湖以西３００ ｋｍ附近，高達１ ０３８ hＰａ。冷高壓的范圍覆蓋了整個東亞地區(qū)，相當穩(wěn)定，觀察逐日場發(fā)現(xiàn)，即使它在短時間內(nèi)受到破壞，也會迎來新一次的強冷高壓入侵。因此，這個季節(jié)里冷高壓的氣流（冬季風）也十分穩(wěn)定，我國北部盛行西北－北氣流。另外，在東亞大槽東邊的阿留申群島有一個強大的低壓，其中心位于１７０°—１８０°Ｅ，４５°—５５°Ｎ，低壓值達９９５ hＰａ。由于能夠得到能量和渦度不斷地補充，這個大低壓在整個冬季基本維持穩(wěn)定不變，它與蒙古冷高壓共同影響著我國北部乃至整個亞洲的天氣形勢。

３模式檢驗

選取最外重網(wǎng)格的計算結(jié)果，與再分析資料的８５０ hＰａ 1月份平均高度場的比較可以發(fā)現(xiàn)，由于地形的影響，模擬的等值線不平滑，且出現(xiàn)了若干個小的高、低值中心?？傮w來說，模擬得到的高度等值線的數(shù)值和分布與實況非常相近（圖３），模式模擬性能較好，具有一定的可信度。

４污染實例的模擬分析

下面結(jié)合污染實例，具體分析２００５年１月６日影響遼寧地區(qū)的環(huán)流形勢和天氣系統(tǒng)。

４．１大尺度環(huán)流特征

圖４ａ為外網(wǎng)格區(qū)域2005年１月６日０５∶００的地面流場，同時刻對應的海平面氣壓場分布（圖４ｂ）為西高東低，遼寧處于蒙古高壓前緣大范圍的均壓帶控制區(qū)域，等壓線稀疏，風速較??；內(nèi)蒙古、河北境內(nèi)等壓線密集，氣壓較高。受高層西風帶的影響，形成了輸入遼寧的西北氣流，而氣流在遼寧西北角轉(zhuǎn)向東北方向則是地形影響造成的。

４．２天氣尺度大氣輸送特征

遼寧省地勢從東北向西南傾斜，東、西兩側(cè)是丘陵，中部是長方形遼河平原，呈馬蹄形向渤海傾斜。遼寧西北角地形平均高度超過６００ ｍ，與其相鄰的河北東部及內(nèi)蒙古東北部地勢高度為８００～１ ４００ ｍ，是典型的山地、丘陵地形（圖５ａ）。在這個過程中，由于夜間河北境內(nèi)的山坡輻射冷卻，其降溫速度比同高度的空氣要快，冷氣流沿地勢高的位置流向地勢較低的遼西山區(qū)，進而繼續(xù)流向東部地勢更低的平原地區(qū)，這一現(xiàn)象正是由于熱力原因引起的山風環(huán)流。來自西北方向的另一股氣流從遼寧正北方輸入，與前者匯聚成輻合帶，即輸送匯（圖４ａ）。輸送匯的形成常常伴隨有較低的大氣壓［２］，這在圖５ｂ中有所體現(xiàn)，圖中虛線所包括的范圍是１ ０２１ ｈＰａ的低壓帶，也正對應著此次過程的輸送匯。

４．３局地性風場特征及其對邊界層的影響

由最內(nèi)重網(wǎng)格的計算結(jié)果得到２００５－０１－０６T０５：００輸送匯中包含遼寧中部城市群在內(nèi)的一部分區(qū)域的風場形勢（圖６）。由圖６可知，風矢量在沈陽發(fā)生氣旋性切變，輸送到撫順并繼續(xù)向東北方向傳播，沈陽、撫順風速不超過３ ｍ／ｓ。由于風速較小，污染物不易疏散而產(chǎn)生堆積造成局地污染。此時大氣層結(jié)穩(wěn)定，邊界層中各種特征量都有典型的分布。沿４１．５°Ｎ作覆蓋遼寧中部城市群的垂直剖面，輸出同一時刻的物理量（圖７）。邊界層內(nèi)存在較強的逆溫和逆濕梯度，穩(wěn)定邊界層厚度可達到１ ０００ ｍ以上，２００ ｍ高度以內(nèi)幾個主要城市逆溫強度為１．２５ ℃／（１００ ｍ），逆濕強度達５×１０－３ｇ／（１０ｍ·ｋｇ）。而水平風速隨高度分布最明顯的特征是邊界層頂附近出現(xiàn)風速最大值（１５ ｍ／ｓ），即形成了低空急流。

５結(jié)論

結(jié)合客觀場分析資料研究發(fā)現(xiàn)，冬季蒙古冷高壓和阿留申低壓共同影響著我國北部乃至整個亞洲的天氣形勢，強度很大而且相當穩(wěn)定，形成穩(wěn)定的冬季風，我國北部盛行西北－北氣流。由于再分析資料時間、空間上的精度都比較差，因此利用中尺度數(shù)值模式ＭＭ５對不同尺度大氣邊界層進行數(shù)值模擬。結(jié)果顯示，ＭＭ５模式能夠較好地再現(xiàn)冬季的氣象場形勢，模擬結(jié)果具有一定的說服力。實例分析發(fā)現(xiàn)，地形影響造成的氣流輸送匯是造成遼寧重污染過程的主要原因。遼寧中部城市群出現(xiàn)了較強的逆溫、逆濕等典型的穩(wěn)定邊界層的特征，邊界層頂附近出現(xiàn)了低空急流。

參考文獻：

［１］ 朱乾根，林錦瑞，壽紹文，等． 天氣學原理和方法［Ｍ］． 第二版.北京：氣象出版社，１９９２.２０８－２８７．

［２］ 任陣海，萬本太，虞統(tǒng)，等． 不同尺度大氣系統(tǒng)對污染物邊界層的影響及其水平流場輸送［Ｊ］． 環(huán)境科學研究，２００４，１７（１）：７－１３．

［３］ 江燕如，梁漢明，高慶先，等． 我國大氣污染物在排污層的輸送特征［Ｊ］． 中國環(huán)境科學，１９９８，１８（６）：４８６－４８８．

［４］ 蔡旭暉，張禎璽，康凌，等． 珠江三角洲東部低層大氣流動與中尺度擴散特性［Ｊ］． 熱帶氣象學報，２００５，２１（５）：４９７－５０５．

［５］ 吳澗，蔣維楣，王衛(wèi)國． 山谷地形流場和擴散的數(shù)值研究［Ｊ］． 高原氣象，２００１，２０（２）：１４０－１４７．

［６］ ＴＥＲＭＯＮＩＡ Ｐ， ＱＵＩＮＥＴ Ａ． Ａ ｎｅｗ ｔｒａｎｓｐｏｒｔ ｉｎｄｅｘ ｆｏｒ ｐｒｅｄｉｃｔｉｎｇ ｅｐｉｓｏｄｅｓ ｏｆ ｅｘｔｒｅｍｅ ａｉｒ ｐｏｌｌｕｔｉｏｎ［Ｊ］． Ｊｏｕｒｎａｌ ｏｆ Ａｐｐｌｉｅｄ Ｍｅｔｅｏｒｏｌｏｇｙ，２００４，４３（４）：６３１－６４０．

［７］ ＥＧＡＮ Ｂ Ａ， ＭＡＨＯＮＥＹ Ｊ Ｒ． Ｎｕｍｅｒｉｃａｌ ｍｏｄｅｌｉｎｇ ｏｆ ａｄｖｅｃｔｉｏｎ ａｎｄ ｄｉｆｆｕｓｉｏｎ ｏｆ ｕｒｂａｎ ａｒｅａ ｓｏｕｒｃｅ ｐｏｌｌｕｔａｎｔｓ［Ｊ］． Ｊｏｕｒｎａｌ ｏｆ Ａｐｐｌｉｅｄ Ｍｅｔｅｏｒｏｌｏｇｙ，１９７２，１１（2）：312－322．