次超级单体风暴多普勒天气雷达资料-气象科学.PDF

维普资讯 第 27卷 第 2期 气 象 科 学 Vo1．27。No．2 2007年 4月 SCIENTIA M ETEOROLOGICA SINICA Apr．。2007 一 次超级单体风暴多普勒天气雷达资料 的低层风场反演研究 张少波 陆汉城 康建伟 (解放军理工大学气象学院大气科学系 ，南京211101) 摘 要 本文使用二维变分风场反演方法对浙江 台州 2003年 4月 12日一次超级单体风暴 的多普勒天气雷达探测资料进行低层风场反演研究 。二维变分方法 以均匀风场反演结果作为初 猜场 ，在 目标泛函中加入散度和涡度约束作为正则项来降低不适定性的影响。反演得到的低层水 平风场 比较好 的符合超级单体风暴 的组织化 的气流结构和强 的风切变 ，通过对散度场垂直积分获 得垂直速度 ，上升气流主要集 中在强回波区域 ，下沉气流一个位于风暴前进 的左侧 的弱回波区，一 个位于钩状 回波的内侧 。 关键词 多普勒天气雷达 超级单体风暴 二维变分 分类号 P415．2 文献标识码 A 引 言 多普勒天气雷达探测信息及其所提供的丰富的物理量产品在危险天气监测和预警方面发挥着重要的作 用。虽然如此 ，目前对多普勒雷达资料的使用还远远没有发挥它应有的作用。怎样在数值模式 中使用丰富 的雷达资料来提高预报水平 ，这是很多气象工作者 目前所关心的问题 ，在这个 问题 中主要面临两个主要 困 难 ，一个是多普勒雷达资料 中风场信息的提取 ，一个是雷达资料的同化。 使用多普勒雷达资料进行风场反演 的研究从上世纪 6o年代就已经开始，从最早 的VAD(Lhermitte、 Atlas1961)、VARD(Easterbrook1975)、VVP(Waldteufel、Corbin1979)到 UW (Persson1987)、VAP(陶祖 钰 1992) 、VPP(郎需兴 2001) 等方法都是通过在某 区域内风场呈均匀或者线性分布的假设条件下直接 从径 向速度来反演风场 。VAD和EVAD方法都是反演整层 的平均风速风 向；VARD、Uw 和 VAP方法在 某-d,的区域 内假定风场均匀或线性分布，但它们都对径 向速度资料误差非常敏感 ；VVP方法在一空间体 积 内假定风场呈线性分布 ，可以反演包括三维风场 、散度 、涡度在内的多个参数，但是该方法在计算中遇到的 病态矩阵求解问题一直没有得到很好的解决 3『]。Rinehart(1979)、Tuttle和 Foote(1989)提 出的示踪相关反 演方法 ，还有胡明宝和汤达章(1998)]提出的矩不变量反演方法对强度场进行跟踪 ，进而获得风场信息的方 法 ，这种方法虽然能够在低仰角得到 比较好的二维风场 ，但是方法本身并没有考虑追踪过程 中云物理过程对 强度场的改变。此外 ，也有结合动力学方程进行反演如涡度 一散度方法 (姜海燕 、葛润生 1997)]，这种方法 使用涡度 一散度方程进行反演 ，而一些不太强 的天气过程中的涡度散度往往会被夸大 。因此 ，这些方法大多 数都只能适用于某一类的风场反演问题 ，对其它情况下的风场反演则存在局限，但是这些方法为真实风场的 反演提供了较多的可选择项。 随着变分方法的发展和计算机技术的进步 ，变分方法被 引入到风场反演 中来。最初 Liouetal(1991) 提出完全模式逼近 (Nudging)方法 ，使用完全数值模式 (FullNumericalMode1)通过数值迭代反演风场和温 度场等参数 ；随后 Suneta1．(1991)利用共轭方程理论和最优控制技术 ，通过伴随方法反演风场和温度场等 参数。这两种方都受到模式计算量大 、运算要求高 的限制 ，为 了解决这些 困难 ，Qiueta1．、xuetal等 收稿 日期 ：2005-10—10；修改稿 日期 ：2006—03—07 基金项 目：本文由国家重点基础研究专项经费(编号：2004CB418300)资助