Using a Sodar to Measure Optical Turbulence and Wind Speed for the Thirty Meter Telescope Site Testing. Part II: Comparison with Independent Instruments

In the second part of this study, we compare both the wind speed and turbulence given by the Sodars with independent sets of measurements. In the case of the wind speed we compare the lowest Sodar data bin with a sonic anemometer located on a 7-m tower. The agreement between the two instruments was convincing with a regression slope near unity. The integrated turbulence measurements of the Sodars are compared with those obtained with a combined multi-aperture scintillation sensor and differential image motion monitor (MASS/DIMM) unit. It was found that the Sodars are indeed capable of quantitatively measuring optical turbulence, and agree with the MASS/DIMM measurements with a correlation coefficient of approximately 80% and a regression slope within 10% of unity. Additional acoustic noise in the Sodar data was identified using this comparison and removed from the data.     

Variability and Maintenance of Turbulence in the Very Stable Boundary Layer

The relationship of turbulence quantities to mean flow quantities, such as the Richardson number, degenerates substantially for strong stability, at least in those studies that do not place restrictions on minimum turbulence or non-stationarity. This study examines the large variability of the turbulence for very stable conditions by analyzing four months of turbulence data from a site with short grass. Brief comparisons are made with three additional sites, one over short grass on flat terrain and two with tall vegetation in complex terrain. For very stable conditions, any dependence of the turbulence quantities on the mean wind speed or bulk Richardson number becomes masked by large scatter, as found in some previous studies. The large variability of the turbulence quantities is due to random variations and other physical influences not represented by the bulk Richardson number. There is no critical Richardson number above which the turbulence vanishes. For very stable conditions, the record-averaged vertical velocity variance and the drag coefficient increase with the strength of the submeso motions (wave motions, solitary waves, horizontal modes and numerous more complex signatures). The submeso motions are on time scales of minutes and not normally considered part of the mean flow. The generation of turbulence by such unpredictable motions appears to preclude universal similarity theory for predicting the surface stress for very stable conditions. Large variation of the stress direction with respect to the wind direction for the very stable regime is also examined. Needed additional work is noted.     

Two-dimensional field of thermal turbulence at the edge of an escarpment

Acoustic sounder measurements of the temperature structure parameter were obtained at the edge of an escarpment which is part of a ridge of mountains. These measurements indicate that in mountainous terrain, the daytime two-dimensional field of thermal turbulence is strongly affected by relative sun-slope orientation and wind direction out to ranges of at least 200–300 m. For the geometry of this site, westerly flow results in a field which tends to decrease rapidly to the west in the morning with a much less rapid decrease in the afternoon. At night, easterly flow results in significantly higher thermal turbulence compared to that obtained during westerly flow.These measurements show an increase in thermal turbulence at horizontal ranges of 100–200 m to the west of the escarpment during early afternoon on days with deep mixed layers. It is conjectured that this is due to the mountain upslope wind.     

Characterizing the dispersive state of convective boundary layers for applied dispersion modeling

Estimates from semiempirical models that characterize surface heat flux, mixing depth, and profiles of temperature, wind, and turbulence are compared with observations from atmospheric field studies conducted in Colorado, Illinois, Indiana, and Minnesota. Sodar observations are compared with tower measurements at the Colorado site, for wind and turbulence profiles. The median surface heat flux, as calculated using surface-layer flux-profile relationships and an energy budget model, was consistently overestimated by 20 to 80%. Several mixing-depth models were evaluated: (1) integration of the hourly surface heat flux and friction velocity, (2) solving for the time rate of change of profiles of virtual potential temperature, and (3) an interpolation scheme used by the U.S. Environmental Protection Agency in regulatory dispersion models. For the late afternoon, 80 to 90% of the estimates from the first and third models were within 40% of the observed values. For the morning hours after sunrise, all were less accurate. Temperature estimates from surface-layer flux-profile relationships compared well with observations within the mixed layer, but were too low for the inversion layer aloft. Wind profiles were derived using surface-layer flux-profile relationships, a windprofile power-law based on Pasquill stability category, and sodar measurements. The sodar measurements were superior to both types of model estimates. Turbulence profiles were derived from sodar measurements and from semiempirical similarity relationships based on mixing depth and Obukhov length. The scatter in the comparisons with the sodar observations is twice that seen in the comparisons with empirical profile relationships. Overall, it appears that uncertainty of as low as 20 to 30% in the characterization of the diffusion meteorology is the exception rather than the rule.On assignment from the National Oceanic and Atmospheric Administration, U. S. Department of Commerce.Disclaimer: Although the research described in this article has been supported by the United States Environmental Protection Agency, it has not been subjected to Agency review and therefore does not necessarily reflect the views of the Agency and no official endorsement should be inferred. Mention of trade names or commercial products does not constitute endorsement or recommendation for use.     

An Investigation into Ridge-Top Turbulence Characteristics During Neutral and Weakly Stable Conditions: Velocity Spectra and Isotropy

An investigation into high Reynolds number turbulent flow over a ridge top in New Zealand is described based on high-resolution in-situ measurements, using ultrasonic anemometers for two separate locations on the same ridge with differing upwind terrain complexity. Twelve 5-h periods during neutrally stratified and weakly stable atmospheric conditions with strong wind speeds were sampled at 20 Hz. Large (and small) turbulent length scales were recorded for both vertical and longitudinal velocity components in the range of 7–23 m (0.7–3.3 m) for the vertical direction and 628–1111 m (10.5–14.5 m) for the longitudinal direction. Large-scale eddy sizes scaled to the WRF (Weather Research and Forecasting) numerical model simulated boundary-layer thickness for both sites, while small-scale turbulent features were a function of the complexity of the upwind terrain. Evidence of a multi-scale turbulent structure was obtained at the more complex terrain site, while an assessment of the three-dimensional isotropy assumption in the inertial subrange of the spectrum showed anisotropic turbulence at the less complex site and evidence of isotropic turbulence at the more complex site, with a spectral ratio convergence deviating from the 4/3 or unity values suggested by previous theory and practice. Existing neutral spectral models can represent locations along the ridge top with simple upwind complexity, especially for the vertical wind spectra, but sites with more orographic complexity and strong vertical wind speeds are often poorly represented using these models. Measured spectra for the two sites exhibited no significant diurnal variation and very similar large-scale and small-scale turbulent length scales for each site, but the turbulence energy measured by the variances revealed a strong diurnal difference.     

A Wind Tunnel and Numerical Investigation of Turbulent Dispersion in Coastal Atmospheric Boundary Layers

Thermal internal boundary layers in onshore air flows have a significant influence on pollutant diffusion in coastal areas. Although several models for this diffusion problem exist, measurements for model verification are scarce. In this paper, we present a set of wind tunnel observations and examine the performance of a Lagrangian stochastic model. The good agreement between the model simulation and the tunnel measurements confirms the usefulness of the Lagrangian stochastic model for practical purposes. Sensitivity tests of the model to turbulence statistics show that uncertainty in velocity skewness to the extent of observational scatter does not seem to have a significant influence on pollutant dispersion, while uncertainties in turbulence intensity (variance) significantly influence the dispersion pattern.     

Atmospheric Turbulence Characteristics Over a Temperature-Inhomogeneous Land Surface. Part II: The Effect of Small-Scale Inhomogeneities of Surface Temperature on Some Characteristics of the Atmospheric Surface Layer

Turbulence characteristics, vertical profiles of wind velocity u(z) and air temperature T(z), and also spatial variations in steppe surface radiation temperature T_r(x) are measured simultaneously. A marked effect of T_r(x) characteristics on the turbulence characteristics and T(z) profiles is observed in the lower part of the atmospheric surface layer. We suggest that variability in T_r(x) noticeably influences the surface-layer temperature field and leads to scatter in the values of the universal functions obtained by different authors; effects of T_r(x) are not accounted for in similarity theory. The introduction of the value of temperature zero-plane displacement d_T in the calculation formulae (to determine temperature flux) noticeably improves the agreement between calculated and measured (by eddy-correlation method) results. The influence of footprint (or Source Area) on the obtained results leads to noticeable scatter in the data obtained from measurements of atmospheric turbulence.     

Influence of mesoscale winds on the turbulent structure of the urban boundary layer over St. Louis

Two fair weather afternoons have been examined, where the urban boundary layer over St. Louis, though exhibiting similar thermal characteristics, had a markedly different kinematic structure. The turbulent nature of the boundary layer was examined through analysis of double theodolite wind profiles at an urban and at a rural site on each day. On 14 July 1975, the winds increased with height above the inversion at both sites and on the following day, the winds decreased above the boundary layer in the same region. While the mean wind speed in the lowest 0.8 km agl was similar on both days, the turbulence characteristics of the urban boundary-layer winds were distinctly different on these two afternoons. This was evidenced by the variance of the wind and is in agreement with simultaneous aircraft measurements reported by Hildebrand and Ackerman (1984). A similar difference in turbulence was not found over the rural site. It is suggested that the enhanced turbulence at the urban site on 14 July is likely associated with the wind profile immediately above the boundary layer, where the downward flux of high momentum air from above the inversion may have resulted in stronger mechanical mixing within the boundary layer.     

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.     

Horizontal, Two Point Coherence for Separations Greater Than the Measurement Height

Wind speed measurements from the test site at Rutherford Appleton Laboratory have been evaluated with respect to the spatial coherence function. The experimental arrangement provides coherence information for separation distances of 62, 80 and 102 m. These are at least three times greater than the measurement heights of 18 m and 18.7 m. Based on these experimental data and data published in the literature, different theoretical formulations are compared and a new, but simple, model for longitudinal and lateral coherence is proposed. At large separations the turbulent wind field is not isotropic, theoretical models to describe the coherence function for such distances are not available. The new model we propose builds on the classical exponential approach. It takes into account the influence of turbulence intensity and models the angular dependence of horizontal coherence. It is found that, for constant turbulence intensity, the lateral coherence decay becomes independent of the mean wind speed.     

Effects of Turbulence and Flow Inclination on the Performance of Cup Anemometers in the Field

Four commercial and one research cup anemometers were comparatively tested in a complex terrain site to quantify the effects of turbulence and flow inclination on the wind speed measurements. The difference of the mean windspeed reading between the anemometers was as much as 2% for wind directions where the mean flow was horizontal. This difference was large enough to be attributed to the well-known overspeeding effect related to the differing distance constant (ranging from 1.7 to 5 m) of the cup anemometers. The application of a theoretical model of the cup-anemometer behaviour in athree-dimensional turbulent wind field proved successful in explaining theobserved differences.Additional measurements were taken with the anemometers tilted at known angles into and out of the incident wind flow. Thus, a field-derived angular response curve is constructed for each anemometer and the deviations from publishedwind-tunnel results are discussed.The uncertainties of, or false assumptions about, the angular response characteristics of the anemometers contribute the largest amount inthe observed errors of mean wind speed even for a horizontal mean flow. The angular response curves are finally used to correct the 10-min mean windspeed. The necessary information for the correction is the turbulent intensity (preferably in the vertical direction) and the mean flow inclination.For demanding applications, the angular response parameters of cup anemometers should be taken into account. The incorporation of the angular response parameters in a correction scheme would be most robustly applied if their variation with inclination and wind speed was smooth.     

Heterogeneity-Induced Heat-Flux Patterns in the Convective Boundary Layer: Can they be Detected from Observations and is There a Blending Height?—A Large-Eddy Simulation Study for the LITFASS-2003 Experiment

An understanding of how the convective boundary layer (CBL) is mixed under heterogeneous surface forcing is crucial for the interpretation of area-averaged turbulence measurements. To determine the height and degree to which a complex heterogeneous surface affects the CBL, large-eddy simulations (LES) for two days of the LITFASS-2003 experiment representing two different wind regimes were undertaken. Spatially-lagged correlation analysis revealed the turbulent heat fluxes to be dependent on the prescribed surface flux pattern throughout the entire CBL including the entrainment layer. These findings prompted the question of whether signals induced by surface heterogeneity can be measured by airborne systems. To examine this question, an ensemble of virtual flights was conducted using LES, according to Helipod flight measurements made during LITFASS-2003. The resulting ensemble-averaged heat fluxes indicated a clear dependence on the underlying surface up to the top of the CBL. However, a large scatter between the flux measurements in different ensemble runs was observed, which was the result of insufficient sampling of the largest turbulent eddies. The random and systematic errors based on the integral length scale did not indicate such a large scatter. For the given flight leg lengths, at least 10–15 statistically independent flight measurements were necessary to give a significant estimate of heterogeneity-induced signals in the CBL. The need for ensemble averaging suggests that the observed blending of heterogeneity-induced signals in the CBL can be partly attributed to insufficient averaging.     

Long-Term Profiles of Wind and Weibull Distribution Parameters up to 600 m in a Rural Coastal and an Inland Suburban Area

An investigation of the long-term variability of wind profiles for wind energy applications is presented. The observations consists of wind measurements obtained from a ground-based wind lidar at heights between 100 and 600 m, in combination with measurements from tall meteorological towers at a flat rural coastal site in western Denmark and at an inland suburban area near Hamburg in Germany. Simulations with the weather research and forecasting numerical model were carried out in both forecast and analysis configurations. The scatter between measured and modelled wind speeds expressed by the root-mean-square error was about 10 % lower for the analysis compared to the forecast simulations. At the rural coastal site, the observed mean wind speeds above 60 m were underestimated by both the analysis and forecast model runs. For the inland suburban area, the mean wind speed is overestimated by both types of the simulations below 500 m. When studying the wind-speed variability with the Weibull distribution, the shape parameter was always underestimated by the forecast compared to both analysis simulations and measurements. At the rural coastal site although the measured and modelled Weibull distributions are different their variances are nearly the same. It is suggested to use the shape parameter for climatological mesoscale model evaluation. Based on the new measurements, a parametrization of the shape parameter for practical applications is suggested.     

Analysis and numerical simulation of the meteorological observations at a tropical coastal site Chennai in India

Meteorological measurements were carried out at North Chennai semi rural area during pre-monsoon period as a part of an air quality study program. Analysis of the data showed the effects of coastal terrain namely the land-sea breeze circulation, temperature cooling during the sea breeze, difference in onset times at these sites etc. Sea breeze onset was observed with a sharp turning of the wind from westerly to south easterly associated with rise in wind speed. Advection speed of the front was about 2.0 m s^− 1. A simple mesoscale meteorological model (MAM-I) developed at Kalpakkam for coastal atmospheric dispersion estimation was used to simulate the observed characteristics. All the major features observed could be simulated by the model while significant difference was noticed in sea breeze frontal movement. MAM results were also inter-compared with MM5. There were no significant differences in the estimate of mean parameters by both the models. It is concluded that the simple model, which takes less run time in a desktop PC, is adequate enough for practical application of providing wind field for plume dispersion models at coastal sites.     

étude expérimentale de la couche limite au-dessus d'un relief modéré proche d'une chaîne de montagne

For the heterogeneous site described in the first part, some aspects of the turbulent structure of the planetary boundary layer are studied. Using mixed-layer scaling, the normalized profiles are compared with those obtained over flat terrain during convective conditions. The measurements were made with the same instrumented aircraft at both sites. The dissipative and spectral length scales are smaller over complex terrain within the whole boundary layer. This is due to the shifting of the wavelength peak toward the high frequencies by dynamic turbulence.This last effect can also explain the increase of the dissipation rate over the heterogeneous site during strong wind conditions. The vertical profiles of sensible heat flux and temperature-water vapor correlation show a lack of entrainment process at the top of the boundary layer. This fact suggests that the investigated boundary layer is advected from the neighbouring plain over the complex site (plateau de Lannemezan).     

Internal equilibrium layer growth over forest

Summary  Friction velocity data from different heights above a forest are used to evaluate the influence from the surrounding landscape on forest micro-meteorological measurements under near-neutral conditions. Data are used from one field site and two forest sites. The field site data are used to estimate the magnitude of the scatter. Different theoretical friction velocity profiles for the Internal Boundary Layer (IBL) are tested against the forest data. The results yield information on the Internal Equilibrium Layer (IEL) growth and an equation for the IEL height for neutral conditions is derived. For stable conditions the results indicate that very long fetches are required in order to measure parameters in equilibrium with the actual surface. Received November 22, 1999 Revised February 22, 2000     

A Comparative Study of Stationary and Non-stationary Wind Models Using Field Measurements

We present a comparative study of the conventional stationary wind speed model and a newly proposed non-stationary wind speed model using field measurements. The concept of, and the differences between, the two wind models are briefly reviewed. Wind data recorded by a field measurement system for wind turbulence parameters (FMS-WTP) of 1-year duration are analyzed using the two wind models. Comparisons were made between the wind characteristics obtained from the two models, including hourly mean wind speed, turbulence intensity, the wind spectrum, integral length scale, root coherence function and probability density function. The effects of wind types (monsoon or typhoon), statistical properties (stationary or non-stationary), and surface roughness (open-sea fetch or overland fetch) on wind characteristics are discussed. The comparative study demonstrates that the non-stationary wind model appears to be more appropriate than the conventional stationary wind speed model for characterizing turbulent winds of one-hour duration over complex terrain.     

Impact of waves on air-sea exchange of sensible heat and momentum

The impact of sea waves on sensible heat and momentum fluxes is described. The approach is based on the conservation of heat and momentum in the marine atmospheric surface layer. The experimental fact that the drag coefficient above the sea increases considerably with increasing wind speed, while the exchange coefficient for sensible heat (Stanton number) remains virtually independent of wind speed, is explained by a different balance of the turbulent and the wave-induced parts in the total fluxes of momentum and sensible heat.Organised motions induced by waves support the wave-induced stress which dominates the surface momentum flux. These organised motions do not contribute to the vertical flux of heat. The heat flux above waves is determined, in part, by the influence of waves upon the turbulence diffusivity.The turbulence diffusivity is altered by waves in an indirect way. The wave-induced stress dominates the surface flux and decays rapidly with height. Therefore the turbulent stress above waves is no longer constant with height. That changes the balance of the turbulent kinetic energy and of the dissipation rate and, hence the diffusivity.The dependence of the exchange coefficient for heat on wind speed is usually parameterized in terms of a constant Stanton number. However, an increase of the exchange coefficient with wind speed is not ruled out by field measurements and could be parametrized in terms of a constant temperature roughness length. Because of the large scatter, field data do not allow us to establish the actual dependence. The exchange coefficient for sensible heat, calculated from the model, is virtually independent of wind speed in the range of 3–10 ms^-1. For wind speeds above 10 ms^-1 an increase of 10% is obtained, which is smaller than that following from the constant roughness length parameterization.The investigation was in part supported by the Netherlands Geosciences Foundation (GOA) with financial aid from the Netherlands Organization for Scientific Research (NWO).     

Foehn signals detected by sodar wind and turbulence measurements in the Rhine Valley, Austria, during the MAP field phase

Summary Within the Mesoscale Alpine Programme MAP conducted in autumn 1999, the vertical structure and the temporal evolution of the planetary boundary layer (PBL) in the Rhine Valley 2km south of Lake Constance were observed with a Remtech PA2 sodar (sound-detection-and-ranging instrument) rendering half-hour averages of the three-dimensional wind profile within the lowest kilometre above ground. During Foehn events, tethered balloon soundings and wind profiler measurements were conducted in addition to the rawinsonde network which was built up for the MAP field campaign.The remote sensing instrument renders a surprisingly high number of valid data during south Foehn. Due to the frequent formation of a cold air pool with stable conditions below the Foehn flow with near-neutral static stability, even more sodar data is valid during Foehn periods than during no Foehn periods. A significant reduction of the sodar data quality is only observed during Foehn events with grounding of the Foehn at the sodar site due to high background noise. At higher levels, a Foehn signal can be detected from the sodar wind and turbulence intensitiy information. With Foehn, higher wind speeds and larger turbulence intensities occur than without Foehn. Comparisons to rawinsonde and tethersonde soundings and wind profiler measurements at sites nearby reveal the spatial inhomogeneity of the Foehn flow within this part of the valley as well as instrumental short-comings. Different methods to determine the mixing height using the vertical sounding devices lead to some uncertainty of mixing height estimates which however can reasonably be explained.