P-σ坐标系数值模式中气压梯度力的递推算法

在静力关系成立的条件下,推导出气压梯度力的递推算法,用以计算有地形P-σ坐标系模式中的气压梯度力.将递推算法与对应的6种直接计算方案相比后发现,在有气压梯度力的参考大气中,递推算法的计算精度相对于直接计算方案的精度有较大改进,特别是采用递推后的两种经典修正格式和Corby格式,能较大提高水平气压梯度力的计算精度,从而改善数值模式的预报效果和稳定性.     

有地形模式中气压梯度力误差扣除法

：在现有有地形的数值模式中，陡峭地形区气压梯度力的计算存在一个普遍问题，即计算精度较高的格式比较繁琐而费时，简单省时的格式又不精确和不稳定。为解决这个问题，作者等在最近提出了一种方法，称为气压梯度力的误差扣除法。该方法假定：气压梯度力的计算误差主要来自地形而与气压形势关系较小。用理想场对该方法进行检验后表明，这一方法是成功的。本文则用气候模式对作者提出的有地形数值模式中气压梯度力的误差扣除法进行了模拟检验。所用模式是作者等使用多年的Ｐ－混合坐标系５层模式，选用了四种气压梯度力的计算格式，即ＤＤＤ格式、Ｃｏｒｂｙ格式、平均温度格式和经典中央差格式。比较了这四种格式在有无误差扣除时的模拟结果，发现：对于计算精度较高的格式，如ＤＤＤ格式、Ｃｏｒｂｙ格式及平均温度格式，有无误差扣除的结果相差不大，但误差扣除法仍可在一定程度上改善模拟效果。对于计算精度差的格式，如经典中央差格式，在无误差扣除时计算不稳定，得不到模拟结果，进行误差扣除后，从根本上提高了其计算精度，因而也提高了计算的稳定性，达到了较满意的模拟效果。而且与其它格式的模拟结果相当接近。本文提出的误差扣除法可同时用于格点模式和谱模式。     

误差扣除法在云南中尺度数值模式中对气压梯度力的计算试验

使用云南中尺度数值预报试验模式^[1]，引入实际地形，选用静力扣除格式、回插格式、局地等温格式、经典修正格式以及模式中原计算格式（Corby格式）等5种计算格式计算气压梯度力，对气压梯度力的误差和扣除法进行了模拟试验。比较了这5种格式在有、无误差扣除时的模拟结果之后，发现误差扣除法与上述格式结合使用时，均能够不同程度地提高预报效果，尤其是在与经典修正格式结合使用以后，预报效果进一步改善，得到了较满意的模拟效果。     

数值模式中气压梯度力的算法试验

本文在差微差一致性二阶精度的坐标变换公式的基础上,推导出了气压梯度力计算的新公式。同时,根据静力扣除法原理,提出了初值静力扣除法,用以计算有地形数值模式中的气压梯度力。将这两种新的方法与二阶计算精度的回插法和局地等温法(Corby格式)相比后发现:新方法(文中称为经典修正法)与回插法精度相近,但可节省计算时间,比局地等温法精度高。静力扣除法则效果较差,但其原因可能在于地转风垂直插值方法有误差。本文的试验模式是上海台风研究所的业务数值预报模式,文中对模式及若干技术处理也作了简要介绍。     

温度初值确定之合理性与气压梯度力项的计算

本文讨论数值试验中两个基本问题:温度初值确定和气压梯度力的计算。首先比较几种利用静力方程差分解求温度初值方案所引起的误差,提出一种把等压面高度场插值与静力方程微分解结合起来求取等σ面温度初值的方法,结果表明用这种方法确定的温度初值及其水平分布和垂直递减率接近实测值。文章又讨论山脉地区气压梯度力计算问题,表明气压梯度力的计算精度不仅与计算方案有关,而且与温度初值有很大的关系。利用本文提出的方案得到的温度初值来计算气压梯度力,误差达到比较满意的精度。     

有地形模式中气压梯度力误差和扣除法的模拟试验

在现有有地形的数值模式中，陡峭地形区气压梯度力的计算存在一个普遍问题，即计算精度较高的格式比较繁琐而费时，简单省时的格式又精确和不稳定，为解决这个问题，作者等在最近提出了一种方法，称为气压梯度力的误差扣除法，该方法假定：气压梯度力的计算误差主要来自地形而与气压形势关系较小，用理想场对该方法进行检验后表明，这一方法是成功的。     

一个六层亚洲有限区域模式及对一次西南涡过程的数值模拟

本文介绍的六层亚洲有限区域模式包含了地形作用和地面摩擦、动量、热量和水汽的水平和垂直方向扩散过程（其中包括了下垫面感热作用）、水汽凝结潜热、云和辐射作用以及积云对流等物理过程。并运用新的有效的计算气压梯度力的方法，为减少青藏高原引起的σ面上的初值（由p面到σ面）的插值误差，用了初始时刻的迭代插值方法。另外，利用上述模式，以1981年7月11日12时客观分析为初始场，对青藏高原背风坡地区的西南涡形成过程进行了数值模拟试验，较成功地研究了青藏高原动力作用对西南涡形成和发展过程的影响，模式还较成功地模拟出在西南涡形成过程中在我国西南地区低层西南风急流的迅速形成和维持过程。如果降低模式地形高度（如将青藏高原降至3000 m以下），地形的坡度也变小了，西南涡的模拟就不能成功。本模式所用的新的气压梯度力的计算方法，使预报效果得到了改进。     

复杂地形强降雪过程中垂直运动诊断分析

利用WRF模式对2018年11月30日伊犁河谷和天山北坡强降雪过程进行数值模拟,并分析复杂地形强降雪过程垂直速度和垂直动能变化机制。研究表明,冷锋过境造成地表气压升高,干空气气柱质量增大,从而导致垂直气压梯度力和干空气气柱浮力发生变化,进而引起垂直运动发生发展。垂直速度局地时间变化主要取决于扰动垂直气压梯度力、水物质拖曳力和扰动干空气浮力。在天山北坡,气流过山时,迎风坡的扰动垂直气压梯度力较大,扰动干空气浮力较小,二者合力促进上升运动;在背风坡,扰动垂直气压梯度力和扰动空气浮力形成向下的合力,产生下沉加速度,导致背风坡下沉大风。扰动垂直气压梯度力做功和扰动干空气浮力做功情况基本相反,背风坡扰动垂直气压梯度力和综合强迫做功项抑制垂直动能,扰动干空气浮力和水物质拖曳力做功项增强垂直动能。此外,扰动垂直气压梯度力和扰动干空气浮力做功项主要出现在中低层,水物质拖曳力做功项主要位于低层,平缓地形处的综合强迫做功明显小于地形复杂处。     

一个用阶梯山脉坐标表示的大气环流模式：模式的描述及其对中期预报的有效性

阶梯山脉或η垂直坐标已成为消除在沿陡坡面计算气压梯度力时出现的数值误差而建议的一种解决方法。这篇文章的主要目的是描述采用η坐标而设计的全球大气环流模式的发展，并检验该模式对中期预报的能力。首先，给出极地边界和极地滤波的处理。为了检验极地处理，给出了用浅水波方程得到的数值结果；其次，还把各种物理过程参数化引入多层η坐标模式。为了证实该模式有效，给出了几个１月份个例的数值积分。η坐标所列方程与随地形起伏的σ坐标所列方程很相似，允许用其中一种垂直坐标试验模式。于是，在同样的物理过程参数化下，我们用η坐标大气环流模式和σ坐标的大气环流模式进行了比较。另外，还与σ坐标的谱模式作了比较。还就模式的有效性，用几种不同的水平分辨率，４种观测的初始状态作了１０天积分。对相对低的分辨率的模式，预报结果稍微有利于谱模式和σ坐标模式。然而，采用较高分辨率，η坐标模式预报技巧评分与σ坐标模式结果不再有明显的差别。给出的其它结果还证实了η坐标在陡峭地形附近的优势以及与地表边界层处理有关的η坐标的潜在缺陷。     

不同垂直坐标系对垂直速度计算的影响

近年来，随着高精度中尺度数值预报模式和多尺度统一模式的发展，与之相关的一些非静力中尺度模式技术也随之提了出来，垂直坐标系的影响问题就是其中之一。本文借助于高精度美国新一代数值预报模式WRF(Weather Reseach and Forecast)，模拟分析了两种不同的坐标系对模式大气中垂直速度计算的影响。试验结果表明，不同的坐标系对大气垂直速度大小计算存在着差异，特别是当分辨率提高时，不同垂直坐标系对大气垂直结构的刻画的差别更加明显。由于垂直速度与降雨过程密切相关，垂直速度的计算差异必然影响模式对降雨过程的描述。理想试验结果表明，在气压地形追随坐标中气压梯度力的计算对地表气压计算很敏感。