蘇 翔 康志明 莊瀟然 陳圣劼

1 江蘇省氣象臺(tái)，南京 210008 2 中國(guó)氣象局交通氣象重點(diǎn)開(kāi)放實(shí)驗(yàn)室，南京 210009

提 要： 針對(duì)ECMWF模式在2020年我國(guó)江淮流域超長(zhǎng)梅雨期暴雨雨帶預(yù)報(bào)中的不確定性，選取10個(gè)典型的狹長(zhǎng)暴雨雨帶，基于對(duì)象診斷方法(MODE)對(duì)雨帶東西段的位置預(yù)報(bào)偏差、穩(wěn)定性以及偏差的連續(xù)性等進(jìn)行分析，并基于集合敏感性方法(ESA)分析揭示雨帶典型偏差特征的關(guān)鍵天氣系統(tǒng)。結(jié)果表明：雨帶東段的預(yù)報(bào)不確定性總體高于西段，尤其是到了中期預(yù)報(bào)時(shí)效，預(yù)報(bào)不確定性進(jìn)一步加強(qiáng)。短期時(shí)效內(nèi)，雨帶東西段均存在偏北的系統(tǒng)偏差，雨帶西段的漏報(bào)率低于東段，穩(wěn)定性和連續(xù)性均好于東段；中期時(shí)效內(nèi)，雨帶東段的預(yù)報(bào)誤差增長(zhǎng)較為明顯，且緯度預(yù)報(bào)平均變幅也高于西段。集合敏感性分析結(jié)果揭示，500 hPa高空槽、副熱帶高壓、850 hPa低空急流和切變線(xiàn)的預(yù)報(bào)偏差及它們與降水潛熱之間存在的非自然正反饋共同作用造成雨帶東段預(yù)報(bào)偏北，其中后半夜高空槽和低空急流的發(fā)展加強(qiáng)是出現(xiàn)梅雨鋒暴雨雨帶東段北偏的背景場(chǎng)特征之一。

引 言

受東亞夏季風(fēng)影響，我國(guó)江淮流域地區(qū)夏季暴雨災(zāi)害頻發(fā)(丁一匯，2019)。2020年我國(guó)該地區(qū)出現(xiàn)了持續(xù)時(shí)間長(zhǎng)達(dá)62天的超長(zhǎng)梅雨期，平均降水量達(dá)759.2 mm，為1961年以來(lái)最多。這次梅雨過(guò)程具有入梅時(shí)間早(6月1日)、出梅時(shí)間晚(8月2日)、覆蓋范圍大、極端降水事件頻發(fā)等特點(diǎn)，為各地的氣象預(yù)報(bào)服務(wù)帶來(lái)不小的壓力和挑戰(zhàn)。

針對(duì)梅雨鋒暴雨過(guò)程，許多學(xué)者已通過(guò)數(shù)值模擬和診斷分析等方式研究了其動(dòng)力學(xué)特征和發(fā)生發(fā)展機(jī)制。劉建勇等(2011)通過(guò)數(shù)值試驗(yàn)，從初始誤差增長(zhǎng)的角度闡釋了梅雨期暴雨系統(tǒng)的可預(yù)報(bào)性。尹韓笑等(2019)利用WRF模式和NCEP再分析資料研究了一次梅雨暴雨過(guò)程的動(dòng)能譜特征并診斷了動(dòng)能收支方程。閔錦忠等(2018)使用NCEP/NCAR再分析資料、衛(wèi)星和自動(dòng)站等觀測(cè)資料，結(jié)合WRF模式，對(duì)一次梅雨鋒暴雨過(guò)程的發(fā)生發(fā)展機(jī)制進(jìn)行了診斷和模擬。蔣薇和高輝(2013)利用觀測(cè)資料統(tǒng)計(jì)了1961—2011年長(zhǎng)江中下游梅雨各特征量的時(shí)空變化特征。趙曉琳(2019)基于NCEP/NCAR再分析資料、NOAA OLR資料和降水站點(diǎn)觀測(cè)資料分析揭示了2018年我國(guó)梅雨期降水異?，F(xiàn)象及其成因。金瓊等(2020)則從環(huán)流因子特征的角度對(duì)該現(xiàn)象展開(kāi)了解釋。上述研究多偏向于機(jī)理分析，而在實(shí)際業(yè)務(wù)中，預(yù)報(bào)員面臨的最大挑戰(zhàn)是梅雨期暴雨雨帶的位置預(yù)報(bào)把握不準(zhǔn)的問(wèn)題，相關(guān)研究仍相對(duì)缺乏。

在2020年梅雨期暴雨預(yù)報(bào)實(shí)踐中，預(yù)報(bào)員們發(fā)現(xiàn)ECMWF模式對(duì)一些狹長(zhǎng)暴雨雨帶的預(yù)報(bào)出現(xiàn)了東段明顯偏北而西段預(yù)報(bào)卻與實(shí)況較為一致的現(xiàn)象。以往的研究大多將梅雨雨帶作為整體進(jìn)行分析(王勇等，2012；陳艷麗等，2016)，沒(méi)有考慮雨帶東、西段位置預(yù)報(bào)偏差的差異性。梅雨鋒東西跨度大，在實(shí)際預(yù)報(bào)中不同位置梅雨暴雨常呈現(xiàn)出不同的預(yù)報(bào)誤差特征，尤其是雨帶東段位于亞歐大陸東側(cè)，同時(shí)受到西風(fēng)帶冷空氣、西南季風(fēng)以及副熱帶高壓(以下簡(jiǎn)稱(chēng)副高)等多個(gè)天氣系統(tǒng)綜合影響，經(jīng)常呈現(xiàn)出預(yù)報(bào)誤差大、預(yù)報(bào)不穩(wěn)定等特點(diǎn)，因此有必要將梅雨鋒暴雨的狹長(zhǎng)雨帶分割為東西兩段分別進(jìn)行研究。近年來(lái)，隨著空間檢驗(yàn)技術(shù)的不斷發(fā)展，基于對(duì)象的診斷評(píng)估方法(MODE)(Davis et al，2006;2009)已成為雨帶形態(tài)屬性和位置預(yù)報(bào)分析的較好工具(尤鳳春等，2011；劉湊華和牛若蕓，2013；茅懋等，2016；陳笑等，2018；蘇翔和康志明，2020)。另一方面，集合預(yù)報(bào)產(chǎn)品包含了預(yù)報(bào)誤差隨天氣形勢(shì)的變化(陳偉斌等，2019)，可用于對(duì)梅雨鋒暴雨預(yù)報(bào)中存在的不確定性信息進(jìn)行分析。集合敏感性分析(ESA)方法早期由Ancell and Hakim(2007)提出，是一種基于線(xiàn)性相關(guān)定量評(píng)估預(yù)報(bào)量對(duì)前期狀態(tài)量敏感性的方法，不僅有助于揭示影響預(yù)報(bào)對(duì)象的動(dòng)力學(xué)特征，還能夠增加對(duì)于預(yù)報(bào)誤差及傳播機(jī)制的認(rèn)知。近年來(lái)，ESA技術(shù)已被應(yīng)用于諸如熱帶氣旋(Xie et al，2013；Keller，2017)、暴雨(Lynch and Schumacher，2014；Yu and Meng，2016；Du and Chen，2018)、強(qiáng)對(duì)流(Hill et al，2016；Limpert and Houston，2018)等不同尺度天氣系統(tǒng)的可預(yù)報(bào)性研究中。

針對(duì)ECMWF模式在2020年梅雨期暴雨狹長(zhǎng)雨帶位置預(yù)報(bào)上存在偏差的問(wèn)題，本文首先基于MODE技術(shù)分別對(duì)雨帶東、西段位置預(yù)報(bào)的偏差特征展開(kāi)統(tǒng)計(jì)分析，并就不同時(shí)效上雨帶東、西段位置預(yù)報(bào)的穩(wěn)定性和持續(xù)性降水過(guò)程中雨帶東、西段位置預(yù)報(bào)偏差的連續(xù)性進(jìn)行評(píng)估，在此基礎(chǔ)上，結(jié)合ESA技術(shù)研究揭示雨帶典型偏差特征的預(yù)報(bào)不確定性來(lái)源，為提升梅雨期暴雨預(yù)報(bào)準(zhǔn)確度提供科學(xué)參考。

1 資料和方法

1.1 資 料

所使用的預(yù)報(bào)數(shù)據(jù)為ECMWF模式確定性預(yù)報(bào)和集合預(yù)報(bào)數(shù)據(jù)，分別用于2020年梅雨期暴雨雨帶的預(yù)報(bào)偏差分析和集合敏感性分析。ECMWF模式確定性預(yù)報(bào)的水平分辨率為0.125°，選取12 時(shí)(世界時(shí))起報(bào)的36～228 h(間隔24 h)時(shí)效的24 h累積降水預(yù)報(bào)進(jìn)行分析。ECMWF集合預(yù)報(bào)資料的水平分辨率為0.5°，包含1個(gè)控制預(yù)報(bào)和50個(gè)擾動(dòng)成員。為保障研究中暴雨雨帶的完整性，對(duì)中國(guó)氣象局2014年印發(fā)的梅雨監(jiān)測(cè)有關(guān)業(yè)務(wù)規(guī)定中劃分的梅雨區(qū)域范圍略做擴(kuò)大調(diào)整，選取25°～38°N、105°～123°E的2020年梅雨期暴雨雨帶進(jìn)行研究。所使用的實(shí)況資料為全國(guó)2 411個(gè)氣象站的逐日(08時(shí)至次日08時(shí)，北京時(shí)，下同)降水?dāng)?shù)據(jù)，通過(guò)反距離權(quán)重插值到與預(yù)報(bào)數(shù)據(jù)相同的0.125°網(wǎng)格上。為了研究不同格點(diǎn)降水實(shí)況資料對(duì)于本文的研究結(jié)果是否敏感，還測(cè)試了逐小時(shí)0.1°分辨率的全國(guó)自動(dòng)氣象站與CMORPH融合降水產(chǎn)品和逐日0.05°分辨率的中國(guó)氣象局下發(fā)格點(diǎn)降水實(shí)況分析產(chǎn)品，發(fā)現(xiàn)本文插值所得的格點(diǎn)降水實(shí)況與這些格點(diǎn)降水分析產(chǎn)品雖存在細(xì)微的差別，但并不影響本文的研究結(jié)果和結(jié)論，這是由于本文的研究對(duì)象都是大尺度暴雨雨帶的緣故(見(jiàn)1.2節(jié))。

1.2 基于對(duì)象的診斷評(píng)估方法(MODE)

MODE可根據(jù)設(shè)定的目標(biāo)閾值識(shí)別出預(yù)報(bào)和觀測(cè)降水對(duì)象，并計(jì)算對(duì)象的面積、中心位置、長(zhǎng)軸、短軸、軸角、曲率等屬性，具體內(nèi)容可參考Davis et al(2006;2009)的研究。MODE分析出的降水對(duì)象的位移偏差和強(qiáng)度偏差可在一定程度上分別代表來(lái)源于天氣系統(tǒng)的空間偏差和剔除天氣系統(tǒng)偏差后的降水參數(shù)化或物理方案導(dǎo)致的偏差。設(shè)定的暴雨目標(biāo)閾值為50 mm·d-1?？紤]到觀測(cè)對(duì)象不隨預(yù)報(bào)時(shí)效變化，故以觀測(cè)對(duì)象為基準(zhǔn)，在2020年江淮流域梅雨期(6月1日至8月1日)選取滿(mǎn)足如下條件的狹長(zhǎng)暴雨樣本：

(1)對(duì)象面積≥60 000 km2；

(2)對(duì)象中心經(jīng)緯度范圍：28°～34°N、110°～123°E；

(3)對(duì)象縱橫比(長(zhǎng)軸與短軸的比例)≥3.0。

通過(guò)以上條件篩選出10個(gè)狹長(zhǎng)暴雨樣本(表1)，基本涵蓋2020年梅雨期所有典型的狹長(zhǎng)暴雨過(guò)程。梅雨期暴雨雨帶大多呈東西帶狀分布，預(yù)報(bào)的重點(diǎn)是雨帶的南北位置(緯度)。將篩選的10個(gè)狹長(zhǎng)暴雨樣本和與之匹配的不同預(yù)報(bào)時(shí)效的暴雨預(yù)報(bào)對(duì)象，按對(duì)象的中心經(jīng)度進(jìn)行經(jīng)向切割，并重新進(jìn)行MODE計(jì)算，獲得雨帶東段和西段的中心位置信息。切割后，對(duì)于預(yù)報(bào)和觀測(cè)對(duì)象面積比(預(yù)報(bào)和觀測(cè)對(duì)象面積較小者與較大者的比值)低于1/3的樣本，做漏報(bào)處理。這里將基于MODE的暴雨雨帶漏報(bào)率定義為：不存在與之匹配的預(yù)報(bào)暴雨雨帶的觀測(cè)暴雨雨帶對(duì)象數(shù)占總觀測(cè)暴雨雨帶對(duì)象數(shù)的比例。

表1 根據(jù)設(shè)定條件篩選出的2020年 梅雨期10個(gè)狹長(zhǎng)暴雨樣本Table 1 10 long and narrow heavy rain belts selected based on the set criteria during 2020 Meiyu period

1.3 預(yù)報(bào)的穩(wěn)定性

預(yù)報(bào)的穩(wěn)定性是用于衡量模式預(yù)報(bào)量隨預(yù)報(bào)時(shí)效變化是否一致的指標(biāo)，若變化不一致，說(shuō)明預(yù)報(bào)存在較大的不確定性。若相鄰兩個(gè)預(yù)報(bào)時(shí)效的預(yù)報(bào)量變幅不超過(guò)特定閾值，則稱(chēng)這兩個(gè)預(yù)報(bào)時(shí)效之間的預(yù)報(bào)是穩(wěn)定的，反之則是不穩(wěn)定的(其中一個(gè)預(yù)報(bào)時(shí)效出現(xiàn)漏報(bào)的情況也算不穩(wěn)定)。本研究中預(yù)報(bào)量選取為梅雨期暴雨雨帶東、西段預(yù)報(bào)對(duì)象中心的緯度值，變幅閾值為0.5°。圖1a、1b、1c分別給出了36 h、60 h、84 h時(shí)效對(duì)應(yīng)的7月6日08時(shí)至7日08時(shí)暴雨預(yù)報(bào)對(duì)象，用于示例預(yù)報(bào)的穩(wěn)定性。

1.4 預(yù)報(bào)偏差的連續(xù)性

預(yù)報(bào)偏差的連續(xù)性是用于評(píng)價(jià)持續(xù)性降水過(guò)程中預(yù)報(bào)偏差隨時(shí)間變化連續(xù)性的指標(biāo)，若變化連續(xù)，說(shuō)明前期的預(yù)報(bào)偏差信息能夠?yàn)楹笃诘念A(yù)報(bào)訂正提供參考。若相鄰兩個(gè)時(shí)次的預(yù)報(bào)偏差變幅不超過(guò)特定閾值，則稱(chēng)這兩個(gè)時(shí)次之間的預(yù)報(bào)偏差是連續(xù)的，反之則是不連續(xù)的(其中一個(gè)時(shí)次出現(xiàn)漏報(bào)的情況也算不連續(xù))。選取的預(yù)報(bào)量是梅雨期暴雨雨帶東、西段對(duì)象中心的緯度預(yù)報(bào)偏差，設(shè)定的緯度預(yù)報(bào)偏差變幅閾值為0.5°。圖1d、1e、1f分別給出了60 h時(shí)效的7月5日08時(shí)至6日08時(shí)、7月6日08時(shí)至7日08時(shí)、7月7日08時(shí)至8日08時(shí)暴雨預(yù)報(bào)和觀測(cè)對(duì)象重疊圖，用于示例預(yù)報(bào)偏差的連續(xù)性。

圖1 以(a)36 h、(b)60 h、(c)84 h預(yù)報(bào)時(shí)效對(duì)應(yīng)的2020年7月6日08時(shí)至7日08時(shí)暴雨預(yù)報(bào)對(duì)象(藍(lán)色) 示例預(yù)報(bào)的穩(wěn)定性概念，和以60 h預(yù)報(bào)時(shí)效對(duì)應(yīng)的(d)7月5日08時(shí)至6日08時(shí)、(e)7月6日08時(shí) 至7日08時(shí)、(f)7月7日08時(shí)至8日08時(shí)暴雨預(yù)報(bào)(藍(lán)色)和觀測(cè)(橙色)對(duì)象示例預(yù)報(bào)偏差的連續(xù)性概念Fig.1 The concept of forecast ability illustrated by using the (a) 36 h, (b) 60 h, (c) 84 h forecast of the heavy rain (blue) from 08:00 BT 6 to 08:00 BT 7 July 2020, and the concept of continuity of forecast bias illustrated by using the 60 h forecast of the heavy rain (blue) and corresponding observation (orange) from (d) 08:00 BT 5 to 08:00 BT 6, (e) 08:00 BT 6 to 08:00 BT 7 and (f) 08:00 BT 7 to 08:00 BT 8 July 2020

1.5 集合敏感性分析方法(ESA)

根據(jù)ESA的定義，假設(shè)集合預(yù)報(bào)在某個(gè)預(yù)報(bào)時(shí)刻的預(yù)報(bào)狀態(tài)量與后期某預(yù)報(bào)量之間存在線(xiàn)性關(guān)聯(lián)，則可通過(guò)建立預(yù)報(bào)量(J)與狀態(tài)量(xt)之間的線(xiàn)性回歸方程用以表征二者的關(guān)聯(lián)性：

(1)

2 預(yù)報(bào)偏差分析

2.1 雨帶東段和西段位置預(yù)報(bào)偏差分析

圖2展示了2020年梅雨期暴雨雨帶東、西段緯度預(yù)報(bào)平均誤差和漏報(bào)率。從緯度預(yù)報(bào)平均誤差來(lái)看，除了雨帶西段在中期時(shí)效前半段(108～156 h)平均誤差有所降低以外，平均誤差整體上隨預(yù)報(bào)時(shí)效的增加不斷上升。在短期時(shí)效內(nèi)(36～84 h)，雨帶東、西段的緯度預(yù)報(bào)誤差大小相當(dāng)(平均為0.6°)，而在36 h時(shí)效下僅為0.4°左右。而到了中期時(shí)效(108～228 h)，雨帶東段的緯度預(yù)報(bào)誤差明顯高于西段，尤其是中期時(shí)效前半段(108～156 h)，東、西段的緯度預(yù)報(bào)誤差分別為1.2°和0.7°。從漏報(bào)率上看，短期時(shí)效內(nèi)雨帶西段的漏報(bào)率(10%以?xún)?nèi))相比于東段更低(20%以?xún)?nèi))。到了中期時(shí)效后半段(108～228 h)，雨帶東段和西段的漏報(bào)率相對(duì)較高，分別在30%以上和40%以上。

圖3展示了2020年梅雨期暴雨雨帶東、西段預(yù)報(bào)偏北和偏南的比例(比例的分母已去除漏報(bào)樣本)。這里僅統(tǒng)計(jì)緯度預(yù)報(bào)偏差達(dá)0.25°(2倍格距)以上的樣本，并認(rèn)為緯度預(yù)報(bào)偏差在0.25°以下(不超過(guò)1倍格距)的樣本沒(méi)有南北預(yù)報(bào)偏差。可以看出，在短期時(shí)效內(nèi)，雨帶東、西段偏北的比例(60%～90%)遠(yuǎn)高于偏南的比例(0%～20%)，偏北的系統(tǒng)偏差非常明顯。為了進(jìn)一步考察預(yù)報(bào)偏北和偏南的雨帶在環(huán)流形勢(shì)上的差異，選取一天四次的FNL再分析數(shù)據(jù)(分辨率為0.25°)對(duì)偏北個(gè)例與偏南個(gè)例分別進(jìn)行合成分析。由圖4可見(jiàn)，梅雨帶預(yù)報(bào)位置偏北和偏南時(shí)的環(huán)流背景存在著顯著差異。在偏北型中(圖4a)，對(duì)流層中低層西南暖濕氣流較強(qiáng)，850 hPa 低空急流出口位于江淮中部地區(qū)，副高北抬至江淮地區(qū)南部，強(qiáng)度較強(qiáng)，受此形勢(shì)場(chǎng)控制，ECMWF模式對(duì)梅雨帶東段預(yù)報(bào)易發(fā)生北偏。反之，在偏南型中(圖4b)，對(duì)流層中低層暖濕氣流相對(duì)偏弱，副高偏南、偏東，低空西南暖濕氣流偏弱、偏南，我國(guó)東部總體處于西風(fēng)槽控制下，冷空氣偏強(qiáng)，此時(shí)ECMWF模式對(duì)雨帶東段預(yù)報(bào)易發(fā)生偏南。到了中期時(shí)效，雨帶東段偏北的比例開(kāi)始減少，偏南的比例開(kāi)始增多，偏南偏北的系統(tǒng)性偏差逐步消失。而雨帶西段在中期時(shí)效前半段出現(xiàn)了偏南比例上升并超過(guò)偏北比例的現(xiàn)象，而這可能是導(dǎo)致其緯度預(yù)報(bào)平均誤差降低的原因(圖2)。

圖3 2020年梅雨期暴雨雨帶東段和西段 預(yù)報(bào)(a)偏北和(b)偏南的比例Fig.3 Proportion of the (a) northward and (b) southward forecast biases for the western and eastern parts of heavy rain belts during the 2020 Meiyu period

圖4 2020年梅雨期梅雨帶預(yù)報(bào)(a)偏北型和(b)偏南型500 hPa位勢(shì)高度(黑色等值線(xiàn))和850 hPa低空急流的合成分析 (填色為850 hPa風(fēng)速≥12 m·s-1區(qū)域，紅線(xiàn)為588 dagpm等值線(xiàn))Fig.4 Composite analysis of 500 hPa geopotential height (black isoline) and 850 hPa low-level jet for the (a) northward and (b) southward patterns of Meiyu rain belt forecast bias during the 2020 Meiyu period (colored: 850 hPa wind speed ≥12 m·s-1， red line: 588 dagpm isoline)

2.2 雨帶東段和西段位置預(yù)報(bào)的穩(wěn)定性

圖5展示了2020年梅雨期暴雨雨帶東、西段緯度預(yù)報(bào)平均變幅和穩(wěn)定比例。在短期時(shí)效內(nèi)，雨帶東段穩(wěn)定的比例為40%，西段穩(wěn)定的比例略高，為50%～60%，東、西段緯度預(yù)報(bào)平均變幅差距不大，均在0.6°左右。到了中期時(shí)效，雨帶東、西段緯度預(yù)報(bào)穩(wěn)定的比例總體呈下降趨勢(shì)，在中期時(shí)效后半段穩(wěn)定的比例不超過(guò)20%。從中期時(shí)效的緯度預(yù)報(bào)平均變幅來(lái)看，西段沒(méi)有明顯的上升趨勢(shì)，而東段的上升趨勢(shì)相對(duì)明顯，且緯度預(yù)報(bào)的平均變幅(0.8°)總體高于西段(0.6°)。

圖5 2020年梅雨期暴雨雨帶東段和西段 緯度預(yù)報(bào)的平均變幅和穩(wěn)定比例Fig.5 The average latitudinal amplitude of variation and proportion of steability of the western and eastern parts of heavy rain belts during the 2020 Meiyu period

2.3 雨帶東段和西段位置預(yù)報(bào)偏差的連續(xù)性

選取2020年7月5日08時(shí)至8日08時(shí)持續(xù)性暴雨過(guò)程中逐日(08時(shí)至次日08時(shí))暴雨雨帶對(duì)象，用于分析雨帶東、西段位置預(yù)報(bào)偏差的連續(xù)性，考察雨帶前期的位置預(yù)報(bào)偏差信息能否為后一天的預(yù)報(bào)訂正提供參考。圖6展示了該暴雨過(guò)程36 h、60 h和84 h時(shí)效雨帶東、西段緯度預(yù)報(bào)偏差的逐日變化情況。在36 h時(shí)效內(nèi)，雨帶東、西段的預(yù)報(bào)連續(xù)性均表現(xiàn)較好，其中西段表現(xiàn)為5—7日都沒(méi)有明顯的南北預(yù)報(bào)偏差，東段表現(xiàn)為5—7日都有小幅的偏北預(yù)報(bào)偏差；在60 h時(shí)效內(nèi)，雨帶西段的連續(xù)性依舊表現(xiàn)較好，表現(xiàn)為連續(xù)的偏北預(yù)報(bào)偏差，而東段先后出現(xiàn)了漏報(bào)和不連續(xù)的情況；到了84 h時(shí)效，雨帶東、西段均出現(xiàn)一半比例的不連續(xù)。

圖6 2020年7月5日08時(shí)至8日08時(shí)(a,d)36 h、(b,e)60 h、(c,f)84 h 時(shí)效 暴雨雨帶西段(a～c)和東段(d～f)緯度預(yù)報(bào)偏差的逐日變化 (實(shí)線(xiàn)代表預(yù)報(bào)偏差連續(xù)，虛線(xiàn)代表預(yù)報(bào)偏差不連續(xù)，空白代表漏報(bào))Fig.6 Daily variation of latitudinal forecast bias of the (a-c) western, (d-f) eastern parts and the (a, d) 36 h, (b, e) 60 h, (c, f) 84 h forecast of heavy rain belts from 08:00 BT 5 to 08:00 BT 8 July 2020 (Solid lines represent continuous forecast bias, dashed lines represent discontinuous forecast bias, blank represents miss forecast)

2.4 預(yù)報(bào)偏差分析小結(jié)

通過(guò)上述分析可知：短期時(shí)效內(nèi)，雨帶東、西段均存在偏北的系統(tǒng)偏差，雨帶西段的漏報(bào)率低于東段，穩(wěn)定性和連續(xù)性均好于東段；中期時(shí)效內(nèi)，雨帶東段的預(yù)報(bào)誤差增長(zhǎng)較為明顯，且緯度預(yù)報(bào)平均變幅也高于西段。由此可見(jiàn)，雨帶東段的預(yù)報(bào)不確定性總體高于西段，尤其是到了中期預(yù)報(bào)時(shí)效，預(yù)報(bào)不確定性進(jìn)一步加強(qiáng)。在盛夏季節(jié)，降水系統(tǒng)在自西向東發(fā)展的過(guò)程中往往伴隨著地面鋒區(qū)和低空急流的增強(qiáng)，而預(yù)報(bào)員們?cè)陬A(yù)報(bào)實(shí)踐中發(fā)現(xiàn)ECMWF模式常常對(duì)這種增強(qiáng)有過(guò)度預(yù)報(bào)的趨勢(shì)。為了進(jìn)一步分析造成雨帶東段預(yù)報(bào)偏差較大的關(guān)鍵天氣系統(tǒng)，下面將基于ESA方法分析揭示雨帶典型偏差特征的不確定性來(lái)源。

3 集合敏感性分析

3.1 梅雨雨帶東段預(yù)報(bào)明顯偏北的個(gè)例

選取2020年7月18日發(fā)生的一次典型梅雨鋒鋒面雨帶東段預(yù)報(bào)偏北的個(gè)例，采用ESA方法對(duì)其不確定性來(lái)源及相應(yīng)的動(dòng)力學(xué)特征進(jìn)行分析。圖7分別給出了7月18日08時(shí)至19日08時(shí)控制預(yù)報(bào)和觀測(cè)的24 h累計(jì)降水分布(圖7a、7b)及二者在暴雨量級(jí)上的雨帶對(duì)象分布(圖7c)，可以看出預(yù)報(bào)和觀測(cè)雨帶整體均呈現(xiàn)東北—西南走向，二者在西段(114°E以西)的空間位置上呈現(xiàn)出較高的匹配度。但對(duì)于東段而言(114°E以東)，觀測(cè)雨帶呈東西走向影響河南南部到安徽中部一帶，而預(yù)報(bào)雨帶則較觀測(cè)出現(xiàn)明顯偏北，對(duì)河南南部到安徽北部地區(qū)產(chǎn)生影響。

圖7 (a)36 h時(shí)效的控制預(yù)報(bào)和(b)觀測(cè)在2020年7月18日08時(shí)至19日08時(shí)的 累計(jì)降水分布以及(c)預(yù)報(bào)和觀測(cè)的暴雨雨帶對(duì)象 (圖7c中，藍(lán)色為預(yù)報(bào)暴雨對(duì)象，橙色為觀測(cè)暴雨對(duì)象)Fig.7 (a) The 36 h control forecast and (b) observation of the accumulated precipitation distribution from 08:00 BT 18 to 08:00 BT 19 July 2020, and (c) heavy rain belt objects of forecast and observation (inFig.7c, blue shaded: forecast heavy rain object, orange shaded: observed heavy rain object)

圖8給出了7月18日08時(shí)到19日08時(shí)ECMWF集合預(yù)報(bào)逐6 h累積降水的離散度和集合平均分布。通過(guò)繪制相應(yīng)的集合預(yù)報(bào)郵票圖(圖略)發(fā)現(xiàn)，集合離散度大值區(qū)處的強(qiáng)降水雨帶位置存在較大的不確定性。由圖8可以看出降水預(yù)報(bào)不確定性較大(離散度較高)的區(qū)域基本與強(qiáng)降水區(qū)相一致，均首先在西段快速發(fā)展，此后在7月18日20時(shí)后雨帶東段的集合平均降水中心逐漸發(fā)展至紅框所標(biāo)注的預(yù)報(bào)北偏區(qū)域。與之相對(duì)應(yīng)，在19日02—08時(shí)紅框附近的降水離散度呈明顯增加并達(dá)到與西段相近的水平，最大值超過(guò)20 mm，說(shuō)明在該時(shí)段的預(yù)報(bào)不確定性對(duì)雨帶東段預(yù)報(bào)偏北的貢獻(xiàn)最大，可作為此次預(yù)報(bào)“轉(zhuǎn)折點(diǎn)”。

圖8 2020年7月18日08時(shí)至19日08時(shí)ECMWF集合預(yù)報(bào)逐6 h累積降水離散度(填色) 和集合平均(等值線(xiàn)，單位：mm) (a)18日08—14時(shí),(b)18日14—20時(shí),(c)18日20時(shí)至19日02時(shí),(d)19日02—08時(shí) (黑色打點(diǎn)區(qū)域?yàn)榧掀骄邓笥?0 mm的區(qū)域，紅色方框標(biāo)出了ECMWF確定性預(yù)報(bào)雨帶偏北的區(qū)域，下同)Fig.8 The every 6 h ECMWF accumulated precipitation ensemble forecast spread (colored) and ensemble mean (isoline, unit: mm) from (a) 08:00 BT to 14:00 BT 18, (b) 14:00 BT to 20:00 BT 18， (c) 20:00 BT 18 to 02:00 BT 19,(d) 02:00 BT to 08:00 BT 19 July 2020 (black dotted: area of ensemble mean precipitation >10 mm， red box: area of northward rain belt bias of ECMWF determinictic forecast， the same below)

3.2 與500 hPa高度場(chǎng)的集合敏感性

對(duì)預(yù)報(bào)“轉(zhuǎn)折點(diǎn)”進(jìn)行敏感性分析，用以解釋雨帶東段偏北的成因，將7月19日02—08時(shí)紅框區(qū)域內(nèi)的區(qū)域平均6 h降水量作為預(yù)報(bào)量，分別計(jì)算其與7月19日02時(shí)(-6 h)、18日20時(shí)(-12 h)、18日14時(shí)(-18 h)環(huán)境場(chǎng)的ESA。圖9給出了預(yù)報(bào)量與前期500 hPa高度場(chǎng)的ESA，從集合平均場(chǎng)上可以看出在預(yù)報(bào)發(fā)生“轉(zhuǎn)折”前500 hPa高空槽處于發(fā)展加深狀態(tài)，且降水偏差的發(fā)生區(qū)域均位于高空槽前(圖9紅框)，附近存在正負(fù)偶極型分布且通過(guò)顯著性水平檢驗(yàn)的敏感區(qū)，說(shuō)明該高空槽越深(或向下游移速度越慢)，則集合成員在響應(yīng)區(qū)域內(nèi)預(yù)報(bào)的6 h累計(jì)降水量越大，這意味著數(shù)值模式對(duì)預(yù)報(bào)偏北發(fā)生前期高空槽強(qiáng)度的過(guò)報(bào)是導(dǎo)致此次預(yù)報(bào)降水雨帶東段偏北的重要因素之一。注意到越臨近降水發(fā)生偏差的時(shí)刻，槽線(xiàn)位置的負(fù)敏感區(qū)和槽前的正敏感區(qū)均出現(xiàn)緩慢東移且強(qiáng)度增加，呈現(xiàn)出“上游消、下游長(zhǎng)”的發(fā)展效應(yīng)，說(shuō)明ESA在該案例中的應(yīng)用具有清晰的物理意義。另一方面，注意到南側(cè)副高脊線(xiàn)附近存在顯著的正敏感區(qū)域，說(shuō)明副高預(yù)報(bào)的偏強(qiáng)也是導(dǎo)致雨帶出現(xiàn)預(yù)報(bào)偏北的原因之一。

3.3 與850 hPa風(fēng)場(chǎng)的集合敏感性

圖10給出了預(yù)報(bào)量與低層850 hPa水平風(fēng)場(chǎng)的ESA及對(duì)應(yīng)時(shí)次的集合平均場(chǎng)，從集合平均風(fēng)場(chǎng)上可以看出從18日14時(shí)至19日02時(shí)低層存在穩(wěn)定的低空急流向降雨區(qū)輸送暖濕氣流，且急流強(qiáng)度入夜后逐漸加強(qiáng)。在緯向風(fēng)和經(jīng)向風(fēng)的ESA場(chǎng)上，預(yù)報(bào)偏差區(qū)域西南側(cè)均存在明顯的正敏感區(qū)，說(shuō)明模式對(duì)西南低空急流預(yù)報(bào)越強(qiáng)則偏差區(qū)域的降水量越大，與500 hPa高度場(chǎng)的ESA相比，此處敏感區(qū)數(shù)值的絕對(duì)值更大，即850 hPa西南低空急流強(qiáng)度預(yù)報(bào)偏強(qiáng)是導(dǎo)致此次降水雨帶西段偏北的最主要因素。此外，對(duì)比緯向風(fēng)場(chǎng)(圖10a～10c)和經(jīng)向風(fēng)場(chǎng)(圖10d～10f)二者在急流區(qū)正ESA的數(shù)值可以發(fā)現(xiàn)，經(jīng)向風(fēng)場(chǎng)的影響更為突出。另一方面，與500 hPa高度場(chǎng)的ESA不同(圖9)，盡管低空急流在凌晨(19日02時(shí))發(fā)展最為強(qiáng)盛(圖10c、10f)，但此時(shí)的ESA較前兩個(gè)時(shí)次反而有所減弱，考慮到鋒面降水往往具有凌晨加強(qiáng)的日變化特征(Xue et al，2018)，因此在實(shí)際預(yù)報(bào)業(yè)務(wù)中考慮低空急流對(duì)于梅雨鋒暴雨預(yù)報(bào)偏差的判斷時(shí)，需著重對(duì)其前期(6 h前)的急流預(yù)報(bào)進(jìn)行研判。此外，注意到偏差區(qū)域西側(cè)和東南側(cè)存在緯向風(fēng)場(chǎng)的負(fù)敏感區(qū)而西南側(cè)有正敏感區(qū)，說(shuō)明鋒面附近氣旋性風(fēng)切變的預(yù)報(bào)偏高也會(huì)導(dǎo)致降水雨帶預(yù)報(bào)偏北。陳濤等(2020)也對(duì)類(lèi)似的氣旋鋒生梅雨暴雨個(gè)例進(jìn)行了研究，發(fā)現(xiàn)該類(lèi)過(guò)程中ECMWF數(shù)值預(yù)報(bào)常出現(xiàn)低渦急流系統(tǒng)和降水潛熱之間的非自然正反饋，致使天氣系統(tǒng)和雨帶預(yù)報(bào)偏差明顯，與本文結(jié)論一致。

圖9 2020年7月19日02—08時(shí)降水預(yù)報(bào)偏北區(qū)域(紅框)的平均6 h累計(jì)降水與前期500 hPa位勢(shì)高度場(chǎng)的ESA(填色) (a)18日14時(shí)(18 h前)，(b)18日20時(shí)(12 h前)，(c)19日02時(shí)(6 h前) (等值線(xiàn)為對(duì)應(yīng)時(shí)次的集合平均位勢(shì)高度場(chǎng)，單位：gpm;黑色打點(diǎn)區(qū)域表示通過(guò)0.05顯著性水平檢驗(yàn))Fig.9 ESA (colored) of the 6 h average accumulated precipitation with northward forecast bias (red box) from 02:00 BT to 08:00 BT 19 July and the earlier 500 hPa geopotential heights at (a) 14:00 BT 18 (18 h before), (b) 20:00 BT 18 (12 h before), (c) 02:00 BT 19 (6 h before) in July 2020 (Isoline is the ensemble mean geopotential height of corresponding forecast lead time, unit: gpm; black dotted area represents the area having passed the 0.05 significant level test)

圖10 2020年7月19日02—08時(shí)降水預(yù)報(bào)偏北區(qū)域(紅框)的平均6 h累計(jì)降水與 前期850 hPa(a～c)緯向、(d～f)經(jīng)向水平風(fēng)場(chǎng)的ESA(填色) (a，d)18日14時(shí)(18 h前)，(b，e)18日20時(shí)(12 h前)，(c，f)19日02時(shí)(6 h前) (風(fēng)矢量為對(duì)應(yīng)時(shí)次的850 hPa集合平均風(fēng)場(chǎng)，單位：m·s-1；黑色打點(diǎn)區(qū)域表示通過(guò)0.05顯著性水平檢驗(yàn))Fig.10 ESA (colored) of the 6 h average accumulated precipitation with northward forecast bias (red box) from 02:00 BT to 08:00 BT 19 July and the earlier 850 hPa (a-c) zonal and (d-f) meridional horizontal wind fields at (a, d) 14:00 BT 18 (18 h before), (b, e) 20:00 BT 18 (12 h before), (c, f) 02:00 BT 19 (6 h before) in July 2020 (Wind vector is the 850 hPa ensemble mean wind field of corresponding forecast lead time, unit: m·s-1; black dotted area represents the area having passed the 0.05 significant level test)

3.4 集合敏感性分析小結(jié)

通過(guò)上述分析可知，造成鋒面暴雨?yáng)|段預(yù)報(bào)偏北的主要天氣尺度系統(tǒng)包括500 hPa高空低槽、副高以及850 hPa上的低空急流和切變線(xiàn)，數(shù)值模式中高空槽預(yù)報(bào)偏深(曲率增加)、副高偏強(qiáng)、低空急流和鋒面氣旋性切變的偏強(qiáng)均是構(gòu)成此次降水雨帶預(yù)報(bào)偏差的因素，其中低空急流的影響最為明顯。低空急流南風(fēng)分量的預(yù)報(bào)偏強(qiáng)配合副高偏強(qiáng)，導(dǎo)致雨帶北推速度加快，并加強(qiáng)了向降水區(qū)的水汽和能量輸送，和降水潛熱之間存在非自然正反饋，造成虛假?gòu)?qiáng)降水，最終導(dǎo)致了雨帶東段整體偏北。此外，注意到集合平均場(chǎng)上高空槽和低空急流的過(guò)度預(yù)報(bào)均發(fā)生于二者發(fā)展加強(qiáng)的時(shí)段中，因此就這個(gè)典型個(gè)例可認(rèn)為高空槽和低空急流在后半夜發(fā)展加強(qiáng)的時(shí)段是出現(xiàn)梅雨鋒暴雨雨帶東段北偏的背景場(chǎng)特征之一。為了對(duì)該情形加以驗(yàn)證，對(duì)另外4次雨帶東段偏北案例(表1中序號(hào)1、4、7、8)的環(huán)流形勢(shì)和逐6 h 累計(jì)降水進(jìn)行分析(圖略)，得到了與此次個(gè)例一致的結(jié)論。

4 結(jié)論和討論

針對(duì)ECMWF模式在2020年我國(guó)江淮流域超長(zhǎng)梅雨期暴雨雨帶預(yù)報(bào)中的不確定性，選取10個(gè)典型的狹長(zhǎng)暴雨雨帶，基于MODE技術(shù)分析了雨帶東、西段的位置預(yù)報(bào)偏差、穩(wěn)定性以及預(yù)報(bào)偏差的連續(xù)性，并基于ESA方法分析揭示了雨帶典型偏差特征的關(guān)鍵天氣系統(tǒng)。結(jié)果表明：

(1)雨帶東段的預(yù)報(bào)不確定性總體高于西段，尤其是到了中期預(yù)報(bào)時(shí)效，預(yù)報(bào)不確定性進(jìn)一步加強(qiáng)。

(2)短期時(shí)效內(nèi)，雨帶東、西段均存在偏北的系統(tǒng)偏差，雨帶西段的漏報(bào)率低于東段，穩(wěn)定性和連續(xù)性均好于東段；中期時(shí)效內(nèi)，雨帶東段的預(yù)報(bào)誤差增長(zhǎng)較為明顯，且緯度預(yù)報(bào)平均變幅也高于西段。

(3)500 hPa高空槽、副高、850 hPa低空急流和切變線(xiàn)的預(yù)報(bào)偏差及其和降水潛熱之間存在的非自然正反饋共同作用造成雨帶東段預(yù)報(bào)偏北，其中后半夜高空槽和低空急流的發(fā)展加強(qiáng)是出現(xiàn)梅雨鋒暴雨雨帶東段北偏的背景場(chǎng)特征之一。

梅雨期暴雨雨帶的位置預(yù)報(bào)是預(yù)報(bào)業(yè)務(wù)的關(guān)鍵著眼點(diǎn)，掌握模式雨帶預(yù)報(bào)偏差的特性對(duì)于提升預(yù)報(bào)準(zhǔn)確率和防災(zāi)減災(zāi)具有重要意義。在今后的工作中，還需要進(jìn)一步增加更多年份的梅雨期暴雨樣本以提升研究結(jié)論的說(shuō)服力，并針對(duì)更多常用數(shù)值預(yù)報(bào)產(chǎn)品(例如GFS,GRAPES-GFS,GRAPES-3km等)進(jìn)行分析研究。