三维变量配置对惯性重力波频散性模拟的影响

在线性斜压原始方程组的基础上，从频率和群速方面讨论了由水平网格 (C、Z网格) 和垂直网格 (L、CP、LZ、LY网格) 组合而成的几种三维网格 (C/L、C/CP、Z/LZ、Z/LY) 的计算频散性并分析了各种网格出现偏差的原因，结果表明三维网格C/CP (水平网格为C网格垂直网格为Charney-Phillips网格) 与Z/LZ（水平网格为Z网格垂直网格为LZ网格）计算频散性能较好。从而为原始方程大气模式选取三维网格提供指导。     

几种垂直跳点网格计算频散性的比较研究

为证明两种新的垂直跳点网格--LZ网格和LY网格的有效性,从频率、群速和出现错误群速的垂直尺度范围等方面与现有的几种垂直跳点网格和时间--垂直跳点网格进行了比较.结果表明:LZ网格计算频散性与目前使用的最广泛的并且计算频散性也是最佳的CP网格相当,而LY网格次之,从而为模式设计者选取垂直网格提供了又一种选择.     

非静力滞弹性模式中水平分辨率与垂直分辨率的协调性初探

为了理论上研究非静力滞弹性模式中水平分辨率与垂直分辨率之间的搭配问题，以线性斜压滞弹性方程组为出发方程，导出解析的频散方程、水平群速方程和垂直群速方程，同时在最优三维网格C／CP下进行离散，也导出频散方程、水平群速方程和垂直群速方程，然后用图示的方法，比较了在不同的水平和垂直格距下模拟解析解所产生的误差。结果表明：当垂直格距和水平格距相当时，在C／CP网格上对惯性重力波的频率、水平群速和垂直群速都能较好地模拟，产生的相对误差均在5％以下。     

非静力模式中几种垂直网格计算特性的对比研究

从非静力滞弹性方程组出发，将惯性重力波的水平波长分几百公里、几十公里和几公里三种情况，从频率、群速和各速的垂直尺度范围等方面对现有的几种垂直网格进行比较讨论。结果表明：在非静力模式中，CP网格和LZ网格适合于各种不同水平尺度惯性重力波的模拟；L网格和LY网格适合于水平尺度为几十公里以上的惯性重力波的模拟；LTS网格和CPTS网格比较适合于水平尺度的几十公里以下的惯性重力波的模拟。     

垂直变量配置对静力适应过程的影响

为了研究垂直变量配置对静力适应过程的影响,本文从描写静力适应过程的方程组出发,分别在将所有变量置于整层上的非跳点N网格;将垂直速度和温度放置在整层,水平速度、气压和密度等变量放置在半层的Charney Phillips跳点网格(CP网格);将水平速度、气压和温度放置在整层,将垂直速度和密度放置在半层的Lorenz跳点网格(L网格);将密度变量放置在整层的Charney Phillips跳点网格(CP_N网格);将密度放置在整层的Lorenz跳点网格(L_N网格)上进行离散,垂直格距分1 km、0.5 km、0.2 km和0.01 km,研究了在这5种网格上产生的频率和垂直群速的相对误差。结果表明:(1)L_N网格和CP网格是完全等效的两种网格。(2)不论垂直格距为多少,CP网格和L网格的误差都是最小,N网格次之,CP_N网格的误差最大。(3)随着垂直格距的减少,在这几种网格上产生的误差都在减小。对于CP网格、L网格和N网格,在水平长波和垂直短波处产生的误差较大。而CP_N网格对水平波长变化不敏感,垂直波长越短,误差越大。(4)当垂直格距为0.01 km时,这几种网格都对水平波长的变化不敏感了,仅对垂直波长敏感。(5)CP网格、L_N网格和L网格在描写静力适应过程和斜压地转适应过程都是误差最小的垂直变量配置方案,因此在非静力完全可压缩深层大气数值预报模式中应优先选择这3种方案。     

水平网格计算频散性的研究

从线性浅水方程组出发, 在Arakawa A—E网格、Z网格和Eliassen网格上, 分可分辨和不可分辨两种情况, 从频率和群速方面讨论了这7种网格的计算频散性, 并指出了每种网格出现错误群速的水平尺度范围。结果表明:在可分辨的情况下, Z、C和Eliassen网格较其它网格效果好; 在不可分辨的情况下, B网格较其它网格效果好。     

水平网格对Rossby波的影响

文中从描述准地转、β平面近似下的Rossby波方程组出发,在Arakawa A-E网格和Z网格上,分可分辨和不可分辨两种情况,从频率和群速方面讨论了这6种网格的计算频散性。结果表明:在可分辨的惰况下,Z,D,C网格比其它网格对Rossby波产生的歪曲都要小;而在不可分辨的情况下,这6种网格的描述能力都比较差。     

水平变量配置对地转演变模拟的影响

从描述准地转、β-平面近似下的Rossby波方程组出发，在Arakawa A-E网格上，从频率和群速方面讨论了这五种网格的计算频散性。结果表明：D、C网格比其它网格对Rossby波产生的歪曲都要小。     

静力模式中水平分辨率与垂直分辨率的协调性问题

本文将线性斜压适应方程组在最优三维网格C/CP下进行离散,讨论了不同的水平和垂直格距大小对解析解模拟的影响。结果表明:垂直格距为1 km量级时,水平格距应为100km量级;垂直格距为100 m量级时,水平格距应为10 km量级较好。     

水平离散网格与大气运动

从线性浅水方程组出发,在Arakawa A-E网格和z网格上,从频率和群速方面讨论了这六种网格对大气中的两种重要的波动：惯性重力波和Rossby波的描述能力。结果表明：C、z、E网格比其它网格产生的歪曲都要小。     

THREE-DIMENSIONAL VARIABLES ALLOCATION IN MESOSCALE MODELS

Forecasts and simulations are varied owing to different allocation of 3-dimensional variables in mesoscale models. No attempts have been made to address the issue of optimizing the simulation with a 3-dimensional variables distribution that should come …     

食指-拇指夹捏运动的动力学建模及控制

本文研究了具有冗余肌肉的双自由度手指三维夹捏运动的运动学分析和动力学建模.拇指和食指均采用单关节肌和双关节肌模型,食指-拇指通过对向作用实现夹捏动作.在夹捏运动中,物体可以沿两个轴旋转,从而使物体实现时变的三维夹捏运动.考虑到屈肌和伸肌的共同收缩、非线性约束及重力的影响建立了动力学模型,并通过简单的控制对所提出的动力学模型的稳定性进行了数值仿真实验验证.     

Boundary-layer pollutant concentrations over complex terrain

Given the incident profiles of wind velocity and pollutant concentration, we seek to determine the 3-dimensional concentration field of a pollutant upon a region with complex terrain. The analytic solution of the wind velocity in a 3-dimensional boundary-layer model by Walmsley et al. (1980) is utilized as a forcing function in the simplified concentration perturbation equation for a pollutant. The resulting solution applies to an isolated cosine-squared hill in a neutrally stratified boundary-layer flow with a surface type which absorbs the pollutant totally. The solution shows that the concentration perturbation field is organized in accordance with the wind field. In particular, the east-west cross-section is 180° out of phase with the velocity perturbation field. The vertical profiles of the concentration perturbations for selected grid points approach the value of the upper boundary condition very rapidly.     

Flux agreement above a Scots pine plantation

Summary The surface energy exchange of 12m high Scots pine plantation at Hartheim, Germany, was measured with a variety of methods during a 11-day period of fine weather in mid-May 1992. Net radiation and rate of thermal storage were measured with conventional net radiometers, soil heat flux discs and temperature-based storage models. The turbulent fluxes discussed in this report were obtained with an interchanging Bowen ratio energy budget system (BREB, at 14 m), two one-propeller eddy correlation systems (OPEC systems 1 and 2 at 17m), a 1-dimensional sonic eddy correlation system (SEC system 3) at 15 m, all on one low tower, and a 3-dimensional sonic eddy correlation system (SEC system 22) at 22 m on the high tower that was about 46 m distant. All systems measured sensible and latent heat (H and LE) directly, except for OPEC systems 1 and 2 which estimated LE as a residual term in the surface energy balance. Closure of turbulent fluxes from the two SEC systems was around 80% for daytime and 30% for night, with closure of 1-dimensional SEC system 3 exceeding that of 3-dimensional SEC system 22. The night measurements of turbulent fluxes contained considerable uncertainty, especially with the BREB system where measured gradients often yielded erroneous fluxes due to problems inherent in the method (i.e., computational instability as Bowen's ratio approaches –1). Also, both eddy correlation system designs (OPEC and SEC) appeared to underestimate |H| during stable conditions at night. In addition, both sonic systems (1- and 3-dimensional) underestimated |LE| during stable conditions. The underestimate of |H| at night generated residual estimates of OPEC LE containing a phantom dew error that erroneously decreased daily LE totals by about 10 percent. These special night problems are circumvented here by comparing results for daytime periods only, rather than for full days. To summarize, turbulent fluxes on the low tower from OPEC system 2 and the adjacent SEC system 3 were in reasonable agreement, while the BREB system appeared to overestimate H and underestimate LE; H and LE measured by SEC system 22 on the high tower were lower than from OPEC and SEC3 on the low tower. The turbulent flux measurements tended to converge, but the data exhibit unexplained differences between days, between systems, and between locations.With 7 Figures     

一次对流层异常臭氧次峰的观测研究及其动力输送过程的分析

文中分析了 1996年 8月 1日发生在西宁 (36 .4 3°N ,10 1.4 5°E ,海拔 :2 2 96m)地区对流层异常臭氧次峰现象。观测资料揭示了高空低压槽东移是臭氧次峰的主要天气特征。三维后向轨迹计算表明 ,尽管代表臭氧次峰的气团可以追溯到中亚地区 ,但是明显的气团向下输送则发生在新疆、青海间的高空低压槽内。中尺度模拟进一步确认了对流层顶折叠和平流层向下输送是臭氧次峰出现的动力机制。臭氧次峰在对流层高度位置与准无辐散层有关     

两种样本生成方法对DRP-4DVar同化方法性能影响的比较及分析

降维投影四维变分同化方法(DRP-4DVar)利用历史预报的集合来统计背景误差协方差,并将分析变量投影到样本空间下求解代价函数,因而集合样本的质量对DRP-4DVar同化方法的性能有着重要影响.文中尝试使用三维变分(3DVar)控制变量的扰动方法来产生集合样本,并与原来的历史预报扰动方法做比较.历史预报扰动样本具有随流...     

基于模式约束三维变分技术的连续循环同化试验研究

由于模式约束三维变分同化技术中考虑了模式的动力和物理过程,因此能保证各物理量间的平衡关系,从而滤除由于观测资料引入导致的高频波动,减小模式与初始场的协调时间.由于能在较短时间调整到稳定状态,采用模式约束三维变分同化进行连续循环同化可用较少的计算量达到同化多时次的多种观测资料的目的.该研究利用模式约束三维变分技术,针对2006年"桑美"台风个例,进行了连续循环同化卫星云导风、QuikSCAT海面风、Bogus海平面气压的试验.在台风数值预报中往往需要使用经验构造的台风信息(如Bogus海平面气压,Bogus风场等),该研究采用了模式约束三维变分同化技术同化Bogus海平面气压.由于模式约束三维变分同化技术充分考虑了各物理量间的约束,因此通过同化Bogus海平面气压也调整了初始场中相应的高度场、温度场、风场等变量,使得初始场中的台风涡旋具有较强的协调性,提高了对台风的模拟能力.采用AVN模式6小时间隔的分析场作为侧边界,2006年8月8日20时的分析场作为初估场,文中对8月8日20时到9日05时"桑美"台风的观测资料进行了连续循环同化.采用连续循环同化后台风路径的模拟精度得到了显著提高,对台风降水结构等的模拟也得到了改善.     

2008年初冰雪灾害和华南准静止锋

2008年1月中、下旬中国南方广大地区受到低温雨雪冰冻的极端异常气象灾害(简称"0801冰雪灾害").其降水性质属于隆冬季节华南准静止锋产生的持续性降水.以1月28日的锋面为例,用经典天气学分析方法从等压面图、垂直剖面图和单站探空曲线,讨论了华南准静止锋的三维空间配合.     

The importance of mesoscale phenomena on the dispersion of pollutants in the vicinity of a point-source emission

The effects of mesoscale meteorological phenomena on the transport of pollutants were analyzed by 3-dimensional dispersion calculations in case of a hypothetical nuclear accident. The results were in agreement with the earlier results obtained by Kestiet al. (1988) clearly indicating the importance of mesocale meteorological phenomena in short-range dispersion of pollutants.     

Initialization and simulation of a landfalling typhoon using a variational bogus mapped data assimilation (BMDA)

Summary Recently, a new data assimilation method called “3-dimensional variational data assimilation of mapped observation (3DVM)” has been developed by the authors. We have shown that the new method is very efficient and inexpensive compared with its counterpart 4-dimensional variational data assimilation (4DVar). The new method has been implemented into the Penn State/NCAR mesoscale model MM5V1 (MM5_3DVM). In this study, we apply the new method to the bogus data assimilation (BDA) available in the original MM5 with the 4DVar. By the new approach, a specified sea-level pressure (SLP) field (bogus data) is incorporated into MM5 through the 3DVM (for convenient, we call it variational bogus mapped data assimilation – BMDA) instead of the original 4DVar data assimilation. To demonstrate the effectiveness of the new 3DVM method, initialization and simulation of a landfalling typhoon – typhoon Dan (1999) over the western North Pacific with the new method are compared with that with its counterpart 4DVar in MM5. Results show that the initial structure and the simulated intensity and track are improved more significantly using 3DVM than 4DVar. Sensitivity experiments also show that the simulated typhoon track and intensity are more sensitive to the size of the assimilation window in the 4DVar than that in the 3DVM. Meanwhile, 3DVM takes much less computing cost than its counterpart 4DVar for a given time window.