森林冠层和森林边界层大涡模拟

在采用各向异性湍流动能闭合方案和3阶Runge-Kutta时间积分方案的大涡模式中,引入由森林冠层粗糙元造成的动量拖曳项、热量输入项和TKE耗散项,以模拟森林冠层和森林边界层的气象场. 通过中性和不稳定层结条件下不同叶面积指数算例的模拟分析及与已有观测结果的比较表明,本文所建大涡模式对森林冠层和森林边界层有较好的模拟效果. 进一步研究表明：不稳定层结条件下较稠密的森林冠层中特有的Kinking & Pairing湍涡结构与森林边界层中湍流的大涡运动相互作用,形成了森林冠层附近的温度斜坡型结构.     

Modeling high-order statistics in the turbulent Ekman layer

Results from a direct numerical simulation (DNS) of the neutral and unstable turbulent Ekman layer at a Reynolds number of 1000 were used to evaluate turbulence closure models. For the neutrally stratified Ekman layer, the higher-order moments of velocity were examined and the accuracy of a kurtosis model was assessed. For the unstable Ekman layer, the analysis of higher-order moments was extended to temperature-velocity correlations. Model coefficients were optimised using DNS data and it was shown that the optimised models accurately captured the distributions of all fourth-order moments. These low-Reynolds number results can be extrapolated to higher Reynolds numbers to parameterise turbulence in other flow fields with rotational effects such as the atmospheric boundary layer.     

大气边界层研究进展

大气边界层对云和对流的发展、演变有重要作用.本文回顾了在大气边界层高度计算方法,边界层的时空分布特征、结构和发展机理,以及边界层参数化方案等方面的主要研究进展.大气边界层高度计算方法主要分为基于大气廓线观测数据计算和基于模式参数化方案计算两大类;大气边界层高度频率分布形态具有明显的日变化特征,并且稳定、中性和对流边界层高度的频率分布呈现出不同的Gamma分布特征;地面湿度状况对边界层发展影响明显,对于不同的下垫面热力性质和地形状况,大气边界层高度呈现出明显的空间差异,青藏高原边界层高度明显高于一般平原地区;在强烈的地面加热驱动下,对流边界层与残余层通过正反馈机制循环增长可以形成4000 m以上的超高大气边界层;研制大气边界层、浅对流以及云物理方案的统一参数化框架是未来数值预报模式的发展趋势.

     

大气边界层湍流的动力非平稳性的验证

首次用验证时间序列中是否存在动力非平稳性的一种简单图示方法——space time-index法来分析大气边界层湍流的动力平稳性特征.本文以取自淮河流域和威斯康星森林下垫面条件下的三维高精度风速和温度、湿度湍流脉动资料对大气边界层湍流的平稳性特征进行了分析.结果表明space time-index方法能有效地检验大气边界层湍流信号中是否存在动力平稳性.另外,均匀下垫面条件（水稻田）及复杂下垫面条件（森林）下的大气边界层湍流信号中几乎都存在动力非平稳性,大气湍流动力学非平稳性可能是边界层湍流信号相当普遍具有的一种特性.大气边界层湍流中的间歇性和相干结构使得其非平稳性图形的特征不同于一般时间序列非平稳性图形的“V”型特征;森林下垫面条件下的湍流信号比相对均匀下垫面（水稻田）下的湍流信号更有组织性,相干结构更强.     

Wake effects on turbulent transport in the convective boundary layer

The effect of the downstream propagation of a wake on the transport of momentum, energy and scalars (such as humidity) in the convective boundary layer (CBL) is studied using a direct numerical simulation. The incompressible Navier–Stokes and energy equations are integrated under neutral and unstable thermal stratification conditions in a rotating coordinate frame with the Ekman layer approximation. Wake effects are introduced by modifying the mean velocity field as an initial condition on a converged turbulent Ekman layer flow. With this initial velocity distribution, the governing equations are integrated in time to determine how turbulent transport in the CBL is affected by the wake. Through the use of Taylor’s hypothesis, temporal evolution of the flow field in a doubly periodic computational domain is transformed into a spatial evolution. The results clearly indicate an increase in the scalar flux at the surface for the neutrally stratified case. An increase in wall scalar and heat flux is also noted for the CBL under unstable stratification, though the effects are diminished given the enhanced buoyant mixing associated with the hot wall.     

民航气象观测的中国中东部大气边界层廓线数据的评估与应用

本文从民航气象报告(AMDAR)数据库中提取了2008—2019年期间我国中东部18个机场飞机起飞和降落期间的大气边界层(ABL)廓线数据,并利用常规探空数据评估了不同季节和高度范围的AMDAR数据准确性及其随探空站与机场距离和飞行状态的变化.与常规探空数据相比,低层(P>850 hPa)的AMDAR温度和风分量的均方根误差(RMSE)均比较高层(850~500 hPa)小,且各高度层上的温度和风分量的RMSE都随间隔距离的增加而增大.相对而言,飞行状态对AMDAR数据准确性的影响很小.在间隔距离小于50 km且P>850 hPa高度层范围内,温度、纬向风和经向风的平均RMSE分别为0.68~0.83 K、1.48~1.62 m·s^-1和1.73~1.81 m·s^-1.相比于RMSE,对应的平均偏差都显著更小,且随间隔距离无明显变化.数据应用示例表明,AMDAR廓线数据由于其较高的时间分辨率,不仅可以反映ABL结构的日变化特征,还可用于研究ABL局地环流特征.基于18个机场AMDAR廓线数据分析发现,我国远内陆、干燥地区大气边界层高度(ABLH)具有较强的日变化,尤其是在春夏两季,而沿海和近内陆地区ABLH的日变化强度和季节差异都较小.

     

Differences in the concentrations of atmospheric trace gases in and above the tropical boundary layer

Weekly air samples were collected at Cape Kumakahi (0 km) and at nearby Mauna Loa Observatory (3.4 km) which is above the boundary layer. EC/GC and GC/FID techniques were used to measure CH₃I, CHCl₃, CO and CH₄, which are largely natural in origin, and C₂Cl₄, CCl₄, CH₃, CCl₃, CCl₃F (F-11), CCl₂F₂ (F-12), CHClF₂ (F-22) and C₂F₃Cl₃ (F-113), which are due to anthropogenic activities. It was found that all these gases are significantly (0.05) more abundant in the boundary layer than above it.     

Wave energy input into the Ekman layer

This paper is concerned with the wave energy input into the Ekman layer, based on 3 observational facts that surface waves could significantly affect the profile of the Ekman layer. Under the assumption of constant vertical diffusivity, the analytical form of wave energy input into the Ekman layer is derived. Analysis of the energy balance shows that the energy input to the Ekman layer through the wind stress and the interaction of the Stokes-drift with planetary vorticity can be divided into two kinds. One is the wind energy input, and the other is the wave energy input which is dependent on wind speed, wave characteristics and the wind direction relative to the wave direction. Estimates of wave energy input show that wave energy input can be up to 10% in high-latitude and high-wind speed areas and higher than 20% in the Antarctic Circumpolar Current, compared with the wind energy input into the classical Ekman layer. Results of this paper are of significance to the study of wave-induced large scale effects.     

The bottom Ekman layer and the apparent violation of the maximum principle

The solution for the bottom Ekman layer has a somewhat counter intuitive character, which seems to violate the maximum principle: at a certain level the velocity within the Ekman layer is higher than the velocity in the geostrophic layer above. I explain this character by looking at an analogous problem in an inertial frame of reference and show that it is the result of observing the flow from a rotating frame of reference (i.e. within a system that is not in steady state). The flow in the bottom Ekman layer is a superposition of the flow that results from the force exerted on the fluid by the rotating Earth and of the flow that results from the pressure-gradient term. Therefore, at a certain level the speed is higher than the speed of the geostrophic layer above which results from the pressure gradient alone.     

Variations of thermodynamical parameters in the atmospheric boundary layer over the Deccan plateau region,India

The thermodynamical structure of the atmospheric boundary layer over the Deccan plateau region has been studied with aerological data for 1980 and 1981.The temperatures in the sub-cloud layer were lower on active monsoon days than dose on weak monsoon days. An opposite trend was noticed in the layer above the 900 mbar level. The moisture content on active monsoon days was higher than that on weak monsoon days. The profile of relative humidity above the 850 mbar level showed large deviations between the active and weak monsoon conditions. On active monsoon days the values of dry static stability were higher than those on weak monsoon days. An opposite trend was noticed in the case of moist static stability. On active monsoon days the magnitudes of thev components of wind were small compared with those of theu components, and the latter showed a gradual decrease with height.     

Atmospheric stability effect on subgrid-scale physics for large-eddy simulation

Field measurements in the atmospheric boundary layer were carried out to identify the effect of atmospheric stability on subgrid-scale physics for large-eddy simulation. The basic instrumentation setup consisted of 12 three-dimensional sonic anemometers arranged in two parallel horizontal arrays (seven sensors in the lower array and five sensors in the upper array). Data from this setup are used to compute the subgrid-scale (SGS) heat fluxes and SGS dissipation of the temperature variance under stable and unstable stability conditions. The relative contribution of the SGS vertical flux to the total turbulent flux increases when going from unstable to stable conditions. The relative importance of negative SGS dissipation (backscatter) events becomes larger under stable conditions. The model coefficients for two well-known SGS models (eddy-viscosity and non-linear) are computed. Model coefficients are found to depend strongly on stability. Under both stable and unstable conditions, large negative SGS dissipation is associated with the onset of ejection events while large positive SGS dissipation tends to occur during the onset of sweep events. These findings are also supported by conditionally sampled 2D velocity and temperature fields obtained using the 12 anemometers placed in a vertical array.     

北京城市化对夏季大气边界层结构及 降水的月平均影响

为克服针对一次或几次天气过程研究城市化对边界层结构及降水影响的局限性,尝试研究北京城市化对夏季大气边界层结构及降水的月平均影响,本文首先总结了2006年8月份的主要天气过程,分析了气象站观测的10 m高度风速、2 m高度气温、2 m高度比湿和24 h降水的月平均分布特征,然后利用WRF/Noah/UCM模拟系统,进行了该月30个个例的高分辨率数值模拟及检验分析,并通过多组不同城市化情景的敏感性试验对比分析了城市化对夏季大气边界层结构及降水的月平均影响.研究表明:本文所用对高分辨率数值模拟结果进行月平均的方法可以较合理地模拟出城市化对大气边界层结构及降水的影响,并再现观测到的各站风频差异.8月份,北京城市化对气温的影响高度白天约为800 m,近地面气温升高1℃以上;夜间约为200 m,对近地面气温的影响达到最大(1.4℃以上).白天,城市化使城市及下风向的一些区域风速略有减小;夜间,城市及周边区域200 m以下风速明显减小,且在100 m左右高度处风速减小最明显,减小达0.8 m/s以上.城市化白天使700 m以下比湿减小,近地面处减小达1.2g/kg以上,夜间使近地面空气比湿略有减小.城市化对城市区域平均降水量的影响随城市发展的不同阶段而不同.初步模拟分析表明, 北京城市化已使上风向区域以及城区三环以内降水量减少,海淀和昌平降水明显增加.     

Effects of wind shear on the atmospheric convective boundary layer structure and evolution

The article reviews past accomplishments and recent advances in conceptual understanding, numerical simulation, and physical interpretation of the wind shear phenomena in the atmospheric convective boundary layer.     

Gradient-based similarity in the atmospheric boundary layer

The “flux-based” and “gradient-based” similarity in the stable boundary layer and also in the interfacial part of the convective boundary layer is discussed. The stable case is examined on the basis of data collected during the CASES-99 experiment. Its interfacial counterpart is considered in both the quasi-steady (mid-day) and non-steady states, utilizing the results of large-eddy simulations. In the stable regime, the “gradient-based” approach is not unique and can be based on various master length scales. Three local master length scales are considered: the local Monin-Obukhov scale, the buoyancy scale, and the Ellison scale. In the convective “quasi-steady” (mid-day) case, the “mixed layer” scaling is shown to be valid in the mixed layer and invalid in the interfacial layer. The temperature variance profile in non-steady conditions can be expressed in terms of the convective temperature scale in the mixed layer. The analogous prediction for velocity variances is not valid under non-steady conditions.     

多孔弹性波方程的多尺度波场模拟

本文研究了二维多孔弹性波方程的多尺度波场数值模拟方法.该多尺度方法可采用较粗的网格计算,同时又能反映细尺度上物性参数的变化信息.文中详细阐述了多尺度模拟方法与算法,并推导了相应的计算格式.基本思想是建立粗细两套网格,在粗网格上,基于有限体积方法计算更新波场;在细网格上,计算多尺度基函数,这基于有限元方法通过求解一个局部化问题得到.对含有随机分布散射体的多孔介质模型进行了数值计算,计算中应用了完全匹配层（PML）吸收边界条件,数值结果验证了本文方法和算法的正确性和有效性.

     

一个大气边界层非线性三维数值模式及其在太湖流域的应用

本模式用以K半经验理论为基础的闭合方程组求解,引进了地形作用、动量、热量和水汽的侧向扩散和垂直湍流交换以及凝结加热等物理过程,其中湍流交换系数K是地表粗糙度、几何高度、大气稳定度及风速切变等因子的函数。本模式还考虑了山地、平原和水面在地形高度、地表粗糙度、辐射、蒸发及热交换等方面的差异,在下垫面上建立热平衡方程与边界层内控制方程耦合。将模式应用于太湖流域,计算得到的该流域边界层内温、压、湿、风的分布特征与实际情况相似,模式具有一定的实用性。     

Some effects of intra-red radiation on the dynamics of the atmospheric boundary layer

An approximate infra-red radiation scheme, employing essentially the cooling to space approximation, is included in a one-dimensional model of the atmospheric boundary layer. The approximate scheme is found not to produce significant errors in the behaviour of the dynamical model when integrated over a few hours. Radiative cooling is shown to be important in the development of a clear well-mixed layer which is capped by an essentially dry region; in particular, convetive instability is induced and this enhances the rate of entrainment of dry air. The development of fog is found to require sound models of both radiative transfer and turbulent diffusion.     

完全匹配层吸收边界在孔隙介质弹性波模拟中的应用

模拟弹性波在孔隙介质中传播,需要稳定有效的吸收边界来消除或尽可能的减小由人工边界引起的虚假反射. 本文在前人工作基础上,首次建立了弹性孔隙介质情况下完全匹配层吸收边界的高阶速度-应力交错网格有限差分算法,并详细讨论了完全匹配层的构建及其有限差分算法实现. 首先,本文通过均匀孔隙模型的数值解与解析解的对比,验证所提出的数值方法的正确性;然后,本文考察了完全匹配层对不同入射角度入射波和自由表面上的瑞利波的吸收性能,将完全匹配层与廖氏和阻尼吸收边界进行了对比,研究了这三种吸收边界在不同吸收厚度情况下对弹性波吸收能力. 数值结果表明,在孔隙介质中,完全匹配层作为吸收边界能十分有效地吸收衰减外行波,无论对体波还是面波,是一种高效边界吸收算法.     

人为热源对城市边界层结构影响的数值模拟研究

用南京大学多尺度模式系统在不同区域进行了多种人为热源引入方案的研究,结果表明：考虑时空变化的人为热源按比例分别引入到地表能量平衡方程和大气热量守恒方程是将人为热源引入模式的最优方案.人为热通量密度变化的敏感性试验结果发现：人为热源的存在对城市热岛的生成有重要作用.计算结果表明,南京现有的人为热源排放量对该地区的城市热岛贡献率约为29.6％,若人为热通量密度在现有量值的基础上增大1倍,则其热岛贡献率可达42.9％;此外,人为热的排放对清晨城市边界层逆温结构有一定程度的破坏作用,能明显升高夜间近地层气温达0.5～1.0℃,并能使白天湍流活动的影响范围增大,混合层高度抬高,使夜间城市热岛环流的影响范围扩大.