基于CMA-TRAMS集合預(yù)報(bào)的“5·22”極端降水事件可預(yù)報(bào)性分析*

肖柳斯 張華龍 張旭斌 馮 璐 諶志剛 戴光豐

1.廣州市氣象臺(tái)，廣州，511430 2.廣東省氣象臺(tái)，廣州，510641 3.中國(guó)氣象局廣州熱帶海洋氣象研究所，廣州，510641 4.廣東省區(qū)域數(shù)值天氣預(yù)報(bào)重點(diǎn)實(shí)驗(yàn)室，廣州，510641 5.廣州市氣象局，廣州，511430

1 引言

近30年來，受全球變暖及城市化進(jìn)程的影響，極端強(qiáng)降水有向超大城市群集中并增強(qiáng)的趨勢(shì)（Liang，et al，2017；唐永蘭等，2019）。在強(qiáng)熱島效應(yīng)影響下，珠江三角洲城市群突發(fā)的極端小時(shí)降水事件顯著增加（Wu，et al，2019），破紀(jì)錄的降水事件屢屢發(fā)生。2020年5月21日夜間至22日早晨，廣東省出現(xiàn)了一次暴雨到大暴雨、局地特大暴雨的過程（“5·22”降水事件）。降水中心位于廣州市、東莞市等經(jīng)濟(jì)繁榮、人口密集區(qū)域，1、3 h降水超過當(dāng)?shù)貧v史極值，累積強(qiáng)度大，強(qiáng)降水造成了嚴(yán)重的城市內(nèi)澇、村莊受浸、河堤塌方及山體滑坡等次生災(zāi)害，并導(dǎo)致4人死亡，公交、地鐵、廣深動(dòng)車停運(yùn)，經(jīng)濟(jì)損失非常嚴(yán)重。

“5·22”降水事件為一次季風(fēng)爆發(fā)背景下，低渦和低空急流共同影響的極端強(qiáng)天氣過程。事后，廣東省氣象局組織了專題分析，預(yù)報(bào)員從環(huán)流背景、多源觀測(cè)資料分析和歷史個(gè)例對(duì)比等傳統(tǒng)角度開展了研究，有助于了解該降水事件的天氣學(xué)成因。動(dòng)、熱力和水汽條件均有利于暴雨事件的發(fā)生，短期預(yù)報(bào)能夠關(guān)注到此次暴雨過程，但無(wú)論是主觀預(yù)報(bào)或是確定性模式預(yù)報(bào)對(duì)其極端性估計(jì)均存在明顯不足。由于極端降水是小概率事件，其發(fā)展具有很大的不確定性，用數(shù)值天氣預(yù)報(bào)模式對(duì)其進(jìn)行單一初值的確定性預(yù)報(bào)遠(yuǎn)遠(yuǎn)不夠（高麗等，2019；Wu，et al，2020）。

集合預(yù)報(bào)能夠包容部分偏差，并反映未來天氣的多種可能變化狀況，對(duì)小概率天氣事件具有更好的預(yù)報(bào)效果，是改善數(shù)值模式的不確定性及提高極端事件可預(yù)報(bào)性的科學(xué)有效方法（杜鈞等，2014；高麗等，2019）。目前，世界上多個(gè)氣象中心正在著力發(fā)展集合預(yù)報(bào)業(yè)務(wù)（陳靜等，2020；陳超輝等，2020；Zhang，2018a，2018b； Medina，et al，2019），其中歐洲中期天氣預(yù)報(bào)中心集合預(yù)報(bào)系統(tǒng)（ECMWF-EPS）的業(yè)務(wù)產(chǎn)品在中國(guó)應(yīng)用最為廣泛（Yu，et al，2016；Du，et al，2018；伍志方等，2018；陳濤等，2019）。ECMWF-EPS是全球集合預(yù)報(bào)系統(tǒng)，降水預(yù)報(bào)強(qiáng)度普遍偏弱，對(duì)中尺度對(duì)流系統(tǒng)發(fā)展演變的描述能力不足（Yu，et al，2016；王毅等，2020），因此區(qū)域甚至對(duì)流尺度集合預(yù)報(bào)系統(tǒng)的研發(fā)備受重視。中國(guó)氣象局廣州熱帶海洋氣象研究所基于CMA-MESO（原GRAPES_MESO）模式開發(fā)了中國(guó)氣象局南海臺(tái)風(fēng)數(shù)值預(yù)報(bào)系統(tǒng)（Tropical Regional Atmosphere Model for the South China Sea，CMA-TRAMS；Chen，et al，2014）（原南海臺(tái)風(fēng)數(shù)值預(yù)報(bào)系統(tǒng)，TRAMS）。Zhang（2018a）基于CMA-TRAMS建立了熱帶中尺度集合預(yù)報(bào)系統(tǒng)（CMA-TRAMS（EPS）），并通過3 a的批量試驗(yàn)證明CMA-TRAMS（EPS）在臺(tái)風(fēng)強(qiáng)度、強(qiáng)降水及大風(fēng)預(yù)報(bào)方面比ECMWFEPS更具優(yōu)勢(shì)。2020年針對(duì)CMA-TRAMS進(jìn)行了版本升級(jí)，CMA-TRAMS（EPS）已投入業(yè)務(wù)試運(yùn)行（陳子通等，2020）。

集合預(yù)報(bào)系統(tǒng)中各成員具有相似的天氣、氣候背景，經(jīng)過積分計(jì)算后，預(yù)報(bào)結(jié)果存在差異，從誤差演變的角度能夠衡量預(yù)報(bào)變量與集合成員模式初值及物理過程的耦合關(guān)系，近年來被廣泛應(yīng)用于天氣系統(tǒng)動(dòng)力特征和實(shí)際可預(yù)報(bào)性的研究（Torn，et al，2008；Yu，et al，2016；王毅等，2017；閔錦忠等，2020）。常用的研究方法包括組間差異分析和集合敏感性分析（ESA）（Schumacher，2011；代刊等，2018）。相較于歷史個(gè)例的對(duì)比，集合敏感性分析方法更適用于挖掘影響模式降水的關(guān)鍵因子，有助于加深對(duì)典型天氣發(fā)生演變的機(jī)理研究及模式可預(yù)報(bào)性能研究。作為研究預(yù)報(bào)誤差和可預(yù)報(bào)性的新工具，集合敏感性分析方法在溫帶和熱帶氣旋（Torn，et al，2009；Brown，et al，2015；代 刊 等，2018；陳偉斌等，2019）、阻塞高壓（Matsueda，et al，2011；Quandt，et al，2019）、東風(fēng)波（Torn，2010）、強(qiáng)對(duì)流（Hill，et al，2016； Limpert，et al，2018）、暴雨（Yu，et al，2016；王毅等，2017；Du，et al，2018；陳濤等，2019；Wu，et al，2020）等多種類型天氣系統(tǒng)和高影響天氣的分析和機(jī)理研究中得到應(yīng)用。

CMA-TRAMS（EPS）上線運(yùn)行以來，目前仍少見針對(duì)該集合預(yù)報(bào)系統(tǒng)在極端暴雨天氣中的業(yè)務(wù)性能檢驗(yàn)和實(shí)例分析方面的研究。本研究針對(duì)CMA-TRAMS（EPS）在“5·22”降水事件中反映的集合成員降水與天氣尺度系統(tǒng)、中尺度系統(tǒng)的因果聯(lián)系展開討論，探討CMA-TRAMS（EPS）集合預(yù)報(bào)產(chǎn)品對(duì)該類極端暴雨個(gè)例的預(yù)報(bào)能力，以期為提高該類極端天氣的預(yù)報(bào)水平提供參考，增強(qiáng)對(duì)過程機(jī)理和CMA-TRAMS（EPS）性能的理解，并為后續(xù)深入研究集合預(yù)報(bào)系統(tǒng)對(duì)于華南極端天氣預(yù)報(bào)的適用和改進(jìn)方向提供基礎(chǔ)。

2 資料和方法

2.1 資 料

為了了解CMA-TRAMS（EPS）的基本性能，用ECMWF-EPS做參考，分別對(duì)兩套集合預(yù)報(bào)產(chǎn)品進(jìn)行主、客觀檢驗(yàn)。起報(bào)時(shí)間為2020年5月21日08時(shí)和20時(shí)（北京時(shí)，下同）。ECMWF-EPS共50個(gè)擾動(dòng)集合成員和1個(gè)對(duì)照預(yù)報(bào)成員，預(yù)報(bào)產(chǎn)品的水平網(wǎng)格分辨率為0.5°×0.5°，預(yù)報(bào)間隔為6 h。CMA-TRAMS（EPS）以CMA-TRAMS（陳子通等，2020）為基礎(chǔ)，通過對(duì)ECMWF-EPS分析場(chǎng)與預(yù)報(bào)場(chǎng)進(jìn)行降尺度處理獲得各個(gè)集合成員的擾動(dòng)初值與側(cè)邊界條件。CMA-TRAMS垂直方向上采用高度地形追隨坐標(biāo)，并采用Charney-Philips躍層。積云對(duì)流參數(shù)化方案采用SAS（Han，et al，2011），微物理過程采用WSM6方案（Hong，et al，2006），邊界層物理方案使用MRF（Hong，et al，1996），長(zhǎng)波輻射過程采用RRTM方案，短波輻射過程采用Dudhia方案。CMA-TRAMS（EPS）共30個(gè)擾動(dòng)集合成員和1個(gè)對(duì)照預(yù)報(bào)成員，范圍覆蓋（0.8°—54.8°N，70°—160°E），水平網(wǎng)格數(shù)為1001（經(jīng)向）、601（緯向），垂直方向共65層，輸出等壓面標(biāo)準(zhǔn)為17層（1000—10 hPa），水平網(wǎng)格分辨率為9 km，預(yù)報(bào)間隔為1 h。本次過程分析重點(diǎn)關(guān)注預(yù)報(bào)時(shí)效36 h以內(nèi)。

實(shí)況資料包括廣東省內(nèi)共2895個(gè)自動(dòng)雨量站（含國(guó)家級(jí)氣象站及區(qū)域氣象站）逐5 min觀測(cè)數(shù)據(jù)、中國(guó)氣象局陸面融合實(shí)況分析產(chǎn)品（CLDAS）及ECMWF的再分析資料（ERA5），具體參數(shù)如表1所示。由于各數(shù)據(jù)時(shí)、空分辨率不一致，首先對(duì)數(shù)據(jù)進(jìn)行預(yù)處理。

表1 數(shù)據(jù)集說明Table 1 Description of datasets

2.2 集合敏感性分析（ESA）方法

集合敏感性分析（ESA）方法假設(shè)集合成員的預(yù)報(bào)變量與某個(gè)因子預(yù)報(bào)場(chǎng)存在線性關(guān)系，分析輸入變量或模型的變化可了解所關(guān)心變量的受影響程度。CMA-TRAMS（EPS）有30個(gè)擾動(dòng)成員，預(yù)報(bào)敏感度（S）可定義為（Torn，et al，2008）

式中，J為預(yù)報(bào)變量值，Xi為i成員的狀態(tài)變量，cov為協(xié)方差算子，var為方差算子。經(jīng)過數(shù)學(xué)轉(zhuǎn)換可得S與J、Xi的相關(guān)系數(shù)cor成正比（王毅等，2020）。即

因此，相關(guān)系數(shù)的絕對(duì)值越大，初始狀態(tài)變量的擾動(dòng)引起預(yù)報(bào)結(jié)果的變動(dòng)越大，即預(yù)報(bào)量對(duì)狀態(tài)變量較為敏感。

2.3 集合組間差異分析方法

集合組間差異分析法指設(shè)定分類標(biāo)準(zhǔn)將各成員分組，再合成分析各組成員的環(huán)境參量并對(duì)比差異。以降水為例，降水預(yù)報(bào)的落區(qū)和強(qiáng)度均可作為分組依據(jù)。如果把落區(qū)準(zhǔn)確而強(qiáng)度偏弱或偏強(qiáng)的成員歸類，也可作為一個(gè)成員組。合成場(chǎng)反映各組成員的共有特征，通過對(duì)比分組合成場(chǎng)的差異，能夠挖掘兩組成員共同的差別（Lamberson，et al，2016；Wu，et al，2020）。例如，伍志方等（2018）通過對(duì)比“5·7”過程中ECMWF-EPS分組成員各層要素的分布差異，指出地面系統(tǒng)的異常是中尺度系統(tǒng)觸發(fā)的主因。

文中將預(yù)報(bào)落區(qū)和強(qiáng)度均相對(duì)較好的成員歸類為“好成員（GOOD）”。除GOOD組以外，其他成員的結(jié)果會(huì)有多種表現(xiàn)，可能存在顯著的預(yù)報(bào)位置或強(qiáng)度偏差，甚至部分成員強(qiáng)度和位置的偏差均較大。由于本研究的目標(biāo)在于研究GOOD成員跟大部分成員的差異，因此，除“GOOD”成員之外，其余成員均歸入“其他成員（OTHERS）”組。

3 “5·22”降水事件概況

3.1 降水分布情況

“5·22”降水事件中珠江三角洲和蓮花山南部的惠州市至汕頭市一帶出現(xiàn)了暴雨以上降水，最強(qiáng)降水中心位于珠江三角洲中部，東莞市高埗鎮(zhèn)中心村測(cè)得全省最大累計(jì)雨量399.50 mm（圖1a）。該過程短時(shí)雨量大，局地強(qiáng)降水持續(xù)時(shí)間長(zhǎng)，廣州市東部和東莞市中北部地區(qū)部分站點(diǎn)（圖1b）小時(shí)雨強(qiáng)非常極端（120 mm/h，曾智琳等，2020），廣州市黃埔大橋站測(cè)得最大整點(diǎn)雨強(qiáng)167.8 mm/h，滑動(dòng)小時(shí)雨量達(dá)到201.8 mm（22日02時(shí)15分到03時(shí)15分）；東莞市東城街道的最大3 h雨量達(dá)到351 mm（22日00時(shí)30分到03時(shí)30分），打破東莞市本地歷史記錄。超過5 mm/（5 min）的強(qiáng)降水持續(xù)了4—5 h（圖1c），以黃埔大橋站（圖1d）為例：強(qiáng)降水從01時(shí)40分開始，02時(shí)和03時(shí)前后分別達(dá)到次峰值和峰值，03時(shí)前后，連續(xù)3個(gè)時(shí)次的最大5 min降水量超過30 mm，0.5 h降水量已超過100 mm，降水效率極高。

圖1 2020 年5 月21 日14 時(shí)至22 日14 時(shí)（a）廣東省累計(jì)降水量（黑框?yàn)榻邓行膮^(qū)，與b、c范圍一致）、（b）最大小時(shí)降水量、（c）超過5 mm/（5 min）降水的累計(jì)時(shí)間和（d）21 日20 時(shí)至22 日08 時(shí)廣州市黃埔區(qū)黃埔大橋逐5 min降水量Fig.1 （a）Accumulative precipitation（the black box is the central area of precipitation，which is the same as the area shown in b and c）in Guangdong，（b）maximum hourly precipitation，（c）amount of time with 5 min precipitation exceeding 5 mm from 14：00 BT 21 to 14：00 BT 22 ，and（d）time series of 5 min precipitation at Huangpu Bridge Station from station from 20：00 BT 21 to 08：00 BT

3.2 天氣形勢(shì)分析

以ERA5資料為基礎(chǔ)，對(duì)2020年5月21日20時(shí)至22日02時(shí)的環(huán)流形勢(shì)進(jìn)行分析。200 hPa，華南地區(qū)位于南亞高壓東側(cè)及高空西風(fēng)急流右側(cè)的扇形輻散區(qū)（圖2a），夜間到凌晨良好的高空輻散條件維持并略有加強(qiáng)（圖2e）。500 hPa華南地區(qū)位于東亞大槽后部，受槽后脊前西北氣流控制影響，與同一時(shí)刻多年（2001—2018年）的平均態(tài)相比，中國(guó)東部地區(qū)高度場(chǎng)顯著偏低，而大陸西部地區(qū)偏高（圖2b）。02時(shí)西部副熱帶高壓的強(qiáng)度維持，114°E以東副熱帶高壓的強(qiáng)度有所減弱，“西高東低”的特征更加明顯，西偏北氣流引導(dǎo)低層系統(tǒng)東移南壓，珠江三角洲地區(qū)東側(cè)可見短波槽活動(dòng)，中層的動(dòng)力條件進(jìn)一步轉(zhuǎn)好（圖2f）。925 hPa，華南地區(qū)的暖濕區(qū)與低壓槽相伴發(fā)展東移（圖2c），02時(shí)廣東省東部地區(qū)渦度增大并可見閉合等高線，低渦南側(cè)的偏西和西南兩支超低空急流顯著增強(qiáng)，低層輻合和水汽輸送條件進(jìn)一步改善（圖2g）。地面，中尺度輻合線逐漸東移南壓至珠江三角洲東部地區(qū)，有利于局地強(qiáng)降水系統(tǒng)的觸發(fā)。整層可降水量高值區(qū)域范圍顯著擴(kuò)大，強(qiáng)度增強(qiáng)，中心值超過85 kg/m2（圖2d、h），遠(yuǎn)大于當(dāng)日多年（2001—2018年）平均值（58 kg/m2），充沛的水汽條件異常有利于極端性暴雨過程的出現(xiàn)?？偠灾?，夜間系統(tǒng)配置向更有利于強(qiáng)降水的趨勢(shì)轉(zhuǎn)變。

圖2 2020年5月 （a、b、c、d） 21日20時(shí)和 （e、f、g、h） 22日02時(shí) （a、e） 200 hPa風(fēng)場(chǎng) （色階），（b、f） 500 hPa高度場(chǎng) （等值線，單位：gpm） 及其與多年平均差值場(chǎng) （色階），（c、g） 925 hPa高度場(chǎng)（等值線，單位：gpm）、假相當(dāng)位溫 （色階） 和急流 （大于12 m/s），（d、h） 海平面氣壓 （等值線，單位：hPa）、整層可降水量 （色階） 和地面風(fēng)場(chǎng)Fig.2 （a，e） Wind field at 200 hPa （shaded），（b，f） geopotential height （contour，unit：gpm） and anomalies from multi-year average （shaded） at 500 hPa，（c，g） geopotential height （contour，unit：gpm），pseudo-equivalent potential temperature （shaded）and jet stream （over 12 m/s） at 925 hPa，and （d，h） sea-level pressure （contour，unit：hPa），total precipitable water （shaded） and wind field at 20：00 BT 21 （a，b，c，d） and 02：00 BT 22 （e，f，g，h） May 2020

續(xù)圖2Fig.2 Continued

4 集合預(yù)報(bào)結(jié)果檢驗(yàn)

以當(dāng)今預(yù)報(bào)水平較高的業(yè)務(wù)系統(tǒng)ECMWFEPS為參考，與CMA-TRAMS（EPS）進(jìn)行主、客觀對(duì)比檢驗(yàn)，可了解CMA-TRAMS（EPS）對(duì)本次極端降水過程的預(yù)報(bào)能力。由于過程的主要降水時(shí)段為5月22日凌晨，檢驗(yàn)的目標(biāo)時(shí)段選定為21日20時(shí)至22日08時(shí)的12 h降水。

分別對(duì)比2個(gè)模式在2個(gè)起報(bào)時(shí)（21日08時(shí)和20時(shí)）的具體預(yù)報(bào)。從08時(shí)起報(bào)的降水分布郵票圖（圖3）可見，ECMWF-EPS中主要降水落區(qū)與CMA-TRAMS（EPS）相似，但降水中心強(qiáng)度差異較為明顯； 20時(shí)起報(bào)的ECMWF-EPS降水量調(diào)整增大，但仍只有少數(shù)成員預(yù)報(bào)出局地100 mm量級(jí)的降水（圖4）。08時(shí)起報(bào)的CMA-TRAMS（EPS）降水強(qiáng)度普遍較大，強(qiáng)降水范圍也偏大，強(qiáng)降水中心主要位于粵東地區(qū)；20時(shí)起報(bào)的CMA-TRAMS（EPS）雨量較08時(shí)有所減小，成員之間分歧增大，但部分成員的量級(jí)、落區(qū)位置和形態(tài)預(yù)報(bào)均與實(shí)況更相符。

圖3 2020年5月21日08時(shí)起報(bào)的2020年5月21日20時(shí)至22日08時(shí)12 h累計(jì)降水量 （ECMWF-EPS成員以“E”開頭，CMA-TRAMS（EPS）成員以“C”開頭，單位：mm）Fig.3 Forecasts of 12 h accumulated precipitation from 20：00 BT 21 to 08：00 BT 22 May 2020 initialized at 08：00 BT 21 May 2020 （ECMWF-EPS members begin with "E" and CMA-TRAMS（EPS） members begin with "C"，unit：mm）

為客觀評(píng)價(jià)兩個(gè)模式各成員之間及不同起報(bào)時(shí)刻的結(jié)果差異，利用雙線性插值法將模式預(yù)報(bào)結(jié)果插值到2895個(gè)雨量站上，計(jì)算21日20時(shí)至22日08時(shí)廣東省全部自動(dòng)雨量站的TS（Threat score）評(píng)分和空?qǐng)?bào)率（FAR）以進(jìn)行定量檢驗(yàn)。結(jié)果顯示，ECMWF-EPS的預(yù)報(bào)結(jié)果調(diào)整明顯，08時(shí)起報(bào)的100 mm降水TS評(píng)分均為0，僅1（8）個(gè)成員50 mm降水TS大于0.4（0.2）。部分成員TS較低的同時(shí)，F(xiàn)AR較高（圖5a1、a2），表明預(yù)報(bào)落區(qū)與實(shí)況位置有偏差。20時(shí)起報(bào)的結(jié)果中，TS明顯提升，8個(gè)成員100 mm降水TS大于0，17（6）個(gè)成員50 mm降水TS超過0.2（0.4）。但FAR也普遍上升，部分成員（如25號(hào)）FAR比TS更大（圖5b1、b2），配合落區(qū)分布（圖4）可見落區(qū)的范圍大小及位置與實(shí)況存在明顯偏差。

圖4 同圖3，但起報(bào)時(shí)刻為2020年5月21日20時(shí) （紅框?yàn)椤癎OOD”組成員）Fig.4 Same as Fig.3 except that the initial time is 20：00 BT 21 May 2020 （the red boxes indicate GOOD group members）

CMA-TRAMS（EPS）不同時(shí)刻起報(bào)的結(jié)果也有明顯不同，08時(shí)起報(bào)的降水范圍和強(qiáng)度均較大，50 mm降水TS超過0.2的成員（圖5a3）數(shù)量遠(yuǎn)多于20時(shí)起報(bào)的結(jié)果（圖5b3），19個(gè)成員100 mm降水TS大于0，但差異較大。大多數(shù)成員50 mm降水FAR超過0.2（圖5a4），意味著08時(shí)的高TS是以高FAR為代價(jià)的。而20時(shí)起報(bào)的TS評(píng)分（圖5b3）顯示各成員的差異較大，部分成員（如7、23號(hào)）在50 mm和100 mm降水中均具有較高TS，且FAR小于0.2（圖5b4），表明這些成員對(duì)降水落區(qū)的空間位置把握能力更強(qiáng)。

圖5 2020年5月21日08時(shí) （a1—a4） 和20時(shí) （b1—b4） 起報(bào)的廣東省21日20時(shí)至22日08時(shí)12 h累計(jì)降水量的ECMWF-EPS （a1、b1、a2、b2） 和 CMA-TRAMS（EPS） （a3、b3、a4、b4） 的TS評(píng)分和空?qǐng)?bào)率Fig.5 TS and FAR for ECMWF-EPS （a1，b1，a2，b2） and CMA-TRAMS（EPS） （a3，b3，a4，b4） forecasts of 12 h accumulated precipitation from 20：00 BT 21 to 08：00 BT 22 May 2020 in Guangdong province（the forecasts are initialized at 08：00 （a1—a4） and 20：00 （b1—b4） BT 21 May 2020，respectively）

通過對(duì)比集合預(yù)報(bào)的平均場(chǎng)和離散度可進(jìn)一步了解成員降水的預(yù)報(bào)性能。20時(shí)起報(bào)的12 h集合平均降雨量和離散度分布（圖6）表明，ECMWFEPS平均降水量均在100 mm以下，存在3個(gè)超過50 mm的降水中心，其中珠江三角洲北部和東南部局地性較強(qiáng)，而粵東降水中心沿著海岸線分布且范圍較大（圖6a）。離散度大值中心位于珠江三角洲北部，次中心位于粵東沿海（圖6b），分別與兩個(gè)降水中心相對(duì)應(yīng)。

CMA-TRAMS（EPS）預(yù)報(bào)降水中心位置偏東，雨量超過50 mm的范圍更廣，集中在粵東海岸線附近（圖6c）。CMA-TRAMS（EPS）的集合離散度較ECMWF-EPS顯著偏大，最大值位于珠江三角洲東南部（圖6d紅色框），其中心值超過50 mm。次中心（圖6d白色框）位置與ECMWF-EPS相似，但是量值明顯偏大。

降水量與集合離散度存在一定相關(guān)，降水量較大時(shí)，離散度絕對(duì)值往往較大（蘇翔等，2012；Su，et al，2014），兩個(gè)模式預(yù)報(bào)的粵東沿海（圖6b、d白色框）及ECMWF-EPS預(yù)報(bào)的珠江三角洲北部地區(qū)（圖6b紅色框）離散度高值中心除了與模式預(yù)報(bào)不確定程度有關(guān)外，也與集合平均值較大有關(guān)。珠江三角洲東部，尤其是東南部地區(qū)（圖6d紅色框）CMA-TRAMS（EPS）預(yù)報(bào)的平均降水量25—50 mm，量值較小，但是離散度較周邊均顯著偏大，表明CMA-TRAMS（EPS）對(duì)珠江三角洲東部地區(qū)的降水預(yù)報(bào)具有高度的不確定性，部分成員預(yù)報(bào)的降水量具有一定極端性。實(shí)際出現(xiàn)的極端降水恰是在此區(qū)域附近，反映出CMA-TRAMS（EPS）部分成員對(duì)此次極端性降水具有一定的預(yù)報(bào)能力。該區(qū)域內(nèi)ECMWF-EPS預(yù)報(bào)的平均場(chǎng)量級(jí)與CMA-TRAMS（EPS）類似，但是離散度較小，預(yù)報(bào)不確定性較小，ECMWF-EPS預(yù)測(cè)該區(qū)域出現(xiàn)極端降水的概率較小。

圖6 2020年5月21日20時(shí)起報(bào)的2020年5月21日20時(shí)至22日08時(shí)ECMWF-EPS （a、b） 和CMA-TRAMS（EPS）（c、d） 的 （a、c） 集合平均降水量和 （b、d） 離散度 （單位：mm，紅色虛線框和白色實(shí)線框代表高離散度區(qū)域）Fig.6 （a，c） Ensemble mean forecasts and （b，d） spreads of accumulative precipitation forecasts for the period from 20：00 BT 21 to 08：00 BT 22 May by （a，b） ECMWF-EPS and （c，b） CMA-TRAMS（EPS） initialized at 20：00 BT 21 May 2020（unit：mm，the red dashed box and the white solid box indicate two areas of high spread）

上述分析表明，ECMWF-EPS預(yù)報(bào)的雨量普遍較小，對(duì)極端降水的預(yù)報(bào)明顯不足，而CMA-TRAMS（EPS）預(yù)報(bào)的強(qiáng)度更接近實(shí)況。不同時(shí)刻起報(bào)的集合預(yù)報(bào)結(jié)果均有一定調(diào)整，CMA-TRAMS（EPS）在20時(shí)起報(bào)的結(jié)果中，部分成員無(wú)論對(duì)降水強(qiáng)度或落區(qū)均有較好的預(yù)報(bào)效果，與其他成員存在較大的差異，不同成員組之間具有較高的對(duì)比意義，以后將選用最臨近時(shí)刻（20時(shí)起報(bào)）的結(jié)果進(jìn)行分析。通過對(duì)比主觀及客觀定量檢驗(yàn)結(jié)果，CMA-TRAMS（EPS）的2、7、11、17、23和25號(hào)成員被選入“GOOD”組，其他成員歸入“OTHERS”組。在粵東沿海地區(qū)，兩組成員預(yù)報(bào)均存在一個(gè)集合平均雨量大值區(qū)，范圍及強(qiáng)度相當(dāng)。而“GOOD”組預(yù)報(bào)在珠江三角洲東部地區(qū)還有一個(gè)降水中心，落區(qū)的形態(tài)與實(shí)況相似，但是極值中心較實(shí)況弱。其中23號(hào)是最佳成員，其預(yù)報(bào)落區(qū)和強(qiáng)度與實(shí)況都具有較高的相似度。

5 結(jié)果

5.1 多因子的集合敏感性分析

文中將集合敏感性分析方法應(yīng)用于CMATRAMS（EPS）中，針對(duì)動(dòng)力因子（海平面氣壓、500 hPa位勢(shì)高度、200 hPa位勢(shì)高度、850 hPa緯向風(fēng)、925 hPa緯向風(fēng)、850 hPa經(jīng)向風(fēng)、925 hPa經(jīng)向風(fēng)）、水汽因子（925 hPa比濕、地面比濕）和熱力因子（925 hPa假相當(dāng)位溫、925與500 hPa假相當(dāng)位溫差、地面溫度）計(jì)算集合敏感度，以對(duì)比極端降水事件對(duì)不同因子的敏感程度。

結(jié)果（圖7）表明，敏感度最高的是850 hPa經(jīng)向風(fēng)（0.84），其次為500 hPa位勢(shì)高度場(chǎng)（?0.83），敏感度絕對(duì)值超過0.7的主要為動(dòng)力因子，如850 hPa緯 向 風(fēng)（?0.78）、925 hPa經(jīng) 向 風(fēng)（0.77，?0.73）、925 hPa緯 向 風(fēng)（?0.73）和 海 平 面 氣 壓（0.73）。由此可見，850 hPa和925 hPa風(fēng)場(chǎng)的配置與模式預(yù)報(bào)降水關(guān)系非常密切，其中經(jīng)向風(fēng)尤為重要。模式降水與這兩層經(jīng)向風(fēng)的正相關(guān)顯著高于負(fù)相關(guān)，南風(fēng)正向擾動(dòng)更有利于降水增加，3.2節(jié)中所提及的（超）低空急流增強(qiáng)正是降水增加的關(guān)鍵配置。此外，除了500 hPa高度場(chǎng)外，其他因子正、負(fù)敏感區(qū)多數(shù)成對(duì)出現(xiàn)，相關(guān)均通過了99%顯著性t檢驗(yàn)。850 hPa和925 hPa的風(fēng)場(chǎng)尤為顯著，其中經(jīng)向風(fēng)正、負(fù)敏感區(qū)分別位于東南和西北方向，而緯向風(fēng)大致為南北向分布，成對(duì)風(fēng)場(chǎng)敏感度的配置體現(xiàn)了CMA-TRAMS（EPS）降水與低層渦度場(chǎng)的強(qiáng)度具有密切的正相關(guān)關(guān)系。

圖7 CMA-TRAMS（EPS）集合成員預(yù)報(bào)敏感度極值（黑色：正，灰色：負(fù)；黑色虛線對(duì)應(yīng)相關(guān)系數(shù)通過99%的顯著性t檢驗(yàn)，GH和PPT分別為位勢(shì)高度和假相當(dāng)位溫縮寫）Fig.7 Extremum of forecast sensitivity of CMA-TRAMS（EPS） members （black：positive sensitivity，grey：minus sensitivity；the black dotted line corresponds to the correlation coefficient passing the 99% confidence level；GH and PPT are the abbreviation names of geopotential height and potential pseudo-equivalent temperature，respectively）

低層溫度、濕度和不穩(wěn)定能量條件與降水也存在一定的相關(guān)，敏感度量級(jí)相當(dāng)，絕對(duì)值為0.5—0.6。廣東省中南部陸地的地面溫度敏感度超過0.6，該地區(qū)溫度越高，越有利于研究區(qū)（圖1a黑框所示）出現(xiàn)強(qiáng)降水，但是周邊的溫度負(fù)相關(guān)未能通過99%顯著性t檢驗(yàn)。

綜上所述，集合預(yù)報(bào)對(duì)降水強(qiáng)度的預(yù)報(bào)與不同高度層的動(dòng)、熱力因子均具有一定的聯(lián)系，但相關(guān)存在差別，相關(guān)最高為表征中、低層的動(dòng)力因子（如850 hPa和 925 hPa風(fēng)場(chǎng)、500 hPa高度場(chǎng)），表明此次降水過程的天氣系統(tǒng)（即低渦和低空急流）的動(dòng)力特征與降水強(qiáng)度關(guān)系最密切，需開展進(jìn)一步的分析。

5.2 天氣尺度系統(tǒng)的可預(yù)報(bào)性分析

5.2.1 低渦

本次降水過程的低層切變系統(tǒng)由中國(guó)西南地區(qū)的低渦東移演變而來。南方汛期，西南渦與南風(fēng)急流是華南地區(qū)致洪暴雨重要組成成員（黃明策等，2010），并且強(qiáng)降水主要分布在環(huán)流中心或系統(tǒng)移向的右前側(cè)。由于局地累計(jì)雨量嚴(yán)重依賴于降水系統(tǒng)的移動(dòng)方式，移速慢的氣旋有助于降水維持和積累，因此低渦的移動(dòng)路徑直接影響落區(qū)的位置和降水強(qiáng)度（Li，et al，2014）。

CMA-TRAMS（EPS）各成員低渦的預(yù)報(bào)路徑均大致呈緯向，位于粵中的偏東地區(qū)至粵東偏北地區(qū)，分布跨度約為2.5緯度、3.5經(jīng)度，離散度較大，部分成員移速較快，12 h內(nèi)移動(dòng)路徑更長(zhǎng)（圖8a）；隨著預(yù)報(bào)時(shí)效延長(zhǎng)，離散度更趨明顯。在所有成員中，GOOD組成員的低渦東移為主，移速呈先慢后快的特征（圖8b），22日02時(shí)該組成員的低渦中心位置較其他成員明顯偏西（圖8c），演變趨勢(shì)與實(shí)況較為一致。其中7和23號(hào)成員的移動(dòng)路徑與實(shí)況最接近，兩者TS評(píng)分也最高，F(xiàn)AR較?。欢?1和25號(hào)較實(shí)況略偏南，F(xiàn)AR較高，表明降水落區(qū)仍具有一定偏差（圖5）。

以低渦東移速度作為預(yù)報(bào)量，對(duì)各因子做敏感度計(jì)算，可見低渦東移速度對(duì)經(jīng)向風(fēng)速敏感度較高。集合敏感度正值較小，低渦東移速度與北風(fēng)的相關(guān)較弱。負(fù)值中心粵中偏東地區(qū)，極小值低于?0.8，表明該區(qū)域南風(fēng)越大，低渦東移速度越慢（圖8d），珠江三角洲東部地區(qū)的南風(fēng)匯入有利于低渦在其北部停滯。

圖8 850 hPa層CMA-TRAMS（EPS）2020年5月21日21時(shí)至22日08時(shí)（a）所有成員和（b）GOOD組成員預(yù)報(bào)的低渦移動(dòng)路徑 （黑色路徑為再分析場(chǎng)），（c）22日02時(shí)低渦中心位置，（d）低渦中心經(jīng)度與南風(fēng)的敏感度（網(wǎng)格表示通過99%顯著性t檢驗(yàn)區(qū)域）Fig.8 Low vortex tracks from 21：00 BT 21 to 08：00 BT 22 May predicted by （a） all members and （b） GOOD group members of CMA-TRAMS（EPS） , （c） location of low vortex center at 02：00 BT 22 May 2020, （d） sensitivity between longitude of low vortex and southerly wind at 850 hPa （areas that pass the 99% confidence levelt-test are represented by grids）

對(duì)比CMA-TRAMS（EPS）和ERA5的低空風(fēng)分布有助于了解模式對(duì)低層形勢(shì)的預(yù)報(bào)偏差，進(jìn)一步得出模式對(duì)低渦預(yù)報(bào)偏差的成因。CMATRAMS（EPS）的偏東風(fēng)場(chǎng)較再分析場(chǎng)更強(qiáng)，珠江三角洲北部地區(qū)尤為明顯（圖9）。其中GOOD組成員（圖9a、d）對(duì)珠江三角洲東側(cè)的偏南風(fēng)和西側(cè)的偏北風(fēng)均預(yù)報(bào)偏強(qiáng)、南側(cè)的西南風(fēng)預(yù)報(bào)略偏強(qiáng)，各向風(fēng)場(chǎng)的配置致使GOOD組成員在珠江三角洲地區(qū)形成了偏強(qiáng)的氣旋性環(huán)流，低渦預(yù)報(bào)較實(shí)況偏強(qiáng)。OTHERS組成員在廣東省境內(nèi)的850 hPa偏北風(fēng)較實(shí)況顯著偏強(qiáng)，在珠江三角洲北部形成北風(fēng)差值中心。沿海差值西南風(fēng)與陸上差值東北風(fēng)交匯，在珠江口東側(cè)沿海增強(qiáng)了氣旋性環(huán)流（圖9b），造成與實(shí)況偏離的降雨中心。由于925 hPa上OTHERS組成員的偏北風(fēng)偏差更顯著，經(jīng)向風(fēng)的輻合區(qū)南推至海上，偏差風(fēng)的氣旋性特征不明顯（圖9e）。兩組成員的差異表明，廣東省中部GOOD組成員氣旋性環(huán)流更強(qiáng)，低渦強(qiáng)度比OTHERS組成員大得多（圖9c、f）。

上述特征表明，GOOD組成功預(yù)報(bào)了低渦緩慢移動(dòng)的特征，同時(shí)預(yù)報(bào)氣旋性渦度偏強(qiáng)，有利于在珠江三角洲地區(qū)預(yù)報(bào)較強(qiáng)的降水，部分成員24 h降雨量可達(dá)200 mm。即便如此，模式降水的強(qiáng)度仍遠(yuǎn)不及實(shí)際出現(xiàn)的極端降水強(qiáng)度。因此，除低渦以外，還有低空急流等其他系統(tǒng)對(duì)極端降水具有影響。

5.2.2 低空急流

低空急流輸送暖濕空氣，增大風(fēng)垂直切變，從而增強(qiáng)大氣環(huán)境的不穩(wěn)定性，對(duì)鋒面降水或暖區(qū)對(duì)流觸發(fā)和維持均有重要作用（劉鴻波等，2014；Du，et al，2014；Huang，et al，2017）。由上文分析也可見偏南風(fēng)急流有利于低渦的停滯，兩個(gè)天氣系統(tǒng)的相互作用對(duì)局地暴雨具有一定貢獻(xiàn)。

Du等（2018）指出，鋒面和暖區(qū)降水強(qiáng)度分別受到低空急流和超低空急流主導(dǎo)，并且兩者的敏感中心分別位于陸面和海面。就本次過程中850 hPa的低空急流而言，兩組成員差異最明顯的區(qū)域位于珠江三角洲東側(cè)陸面（圖9c），與Du等（2018）分析的鋒面降水的敏感中心一致。與實(shí)況對(duì)比，GOOD組成員的西南急流強(qiáng)度偏強(qiáng)（圖9a），自海面向陸地，經(jīng)向風(fēng)（南風(fēng)）正偏差逐漸增大；OTHERS組成員陸面上經(jīng)向偏差風(fēng)以北風(fēng)為主，偏南風(fēng)較實(shí)況明顯弱（圖9b）。陸面風(fēng)場(chǎng)的差異配置對(duì)低渦強(qiáng)度有重要影響，偏強(qiáng)的低空急流通過夾卷作用在暴雨區(qū)東南側(cè)進(jìn)行動(dòng)量、水汽和能量輸送（李建輝，1982），有助于加劇降水強(qiáng)度。

圖9 850 hPa（a、b、c）和925 hPa（d、e、f）層（a、d）GOOD組成員和再分析場(chǎng)、（b、e）OTHERS組成員和再分析場(chǎng)及（c、f）GOOD組成員和OTHERS組成員的差值風(fēng)場(chǎng)（色階為經(jīng)向風(fēng)偏差，單位：m/s）Fig.9 Differences in wind field between （a，d） GOOD group and reanalysis，（b，e） OTHERS group and reanalysis and（c，f） GOOD group and OTHERS group at 850 hPa （a，b，c） and 925 hPa （d，e，f） levels，respectively（the shadings are for deviations of meridional wind，unit：m/s）

兩組成員陸面上超低空（925 hPa）急流南風(fēng)分量較實(shí)況強(qiáng)的范圍和幅度比低空（850 hPa）急流小。近岸海面上，GOOD組的南風(fēng)預(yù)報(bào)顯著偏強(qiáng)（正偏差），112°E和115°E附近海面存在2個(gè)南風(fēng)偏差核（圖9d），在該核區(qū)位置，OTHERS組同樣存在南風(fēng)偏強(qiáng)的現(xiàn)象，但是幅度較小，2個(gè)正偏差區(qū)有斷裂（圖9e）。換言之，GOOD組超低空偏南急流強(qiáng)度明顯比OTHERS強(qiáng)（圖9f）。超低空急流是季風(fēng)在天氣尺度上的重要載體，其物質(zhì)及能量的輸送方式與低空急流有所不同。本次極端個(gè)例中，模式對(duì)珠江三角洲東部地區(qū)邊界層北風(fēng)的預(yù)報(bào)偏強(qiáng)，GOOD組成員由于偏南急流較實(shí)況強(qiáng)，使得急流前端輻合上升區(qū)維持在珠江三角洲南部，從海上來的暖濕氣流從暴雨底部匯入，提升了低層雨滴濃度，促進(jìn)底層加濕過程，有利于對(duì)流不穩(wěn)定層結(jié)的建立和對(duì)流降水的形成（李建輝，1982；Du，et al，2018）。由于模式同時(shí)存在南、北兩個(gè)方向的風(fēng)速偏差，在偏差耦合作用下加強(qiáng)了低層異常上升運(yùn)動(dòng)，較實(shí)況強(qiáng)的垂直上升運(yùn)動(dòng)導(dǎo)致預(yù)報(bào)降水率的增大，從而使模式預(yù)報(bào)降水量接近極端值的概率增大，降水增幅與偏強(qiáng)的動(dòng)力作用具有一定關(guān)系。而OTHERS組成員在模式偏強(qiáng)北風(fēng)的背景下，輻合區(qū)快速南壓至海上，使得落區(qū)偏移的同時(shí)，降水的累積量也不足。

（超）低空急流有顯著的日變化，通常在夜間增強(qiáng)，早晨達(dá)到最大值，使得華南沿海季風(fēng)降水有夜間多發(fā)的特點(diǎn)（Du，et al，2014）。ERA5再分析場(chǎng)（圖10a、d）顯示，（超）低空急流均有夜間增強(qiáng)的特征，陸上增速中心位于粵東，海上急流加速軸位于海南島東部海面上。925 hPa上，海陸的超低空急流加速中心之間由順暢的增強(qiáng)西南風(fēng)場(chǎng)相連，意味著上游風(fēng)速的日變化增大對(duì)下游的風(fēng)速增大具有一定作用；但下游的急流增速比上游更大，在沿岸出現(xiàn)了加速度輻散區(qū)，以850 hPa層最明顯，這意味著陸地的急流強(qiáng)度變化不完全是由日變化造成的。

從天氣系統(tǒng)的分布配置可見，低空急流與低渦在陸地連接，低空急流的大風(fēng)速區(qū)是低渦環(huán)流的一部分，暖濕氣流的動(dòng)力輻合增強(qiáng)低渦強(qiáng)度，而低渦自身發(fā)展又促進(jìn)了低空急流的加速。在對(duì)應(yīng)位置上，天氣系統(tǒng)之間正反饋?zhàn)饔檬箖烧吖餐l(fā)展。這種作用在850 hPa上更顯著，925 hPa的急流在邊界層的摩擦作用下，加速度有所減小，輻合相對(duì)較弱。

CMA-TRAMS（EPS）兩組成員對(duì)夜間低空急流風(fēng)速變化預(yù)報(bào)的差異較大，導(dǎo)致模式降水的強(qiáng)度和分布出現(xiàn)差異。GOOD組西南急流的加速中心與實(shí)況接近，急流增長(zhǎng)強(qiáng)度（尤其是南風(fēng)分量）略偏強(qiáng)（圖10b、e），體現(xiàn)了CMA-TRAMS在華南地區(qū)常存在的系統(tǒng)性偏強(qiáng)特征（Zhong，et al，2015）。一般而言，環(huán)境風(fēng)偏強(qiáng)使對(duì)流系統(tǒng)移動(dòng)加快，但是本組成員中南北風(fēng)的增速對(duì)比強(qiáng)烈，并且加速區(qū)呈東南—西北向配置，有利于低渦旋轉(zhuǎn)和輻合加強(qiáng)。相比之下，OTHERS組的低空急流風(fēng)速增長(zhǎng)雖也較大，但以偏西分量為主，偏南分量增幅較小，大風(fēng)速軸東移至粵東沿海，陸上的急流增速中心強(qiáng)度顯著偏弱，風(fēng)場(chǎng)輻合強(qiáng)度偏弱，由于急流增長(zhǎng)的動(dòng)力作用導(dǎo)致垂直上升運(yùn)動(dòng)增量較?。▓D10c、f）。

圖10 不同情景下 （a、d.ERA5，b、e.GOOD，c、f.OTHERS） 2020年5月22日02時(shí)和21日20時(shí)風(fēng)速差（a、b、c.850 hPa，d、e、f.925 hPa；單位：m/s）Fig.10 Wind variations （a，b，c.850 hPa；d，e，f.925 hPa；unit：m/s） for 02：00 BT 22 and 20：00 BT 21 May 2020 in different conditions （a，d.ERA5；b，e.GOOD；c，f.OTHERS）

低空急流的強(qiáng)度和輻合程度既受自身日變化的影響，還因增強(qiáng)的低渦環(huán)流反饋而小范圍加速，而低空急流的動(dòng)力輻合作用促進(jìn)低渦的發(fā)展。模式只有能夠較好捕捉到兩系統(tǒng)位置及強(qiáng)度的耦合關(guān)系（如GOOD組），才能較好地預(yù)報(bào)降水落區(qū)及強(qiáng)度。

5.3 中尺度特征的可預(yù)報(bào)性分析

多數(shù)極端降水事件都伴有后向傳播、“列車效應(yīng)”等中尺度對(duì)流組織化特征（徐珺等，2018；曾智琳等，2020；盛杰等，2020）。在本次極端降水過程中，中尺度輻合線在珠江三角洲東北部的長(zhǎng)期維持、強(qiáng)降水單體沿中尺度輻合線后向傳播是導(dǎo)致極端降水發(fā)生的關(guān)鍵中尺度特征。中尺度輻合線的維持依賴兩方面條件：（1）輻合線南側(cè)暖區(qū)的維持利于水平溫度梯度增大。從22日02時(shí)地面溫度場(chǎng)（圖11a）可見，粵東和粵北的平均溫度僅為24—25℃，而珠江口西側(cè)沿海及粵西地區(qū)的溫度超過29℃，珠海至陽(yáng)江一帶的暖中心溫度略高于海面溫度，23°N附近形成一條接近與緯向平行的溫度鋒線。（2）地面冷池的形成及其與暖氣團(tuán)的對(duì)峙促使中尺度鋒生。降水發(fā)生后，珠江三角洲中部地區(qū)的冷舌伸向暖中心。22日03時(shí)廣州市珠江街和東莞市沿江東路形成冷池，兩站與西南或偏南方向站點(diǎn)的水平溫度梯度為0.23—0.76℃/km。其中東莞市沿江東路與官橋滘村局地的強(qiáng)溫度梯度（0.76℃/km）超過2018年華南超歷史極值降水事件中冷池強(qiáng)度（溫度梯度為0.33—0.35℃/km，曾智琳等，2020）的2倍。冷、暖氣團(tuán)的激烈對(duì)峙，有利于鋒區(qū)新生和對(duì)流后向傳播。

圖11 2020年5月22日02時(shí) （a） 地面溫度場(chǎng)、（b） GOOD組溫度場(chǎng)、（c） OTHERS組溫度場(chǎng)和（d） 降水對(duì)地面溫度的敏感度Fig.11 Surface temperature from （a） observations，（b） GOOD group，（c） OTHERS group and（d） precipitation sensitivity to surface temperature at 02：00 BT 22 May 2020

集合敏感性分析結(jié)果表明，CMA-TRAMS（EPS）的降水對(duì)珠江口兩岸的氣溫預(yù)報(bào)有較高的敏感度（圖11d），氣溫越高，越有利于降水增幅，體現(xiàn)了模式降水預(yù)報(bào)對(duì)熱力作用的響應(yīng)。CMA-TRAMS（EPS）所有成員均預(yù)測(cè)粵西冷空氣南下速度偏快（圖9d、e），導(dǎo)致暖脊范圍較實(shí)況明顯偏小，強(qiáng)度偏弱。雖然相比于OTHERS組，GOOD組預(yù)測(cè)珠江口兩岸的氣溫更高，熱力條件更有利，但是仍大幅度低于實(shí)況（圖11b、c）。從海面向陸地伸展的暖舌不僅輸送水汽和能量，而且增強(qiáng)了中尺度鋒區(qū)的斜壓性，對(duì)中尺度天氣系統(tǒng)的維持和降水的增強(qiáng)均有重要作用。而各成員對(duì)珠江口兩岸氣溫預(yù)報(bào)的冷偏差一定程度上導(dǎo)致了該集合預(yù)報(bào)系統(tǒng)對(duì)于極端降水效率的低估。

21日入夜后粵北以平流降溫為主，強(qiáng)降水區(qū)的降溫幅度顯著高于冷平流造成的降溫幅度（圖12a），形成一片孤立于天氣尺度鋒區(qū)的對(duì)流冷池。大氣中的水滴蒸發(fā)吸熱過程是對(duì)流冷池形成的主要機(jī)制，高空干冷空氣夾卷有助于加速該熱力過程。冷池致使地面高壓形成，其風(fēng)場(chǎng)具有多向擴(kuò)散特征，南部和偏西南部均為偏北風(fēng)，而同一時(shí)刻，珠江口西部暖區(qū)受兩支邊界層急流主導(dǎo)，存在西南風(fēng)和東南風(fēng)的輻合，因此在冷池的西側(cè)與南側(cè)，與暖區(qū)的暖濕平流形成風(fēng)向輻合，輻合強(qiáng)度較強(qiáng)（圖12b），對(duì)對(duì)流的觸發(fā)和向后傳播起關(guān)鍵作用。

圖12 2020 年5 月22 日02 時(shí)與21 日20 時(shí)地面溫差（單位：℃）（a中黑框指示冷池所在區(qū)域，范圍與b相同，b疊加22 日02 時(shí)自動(dòng)站點(diǎn)風(fēng)速大于2.5 m/s的風(fēng)場(chǎng)，紅色線為風(fēng)向輻合線；c為GOOD組22 日02 時(shí)的風(fēng)場(chǎng)和地面溫差）Fig.12 Temperature difference between 02：00 BT 22 and 20：00 BT 21 May 2020（unit：℃）（a.the black box indicates the area where the cold pool is located，and the area is the same as that in figure b；b.automatic stations with wind speed greater than 2.5 m/s at 02：00 BT 22 May are superimposed，the red line is the wind direction convergence line；c.wind field and surface temperature differences of the GOOD group at 02：00 BT 22 May 2020）

CMA-TRAMS（EPS）在反饋強(qiáng)降水帶來的冷池特征方面仍存在較大不足。從GOOD組成員的溫度場(chǎng)與風(fēng)場(chǎng)配置（圖12c）可見，模式模擬的暖脊和冷區(qū)主要為環(huán)境冷、暖平流的產(chǎn)物，低渦中心負(fù)變溫幅度小，風(fēng)場(chǎng)也未見中尺度輻合線特征。弱變溫可能是由于冷平流而非中尺度動(dòng)力效應(yīng)造成的。集合敏感度分析（圖11d）顯示，降水對(duì)珠江三角洲北側(cè)的氣溫有一定的負(fù)敏感度，但未能通過顯著性t檢驗(yàn)，意味著大部分成員預(yù)報(bào)的降水與負(fù)變溫并不高度相關(guān)，反映模式模擬對(duì)流系統(tǒng)組織性和中尺度傳播動(dòng)力機(jī)制仍有不足，限制了CMA-TRAMS（EPS）模式對(duì)極端降水的可預(yù)報(bào)能力。

6 結(jié)論和討論

本研究著眼于集合預(yù)報(bào)在極端天氣過程中的應(yīng)用，利用集合敏感性分析和集合組間分析方法對(duì)熱帶中尺度集合預(yù)報(bào)系統(tǒng)（CMA-TRAMS（EPS））在2020年5月22日珠江三角洲地區(qū)一次極端降水天氣過程中的可預(yù)報(bào)性進(jìn)行研究。這次過程發(fā)生在華南前汛期季風(fēng)爆發(fā)的背景下，西南渦向東南方向移至珠江三角洲地區(qū)，與夜間增強(qiáng)的低空急流耦合產(chǎn)生了一次極端強(qiáng)降水天氣。降水過程具有降水效率極端性強(qiáng)、累計(jì)雨量大，造成人員傷亡及經(jīng)濟(jì)損失嚴(yán)重的特點(diǎn)。

CMA-TRAMS（EPS）部分成員臨近時(shí)刻起報(bào)的12 h降水強(qiáng)度和落區(qū)預(yù)報(bào)結(jié)果與實(shí)況較為接近，對(duì)此次極端降水具有一定的預(yù)報(bào)能力。集合敏感性分析表明，低渦和低空急流的發(fā)展是模式降水的關(guān)鍵系統(tǒng)。本次降水與低層溫度和濕度也存在一定的相關(guān)。

基于以上的認(rèn)識(shí)，著眼于CMA-TRAMS（EPS）在天氣尺度系統(tǒng)（低渦和低空急流）和中尺度系統(tǒng)（地面中尺度輻合線、暖脊和冷池）的可預(yù)報(bào)性能進(jìn)行深入討論。基于對(duì)比研究?jī)r(jià)值考慮，根據(jù)主觀對(duì)比和定量檢驗(yàn)結(jié)果，選出“好成員（GOOD）”組，其余成員均歸入“OTHERS”組。其中，GOOD組對(duì)本次過程降水的強(qiáng)度及空間位置的預(yù)報(bào)結(jié)果與實(shí)況更接近，但極端性方面仍有欠缺。該組成員較好的預(yù)報(bào)能力來自于對(duì)低渦和（超）低空急流的演變特征，以及兩者強(qiáng)度和位置耦合的有效預(yù)測(cè)。集合敏感性分析結(jié)果表明南風(fēng)動(dòng)量匯入有利于低渦的緩慢移動(dòng)和增強(qiáng)。利用組間差異分析方法進(jìn)一步探討可知，GOOD組預(yù)報(bào)珠江三角洲東部（超）低空急流南風(fēng)分量較強(qiáng)，不僅有利于低渦緩慢移動(dòng)，并增強(qiáng)了氣旋性輻合，致使降水持續(xù)時(shí)間長(zhǎng)、效率高。而低渦自身發(fā)展又反過來促進(jìn)急流強(qiáng)度的變化。GOOD組既反映出急流強(qiáng)度日變化狀況，并較準(zhǔn)確地刻畫了增強(qiáng)的低渦環(huán)流反饋導(dǎo)致急流小范圍加速的耦合特征。而OTHERS組預(yù)報(bào)的低渦和（超）低空急流的耦合位置偏東偏南，輻合強(qiáng)度偏弱，導(dǎo)致降水強(qiáng)度或落區(qū)出現(xiàn)偏差?？偠灾?dāng)?shù)蜏u和（超）低空急流的位置和強(qiáng)度耦合程度較高時(shí)，模式對(duì)極端天氣的預(yù)報(bào)能力才會(huì)提高。

由于中尺度對(duì)流組織化的模擬能力是能否預(yù)報(bào)極端降水量的關(guān)鍵因素。本個(gè)例中，冷池與其西南或偏南方向上暖區(qū)之間的水平溫度梯度達(dá)0.23—0.76℃/km，激烈的冷、暖氣團(tuán)對(duì)峙加強(qiáng)了中尺度鋒區(qū)和對(duì)流后向傳播。而CMA-TRAMS（EPS）對(duì)珠江口西側(cè)的暖脊強(qiáng)度預(yù)報(bào)偏弱，并且反映的冷區(qū)主要為溫度平流的產(chǎn)物，大部分成員的降水對(duì)負(fù)變溫并未能形成高相關(guān)，CMA-TRAMS（EPS）兩組成員在預(yù)報(bào)中尺度系統(tǒng)的組織性和傳播特征方面均存在明顯不足，從而影響集合預(yù)報(bào)系統(tǒng)對(duì)極端降水的預(yù)報(bào)能力。

提高分辨率可以增強(qiáng)模式對(duì)中尺度系統(tǒng)的預(yù)報(bào)能力，以分辨率達(dá)到3 km的中國(guó)氣象局廣東快速更新同化數(shù)值預(yù)報(bào)系統(tǒng)（CMA-GD）為例，在臨近時(shí)刻，模式地面要素預(yù)報(bào)較好地反映了與強(qiáng)降水區(qū)匹配的冷池、地面輻合線和中尺度低渦等系統(tǒng)的演變特征，該模式預(yù)報(bào)降水效率也較高，達(dá)到140 mm/（3 h）。未來的工作將考慮采用更高分辨率的集合預(yù)報(bào)系統(tǒng)來分析此次強(qiáng)降水過程的中小尺度天氣系統(tǒng)的可預(yù)報(bào)性。此外，文中僅針對(duì)單個(gè)極端降水過程做了集合敏感性及成員對(duì)比分析，而面向不同的天氣流型配置，CMA-TRAMS（EPS）的分析結(jié)果可能出現(xiàn)較大差異，為了得到更普適的結(jié)果，仍需對(duì)更多的個(gè)例進(jìn)行分析。