楊秀莊，彭 芳，童碧慶

(貴州省氣象臺，貴州 貴陽 550002)

1 引言

我國是典型的季風(fēng)氣候區(qū)，夏半年受東亞季風(fēng)環(huán)流影響，強(qiáng)對流天氣頻繁，雷電是強(qiáng)對流天氣最明顯的特征。雷暴一般具有突發(fā)性強(qiáng)、范圍小、局地災(zāi)害嚴(yán)重等中小尺度特征，因此對雷暴的預(yù)報難度較大。隨著大氣探測設(shè)備的更新，如閃電定位儀、大氣電場儀、雷達(dá)衛(wèi)星等的廣泛應(yīng)用，雷電特別是強(qiáng)雷電的研究取得了很大的進(jìn)展。張義軍［1］、張延龍［2］等對雷電的電荷結(jié)構(gòu)及其成因做了很多研究，得出一些很有用的結(jié)論。郝瑩［3］等做了關(guān)于判別分析法和指標(biāo)疊加法應(yīng)用于雷暴潛勢預(yù)報方面的研究;鄭傳烈［4］等將獨立參數(shù)條件概率法應(yīng)用于雷暴的預(yù)報;AgostinoManzato 等［5］利用探空資料以及閃電定位資料對神經(jīng)網(wǎng)絡(luò)進(jìn)行了訓(xùn)練，并利用神經(jīng)網(wǎng)絡(luò)法對雷暴的發(fā)生以及閃電的密度進(jìn)行了預(yù)報;金米娜［6］等利用衛(wèi)星資料分析了中尺度對流系統(tǒng)(MCS)的雷電演變特征。

貴州地處青藏高原東南部，橫斷山脈的下坡段，全省多山，對流天氣頻繁，雷電天氣較多。楊秀莊等［7］通過分析貴州雷電時空分布特征得到貴州雷電主要發(fā)生時間為4—10月份，呈不均勻分布，西部、南部偏多，北部、東部偏少，特別是省的西南部雷電天氣頻繁。因此，加強(qiáng)對貴州雷電天氣的分析與研究，提高雷電天氣的監(jiān)測預(yù)報能力，對于防災(zāi)減災(zāi)有十分重要的意義。本文通過中尺度模式輸出的時空分辨率較高的物理量場資料，開展對貴州雷電潛勢預(yù)報研究。

2 研究資料和MOS 方法介紹

選取2011年5—8月份ADTA 閃電定位儀地閃資料，范圍為貴州區(qū)域(22～30°N，103～110°E)。利用ADTA 地閃資料通過編程實現(xiàn)經(jīng)緯網(wǎng)格(0.5°×0.5°)插值，獲取經(jīng)緯網(wǎng)格上地閃頻數(shù)。物理量場資料為以T639 數(shù)值模式為背景場的中尺度模式WRF 輸出的物理量場，時間分辨率為3 h，空間分辨率為5 km。

MOS 預(yù)報方法是統(tǒng)計方法的一種。選取一定數(shù)量的樣本，以數(shù)值模式輸出的數(shù)據(jù)作為預(yù)報因子，以同時刻實際觀測的數(shù)據(jù)為預(yù)報量，二者建立統(tǒng)計關(guān)系。其推導(dǎo)方程的函數(shù)關(guān)系式如下:

3 WRF 模式物理量選取及物理意義

以T639 數(shù)值模式為背景場的中尺度模式WRF輸出的物理量很多，包括高度場、風(fēng)場、溫度場等等，在分析雷暴發(fā)生基本條件基礎(chǔ)上，選取10 個能有效反映雷暴發(fā)生發(fā)展的物理量參數(shù)。

根據(jù)雷暴發(fā)生的基本條件，要出現(xiàn)明顯的雷暴天氣，大氣必須為層結(jié)不穩(wěn)定，通常來說，主要表現(xiàn)為大氣的熱力不穩(wěn)定。結(jié)合天氣學(xué)知識，選取了常用的K 指數(shù)、對流有效不穩(wěn)定能量(CPAE)、對流抑制能量(CIN)、地面總溫度(TT)，這4 項物理量參數(shù)主要反映大氣聚集的有效能量，選取沙氏指數(shù)(SI)、850 hPa 與300 hPa 溫度之差(T850－T300)，主要反映大氣的不穩(wěn)定指標(biāo)。同時，雷暴的發(fā)生必須有合適的水汽條件，所以選取了物理量850 hPa、700 hPa、500 hPa 溫度露點差之和(TTd)。不穩(wěn)定條件和水汽條件都滿足的情況下，不一定會有雷暴的發(fā)生，必須有一定的觸發(fā)機(jī)制，才有雷暴天氣的發(fā)生，反映雷暴觸發(fā)機(jī)制的動力條件選取了850 hPa 渦度(vort850)、850 hPa 散度(div850)、300 hPa 散度(div300)。最后，選取了850 hPa 與500 hPa 假相當(dāng)位溫之差(θse850－θse500)、850 hPa 水汽通量散度(fluxdiv850)等綜合指標(biāo)來加強(qiáng)對雷暴發(fā)生發(fā)展的反映，更好的開展雷暴潛勢預(yù)報研究。

4 物理量參數(shù)和雷電相關(guān)性研究及預(yù)報因子的篩選

上述選取的WRF 模式輸出的物理量參數(shù)是結(jié)合天氣學(xué)知識人為選取的，在貴州地區(qū)雷電特別是強(qiáng)雷電的發(fā)生發(fā)展是否與這些所有的物理量有關(guān)還有待科學(xué)的檢驗。收集整理5—8月份ADTA 閃電定位儀的地閃資料和WRF 模式輸出的物理量參數(shù)資料，選取了符合條件的188 個樣本。以經(jīng)緯網(wǎng)格內(nèi)的閃電頻數(shù)作為雷暴強(qiáng)度來表示，WRF 模式輸出物理量場選取上一日20 時起報的18 h 預(yù)報產(chǎn)品和當(dāng)日08 時起報的18 h 預(yù)報產(chǎn)品。二者之間建立相關(guān)性統(tǒng)計關(guān)系。其關(guān)系式如下:

其中x 為經(jīng)緯網(wǎng)格內(nèi)的閃電頻數(shù)，y 為WRF 模式輸出的18 h 預(yù)報產(chǎn)品。通過計算，上述選取的10個物理量與雷電的相關(guān)系數(shù)如表1。

表1 中NaN 表示數(shù)據(jù)缺失，本文不予討論，其余數(shù)據(jù)為物理量產(chǎn)品與雷電的相關(guān)系數(shù)。由于所選取樣本數(shù)較少，所有必須對相關(guān)系數(shù)數(shù)據(jù)進(jìn)行檢驗。在顯著水平α =0.05 情況下，樣本容量為188個時臨界相關(guān)系數(shù)為0.138 4，即如果相關(guān)系數(shù)在0.138 4 以上，則呈顯著相關(guān)。通過分析，選取了呈顯著相關(guān)的6 個物理量產(chǎn)品。包括K 指數(shù)、CAPE(對流不穩(wěn)定能量)、SI(沙氏指數(shù))、TT(地面總溫度)、θse850－θse500(850 hPa 與500 hPa 假相當(dāng)位溫之差)、TTd(850 hPa、700 hPa、500 hPa 溫度露點差之和)。其中SI、TTd 和雷電為明顯的負(fù)相關(guān)。

表1 WRF 模式輸出的物理量18 h預(yù)報產(chǎn)品與雷電相關(guān)系數(shù)

5 利用判別分析法進(jìn)行雷電潛勢預(yù)報

利用WRF 模式輸出的上一日20 時起報的18 h預(yù)報產(chǎn)品和當(dāng)日08 時起報的18 h 預(yù)報產(chǎn)品，選取呈顯著相關(guān)6 個物理量作為預(yù)報因子，以ADTA 閃電定位儀監(jiān)測到的單位經(jīng)緯網(wǎng)格的閃電頻數(shù)為預(yù)報對象，在樣本容量為188的情況下用判別分析法使二者之間建立統(tǒng)計關(guān)系，建立判別分析方程。計算得到的雷電潛勢預(yù)報方程如下:

其中X1～X6為預(yù)報因子，Y 為預(yù)報量。X1為K指數(shù)，X2為CAPE(對流不穩(wěn)定能量)，X3為SI(沙氏指數(shù))，X4為TT(地面總溫度)，X5為θse850－θse500，x6為TTd。

YC=－215.1 為預(yù)報量閾值。

當(dāng)計算的Y ＞YC=－215.1 時，則無雷電。

當(dāng)計算的Y ＜YC=－215.1 時，則有雷電。

在實際預(yù)報中，利用WRF 模式輸出的上述6 個物理量18 h 預(yù)報產(chǎn)品，代入上述雷電潛勢預(yù)報方程，根據(jù)計算出的預(yù)報量Y 值與YC對比，就可以判定未來12 h(當(dāng)日08—20 時或20 時—次日08 時)單位經(jīng)緯網(wǎng)格有無閃電，從而建立12 h 雷電潛勢預(yù)報。

圖1 是以9月3日20 時為起始場的WRF 模式輸出的18 h 預(yù)報產(chǎn)品，利用上述的雷電潛勢預(yù)報方程計算得到的9月4日08—20 時雷電潛勢預(yù)報產(chǎn)品。產(chǎn)品以Micaps 第4 類格點數(shù)據(jù)顯示，100 表示有雷電，0 表示無雷電，999 表示數(shù)據(jù)缺失。圖中可以看到，除東北部地區(qū)外，其余大部分地區(qū)都有明顯的雷電潛勢，西部邊緣地區(qū)數(shù)據(jù)缺失。由于西部邊緣地區(qū)海拔較高，某些物理量參數(shù)無法計算出來，如850 hPa 與500 hPa 假相當(dāng)位溫之差，因此得到的雷電潛勢預(yù)報產(chǎn)品數(shù)據(jù)缺失。圖2 是9月4日08 時—20 時的ADTA 閃電定位儀閃電監(jiān)測情況，圖中顯示，貴州北部地區(qū)雷電少，其余大部分地區(qū)都有雷電，其中省的西南部雷電比較密集，雷電頻數(shù)高，雷電強(qiáng)。對比分析圖1、圖2 可以看到，利用上述雷電潛勢預(yù)報方程得到的結(jié)果和實況基本吻合，但省北部地區(qū)有一定程度的空報，可能的原因是其受數(shù)值預(yù)報影響較大，WRF 模式當(dāng)日對北部的整體預(yù)報效果較差，且由于雷電潛勢預(yù)報方程本身的誤差以及閃電離散化程度較高的原因，造成了一定程度空報。但總的來說，本次個例雷電潛勢預(yù)報產(chǎn)品在實際的預(yù)報業(yè)務(wù)中具有一定的參考作用。通過選取其他幾次出現(xiàn)雷電天氣明顯的個例，雷電潛勢預(yù)報產(chǎn)品大趨勢預(yù)報較為準(zhǔn)確，但局部地區(qū)有一定程度空報或漏報。由于閃電數(shù)據(jù)離散程度較高，且ADTA 閃電定位儀探測的數(shù)據(jù)本身誤差較大，對閃電的預(yù)報預(yù)測一直是雷電預(yù)報的難點。本文中雷電潛勢預(yù)報產(chǎn)品對雷電有不同程度的空報或漏報，這也屬正?，F(xiàn)象，希望在今后的工作中對其進(jìn)一步改進(jìn)和提高，使其在實際的業(yè)務(wù)應(yīng)用中能發(fā)揮較大作用。

6 小結(jié)

①本文采用ADTA 閃電定位儀資料，通過編程實現(xiàn)獲取經(jīng)緯網(wǎng)格內(nèi)的閃電頻數(shù)，采用WRF 模式輸出的高分辨率的物理場資料，二者之間建立統(tǒng)計關(guān)系，利用MOS 方法開展雷電潛勢預(yù)報研究。

②通過研究得到對貴州雷電潛勢有明顯影響的物理量預(yù)報因子有K 指數(shù)、對流不穩(wěn)定能量、沙氏指數(shù)、地面總溫度、850 hPa 與500 hPa 假相當(dāng)位溫之差、850 hPa、700 hPa、500 hPa 溫度露點差之和，利用這些因子采用判別分析法建立了貴州地區(qū)雷電潛勢預(yù)報方程，通過檢驗，此雷電潛勢預(yù)報方程在實際的雷電潛勢業(yè)務(wù)中有一定的參考作用。

［1］張義軍，劉欣生，Krehbie P R..雷暴中的反極性放電和電荷結(jié)構(gòu)［J］.科學(xué)通報，2002，47(15):1192－1195.

［2］張廷龍，郄秀書，言穆弘.青藏高原雷暴的閃電特征及其成因探討［J］.高原氣象，2007，26(4):774－782.

［3］郝瑩，姚葉青，陳焱，等.基于對流參數(shù)的雷暴潛勢預(yù)報研究［J］.氣象，2007，33 (1):51－56.

［4］鄭傳烈.用獨立參數(shù)條件概率法做汕頭機(jī)場雷暴預(yù)報［J］.廣東氣象，2000 (S1):38－39.

［5］AGOSTINOMANZATO.Sounding－derived indices for neural network based Short－term thunderstorm and rainfall forecasts［J］.Atmospheric Research，2007，83:49－365.

［6］金米娜，陳云輝，徐星生.2007年夏季兩次強(qiáng)雷電過程的氣象衛(wèi)星資料分析［J］.氣象與減災(zāi)研究，2009，32(3):31－35.

［7］楊秀莊，盧璐，王宇.貴州雷電時空分布及強(qiáng)度特征分析［J］.貴州氣象，2008，32(5):4－7.