WRF模式對烏魯木齊機(jī)場一次凍霧天氣的數(shù)值預(yù)報對比試驗分析

陳陽權(quán)，杜安妮，竇新英，王清平，張利平

（1.民航新疆空中交通管理局氣象中心，新疆 烏魯木齊830016；2.新疆氣象臺，新疆 烏魯木齊830002）

霧的形成受到天氣系統(tǒng)、本地溫度、相對濕度、風(fēng)向、風(fēng)速、大氣穩(wěn)定度、大氣成分、地形等諸多因素的影響，其預(yù)報是一個非常復(fù)雜的難題。目前對大霧的預(yù)報主要依靠天氣分型法、指標(biāo)法、傳統(tǒng)的統(tǒng)計學(xué)方法等，在實際業(yè)務(wù)應(yīng)用中已經(jīng)取得了一定效果，但這些方法嚴(yán)重依賴觀測資料的時間和空間分辨率，在一些資料匱乏的地區(qū)難以應(yīng)用，對于一些局地性的大霧預(yù)報較難。近年來，國內(nèi)學(xué)者已經(jīng)開始利用中尺度數(shù)值預(yù)報模式對霧進(jìn)行試驗研究，得到了有意義的研究成果[1-3]。已有的這些研究成果表明，霧的數(shù)值預(yù)報水平與模式的各物理過程參數(shù)化方案密切相關(guān)，這些參數(shù)化方案包括微物理過程、邊界層過程、陸面過程、垂直分辨率等。

WRF（Weather Research and Forecasting）模式作為新一代中尺度數(shù)值預(yù)報模式，提供了多種完善的參數(shù)化方案可供選擇，模式分辨率在初始場分辨率較高的條件下可以精確到3 km或者更細(xì)，垂直分辨率甚至可以達(dá)到幾十米，被廣泛應(yīng)用于降水、大風(fēng)等業(yè)務(wù)系統(tǒng)，也被廣泛應(yīng)用于霧的診斷與預(yù)報試驗[1]。在民航氣象領(lǐng)域，部分地區(qū)空管局氣象中心已經(jīng)開始利用WRF進(jìn)行低云、低能見度預(yù)報業(yè)務(wù)化系統(tǒng)的試驗。

在大霧的模擬研究中，常將特定的液態(tài)含水量臨界值作為判斷霧區(qū)的條件之一[4-6]，有時也選定相對濕度的閾值作為判斷大霧出現(xiàn)和消散的條件[7]。不同的地域特點及不同的出霧條件決定了各地霧中液態(tài)含水量和相對濕度條件必然各不相同。據(jù)統(tǒng)計，新疆80豫耀90豫的大霧天氣出現(xiàn)在相對濕度>90豫的環(huán)境，5豫耀20豫的大霧出現(xiàn)在相對濕度80豫耀90豫的環(huán)境下，低于80豫的相對濕度環(huán)境下很少出現(xiàn)大霧天氣[8]。王春紅等[9]對烏魯木齊機(jī)場一次初霧天氣的模擬分析中選用地面濕度逸80%作為出現(xiàn)大霧的潛在區(qū)域，并取得了較好的效果。

航空氣象上霧是指近地面中水汽凝結(jié)或凝華而使能見度低于1000 m的視程障礙現(xiàn)象，凍霧是霧的一種形態(tài)，包括凍結(jié)和過冷卻兩種，凍結(jié)指在地面上形成了霧凇，過冷卻指地面上未凍結(jié)，但霧中由過冷卻水滴組成。烏魯木齊機(jī)場出現(xiàn)能見度低于800 m的凍霧時就會影響航空器運行正常，而每年11月—次年3月，都會頻繁出現(xiàn)能見度低于800 m的凍霧天氣，王楠等[10]通過對烏魯木齊機(jī)場37 a霧的特征進(jìn)行分析發(fā)現(xiàn)，機(jī)場霧累積出現(xiàn)時間呈上升趨勢。據(jù)機(jī)場氣象部門統(tǒng)計，在2015年12月、2016年1月分別出現(xiàn)了27 d和29 d的凍霧，凍霧日數(shù)突破歷史同期的極值，凍霧成為冬季影響機(jī)場航班安全和正常的重要復(fù)雜天氣。本文利用WRF模式對2016年11月28—29日烏魯木齊機(jī)場長時間凍霧天氣過程進(jìn)行數(shù)值預(yù)報，以2 m相對濕度逸80%作為評估參數(shù)，機(jī)場實況進(jìn)行對比，研究WRF模式在凍霧預(yù)報過程中對各種參數(shù)化方案的敏感性，選擇出最優(yōu)的參數(shù)化方案組合，并在此基礎(chǔ)上分析預(yù)報結(jié)果，期望直接利用數(shù)值預(yù)報診斷出此次凍霧過程，為凍霧預(yù)報積累經(jīng)驗。

1 天氣過程

2016年11月28日10時—29日06時（文中時間均為世界協(xié)調(diào)時UTC），北疆沿天山一帶的烏蘭烏蘇至烏魯木齊機(jī)場一帶出現(xiàn)了大范圍的凍霧天氣，能見度普遍在100 m左右。烏魯木齊機(jī)場28日10時，開始出現(xiàn)凍霧天氣，至29日06時能見度迅速好轉(zhuǎn)到900 m，凍霧持續(xù)時間達(dá)20 h以上，能見度最低100 m，對航空運行造成了極大影響，凍霧期間主導(dǎo)能見度變化如圖1所示。凍霧在28日10時出現(xiàn)，實屬罕見。烏魯木齊市區(qū)從29日03時開始能見度迅速降低至800~1000 m，在07時后能見度迅速下降至500 m及以下，最低100 m。

圖1 2016年11月28—29日凍霧期間主導(dǎo)能見度變化

2 資料與方法

2.1 所用資料包括

2016年11月28—29日機(jī)場METAR報文資料及新疆地區(qū)逐小時自動站資料；28日00時UTC GFS 1毅伊 1毅預(yù)報場資料。

2.2 預(yù)報試驗設(shè)計

采用WRF3.7.1版本，使用三重雙向嵌套網(wǎng)格，中心點位于烏魯木齊附近（44毅N，87毅E），網(wǎng)格分辨率分別為27、9、3 km，區(qū)域設(shè)計見圖2a。垂直層次為45層，對底層進(jìn)行加密，3 km以下15層，垂直層次見圖2b。以28日00時GFS全球預(yù)報場作為初始場，積分步長90 s，積分48 h，大區(qū)域每3 h輸出一次結(jié)果，小區(qū)域每1 h輸出一次結(jié)果，分析中以前6 h作為SPIN-UP時間。

圖2 區(qū)域嵌套設(shè)計及垂直層次的高度變化

目前，WRF模式參數(shù)化方案主要包括微物理方案、長波輻射方案、短波輻射方案、近地層方案、陸面過程方案、邊界層方案及積云對流方案等。由于積云對流參數(shù)化方案對新疆區(qū)域冬季影響較小，本文選擇KF方案（Kain-Fritsch（new Eta）scheme）。長波輻射方案選用RRTM方案，短波輻射方案選用Dudhia方案[11-15]。由于霧是近地層的視程障礙現(xiàn)象，且其主要發(fā)生在穩(wěn)定的邊界層中，其對下墊面性質(zhì)極為敏感[16]，因此，本文將重點研究近地層方案、邊界層方案、陸面過程方案的組合對霧的影響。在利用模式模擬霧的研究中，經(jīng)常用到的陸面過程方案有SLAB方案（5-layer thermal diffusion scheme）和Noah方案（Unified Noah land-surface model）；邊界層方案主要有YSU方案 （Yonsei University scheme）、MYJ方案（Mellor-Yamada-Janjic TKE scheme）和QNSE方案（Quasi-Normal Scale Elimination PBL），與之對應(yīng)的近地層方案分別是MM5 similarity方案、Eta similarity方案和QNSE surface layer方案；微物理過程采用Lin方案（Lin et al.scheme）和WDM6方案（WDM 6-class scheme）。為了找出不同方案對霧的預(yù)報效果，參考黃政[3]等在高速公路大霧過程研究中的參數(shù)化方案組合方法，將不同方案進(jìn)行組合得到12套物理方案選擇。另外，榮昕等[14]在其研究中提到微物理方案選擇Morrison方案（Morrison double-moment scheme）、積云對流方案選擇G3方案（Grell 3D ensemblescheme）、陸面過程方案選擇SLAB方案、邊界層方案選擇QNSE方案時，對霧的模擬效果較好，因此本文也考慮這套方案的預(yù)報效果，方案組合見表1。

表1 各參數(shù)化方案組合試驗設(shè)計[3]

3 不同參數(shù)化方案效果評估

對13個試驗方案組合的500 hPa環(huán)流形勢及海平面氣壓場進(jìn)行對比分析發(fā)現(xiàn)，在28日12時UTC 13個試驗方案預(yù)報新疆西部地區(qū)為脊區(qū)控制，新疆其它大部地區(qū)為脊前西北氣流控制，海平面氣壓場上北疆大部為高壓控制，13個方案預(yù)報的高空和地面環(huán)流形勢完全一致，且預(yù)報高空環(huán)流形勢與實況基本吻合，北疆地面高壓位置與實況也一致，強(qiáng)度略偏小。通過對高空地面形勢分析可排除降水對2 m相對濕度的影響，在13個試驗方案形勢預(yù)報一致的條件下，本文將通過模式預(yù)報的2 m相對濕度與機(jī)場2 m相對濕度實況進(jìn)行對比分析，選擇最優(yōu)的參數(shù)化方案組合。

由表2可知，機(jī)場實況顯示09時濕度超過80%，對比圖1中能見度可知，能見度從1600 m迅速下降至1000 m，至10時能見度下降至200 m，之后長時間維持在100耀200 m。對比13個方案中2 m相對濕度變化來看，組合方案3、9、13預(yù)報的2 m相對濕度在10時均超過80%，且之后一直維持，而其余方案2 m相對濕度在10時均低于80%，組合方案4、5、6、10、11、12等方案整個預(yù)報過程中2 m相對濕度預(yù)報一直低于70%，與實際差異較大。以2 m相對濕度的變化來代表能見度變化，組合方案3、9、13預(yù)報效果相對于其余方案較好。

單獨分析方案3、9、13方案（圖3）在28日19時前預(yù)報基本一致，與實況趨勢也基本一致，但預(yù)報的相對濕度比實況小。19時后，實況出現(xiàn)了上下波動，方案3和9能預(yù)報出這種波動，方案13預(yù)報結(jié)果仍然持續(xù)上升，3個方案在20時后預(yù)報相對濕度均比實況大。與方案3相比，方案9的預(yù)報更接近于實況。由此可見，方案9在這次試驗中是最優(yōu)的。

對比（表2）可知，方案3、9、13近地層方案均為QNSE、陸面方案均為SLAB、邊界層方案均為QNSE，由此可見，WRF模式2 m相對濕度的預(yù)報對近地層方案、陸面過程方案、邊界層方案組合較為敏感，而微物理過程起到微調(diào)作用，Lin方案和WDM6兩者預(yù)報趨于一致，在數(shù)值上略有差異。Morrison方案在預(yù)報后期表現(xiàn)略差，尤其是在實況出現(xiàn)波動時，預(yù)報效果沒有Lin方案和WDM6好。

圖3 方案3、9、13預(yù)報2 m相對濕度與實況對比

表2 機(jī)場28日06時—29日06時實況相對濕度及WRF各參數(shù)化方案預(yù)報2 m相對濕度 %

4 預(yù)報結(jié)果分析

通過對13個方案的對比分析，方案9（即WDM6、QNSE、SLAB、QNSE、KF）在此次試驗中表現(xiàn)最好，利用方案9的預(yù)報結(jié)果，嘗試對天氣形勢、穩(wěn)定層、溫濕度、風(fēng)場混合等條件進(jìn)行分析，期望能直接通過模式預(yù)報診斷出此次凍霧天氣，對將來分析和預(yù)報凍霧提供一定的幫助。以下分析中所參考圖形高度均為sigma層高度，該高度與實際高度關(guān)系見圖2b，具體對應(yīng)數(shù)值見表3。

表3 sigma層高度與實際高度對應(yīng)數(shù)值 m

4.1 形勢預(yù)報

28日09時預(yù)報的500 hPa高空環(huán)流形勢（圖4a）上，北疆處于槽后西北氣流控制，中亞至新疆西部地區(qū)為淺脊區(qū)，500 hPa環(huán)流形勢逐漸趨于穩(wěn)定，且蒙古至新疆東部以東為高空槽區(qū)，高空環(huán)流形勢非常有利于中低層形成穩(wěn)定的東南風(fēng)層，使得逆溫層穩(wěn)定且深厚，為凍霧出現(xiàn)和長時間維持提供非常好的環(huán)流條件。海平面氣壓場的預(yù)報上（圖4b），北疆大部處于高壓控制，北疆沿天山正好處于高壓底部區(qū)域，氣壓梯度小，地面風(fēng)場弱，且在高壓控制下天空放晴，利于輻射降溫而形成輻射霧。

4.2 穩(wěn)定層預(yù)報

穩(wěn)定層結(jié)是凍霧出現(xiàn)和持續(xù)的必要條件之一。從弱減溫層、等溫層至出現(xiàn)逆溫層，層結(jié)越來越穩(wěn)定，凍霧出現(xiàn)和持續(xù)的可能性越來越大。從06時開始，溫度層結(jié)在615 m以下由等溫層逐漸發(fā)展為逆溫層，在17時以后，268耀652 m出現(xiàn)東南風(fēng)層，東南風(fēng)增溫作用使得該層次溫度升高，而東南風(fēng)層以下高度溫度變化不大，逆溫強(qiáng)度加強(qiáng)（圖5a）。由圖5b可知，在09時107耀450 m高度已經(jīng)出現(xiàn)1益的弱逆溫層，與10時出現(xiàn)的凍霧天氣有很好的對應(yīng)關(guān)系，至18時，逆溫層加厚至652 m，逆溫層底部貼近地面，逆溫強(qiáng)度明顯加強(qiáng)超過6益，持續(xù)且強(qiáng)盛的逆溫層對凍霧強(qiáng)度加強(qiáng)和持續(xù)非常有利。29日03時預(yù)報逆溫層明顯抬高（圖5c），逆溫層底抬高至107 m，逆溫層頂抬高至827 m，至06時逆溫層底抬高至180 m，逆溫層頂抬高至1305 m，逆溫層的抬高表明機(jī)場凍霧將逐漸減弱抬高，對06時凍霧減弱有較好的指示作用，同時逆溫層抬高與烏魯木齊市氣象站的能見度變差的時間吻合較好。

圖4 方案9預(yù)報28日09時500 hPa高度場（a）和海平面氣壓場（b）

圖5 方案9預(yù)報機(jī)場上空風(fēng)溫時間剖面及溫度垂直廓線

由此可見，本方案中逆溫層的變化與實況的能見度的變化有很好的對應(yīng)關(guān)系，逆溫層的變化對凍霧生消、加強(qiáng)及減弱有明顯的參考價值。

4.3 溫濕條件

在干旱和半干旱地區(qū)，水汽條件在形成霧的過程中起著十分關(guān)鍵的作用[9]。從方案9預(yù)報的2 m對相對濕度上分析機(jī)場上空相對濕度變化情況（圖6a），在2 m溫度預(yù)報上，從28日08時開始，機(jī)場溫度迅速下降，導(dǎo)致相對濕度快速上升，2 m濕度維持超過80%的時段，與機(jī)場出霧時間基本一致，在17時后維持超過90%的濕度。從28日10—17時，機(jī)場上空相對濕度幾乎小于80%（圖6b），17時后相對濕度大于80%的厚度開始增加，29日07時達(dá)到了268 m，最厚時高度達(dá)372 m，相對濕度強(qiáng)度加強(qiáng)和厚度增加對霧的加強(qiáng)和維持是有利的。從濕度垂直結(jié)構(gòu)預(yù)報來看，10—17時霧的厚度比較小，17時后霧區(qū)厚度加大，但是由于沒有資料佐證，并不能確定霧區(qū)真正達(dá)到的高度，而烏魯木齊市區(qū)（較機(jī)場高約282 m）29日07時后能見度迅速下降至500 m或以下，能從側(cè)面反映濕度逸80%的區(qū)域確實與霧區(qū)有一定的對應(yīng)關(guān)系，可作為霧出現(xiàn)和霧厚度的判據(jù)之一。

4.4 混合條件

圖6 方案9預(yù)報機(jī)場28日06時—29日15時溫度、2 m濕度（a）及濕度隨高度的時間剖面（b）

圖7 風(fēng)場及機(jī)場上空風(fēng)、垂直速度的時空剖面（b圖中等值線為垂直速度，單位為100 m/s）

28日09時10 m風(fēng)場（圖7a）上，烏魯木齊機(jī)場附近有偏北風(fēng)和偏南風(fēng)的弱輻合區(qū)，風(fēng)速維持在1~3 m/s，弱風(fēng)場輻合區(qū)有利于污染物在烏魯木齊機(jī)場附近聚集，在其它條件具備情況下，凍霧易在機(jī)場附近生成且維持。在垂直方向上（圖7b），28日06時—29日06時，機(jī)場上空180 m下為弱風(fēng)場，風(fēng)速1耀2 m/s，垂直速度較弱，但大部分時間為弱的上升運動，這種弱上升運動有利于水汽和污染物在垂直方向上混合，使得霧區(qū)厚度加大，持續(xù)時間更長。

5 結(jié)論與討論

（1）此次個例中，通過對參數(shù)化方案進(jìn)行組合試驗，以2 m相對濕度逸80%來表示凍霧進(jìn)行分析，發(fā)現(xiàn)模式預(yù)報對近地層、陸面過程、邊界層等參數(shù)化方案組合較為敏感，微物理過程有微調(diào)作用，最好的方案為方案9，即微物理過程為WDM6、近地層方案為QNSE，陸面過程方案為Noah，邊界層方案QNSE。

（2）在合適的方案下，通過對此次凍霧過程的要素預(yù)報進(jìn)行分析發(fā)現(xiàn)，模式預(yù)報的環(huán)流形勢、層結(jié)條件、溫濕條件、混合條件能夠幫助預(yù)報員判別出此次凍霧的生消。針對本次凍霧過程，模式預(yù)報中機(jī)場上空穩(wěn)定層的變化，濕層結(jié)構(gòu)、風(fēng)場的水平和垂直結(jié)構(gòu)等，對凍霧的生消以及大霧過程的能見度變化有一定的指示意義。

（3）本文對凍霧天氣的預(yù)報分析說明，利用WRF模式對烏魯木齊機(jī)場凍霧進(jìn)行預(yù)報具有一定的可行性，通過對出霧條件進(jìn)行診斷，可以診斷出凍霧過程，為預(yù)報提供高時空分辨率的資料，但仍需要進(jìn)一步驗證，在驗證完成后，建立診斷預(yù)報模型，直接輸出相應(yīng)的產(chǎn)品，為將來機(jī)場凍霧預(yù)報提供參考資料。

參考文獻(xiàn)：

[1]王益栢，梅娜，范磊，等.WRF模式對2013年1月華北一次大霧的數(shù)值對比試驗[J].氣象，2014，40（12）：1522-1529.

[2]陸雪，高山紅，饒莉娟，等.春季黃海海霧宰砸云參數(shù)化方案敏感性研究[J].應(yīng)用氣象學(xué)報，2014，25（3）：312-320.

[3]黃政，袁成松，包云軒，等.基于不同參數(shù)化方案的高速公路大霧過程的數(shù)值模擬試驗[J].氣象，2016，42（8）：944-953.

[4]周梅，銀燕，王巍巍.2006年12月24—27日大范圍大霧過程數(shù)值模擬[J].應(yīng)用氣象學(xué)報，2008，19（5）：602-610.

[5]梁愛民，張慶紅，劉開宇，等.華北地區(qū)一次大霧過程的三維變分同化試驗[J].氣象學(xué)報，2007，65（5）：792-803.

[6]劉開宇，梁愛民，蔡軍，等.貴陽機(jī)場一次輻射大霧的數(shù)值模擬研究[J].貴州氣象，2007，31（3）：10-13.

[7]王金蘭，壽紹文，劉澤軍，等.河南省一次大霧的數(shù)值模擬及生消機(jī)制分析[J].氣象與環(huán)境科學(xué).2008，31（1）：39-44.

[8]《新疆短期天氣預(yù)報指導(dǎo)手冊》編寫組.新疆短期天氣預(yù)報指導(dǎo)手冊[M].烏魯木齊：新疆人民出版社，1986.

[9]王春紅.烏魯木齊國際機(jī)場2009年秋冬初霧天氣的數(shù)值模擬分析[J].沙漠與綠洲氣象，2011，5（5）：15-19.

[10]王楠，楊洪儒，梁艷.近37 a烏魯木齊機(jī)場霧的不均勻性特征分析[J].沙漠與綠洲氣象，2015，9（4）：43-49.

[11]萬小雁，包云軒，嚴(yán)明良，等.不同陸面方案對滬寧高速公路團(tuán)霧的模擬[J].氣象科學(xué)，2010，30（4）.487-494.

[12]彭雙姿，劉從省，屈右銘，等.一次大范圍輻射霧天氣過程的觀測和數(shù)值模擬分析[J].氣象，2012，38（6）：679-687.

[13]包云軒，丁秋冀，袁成松，等.滬寧高速公路一次復(fù)雜性大霧過程的數(shù)值模擬試驗[J].大氣科學(xué)，2013，37（1）：124-136.

[14]榮昕.基于WRF-CHEM模式的連續(xù)霧霾過程數(shù)值模擬及其能見度參數(shù)化[D].南京：南京信息工程大學(xué)，2013，25-26.

[15]石春娥，李耀孫，楊軍，等.MM5和WRF對中國東部地區(qū)冬季邊界層模擬效果比較[J].高原氣象，2015，34（2）：389-399.

[16]周旺，劉昊，樊斌.重慶機(jī)場輻射霧消散機(jī)制[J].氣象，2008，34（T1）：231-237.