劉紅武 , 黃楚惠 , 王彬雁 , 牛金龍 , 尹忠海 , 唐 林

（1.湖南省氣象臺, 長沙 410118；2.四川省氣象臺/高原與盆地暴雨旱澇災(zāi)害四川省重點實驗室, 成都 610072；3.成都市氣象局, 成都 610071；4.湖南省人工影響天氣領(lǐng)導小組辦公室, 長沙 410118）

引言

隨著社會經(jīng)濟的發(fā)展和城市化進程的加快，特別是交通、通訊的迅猛發(fā)展，低能見度天氣對人們?nèi)粘Ｉ畹挠绊懭遮吤黠@[1]。近年來，國內(nèi)學者對能見度問題開展了諸多研究，主要集中于大氣能見度趨勢、影響機理、數(shù)值模擬和預報方法等方面[2-10]。對于能見度預報訂正方法研究，常見的有利用多元回歸、相關(guān)系數(shù)、均方根誤差等對自動站數(shù)據(jù)進行訂正[11-13]。胡海川等[14]采用BP神經(jīng)網(wǎng)絡(luò)方法構(gòu)建環(huán)渤海沿海城市能見度預報模型，預報效果明顯高于ECMWF模式；黃楚惠等[15]在四川盆地低能見度的氣象要素訂正閾值基礎(chǔ)上結(jié)合概率密度匹配的方法，對四川盆地低能見度天氣進行訂正并取得了較好的效果；王媛媛等[16]對北京區(qū)域各站點能見度觀測值與模式預報結(jié)果進行比較和分時段逐級偏差訂正，改善了模式對高海拔地區(qū)的低能見度高估現(xiàn)象，也能更好地預報出低能見度現(xiàn)象。

總之，我國對于能見度的研究已取得了顯著成果，但大多集中于東部大中城市。湖南地處云貴高原向江南丘陵和南嶺山脈向江漢平原過渡的地帶，地勢呈三面環(huán)山、朝北開口的馬蹄形地貌，低能見度天氣多發(fā)。目前，針對湖南能見度長期演變特征的成果并不多見，對其預報訂正方法的研究更是匱乏。因此，在了解湖南低能見度天氣時空分布特征的基礎(chǔ)上，開展能見度預報訂正的相關(guān)研究，有利于提高低能見度預警預報能力，對減輕低能見度天氣所造成的災(zāi)害損失具有十分重要的意義。本文利用2011～2020年國家基本站觀測資料，研究湖南省低能見度時空分布特征及低能見度與地面氣象要素的關(guān)系，獲取低能見度與要素之間的預報訂正閾值，進而基于ECMWF模式（簡稱EC模式，下文同）采用概率密度匹配結(jié)合要素的雙重試驗訂正方法，對湖南省兩次低能見度天氣個例預報進行訂正，評估模式訂正前后的預報效果，以期為湖南乃至其他地區(qū)能見度預報提供參考。

1 資料和方法

1.1 資料

本文選取2011～2020年湖南省國家基本站（共計97站）逐時觀測資料，包括能見度、地面風速、24 h變壓、相對濕度、地面溫度等要素。個例預報訂正選取EC模式2019年12月23日20時（北京時，下同）起報和2020年2月10日20時起報的逐3 h能見度預報資料。

1.2 方法

1.2.1 低能見度天氣定義

對湖南省97站能見度資料進行統(tǒng)計，得到湖南地區(qū)能見度變化特征。以≤1km為標準，統(tǒng)計了湖南地區(qū)低能見度的分布情況。此外，結(jié)合地面風速、24 h變壓等要素，分級統(tǒng)計不同能見度下的要素分布比例，進一步確定低能見度的要素訂正閾值。

1.2.2 能見度訂正方法

訂正方法采用概率密度匹配結(jié)合要素的雙重試驗訂正方法。由于概率密度匹配方法曾被廣泛用于具有偏差非獨立特性（即誤差隨著資料的大小、時間及區(qū)域的不同而變化）的資料誤差訂正上，可以很好地消除多元線性回歸值與實況值的非獨立系統(tǒng)誤差[17]。因此，本文利用該方法結(jié)合氣象要素進行雙重訂正試驗，以期增加低能見度事件的命中率，提高能見度預報水平。

由于Gamma分布可以很好地對各種物理量進行擬合，并且適用于各地任一累積時段物理量，具有一定的普適性。因此，本文選取Gamma分布對能見度進行訂正。Gamma分布的概率函數(shù)[18]如下：

式（1）中： Γ (α)是 Gamma函數(shù)，其函數(shù)公式如式（2）所示；α為Gamma函數(shù)的形狀參數(shù)，表示分布曲線的形狀，值越大（小），曲線越平緩（陡峭），小值出現(xiàn)的概率越?。ù螅?；β為Gamma函數(shù)的尺度參數(shù)，值越大（?。怕屎瘮?shù)取值范圍越大（?。植荚椒稚ⅲ校?。

α和β的估算采用最大似然法，具體公式如下：

式（3）和（4）中，u為序列的均值，σ2為序列的方差。

1.2.3 能見度 TS 評分

根據(jù)中國氣象局制訂的“重要天氣預報質(zhì)量評定辦法”，TS（Threat Score）評分計算公式為:

式中：TS為預報準確率，Na、Nb、Nc分別為預報正確、漏報、空報的站(次)數(shù)。TS評分介于[0，1]，體現(xiàn)了預報的準確程度[19]。本文將能見度分級檢驗區(qū)間分為 0～1 km、1～3 km、3～10 km 及10 km 以上。

2 能見度時空分布特征

2.1 低能見度天氣空間分布特征

由2011～2020年湖南地區(qū)各季能見度值空間分布（圖1）可知，湖南北部能見度值顯著低于南部，春季及秋冬季低能見度分布范圍大于夏季，冬季低能見度分布范圍最廣?！? km低能見度區(qū)域主要位于湖南中北部、通道縣及衡陽市附近，并且從秋季開始整個低能見度區(qū)域有明顯向南擴大的趨勢，低能見度中心分布不均，湘東北的洞庭湖區(qū)最低，山區(qū)至平原的過渡地帶次之。

圖1 2011～2020 年湖南省各季能見度空間分布（a.春季 3～5 月，b.夏季 6～8 月，c.秋季 9～11 月，d.冬季 12 月～次年 2 月，單位：km）

2.2 低能見度天氣出現(xiàn)時數(shù)空間分布

通過對2011～2020年湖南省各季低能見度（≤1 km）出現(xiàn)時數(shù)進行統(tǒng)計，發(fā)現(xiàn)各季南岳（57776）和郴州（57972）是出現(xiàn)時數(shù)最多的地區(qū)，南岳站冬季低能見度出現(xiàn)時數(shù)達350 h，即使是在夏季也有240 h，郴州站次之，春、秋和冬季低能見度出現(xiàn)時數(shù)分別為117 h、123 h和193 h。春季低能見度時數(shù)分布范圍較秋季稍大，多集中在8～25 h，主要位于湖南東北至南部一帶（圖2a）。夏季低能見度時數(shù)分布范圍最小，主要分布在中北部，多集中在 3～12 h，常德（57662）低能見度時數(shù)（約55 h）僅次于南岳（圖2b）。秋季低能見度時數(shù)分布范圍較夏季明顯增多，平均為5～15 h，其中桃江超過55 h（圖2c）。冬季低能見度時數(shù)分布范圍明顯擴大，與春季分布相似，中東部低能見度時數(shù)更為集中，平均 10～30 h，馬坡嶺達 108 h（圖2d）。

圖2 2011～2020年湖南省各季低能見度（≤1km）出現(xiàn)時數(shù)空間分布（a.春季3～5月，b.夏季6～8月，c.秋季9～11月，d.冬季12月～次年2月，單位：h）

2.3 低能見度天氣過程日變化

為了更好地揭示低能見度天氣日變化特點，了解其生消情況，選取探測環(huán)境較好且發(fā)生低能見度天氣概率相對較多的常德站（57662）為例，給出了2011～2020年常德站逐月能見度的日變化特征（圖3）。如圖所示，能見度逐月日變化均呈現(xiàn)單峰型，能見度低值時段主要集中在20時～次日09時，18時之后能見度逐漸降低，次日09時后能見度逐漸轉(zhuǎn)好。能見度分布存在明顯的季節(jié)變化特征，冬季的12月、1月為能見度低值月份。從能見度逐時時數(shù)分布（圖4）來看，≤1 km低能見度出現(xiàn)時數(shù)10次以上的情況主要發(fā)生在22時～次日12時。

圖3 2011～2020年常德站逐月能見度的日變化特征（a.3～5月，b.6～8月，c.9～11月，d.12～2月，單位：km）

圖4 常德站能見度分等級逐時時數(shù)分布

3 能見度與地面要素的關(guān)系

選取地面要素（風速、24 h變壓、相對濕度、地面溫度）作為研究對象，分別對不同能見度閾值出現(xiàn)時，地面要素在不同區(qū)間的變化情況進行統(tǒng)計分析。如圖所示，＜3 km能見度主要出現(xiàn)在風速為1～2 m/s，≤1 km低能見度大都在風速＜2 m/s下出現(xiàn)（圖5a）；低能見度主要出現(xiàn)在20℃以下，并且溫度越高出現(xiàn)越低等級能見度可能性越?。▓D5b）；≤1 km 低能見度24 h變壓范圍＜ 2 hPa（圖5c）；相對濕度 80%～90% 易出現(xiàn)在1～3 km能見度，能見度≤1 km的相對濕度大多處于90%以上（圖5d）。綜合各要素情況來看，低能見度與地面風速、相對濕度及溫度相關(guān)性較好，與地面氣壓相關(guān)性一般。因此，進一步統(tǒng)計得到能見度≤1 km的訂正閾值如表1所示。當最大風速在2 m/s以下，低能見度事件出現(xiàn)頻率是87.47%；溫度在20℃以下時，低能見度事件出現(xiàn)頻率是82.88%；相對濕度在90%以上時，低能見度事件出現(xiàn)頻率是95.79%；24 h變壓在2 hPa以下時，低能見度事件出現(xiàn)頻率是81.49%。

表1 能見度≤1 km訂正閾值

圖5 不同能見度等級下地面要素在不同區(qū)間占比（a.風速，b.地面溫度，c.24 h變壓，d.相對濕度)

4 基于 EC 模式的能見度預報個例訂正

由上述統(tǒng)計可知，低能見度主要出現(xiàn)在20時～次日09時，因此，本節(jié)重點關(guān)注23時～次日11時的訂正能力。先采用概率匹配訂正方案進行訂正，具體過程如下：首先對能見度閾值進行分區(qū)，這里取 0、1、3、10和30 km，分別算出EC模式及觀測對應(yīng)的累積概率值。對某一區(qū)間，假設(shè)預報概率值與實況概率值應(yīng)該一致，反求得 0、1、3、10和30 km 對應(yīng)的訂正閾值，

再利用內(nèi)插法得到任意值的訂正值。在此基礎(chǔ)上，結(jié)合第 3 節(jié)對要素閾值的分析，加入云量、風速、相對濕度、邊界層高度等進行訂正，采用概率匹配法和要素訂正相結(jié)合的方法對能見度預報進行再次訂正。

4.1 2019 年 12 月 25 日湖南中東部低能見度天氣過程

從實況上看，24 日 23 時（圖6a）湖南中東部出現(xiàn)了分散的低能見度＜1 km區(qū)域，到25日05時（圖6i），低能見度站數(shù)略有增加，仍位于湖南中東部區(qū)域，08時（圖6m）開始低能見度范圍已減小。EC預報24 日23 時～25 日 11 時（圖6b、f、j、n、r）湖南僅有個別站點出現(xiàn)了＜1 km低能見度天氣，湖南中西部及東北部能見度約為1～3 km，對低能見度天氣出現(xiàn)了明顯的漏報。因此，EC模式對這次低能見度過程預報范圍較實況顯著偏小，漏報了東北部和中南部的低能見度區(qū)域，效果并不理想。通過概率匹配法訂正后（圖6c、g、k、o、s），量級得到了很好的調(diào)整；但在24日23時低能見度范圍較實況明顯偏大，空報了湖南西部低能見度；25日02時和05時，對于模式未預報的低能見度區(qū)域訂正效果也不佳；08時和11時，量級雖有調(diào)整，仍存在對湖南西部的空報。通過概率匹配加要素雙重訂正后（圖6d、h、l、p、t），對湖南中東部（紅圈所示）的低能見度有較好反映，各時次低能見度預報更接近實況，23時和02時對湖南西部（藍圈所示）低能見度區(qū)域空報范圍較單純使用概率匹配的方法明顯減小，表明概率匹配結(jié)合要素的訂正方法對部分低能見度區(qū)域有較好的訂正效果。

圖6 2019 年 12 月低能見度實況（a、e、i、m、q）、EC 模式 23 日 20 時預報（b、f、j、n、r）、概率匹配訂正（c、g、k、o、s）和概率匹配+物理量訂正（d、h、l、p、t）結(jié)果對比（a～d.24 日 23 時，e～h.25 日 02 時，i～l.25 日 05 時，m～p.25 日 08 時，q～t.25 日 11 時）

圖7給出了過程訂正后能見度預報的TS評分。從評分結(jié)果來看，各時次＜1 km能見度TS評分在兩種方案中均為訂正后明顯高于訂正前，其中02時（圖7b，紅圈所示）評分最高約為 0.24。 此外，結(jié)合空報率（圖略）來看，概率匹配結(jié)合要素的訂正方案比僅使用概率匹配的訂正方案效果更好。

圖7 EC模式預報、概率匹配訂正及概率匹配+物理量訂正TS評分對比（a.24日23時，b.25日02時，c.25日05時，d.25日08時，e.25日11時）

4.2 2020 年 2 月 12 日湖南大部低能見度天氣過程

由11日23時能見度實況（圖8a）可知，湖南中部和東北部出現(xiàn)了＜1 km低能見度天氣，且范圍逐漸擴大，至12日08時（圖8m），湖南大部地區(qū)出現(xiàn)了低能見度天氣，11時（圖8q）明顯減小，而EC模式10日20時起報相應(yīng)時次低能見度區(qū)域較實況顯著偏小，僅12日02時和05時在湖南北部有所體現(xiàn)（圖8b、f、j、n、r）。單純應(yīng)用概率匹配方案訂正后（圖8c、g、k、o、s），11日23時～12日05時湖南省能見度量級降低，＜1 km范圍低能見度區(qū)和EC模式預報差別不大，仍然漏報了湖南南部低能見度區(qū)。通過應(yīng)用概率匹配加要素雙重訂正后（圖8d、h、l、p、t），湖南南部低能見度區(qū)體現(xiàn)出來，雖然低能見度范圍較實況偏大（特別是12日02和05時），但也較好地指示了湖南大范圍低能見度天氣，且12日08時在湖南中部一帶預報了＜3 km能見度區(qū)域。此外，從各時次TS評分（圖9）來看，兩種方案訂正后TS評分均比模式預報要高，尤其是12日05時（圖9c），概率匹配結(jié)合要素的訂正方案TS評分達到0.55。綜合兩次個例訂正效果來看，概率匹配結(jié)合要素的訂正方案優(yōu)于僅應(yīng)用概率匹配的訂正方案。

圖8 2020 年 2 月低能見度實況（a、e、i、m、q），EC 模式 10 日 20 時預報（b、f、j、n、r），概率匹配訂正（c、g、k、o、s）和概率匹配+物理量訂正（d、h、l、p、t）結(jié)果對比（a～d.11 日 23 時，e～h.12 日 02 時，i～l.12 日 05 時，m～p.12 日 08 時，q～t.12 日 11 時）

圖9 EC模式預報、概率匹配訂正及概率匹配+物理量訂正TS評分對比（a.10日23時，b.11日02時，c.11日05時，d.11日08時，e.11日11時）

5 結(jié)論與討論

本文利用2011～2020年湖南省國家基本站觀測資料，研究了湖南低能見度天氣時空分布特征，并結(jié)合低能見度與地面氣象要素的關(guān)系獲取預報訂正閾值，在此基礎(chǔ)上對比分析了概率匹配和概率匹配結(jié)合要素的訂正方法對EC模式能見度預報的訂正效果，得到如下主要結(jié)論：

（1）低能見度區(qū)域主要位于湖南中北部、通道縣及衡陽市附近，并且從秋季開始整個低能見度區(qū)域有明顯向南擴大的趨勢；低能見度出現(xiàn)時數(shù)最多的是冬、春和秋季，冬季低能見度時數(shù)分布范圍與春季類似，范圍最廣，集中于湖南中東部。低能見度天氣日變化呈單峰型，主要集中在20時～次日09時。

（2）≤1 km低能見度主要出現(xiàn)在風速＜2 m/s、相對濕度基本高于90%、地面溫度低于20℃及24 h變壓＜2 hPa的氣象環(huán)境。

（3）基于EC模式的能見度預報產(chǎn)品，采用概率密度匹配結(jié)合要素的訂正方案優(yōu)于僅應(yīng)用概率匹配的訂正方案，可以很好地對湖南大部分能見度預報進行訂正，訂正后TS評分顯著提高。

本文僅對湖南兩次低能見度天氣過程進行了訂正試驗，后期需挑選更多個例進行長時間序列的訂正，以更好評估所用方法對湖南能見度的訂正效果。此外，由于影響能見度變化的要素較多，如降水強度、近地面層溫差等，而本文訂正方法僅使用了6個要素，存在一定的局限性，且未分季節(jié)進行要素關(guān)系的分析，這些問題將在下一步工作中深入探討。