细网格非静力二维高阶闭合模式对山体流场特征的模拟分析

本文采用二阶矩湍流闭合方案，分别就非静力和准静力两种条件建立ＰＢＬ数值模式，并模拟了二维山脊流场的结构，计算了二维山体对气流平均场和湍流场的影响。计算结果表明，无论是平均场还是湍流场，采用准静力假设都会造成模拟结果的较大偏差，且对平均场的影响比对湍流场的影响更为显著。两种模式计算结果的差别还随地形坡度和背景风增大而增大。     

三维非静力二阶闭合PBL模式的初步研究

采用二阶距湍流闭合方案，建立了一个复杂地形下的三维非静力边界层预报模式。通过该模式，细致了模拟了北京地区的风、温时空变化和湍流场变化。结合地面、探空和气象塔的实测资料，进一步研究了该方案在三维空间上对平均场和湍流场的预报模拟效能。结果表明该模式能较好地进行平均量和湍流量的模拟和预报。     

天津市郊大气边界层湍谱特征分析

用１９９０年７月在天津市郊的观测资料，计算了风速谱、温度谱及动量和热量通量的协谱。结果表明，市郊下垫面上的大气湍流谱特征与平坦地形上得到的典型结果基本一致，但满足各向同性的湍流尺度似乎比在平坦地形上的略大，谱的峰值频率区域略窄，谱的低频区（含能区）特征与平坦地形上的情况也不同。     

细网络非静力二维高阶闭合模式对山体流场特征的模拟分析

本文采用二阶矩湍流闭合方案，分别就非静力和准静力两种条件建立ＰＢＬ数值模糊，并模拟了二维山脊流场的结构，计算二维山体对气流平均场和湍流场的影响，计算结果表明，无论是平均场还是湍流场，采用准静力假设都会造成模拟结果的较大偏差，且对平均场的影响比对湍流场的影响更为显著，两种模式计算结果的差别还随地形坡度和背景风增大而增大。     

大地形上边界层流场的动力学研究

本文应用边界层气象学中Estoque数值模式关于湍流交换系数及分层的处理方法,求得了大地形存在时定常边界展方程的零、一级解析解,并用来得到大地形存在时边界层的散度场、垂直速度场,改进了前人的结果。     

地形湍流拖曳力参数化及在GRAPES中的应用

分辨率的限制使得不能被模式识别的地形称为次网格尺度地形,次网格尺度地形在热力和动力方面对实际大气有着不可忽略的作用,其效应只能通过参数化的形式回馈给模式。分辨率的提高使得与较小尺度地形相联系的地形湍流拖曳力凸显其重要性。数值模式中地形湍流拖曳力的参数化对完善模式物理过程和改善模式近地层预报效果具有积极意义,其方法包括有效粗糙度法和直接参数化法,而GRAPES模式中并未以任何方法考虑次网格尺度地形的影响。该文通过单柱模式比较了有效粗糙度法和直接参数化法的优劣, 发现后者在有些方面优于前者。最后,将应用于实际的一个直接参数化方案接入GRAPES中尺度模式中,进行个例模拟,并与NCEP再分析资料进行对比,结果表明:考虑地形湍流拖曳力方案对模式预报具有改进作用,尤其对局地低层风场具有积极影响。     

西双版纳地区雾的数值模拟研究

建立了一个适用于复杂地形上的三维非定常雾模式,用于研究西双版纳地区雾的生消规律。模式详细考虑了湍流、长短波辐射、凝结、蒸发、重力沉降等影响因子,特别注意了植被和气溶胶辐射效应对雾的生消的影响。利用本模式对山谷风及雾的生消过程进行了实况模拟,结果表明,模式在一定程度上反映了实际情况。     

三维非静力E—ε闭合模式对山体流场及浓度分布的模拟

本文采用三维非静力E-ε闭合模式,在对流层边界层方程组中,加入湍能和耗散率方程描述边界层中的湍流运动。求解孤立三维山地地形上的流场、湍流场和浓度场,并分析地形对气流的影响及对污染物浓度分布的影响。同时做了不同坡度和不同风速下的孤立三维山体的流场、湍流场和浓度场比较。数值计算结果表明,地形能改变气流运动,气流过山时将出现分支现象,山后背坡有反向回流出现,并且回流区大小随着山体坡度的增大而增大。山体对气流的阻挡作用随山体高度的增加而更明显,随风速的增加而减弱。湍能及耗散率随高度衰减较快。随着山体坡度增加,山前气流的垂直速度也增加,湍流运动更强。风速大时湍流运动越显著。污染物在山前出现辐散,迎风坡是浓度高值区,山后背风坡也是相对高值区。山坡度越大,山前迎风坡污物浓度越高。     

三维山体过山气流流场特征的数值模拟

采用高阶矩湍流闭合方案，建立了一个细网格、高分辨率、三维非静力ＰＢＬ数值模式，并由其模拟了三维山体条件下的流场结构和湍流场特征。为反映数值模拟结果的可靠性，对中性条件下三维山体流场进行了风洞试验。与数值模拟结果对比分析表明，数值计算与风洞试验结果有较好的一致性；使用该模式模拟山体条件下的流场结构能取得较好的结果；将模拟结果作为随机游动扩散模式的三维风场及湍流场资料输入，为复杂地形条件下大气污染的研     

采用隐式立方样条计算平流过程的数值模式及理想实验

本文介绍了一种用隐式立方样条求解平流方程的数值方法,并从理论上对其无条件稳定性进行了证明,在此基础上建立了一个在地形坐标系下的两维原始方程模式,模式在行星边界层参数化中引入了湍流动能方程,在模式顶部引入了吸收层.数值实验表明:模式有较好的计算稳定性,对较高的模式水平分辨率和复杂地形均有较强的适应能力;对复杂地形和下垫面非均匀热源条件下中尺度系统的模拟能获得合理的结果,并具有较高的精确度.     

A global spectral element model for poisson equations and advective flow over a sphere

A global spherical Fourier–Legendre spectral element method is proposed to solve Poisson equations and advective flow over a sphere. In the meridional direction, Legendre polynomials are used and the region is divided into several elements. In order to avoid coordinate singularities at the north and south poles in the meridional direction, Legendre–Gauss–Radau points are chosen at the elements involving the two poles. Fourier polynomials are applied in the zonal direction for its periodicity,with only one element. Then, the partial differential equations are solved on the longitude–latitude meshes without coordinate transformation between spherical and Cartesian coordinates. For verification of the proposed method, a few Poisson equations and advective flows are tested. Firstly, the method is found to be valid for test cases with smooth solution. The results of the Poisson equations demonstrate that the present method exhibits high accuracy and exponential convergence. Highprecision solutions are also obtained with near negligible numerical diffusion during the time evolution for advective flow with smooth shape. Secondly, the results of advective flow with non-smooth shape and deformational flow are also shown to be reasonable and effective. As a result, the present method is proved to be capable of solving flow through different types of elements, and thereby a desirable method with reliability and high accuracy for solving partial differential equations over a sphere.     

重庆山地区域气象要素空间插值方法对比

利用重庆地区1999年和2018年气象数据, 分别采用薄盘光滑样条、协同克里金、普通克里金、反距离加权4种方法, 从年和月两种尺度对气温、降水、太阳总辐射三个要素进行空间插值; 采取交叉验证方法, 用MAE、MRE、RMSE评估插值精度, 确定各要素最优插值方法。结果表明: 气温和太阳总辐射最优插值方法为薄盘光滑样条, 降水为反距离加权; 插值精度上气温、太阳总辐射高值月份优于低值月份, 降水则相反, 但三个要素均表现出年尺度优于月尺度。MRE检验表明, 插值精度为气温>太阳总辐射>降水, 1999年年尺度插值精度分别为1.86%、4.60%、6.87%, 月尺度插值精度分别为2.79%、5.82%、17.42%;2018年太阳总辐射年、月尺度插值精度分别为3.03%、4.88%, 区域站加密后气温、降水年尺度插值精度分别为2.03%、11.20%, 月尺度对应插值精度分别为3.20%、23.14%。     

OBJECTIVE ANALYSIS OVER INDIAN REGION BY BRATSETH''S SUCCESSIVE CORRECTION SCHEME

Bratseth (1986) has proposed a successive correction scheme for numerical analysis where the solution converges towards the solution obtained by optimum interpolation. This paper presents the results of some experiments carried out over indian reston during summer monsoon period using this Bratseth's method. The convergence of the method is investigated.     

Investigation of a Complex Nocturnal Flow in Owens Valley,California Using Coherent Doppler Lidar

A study of an interesting meteorological episode over the Owens Valley, California, USA during the Terrain-Induced Rotor EXperiment was conducted using a recently adapted statistical interpolation method to retrieve wind-velocity vectors from Doppler lidar data. This vector retrieval method has been adapted from radar data assimilation techniques. Results show that the method allows better preservation of local variations in the flow field than other techniques. In addition, a high resolution Coupled Ocean/Atmosphere Mesoscale Prediction System (COAMPS^®) run is used to understand the large-scale flow within the valley and compared with lidar retrievals. Observations from 1030 UTC to 1230 UTC (0230 local time to 0430 local time) on March 27, 2006 are presented. Lidar observations show complex and uncharacteristic flows such as sudden bursts of westerly cross-valley wind mixing with the dominant up-valley wind. Model results from COAMPS and other in-situ instrumentation are used to corroborate and complement these observations. The optimal interpolation technique for Doppler lidar data vector retrieval appears well suited for scenarios with complex spatial variations in the flow field.     

基于DEM修订的新疆区域气温空间插值精度优化研究

为研究新疆地区气温的空间变异性,以新疆66个国家气象台站1981—2010年月平均气温和30 m空间分辨率DEM数据为基础,采用传统插值法、基于DEM多元线性回归插值和基于DEM修正的空间插值方法对新疆区域气温数据进行栅格化,并分析年平均气温与海拔的相关关系。通过采用反距离权重法(IDW),普通克里格法(Kriging),样条函数法(Spline)和趋势面分析法(Trend)4种空间插值方法对气象要素进行直接插值、气温多元回归模型残差结果插值、基于DEM修正插值对比分析。通过针对插值方法进行基于MAE和RMSIE的交叉验证,结果表明传统插值方法、基于多元线性回归和基于DEM修正4种空间插值精度均为IDWKrigingSplineTrend。反距离权重(IDW)空间插值方法最优,基于DEM修正IDW插值、基于多元线性回归IDW插值与传统IDW插值精度分别是0.039、0.477、1.038,插值结果客观的表达了新疆区域气温随空间梯度的变化趋势。     

An unsaturated soil water flow problem and its numerical simulation

1.IntroductionUnsaturatedsoilwaterflowisaflowwherewaterisnotfullofsoilhole,whichisanimportantformofflowinporousmedia.Predictionofanunsaturatedflowisprovidedwithsignificanceinmanybranchesofscienceandengineering.Theseincludeatmosphericscience,soilscien...     

基于台站降水资料对不同空间内插方法的比较

如何把离散的气象台站资料通过合适的空间内插方法转变成规则的网格数据,对于气候变化分析和模拟研究具有重要的意义.作者利用中国区域160个常用台站10年降水观测资料为例,分别采用克里格(Kriging)插值、反距离加权、Delaunay三角剖分线性插值、双谐样条(Biharmonic Spline)插值和Cressman客观分析等几种常见的空间内插方法,较为系统地分析和比较了这几种内插方法插值结果之间的异同,对其优缺点和适用范围进行了适当讨论,并采用不同的网格分辨率和影响半径对Cressman客观分析方法做了进一步分析比较.结果表明:在台站分布密集的区域不同空间内插方法之间差异较小,在台站分布稀疏的区域则差异较大;与更加密集的台站观测资料的比较显示,自动调节影响半径的Cressman客观分析方法与双谐样条插值方法误差相对较小.     

Optimization of Upstream Profiles in Modelled Flow Over Complex Terrain

This paper describes an adjoint method for data assimilation intoupstream boundary conditions of numerical modelsusing optimal control theory. Mathematical formalisms are given along with the numerical implementation of the schemein a column model of the atmospheric boundary layer. The optimized mean and turbulence profiles are used as an upstream solutionin a model of turbulent flow in complex terrain. To contrast thiswith other methods, two solutions for flow over an isolatedhill are calculated, one with an optimized upstream solution andone with a simple surface-layer formulation for the upstream solution.These two solutions are compared to observations and analytical theory. The adjoint optimization method is shown to producesolutions of flow in complex terrain that are substantively differentat the two solutions, with the optimized solution giving more accurate results.     

Spatial analysis of monthly precipitation in Turkey

Summary  The principal objectives of this paper are to develop and validate an optimum interpolation method for the spatial analysis of monthly precipitation in Turkey. A two-dimensional optimum interpolation objective analysis scheme has been developed for the spatial analysis of precipitation. The model is developed for generating statistically optimum interpolation based on the irregular distribution of meteorological stations. One question that affects the optimum interpolation method and, indeed, all such techniques, is how many observations should be allowed to influence a given grid point? The method developed in this paper addresses this question. For the implementation of the method, 52 stations are considered for Turkey, with 30 years of monthly data at each point. It is observed that each monthly average spatial correlation function shows a monotonically decreasing pattern based on 15 km interval averages. The method provides high estimation accuracy in dense station locations such as in northwestern Turkey. Precipitation contour maps obtained by the optimum interpolation method indicate two spatial trends over Turkey which have not been identified in any previous study. Received June 24, 1999/Revised April 26, 2000