王姝蘇　李敏　丁偉　顧浩

摘要? ? 利用ECMWF、GRAPES_GFS、GRAPES_MESO、JAPAN 4個模式1～3 d日降水量預(yù)報結(jié)果，并以2 402個國家地面氣象觀測站24 h累計降水量作為觀測資料，基于全國站點進(jìn)行降水預(yù)報的消空處理，并與直接插值到站點的模式降水預(yù)報值進(jìn)行對比。結(jié)果表明，經(jīng)過消空處理后，各個模式不同時間起報、不同預(yù)報時效無雨空報的現(xiàn)象相比直接插值都得到明顯改善，因為晴雨預(yù)報準(zhǔn)確率更高，空報率更低;20：00起報的直接插值空報現(xiàn)象更嚴(yán)重，消空處理后空報現(xiàn)象少于8：00起報的;GRAPES_MESO是4個模式中直接插值空報現(xiàn)象最少的，直接插值空報現(xiàn)象相對較多的模式經(jīng)過消空處理后空報現(xiàn)象改善也更明顯。

關(guān)鍵詞? ? 日降水;站點;晴雨預(yù)報準(zhǔn)確率;空報率;消空

中圖分類號? ? P457.6? ? ? ? 文獻(xiàn)標(biāo)識碼? ? A

文章編號? ?1007-5739（2020）09-0201-02? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ?開放科學(xué)（資源服務(wù)）標(biāo)識碼（OSID）

數(shù)值預(yù)報在氣象預(yù)報業(yè)務(wù)中發(fā)揮著越來越重要的作用，成為預(yù)報員日常預(yù)報不可或缺的參考產(chǎn)品。相比其他要素，降水存在不連續(xù)性、非正態(tài)分布的特點，因而一直是預(yù)報中的重點和難點，同時也是生產(chǎn)生活關(guān)注的熱點[1-3]。數(shù)值模式的降水預(yù)報結(jié)果常存在空報、漏報的現(xiàn)象，影響了預(yù)報的準(zhǔn)確率[4-6]。

對數(shù)值預(yù)報結(jié)果進(jìn)行統(tǒng)計后處理，可以較大程度地提高預(yù)報技巧[7-8]。王亞男等[9]利用統(tǒng)計降水尺度方法有效改進(jìn)了TIGGE各中心模式直接插值的日降水量預(yù)報效果。李? 俊等[10]基于AREM模式預(yù)報產(chǎn)品，研究發(fā)現(xiàn)頻率匹配法能顯著改善模式日降水預(yù)報中雨量和雨區(qū)范圍的系統(tǒng)性偏差，訂正后降水預(yù)報的范圍和平均強(qiáng)度與實況更加接近。劉綠柳等[11]采用百分位訂正方法對全球氣候模式（GCM）輸出的日降水量進(jìn)行了處理，減小了降水量偏差。王佳津等[12]基于T639模式建立最優(yōu)多因子降尺度降水預(yù)報模型，發(fā)現(xiàn)對日降水量的預(yù)報效果優(yōu)于T639模式直接輸出的結(jié)果。

“有雨或無雨”的定性降水預(yù)報在天氣預(yù)報中備受關(guān)注，本文針對模式預(yù)報存在的無雨空報現(xiàn)象進(jìn)行研究，運(yùn)用預(yù)報中常用的4種數(shù)值模式預(yù)報產(chǎn)品，通過消空處理，提高“有雨或無雨”的定性降水預(yù)報技巧。

1? ? 資料與方法

1.1? ? 資料來源

本文采用的預(yù)報資料是ECMWF、GRAPES_GFS、GRAPES_ MESO、JAPAN 4個模式24 h累計降水量預(yù)報，預(yù)報時段取2018年5月1日至9月30日，起報時間為8：00和20：00（北京時間），預(yù)報時效為24～72 h，預(yù)報間隔24 h，空間范圍取東經(jīng)70°～140°、北緯15°～55°，空間分辨率依次分別為0.125°×0.125°、0.25°×0.25°、0.1°×0.1°、0.5°×0.5°。

觀測資料是全國范圍內(nèi)2 402個國家地面氣象觀測站24 h累計降水量，觀測時段取與預(yù)報資料相對應(yīng)的2018年5月2日至10月3日的8：00和20：00。

1.2? ? 研究方法

1.2.1? ? 消空方法。對于每個站點，將模式預(yù)報值與實況值進(jìn)行比較，統(tǒng)計無雨日的模式雨量，按大小排序，選取一個臨界值，雨量低于臨界值設(shè)為無雨。為了有效減少無雨空報的現(xiàn)象，并使處理后的漏報也盡量減少，本文選取雨量臨界值使無雨空報現(xiàn)象減少90%。消空處理方法主要是針對小雨空報問題，并不會影響其他量級降水的預(yù)報結(jié)果。

1.2.2? ? 檢驗方法。

式中：PC為晴雨預(yù)報準(zhǔn)確率，NA為降水預(yù)報正確的站（次）數(shù)，NB為空報的站（次）數(shù)，NC為漏報的站（次）數(shù)，ND為無降水預(yù)報正確的站（次）數(shù)。

式中：FAR是降水空報率，NA是預(yù)報正確的站（次）數(shù)，NB是空報的站（次）數(shù)。

2? ? 結(jié)果與分析

將4個模式的降水預(yù)報值利用雙線性插值的方法統(tǒng)一插值到全國2 402個國家地面氣象觀測站點上，然后利用消空方法對模式預(yù)報值進(jìn)行處理，并與直接雙線性插值的預(yù)報結(jié)果進(jìn)行對比，計算晴雨預(yù)報準(zhǔn)確率，結(jié)果如圖1所示。由圖1可以看出，4個模式直接插值到站點的預(yù)報結(jié)果與觀測資料相比，晴雨預(yù)報準(zhǔn)確率各有差異，GRAPES_MESO在各個預(yù)報時效的晴雨預(yù)報準(zhǔn)確率均相對較高，ECMWF和GRAPES_GFS的晴雨預(yù)報準(zhǔn)確率比較相近;8：00起報的4個模式晴雨預(yù)報準(zhǔn)確率在1～3 d預(yù)報時效均高于20：00起報的;GRAPES_MESO和JAPAN模式隨著預(yù)報時效的延長，晴雨預(yù)報準(zhǔn)確率不斷減小，而且JAPAN在24 h預(yù)報時效的晴雨預(yù)報準(zhǔn)確率明顯高于48 h和72 h，而ECMWF和GRAPES_GFS 2個模式晴雨預(yù)報準(zhǔn)確率在各個預(yù)報時效相差不大。經(jīng)過消空處理后，4個模式8：00和20：00起報的晴雨預(yù)報準(zhǔn)確率在各個預(yù)報時效都有明顯的提高，20：00起報的晴雨預(yù)報準(zhǔn)確率提高程度均大于8：00起報的。GRAPES_MESO經(jīng)過消空處理后，晴雨預(yù)報準(zhǔn)確率的提高程度是4個模式中最小的，但晴雨預(yù)報準(zhǔn)確率仍然是4個模式中最高的。

為了進(jìn)一步研究消空處理對降水空報現(xiàn)象的改進(jìn)情況，圖2分別計算了4個模式8：00、20：00起報的直接插值和消空處理后無雨的空報率。將直接插值的預(yù)報結(jié)果與觀測資料進(jìn)行對比，4個模式8：00起報的降水空報率均小于20：00起報的;GRAPES_MESO是4個模式中空報率在各個預(yù)報時效都相對最小的，ECMWF和GRAPES_GFS的空報率在各個預(yù)報時效都比較接近，JAPAN在24 h預(yù)報時效的空報率明顯小于48 h和72 h預(yù)報時效的。經(jīng)過消空處理，4個模式8：00、20：00起報的降水空報率都在20%上下，相比直接插值在各個預(yù)報時效也都有了明顯的減小，空報現(xiàn)象得到了較大的改善;直接插值空報率相對較高的模式經(jīng)過消空處理，空報率的減小幅度也相對更大。

3? ? 結(jié)論

本文利用消空方法對ECMWF、GRAPES_GFS、GRAPES_ MESO、JAPAN 4個模式8：00和20：00起報的日累計降水量進(jìn)行了處理，并與直接插值到站點的模式預(yù)報值進(jìn)行對比，得出以下結(jié)論。

（1）經(jīng)過消空處理后，各個模式相比直接插值，不同時間起報、不同預(yù)報時效的晴雨預(yù)報準(zhǔn)確率都有了明顯的提高，無雨的空報率也明顯減小，無雨空報的現(xiàn)象得到了明顯的改善。

（2）總體來看，各模式在1～3 d預(yù)報時效20：00起報直接插值的空報現(xiàn)象相比8：00更嚴(yán)重，因為20：00起報的晴雨預(yù)報準(zhǔn)確率更低，空報率更高;經(jīng)過消空處理后，各個模式在大多數(shù)預(yù)報時效20：00起報的空報現(xiàn)象少于8：00起報的，20：00起報的空報改善情況比8：00更明顯。

（3）就直接插值結(jié)果看，GRAPES_MESO是4個模式中空報現(xiàn)象最少的，經(jīng)過消空處理后，ECMWF是4個模式中空報現(xiàn)象最少的。直接插值空報現(xiàn)象相對嚴(yán)重的模式經(jīng)過消空處理后的空報改善情況更明顯。

4? ? 參考文獻(xiàn)

[1] 智協(xié)飛，王姝蘇，周紅梅，等.我國地面降水的分級回歸統(tǒng)計降尺度預(yù)報研究[J].大氣科學(xué)學(xué)報，2016，39（3）：329-338.

[2] 王姝蘇，智協(xié)飛，俞劍蔚，等.降水的統(tǒng)計降尺度預(yù)報及其空間相關(guān)性和時間連續(xù)性重建[J].大氣科學(xué)學(xué)報，2018，41（1）：36-45.

[3] 王海霞，智協(xié)飛.基于TIGGE多模式降水量預(yù)報的統(tǒng)計降尺度研究[J].氣象科學(xué)，2015，35（4）：430-437.

[4] 孫軍，馬學(xué)款，代刊，等.北京幾次弱降水過程預(yù)報失誤分析[J].氣象，2018，44（1）：26-41.

[5] 翟振芳，魏春璇，鄧斌，等.安徽省ECMWF數(shù)值模式降水預(yù)報性能的檢驗[J].氣象與環(huán)境學(xué)報，2017，33（5）：1-9.

[6] 高榮珍，李欣.青島短期降水空漏報及小量降水預(yù)報指標(biāo)分析[J].山東氣象，2015，35（4）：18-22.

[7] 代刊，朱躍建，畢寶貴.集合模式定量降水預(yù)報的統(tǒng)計后處理技術(shù)研究綜述[J].氣象學(xué)報，2018，76（4）：493-510.

[8] 畢寶貴，代刊，王毅，等.定量降水預(yù)報技術(shù)進(jìn)展[J].應(yīng)用氣象學(xué)報，2016，27（5）：534-549.

[9] 王亞男，智協(xié)飛.多模式降水集合預(yù)報的統(tǒng)計降尺度研究[J].暴雨災(zāi)害，2012，31（1）：1-7.

[10] 李俊，杜鈞，陳超君.降水偏差訂正的頻率（或面積）匹配方法介紹和分析[J].氣象，2014，40（5）：580-588.

[11] 劉綠柳，任國玉.百分位統(tǒng)計降尺度方法及在GCMs日降水訂正中的應(yīng)用[J].高原氣象，2012，31（3）：715-722.

[12] 王佳津，陳朝平，王春學(xué).BP-CCA方法用于四川盆地夏季日降水量的可預(yù)報性研究[J].高原山地氣象研究，2013，33（4）：35-41.