魯中地區霾期間氣溶膠的垂直分布及其與氣象條件的關系*

張 琴 姚秀萍

1 山東省淄博市氣象局，淄博 255048 2 中國氣象局氣象干部培訓學院，北京 100081

提 要： 利用CALIPSO氣溶膠垂直廓線產品和NCEP提供的FNL 1°×1°再分析數據，以清潔日為對比，統計分析了2007—2017年霾和重度霾期間魯中地區氣溶膠的垂直分布規律，并探討了兩者之間的關系。結果表明：霾日、重度霾日的氣溶膠粒子主要聚集在對流層低層2.7 km、1.5 km以下。以0.9、1.66和1.34 km為臨界高度，清潔日、霾日、重度霾日的消光系數(EC)在臨界高度以上呈指數遞減，以下呈對數遞減。霾日和重度霾日年平均體積退偏比(PDR)和色比(CR)范圍分別為0.1～0.3和0.5～0.9。低層2 km以下，重度霾日以相對規則的細粒子為主；2～4 km，不規則度和尺寸都急劇增大。PDR在7 km 及以下隨高度逐步增加，7 km以上逐步減小。與清潔日不同，霾日CR基本隨高度逐漸增大,重度霾日則在6 km以下隨高度逐漸增大。污染沙塵型氣溶膠對重度霾日和霾日近地面的EC貢獻最大，分別為0.58 km-1和0.36 km-1,其次是大陸污染型，分別為0.27 km-1和0.20 km-1。重度霾日在1.5 km及以下維持了較小的平均風速和較大的相對濕度，有利于EC的增大。霾日在較高的高度上維持了相對較小的平均風速和相對較大的相對濕度，使得氣溶膠高濃度層達2.7 km。霾時近地面污染物來源為蒙古、內蒙古等地的長距離傳輸和周邊地區的傳輸。重度霾時蒙古方向和河南方向的污染物傳輸占比更大，分別達25.26%和31.58%。重度霾對應的多層逆溫比例接近50%，貼地+低懸型和貼地+高懸型逆溫在其中占比最大。由于這兩種逆溫型在逆溫層以下對應更小的出流量，大氣更穩定，使得重度霾時氣溶膠粒子聚集在低層。

引 言

霾是大量粒徑為幾微米以下的大氣氣溶膠粒子使水平能見度小于10.0 km、空氣普遍渾濁的天氣現象(國家市場監督管理總局和國家標準化管理委員會，2018)，其本質是細粒子污染(吳兌，2013)。霾發生時，大量的有害氣體和顆粒物使空氣質量下降，嚴重危害人體健康。《大氣科學辭典》中提及“霾可看做氣溶膠”(《大氣科學辭典》編委會，1994)。近些年來我國大面積的霾天氣就與細顆粒氣溶膠有關，氣溶膠顆粒是霾的主要成分(穆穆和張人禾，2014)。

氣溶膠垂直分布是氣溶膠輻射強迫研究(Haywood and Ramaswamy, 1998)和氣溶膠光學厚度衛星反演(Rozwadowska，2007)的重要影響因素。目前，對于氣溶膠的垂直分布的研究多集中于京津冀、淮河流域、珠江三角洲地區(Gui et al,2016;范偉等,2020;Cao et al,2018;劉瓊等,2012)，也有對華北地區和中國大范圍地區氣溶膠垂直分布的研究(Cao et al,2018;Liu et al,2017;許瀟鋒等,2018)。但對于霾時，尤其是重度霾時氣溶膠垂直分布的統計分析不是很多。近年來，關于霾時氣溶膠垂直分布與不同的氣象條件的相關性也有不少研究。如不同地區溫度垂直廓線的日變化所導致的邊界層高度發生變化，從而影響近地層污染物濃度；秋、冬季重污染過程頻發與邊界層的季節變化的關系(Miao et al,2016;2018；Xu et al,2019)；在京津冀特殊地形的影響下，局地熱力環流的晝夜變化以及地形通過影響邊界層來影響污染物的輸送(Miao et al,2020),消光系數演變與PM2.5濃度值變化以及天空狀況對氣溶膠垂直分布影響顯著(徐棟夫等，2020)。由于氣溶膠濃度不僅是人為排放所導致的，也與氣象條件息息相關，因此不同地域間的氣溶膠垂直分布有很大差異。

魯中地區位于山東中部和南部，是一個孤立的山區，面積為6.5萬km2。魯中地區西面250 km處是平均海拔高度超過1 000 m的太行山脈，東部瀕臨渤海，猶如一座孤島。魯中地區有著獨特的氣候條件，對該地區氣溶膠的垂直分布特性及其與氣象條件的關系進行統計分析，有助于對魯中地區霾污染事件更深入地理解，為魯中地區大氣氣溶膠垂直分布特征及其成因研究提供參考。

1 數據與方法

1.1 霾的判別

選取魯中地區5個國家級地面氣象觀測站(濟南、淄博、萊蕪、泰安、濰坊；圖1)的氣象要素觀測資料對霾日進行判別。根據吳兌(2008)文中的方法選取霾日，只要有3個及以上的觀測站滿足下列條件：排除降水、吹雪、揚沙、沙塵暴和浮塵等影響能見度的天氣現象，相對濕度<90%，能見度(VS)<10 km時，則認為有霾發生。當VS<2.0 km時為重度霾，2.0 km≤VS<3.0 km時為中度霾，3.0 km≤VS<5.0 km 時為輕度霾。另外，本文將相對濕度<90%，且VS≥15 km時定義為清潔日。經統計，魯中地區2007—2017年間，霾日為1 655 d，重度霾日為95 d，中度霾日為283 d，輕度霾日為552 d，清潔日為235 d。

圖1 魯中地區地形圖(紅色圓點為氣象觀測站點的位置)Fig.1 Topographic map of central part of Shandong Province (Luzhong Area)(Red dots are the locations of meteorological observation stations)

1.2 數據來源

本文使用的衛星數據來自云-氣溶膠激光雷達和紅外探測者觀測衛星(CALIPSO)。CALIPSO首次實現了云和氣溶膠的全球觀測，可提供云和氣溶膠的物理特性及其時空分布，有利于更好地研究云和氣溶膠的形成、發展過程及其相互作用，通過提供氣溶膠的垂直結構區分其種類(張鵬等，2018)。

本文采用了2007—2017年霾時CALIIPSO衛星資料中的L2產品途經魯中地區的氣溶膠廓線數據以及垂直特征層VFM(vertical feature mask)產品。 CALIPSO衛星于13時和02時(世界時，下同)左右經過魯中地區上空。CALIPSO的廓線在魯中地區(36°～37°N、116°～120°E)篩選的清潔日、霾日和重度霾日分別為49、367、23 d。本文選取霾期間的夜間數據，有利于降低白天噪音的影響造成的偏差；并對廓線數據進行了嚴格的質量控制。利用VFM產品中的氣溶膠類型數據，對L2中數據進行篩選，剔除氣溶膠以外的數據。選取的L2廓線數據包括4個參量：532 nm消光系數(extinction coefficient，EC)、532 nm平行后向散射系數(β532，∥)、532 nm垂直后向散射系數(β532，⊥)和1 064 nm后向散射系數(β1 064)。CALIPSO VFM產品將氣溶膠分為潔凈海洋型、沙塵型、污染大陸型、潔凈大陸型、污染沙塵型、煙塵型、海洋沙塵型7種類型。

式中：z為高度；532 nm總后向散射系數(β532)為垂直(β532，⊥)與平行(β532，∥)后向散射系數之和；RPDR為體積退偏比；色比(RCR)為1 064 nm總后向散射系數(β1 064)β532的比值。

NCEP FNL再分析數據為2007—2017年1°×1°數據，包括1000～500 hPa濕度、風向、風速。為了與CALIPSO衛星數據相匹配，逆溫計算過程中溫度數據選用L2產品中的溫度反演量。

1.3 分析方法

為追蹤污染來源，利用NOAA開發的HYSPLIT4.8軌跡模式定性分析了霾天污染物的傳輸途徑。HYSPLIT4.8是一種歐拉和拉格朗日混合型的計算模式，包含多種物理過程，可針對不同類型排放源進行較完整的輸送擴散和沉降過程模擬，廣泛應用于大氣污染物輸送研究(Draxler and Hess,1997)。數據來源于NCEP提供的全球數據同化系統GDAS。設置計算方式為后向軌跡，設置高度為500 m，計算時長為24 h，聚類計算選取霾日和重度霾日18時的軌跡。

2 結果與討論

2.1 氣溶膠垂直分布特征

EC體現了氣溶膠吸收和散射截面隨高度的分布，反映了顆粒物的含量(Wang et al，2014)。圖2為2007—2017年魯中地區清潔日、霾日、重度霾日532 nm的EC年平均廓線。可以看出，重度霾日、霾日、清潔日的EC隨高度均呈遞減趨勢，在地面附近，分別是1.2、0.75、0.5 km-1。2 km以下遞減迅速。重度霾日EC在0.7 km以下最強，隨后迅速遞減，在1.5 km以上遞減至三種天氣中的最小值。在高層，霾日EC在2.7 km以下大于清潔日，說明霾時氣溶膠大值層頂高約為2.7 km，氣溶膠粒子向高層擴散能力強于重度霾日。3 km以上，EC很小，廓線也基本重合，大氣較為清潔。

圖2 2007—2017年魯中地區清潔日、霾日、重度霾日532 nm的EC年平均廓線Fig.2 Annual average profile of 532 nm EC on clear days, haze days and heavy haze days in central part of Shandong Province during 2007-2017

有研究表明，隨著高度增加，氣溶膠濃度一般呈指數關系下降(韓永等，2009)。標高(Ha)是假定大氣氣溶膠濃度隨高度指數衰減到地面濃度的1/e時氣溶膠層的高度。通過計算，清潔日、霾日、重度霾日氣溶膠的Ha分別為2.00、1.58、0.84 km。為了更清楚地表示氣溶膠的垂直變化規律，以z/Ha為自變量，對氣溶膠濃度σ(z)進行曲線擬合(圖3)。發現指數方程和對數方程的擬合優度R2最大。在每層高度上采用絕對估計偏差(DE)的方式來判定兩種擬合效果的好壞，定義DE=|VT-VM|,VT為觀測值，VM為擬合值。當DE最小時，其擬合效果最好。指數曲線和對數曲線上的DE相等時，得到z/Ha的臨界值。可以發現，三種不同天氣狀態下，z/Ha分別在大于臨界值0.45、1.05、1.60，即z分別為0.9、1.66、1.34 km時，指數方程擬合效果較好。低于此臨界值，對數方程擬合效果較好。

圖3 2007—2017年魯中地區清潔日(a)、霾日(b)、重度霾日(c)532 nm的EC垂直廓線擬合Fig.3 Fitting of 532 nm EC vertical profile on clear days (a), haze days (b) and heavy haze days (c) in central part of Shandong Province during 2007-2017

體積退偏比(particulate deploarization ratio,PDR)可以反映氣溶膠粒子的不規則程度，色比(color ratio,CR)反映氣溶膠粒子的大小。圖4是2007—2017年清潔日、霾日、重度霾日期間PDR和CR的垂直分布。可以看出，三條PDR廓線范圍在0.1～0.3浮動。2 km以下，霾和重度霾日PDR值小于清潔日。2～4 km，重度霾日PDR急劇增大，霾日最小。4 km以上，清潔日PDR最大，霾日PDR較小。三條CR廓線范圍在0.5～0.9浮動。2 km以下，重度霾日CR最小，其次是霾日，清潔日最大。2～3 km，重度霾日CR隨高度急劇增大至0.87。3 km以上，重度霾期間CR隨高度變動較大。2～5 km，清潔日CR大于霾日。5 km以上，清潔日與霾日CR基本一致。由于魯中地區高層長期受到外來沙塵粒子的影響，而地面又多以細粒子為主，造成了PDR與CR隨高度增大的特征。低層2 km以下，重度霾時以大量的相對規則的細粒子為主，霾時以相對規則的粗粒子為主，清潔日則以不規則的粗粒子為主。2～4 km，重度霾不規則度和尺寸都急劇增大，可能與高空的不規則沙塵粒子有關。如2016年11月13日出現的重度霾天氣(圖5)，CALIPSO衛星監測到魯中地區2～4 km分布著沙塵型氣溶膠。利用后向軌跡模式模擬24 h的氣溶膠輸送路徑后發現，12日，來自甘肅騰格里沙漠附近4～5 km高度的沙塵型氣溶膠，經過黃土高原向東輸送至魯中地區的3 km高度附近，造成13日魯中地區2～4 km氣溶膠粒子尺寸和不規則度的急劇增大。

圖4 2007—2017年魯中地區清潔日、霾日、重度霾日532 nm PDR(a)和CR(b)年平均廓線Fig.4 Annual average profile of 532 nm PDR (a) and CR (b) on clear days, haze days and heavy haze days in central part of Shandong Province during 2007-2017

圖5 2016年11月12日19：13(a)和13日18：34(b)氣溶膠類型垂直剖面，13日18：34 (c)魯中地區氣團24 h HYSPLIT后向軌跡(黑線)、CALIPSO衛星軌跡(藍線、紅線)Fig.5 Vertical profile of aerosol types at 19：13 UTC 12 (a) and 18：34 UTC 13 (b)，24 h HYSPLIT backward trajectory (black line), CALIPSO satellite track (blue and red solid lines) of air mass in central part of Shandong Province at 18:34 UTC 13 (c) November 2016

為了清楚地說明PDR和CR在垂直方向上的變化趨勢，將PDR和CR與高度進行擬合，最優擬合結果如圖6所示。三種天氣下，PDR和CR的最優擬合方程均為三次方擬合方程。PDR擬合優度R2分別為0.868、0.954、0.310。CR擬合優度R2分別為0.571、0.922、0.420(表略)。霾日的擬合優度最好，重度霾日最差，可能與重度霾日個例偏少有關。從圖6中可以看出，地面上三種天氣下PDR的值分別為0.131、0.102、0.105。三條PDR廓線在7 km 以下隨高度逐步增加；7 km以上，逐步減小。三種天氣下，PDR最大值分別為0.25、0.22、0.23。地面上CR值在三種天氣下分別為0.602、0.602、0.520。清潔日時，CR在3 km以下隨高度增大，3～8 km保持不變，8 km以上隨高度增大。與清潔日不同，霾時CR基本隨高度逐漸增大。重度霾則在6 km以下隨高度逐漸增大，6 km以上逐漸減小。

圖6 2007—2017年魯中地區清潔日(a,b)、霾日(c,d)、重度霾日(e,f)的PDR(a,c,e)和CR(b,d,f)的垂直廓線擬合Fig.6 Vertical profile fitting of PDR (a, c, e) and CR (b, d, f) on clear days (a, b), haze days (c, d) and heavy haze days (e, f) in central part of Shandong Province during 2007-2017

2.2 氣溶膠類型的垂直分布特征

統計2007—2017年魯中地區清潔日、霾日、重度霾日氣溶膠類型出現的次數，計算相應占比(圖7)。清潔日、霾日、重度霾日氣溶膠類型總體上以沙塵型、污染沙塵型為主，分別占比30%和50%左右。清潔日沙塵型氣溶膠明顯偏多，且重度霾日的沙塵型氣溶膠多于霾日。在霾日和重度霾日，由于本地污染物的生成和外來污染物的傳輸，污染大陸型氣溶膠多于清潔日。霾時煙塵型氣溶膠最多，占比在10%左右。由于魯中地區距離東部沿海較遠，清潔日、霾日、重度霾日潔凈海洋型與潔凈大陸型氣溶膠最少。

圖7 2007—2017年魯中地區清潔日、霾日、重度霾日氣溶膠的成分占比(1：潔凈海洋型，2：沙塵型，3：污染大陸型，4：潔凈大陸型，5：污染沙塵型，6：煙塵型，7：海洋沙塵型)Fig.7 Aerosol composition proportion on clear days, haze days and heavy haze days in central part of Shandong Province during 2007-2017(1： clear marine, 2： dusty, 3： polluted continental, 4： clear continental, 5： polluted dusty, 6： elevated smoke， 7： dusty marine)

圖8是2007—2017年魯中地區清潔日、霾日、重度霾日不同高度上的氣溶膠成分比例分布。由圖可知，清潔日、霾日、重度霾日均以污染沙塵型和沙塵型為主。隨著高度的增加，污染沙塵比例逐漸減小，沙塵逐漸增多。在0～4 km的低層大氣，主要以污染沙塵型為主。4 km以上，則主要以沙塵型為主。在高層大氣中，由于沙塵粒子更大且不規則，氣溶膠粒子的PDR和CR更大，這與圖4的分析一致。不同的是，2～4 km，霾日期間污染沙塵氣溶膠比例小于清潔日和重度霾日，4 km以上相反。整層大氣中，霾日期間沙塵氣溶膠比例均小于清潔日和重度霾日。污染大陸型氣溶膠主要分布在低層大氣中。0～2 km，霾日、重度霾日、清潔日污染大陸型氣溶膠比例依次減少，分別為0.23、0.20、0.11。 煙塵型氣溶膠主要分布在2～6 km的中層大氣中。霾日、清潔日、重度霾日煙塵型氣溶膠比例依次減少，最大分別為0.20、0.10、0.05。海洋沙塵型氣溶膠主要分布在低層大氣，清潔日偏多。清潔大陸型氣溶膠比例較低，主要分布在中層大氣中。由于魯中地區位于山東省內陸，清潔海洋型氣溶膠分布很少。

圖8 2007—2017年魯中地區清潔日(a)、霾日(b)、重度霾日(c)不同高度上氣溶膠成分比例Fig.8 Aerosol composition ratio at different heights on clear days (a), haze days (b) and heavy haze days (c) in central part of Shandong Province during 2007-2017

2.3 不同氣溶膠類型對消光的貢獻

通過前面的分析，清潔日、霾日、重度霾日的EC和氣溶膠類型在不同高度上分布不同。為了給出各類型氣溶膠對EC的定量貢獻，本文計算了不同高度上各氣溶膠類型對應的EC年平均值(圖9)。可以看出，與EC的高度廓線相對應(圖2)，各類型氣溶膠對EC的貢獻逐漸減小。清潔日時，2 km以下污染沙塵型的貢獻最大，其EC年平均最大達0.28 km-1，其次是沙塵型，再次是污染大陸型。2～4 km貢獻大小依次為沙塵型、污染沙塵型、煙塵型。霾日時，2 km以下也是污染沙塵型的貢獻最大，EC年平均值為0.36 km-1，大于清潔日；污染大陸型的貢獻是0.21 km-1；沙塵型的貢獻約為0.05 km-1，小于清潔日。重度霾時，0～2 km，污染沙塵型的貢獻最大，其次是污染大陸型，分別達0.58 km-1和0.27 km-1，遠大于霾日和清潔日。由于這兩種氣溶膠對應的顆粒物在貼近地面處分布多，濃度大，且粒子較小，吸濕性強，導致重度霾日近地面顆粒物EC偏大且規則度和尺寸偏小。2～4 km，沙塵型貢獻最大，其次是污染沙塵型。不規則的沙塵粒子尺寸大，導致了PDR和CR的急劇增大。從圖9中還可以看出，在接近地面的低層大氣，污染沙塵型對清潔日、霾日、重度霾日EC的貢獻最大，其中貢獻大小依次為重度霾日、霾日、清潔日。與清潔日不同，霾日和重度霾日期間，污染大陸型對EC的貢獻僅次于污染沙塵型。沙塵型在清潔日和重度霾日中對EC的貢獻大于霾日，尤其是在地面附近區域。霾日時，煙塵型在2～4 km的中層大氣對EC有一定的貢獻，是造成霾日污染物層比重度霾日污染物層高的原因。

圖9 2007—2017年魯中地區清潔日(a)、霾日(b)、重度霾日(c)各類型氣溶膠對EC貢獻的垂直廓線Fig.9 Vertical profile of the contribution of various aerosols to EC on clear days (a), haze days (b) and heavy haze days (c) in central part of Shandong Province during 2007-2017

2.4 氣象因素對霾時氣溶膠垂直分布的影響

2.4.1 風速

圖10給出了魯中地區清潔日、霾日、重度霾日近地面風速分布和1 000～500 hPa垂直平均風速廓線。從圖中可知，與清潔日相比，霾日和重度霾日的地面風速主要集中在0～4 m·s-1，在0～2 m·s-1的占比最大。其中，重度霾日地面風速在0～2 m·s-1區間的占比大于霾日。清潔日在地面風速2～6 m·s-1的占比最大。從垂直廓線可以看出，霾日在垂直方向的平均風速始終小于清潔日。近地面925 hPa(約0.8 km)以下重度霾日平均風速最小，穩定的風場使污染物難以擴散，從而有利于污染物的積聚。925 hPa以上，霾日的平均風速維持在6 m·s-1左右，遠小于重度霾日和清潔日，使得霾日污染物層頂高達2.7 km，高于重度霾日。850 hPa(1.5 km)以上，重度霾日與清潔日的平均風速重合，使得重度霾在1.5 km以上污染物得到有效擴散，污染物層頂維持在1.5 km附近。可見，地面風速越小，越有利于霾日顆粒物的聚集。垂直方向上風速的大小與氣溶膠層的垂直分布有很好的對應關系。

圖10 2007—2017年魯中地區清潔日、霾日、重度霾日近地面風速分布(a)和垂直平均風速廓線(b)Fig.10 Distribution of near-surface wind speed (a) and vertical mean wind speed profile (b) on clear days, haze days and heavy haze days in central part of Shandong Province during 2007-2017

2.4.2 濕度

許多氣溶膠成分具有吸濕性，大小與相對濕度有關(韓永等,2009；Clarke et al,2004;Covert et al,1972)。有研究表明，當相對濕度超過70%～80%時，EC可增加100%以上(Swietlicki et al，2008)。由于清潔日、霾日、重度霾日氣溶膠的EC在500 hPa及其以上很小，本文只討論500 hPa以下的大氣相對濕度。

計算魯中地區近地面相對濕度分布發現(圖11a)，清潔日、霾日及重度霾日的平均地面相對濕度分別為31.38%、45.20%、55.13%。清潔日地面相對濕度分布較為集中，主要介于10%～40%，霾日則主要集中在20%～40%。重度霾日地面相對濕度主要集中在70%～80%與30%～40%。研究表明，霾時大氣中人為排放的硫酸鹽、硝酸鹽和碳氫化物等氣溶膠粒子有較強的吸濕性，當風速小、擴散條件不利時，氣溶膠粒子快速吸濕增長，易形成霾。可以看出，大部分重度霾日對應著地面上較大的相對濕度，但仍有少部分重度霾日對應較小的相對濕度，說明大的地面相對濕度可以導致重度霾的發生，但濕度不是重度霾發生的必要條件。

從垂直方向上來看(圖11b)，地面到5.6 km(約500 hPa)，霾日的相對濕度在40%～45%，清潔日則集中在30%～38%。重度霾時，1 km(約800 hPa)以下的濕度在50%以上，明顯大于霾日，雖然重度霾和霾日均以污染大陸型和污染沙塵型為主，但重度霾氣溶膠粒子更小。大量的細粒子在吸濕增長后，有利于重度霾的維持和加強。1 km以上，重度霾的相對濕度隨高度逐漸減小到與霾日相同。由氣溶膠的類型來看，霾日在1 km以上對應的煙塵型逐漸增多，較小的氣溶膠粒子吸濕后，使得EC大于重度霾，且粒子更加規則。對于清潔日，3 km(約700 hPa)以下，相對濕度遠小于霾日，加上整層風速偏大，沙塵型氣溶膠較多，其EC小于霾日。3 km以上，相對濕度開始增大，清潔日、霾日、重度霾日相對濕度基本相同，氣溶膠的EC也基本一致。

圖11 同圖10，但為近地面相對濕度分布(a)和垂直平均相對濕度廓線(b)Fig.11 Same as Fig.10, but for distribution of near-surface relative humidity (a) and vertical average relative humidity profile (b)

2.4.3 風向

為統計清潔日、霾日、重度霾日不同高度層的風向分布，分別選取了1 000、850、500 hPa三個氣壓層進行統計(圖12)。可以看出，清潔日，500 hPa上，魯中地區以西偏西北風為主導風向，占43.00%；850 hPa以偏西北風為主，西北風占18.30%，西北偏北風占17.00%；1 000 hPa以南風為主，占20.43%，其次為西偏西南風，占16.60%。霾時500 hPa風向以西風和西北偏西風為主，分別占31.52% 和27.45%，可見霾時高空氣流更為平直。與清潔日不同，霾時850 hPa風向以偏西風為主，西南風、西南偏西風及西風分別占14.13%、11.96%、11.96%。1 000 hPa風向以偏南風為主，南風、西南偏南風分別占比18.75%與19.00%。可見，霾時風向隨高度的增加呈順時針旋轉，易出現暖平流。重度霾時，500 hPa風向與霾時類似，也是以偏西風為主。不同的是，西風出現的比例略大，說明高空氣流更加平直，同時西南風出現的比例比霾時更大，達13.00%，這與重度霾時易出現的高空短波槽脊有關。重度霾時850 hPa風向主要為偏西風和偏北風，偏北風出現的比例明顯大于霾日，達17.39%，是重度霾日污染沙塵型氣溶膠偏多的原因。1 000 hPa風向以西南風為主，西南風的出現有利于河南、安徽等地污染物向山東傳輸，從而有利于重度霾的發生。

圖12 2007—2017年魯中地區清潔日(a)、霾日(b)、重度霾日(c)風向分布Fig.12 Distribution of wind direction distribution on clear days (a), haze days (b) and heavy haze days (c) in central part of Shandong Province during 2007-2017

2.4.4 污染物傳輸路徑

魯中地區霾和重度霾時污染物后向軌跡聚類分析見圖13。霾日，來自蒙古高原的軌跡占總量的7.07%，軌跡在高空700 hPa(3 km)向東南方向傳輸，越過內蒙古后逐漸下沉，經河北到達魯中地區；來自內蒙古東北部的軌跡占總量的6.52%，軌跡自800 hPa(1.5 km)逐漸下沉至魯中地區；自東南方向的江蘇和西南方向的河南、安徽占比分別為23.37%和12.50%；河北方向為26.90%；本地傳輸占比為23.64%。重度霾時，來自蒙古高原的遠距離高空傳輸軌跡占比為4.21%，少于霾日，傳輸距離和高度均大于霾日；來自內蒙古東部的軌跡占比合計為25.26%，垂直方向上自850～700 hPa(1.5～3 km)向南傳輸后下沉至魯中；來自河南和江蘇方向的軌跡占比分別為31.58%和13.68%；本地傳輸為25.26%。

圖13 2007—2017年魯中地區霾日(a)、重度霾日(b)后向軌跡聚類分析Fig.13 Backward trajectory clustering analysis of haze days (a) and severe haze days (b) in central part of Shandong Province during 2007-2017

2.4.5 逆溫層

大氣邊界層內逆溫抑制近地層污染物的擴散(Zhang and McMurry，1992)。根據文中2.1節的分析，污染物主要集中在2 000 m以下，因此本文主要統計底高在2 000 m以下的逆溫。用上層溫度減去下層溫度，當ΔTi>0時，該層判定為逆溫層。當兩層逆溫層之間小于100 m時，作為一個逆溫層統計。本文將逆溫分為三類：底高小于300 m的貼地逆溫，底高位于300～1 000 m的低懸逆溫，底高在1 000～2 000 m的高懸逆溫。

對霾和重度霾進行逆溫計算，同時利用中度和輕度霾來進行比較分析，以期更好地揭示逆溫在霾和重度霾中的作用。可以發現，霾時逆溫出現的比例很高，為97.80%。其中，逆溫出現的比例隨霾的加重而增大，重度霾為100%，中度霾為98.20%，輕度霾為97.20%。以出現逆溫層的霾日作為統計對象，統計不同強度霾時逆溫層數分布(表1)發現，霾時出現一層逆溫的比例達83.8%，重度、中度及輕度霾比例依次增大。霾時出現二層逆溫的比例為14.3%，與中度、輕度霾相差不大；重度霾比例達42.9%，明顯大于其他強度霾。霾時三層逆溫的比例很小，為2.0%；重度、中度霾比例分別為4.6%與5.5%，輕度霾未出現三層逆溫。可以看出，大多數霾日出現一層逆溫，其次是二層逆溫。霾的強度越大，多層逆溫出現的比例越高，重度霾對應的多層逆溫比例接近50%。

表1 2007—2017年魯中地區不同強度霾時逆溫層數分布(單位：%)Table 1 Distribution of temperature inversion layers on haze days with different intensities in central part of Shandong Province during 2007-2017 (unit: %)

經過統計，魯中地區霾時共有七種逆溫層類型(表2)。從逆溫類型比例的分布來看，霾日期間，貼地逆溫出現的比例最大，達84.96%。其次是貼地+低懸型和貼地+高懸型，占比分別為4.42%和6.19%。隨著霾強度的增大，貼地型逆溫占比越小，而其他類型占比則越大。其中，重度霾在不同強度霾中貼地型比例最小(52.63%)，貼地+低懸型、貼地+高懸型、高懸型逆溫比例最大(分別為15.79%，15.79%，10.53%)。可見，除了貼地型，貼地+低懸型和貼地+高懸型逆溫也對魯中地區重度霾的發展維持起到了很重要的作用。

表2 2007—2017年魯中地區不同強度霾時逆溫類型分布(單位：%)Table 2 Distribution of temperature inversion types on haze days with different intensities in central part of Shandong Province during 2007-2017 (unit: %)

2.4.6 不同類型逆溫對大氣擴散能力的影響

逆溫使大氣穩定地分層，導致低層擴散能力較差。風速被認為是穩定邊界層的一個關鍵參數(廖曉農等，2014)。弱風有利于逆溫的發展和維持;強風則會導致湍流加強，導致地表與高空的冷暖空氣混合，從而導致逆溫減弱(Guo et al,2015;Pietroni et al,2014)。另外，在穩定的條件下，逆溫高度以下湍流交換減弱，易造成污染物的積累(van Hooijdonk et al,2017)。因此通過計算逆溫層下空氣出流量(F)來考察不同逆溫型時氣溶膠的擴散能力。其中:

F=WS×IH

式中WS為逆溫層以下的平均全風速。

逆溫層高度(IH):

IH=(頂高-底高)/2

通過計算，貼地+低懸型和貼地+高懸型逆溫層高度、平均全風速、出流量相對較低。其次是貼地型、貼地+低懸+高懸型。低懸型和高懸型的逆溫層高度、平均全風速、出流量相對較大。可見，當逆溫為貼地+低懸型和貼地+高懸型逆溫時，逆溫層下的大氣最穩定，有利于重度霾時污染物的積聚；其次是貼地型和貼地+低懸+高懸型逆溫。三種懸浮型逆溫時，逆溫層下大氣穩定性最差。

表3 2007—2017年魯中地區霾時不同類型的逆溫層高、逆溫層下的平均全風速和大氣出流量Table 3 Height of different types of inversion layers, average total wind speed under the inversion layer and atmospheric outflow flow on haze days in central part of Shandong Province during 2007-2017

3 結 論

基于CALIPSO衛星數據和NCEP提供的FNL再分析數據，對2007—2017年魯中地區霾期間的氣溶膠垂直分布及其氣象條件進行了統計分析，結果表明：

(1)魯中地區氣溶膠EC隨高度總體上呈遞減趨勢。霾時氣溶膠粒子主要聚集在對流層低層2.7 km 以下，重度霾時氣溶膠粒子多集中在1.5 km以下。以0.9、1.66、1.34 km為臨界高度，清潔日、霾日、重度霾日的EC在臨界高度以上呈指數遞減，以下呈對數遞減。年平均PDR和CR總體上隨高度呈遞增趨勢，范圍為0.1～0.3和0.5～0.9。2 km以下，重度霾日氣溶膠以規則的小粒子為主；霾日氣溶膠粒子則以規則的大粒子為主。2～4 km，重度霾日氣溶膠粒子的規則度和大小均為最大，霾日最小。

(2)霾日、重度霾日氣溶膠類型以沙塵型、污染沙塵型為主。清潔日沙塵型氣溶膠最多。清潔日、霾日、重度霾日隨著高度的增加，污染沙塵型比例逐漸減小，沙塵型逐漸增多。霾日和重度霾日，近地面2 km以下存在一部分污染大陸型氣溶膠。霾日2～6 km煙塵型氣溶膠最多。

(3)污染沙塵型對霾日和重度霾日消光的貢獻最大，其次是大陸污染型。沙塵型在近地面1 km以下對重度霾有一定的貢獻。煙塵型在2～4 km對霾日的消光有一定的作用。

(4)垂直方向上風速和相對濕度的分布與氣溶膠層的垂直分布有很好的對應關系。霾時風向隨高度的增加呈順時針旋轉，易出現暖平流。500 hPa霾日與重度霾日高空以西風為主，重度霾日的高空氣流更為平直；850 hPa霾日以西南風為主，重度霾日的北風分量有利于北方的沙塵輸送；1 000 hPa西南風有利于外地污染物的輸送，促進重度霾的發生。

(5)后向軌跡聚類分析結果表明，魯中地區霾時近地面的污染物不僅來源于本地污染排放和河北、江蘇等省份的傳輸，也存在自蒙古、內蒙古等地的長距離傳輸。重度霾期間，蒙古方向的傳輸距離長且高度高，來自內蒙古的傳輸量占比更大，來自河南方向的污染物傳輸量占比最大。

(6)霾的強度越大，多層逆溫出現的比例越高。重度霾對應的多層逆溫比例接近一半。當逆溫型為貼地+低懸型和貼地+高懸型逆溫時，逆溫層下的大氣最穩定，有利于污染物的聚集。