Study of near-surface models for large-eddy simulations of a neutrally stratified atmospheric boundary layer

In large-eddy simulations (LES) of the atmospheric boundary layer (ABL), near-surface models are often used to supplement subgrid-scale (SGS) turbulent stresses when a major fraction of the energetic scales within the surface layer cannot be resolved with the temporal and spatial resolution at hand. In this study, we investigate the performance of both dynamic and non-dynamic eddy viscosity models coupled with near-surface models in simulations of a neutrally stratified ABL. Two near-surface models that are commonly used in LES of the atmospheric boundary layer are considered. Additionally, a hybrid Reynolds- averaged/LES eddy viscosity model is presented, which uses Prandtl’s mixing length model in the vicinity of the surface, and blends in with the dynamic Smagorinsky model away from the surface. Present simulations show that significant portions of the modelled turbulent stresses are generated by the near-surface models, and they play a dominant role in capturing the expected logarithmic wind profile. Visualizations of the instantaneous vorticity field reveal that flow structures in the vicinity of the surface depend on the choice of the near-surface model. Among the three near-surface models studied, the hybrid eddy viscosity model gives the closest agreement with the logarithmic wind profile in the surface layer. It is also observed that high levels of resolved turbulence stresses can be maintained with the so-called canopy stress model while producing good agreement with the logarithmic wind profile.     

基于气象站资料的中国地区太阳日辐射量算法研究

现行计算水平面太阳日辐射主要有两种方式:一种是利用影响辐射的相关要素建立模型,另一种是依据实测资料进行空间插值。但后一种方法若要保证精度则需有足够多的样本。针对上述问题,利用我国不同区域67个站点的数据,在VP-RAD模型的基础上,建立了一个适用于中国地区的逐日太阳辐射算法CNR,该算法仅需要输入站点基本信息、最高最低温度和降水量。模拟结果与实际观测结果比较吻合。     

Turbulence Structure of the Unstable Atmospheric Surface Layer and Transition to the Outer Layer

We present a new model of the structure of turbulence in the unstable atmospheric surface layer, and of the structural transition between this and the outer layer. The archetypal element of wall-bounded shear turbulence is the Theodorsen ejection amplifier (TEA) structure, in which an initial ejection of air from near the ground into an ideal laminar and logarithmic flow induces vortical motion about a hairpin-shaped core, which then creates a second ejection that is similar to, but larger than, the first. A series of TEA structures form a TEA cascade. In real turbulent flows TEA structures occur in distorted forms as TEA-like (TEAL) structures. Distortion terminates many TEAL cascades and only the best-formed TEAL structures initiate new cycles. In an extended log layer the resulting shear turbulence is a complex, self-organizing, dissipative system exhibiting self-similar behaviour under inner scaling. Spectral results show that this structure is insensitive to instability. This is contrary to the fundamental hypothesis of Monin--Obukhov similarity theory. All TEAL cascades terminate at the top of the surface layer where they encounter, and are severely distorted by, powerful eddies of similar size from the outer layer. These eddies are products of the breakdown of the large eddies produced by buoyancy in the outer layer. When the outer layer is much deeper than the surface layer the interacting eddies are from the inertial subrange of the outer Richardson cascade. The scale height of the surface layer, z _s, is then found by matching the powers delivered to the creation of emerging TEAL structures to the power passing down the Richardson cascade in the outer layer. It is z _s = u _* ³ /k_s, where u _* is friction velocity, k is the von Kármán constant and _s is the rate of dissipation of turbulence kinetic energy in the outer layer immediately above the surface layer. This height is comparable to the Obukhov length in the fully convective boundary layer. Aircraft and tower observations confirm a strong qualitative change in the structure of the turbulence at about that height. The tallest eddies within the surface layer have height z _s, so z _s is a new basis parameter for similarity models of the surface layer.     

Wind-Tunnel and Numerical Simulations of the Coastal Thermal Internal Boundary Layer

Wind-tunnel experiments in a thermally stratified wind tunnel and direct numerical simulations were performed to simulate the thermal internal boundary layer (TIBL) that developed over a coastal area in a sea-breeze flow. The results of the simulations were analyzed to investigate turbulence structure in the TIBL. To study the effects of the atmospheric stability over the sea on the TIBL, two vertical profiles of temperature were created in the upstream portion of the wind-tunnel experiment and the direct numerical simulation. Turbulence statistics of the TIBL changed significantly according to the temperature profile over the sea, indicating that the stability of the flow over the sea has a significant effect on the structure and turbulence characteristics of the TIBL. Furthermore, the TIBL heights were estimated from the vertical profiles of the local Richardson number. The estimated TIBL heights agreed with those predicted by a pre-existing relation, suggesting that both the wind-tunnel experiment and the direct numerical simulation accurately reproduced the growth of the TIBL.     

亚欧及北太平洋相当正压层高度的分析及数值模拟

本文利用相当正压层方法,分析了亚欧及北太平洋中纬度地区47个站点的相当正压层高度,其结果与文献[1]类似。计算表明:亚欧及北太平洋的相当正压层高度平均位于359hPa,略去高原站点后,平原地区的平均高度为374hPa。相当正压层是一个波动很大的不同气压的曲面,并随地区和季节变化较大。高原和西风急流的强弱可能都是其影响因子,尤其是高原的作用值得进一步探讨。数值模拟也证实了300hPa上Rossby波比500hPa效果好。     

鄂尔多斯市城区发展对局地大气环境影响的数值模拟

选取鄂尔多斯市东胜区为代表区域。基于该区域1961—2006年定时地面气象观测资料和2002—2006年08时、20时各标准层探空以及2001—2006年主要污染源源强、1961—2006年社会经济数据、2002年卫星遥感反演资料、2003—2020年城市总体规划资料等,利用区域边界层模式模拟方法,分析鄂尔多斯市城市化进程对局地大气环境的影响。结果表明：东胜区发展后（2020年）气温低于现状（2002年）,使城区中心气温与周边地区气温差异变小,对城区大气环境的改善较为有利。东胜城区扩展后,在污染严重的冬季,气流辐合区域减少,辐散能力增强,对减轻城区污染和改善大气环境有利。城区拓展后由于下垫面等布局发生变化,气流场变化较大,但从污染物的扩散能力来看,按规划发展后城区污染轻于现状。     

乌鲁木齐重污染日的天气分型和边界层结构特征研究

利用2004年1月-2009年4月高空、地面气象观测资料和逐日空气污染指数,对乌鲁木齐空气污染≥Ⅳ级的重污染日持续时间特征、500hPa高空环流形势、地面气压场及相应的边界层结构特征进行了统计分析。结果表明,乌鲁木齐重污染过程发生1天和持续2,3和4天的比例分别是32.2%,23.3%,18.5%和11.0%;发生重污染时500hPa以纬向环流型居多,占重污染总日数的84.2%,经向环流型为15.8%。从地面气压场来看,高压后部型出现重污染的频率最高,达86.3%;高压底部型次之,为9.6%;高压前部型和南高北低型出现重污染的几率较小。乌鲁木齐冬季Ⅲ级污染日对应的温度、湿度及风等要素廓线的垂直结构与冬季平均状况几乎一致,而重污染出现时,边界层逆温较强、风速较低且低空伴随有较厚的偏东风或东南风气流;重污染日和雾的关系密切,伴随有雾或轻雾的频率高达81.3%;前一日20:00(北京时)上干下湿并伴有逆温的边界层结构极易导致空气质量恶化。     

全球大气角动量变化的比较分析

利用日本气象局ＧＳＭ９６０３大气数值模式和客观分析资料以及美国环境预报中心大气再分析数据研究比较全球大气角动量变化。分析表明，三者大气质量再分布角动量东经９０°Ｅ分量的符合程度要好于格林威治方向的分量；而且ＧＳＭ９６０３对大气轴向相对角动量季节变化的强度模拟与其他两者一致，并明显优于第一阶段 ＡＭＩＰ ２３个大气环流模式的平均模拟结果。另外，ＧＳＭ９６０３大气相对角动量的年际变化基本显示了厄尔尼诺变化的主要历程。     

南京市城市不同功能区PM10和PM2.1质量浓度的季节变化特征

使用Anderson-Ⅱ型9级撞击采样器测量了南京市鼓楼商业区、江北工业区、钟山风景区和宁六高速公路交通源春、夏、秋三季的大气气溶胶质量浓度。分析结果表明:南京市PM_2.1和PM₁₀的质量浓度存在明显的季节变化,秋季>春季>夏季;ρ_PM10春季为167.47 μg/m³,夏季为 85.99 μg/m³,秋季为238.99 μg/m³;ρ_PM2.1春季为59.66 μg/m³,夏季为42.80 μg/m³,秋季为100.15 μg/m³。不同季节中ρ_PM10和ρ_PM2.1均存在较好的相关性,夏季相关性最好,相关系数为0.952;秋季次之,相关系数为0.783;春季相对较差,相关系数为0.613。城市不同功能区之间ρ_PM2.1和ρ_PM10的质量浓度值差异很大,交通源>工业区>商业区>风景区。城市不同功能区的质量浓度谱分布基本一致,均为双峰型分布,峰值分别位于0.43~0.65 μm/m³和9.0~10.0 μm/m³。南京市春、夏、秋三个季节大气粒子质量浓度谱为双峰分布,粒子主要集中在0.43~3.3 μm/m³的粒径段。江北工业区ρ_PM10和ρ_PM2.1质量浓度的相关系数为0.814,略高于鼓楼商业区的0.797。