2 中国科学院大气物理研究所云降水物理与强风暴实验室, 北京 100029;
3 南京信息工程大学气象灾害预报预警与评估协同创新中心, 南京 210044
2 Key Laboratory of Cloud-Precipitation Physics and Severe Storms, Institute of Atmospheric Physics, Chinese Academy of Sciences, Beijing 100029;
3 Collaborative Innovation Center on Forecast and Evaluation of Meteorological Disasters, Nanjing University of Information Science and Technology, Nanjing 210044
风切变是影响对流发展和组织的重要动力因素,很多学者通过观测和数值模式研究了风切变对对流系统的影响机制。Lilly(1979)认为超级单体排列成线状且与风切变成一夹角,从而使得各个单体之间的环流不相互干预,他们能够组成一个稳定的飑线,这一观测发现也得到了数值模拟结果的证实(Moncrieff,1978)。Rotunno et al.(1988)利用数值模拟进一步证实了冷池和低层风切变在非超级单体飑线维持中的作用,提出了强飑线维持的RKW理论,当冷池产生的负涡度与环境风切变产生的正涡度达到平衡时是强飑线维持的“最佳状态”,并提出了长生命史飑线的两种类型:一种发生在低层强切变、高层弱切变下,另一种发生在强的深厚切变下,这种环流特征使得单体之间的三维环流场互不影响。从20世纪80年代至今,风切变的研究越来越受到学者的关注,很多研究采用观测资料揭示了对流组织化与环境风垂直切变的密切关系(Bluestein and Jain,1985; Parker and Johnson,2000; Schumacher and Johnson,2005; Johnson et al.,2005),而且数值模拟研究对此给出了定量表述(Rotunno et al.,1988; Lafore and Moncrieff,1989; Weisman and Rotunno,2004),这些研究表明低层切变在对流发展和组织中有重要作用。Johnson et al.(2005)揭示了在热带与副热带地区,低层和中层的风切变对对流的组织作用是相同的,他们采用观测证实了低层切变垂直于对流 线时,由于其阻止冷出流快速远离风暴从而使风暴长时间维持。如果切变集中在高层(低层无切变 或弱切变的情况下),与系统在无切变环境下的相似,即对流发展和组织性弱(Weisman and Rotunno,2004)。虽然冷池与低层切变的相互作用(RKW理论)对飑线的维持有重要作用,Lafore and Moncrieff(1989)认为对于有层状云的飑线,层状云与对流之间的相互作用对飑线的影响作用也不可忽视。
孙建华等(2014)针对水汽整层含量和垂直 分布的数值试验表明:整层水汽增加(减少),对流增强(减弱),冷池和雷暴高压增强(减弱)导致大风增强(减弱);在保持整层水汽含量不变的情况下,线状对流和雷暴大风易发生在中层干、下层湿的环境中,这种层结条件对雷暴高压的增强有重要作用,但不利于整个对流系统的长时间维持。水汽试验主要反映的是热力作用对飑线的影响,而过去很多研究发现中层特别是低层风切变对中尺度对流系统组织的作用,高层的风切变对对流影响较弱。Zheng et al.(2013)统计了中国江淮地区的线状对流,由于观测资料的时间分辨率问题,风切变对不同组织类型的线状对流的影响差异不大,因此,需要通过数值试验来研究东亚季风区风切变对对流系统的影响。
过去的研究强调了低层风切变在对流发展和组织中的重要作用,特别是低层切变与冷池达到平衡时对于飑线长时间维持的作用。很少有研究涉及到中层切变和整层切变对于对流强度和组织的影响,另外,前人的理想数值模拟研究多采用北美的探空特征,而孙建华等(2014)的工作又只针对了水汽特征开展了试验。本文的理想数值模拟将采用我国东部江淮地区的观测探空,研究不同层次和不同强度的风切变对对流的强度和组织结构的影响。主要研究低层(0~3 km)、中层(3~6 km)和整层(地面到18 km)风切变对中尺度对流系统强度和组织类型的影响,哪种切变分布最有利线状中尺度对流系统的形成及其对对流强度的可能影响。
2 试验方案采用WRF模式进行中层、低层和整层的风切变对对流的强度和组织作用的理想试验,理想试验是在均一场中加温度扰动场,均一场的温度、气压、湿度和风场是由观测探空决定的。模式区域为360 km×360 km×18 km,水平分辨率为2 km,垂直方向为41层,侧边界条件采用开放边界。微物理方案采用Lin方案,不考虑辐射、边界层过程。背景场的探空是山东德州陵县2010年8月8日14时(北京时间,下同)修正的探空资料(图 1),即采用14时的地面观测修正08 时的探空,具体的修正方法见(Johnson and Bresch,1991; 潘玉洁等,2008)。该探空从最不稳定层抬升的对流有效位能是4471 J kg-1,整层可降水量达到68.52 mm,中低层水汽充沛(图 1),为郑淋淋和孙建华(2013)定义的湿环境,即整层可降水量大于50 mm。与该探空对应,8月8日20时山东西北部出现了降水回波,主要以短时强降水为主,其中观测到9日02~03时最强降水达到了102 mm。模式启动扰动是WRF理想试验模块默认的扰动,即初始扰动采用水平方向半径10 km,垂直方向半径1.5 km的温度扰动结构,简称扰动泡,扰动泡的中心位于地面1.5 km高度上,扰动泡位于模式模拟区域水平方向的中心位置。扰动泡的温度扰动分布为
$\Delta T={{T}_{0}}\times {{\cos }^{2}}(\beta \text{ }\!\!\pi\!\!\text{ }/2)(\beta \le \text{1})$ | (1) |
![]() |
图 1 山东德州陵县(37.32°N, 116.57°E)2010年8月8日14时(北京时间,下同)修正的探空图(其中实线为温度曲线,虚线为露点曲线)Fig.1 A skew-t diagram of the modified sounding data at Ling County (37.32°N, 116.57°E), Dezhou City, Shandong Province at 1400 BT (Beijing time) 8 August 2010 (the solid line represents temperature profile and the dashed line represent dew point profile) |
由于我国处于西风带,一般情况下观测风是纬向风的分量较大,这个探空下发生的真实个例和本研究中的模拟结果都是南北向分布的对流系统,过去的研究也认为垂直系统的风切变影响较大(Parker and Johnson,2000),此外,只改变纬向风使得试验较为简单,但也能反映风切变的影响,因此,本研究的试验中只改变纬向风的垂直切变。纬向风曲线的计算公式如下(Weisman and Klemp,1982):
$u={{u}_{0}}+{{u}_{\text{s}}}\tanh z/{{z}_{\text{s}}}$ | (2) |
![]() |
表 1 切变试验方案及其各层切变值 Table 1 The configuration and shear parameters in different levels of the experiments |
![]() |
图 2 CTRL、US-5、US5、L_US-5、L_US5、M_US-5和M_US5试验的纬向风廓线(单位: m s-1):(a)CTRL、US-5和US5试验;(b)CTRL、L_US-5和L_US5试验;(c)CTRL、M_US-5和M_US5试验Fig.2 The zonal wind profiles(units: m s-1)of different experiments:(a)CTRL,US-5,and US5;(b)CTRL,L_US-5,and L_US5;(c)CTRL,M_US-5,and M_US5 |
从回波的演变情况看,CTRL试验的触发、发展、成熟和衰减阶段分别为积分0~2 h、2~4 h、4~6 h和6 h以后。图 3为CTRL试验触发、发展和衰减阶段代表时刻的雷达组合反射率。触发阶段,回波为单体形式,最强回波强度达到50 dBZ以上。发展阶段,回波范围扩大,开始有西北—东南向的相对窄回波区出现,强回波带开始呈线状排列。成熟阶段,回波的线状特征明显,对流线向东南方向延长。从发展到成熟阶段,不断有新单体在对流线的东南方向形成,造成对流线延长。衰减阶段,对流尺度继续增大,强回波区减少,线状组织的对流区开始分散。从回波的整个发展过程看(图 3),回波有从中间向四周发展的趋势,在整个对流区的东部边缘形成了线状回波。下面的分析中分别采用积分4 h和7 h的模拟结果研究系统成熟和衰减阶段的特征。
![]() |
图 3 CTRL试验积分1 h、3 h、4 h、5 h和7 h的雷达组合反射率(单位:dBZ)。图上的坐标是格点数,格距为2 km(下同)Fig.3 The composited radar reflectivity(units: dBZ)after integrating 1 h,3 h,4 h,5 h and 7 h in the CTRL experiment. The coordinates are grid numbers and the grid interval is 2 km(the same below) |
对比分析几个垂直风切变试验,积分到2 h,7个试验的回波形状和强度差异不是很大,只是US-5试验的回波范围小一些(图略)。积分4 h,即系统的发展阶段,回波渐渐组织成线状,US-5和M_US-5试验强回波区没有组织成明显的线状结构,但L_US-5试验的结果仍然出现了线状对流,只是强度比L_US5弱一些(图 4)。积分7 h,即系统的衰减阶段,US-5和M_US-5试验的对流系统前部的强对流强度比其他试验的要弱,组织结构也更松散,但M_US-5试验对流的组织结构比US-5稍强(图 5),而L_US-5比L_US5的线状结构更清楚,只是对流的范围减小。从移动速度看,增加切变的试验均比减少切变的试验系统的移动速度快(图 4和图 5)。这些试验的结果表明,改变整层风切变对对流系统的强度和组织结构影响最大,其次是改变中层的切变。减小整层垂直风切变对流强度减弱,不利于对流组织成线状,增加整层垂直风切变,对流强度强,有利于对流组织成线状。垂直风切变减小,回波移动慢,垂直风切变增加,回波移动快。过去研究认为飑线的强度和组织类型主要依赖于低层切变的强度(Fovell and Ogura,1989),低层切变和冷池的相互作用影响飑线维持时间和强度(Weisman et al.,1988),除了低层切变,中层切变对飑线的结构和演变也有重要影响(Alexander and Young,1992; Ferrier et al.,1996; LeMone et al.,1998; Montmerle et al.,2000; Parker and Johnson,2000; Xue,2000)。本研究的试验结果表明整层风切变对对流强度和组织类型影响最大,与前人的结论有差异,有必要深入研究整层风切变如何影响中尺度对流系统强度和组织类型。由于US5与US-5的试验结果差别较大,下面主要针对这两组试验的来分析垂直风切变对对流系统的影响。
![]() |
图 4 US-5、US5、L_US-5、L_US5、M_US-5和M_US5试验在积分4 h的组合雷达回波反射率(单位:dBZ)Fig.4 The composited radar reflectivity (units: dBZ) after integrating 4 h in experiments US-5, US5, L_US-5, L_US5, M_US-5, and M_US5 |
![]() |
图 5 同图 4,但为积分7 h的试验结果Fig.5 As in Fig. 4,but for the simulation results after integrating 7 h |
积分1 h后,US5试验的近地面最大风速就明显比CTRL试验大,US-5试验最大近地面风速比CTRL试验小(图 6a),尤其在积分3 h后,US-5试验的近地面最大风速比CTRL试验小3~5 m s-1。但是US5试验近地面最大风速只在2~4 h明显比CTRL试验大,之后,与CTRL试验基本相似。虽然增加风切变,近地面大风增强,减小垂直风切变,不利于近地面大风的增强,但近地面大风增强主要在发展阶段。孙建华等(2014)的研究中发现增加水汽越多发展阶段冷池强度越强,最大风速越强,但成熟阶段后期冷池减弱的越快,层状云区的后部入流减弱,不利于雷暴大风的出现和维持。US-5试验最大降水量比CTRL试验小(图 6b),US5试验每小时最大降水量比CTRL试验大,但明显增强出现在4 h之后,即在成熟阶段,最强超过了50 mm h-1(图 6b)。表明增加垂直风切变,对流强度增强,降水量增大,减小垂直风切变,对流强度减弱,降水量减小,但降水主要是在成熟阶段增加。积分1 h后,US5试验最大上升速度比CTRL试验大,US-5试验最大上升速度比CTRL试验小(图 6c)。最大下沉速度也是积分1 h后开始出现明显的差异,US-5试验最大下沉速度比CTRL试验小(图 6d),US5试验比CTRL试验大,在发展和成熟阶段出现了5次大于20 m s-1的下沉风速,而CTRL试验只出现了一次。这与Schlesinger(1978)在数值模拟中得出的结论一致:切变环境中的上升运动与下沉运动比非切变环境中的强。垂直风切变的大小对近地面风速、降水量和垂直速度的影响结果表明,增加垂直风切变,上升速度增强,降水量增强,下沉速度增强,近地面风增强。试验结果表明垂直风切变的增加有助于对流组织成线状和增强下沉运动、近地面风速。此外,地面最强降雨出现在近地面最强大风之后,其原因何在,也需要进一步分析。
![]() |
图 6 US-5和US5试验(a)1 km高度上最大水平风速(单位:m s-1)、(b)每小时最大降水量(单位:mm)、(c)最大上升速度(单位:m s-1)和(d)下沉速度(单位:m s-1)的时间演变Fig.6 Temporal variation of the(a)maximum horizontal wind velocity(units: m s-1)at the 1-km height,(b)hourly precipitation amount(units: mm),(c)maximum ascending velocity(units: m s-1),and (d)maximum descending velocity(units: m s-1)in experiments US-5 and US5 |
在第3部分中对比分析了US-5与US5试验的组织结构、降雨量以及近地面风速等结果,分析表明不同的垂直风切变强度对对流系统的组织特征和强度等影响较大。积分3 h和5 h分别对应各试验的发展阶段和成熟阶段,且US5试验在积分3 h和5 h分别出现了近地面大风和降水量的大值,选定这两个时次进行风切变对流强度和组织类型的影响机理分析可以深入了解大风和强降水出现的原因。
4.1 垂直速度与水平风场1 km高度上,积分3 h 时US-5试验下沉区位于距离上升区约20 km左侧,大风区主要位于下沉区的北部,积分5 h上升区与下沉区的距离约为45km,上升运动区在下沉运动区的北部和东部,大风区位于下沉区的东部。积分3 h和5 h降水量大值区均位于下沉运动区的大值区(图 7),即为层状云降水区。US5试验上升运动区与下沉运动区基本是成对出现的,从积分3 h到积分5 h,沿着系统的东南方向不断有新单体出现,积分5 h形成线状组织类型。US5和US-5试验北部的扰动区形成的回波对应的上升运动、下沉运动、水平风速和降水量均很强(图 7和图 8),因此,下面的剖面图选择这一区域进行研究。
![]() |
图 7 US-5与US5试验的结果,上图为1 km高度上的垂直速度(等值线:单位:m s-1,等值线标值从小到大依次为-2、-1.5、-1、1、3、5 m s-1)、大于和等于17 m s-1的水平风速(风向杆,全风向杆代表 4 m s-1)以及大于和等于30 mm h-1的降水量(彩色阴影),AB为图 11的剖面位置,中图和下图分别4 km和8 km高度上的垂直速度(单位:m s-1,等值线标值从小到大依次为-2、-1.5、-1、1、3、5 m s-1)和相对系统的水平风场(单位:m s-1)Fig.7 The results of experiments US-5 and US5. The top row of figures are the vertical velocity(contours; units: m s-1; values of contours: -2,-1.5,-1,1,3 and 5 m s-1),horizontal wind speed ≥17 m s-1(barbs; full barb represents 4 m s-1),and hourly precipitation amount ≥30 mm h-1(shaded)at the 1 km level(line AB represents the cross section in Fig. 11). The middle and bottom rows of figures are the vertical velocity(units: m s-1; values of contours: -2,-1.5,-1,1,3 and 5 m s-1) and the system relative horizontal velocity(units: m s-1)at 4 and 8 km,respectively |
![]() |
图 8 图 7中US-5和US5试验南部对流区的放大图,1 km、4 km和8 km高度上的垂直速度(单位:m s-1,等值线标值从小到大依次为-5、-4、-3、-2、-1、-0.5、0.5、1、2、3、4、5 m s-1)和相对系统的水平风场(单位:m s-1)Fig.8 The enlarged southern convective areas in Fig. 7 in experiments US-5 and US5,the vertical velocity(units: m s-1,values of contours: -5,-4,-3,-2,-1,-0.5,0.5,1,2,3,4,and 5 m s-1) and the system relative horizontal velocity(units: m s-1)at 1,4,and 8 km |
US-5试验中低层(1和4 km)的下沉气流正好位于8 km高度上升气流的正下方(图 7和图 8)。US5试验中低层的下沉气流不仅与对应层次的上升气流成对出现,形成完整的对流单体,高层的上升气流也与中低层的下沉气流有一定的距离,这样下沉气流对上升气流的削弱作用减弱。这与Lilly(1979)的结论一致:垂直风切变使得单体的环流互不干扰,系统稳定发展。
积分3 h和5 h上升运动随着高度的变化趋势基本一致,但积分3 h时各个高度上的上升运动比积分5 h时稍强(图 9a,c)。下沉运动随着高度的变化在积分3 h和5 h差别较大,积分3 h,最大下沉速度位于8 km高度上,且低层3和5 km高度上的下沉运动也较大,积分5 h,8 km以下的下沉运动明显比积分3 h弱,10 km以下的最大下沉运动位于3 km和7 km,3 km和8 km高度上的下沉运动比积分3 h分别小1 m s-1和4m s-1(图 9b,d)。积分3 h和5 h,3 km以下下沉运动迅速减小,到地面下沉运动为3 m s-1,这是由于3 km以下已经转变为水平运动,且积分3 h,US5试验3 km高度上强下沉运动对应的近地面风速大(图 10)。表明系统发展阶段(积分3 h)的下沉运动比成熟阶段(积分5 h)强,有利于出现近地面大风,这与US5试验在积分3 h的近地面风速比积分5 h强的结论一致(图 6a)。另外,US5试验比US-5试验的垂直运动强,与图 5c,d的结论一致。
![]() |
图 9 US-5和US5试验积分3 h(上)和5 h(下)对流区(a、c)最大上升速度(单位:m s-1)和(b、d)最大下沉速度(单位:m s-1)的随高度的变化Fig.9 The height variation of the(a,c)maximum ascending velocity(units: m s-1) and (b,d)maximum descending velocity(units: m s-1)in convective areas in experiments US-5 and US5 at 3 h(top) and 5 h(bottom) |
![]() |
图 10 US-5和US5试验结果,图 7中AB剖面的垂直速度(等值线,红色实线为上升运动,蓝色虚线为下沉运动,单位:10-1 m s-1,等值线间隔1 m s-1)、水平与垂直风场(单位:m s-1,风速矢量的垂直分量表示实际风速的10倍)以及水平风速大于和等于17 m s-1的区域(黄色阴影)Fig.10 The results of experiments US-5 and US5. The cross sections of vertical wind velocity(contours; the red solid lines represent ascending motion,the blue dashed lines represent descending motion; units: 10-1 m s-1,contoured at 1 m s-1 interval),horizontal and vertical wind field(units: m s-1; vertical component is 10 times the real wind velocity),and the areas with horizontal wind speed ≥17 m s-1(shaded in yellow)along AB in Fig. 7 |
从不同高度气流的配置看,US5试验各个高度上的上升气流比US-5试验强(图 9a,c),且1 km高度上南部回波的上升和下沉气流呈现南北向线状分布(图 8)。 US5试验的最大上升速度在各个高度上均比US-5试验大,其中6~12 km高度二者差异达到最大(图 9a,c)。US5比US-5试验的下沉速度大,低层在3 km高度二者差异达到最大(图 9b,d),在3 km以下的差异较小,这是由于3 km以下已经转变为水平运动。积分3 h,US5试验3 km高度上下沉运动强对应的近地面风速大(图 10)。从系统南部的放大图看(图 8),在中高层(8 km),US5试验比US-5试验的上升运动和下沉运动区东西向分布更明显。这样US5试验低层的上升运动与下沉运动区为南北向分布,中高层为东西向分布,表明气流边上升边呈现顺时针旋转、边下沉边逆时针旋转,上升气流与下沉气流相互干扰较弱,有利系统的长时间维持和增强,而US-5试验气流在上升过程中的顺时针旋转较弱、下沉过程中的逆时针旋转也较弱,上升气流与下沉气流相互干扰较强,使系统难以继续生成或者有利已有系统的消散。
4.2 散度场积分3 h和5 h时,US-5试验,1 km高度以辐散为主,辐散区位于对流区的后部,对流区对应辐合气流,4 km高度辐合区对应1 km高度上的辐散区,8 km高度以辐散为主,主要辐散区对应1 km高度的辐合区(对流区)(图 11)。US5试验,1 km高度对流区对应气流辐合区,而4 km高度为气流辐散区,低层的上升运动强,有利于对流的发展和增强,1 km高度气流辐散区在辐合区的后部,4 km高度辐散区在辐合区的前部,8 km高度气流以辐合为主,且辐合辐散区分别呈南北向排列成线状(图 11)。积分3 h最北边的辐散比积分5 h强。表明在1~4 km的高度系统前部以上升气流为主、后部以下沉气流为主,上升与下沉气流互不干扰分别排列成南北向带状,积分3 h北部辐散区对应的下沉气流比积分5 h强,这与图 9b,d的结论一致。
![]() |
图 11 US-5与US5试验积分3 h和5 h,1 km、4 km和8 km高度上的散度(实线为正值,代表辐散区,虚线为负值,代表辐合区,等值线标值从小到大依次为-100、-70、-40、40、60、80,单位:10-5 s-1)和垂直速度(阴影,单位: 10-1 m s-1,等值线间隔0.5 m s-1)Fig.11 Divergence(solid lines are positive values and represent divergence areas; dashed lines are negative values and represent convergence areas,values of contours: -100,-70,-40,40,60,80,units: 10-5 s-1) and vertical velocity(shaded; units: 10-1 m s-1,contoured at 0.5 m s-1 interval)in experiments US-5 and US5 at 3 h and 5 h,at the l,4 and 8 km levels,respectively |
从垂直速度、风场、散度场、雨水比湿和涡度场的剖面图(图 10,12)看,US5试验在系统前部,积分3 h,0~3 km高度为辐合气流,3~6 km为辐散气流,这样在0~6 km形成上升运动区,在系统的后部(AB剖面上50~70 km处的低层1~2 km为水平风速大值区),在风速大值区的后部0~8 km为下沉运动区,对应这里为负涡度区,此时系统前部的正涡度区和后部的负涡度区的界限明显,互相独立。积分5 h,系统前部的上升运动加强,后部0~3 km为下沉运动,对应地面辐散和负涡度区(y方向的涡度分量)(图 12),但下沉运动比积分3 h弱,负涡度区(y方向的涡度分量)集中在上升运动后部的地面到2 km的高度上,没有像积分3 h与中高层的负涡度区连成一体,造成积分5 h近地面风速减小,而下沉运动正好位于上升运动的下方,对应这里雨水比湿较大,下沉运动可能是由于降水粒子的拖曳作用造成的,此时降水量比积分3 h时强。US-5试验系统上升运动主要集中在中高层,中低层以下沉气流为主,这与图 11的结论一致,涡度场上5 km以下以负涡度为主,5 km以上为正涡度区。US5试验系统前部的垂直环流比US-5试验更加明显,造成US5试验前部的上升运动相对集中且强度大。在系统前部,US5试验由于垂直风切变加大,y方向的正涡度区(${{\varsigma }_{y}}={\partial u}/{\partial z-{\partial w}/{\partial x}\;}\;$)比US-5试验强(图 12),特别是近地面层的正涡度区较强,同时,在系统的前部,上升运动在对流的北部是减小的[$\left( {\partial w}/{\partial y}\; \right)<0$],扭转项的作用为减弱水平涡度,对流减弱,而上升运动在对流南部是增强的$({\partial w}/{\partial y}\;>0)$,扭转项的作用是增加水平涡度,上升运动增强,这样在扭转项的作用下南部对流增强,北部对流减弱,对流向南发展。
![]() |
图 12 US-5和US5试验结果,图 7中AB剖面的散度(上,单位:10-4 s-1),y方向涡度(中,单位:10-4 s-1)和雨水比湿(下,单位:g kg-1)。Fig.12 The results of experiments US-5 and US5. The cross sections of divergence(top,units: 10-4 s-1),y-component vorticity(middle,units: 10-4 s-1),specific humidity of rainwater(bottom,units: g kg-1)along AB in Fig. 7 |
假相当位温的扰动量可以用来表征冷池,其计算方法为扰动量${{{\theta }'}_{v}}={{\theta }_{v}}-{{\bar{\theta }}_{v}}$(Coniglio and Stensrud,2001),平均量(${{\bar{\theta }}_{v}}$)是用飑线前部100 km处的$\text{10 km}\times \text{20 km}$区域的平均值表示。积分3 h,US5试验系统存在一支向后的斜升气流和一支向前的下沉气流,与下沉气流对应,在0~2 km高度上为位温的负扰动,且梯度大值区集中在水平运动的辐合区,对应这里出现水平风速的大值区,US-5试验也存在两支气流(斜升气流和下沉气流),但下沉气流位于上升气流的下方,对上升气流有削弱作用,这与图 7的结论一致,下沉运动区在0~2 km高度上也存在位温负扰动,但比US5试验扰动强度小,造成的近地面大风较弱。积分5 h时,US5试验比积分3 h系统的斜升气流增强,后部的下沉气流减弱,下沉运动区主要位于斜升气流的正下方,可能是降水增强造成的,下沉运动造成位温负扰动比积分3 h弱,此时降水量增加而近地面风速小,US-5试验,斜升运动和位温负扰动比积分3 h减弱,比US5试验也弱。增加垂直风切变,对流发展高度高,旧单体消散快,假相当位温负扰动集中在新单体后部,冷池高度高且相对集中在新单体后部(图 13)。减小垂直风切变,对流发展高度低,旧单体消散慢,假相当位温负扰动分布相对均匀,冷池高度低且分布较分散,不利于新单体的生成(图 13)。计算冷池的强度:
${{c}^{2}}=-2\int_{z=0}^{z=H}{B\text{d}z}$ | (3) |
![]() |
图 13 同图 10,但为组合雷达反射率(黑色实线,单位:dBZ)、假相当位温负扰动(紫色虚线,单位:K,间隔2 K)、纬向风(阴影,单位:m s-1)和风场矢量(箭头,单位:m s-1,风速矢量的垂直分量表示实际风速的10倍)Fig.13 As in Fig. 10,but for the composited radar reflectivity(black line; units: dBZ),negative pseudoequivalent potential temperature turbulence(purple dashed isolines with 2 K interval; units: K),zonal wind(color shading; units: m s-1),and wind field(vectors; units: m s-1; vertical component is 10 times the real wind velocity) |
![]() |
图 14 US-5和US5试验(a)冷池强度(单位:m s-1)和(b)高度(单位:m)的时间变化Fig.14 Temporal variation of the(a)intensity(units: m s-1) and (b)height(units: m)of the cold pool in experiments US-5 and US5 |
本文利用我国江淮地区一次发生强对流的探空和WRF模式,进行了垂直风切变对对流强度、组织形态影响的理想数值试验。通过增加和减小整层、低层和中层的风切变的试验,发现:改变整层垂直风切变对对流的强度和组织类型影响最大,其次是改变中层的风切变强度,而改变低层的风切变影响最小;增加整层垂直风切变,对流发展的强度强,且易组织成线状,减少垂直风切变,对流发展强度弱,对流呈分散的状态;增加垂直风切变,上升速度增强,最大小时降水量增强,下沉速度增强,近地面风增强。
从垂直速度、风场、散度场和垂直结构以及冷池分析了垂直风切变对系统强度和组织影响的可能原因。增加垂直风切变,系统的上升运动与下沉运动的旋转增强,二者之间的相互干扰减弱,使系统能够长时间维持,且系统的上升运动与下沉运动增强。从涡度场看,垂直风切变增强,涡度增强,这样在系统的前部,由于强上升气流的作用,在扭转项的作用下,系统的垂直涡度增强,加强原来的上升运动,相反在系统后部,由于下沉气流的作用,在扭转项的作用下,系统的垂直涡度减小,加强原来的下沉运动,这样造成系统强度增强。垂直风切变增加,冷池高度高且相对集中,垂直风切变减少,冷池高度与冷池强度均降低。因此,垂直风切变增加,加强了系统后部的下沉运动,此外冷池在系统后部的集中度高表明系统的组织性强,不易呈现分散的状态。
研究还发现增加整层垂直风切变,近地面风速和降水量增加,且降水量的大值出现在近地面大风大值之后。在系统的发展阶段,出现近地面大风的原因是系统的正负涡度区相互独立,上升运动与下沉运动互不干扰,后部的下沉运动强,下沉运动到达近地面后形成强辐散造成近地面大风,成熟阶段,系统后部的正涡度区减弱,没有与中高层的负涡度大值区连成一体,下沉运动减弱,近地面大风减弱。在系统的成熟阶段,出现降水量大值的原因是,随着上升运动的不断增强,雨水比湿增大,下降到地面形成强降水,因此强降水产生是一个上升运动不断加强或维持造成雨水比湿不断增强的过程。
本文的研究虽然采用的是实际观测的探空作为背景场,但采用的是理想模式,只考虑了垂直风切变一个因素对对流的影响,实际中影响对流发展和组织的因素相对较多且复杂,今后应考虑更多环境因子对对流的影响。
[1] | Alexander G P, Young G S. 1992. The relationship between EMEX mesoscale precipitation feature properties and their environmental characteristics[J]. Mon. Wea. Rev., 120:554-564. |
[2] | Bluestein H B, Jain M H. 1985. Formation of mesoscale lines of precipitation:Severe squall lines in Oklahoma during the spring[J]. J.Atmos. Sci., 42(16):1711-1732. |
[3] | Coniglio M C, Stensrud D J. 2001. Simulation of a progressive derecho using composite initial conditions[J]. Mon. Wea. Rev., 129:1593-1616. |
[4] | Ferrier B S, Simpson J, Tao W K. 1996. Factors responsible for precipitation efficiencies in midlatitude and tropical squall simulations[J]. Mon. Wea.Rev., 124:2100-2125. |
[5] | Fovell R G, Ogura Y. 1989. Effect of vertical wind shear on numerically simulated multicell storm structure[J]. J. Atmos. Sci., 46:3144-3176. |
[6] | Johnson R H, Bresch J F. 1991. Diagnosed characteristics of precipitation systems over Taiwan during the May-June 1987 TAMEX[J]. Mon. Wea.Rev., 119(11):2540-2557. |
[7] | Johnson R H, Aves S L, Ciesielski P E, et al. 2005. Organization of oceanic convection during the onset of the 1998 East Asian summer monsoon[J]. Mon. Wea. Rev., 133(1):131-148. |
[8] | Lafore J P, Moncrieff M W. 1989. A numerical investigation of the organization and interaction of the convective and stratiform regions of tropical squall lines[J]. J. Atmos. Sci., 46(4):521-544. |
[9] | LeMone M A, Zipser E J, Trier S B. 1998. The role of environmental shear and thermodynamic conditions in determining the structure and evolution of mesoscale convective systems during TOGA COARE[J]. J. Atmos.Sci., 55:3493-3518. |
[10] | Lilly D K. 1979. The dynamical structure and evolution of thunderstorms and squall lines[J]. Annual Review of Earth and Planetary Sciences, 7:117-161. |
[11] | Moncrieff M W. 1978. The dynamical structure of two-dimensional steady convection in constant vertical shear[J]. Quart. J. Roy. Meteor. Soc., 104:543-567. |
[12] | Montmerle T, Lafore J P, Redelsperger J L. 2000. A tropical squall line observed during TOGA COARE:Extended comparisons between simulations and Doppler radar data and the role of midlevel wind shear[J]. Mon. Wea. Rev., 128:3709-3730. |
[13] | 潘玉洁, 赵坤, 攀益农. 2008. 一次强飑线内强降水超级单体风暴的单多普勒雷达分析[J]. 气象学报, 66(4):621-636. Pan Yujie, Zhao Kun, Pan Yinong. 2008. Single-Doppler radar observation of a heavy precipitation supercell on a severe squall line[J]. Acta Meteor. Sinica(in Chinese), 66(4):621-636. |
[14] | Parker M D, Johnson R H. 2000. Organizational modes of midlatitude mesoscale convective systems[J]. Mon. Wea. Rev., 128(10):3413-3436. |
[15] | Rotunno R, Klemp J B, Weisman M L. 1988. A theory for strong, long-lived squall lines[J]. J. Atmos. Sci., 45:463-485. |
[16] | Schlesinger R E. 1978. A three-dimensional numerical model of an isolated thunderstorm. Part I:Comparative experiments for variable ambient wind shear[J]. J. Atmos. Sci., 35:690-713. |
[17] | Schumacher R S, Johnson R H. 2005. Organization and environmental properties of extreme-rain-producing mesoscale convective systems[J]. Mon. Wea. Rev., 133:961-976. |
[18] | 孙建华, 郑淋淋, 赵思雄. 2014. 水汽含量对飑线组织结构和强度影响的数值试验[J]. 大气科学, 38(4):742-755. Sun Jianhua, Zheng Linlin and Zhao Sixiong. 2014. Impact of moisture on the organizational mode and intensity of squall lines determined through numerical experiments[J].Chinese J. Atmos. Sci.(in Chinese), 38(4):742-755, doi:10.3878/j.issn.1006-9895. 2013.13187. |
[19] | Weisman M L, Klemp J B. 1982. The dependence of numerically simulated convective storms on vertical wind shear and buoyancy[J]. Mon. Wea.Rev., 110:504-520. |
[20] | Weisman M L, Rotunno R. 2004. "A theory for strong long-lived squall lines" revisited[J]. J. Atmos. Sci., 61:361-382. |
[21] | Weisman M L, Klemp J B, Rotunno R. 1988. Structure and evolution of numerically simulated squall lines[J]. J. Atmos. Sci., 45:1990-2013. |
[22] | Xue M. 2000. Density currents in two-layer shear flows[J]. Quart. J. Roy.Meteor. Soc., 126:1301-1320. |
[23] | 郑淋淋, 孙建华. 2013. 干、湿环境下中尺度对流系统发生的环流背景和地面特征分析[J]. 大气科学, 37(4):891-904. Zheng Linlin, Sun Jianhua. 2013. Characteristics of synoptic and surface circulation of mesoscale convective systems in dry and moist environmental conditions[J]. Chinese J. Atmos. Sci.(in Chinese), 37(4):891-904, doi:10.3878/j.issn.1006-9895.2012.12090. |
[24] | Zheng L L, Sun J H, Zhang X L, et al. 2013. Organizational modes of mesoscale convective systems over central East China[J]. Wea.Forecasting, 28(5):1081-1098. |