楊軍,高雅,吳彬貴,董琪如,王兆宇,胡漢峰

摘要 2016/2017冬季在天津開展了平流霧微物理特征觀測試驗，結(jié)合距地66 m高度處霧滴譜和255 m氣象塔大氣邊界層資料，借助突變和趨勢一致性非參數(shù)檢驗方法對重度霾后接連發(fā)生的兩次平流霧過程發(fā)展階段進行客觀劃分，揭示霧體內(nèi)部一定高度處霧滴微物理特征和尺度分布特征的觀測事實，討論其生消演變規(guī)律。結(jié)果表明，伴隨西南暖濕平流，飽和層首先在空中出現(xiàn)并向地面擴展，霧過程中成熟階段觀測高度范圍內(nèi)升溫，霧層處于中性或弱不穩(wěn)定層結(jié)狀態(tài)。66 m高度處大霧滴持續(xù)存在，微物理特征與地面能見度準同步變化，數(shù)濃度高值出現(xiàn)在成熟階段初期，而含水量、特征直徑高值出現(xiàn)在成熟階段后期，對應成熟階段后期霧滴數(shù)濃度減少、地面能見度小幅躍升。消散階段各尺度數(shù)濃度因霧滴蒸發(fā)同步減小。

關(guān)鍵詞 平流霧;微物理特征;尺度譜分布;霾

霧是具有垂直結(jié)構(gòu)特征的大氣水物質(zhì)演變現(xiàn)象，霧滴尺度譜分布、含水量等相關(guān)微物理特性是霧的最基本物理屬性，其垂直分布決定了整個霧層的生消演變規(guī)律，是霧數(shù)值模型開發(fā)、預報、遙感等應用研究的科學基礎，也是當前面臨的瓶頸（Gultepe et al.，2007）。

長期以來，霧滴微物理特征的直接觀測主要在地面（即霧層的底部）進行（李子華等，2011;陸春松等，2011;趙麗娟和牛生杰，2012;呂晶晶等，2014;費冬冬等，2016;董瑩等，2021），對霧層內(nèi)部微物理特征的認識十分有限。到目前為止，限于觀測條件，對空中霧體內(nèi)部霧滴尺度譜分布及其微物理特征量的認識僅有以下零星結(jié)果。1959年秋季，Okita （1962）利用氣球（升限約300 m）在日本北海道測量了山地濃霧的液水含量、霧滴數(shù)濃度隨高度的分布，多數(shù)個例顯示含水量和數(shù)濃度在霧層中上部較大，這種垂直結(jié)構(gòu)有利于地面毛毛雨的形成。NASA支持Cornell航空實驗室（CAL）于1970年秋季在紐約州開展了以人工消霧為目的的外場觀測實驗，利用17 m氣象塔和飛機獲得了霧微物理特征的垂直結(jié)構(gòu)，顯示局地生成的輻射霧隨高度和持續(xù)時間增加其譜寬變窄、平均半徑減小（Pilié et al.，1975）。與上述結(jié)果不同的是，1974和1975年夏季在美國舊金山西海岸利用250 m電視塔觀測了平流霧譜分布隨高度的變化特征，發(fā)現(xiàn)無論是平均尺度還是液水含量均隨高度升高而增大（Goodman，1977）。隨后Pinnick et al.（1978）于1976年冬季在德國Grafenwhr利用氣球（升限約250 m，濃霧期間實際約150 m）攜帶前向散射粒譜儀也測得從地面至150 m高度（未能達到霧頂）霧滴譜明顯變寬、數(shù)濃度增大。1992年在美國數(shù)地進行的飛機探測表明，平流霧峰值直徑在低層為4～8 μm，向上往霧體高層則增大到12～20 μm，含水量最大值亦出現(xiàn)在霧體中上層，而輻射霧大滴和高數(shù)濃度出現(xiàn)在近地面50 m高度范圍內(nèi)（Zak，1994）。由于霧微物理垂直結(jié)構(gòu)在遙感探測中的重要性，加之與云微物理垂直結(jié)構(gòu)的差異，Egli et al.（2015）利用氣球攜帶改進的云滴譜探頭，并結(jié)合調(diào)頻連續(xù)波云雷達、云高儀、能見度儀等，在德國Linden進行霧滴譜的垂直探測，通過霧滴譜確定液水含量的垂直廓線。結(jié)果顯示，含水量的垂直結(jié)構(gòu)在部分霧過程的成熟階段與理論廓線較為一致，而在發(fā)展和消散階段則差異較大。類似目的的觀測試驗還在德國波河谷地（Fuzzi et al.，1992;Heintzenberg，1992）開展，得到兩個高度的微物理特征。除此之外，鮮見相似研究報道。

以上零星研究得到了空中霧滴譜及其積分特征量垂直結(jié)構(gòu)的初步結(jié)果，但結(jié)論不完全一致，且有限的認識尚無法滿足相關(guān)科學問題和業(yè)務應用的需要。例如在實際應用中，多采用均勻的垂直分布（Hutchison，2002）、或者可逆絕熱含水量廓線、或者直接利用云厚和光學厚度建立垂直廓線的簡單經(jīng)驗關(guān)系（Brenguier et al.，2000）。而且不同地理環(huán)境、霧的不同發(fā)展階段、不同物理機制組合造成了霧垂直結(jié)構(gòu)更具有多樣性，針對不同霧過程及生消階段進一步認識空中霧微物理垂直結(jié)構(gòu)的基本規(guī)律是相關(guān)研究的基礎。

針對空中霧微物理特征探測的迫切性（Gultepe et al.，2007;Niu et al.，2010），2016/2017年冬季在天津開展了霧過程外場觀測試驗，本文選取2016年12月31日至2017年1月1日重度霾后接連發(fā)生的兩次平流霧過程，將距地面66 m高度霧微物理探測與氣象塔邊界層探測相結(jié)合，在對霧過程發(fā)展階段客觀劃分的基礎上，討論霧過程的發(fā)展規(guī)律，揭示關(guān)于霧層內(nèi)部一定高度霧滴微物理特征的觀測事實，增加對霧微物理演變規(guī)律的科學認識。

1 觀測與資料處理

1.1 觀測概況

觀測于2016/2017年冬季進行，F(xiàn)M-120型霧滴譜儀位于天津市氣象局預警中心大樓樓頂平臺，距地高度66 m。大氣邊界層資料來自天津城市氣候監(jiān)測站255 m大氣邊界層氣象塔（117°12′E，39°04′N），塔上120、200 m高度處設置大氣顆粒物監(jiān)測儀，常規(guī)氣象要素15層。地面設備包括前向散射能見度儀和常規(guī)要素（表1）。

觀測區(qū)域周邊為城市建成區(qū)下墊面，東邊距海岸線約58 km。樓頂平臺與氣象塔相距2.9 km，考慮本文研究時段背景大氣條件和平流霧的水平均勻性（如表2），認為兩位置66 m高度處大氣狀態(tài)具有一致性，下文分析中地面能見度與空中霧微物理參數(shù)的時間演變規(guī)律也表明了該一致性。

霧滴譜資料采用切尾均值（Rosenberger and Gasko，1983）得到分鐘數(shù)濃度尺度譜分布，在此基礎上計算液水含量等統(tǒng)計特征量。天氣形勢分析使用NCEP全球再分析資料和MICAPS地面常規(guī)氣象資料。

1.2 霧發(fā)展階段客觀劃分方法

對霧過程發(fā)展階段的客觀劃分，有助于準確揭示其階段性微物理特征和演變機制。但長期以來，霧過程發(fā)展階段多為主觀劃分。Maier et al.（2013）采取氣象、水文等領(lǐng)域時間序列突變分析中廣泛應用的非參數(shù)統(tǒng)計方法（Kundzewicz and Robson，2004;Li et al.，2005），利用霧過程的能見度和霧滴微物理參數(shù)觀測結(jié)果，提出了霧過程發(fā)展階段的客觀劃分方法。

該方法分3步進行：1）可能變點識別：要素時間序列發(fā)生趨勢性改變的變點是過程發(fā)展階段劃分的基礎，可采用Mann and Whitney（1947）提出并經(jīng)Pettitt（1979）完善的Mann-Whitney同質(zhì)性檢驗方法識別可能變點;2）變化趨勢顯著性檢驗：從可能變點位置將時間序列分割為兩段序列，使用Mann-Kendall趨勢檢驗分別確認兩段曲線是否存在顯著的變化趨勢，若存在顯著變化趨勢則進行下一步;3）變化趨勢一致性檢驗：對于具有顯著變化趨勢的兩段序列，對要素-時間雙累積曲線進行雙側(cè)t檢驗，確認兩段序列變化趨勢是否具有顯著差異。若具有顯著差異，則認為第1步識別的可能變點為過程發(fā)展階段的轉(zhuǎn)折點。以上檢驗顯著性水平均取0.05。

對于通過檢驗的變點分割的兩段序列，再次分別進行變點識別，直至序列段中無符合條件的變點為止。利用得到的可能變點，結(jié)合要素序列及其雙累積曲線，可對要素序列發(fā)展階段進行客觀劃分，從而更準確地得到霧過程階段性特征。

2 霧過程概況

2016年12月29日至2017年1月1日華北、黃淮出現(xiàn)大范圍霧和霾過程，前期以霾為主，12月31日至次年1月1日天津到河北中南部大范圍濃霧發(fā)展，局地能見度低于50 m，天津受到重度霾和特強濃霧影響。

12月31日華北地區(qū)500 hPa受平直西風氣流控制，脊前負渦度平流對地面的加壓作用促進了地面均壓場的維持。天津位于地面高壓西北側(cè)（圖1a），低層反氣旋環(huán)流外圍西南氣流引導暖濕空氣向天津地區(qū)輸送，并在夜里加強（圖1b），為霧的形成與維持提供了水汽條件。31日全天，從北京、天津往西南方向一直到邯鄲均被大范圍霧區(qū)持續(xù)覆蓋，天津城市效應使得該大范圍霧區(qū)東北部形成無霧缺口（圖2a）。

31日下午隨著低層西南氣流的平流輸送作用，霧區(qū)進一步向東北方向發(fā)展。從圖2b可以看到，1月1日09：00（北京時，下同）均勻霧區(qū)覆蓋了整個天津市范圍。

本文觀測點地面水平能見度低于1 000 m的時段從31日15：23開始，至1日12：41結(jié)束（圖3），以該時段作為本文研究的時間范圍?？紤]1日01：00前后能見度兩次升高到893 m（1日00：13及01：27），該時間前后霧過程均經(jīng)歷了發(fā)展、成熟、減弱的過程，因此將研究時間范圍進一步劃分為兩次霧過程：F1和F2，時間范圍分別為31日15：23—1日00：13、1日01：27—12：41。

從氣象塔觀測結(jié)果（圖4，數(shù)據(jù)采集器故障導致16：00—00：01缺測）來看，霧過程及其前后觀測點近地層均以西南氣流為主。31日霧前05：00—08：00底層明顯降溫、穩(wěn)定度增大并出現(xiàn)等溫、逆溫現(xiàn)象。隨后塔層上部降溫，到14：15塔層高度范圍內(nèi)氣溫降至本過程最低狀態(tài)，為平流霧的發(fā)展提供了冷的近地層條件。14：15以后，塔層范圍內(nèi)持續(xù)升溫，能見度于15：23降到1 000 m以下，F(xiàn)1過程結(jié)束前的00：10與開始時刻15：23相比，除250 m高度升溫1.4 ℃外，其下各層均升溫2.1～2.4 ℃。吳彬貴等（2008）和蔡子穎等（2012）研究天津平流霧的生消及不同發(fā)展階段的溫、濕、風垂直結(jié)構(gòu)時指出，西南暖濕平流和西北干冷平流對平流霧的生消起重要作用。

從相對濕度來看，霧前塔層上部先存在飽和層，伴隨降溫過程自上而下相對濕度逐漸增大，30、20、10、5 m高度先后于07：55、08：02、08：43、09：26達到飽和。雖然濕度達到飽和，但由于處于重度霾狀態(tài)，霧前能見度一直處于1 000～2 000 m。

F2過程03：00前后120 m以下有微弱降溫，但與F1過程相同的是，隨后濃霧維持期間同樣伴隨溫度上升，整個F2過程塔層高度范圍內(nèi)均處于飽和狀態(tài)。從06：30開始，個別高度開始有弱西北風出現(xiàn)，溫度小幅下降，上層降溫幅度略大于近地層。溫度廓線及其遞減率在濃霧階段及其前后較為一致，推斷霧頂應高于200 m，即微物理觀測所在高度處于霧層中下部。

研究時段當?shù)厝粘?、日落時間分別約為07：30、17：00，以上近地層溫度變化基本沒有顯示出太陽短波輻射加熱的日變化特點。濃霧階段伴隨西南暖濕氣流，觀測高度范圍內(nèi)升溫，且溫度層結(jié)接近或略大于對流層平均減溫率，處于中性或弱不穩(wěn)定層結(jié)狀態(tài)。以上特點結(jié)合天氣形勢來看，本文研究對象主要表現(xiàn)為平流霧特征。

3 微物理特征

利用2.2節(jié)階段劃分方法，對F1、F2兩次霧過程2 m能見度和66 m霧滴液水含量、數(shù)濃度和平均半徑（分鐘值）進行變點分析，得到霧過程發(fā)展階段轉(zhuǎn)換時刻（圖5）。因此將兩次霧過程按地面能見度分別劃分為3個階段，以F1A、F1B、F1C和F2A、F2B、F2C表示，A、B、C分別對應形成、成熟、消散階段。

兩次霧過程能見度時間變化的共同特點有：F1過程21：02—22：49、F2過程08：04—10：07兩時段能見度均較前一階段有小幅度升高但明顯低于消散階段，對應空中含水量處于整個過程的高值階段。分析F2過程的塔層廓線，從06：30開始，個別高度開始有弱偏北風出現(xiàn)，溫度小幅下降，在水汽充足條件下這種降溫對應了低能見度的維持以及08：00—10：00霧滴凝結(jié)增大、含水量增加。

圖5顯示，F(xiàn)1過程含水量、數(shù)濃度開始顯著增大的時刻（17：45、17：39）與地面能見度降到最低階段的時刻（17：54）基本吻合、僅略早10 min左右。含水量、數(shù)濃度最后一個變點（22：57、22：56）與能見度快速升高時刻偏晚僅在10 min之內(nèi)。F2過程含水量、數(shù)濃度、表面積平均直徑最后變點與霧過程結(jié)束前能見度快速增大時刻（10：07）一致。該過程能見度第1個變點（02：20）與數(shù)濃度和含水量（02：36）也基本匹配。

上述地面能見度與空中微物理參數(shù)的變點對應關(guān)系表明本文霧過程在地面至66 m甚至更高范圍內(nèi)具有同步生消的特點，再次說明研究時段大氣狀態(tài)水平較均勻，同時也表明該變點分析方法對霧過程不同要素具有適用性。

從圖5還可進一步看到，兩次霧過程含水量與表面積平均直徑的高值均出現(xiàn)在成熟階段的后期，與之相反數(shù)濃度高值出現(xiàn)在成熟階段初始時段，這是重污染大氣背景大量顆粒物競爭水分與霧過程濕清除的共同結(jié)果。Goodman（1977）在大陸西海岸夏季平流霧研究中也同樣觀測到霧滴平均直徑隨時間持續(xù)增長至霧消散前一階段的現(xiàn)象。

表3給出了由分鐘數(shù)據(jù)統(tǒng)計的霧過程不同階段能見度與微物理特征量，每一過程的時間范圍基于圖5中微物理特征量的變點，并結(jié)合成熟階段的要素特征確定，與地面能見度階段劃分有數(shù)分鐘至十幾分鐘的差異。從表中可以看到，兩次霧過程在3個不同階段能見度平均值相近，其中F2過程成熟階段平均能見度略低于F1過程成熟階段，但F2過程最低能見度略大于F1過程。F2過程略低的平均能見度對應了相比F1過程略高的含水量和特征直徑，但F2過程霧滴數(shù)濃度較F1過程小，這可能是由于F1過程顆粒物濃度高（圖6），導致霧滴數(shù)濃度大，但尺度偏小。另從圖6還可看到，隨著F1過程開始，PM2.5和PM10均迅速下降，表明重度霾背景條件下平流霧的有效清除作用。這種清除效應在F2過程消散前亦表現(xiàn)明顯，F(xiàn)2過程結(jié)束后顆粒物濃度又有所上升。

4 霧滴尺度譜分布特征

為揭示空中霧滴的尺度特征，圖7給出了基于分鐘數(shù)據(jù)得到的兩次霧過程不同發(fā)展階段的平均數(shù)濃度尺度譜分布，各階段時間范圍如表3。由圖可見，小粒子端存在4～5 μm的數(shù)濃度峰值，特別是受前期重度霾影響的F1A階段也是如此。這不同于大量局地輻射霧形成階段的譜分布規(guī)律，即核化過程形成霧滴，其后逐步增長，開始階段霧滴數(shù)濃度譜傾向于隨尺度增大而遞減（Price，2011）。這說明本文霧過程從開始便受到上游霧滴平流輸入的影響。兩次過程3個階段均在18～20 μm附近出現(xiàn)數(shù)濃度峰值，這一現(xiàn)象還有待進一步研究。

雖然兩次霧過程66 m高度處微物理特征量與其他地面相對清潔環(huán)境的觀測結(jié)果相比偏低，但從數(shù)濃度尺度譜分布（圖7）來看，大滴依然存在，即成熟階段存在直徑40 μm以上的大滴，

這與大量平流霧地面觀測結(jié)果一致。Goodman （1977）觀測發(fā)現(xiàn)毛毛雨滴在海岸平流霧（低云）上部形成并向下沉降，使得各高度均有大滴存在。而對輻射霧微物理特性垂直分布觀測結(jié)果，Egli et al.（2015）認為由于重力沉降作用，霧層下部霧滴譜分布更寬。本文霧過程與Goodman （1977）不同的是，背景污染程度大，霧區(qū)上游位于河北南部，水汽條件較之夏季向岸海洋濕團弱，因此大霧滴存在但含量偏低，沒有毛毛雨出現(xiàn)。

消散階段與成熟階段譜型相似，各尺度霧滴數(shù)濃度同步增減，表明從成熟階段到消散階段各尺度霧滴均蒸發(fā)變小變少。F1消散階段各尺度霧滴數(shù)濃度大于形成階段，這與兩次過程連續(xù)發(fā)生有關(guān)，F(xiàn)1C結(jié)束和F2A開始的能見度均在1 000 m以下。兩次過程消散階段大滴均多于形成階段，與圖5中霧滴尺度持續(xù)增長至消散階段前相對應。消散階段與形成階段相比大滴多、小滴少，說明消散階段能見度從數(shù)百米回升到1 000 m以上的過程中，小滴減少的比大滴快，同樣反映的是霧滴蒸發(fā)的結(jié)果。

5 結(jié)論與討論

通過對天津兩次連續(xù)平流霧過程距地66 m高度霧滴譜及255 m氣象塔外場試驗結(jié)果分析，得到如下觀測事實：

1）飽和層首先自空中向地面擴展，隨后66 m高度霧微物理特征與地面能見度同步演變，表明平流霧過程可在一定高度范圍內(nèi)準同步發(fā)展。

2）濃霧階段伴隨西南暖濕氣流，觀測高度范圍內(nèi)升溫，溫度層結(jié)接近或略大于對流層平均減溫率，霧層處于中性或弱不穩(wěn)定層結(jié)狀態(tài)。

3）數(shù)濃度高值出現(xiàn)在霧過程成熟階段初期，成熟階段后期含水量、特征直徑持續(xù)增長到過程最大值，但霧滴數(shù)濃度減少，對應地面能見度小幅升高。

4）66 m高度霧滴部分來自平流輸入，雖然微物理特征參數(shù)值偏低，但大霧滴持續(xù)存在。消散階段與成熟階段譜型相似，即消散過程各尺度霧滴數(shù)濃度同步蒸發(fā)減少。

參考文獻（References）

Brenguier J L，Pawlowska H，Schüller L，et al.，2000.Radiative properties of boundary layer clouds：droplet effective radius versus number concentration[J].J Atmos Sci，57（6）：803-821.doi：10.1175/1520-0469（2000）057<0803：rpoblc>2.0.co;2.

蔡子穎，韓素芹，吳彬貴，等，2012.天津一次霧過程的邊界層特征研究[J].氣象，38（9）：1103-1109. Cai Z Y，Han S Q，Wu B G，et al.，2012.Analysis on characteristics of atmospheric boundary layer during a fog process in Tianjin[J].Meteor Mon，38（9）：1103-1109.（in Chinese）.

董瑩，尹志聰，段明鏗，2021.長三角地區(qū)冬季霾污染日數(shù)的季節(jié)預測[J].大氣科學學報，44（2）：290-301. Dong Y，Yin Z C，Duan M K，2021.Seasonal prediction of winter haze days in the Yangtze River Delta[J].Trans Atmos Sci，44（2）：290-301.doi：10.13878/j.cnki.dqkxxb.20200525001.（in Chinese）.

Egli S，Maier F，Bendix J，et al.，2015.Vertical distribution of microphysical properties in radiation fogs：a case study[J].Atmos Res，151：130-145.doi：10.1016/j.atmosres.2014.05.027.

費冬冬，牛生杰，楊軍，等，2016.長江中上游冬季山地霧邊界層特征及生消過程分析[J].大氣科學學報，39（2）：221-231. Fei D D，Niu S J，Yang J，et al.，2016.Boundary layer characteristics and formation processes of winter valley fog in the upper and middle reaches of the Yangtze River[J].Trans Atmos Sci，39（2）：221-231.doi：10.13878/j.cnki.dqkxxb.20140519001.（in Chinese）.

Fuzzi S，F(xiàn)acchini M C，Orsi G，et al.，1992.The Po valley fog experiment 1989[J].Tellus B：Chem Phys Meteor，44（5）：448-468.doi：10.3402/tellusb.v44i5.15561.

Goodman J，1977.The microstructure of California coastal fog and stratus[J].J Appl Meteor，16（10）：1056-1067.doi：10.1175/1520-0450（1977）016<1056：tmoccf>2.0.co;2.

Gultepe I，Tardif R，Michaelides S C，et al.，2007.Fog research：a review of past achievements and future perspectives[J].Pure Appl Geophys，164（6/7）：1121-1159.doi：10.1007/s00024-007-0211-x.

Heintzenberg J，1992.The Po Valley Fog Experiment 1989：what have we learned，where do we go from here？[J].Tellus B：Chem Phys Meteorol，44（5）：443-447.doi：10.3402/tellusb.v44i5.15558.

Hutchison K D，2002.The retrieval of cloud base heights from MODIS and three-dimensional cloud fields from NASAs EOS Aqua mission[J].Int J Remote Sens，23（24）：5249-5265.doi：10.1080/01431160110117391.

Kundzewicz Z W，Robson A J，2004.Change detection in hydrological records：a review of the methodology[J].Hydrolog Sci J，49（1）：7-19.

Li Y，Cai W，Campbell E P，2005.Statistical modeling of extreme rainfall in southwest western Australia[J].J Climate，18（6）：852-863.doi：10.1175/jcli-3296.1.

李子華，劉端陽，楊軍，等，2011.南京市冬季霧的物理化學特征[J].氣象學報，69（4）：706-718. Li Z H，Liu D Y，Yang J，et al.，2011.Physical and chemical characteristics of winter fogs in Nanjing[J].Acta Meteorol Sin，69（4）：706-718.（in Chinese）.

陸春松，牛生杰，岳平，等，2011.南京冬季霧多發(fā)期邊界層結(jié)構(gòu)觀測分析[J].大氣科學學報，34（1）：58-65. Lu C S，Niu S J，Yue P，et al.，2011.Observational research on boundary layer structure during high incidence period of winter fog in Nanjing[J].Trans Atmos Sci，34（1）：58-65.doi：10.13878/j.cnki.dqkxxb.2011.01.008.（in Chinese）.

呂晶晶，牛生杰，張羽，等，2014.湛江東海島一次春季海霧的宏微觀結(jié)構(gòu)及邊界層演變特征[J].氣象學報，72（2）：350-365. Lü J J，Niu S J，Zhang Y，et al.，2014.Evolution characteristics of the macro-/micro-structure and the boundary layer during a spring heavy sea fog episode in Donghai Island in Zhanjiang[J].Acta Meteorol Sin，72（2）：350-365.doi：10.11676/qxxb2014.038.（in Chinese）.

Maier F，Bendix J，Thies B，2013.Development and application of a method for the objective differentiation of fog life cycle phases[J].Tellus B：Chem Phys Meteorol，65（1）：19971.doi：10.3402/tellusb.v65i0.19971.

Mann H B，Whitney D R，1947.On a test of whether one of two random variables is stochastically larger than the other[J].Ann Math Statist，18（1）：50-60.doi：10.1214/aoms/1177730491.

Niu S J，Lu C S，Yu H Y，et al.，2010.Fog research in China：an overview[J].Adv Atmos Sci，27（3）：639-662.doi：10.1007/s00376-009-8174-8.

Okita T，1962.Observations of the vertical structure of a stratus cloud and radiation fogs in relation to the mechanism of drizzle formation[J].Tellus，14（3）：310-322.doi：10.3402/tellusa.v14i3.9556.

Pettitt A N，1979.A non-parametric approach to the change-point problem[J].J Royal Stat Soc：Ser C Appl Stat，28（2）：126-135.doi：10.2307/2346729.

Pilié R J，Mack E J，Kocmond W C，et al.，1975.The life cycle of valley fog.part II：fog microphysics[J].J Appl Meteor，14（3）：364-374.doi：10.1175/1520-0450（1975）014<0364：tlcovf>2.0.co;2.

Pinnick R G，Hoihjelle D L，F(xiàn)ernandez G，et al.，1978.Vertical structure in atmospheric fog and haze and its effects on visible and infrared extinction[J].J Atmos Sci，35（10）：2020-2032.doi：10.1175/1520-0469（1978）035<2020：vsiafa>2.0.co;2.

Price J，2011.Radiation fog.part I：observations of stability and drop size distributions[J].Bound Layer Meteor，139（2）：167-191.doi：10.1007/s10546-010-9580-2.

Rosenberger J L，Gasko M，1983.Comparing location estimators：trimmed means，medians，and trimean [C]//Understanding Robust and Exploratory Data Analysis.Wiley：297-336.

吳彬貴，張宏升，張長春，等，2008.天津南部地區(qū)平流霧過程塔層氣象要素特征分析[J].北京大學學報（自然科學版），44（5）：744-750. Wu B G，Zhang H S，Zhang C C，et al.，2008.Analysis of the micro-meteorologic element during the advection fog period in the south of Tianjin city[J].Acta Sci Nat Univ Pekin，44（5）：744-750.doi：10.13209/j.0479-8023.2008.115.（in Chinese）.

Zak J A，1994.Drop size distributions and related properties of fog for five locations measured from aircraft[R].Contractor Report NASA CR-4585，DOT/FAA/CT-94/02：136.

趙麗娟，牛生杰，2012.近地層輻射過程與霧微結(jié)構(gòu)的相互作用特征[J].大氣科學學報，35（6）：673-679. Zhao L J，Niu S J，2012.Characteristics of interactions between radiation processes and fog microphysical structure[J].Trans Atmos Sci，35（6）：673-679.doi：10.13878/j.cnki.dqkxxb.2012.06.008.（in Chinese）.

Advection fog process and its microphysical properties：a case study in Tianjin

YANG Jun1，GAO Ya1，WU Bingui2，DONG Qiru2，WANG Zhaoyu3，HU Hanfeng1

1Key Laboratory for Aerosol-Cloud-Precipitation of China Meteorological Administration，Nanjing University of Information Science & Technology，Nanjing 210044，China;

2Tianjin Institute of Meteorological Science，Tianjin 300074，China;

3Tianjin Weather Modification Office，Tianjin 300074，China

The microphysical characteristics observation test of advection fog was carried out in Tianjin in the winter of 2016/2017.Combined with the droplet spectrum at a height of 66 m above the ground and the atmospheric boundary layer data of 255 m meteorological tower，the development stages of two consecutive advection fog processes after heavy haze were objectively divided by means of change-point detection and trend consistency nonparametric test，which reveals the observation facts of droplet microphysical characteristics and size distribution characteristics at a certain height in the fog，and discusses the evolution law of its generation and disappearance.The results show that with the southwest warm and humid advection，the saturated layer first appears in the air and extends to the ground.During the fog process，the temperature rises within the observation height range in the mature stage，and the fog layer is in a neutral or weakly unstable stratification state.At the height of 66 m，large fog droplets persist，and the microphysical characteristics change quasi synchronously with the ground visibility.The high value of number concentration appears in the early stage of maturity，while the high value of water content and characteristic diameter appears in the late stage of maturity，corresponding to the decrease of fog droplet number concentration and a small jump in ground visibility in the late stage of maturity.In the dissipation stage，the number concentration of each size decreases synchronously due to droplet evaporation.

advection fog;microphysical property;size distribution;haze

doi：10.13878/j.cnki.dqkxxb.20200516012

（責任編輯：劉菲）