三维云模式解析课件

1、第三章 对流云数值模拟3.2 三维云模式概述一、国内著名三维云模式二、大气所三维冰雹云模式3.1 对流云的宏观特征 空间尺度、生命期、气流、温度、 含水量3.3 大气所三维雹云数值模式的 动力学框架 一、模式控制方程组二、湍流参数化次网格尺度混合项D三、边界条件四、初始条件3.4 大气所三维雹云数值模式的 微物理过程参数化方案（略）3.5 该模式程序结构和运行方法3.6 实习内容和安排3.1 对流云的宏观特征一.对流云的空间尺度对流云高度因纬度, 气团属性,发展阶段而异一般而言,处于相同发展阶段的对流云,中低纬高,高纬低;海洋性气团中的高于大陆性气团中的1. 高度3.1 对流云的宏观特征3.

2、垂直厚度H与水平宽度L具有同一量级2. 水平尺度一般, 在淡积云阶段为数百米到1km量级;气团性浓积云或积雨云单体为几km至十几km二. 对流云的生命期不同尺度的对流云,生命期是不一样的.同一尺度的对流云,其生命期也因地有很大差异.根据我国九江地区的观测结果,频数中值生命期,无降水的浓积云为12min,降水性浓积云为25min,雷雨云为65min海洋性晴天积云同一高度多次穿云飞行结果三. 对流云中的垂直气流小尺度湍流: 水平尺度L小于100-200米(在一些积雨云研究中定为400-500米)的垂直运动系统性垂直气流: 将L大于上述临界尺度的垂直运动称为系统性垂直气流.分类:飓风系统积雨云中上升

3、和下沉气流合计的水平尺度频数分布图众数直径为1-2km1. 系统性垂直气流的尺度热带浓积云中垂直气流的水平尺度可达1-2 km. 积雨云中上升气流的平均”直径”为1.5 km,下沉气流的平均”直径”为1.2 km2. 系统性垂直气流区的相对面积发展中的晴天积云的中下部,上升气流占据云区的70%-80%,只在边缘部位有下沉气流发展旺盛的浓积云中,上升气流区的体积约占云体积的40%.在各类对流云的上部,上升气流速度较小,占据的范围不超过云体积的30-40%,其余部分是高湍流性的3.垂直气流的大小澳大利亚东海岸活跃和成熟期积云垂直气流的平均值一般为几米/秒,最大垂直速度可达30-40米/秒, 甚至也

4、出现过63米/秒 的最大上升气流速度了解对流云中的垂直气流的空间分布和时间演变特征很重要4. 垂直气流随高度的分布一般来说,上升速度随高度的分布曲线近似为抛物线型上升速度随高度的分布(苏联)但不同部位上升气流的垂直分布有区别5. 垂直气流的水平分布浓积云中飞机穿云观测事实上,对流云中垂直气流的分布非常复杂,与云发展阶段具有密切联系.三. 对流云中的湍流湍流可伸展到对流云外,当然,云内湍流还是强于云外.当环境盛行风存在垂直切变时,对流云周围的湍流分布就不再对称,而是顺着风向向下延伸,但湍流强度却是迎风面的大.四. 对流云中的温度一般来说, 温度的正距平和上升运动配置在一起,负距平则和下沉运动相联

5、系.1-2 4(甚至6-10) 0-1 正距平随对流发展阶段的演变:一般情况对流云中温度的分布存在着起伏水平分布垂直分布对流云内温度的垂直递减率比湿绝热递减率要大,它与云外的空气温度递减率相近五. 对流云的含水量对流云中的含水量随时间,地点和云的发展阶段而异1. 含水量数值与云底温度和环境空气湿度以及云内外空气的混合程度有关.苏联,雷雨云中一般为4-7g.m-3, 最大值达20g. m-3 美国,最大值也达20g. m-3 2. 含水量影响因子云底温度云顶部位最大含水量湖南地区浓积云飞机观测结果(1962-1963)3. 含水量的垂直分布极大值位于云的中上部4. 含水量的水平分布含水量高值区往

6、往与升速较大的云泡相对应Warner概括出”大礼帽”式的水平分布形式.含水量在云边急剧上升,在云体内分布平缓在不同高度上,含水量的水平起伏不一样, 云的下部变化小,而云的中上部变化大值得一提的是对流云中含水量的空间分布是十分复杂的.发展旺盛的对流云中常有多个极大值中心,最大值中心未必位于云中心,有时在云边缘含水量可急剧增加到很大.七. 对流云宏观特征小结垂直厚度H与水平宽度L具有同一量级 尺度从几百米到十几km2. 生命期几min到几h3. 气流 包含有系统性上升气流区,系统性上升气流区,湍流区. 上升气流区和下沉气流区的水平尺度众数为1-2km 垂直气流的平均值一般为几米/秒,最大垂直速度可

7、达30-40米/秒, 甚至也出现过63米/秒 的最大上升气流速度一般来说,上升速度随高度的分布曲线近似为抛物线型积云中的湍流交换系数达102m2/s一般来说, 温度的正距平和上升运动配置在一起,负距平则和下沉运动相联系.正,负距平可达甚至超过44. 温度对流云内温度的垂直递减率比湿绝热递减率要大,它与云外的空气温度递减率相近5. 含水量平均值约为几g.m-3, 最大可超过10 g.m-33.2 概述云数值模拟最早可追溯到五十年代中期, 建立了各种一维、二维云模式。直到70年代以后随着计算机性能的提高，三维云模式才得以建立。一、国内著名三维云模式:1、中科院大气所孔凡铀1990：完全弹性三维冷云

8、单参数模式洪延超1998：三维弹性冰雹云催化模式 在以上模式基础上，对云微物理过程的处理进行了改进，并结合进催化剂的核化过程，可以用来研究催化防雹技术。郭学良2001：三维冰雹分档强对流云模式2、中国气象科学研究院刘玉宝1992：双向套网格、变网格距准弹性中-大气模式胡志晋2000：三维层状云模式胡志晋2000：三维对流云模式3、总参大气所许焕斌等1990：三维冰雹云模式这门课里主要以大气所的三维冰雹云模式为例， 学习三维云模式的运转和结果分析二、大气所三维冰雹云模式 1998年洪延超研究员保留孔的对流云模式中的动力学框架（包括初始条件、边界条件、启动云发展方法在内），对云微物理过程参数化方案

9、进行了比较大的改进，主要是对水物质的类别分得更细（5类-8类），数浓度谱函数假定为GAMMA函数（孔假定为M-P分布函数），采用双参数方案（孔的为单参数方案）。 据不完全统计，到目前为止，利用该模式发表研究论文11篇，而且一直有人对该模式进行改进当中，除了大气所还有南大都有人在不断改进该模式。改进主要体现在：1、对流启动方式；2、参数化方案；3、与化学模式耦合3.3 三维雹云数值模式的 动力学框架一、模式控制方程组1、假设1）大气为无粘性可压缩流体2）忽略地转偏向力3）任一时刻预报变量场A可以分解为：4）云滴完全随气流运动，其他水成物质粒均具有相对空气的沉降运动。 将Navier-Stokes

10、运动方程对网格元进行雷诺平均，并运用对流近似处理，再用无量纲气压扰动量和位温分别代替气压P和温度T，即可导出模式运动方程组。此外，还有位温预报方程、质量连续方程、水物质预报方程。2、控制方程qt 是液、固态水总比质量其中，M= qx,nx，下标 x= v,c,r,i,s,g,f,h,a, 8类水物质变量和催化剂，除水汽和云滴只用比质量外，其他都包含比质量和比数浓度两类变量。另外，上式中右边第一项为粒子降落项。上式中，共20个预报方程，共20个预报量二、湍流参数化次网格尺度混合项D有如下表达式：为了闭合模式方程组，必须用网格元平均场量来表示上式中的次网格湍流通量。该模式采用一阶湍流粘性闭合参数化

11、方案，并考虑了大气稳定度的贡献以及水成物拖曳影响，湍流脉动相关量可以表示为其中k为湍流交换系数，该模式分热量和动量两种，具体计算方法略。三、边界条件1、侧边界条件 雹云模式属于有限区域模式，必须人为建立模拟域的侧边界和给出相应的侧边界条件。较好的侧边界条件应当是能够允许模拟域内外的各种波动和平均流自由穿过侧边界。就小尺度雹云模式而言，最主要的是要让区域内对流活动等激发的各种扰动通过侧边界向外传出，而只产生最低限度的反射。本模式对侧边界的法向速度un采用如下辐射边界条件上式中n为侧边界法向，C*相当于域内优势重力内波指向域外的水平固有相速。该模式取30m/s.2、上下边界条件模式的上、下边界首先

12、限定ZT为模拟域顶高。只要模式上界取得距云顶有足够距离，以上限定还是合理的。在上下边界以内z/2处，方程中的垂直湍流混合项取零值，这个条件可使初始平均状态廓线在边界附近的垂直梯度不受到湍流混合作用的歪曲。四、初始条件1、模式方程组初始条件在初始时刻，令（不包括比湿）是模式水平均匀初始环境场由探空资料确定，并在整个模拟时段内保持不变。是位温初始扰动场，用以启动对流，其具体形式与对流启动方式有关。2、对流启动方式热泡、湿泡、湿热泡低层冷空气出流辐合四种启动方式设 (,)是扰动场中心坐标, (r,r,r)分别是扰动区在三个坐标方向的半径,=(,)是中心最大扰动位温值,、0分别为饱和混合比和初始水汽混

13、合比。1）热泡启动方式湿热泡内位温偏差由下式确定:2) 湿泡启动方式 在湿泡内，水汽混合比由下式给出:3.4 大气所三维雹云数值模式的 微物理过程参数化方案参见下列文献3) 湿热泡启动方式洪延超, 积层混合云的数值模拟研究() 模式及其微物理过程参数化,气象学报,1996, 54,544-557洪延超,三维冰雹云催化数值模式,气象学报,1998, 56（6）, 641-6533.5 该模式程序结构和运行方法FM33wc.F90 二、运行方法第一步：准备好探空资料（indata,wind），设 定格点数（domain.DEF） 对流扰动等参数（model.input）第二步：把include文件

14、、主程序和所有有关子程序编译、链接、运行第三步：处理输出结果dname所对应的文件Tiwi岛(130.0-131.55E，11.16-11.95S )由两块相对平坦的小岛组成，海拔最高为120米。Tiwi岛东西长约150km，南北约50km。岛上大多数区域覆盖着桉树林和红树林。 湿季里，对流频繁 3.6 所模拟强风暴的大体情况2.2 模式控制试验强雷暴在地方时14:20达到最强,中心54-58dBz,右边的单体很快消散11:30 , 出现初始回波见CooK(2004)17时对流云主体消失卷云持续到19时以后2003年11月23日卷云伸展18时的红外卫星云图.伸展约300km图2 对流云主体移动

15、方向和云砧伸展方向Fig.2 Directions of the storm motion and cirrus outflowAustralia图1 2003年11月23日飞机观测飞行路线图Fig. 1 Location of Tiwi Islands and the flight tracks on 23rd November 2002 2.2 模式控制试验中科院大气所三维完全弹性强风暴模式 (IAP-CSM3D)包含了水汽、云水、雨水、冰晶、雪、霰、冻滴和冰雹等8种水物质的微物理过程, 各类粒子的数浓度分布密度函数用Gamma函数表达 双参数，既预报水质粒的比质量，也预报比数浓度 本文修

16、改了初始扰动方案 雷达反射率因子的计算方法2.2 模式控制试验-模拟方案设计x= y=1.0km, z=0.5km, t=5s, t=1sDomain: 200km200km22kmT=360min最大云滴数浓度 (unit:cm-3) Ctl 600 (M.O.Andreae等, science,2004,亚马逊地区)2.2 模式控制试验模拟总的情况模拟的第250分钟总含水量Y-Z 垂直剖面图(单位：g/m3) (units：g/m3) .等值线分别代表 0.001, 0.01, 0.1 and 0. 5 g/m3250min对流结束云体开始演变为层状云模拟到90分钟时, 已经形成明显的云砧. 3.7 实习目的、内容一、目的1、大致了解大气所三维强风暴模式（也叫三维冰雹云模式）2、利用雷达、卫星图像资料分析强风暴演变过程、云砧伸展情况及高度范围、并利用探空资料解释强风暴的发生及其特征与环境条件间的联系（例如，静力稳定度与强对流的产生、低空风与雷暴移动方向、高空风与云砧伸展方向）3、模拟结果与实测对比 4、对流云主体宏微观特征模拟二、内容1、分析2002年11月23日发生在澳大利亚Tiwi岛上的一次强风暴的宏观演变过程及环境