劉 靜 任 川 趙梓淇 陳傳雷 王 瀛 才奎志
1 遼寧省氣象災(zāi)害監(jiān)測預(yù)警中心,沈陽 110166 2 中國氣象局沈陽大氣環(huán)境研究所,沈陽 110166 3 沈陽市氣象局,沈陽 110166 4 遼寧省氣象信息中心,沈陽 110166
提 要:利用氣象大數(shù)據(jù)云平臺中逐小時降水資料,基于目標對象檢驗法和鄰域法,評估2019—2020年遼寧主汛期降水過程中中國氣象局上海數(shù)值預(yù)報模式系統(tǒng)(CMA-SH9)、中國氣象局中尺度天氣數(shù)值預(yù)報系統(tǒng)(CMA-MESO)、中國氣象局睿圖東北數(shù)值預(yù)報模式系統(tǒng)(CMA-DB)的預(yù)報性能。結(jié)果表明:千米尺度或接近千米尺度的上述三個模式,在36 h時效內(nèi),對于累積強降水(12 h 降水量≥50 mm)落區(qū)形態(tài)預(yù)報與實況有相似性,落區(qū)質(zhì)心預(yù)報偏差一般在20 km左右。然而,預(yù)報落區(qū)與實況重疊的面積一般都在10%以下,個別情形下(如CMA-MESO對于氣旋型降水過程)累積強降水落區(qū)預(yù)報與實況重疊度能夠接近20%;位置偏離的直接結(jié)果是導(dǎo)致漏報率高(一般在75%左右,CMA-MESO模式漏報率略低,為10%~20%),其中高壓后部型降水過程中累積強降水的漏報率超過80%,位置偏離也造成較高空報率。對于短時強降水(1 h降水量≥20 mm)預(yù)報,在方圓40 km內(nèi)不計偏差情況下,各模式預(yù)報命中率平均在10%以下(最大值為9.2%),空報率平均為58.7%;三種降水類型中,模式對臺風(fēng)型降水過程的短時強降水預(yù)報性能最低。
遼寧省暴雨具有頻發(fā)性、歷時短、致災(zāi)重等特點(孫欣等,2016)。受穩(wěn)定的副熱帶高壓(以下簡稱副高)影響,低層暖濕空氣和中高緯西風(fēng)帶高空槽為重要的暴雨影響系統(tǒng)(陳傳雷和楊青,2018)。暴雨多發(fā)生在中尺度系統(tǒng)觸發(fā)下(劉靜等,2016),區(qū)域模式預(yù)報和全球模式相比,對中尺度系統(tǒng)誘發(fā)的暴雨表現(xiàn)出較好的預(yù)報能力(漆梁波和徐珺,2018),因此加深對區(qū)域模式的了解,有助于制作暴雨預(yù)報和發(fā)布暴雨紅色預(yù)警信號(劉靜等,2021),檢驗評估是提升對區(qū)域模式的認知水平和應(yīng)用能力的有效途徑(Gofa et al,2018;Chakraborty,2010)。
在檢驗技術(shù)方面,全國防災(zāi)減災(zāi)標準化技術(shù)委員會(SAC/TC 345)(2017)給出了誤差檢驗等常規(guī)方法,這些方法雖然可以評價預(yù)報正確與否或準確程度,卻往往忽略了預(yù)報中一些積極的信息(戴建華等,2013)。目前對高分辨率區(qū)域模式的檢驗評估,多應(yīng)用空間檢驗技術(shù)。Gilleland et al(2009)將現(xiàn)有的空間檢驗技術(shù)歸納為兩類:一類是通過空間濾波的方式對原始場進行處理進而統(tǒng)計各個檢驗物理量,通過對比模式與實況降水場的統(tǒng)計量差異對模式預(yù)報性能進行評估,以鄰域法(模糊檢驗法)為代表;另一類是通過從預(yù)報場和實況場中識別目標個體,在相互匹配的過程中計算兩者之間的空間位移形變統(tǒng)計量,以基于目標對象檢驗法(method of object-based diagnostic evaluation,MODE)為代表。李佰平等(2016)認為鄰域法能夠在不同尺度和評價策略上給出與預(yù)報有關(guān)的更多信息。鄰域法通常采用分數(shù)技巧評分(fractional skill score,F(xiàn)SS)作為評價指標,王新敏和栗晗(2020)研究指出,F(xiàn)SS評分相對于傳統(tǒng)的TS評分能夠更好地通過量化方式反映出不同模式的預(yù)報能力差別。劉靜等(2019)應(yīng)用鄰域法對遼寧地區(qū)的雷達回波組合反射率因子進行檢驗,得出在11 km的鄰域半徑下FSS評分較高。在鄰域法應(yīng)用技術(shù)上,潘留杰等(2017)根據(jù)不同空間尺度和鄰域半徑評估數(shù)值模式降水預(yù)報的CSI 評分表現(xiàn),唐文苑等(2017)應(yīng)用固定鄰域半徑(40 km)對強對流天氣進行檢驗,得出空報率、漏報率等檢驗評分。馬申佳等(2018)認為,改良后的鄰域概率法可以在時空尺度上綜合評估對流尺度降雨的特征。美國NCAR中心研發(fā)了MODE,可以客觀地對比分析預(yù)報和實況的位置、形狀、方向和尺寸等相關(guān)屬性(Davis et al,2006a;2006b)。茅懋等(2016)設(shè)計了和MODE類似的目標對象檢驗方法,量化評價預(yù)報對象的面積、強度等要素的吻合程度及對應(yīng)的偏差,給出常規(guī)檢驗TS評分偏低的原因。在預(yù)報檢驗的時間尺度方面,Duc et al(2013)在高分辨率模式檢驗中引入了時間維度并進行了初步分析。
在檢驗結(jié)果方面,全球模式對中到大雨的預(yù)報能力有限,空報和漏報都很明顯(王曉峰和周榮衛(wèi),2021)。多區(qū)域中心數(shù)值模式的對比結(jié)果顯示,北京、沈陽、上海、廣州區(qū)域模式的預(yù)報穩(wěn)定性較好(Chen et al,2013)。華東區(qū)域數(shù)值模式預(yù)報系統(tǒng)SMS-WARR(9 km),(現(xiàn)為CMA-SH9)經(jīng)過版本的更新,提升了對降水的預(yù)報性能(徐同等,2016),對西南渦誘發(fā)的暴雨過程能夠指示出強降水空間落區(qū),但強度有偏差(徐同等,2019)。CMA-SH9對臺風(fēng)快速移動過程的預(yù)報能力偏弱(薛文博等,2020)。羅月琳等(2021)對GRAPES-GEPS(現(xiàn)為CMA-GEPS)環(huán)流集合預(yù)報的分類釋用方法進行研究與檢驗,劃分環(huán)流形勢類型,從而為天氣預(yù)報提供參考依據(jù),GRAPES區(qū)域集合預(yù)報系統(tǒng)對中國汛期降水具有很好的預(yù)報能力(王婧卓等,2021)。與華東模式、華南模式相比,中國氣象局中尺度天氣數(shù)值預(yù)報系統(tǒng)(CMA-MESO)對于云南地區(qū)24 h降水量≥50 mm 的降水預(yù)報效果較好(趙寧坤等,2021)。何爽爽等(2018)通過研究指出,中國氣象局睿圖東北數(shù)值預(yù)報模式系統(tǒng)(CMA-DB)對降水的分布特征有所體現(xiàn),但對于降水中心位置和強度預(yù)報偏差較大。區(qū)域模式對暴雨的預(yù)報空報率較高(蘇翔等,2022)。
多方位評估可以挖掘多區(qū)域模式的附加值,為預(yù)報員提供更加全面和客觀的參考信息。目前,應(yīng)用空間檢驗法對多區(qū)域高分辨率模式進行檢驗的研究尚不多見,對遼寧省12 h累積強降水預(yù)報檢驗研究較少,尤其是對短時強降水預(yù)報的檢驗分析。為綜合評估多區(qū)域高分辨率模式在遼寧汛期暴雨過程中的預(yù)報表現(xiàn),對2019—2020年發(fā)生在遼寧省內(nèi)的16次天氣過程(共計408個時次的暴雨實況)進行累積強降水目標對象識別和短時強降水鄰域法檢驗,從短時預(yù)報效果分析CMA-SH9、CMA-MESO和CMA-DB的表現(xiàn),對高壓后部型、氣旋型、臺風(fēng)型天氣過程的降水量進行評估,得出不同模式預(yù)報性能優(yōu)勢和偏差特征。
1.1.1 實況資料
實況資料來自氣象大數(shù)據(jù)云平臺“天擎”。采用遼寧省1 282 個觀測站逐小時降水資料,去除異常值確保降水?dāng)?shù)據(jù)連續(xù)準確。
1.1.2 預(yù)報資料
預(yù)報數(shù)據(jù)來自業(yè)務(wù)化的多區(qū)域高分辨率模式,包括CMA-SH9、CMA-MESO及CMA-DB三個區(qū)域模式的初始場均為NCEP-GFS模式,其中CMA-SH9的空間分辨率為9 km(沒有選用3 km空間分辨率的CMA-SH3模式是因為其預(yù)報范圍不能涵蓋所有遼寧站點),其他兩種區(qū)域模式分辨率均為3 km。
1.2.1 天氣學(xué)分型
孫欣等(2016)分析了1960—2013年發(fā)生在遼寧地區(qū)的60次大暴雨過程,指出遼寧強降水的主要影響系統(tǒng)包括副高后部、氣旋(河套氣旋、江淮氣旋、蒙古氣旋等)影響及臺風(fēng)北上等。本文挑選2019—2020年發(fā)生在遼寧省內(nèi)的16次強降水過程,共計408個時次進行多區(qū)域高分辨率模式降水預(yù)報檢驗分析,降水分型分別簡稱為高壓后部型、氣旋型和臺風(fēng)型降水。具體分型時段詳見表1。
表1 不同影響天氣系統(tǒng)下的個例時段
1.2.2 不同時間定義
文中時間表述定義如下:
“08時和20時起報”表示模式起報時刻,均為北京時?!?4 h和36 h”表示預(yù)報時效。由于區(qū)域模式數(shù)據(jù)延時、本地化解析處理等問題,區(qū)域模式資料的前12 h業(yè)務(wù)應(yīng)用價值不大,因此本文選用預(yù)報時效13~36 h的結(jié)果開展工作,針對不同檢驗方法,分別對13~24 h(以下簡稱24 h)累積降水量、25~36 h(以下簡稱36 h)累積降水量及逐小時降水量進行檢驗?!? h和12 h”表示累計時間。
1.2.3 強降水定義
1.2.3.1 累積強降水
在累積強降水預(yù)報檢驗中,對多區(qū)域高分辨率模式的12 h累積降水量預(yù)報進行檢驗評估。本研究針對12 h累積降水量≥50 mm的降水展開檢驗,表述為“12 h累積強降水”(以下簡稱為累積強降水)。
1.2.3.2 短時強降水
在短時強降水預(yù)報檢驗中,對多區(qū)域高分辨率模式的1 h降水量預(yù)報進行檢驗評估。根據(jù)中國氣象局2020年6月印發(fā)的《智能預(yù)報技術(shù)方法競賽檢驗方案》(以下簡稱《檢驗方案》)規(guī)定,在遼寧地區(qū)定義1 h降水量≥20 mm的降水為短時強降水。
1.2.4 檢驗方法
區(qū)域模式網(wǎng)格預(yù)報采用鄰近點插值方法插值到檢驗站點后和對應(yīng)的實況檢驗站點觀測進行對比檢驗。
有研究表明,模式空間分辨率改變并不是對所有的降水個例都有顯著影響(于翡等,2018),因此本文嘗試對不同空間分辨率的多區(qū)域模式進行對比分析,指出各區(qū)域模式性能優(yōu)勢和偏差特征。
1.2.4.1 累積強降水預(yù)報目標對象檢驗
基于MODE,對表1中所有個例進行檢驗,評估多區(qū)域高分辨率模式對累積強降水的預(yù)報效果,其結(jié)果是對各模式08時和20時起報的產(chǎn)品進行的整合檢驗。
MODE(Davis et al,2009)可以辨識預(yù)報二維場的空間特征,它除了給出臨界成功指數(shù)(critical success index,CSI)評分等傳統(tǒng)的檢驗指標外,還可以確定空間特征在不同場中的相似程度,進而開展相似性統(tǒng)計。MODE主要步驟包括:卷積處理、去背景、確定對象及對象屬性等(尤鳳春等,2011;茅懋等,2016)。
本文應(yīng)用MODE檢驗方法對累積強降水預(yù)報進行空間分布特征分析,檢驗區(qū)域選為38°~44°N、118°~126°E,降水量閾值設(shè)定為12 h降水量≥50 mm,采用預(yù)報和實況目標對的位置評價、縱橫比、面積比、質(zhì)心距離、重疊面積與聯(lián)合面積比、CSI評分、漏報率、總評分等屬性來評價多區(qū)域模式的預(yù)報技巧。
MODE包涵的六類評分指標的具體含義及算法:
(1)位置評價。預(yù)報目標相對于實況目標的位置信息、質(zhì)心距離。
(2)縱橫比。對于某個檢驗?zāi)繕藢ο?,提取邊界點信息,通過中心的最長兩個邊界點連線即為長軸,通過重心且垂直于長軸的最短兩個邊界點連線即為短軸,短軸和長軸的比值為縱橫比。縱橫比在0~1,值越小,說明強降水預(yù)報落區(qū)帶狀分布明顯。
(3)面積比。實況面積和預(yù)報面積的比值。
(4)重疊面積/聯(lián)合面積。在格點場中,實況場和預(yù)報場交叉的面積與二者聯(lián)合面積的比值。
(5)評分指標。包括CSI評分、漏報率、總評分三個指標。
CSI、漏報率計算方法如式(1)、式(2)所示。
(1)
(2)
式中:n為參與平均的樣本數(shù),A、B、C的定義參見表2。
表2 降雨檢驗分類表
總評分是通過對所有屬性的計算,得出的總體評價指標。
1.2.4.2 短時強降水預(yù)報鄰域檢驗
本文采用鄰域法對短時強降水預(yù)報進行空間鄰域檢驗。短時強降水時空尺度小、局地性強,常規(guī)的地面觀測站很難完全觀測到,采用傳統(tǒng)的“點對點”檢驗方法難以客觀準確地反映區(qū)域模式預(yù)報質(zhì)量??臻g檢驗法在強對流天氣預(yù)報檢驗中應(yīng)用較為廣泛(鄭永光等,2015)。短時強降水檢驗多采用美國SPC的“點對面”空間檢驗辦法,即對于每一個區(qū)域站點上的短時強降水預(yù)報正確與否,是用以該站點為圓心,40 km為半徑的圓面上是否出現(xiàn)了短時強降水來判別(唐文苑等,2017),本文參考該方法對短時強降水進行檢驗。
《檢驗方案》規(guī)定短時強降水檢驗指標為命中率和空報率,本文采用這兩個指標對多區(qū)域模式的預(yù)報效果進行分析。
命中率(POD)、空報率(FAR)計算方法分別如式(3)、式(4)所示。
(3)
(4)
2.1.1 累積強降水預(yù)報目標對象檢驗
在高壓后部型降水過程中,使用區(qū)域模式和觀測數(shù)據(jù)分別做為預(yù)報場和實況場,通過MODE進行分析。
由表3可見重疊面積/聯(lián)合面積為0.040~0.052(即4.0%~5.2%),結(jié)合面積比(實況面積/預(yù)報面積)的0.274~0.466(即27.4%~46.6%),可以分析出,三個模式預(yù)報累積強降水的落區(qū)大部分(80%~86%)落在實際強降水發(fā)生區(qū)之外,這也就不難理解為何模式預(yù)報的漏報率那么高,都在83%左右(模式空報率也很大,因為面積比顯示實況面積比預(yù)報面積小得多),CSI評分只有11.6%~13.1%。
表3 高壓后部型降水過程中累積強降水預(yù)報目標對象檢驗結(jié)果
綜合得到評價結(jié)果:在高壓后部型降水過程中,三個模式的位置偏差接近,質(zhì)心距離在18~21 km,屬帶狀目標;CMA-MESO預(yù)報形態(tài)偏差最大,過高估計了累積強降水落區(qū)的帶狀分布;CMA-DB預(yù)報累積強降水落區(qū)偏大,漏報率高;CMA-SH9的面積比和縱橫比均表現(xiàn)較好,漏報率低,因此獲得較高的CSI和總評分。MODE能較好地判斷累積強降水目標的形態(tài),為預(yù)報員提供直觀的參考信息。
2.1.2 短時強降水預(yù)報鄰域檢驗
在高壓后部型降雨過程中,基于鄰域法,對短時強降水預(yù)報進行命中率、空報率評分檢驗。
為了分析不同預(yù)報時效多區(qū)域高分辨率模式短時強降水預(yù)報性能差異,將各模式08時和20時起報的產(chǎn)品進行整合檢驗(圖1),對24 h、36 h時效內(nèi)的短時強降水預(yù)報命中率、空報率評分進行分析,結(jié)果顯示:在高壓后部型降水過程中,三個模式所有時次的預(yù)報命中率評分都低于12%,最高值為11.6%;分析12 h短時強降水預(yù)報檢驗平均值(Avg_12)得出,在方圓40 km內(nèi)不計偏差情況下,模式預(yù)報命中率平均在4%以下,最大值達3.8%,CMA-MESO在36 h時效內(nèi)的預(yù)報命中率高于24 h,CMA-SH9則相反;空報率平均值為61%,三個模式均表現(xiàn)為臨近時刻短時強降水預(yù)報空報率高(24 h 預(yù)報空報率大于36 h)。
圖1 高壓后部型降水過程中多區(qū)域模式(a,c)24 h和(b,d)36 h的短時強降水預(yù)報檢驗評分(數(shù)字,陰影)(a,b)命中率,(c,d)空報率
以上結(jié)果表明:CMA-MESO在36 h的預(yù)報效果好于其他兩個區(qū)域模式(命中率高、空報率低);CMA-SH9在24 h的預(yù)報性能高于36 h;CMA-DB的24 h、36 h預(yù)報效果持平。多區(qū)域高分辨率模式短時強降水預(yù)報命中率較低、空報率較高的可能原因包括:遼寧(北方)地區(qū)短時強降水發(fā)生概率和頻次低于南方地區(qū),因此預(yù)報檢驗評分低于全國(唐文苑等,2017);模式本身存在一定誤差;空間分辨率決定三個模式均對產(chǎn)生短時強降水的中小尺度系統(tǒng)預(yù)報能力有限。
為了分析不同起報時刻多區(qū)域模式預(yù)報性能差異,分別對08時和20時起報的短時強降水預(yù)報的命中率、空報率評分進行單獨分析(圖略),結(jié)果顯示:與20時相比,CMA-SH9在08時起報的短時強降水命中率低、空報率高;CMA-MESO表現(xiàn)為08時的命中率和空報率均較高;CMA-DB與CMA-MESO相反。
CMA-DB模式的短時強降水預(yù)報檢驗結(jié)果顯示:預(yù)報命中率在起報時間后的13~24 h,其中間預(yù)報效果好于兩端,這與區(qū)域模式雷達回波逐小時預(yù)報效果隨時間的變化表現(xiàn)一致(劉靜等,2019)。
當(dāng)遼寧省受氣旋影響時,通常為偏北風(fēng)和東南風(fēng)輻合,冷空氣和暖濕空氣結(jié)合后,促進中尺度對流系統(tǒng)在地面鋒前暖區(qū)發(fā)展,產(chǎn)生累積強降水(孫欣等,2016)。
2.2.1 累積強降水預(yù)報目標對象檢驗
由表4可見重疊面積/聯(lián)合面積值為0.002~0.186(即0.2%~18.6%),結(jié)合面積比(實況面積/預(yù)報面積)值的0.146~0.301(即14.6%~30.1%),可以計算出,CMA-SH9、CMA-DB兩個模式預(yù)報的累積強降水落區(qū)大部分(81%~98%,而CMA-MESO僅有25%)落在實際強降水發(fā)生區(qū)之外,因此CMA-SH9和CMA-DB兩個模式預(yù)報的漏報率均很高,在66.8%~72.1%,而CMA-MESO模式的漏報率比其他兩個模式小得多,漏報率低于37%;模式空報率也很大,因為面積比指標顯示預(yù)報面積比實況面積大得多,因此CMA-SH9和CMA-DB兩個模式的CSI評分只有21.7%~30.6%。
表4 同表3,但為氣旋型降水過程
綜合得到評價結(jié)果:在氣旋型降水過程中,三個區(qū)域模式的質(zhì)心距離均在17~22 km,屬帶狀目標;CMA-MESO預(yù)報的質(zhì)心距離最近,但預(yù)報縱橫比和實況相差較多,對累積強降水落區(qū)的形態(tài)描述狹長,優(yōu)勢在于預(yù)報與實況重疊面積大,因此獲得較高的CSI評分和總評分。
2.2.2 短時強降水預(yù)報鄰域檢驗
在氣旋型降水過程中,將各模式08時和20時起報的產(chǎn)品進行整合檢驗(圖2),對24 h、36 h 的短時強降水預(yù)報命中率、空報率評分進行分析,結(jié)果顯示:在氣旋型降水過程中,三個模式所有時次的預(yù)報命中率評分都低于20%,最高值為19.7%;分析12 h 短時強降水預(yù)報檢驗平均值(Avg_12)得出,在方圓40 km內(nèi)不計偏差情況下,模式預(yù)報命中率平均在10%以下,最大值達9.2%;在24 h預(yù)報時效內(nèi),CMA-MESO預(yù)報命中率為CMA-SH9的兩倍;CMA-SH9在36 h時效內(nèi)的預(yù)報命中率高于24 h,CMA-DB則相反;空報率平均值為41%,CMA-SH9在24 h時效內(nèi)的預(yù)報空報率大于36 h,其他兩個模式則相反。
圖2 同圖1,但降水類型為氣旋型降水
以上結(jié)果表明:CMA-MESO的命中率大于其他兩個區(qū)域模式;CMA-SH9在24 h的預(yù)報效果差,即命中率(空報率)低(高),CMA-DB則相反。
在起報時刻后的第32個時次,CMA-SH9呈現(xiàn)出預(yù)報命中率接近空報率的特征,這與CMA-SH9對西南地區(qū)的累積強降水預(yù)報檢驗結(jié)論一致(徐同等,2019)。
三個區(qū)域高分辨率模式在氣旋型降水中的短時強降水預(yù)報命中率高于高壓后部型降水,客觀上表明多區(qū)域模式對氣旋型降水過程中的短時強降水預(yù)報效果較好。
為了分析不同起報時刻多區(qū)域高分辨率模式短時強降水預(yù)報性能差異,分別對08時和20時起報的短時強降水預(yù)報命中率、空報率評分進行單獨分析(圖略),結(jié)果顯示:與20時相比,CMA-SH9表現(xiàn)為08時起報的短時強降水命中率和空報率均較低,CMA-MESO和CMA-SH9預(yù)報表現(xiàn)相反;CMA-DB表現(xiàn)命中率低,空報率高。
遼寧省臺風(fēng)型降水多受臺風(fēng)倒槽影響,鋒生次級環(huán)流與風(fēng)場輻合作用增強了動力抬升,水汽達到暴雨閾值,不穩(wěn)定能量釋放產(chǎn)生大暴雨(閻琦等,2019)。
2.3.1 累積強降水預(yù)報目標對象檢驗
由表5可見重疊面積/聯(lián)合面積值為0.084~0.122(即8.4%~12.2%),結(jié)合面積比(實況面積/預(yù)報面積)的0.369~0.581(即36.9%~58.1%),可以總結(jié)出,三個模式預(yù)報的累積強降水落區(qū)幾乎(60%~76%)落在實際強降水發(fā)生區(qū)之外,這就可以合理解釋模式預(yù)報的漏報率為何那么高,均值在67%左右(模式空報率也很大,因為面積比指標顯示實況面積比預(yù)報面積小得多),CSI評分只有9.9%~12.4%。
表5 同表3,但為臺風(fēng)型降水過程
由于預(yù)報面積偏小、落區(qū)形態(tài)差異等問題,三個區(qū)域模式在臺風(fēng)型累積強降水過程中的CSI評分偏低,為預(yù)報員提供參考信息的能力不足。
2.3.2 短時強降水預(yù)報檢驗
在臺風(fēng)型降水類型中,為了分析不同預(yù)報時效多區(qū)域高分辨率模式短時強降水預(yù)報性能差異,將各模式08時和20時起報的產(chǎn)品進行整合檢驗(圖3),對24 h、36 h時效內(nèi)的短時強降水預(yù)報命中率、空報率評分進行分析,結(jié)果顯示:三個模式所有時次的預(yù)報命中率評分都低于0.11(11%);分析12 h 短時強降水預(yù)報檢驗平均值(Avg_12)得出,在方圓40 km內(nèi)不計偏差情況下,模式預(yù)報命中率最大值僅為3%,CMA-MESO在36 h時效內(nèi)的預(yù)報命中率好于24 h,CMA-SH9和CMA-DB則相反;空報率平均值為75%,CMA-DB在24 h時效內(nèi)的預(yù)報空報率高于36 h,其他兩個模式相反。
圖3 同圖1,但降水類型為臺風(fēng)型降水
以上結(jié)果表明:在24 h和36 h的整體時效內(nèi),CMA-SH9的命中率大于其他兩個區(qū)域模式,但空報率表現(xiàn)為三個區(qū)域模式最高值,平均值達75%,這與王新敏和栗晗(2020)研究結(jié)論一致,CMA-SH9模式對臺風(fēng)暴雨的降水范圍和強度比實況大得多,空報也較多;所有模式12 h短時強降水預(yù)報命中率均值小于2.1%。綜合考慮命中率和空報率,多區(qū)域高分辨率模式在臺風(fēng)型降水過程中的預(yù)報效果顯著低于高壓后部型和氣旋型降雨,這可能是因為,區(qū)域模式在長時效內(nèi)預(yù)報臺風(fēng)路徑的能力有限,對于臺風(fēng)系統(tǒng)影響陸上降水預(yù)報的穩(wěn)定性和準確性均表現(xiàn)欠佳(宮宇等,2018),也有可能因為區(qū)域模式對臺風(fēng)快速增強的預(yù)報能力較弱(薛文博等,2020)。
為了分析不同起報時刻多區(qū)域模式短時強降水預(yù)報性能差異,分別對08時和20時起報的短時強降水預(yù)報命中率、空報率評分進行單獨分析(圖略),結(jié)果顯示:與20時相比,CMA-SH9、CMA-MESO和CMA-DB均表現(xiàn)為08時起報的短時強降水預(yù)報效果好,即命中率(空報率)高(低)。
目前,天氣預(yù)報業(yè)務(wù)和多區(qū)域高分辨率模式研究常用長時效基于點對點的檢驗方法,面向多區(qū)域高分辨率模式,基于短時效和逐小時的強降水檢驗尚不多見。本文綜合已有的研究,對東北地區(qū)業(yè)務(wù)常用多區(qū)域高分辨率模式進行累積強降水目標對象檢驗和短時強降水鄰域法檢驗,從不同角度挖掘區(qū)域模式附加信息,為預(yù)報員決策提供參考。
(1)千米尺度或接近千米尺度的三個模式,在36 h時效內(nèi),對于累積強降水落區(qū)的形態(tài)預(yù)報與實況有相似性,落區(qū)質(zhì)心預(yù)報偏差一般在20 km左右。然而,強降水預(yù)報落區(qū)與實況重疊的面積一般都在10%以下,個別情形下(如CMA-MESO對于氣旋型降水過程),累積強降水落區(qū)預(yù)報與實況重疊度能夠接近20%;位置偏離的直接結(jié)果是導(dǎo)致漏報率高(一般在75%左右,CMA-MESO漏報率略低,為10%~20%),其中高壓后部型降水過程中累積強降水的漏報率超過80%,位置偏離也造成較高空報率。
(2)對短時強降水預(yù)報,在方圓40 km內(nèi)不計偏差情況下,模式預(yù)報命中率平均在10%以下、最大值為9.2%,氣旋型降水過程的命中率評分高于高壓后部型和臺風(fēng)型;空報率平均值為58.7%,臺風(fēng)型降水過程中空報率最高,達到74.6%;三種天氣型中,區(qū)域模式對臺風(fēng)型降水過程的短時強降水預(yù)報性能最低。
(3)從累積強降水目標對象檢驗得出:重疊面積與聯(lián)合面積比與CSI評分有很好的正相關(guān),預(yù)報的縱橫比小(高估累積強降水落區(qū)帶狀分布)、重疊面積與聯(lián)合面積比小、位置偏差大、漏報率高等因素均易造成CSI評分偏低;CMA-SH9對累積強降水落區(qū)形態(tài)帶狀分布估計過高,CMA-MESO累積強降水落區(qū)預(yù)報偏北,CMA-DB與實況目標匹配的預(yù)報目標個數(shù)最多。
(4)在高壓后部型降水過程中,應(yīng)充分考慮南側(cè)的累積強降水落區(qū);在氣旋型降水過程中,可以結(jié)合氣旋系統(tǒng)強弱,調(diào)整累積強降水落區(qū)的帶狀分布形態(tài);在臺風(fēng)型降水過程中,應(yīng)結(jié)合全球模式預(yù)報結(jié)論,預(yù)估多區(qū)域模式短時強降水預(yù)報面積,判斷累積強降水預(yù)報的形態(tài)和落區(qū)。
受時間尺度和降水日期數(shù)量的限制,檢驗結(jié)果不能完全準確地描述區(qū)域高分辨率模式預(yù)報的能力。不同區(qū)域模式在遼寧地區(qū)降水過程中預(yù)報表現(xiàn)差異較大,具體是由哪些因素導(dǎo)致了不同的預(yù)報效果有待進一步研究。
致謝:本文在撰寫過程中,得到上海市氣象臺戴建華研究員、陜西省氣象臺潘留杰研究員的悉心指導(dǎo)和幫助,在此表示衷心感謝!