一个诊断非平坦地形上边界层风的数值模式

根据半地转大气边界层模式，由大尺度数值模式并考虑了下垫面地形及粗糙度的水平非均匀性及大尺度气压场的时空变化，给出了一个诊断边界层风的数值模式。对低纬度运用塔层风模式进行诊断。诊断结果与实测资料比较，风向风速均达到了一定的精确度。     

The numerical experiments on applicating geostrophic momentum approximation to the baroclinic and non-neutral PBL

In this paper, Wu and Blumen’s boundary layer geostrophic momentum approximation model (Wu and Blumen, 1982) is applied to baroclinic and non-neutral PBL, the motion equations for the PBL under the geostrophic momentum approximation are solved, in which the eddy transfer coefficient is a function of the distributions of the wind and temperature. The results are compared with those in barotropic and neutral conditions with the geostrophic momentum approximation. It is found that in the baroclinic condition, the wind distribution has both the characteristics of a steady, homogeneous and baroclinic PBL and those caused by the geostrophic momentum approximation. Those in non-neutral conditions show that they retain the intrinsic characteristics for the wind in non-neutral PBL, at the same time, the effects of the large-scale advection and local variation are also included. We can predict the wind in the non-neutral and baroclinic PBL by use of the geostrophic mo-mentum approximation when the temporal and spatial distributions of the geostrophic wind, as well as the po-tential temperatures and their variation rates at the upper and lower boundary of the PBL are given by large-scale model. Finally, the model is extended to the case over sea surface.     

A barotropic planetary boundary layer

The temperature and wind profiles in the planetary boundary layer (PBL) are investigated. Assuming stationary and homogeneous conditions, the turbulent state in the PBL is uniquely determined by the external Rossby number and the stratification parameters. In this study, a simple two-layer barotropic model is proposed. It consists of a surface (SL) and overlying Ekman-type layer. The system of dynamic and heat transfer equations is closed usingK theory. In the SL, the turbulent exchange coefficient is consistent with the results of similarity theory while in the Ekman layer, it is constant. Analytical solutions for the wind and temperature profiles in the PBL are obtained. The SL and thermal PBL heights are properly chosen functions of the stratification so that from the solutions for wind and temperature, the PBL resistance laws can be easily deduced. The internal PBL characteristics necessary for the calculation (friction velocity, angle between surface and geostrophic winds and internal stratification parameter) are presented in terms of the external parameters. Favorable agreement with experimental data and model results is demonstrated. The simplicity of the model allows it to be incorporated in large-scale weather prediction models as well as in the solution of various other meteorological problems.     

地转动量近似下边界层中风场的调整

本文利用地转动量近似,并假设气压场为定常的圆形涡旋和初始风场不满足四力平衡(气压梯度力、科里奥利力、湍流粘性力和半地转惯性力)的条件下,求解了正压边界层中风场向气压场调整的初边值问题,得到了一些初步结论。本工作为利用四力平衡下的风速分布来诊断预报边界层风场提供了理论依据。     

Frequency of Boundary-Layer-Top Fluctuations in Convective and Stable Conditions Using Laser Remote Sensing

The planetary boundary-layer (PBL) height is determined with high temporal and altitude resolution from lidar backscatter profiles. Then, the frequencies of daytime thermal updrafts and downdrafts and of nighttime gravity waves are obtained applying a fast Fourier transform on the temporal fluctuation of the PBL height. The principal frequency components of each spectrum are related to the dominant processes occurring at the daytime and nighttime PBL top. Two groups of cases are selected for the study: one group combines daytime cases, measured in weak horizontal wind conditions and dominated by convection. The cases show higher updraft and downdraft frequencies for the shallow, convective boundary layer and lower frequencies for a deep PBL. For cases characterized by strong horizontal winds, the frequencies directly depend on the wind speed. The temporal variation of the PBL height is determined also in the likely presence of lee waves. For nighttime cases, the main frequency components in the spectra do not show a real correlation with the nocturnal PBL height. Altitude fluctuations of the top of the nocturnal boundary layer are observed even though the boundary layer is statically stable. These oscillations are associated with the wind shear effect and with buoyancy waves at the PBL top.     

The wind field in the nonlinear multilayer planetary boundary layer

By use of the small parameter expansion method, the nonlinear planetary boundary layer (PBL) is studied in this paper. The PBL is divided into the surface layer and the Ekman layer, which is divided into several sublayers. In the surface-layer, the eddy coefficient K is taken as a linear function of height; in the Ekman layer, different constant K values are taken within different sublayers: these values are determined from O'Brien's formula (O'Brien, 1970) approximately. Under the upper and lower boundary conditions and the continuity conditions of the wind velocities and turbulent stresses at each boundary between sublayers, analytical expressions for wind velocity in all sublayers and the vertical velocity at the top of the PBL are obtained. A specific example of steady axisymmetrical circular high and low pressure areas is analysed, and some new conclusions are obtained. The results are in better agreement with reality than previous results. This example also shows that the vertical velocity at the top of the PBL caused by friction approaches zero near the center of a high or low pressure system for this model, but attains its maximum absolute values near the center of the high or low pressure area for Wu's (1984) model. This is due to the fact that in our model, the geostrophic wind speed near the center of this specific vortex approaches zero, which causes the wind shear and the friction effect to be very weak. Therefore the wind distribution in the PBL is very sensitive to the type of eddy coefficient.     

An ice breeze mechanism for boundary-layer jets

The existence of a low-level (z=~1000 m) jet adjacent to a sea-ice boundary is investigated with a two-dimensional numerical model. A thermally-direct ice breeze circulation is induced by specifying an ice-sea surface temperature gradient, with the mean geostrophic wind parallel to the ice edge. Pressure changes associated with over-water mixed-layer development create an increase in geostrophic velocity that accounts for most of the increase in wind speed. A change in initial geostrophic wind direction has significant effects on location and intensity of the low-level jet; geostrophic winds parallel to the ice edge result in stronger jets than occur with cross-ice geostrophic winds. An inertial oscillation simulated by the model in 1-D makes a negligible contribution to the low-level jet.     

Effects of local accelerations and baroclinity on the mean structure of the atmospheric boundary layer over the sea

Substitution of the geostrophic wind by the actual upper wind in the equations of motion for the boundary layer implies less sensitivity of the mean wind to inertial effects. This is confirmed by observations, although the problem of computing time or spatial derivatives from scattered data reduces the accuracy and the clarity of the results. It is found that acceleration (deceleration) increases (decreases) the cross-isobar angle whereas the geostrophic drag coefficient is a minimum (maximum) for crosswind acceleration (deceleration). On the other hand, cold air advection increases the cross-isobar angle whereas the geostrophic drag coefficient is a maximum when the thermal wind is parallel to the surface wind. The universal functions A _m and B _m based on vertically averaged winds are also rather insensitive to inertial influences.     

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.     

边界层动力学中的Ekman动量近似

自由大气中,大气运动的基本状态是地转风,近年来发展的地转动量近似,是为了进一步研究非均匀地转流的动力学问题,然而,在边界层大气中,运动的基本状态是经典的Ekman流,所以对边界层运动来说,地转动量近似是不合适的,需作一推广。本文提出了一种所谓Ekman动量近似,它相似于自由大气中的地转动量近似,并讨论了Ekman动量近似的物理基础,对边界层的风场结构及边界屋顶部的垂直速度也作了详细分析。     

Remote Sensing Of Three-Dimensional Winds With Elastic Lidar: Explanation Of Maximum Cross-Correlation Method

Maximum cross-correlation provides a method to remotely determine highly resolved three-dimensional fields of horizontal winds with elastic lidar throughout large volumes of the planetary boundary layer (PBL). This paper details the technique and shows comparisons between elastic lidar winds, remotely sensed laser Doppler velocimeter (LDV) wind profiles, and radiosonde winds. Radiosonde wind data were acquired at Barcelona, Spain, during the Barcelona Air-Quality Initiative (1992), and the LDV wind data were acquired at Sunland Park, New Mexico (N.M.), during the Border Area Air-Quality Study (1994). Comparisons show good agreement between the different instruments, and demonstrate the method useful for air pollution management at the local/regional scale. Elastic lidar winds could thus offer insight into aerosol and pollution transport within the PBL. Lidar wind fields might also be used to nudge or improve initialization and evaluation of atmospheric meteorological models.     

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

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

On the influence of variations in Monin-Obukhov length on horizontal roll vortices in an inversion-capped planetary boundary layer

Observations of the planetary boundary layer under conditions with strong or moderate winds often show the large-scale boundary-layer motions to be highly elongated in a direction close to that of the geostrophic wind. The properties of such large-scale motions are examined by means of a two-dimensional numerical model. The small-scale turbulence is parameterized using a buoyancy-dependent mixing-length hypothesis. The objective is to understand how the properties and dynamics of such rolls depend upon the relative importance of shear and buoyancy forces.     

Steady state hodographs in a baroclinic boundary layer

If the steady-state geostrophic wind vector varies exponentially with height in the planetary boundary layer, calculated hodographs of the 24-hour mean wind resemble the classic Ekman spiral distorted by thermal effects. For such an assumed distribution, Lettau's (1967) conclusion that the effects of thermal influence on the steady-state boundary-layer winds can be linearly superimposed on the effects of internal friction is justified.The minimum value of the cross-isobar angle of the surface wind for a given magnitude of the thermal wind vector occurs when that vector points about 345° to the right of the surface geostrophic wind vector and the maximum value occurs when it points at an azimuth of about 120° relative to the surface geostrophic wind vector. The range of values of the cross-isobar angle is almost directly proportional to the magnitude of the thermal wind vector.Hodographs resulting from the assumed variation of the geostrophic wind have approximately the same shape as 24-hour mean hodographs at two locations over the Great Plains.This work is part of a thesis submitted to the University of Wisconsin in partial fulfillment of the requirements for the Ph.D. Degree, written under the supervision of Professor H. Lettau, Department of Meteorology.     

Variability of surface wind directions over Finnmark, Norway, and coupling to the larger-scale atmospheric circulation

This study documents the temporal and spatial variability of surface wind conditions over the Norwegian county of Finnmark and the coupling of local surface winds to the larger-scale atmospheric circulation, represented by the mean sea level pressure field. At locations along the northern coast, thermally driven offshore winds from the south dominate, especially during the cold season. During the warm season, downward mixing of westerly overlying winds becomes more important as the stability of the boundary layer stratification decreases. In the western part, locations are situated in valleys, resulting in two opposing along-channel dominant wind directions. Seasonal changes in atmospheric conditions are reflected in a reversal of the dominant wind direction, with a component towards the coast during the cold season. At all locations, a clear separation between different prevailing surface wind directions in each season can be achieved based exclusively on local mean sea level geostrophic wind direction. This allows statistical downscaling of the prevailing surface wind conditions from lower-resolution simulations of the surface pressure field and may improve local wind forecasts over complex terrain.     

大气边界层的初值形成

采用850毫巴和地面的大尺度天气状况和物理特性作为外参数,将大气边界层下层的地表层进行特殊处理,考虑了层结的作用,其上分成等距15层,考虑地转民随高度线性变化。上下层耦合,建立一个定常平衡初值模型,并将所得的初值作为第一级近似,再利用地面、300、600、900、1500米五个高度测风资料进行动力调整,试图建立一个更加接近实际状态的边界层初值。 文中利用所建立的模型,对大气边界层的一些物理参数的作用进行了分析,并对动力调整过程以及垂直分层的选取等作了讨论。 最后例举华北地区1966年8月23日08时边界层的各层次风场和垂直运动,清楚地显示出边界层内的细致结构。     

边界层风场调整的分析

分析任意气压场在中性、非线性正压边界层中风场的调整问题。结果表明:任意初始风场,将向给定的气压调整。具休调整过程是通过绝对涡度振荡及粘性系数的耗散,其调整的时间尺度为[O(t^(-3/2))],比自由大气地转适应时间尺度[O(t^(-1))]快。     

Wind energy estimates by use of a diagnostic model

A diagnostic model is a relatively simple and practical tool for modeling the wind flow of the boundary layer in complex terrain. The model begins with a wind analysis based on available surface wind reports and geostrophic winds (computed from pressure data). The height of the boundary layer top (upper surface of the computational domain) is prescribed to fit local conditions. Using the continuity equation in terrain-following coordinates, the winds at mesh points are adjusted to produce nondivergence while maintaining the original vertical component of vorticity. The method of computing the nondivergent winds uses direct alterations. This method may be useful for other modeling purposes and will be described. Data for a long period (usually a year) are analyzed to obtain eigenvectors and the associated time series of their coefficients at each observation time. The model is run only for the five or six eigenvectors that explain most of the variance. The wind field at any particular time is reconstructed from the eigenvector solutions and their appropriate coefficients. Comparisons of model results with measured winds at sites representing different types of terrain will be shown. The accuracy and economy of the model make it a useful tool for estimating wind energy and also for giving wind fields for low-level diffusion models.     

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.