馬杰 尹姍 金榮花 李勇 任宏昌 梅雙麗
1 國(guó)家氣象中心,北京 100081
2 中國(guó)氣象局—河海大學(xué)水文氣象研究聯(lián)合實(shí)驗(yàn)室,北京 100081
眾所周知,隨著數(shù)值預(yù)報(bào)模式性能和分辨率的不斷提升,現(xiàn)代天氣預(yù)報(bào)越來(lái)越依賴(lài)數(shù)值模式的輸出結(jié)果。目前,歐洲中期天氣預(yù)報(bào)中心(European Centre for Medium Range Weather Forecasting,簡(jiǎn)稱(chēng)ECMWF)和美國(guó)環(huán)境預(yù)報(bào)中心(National Centers for Environmental Prediction,簡(jiǎn)稱(chēng)NCEP)的全球確定性和集合預(yù)報(bào)模式是最先進(jìn)的兩個(gè)主流業(yè)務(wù)模式。前者的確定性預(yù)報(bào)時(shí)效有240小時(shí),水平分辨率0.1°(約10公里)。后者預(yù)報(bào)時(shí)效16天,前10天分辨率13公里,11~16天時(shí)效分辨率34公里。2019年開(kāi)始二者均實(shí)現(xiàn)了每日發(fā)布4次預(yù)報(bào)產(chǎn)品(金榮花等, 2019)。當(dāng)前在我國(guó)國(guó)家級(jí)和省級(jí)業(yè)務(wù)預(yù)報(bào)臺(tái)站中,二者應(yīng)用非常廣泛(李澤椿等, 2014; 薛諶彬等, 2019)。因此,分析總結(jié)其預(yù)報(bào)能力,尤其針對(duì)重大強(qiáng)降雨過(guò)程的檢驗(yàn)評(píng)估,無(wú)疑將有助于提升未來(lái)預(yù)報(bào)準(zhǔn)確率,對(duì)防災(zāi)減災(zāi)和氣象決策服務(wù)均具有十分重要的參考意義。眾多業(yè)務(wù)和科研工作者圍繞實(shí)際業(yè)務(wù)需求,針對(duì)數(shù)值模式檢驗(yàn)評(píng)估也進(jìn)行了大量的研究(吳志鵬等, 2014; 潘留杰等, 2014)。如公穎(2010)從雨帶結(jié)構(gòu)、強(qiáng)度和位置等方面對(duì)多個(gè)模式的定量降水效果進(jìn)行了檢驗(yàn)評(píng)估。王雨和閆之輝(2007)對(duì)比了多個(gè)檢驗(yàn)方案對(duì)降水評(píng)估效果的影響。孫軍等(2018)圍繞北京地區(qū)多次弱降水過(guò)程,利用多個(gè)業(yè)務(wù)數(shù)值模式分析了預(yù)報(bào)偏差及可能原因。
暴雨及其造成的洪澇災(zāi)害是我國(guó)主要自然災(zāi)害之一。每年6月初到7月間的江淮梅汛期,是我國(guó)防汛關(guān)鍵時(shí)段,暴雨頻繁、降水量大,雨日集中(Tao and Chen, 1987; 張德二和王寶貫, 1990; 魏鳳英和張京江, 2004; 王遵婭和丁一匯, 2008; 李勇等,2017)。據(jù)統(tǒng)計(jì),江淮梅雨降水占到了夏季總降雨量的45%左右。研究表明,梅雨雨帶在30°N附近通常呈準(zhǔn)緯向的分布,并具有狹長(zhǎng)的帶狀特點(diǎn)(徐衛(wèi)國(guó)和江靜, 2004; Ding and Chan, 2005; 司東等,2010),而在朝鮮半島和日本地區(qū)雨帶呈東—東北向(丁一匯等, 2007; 蘇同華和薛峰, 2010)。梅汛期雨帶位置的變化和停滯時(shí)間往往決定著江淮地區(qū)的降水空間分布,如當(dāng)雨帶偏南時(shí),即“南澇北旱”,江淮流域會(huì)出現(xiàn)旱情;相反,當(dāng)雨帶北移時(shí),江淮流域會(huì)出現(xiàn)嚴(yán)重洪澇(陳麗娟等, 2019)。因此,準(zhǔn)確預(yù)報(bào)梅雨雨帶的位置在業(yè)務(wù)預(yù)報(bào)中具有重要意義。
2016年6月30日至7月4日,出現(xiàn)了當(dāng)年入汛以來(lái)最強(qiáng)的降雨過(guò)程。自西南地區(qū)東部至長(zhǎng)江中下游地區(qū)連續(xù)數(shù)日普降大到暴雨。以武漢為例,持續(xù)5天的強(qiáng)降雨,使得累計(jì)降雨量達(dá)到了560.5毫米,造成嚴(yán)重的人員傷亡和經(jīng)濟(jì)損失。此次降雨過(guò)程受到了國(guó)務(wù)院防汛部門(mén)以及氣象部門(mén)的高度重視,然而多個(gè)業(yè)務(wù)數(shù)值模式,包括ECMWF和NCEP模式,對(duì)主雨帶位置和降雨強(qiáng)度的預(yù)報(bào),在中短期時(shí)效均出現(xiàn)了明顯偏差。ECMWF連續(xù)多日的確定性和集合預(yù)報(bào),持續(xù)穩(wěn)定表明主雨帶的中心將位于江淮至黃淮一帶,而NCEP模式預(yù)報(bào)卻指出雨帶位置將位于長(zhǎng)江以南地區(qū)。什么原因造成了兩家模式出現(xiàn)如此明顯的差異?本文將在對(duì)比分析兩家模式性能差異基礎(chǔ)上,分析造成這一現(xiàn)象的原因,以期為未來(lái)梅汛期降雨預(yù)報(bào)提供參考。
本文采用了2016年6月24~29日20:00(北京時(shí),下同)起報(bào)的ECMWF和NCEP模式的集合和確定性預(yù)報(bào)資料。集合預(yù)報(bào)數(shù)據(jù)來(lái)自于THORPEX計(jì)劃的TIGGE數(shù)據(jù)集(Bougeault et al.,2010)。前者每個(gè)預(yù)報(bào)時(shí)次對(duì)應(yīng)51個(gè)預(yù)報(bào)成員,后者成員數(shù)為21個(gè)。變量包括分辨率為2°的500 hPa位勢(shì)高度場(chǎng),1°的850 hPa風(fēng)場(chǎng)以及0.5°的降雨場(chǎng)。時(shí)間分辨率為12小時(shí),預(yù)報(bào)時(shí)長(zhǎng)360小時(shí)。此外,還選取ECMWF確定性模式的500 hPa高度場(chǎng)和850 hPa風(fēng)場(chǎng)的初始場(chǎng)作為分析場(chǎng),來(lái)揭示降雨過(guò)程期間環(huán)流形勢(shì)演變特征。實(shí)況降雨采用了中國(guó)氣象局國(guó)家氣象信息中心提供的2016年6月30日至7月4日2472個(gè)臺(tái)站的逐日降雨資料,這里日降水量指的是08:00至第二天08:00的24小時(shí)累計(jì)降雨量。
為定量表征雨帶位置,本文定義了一種雨帶位置指數(shù)。具體計(jì)算方法為:(1)計(jì)算110°~120°E降雨量的緯向平均值;(2)選取25°~35°N范圍內(nèi)降雨量緯向平均極大值所在的緯度作為整體雨帶位置(Yin et al., 2019)。由于ECMWF(NCEP)集合預(yù)報(bào)每個(gè)時(shí)次對(duì)應(yīng)51(21)個(gè)預(yù)報(bào)成員,每個(gè)成員均能計(jì)算得到一個(gè)雨帶位置值,因而所有預(yù)報(bào)成員就組成了一個(gè)雨帶位置序列,即位置指數(shù)。此外,為便于和預(yù)報(bào)資料比較,本文還將站點(diǎn)降雨觀測(cè)資料利用插值方法,插值到了0.5°×0.5°網(wǎng)格上?;谏鲜鲇陰恢玫亩x,計(jì)算得到2016年6月30日至7月4日降雨過(guò)程的實(shí)況雨帶位置為31.5°N。
由圖1可見(jiàn),2016年6月29日至7月4日降雨過(guò)程期間,在500 hPa高度層上的西太平洋副熱帶高壓(簡(jiǎn)稱(chēng)副高)位置較為穩(wěn)定少動(dòng),代表副高的588位勢(shì)什米線持續(xù)分布于江南、華南一帶。而在中高緯度地區(qū),不斷有西風(fēng)槽自西向東傳播,其與東北冷渦相互配合,引導(dǎo)冷空氣分裂南下。在850 hPa層,西南季風(fēng)沿副高西北側(cè)不斷將暖濕空氣向北輸送,與北方持續(xù)南下的偏北氣流交匯形成一穩(wěn)定切變線,并維持在沿長(zhǎng)江一線,使得強(qiáng)降雨不斷產(chǎn)生,最終造成了長(zhǎng)江中下游的強(qiáng)降雨過(guò)程。
圖1 (a-f)2016年6月29日至7月4日20:00 500 hPa高度場(chǎng)(等值線,單位:dagpm)和850 hPa風(fēng)場(chǎng)(風(fēng)矢量,單位:m s-1)分布Fig. 1 Distribution of 500-hPa geopotential height field (isolines, units: dagpm) and 850-hPa wind field (vectors, units: m s-1) at 2000 BJT (Beijing time) from (a-f) June 29 to July 4, 2016
圖2為此次降雨過(guò)程期間的逐日降水量分布。由圖可見(jiàn),主雨帶位于江漢、江淮及江南北部一帶,其中沿長(zhǎng)江地區(qū)出現(xiàn)了大到暴雨,部分地區(qū)大暴雨。不難發(fā)現(xiàn),強(qiáng)降雨中心雖然緩慢東移,但主雨帶位置在南北方向較為穩(wěn)定,造成了多地連續(xù)出現(xiàn)了強(qiáng)降雨天氣,如武漢,在這次降水過(guò)程中連續(xù)4天出現(xiàn)了大到暴雨,引發(fā)了城市內(nèi)澇等災(zāi)害。
圖2 (a-d)2016年7月1~4日逐日降水量(填色,單位:mm)分布及對(duì)應(yīng)緯向平均值(右側(cè)紅色實(shí)線,單位:mm;圖中最右側(cè)標(biāo)注出了緯向平均極大值的位置)Fig. 2 Daily precipitation (shaded; units: mm) and the zonal mean precipitation (red line, units: mm; the postion of zonal mean maximum precipitation is marked on the right ordinate) from (a-d) July 1 to July 4, 2016
2016年梅汛期從6月18日開(kāi)始,截止于7月16日,前后共持續(xù)了29天。圖3給出了整個(gè)梅汛期中期預(yù)報(bào)時(shí)效的主雨帶位置預(yù)報(bào)性能(為更清楚地展現(xiàn)平均差異,該圖使用了9點(diǎn)平滑處理)。由圖可見(jiàn),ECMWF模式預(yù)報(bào)的雨帶位置呈現(xiàn)出穩(wěn)定偏北的特征(周寧芳等, 2018),而NCEP模式則表現(xiàn)出偏北、偏南誤差交替出現(xiàn)(圖3)。即前者誤差較為穩(wěn)定,這一特征在業(yè)務(wù)預(yù)報(bào)中也有利于預(yù)報(bào)員根據(jù)前期模式誤差進(jìn)行下一時(shí)刻的訂正;而后者則無(wú)規(guī)律可言。此外,圖3還揭示出隨著預(yù)報(bào)時(shí)效臨近,ECMWF模式預(yù)報(bào)誤差不斷縮小,雨帶位置不斷接近實(shí)況。而NCEP模式在某些時(shí)段(如6月30日至7月4日)反而隨預(yù)報(bào)時(shí)效臨近誤差呈擴(kuò)大趨勢(shì)。
圖3 梅汛期主雨帶位置的實(shí)況觀測(cè)結(jié)果和(a)ECMWF模式、(b)NCEP模式不同時(shí)效預(yù)報(bào)結(jié)果的對(duì)比Fig. 3 Observed and forecasted positions of (a) ECMWF and (b) NCEP model of rain bands at different leading time in Meiyu period
對(duì)于6月30日至7月4日極端降雨過(guò)程,圖4進(jìn)一步分析了不同預(yù)報(bào)時(shí)效的110°~120°E緯向平均的降雨量分布,從中可以看到降雨強(qiáng)度及雨帶位置的預(yù)報(bào)性能變化特征。由圖4a同樣可見(jiàn),預(yù)報(bào)時(shí)效越長(zhǎng),ECMWF模式預(yù)報(bào)雨帶位置越偏北、雨強(qiáng)越弱。6月24日(提前時(shí)效144小時(shí))預(yù)報(bào)的主雨帶位置北界達(dá)到了34°N附近,即強(qiáng)降雨涵蓋了黃淮地區(qū),較實(shí)況的長(zhǎng)江沿岸(31.5°N)偏北了2~3個(gè)緯距。之后隨著預(yù)報(bào)時(shí)效臨近,雨帶位置雖不斷向南調(diào)整,但直至臨近的24小時(shí)時(shí)效(29日的預(yù)報(bào)),雨帶位置依然偏北了1個(gè)緯距。一方面說(shuō)明ECMWF模式對(duì)此次降雨過(guò)程的可預(yù)報(bào)性較低,另一方面揭示出討論此次降雨過(guò)程誤差的產(chǎn)生,將有助于解釋整個(gè)梅汛期內(nèi)模式預(yù)報(bào)性能較低原因。對(duì)比而言,NCEP模式預(yù)報(bào)的主雨帶位置(圖4b)誤差更小。6月24日給出預(yù)報(bào),雨帶在經(jīng)向上的分布雖較寬,但位置卻較ECMWF偏南,更接近于實(shí)況。自26日后,預(yù)報(bào)的雨帶位置更是持續(xù)位于30°N附近。而對(duì)降雨強(qiáng)度的預(yù)報(bào),NCEP模式雖依然較實(shí)況偏弱,但卻較ECMWF偏強(qiáng),且27日后強(qiáng)度較為穩(wěn)定。綜上所述,針對(duì)此次降雨過(guò)程的雨帶位置和降雨強(qiáng)度的確定性預(yù)報(bào),NCEP模式優(yōu)于ECMWF模式。
圖4 不同時(shí)刻預(yù)報(bào)的緯向平均(110°~120°E)的2016年6月30日至7月4日累計(jì)降水量(單位:mm,陰影)分布:(a)ECMWF模式結(jié)果;(b)NCEP模式結(jié)果。橫坐標(biāo)為不同起報(bào)日期(從左到右對(duì)應(yīng)2016年6月24日20:00至29日20:00)以及實(shí)況(最右側(cè))Fig. 4 Zonal average (110°-120°E) of the accumulated precipitation from June 30 to July 4, 2016, forecasted by (a) the ECMWF model and (b) the NCEP model initialized at different time. The x-axis represents consecutive runs from 2000 BJT 24 to 2000 BJT 29 June 2016 and the observed zonal mean precipitation totals in the right edge
圖5a-b分別為ECMWF和NCEP模式給出的集合平均的6月30日至7月4日過(guò)程累計(jì)降雨量水平分布。比較發(fā)現(xiàn)二者差異明顯。首先,ECMWF模式的雨帶位置較NCEP(圖5b)和實(shí)況(圖5c)均顯著偏北,尤其在115°E以東地區(qū),大于100 mm的雨區(qū)位于黃淮地區(qū),雨區(qū)北界更是伸展到了沿黃河沿岸。而NCEP模式的雨區(qū)偏南,雨帶位于江淮地區(qū)。二者預(yù)報(bào)結(jié)果相差了2~3個(gè)緯距。若與實(shí)況觀測(cè)(圖5c)相比,最大偏差更是達(dá)到了4個(gè)緯距左右。其次,由圖5c的實(shí)況觀測(cè)可見(jiàn),沿長(zhǎng)江中下游地區(qū)的累計(jì)降水量超過(guò)了250 mm,然而兩模式的集合平均強(qiáng)度均明顯偏弱。最后,從雨帶水平分布特征看,ECMWF模式的雨帶更為緊湊,其經(jīng)向范圍更小。同圖4一樣,集合平均的結(jié)果表明了NCEP模式更接近于實(shí)況。
圖5 2016年6月29日20:00起報(bào)的(a)ECMWF和(b)NCEP模式集合平均的6月30日至7月4日過(guò)程累計(jì)降雨量分布與(c)過(guò)程累計(jì)降雨量實(shí)況Fig. 5 Accumulated precipitation forecasted by (a) the ECMWF model, (b) the NCEP model initialized at 2000 BJT June 29, and (c) the observations from June 30 to July 4, 2016
圖6a為ECMWF模式的51個(gè)集合成員預(yù)報(bào)此次降水過(guò)程的雨帶位置。由圖可見(jiàn),大部分成員給出的雨帶位置較實(shí)況依然明顯偏北,只有個(gè)別成員偏南。與之相反,NCEP模式的絕大部分集合成員預(yù)報(bào)的雨帶則偏南(圖6b)。
在圖6a的基礎(chǔ)上,本文對(duì)雨帶位置指數(shù)進(jìn)行了標(biāo)準(zhǔn)化和排序,選取大于兩個(gè)標(biāo)準(zhǔn)差的成員作為偏差組;小于0.24個(gè)標(biāo)準(zhǔn)差的成員作為準(zhǔn)確組。二者分別涵蓋了15個(gè)成員,準(zhǔn)確組為14、38、39、51、4、5、7、12、15、20、27、28、29、43、48,偏差組對(duì)應(yīng)11、24、36、40、41、42、44、47、49、2、13、18、22、32和37號(hào)成員。每日最準(zhǔn)確成員與圖6中整體降雨過(guò)程的準(zhǔn)確成員是不固定的。進(jìn)一步分析了兩組成員的逐日演變差異(圖7a)。不難看到,在逐日變化中ECMWF模式預(yù)報(bào)的雨帶位置在大多數(shù)時(shí)段依然表現(xiàn)為偏北誤差。而對(duì)于NCEP模式,圖7b表明其預(yù)報(bào)的雨帶位置在31°N附近擺動(dòng)。后者與實(shí)況觀測(cè)(圖2)的雨帶位置更為一致。此外,在降雨過(guò)程的大多數(shù)時(shí)段ECMWF模式的偏差組預(yù)報(bào)的雨帶位置較準(zhǔn)確組偏北,且最大偏差出現(xiàn)在7月1~2日,說(shuō)明ECMWF對(duì)整個(gè)降雨過(guò)程的預(yù)報(bào)誤差主要來(lái)源于這兩日(圖7a)。而NCEP模式兩組成員之間的差異不穩(wěn)定,7月1日偏差組預(yù)報(bào)的雨帶位置偏北,而在7月2日和4日則以偏南誤差為主(圖7b)。不難看出,圖7與圖3均揭示出ECMWF模式對(duì)雨帶位置的預(yù)報(bào)呈持續(xù)穩(wěn)定的偏北。而NCEP模式誤差不穩(wěn)定,呈現(xiàn)偏北、偏南交替出現(xiàn)。
圖6 2016年6月29日20:00起報(bào)的(a)ECMWF模式51個(gè)成員和(b)NCEP模式21個(gè)成員集合預(yù)報(bào)的此次降雨過(guò)程的雨帶位置序列Fig. 6 Rain belt location serials of this precipitation event forecasted by (a) 51 members of the ECMWF model and (b) 21 members of the NCEP model initialized at 2000 BJT June 29, 2016
圖8為兩組成員在7月1日的降水量水平分布。由圖8b可見(jiàn),ECMWF偏差組的主雨帶沿黃河下游分布,而準(zhǔn)確組主雨帶分布于江淮地區(qū)(圖8a),二者雨帶位置相差了2~3個(gè)緯距(200~300 km)。此外,偏差組降雨強(qiáng)度更大,圖8b中顯示出成片降雨強(qiáng)度超過(guò)100 mm的大暴雨區(qū)。NCEP模式兩組成員的雨帶位置差異較?。▓D8c-d)。同樣,它的偏差組降雨較準(zhǔn)確組更強(qiáng)。由此,兩集合模式均表明當(dāng)預(yù)報(bào)成員給出的雨帶位置預(yù)報(bào)越偏北,對(duì)應(yīng)的降雨強(qiáng)度將更大。
圖8 2016年7月1日(a、b)ECMWF模式和(c、d)NCEP模式準(zhǔn)確組(左列)與偏差組(右列)成員模擬的降雨量分布。起報(bào)時(shí)間為6月29日20:00Fig. 8 Precipitation on July 1, 2016 forecasted by (a, b) the ECMWF model and (c, d) the NCEP model initialized at 2000 BJT June 29, 2016:(a, c) Correct group, (b, d) error group
在7月2日兩個(gè)模式的預(yù)報(bào)誤差剛好相反,ECMWF模式偏北,而NCEP模式偏南(圖7b)。圖9分別顯示了7月2日兩個(gè)模式的準(zhǔn)確組和偏差組降雨量水平分布。由圖可見(jiàn),相較7月1日,兩個(gè)模式均揭示出了雨帶將緩慢南移趨勢(shì)。同樣,ECMWF模式的偏差組雨強(qiáng)較準(zhǔn)確組略偏大(圖9ab)。而NCEP模式則表明,當(dāng)雨帶偏南時(shí)降雨偏弱。
圖7 2016年6月29日20:00起報(bào)的(a)ECMWF模式和(b)NCEP模式準(zhǔn)確組(點(diǎn)線)、偏差組(短虛線)成員及實(shí)況觀測(cè)(實(shí)線)的雨帶位置逐日變化Fig. 7 Daily rain belt location variation in the correct group (dotted line), error group (dashed line), and observation (solid line), forecasted by (a) the ECMWF model and (b) the NCEP model initialized at 2000 BJT June 29, 2016
圖9 同圖8,但為7月2日Fig. 9 Same as Fig. 8, but for July 2
綜上所述,對(duì)雨帶位置的預(yù)報(bào),雖然NCEP模式總體優(yōu)于ECMWF,但這是由于NCEP模式的誤差呈現(xiàn)偏北、偏南交替出現(xiàn),使得平均的位置誤差較小造成的。在雨帶逐日演變以及誤差穩(wěn)定性方面,ECMWF模式卻明顯好于NCEP模式,其雨帶預(yù)報(bào)呈穩(wěn)定偏北,且降雨強(qiáng)度更接近于實(shí)況。此外,兩家模式的集合成員均揭示出,當(dāng)預(yù)報(bào)成員預(yù)報(bào)的雨帶位置偏北時(shí),降雨強(qiáng)度更大,反之亦然。
圖10計(jì)算了7月1日的位置指數(shù)與多個(gè)時(shí)次500 hPa高度場(chǎng)之間的相關(guān)。圖10d給出了7月1日的雨帶位置指數(shù)與同期500 hPa位勢(shì)高度場(chǎng)所計(jì)算的相關(guān)系數(shù)分布。由圖可知,500 hPa層上的西風(fēng)槽強(qiáng)度與降雨雨帶位置關(guān)系密切。即當(dāng)預(yù)報(bào)的西風(fēng)槽偏強(qiáng)時(shí),雨帶位置偏北,反之亦然(Ma et al.,2019)。圖10a-c為超前時(shí)刻的高度場(chǎng)與7月1日雨帶位置之間的相關(guān)系數(shù),其中圖10a對(duì)應(yīng)初始時(shí)刻。由圖可見(jiàn)兩組成員在青藏高原東北側(cè)(青海到四川西部)的西風(fēng)槽區(qū)域?qū)?yīng)著明顯負(fù)相關(guān)系數(shù)。說(shuō)明初始時(shí)刻(6月29日20:00)高空槽強(qiáng)度的微小偏差與7月1日江淮地區(qū)降雨位置存在密切聯(lián)系。之后,隨著西風(fēng)槽東移出青藏高原(圖10b-c),西風(fēng)槽明顯增強(qiáng)(兩組成員預(yù)報(bào)的西風(fēng)槽區(qū)強(qiáng)度差異以及對(duì)應(yīng)的相關(guān)系數(shù)更加明顯)。當(dāng)西風(fēng)槽移動(dòng)至江淮流域時(shí),造成了降雨落區(qū)與強(qiáng)度之間的差異(圖10d)。
圖10 ECMWF模式預(yù)報(bào)的2016年7月1日雨帶位置指數(shù)與2016年6月(a)29日20:00、(b)30日08:00、(c)30日20:00和(d)7月1日08:00 500 hPa高度場(chǎng)(單位:dagpm)的相關(guān)系數(shù)(填色)分布。藍(lán)色實(shí)線(藍(lán)色虛線)為準(zhǔn)確組(偏差組)成員的500 hPa高度場(chǎng),黑色線為兩組高度差值場(chǎng),起報(bào)時(shí)間為6月29日20:00Fig. 10 Lag correlation coefficients (shaded) between the rain belt location index (RBI) on July 1, 2016 and the 500-hPa geopotential height at different times forecasted by the ECMWF model: (a) 2000 BJT June 29, (b) 0800 BJT June 30, (c) 2000 BJT June 30, and (d) 0800 BJT July 1, 2016.Blue solid (dashed) lines: the 500-hPa geopotential height for the correct (error) group, black contours: difference field between the two groups.Forecast initial time: 2000 BJT 29 June
圖11進(jìn)一步分析了兩組成員預(yù)報(bào)的850 hPa風(fēng)場(chǎng)差異。由圖11a和b可見(jiàn),當(dāng)ECMWF預(yù)報(bào)的西風(fēng)槽偏深時(shí),對(duì)應(yīng)850 hPa層低渦強(qiáng)度偏強(qiáng),表現(xiàn)為低渦東側(cè)切變線南側(cè)(北側(cè))的偏南季風(fēng)(偏北風(fēng))偏強(qiáng),同時(shí)切變線位置也偏北,因而造成了降雨偏強(qiáng)、雨區(qū)位置偏北(圖8a和b)。這一結(jié)論同樣存在于NCEP模式7月1日的預(yù)報(bào)中(圖11c-d和圖8c-d)。反之,當(dāng)西風(fēng)槽預(yù)報(bào)偏弱時(shí),850 hPa層偏南季風(fēng)偏弱、切變線偏南(圖11ef),造成降雨偏弱、雨區(qū)偏南(圖9c-d)。
圖11 (a、b)ECMWF模式預(yù)報(bào)的2016年7月1日08時(shí)500 hPa高度場(chǎng),以及(a)準(zhǔn)確組和(b)偏差組850 hPa風(fēng)場(chǎng)(黑色風(fēng)矢量)分布。藍(lán)色實(shí)線(虛線)對(duì)應(yīng)準(zhǔn)確組(偏差組)高度場(chǎng),紅色實(shí)線(虛線)分別為準(zhǔn)確組(偏差組)副熱帶高壓位置,藍(lán)色粗實(shí)線為其東側(cè)切變線位置,(b)中陰影為偏差組與準(zhǔn)確組經(jīng)向風(fēng)速的差值,起報(bào)時(shí)間為6月29日20:00。(c、d)和(e、f)同(a、b),但分別為NCEP模式7月1日和2日的情形Fig. 11 (a, b) 500-hPa geopotential height (contour) and 850-hPa wind [black wind vector, (a) correct group, (b) error group] on July 1, 2016 forecasted by the ECMWF model initialized at 20000 BJT June 29. The solid (dashed) blue lines: geopotential height for the correct (error) group, the solid (dashed) red line: the position of subtropical high in the correct (error) group, the solid bold blue lines: the position of shear line, shadings in (b):the difference of 850-hPa meridional winds between the correct group and error group. (c) and (d) are the same as (a) and (b), but for the NCEP model on July 1. (e) and (f) are the same as (a) and (b), but for the NCEP model on July 2
綜上所述,此次降雨過(guò)程中西風(fēng)槽預(yù)報(bào)誤差是造成長(zhǎng)江中下游地區(qū)雨帶位置出現(xiàn)誤差的主要原因。當(dāng)西風(fēng)槽偏強(qiáng)時(shí),對(duì)應(yīng)雨帶偏北、降雨強(qiáng)度偏強(qiáng),反之亦然。具體到ECMWF模式,其西風(fēng)槽強(qiáng)度預(yù)報(bào)持續(xù)偏強(qiáng)(高度場(chǎng)誤差表現(xiàn)為持續(xù)的負(fù)值),對(duì)應(yīng)的850 hPa層低渦強(qiáng)度穩(wěn)定偏強(qiáng)、切變線位置偏北,造成了雨區(qū)位置偏北、降雨偏強(qiáng)。而NCEP模式,其7月1~2日預(yù)報(bào)的雨帶位置呈偏北、偏南差異,源于對(duì)西風(fēng)槽的預(yù)報(bào)偏強(qiáng)、偏弱,造成了850 hPa層切變線位置偏北、偏南以及切變線南側(cè)季風(fēng)偏強(qiáng)和偏弱。
2016年6月30日至7月4日降雨是當(dāng)年入汛后最強(qiáng)的一次極端降雨過(guò)程。在實(shí)際業(yè)務(wù)預(yù)報(bào)中,針對(duì)其雨帶位置的預(yù)報(bào),無(wú)論中期還是臨近時(shí)效均出現(xiàn)了較明顯的預(yù)報(bào)誤差。為此,本文基于ECMWF和NCEP模式的確定和集合預(yù)報(bào)資料,比較二者的預(yù)報(bào)性能差異及其出現(xiàn)的原因。得出以下結(jié)論:
(1)針對(duì)此次降雨過(guò)程整體雨帶位置的預(yù)報(bào),雖然NCEP的確定性和集合平均結(jié)果較ECMWF更為準(zhǔn)確。但進(jìn)一步分析卻表明,NCEP預(yù)報(bào)誤差在逐日變化中呈現(xiàn)偏北、偏南交替出現(xiàn),因而造成了平均的預(yù)報(bào)誤差較小。而ECMWF模式無(wú)論確定性模式還是大部分集合成員,雨帶預(yù)報(bào)呈現(xiàn)穩(wěn)定偏北誤差。
(2)基于集合預(yù)報(bào)成員的分組對(duì)比分析表明,當(dāng)雨帶位置預(yù)報(bào)偏北時(shí),對(duì)應(yīng)降雨強(qiáng)度偏強(qiáng),反之亦然。這一結(jié)論在NCEP模式和ECMWF模式中均存在。
(3)500 hPa高度場(chǎng)的西風(fēng)槽預(yù)報(bào)差異是造成雨帶位置誤差的重要原因。當(dāng)西風(fēng)槽預(yù)報(bào)偏強(qiáng)時(shí),雨帶位置偏北、降雨強(qiáng)度偏強(qiáng)。反之,當(dāng)西風(fēng)槽強(qiáng)度偏弱時(shí),雨帶位置偏南、降雨強(qiáng)度偏弱。NCEP模式西風(fēng)槽強(qiáng)度預(yù)報(bào)偏強(qiáng)、偏弱交替出現(xiàn),對(duì)應(yīng)了其雨帶位置偏北、偏南,而ECMWF模式西風(fēng)槽預(yù)報(bào)持續(xù)偏強(qiáng),造成了雨帶位置的預(yù)報(bào)穩(wěn)定偏北。
雖然本工作只討論了2016年6月30日至7月4日的降雨過(guò)程,但分析中利用了ECMWF和NCEP集合模式的連續(xù)多個(gè)預(yù)報(bào)時(shí)次資料,且前者有51個(gè)成員、后者有21個(gè)成員,因此結(jié)論具有一定的代表性。ECMWF模式預(yù)報(bào)的穩(wěn)定性,使得其在業(yè)務(wù)預(yù)報(bào)中更受歡迎。而NCEP模式的預(yù)報(bào)雖存在明顯不穩(wěn)定現(xiàn)象,使得基于其進(jìn)行模式檢評(píng)估和訂正困難較大。但另一方面,其集合平均結(jié)果卻明顯優(yōu)于ECMWF模式預(yù)報(bào),更接近于實(shí)況觀測(cè),揭示在未來(lái)的預(yù)報(bào)實(shí)踐中,還需要進(jìn)一步研究和討論二者的具體差異。此外,本工作還表明數(shù)值模式初始場(chǎng)上500 hPa層西風(fēng)槽微小的差異在造成雨帶位置和強(qiáng)度誤差的重要作用。當(dāng)其東移出青藏高原時(shí),西風(fēng)槽強(qiáng)度及對(duì)應(yīng)的誤差同時(shí)突然增強(qiáng),是什么原因(如對(duì)流潛熱或地形影響)造成了這一現(xiàn)象?除了初始場(chǎng)誤差外,數(shù)值模式自身預(yù)報(bào)誤差,對(duì)流參數(shù)化、潛熱釋放等因素,造成了什么樣的影響,這些均有待于進(jìn)一步研究。