李青春 程丛兰 全继萍 陈 敏 窦有俊 仲跻芹
北京城市气象研究院,北京 100089
提 要: 为了揭示北京地区夏季暴雨形成过程中的低空急流成因、结构特征及其作用,利用风廓线雷达测风数据、NCEP/GFS再分析资料(0.5°×0.5°)、国家级地面气象观测站和区域自动气象站小时降水量观测数据,对两次不同天气系统类型和强度等级的暴雨过程开展对比分析,得到结论如下:2018年7月16日(简称“7·16”)北京特大暴雨是在副热带高压边缘低空西南急流的影响下形成的,低涡低槽东移发展、低空急流不断加强,加之地形辐合线的影响,为大暴雨的发生提供了极为有利的低空水汽输送、辐合和抬升条件。2019年7月22日(简称“7·22”)北京暴雨是在高空槽前低空急流的影响下形成的,低空急流是由高空急流东移、下传所产生,低空水汽输送和辐合作用相对较小。在两次暴雨过程中的强降水发生前3 h,均出现低空急流强度增大、低空急流最低高度下降、低空急流指数激增、1500 m以下出现明显的垂直风切变并且逐渐增大等特征。在“7·16”暴雨过程中,边界层急流遇到山脉地形触发了对流新生,低空急流、地形辐合线共同作用使对流系统不断发展和组织化,是沿山地区极端强降水形成的关键原因。低空急流特征量(低空最大风速、急流最低高度和低空急流指数)和1500 m以下的垂直风切变大小对强降水的发生具有重要的指示意义。
自20世纪70、80年代以后,我国学者针对低空急流对华南前汛期暴雨、江淮气旋暖区暴雨形成的作用开展了大量研究,结果表明:低空急流轴左前侧存在较强辐合区、理查森数Ri负值区对应的低层最不稳定区,是产生暴雨的有利动力条件(孙淑清和翟国庆,1980;高守亭和孙淑清,1984;翟国庆等,1999)。陶祖钰(1980)分析得出,低空急流总是与暖舌相结合,使低层的假相当位温θse增大,导致500 hPa 以下对流层中下层条件不稳定和对流有效位能CAPE增大,对流抑制能量CIN减小,有利于对流的发展。同时,急流轴上的风速不连续性以及水汽输送作用是引起水汽聚集的主要原因。副热带高压(以下简称副高)型和强西南急流型的暴雨落区分别位于强低空急流出口区左前侧和副高内的水汽充沛且大气层结不稳定区内(汪玲瑶等,2018)。
以往针对低空急流、边界层急流的研究多数是利用常规气象探空资料、极少的高山气象站测风数据、NCEP全球再分析资料,或者外场试验期间获取的加密探空数据开展的。Mitchell et al(1995)利用风廓线雷达观测资料研究美国夏季大草原上的低空急流时发现,一天两次的常规气象探空不能完整地监测低空急流,常漏掉最强急流的出现时刻;Zhong et al (1996)利用美国风廓线雷达组网资料开展个例研究,揭示了夏季美国大平原的低空急流的演变过程。对布设在美国中部的风廓线雷达网做出的评估结果显示(The National Weather Service and the Office of Oceanic and Atmospheric Research,1994),风廓线雷达的时空分辨率测风数据超过任何一种高空风探测系统,可显示出短波天气系统的连续变化过程;Joshi et al(2006)利用风廓线雷达资料研究了印度普纳季风期的低空急流的高发时间,白天与夜间的低空急流高度的差异等。
近年来,随着中国气象部门风廓线雷达站网的建设,其高时间分辨率的测风数据成为研究低空急流特征的重要手段。Du et al(2012)利用风廓线雷达数据开展上海梅雨期前后的低空急流气候特征分析发现,存在两类低空急流,分别是边界层急流和与天气系统相关的低空急流。Du and Chen(2019)研究认为,华南沿海暖区暴雨的发生发展与边界层急流、天气尺度急流的相互作用和日变化特征有密切关系。周芯玉等(2015;2019)研究了广东前汛期的短时强降水(降水强度,简称雨强,≥20 mm·h-1)临近时次的低空急流特征,得出:在广东前汛期,有60%以上的短时强降水在发生前3 h出现低空急流,低空急流与低层垂直风切变的变化对于短时强降水类暴雨的发生具有重要的影响。一些学者研究了低空急流对强降水等天气的指示意义,刘淑媛等(2003)、张京英等(2003)和曹春燕等(2006)研究指出,风廓线雷达能够很好地揭示低空急流脉动的强度及向下扩展的程度与中小尺度强降水之间的关系,暴雨的产生主要是由低空急流的下传和加强引起。金巍等(2007)研究也表明强降水的强弱和低空急流的强度和向下伸展高度直接相关;2012年北京“7·21”极端暴雨发生后,在京津冀地区的暖区极端暴雨形成过程中低空急流的重要作用受到广泛关注(孙军等,2012;孙继松等,2012;姜晓曼等,2014;谌芸等,2018;雷蕾等,2020)。一些针对强暴雨过程中边界层急流的作用研究发现(朱乾根等,2001;魏绍远和林锡怀,1995),边界层急流是低层水汽的最大输送者和暴雨区对流不稳定能量释放的触发者。吴海英等(2010)、马旭林等(2004)的研究指出,边界层急流对暴雨天气过程中的中小尺度系统活动有重要的指示意义。
廖晓农等(2011)将北京暴雨日的环流型分为两大类:一类是热带、副热带系统和来自35°N以南的低涡及暖切变等系统。其中,副高边缘偏南暖湿气流中的对流云团带来的暴雨占比最高(占16%)。另一类是与北方冷空气活动有关的低涡低槽等系统。其中,蒙古低涡低槽的影响占比最高(占32%),副高边缘和西来槽的影响占比分别位列第二和第三位(占16%和14%)。另外,对近5年(2016—2020年)北京地区夏季暴雨日的主要影响系统的统计分类结果也证明了西来槽型、副高边缘型是北京夏季暴雨的最主要影响系统类型。在暴雨日中(共计42 d)这种两种类型的暴雨占比依次为26.2%、21.4%,其次为蒙古低涡低槽型占19.0%,其他类型占比较小。同时,在预报中发现,低空急流主要发生在西来槽前暖湿气流、副高边缘天气类型中,低空急流的成因和结构特征的不同,会导致降水量级和降水强度、落区范围存在明显的差异。因此,需要通过个例分析,揭示不同天气形势下的低空急流结构特征差异所引起的水汽输送、水汽辐合、抬升条件的差异,及其对降水量级和降水强度所造成的影响。本研究利用北京地区风廓线雷达30 min间隔的多层测风数据、NCEP/GFS再分析资料(0.5°×0.5°)、地面加密自动气象站小时降水量观测数据,并引入了低空急流特征量(低空最大风速、低空急流最低高度、低空急流指数)和低空垂直风切变,分别对副高边缘型和西来槽型暴雨过程(2018年7月15日20时至18日02时副高边缘型大暴雨过程,简称“7·16”过程;2019年7月22日15时至23日08时高空槽型暴雨过程,简称“7·22”过程;如无特殊说明,文中时间指北京时)的低空急流成因、低空急流特征量和垂直风切变演变特征及其对暴雨形成、降水强度的影响机制进行对比分析,进一步提高对北京地区极端暴雨形成机理的认识,归纳出引发两种天气类型暴雨的低空急流特征的预报要点,提高暴雨预报准确率。
北京地区位于华北平原北部,西部、西北部和东北部为山区,中部、东部和南部为平原(如图1),山区海拔高度为1000~1500 m。其中,北部属于燕山山脉,怀柔区、密云区和平谷区位于燕山南麓。房山区位于北京西南部的百花山东南坡和山前平原地区。
为了弥补常规气象探空数据(每天两次)、NCEP/GFS再分析数据(0.5°×0.5°,每天四次:02、08、14和20时)时间密度的不足,利用位于东北部降水量中心较近的怀柔风廓线雷达站测风分钟数据(图1),开展本站低空急流精细时间尺度的变化特征分析。
利用地面加密自动气象站小时降水量观测数据分析两次暴雨过程的降水量分布和时间演变,站点分布见图1。
图1 北京地形高度(填色)、行政区划(黑色虚线)和站点位置示意图(圆锥表示风廓线雷达站,红色圆圈表示自动气象站点;CP:昌平,YQ:延庆,HR:怀柔,MY:密云,PG:平谷,SY:顺义,TZ:通州,GXT:观象台,DX:大兴,FS:房山,MTG:门头沟)
对风廓线雷达0~5000 m各层6 min间隔的水平风观测数据进行质量控制和缺测插值。步骤为:(1)去除奇异值(≥40 m·s-1);(2)当某时次的水平风向风速出现缺测(或为奇异值)时,缺测值(或奇异值)由最近两个时次(-6 min、+6 min)的同层数据进行线性插值得出。若最近两个时次(-6 min、+6 min)有一个时次缺测(或为奇异值),缺测值(或奇异值)则由另一个时次的同层数据进行代替。若最近两个时次同时缺测(或为奇异值),则认为该时次数据缺测。
插值后获得2018年7月15—18日、2019年7月22—23日各层水平风向风速的整点和半点数据(30 min间隔)的时间序列。
在本研究中,将3000 m(700 hPa)及以下出现偏南风且风速≥12.0 m·s-1,风向为91.0°~269.0°视为出现了低空急流(LLJ)。将1000 m(900 hPa)以下出现偏南风且风速≥12.0 m·s-1,风向为91.0°~269.0°视为出现了边界层急流(BLJ)。
为了定量地反映低空急流强度和高度的变化,引入低空急流指数Ij,进一步分析暴雨形成与Ij的关系,计算公式为(刘淑媛等,2003;杨毅和邱崇践,2006):
(1)
式中:Wmax为3000 m高度以下低空急流中心的最大风速,Hmin为风速达到12.0 m·s-1及以上的最低高度。低空风速越大、急流所在高度越低,低空急流指数Ij越大(单位:10-3s-1)。
在分析低空急流引发的低空垂直风切变时,利用以下计算公式(周芯玉等,2019):
(2)
(3)
式中:M为上下层的垂直风切变,ΔZ为上下层的高度差值,V1、V2为上下层的风速,D为上下层的风向差值。由于受地物杂波影响,风廓线雷达低层数据质量较差。因此,计算150~870、150~1470、150~3030 m风的垂直切变以代表地面与925、850以及700 hPa之间的垂直风切变(单位:10-3s-1)。
低空急流的重要作用之一是对水汽和热量的输送作用。为了分析低空急流的水汽输送和辐合作用,利用GFS再分析数据的各层温度、水平风分量,以及相对湿度,计算水平水汽通量和水汽通量散度(董晓敏和田盛培,1986)。
“7·16”北京特大暴雨从2018年7月15日20时开始,18日02时以后降水逐渐结束,总历时超过52 h,整个过程达到大暴雨量级。在北京地区436个气象站中(63个站缺测),有207个站(约占55.5%)过程累计降水量超过100.0 mm,16个站超过200.0 mm, 3个站超过300.0 mm(其中,密云区西白莲峪最大达351.3 mm)。分析整个过程累计降水量分布(图2a),降雨落区呈东北—西南走向,与北部、西部和西南部地形走向一致,暴雨中心位于北京北部、东北部山区(密云、怀柔)和西南部山区(房山)沿山地区。由北京地区自动气象站的最大雨强(单位:mm·h-1)和短时强降水(雨强≥20 mm·h-1)出现站数百分比的时序图(图2c)可知,此次大暴雨过程主要由三轮强降水构成。其中,第一轮强降水(7月15日20时至16日14时,图2c中红色方框所示)为整个“7·16”过程的最强降雨时段。“7·16”第一轮强降水时长约为18 h,累计降水量达到大暴雨量级,有125个站(占总站数的33.5%)累计降水量超过50.0 mm,28个站超过100.0 mm,最大降水量中心出现在密云西白莲峪(301.9 mm)、捧河岩村(271.3 mm),均达到特大暴雨量级。最大降水量的贡献率为85.5%(该时段最大降水量/过程总累计降水量)。因此,本研究将第一轮强降水作为重点研究时段。
图2b为第一轮强降水时段累计降水量分布,暴雨中心位于北京东北部山区(密云、怀柔)和西南部山区(房山)沿山地区。在该时段先后出现两个强降水峰值时段,分别为16日01—05时和07—11时。期间的16日02—03时在密云区西白莲峪、捧河岩村出现最大雨强,分别达到117.0 mm·h-1和101.6 mm·h-1,超过2012年北京“7·21”特大暴雨过程的最大雨强(平谷挂甲峪为100.3 mm·h-1)。另外,16日02—03时出现短时强降水的站点数最多、影响范围最大,占32.4%。
图2 2018年7月(a)15日20时至18日08时,(b)15日20时至16日14时,北京地区累计降水量(紫色实线,单位:mm)、地形高度(填色)、风廓线雷达站点(空心三角)分布,以及(c)15日17时至18日08时最大雨强和短时强降水出现站数百分比时序
第二轮强降水从7月16日18时开始,最大雨强(46.2 mm·h-1,18—19时出现在昌平区沙河水库),20时以后雨势减小;17日07—11时出现第三轮较强降水,最大雨强进一步减小为44.0 mm·h-1(08—09时出现在密云区河南寨),17日12时以后雨势逐渐减小,18日02时以后降水逐渐结束。
“7·22”暴雨过程从2019年7月22日15时开始,23日08时结束,总时长近18 h,过程累计降水量达到暴雨量级。在445个自动站中,有71个站(占总站数的16.0%)累计降水量超过50.0 mm、4个站超过100.0 mm,最大降水量出现在密云石城(109.7 mm),过程出现暴雨的站数相对较少。降水量分布与“7·16”过程类似,降雨落区分布呈东北—西南走向(图3a),降水中心位于北京北部、东北部山区(密云、怀柔)和西南部山区(房山)的沿山地区。从北京地区自动气象站的最大雨强和短时强降水出现站数的百分比时序图可知(图3b),22日20时至23日04时出现短时强降水,最大雨强出现在22日22—23时,为40.3 mm·h-1(出现在密云区三峪培训基地),远小于“7·16”过程中在16日02—03时(117.0 mm·h-1)出现的雨强。23日00—01时短时强降水出现站数最多,仅占比为13.0%。
图3 2019年7月(a)22日14时至23日08时北京地区累计降水量(紫色实线,单位:mm)、地形高度(色标)、风廓线雷达站点(空心三角)分布,(b)22日08时至23日08时最大雨强和短时强降水出现站数百分比时序
分析高空环流形势得出,2018年“7·16”特大暴雨过程是在西太平洋强盛的副高西伸北进,北京处于副高边缘不断加强的西南气流控制,在500 hPa上游河套地区有一明显的高空槽东移(图4a),以及850 hPa北京附近有一条暖切变线的形势下发生的。地面图上(图4b),7月15日20时前后在上游河套地区北部有低压发展,我国华北地区处于东南高、西北低的气压场控制,风向由偏南转为东南风。导致在北京西北部山麓地区和太行山东部山麓地区沿地形形成一条辐合线,有利于山麓地区降水天气的形成;16日08时以后,由于850 hPa、700 hPa在上游河套东部地区有低涡发展(图略),导致副高边缘西南急流不断加强,地面东高、西低气压场长时间稳定维持,造成16日傍晚到前半夜、17日早晨到上午强降水的出现并长时间持续。
与“7.16暴雨”不同的是,在2019年“7·22”暴雨过程中副高位置偏南,副高588 dagpm线位于华北南部。暴雨是在500 hPa上游河套地区附近的高空槽东移、槽前伴有西南急流,以及850 hPa北京南部出现暖切变线的形势下发生的(图4c)。分析7月22日20时地面图(图4d),河北中南部地区(包括北京)处于低压辐合区北部,在北京东部—南部有一条偏东风与偏北风辐合线,对北京东部和南部地区降水天气的形成较为有利。
图4 (a,b)2018年7月15日20时,(c,d)2019年7月22日20时(a,c)500 hPa风场(黑色风羽)、位势高度场(黑色实线,单位:dagpm;棕色粗实线:高空槽)、850 hPa风场(红色风羽,红色实线:暖切变线),(b,d)地面海平面气压场(实线,单位:hPa)、风场(风羽,棕色粗实线:辐合线)
对两次暴雨过程中北京上空的垂直运动和辐合条件进行对比分析(图5),在“7·16”过程第一轮强降水发生之前,15日18时以后,随着地形辐合线和低空西南急流的先后出现,在900 hPa以下出现明显辐合(散度D<-12×10-5s-1)(图5a),900 hPa以上为较明显的辐散。与其相对应,在低层(800 hPa以下)为明显的上升气流(上升速度ω<-0.6 Pa·s-1),为第一轮极端强降水提供极为有利的低空辐合和抬升条件;受中低空河套低涡的东移发展影响,从16日17时开始出现整层上升气流,以700~300 hPa高度上升速度最大(上升速度ω<-1.2 Pa·s-1),900 hPa以下为明显辐合(散度D<-8×10-5s-1),一直维持到17日10时,为7月16日傍晚和17日早晨的第二轮、第三轮强降水的形成提供了有利的辐合条件。但近地层上升速度相对较小,抬升条件略差;17日17—24时,900 hPa以下再次出现辐合(散度D<-6×10-5s-1),但近地层上升速度较小且抬升条件较差。对应时间虽然出现了降水,但是降水强度远小于前三轮降水。
由图5b可以看出,在“7·22”过程中,随着高空槽、槽前西南急流的东移发展,7月22日17时至23日05时北京上空出现强烈的上升气流,850~700 hPa上升速度最大(上升速度ω<-1.8 Pa·s-1)。从22日15时开始近地层900 hPa以下出现较为明显的辐散,22日20时以后逐渐转为辐合。天气形势为强降水天气的形成提供了有利的抬升条件,近地层辐合条件由差逐渐转好。
图5 (a)2018年7月15日08时至18日08时和(b)2019年7月22日08时至23日20时两次暴雨过程中北京观象台上空垂直速度和散度的时间-高度分布(黑色等值线:垂直速度,单位:Pa·s-1;填色:散度,单位:10-5 s-1)
利用北京观象台探空数据分析两次暴雨发生之前北京的温度、湿度,以及对流不稳定条件。在2018年“7·16”副高边缘型大暴雨发生之前(7月15日14时、20时),对流层中低层比湿明显增大、湿层明显增厚,700 hPa增湿最为明显,比湿由3.3 g·kg-1增至10.3 g·kg-1。明显大于孙继松等(2012)对北京4次暴雨过程的700 hPa最大比湿为9.0 g·kg-1的统计结果;相比之下,2019年“7·22”降水发生前至短时强降水开始时刻(7月22日14时、20时)850 hPa 以下增湿明显。850 hPa比湿由7.4 g·kg-1增至15.2 g·kg-1,700 hPa比湿由3.0 g·kg-1增至8.4 g·kg-1,湿层厚度不及“7·16”暴雨之前。低层显著增湿与低空急流的水汽输送作用有关,两种天气类型的低空急流结构特征及其水汽输送作用、辐合作用存在差异。
3.1.1 “7·16”过程
分析“7·16”过程的850 hPa风场和全风速分布图发现:在北京第一轮强降水出现前,从7月15日20时开始河北南部地区出现较大风速区(风速≥8 m·s-1)(图6a)。15日20时至16日02时京津冀中东部地区出现西南低空急流,北京位于低空急流核的下游左前方(图6b)。因低空急流核的下游对应辐合区和不稳定区,有利于对流系统在该区域的发生和加强(Bonner,1968;孙淑清和翟国庆,1980;翟国庆等,1999;汪玲瑶等,2018),由此引发第一个强降水峰值出现。在北京所在40 °N做的水平风纬向剖面图上(图6e),相比于15日20时,16日02时低层风速加大,在900~800 hPa上空出现了低空急流。16日08时京津冀中东部地区的西南低空急流进一步增强,北京处于低空急流轴左侧(图6c),低空急流最低高度Hmin下降(图6f),对应时间出现第二个强降水峰值。由此可见,在2018年7月16日凌晨和上午的两个强降水峰值发生之前和发生时段,北京位于低空急流核区的下游左前方或者急流轴左侧,并且出现了低空偏南急流强度增大、低空急流高度下降等特征。
图6 2018年7月(a,d)15日20时,(b,e)16日02时,(c,f)16日08时(a~c)850 hPa风场分布和(d~f)沿40°N水平风的纬向剖面(填色为全风速,细实线为≥12 m·s-1全风速等值线;横坐标轴上黑色粗短线表示北京市域东西范围)
由“7·16”过程北京上空的水汽通量和水汽通量散度时间-高度(图7e)可以看出,在第一轮强降水发生前7~8 h,从7月15日17时开始,北京低空800 hPa以下出现水汽辐合并逐渐增强。近地层水汽辐合增强的时间与2.2 节中提到的东南风出现时间吻合。说明近地层水汽辐合的加强与东南风出现、地形辐合线的影响有关。在第一轮强降水发生前3 h (15日22时前后)低空急流开始出现。与之相对应,在900~800 hPa的水汽通量(填色所示)突然增大,16日凌晨至上午水汽通量达到18.0 g·s-1·cm-1·hPa-1以上。北京低空950 hPa以下水汽辐合进一步增强(散度D<-10×10-5s-1);从低空水汽通量和水汽通量散度分布可以看出,在16日凌晨、早晨至上午的两个强降水峰值时段,在850 hPa(图7c,7d),京津冀中东部地区(包括北京)位于低空急流轴左侧和水汽通量大值区左侧(最大水汽通量>18.0 g·s-1·cm-1·hPa-1)。在950 hPa,北京大部分地区及以东地区为水汽辐合区控制(图7a,7b),这些地区的最大水汽通量散度达-12×10-6g·s-1·cm-2·hPa-1,有利于对应区域的强降水形成。由以上分析得出,15日22时开始出现的低空急流,在900~800 hPa形成明显的水汽输送,并使地形辐合线影响的低空950 hPa以下水汽辐合进一步加强。两者共同作用是引发16日凌晨至上午北京出现第一轮强降水的关键原因。
图7 2018年7月(a,c)16日02时,(b,d)16日08时(a,b)950 hPa,(c,d)850 hPa水汽通量和水汽通量散度及风速(风羽)分布以及(e)15日08时至18日08时北京观象台的水汽通量、水汽通量散度时间-高度分布(填色为水汽通量,单位: g·s-1·cm-1·hPa-1;黑色实线为水汽通量散度,单位:10-6 g·s-1·cm-2·hPa-1;红色粗虚线表示全风速≥12 m·s-1的范围)
由图7e还可以看出,7月16日18—20时出现的第二轮强降水、17日早晨至上午(07—11时)出现的第三轮强降水均与北京近地层明显的水汽辐合有密切关系。主要原因是由于16日14—20时和17日02—08低空急流明显加强(图略),北京先后位于低空急流轴左侧、急流核的下游左前方,使近地层水汽辐合随着低空急流的变化而明显增强所造成;7月17日14时以后,尽管低空(850~700 hPa)仍存在明显的水汽输送(图7e),但是近地层水汽辐合明显较弱,入夜后出现降水的强度远小于前三轮降水强度。
总之,低空急流的水汽输送以及与地形辐合线、低空急流有关的近地层水汽辐合共同作用对强降水的发生起关键作用。近地层水汽辐合作用对雨强、累计降水量量级大小的影响更为重要。低空急流增强和维持,使低空水汽输送作用和近地层辐合作用加强和持续,使强降水反复出现,从而造成影响区域出现极端降水。
3.1.2 “7·22”过程
由“7·22”过程沿40°N的水平风纬向剖面可知,在暴雨过程发生之前,从7月22日08时开始高空急流下传,在北京西部上游地区中低空出现急流(图略)。在北京降水出现之前1 h (22日14时),850 hPa河北西南部出现偏南大风速区(图8a)。在对应时间的北京所在纬度40°N的水平风纬向剖面图上(图8d),低空急流出现在北京西部800~600 hPa上空。22日20时大风速区东移北上、增强(图8b),在北京南部出现低空偏南急流,北京位于低空急流核下游。分析对应的40°N水平风纬向剖面(图8e),急流中心位于850~550 hPa上空,急流强度加强,北京开始出现短时强降水。22日20时至23日02时低空急流继续东移加强、范围扩大,北京位于低空急流轴左侧(图8c)。急流中心和急流最低高度Hmin明显下降至900 hPa(图8f),22日22—23时北京出现最大雨强降水。以上分析说明,影响此次暴雨过程的低空急流是由东移的高空槽前西南气流加强所产生。随着强降水的临近,急流强度加强、急流高度明显下降,在时间上与强降水、最大雨强的出现时间有对应关系。
图8 同图6,但为2019年“7·22”过程
进一步分析2019年7月22日傍晚至23日凌晨的低空水汽通量和水汽通量散度分布(图9c),可以看出:22日20时位于850 hPa高空槽前的偏南急流控制了京津冀大部分地区,水汽通量大值区位于北京上游地区(河北西北部)和河北南部地区,并且沿急流轴有多个水汽辐合区。尽管950 hPa水汽通量较小,但是北京北部、南部仍处于水汽辐合区(图9a),特别是渤海湾附近受东南风的影响,河北中东部地区为水汽通量大值区和水汽辐合区。由北京上空的水汽通量和水汽通量散度时间-高度图可以看出(图9e),在7月22日暴雨发生前(22日08—12时),水汽通量大值区所在高度随着急流高度呈下降趋势。在22日傍晚前后,北京上空650~450 hPa的水汽通量突然增大(>18.0 g·s-1·cm-1·hPa-1)。另外,与渤海湾地区低空东南风相对应,在850 hPa附近出现了水汽通量大值区(>6.0 g·s-1·cm-1·hPa-1)。这支东南风的作用是在低空打通了一条由海上伸向内陆的水汽输送通道(侯淑梅等,2018)。与此同时,近地层出现水汽辐合明显(<-4×10-6g·s-1·cm-2·hPa-1),20时开始出现短时强降水。从22日20时至23日02时,低空急流移至北京下游地区(图9c,9d),北京位于低空急流轴左侧以及水汽通量大值区,并且850 hPa北京中部地区存在较明显的水汽辐合区。期间,22—23时出现了最大雨强降水;23日04时以后,中低层700 hPa以上的水汽通量快速减小(图9e),降水结束。由于在此次暴雨过程中,急流对北京上空的主要水汽输送集中在对流层中层650~450 hPa,对低层的水汽输送作用、近地层水汽辐合作用较弱,共同作用的时段较短,因此形成的降水量和降水强度远小于“7·16”过程的第一轮强降水时段。
图9 同图7,但为2019年“7·22”过程
3.2.1 低空急流特征量变化与作用
利用风廓线雷达测风数据可以得到常规气象观测无法反映的本地局地风场精细变化。分析2018年“7·16”过程第一轮强降水期间各层水平风的时间变化可以看出(图10a,10b):在第一轮强降水的第一个强降水峰值发生前3 h(7月15日22时前后),2500~3500 m高度出现低空急流(低空最大风速Wmax为12.6 m·s-1)。23时急流最低高度Hmin下降至1470 m(Wmax加大到15.1 m·s-1)。16日00时Hmin迅速下降到390 m(Wmax为16.5 m·s-1),随后出现强降水峰值。在强降水峰值期间,边界层风速继续增大,16日04:30的Wmax达到21.9 m·s-1,一直持续到16日05时前后;在第二个强降水峰值发生前2 h(7月16日06时),Wmax再次增大(由05时的14.4 m·s-1增至19.4 m·s-1)。Hmin再次下降,由06时的2310 m降至07时的870 m(Wmax达28.0 m·s-1)。08时的Hmin下降至750 m(Wmax为15.0 m·s-1)。风廓线雷达数据反映出第一个、第二个强降水峰值的出现与西南急流强度加强、所在高度不断下降,以及边界层急流的出现有关。
图10 2018 年7月15—16日怀柔站风廓线雷达观测的(a)各层水平风(填色为风速,黑色实线为≥12 m·s-1的等风速线)时间-高度分布,(b)Wmax和Hmin,以及(c)低空急流指数Ij的时序(横坐标轴上黑色短粗线表示降水量峰值时段)
低空急流的作用除了在900~800 hPa带来的明显水汽输送外,另一个重要作用是遇到山体后产生强烈的动力抬升作用。徐珺等(2014)对我国华北暖区暴雨的研究指出,超低空急流和低层垂直风切变在触发对流和对流维持中有重要作用。2018年7月15日23时Hmin由2550 m降至1470 m,在北京的西部和北部山区出现一些弱对流云团。7月16日00时的Hmin迅速下降至390 m(Wmax为16.5 m·s-1),即出现边界层急流。由于边界层急流所在高度明显低于北京西南部山区和东北山区的地形高度(山地海拔高度为1000~1500 m),由此引发山体对边界层急流的强烈抬升作用,这种抬升作用加剧了天气系统造成的垂直上升运动,进而触发了西南部山区、北部山区对流单体的形成(图11a),并在地形辐合线西侧的山前地区形成并强烈发展成中尺度对流雨带(图11b),造成了第一个强降水峰值的出现;同样,7月16日07时当边界层急流再次出现(Hmin为870 m,Wmax达28.0 m·s-1),首先在西南部山区(房山区)触发了对流单体形成(图11c)。随后对流单体沿着山前辐合线附近发展成中尺度对流雨带(图11d),造成第二个降水量峰值的出现。由此可见,边界层急流对对流系统具有重要的触发和组织作用。
图11 2018年7月16日(a)00时,(b)02时,(c)06时,(d)08时的地面风场(风矢)和雷达反射率因子(填色)分布(棕色双虚线表示辐合线;黑色粗实线表示300 m地形高度线)
低空急流指数可以直观地反映出低空急流强度和高度的变化。曹春燕等(2006)、金巍等(2007)、周芯玉等(2015)研究指出,强降水的出现与低空急流指数的增大有密切关系。计算分析2018 年“7·16”过程的低空急流指数Ij(图10c),在第一个强降水峰值发生前3 h、2 h(7月15日22时、23时)的Ij分别为4.9×10-3s-1、10.3×10-3s-1,发生前1 h(16日00时)激增至42.3×10-3s-1。同样,在第二个强降水峰值发生前2 h(16日06时)Ij为8.4×10-3s-1,发生前1 h(07时)激增至32.2×10-3s-1。由此可见,在两个降水峰值发生之前,Ij具有明显增大特征。
同样,由2019 年“7·22”过程的风廓线雷达观测的各层水平风时间-高度分析可知(图12a,12b),7月22日小雨出现前4 h(22日11时),在北京上空2500~3000 m出现低空西南急流,14时Wmax增至16.9 m·s-1,Hmin下降至1110 m,15时Hmin突然降至630 m以下,Wmax达20.0 m·s-1,15时出现小雨。7月22日17—19时的Hmin基本维持在2000 m高度上下,Wmax分别为17.6 m·s-1(17时)、16.0 m·s-1(18时),19时增至22.9 m·s-1,20时短时强降水开始出现。21时开始Hmin再次下降,22时下降至1470 m,23时下降至1230 m,Wmax由12.3 m·s-1增至17.4 m·s-1。22—23时出现最大雨强降水。由此可见,低空急流强度增大、低空急流最低高度下降,与强降水、最大雨强降水出现时间有对应关系。另外, 22日下午Hmin下降至1000 m以下,但低空水汽输送弱、近地层为辐散(图9e),仅出现了小雨天气。在22日21时以后低空急流高度下降高度在1000 m以上,引起的地形抬升作用并不显著(图略)。
分析“7·22”过程的Ij时间变化(图12c),在7月22日下午小雨出现前1 h(14—15时)Ij激增,由3.8×10-3s-1增至51.3×10-3s-1。在22日20时开始的强降水之前,Ij增大并不明显。但在强降水期间,Ij增大,20时、21时分别为7.2×10-3s-1、7.7×10-3s-1,22时增至12.3×10-3s-1,22—23时出现最大雨强降水。可见,Ij对最大雨强出现时间具有指示性。
图12 同图10,但为2019年“7·22”过程
3.2.2 低空垂直风切变演变特征与作用的比较
在低空急流覆盖范围内,除了等压面上存在着急流核区并伴有强烈的水平风切变外,垂直方向的风速廓线也存在着极大风速值和显著的垂直风切变。垂直风切变即指水平风(包括大小和方向)随高度的变化,其大小往往和形成风暴的强弱有密切关系(俞小鼎等,2006)。边界层内较大的风速垂直切变使低层大气处于极不稳定状态,有利于本地对流活动发展和组织化(《华北暴雨》编写组,1992)。
图13a为“7·16”过程第一轮强降水期间各层垂直风切变的时间-高度,可以看出:在第一个强降水峰值出现前6~7 h(从7月15日18时开始),低空1200 m以下开始出现垂直风切变。随着强降水的临近,低空风垂直切变明显增大、大值区明显增厚。16日06时低空1500 m以下再次出现垂直风切变增大、大值区增厚现象,08时开始出现第二个强降水峰值;同样,在“7·22”过程中(图13b),从7月22日上午小雨出现前6 h(从08时开始)900 m以下出现垂直风切变,小雨出现前1 h(14—15时)再次出现明显的垂直风切变。强降水出现以后(22日20时开始),1500 m以下再次出现垂直风切变,其大小和大值区所达高度远不及“7·16”过程,22日22时至23日00时出现的最大雨强也相对较小。
图13 (a)2018 年7月15日14时至16日14时,(b)2019年7月22日08时至23日08时两次暴雨过程的低层垂直风切变(填色)时间-高度变化(细实线表示垂直风切变≥120×10-3 s-1的等值线;横坐标轴上黑色短粗线段表示降水量峰值时段)
经与本站低空急流特征量的时间变化比较发现,垂直风切变增大与急流强度增大、急流高度下降在时间上具有对应关系。这说明低空急流高度、强度的变化,导致低层风场垂直切变增大。分析涡度方程的涡度倾侧项(朱乾根等,1991):
(4)
低空急流的出现导致低空垂直风切变,即出现∂v/∂p<0,或者∂u/∂p<0。根据式(4),垂直风切变在∂ω/∂x>0急流左侧区域,或者∂ω/∂y<0急流前部区域造成低层垂直涡度增大∂ζ/∂t>0。由此得出:低空急流强度增大导致低空垂直风切变增大,使急流左侧区域、急流前部区域的低层垂直涡度∂ζ/∂t>0进一步增大,从而加剧上述区域的辐合作用。另外,垂直风切变增大能够使上升气流倾斜,使得在上升气流中形成的降水质点脱离上升气流,而不会因拖曳作用减弱上升气流的浮力(俞小鼎等,2006)。
对两次过程强降水出现前3 h的低空三层垂直风切变(地面与925 hPa 、地面与850 hPa、地面与700 hPa)进行分析比较发现(表1):三层垂直风切变随高度依次减小。地面与 925 hPa的垂直风切变最大,地面与850 hPa其次,地面与700 hPa最小;在“7·16”过程第一轮强降水的两个强降水峰值发生前,低空风垂直切变明显增大,尤其是以边界层(地面与925 hPa)增大最为明显。在第一个强降水峰值出现前3 h和发生时刻(15日22时至16日01时)地面与925 hPa的垂直风切变依次为6.5×10-3、19.8×10-3、24.4×10-3和20.3×10-3s-1),地面与850 hPa的依次为8.8×10-3、11.9×10-3、13.4×10-3和8.5×10-3s-1,在强降水发生前1 h达最大。在第二个强降水峰值出现前3 h和发生时刻(7月16日05—08时),地面与925 hPa、地面与850 hPa的垂直风切变随时间增大,分别为7.5×10-3、2.8×10-3、6.0×10-3和18.6×10-3,以及4.2×10-3、4.5×10-3、5.5×10-3和6.6×10-3s-1),在强降水发生时刻达到最大;在“7·22”过程的强降水发生前3 h和发生时刻(7月22日17—20时),地面与850 hP的垂直风切变随时间增大并不明显,依次为6.6×10-3、7.1×10-3、7.2×10-3和7.9×10-3s-1,在强降水发生时刻达最大。
表1 强降水发生前3 h低层垂直风切变统计表(单位:10-3 s-1)
由此得出:两次暴雨过程发生前低空均出现垂直风切变。“7·16”过程第一轮强降水期间,低空1500 m以下垂直风切变的大小、大值区所达高度远大于“7·22”过程的强降水时段。随着强降水的临近,低空1500 m以下(地面与925 hPa、地面与850 hPa)的垂直风切变增大,在发生时刻达到最大,1500 m以下垂直风切变的大小与降水强度有密切关系,对强降水预报有明显的指示意义。
本研究利用北京地区风廓线雷达30 min间隔的多层测风数据、美国NCEP/GFS再分析资料(0.5°×0.5°)、地面加密自动气象站小时降水量观测数据,并引入了低空急流特征量(低空急流最大风速Wmax、急流最低高度Hmin、低空急流指数Ij)和垂直风切变,对两次不同天气类型和不同强度的典型暴雨过程中的低空急流成因、结构特征,以及边界层急流及其作用差异开展对比分析,得出以下主要结论:
(1)“7·16”北京特大暴雨是在北京处于副高边缘的西南急流控制、受高空低涡低槽东移和850 hPa暖切变线影响,以及地面东高西低的形势下发生的。副高边缘的低空西南急流加强、低涡低槽东移,以及近地层辐合线的共同作用,为第一轮极端强降水的发生提供了极为有利的水汽输送、辐合和抬升条件;副高、地面东高西低气压场的稳定维持,以及低空西南急流不断加强和维持,是造成强降水重复出现、降水天气持续的主要原因;“7·22”北京暴雨是在高空槽前偏南急流东移,850 hPa河套东部暖切变线,以及北京地面处于倒槽北部辐合区的形势下发生的。高空槽前偏南急流附近的较强的、深厚上升运动,为强降水天气的形成提供了有利的抬升条件和中层水汽输送条件,但是低空水汽输送条件和近地层的辐合条件较差。
(2)影响“7·16”过程的副高边缘西南低空急流在900~800 hPa具有明显的水汽输送作用。同时,由于地形辐合线的持续作用以及始终位于低空急流核区的下游左前方或者低空急流轴左侧,使北京低层(950 hPa以下)的水汽辐合不断加强。低空急流的明显的水汽输送作用及其与近地层不断加强的水汽辐合作用相互配合,是引发7月16日凌晨至上午出现极端强降水的重要原因;相比之下,影响“7·22”过程的西来槽前偏南急流的水汽输送主要在对流层中层650~450 hPa,低空(850 hPa以下)的水汽输送作用较小。另外,西来槽前偏南急流随高空槽在自西向东移动过程中,前期北京处于急流轴右侧、后期处于急流前方或急流轴左侧,到后期低空水汽辐合作用才逐渐转好。低空水汽输送和水汽辐合作用相互配合的时间较短。
(3)对利用风廓线雷达测风数据计算低空急流特征量的分析得出:两次暴雨过程中在强降水发生前3 h,均出现Wmax增大、Hmin下降,Ij增大,以及低空1500 m以下均出现明显的垂直风切变,且随着短时强降水临近而数值增大等特征。
(4)在“7·16”过程第一轮强降水的两个降水量峰值发生之前,急流最低高度下降至几百米,即出现边界层急流。边界层急流遇地形产生的强烈抬升作用,触发对流单体的形成。对流单体沿地形排列,在地形辐合线附近或者沿山一侧强烈发展成中尺度对流雨带,造成沿山地区极端强降水的发生。在“7·22”过程强降水发生前,低空急流最低高度在 1000 m以上,地形抬升作用不明显。
综上所述,在北京暴雨形成过程中,副高边缘低空急流比西来槽前低空急流的水汽输送作用、低空辐合和抬升作用更为显著。边界层急流对中尺度系统的触发和组织作用是极端强降水形成的重要原因。低空急流特征量(Wmax、Hmin、Ij),以及1500 m以下的垂直风切变的大小对强降水的发生具有重要的指示意义。
致谢:本研究得到北京城市气象研究院肖现研究员提供的雷达数据支持,在此表示衷心感谢!