Structure of the planetary boundary layer over southeast England: modeling and measurements

The Weather Research and Forecasting model was applied to analyze variations in the planetary boundary layer (PBL) structure over Southeast England including central and suburban London. The parameterizations and predictive skills of two nonlocal mixing PBL schemes, YSU and ACM2, and two local mixing PBL schemes, MYJ and MYNN2, were evaluated over a variety of stability conditions, with model predictions at a 3 km grid spacing. The PBL height predictions, which are critical for scaling turbulence and diffusion in meteorological and air quality models, show significant intra-scheme variance (> 20%), and the reasons are presented. ACM2 diagnoses the PBL height thermodynamically using the bulk Richardson number method, which leads to a good agreement with the lidar data for both unstable and stable conditions. The modeled vertical profiles in the PBL, such as wind speed, turbulent kinetic energy (TKE), and heat flux, exhibit large spreads across the PBL schemes. The TKE predicted by MYJ were found to be too small and show much less diurnal variation as compared with observations over London. MYNN2 produces better TKE predictions at low levels than MYJ, but its turbulent length scale increases with height in the upper part of the strongly convective PBL, where it should decrease. The local PBL schemes considerably underestimate the entrainment heat fluxes for convective cases. The nonlocal PBL schemes exhibit stronger mixing in the mean wind fields under convective conditions than the local PBL schemes and agree better with large-eddy simulation (LES) studies.

Study of mineral dust entrainment in the planetary boundary layer by lidar depolarisation technique

Measurements on 27 June 2011 were performed over the Southern Iberian Peninsula at Granada EARLINET station, using active and passive remote sensing and airborne and surface in-situ data in order to study the entrainment processes between aerosols in the free troposphere and those in the planetary boundary layer (PBL). To this aim the temporal evolution of the lidar depolarisation, backscatter-related Angström exponent and potential temperature profiles were used in combination with the PBL contribution to the aerosol optical depth (AOD). Our results show that the mineral dust entrainment in the PBL was caused by the convective processes which ‘trapped’ the lofted mineral dust layer, distributing the mineral dust particles within the PBL. The temporal evolution of ground-based in-situ data evidenced the impact of this process at surface level. Finally, the amount of mineral dust in the atmospheric column available to be dispersed into the PBL was estimated by means of POLIPHON (Polarizing Lidar Photometer Networking). The dust mass concentration derived from POLIPHON was compared with the coarse-mode mass concentration retrieved with airborne in-situ measurements. Comparison shows differences below 50 µg/m³ (30% relative difference) indicating a relative good agreement between both techniques.

Boundary layer structure and decoupling from synoptic scale flow during NAMBLEX

This paper presents an overview of the meteorology and planetary boundary layer structure observed during the NAMBLEX field campaign to aid interpretation of the chemical and aerosol measurements. The campaign has been separated into five periods corresponding to the prevailing synoptic condition. Comparisons between meteorological measurements ( UHF wind profiler, Doppler sodar, sonic aneometers mounted on a tower at varying heights and a standard anemometer) and the ECMWF analysis at 10 m and 1100 m identified days when the internal boundary layer was decoupled from the synoptic flow aloft. Generally the agreement was remarkably good apart from during period one and on a few days during period four when the diurnal swing in wind direction implies a sea/land breeze circulation near the surface. During these periods the origin of air sampled at Mace Head would not be accurately represented by back trajectories following the winds resolved in ECMWF analyses. The wind profiler observations give a detailed record of boundary layer structure including an indication of its depth, average wind speed and direction. Turbulence statistics have been used to assess the height to which the developing internal boundary layer, caused by the increased surface drag at the coast, reaches the sampling location under a wide range of marine conditions. Sampling conducted below 10 m will be impacted by emission sources at the shoreline in all wind directions and tidal conditions, whereas sampling above 15 m is unlikely to be affected in any of the wind directions and tidal heights sampled during the experiment.

Implementation of non-local boundary layer schemes in the Regional Atmospheric Modeling System and its impact on simulated mesoscale circulations

This paper proposes the implementation of different non-local Planetary Boundary Layer schemes within the Regional Atmospheric Modeling System (RAMS) model. The two selected PBL parameterizations are the Medium-Range Forecast (MRF) PBL and its updated version, known as the Yonsei University (YSU) PBL. YSU is a first-order scheme that uses non-local eddy diffusivity coefficients to compute turbulent fluxes. It is based on the MRF, and improves it with an explicit treatment of the entrainment. With the aim of evaluating the RAMS results for these PBL parameterizations, a series of numerical simulations have been performed and contrasted with the results obtained using the Mellor and Yamada (MY) scheme, also widely used, and the standard PBL scheme in the RAMS model. The numerical study carried out here is focused on mesoscale circulation events during the summer, as these meteorological situations dominate this season of the year in the Western Mediterranean coast. In addition, the sensitivity of these PBL parameterizations to the initial soil moisture content is also evaluated. The results show a warmer and moister PBL for the YSU scheme compared to both MRF and MY. The model presents as well a tendency to overestimate the observed temperature and to underestimate the observed humidity, considering all PBL schemes and a low initial soil moisture content. In addition, the bias between the model and the observations is significantly reduced moistening the initial soil moisture of the corresponding run. Thus, varying this parameter has a positive effect and improves the simulated results in relation to the observations. However, there is still a significant overestimation of the wind speed over flatter terrain, independently of the PBL scheme and the initial soil moisture used, even though a different degree of accuracy is reproduced by RAMS taking into account the different sensitivity tests.

Numerical prediction of gas concentrations and fluctuations above a triangular hill within a turbulent boundary layer

The objective of the present work is to propose a numerical and statistical approach, using computational fluid dynamics, for the study of the atmospheric pollutant dispersion. Modifications in the standard k-epsilon turbulence model and additional equations for the calculation of the variance of concentration are introduced to enhance the prediction of the flow field and scalar quantities. The flow field, the mean concentration and the variance of a flow over a two-dimensional triangular hill, with a finite-size point pollutant source, are calculated by a finite volume code and compared with published experimental results. A modified low Reynolds k-epsilon turbulence model was employed in this work, using the constant of the k-epsilon model C(mu)=0.03 to take into account the inactive atmospheric turbulence. The numerical results for the velocity profiles and the position of the reattachment point are in good agreement with the experimental results. The results for the mean and the variance of the concentration are also in good agreement with experimental results from the literature. (C) 2009 Elsevier Ltd. All rights reserved.

Boundary-layer transition measurements in hypervelocity flows in a shock tunnel

Experiments to investigate the transition process in hypervelocity boundary layers were performed in the T4 free-piston shock tunnel. An array of thin-film heat-transfer gauges was used to detect the location and extent of the transitional region on a 1500 mm long x 120 turn wide flat plate, which formed one of the walls of a duct. The experiments were performed in a Mach 6 flow of air with 6- and 12-MJ/kg nozzle-supply enthalpies at unit Reynolds numbers ranging from 1.6 x 10(6) to 4.9 x 10(6) m(-1). The results show that the characteristics typical of transition taking place through the initiation, growth, and merger of turbulent spots are evident in the heat-transfer signals. A 2-mm-high excrescence located 440 turn from the leading edge was found to be capable of generating a turbulent wedge within an otherwise laminar boundary layer at a unit Reynolds number of 2.6 x 10(6) m(-1) at the 6-MJ/kg condition. A tripping strip, located 100 mm from the leading edge and consisting of a line 37 teeth of 2 rum height equally spaced and spanning the test surface, was also found to be capable of advancing the transition location at the same condition and at the higher enthalpy condition.

The Lake Tekapo experiment (LTEX): An investigation of atmospheric boundary layer processes in complex terrain

A research program on atmospheric boundary layer processes and local wind regimes in complex terrain was conducted in the vicinity of Lake Tekapo in the southern Alps of New Zealand, during two 1-month field campaigns in 1997 and 1999. The effects of the interaction of thermal and dynamic forcing were of specific interest, with a particular focus on the interaction of thermal forcing of differing scales. The rationale and objectives of the field and modeling program are described, along with the methodology used to achieve them. Specific research aims include improved knowledge of the role of surface forcing associated with varying energy balances across heterogeneous terrain, thermal influences on boundary layer and local wind development, and dynamic influences of the terrain through channeling effects. Data were collected using a network of surface meteorological and energy balance stations, radiosonde and pilot balloon soundings, tethered balloon and kite-based systems, sodar, and an instrumented light aircraft. These data are being used to investigate the energetics of surface heat fluxes, the effects of localized heating/cooling and advective processes on atmospheric boundary layer development, and dynamic channeling. A complementary program of numerical modeling includes application of the Regional Atmospheric Modeling System (RAMS) to case studies characterizing typical boundary layer structures and airflow patterns observed around Lake Tekapo. Some initial results derived from the special observation periods are used to illustrate progress made to date. In spite of the difficulties involved in obtaining good data and undertaking modeling experiments in such complex terrain, initial results show that surface thermal heterogeneity has a significant influence on local atmospheric structure and wind fields in the vicinity of the lake. This influence occurs particularly in the morning. However, dynamic channeling effects and the larger-scale thermal effect of the mountain region frequently override these more local features later in the day.

An approximate solution method for boundary layer flow of a power law fluid over a flat plate

The work in this paper deals with the development of momentum and thermal boundary layers when a power law fluid flows over a flat plate. At the plate we impose either constant temperature, constant flux or a Newton cooling condition. The problem is analysed using similarity solutions, integral momentum and energy equations and an approximation technique which is a form of the Heat Balance Integral Method. The fluid properties are assumed to be independent of temperature, hence the momentum equation uncouples from the thermal problem. We first derive the similarity equations for the velocity and present exact solutions for the case where the power law index n = 2. The similarity solutions are used to validate the new approximation method. This new technique is then applied to the thermal boundary layer, where a similarity solution can only be obtained for the case n = 1.

Effect of different surface roughnesses on a turbulent boundary layer

The classical treatment of rough wall turbulent boundary layers consists in determining the effect the roughness has on the mean velocity profile. This effect is usually described in terms of the roughness function delta U+. The general implication is that different roughness geometries with the same delta U+ will have similar turbulence characteristics, at least at a sufficient distance from the roughness elements. Measurements over two different surface geometries (a mesh roughness and spanwise circular rods regularly spaced in the streamwise direction) with nominally the same delta U+ indicate significant differences in the Reynolds stresses, especially those involving the wall-normal velocity fluctuation, over the outer region. The differences are such that the Reynolds stress anisotropy is smaller over the mesh roughness than the rod roughness. The Reynolds stress anisotropy is largest for a smooth wall. The small-scale anisotropy and interniittency exhibit much smaller differences when the Taylor microscale Reynolds number and the Kolmogorov-normalized mean shear are nominally the same. There is nonetheless evidence that the small-scale structure over the three-dimensional mesh roughness conforms more closely with isotropy than that over the rod-roughened and smooth walls.

An experimental/numerical study of the turbulent boundary layer development along a surface with a sudden change in roughness

A theory for the description of turbulent boundary layer flows over surfaces with a sudden change in roughness is considered. The theory resorts to the concept of displacement in origin to specify a wall function boundary condition for a kappa-epsilon model. An approximate algebraic expression for the displacement in origin is obtained from the experimental data by using the chart method of Perry and Joubert(J.F.M., vol. 17, pp. 193-122, 1963). This expression is subsequently included in the near wall logarithmic velocity profile, which is then adopted as a boundary condition for a kappa-epsilon modelling of the external flow. The results are compared with the lower atmospheric observations made by Bradley(Q. J. Roy. Meteo. Soc., vol. 94, pp. 361-379, 1968) as well as with velocity profiles extracted from a set of wind tunnel experiments carried out by Avelino et al.( 7th ENCIT, 1998). The measurements are found to be in good agreement with the theoretical computations. The skin-friction coefficient was calculated according to the chart method of Perry and Joubert(J.F.M., vol. 17, pp. 193-122, 1963) and to a balance of the integral momentum equation. In particular, the growth of the internal boundary layer thickness obtained from the numerical simulation is compared with predictions of the experimental data calculated by two methods, the "knee" point method and the "merge" point method.

Wind tunnel simulation of atmospheric boundary layer flows

The present work shows how thick boundary layers can be produced in a short wind tunnel with a view to simulate atmospheric flows. Several types of thickening devices are analysed. The experimental assessment of the devices was conducted by considering integral properties of the flow and the spectra: skin-friction, mean velocity profiles in inner and outer co-ordinates and longitudinal turbulence. Designs based on screens, elliptic wedge generators, and cylindrical rod generators are analysed. The paper describes in detail the experimental arrangement, including the features of the wind tunnel and of the instrumentation. The results are compared with experimental data published by other authors and with naturally developed flows.

Studies on the Vertical Structure of Horizontal Wind Variability in the Surface Boundary Layer Over Sriharikota

The present study on the vertical structure of horizontal wind variability in the surface boundary layer over Sriharikota. Based on clock wind speed and direction measuring meteorological tower facility from seven levels in the 100 m layer. The study on wind variability and elliptical approximation of wind hodographs investigated for this tropical coastal station established that Sriharikota is of meso-scale weather entity. Wind variability ratio increases from lower levels to upper levels. In South West monsoon months the station is of high ratio values and it gets affected with meso-scale weather features like thunderstorms. Average total shears are observed greater values than scalar shears. Scalar shears are high in the lowest shear levels compared to upper levels. Semi diurnal types of oscillation in average total shears are found in south west monsoon months. During cyclonic storm passage it is observed that there can be significant difference in mean wind speed from 10 m to 100 m level, but it is not so for peak wind speeds. The variations in wind variability ratio in different months is clearly depicted its strong link to define meso-scale or synoptic –scale forcing domination for this station. Meso-scale forcing is characterized by diurnal wind variability and synoptic- scale forcing by interdiurnal wind variability.