层结性对大气边界层风场的影响

本文在地转动量[1]近似下，采用与Agee[2]公式相近的分段常数来描述湍流粘性系数K及其一阶导数（湍流摩擦力所包含的两项和均得到考虑），求得了层结大气中的边界层风场分布和边界层顶垂直速度的解析表达式。计算结果表明；不稳定层结的K比稳定层结的大，且Ekman抽吸作用强；在边界层低层，不稳定层结的全风速比稳定层结的大，而在中层以上则相反。还有，摩擦力中的一阶导数项和二阶导数项量级相同，以往忽略一阶导数项的做法会导致在边界层低层，不稳定层结的全风速反而较小，这与事实不符。因而考虑一阶导数项将进一步完善边界层动力学的理论和应用。     

非线性多层行星边界层风场

本文利用小参数展开法研究非线性行星边界层,将其分为近地层和Ekman层,而Ekman层又分为若干层。在近地层中,取湍流粘性系数K为高度的线性函数;在Ekman层中,取K为逼近O′Brien[1]公式的分段常数。利用上下边界条件和分层界面上由风速连续和粘性应力连续而列出的衔接条件,求得了各层的风速和边界层顶的垂直速度的解析表达式,并计算和作图分析了特定的圆形轴对称气旋和反气旋的个例,得到了一些新的结论。例如,边界层中风速与其顶部风速的最大夹角远小于Wut的结果,只有25°左右,而近地层风速却远大于Wu的结果,50米处即可达地转风的一半左右。这些都与实测更趋于一致。再者,此个例还表明,Wu所得到的关于摩擦作用造成的边界层顶垂直速度在中心处达最大值,而我们得到的却是零,这与近中心处地转风很小,使边界层顶风速很小,因而风速切变很小,摩擦作用便很微弱的结论是一致的。由此可见,边界层中的风分布对湍流粘性系数是敏感的。     

山地边界层中场变量的动力学研究

采用与实测较接近的二次函数来表达Ｅｋｍａｎ层中的湍流粘性系数Ｋ，在圆形气压场条件下，求得了山地上空边界层中的风速，进而求得散度、涡度和垂直速度等场变量随高度的分布。并作图分析了这些场变量的一些动力学特征。改进了以往在求解析解时，略去运动方程中湍流粘性力项中的关于高度的一阶导数项，以及取山坡面上风速为零作下边界条件等欠合理欠精确的做法。所求得的风速、散度、涡度和垂直速度均用简单的初等函数表示出来，有助于边界层参数化和深化对边界层动力学的认识。     

边界层特征参数对边界层顶垂直速度的影响

本文从正斜压及有层结时的边界层相似理论及阻力定律出发，由边界层顶垂直速度与地面湍应力的关系求出了层结、粗糙度、它们的水平梯度及地转风的水平梯度、斜压性对w的影响的解析式，可用于模式计算。计算结果表明层结影响可使w差1-２个量级，不稳定时粗糙度影响也使w差几倍。除地转涡度决定w外，地转风、层结稳定度和粗糙度及其水平梯度也起了重要作用，还讨论了斜压性的影响。     

北京城市化发展对大气边界层特性的影响

利用中国科学院大气物理研究所大气边界层物理与大气化学国家重点实验室的北京325 m气象塔1993年～2003年夏季 (7月～9月) 的观测资料, 统计分析了各年的风速与温度廓线分布特征。统计分析结果表明, 随着城市化的发展, 相对风速有逐年减小的趋势, 并且越靠近地面, 相对风速的减小越明显, 这反映了城市建筑对近地面层空气流动的摩擦作用。对风速廓线进行线性拟合得到风速随高度的垂直递增率, 发现无论是100 m以下的近地面层还是较高层, 风速的垂直递增率都随城市化发展存在逐年增大的趋势, 表明粗糙下垫面的影响已经向高层扩展。根据温度廓线计算了各年的温度垂直递减率, 发现其有增大的趋势, 这表明城市化发展对边界层热力结构同样有显著影响。本文还依据统计整理得到的近中性层结下的风速廓线资料, 利用莫宁-奥布霍夫相似理论计算了下垫面的空气动力学参数, 结果表明, 地表粗糙度、 零平面位移随着城市化发展皆有明显增加的趋势。同时, 分析了各空气动力学参数与平均风速及无量纲风速的关系。其中, 摩擦速度和平均风速二者基本成正相关, 且摩擦速度随平均风速的增大而增大的趋势越发明显。本文研究结果对研究城市化发展对区域大气边界层结构、 气候和环境影响有参考意义, 可为城市大气边界层模式和区域气候模式提供参数化依据。     

一次强飑线边界层气象要素模拟分析

使用WRF中尺度数值模式,分别选用两种不同边界层参数化方案（MYJ,YSU）,对2006—06—25洛阳地区一次强飑线边界层结构进行模拟,对比分析近地面层风场、温度场以及边界层日变化特征,结果发现：WRF模式基本模拟出了强飑线过程边界层变化特征;在边界层方案中,MYJ方案描述的边界层结构较YSU方案合理。这表明,用WRF模式能较好地模拟预报飑线的边界层气象要素特征。     

一次东北冷涡MCS边界层特征数值模拟分析

边界层是中尺度对流系统(MCS)动力和水汽的主要来源.针对东北冷涡中MCS的边界层特征及边界层对MCS的触发和维持机理,应用MM5模式耦合Noar陆面模式对2004年7月5日发生在辽宁中西部的强对流过程进行了数值模拟.重点分析了MCS不同发展阶段边界层三维气流结构、边界层冷丘、边界层层结结构及其在MCS触发和维持中的影响.结果表明,本次MCS在边界层有3股明显的气流汇合于雷暴区,1股是来源于东北部长白山稳定气团的东北气流,1股是西北部的西北下沉气流,1股是西南气流,浅薄的东北部底层冷空气有利于西南气流的上升,形成对流.MCS初期边界层低层的辐合强于边界层高层辐合,边界层气流旋转作用较弱,边界层辐合线是对流触发的重要因素之一.MCS边界层卜层的中尺度辐合涡旋主要由环境场辐合,边界层摩擦抽吸形成,是对流重要的能量、水汽输送系统,它伴随着强对流的发展而出现,同时对对流的维持、发展有重要的反馈作用,是边界层与自由大气进行交换的重要系统,是对流系统维持的人流的主要入口.MCS边界层冷丘内部为潜在稳定层结,特别在近地面层形成了非常稳定的层结,其温湿层结及气流结构改变了边界层人流.     

北京地区夏季边界层结构日变化的高分辨模拟对比

使用WRF中尺度数值模式, 分别选用两种不同的边界层参数化方案 (MYJ, YSU) 和3种陆面参数化方案 (SLAB, Noah, RUC), 对2004年7月1日08:00—7月4日20:00 (北京时) 北京地区夏季边界层结构进行1 km的高分辨模拟。对比分析了近地面层风场、温度场以及边界层的日变化特征, 结果发现:WRF模式基本模拟出了北京夏季边界层的日变化特征; 在边界层方案中, MYJ方案描述的边界层结构较YSU方案合理; Noah陆面模式较好地反映了城市的热岛效应; 无降水时, 风速及边界层高度对于陆面过程不敏感, 而降水发生后, 陆面过程对于边界层结构的影响增大; 各方案模拟的城区风速明显偏大, 这是因为没有充分考虑城市建筑物的阻力作用。     

对流边界层发展受覆盖逆温影响的大涡模拟研究

本文运用大涡模拟方法研究了存在覆盖逆温时对流边界层顶部的夹卷过程特征,并通过敏感性试验着重分析了在此情况之下的夹卷速度参数化问题,以及风切变的作用。模拟结果表明,初始覆盖逆温改变了夹卷层的结构特征,使得夹卷层结构参数明显增大,风切变使这一效应略有增强,以往关于夹卷速度的参数化方案不再适用。分析研究表明,在有覆盖逆温的情况下,夹卷层结构参数是控制夹卷速度的关键因子,应该将夹卷层结构参数作为变量引入夹卷速度参数化方案,该方案能够很好地预报出受覆盖逆温影响的对流边界层的发展过程。     

赤道带边界层结构的数值研究

在气压随纬度的不同分布条件下，数值求解了二维赤道带大气边界层非线性运动方程组，得到了不同气压分布下赤道带边界层运动特征，显示了平流的重要作用。在赤道带存在低压槽或气压单调随纬度增加而减少的情况下，都得出在赤道附近某纬度有一相应于赤道辐合带（ＩＴＣＺ）的垂直运动大值集中区域存在，克服了前人用临界纬度机制解释边界层大抽吸速度时略去重要的非线性项的处理上的困难。     

风沙流中风速廓线的数值模拟与实验验证

如何描述风沙流中被风沙运动改变了的风速廓线是风沙相互作用研究中的关键问题之一.该文中将跃移风沙流视为一种颗粒拟流体,将跃移颗粒对气流产生的阻力用颗粒流的阻力系数来表达,建立了描写两场相互作用的数学模型.颗粒流的阻力系数采用了前人在液态流化床研究中得出的阻力系数表达形式,通过引入一个修正系数,使其适用于风沙流(气-固两相流).将风沙边界层划分为跃移颗粒所产生的阻力不可忽略的内边界层和跃移颗粒阻力可以忽略但受内边界层影响的外边界层,分别建立了内边界层和外边界层的风速廓线表达式.应用所建立的数学模型,根据由风洞实验测定的跃移风沙流的浓度分布和速度分布资料,计算了跃移风沙流中的风速廓线,并与风洞实测结果进行了对比.结果表明,计算风速廓线与实测风速廓线吻合得比较好,在半对数图上均为上凸的曲线,有别于无风沙运动时的直线.跃移边界层外风速分布可较好地用对数函数来描述.对风沙流中风速廓线的进一步分析证实了风沙物理学奠基人Bagnold在其早期观测风沙流中的风速廓线时提出的"结点现象"(Bagnold结),该结点的高度随风速的增大而升高,随颗粒粒径的增大而降低.根据数值模拟和模拟实验,可以认为有风沙运动的动床剪切风速是综合反映风场与跃移层以及地表之间相互作用的物理量.     

A Homogeneous Langevin Equation Model, Part ii: Simulation of Dispersion in the Convective Boundary Layer

We present a Lagrangian stochastic model of vertical dispersion in the convective boundary layer (CBL). This model is based on a generalized Langevin equation that uses the simplifying assumption that the skewed vertical velocity probability distribution is spatially homogeneous. This approach has been shown to account for two key properties of CBL turbulence associated with large-scale coherent turbulent structures: skewed vertical velocity distributions and long velocity correlation time. A 'linear-skewed' form of the generalized Langevin equation is used, which has a linear (in velocity) deterministic acceleration and a skewed random acceleration. 'Reflection' boundary conditions for selecting a new velocity for a particle that encounters a boundary were investigated, including alternatives to the standard assumption that the magnitudes of the particle incident and reflected velocities are positively correlated. Model simulations were tested using cases for which exact, analytic statistical properties of particle velocity and position are known, i.e., well-mixed spatial and velocity distributions. Simulations of laboratory experiments of CBL dispersion show that (1) the homogeneous linear-skewed Langevin equation model (as well as an alternative 'nonlinear-Gaussian' Langevin equation model) can simulate the important aspects of dispersion in the CBL, and (2) a negatively-correlated-speed reflection boundary condition simulates the observed dispersion of material near the surface in the CBL significantly better than alternative reflection boundary conditions. The homogeneous linear-skewed Langevin equation model has the advantage that it is computationally more efficient than the homogeneous nonlinear-Gaussian Langevin equation model, and considerably more efficient than inhomogeneous Langevin equation models.     

Dynamic Impact of the Vertical Shear of Gradient Wind on the Tropical Cyclone Boundary Layer Wind Field

This work studies the impact of the vertical shear of gradient wind （VSGW） in the free atmosphere on the tropical cyclone boundary layer （TCBL）. A new TCBL model is established, which relies on five- force balance including the pressure gradient force, Coriolis force, centrifugal force, turbulent friction, and inertial deviation force. This model is then employed to idealize tropical cyclones （TCs） produced by DeMaria＇s model, under different VSGW conditions （non-VSGW, positive VSGW, negative VSGW, and VSGW increase/decrease along the radial direction）. The results show that the free-atmosphere VSGW is particularly important to the intensity of TC. For negative VSGW, the total horizontal velocity in the TCBL is somewhat suppressed. However, with the maximum radial inflow displaced upward and outward, the radial velocity notably intensifies. Consequently, the convergence is enhanced throughout the TCBL, giving rise to a stronger vertical pumping at the TCBL top. In contrast, for positive VSGW, the radial inflow is significantly suppressed, even with divergent outflow in the middle-upper TCBL. For varying VSGW along the radial direction, the results indicate that the sign and value of VSGW is more important than its radial distribution, and the negative VSGW induces stronger convergence and Ekman pumping in the TCBL. which favors the formation and intensification of TC.     

Characteristics of Submeso Winds in the Stable Boundary Layer

The characteristics of submeso motions in the stable boundary layer are examined using observations from networks of sonic anemometers with network sizes ranging from a few hundred metres to 100 km. This study examines variations on time scales between 1 min and 1 h. The analysis focuses on the behaviour of the spectra of the horizontal kinetic energy, the ratios of the three velocity variances, their kurtosis, the dependence of horizontal variability on time scale, and the inter-relationship between vertical vorticity, horizontal divergence and deformation. Motions on larger time and space scales in the stable boundary layer are found to be nearly two-dimensional horizontal modes although the ratio of the vorticity to the divergence is generally on the order of one and independent of scale. One exception is a small network where stronger horizontal divergence is forced by a decrease in surface roughness. The horizontal variability, averaged over 1 h, appears to be strongly influenced by surface heterogeneity and increases with wind speed. In contrast, the time dependence of the horizontal structure on time scales less than one hour tends to be independent of wind speed for the present datasets. The spectra of the horizontal kinetic energy and the ratio of the crosswind velocity variance to the along-wind variance vary substantially between networks. This study was unable to isolate the cause of such differences. As a result, the basic behaviour of the submeso motions in the stable boundary layer cannot be generalized into a universal theory, at least not from existing data.     

Statistics of Scalar Fields in the Atmospheric Boundary Layer Based on Large-Eddy Simulations. Part II: Forced Convection

Forced convection in a quasi-steady atmospheric boundary layer is investigated based on a large-eddy simulation (LES) model. The performed simulations show that in the upper portion of the mixed layer the dimensionless (in terms of mixed layer scales) vertical gradients of temperature, humidity, and wind velocity depend on the dimensionless height z/z _i and the Reech number Rn. The peak values of variances and covariances at the top of the mixed layer, scaled in terms of the interfacial scales, are functions of the interfacial Richardson number Ri. As a result expressions for the entrainment rates, in the case when the interfacial layer has a finite depth, and a condition for the presence of moistening or drying regimes in the mixed layer, are derived. Profiles of dimensionless scalar moments in the mixed layer are proposed to be expressed in terms of two empirical similarity functions F _m and F _i , dependent on dimensionless height z/z _i , and the interfacial Richardson number Ri. The obtained similarity expressions adequately approximate the LES profiles of scalar statistics, and properly represent the impact of stability, shear, and entrainment. They are also consistent with the parameterization proposed for free convection in the first part of this paper.     

Coherence and Scale of Vertical Velocity in the Convective Boundary Layer from a Doppler Lidar

We utilized a Doppler lidar to measure integral scale and coherence of vertical velocity w in the daytime convective boundary layer (CBL). The high resolution 2 μm wavelength Doppler lidar developed by the NOAA Environmental Technology Laboratory was used to detect the mean radial velocity of aerosol particles. It operated continuously in the zenith-pointing mode for several days in the summer 1996 during the “Lidars in Flat Terrain” experiment over level farmland in central Illinois. We calculated profiles of w integral scales in both the alongwind and vertical directions from about 390 m height to the CBL top. In the middle of the mixed layer we found, from the ratio of the w integral scale in the vertical to that in the horizontal direction, that the w eddies are squashed by a factor of about 0.65 as compared to what would be the case for isotropic turbulence. Furthermore, there is a significant decrease of the vertical integral scale with height. The integral scale profiles and vertical coherence show that vertical velocity fluctuations in the CBL have a predictable anisotropic structure. We found no significant tilt of the thermal structures with height in the middle part of the CBL.The National Center for Atmospheric Research is sponsored by the National Science Foundation.     

Circulation and boundary layers in differentially heated rotating stratified fluid

Recent laboratory experiments with rotating stratified water in a cylinder have revealed many of the predictions of linearized, analytic theory. Earlier measurements of the velocity field generated in a cylinder by top heating compared well with theory. Large stratification clearly suppressed Ekman pumping so that the interior velocity field (primarily azimuthal) responded by satisfying no-slip top and bottom boundary conditions without the need for Ekman layers. This interior flow also occupied a boundary layer of greater thickness than the Ekman layer under some conditions. Theory and experiments have now been conducted for sidewall heating. As before, experiment and theory agree well over some parameter ranges. But for some parameters, the flow is unstable. The exact nature of the instability remains poorly understood. The size of one combination of both vertical and horizontal boundary layers is governed by the Rossby radius of deformation multiplied by the square root of the Prandtl number. Sidewall boundary layers and their scales will be reviewed with the present results in mind.     

Semigeostrophic Frontal Boundary Layer

The viscous semigeostrophic solutions obtained for the baroclinic Eady wave fronts are analyzed for the generation of the cross-frontal temperature gradient in the boundary layer. In the case of free-slip boundaries, the cross-frontal gradient is maximally generated at the surface by meridional temperature advection. In the case of no-slip boundaries, surface friction reduces the meridional temperature advection in the boundary layer: The maximum generation occurs above the surface layer and the temperature gradient at the surface is maintained by vertical diffusion. The no-slip solution is compared with the Ekman-layer model solution. Errors are quantified for the use of the Ekman-layer model in the mature state of frontogenesis.The surface frontogenesis is found to be affected by diffusivity both directly and indirectly. The direct effect of diffusivity is represented explicitly by the diffusion term in the potential temperature equation. The indirect effect of diffusivity is related implicitly to the temperature advection caused by the viscous part of the ageostrophic motion whose horizontal velocity component is defined by the frictional wind deflection (away from the geostrophy). The direct effect of diffusivity is frontolytical, whilst theindirect effect of diffusivity is frontogenetic in the mesoscale vicinity of the front. The indirect effect of diffusivity contributes dominantly to the mesoscale surface frontogenesis for the free-slip case, but it is offset by the divergence of the dynamic part of the ageostrophic motion at the surface level for the non-slip case.     

Improvement of the K-profile Model for the Planetary Boundary Layer based on Large Eddy Simulation Data

Modifications of the widely used K-profile model of the planetary boundary layer (PBL), reported by Troen and Mahrt (TM) in 1986, are proposed and their effects examined by comparison with large eddy simulation (LES) data. The modifications involve three parts. First, the heat flux from the entrainment at the inversion layer is incorporated into the heat and momentum profiles, and it is used to predict the growth of the PBL directly. Second, profiles of the velocity scale and the Prandtl number in the PBL are proposed, in contrast to the constant values used in the TM model. Finally, non-local mixing of momentum was included. The results from the new PBL model and the original TM model are compared with LES data. The TM model was found to give too high PBL heights in the PBL with strong shear, and too low heights for the convection-dominated PBL, which causes unrealistic heat flux profiles. The new PBL model improves the predictability of the PBL height and produces profiles that are more realistic. Moreover, the new PBL model produces more realistic profiles of potential temperature and velocity. We also investigated how each of these three modifications affects the results, and found that explicit representation of the entrainment rate is the most critical.     

Effect of Roughness on Surface Boundary Conditions for Large-Eddy Simulation

An important parameterization in large-eddy simulations (LESs) of high- Reynolds-number boundary layers, such as the atmospheric boundary layer, is the specification of the surface boundary condition. Typical boundary conditions compute the fluctuating surface shear stress as a function of the resolved (filtered) velocity at the lowest grid points based on similarity theory. However, these approaches are questionable because they use instantaneous (filtered) variables, while similarity theory is only valid for mean quantities. Three of these formulations are implemented in simulations of a neutral atmospheric boundary layer with different aerodynamic surface roughness. Our results show unrealistic influence of surface roughness on the mean profile, variance and spectra of the resolved velocity near the ground, in contradiction of similarity theory. In addition to similarity-based surface boundary conditions, a recent model developed from an a priori experimental study is tested and it is shown to yield more realistic independence of the results to changes in surface roughness. The optimum value of the model parameter found in our simulations matches well the value reported in the a priori wind-tunnel study.