适用于非静力大气模式的完全平衡多矩约束有限体积方法

大气数值模式离散垂直动量方程时,一般而言气压梯度力和重力不易维持严格的静力平衡关系。为精确平衡数值离散的垂直气压梯度力和重力,基于高精度多矩约束有限体积方法引入完全平衡数值公式,即以满足静力平衡关系的热力学参考态对重力源项进行数值离散构造,发展了适用于非静力大气的完全平衡多矩约束有限体积方法。一维标准数值试验表明,完全平衡多矩约束有限体积方法能在较粗糙的计算网格点上保持静力平衡参考态的数值计算误差在计算机的单精度（10-6）和双精度（10-14）水平,在具有小量级扰动的初始条件下,完全平衡多矩约束有限体积方法能较好地模拟扰动的传播,二维非静力热泡试验进一步验证了完全平衡多矩约束有限体积方法对非静力大气运动的模拟能力。数值试验结果验证了所发展方法的完全平衡属性和适用性,这为非静力大气模式发展提供了良好参考价值。     

Implementation of CMA-GFS 3D-Var System on the Yin-Yang Grid

The global three-dimensional variational (3D-Var) data assimilation is implemented on a new quasi-uniform overset (Yin-Yang) grid on sphere. As a quasi-uniform spherical grid, it covers the sphere by overlapping two perpendicularly oriented grid components which is nothing but low latitude region of the usual latitude-longitude grid. Based on this characteristic of the Yin-Yang grid, it enables us to implement the regional 3D-Var system efficiently and accurately on the Yin or Yang component grid, respectively. The global analysis could update directly from the regional analysis since they have the same configurations like the precondition of eigenvalue decomposition for vertical direction, recursive filtering for horizontal direction, minimization method and observation operator and so on. However, the balance equation and vector wind are needed to be paid more attention on the Yin grid analysis due to its coordinate transformation. How to spread the observation information near the boundary of Yin and Yang grid is a key to the 3D- Var analysis. Extending double the horizontal correlation length distance in the overset boundary of Yin and Yang grid has successfully solved the problem. The results show that the analysis on the Yin-Yang grid is reasonable and similar to the result on the latitude-longitude (LAT-LON) grid. This paper provides a promising strategy for the development of a 3D-Var global system for overset grids.     

Comparison and Analysis of Two Different Calculation Schemes of SLEVE- Hybrid Coordinate in GRAPES_MESO Model

The calculation scheme of the smoothed-level and hybrid (SLEVE-hybrid for short) coordinates in numerical forecasting model is not limited to one. It is divided into the semi-analytical scheme and the finite differential scheme in terms of the various differential methods of the coordinate deformation variables. Comparing the dynamic equation and the long-time batch simulation results of the two schemes, the present study draws the following conclusions. The first- order finite difference accuracy of the coordinate deformation variables in the finite differential scheme is theoretically lower than that in the semi-analytical scheme. The larger the vertical gradient of the layer thickness is, the larger the relative errors of the finite differential scheme are. The long-time batch simulation test in the GRAPES model dynamic core demonstrates that the bias of the temperature and the geopotential height in the semi-analytical scheme is smaller under the default layering, while the simulation difference of the two schemes is greatly reduced when the layering is more uniform.     

数理统计方法在区域重力调查资料解释中的应用

在研究区域重力场解释中的一些问题的时候,如何客观地描述重力场特征,研究区内的基本构造格架;在无深部资料的情况下,如何获得区内地壳结构的基本模式;通过布格重力异常特征线,重力异常与高程的相关分析、根据研究区内的布格重力异常进行定性分析,结合地质和物性资料做出定量解释。我们用上述方法在近年来的1:20万区域重力调查资料解释中取得了一些有益的地质成果。     

GRAPES全球模式的模式误差估计

现代数值天气模式考虑的物理过程和边界条件越来越复杂, 但是它描述的大气状态和真实的大气流体运动轨迹还有一定的差距, 存在模式误差。在以往的研究中, 模式误差往往被忽略, 在集合卡尔曼滤波同化系统中, 如果忽略模式误差会导致滤波发散现象。本文用不同分辨率的模式预报差异估计了GRAPES全球模式的模式误差, 研究发现模式误差随着分辨率降低而线性增加, 而且模式误差随着预报时效的增加呈现线性增长的趋势。     

高阶正定守恒的中心约束多矩有限体积平流模式

采用新的均匀三点中心约束多矩有限体积方法（3-point Multi-moment Constrained finite-Volume scheme for Uniform Points with Center Constraints, MCV3_UPCC）,发展了一个三阶正定守恒的平流模式。三点多矩有限体积方法在单网格内定义等距的3个自由度,采用多矩约束条件并通过控制方程获得时间演变方程。新的三点中心约束多矩方法能在单网格内采用等距的3个点值及中心一阶、二阶导数作为约束条件进行空间4次多项式数值重构,获得3个自由度的时间演变方程;所构建的新数值方案具有三阶精度,边界通量连续性保证了其数值严格守恒。为了抑制该方法的非物理数值振荡,引入了边界保型限制器技术,它能够把数值解控制在既定物理场最小值（最小值为0时则保持数值正定）与最大值之间。数值试验表明新发展的三阶平流模式具有良好的计算精度,能够严格保持数值解的正定性和守恒性,同其他高精度平流模式相当,在实际大气模式水汽等平流输送应用中具备良好的发展潜力。      

逐层平滑地形的平缓地形追随坐标在高分辨率GRAPES模式中的应用研究

高分辨率GRAPES模式如3?km模式对地形的识别程度更高,模式中各高度坐标面可识别的地形坡度也更大,地形作用带来的气压梯度力计算误差和平流输送误差更突出.平缓地形追随坐标可以通过多种方式衰减坐标面上的地形影响进而减小这些计算误差.选择一种逐层平滑地形的平缓地形追随坐标,基于GRAPES-3km模式进行理想试验和批量模...     

一种改进的平缓-混合地形追随坐标在GRAPES中尺度模式中的应用研究

平缓-混合地形追随坐标（T-F坐标）可以减小坐标面上的地形影响带来的各种计算误差。以余弦三角函数为基函数的平缓-混合坐标（COS坐标）高层坐标面水平,计算误差较小,但是低层坐标面之间的厚度较薄,增大了计算误差,给模式稳定性及模拟效果带来较大的影响。设计一种改进的COS坐标,使低层坐标面垂直分布更加均匀,应用于GRAPES-Meso模式进行理想试验和实际模拟试验。结果表明,改进的COS坐标相对COS坐标,中高层计算误差相当,低层地形作用衰减的垂直变化更加均匀,减小了计算误差,提高了计算稳定性;地形重力波试验结果显示,改进的COS坐标重力波破碎相对COS坐标有一定缓解,更接近解析值;批量模拟试验结果显示,改进的COS坐标各个层次上的月平均模拟偏差比单尺度双曲函数平缓-混合坐标（简称SLEVE1坐标）更小,均方根误差减小,距平相关系数增大。改进的COS坐标有效地解决了COS坐标的计算问题,提高了模式预报效果。      

非静力中尺度高分辨率模式模拟中的垂直坐标影响研究

近几年来,随着高性能计算机技术的高速发展,甚高分辨率的中尺度数值预报模式业务应用已成为可能,与之相关的一系列模式技术的新问题也随之提了出来,垂直坐标系的影响就是其中之一。文中借助于美国新一代数值预报模式WRF(WeatherReseachandForecast),比较了非静力中尺度模式高分辨率模拟应用的垂直坐标影响问题。研究表明,当选用几何高度(z)和气压(p)来构造地形追随坐标时,低层两坐标引起的误差基本一样,中高层高度地形追随坐标引起的误差小于气压地形追随坐标;而且分辨率越高差异愈大。高分辨率模拟结果也表明,这种差异趋势是存在的;此外,本文对天气过程的预报要素场进行了相关的分析。     

PRM标量平流方案在GRAPES全球预报系统中的应用

如何更好地模拟水物质的分布,对于数值天气预报效果的改进,特别是对于更好地模拟降水过程,具有重要的意义。半拉格朗日模式中的标量平流计算要求做到高精度、守恒、正定和保形,但GRAPES_GFS (Global-Regional Assimilation and PrEdiction System, Global Forecast System) 中采用的QMSL(Quasi-Monotone Semi-Lagrangian)平流方案在水汽的强梯度、不连续区域计算精度较低,且不能做到严格守恒。本研究借鉴计算流体力学领域的研究进展,将一个基于分段有理函数的物质平流方案PRM(Piecewise Rational Method)引入GRAPES_GFS中,按照通量形式求解水汽方程,并对极区进行了混合等技术处理。通过一系列理想试验对两种平流方案进行了对比,证明了PRM方案精度较高,特别是在水汽梯度大的区域优势明显,频散、耗散误差较小,守恒、保形性也要好于QMSL方案。通过对GRAPES_GFS中批量预报试验效果的检验,验证了PRM方案可以有效地改进模式对水物质分布的模拟,提高了降水的预报效果,对模式综合预报性能的提升也有明显作用。