第 63 卷第 4 期 2005 年 8 月气象学报 ACTA METEOROLOGICASINICA Vol.63,No.4 August 年 8 月 23 日华北强风暴动力机制的数值研究龚佃利 ( 山东省气象科学研究所, 济南,250031; 中国海洋大学海洋环境学院,

Size: px

Start display at page:

Download "第 63 卷第 4 期 2005 年 8 月气象学报 ACTA METEOROLOGICASINICA Vol.63,No.4 August 年 8 月 23 日华北强风暴动力机制的数值研究龚佃利 ( 山东省气象科学研究所, 济南,250031; 中国海洋大学海洋环境学院,"

靓常
5 years ago
Views:

1 第 63 卷第 4 期 2005 年 8 月气象学报 ACTA METEOROLOGICASINICA Vol.63,No.4 August 年 8 月 23 日华北强风暴动力机制的数值研究龚佃利 ( 山东省气象科学研究所, 济南,250031; 中国海洋大学海洋环境学院, 青岛,266003) 吴增茂傅刚 ( 中国海洋大学海洋环境学院, 青岛,266003) 摘要利用非静力平衡模式 (MM5V3.5) 对 2001 年 8 月 23 日发生在华北地区的一次强对流风暴过程进行了数值模拟, 在取得合理模拟结果的基础上, 着重分析了风暴发生的热力条件太行山地形动力作用和风暴中尺度结构特征, 并依据倾斜涡度发展理论, 分析了强风暴的发展机制结果表明 : 本次风暴是在以西北气流为主导的环流背景下产生的, 低层增温, 中高层降温和整层增湿是造成层结不稳定的重要原因上游地区对流引发的干冷下沉气流沿太行山背风面下滑形成的下坡风是触发强风暴发生的直接动力机制湿位涡分析表明, 风暴发生区具有明显的等熵面倾斜, 对流层中低层狆 m1<0 区域, 同时狆 m2>0, 满足倾斜涡度发展的条件, 对强风暴的发展具有一定指示性风暴发展早期, 其垂直方向次级环流可能与新雷暴的产生和雷暴的跳跃式传播有关太行山地形引起的边界层风场的变化, 包括下坡气流和边界层中尺度辐合线对风暴的触发组织和移动发挥着重要作用关键词 : 强对流风暴, 数值模拟, 动力机制, 地形作用, 倾斜涡度发展 1 引言华北地区西靠黄土高原东临黄渤海中部为华北平原, 太行山燕山山脉位于华北平原西侧和北侧, 具有多样的下垫面地形条件, 是中国夏季飑线 MCC 等强对流天气多发区, 常有大风冰雹等灾害出现早期对华北强对流天气的研究主要从系统发生的环流背景形成机制和演变特征进行天气学分析, 初步给出了该地区强对流天气发生的气候学特 [1,2] 征山东位于华北南部, 许多强对流天气与华北上游地区对流活动和太行山地形作用密切相关 [3] 据统计, 华北飑线的初生区主要集中在华北西北部的阴山太行山北部山区 ; 其中 30% 源于太行山脉的飑线向东南方向移动, 影响河北东南部和山东 [4] 地区朱官忠等通过对典型 MCC 过程的诊断分析, 对盛夏华北南部集中出现 MCC 的原因进行了探讨, 认为华北南部特殊地理条件产生的准定常暖区及相应热低压和鲁中有利地形对暖湿气流的抬升作用, 是华北南部冀鲁豫交界处多有 MCC 发生的重要原因上述研究表明, 华北特定地形在强对流天气发展中具有重要作用, 许多人通过理论分析或数值模拟方法对此进行了研究并对地形强迫产生背风波 [5~6] 与下坡风的认识较为充分, 即在稳定层结大气, 一定强度的气流越过山脊后可形成准静止背风波和下坡风, 山脉越高背风坡度越大低层大气越稳定, [3] [7] 则产生的下坡风越强李鸿洲张群等利用中国科学院大气物理研究所 LASG 建立的 η 坐标有限区域模式, 对华北强飑线生成过程中燕山地形及边界层中尺度辐合线的作用进行了数值研究, 认为初稿时间 :2003 年 10 月 14 日 ; 修改稿时间 :2004 年 8 月 10 日资助课题 : 山东省气象局重点课题中尺度模式业务系统及高性能计算环境应用山东省中尺度数值预报系统的建立及应用和山东省自然科学基金课题 (Y2003E01) 作者简介 : 龚佃利, 男,1969 年生, 硕士, 高级工程师, 主要从事中尺度数值模拟分析研究 gdlj@sohu.com

2 4 期龚佃利等 :2001 年 8 月 23 日华北强风暴动力机制的数值研究 505 在高层西北气流控制, 有较强冷空气侵入时, 地形强迫引起的低层强烈锋生是华北飑线形成和发展的主要动力过程, 边界层辐合线对飑线的发展具有重要 [8] 作用范广洲等利用 NCAR RegCM2 模式就华北西部和北部山脉地形对华北夏季降水的影响进行了数值试验, 结果表明, 当降低地形高度时, 华北夏季降水将明显减小但上述研究要么是基于理想地形条件下的理论研究, 要么受模式性能限制, 尚不能很好地揭示实际强风暴过程中上游地区对流活动太行山地形在下游强对流天气发展演变过程中的作用机制研究气旋性涡度的发展机制是分析强天气产生原因的重要内容位涡是一个将热力和动力作用结 [9] 合在一起的物理量吴国雄等从完整的原始运动方程出发, 导出了精确形式的湿位涡方程, 并证明其守恒性, 提出了倾斜涡度发展 (SVD) 理论, 指出涡旋易于在等熵面陡立的地方发展近几年, 湿位涡理论较多地应用于暴雨诊断气旋的移动发展等方面, [10~12] 也有应用于冰雹天气中的诊断分析, 取得了较好的结果本文利用非静力学模式 MM5V3.5, 对 2001 年 8 月 23 日发生在华北南部的强风暴发展演变的热力动力学机制进行了数值分析, 特别是对太行山地形的作用进行了研究, 并从湿位涡角度分析了强风暴的倾斜涡度发展过程这对于认识华北强对流风暴演变的动力过程, 为更好地监测预报该地区强风暴灾害提供依据 2 天气概述和系统演变 2001 年 8 月 22~24 日, 受西北气流和冷涡影响, 华北地区连续多日出现强对流天气, 其中 23~ 24 日, 山东连续两天出现罕见的大范围降雹和强雷雨大风天气 23 日 05:00( 北京时, 下同 ) 左右, 强对流风暴从鲁西北侵入山东后, 横扫鲁中山区和鲁东南地区,23 日中午前后从鲁东南沿海移出本次强风暴在山东境内持续影响超过 6h, 移向基本为 NW SE 走向, 移速约 65km/h, 期间, 山东境内有 20 多个市县遭受严重大风冰雹灾害, 直接经济损失 13 亿多元 2.1 环流形势图 1 为根据 22 日 20:00T106 分析同化数据绘制的 500,850hPa 环流形势其中,500hPa 图上, 在东北和阿尔泰山以东地区为两个低涡系统, 且后者配合 -16 的冷中心, 此后, 该冷涡向东南方向移动, 成为造成山东连续两天强对流天气的冷空气源分析 850hPa 形势图可见,22 日 08:00, 黄河中下游地区为高于 20 的暖区, 华北南部低层有大范围逆温区, 随着高层冷空气的侵入, 导致大气层结不稳定 22 日 20:00( 图 1b), 原位于蒙古境内低压环流东移与东北冷涡环流相接, 在内蒙古中西部形成弱横槽, 引导冷空气和正涡度平流向华北南部移动同日, 华北地区地面为前期 0111 号台风减弱后形成的弱低压系统 ( 图略 ),22 日 14:00, 河北东南部和鲁西北地区地面气温超过 35 分析表明, 中纬度高空冷涡南侧的西北气流与低层横槽相配合是本次强风暴发生前的主要环流形 [13] 势, 但与以往总结的山东冰雹天气形势有所不同 : 高空冷涡位置偏北无横槽, 低层横槽前期也不明显, 形势并不利于山东降雹因此, 这种高层以西图 1 8 月 22 日 20:00500hPa(a) 和 850hPa(b) 环流形势分析 ( 实线为等高线, 间隔 2dagpm; 虚线为等温线, 间隔 2 ) Fig.1 Analysisat500hPa(a)and850hPa(b)showinggeopotentialheight(solidevery2dagpm) andtemperature(dashedevery2 )for20:00lst22august2001

3 ５０６气象学报６３卷北气流控制为主的环流形势常常使预报员容易忽与上游地区对流活动的关系视强对流天气发生的可能性给业务预报带来了很大困难２３日００３２图２ａ如前分析位于北京呼和浩特一带的强雷暴云团正处于强盛阶段华北南部２２红外云图演变特征随着高空西北冷涡云系东移受地面弱冷锋和为中高云覆盖区２３日０２３２图２ｂ覆盖京津冀地区的对流云开始减弱但在河北东南部和太行山午后热对流的共同作用２２日下午内蒙古中西部地区有对流云发展并局地发展成两条东北西南地区开始有雷暴出现对流云团东移加强并于夜间向的对流云带图２ｃ此后这两块对流云合并迅造成华北北部强对流天气由２３日００３２至１０３２速加强为强风暴并向东南方向移动０５００风暴红外云图演变图２可以清楚地看到本次强风暴侵入山东其云顶亮温ＴＢＢ低于４０云区呈椭图２８月２３日强风暴红外卫星云图ＴＢＢ演变ａ００３２ｂ０２３２ｃ０４３２ｄ０６３２ｅ０８３２ｆ１０３２Ｆｉ２Ｔｈｅｅｖｏｌｕｔｉｏｎａｒｒｏｃｅｓｓｏｆｔｂｂｉｎｇｍｓ５ｉｎｆｒａｒｅｄｉｍａｇｅｓｇｙｐｆｒｏｍ００３２ＬＳＴｔｏ１０３２ＬＳＴ２３Ａｕｇｕｓｔ２００１

4 4 期龚佃利等 :2001 年 8 月 23 日华北强风暴动力机制的数值研究 507 圆型, 长轴为南北向, 强中心 (TBB< -50 ) 位于风暴南部, 大风冰雹等强天气出现在风暴强中心附近, 而风暴北部为雷阵雨天气 08:32( 图 2e), 对流云强中心移到鲁中山区, 不仅造成博山莱芜沂源等地强风雹灾害, 且伴有局地强降水, 其中, 博山站 08:00~09:00 降水量达到 48.7mm 10:32( 图 2f) 后, 对流风暴逐渐减弱, 但地面的大风冰雹不减, 11:32 后, 风暴中心自鲁东南沿海入海纵观本次强对流风暴的演变可见, 它起源于华北北部的对流活动, 发展迅速, 虽然风暴强中心范围较小, 云顶 TBB 最低也只有 -58, 但在其影响山东省的 6h 内, 恰处于对流云体发展的强盛时期, 因而造成大范围风雹灾害 3 数值模拟方案和结果分析 3.1 模拟方案采用双向作用 3 层套网格, 格距分别为 81,27, 9km, 垂直 σ 坐标分 23 层大网格范围覆盖了直接影响风暴产生的天气尺度系统, 这有利于提高对细网格区域内小尺度系统的模拟能力 ; 内层网格定位则主要为了跟踪风暴发生发展的全过程, 并便于对太行山地形作用进行分析采用 NCAR 的 10,5,2 地形及 25 类植被资料形成模式地形和下垫面数据图 3 给出了第 3 层 9km 嵌套区域的模式地形分析剖线和站点分布示意图, 山阴日照剖线与风暴移动路径基本一致, 以下的剖面分析均沿此剖线进行以 2001 年 8 月 22 日 08:00 为初始时刻, 采用北京 T106 同化分析数据作为背景场, 并用风暴期间的高空地面资料对其进行订正, 构成模式的初边值条件采用时变侧边界 MRF(Hong Pan,1996 年 ) 行星边界层参数化方案和云辐射处理方案对于云降水物理过程,3 层区域均考虑霰和冰相过程, 选用 Reisner 霰方案 ; 另外, 第 1 层网格选用 Bets Miler 方案 (1993 年 ), 第 2,3 层选用 Grel 方案 (1993 年 ), 以反映次网格尺度的对流过程海温 SST 利用 NOAA 周平均资料图 3 模式第 3 层网格区域地形 ( 虚线,m) 分析剖线 ( 直线 ) 和站点示意 Fig.3 Theterrain (dashed,m)andobservationstationsin9km MM5domain (ThebeelinefromShanyitoRizhaoshowsthecross sectionalongthestormtrack, thedistancesfromshanyinarelabeledatthecross section)

5 508 气象学报 63 卷 3.2 模拟结果分析降水图 4 给出了强风暴影响山东期间 23 日 06:00 ~12:00 模拟降水和实况降水量对照可见, 两者的雨区分布和中心位置基本一致, 其中模拟最大降水 65 mm, 位于博山附近, 实际观测博山站降水 48.9mm, 为本次过程最大雨量站另外, 模拟的逐时降水范围降水中心及其强度变化与实况基本一致, 只是 09:00~11:00 降水中心值偏大图 4 23 日 06:00~12:00 模拟 (a) 和实况 (b) 降水量 (mm) Fig.4 Spatialdistributionofmodelprecipitation(a)andobservedprecipitation(b) from06:00lstto12:00lst23august2001(unit:mm) 地面气象要素演变图 5 为模拟的济南站 8 月 23 日 02:00~14:00 地面 2m 温度混合比 10 m 风速和地面气压波动演变可见, 上述各量在强风暴飑锋影响时均出现明显的突变, 即地面风速骤增地面温度和水汽混合比迅速下降, 地面气压波动出现明显陡升陡降特征, 与已知的飑线影响时天气现象变化一致云场和风场图 6( 见彩页图 ) 给出了模拟的 8 月 23 日 06:00 雷达综合反射率 ( 由模拟全部水成物含量和下落速 [14] 度计算, 并取各层最大值 ) 10 m 风场和滨州 CINRAD/SC 观测的 06:23 雷达回波强度图可见, 模拟的云场结构位置与实际卫星云图 ( 图 2d) 雷达观测 ( 图 6b) 非常相似 ; 地面雷暴大风也很清楚, 与实况一致总之, 模式较好地再现了强风暴发展演变过程, 为进一步分析风暴形成的动力机制及其中尺度结构奠定了基础 4 强风暴发展机制分析图 5 模拟的济南站 8 月 23 日 02:00~14:00 2m 温度 ( 犜 2m) 地面混合比 ( 犙 s) 10m 风速 ( 犞 10m) 和地面气压波动 ( 犘 Pa) 时间演变 ( 各量采用一致的纵坐标, 单位见图注 ) Fig.5 Simulated2mtemperature( 犜 2m, ), surfacemixingratio( 犙 s,g/kg),10m wind speed( 犞 10m,m/s)andsurfacepressure perturbation( 犘 Pa,hPa)atJinanstation from02:00lstto14:00lst23august 热力条件图 7 给出了济南站 22 日 08:00( 虚线 ) 和 23 日 02:00( 实线 ) 探空曲线及垂直风廓线可见,22 日 08:00, 济南 800hPa 以下空气湿度较小, 层结廓线呈倒 V 型,966 hpa 以下有 2 逆温层 ;23 日 02:00, 整层大气温度露点差 ( 犜 - 犜 d) 较 08:00 明显减小垂直风廓线显示,23 日 02:00 前, 济南 850 hpa 以下基本为西南风, 中高层为西北风,850~700 hpa 之间具有明显的风切变 23 日 08:00, 释放探空气球后, 恰遇风暴袭击, 仅有低层风记录, 可见低层已转偏北风, 且风速明显增大

6 4 期龚佃利等 :2001 年 8 月 23 日华北强风暴动力机制的数值研究 509 图 7 济南站 8 月 22 日 08:00( 虚线 ) 23 日 02:00 ( 实线 ) 探空曲线和垂直风廓线 Fig.7 Atmosphericsounding (skew 犜 log 狆 format)atjinanfor08:00 LST22August(dashed)and02:00LST23August(solid),thewind vectorofeachheightlayerfor08:00lst22,20:00lst22, 02:00LST23and08:00LST23August2001areshown 分析对流风暴发生前济南探空站大气热力参数变化 ( 表 1) 可见,22 日白天, 济南附近对流有效位能 (CAPE) 增大, 对流抑制能量 (CIN) 减小, 自由对流高度 (LFC) 降低, 大气层结积累了大量不稳定能量从 22 日 20:00 至 23 日 02:00, 低层增温高层降温和整层增湿特征明显,850hPa 与 500hPa 间相当位温差 ( θe) 和沙氏指数 ( 犛犻 ) 表明, 大气层结处于强对流不稳定状态, 极易引发强对流天气的发生表 1 济南探空站大气层结热力参数变化 Table1 ThevarietyofthermodynamicindexcalculatedfromJinansoundingdata 时间 θe(k) 犛犻 CAPE(J/kg) CIN(J/kg) LFC(hPa) 温度 ( ) 比湿 (g/kg) 500hPa 850hPa 500hPa 850hPa 22 日 08: 日 20: 日 02: 根据模式和实况资料分析,22 日白天,850hPa 以下干暖的偏南气流 400hPa 以上干冷偏西气流和中层湿冷西北气流在华北南部到山东地区层叠, 是造成该地区层结不稳定的主要原因特别是 22 日白天,650~500hPa 高度层有湿冷平流 ( 犜 - 犜 d< 4,θe<328K) 移向华北南部地区, 致使未来的雷暴区中高层湿度条件好于低层, 这与大多数华北强对流天气发生时, 低层暖湿, 中层干燥的特征不 [4] 同因此, 本过程中层冷云存在的辐射效应是区别于其他强对流天气的重要特征 4.2 太行山的地形作用强对流不稳定层结为对流的发生准备了必要的能量条件, 但对流何时何地发生, 则取决于触发对流不稳定能量释放的机制图 8 给出了风暴发展前期 CAPE 分布和对流层低层 925hPa 风场可见,23 日 03:00( 图 8a), 山东境内大部分地区 CAPE 高于 1000J/kg, 其中鲁中山区鲁北地区超过 2000J/ kg 太行山东侧河北中部有一条东北西南向大风带, 风向与太行山近乎垂直, 最大风力近 20m/s, 较前一时次增加近 10 m/s 根据卫星云图和山阴日照剖面狏, 狑合成风场 ( 图 9a) 演变可见, 该大风带源于太行山区对流产生的干冷下沉气流, 犜 - 犜 d>10, 该气流沿太行山背风坡冲击下来, 从低层侵入其前方的暖湿气团, 是启动本次强对流风暴

7 510 气象学报 63 卷发展的最重要机制 06:00( 图 8b), 来自太行山的干冷气流推进到鲁西北地区, 与其前方的西南暖湿气流形成中尺度辐合线, 促使对流迅速发展同时可以看到, 在太行山东侧, 除了来自太行山的下坡气流外, 来自燕山山脉的偏北气流也沿太行山东侧南下侵入到华北平原南部和鲁西北地区分析低层风场演变可见, 这两支气流具有顺太行山山体右偏呈气旋性弯曲南下的特征 : 冷空气流如河水, 而太行山如河岸造成这种风场的原因, 一方面受河北境内中低压系统的影响, 另一方面受柯氏力作用对历史降雹过程的统计分 [13] 析表明, 自黄土高原太行山脉快速移动的冷锋或飑线可造成山东大范围强对流天气, 低层冷空气的移动路径也具有类似现象本次模拟所揭示的低层中尺度风场特征, 可为研究太行山华北平原特定地形条件下, 冷空气的侵入路径与华北南部冰雹天气发生规律之间的关系提供了新的启示图 8 23 日 03:00(a),06:00(b)925hPa 风场和 CAPE 分布 ( 实线为 CAPE, 单位 :J/kg; 虚线阴影区为地形, 单位 :m) Fig.8 Thedistributionofmodelwindvectorat925hPaandCAPE(solid,unit:J/kg)for03:00LST(a) and06:00lst (b)23august2001(theshadedareaisterainin9km MM5domain,unit:m) 4.3 倾斜湿位涡发展机制根据文献 [9], 在静力近似, 并假定垂直速度的水平变化比水平速度的垂直切变小得多, 狆坐标下的湿位涡 (MPV) 表达式为 : 狆 m =- 犵 [ 狏犳 +ζp ] θe θe 狌 θe + 犵 - 犵狆狆狓狆狔可以证明, 绝热无摩擦的饱和大气具有湿位涡守恒特征, 即狆 m= 常数将其写成分量形式 : 狆 m1 =- 犵 ( ζ p + 犳 ) θe 狆狌 θe 狏 θe 狆 m2 =- 犵 + 犵狆狔狆狓其中 : ζ p 为绝对涡度垂直分量,θe 为相当位温狆 m1 是湿位涡的垂直分量, 为湿正压项, 北半球大气对流不稳定时, 狆 m1<0; 大气对流稳定时, 狆 m1>0 狆 m2 是湿位涡在等压面上的水平分量, 为湿斜压项当等熵面与等压面近似平行时, 狆 m1 狆 m1 在湿位涡守恒的制约下,θe 面倾斜越大, 气旋性涡度发展越剧烈, 这种涡度的增长称为倾斜涡度发展 (SVD),SVD 发展的条件为 : 犆 d = 狆 m2/ θe >0 狆因此, 在湿绝热无摩擦大气中, 系统涡度的发展由大气层结稳定度斜压性和风的垂直切变等因素所决定湿等熵面的倾斜, 大气水平风垂直切变或湿斜压性增加, 将导致垂直涡度的显著性发展图 9 给出的山阴日照大气层结和狆 m1, 狆 m2 剖面较好地揭示了对流风暴倾斜涡度的发展过程 23 日 03:00, 风暴发生前 ( 图 9a), 太行山区附近低层为中性或稳定层结, 与图 8a 一致, 山脉东坡已出现干冷下坡风, 犜 - 犜 d>8, 低层风切变增大 ; 华北平原南部 (250~350km) 低层为暖湿层结, 等熵面具有明显的倾斜特征图 9b 显示, 对应自太行山东移的干冷气团狆 m1>0pvu 1, 气团稳定度较大, 而河 1 PVU=Potential Vorticity Unit

8 4 期龚佃利等 :2001 年 8 月 23 日华北强风暴动力机制的数值研究 511 北东南部山东地区中低层狆 m1<0pvu, 且负值中心在鲁中山区 800hPa 高度附近, 表明上述地区为强对流不稳定层结图 9c 显示, 狆犿 2 在 400hPa 以上高度明显小于狆 m1, 但低层与狆 m1 具有相当的量 θe 级, 对应 >0 的区域, 狆 m2>0, 显示斜压不稳定较狆强, 满足 SVD 发展的条件因此, 随着低层冷暖气团交界面倾斜度加大水平风垂直切变和湿斜压性增加, 并有辐合上升运动时, 可使倾斜涡度强烈发展, 引起对流的发生 06:00( 图 9d), 从太行山冲击下来的干空气团将其前方的暖湿气流抬升起来, 形成高能舌 Ω 型锢囚, 暖气团中强上升运动可达到 400hPa 高度以上图 9e 所示风暴移动前方 ( 山东境内 ) 的中低层狆 m1 仍为负值狆 m2 大多为正值 ( 图 9f), 犆 d>0, 因此, 位于风暴移动前方的鲁中山区鲁东南地区具有很好的湿斜压涡度发展条件, 致使强风暴在移动过程中继续加强, 造成严重大风冰雹灾害风暴影响区域, 由于冷空气的侵入, 底层狆 m1 演变为正值区, 层结趋于稳定, 狆 m2 出现正负中心交替排列特征, 数值与狆 m1 接近, 甚至大于后者狆 m2 的变化表明, 由于强对流风暴的强非地转非静力平衡运动特征, 特别是湿下沉运动出现后造成的质量强迫, 使湿位涡守恒性被破坏, 导致风暴低层的位涡发生异常, 表现出与大尺度天气系统不同的位涡特征, 这与文献 [9,10] 的分析一致图 9 23 日 03:00(a,b,c) 和 06:00(d,e,f) 山阴日照垂直层结和狆 m1, 狆 m2 剖面 ( 图 a,d 实线为 θe, 间隔 3K; 虚线为犜 - 犜 d, 间隔 4, 风场为狏, 狑 ( 放大 10 倍 ) 合成, 图底阴影区为太行山鲁中山区地形 ; 图 b,e 为狆 m1, 单位为 0.1PVU; 图 c,f 为狆 m2, 单位为 0.1PVU) Fig.9 (a)simulatedverticalcross section(seefig.3)ofequivalentpotentialtemperature (θe,solidlinesevery3k),verticalwindfield (m/s)and 犜 - 犜 d(dashedlines every4 )at03:00lst23august2001,theshadedareasalong 狓 axisshowingthe terrainoftaihangmountain(0-200km)andluzhongmountain( km),(b)the 狆 m1 distribution (unitin0.1pvu ),(c)the 狆 m2distribution (unitin0.1pvu), (d)-(f)sameas(a)-(c)butfor06:00lst23august2001

9 512 气象学报 63 卷 4.4 风暴中尺度结构及其与边界层地形关系图 10 给出了对流风暴低层与高层的垂直运动场 23 日 05:00,925 hpa 垂直速度场显示 ( 图 10a), 位于冀鲁交界处的上升下沉区呈东西带状, 与雷达回波观测的多单体群的分布基本相符, 配合温湿度场分布, 充分反映出了风暴 β 中尺度的结构特征地面强风带与上升运动区一致, 而风场中尺度辐合线位于上升区前侧, 这种结构有利于在对流风暴前部激发新的对流单体 08:00, 强对流风暴推进到鲁中山区北侧, 对流风暴呈现出多个对流中心 ( 图 10b), 其中位于鲁中山区南侧的上升运动带与低层中尺度辐合线重合, 预示着该对流带将减弱 09:00 后,925hPa 表现出下击暴流的风场特征, 与下沉运动对应的是地面雷暴高压和强降水, 位于风暴后部的弱下沉区对应层状云区和雷阵雨图日对流风暴低层和高层垂直速度场演变 (a,b 为 05:00,08:00925hPa;c,d 为 05:00,08:00400hPa; 阴影区为狑 >0.1m/s 的上升运动区 ) Fig.10 Theevolutionoftheverticalvelocityandwindfieldat925hPa(a,b)and400hPa (c,d)for05:00lst(a,c)and08:00lst(b,d)(theshadedareasshow 狑 >0.1m/s) 400hPa 反映的高层上升运动强而集中, 下沉运动弱且位于上升运动区边缘 ( 图 10c,d) 08:00 后, 成熟阶段的对流云中高层只维持一个强对流中心, 模拟的中心最大上升速度为 4.5 m/s 高层上升下沉区的分布结构, 反映出高层对流云的运动特点, 即强对流云的上升气流受对流层顶稳定层结的 [15] 阻挡, 向四周辐散下沉所致这与陈泽宇对一次深对流过程的数值模拟结果一致分析风暴垂直运动场演变可见, 对流风暴影响华北平原南部地区时, 下垫面较为平坦, 近地层风场较好地反映出了风暴的风场结构 ; 而风暴移动到鲁中山区后, 受复杂地形的影响, 近地层风场变得较为复杂, 而中高层却可以很好地分析出对流核心的演变 4.5 次级环流作用图 11a 为模拟的 23 日 05:00 沿 E 经向 θe, 雷达回波和垂直环流结构剖面可见, 对流初期, 来自对流云中后部的中层冷气流从低层侵入对流云前暖湿云团, 促使暖湿空气上升, 在对流云前部

10 4 期龚佃利等 :2001 年 8 月 23 日华北强风暴动力机制的数值研究 513 形成一支中低层次级环流 ; 对流云前的上升气流倾斜上升, 在 400hPa 以上形成前后两支上升气流 : 向前的一支上升到 200hPa 后下沉, 在 200~400hPa 间形成另一支次级环流 ; 向后的上升气流与来自对流云后部中层的弱上升气流汇合, 形成对流云后部的高空出流, 与模拟的砧云回波对应图 11 模拟 23 日 05:00(a) 和 08:00(b) 强风暴中心次级环流结构 (a. 沿 E 经线剖面, 风场由狏, 狑合成, 狑放大 10 倍 ; 阴影区为模拟雷达综合反射率大于 1dBz 的云体, 虚线为等 θe 线,b. 沿 E 经线剖面, 狑为实际模拟值, 其他同 a) Fig.11 (a)simulatedverticalcross sectionofθe(dashedlinesevery5k),verticalwindfield (composedof 狏 and 狑, 狑 isenlarged10times,unit:m/s)andradarreflectivity (shaded,>1dbz)along116.6 Efor05:00LST23August2001,(b)sameas(a) expectalong117.6 Eandfor08:00LST23August2001 分析 23 日 05:00~08:00 滨州多普勒雷达速度图可见, 风暴中心附近的多个中气旋是直接造成地面冰雹大风的 β 中尺度系统风暴主体回波移动前方 20~50km 处不断有新对流单体 ( 中气旋 ) 被激发并快速发展为新的强风暴中心, 引起强风暴呈现不连续传播特征这种强雷暴的发展模式可能与该时段中低层次级环流的作用具有密切关系 23 日 08:00 后 ( 图 11b), 来自对流层中层的冷空气已减弱, 高低层的次级环流已不清楚, 对流云表现出低层辐合高层辐散的特征, 表明对流风暴已达到强盛阶段雷达观测此后的对流风暴前部无新生单体产生, 对流风暴云体外形趋于光滑规则强对流风暴的发展将主要依靠自身的动力热力作用, 包括潜热释放对流反馈等机制维持风暴的强盛, 并继续造成鲁中山区局地暴雨冰雹和大风灾害 5 结论和讨论 (1) 本次强风雹灾害天气是在西北气流为主导的环流形势下产生的, 在西北流型和凌晨时段发生如此强的对流风暴, 为历史罕见与典型西北流型华北强对流天气中层干冷的层结特征不同, 本过程对流层中层一直为湿冷西北气流控制 ; 低层增温, 中高层降温和整层增湿是造成层结不稳定的重要原因上游地区对流活动引发的干冷下沉气流沿太行山背风面下滑形成的下坡风是触发强风暴发生的直接动力机制 (2) 应用倾斜湿位涡发展理论分析风暴的发展机制表明, 风暴发生区具有明显的等熵面倾斜, 对流层中低层狆 m1<0 区域, 同时狆 m2>0, 满足倾斜涡度发展的条件, 对强风暴的发展具有一定指示性同时, 由于强对流风暴运动的非地转非静力平衡特征, 特别质量强迫的存在, 位涡守恒性被破坏, 使风暴低层的位涡发生异常, 表现出与大尺度天气系统不同的位涡特征 (3) 太行山地形对华北强对流风暴天气的作用表现为 : 一方面地形可引起边界层风场的变化, 包括风垂直切变下坡气流和中尺度辐合线, 从而对风暴的触发组织和移动发挥作用, 不同的下垫面地形可使风暴低层结构发生变化另一方面, 山脉背风波动和地形产生的不平衡气流, 有利于中尺度重力波的产生和维持, 并对下游地区强对流天气的发展产生影响

11 514 气象学报 63 卷参考文献 [1] 蔡则怡, 李鸿洲, 李焕安. 华北飑线系统的结构与演变特征. 大气科学,1988,12(2):191~199 CaiZeyi,Li Hongzhou,Li Huanan.Thecharacteristicsof structureandevolvementaboutsquallineinnorth China. ChineseJAtmosphericSci(inChinese),1988,12(2): [2] 丁一汇, 李鸿洲, 章名立等. 我国飑线发生条件的研究. 大气科学,1982,6(1):18~27 DingYihui,LiHongzhou,Zhang Mingli,etal.Researchfor thedevelopingconditionofthesquallineinchina.chinesej AtmosphericSci(inChinese),1982,6(1):18-27 [3] 李鸿洲, 蔡则怡, 徐元泰. 华北强飑线生成环境与地形作用的数值试验研究. 大气科学,1999,23(6):713~721 LiHongzhou,CaiZeyi,XuYuantai.Anumericalexperiment oftopographicefectongenesisofthesquallineinnorthchi na.chinesejatmosphericsci(inchinese),1999,23(6): [4] 朱官忠, 刘恭淑. 华北南部产生中尺度对流复合体的环境条件分析. 应用气象学报,1998,9(4):441~448 ZhuGuanzhong,LiuGongshu.Analysisoftheenvironmental conditionsforformingmesoscaleconvectivecomplexesinthe southernpartofnorth China.J Appli MeteorSci (in Chi nese),1998,9(4): [5] 朱民, 余志豪, 陆汉城. 中尺度地形背风波的作用及其应用. 气象学报,1999,57(6):705~714 Zhu Min,YuZhihao,LuHancheng.Theefectofmeso scale leewaveanditsapplication.actameteorsinica(inchinese), 1999,57(6): [6] 肖庆农, 伍荣生. 地形对于气流运动影响的数值研究. 气象学报,1995,53(1):38~49 Xiao Qingnong,Wu Rongsheng.Numericalexperimentsof influenceoftopographyonatmosphericmotion.actameteor Sinica(inChinese),1995,53(1):38-49 [7] 张群, 张维桓, 姜勇强. 边界层辐合线发展成飑线的数值试验. 气象科学,2001,21(3):308~315 ZhangQun,ZhangWeihuan,JiangYongqiang.Numericaltri alofpblconvergencedevelopingtosqualline.scientiame teorsinica(inchinese),2001,21(3): [8] 范广洲, 吕世华. 地形对华北地区夏季降水影响的数值模拟研究. 高原气象,1999,18(4):659~667 FanGuangzhou,LuShihua.Numericalsimulationstudyfor theefectofterrainonnorth Chinasummerprecipitation. Plateau Meteor(inChinese),1999,18(4): [9] 吴国雄, 蔡雅萍, 唐晓箐. 湿位涡和倾斜涡度发展. 气象学报, 1995,53(4):387~405 Wu Guoxiong,Cai Yaping,Tang Xiaojing.Moistpotential vorticityandslantwisevorticitydevelopment.acta Meteor Sinica(inChinese),1995,53(4): [10] 高守亭, 雷霆, 周玉淑. 强暴雨系统中湿位涡异常的诊断分析应用. 应用气象学报,2002,13(6):662~670 GaoShouting,LeiTing,Zhou Yushu.Diagnosticanalysisof moistpotentialvorticityanomalyintorrentialrainsystems.j AppliMeteorScience(inChinese),2002,13(6): [11] 崔晓鹏, 吴国雄, 高守亭. 西大西洋锋面气旋过程的数值模拟和等熵分析. 气象学报,2002,60(4):385~399 CuiXiaopeng,WuGuoxiong,GaoShouting.Numericalsimu lationandisentropicanalysis offrontalcyclones overthe western Atlanticocean.Acta MeteorSinica (in Chinese), 2002,60(4): [12] 李英, 段旭. 湿位涡在云南冰雹天气分析中的应用. 应用气象学报,2000,11(2):242~248 LiYing,DuanXu.Diagnosticanalysisofmoistpotentialvor ticityforhailinsouthernyunnan.japplimeteorologicalsci (inchinese),2000,11(2):242~248 [13] 曹钢锋, 张善君, 朱官忠等. 山东天气分析与预报. 北京 : 气象出版社, pp CaoGangfeng,ZhangShanjun,ZhuGuanzhong,etal.Analy sisandforecastofshandong Weather.Beijing:ChinaMete orologicalpress, pp [14] 楼小凤.MM5 模式的新显式云物理方案的建立和耦合及原微物理方案的对比分析 :[ 博士学位论文 ]. 北京 : 北京大学, 2002 LouXiaofeng.Developmentandimplementationofanewex plicit microphysical scheme and comparisions of original schemesofmm5:[doctorpaper].beijing:pekinguniversi ty,2002 [15] 陈泽宇, 吕达仁, 刘锦丽. 中纬度夏季一次深厚对流过程的数值模拟研究 : 高空热力层结和风切变影响. 大气科学,2002,26 (6):744~750 ChenZeyu,LuDaren,LiuJinli.Anumericalstudyofamid latitudedeepconvectioncase:influenceoftheuppertropo sphericstratificationandwindshear.chinesejournalofat mosphericsciences(inchinese),2002,26(6):

12 4 期龚佃利等 :2001 年 8 月 23 日华北强风暴动力机制的数值研究 515 犖犝犕犈犚犐犆犃犔犛犜犝犇犢犗犉犃犛犈犞犈犚犈犕犈犛犗犛犆犃犔犈犆犗犖犞犈犆犜犐犞犈犛犢犛犜犈犕犐犖犜犎犈犖犗犚犜犎犆犎犐犖犃犗犖 23 犃犝犌犝犛犜 2001 GongDianli ( 犛犺犪狀犱狅狀犵犕犲狋犲狅狉狅犾狅犵犻犮犪犾犐狀狊狋犻狋狌狋犲, 犑犻狀犪狀 ; 犆狅犾犲犵犲犘犺狔狊犻犮犪犾犪狀犱犈狀狏犻狉狅狀犿犲狀狋犪犾犗犮犲犪狀狅犵狉犪狆犺狔, 犗犮犲犪狀犝狀犻狏犲狉狊犻狋狔狅犳犆犺犻狀犪, 犙犻狀犵犱犪狅 ) WuZengmao FuGang ( 犆狅犾犲犵犲犘犺狔狊犻犮犪犾犪狀犱犈狀狏犻狉狅狀犿犲狀狋犪犾犗犮犲犪狀狅犵狉犪狆犺狔, 犗犮犲犪狀犝狀犻狏犲狉狊犻狋狔狅犳犆犺犻狀犪, 犙犻狀犵犱犪狅 ) 犃犫狊狋狉犪犮狋 Byusingalavailableobservationaldata,suchassatelite,Dopplerradarimagesaswelasroutineob servationaldata,andnon hydrostaticmodelmm5v3.5,aseveremcs(mesoscaleconvectivesystem)oc curredon23august2001inthenorthchinawasinvestigated.theevolutionaryprocessanddevelopment mechanismofthismcscasewasanalyzedwithemphasisontheorographicefectoftaihang Mountainon itsdevelopment.theobservationaldatashowthisseverestormoccurredinthepaternofhigh levelnorth westflow,thetemperatureincreaseinthelowertropospheredescendsinhigh midtroposphereandthehu midityenhancedinthewholetroposphereinducedtheconvectiveinstability. Thesimulationresultsindicatedthatthedownhilflowform Taihang Mountainisthekeymechanism toinitiatetheconvection.thecharacteristicsofmoistpotentialvorticity(mpv)showthatthemoistisen tropesslantisthefoundationoftheconvection,theverticalandhorizontalcomponents,i.e. 狆 m1and 狆 m2 respectively, 狆 m1<0and 狆 m2>0intheloweroftropospherearefavorableconditionforthedevelopmentof MSC,itcanserveasapowerfultoolsforthediagnosisandpredictionofMCS.Ontheotherhand,because theconditionofmpvconservationcanbedestroyedpossiblyinthemcsprocess,forexamplethehydro staticbalanceandthemassforcinginducedbytheprecipitationetal.,canleadto MPVanomaly.Inthe earlyperiodofthe MCS,theverticalsecondcirculationcanbeanimportant mechanism totriggernew thunderstorm andconducethestorm spreadingbyjump model.thetaihang Mountaincanchangethe windfieldinboundarylayer;includedownhilflow,verticalshearofwindandmesoscaleconvergenceline. Asaresult,itplaysanimportantroleininitiating,organizationandmovementofthestorm. 犓犲狔狑狅狉犱狊 :Severethunderstorm,Numericalsimulation,Dynamicalmechanism,Orographicefect, Slantwisevorticitydevelopment.

气象第卷第期对流不稳定度

气象第卷第期对流不稳定度张芳华张涛周庆亮杨克明国家气象中心北京中国气象科学研究院灾害天气重点实验室利用上海自动观测站资料的逐小时分析资料水汽云图和多普勒雷达探测等资料对年月日上海飑线天气过程进行天气动力学诊断分析结果表明上海处在上干冷下暖湿的不稳定大气层结中对流层上层干冷空气的侵入和边界层暖湿气流强烈辐合是飑线天气的触发机制飑线以断续线型与后续线型相结合的型式形成在对流单体的右侧不断新生单体

第 63 卷第 4 期 2005 年 8 月 气象学报 ACTA METEOROLOGICASINICA Vol.63,No.4 August 年 8 月 23 日华北强风暴 动力机制的数值研究 龚佃利 ( 山东省气象科学研究所, 济南,250031; 中国海洋大学海洋环境学院,

第 63 卷第 4 期 2005 年 8 月气象学报 ACTA METEOROLOGICASINICA Vol.63,No.4 August 年 8 月 23 日华北强风暴动力机制的数值研究龚佃利 ( 山东省气象科学研究所, 济南,250031; 中国海洋大学海洋环境学院,