1. 引言
云南省地处低纬高原,地形复杂多样,局地暴雨频发且具有突发性强、历时短、雨量集中等特点,常引发山洪、滑坡、泥石流及城市内涝等次生灾害,对人民生命财产安全造成严重威胁[1]。马关县位于云南省东南部,地处哀牢山余脉与滇东南喀斯特地貌过渡带,受复杂地形与季风环流共同影响,暴雨事件频发且局地特征显著。例如,2015年9月27日马关县遭遇特大暴雨,单日降水量突破历史极值,导致交通中断、农田损毁,凸显了该地区暴雨灾害的严峻性。然而,现有研究多聚焦于云南典型区域(如昆明)的暴雨过程,针对马关地区暴雨的系统性研究仍较匮乏,其形成机制、中小尺度特征及地形效应尚未明确,制约了精细化预报能力的提升[2]-[5]。
已有研究表明,云南夏季常受孟加拉湾风暴、西行台风、副热带高压、高原低涡、切变线等天气系统的影响,产生暴雨[1] [6]。例如,副热带高压外围型暴雨常伴随西南或东南水汽输送,而高原低涡切变则通过动力抬升触发强对流[7]。此外,地形对暴雨的增幅作用显著,如昆明主城区的喇叭口地形可增强辐合并延长降水时间[8]。
马关地区地处滇东南迎风坡,地形坡度陡峭,可能通过强迫抬升机制加剧水汽辐合,但其具体作用机制仍需量化分析。同时,数值模拟显示,暴雨过程中低层水汽通量散度与垂直上升运动的时空配置对雨强分布具有关键指示意义[4] [9],而马关地区在此方面的研究尚属空白。鉴于此,本研究以马关地区暴雨为对象,分析三次不同背景下暴雨过程的环流形势、形成原因及雷达特征,为提升该地区暴雨预报精度和防灾减灾能力提供科学依据。
2. 资料与方法
本文采用的资料主要包括:(1) 马关市地面常规观测站和加密自动站逐小时的温、压、湿、风等常规气象要素观测资料;(2) ERA5 0.25˚ × 0.25˚逐小时的再分析资料;(3) 蒙自探空站的探空资料;(4) 文山S波段多普勒雷达的产品资料。
3. 天气实况与背景分析
3.1. 天气实况
个例1,2016年第3号台风“银河”于7月28日00:20在越南北部南定省沿海登陆,登陆时中心附近最大风力达10级,中心最低气压982 hPa。登陆后,“银河”减弱为热带风暴级。28日08时台风中心位于越南北部,此时马关位于台风的第二象限。随着台风西行、受内陆摩擦减弱影响,台风逐渐减弱为热带低压,20时低压在马关打转。受其影响,马关县28日08时至29日20时发生了一次台风型暴雨过程。受台风“银河”残余低压影响,过程累积雨量达137.8 mm,其中28日08时至29日08时24小时累积雨量达80.8 mm (图1(a)),29日08时至30日08时24小时累积雨量达66.5 mm (图1(b)),最大小时雨强为22.9 mm,出现在29日09时~10时(图2(a))。
个例2,2017年6月19日08时至20日20时马关地区出现暴雨过程,过程累积雨量达113.9 mm,19日08时至20日08时24小时累积雨量达100.0 mm (图1(c)),最大小时雨强为33.9 mm,出现在19日18时(图2(b))。
个例3,2018年6月23日08时~25日08时马关过程累积雨量达135.4 mm,其中24日08时至25日08时24小时累积雨量达88.2 mm (图1(d)),最大小时雨强为18.1 mm,出现在24日08时(图2(c))。
Figure 1. 24-hour accumulated precipitation from MICAPS (shaded, unit: mm). (a) Case 1: 08:00 BST 28 July to 08:00 BST 29 July 2016; (b) Case 1: 08:00 BST 29 July to 08:00 BST 30 July 2016; (c) Case 2: 08:00 BST 19 June to 08:00 BST 20 June 2017; (d) Case 3: 08:00 BST 24 June to 08:00 BST 25 June 2018
图1. MICAPS 24小时累积降水(阴影,单位:mm)。(a) 个例1 2016年7月28日08时至29日08时;(b) 个例1 2016年7月29日08时至30日08时;(c) 个例2 2017年6月19日08时至20日08时;(d) 个例3 2018年6月24日08时至25日08时
Figure 2. Hourly precipitation variation (unit: mm) of Maguan observe station. (a) Case 1; (b) Case 2; (c) Case 3
图2. 马关站1小时降水量(单位:mm)随时间的变化。(a) 个例1;(b) 个例2;(c) 个例3
3.2. 大气形势
个例1中,500 hPa (图3(a)、图3(b))28日08时亚欧中高纬两槽一脊,乌拉尔山东部和内蒙古东部为槽,并伴有低涡,贝加尔湖东部为一脊。中纬度地区有两大高压影响,分别是伊朗高压和西太平洋副热带高压,两高之间以均压场为主,引导气流弱,有一台风低压影响越南北部。2908中高纬度两槽一脊东移。中低纬度副高增强,云南受两高之间均压场控制,马关有风场的辐合。700 hPa (图4(a)、图4(b)) 28日08时云南省主要位于台风低压第二象限,随着台风低压继续西行,马关逐渐受第一象限影响。至2908时,低压中心移出文山,马关位于低压东侧切变线附近。
Figure 3. ERA5 reanalysis data at 500 hPa: wind field (vectors, unit: m/s), temperature (red contours, unit: ˚C) and geopotential height field (blue contours, unit: dagpm). 2016 July: (a) 08:00 BST 28 July, (b) 08:00 BST 29 July; 2017 June: (c) 08:00 BST 19 June, (d) 20:00 BST 19 June; 2018 June: (e) 20:00 BST 23 June, (f) 08:00 BST 24 June
图3. ERA5再分析资料500 hPa风场(矢量,单位:m/s)、温度场(红色等值线,单位:℃)和位势高度场(蓝色等值线,单位:dagpm)。2016年7月:(a) 28日08时,(b) 29日08时;2017年6月:(c) 19日08时,(d) 19日20时;2018年6月:(e) 23日20时,(f) 24日08时
个例2中,500 hPa (图3(c)、图3(d)) 19日08时中高纬两槽一脊,乌拉尔山东部和东北北部均有低涡,东北北部低涡中心气压低于560 dagpm,贝加尔湖西部为一脊,引导冷空气南下影响中纬度地区。马关位于588 dagpm和584 dagpm之间,有短波槽影响。19日20时两槽一脊稳定少动,脊前冷空气南下位置更深。低纬短波槽加深,副高略有增强。脊前冷空气向南输送,易在副高北侧高能高湿区激发对流。700 hPa (图4(c)、图4(d)) 19日08时马关上游有弱的切变,19日20时马关附近切变加深。切变线位于低涡底部,有利于弱冷空向南渗透,有利于激发对流[6]。
个例3中,500 hPa (图3(e)、图3(f)) 23日20时高纬两槽两脊,中纬度多短波槽活动。低纬分别受伊朗高压和西太平洋副热带高压影响,副高脊线位置约25˚N,马关处于两高之间,受高原涡影响。24日08时高纬槽脊位置稳定少动,短波槽向南加深与高原涡合并,西太副高西移。后期副高持续加强西伸,阻挡短波槽东移,马关长时间处于短波槽前,为降水持续提供了非常有利的大尺度环流形势。期间700 hPa (图4(e)、图4(f))马关受低涡切变影响,切变线位置稳定,马关处于切变线前部上升气流区,其上不断有高原涡生成后消失。
Figure 4. ERA5 reanalysis data at 700 hPa: wind field (vectors; color shading for wind speeds ≥12 m/s, unit: m/s), temperature (red contours, unit: ˚C) and geopotential height field (blue contours, unit: dagpm). 2016 July: (a) 08:00 BST 28 July, (b) 08:00 BST 29 July; 2017 June: (c) 08:00 BST 19 June, (d) 20:00 BST 19 June; 2018 June: (e) 20:00 BST 23 June, (f) 08:00 BST 24 June
图4. ERA5再分析资料700 hPa风场(矢量,填色为12 m/s以上风速,单位:m/s)、温度场(红色等值线,单位:℃)和位势高度场(蓝色等值线,单位:dagpm)分布。2016年7月:(a) 28日08时,(b) 29日08时;2017年6月:(c) 19日08时,(d) 19日20时;2018年6月:(e) 23日20时,(f) 24日08时
3.3. 探空分析
分析个例1中蒙自探空站29日08时T-lnP图可知(图5(a)),温度露点差小,水汽条件好,大气湿层深厚,500 hPa以下大气近似饱和。抬升凝结高度低接近地面,CAPE值为338.9 J/Kg,CIN为0 J/Kg,说明完全不需要克服外力就可以自由对流,K指数较大为39℃,SI指数<0℃,有利于短时强降水的出现。垂直风切变较小,不利于对流大风出现。
分析个例2中蒙自探空站19日08时T-lnP图可知(图5(b)),大气湿层深厚,600 hPa以下大气近似饱和,400 hPa以上有干层。抬升凝结高度约在820 hPa,CAPE值为375 J/Kg,CIN为18.4 J/Kg,说明需要克服一定外力才能实现自由对流,K指数为41℃,SI指数<0℃。有利于短时强降水和对流大风天气的出现。300 hPa以下风随高度顺转,有暖平流,有利于大气不稳定层结加强,利于对流天气的发生发展。
分析个例3中蒙自探空站23日20时T-lnP图可知(图5(c)),大气湿层深厚,整层大气近似饱。抬升凝结高度低,约850 hPa,CAPE值为998.7 J/Kg,CIN为17.7 J/Kg,说明需要克服一定外力实现自由对流。K指数为43℃,SI指数<0℃,400 hPa以下风随高度顺转,有暖平流,垂直风切变较小。以上条件有利于短时强降水的出现。
Figure 5. Sounding data of Mengzi observe station. (a) 08:00 BST 29 July 2016; (b) 08:00 BST 19 June 2017; (c) 20:00 BST 23 June 2018
图5. 蒙自站探空资料。(a) 2016年7月29日08时;(b) 2017年6月19日08时;(c) 2018年6月23日20时
4. 暴雨成因分析
形成暴雨需要有充足的水汽、强的上升运动和较长的持续时间。水汽条件是暴雨发生的必要条件,如果没有水汽补充,某一局地上空的水汽总量是有限的,要使降水维持,必须有源源不断的水汽输送[8] [10]如图6(a)、图7(a)、图8(a)所示,个例1中,马关28日08时~29日20时水汽条件均较好,850 hPa比湿均超过12 g/kg,700 hPa比湿超过8 g/kg。其中28日21时~29日06时比湿增强,700 hPa比湿明显增加,均超过10 g/kg,28日23时700 hPa以下比湿达到最大值。垂直运动方面,28日18时~29日08时马关地区均为上升运动,且垂直上升运动发展高度高,其中29日03时上升运动最强,绝对值超过3 Pa/s。29日14时~15时再次有垂直上升运动发展。分析马关水汽通量散度随时间的变化可见,28日20时~29日08时约600 hPa以下有水汽通量辐合,最大超过3·10−7 kg/(m2·s·hPa)。29日08时~29日15时低层约800 hPa以下有弱的水汽通量辐合。分析马关地区该次暴雨期间的热力条件发现,28日08时~28日19时700 hPa以下大气为位势不稳定层结,28日20时~29日05时受台风残余低压影响,为稳定层结,29日06时~29日20时约600 hPa以下大气再次转为位势不稳定层结。对应降水发现,29日04时~05时和08时~09时为主要强降水时间段。04时~05时此时受台风残余低压影响,大气稳定层结,但03时起水汽和垂直运动条件增强,故降水增强。08时~09时低层比湿增强,有一定的垂直上升运动触发位势不稳定,故降水增强。29日08时~29日09时降水增强,此时的垂直运动和水汽通量散度并未最强,但由于此时大气为不稳定层结,仅需一定的垂直运动触发不稳定,形成对流降水。
Figure 6. Temporal evolution of specific humidity (contours, unit: g/kg) and vertical velocity (shaded, unit: Pa/s) at different heights of Maguan observe station. (a) Case 1; (b) Case 2; (c) Case 3
图6. 马关站不同高度比湿(等值线,单位:g/kg)、垂直运动(阴影,单位:Pa/s)随时间的演变。(a) 个例1;(b) 个例2;(c) 个例3
如图6(b)、图7(b)、图8(b)所示,个例2中,过程期间850 hPa比湿超过14 g/kg,700 hPa比湿超过10 g/kg,其中19日12时~20时和20日10时~17时850 hPa比湿超过16 g/kg。垂直运动方面,19日13时~15时、19日18时~20时和20日17时~18时有较明显的垂直运动,其中19日18时~20时最强上升运动超过0.9 Pa/s,且低层至高层均为上升运动;20日17时~18时垂直运动也较强,但700 hPa以下为弱的上升运动,不利于触发位势不稳定。分析马关水汽通量散度随时间的变化可见,19日18时~20时、19日21时~24时和20日18时低层有较明显的水汽通量辐合。热力条件方面,19日08时~22日20时马关600 hPa以下均为位势不稳定,有一定的强对流潜势,且午后低层增温增湿后大气位势不稳定性增强。700 hPa以下西南气流为主,有利于向暴雨区输送水汽。由于水汽通量散度和垂直上升运动均在19日18时前后最强,触发位势不稳定,故18时降水最强。
如图6(c)、图7(c)、图8(c)所示,个例3中,过程期间850 hPa比湿超过14 g/kg,700 hPa比湿超过10 g/kg。24日04时~12时、25日00时~07时马关受上升运动影响,且垂直运动层深厚,最大上升速度超过3 Pa/s,出现在25日05时。分析马关水汽通量散度随时间的变化可见,24日04时~14时和25日02时~06时马关地区均有水汽通量辐合,且24日02时和25日04时前后水汽通量辐合强,且高度高,辐合层最大高度超过700 hPa。热力条件方面,700~500 hPa有弱冷空气影响,假相当位温低,23日08时~25日08时马关地区均为位势不稳定层结或者中性层结,其中24日06时和25日05时稳定度略有增强。结合降水随时间的演变可见,24日08时~09时降水增强,对应低层水汽通量辐合增强和垂直运动增强,且水汽和垂直运动增强较降水增强提前2小时左右。
Figure 7. Temporal evolution of moisture flux (vectors, unit: kg/( m·s·hPa)) and moisture flux divergence (shaded, unit: 10−7 kg/(m2·s·hPa)) at different heights of Maguan observe station. (a) Case 1; (b) Case 2; (c) Case 3
图7. 马关站不同高度水汽通量(矢量,单位:kg/(m·s·hPa))和水汽通量散度(阴影,单位:10−7 kg/( m·s·hPa))随时间的演变。(a) 个例1;(b) 个例2;(c) 个例3
Figure 8. Temporal evolution of horizontal wind field (vectors, unit: m/s) and pseudo-equivalent potential temperature (shaded, unit: K) at different heights of Maguan observe station. (a) Case 1; (b) Case 2; (c) Case 3
图8. 马关站不同高度水平风场(矢量,单位:m/s)和假相当位温(阴影,单位:K)随时间的演变。(a) 个例1;(b) 个例2;(c) 个例3
个例1和个例3降水的水汽,上升运动条件较个例2好,但降水最大小时雨强弱于个例2。个例2中19日午后低层升温,中层有弱的冷空气影响,大气不稳定度增强。14时起马关基本受上升运动控制,18时触发位势不稳定,降水增强。综合以上分析可见,降水受水汽、垂直运动加强而加强[4],且水汽和垂直运动增强可较降水增强提前1~2小时。另外,位势不稳定触发也可增强降水。
5. 雷达特征分析
个例1中,0.5˚仰角反射率因子雷达回波显示(图9),28日21:54显示上游大片降水回波将影响马关地区,回波结构紧凑,最大值超过40 dBZ,说明有一定对流性。之后回波后向发展,向西北移动并伴有旋转,持续经过马关,形成列车效应,造成马关地区长时间降水。29日04:48强回波呈块状影响马关–麻栗坡一带,反射率因子最大超过40 dBZ,过雷达中心点和文山做垂直剖面可见,回波质心低。径向速度图显示,速度零线呈“S”形,说明有暖平流,有利于降水的发展和维持。29日8:29文山地区受两条带状回波影响,其中一条经过文山,回波向东北方向移动,伴有后向发展,造成降水持续时间长。过雷达中心点和文山做垂直剖面可见,降水回波质心较低,降水效率高。该时刻径向速度图上显示“S”形,说明有暖平流,暖平流有利于降水的维持。
个例2中(图10),午后文山西北部有分散性的对流回波出现,并向东南方向移动,13:00回波经过文山雷达并有组织加强向文山移动,由于马关对流条件较好,回波东移南下至马关加强。19日17:11回波显示,马关东部回波呈块状,中心强度超过50 dBZ。过强中心和雷达中心点做垂直剖面可见,马关回波发展高度较高,15 dBZ回波伸展至7~8 km高度,且45 dBZ区发展高度高。径向速度图上显示,马关东部有径向速度辐合,垂直剖面图上也可见速度辐合区,且辐合层厚度厚,超过5 km。后期上游仍有分散回波移入,马关降水强度减弱,但持续时间长。
个例3中(图11),24日08:20文山大部分地区受大范围层状降水回波影响,回波强度最大超过40 dBZ,降水回波东移或西移南压,且上游不断有回波新生,形成列车效应,持续影响马关,造成马关地区降水持续,达到暴雨。过马关附近强中心和雷达中心点做垂直剖面可见,回波影响范围大,有多个中心,质心高度低,最大伸展高度约5 km。径向速度图上显示,马关地区有风场的切变,有利于降水回波移入或加强。
Figure 9. Radar reflectivity (unit: dBZ) and radial velocity (unit: m/s) of Wenshan station at 0.5˚ elevation angle. 2016 July: (a) Reflectivity at 21:54 BST 28 July; (b) Reflectivity at 04:48 BST 29 July; (c) Vertical cross-section of reflectivity along the white line in panel (b) (the same below); (d) Radial velocity at 04:48 BST 29 July; (e) Reflectivity at 08:29 BST 29 July; (f) Vertical cross-section of reflectivity at 08:29 BST 29 July; (g) Radial velocity at 08:29 BST 29 July
图9. 文山站雷达反射率因子(单位:dBZ)和径向速度图(单位:m/s)。0.5˚仰角2016年7月(a) 28日21:54反射率因子;(b) 29日04:48反射率因子;(c) 29日04:48过图(b)中白色线条(下同)反射率因子垂直剖面;(d) 29日04:48径向速度;(e) 29日8:29反射率因子;(f) 29日8:29反射率因子垂直剖面;(g) 29日08:29径向速度
Figure 10. Same as Figure 9, but for 0.5˚ elevation angle at 17:11 BST 19 June 2017. (a) Reflectivity; (b) Vertical cross-section of reflectivity; (c) Radial velocity; (d) Vertical cross-section of radial velocity
图10. 同图9,但为0.5˚仰角2017年6月19日17:11。(a) 反射率因子;(b) 反射率因子垂直剖面;(c) 径向速度;(b) 径向速度垂直剖面
Figure 11. Same as Figure 9, but for 0.5˚ elevation angle at 08:20 BST 24 June 2018. (a) Reflectivity; (b) Vertical cross-section of reflectivity; (c) Radial velocity
图11. 同图9,但为0.5˚仰角2018年6月24日08:20。(a) 反射率因子;(b) 反射率因子垂直剖面;(c) 径向速度
综合以上分析可见,个例1和个例3中雷达特征显示降水以层状云降水为主,回波强度在25~35 dBZ之间,质心高度低,低于5 km,径向速度场有暖平流特征或风向切变特征,有利于降水的持续。个例2中降水具有一定的对流性,回波强度最大超过50 dBZ,回波质心高度较高,约达7~8 km,径向速度场显示有中尺度的辐合,垂直剖面可见辐合层厚度较厚,超过5 km,强的中尺度辐合提供了较好的动力条件,有利于触发和加强对流。
6. 结论
本文利用ERA5 0.25˚ × 0.25˚逐小时的再分析资料、Micaps国家站和自动站气象要素资料、文山雷达资料等,对马关地区三次不同背景下的暴雨过程天气背景、成因及雷达特征进行了分析,以提升马关地区暴雨预报准确率。结论如下:
分析形势场发现,个例1中马关500 hPa位于两高之间均压场,引导气流弱,残余台风低压移入后,在马关附近打转,后期减弱,向东北方向移去,影响时间长。个例2中500 hPa位于588 dagpm和584 dagpm之间,有短波槽影响,中高纬冷空气渗透,短波槽加深;700 hPa受低涡底部切变线影响。个例3中500 hPa位于两高之间,有高原涡影响,并伴有短波槽东移;700 hPa有低涡切变影响,由于副高增强,阻挡短波槽东移。以上形势均有利于马关地区长时间受上升气流影响,有利于长时间降水。
分析探空可见,三次个例中马关地区均有一定的对流能量,对流抑制小,抬升凝结高度低,K指数在40℃附近,SI指数<0℃。个例1和个例3湿层深厚,有利于短时强降水的出现,个例2中400 hPa有干冷空气侵入有利于出现短时强降水和大风。
分析暴雨成因,三次过程比湿条件均较好,个例1中,850 hPa超过12 g/kg,700 hPa超过8 g/kg。个例2和个例3中,850 hPa超过14 g/kg,700hPa超过10 g/kg。三次过程期间均伴有水汽通量散度辐合,个例1辐合强度最强,个例2和个例3辐合相对较弱。三次过程均有上升运动,个例1和个例3中上升运动强,且持续时间长,最大超过3 Pa/s;个例2中低层以弱的上升运动为主,部分时刻上升运动增强,超过0.8 Pa/s。降水受水汽、垂直运动加强而加强,且水汽和垂直运动增强可较降水增强提前1~2小时。另外,位势不稳定触发也可增强降水。
分析雷达演变特征可见,个例1和个例2以层状云降水为主,最大反射率因子超过40 dBZ。马关地区持续受回波影响,形成列车效应。从垂直剖面可见,回波质心低。径向速度特征上,个例1中有“S”形特征,说明有暖平流;个例3中有风场的切变。个例2中最强时段降水有对流性,19日17:11回波呈块状,中心强度超过50 dBZ。垂直剖面图上,马关回波发展高度较高,可伸展至7~8 km高度。径向速度图上有速度辐合,辐合层厚度超过5 km。
基金项目
文山州气象局2023年科研项目(WSZL202319)。
NOTES
*通讯作者。