边界层摩擦和地形对斜压波不稳定性的影响

本文用新的边界层顶垂直速度参数化方案研究了当地形和边界层摩擦同时存在时二者对Eady波不稳定性的影响,得到了边界层层结、地面粗糙度、地形坡度的影响,还研究了摩擦和地形对一般化Eady波不稳定性的影响。     

EFFECTS OF BOUNDARY LAYER FRICTION AND TOPOGRAPHY ON THE INSTABILITY OF BAROCLINIC WAVES

In this paper,the simultaneous effects of boundary layer and topography on the instability of Eady waveare investigated by using a new parameterization of the vertical velocity at the top of PBL and the influencesof the stratification of the PBL,roughness and the slope of terrain are shown.Furthermore,the effects of theboundary layer friction and topography on generalized Eady wave are also investigated.     

On the extension of the wind profile over homogeneous terrain beyond the surface boundary layer

Analysis of profiles of meteorological measurements from a 160 m high mast at the National Test Site for wind turbines at Høvsøre (Denmark) and at a 250 m high TV tower at Hamburg (Germany) shows that the wind profile based on surface-layer theory and Monin-Obukhov scaling is valid up to a height of 50–80 m. At higher levels deviations from the measurements progressively occur. For applied use an extension to the wind profile in the surface layer is formulated for the entire boundary layer, with emphasis on the lowest 200–300 m and considering only wind speeds above 3 m s^?1 at 10 m height. The friction velocity is taken to decrease linearly through the boundary layer. The wind profile length scale is composed of three component length scales. In the surface layer the first length scale is taken to increase linearly with height with a stability correction following Monin-Obukhov similarity. Above the surface layer the second length scale (L _MBL ) becomes independent of height but not of stability, and at the top of the boundary layer the third length scale is assumed to be negligible. A simple model for the combined length scale that controls the wind profile and its stability dependence is formulated by inverse summation. Based on these assumptions the wind profile for the entire boundary layer is derived. A parameterization of L _MBL is formulated using the geostrophic drag law, which relates friction velocity and geostrophic wind. The empirical parameterization of the resistance law functions A and B in the geostrophic drag law is uncertain, making it impractical. Therefore an expression for the length scale, L _MBL , for applied use is suggested, based on measurements from the two sites.     

Turbulent diffusivity and diurnal variations in the atmospheric boundary layer

Diurnal variations of the vertical profiles of wind and temperature have been surveyed, and the diffusivity and the dimensionless gradient function in the atmospheric boundary layer have been estimated. Even in the middle of the atmospheric boundary layer (e.g., below a height of 442 m), the vertical wind profile normalized by the surface friction velocity has approximately a universal profile function different from that in the surface boundary layer. Under strong stability conditions, the dimensionless gradient function has a value of about 9.     

Effects of Topographical Slope Angle and Atmospheric Stratification on Surface-Layer Turbulence

The effect of topographical slope angle and atmospheric stratification on turbulence intensities in the unstably stratified surface layer have been parameterized using observations obtained from a three-dimensional sonic anemometer installed at 8 m height above the ground at the Seoul National University (SNU) campus site in Korea for the years 1999–2001. Winds obtained from the sonic anemometer are analyzed according to the mean wind direction, since the topographical slope angle changes significantly along the azimuthal direction. The effects of the topographical slope angle and atmospheric stratification on surface-layer turbulence intensity are examined with these data. It is found that both the friction velocity and the variance for each component of wind normalized by the mean wind speed decrease with increase of the topographical slope angle, having a maximum decreasing rate at very unstable stratification. The decreasing rate of the normalized friction velocity (u _* /U) is found to be much larger than that of the turbulence intensity of each wind component due to the reduction of wind shear with increase in slope angle under unstable stratification. The decreasing rate of the w component of turbulence intensity (σ _w /U) is the smallest over the downslope surface whereas that of the u component (σ _u /U) has a minimum over the upslope surface. Consequently, σ _w /u _* has a maximum increasing rate with increase in slope angle for the downslope wind, whereas σ _u /u _* has its maximum for the upslope wind. The sloping terrain is found to reduce both the friction velocity and turbulence intensity compared with those on a flat surface. However, the reduction of the friction velocity over the sloping terrain is larger than that of the turbulence intensity, thereby enhancing the turbulence intensity normalized by the friction velocity over sloping terrain compared with that over a flat surface.     

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

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

Observations of the Relation Between Upslope Flows and the Convective Boundary Layer in Steep Terrain

Slope flow mechanisms are crucial for the transport of air pollutants in complex terrain. Previous observations in sloping terrain showed upslope flows filling the entire convective boundary layer (CBL) and reducing air pollution concentrations by venting air pollutants out of the CBL into the free atmosphere. During the Pacific 2001 Air Quality Field Study in the Lower Fraser Valley, British Columbia, Canada, we observed slope flows during weak synoptic winds, clear skies, and strong daytime solar heating. With a Doppler sodar we measured the three wind components at the foot of a slope having an average angle of 19° and a ridge height of 780 m. We operated a scanning lidar system and a tethersonde at a nearby site on the adjacent plain to measure backscatter of particulate matter, temperature, wind speed, wind direction, and specific humidity. Strong daytime upslope flows of up to 6 m s⁻¹ through a depth of up to 500 m occurred in the lower CBL, but with often equally strong and deep return flows in the upper part of the CBL. The mass transport of upslope flow and return flow approximately balanced over a 4-h morning period, suggesting a closed slope-flow circulation within the CBL. These observations showed that air pollutants can remain trapped within a CBL rather than being vented from the CBL into the free atmosphere.     

The temperature profile and heat transfer law in a neutrally and stably stratified planetary boundary layer

A logarithmic + polynomial approximation is proposed for the vertical temperature profile in a neutrally or stably stratified planetary boundary layer (PBL) in conditions of quasi-stationarity. Using this approximation with the asymptotic logarithmic + linear law of the Monin-Obukhov similarity theory for the near-surface layer and with the Zilitinkevich formula for the PBL thickness allows one to derive an analytical expression for the function C in the heat transfer law, which permits simple parameterization of the thermal interaction between the atmosphere and the underlying medium in terms of external parameters, such as the geostrophic wind velocity and the temperature difference across the PBL.     

地形上空边界层流中低层锋面结构的理论研究Ⅱ:暖锋、均匀地转流

利用一个有地形、边界层摩擦作用、简化的二层浅水锋面模型,在理论上研究了地形上空边界层流动中地面暖锋的结构及环流分布特征问题。暖锋的坡度主要取决于其暖域地转流、锋面移速,它随锋面移速增大而减小,这与冷锋特征相反。地形对暖锋坡度的影响作用较小。与无地形作用时相比,静止性暖锋冷域中,位于锋面界面附近的闭合正环流系,当暖锋位于地形上游,其伸展范围增大;当暖锋位于迎风坡时,其伸展范围缩小,中心位置上抬;锋面移至背风坡时,其伸展范围重新增大。对于冷域中远离地面暖锋的另一支正环流系来说,当暖锋位于地形上游或迎风坡时,它可被地形完全阻塞于背风侧,地形高度越高,地形阻塞作用越大。在暖锋锋区附近主要存在三支垂直上升运动带：（a）由于边界层摩擦辐合作用,导致在地面暖锋后缘暖区中形成一支水平尺度较小、强度较大的垂直运动带,它随着暖锋移速增大而减弱。该垂直运动带,当暖锋位于地形迎风侧,强度增加;暖锋位于地形背风侧,其强度减弱。（b）在锋区暖域沿锋面存在均匀的上升运动,（c）在冷域远离地面暖锋处,存在一支水平范围较宽,其中心位于边界层顶部附近的垂直运动带,当暖锋位于迎风坡时,这支垂直运动带可被地形阻塞于地形背风侧。     

A ROUTINE FOR THE CALCULATION OF THE TIME-DEPENDENT HEIGHT OF THE ATMOSPHERIC BOUNDARY LAYER FROM SURFACE-LAYER PARAMETERS

A simple routine has been implemented to deduce the 24-hourevolution of the height of the atmospheric boundary layer. This uses a reduced data set of surface-layer parameters, as obtained for examplefrom surface automatic stations.The routine is based on the combination and fitting of the three alreadyexistent models for the evolution of the convectiveboundary layer, the stable boundary layer, and the surface inversionlayer.Hourly values of temperature, friction velocity and potentialtemperature scale (or sensible heat flux) in the surface layer need onlyto be supplied as input data. The lapse rate at the top of the daytime mixed layer is derived fromthe calculated surface inversion profile at sunrise, so that only a roughevaluation of the lapse rate in the free atmosphere remains to be given.The sensitivity of the mixed-layer height is expected to be not verystrong with respect to this last parameter (final part of the growth). The routine has shown satisfactory performances when compared withsodar measurements, working with only a rough average estimate of thefree atmosphere lapse rate.     

Parameterization of Entrainment in a Sheared Convective Boundary Layer Using a First-order Jump Model

Basic entrainment equations applicable to the sheared convective boundary layer (CBL) are derived by assuming an inversion layer with a finite depth, i.e., the first-order jump model. Large-eddy simulation data are used to determine the constants involved in the parameterizations of the entrainment equations. Based on the integrated turbulent kinetic energy budget from surface to the top of the CBL, the resulting entrainment heat flux normalized by surface heat flux is a function of the inversion layer depth, the velocity jumps across the inversion layer, the friction velocity, and the convection velocity. The developed first-order jump model is tested against large-eddy simulation data of two independent cases with different inversion strengths. In both cases, the model reproduces quite reasonably the evolution of the CBL height, virtual potential temperature, and velocity components in the mixed layer and in the inversion layer.The part of this work was done when the first author visited at NCAR.     

Velocity profiles,the resistance law and the dissipation rate of mean flow kinetic energy in a neutrally and stably stratified planetary boundary layer

The logarithmic + polynomial approximation is suggested for vertical profiles of velocity components in a planetary boundary layer (PBL) at neutral and stable stratification. The resistance law functions A and B are determined on the basis of this approximation, using integral relations derived from the momentum equations, the Monin-Obukhov asymptotic formula for the wind profile in a stably stratified near-surface layer and the known expressions for the PBL depth. This result gives a realistic and convenient method for calculating the surface friction velocity and direction and the total dissipation rate of mean flow kinetic energy in terms of geostrophic velocity, buoyancy flux at the surface, the roughness parameter and the Coriolis parameter. In the course of these derivations a review is given of current views on the main problems of the neutral and stable PBL.     

Numerical Modelling of Katabatic Flow Over a Melting Outflow Glacier

A realistic simulation of katabatic flows is not a straightforward task for numerical models. One complicating factor is that katabatic flows develop within a stably stratified boundary layer, which is poorly resolved and described in many numerical models. To capture the jet-shaped shallow flow a model set-up with high vertical resolution is also required. In this study, ‘a state of the art’ mesoscale numerical model is applied in a simulation of katabatic flow over a melting glacier. A basic agreement between observations and model results is found. From scale analysis, it is concluded that the simulated flow can be classified as katabatic. Although the background flow varies in strength and direction, the simulated katabatic flow over Breidamerkurjökull is persistent. Two factors vital for this persistence are identified. First, the melting snow maintains the surface temperature close to 0 °C while the air temperature warms adiabatically as it descends the slope. This provides a ‘self enhanced’ negative buoyancy that drives the flow to a balance with local friction. Second, the jet-like shape of the resulting flow gives rise to a large ‘curvature term’ in the Scorer parameter, which becomes negative in the upper jet. This prevents vertical wave propagation and isolates the katabatic layer of the influence from the free troposphere aloft. Our results suggest that the formation of local microclimates dominated by katabatic flow is a general feature over melting glaciers. The modelled turbulence structure illustrates the importance of non-local processes. Neglecting the vertical transport of turbulence in katabatic flows is not a valid assumption. It is also found that the local friction velocity remains larger than zero through the katabatic jet, due to directional shear where the scalar wind speed approaches its maximum.     

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

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

Turbulence and mixing in the boundary layer generated by shoaling internal waves

Over a range of incidence angles the energy of an internal wave propagating onto a sloping boundary is concentrated in a boundary layer on the slope. As a wave propagates upslope the change in its amplitude and interaction with the downslope flow remaining from previous waves results in the wave breaking and the generation of turbulence and mixing in the boundary layer. Measurements of the overturning and buoyancy scales on the slope show that turbulence is generated and decays during each wave cycle and that much of the energy input to mixing scales is extracted from density inversions generated by the wave-induced mean flow. A comparison with decaying turbulence behind a grid in a stratified water tunnel suggests that the criterion for the extinction of the buoyancy flux is similar in the two cases.     

Influence of soil moisture on simulations of katabatic flow

Summary Simulations of katabatic flow with a two-dimensional dynamic numerical model with a soil parameterization indicate that downslope flow developed over a moist slope is weaker than that over a dry slope. This agrees with earlier findings that daytime anabatic (upslope) flow is weaker over a moister slope. But, whereas the weaker anabatic flow is produced because surface evaporation prevents the moist slope from heating as much as a dry slope, the weaker katabatic flow is produced over moist slopes because (1) the soil thermal conductivity is greater in moist soil, and (2) downward longwave radiation flux from the atmosphere to the surface is greater because of higher humidity in the air near the surface from evaporation. The higher thermal conductivity allows warmer soil temperatures (heat) to diffuse upward to the soil surface and prevents the surface temperature from becoming as cold in the moist run as in the dry run.With 6 Figures     

BAROCLINITY AND NONLINEAR EKMAN PUMPING DYNAMICS

With the Ekman momentum approximation,the influence of atmospheric baroclinity on the dynamics of boundarylayer is studied.Some new results are obtained.These results show that the atmospheric baroclinity plays an importantrole in altering the horizontal velocity of Ekman boundary layer and its angle with the horizontal wind velocity compo-nent near the surface.There are three different physical factors affecting the nonlinear Ekman suction,the vertical mo-tion at the top of boundary layer:first,barotropic geostrophic relative vorticity at the ground;second,the thermal windvorticity induced by the baroclinity;and third,the nonlinear interaction between the barotropic geostrophic relativevorticity and the baroclinic thermal wind vorticity.These results may provide a better physical basis for theparameterization of boundary layer and the interpretation of the numerical modeling results.     

地形辐射效应参数化对海南岛海风环流结构和云水分布模拟的影响

采用非静力中尺度模式WRF研究地形辐射效应参数化对海南岛2012年7月5日多云天气条件下的海风环流结构和云水分布模拟的影响, 并分析了差异产生的可能原因。地形辐射效应是考虑坡地辐射强迫后, 辐射与大气中的各种气体、云以及非均匀下垫面间相互反馈的累积效应。其中地形辐射效应参数化的使用使得温度的模拟更接近实际, 对水汽的模拟也有一定改进能力, 对风速、风向的改进效果不明显。考虑地形辐射效应后, 海风的发生发展演变过程及风场的水平分布无显著变化, 但局地海风以及海风对流云的位置和强度有较明显的改变。山区四周的海风环流结构和对流云的变化与坡地辐射强迫直接相关, 考虑地形辐射效应后, 山坡向阳面的海风有所增强, 背阳面的海风减弱; 向阳坡谷风减弱, 背阳坡谷风增强; 同时紧临海岸的山坡对海风的影响与岛上山坡对谷风的作用类似。平坦地区的海风环流和海风对流云总体上有所减弱, 其变化原因比较复杂。各向海风的强度变化最终会改变海风辐合线的分布, 使海风潜在降水区域发生变化。      

Third-Order Moment Closure Through A Mass-Flux Approach

The parameterization of the third moments, the flux of the heat flux and the flux of the potential temperature variance, is considered. It is shown that present parameterizations of these moments using the mass-flux approach with a `top-hat' profile assumption lead to a significant underestimation, resulting in an inaccurate representation of second moments in the convective boundary layer. It is also shown that the underestimation is a result of the `top-hat' profileassumption in which the sub-plume contributions to the total fluxes are ignored. By including these contributions a new parameterization is proposed, whichsatisfies the physical requirements of symmetry andrealizability, and gives results that are in fair agreement with thelarge-eddy simulation data.