CMA-SH9及ECMWF模式對(duì)龍巖市一次暴雨預(yù)報(bào)檢驗(yàn)分析*

曾凡輝, 倪 悅, 鄭敏煜

（福建省龍巖市氣象局, 福建 龍巖 364000）

0 引言

每年前汛期（4～6 月）暖區(qū)暴雨, 因其具有降水強(qiáng)度大、時(shí)效性短、局地性強(qiáng)等特點(diǎn)[1］, 極易造成洪澇災(zāi)害, 常常給人民生命財(cái)產(chǎn)和社會(huì)經(jīng)濟(jì)帶來(lái)巨大損失。目前, 暖區(qū)暴雨已成為福建省最主要的氣象災(zāi)害之一, 暖區(qū)暴雨由于其受復(fù)雜地形、下墊面的作用, 一直以來(lái)都是福建地區(qū)預(yù)報(bào)中的重點(diǎn)和難點(diǎn)之一。近年來(lái), 隨著中尺度數(shù)值模式的迅速發(fā)展, 人們對(duì)暴雨過(guò)程進(jìn)行了大量數(shù)值模擬研究, 加深了對(duì)暖區(qū)暴雨的認(rèn)識(shí)[2-4］。另外, 伴隨公眾對(duì)氣象服務(wù)產(chǎn)品要求的不斷細(xì)化, 數(shù)值預(yù)報(bào)產(chǎn)品尤其是高分辨率數(shù)值模式產(chǎn)品, 已經(jīng)成為目前也是今后一個(gè)較長(zhǎng)時(shí)期內(nèi)現(xiàn)代天氣預(yù)報(bào)技術(shù)的基礎(chǔ)[5］。目前, 模式降水預(yù)報(bào)產(chǎn)品已成為業(yè)務(wù)預(yù)報(bào)中重要的參考依據(jù), 但不同模式在不同地區(qū)的適用性仍存在著一定的差異性[6-8］。

目前, 諸多學(xué)者對(duì)模式降水的檢驗(yàn)工作進(jìn)行了大量研究, 并取得了豐碩的成果。潘留杰等[9］使用ECMWF 模式對(duì)東北半球氣象要素場(chǎng)預(yù)報(bào)能力的檢驗(yàn)指出, 模式對(duì)高層的預(yù)報(bào)能力優(yōu)于低層, 但模式對(duì)700 hPa風(fēng)場(chǎng)的預(yù)報(bào)存在顯著差異。范水勇等[10］采用客觀檢驗(yàn)方法對(duì)不同分辨率模式進(jìn)行研究表明, 高分辨率區(qū)域模式降水預(yù)報(bào)效果較好, 尤其對(duì)大量級(jí)降水預(yù)報(bào)。王雨[11］應(yīng)用TS（Threat Score）評(píng)分等統(tǒng)計(jì)方法研究了不同模式降水預(yù)報(bào)能力的差異。耿建軍[12］用客觀檢驗(yàn)方法對(duì)ECMWF、T639 數(shù)值模式降水預(yù)報(bào)的檢驗(yàn)表明, 數(shù)值模式對(duì)明顯大尺度系統(tǒng)影響的降水過(guò)程預(yù)報(bào)效果較好, ECMWF對(duì)該次過(guò)程的預(yù)報(bào)能力優(yōu)于T639。

文章基于ECMWF 模式及CMA-SH9 模式, 通過(guò)計(jì)算實(shí)況與各家模式的相關(guān)系數(shù)和均方根誤差等參量, 對(duì)2019年5月27日午后發(fā)生于福建省龍巖市境內(nèi)的一次短時(shí)強(qiáng)降水過(guò)程進(jìn)行模式檢驗(yàn)分析, 綜合評(píng)估適用于龍巖地區(qū)的模式產(chǎn)品, 從而進(jìn)一步提高對(duì)該地區(qū)模式預(yù)報(bào)產(chǎn)品的釋用能力, 為預(yù)報(bào)員的預(yù)報(bào)和決策服務(wù)提供更科學(xué)的參考。

1 天氣實(shí)況

2019 年5 月27 日下午, 福建省龍巖市境內(nèi)出現(xiàn)一次局地短時(shí)強(qiáng)降水過(guò)程, 根據(jù)區(qū)域自動(dòng)站的降水統(tǒng)計(jì), 27日08：00～28日08：00（北京時(shí), 下同）共計(jì)14個(gè)站點(diǎn)24 h累計(jì)雨量超過(guò)50 mm, 24 h累計(jì)雨量以長(zhǎng)汀縣長(zhǎng)汀一中站點(diǎn)97.0 mm為最大, 最大小時(shí)雨量達(dá)48.7 mm（長(zhǎng)汀一中站點(diǎn)）。從雨區(qū)空間分布來(lái)看（圖1）, 此次降水過(guò)程大致可劃分南北兩個(gè)強(qiáng)降水區(qū), 27 日18：00 北部開始受到上游回波移入影響, 出現(xiàn)超過(guò)40 mm/h的短時(shí)強(qiáng)降水, 降水持續(xù)長(zhǎng)達(dá)3 h, 且降水區(qū)域范圍較集中；南部降水則相對(duì)分散, 且持續(xù)時(shí)間較短（1～2 h）。此次短時(shí)強(qiáng)降水主要集中于長(zhǎng)汀縣城區(qū)境內(nèi), 造成長(zhǎng)汀縣城區(qū)多處民房被淹, 城區(qū)部分低洼路段短時(shí)積水, 最深約20 cm；共計(jì)13 名群眾被水圍困；部分鄉(xiāng)鎮(zhèn)民房出現(xiàn)溜方, 過(guò)程未造成人員傷亡。

圖1 2019年5月27日08：00～28日08：00累計(jì)降水雨量（單位：mm）分布

2 環(huán)流背景及物理量場(chǎng)分析

曾凡輝[13］對(duì)此次降水過(guò)程環(huán)流背景場(chǎng)進(jìn)行詳細(xì)分析, 此處僅作簡(jiǎn)要概述。此次過(guò)程中低緯多短槽波動(dòng), 東移高原槽與強(qiáng)盛南之槽在雨區(qū)上游相結(jié)合, 是誘發(fā)此次強(qiáng)降水過(guò)程的主要天氣因素之一。暴雨發(fā)生前, 副高穩(wěn)定維持于我國(guó)南海沿岸, 龍巖處副高北側(cè)邊緣, 受西南氣流影響, 穩(wěn)定的環(huán)流形勢(shì)有利于暖濕水汽向暴雨區(qū)上空的持續(xù)輸送, 為強(qiáng)降水的發(fā)生提供充足的水汽。27 日午后雨區(qū)上空大氣層結(jié)處于不穩(wěn)定狀態(tài), 而中低層受暖脊控制, 近地面有弱冷空氣擴(kuò)散配合, 有利于對(duì)流降水進(jìn)一步的觸發(fā)與維持。

對(duì)比分析地面區(qū)域自動(dòng)氣象站的風(fēng)場(chǎng)與降水隨時(shí)間的分布特征（圖2）, 發(fā)現(xiàn)近地面層輻合線與降水有較好的匹配關(guān)系, 輻合線的生成、發(fā)展和消散過(guò)程與此次降水過(guò)程有著密切的關(guān)系。27日14：00, 地面輻合線主要位于龍巖偏北位置, 受其影響長(zhǎng)汀境內(nèi)局部近地面風(fēng)場(chǎng)開始轉(zhuǎn)北風(fēng), 但大部分仍以南風(fēng)為主, 輻合線附近產(chǎn)生分散弱降水；15：00 龍巖市西南部地面風(fēng)速加大（局部超過(guò)8 m·s-1）, 地面輻合線在南部也有生成, 有利南部對(duì)流觸發(fā)并導(dǎo)致強(qiáng)降水的發(fā)生。與此同時(shí), 北部輻合線持續(xù)維持于長(zhǎng)汀境內(nèi)。南、北輻合線附近仍有零星弱降水的產(chǎn)生；17：00 西南部的地面輻合線東移, 對(duì)應(yīng)南部雨區(qū)整體有所東移。北部長(zhǎng)汀境內(nèi)的輻合線略有南壓（大致位于長(zhǎng)汀、連城、武平三縣交界處）, 輻合線北側(cè)北風(fēng)有所加大（局部超過(guò)8 m·s-1）, 長(zhǎng)汀主要降雨區(qū)處于風(fēng)速輻合區(qū)內(nèi), 有利于北部強(qiáng)降水的觸發(fā)與維持, 長(zhǎng)汀縣局部（城區(qū)）17:00 降雨達(dá)大雨量級(jí)；至27日22：00, 南、北輻合線完全東移出龍巖市境內(nèi), 降水也趨于減弱并結(jié)束。

圖2 龍巖市區(qū)域自動(dòng)站極大風(fēng)向風(fēng)速（單位：m·s-1）和小時(shí)雨量演變

3 數(shù)值預(yù)報(bào)應(yīng)用情況與檢驗(yàn)分析

3.1 檢驗(yàn)要素

環(huán)流形勢(shì)包括U和V風(fēng)場(chǎng)（925、850、700、500、200 hPa, 簡(jiǎn)稱U925、U850、U700、U500、U200, V 同理）、500 hPa高度場(chǎng)（簡(jiǎn)稱HGT500）；物理量要素包括2 m氣溫（簡(jiǎn)稱2t）、溫 度（925、850、700、500、200 hPa, 簡(jiǎn) 稱tt925、tt850、tt700、tt500、tt200）、K指數(shù)、比濕（925、850、700 hPa, 簡(jiǎn)稱q925、q850、q700）和相對(duì)濕度（925、850、700 hPa, 簡(jiǎn)稱rh925、rh850、rh700）。要素檢驗(yàn)范圍為（24°～26.5°N, 115.5°～118°E）。環(huán)流形勢(shì)和物理量要素檢驗(yàn)時(shí)次為2019 年5 月26 日～27 日08：00、20：00 起報(bào)的27～28 日的要素場(chǎng), 降水檢驗(yàn)時(shí)段為2019 年5 月27 日08：00～28日08：00（BT）。

3.2 檢驗(yàn)結(jié)果

3.2.1 降水檢驗(yàn)

從實(shí)況降水空間分布來(lái)看（圖1）, 2019年5月27日08：00～28 日08：00, 龍巖市西北、西南及中東部實(shí)況均出現(xiàn)大雨以上量級(jí)降水, 長(zhǎng)汀與新羅的北部、上杭南部等部分區(qū)域達(dá)暴雨量級(jí)；全球模式ECWMF在5月27日08：00起報(bào)的24 h預(yù)報(bào)降水強(qiáng)度偏弱, 以中等強(qiáng)度降水為主, 西北部長(zhǎng)汀縣局部有大雨, 對(duì)長(zhǎng)汀的暴雨降水存在漏報(bào), 且未準(zhǔn)確預(yù)報(bào)出南部大雨及以上量級(jí)的強(qiáng)降水, 強(qiáng)降水中心整體偏西, 降水落區(qū)形態(tài)與實(shí)況相差較大；與ECWMF 模式相比, 國(guó)產(chǎn)中尺度數(shù)值模式CMASH9 在相同起報(bào)時(shí)間內(nèi), 對(duì)24 h 降水量級(jí)和降水落區(qū)形態(tài)的預(yù)報(bào)較為準(zhǔn)確。與降水實(shí)況相比較, CMA-SH9對(duì)于強(qiáng)降水落區(qū)及強(qiáng)度的預(yù)報(bào)上仍存在一定偏差, 長(zhǎng)汀暴雨中心預(yù)報(bào)位置偏北、降水強(qiáng)度偏強(qiáng)（24 h預(yù)報(bào)降水超過(guò)100 mm, 達(dá)大暴雨量級(jí), 但實(shí)況為暴雨量級(jí), 累計(jì)雨量最大為長(zhǎng)汀一中站點(diǎn)99 mm）。同時(shí), CMA-SH9模式對(duì)較分散的南部降水區(qū)預(yù)報(bào)強(qiáng)降水范圍偏大、強(qiáng)度偏強(qiáng)。綜上分析, ECWMF 與CMA-SH9 模式在對(duì)強(qiáng)降水中心位置及形態(tài)的預(yù)報(bào)上均存在偏差, 相較而言, CMA-SH9 模式的降水預(yù)報(bào)產(chǎn)品對(duì)此次強(qiáng)降水落區(qū)和強(qiáng)度預(yù)報(bào)具有較好指導(dǎo)意義。

3.2.2 環(huán)流形勢(shì)場(chǎng)檢驗(yàn)

從各預(yù)報(bào)時(shí)效檢驗(yàn)來(lái)看, ECWMF 模式預(yù)報(bào)產(chǎn)品500 hPa高度場(chǎng)的均方根誤差明顯低于CMA-SH9模式（圖3a）, 表明ECWMF對(duì)高層的環(huán)流形勢(shì), 尤其是500 hPa高度場(chǎng)預(yù)報(bào)效果要優(yōu)于CMA-SH9模式產(chǎn)品, 其對(duì)于大尺度的高空環(huán)流形勢(shì)具有較好的預(yù)報(bào)參考性。對(duì)于不同高度（200、500、700、850 和925 hPa, 下同）風(fēng)場(chǎng)U 各時(shí)效均方根誤差, ECWMF 模式明顯大于CMA-SH9 模式, 且ECWMF模式產(chǎn)品的預(yù)報(bào)誤差隨高度增加而呈現(xiàn)處出非線性的變化, 925 hPa與500 hPa 的預(yù)報(bào)誤差最小, 850 hPa 的預(yù)報(bào)誤差最大；CMA-SH9模式對(duì)于不同高度（200、500、700、850和925 hPa）上的U風(fēng)場(chǎng)預(yù)報(bào)誤差則相對(duì)比較接近, 且誤差明顯低于ECMWF。ECWMF 各預(yù)報(bào)時(shí)效200 hPa、500 hPa和700 hPa 高度風(fēng)場(chǎng)V 的預(yù)報(bào)誤差大于CMA-SH9, 其余高度風(fēng)場(chǎng)V 預(yù)報(bào)誤差與CMA-SH9 持平, ECWMF 對(duì)于V 風(fēng)場(chǎng)的預(yù)報(bào)誤差隨高度增加而明顯增大, 但是CMA-SH9模式的預(yù)報(bào)誤差隨高度的變化不明顯。

ECWMF模式對(duì)于各時(shí)效500 hPa高度預(yù)報(bào)的相關(guān)系數(shù)（圖3b）要明顯高于CMA-SH9, 200 hPa 層結(jié)高度情況相反, 其余高度場(chǎng)上ECWMF 為負(fù)相關(guān), 可參考性較差。ECWMF模式在各預(yù)報(bào)時(shí)效內(nèi)200 hPa和500 hPa高度場(chǎng)上對(duì)于U 參量預(yù)報(bào)風(fēng)場(chǎng)均為正相關(guān)（圖3b）, 且500 hPa的U風(fēng)場(chǎng)預(yù)報(bào)的相關(guān)性要明顯高于CMA-SH9, 但200 hPa 的情況則相反, 其他高度上（700、850 和925 hPa）U風(fēng)場(chǎng)預(yù)報(bào)均為負(fù)相關(guān), 表明在該范圍高度內(nèi)ECMWF 模式對(duì)U 風(fēng)場(chǎng)預(yù)報(bào)的可參考性較差, 而CMASH9 模式除了700 hPa 高度的U 風(fēng)場(chǎng)預(yù)報(bào)為負(fù)相關(guān)外, 其余高度均為正相關(guān), 且相關(guān)系數(shù)絕對(duì)值要明顯高于ECWMF產(chǎn)品。

圖3 模式對(duì)環(huán)流形勢(shì)場(chǎng)12～36 h時(shí)效預(yù)報(bào)檢驗(yàn)

由圖2可知, 近地層地面輻合線是此次降水過(guò)程的主要影響系統(tǒng)之一, 因此對(duì)近地層風(fēng)場(chǎng)的預(yù)報(bào)效果將直接對(duì)降水預(yù)報(bào)效果產(chǎn)生明顯的影響。參考上述對(duì)ECMWF和CMA-SH9兩種模式對(duì)風(fēng)場(chǎng)預(yù)報(bào)的檢驗(yàn)分析可知, 基于預(yù)報(bào)誤差和相關(guān)性等參量, 國(guó)產(chǎn)CMA-SH9模式對(duì)于近地面層風(fēng)場(chǎng)（925 hPa）的預(yù)報(bào)效果明顯優(yōu)于ECMWF, 因此, 此次過(guò)程CMA-SH9模式對(duì)于對(duì)24 h地面降水預(yù)報(bào)效果優(yōu)于ECMWF, 這與上文降水檢驗(yàn)部分的結(jié)論相一致。

3.2.3 物理量檢驗(yàn)

ECWMF 模式對(duì)于850 hPa和925 hPa 高度上溫度場(chǎng)的預(yù)報(bào)均方根誤差（圖4a）明顯小于CMA-SH9（其中, 24 h預(yù)報(bào)時(shí)效內(nèi), 二者均方根誤差的差值在0.7 ℃～1 ℃）, 其余高度CMA-SH9 的均方根誤差更小, 且兩個(gè)模式預(yù)報(bào)的均方根誤差差值小于0.2 ℃, 差距不明顯。ECWMF對(duì)24 h時(shí)效, 溫度場(chǎng)預(yù)報(bào)的相關(guān)性（圖5b）略高于CMA-SH9, 且均為正相關(guān), 除700 hPa高度上溫度的相關(guān)性較低之外, 其余高度上的相關(guān)性均較大, 超過(guò)0.6。相較而言, 此次降水過(guò)程CMA-SH9 對(duì)于各高度場(chǎng)溫度的預(yù)報(bào)效果則較差。此外, 從濕度條件檢驗(yàn)情況來(lái)看, ECWMF各預(yù)報(bào)時(shí)效比濕的均方根誤差隨高度呈非線性變化, 先增加后減小, 700 hPa和925 hPa預(yù)報(bào)誤差大小接近, 相對(duì)ECMWF而言, 850 hPa的誤差較大。

由圖4 可知, 除700 hPa 外, ECMWF 預(yù)報(bào)比濕的相關(guān)性隨高度增加而有所提高；ECWMF模式各時(shí)效預(yù)報(bào)925 hPa 和850 hPa 相對(duì)濕度的相關(guān)系數(shù)要明顯高于CMA-SH9 模式預(yù)報(bào)產(chǎn)品。24 h 時(shí)效內(nèi), ECMWF 預(yù)報(bào)低層（850 hPa 和925 hPa）的相對(duì)濕度相關(guān)性明顯大于中高層（700 hPa）, CMA-SH9 相對(duì)濕度則以925 hPa 的相關(guān)性為最大, 中上層的相關(guān)性則較低。這表明兩家模式對(duì)于近地層濕度條件的預(yù)報(bào)效果均要好于中高層。此外, ECWMF 模式對(duì)k 指數(shù)預(yù)報(bào)整體效果較好, 均方根誤差及相關(guān)性方面均占優(yōu), 表明ECMWF對(duì)于熱力參數(shù)的預(yù)報(bào)效果要優(yōu)于CMA-SH9, 但二者差距不顯著。

圖4 模式對(duì)物理量12～36 h時(shí)效預(yù)報(bào)檢驗(yàn)

4 結(jié)論與討論

文章基于全球ECWMF 模式及過(guò)程中尺度CMASH9 模式0～24 h 時(shí)效的24 h 降水預(yù)報(bào)資料及同時(shí)段自動(dòng)站觀測(cè)降水資料, 對(duì)2019 年5 月27 日發(fā)生于龍巖境內(nèi)的一次短時(shí)強(qiáng)降水過(guò)程進(jìn)行檢驗(yàn)分析, 評(píng)估兩家模式對(duì)此次過(guò)程的預(yù)報(bào)能力, 得出主要結(jié)論如下：

（1）此次短時(shí)強(qiáng)降水過(guò)程主要影響系統(tǒng)為低空急流和近地層輻合線, 暴雨區(qū)主要位于低空急流右前側(cè)及地面輻合線附近。有利的環(huán)流背景條件為暴雨區(qū)上空水汽的持續(xù)輸送提供了條件, 與此同時(shí), 低層輻合和高層輻散的動(dòng)力條件配置, 有利于水汽和熱量的垂直輸送與不穩(wěn)定能量的釋放, 是觸發(fā)此次短時(shí)強(qiáng)降水非常重要的因素之一。

（2）此次強(qiáng)降水過(guò)程中, ECWMF模式和CMA-SH9模式在對(duì)于強(qiáng)降水中心位置及其形態(tài)的預(yù)報(bào)上存在著明顯的偏差, 總體而言, ECWMF 模式對(duì)于中等強(qiáng)度量級(jí)降水空間分布的預(yù)報(bào)效果要優(yōu)于CMA-SH9模式, 但對(duì)于大雨及以上量級(jí)的降水, 其預(yù)報(bào)效果則要明顯低于CMA-SH9 模式。相較而言, CMA-SH9 模式的降水預(yù)報(bào)產(chǎn)品對(duì)此次強(qiáng)降水落區(qū)和強(qiáng)度預(yù)報(bào)具有較好指導(dǎo)意義。

（3）ECMWF 和CMA-SH9 兩種模式在此次降水過(guò)程中, 二者對(duì)于物理量場(chǎng)的預(yù)報(bào)偏差不是特別顯著, 整體而言, ECWMFF 模式對(duì)于溫度、濕度、熱力等物理量場(chǎng)的預(yù)報(bào)效果要略優(yōu)于CMA-SH9模式。

（4）此次降水過(guò)程中, CMA-SH9 模式對(duì)于近地面層, 尤其是925 hPa風(fēng)場(chǎng)預(yù)報(bào)效果要明顯優(yōu)于ECMWF, 但高層環(huán)流形勢(shì)場(chǎng)則以ECMWF 預(yù)報(bào)更具有參考性。由于此次過(guò)程, 主要天氣系統(tǒng)為近地層輻合線, CMASH9 模式對(duì)于低層輻合線的預(yù)報(bào)效果要優(yōu)于ECMWF。綜合評(píng)估此次降水過(guò)程, 相較于ECMWF, CMA-SH9 模式對(duì)此次強(qiáng)降水整體落區(qū)和強(qiáng)度預(yù)報(bào)具有較好指導(dǎo)意義。