基于随机扰动的积云团动力学分析.docxVIP

基于随机扰动的积云团动力学分析 1 种方法的对比研究 云对称与大环境的相互作用已成为近年来的研究主题之一。在大气环模型和中尺度值监测模式中，云对称与大环境之间的关系尚未得到很好的解决。目前，云参数化参数化方案的检测有三种方法。第一种方法是在大体积模型中测试此方案，并将该模型的预测量与观察到的未来时间的气结构进行比较。这种完全预测方法是最有效的方法，但由于云量散列误差（如空间差分和时间差异的段性误差、初始值化、云对称等物理过程的参数化等），因此云参数化的评估方法具有一定的局限性。 第二种方法是用半预报检验.半预报检验是利用积云参数化方案由大尺度温湿平流、地面通量和辐射增热的观测值作单步预报(积分一步)积云的属性,包括连续观测的降水以及视热源Q1和视湿汇Q2.半预报检验比完全预报检验的优越性在于它把由于参数引起的误差孤立出来,而不会与大尺度模式的其他物理过程相混. 由于积云对流的总体效应不能直接测量很难得到有积云对流的完整资料而积云团数值模式能产生合适的完整的资料,因而可以用来检验积云参数化方案,只不过所用的资料是模拟的观测资料.因此,第三种方法本质上也属于半预报检验.积云团模式能产生非常高分辨率的资料,这样的资料能给积云参数化方案非常详细的检验.虽然模拟的资料依赖于模式,但这方面的工作是令人鼓舞的. 本研究采用第三种方法来进行.首先设计一个二维积云团模式,再用模拟的资料来检查积云团活动对大气层结的影响. 已有学者进行了积云团数值模拟的研究.为了真实地模拟实际积云团的演变和个别积云的行为,模式必须是三维的.但是积云团属于中尺度系统,水平尺度为50～500km,时间尺度为十几小时.如果我们用1 km的水平网格距,10 s的时间步长来模拟三维积云团的演变过程,以现有的计算机能力难以实现的.因此,虽然二维积云团模式不能真实地模拟一些长寿命积云团和超级单体的结构和演变特征,但对于了解积云团的宏观统计特征还是有益的.Lipps和Hemler,Tao等的数值模拟结果表明,三维模式的平均统计结果与二维模式的基本相同.一些观测结果也表明很多中尺度对流系统是二维性质的.特别是近年来,Xu等的积云团数值模式也采用了二维模式. 2 大尺度升降速度的影响 现有积云团模式是从笛卡儿坐标系下的二维平面对称的完全弹性冰雹云模式发展而成的.从模式的基本框架上讲,与以往的单块积云模式基本类似.不同的是,积云团模式允许模式区域内几块对流云同时存在和发展,而且在模式中还包括了水平均一而垂直变化的大尺度抬升速度,就好象模式是嵌套在大尺度数值模式中一样.大尺度抬升速度在对流层中层的量级为10 cm·s-1,它的影响仅仅是在位温方程和水汽守恒方程中体现出来,而在所有的其他方程中不考虑大尺度抬升速度的影响.这样的近似与Krueger近似是一致的. 由于考虑了大尺度抬升速度的影响,热力学方程和水汽守恒方程为: 其中Sθ和Sqv为微物理过程调整项,Dθ和Dqv为次网格尺度湍流通量项. 除了热力学方程和水汽守恒方程,模式的动力框架、微物理过程参数化、次网格尺度湍流参数化、空间网格分布、数值方法和边界条件等均与文献一样. 积云团的启动采用冷池和随机扰动相结合的方法.加入到位温场的冷池随水平位置X的变化写为: 其中为最大冷却率. 随机扰动是每5间隔在模式低层(行星边界层)某一高度上所有的网格点加上位温场的随机扰动,其目的是触发生成积云团.具体的做法是在模式中把一系列均匀分布的范围为-0.5 K到0.5 K的随机数,加在所在高度上位温场的所有水平网格点上.每次加入的随机数系列都是不同的. 3 试验wfns的冷池参数 为了得到足够的关于积云团的宏观特征,共进行了11个试验,其中试验WEFNOS作为控制试验. 试验WEFNOS中,初始温湿层结采用GATE第Ⅲ阶段平均探空资料(图1).大尺度抬升速度采用的是Marshall岛上的典型天气扰动垂直速度.模式区域在水平方向取120个网格点,水平网格距为1 km,垂直方向取37个网格点,垂直网格距为0.5 km,这样相应的水平范围为120 km,垂直高度为18.5 km.由于模式是可压缩系统,采用时间分裂算法,小时间步长为2 s,大时间步长为10 s.实际积分时间取为8 h.试验WEFNOS的冷池启动是从积分开始第5 min后加入到第35 min结束.冷池的参数为X1=65 km,X2=85 km,也即冷池的宽度为20 km,最大冷却率为每秒0.015 K,冷池的高度为从250m到2250 m,而且强度不随高度变化.随机扰动所加入的高度为250 m. 其余10个试验的情况如下: (1)试验WEFNOS1,同试验WEFNOS,但随机扰动的随机数系列不同. (2)试验WEFHAS,同试验WEFNOS,但有初始水平风切变.初始水平风切变廓线采用GATE第