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.     

Surface energy balance,parameterizations of boundary-layer heights and the application of resistance laws near an Antarctic Ice Shelf front

A study of the surface energy balance with turbulent fluxes obtained by the Monin-Obukhov similarity theory and a comparison with results for resistance laws are presented for the strong baroclinic conditions in the vicinity of the Filchner/Ronne Ice Shelf front. The data are taken from a field experiment in the Antarctic summer season 1983/84. For the first time in the coastal Antarctic region, this data set comprises synchronous energy balance measurements over the polynya and the ice shelf together with soundings of the boundary layer, yielding vertical profiles of the wind velocity and temperature over the ice shelf, at the ice shelf front and over the polynya.Over the ice shelf, the radiation balance is the largest component of the energy fluxes and is mainly compensated by the subsurface energy flux and the turbulent heat flux in the daily mean. Over the polynya, turbulent fluxes of sensible and latent heat lead to large energy losses of the water surface in the night-time and in situations of very low air temperatures.Different parameterizations for boundary-layer height are compared using tethered sonde and energy balance measurements. With the height of the inversion base over the polynya and the height of the critical bulk Richardson number over the ice shelf, external parameters for the application of resistance laws were determined. The comparison of turbulent surface fluxes obtained by the energy balance measurements and by the resistance laws shows good agreement for the convective conditions over the polynya. For the stably stratified boundary layer over the ice shelf with small amounts of the turbulent heat flux, the deviation is large for the case of a cold air outflow with a superposed inertial oscillation.     

农林复合带一维非静力大气边界层模式及其模拟分析

建立了一个农林复合带地区一维非静力大气边界层能量闭合模式，对1000m以下的大气边界层内的风、温、湿作了24h的预报，并对下垫面3种不同参数化方案(农作物、森林、无植被)的输出结果与实测值进行了分析和比较；同时通过敏感性试验，突出比较了农作物和森林下垫面对大气边界层垂直流场，湍流垂直交换和湍能的影响。结果表明，本模式能改善边界层风速、位温和湿度预报的模拟效果，下垫面植被对边界层气象要素大小和分布有显著的作用，对湍能垂直分布有一定影响。     

Vertical Heat Transfer in the Lower Atmosphere over the Arctic Ocean During Clear-sky Periods

A diagnostic study of heat transfer within the lower atmosphere and between the atmosphere and the surface of the Arctic Ocean snow/ice pack during clear-sky conditions is conducted using data from the Surface Heat Budget of the Arctic Ocean (SHEBA) field experiment. Surface heat budgets computed for four cloudy and four clear periods show that, while the net turbulent heat fluxes at the surface are small during the cloudy periods, during the clear-sky periods they are a considerable source of surface heating, balancing significant portions of the conductive heat fluxes from within the snow/ice pack. Analysis of the dynamics and thermodynamics of the lower atmosphere during the clear-sky periods reveals that a considerable portion of the heat lost to the surface by turbulent heat fluxes is balanced by locally strong heating near the atmospheric boundary-layer (ABL) top due to the interaction of subsiding motions with the strong overlying temperature inversions surmounting the ABL. This heat is then entrained into the ABL and transported to the surface by turbulent mixing, maintained by a combination of vertical wind shear and wave-turbulence interactions. The frequency of stable, clear-sky periods, particularly during the winter, combined with these results, suggests that the downward transfer of heat through the lower atmosphere and into the surface represents an important component of the heat budgets of the lower atmosphere and snow/ice pack over the annual cycle     

E − ε − 〈θ 2〉 turbulence closure model for an atmospheric boundary layer including the urban canopy

A modified three-parameter model of turbulence for a thermally stratified atmospheric boundary layer (ABL) is presented. The model is based on tensor-invariant parameterizations for the pressure–strain and pressure–temperature correlations that are more complete than the parameterizations used in the Mellor–Yamada model of level 3.0. The turbulent momentum and heat fluxes are calculated with explicit algebraic models obtained with the aid of symbol algebra from the transport equations for momentum and heat fluxes in the approximation of weakly equilibrium turbulence. The turbulent transport of heat and momentum fluxes is assumed to be negligibly small in this approximation. The three-parameter $E - \varepsilon - {\left\langle {\theta ^{2} } \right\rangle }$ model of thermally stratified turbulence is employed to obtain closed-form algebraic expressions for the fluxes. A computational test of a 24-h ABL evolution is implemented for an idealized two-dimensional region. Comparison of the computed results with the available observational data and other numerical models shows that the proposed model is able to reproduce both the most important structural features of the turbulence in an urban canopy layer near the urbanized ABL surface and the effect of urban roughness on a global structure of the fields of wind and temperature over a city. The results of the computational test for the new model indicate that the motion of air in the urban canopy layer is strongly influenced by mechanical factors (buildings) and thermal stratification.     

Dynamical Initialization of Boundary-Layer Data

A procedure for the dynamic initialization of wind and temperature fields within the atmospheric boundary layer (ABL) is presented. The procedure uses transilient turbulence theory to assess the turbulent fluxes from observed wind and temperature data, which are used to estimate the eddy transfer coefficients, so enabling a local closure ABL model to be integrated forward in time to a new steady state. The method has been applied to initialize kytoon data taken at Kharagpur (22.3° N, 87.2° E) during the MONTBLEX field programme. Results of a case study for 17 June 1990 are discussed.     

An Actively Controlled Wind Tunnel And Its Application To The Reproduction Of The Atmospheric Boundary Layer

An actively controlled wind tunnel equipped with multiple fansand airfoils has been developed, mainly for the purpose of reproducing the atmospheric boundary layer (ABL) for wind engineering applications. Various fluctuating flows can be achieved in this wind tunnel by altering the input data of the fans and airfoils through computer control. In this study, the ABL is physically simulated in this wind tunnel, and particular attention ispaid to the simulation of the profile of Reynolds stress. The method of generating the fluctuating flow and the experimental results of reproducing the ABL are presented. As the results show, the spatial distribution of Reynolds stress is satisfactorily simulated, and the profiles of other statisticalturbulent parameters, such as mean velocity, turbulent intensity, integral scale and power spectrum are successfully reproduced simultaneously.     

On the consideration of mesoscale transports in climate modelling

Summary The dynamical effect of land surface heterogeneity on heat fluxes in the atmospheric boundary layer (ABL) is investigated using numerical simulations with a non-hydrostatic model over a wide range of grid resolutions. It is commonly assumed that mesoscale or dynamical fluxes associated with mesoscale and convective circulations simulated by a high-resolution model (subgrid (SG) model) on the subgrid scale of a climate model (large-scale (LS) model) represent additional processes in the ABL, which are not considered by the turbulence scheme of the LS-model, and which can be parameterized using the SG-model. The present study investigates the usefulness of this methodology for small-scale and large-scale idealized heterogeneities using a SG-model resolving mesoscale or even microscale circulations to compute the mesoscale fluxes on the scale of the LS-model. It is shown that the dynamical transports as derived from the SG-model should not be used to correct the parameterized turbulent fluxes of the LS-model. The reason is that the subgrid circulations simulated by the SG-model interact with the fields of wind and scalars in the ABL, which results in reduced turbulent fluxes in the ABL. Thus the methodology of previous studies to use mesoscale/dynamical fluxes for the correction of flux profiles simulated by climate models seems to be questionable.     

A bulk model for the atmospheric planetary boundary layer

The integrated momentum and thermodynamic equations through the planetary boundary layer (PBL) are solved numerically to predict the mean changes of wind and potential temperature from which surface fluxes are computed using bulk transfer coefficients of momentum and heat. The second part of the study involves a formulation and testing of a PBL height model based on the turbulent energy budget equation where turbulent fluxes of wind and heat are considered as the source of energy. The model exhibits capability of predicting the PBL height development for both stable and unstable regimes of observed conditions. Results of the model agree favourably with those of Deardorff's (1974a) and Tennekes' (1973) models in convective conditions.Contribution number 396.     

北京北郊冬季大风过程湍流通量演变特征的分析研究

利用中国科学院大气物理研究所325 m气象观测塔1993年12月～1994年1月大气边界层实验资料, 计算分析了大风过境过程中47 m和120 m高度湍流通量演变特征及其影响因子, 以及与风速、 稳定度等参数的关系。结果表明: 大风过程对近地面层的物质能量输送有着重要影响, 大风之前出现短时间动量上传和热量下传; 大风过程中的湍流通量数值明显高于过境后, 水平方向湍流通量数值和能量增加幅度大于垂直方向; 当风速大于临界值5 m/s时, 湍流通量与风速、 湍流动能的相关迅速增大; 湍流谱特征表现为湍流能量的低频部分增加、 湍流谱曲线变宽; 大风能强烈影响近地面层的能量收支。     

Current problems in the stratocumulus-topped atmospheric boundary layer

Extended sheets of stratocumulus (Sc) in the upper part of the atmospheric boundary layer (ABL) often occur under appropriate meteorological conditions. These cloud decks are important both in climate studies and in weather forecasting. We review the current knowledge of the turbulent structure of the ABL capped by a cloud deck, in the light of recent observations and model studies. The most important physical processes determining this structure are longwave radiative cooling at cloud top, shortwave radiative wanning by absorption in the cloud, surface buoyancy flux, and wind shear in the ABL. As a result, turbulence can cause entrainment against the buoyancy jump at cloud top. In cases where only longwave radiative fluxes and surface buoyancy fluxes are important, the turbulent structure is relatively well understood. When shortwave radiative fluxes and/or wind shear are also important, the resulting turbulent structure may change considerably. A decoupling of the cloud from the sub-cloud layer or of the top of the cloud from the rest of the ABL is then regularly observed. In no cases are the details of the entrainment at cloud top understood well enough to derive a relatively simple formulation that is consistent with observations. Cloud-top entrainment instability may lead to the break-up of a cloud deck (but also to cloud deepening). The role of mesoscale circulations in determining fractional cloudiness is not yet well understood.     

北京地区一次强沙尘暴过程的大气边界层结构和湍流通量输送特征

利用中国科学院大气物理研究所北京325 m气象塔的风速和温度平均场观测资料和湍流资料,以及北京市气象台地面常规气象资料和逐日08:00和20:00(北京时间)的探空资料,分析了2002年3月18～22日沙尘暴过境前后北京城市边界层结构特征和湍流输送特征,结果表明:1)在沙尘暴爆发前,边界层中水平风速一直较小;气温较高,大气层结稳定,在边界层上部有强大的逆温层.随着冷锋过境,沙尘暴爆发,边界层中水平风速和平均湍流速度急剧增强;温度也突然变化,先迅速增强后又持续下降,逆温层迅速被破坏.2)沙尘暴初期,280 m上为系统性上升气流,而47和120 m则为系统性的下沉气流.随着沙尘暴爆发,湍流动能、向下传输的动量以及向上传输的感热也迅速增大,并且120 m高度的湍能、动量通量以及感热通量明显高于47和280 m,这与北京的局地环流有关.3)本次沙尘暴过程中,120和47 m层的摩擦速度都明显超过了北京的临界摩擦速度,表明局地起沙也是本次沙尘暴过程中北京沙尘的一个重要沙源.     

Application of transilient turbulent theory to study interactions between the atmospheric boundary layer and forest canopies

The new Forest-Land-Atmosphere ModEl called FLAME is presented. The first-order, nonlocal turbulence closure called transilient turbulence theory (Stull, 1993) is applied to study the interactions between a forested land-surface and the atmospheric boundary layer (ABL). The transilient scheme is used for unequal vertical grid spacing and includes the effects of drag, wake turbulence, and interference to vertical mixing by plant elements. Radiation transfer within the vegetation and the equations for the energy balance at the leaf surface have been taken from Norman (1979). Among others, the model predicts profiles of air temperature, humidity and wind velocity within the ABL, sensible and latent heat fluxes from the soil and the vegetation, the stomata and aerodynamic resistances, as well as profiles of temperature and water content in the soil. Preliminary studies carried out for a cloud free day and idealized initial conditions are presented. The canopy height is 30 m within a vertical domain of 3 km. The model is able to capture some of the effects usually observed within and above forested areas, including the relative wind speed maximum in the trunk space and the counter gradient-fluxes in the lower part of the plant stand. Of special interest is the determination of the location and magnitude of the turbulent mixing between model layers, which permits one to identify the effects of large eddies transporting momentum and scalar quantities into the canopy. A comparison between model simulations and field measurements will be presented in a future paper.     

Modelling the One-Dimensional Stable Boundary Layer with an E−ℓ Turbulence Closure Scheme

The atmospheric boundary layer (ABL) model of Weng and Taylor with E−ℓ turbulence closure is applied to simulate the one-dimensional stably stratified ABL. The model has been run for nine hours from specified initial wind, potential temperature and turbulent kinetic energy profiles, and with a specified cooling rate applied at the surface. Different runs are conducted for different cooling rates, geostrophic winds and surface roughnesses. The results are discussed and compared with other models, large-eddy simulations and published field data.     

Observations of atmospheric boundary layer evolution above the Gulf Stream frontal zone

Ship borne measurements of atmospheric boundary layer (ABL) parameters, sea-surface temperature and radar signals are analyzed to reveal the effects of the ABL transformation above the Gulf Stream temperature frontal zone. It was found that local changes in vertical gradients of wind speed and air temperature are well correlated with sub-mesoscale (~ 10 km) sea surface temperature variations. These effects are accompanied by appropriate variations in surface wind stresses that were identified from microwave backscatter.For steady atmospheric conditions the same effects were observed on spatial scales of 100 km, demonstrating positive radar signal contrast of the Gulf Stream warm waters with respect to surrounding Sargasso sea and shelf water areas. A simplified model of the ABL, accounting for an effect of spatial inhomogeneity by introducing an internal boundary layer, is used to analyze field observations. The model is able to reproduce both sub-mesoscale and mesoscale ABL evolution.     

UNCERTAINTY IN THE SPECIFICATION OF SURFACE CHARACTERISTICS: A STUDY OF PREDICTION ERRORS IN THE BOUNDARY LAYER

The effects of uncertainty in the specification of surface characteristics on simulated atmospheric boundary layer (ABL) processes and structure were investigated using a one-dimensional soil-vegetation-boundary layer model. Observational data from the First International Satellite Land Surface Climatology Project Field Experiment were selected to quantify prediction errors in simulated boundary-layer parameters. Several numerical 12-hour simulations were performed to simulate the convective boundary-layer structure, starting at 0700 LT 6 June 1987.In the control simulation, measured surface parameters and atmospheric data were used to simulate observed boundary-layer processes. In the remaining simulations, five surface parameters – soil texture, initial soil moisture, minimum stomatal resistance, leaf area index, and vegetation cover – were varied systematically to study how uncertainty in the specification of these surface parameters affects simulated boundary-layer processes.The simulated uncertainty in the specification of these five surface parameters resulted in a wide range of errors in the prediction of turbulent fluxes, mean thermodynamic structure, and the depth of the ABL. Under certain conditions uncertainty in the specifications of soil texture and minimum stomatal resistance had the greatest influence on the boundary-layer structure. A lesser but still moderately strong effect on the simulated ABL resulted from (1) a small decrease (4%) in the observed initial soil moisture (although a large increase [40%] had only a marginal effect), and (2) a large reduction (66%) in the observed vegetation cover. High uncertainty in the specification of leaf area index had only a marginal impact on the simulated ABL. It was also found that the variations in these five surface parameters had a negligible effect on the simulated horizontal wind fields. On the other hand, these variations had a significant effect on the vertical distribution of turbulent heat fluxes, and on the predicted maximum boundary-layer depth, which varied from about 1400–2300 m across the 11 simulations. Thus, uncertainties in the specification of surface parameters can significantly affect the simulated boundary-layer structure in terms of meteorological and air quality model predictions.     

HUBEX试验区近地面层的湍流输送

1998年淮河流域能量和水循环试验(HUBEX)期间曾进行了1个月的近地面层湍流观测.分析不稳定条件下湍流的统计特性和谱特征,并与Monin-Obukhov相似理论进行了比较.结果表明,不稳定的时候各湍流量的统计特征与相似理论的预期相符.虽然不稳定条件下温度和湿度涨落的相关系数很高,谱的式样也相近,但湿度谱的峰值频率高于温度谱.协谱曲线的形状显示感热通量的谱峰较宽,表现出w和T在较宽范围的强相关性,而水汽通量谱在高频段下降很快,说明水汽的输送更多地出现在低频部分.从谱相关系数可见,在近中性的时候,各尺度湍流涡的热量输送效率普遍较低,随着不稳定性增强而显著提高.分析还发现,不论不稳定性的程度如何,小尺度湍流的水汽输送效率都较低.水汽通量谱的相关系数随稳定度的变化不如热通量的谱相关系数大,表明近中性时除小尺度湍流外其他湍流涡的水汽输送效率高于热量输送.     

Equilibrium Atmospheric Boundary-Layer Flow: Computational Fluid Dynamics Simulation with Balanced Forces

Forcing relationships in steady, neutrally stratified atmospheric boundary-layer (ABL) flow are thoroughly analyzed. The ABL flow can be viewed as balanced between a forcing and a drag term. The drag term results from turbulent stress divergence, and above the ABL, both the drag and the forcing terms vanish. In computational wind engineering applications, the ABL flow is simulated not by directly specifying a forcing term in the ABL but by specifying boundary conditions for the simulation domain. Usually, these include the inflow boundary and the top boundary conditions. This ‘boundary-driven’ ABL flow is dynamically different from its real counterpart, and this is the major reason that the simulated boundary-driven ABL flow does not maintain horizontal homogeneity. Here, first a dynamical approach is proposed to develop a neutrally stratified equilibrium ABL flow. Computational fluid dynamics (CFD) software (Fluent 6.3) with the standard \(k\) – \(\varepsilon \) turbulence model is employed, and by applying a driving force profile, steady equilibrium ABL flows are simulated by the model. Profiles of wind speed and turbulent kinetic energy (TKE) derived using this approach are reasonable in comparison with the conventional logarithmic law and with observational data respectively. Secondly, the equilibrium ABL profiles apply as inflow conditions to simulate the boundary-driven ABL flow. Simulated properties between the inlet and the outlet sections across a fetch of 10 km are compared. Although profiles of wind speed, TKE, and its dissipation rate are consistently satisfactory under higher wind conditions, a deviation of TKE and its dissipation rate between the inlet and outlet are apparent (7–8 %) under lower wind-speed conditions (2 m s \(^{-1}\) at 10 m). Furthermore, the simulated surface stress systematically decreases in the downwind direction. A redistribution of the pressure field is also found in the simulation domain, which provides a different driving pattern from the realistic case in the ABL.     

南沙海域近海层大气湍流结构及输送特征研究

本文根据南沙群岛渚碧礁的一次大气湍流观测资料，分析研究南沙海域的湍流热通量输送、方差相似性、风速各分量谱、温度谱及各湍流通量协谱等大气湍流结构和输送特征，并与陆地下垫面上的结果作了比较，得到了一些有关南沙海域大气湍流输送及湍流结构的新认识。     

Observations of Coherent Turbulence Structures in the Near-Neutral Atmospheric Boundary Layer

Turbulence structures of high Reynolds number flow in the near-neutral atmospheric boundary layer (ABL) are investigated based on observations at Shionomisaki and Shigaraki, Japan. A Doppler sodar measured the vertical profiles of winds in the ABL. Using the integral wavelet transform for the time series of surface wind data, the pattern of a descending high-speed structure with large vertical extent (from the surface to more than 200-m level) is depicted from the Doppler sodar data. Essentially this structure is a specific type of coherent structure that has been previously shown in experiments on turbulent boundary-layer flows. Large-scale high-speed structures in the ABL are extracted using a long time scale (240 s) for the wavelet transform. The non-dimensional interval of time between structures is evaluated as 3.0–6.2 in most cases. These structures make a large contribution to downward momentum transfer in the surface layer. Quadrant analyses of the turbulent motion measured by the sonic anemometer (20-m height) suggest that the sweep motion (high-speed downward motion) plays a substantial role in the downward momentum transfer. In general, the contribution of sweep motions to the momentum flux is nearly equal to that of ejection motions (low-speed upward motions). This contribution of sweep motions is related to the large-scale high-speed structures.