陳 健,周后福,周文鱗,方 茸,趙 倩
(1.合肥市氣象局,安徽 合肥230041;2.安徽省氣象科學(xué)研究所,安徽 合肥230031;3.江蘇理工學(xué)院化學(xué)與環(huán)境工程學(xué)院,江蘇 常州213001)
能見度是反映大氣透明度的一個(gè)指標(biāo), 大霧是影響能見度的主要天氣現(xiàn)象,根據(jù)中國氣象局《地面氣象觀測(cè)規(guī)范》[1]規(guī)定:大霧是指大量微小水滴浮游在空中,使水平能見度<1.0 km 的天氣現(xiàn)象。 大霧對(duì)交通安全影響非常大, 導(dǎo)致地方政府對(duì)大霧預(yù)報(bào)的氣象服務(wù)需求不斷提升。 經(jīng)過前期對(duì)合肥市大霧時(shí)空特征統(tǒng)計(jì)發(fā)現(xiàn), 大霧是每年秋冬季節(jié)合肥市最常見的災(zāi)害性天氣。
大霧預(yù)報(bào)常用的方法主要有天氣學(xué)和數(shù)值模式產(chǎn)品釋用方法[2]。 早期以主觀判斷為主,預(yù)報(bào)準(zhǔn)確率建立在生消機(jī)制及天氣氣候形勢(shì)的充分認(rèn)識(shí)之上,查找與大霧關(guān)聯(lián)的宏觀因子, 利用統(tǒng)計(jì)方法開展預(yù)報(bào)[3-4]。 后來數(shù)值模式模擬是大霧的主要研究手段之一,體現(xiàn)出模式的高時(shí)空分辨率的優(yōu)勢(shì),近年來則嘗試使用集合預(yù)報(bào)的方式進(jìn)行霧的預(yù)報(bào)[5-6]。 業(yè)務(wù)數(shù)值預(yù)報(bào)模式直接輸出的產(chǎn)品中沒有大霧的預(yù)報(bào), 有些中尺度數(shù)值天氣模式依托云水、云冰、霰等水凝物的消光性質(zhì)來反算能見度, 類似于對(duì)大霧進(jìn)行預(yù)報(bào)。 這種方式的缺點(diǎn)是沒考慮氣溶膠變化的影響,對(duì)非過飽和條件下的低能見度天氣幾乎無預(yù)報(bào)能力[7-8],故現(xiàn)階段大霧預(yù)報(bào)大多數(shù)還是采用數(shù)值模式的釋用技術(shù)[8-11]。
目前對(duì)影響合肥市城市交通的大霧客觀預(yù)報(bào)方法和數(shù)值預(yù)報(bào)產(chǎn)品釋用技術(shù)研究相對(duì)較少, 預(yù)報(bào)員大多通過天氣分析后的主觀方法來判別大霧, 大霧預(yù)報(bào)還沒有特定的客觀預(yù)報(bào)產(chǎn)品。 本文首先利用主成分分析方法對(duì)合肥市縣觀測(cè)站點(diǎn)進(jìn)行大霧區(qū)域劃分,再基于PP 法建立預(yù)報(bào)模型,形成合肥市大霧預(yù)報(bào)業(yè)務(wù)系統(tǒng),實(shí)時(shí)輸出大霧格點(diǎn)產(chǎn)品。預(yù)報(bào)系統(tǒng)的空間格距約在13 km, 于2015 年1 月投入業(yè)務(wù)運(yùn)行,解決了基層臺(tái)站大霧預(yù)報(bào)的實(shí)際業(yè)務(wù)需求。
區(qū)域劃分:利用主成分分析方法[12-13],將m 個(gè)變量組合成P 個(gè)公因子, 并研究公因子的荷載陣,根據(jù)m 個(gè)點(diǎn)的空間位置了解m 個(gè)變量的相互關(guān)系或?qū) 個(gè)變量進(jìn)行分類;其主要思路是利用從因子群中提取共性因子的方法,繪制因子荷載點(diǎn)聚圖,根據(jù)點(diǎn)聚圖上觀測(cè)站點(diǎn)的分布情況進(jìn)行區(qū)域劃分。
預(yù)報(bào)模型思路:PP 法[14-15]是根據(jù)預(yù)報(bào)量和預(yù)報(bào)因子的同時(shí)性(或近于同時(shí)性)的加權(quán)組合,利用歷史觀測(cè)資料來確定局地天氣要素, 其推導(dǎo)方程的函數(shù)關(guān)系為:Y0=F(X0),其中X0、Y0分別表示T0時(shí)刻的因子向量和預(yù)報(bào)量;為了用導(dǎo)出的方程制作預(yù)報(bào),用模擬的數(shù)值預(yù)報(bào)模式的輸出結(jié)果XT代入方程而求得YT,即YT=F(XT)。 本文基于PP 法統(tǒng)計(jì)模式輸出產(chǎn)品與大霧之間的相關(guān)性, 找出有較好指示意義的預(yù)報(bào)因子,建立預(yù)報(bào)模型,再利用Fortran 編制計(jì)算代碼,基于Visual BASIC 開發(fā)大霧天氣精細(xì)化預(yù)報(bào)系統(tǒng),形成合肥市大霧預(yù)報(bào)格點(diǎn)產(chǎn)品??傮w技術(shù)流程見圖1。
圖1 總體技術(shù)流程
合肥市轄范圍涉及到主城區(qū)和5 個(gè)縣市, 所轄區(qū)域南北距離相對(duì)較長(zhǎng), 大霧發(fā)生的區(qū)域特點(diǎn)有所不同,因此有必要進(jìn)行區(qū)域劃分[16-18]。 ECMWF 細(xì)網(wǎng)格資料于2011 年下半年開始下發(fā)到省級(jí)氣象部門,故本文選取2012—2014 年ECMWF 細(xì)網(wǎng)格資料;再根據(jù)合肥主城區(qū)及5 個(gè)縣市65 個(gè)交通氣象站3 a日最低能見度資料,得到87 個(gè)大霧觀測(cè)樣本。 采用大霧等級(jí)分類法進(jìn)行大霧強(qiáng)度的劃分, 若某個(gè)區(qū)域出現(xiàn)水平能見度V<500 m 的站點(diǎn),則其大霧強(qiáng)度取值為2;若500 m≤V<1000 m,則取值為1;若V≥1000 m,則取值為0。
合肥主城區(qū)及5 個(gè)縣市大霧觀測(cè)樣本組成矩陣mXn,其中m=6, n=87,對(duì)其進(jìn)行因子處理。 表1 說明前5 個(gè)特征值的累積百分率, 可以看出取前5 個(gè)特征值時(shí)精確度達(dá)93.35%。
表1 特征值和累積百分率(前5 個(gè))
因子荷載陣6A5中第i 列元素的平方和為λi,它表示第i 個(gè)公因子對(duì)6 個(gè)原變量所提供的方差之和。 而6A5中第i 行元素的平方和稱為共同度,它表示這5 個(gè)公因子對(duì)原變量Xi的總方差作的貢獻(xiàn)。取2 個(gè)公因子, 則可以用6A5中的左邊兩列的荷載值點(diǎn)圖,這兩列分別表示6 個(gè)因子在第一、第二個(gè)公因子上的投影。
因子荷載點(diǎn)聚圖:將第一、第二個(gè)公因子取為坐標(biāo)軸F1、F2。 將荷載值(Ai1,Ai2),(i=1,2,3,4,5,6)同時(shí)擴(kuò)大10 倍。則變量可分成2 類:X1(6.99,-2.83),X2(7.05,-1.25),X3(5.53,-4.85),X5(7.83, 2.59),X6(7.26,1.78)成一類;X4(2.68,8.31)成一類。這樣可將合肥市縣大霧劃分為2 個(gè)區(qū),分別稱之為中南區(qū)、北區(qū)(表2)。
表2 合肥市大霧區(qū)劃
在因子荷載陣6A5中, 用5 個(gè)主因子計(jì)算6 個(gè)變量的共同度。 可以看出,對(duì)X1,X2,X3,X5,X6這一類,=0.9997 為最大, 表示用X2代表X1,X2,X3,X5,X6比較好, 也就是中南區(qū)中, 肥東站代表性最好;北區(qū)中長(zhǎng)豐=0.993 7。
一般04—08 時(shí)(北京時(shí),下同)為安徽省出霧的最佳時(shí)段[16-18],故將大霧日05 時(shí)ECMWF 細(xì)網(wǎng)格物理量預(yù)報(bào)資料用距離插值法插值到全市和分片的代表區(qū)域站上, 選取和大霧密切相關(guān)的氣象因子與分區(qū)時(shí)采用的大霧等級(jí)分類進(jìn)行相關(guān)分析,確定8、8、7 個(gè)相關(guān)系數(shù)絕對(duì)值≥0.3 為高影響因子,它們均通過α=0.05 顯著性檢驗(yàn)(表3)。 所選取的高影響因子以大氣低層的物理量為主,包括近地層逆溫條件、濕度條件和風(fēng)速大小,上干下濕配置以及云量等[19-20]。
表3 對(duì)全市和分片大霧有高影響的氣象因子及相關(guān)系數(shù)
把全市和分片所選取的高影響因子, 利用大霧等級(jí)分類和逐步回歸分別建立全市、 中南區(qū)和北區(qū)24 h 大霧預(yù)報(bào)模型,它們都通過α=0.05 顯著性水平檢驗(yàn)。 如對(duì)全市大霧有8 個(gè)高影響的氣象因子,在0.05 顯著性水平檢驗(yàn)標(biāo)準(zhǔn)約束時(shí)進(jìn)入模型的氣象因子有4 個(gè), 分別是T3:925 hPa 和2 m 溫度差,R2:925 hPa 相對(duì)濕度,R6:1000 hPa 和850 hPa 相對(duì)濕度差,V1:10 m 風(fēng)速。 故24 h 大霧預(yù)報(bào)方程為Y=2.301+0.028T3-0.376R2+0.313R6-0.203V1。
當(dāng)水平能見度V<500 m 時(shí),即會(huì)對(duì)交通造成影響, 氣象部門就需要及時(shí)向公眾和相關(guān)部門發(fā)布大霧預(yù)警信號(hào), 因此本文主要針對(duì)水平能見度V<500 m的大霧天氣開展預(yù)報(bào)[21]。利用2012—2014 年大霧日05 時(shí)有高影響的氣象因子值代入模型,其結(jié)果與實(shí)測(cè)值進(jìn)行擬合, 確定全市和分片出現(xiàn)大霧時(shí)回歸結(jié)果的閾值均為1.5,擬合準(zhǔn)確率分別為100%、93%和92%; 出現(xiàn)能見度V<1000 m 的霧時(shí)回歸結(jié)果的閾值均為1.2,擬合準(zhǔn)確率分別為98%、92%和90%。
安徽省的霧形成主要以輻射冷卻為主因, 但有時(shí)也有平流的作用,就合肥地區(qū)而言,輻射霧和平流輻射霧占九成以上[16-18]。 吳洪等[22]認(rèn)為若有大霧生成,須在前期(前一日20:00 或當(dāng)日08:00)滿足大氣層結(jié)穩(wěn)定、850 hPa 以下存在逆溫層或等溫層、地面較濕和地面風(fēng)速<4 m/s 等條件; 結(jié)合前期對(duì)大量大霧個(gè)例的分析及高影響因子的選取表明: 近地層濕度、風(fēng)速和逆溫條件可以作為大霧的判定條件,故通過對(duì)87 個(gè)大霧觀測(cè)樣本與代表區(qū)域站大霧日05時(shí)1000 hPa 相對(duì)濕度、10 m 風(fēng)速及925 hPa 和2 m溫差進(jìn)行對(duì)比分析,得到1000 hPa 相對(duì)濕度<50%、10 m 風(fēng)速≥4 m/s 及925 hPa 和2 m 溫差<-1 ℃是不考慮出大霧的消空指標(biāo),其中全市、中南區(qū)和北區(qū)的漏報(bào)率僅為1.8%、1.8%和5.3%;霧的消空指標(biāo)與大霧相同,漏報(bào)率分別為3.5%、5.8%和9.3%。 經(jīng)過以上的篩選,基本上可以提前將無霧日濾掉,借以提高預(yù)報(bào)準(zhǔn)確率。 所選取的消空指標(biāo)雖然濾掉了大部分無霧日,但由于大霧形成的原因比較復(fù)雜,在沒有考慮到暖濕平流、地面有效輻射、空氣質(zhì)量狀況以及局地地理?xiàng)l件等因素的情況下, 可能會(huì)造成漏報(bào)和空?qǐng)?bào)。
大霧預(yù)報(bào)系統(tǒng)(圖2) 主要通過把全市范圍內(nèi)(116.5°~118°E,31°~32.75°N)、中南區(qū)(116.5°~118°E,31°~32.25°N)和北區(qū)(116.5°~118°E,32°~32.75°N)ECMWF 細(xì)網(wǎng)格物理量格點(diǎn)輸出產(chǎn)品與預(yù)報(bào)模型對(duì)接,在考慮到消空指標(biāo)和回歸結(jié)果的閾值情況下,每日定時(shí)輸出合肥市大霧預(yù)報(bào)格點(diǎn)產(chǎn)品。若能見度V<500 m,則格點(diǎn)值取2;若500 m≤V<1000 m,則取1;若V≥1000 m,則取0,方便預(yù)報(bào)員調(diào)用。
圖2 大霧天氣預(yù)報(bào)系統(tǒng)
分為2015—2017 年冬半年大霧預(yù)報(bào)統(tǒng)計(jì)檢驗(yàn)和2 類典型大霧天氣過程檢驗(yàn)。
針對(duì)2015—2017 年冬半年的日最低能見度,以500 m 為分界線對(duì)全市和分片大霧預(yù)報(bào)進(jìn)行檢驗(yàn),為期18 個(gè)月547 d。 大霧的預(yù)報(bào)效果檢驗(yàn)采用預(yù)報(bào)值與日最低能見度觀測(cè)資料的比較[7],選取距離代表區(qū)域站最近的格點(diǎn)預(yù)報(bào)值與該區(qū)域交通氣象站日最低能見度資料進(jìn)行對(duì)比檢驗(yàn)。TS 技巧評(píng)分結(jié)果見表4。 對(duì)于<500 m 的日最低能見度預(yù)報(bào),全市和分片數(shù)值模式產(chǎn)品釋用方法的技巧評(píng)分分別為92%、83%和76%,參照WRF 模式每日輸出的輻射霧能見度預(yù)報(bào)技巧評(píng)分68%、59%和51%,表明數(shù)值模式產(chǎn)品釋用方法在有無大霧預(yù)報(bào)技巧方面明顯占優(yōu),技巧評(píng)分大幅提升。
表4 2015—2017 年數(shù)值模式產(chǎn)品釋用方法、WRF 模式預(yù)報(bào)TS 評(píng)分
本文根據(jù)不同環(huán)流形勢(shì)將影響安徽的大霧天氣主要分為兩類, 分別為Ⅰ類槽后西北氣流型和Ⅱ類西風(fēng)氣流型[16,23]。 選取2017 年11 月15 日和2018年9 月22 日分屬Ⅰ類、Ⅱ類環(huán)流形勢(shì)下典型大霧天氣過程進(jìn)行預(yù)報(bào)效果檢驗(yàn), 可以更直觀地看出數(shù)值模式產(chǎn)品釋用的預(yù)報(bào)方法效果更好。
5.2.1 槽后西北氣流型(I 類)
2017 年11 月14 日20 時(shí)500 hPa 亞歐大陸為兩槽一脊,低槽快速經(jīng)過安徽,安徽大部由槽前西南氣流影響轉(zhuǎn)為槽后西北氣流控制,天氣形勢(shì)由不穩(wěn)定轉(zhuǎn)為穩(wěn)定,天空狀況轉(zhuǎn)好;同時(shí)伴隨低槽過境,降水過程結(jié)束初期底層相對(duì)濕度仍然較大(圖3a),符合Ⅰ類環(huán)流形勢(shì)特征。此次大霧過程發(fā)生時(shí)段合肥市晴天少云,輻射降溫明顯,且近地層相對(duì)濕度大、風(fēng)速較小、存在逆溫層,滿足輻射霧的形成條件。
11 月15 日早晨安徽省沿淮到沿江一帶出現(xiàn)分布不均的大霧天氣, 合肥市中北部大部分地區(qū)最低能見度不足200 m, 而南部和東南部部分地區(qū)則>1000 m(圖4a),由于全省受槽后西北氣流控制,天氣較好,氣溫回升迅速,9:00 后能見度已經(jīng)轉(zhuǎn)好。11 月14 日數(shù)值模式產(chǎn)品釋用的預(yù)報(bào)方法顯示(圖5a),除西南部局部地區(qū)以外,合肥市其他地區(qū)均有大霧; 和實(shí)況對(duì)比發(fā)現(xiàn)僅東南部和南部部分地區(qū)預(yù)報(bào)偏強(qiáng)外,其他地區(qū)與實(shí)況一致。模式預(yù)報(bào)15 日05時(shí)全市925 hPa 和2 m 溫度差為0~2 ℃,1000 hPa相對(duì)濕度為70%~80%,1000 hPa 和850 hPa 相對(duì)濕度差為30%~40%,10 m 風(fēng)速2 m/s 左右, 無低云,進(jìn)入預(yù)報(bào)模型的近地層逆溫條件、 濕度條件和風(fēng)速大小,上干下濕以及云量等氣象因子基本符合Ⅰ類大霧中的物理量配置特征。WRF 模式輸出的能見度預(yù)報(bào)依托對(duì)影響能見度的各因子(云水、云冰、霰和氣溶膠等)消光特征的反演,模式預(yù)報(bào)(圖5b)合肥市有零散的大霧天氣,大部分區(qū)域能見度>1000 m,和實(shí)況對(duì)比發(fā)現(xiàn)中北部大片區(qū)域被漏報(bào)。
5.2.2 西風(fēng)氣流型(II 類)
2018 年9 月21 日20 時(shí)500 hPa 亞歐大陸亦為兩槽一脊, 受東北冷渦影響, 中高緯為寬廣的槽區(qū),安徽位于槽區(qū)底部西風(fēng)氣流中(圖3b),符合Ⅱ類環(huán)流形勢(shì)特征。 此次大霧過程滿足平流輻射霧的形成條件,21 日白天受切變線影響, 江淮之間南部到江南有小到中雨,夜里切變線南壓,江南部分地區(qū)仍有小雨, 而江北轉(zhuǎn)為多云天氣為主, 天空逐漸打開,有一定的輻射降溫,而且降水過程結(jié)束初期底層相對(duì)濕度較大; 從22 日08 時(shí)阜陽站和南京站溫度對(duì)數(shù)壓力圖來看,1000~925 hPa 有明顯逆溫存在,且風(fēng)向順轉(zhuǎn),有弱暖平流。
圖3 2017 年11 月14 日20 時(shí)(a)和2018 年9 月21 日20 時(shí)(b)500 hPa 環(huán)流形勢(shì)
圖4 2017 年11 月15 日(a)、2018 年9 月22 日(b)合肥市日最低能見度(km)實(shí)況
圖5 數(shù)值模式產(chǎn)品釋用方法(a、c)和WRF 模式(b、d)分別輸出的2017 年11 月15 日05 時(shí)(a、b)、2018 年9 月22 日05 時(shí)(c、d)能見度預(yù)報(bào)結(jié)果
9 月22 日早晨安徽省江淮之間部分地區(qū)、沿淮淮北及江南局部地區(qū)出現(xiàn)大霧天氣,合肥市東北部最低能見度普遍不足500 m, 部分地區(qū)<200 m,其他地區(qū)大部分在1000 m 以下, 局部不足500 m(圖4b),全省位于槽區(qū)底部西風(fēng)氣流中,以多云到晴天氣為主,輻射降溫有利于大霧的形成,但同時(shí)也存在暖濕平流作用, 因此直到11:00 后能見度才明顯轉(zhuǎn)好。 9 月21 日數(shù)值模式產(chǎn)品釋用的預(yù)報(bào)方法表明(圖5c) 合肥市中北部大部分地區(qū)能見度<500 m,而南部大部分地區(qū)>1000 m;和實(shí)況對(duì)比發(fā)現(xiàn),除市區(qū)及西北部預(yù)報(bào)偏強(qiáng)外,其他地區(qū)與實(shí)況一致; 模式預(yù)報(bào)22 日05 時(shí)全市925 hPa 和2 m溫度差為0~1 ℃,1000 hPa 相對(duì)濕度為75%~85%,1000 hPa 和850 hPa 相對(duì)濕度差為20%~30%,10 m風(fēng)速為2~3 m/s,低云量20%左右,進(jìn)入預(yù)報(bào)模型的氣象因子沒有考慮到暖濕平流的作用,其他基本符合Ⅱ類大霧中的物理量配置特征。 WRF 模式預(yù)報(bào)(圖5d)合肥市能見度普遍在1000 m 以上,大霧過程完全被漏報(bào)。
本文利用主成分分析方法對(duì)合肥市縣觀測(cè)站點(diǎn)進(jìn)行大霧區(qū)域劃分,再基于PP 法建立預(yù)報(bào)模型,形成合肥市大霧預(yù)報(bào)業(yè)務(wù)系統(tǒng),結(jié)論如下:
(1)合肥市縣大霧劃分可為兩個(gè)區(qū),分別稱之為中南區(qū)、北區(qū),其中肥東、長(zhǎng)豐可以作為它們的代表區(qū)域站。
(2)基于PP 法統(tǒng)計(jì)ECMWF 模式輸出資料與大霧之間的相關(guān)性, 全市和分片分別確定8、8、7 個(gè)與大霧密切相關(guān)的高影響因子, 再利用等級(jí)分類和逐步回歸建立24 h 大霧預(yù)報(bào)模型,其中全市和分片出現(xiàn)大霧時(shí)回歸結(jié)果的判定閾值為1.5,出現(xiàn)霧時(shí)為1.2。 同時(shí)在1000 hPa 相對(duì)濕度<50%、10 m 風(fēng)速≥4 m/s 及925 hPa 和2 m 溫差<-1 ℃等消空指標(biāo)選定的情況下,基本上可以提前將無霧日濾掉,以提高預(yù)報(bào)準(zhǔn)確率。
(3)研發(fā)的大霧天氣精細(xì)化預(yù)報(bào)系統(tǒng)把ECMWF細(xì)網(wǎng)格物理量格點(diǎn)輸出產(chǎn)品與全市和分片大霧預(yù)報(bào)模型對(duì)接,每日定時(shí)輸出合肥市大霧預(yù)報(bào)格點(diǎn)產(chǎn)品。前期業(yè)務(wù)化運(yùn)行和預(yù)報(bào)效果檢驗(yàn)表明: 數(shù)值模式產(chǎn)品釋用方法在有無大霧預(yù)報(bào)技巧方面較WRF 模式明顯占優(yōu),技巧評(píng)分大幅提升,而兩類典型大霧天氣過程預(yù)報(bào)效果檢驗(yàn)則可以更直觀地看出數(shù)值模式產(chǎn)品釋用的預(yù)報(bào)方法效果要更好。
大霧具有突發(fā)性、局地性、頻次少和持續(xù)時(shí)間短等特點(diǎn), 在實(shí)際業(yè)務(wù)中強(qiáng)度和落區(qū)預(yù)報(bào)方面難度較大。 雖然數(shù)值模式產(chǎn)品釋用方法已經(jīng)應(yīng)用在合肥市大霧預(yù)報(bào)業(yè)務(wù)中, 但模式輸出產(chǎn)品相對(duì)于實(shí)況存在誤差,釋用方法無法擺脫統(tǒng)計(jì)本身的缺陷,缺乏對(duì)大霧過程物理規(guī)律的描述, 預(yù)報(bào)模型的建立需要長(zhǎng)期的、穩(wěn)定的歷史數(shù)據(jù),很多原因都可能影響大霧預(yù)報(bào)質(zhì)量。 因此, 隨著模式分辨率和預(yù)報(bào)性能的不斷提高,能夠挑選出反映大霧微物理過程、熱力和動(dòng)力過程的預(yù)報(bào)因子,且隨著模式參數(shù)化方案、預(yù)報(bào)模型的不斷改進(jìn), 大霧的預(yù)報(bào)預(yù)警能力將會(huì)上一個(gè)新的臺(tái)階。