wps大气运动的稳定性理论.pptVIP

两层模式中把大气分成由编号为0, 2,和4的面隔开的两层，如下图所示。 将涡度方程(8.18)写在面1、面3两个中间层上，但要先估计散度项?ω/?p在每一面上的值。将有限差分(finite difference)近似用与垂直导数项可得： Δp为0~2和2~4面之间的气压差，下标表示在(0~4)某一面上。 1、3面上的涡度方程分别为 (8.20) (8.21) 为简单起见，已假定ω4=0，这对于水平地表是近似成立的。 再将热力学能量方程(8.19)写在面2上，用以下差分形式估算?ψ/?p： 能量方程变成 (8.22) (8.22)右边第一项是500hPa风场对250-750hPa厚度场的平流。然而，在此模式中，500hPa的流函数ψ2不是预报值。 所以ψ2必须由250和750hPa之间的线性内插而得到： (8.23) 各变量的排列如图9.1所示。假定流函数ψ1和ψ3由只依赖于y的基本部分和随 x和t变化的扰动部分组成，即 (8.24) 在面1和面3上的纬向风为U1和U3，于是扰动场只与经向和垂直速度场有关。把(8.24)代入(8.20)~(8.22)再线性化，便得扰动方程组： (8.25) (8.26) (8.27) 这里使用了β-平面近似，V2由ψ1和ψ3线性内插得到。 令 , , (8.28) 代入(8.25)~(8.27)中得到关于振幅A，B和C的线性代数方程组： (8.29) (8.30) (8.31) 在系数行列式为零时，才有非平凡解(nontrivial solutions)。因此，相速c满足 可解出c的二次方程： (8.32) 推导中已经取 ，长度尺度λ-1称作变形半径(radius of deformation)。 或从(8.29)~(8.31)消去A、B或C任意两个也可得到(8.32)。(8.32)关于相速的解为 (8.33) 其中 Um和UT分别为垂直平均纬向风和间隔2Δp的背景热成风，即 我们已经证明了只有相速度满足(8.33)时，(8.28)是 (8.25)~(8.27)的一个解，而方程(8.33)的形式非常复杂。 如果注意到当δ0时，相速度存在虚部，扰动将指数增长。在详细讨论(8.33)式的一般性质之前，先考虑两个特例。 一是基流为正压，UT=0，背景流场没有热成风，此时相速度为： (8.34) (8.35) 它们都是方程(8.32)的解。正压基流的两层模式自由扰动有对应的实数值。相速c1就是一维正压Rossby波的频散关系。把(8.34)代入(8.29)~(8.31)可求出A=B，C=0，所以扰动是正压的。另一方面，(8.35)可以解释为内斜压Rossby波的相速度。注意，c2是具有自由面的均值海洋中类似于Rossby波速的频散关系。但在两层模式中分母中用2λ2代替了海洋中的f02/gH。两种情形中都有垂直运动伴随着Rossby波，所以静力稳定性对波速有影响。可以证明，把c2代入(8.29)~(8.31)可看出ψ1和ψ3场的位相相差180o，所以扰动是斜压的，虽然基态是正压的。同时，ω2场超前250 hPa位势场90o(于是最大上升运动出现在250 hPa槽的西侧)。 垂直运动这种型可以这样理解，因c2–Um0，所以扰动型相对于平均气流是向西移动的。在一个随平均风移动的坐标系中，涡度的变化仅由行星涡度平流和辐散项引起，同时，厚度的变化必须唯一地依赖于垂直运动引起的绝热加热或冷却。因此，为使系统的西