一次燕山南麓局地大暴雨的成因分析及數(shù)值預(yù)報檢驗

崔粉娥 暢巨崢2 燕成玉 張晨宇

（1 河北省秦皇島市氣象局，秦皇島 066001；2 河北省青龍滿族自治縣氣象局,秦皇島 066500）

0 引言

致災(zāi)暴雨一直是廣大氣象工作者關(guān)注的重點[1-3]。大范圍對流性暴雨多數(shù)是由中尺度對流系統(tǒng)或中尺度對流復(fù)合體而產(chǎn)生[4-5]，由對流風暴引起的強天氣很容易漏報而成為眾多學(xué)者關(guān)注和研究的重點[6-16]。而對數(shù)值預(yù)報產(chǎn)品進行評估、檢驗，是數(shù)值預(yù)報產(chǎn)品在實際預(yù)報業(yè)務(wù)中有效應(yīng)用，提高天氣預(yù)報準確率的途徑之一。許多學(xué)者從不同角度做了大量的檢驗工作[17-27]，孟英杰等[17]認為，在武漢區(qū)域降水預(yù)報中，JAPAN模式的綜合預(yù)報性能好于T213、AREM、GRE模式。肖明靜等[19]對汛期山東降水預(yù)報進行了多模式檢驗，結(jié)果認為一般性降水和小雨預(yù)報，24 h、48 h降水MM5模式評分結(jié)果均為最差，T639模式預(yù)報效果最好；中雨及以上量級，兩個區(qū)域模式均不如T639模式預(yù)報效果好。蔡薌寧[26]認為ECMWF模式對各影響系統(tǒng)及要素的預(yù)報最接近實況JAPAN模式和T639模式次之。程錦霞[27]對主要物理量因子的T639預(yù)報場進行檢驗，結(jié)果認為垂直速度預(yù)報最優(yōu)，其次是渦度、散度、渦度平流預(yù)報，水汽通量、假相當位溫24～72 h預(yù)報產(chǎn)品參考價值較大。本文通過對基層氣象臺（站）常用的幾家數(shù)值模式在一次局地大暴雨過程中的預(yù)報情況進行檢驗，為預(yù)報工作者正確使用數(shù)值模式產(chǎn)品提供依據(jù)。

2013年6月28—29日，秦皇島市各縣區(qū)出現(xiàn)強降雨天氣，降水分布不均勻，北部青龍為大暴雨，盧龍為暴雨，昌黎、撫寧為大雨，其他地區(qū)為中雨，并伴有較強雷電。5個區(qū)域氣象站降水量超過100 mm，26個區(qū)域氣象站超過50 mm，最大出現(xiàn)在青龍土門子為132.6 mm。強降水時段出現(xiàn)在28日22時至29日00時，小時雨強最大為59.7 mm，3 h累計降水量超過100 mm，達短時暴雨標準。本次過程降水時段集中，降水強度大，且強降水出現(xiàn)在北部山區(qū)，跟燕山南麓地形關(guān)系較大。

1 環(huán)流背景及影響系統(tǒng)分析

500 hPa環(huán)流形勢中高緯度呈一槽一脊型，2013年6月27日08時貝加爾湖附近形成一低渦，在鄂霍次克海存在一阻塞高壓，該低渦系統(tǒng)穩(wěn)定少動，強度不斷加強，到28日20時低渦中心強度從544 dagpm增強至542 dagpm，并配合有-20℃的冷溫度中心。由于渦后部冷空氣不斷分裂南下，在28日20時形成階梯槽，一支位于貝加爾湖南部，另一支位于華北北部，這種形勢有利于冷空氣補充南下影響秦皇島地區(qū)。在相應(yīng)的200 hPa上，秦皇島地區(qū)受高壓脊控制，脊前呈疏散狀，表明高空存在輻散，有利于中低層氣流輻合上升。

在28日20時的700 hPa和850 hPa上，秦皇島受西南氣流和溫度脊控制，700 hPa在河北西北部至內(nèi)蒙古一帶溫度露點差小于2℃，說明水汽充沛。850 hPa在河北、山西、內(nèi)蒙古交界處存在一切變。

本次強降水過程的影響系統(tǒng)是高空短波槽和東蒙低渦。由中尺度環(huán)境場（圖1）可知：由于東蒙冷渦穩(wěn)定少動，且不斷發(fā)展加強，冷空氣不斷分裂南下影響河北東北部地區(qū)。秦皇島位于高空冷溫度槽和低空暖脊之間，且為K指數(shù)大值區(qū)，該區(qū)域為不穩(wěn)定區(qū)域，有利于產(chǎn)生對流活動。強降水落區(qū)位于高空急流出口右側(cè)和低空急流入口左側(cè)，850 hPa顯著濕區(qū)，且K指數(shù)大于32 ℃。切變線和地面輻合線觸發(fā)了秦皇島北部短時強降水。

圖1 2013年6月28日20時高低空系統(tǒng)配置圖（其中紫色線條代表K指數(shù)場，綠色區(qū)域為暴雨區(qū)）Fig. 1 Configuration of high and low weather system at 2000 BT on 28 June 2013(The purple line represents K exponential field, and the green area represents the rainstorm area)

由圖2可知，28日20時秦皇島受低壓控制，輻合線位于河北和內(nèi)蒙交界，青龍站小時降水量為1.9 mm。到23時隨著系統(tǒng)東移和冷空氣補充，輻合區(qū)位于河北北部，并在河北東部產(chǎn)生新的輻合區(qū)，河北西北部處于氣旋性曲率大值區(qū)，青龍站小時雨量達到32.2 mm，說明該區(qū)域輻合強，有利于對流活動的發(fā)展。到29日02時河北東北部轉(zhuǎn)為反氣旋環(huán)流，秦皇島北部受東北氣流控制，降水明顯減弱，04時降水基本結(jié)束。

2 中尺度特征分析

2.1 衛(wèi)星云圖演變分析

從紅外云圖動態(tài)分析可知：28日08:30在河北西北部有一低渦云系，該云系范圍較寬廣（圖3a）。由于高空冷空氣補充南下及日本海阻塞高壓的存在，有利于對流云團加強發(fā)展和維持，在22:30低渦云系發(fā)展為中尺度對流復(fù)合體MCC（圖3b），該云系結(jié)構(gòu)密實，青龍縣區(qū)位于TBB梯度最大處，該時刻與最大小時降水量相吻合（圖3e）。隨著系統(tǒng)東移和低層輻合減弱，從00:30開始，對流云團東南部的云系逐漸減弱與主體分離，云團結(jié)構(gòu)開始松散（圖3c），至03:30對流云團基本移出秦皇島（圖3d），強降水明顯減弱，后期轉(zhuǎn)為分散的降水云系影響，降水趨于結(jié)束（圖3e）。

圖2 2013年6月28日20時和23時加密站風場疊加小時雨量（紅色區(qū)域代表青龍站）Fig. 2 Wind field and hourly rainfall (The red area represents Qinglong station)(a) 2000 BT on 28 June; (b) 2300 BT on 28 June 2013

2.2 多普勒雷達組合反射率演變分析

本次過程的強降水時段出現(xiàn)在28日22時至29日04時，由組合反射率因子的動態(tài)變化（圖4）可知：28日20時青龍上空開始有弱降水回波，在興隆到寬城一帶存在一窄帶回波，大片降水回波中鑲嵌著超級單體，中心強度值為50～60 dBz，該回波帶緩慢向東南方向移動。位于秦皇島上空的回波帶呈線狀，線狀回波中不斷有超級單體影響秦皇島市北部縣區(qū)，類似“列車效應(yīng)”。21時青龍上空降水回波中心強度最大值為60 dBz，由于秦皇島市西北部位于切變線附近，且低層濕度條件較好，并且在北京—天津一帶不斷有帶狀降水回波生成，使該回波帶經(jīng)過降水區(qū)時回波中心強度值穩(wěn)定維持在55～65 dBz，青龍最大小時降水量達到44.3 mm。29日02時降水回波的線狀結(jié)構(gòu)開始松散，秦皇島市轉(zhuǎn)為零散的降水回波影響，該降水回波影響秦皇島市長達4～5 h，從而造成北部山區(qū)暴雨到大暴雨。04時以后降水明顯減弱，降水趨于結(jié)束。

2.3 閃電資料及自動站資料分析

分析強降水時段基本徑向速度圖和閃電實況圖（圖略）可知：秦皇島市降水出現(xiàn)在28日20時至29日04時，北部山區(qū)為雷電極高密度區(qū)域。速度圖在低空有明顯的東南氣流輸送水汽，青龍縣區(qū)附近低層為輻合，高層為輻散，有利于上升運動的維持，且中低層為大范圍正速度區(qū)包圍有負速度區(qū)，有“逆風區(qū)”存在，該區(qū)域與強降水區(qū)域吻合。5 min加密自動站風場顯示（圖略）：28日21：30青龍縣區(qū)有明顯的風場輻合，29日01：30青龍轉(zhuǎn)為南到東南風，中尺度風場輻合減弱、消失，降水開始明顯減弱。

3 物理量診斷分析

3.1 水汽條件分析

分析可知：水汽通量散度輻合區(qū)位于河北東部，最大值為-3×10-6g·s-1，同時次的水汽通量場顯示有東南和西南兩支氣流輸送水汽，水汽主要來源于渤海和黃海（圖5a）。850 hPa比濕大于12 g·kg-1，這充分說明中低層充沛的水汽條件，對降水起到增幅作用，有利于產(chǎn)生短時強降水（圖5b）。

3.2 熱力條件分析

相當位溫θe是一個重要的溫濕特征參數(shù)，由相當位溫的徑向剖面圖（圖6a）可知：秦皇島上空在900 hPa附近存在一閉合暖區(qū)，中心值達342 K，且850 hPa以下為位勢不穩(wěn)定區(qū)域，與徑向分量大正值區(qū)相對應(yīng)，偏南風有利于暖濕氣流輸送。

在T-lnP圖上（圖6b），濕層較厚，高濕區(qū)達700 hPa高度，925 hPa 存在一淺薄的逆溫層，并且存在風切變。K指數(shù)為34，SI指數(shù)為-9.12，CAPE值為2766.1，整層比濕積分為4682.7，這說明大氣為層結(jié)不穩(wěn)定狀態(tài)，蘊含著不穩(wěn)定能量，為強降水的產(chǎn)生提供了不穩(wěn)定條件。同時，-20 ℃層高度為5.8 km，0 ℃層高度為4.3 km，但由于中低層濕層深厚，且風切變小，不滿足大風、冰雹等強對流天氣的條件。

圖3 2013年6月28—29日紅外云圖（其中白圓圈代表秦皇島）：（a）28日08:30，（b）28日22:30，（c）29日00:30,（d）29日03:30和（e）青龍站逐小時降水量Fig. 3 Infrared cloud image (a) 0830 BT on 28 June, (b) 2230 BT on 28 June, (c) 0030 BT on 29 June and (d) 0330 BT on 29 June 2013 and hourly precipitation (e) at Qinglong station (The white circle represents Qinhuangdao)

3.3 動力條件分析

從垂直速度剖面可以看出秦皇島北部地區(qū)600 hPa以下為上升運動區(qū)，最大上升速度中心值為-40×10-3m·s-1（圖7a），對應(yīng)散度場在300 hPa以下為輻合區(qū)（圖7b），中心值為-16×10-6s-1，高層為輻散區(qū)，中心值為24×10-6s-1，正散度明顯大于負散度，說明高空急流抽吸作用明顯，有利于降水增強。

4 數(shù)值預(yù)報產(chǎn)品檢驗分析

4.1 環(huán)流形勢檢驗（T639模式/EC模式）

圖4 2013年6月28日21時—29日02時組合反射率（單位：dBz）與小時降水量（單位：mm）配置圖（a）28日2100 BT；（b）28日2200 BT；（c）28日2300 BT；（d）29日0000 BT；（e）29日0100 BT；（f）29日0200 BTFig. 4 Combined reflectance (unit: dBZ) and hourly precipitation (unit: mm) from 2100BST 28th to 0200BST 29th(a) 2100BST on 28 June; (b) 2200BST on 28 June; (c) 2300BST on 28 June; (d) 0000BST on 29 June; (e) 0100BST on 29 June; (f) 0200BST on 29 June

圖5 2013年6月28日20時850 hPa（a）水汽通量散度（單位：10-6 g·s-1）和（b）比濕場（單位：g·kg-1）Fig. 5 Moisture flux divergence field (units: 10-6 g·s-1) and specific humidity (units: g·kg-1) at 850 hPa at 2000 BT on 28 Jun 2013

對比EC、T639模式500 hPa高度場預(yù)報穩(wěn)定性（圖8），可以看出：對東蒙冷渦系統(tǒng)，EC模式和T639模式各時效均預(yù)報出544 dagpm低渦中心，EC模式和T639模式24 h、48 h預(yù)報接近實況場，72 h預(yù)報EC模式較實況場略偏強，T639模式較實況場略偏弱；對副熱帶高壓584 dagpm特征線時次的預(yù)報，EC模式的預(yù)報比較穩(wěn)定，跟實況基本一致，T639模式誤差較大。對日本海阻塞高壓T639模式各時效均預(yù)報偏弱，且時效越長穩(wěn)定性越差，EC模式各時效預(yù)報效果較好。

4.2 850 hPa風場檢驗

對影響華北東北部的低空急流系統(tǒng)（圖9），兩家模式24 h時效預(yù)報的風速均較實況偏弱，對850 hPa切變線預(yù)報均偏弱、偏北，EC模式對秦皇島北部的風向預(yù)報偏差較大。

4.3 降水強度及落區(qū)檢驗

對降水落區(qū)進行檢驗（圖10），本次強降水過程，各家數(shù)值預(yù)報對雨帶的走向把握較好，從預(yù)報量級上看，均存在預(yù)報偏小的情況，特別是對大暴雨的落區(qū)預(yù)報效果均不理想，穩(wěn)定性差。對華北地區(qū)東北部中雨以上量級的降水，各家模式預(yù)報均偏北，強度偏弱。各家指導(dǎo)產(chǎn)品強降水中心為大雨到暴雨，比較一致，但強降水中心都偏北，且未能體現(xiàn)出大于100 mm的大暴雨區(qū)域。T639預(yù)報偏差較大，量級為大雨并且中心在河北張家口一帶，偏差較大。日本細網(wǎng)格和省臺指導(dǎo)預(yù)報對強降水落區(qū)偏差較小，中心略偏北，更接近實況。

圖6 2013年6月28日20時（a）相當位溫沿119°E徑向剖面（單位：K）疊加風速徑向分量和（b）北京探空圖（黑色：相當位溫；紅色：風速徑向分量）Fig. 6 The vertical cross-section of potential equivalent temperature (units: K) along 119°E with the radial component of the wind velocity (a) and sounding figure from Beijing station at 2000 BT on 28 June 2013 (b)(The black represents potential equivalent temperature; the red represents radial component of the wind velocity)

圖7 2013年6月28日20時沿119°E徑向剖面的（a）垂直速度（單位：m·s-1）和（b）散度（單位：s-1）Fig. 7 The vertical cross-section of vertical velocity (a, units: m/s) and divergence (b, units: s-1) along 119°E at 2000 BT on 28 June 2013

圖8 2013年6月6月28日20時500 hPa高度實況場及對應(yīng)24 h、48 h、72 h預(yù)報（a）EC模式；（b）T639模式（其中黑色為實況場，紅色表示24 h預(yù)報，黃色表示48 h預(yù)報，藍色表示72 h預(yù)報）Fig. 8 The actual field at 500 hPa and the corresponding forecasts at 24 h, 48 h and 72 h: (a) EC model;(b) T639 model (black: actual field, red: 24 h forecast, yellow: 48 h forecast, blue: 72 h forecast)

圖9 2013年6月28日20時850 hPa風場實況及EC（a）、T639（b）模式24 h時效預(yù)報圖（其中黑色表示觀測值，紅色表示預(yù)報值大于觀測值，藍色表示預(yù)報值小于觀測值）Fig. 9 Wind field at 850 hPa and 24 h forecast from EC (a) and T639 (b) models (black: observed value, red:Larger than observed values, blue: Less than observed values)

圖10 2013年6月28日08時至6月29月08時累積降水量預(yù)報及實況（單位：mm）（a）T639；（b）EC；（c）日本細網(wǎng)格；（d）日本傳真圖FAX；（e）GRAPES；（f）省臺預(yù)報；（g）實況Fig. 10 Forecast and reality of cumulative precipitation (unit: mm) from 0800 BT on 28 June to 0800 BT on 29 June 2013(a) T639；(b) EC；(c) Jap-thin；(d) FAX；(e) GRAPES；(f) Provincial forecast；(g) Actual

上圖的基本結(jié)論是：各家預(yù)報均預(yù)報了一次大到暴雨過程（即雨量在25～50 mm），只有省臺綜合預(yù)報為暴雨到大暴雨過程（即雨量為50～100 mm），但位置普遍偏北?？傊骷夷Ｊ綄@次過程的預(yù)報能力基本相同。造成位置偏北的原因可能與大尺度流場的預(yù)報存在缺陷有關(guān)。

整體而言，EC模式的預(yù)報對亞洲地區(qū)中高空環(huán)流預(yù)報的結(jié)果更為準確，對弱高空槽系統(tǒng)預(yù)報偏差較大。日本細網(wǎng)格和省臺指導(dǎo)預(yù)報對暴雨落區(qū)的預(yù)報偏差較小，而降水極值預(yù)報需結(jié)合強對流工具、實況加密資料及時更新預(yù)報，發(fā)布短臨預(yù)報。

冷渦天氣形勢下，各種數(shù)值預(yù)報模式均預(yù)報有降水，且降水量級一致或相差不多時，考慮地形作用，秦皇島北部山區(qū)降水量級可增加一個量級或稍大一些；當細網(wǎng)格預(yù)報有大雨，數(shù)值接近暴雨時，應(yīng)考慮有暴雨或局地暴雨出現(xiàn)。

5 結(jié)論

（1）本次暴雨過程是由于東蒙冷渦穩(wěn)定少動，影響系統(tǒng)是高空低渦和高空槽，觸發(fā)機制是切變線和中尺度風場輻合。強降水落區(qū)位于水汽通量大值區(qū)，600 hPa以下強上升運動區(qū)和850 hPa以下位勢不穩(wěn)定區(qū)域。

（2）對500 hPa中高緯環(huán)流場、槽線進行檢驗，對東蒙冷渦預(yù)報T639和EC兩模式預(yù)報都比較穩(wěn)定，但對日本海阻高預(yù)報T639預(yù)報偏弱，且時效越長穩(wěn)定性越差。EC模式對弱高空槽系統(tǒng)預(yù)報偏差較大，需重點關(guān)注。

（3）降水強度和落區(qū)檢驗：各家數(shù)值模式預(yù)報產(chǎn)品對雨帶走向把握較好，但對降水的強度把握均較實況偏弱、偏北。在降水量級上，日本細網(wǎng)格和省臺指導(dǎo)預(yù)報偏差較小。

綜上所述，檢驗結(jié)果由于其局限性并不能全面精確地評述每種模式的預(yù)報性能，只能從典型個例給相似天氣過程的模式預(yù)報解釋應(yīng)用提供一個參考。