Determination of the Atmospheric Boundary Layer Height from Radiosonde and Lidar Backscatter

The height of the atmospheric boundary layer is derived with the help of two different measuring systems and methods. From radiosoundings the boundary layer height is determined by the parcel method and by temperature and humidity gradients. From lidar backscatter measurements a combination of the averaging variance method and the high-resolution gradient method is used to determine boundary layer heights. In this paper lidar-derived boundary layer heights on a 10 min basis are presented. Datasets from four experiments – two over land and two over the sea – are used to compare boundary layer heights from both methods. Only the daytime boundary layer is investigated because the height of the nighttime stable boundary layer is below the range of the lidar. In many situations the boundary layer heights from both systems coincide within ±200 m. This corresponds to the standard deviation of lidar-derived 10-min values within a 1-h interval and is due to the time and space variability of the boundary layer height. Deviations appear for certain situations and depend on which radiosonde method is applied. The parcel method fails over land surfaces in the afternoon when the boundary layer stabilizes and over the ocean when the boundary layer is slightly stable. An automatic radiosonde gradient method sometimes fails when multiple layers are present, e.g. a residual layer above the growing convective boundary layer. The lidar method has the advantage of continuous tracing and thus avoids confusion with elevated layers. On the other hand, it mostly fails in situations with boundary layer clouds     

Comparing mixing-length models of the diabatic wind profile over homogeneous terrain

Models of the diabatic wind profile over homogeneous terrain for the entire atmospheric boundary layer are developed using mixing-length theory and are compared to wind speed observations up to 300 m at the National Test Station for Wind Turbines at Høvsøre, Denmark. The measurements are performed within a wide range of atmospheric stability conditions, which allows a comparison of the models with the average wind profile computed in seven stability classes, showing a better agreement than compared to the traditional surface-layer wind profile. The wind profile is measured by combining cup anemometer and lidar observations, showing good agreement at the overlapping heights. The height of the boundary layer, a parameter required for the wind profile models, is estimated under neutral and stable conditions using surface-layer turbulence measurements, and under unstable conditions based on the aerosol backscatter profile from ceilometer observations.     

A Comparative Study and Evaluation of Mixing-Height Estimation Based on Sodar-RASS, Ceilometer Data and Numerical Model Simulations

The unmanned aerial system SUMO (Small Unmanned Meteorological Observer) has been used for the observation of the structure and behaviour of the atmospheric boundary layer above the Advent Valley, Svalbard during a two-week period in early spring 2009. Temperature, humidity and wind profiles measured by the SUMO system have been compared with measurements of a small tethered balloon system that was operated simultaneously. It is shown that both systems complement each other. Above 200?m, the SUMO system outperforms the tethered balloon in terms of flexibility and the ability to penetrate strong inversion layers of the Arctic boundary layer. Below that level, the tethered balloon system provides atmospheric profiles with higher accuracy, mainly due to its ability to operate at very low vertical velocities. For the observational period, a numerical mesoscale model has been run at high resolution and evaluated with SUMO profiles reaching up to a height of 1500?m above the ground. The sensitivity to the choice of atmospheric boundary-layer schemes and horizontal resolution has been investigated. A new scheme especially suited for stable conditions slightly improves the temperature forecast in stable conditions, although all schemes show a warm bias close to the surface and a cold bias above the atmospheric boundary layer. During one cold and cloudless night, the SUMO system could be operated nearly continuously (every 30?C45?minutes). This allowed for a detailed case study addressing the structure and behaviour of the air column within and above Advent Valley and its interaction with the local topography. The SUMO measurements in conjunction with a 10-m meteorological mast enabled the identification of a very stable nocturnal surface layer adjacent to the valley bottom, a stable air column in the valley and a strong inversion layer above the summit height. The results indicate the presence of inertial-gravity waves during the night, a feature not captured by the model.     

Turbulence Characteristics Of The Stable Boundary Layer Over A Mid-Latitude Glacier. Part Ii: Pure Katabatic Forcing Conditions

Observations obtained over a glacier surface in a predominantlykatabatic flow and with a distinctwind maximum below 13-m height are presented. The data werecollected using a 13-m high profilemast and two sonic anemometers (at about 2.5-m and 10-m heights).The spectra at frequencies belowthat of the turbulence range appear to deviate considerably fromthe curves obtained by Kaimal andco-workers during the 1968 Kansas experiment. The characteristicsof these deviations are compared tothe observations of others in surface-layers disturbed by anykind of large-scale outer-layer (orinactive) turbulence. In our case the disturbances arelikely to be induced by the highmountain ridges that surround the glacier. Moreover, the deviationsobserved in the cospectra seemto result from an, as yet, unspecified interaction between theinactive outer-layer turbulenceand the local surface-layer turbulence. Near the distinctwind maximum turbulence production ceasedwhile turbulence itself did not, probably the result ofturbulence transport from other levels. Consequently, we studied thelocal similarity relations using _w instead of u_* as an alternative velocity scale. Wellbelow the wind maximum, and for relatively low stability(0< Ri_g <0.2), the flow behaves accordingto well established local-scaling similarity relationshipsin the stable boundary layer. For higherstability (Ri_g > 0.2), and near or above the wind maximum, the boundary-layer structure conforms tothat of z-less stratification suggesting that the eddy sizeis restricted by the local stability ofthe flow. In line with this we observed that the sensibleheat fluxes relate remarkably well to thelocal flow parameters.     

On the Turbulence in the Upper Part of the Low-Level Jet: An Experimental and Numerical Study

The characteristics of low-level jets (LLJ) observed at the “Centro de Investigacion de la Baja Atmósfera” (CIBA) site in Spain are analysed, focussing on the turbulence generated in the upper part of the jet, a feature that is still to be thoroughly understood. During the Stable Boundary Layer Experiment in Spain (SABLES) 1998, captive balloon soundings were taken intensively, and their analyses have highlighted the main characteristics of the jet’s wind and temperature structure, leading to a composite profile. There are indications that the turbulence has a minimum at the level of the wind maximum, with elevated turbulence in a layer at a height between two and three times that of the LLJ maximum, but no direct measurements of turbulence were available at these heights. In September 2001, a 100-m tower at the same site was re-instrumented to give turbulence measurements up to 96.6 m above ground level. All occurrences of LLJ below this height between September 2002 and June 2003 have been selected and significant turbulence above the LLJ has been found. Simulations with a single-column turbulence kinetic energy model have been made in order to further investigate the generation of elevated turbulence. The results correlate well with the measurements, showing that in the layer above the LLJ, where there is significant shear and weakly stable stratification, conditions are conducive to the development of turbulence.     

BOUNDARY LAYER GROWTH AND LAPSE RATE CHANGES DETERMINED BY LIDAR AND SURFACE HEAT FLUX IN SOFIA*

In this study the results from a boundary layer experiment,conducted in autumn 1991 over a flat,build-up urban area in Southeast Sofia,together with some models for mixed layer growth rates are used to investigate the layered structure of the vertical atmospheric stability distribution in the Sofia Valley.Lidar measurements of aerosol layer heights and morning boundary layer development are combined with surface eddy correlation measurements of kinematic heat and moisture fluxes,profiles of temperature and humidity,wind speed and wind direction.A diagnostic method is presented for determining vertical lapse rates using surface meteorological measurements and lidar returns observed during the transition from nighttime stable stratification to daytime convective boundary layer after the sunrise.     

BOUNDARY LAYER GROWTH AND LAPSE RATE CHANGES DETERMINED BY LIDAR AND SURFACE HEAT FLUX IN SOFIA

In this study the results from a boundary layer experiment,conducted in autumn 1991 over a flat,build-up urbanarea in Southeast Sofia,together with some models for mixed layer growth rates are used to investigate the layered struc-ture of the vertical atmospheric stability distribution in the Sofia Valley.Lidar measurements of aerosol layer heightsand morning boundary layer development are combined with surface eddy correlation measurements of kinematic heatand moisture fluxes,profiles of temperature and humidity,wind speed and wind direction.A diagnostic method is pres-ented for determining vertical lapse rates using surface meteorological measurements and lidar returns observed duringthe transition from nighttime stable stratification to daytime convective boundary layer after the sunrise.     

Measurements and Modelling of the Wind Speed Profile in the Marine Atmospheric Boundary Layer

We present measurements from 2006 of the marine wind speed profile at a site located 18 km from the west coast of Denmark in the North Sea. Measurements from mast-mounted cup anemometers up to a height of 45 m are extended to 161 m using LiDAR observations. Atmospheric turbulent flux measurements performed in 2004 with a sonic anemometer are compared to a bulk Richardson number formulation of the atmospheric stability. This is used to classify the LiDAR/cup wind speed profiles into atmospheric stability classes. The observations are compared to a simplified model for the wind speed profile that accounts for the effect of the boundary-layer height. For unstable and neutral atmospheric conditions the boundary-layer height could be neglected, whereas for stable conditions it is comparable to the measuring heights and therefore essential to include. It is interesting to note that, although it is derived from a different physical approach, the simplified wind speed profile conforms to the traditional expressions of the surface layer when the effect of the boundary-layer height is neglected.     

Diabatic wind speed profiles in coastal regions: Comparison of an internal boundary layer (IBL) model with observations

A model is presented to transform wind speed observations at a single height over sea or near the coast to any possible location and height in a topographic flat coastal region (up to distances of about 5 km from the coast and up to heights of 100 m). Only moderate and strong winds from the sea are considered, which are particularly important for wind energy applications. The model, called diabatic coast model, which is based on the well known internal boundary layer (IBL) concept and Monin-Obukhov similarity theory, describes the effects of the roughness transition from sea to land as well as the effect of stability on the shape of the profiles and the IBL growth. The predicted IBL heights are compared with published data.In the second part of this paper, the model is compared with measurements taken at the Maasvlakte location near the Dutch coast. It is shown that a neutral formulation of the IBL height is sufficient to model the overall mean wind speed with height, but that stability corrections are needed to describe the diurnal variations in wind speed properly. Finally, an application is given, where a single routine wind speed observation at the coast, combined with air-water temperature differences is used to predict the wind speed at 500m from the coast at heights of 10 and 53 m. The results are in good agreement with the measurements.     

A multi-limit formulation for the equilibrium depth of a stably stratified boundary layer

Currently no expression for the equilibrium depth of the turbulent stably-stratified boundary layer is available that accounts for the combined effects of rotation, surface buoyancy flux and static stability in the free flow. Various expressions proposed to date are reviewed in the light of what is meant by the stable boundary layer. Two major definitions are thoroughly discussed. The first emphasises turbulence and specifies the boundary layer as a continuously and vigorously turbulent layer adjacent to the surface. The second specifies the boundary layer in terms of the mean velocity profile, e.g. by the proximity of the actual velocity to the geostrophic velocity. It is shown that the expressions based on the second definition are relevant to the Ekman layer and portray the depth of the turbulence in the intermediate regimes, when the effects of static stability and rotation essentially interfere. Limiting asymptotic regimes dominated by either stratification or rotation are examined using the energy considerations. As a result, a simple equation for the depth of the equilibrium stable boundary layer is developed. It is valid throughout the range of stability conditions and remains in force in the limits of a perfectly neutral layer subjected to rotation and a rotation-free boundary layer dominated by surface buoyancy flux or stable density stratification at its outer edge. Dimensionless coefficients are estimated using data from observations and large-eddy simulations. Well-known and widely used formulae proposed earlier by Zilitinkevich and by Pollard, Rhines and Thompson are shown to be characteristic of the above interference regimes, when the effects of rotation and static stability (due to either surface buoyancy flux, or stratification at the outer edge of the boundary layer) are roughly equally important.     

Determining Boundary-Layer Height from Aircraft Measurements

The height of the atmospheric boundary layer (ABL) is an important variable in both observational studies and model simulations. The most commonly used measurement for obtaining ABL height is a rawinsonde profile. Mesoscale or regional scale models use a bulk Richardson number based on profiles of the forecast variables. Here we evaluate the limitations of several frequently-used approaches for defining ABL height from a single profile, and identify the optimal threshold value for each method if profiles are the only available measurements. Aircraft measurements from five field projects are used, representing a variety of ABL conditions including stable, convective, and cloud-topped boundary layers over different underlying surfaces. ABL heights detected from these methods were validated against the ‘true’ value determined from aircraft soundings, where ABL height is defined as the top of the layer with significant turbulence. A detection rate was defined to denote how often the ABL height was correctly diagnosed with a particular method. The results suggest that the temperature gradient method provides the most reasonable estimates, although the detection rate and suitable detection criteria vary for different types of ABL. The Richardson number method, on the other hand, is in most cases inadequate or inferior to the other methods that were tried. The optimal range of the detection criteria is given for all ABL types examined in this study.     

Contrasting vertical structures of nocturnal boundary layers

This study analyzes eight levels of sonic anemometerdata collected on a 60-m towerduring CASES-99, toward the goal of understanding thevertical structure of thenocturnal boundary layer. Several different regimesare found. Thin boundarylayers are often observed where fluxes decrease with height and approximately vanish between 20 and 30 m aboveground. The flow above the thin boundary layeraccelerates and increasing shear oftengenerates significant turbulence in the middle ofthe night. Thisshear-generated turbulence is often stronger thanthat near the surface corresponding to an upside-downboundary layer. During these conditions,the turbulent transport of turbulence is downwardtoward the surface. The turbulence in this regimeshows features of z-less turbulence to the extentthat neither the height above groundnor the boundary-layer depth are primary scalingvariables. This layer isdifferent from a `residual layer' in thatturbulence is actively generated byshear associated with nocturnal accelerationsand often is stronger than that inthe surface-based boundary layer.In many cases, the turbulence does not varysignificantly across the towerlayer, implying that the boundary layer ismuch deeper than the 60-m towerlayer. Several case studies are presentedto illustrate the largevariation of vertical structure betweennights.     

基于风廓线雷达的广东登陆台风边界层高度特征研究

针对8个登陆广东省的热带气旋,利用经过数据质量控制的风廓线雷达连续、高时空分辨率的风场观测数据,对热带气旋边界层特征进行了分析。研究结果表明:热带气旋边界层中切向风速大值区垂直范围越大、风速越强、持续时间越久,则热带气旋强度越大、登陆后强度维持时间越久。眼区外入流层厚度越大,入流层气流越强,热带气旋登陆后强度维持时间则越久。风廓线雷达信噪比垂直梯度对大气湍流信息有一定的指示作用,对于入流层高度在2000 m以下的热带气旋,其入流层顶所在高度与信噪比梯度最大值所在高度相近,对于入流层较为深厚的热带气旋,用信噪比垂直梯度确定的边界层高度虽接近入流层顶高,但仍有一定差距。不同特点的热带气旋其边界层高度并不相同,对于登陆后强度迅速减弱的热带气旋边界层高度在500~1000 m;登陆后强度持续时间短的热带气旋,其边界层高度约1000~2000 m;登陆后强度持续时间长的热带气旋,其边界层高度在2000 m之上,最高可达5000~7000 m。这些结果加深了对登陆台风边界层高度演变特征的认识。     

Formation and Characteristics of Coastal Internal Boundary Layers During Onshore Flows

The development and characteristics of coastal internal boundary layers were investigated in 28 tests. These were made at all seasons and in both gradient and sea-breeze flows but only during mid-day periods. Measurements of turbulence and temperature were taken from a light aircraft which flew traverses across Long Island at successive altitudes parallel to the wind direction. These were used to locate the boundary between modified and unmodified air as a function of height and distance from the coast. The same measurements plus tower measurements of wind, turbulence and temperature, pilot balloon soundings and measurements of land and water surface temperatures by a remote sensing IR thermometer were used to quantify the characteristics of the modified and unmodified air. The boundary layer slope was steep close to the land-water interface and became shallower with downwind distance. Growth of the boundary layer was initially slower with stable lapse rates upwind than with neutral or unstable conditions over the water. An equilibrium height was found in many tests except under conditions of free convection when the internal boundary layer merged into the mixed layer inland and with sea-breeze conditions. The equilibrium height depended on downwind conditions and was greater with low wind speeds and strong land surface heating than with stronger winds and small land-water temperature differences. Current theoretical models are not adequate to predict the height of the boundary layer at the altitudes and distances studied but reasonably good predictions were given by an empirical model developed earlier. Wind speed in the modified air averaged about 70% of that at the coast but turbulence levels were several times higher both near the surface and aloft. These findings have important implications for diffusion from coastal sites.     

Formation and characteristics of coastal internal boundary layers during onshore flows

The development and characteristics of coastal internal boundary layers were investigated in 28 tests. These were made at all seasons and in both gradient and sea-breeze flows but only during mid-day periods. Measurements of turbulence and temperature were taken from a light aircraft which flew traverses across Long Island at successive altitudes parallel to the wind direction. These were used to locate the boundary between modified and unmodified air as a function of height and distance from the coast. The same measurements plus tower measurements of wind, turbulence and temperature, pilot balloon soundings and measurements of land and water surface temperatures by a remote sensing IR thermometer were used to quantify the characteristics of the modified and unmodified air. The boundary layer slope was steep close to the land-water interface and became shallower with downwind distance. Growth of the boundary layer was initially slower with stable lapse rates upwind than with neutral or unstable conditions over the water. An equilibrium height was found in many tests except under conditions of free convection when the internal boundary layer merged into the mixed layer inland and with sea-breeze conditions. The equilibrium height depended on downwind conditions and was greater with low wind speeds and strong land surface heating than with stronger winds and small land-water temperature differences. Current theoretical models are not adequate to predict the height of the boundary layer at the altitudes and distances studied but reasonably good predictions were given by an empirical model developed earlier. Wind speed in the modified air averaged about 70% of that at the coast but turbulence levels were several times higher both near the surface and aloft. These findings have important implications for diffusion from coastal sites.     

Boundary-layer effects on mountain waves: a new look at some historical studies

Early studies of mountain waves reported various results that have rarely been investigated since. These include: large-amplitude mountain waves above an unstable boundary layer much higher than the mountains; a repeated downwind drift and upwind jump of mountain waves; and larger vertical wind magnitude near sunrise and/or sunset. These are investigated using over 3,000 radiosondes and meso-strato-troposphere (MST) radar. Superadiabatic temperature gradients are found beneath mountain waves, explainable by convection which appears to raise the mountain-wave launching height. Movement of mountain-wave patterns is studied by a new method using height–time vertical wind data. A swaying motion of mountain waves, with period of a few minutes, appears to be equally upwind and downwind, rather than asymmetric at the heights measurable. Also, vertical wind shows no change in mean, variance or extreme values near sunrise and sunset, despite the expected diurnal changes of boundary-layer structure. An explanation for differences between MST radar and other measurements and models of mountain waves is suggested in terms of more than one variety of mountain wave. Type 1 has stable air near the ground; type 2 is above a convective/turbulent boundary layer of significant height as compared to the mountains.     

Wind profiles,momentum fluxes and roughness lengths at Cabauw revisited

We describe the results of an experiment focusing on wind speed and momentum fluxes in the atmospheric boundary layer up to 200 m. The measurements were conducted in 1996 at the Cabauw site in the Netherlands. Momentum fluxes are measured using the K-Gill Propeller Vane. Estimates of the roughness length are derived using various techniques from the wind speed and flux measurements, and the observed differences are explained by considering the source area of the meteorological parameters. A clear rough-to-smooth transition is found in the wind speed profiles at Cabauw. The internal boundary layer reaches the lowest k-vane (20 m) only in the south-west direction where the obstacle-free fetch is about 2 km. The internal boundary layer is also reflected in the roughness lengths derived from the wind speed profiles. The lower part of the profile (< 40 m) is not in equilibrium and no reliable roughness analysis can be given. The upper part of the profile can be linked to a large-scale roughness length. Roughness lengths derived from the horizontal wind speed variance and gustiness have large footprints and therefore represent a large-scale average roughness. The drag coefficient is more locally determined but still represents a large-scale roughness length when it is measured above the local internal boundary layer. The roughness length at inhomogeneous sites can therefore be determined best from drag coefficient measurements just above the local internal boundary layers directly, or indirectly from horizontal wind speed variance or gustiness. In addition, the momentum and heat fluxes along the tower are analysed and these show significant variation with height related to stability and possibly surface heterogeneity. It appears that the dimensionless wind speed gradients scale well with local fluxes for the variety of conditions considered, including the unstable cases.     

An Auxiliary Tool to Determine the Height of the Boundary Layer

Results from radiosoundings, performed both over land and over sea, show that the ascent rate of a radiosounding balloon, the vertical velocity of the balloon, can be used to determine the height of the boundary layer. In many cases the balloon has a higher ascent rate in the boundary layer and a lower, less variable, ascent rate above. The decrease in ascending velocity appears as a jump at the top of the boundary layer. Two examples of potential temperature profiles for unstable stratification and one profile for stable conditions are shown with the corresponding ascent rates. A comparison between the boundary-layer height determined from potential temperature profiles and from ascent rates is presented for a larger dataset. The different ascent rates of the balloon in the boundary layer and above can be explained by a decrease in drag on the balloon in combination with a lowering of the critical Reynolds number in the boundary layer caused by turbulence. Hence, by simply logging the time from release of a radiosonde, it is possible to obtain additional information that can be used to estimate the height of both the unstable and stable boundary layers.     

Quantifying the local influence at a tall tower site in nocturnal conditions

The influence of the local terrestrial environment on nocturnal atmospheric CO₂ measurements at a 329-m television transmitter tower (and a component of a CO₂ monitoring network) was estimated with a tracer release experiment and a subsequent simulation of the releases. This was done to characterize the vertical transport of emissions from the surface to the uppermost tower level and how it is affected by atmospheric stability. The tracer release experiment was conducted over two nights in May of 2009 near the Department of Energy’s Savannah River Site (SRS) in South Carolina. Tracer was released on two contrasting nights—slightly stable and moderately stable—from several upwind surface locations. Measurements at the 329-m level on both nights indicate that tracer was able to mix vertically within a relatively short (～24 km) distance, implying that nocturnal stable conditions do not necessarily prevent vertical dispersion in the boundary layer and that CO₂ measurements at the tower are at least partly influenced by nearby emissions. A simulation of the tracer release is used to calculate the tower footprint on the two nights to estimate the degree to which the local domain affects the tower readings. The effect of the nocturnal boundary layer on the area sampled by the tower can be seen clearly, as the footprints were affected by changes in stability. The contribution of local sources to the measurements at the tower was minimal, however, suggesting that nocturnal concentrations at upper levels are contributed mostly by regional sources.     

基于高塔观测的登陆台风边界层风切变指数拟合

利用深圳气象梯度观测塔观测数据,以2017年以来进入深圳150 km范围的7个台风个例为研究对象,基于幂指数律拟合讨论台风边界层风切变指数的变化规律。结果表明:幂指数能较好地拟合台风影响下350 m高度以下风廓线,随着拟合高度范围增加,风切变指数增大,拟合精度基本维持;用深圳气象梯度观测塔等差层数据拟合台风风速效果好于全层次数据和等比层数据拟合;7个台风影响期间拟合风切变指数平均值为0.268,明显高于以往研究(0.1～0.177),主要原因是拟合的高度范围较以往研究明显增大,此外还与强风样本较少以及下垫面更粗糙有关。利用幂函数拟合台风不同风速段最大风切变指数,可在台风过程中预估不同高度极端大风风险。研究还表明:台风眼经过铁塔前后风切变指数明显升高,在抗风设计以及台风防御过程中应充分考虑这一变化。